Reclassification of a likely pathogenic Dutch founder variant in KCNH2; implications of reduced penetrance.

BACKGROUND: Variants in KCNH2, encoding the human ether a-go-go (hERG) channel that is responsible for the rapid component of the cardiac delayed rectifier K+ current (IKr), are causal to long QT syndrome type 2 (LQTS2). We identified eight index patients with a new variant of unknown significance (VUS), KCNH2:c.2717C > T:p.(Ser906Leu). We aimed to elucidate the biophysiological effect of this variant, to enable reclassification and consequent clinical decision-making. METHODS: A genotype-phenotype overview of the patients and relatives was created. The biophysiological effects were assessed independently by manual-, and automated calibrated patch clamp. HEK293a cells expressing (i) wild-type (WT) KCNH2, (ii) KCNH2-p.S906L alone (homozygous, Hm) or (iii) KCNH2-p.S906L in combination with WT (1:1) (heterozygous, Hz) were used for manual patching. Automated patch clamp measured the variants function against known benign and pathogenic variants, using Flp-In T-rex HEK293 KCNH2-variant cell lines. RESULTS: Incomplete penetrance of LQTS2 in KCNH2:p.(Ser906Leu) carriers was observed. In addition, some patients were heterozygous for other VUSs in CACNA1C, PKP2, RYR2 or AKAP9. The phenotype of carriers of KCNH2:p.(Ser906Leu) ranged from asymptomatic to life-threatening arrhythmic events. Manual patch clamp showed a reduced current density by 69.8 and 60.4% in KCNH2-p.S906L-Hm and KCNH2-p.S906L-Hz, respectively. The time constant of activation was significantly increased with 80.1% in KCNH2-p.S906L-Hm compared with KCNH2-WT. Assessment of KCNH2-p.S906L-Hz by calibrated automatic patch clamp assay showed a reduction in current density by 35.6%. CONCLUSION: The reduced current density in the KCNH2-p.S906L-Hz indicates a moderate loss-of-function. Combined with the reduced penetrance and variable phenotype, we conclude that KCNH2:p.(Ser906Leu) is a low penetrant likely pathogenic variant for LQTS2.

A massively parallel assay accurately discriminates between functionally normal and abnormal variants in a hotspot domain of KCNH2.

Many genes, including KCNH2, contain "hotspot" domains associated with a high density of variants associated with disease. This has led to the suggestion that variant location can be used as evidence supporting classification of clinical variants. However, it is not known what proportion of all potential variants in hotspot domains cause loss of function. Here, we have used a massively parallel trafficking assay to characterize all single-nucleotide variants in exon 2 of KCNH2, a known hotspot for variants that cause long QT syndrome type 2 and an increased risk of sudden cardiac death. Forty-two percent of KCNH2 exon 2 variants caused at least 50% reduction in protein trafficking, and 65% of these trafficking-defective variants exerted a dominant-negative effect when co-expressed with a WT KCNH2 allele as assessed using a calibrated patch-clamp electrophysiology assay. The massively parallel trafficking assay was more accurate (AUC of 0.94) than bioinformatic prediction tools (REVEL and CardioBoost, AUC of 0.81) in discriminating between functionally normal and abnormal variants. Interestingly, over half of variants in exon 2 were found to be functionally normal, suggesting a nuanced interpretation of variants in this "hotspot" domain is necessary. Our massively parallel trafficking assay can provide this information prospectively.

A calibrated functional patch-clamp assay to enhance clinical variant interpretation in KCNH2-related long QT syndrome.

Modern sequencing technologies have revolutionized our detection of gene variants. However, in most genes, including KCNH2, the majority of missense variants are currently classified as variants of uncertain significance (VUSs). The aim of this study was to investigate the utility of an automated patch-clamp assay for aiding clinical variant classification in KCNH2. The assay was designed according to recommendations proposed by the Clinical Genome Sequence Variant Interpretation Working Group. Thirty-one variants (17 pathogenic/likely pathogenic, 14 benign/likely benign) were classified internally as variant controls. They were heterozygously expressed in Flp-In HEK293 cells for assessing the effects of variants on current density and channel gating in order to determine the sensitivity and specificity of the assay. All 17 pathogenic variant controls had reduced current density, and 13 of 14 benign variant controls had normal current density, which enabled determination of normal and abnormal ranges for applying evidence of moderate or supporting strength for VUS reclassification. Inclusion of functional assay evidence enabled us to reclassify 6 out of 44 KCNH2 VUSs as likely pathogenic. The high-throughput patch-clamp assay can provide moderate-strength evidence for clinical interpretation of clinical KCNH2 variants and demonstrates the value of developing automated patch-clamp assays for functional characterization of ion channel gene variants.

Domain reorientation and rotation of an intracellular assembly regulate conduction in Kir potassium channels.

Potassium channels embedded in cell membranes employ gates to regulate K+ current. While a specific constriction in the permeation pathway has historically been implicated in gating, recent reports suggest that the signature ion selectivity filter located in the outer membrane leaflet may be equally important. Inwardly rectifying K+ channels also control the directionality of flow, using intracellular polyamines to stem ion efflux by a valve-like action. This study presents crystallographic evidence of interdependent gates in the conduction pathway and reveals the mechanism of polyamine block. Reorientation of the intracellular domains, concomitant with activation, instigates polyamine release from intracellular binding sites to block the permeation pathway. Conformational adjustments of the slide helices, achieved by rotation of the cytoplasmic assembly relative to the pore, are directly correlated to the ion configuration in the selectivity filter. Ion redistribution occurs irrespective of the constriction, suggesting a more expansive role of the selectivity filter in gating than previously appreciated.

Altered left atrial metrics in patients with cryptogenic stroke: A systematic review and meta-analysis.

BACKGROUND: There is no defined cause for cryptogenic stroke/embolic stroke of undetermined source (CS-ESUS). As atrial fibrillation (AF) develops in a significant proportion of these patients, it has been suggested that left atrial (LA) myopathy may predispose to CS-ESUS. We investigated alterations in echocardiographic measures of LA size and function in patients with CS-ESUS. METHODS: A systematic literature review and meta-analysis was performed. PubMed, EMBASE, Cochrane Library, Web of Science and SCOPUS were searched for articles published between 1 January 1990 and 10 February 2023. All observational studies of adult CS-ESUS patients with LA volume or function measurements performed by transthoracic echocardiogram were included. Individual random effects meta-analyses were performed on LA measurements in the CS-ESUS patients using subgroup analysis of comparator groups. RESULTS: We included 29 articles with 3927 CS-ESUS patients. Analysis of weighted mean differences showed CS-ESUS patients had altered LA structure and function parameters, with a larger maximum indexed LA volume, reduced LA emptying fraction and/or LA reservoir strain, compared to healthy controls and noncardioembolic stroke patients. Conversely, CS-ESUS patients had a smaller left atrium with better function, compared to cardioembolic stroke patients and CS-ESUS patients who subsequently developed atrial fibrillation. CONCLUSIONS: LA volume and function are altered in CS-ESUS patients compared to healthy controls and other stroke aetiologies. An underlying atrial myopathy in a subset of CS-ESUS patients may be involved in both thrombogenesis and dysrhythmia (specifically AF). While LA functional assessment is not currently recommended following stroke, it may offer an opportunity for recurrent stroke risk stratification.

The diagnostic role of T wave morphology biomarkers in congenital and acquired long QT syndrome: A systematic review.

INTRODUCTION: QTc prolongation is key in diagnosing long QT syndrome (LQTS), however 25%-50% with congenital LQTS (cLQTS) demonstrate a normal resting QTc. T wave morphology (TWM) can distinguish cLQTS subtypes but its role in acquired LQTS (aLQTS) is unclear. METHODS: Electronic databases were searched using the terms "LQTS," "long QT syndrome," "QTc prolongation," "prolonged QT," and "T wave," "T wave morphology," "T wave pattern," "T wave biomarkers." Whole text articles assessing TWM, independent of QTc, were included. RESULTS: Seventeen studies met criteria. TWM measurements included T-wave amplitude, duration, magnitude, Tpeak-Tend, QTpeak, left and right slope, center of gravity (COG), sigmoidal and polynomial classifiers, repolarizing integral, morphology combination score (MCS) and principal component analysis (PCA); and vectorcardiographic biomarkers. cLQTS were distinguished from controls by sigmoidal and polynomial classifiers, MCS, QTpeak, Tpeak-Tend, left slope; and COG x axis. MCS detected aLQTS more significantly than QTc. Flatness, asymmetry and notching, J-Tpeak; and Tpeak-Tend correlated with QTc in aLQTS. Multichannel block in aLQTS was identified by early repolarization (ERD30% ) and late repolarization (LRD30% ), with ERD reflecting hERG-specific blockade. Cardiac events were predicted in cLQTS by T wave flatness, notching, and inversion in leads II and V5 , left slope in lead V6 ; and COG last 25% in lead I. T wave right slope in lead I and T-roundness achieved this in aLQTS. CONCLUSION: Numerous TWM biomarkers which supplement QTc assessment were identified. Their diagnostic capabilities include differentiation of genotypes, identification of concealed LQTS, differentiating aLQTS from cLQTS; and determining multichannel versus hERG channel blockade.

Arrhythmogenic effects of ultra-long and bistable cardiac action potentials.

Contemporary accounts of the initiation of cardiac arrhythmias typically rely on after-depolarizations as the trigger for reentrant activity. The after-depolarizations are usually triggered by calcium entry or spontaneous release within the cells of the myocardium or the conduction system. Here we propose an alternative mechanism whereby arrhythmias are triggered autonomously by cardiac cells that fail to repolarize after a normal heartbeat. We investigated the proposal by representing the heart as an excitable medium of FitzHugh-Nagumo cells where a proportion of cells were capable of remaining depolarized indefinitely. As such, those cells exhibit bistable membrane dynamics. We found that heterogeneous media can tolerate a surprisingly large number of bistable cells and still support normal rhythmic activity. Yet there is a critical limit beyond which the medium is persistently arrhythmogenic. Numerical analysis revealed that the critical threshold for arrhythmogenesis depends on both the strength of the coupling between cells and the extent to which the abnormal cells resist repolarization. Moreover, arrhythmogenesis was found to emerge preferentially at tissue boundaries where cells naturally have fewer neighbors to influence their behavior. These findings may explain why atrial fibrillation typically originates from tissue boundaries such as the cuff of the pulmonary vein.

Don't Turn Off the Tap! The Importance of Discovery Science to the Australian Cardiovascular Sector and Improving Clinical Outcomes Into the Future.

Despite significant advances in interventional and therapeutic approaches, cardiovascular disease (CVD) remains the leading cause of death and mortality. To lower this health burden, cardiovascular discovery scientists need to play an integral part in the solution. Successful clinical translation is achieved when built upon a strong foundational understanding of the disease mechanisms involved. Changes in the Australian funding landscape, to place greater emphasis on translation, however, have increased job insecurity for discovery science researchers and especially early-mid career researchers. To highlight the importance of discovery science in cardiovascular research, this review compiles six science stories in which fundamental discoveries, often involving Australian researchers, has led to or is advancing to clinical translation. These stories demonstrate the importance of the role of discovery scientists and the need for their work to be prioritised now and in the future. Australia needs to keep discovery scientists supported and fully engaged within the broader cardiovascular research ecosystem so they can help realise the next game-changing therapy or diagnostic approach that diminishes the burden of CVD on society.

Protocol-Dependent Differences in IC50 Values Measured in Human Ether-Á-Go-Go-Related Gene Assays Occur in a Predictable Way and Can Be Used to Quantify State Preference of Drug Binding.

Current guidelines around preclinical screening for drug-induced arrhythmias require the measurement of the potency of block of voltage-gated potassium channel subtype 11.1 (Kv11.1) as a surrogate for risk. A shortcoming of this approach is that the measured IC50 of Kv11.1 block varies widely depending on the voltage protocol used in electrophysiological assays. In this study, we aimed to investigate the factors that contribute to these differences and to identify whether it is possible to make predictions about protocol-dependent block that might facilitate the comparison of potencies measured using different assays. Our data demonstrate that state preferential binding, together with drug-binding kinetics and trapping, is an important determinant of the protocol dependence of Kv11.1 block. We show for the first time that differences in IC50 measured between protocols occurs in a predictable way, such that machine-learning algorithms trained using a selection of simple voltage protocols can indeed predict protocol-dependent potency. Furthermore, we also show that the preference of a drug for binding to the open versus the inactivated state of Kv11.1 can also be inferred from differences in IC50 values measured between protocols. Our work therefore identifies how state preferential drug binding is a major determinant of the protocol dependence of IC50 values measured in preclinical Kv11.1 assays. It also provides a novel method for quantifying the state dependence of Kv11.1 drug binding that will facilitate the development of more complete models of drug binding to Kv11.1 and improve our understanding of proarrhythmic risk associated with compounds that block Kv11.1.

hERG K(+) channels: structure, function, and clinical significance.

The human ether-a-go-go related gene (hERG) encodes the pore-forming subunit of the rapid component of the delayed rectifier K(+) channel, Kv11.1, which are expressed in the heart, various brain regions, smooth muscle cells, endocrine cells, and a wide range of tumor cell lines. However, it is the role that Kv11.1 channels play in the heart that has been best characterized, for two main reasons. First, it is the gene product involved in chromosome 7-associated long QT syndrome (LQTS), an inherited disorder associated with a markedly increased risk of ventricular arrhythmias and sudden cardiac death. Second, blockade of Kv11.1, by a wide range of prescription medications, causes drug-induced QT prolongation with an increase in risk of sudden cardiac arrest. In the first part of this review, the properties of Kv11.1 channels, including biogenesis, trafficking, gating, and pharmacology are discussed, while the second part focuses on the pathophysiology of Kv11.1 channels.

Molecular Docking Guided Grid-Independent Descriptor Analysis to Probe the Impact of Water Molecules on Conformational Changes of hERG Inhibitors in Drug Trapping Phenomenon.

Human ether a-go-go related gene (hERG) or KV11.1 potassium channels mediate the rapid delayed rectifier current (IKr) in cardiac myocytes. Drug-induced inhibition of hERG channels has been implicated in the development of acquired long QT syndrome type (aLQTS) and fatal arrhythmias. Several marketed drugs have been withdrawn for this reason. Therefore, there is considerable interest in developing better tests for predicting drugs which can block the hERG channel. The drug-binding pocket in hERG channels, which lies below the selectivity filter, normally contains K+ ions and water molecules. In this study, we test the hypothesis that these water molecules impact drug binding to hERG. We developed 3D QSAR models based on alignment independent descriptors (GRIND) using docked ligands in open and closed conformations of hERG in the presence (solvated) and absence (non-solvated) of water molecules. The ligand-protein interaction fingerprints (PLIF) scheme was used to summarize and compare the interactions. All models delineated similar 3D hERG binding features, however, small deviations of about ~0.4 Å were observed between important hotspots of molecular interaction fields (MIFs) between solvated and non-solvated hERG models. These small changes in conformations do not affect the performance and predictive power of the model to any significant extent. The model that exhibits the best statistical values was attained with a cryo_EM structure of the hERG channel in open state without water. This model also showed the best R2 of 0.58 and 0.51 for the internal and external validation test sets respectively. Our results suggest that the inclusion of water molecules during the docking process has little effect on conformations and this conformational change does not impact the predictive ability of the 3D QSAR models.

The Temperature Dependence of Kinetics Associated with Drug Block of hERG Channels Is Compound-Specific and an Important Factor for Proarrhythmic Risk Prediction.

Current mandated preclinical tests for drug-induced proarrhythmia are very sensitive, but not sufficiently specific. This has led to concern that there is a high attrition rate of potentially safe drugs that could have been beneficial to patients. The comprehensive in vitro proarrhythmia initiative has proposed new metrics based around in silico risk predictions, which are informed, among other things, by measures of human ether-à-go-go-related gene channel (hERG) block kinetics. However, high-throughput patch-clamp systems set to collect these data largely operate at ambient temperature, whereas the simulations for risk prediction are carried out at physiologic temperature. The aims of this study were to: 1) determine to what degree kinetics of drug block of hERG are temperature-dependent, 2) assess the impact of any temperature dependence of drug binding kinetics on repolarization in silico, and 3) identify whether a common set of Q10 scalars can be used to extrapolate kinetic data gathered at ambient to physiologic temperatures for use in in silico proarrhythmic risk prediction. We show that, for a range of drugs, kinetics of block are temperature-dependent and, furthermore, that the degree of temperature dependence is different for each drug. As a result, no common set of Q10 scalars could describe the observed range of temperature dependencies. These results suggest that if accurate physiologic temperature models of the kinetics of drug binding are important for in silico risk prediction, the in vitro data should be acquired at physiologic temperature.

Never at rest: insights into the conformational dynamics of ion channels from cryo-electron microscopy.

The tightly regulated opening and closure of ion channels underlies the electrical signals that are vital for a wide range of physiological processes. Two decades ago the first atomic level view of ion channel structures led to a detailed understanding of ion selectivity and conduction. In recent years, spectacular developments in the field of cryo-electron microscopy have resulted in cryo-EM superseding crystallography as the technique of choice for determining near-atomic resolution structures of ion channels. Here, we will review the recent developments in cryo-EM and its specific application to the study of ion channel gating. We will highlight the advantages and disadvantages of the current technology and where the field is likely to head in the next few years.

Sinusoidal voltage protocols for rapid characterisation of ion channel kinetics.

KEY POINTS: Ion current kinetics are commonly represented by current-voltage relationships, time constant-voltage relationships and subsequently mathematical models fitted to these. These experiments take substantial time, which means they are rarely performed in the same cell. Rather than traditional square-wave voltage clamps, we fitted a model to the current evoked by a novel sum-of-sinusoids voltage clamp that was only 8 s long. Short protocols that can be performed multiple times within a single cell will offer many new opportunities to measure how ion current kinetics are affected by changing conditions. The new model predicts the current under traditional square-wave protocols well, with better predictions of underlying currents than literature models. The current under a novel physiologically relevant series of action potential clamps is predicted extremely well. The short sinusoidal protocols allow a model to be fully fitted to individual cells, allowing us to examine cell-cell variability in current kinetics for the first time. ABSTRACT: Understanding the roles of ion currents is crucial to predict the action of pharmaceuticals and mutations in different scenarios, and thereby to guide clinical interventions in the heart, brain and other electrophysiological systems. Our ability to predict how ion currents contribute to cellular electrophysiology is in turn critically dependent on our characterisation of ion channel kinetics - the voltage-dependent rates of transition between open, closed and inactivated channel states. We present a new method for rapidly exploring and characterising ion channel kinetics, applying it to the hERG potassium channel as an example, with the aim of generating a quantitatively predictive representation of the ion current. We fitted a mathematical model to currents evoked by a novel 8 second sinusoidal voltage clamp in CHO cells overexpressing hERG1a. The model was then used to predict over 5 minutes of recordings in the same cell in response to further protocols: a series of traditional square step voltage clamps, and also a novel voltage clamp comprising a collection of physiologically relevant action potentials. We demonstrate that we can make predictive cell-specific models that outperform the use of averaged data from a number of different cells, and thereby examine which changes in gating are responsible for cell-cell variability in current kinetics. Our technique allows rapid collection of consistent and high quality data, from single cells, and produces more predictive mathematical ion channel models than traditional approaches.

The S1 helix critically regulates the finely tuned gating of Kv11.1 channels.

Congenital mutations in the cardiac Kv11.1 channel can cause long QT syndrome type 2 (LQTS2), a heart rhythm disorder associated with sudden cardiac death. Mutations act either by reducing protein expression at the membrane and/or by perturbing the intricate gating properties of Kv11.1 channels. A number of clinical LQTS2-associated mutations have been reported in the first transmembrane segment (S1) of Kv11.1 channels, but the role of this region of the channel is largely unexplored. In part, this is due to problems defining the extent of the S1 helix, as a consequence of its low sequence homology with other Kv family members. Here, we used NMR spectroscopy and electrophysiological characterization to show that the S1 of Kv11.1 channels extends seven helical turns, from Pro-405 to Phe-431, and is flanked by unstructured loops. Functional analysis suggests that pre-S1 loop residues His-402 and Tyr-403 play an important role in regulating the kinetics and voltage dependence of channel activation and deactivation. Multiple residues within the S1 helix also play an important role in fine-tuning the voltage dependence of activation, regulating slow deactivation, and modulating C-type inactivation of Kv11.1 channels. Analyses of LQTS2-associated mutations in the pre-S1 loop or S1 helix of Kv11.1 channels demonstrate perturbations to both protein expression and most gating transitions. Thus, S1 region mutations would reduce both the action potential repolarizing current passed by Kv11.1 channels in cardiac myocytes, as well as the current passed in response to premature depolarizations that normally helps protect against the formation of ectopic beats.

Potassium currents in the heart: functional roles in repolarization, arrhythmia and therapeutics.

This is the second of the two White Papers from the fourth UC Davis Cardiovascular Symposium Systems Approach to Understanding Cardiac Excitation-Contraction Coupling and Arrhythmias (3-4 March 2016), a biennial event that brings together leading experts in different fields of cardiovascular research. The theme of the 2016 symposium was 'K+ channels and regulation', and the objectives of the conference were severalfold: (1) to identify current knowledge gaps; (2) to understand what may go wrong in the diseased heart and why; (3) to identify possible novel therapeutic targets; and (4) to further the development of systems biology approaches to decipher the molecular mechanisms and treatment of cardiac arrhythmias. The sessions of the Symposium focusing on the functional roles of the cardiac K+ channel in health and disease, as well as K+ channels as therapeutic targets, were contributed by Ye Chen-Izu, Gideon Koren, James Weiss, David Paterson, David Christini, Dobromir Dobrev, Jordi Heijman, Thomas O'Hara, Crystal Ripplinger, Zhilin Qu, Jamie Vandenberg, Colleen Clancy, Isabelle Deschenes, Leighton Izu, Tamas Banyasz, Andras Varro, Heike Wulff, Eleonora Grandi, Michael Sanguinetti, Donald Bers, Jeanne Nerbonne and Nipavan Chiamvimonvat as speakers and panel discussants. This article summarizes state-of-the-art knowledge and controversies on the functional roles of cardiac K+ channels in normal and diseased heart. We endeavour to integrate current knowledge at multiple scales, from the single cell to the whole organ levels, and from both experimental and computational studies.

Potassium channels in the heart: structure, function and regulation.

This paper is the outcome of the fourth UC Davis Systems Approach to Understanding Cardiac Excitation-Contraction Coupling and Arrhythmias Symposium, a biannual event that aims to bring together leading experts in subfields of cardiovascular biomedicine to focus on topics of importance to the field. The theme of the 2016 symposium was 'K+ Channels and Regulation'. Experts in the field contributed their experimental and mathematical modelling perspectives and discussed emerging questions, controversies and challenges on the topic of cardiac K+ channels. This paper summarizes the topics of formal presentations and informal discussions from the symposium on the structural basis of voltage-gated K+ channel function, as well as the mechanisms involved in regulation of K+ channel gating, expression and membrane localization. Given the critical role for K+ channels in determining the rate of cardiac repolarization, it is hardly surprising that essentially every aspect of K+ channel function is exquisitely regulated in cardiac myocytes. This regulation is complex and highly interrelated to other aspects of myocardial function. K+ channel regulatory mechanisms alter, and are altered by, physiological challenges, pathophysiological conditions, and pharmacological agents. An accompanying paper focuses on the integrative role of K+ channels in cardiac electrophysiology, i.e. how K+ currents shape the cardiac action potential, and how their dysfunction can lead to arrhythmias, and discusses K+ channel-based therapeutics. A fundamental understanding of K+ channel regulatory mechanisms and disease processes is fundamental to reveal new targets for human therapy.

Tyrosine Residues from the S4-S5 Linker of Kv11.1 Channels Are Critical for Slow Deactivation.

Slow deactivation of Kv11.1 channels is critical for its function in the heart. The S4-S5 linker, which joins the voltage sensor and pore domains, plays a critical role in this slow deactivation gating. Here, we use NMR spectroscopy to identify the membrane-bound surface of the S4S5 linker, and we show that two highly conserved tyrosine residues within the KCNH subfamily of channels are membrane-associated. Site-directed mutagenesis and electrophysiological analysis indicates that Tyr-542 interacts with both the pore domain and voltage sensor residues to stabilize activated conformations of the channel, whereas Tyr-545 contributes to the slow kinetics of deactivation by primarily stabilizing the transition state between the activated and closed states. Thus, the two tyrosine residues in the Kv11.1 S4S5 linker play critical but distinct roles in the slow deactivation phenotype, which is a hallmark of Kv11.1 channels.

Convergence of models of human ventricular myocyte electrophysiology after global optimization to recapitulate clinical long QT phenotypes.

In-silico models of human cardiac electrophysiology are now being considered for prediction of cardiotoxicity as part of the preclinical assessment phase of all new drugs. We ask the question whether any of the available models are actually fit for this purpose. We tested three models of the human ventricular action potential, the O'hara-Rudy (ORD11), the Grandi-Bers (GB10) and the Ten Tusscher (TT06) models. We extracted clinical QT data for LQTS1 and LQTS2 patients with nonsense mutations that would be predicted to cause 50% loss of function in IKs and IKr respectively. We also obtained clinical QT data for LQTS3 patients. We then used a global optimization approach to improve the existing in silico models so that they reproduced all three clinical data sets more closely. We also examined the effects of adrenergic stimulation in the different LQTS subsets. All models, in their original form, produce markedly different and unrealistic predictions of QT prolongation for LQTS1, 2 and 3. After global optimization of the maximum conductances for membrane channels, all models have similar current densities during the action potential, despite differences in kinetic properties of the channels in the different models, and more closely reproduce the prolongation of repolarization seen in all LQTS subtypes. In-silico models of cardiac electrophysiology have the potential to be tremendously useful in complementing traditional preclinical drug testing studies. However, our results demonstrate they should be carefully validated and optimized to clinical data before they can be used for this purpose.

The two-pore domain potassium channel, TWIK-1, has a role in the regulation of heart rate and atrial size.

The two-pore domain potassium (K(+)) channel TWIK-1 (or K2P1.1) contributes to background K(+) conductance in diverse cell types. TWIK-1, encoded by the KCNK1 gene, is present in the human heart with robust expression in the atria, however its physiological significance is unknown. To evaluate the cardiac effects of TWIK-1 deficiency, we studied zebrafish embryos after knockdown of the two KCNK1 orthologues, kcnk1a and kcnk1b. Knockdown of kcnk1a or kcnk1b individually caused bradycardia and atrial dilation (p<0.001 vs. controls), while ventricular stroke volume was preserved. Combined knockdown of both kcnk1a and kcnk1b resulted in a more severe phenotype, which was partially reversed by co-injection of wild-type human KCNK1 mRNA, but not by a dominant negative variant of human KCNK1 mRNA. To determine whether genetic variants in KCNK1 might cause atrial fibrillation (AF), we sequenced protein-coding regions in two independent cohorts of patients (373 subjects) and identified three non-synonymous variants, p.R171H, p.I198M and p.G236S, that were all located in highly conserved amino acid residues. In transfected mammalian cells, zebrafish and wild-type human TWIK-1 channels had a similar cellular distribution with predominant localization in the endosomal compartment. Two-electrode voltage-clamp experiments using Xenopus oocytes showed that both zebrafish and wild-type human TWIK-1 channels produced K(+) currents that are sensitive to external K(+) concentration as well as acidic pH. There were no effects of the three KCNK1 variants on cellular localization, current amplitude or reversal potential at pH7.4 or pH6. Our data indicate that TWIK-1 has a highly conserved role in cardiac function and is required for normal heart rate and atrial morphology. Despite the functional importance of TWIK-1 in the atrium, genetic variation in KCNK1 is not a common primary cause of human AF.