Reduced connexin-43 expression, slow conduction and repolarisation dispersion in a model of hypertrophic cardiomyopathy.

Hypertrophic cardiomyopathy (HCM) is an inherited heart muscle disease that is characterised by left ventricular wall thickening, cardiomyocyte disarray and fibrosis, and is associated with arrhythmias, heart failure and sudden death. However, it is unclear to what extent the electrophysiological disturbances that lead to sudden death occur secondary to structural changes in the myocardium or as a result of HCM cardiomyocyte electrophysiology. In this study, we used an induced pluripotent stem cell model of the R403Q variant in myosin heavy chain 7 (MYH7) to study the electrophysiology of HCM cardiomyocytes in electrically coupled syncytia, revealing significant conduction slowing and increased spatial dispersion of repolarisation - both well-established substrates for arrhythmia. Analysis of rhythmonome protein expression in MYH7 R403Q cardiomyocytes showed reduced expression of connexin-43 (also known as GJA1), sodium channels and inward rectifier potassium channels - a three-way hit that reduces electrotonic coupling and slows cardiac conduction. Our data represent a previously unreported, biophysical basis for arrhythmia in HCM that is intrinsic to cardiomyocyte electrophysiology. Later in the progression of the disease, these proarrhythmic phenotypes may be accentuated by myocyte disarray and fibrosis to contribute to sudden death.

Development of a pharmaceutical database as an aid to the nonclinical detection of drug-induced cardiac toxicity.

The Health and Environmental Sciences Institute (HESI) Cardiac Safety Committee designed and created a publicly accessible database with an initial set of 128 pharmacologically defined pharmaceutical agents, many with known cardiotoxic properties. The database includes specific information about each compound that could be useful in evaluating hypotheses around mechanisms of drug-induced cardiac toxicity or for development of novel cardiovascular safety assays. Data on each of the compounds was obtained from published literature and online sources (e.g., DrugBank.ca and International Union of Basic and Clinical Pharmacology (IUPHAR) / British Pharmacological Society (BPS) Guide to PHARMACOLOGY) and was curated by 10 subject matter experts. The database includes information such as compound name, pharmacological mode of action, characterized cardiac mode of action, type of cardiac toxicity, known clinical cardiac toxicity profile, animal models used to evaluate the cardiotoxicity profile, routes of administration, and toxicokinetic parameters (i.e., Cmax). Data from both nonclinical and clinical studies are included for each compound. The user-friendly web interface allows for multiple approaches to search the database and is also intended to provide a means for the submission of new data/compounds from relevant users. This will ensure that the database is constantly updated and remains current. Such a data repository will not only aid the HESI working groups in defining drugs for use in any future studies, but safety scientists can also use the database as a vehicle of support for broader cardiovascular safety studies or exploring mechanisms of toxicity associated with certain pharmacological modes of action.

Assessing drug safety by identifying the axis of arrhythmia in cardiomyocyte electrophysiology.

Many classes of drugs can induce fatal cardiac arrhythmias by disrupting the electrophysiology of cardiomyocytes. Safety guidelines thus require all new drugs to be assessed for pro-arrhythmic risk prior to conducting human trials. The standard safety protocols primarily focus on drug blockade of the delayed-rectifier potassium current (IKr). Yet the risk is better assessed using four key ion currents (IKr, ICaL, INaL, IKs). We simulated 100,000 phenotypically diverse cardiomyocytes to identify the underlying relationship between the blockade of those currents and the emergence of ectopic beats in the action potential. We call that relationship the axis of arrhythmia. It serves as a yardstick for quantifying the arrhythmogenic risk of any drug from its profile of multi-channel block alone. We tested it on 109 drugs and found that it predicted the clinical risk labels with an accuracy of 88.1-90.8%. Pharmacologists can use our method to assess the safety of novel drugs without resorting to animal testing or unwieldy computer simulations.

Functional coupling between Piezo1 and TRPM4 influences the electrical activity of HL-1 atrial myocytes.

The transient receptor potential melastatin 4 (TRPM4) channel contributes extensively to cardiac electrical activity, especially cardiomyocyte action potential formation. Mechanical stretch can induce changes in heart rate and rhythm, and the mechanosensitive channel Piezo1 is expressed in many cell types within the myocardium. Our previous study showed that TRPM4 and Piezo1 are closely co-localized in the t-tubules of ventricular cardiomyocytes and contribute to the Ca2+ -dependent signalling cascade that underlies hypertrophy in response to mechanical pressure overload. However, there was no direct evidence showing that Piezo1 activation was related to TRPM4 activation in situ. In the present study, we employed the HL-1 mouse atrial myocyte-like cell line as an in vitro model to investigate whether Piezo1-TRPM4 coupling can affect action potential properties. We used the small molecule Piezo1 agonist, Yoda1, as a surrogate for mechanical stretch to activate Piezo1 and detected the action potential changes in HL-1 cells using FluoVolt, a fluorescent voltage sensitive dye. Our results demonstrate that Yoda1-induced activation of Piezo1 changes the action potential frequency in HL-1 cells. This change in action potential frequency is reduced by Piezo1 knockdown using small intefering RNA. Importantly knockdown or pharmacological inhibition of TRPM4 significantly affected the degree to which Yoda1-evoked Piezo1 activation influenced action potential frequency. Thus, the present study provides in vitro evidence of a functional coupling between Piezo1 and TRPM4 in a cardiomyocyte-like cell line. The coupling of a mechanosensitive Ca2+ permeable channel and a Ca2+ -activated TRP channel probably represents a ubiquitous model for the role of TRP channels in mechanosensory transduction. KEY POINTS: The transient receptor potential melastatin 4 (TRPM4) and Piezo1 channels have been confirmed to contribute to the Ca2+ -dependent signalling cascade that underlies cardiac hypertrophy in response to mechanical pressure overload. However, there was no direct evidence showing that Piezo1 activation was related to TRPM4 activation in situ. We employed the HL-1 mouse atrial myocyte-like cell line as an in vitro model to investigate the effect of Piezo1-TRPM4 coupling on cardiac electrical properties. The results show that both pharmacological and genetic inhibition of TRPM4 significantly affected the degree to which Piezo1 activation influenced action potential frequency in HL-1 cells. Our findings provide in vitro evidence of a functional coupling between Piezo1 and TRPM4 in a cardiomyocyte-like cell line. The coupling of a mechanosensitive Ca2+ permeable channel and a Ca2+ -activated TRP channel probably represents a ubiquitous model for the role of TRP channels in mechanosensory transduction in various (patho)physiological processes.

Empirical Quantification of Predictive Uncertainty Due to Model Discrepancy by Training with an Ensemble of Experimental Designs: An Application to Ion Channel Kinetics.

When using mathematical models to make quantitative predictions for clinical or industrial use, it is important that predictions come with a reliable estimate of their accuracy (uncertainty quantification). Because models of complex biological systems are always large simplifications, model discrepancy arises-models fail to perfectly recapitulate the true data generating process. This presents a particular challenge for making accurate predictions, and especially for accurately quantifying uncertainty in these predictions. Experimentalists and modellers must choose which experimental procedures (protocols) are used to produce data used to train models. We propose to characterise uncertainty owing to model discrepancy with an ensemble of parameter sets, each of which results from training to data from a different protocol. The variability in predictions from this ensemble provides an empirical estimate of predictive uncertainty owing to model discrepancy, even for unseen protocols. We use the example of electrophysiology experiments that investigate the properties of hERG potassium channels. Here, 'information-rich' protocols allow mathematical models to be trained using numerous short experiments performed on the same cell. In this case, we simulate data with one model and fit it with a different (discrepant) one. For any individual experimental protocol, parameter estimates vary little under repeated samples from the assumed additive independent Gaussian noise model. Yet parameter sets arising from the same model applied to different experiments conflict-highlighting model discrepancy. Our methods will help select more suitable ion channel models for future studies, and will be widely applicable to a range of biological modelling problems.

Development of new adeno-associated virus capsid variants for targeted gene delivery to human cardiomyocytes.

Recombinant adeno-associated viruses (rAAVs) have emerged as one of the most promising gene therapy vectors that have been successfully used in pre-clinical models of heart disease. However, this has not translated well to humans due to species differences in rAAV transduction efficiency. As a result, the search for human cardiotropic capsids is a major contemporary challenge. We used a capsid-shuffled rAAV library to perform directed evolution in human iPSC-derived cardiomyocytes (hiPSC-CMs). Five candidates emerged, with four presenting high sequence identity to AAV6, while a fifth divergent variant was related to AAV3b. Functional analysis of the variants was performed in vitro using hiPSC-CMs, cardiac organoids, human cardiac slices, non-human primate and porcine cardiac slices, as well as mouse heart and liver in vivo. We showed that cell entry was not the best predictor of transgene expression efficiency. The novel variant rAAV.KK04 was the best-performing vector in human-based screening platforms, exceeding the benchmark rAAV6. None of the novel capsids demonstrate a significant transduction of liver in vivo. The range of experimental models used revealed the value of testing for tropism differences under the conditions of human specificity, bona fide, myocardium and cell type of interest.

Development of a robust induced pluripotent stem cell atrial cardiomyocyte differentiation protocol to model atrial arrhythmia.

BACKGROUND: Atrial fibrillation is the most common arrhythmia syndrome and causes significant morbidity and mortality. Current therapeutics, however, have limited efficacy. Notably, many therapeutics shown to be efficacious in animal models have not proved effective in humans. Thus, there is a need for a drug screening platform based on human tissue. The aim of this study was to develop a robust protocol for generating atrial cardiomyocytes from human-induced pluripotent stem cells. METHODS: A novel protocol for atrial differentiation, with optimized timing of retinoic acid during mesoderm formation, was compared to two previously published methods. Each differentiation method was assessed for successful formation of a contractile syncytium, electrical properties assayed by optical action potential recordings and multi-electrode array electrophysiology, and response to the G-protein-gated potassium channel activator, carbamylcholine. Atrial myocyte monolayers, derived using the new differentiation protocol, were further assessed for cardiomyocyte purity, gene expression, and the ability to form arrhythmic rotors in response to burst pacing. RESULTS: Application of retinoic acid at day 1 of mesoderm formation resulted in a robust differentiation of atrial myocytes with contractile syncytium forming in 16/18 differentiations across two cell lines. Atrial-like myocytes produced have shortened action potentials and field potentials, when compared to standard application of retinoic acid at the cardiac mesoderm stage. Day 1 retinoic acid produced atrial cardiomyocytes are also carbamylcholine sensitive, indicative of active Ikach currents, which was distinct from ventricular myocytes and standard retinoic addition in matched differentiations. A current protocol utilizing reduced Activin A and BMP4 can produce atrial cardiomyocytes with equivalent functionality but with reduced robustness of differentiation; only 8/17 differentiations produced a contractile syncytium. The day 1 retinoic acid protocol was successfully applied to 6 iPSC lines (3 male and 3 female) without additional optimization or modification. Atrial myocytes produced could also generate syncytia with rapid conduction velocities, > 40 cm s-1, and form rotor style arrhythmia in response to burst pacing. CONCLUSIONS: This method combines an enhanced atrial-like phenotype with robustness of differentiation, which will facilitate further research in human atrial arrhythmia and myopathies, while being economically viable for larger anti-arrhythmic drug screens.

Action Potential Morphology Accurately Predicts Proarrhythmic Risk for Drugs With Potential to Prolong Cardiac Repolarization.

BACKGROUND: Drug-induced or acquired long QT syndrome occurs as a result of the unintended disruption of cardiac repolarization due to drugs that block cardiac ion channels. These side effects have been responsible for the withdrawal of a range of drugs from market and are a common reason for termination of the development of new drugs in the preclinical stage. Existing approaches to risk prediction are expensive and overly sensitive meaning that recently there have been renewed efforts, largely driven by the comprehensive proarrhythmic assay initiative, to develop more accurate methods for allocation of proarrhythmic risk. METHODS: In this study, we aimed to quantify changes in the morphology of the repolarization phase of the cardiac action potential as an indicator of proarrhythmia, supposing that these shape changes might precede the emergence of ectopic depolarizations that trigger arrhythmia. To do this, we describe a new method of quantifying action potential morphology by measuring the radius of curvature of the repolarization phase both in simulated action potentials, as well as in action potentials measured from induced pluripotent stem cell-derived cardiomyocytes. Features derived from the curvature signal were used as inputs for logistic regressions to predict proarrhythmic risk. RESULTS: Optimal risk classifiers based on morphology were able to correctly classify risk to drugs in the comprehensive proarrhythmic assay initiative panels with very high accuracy (0.9375) and outperformed conventional metrics based on action potential duration at 90% repolarization, triangulation, and charge movement (qNet). CONCLUSIONS: Analysis of action potential morphology in response to proarrhythmic drugs improves prediction of torsadogenic risk. Furthermore, morphology metrics can be measured directly from the action potential, potentially eliminating the burden of undertaking complex screens of potency and drug-binding kinetics against multiple cardiac ion channels. As such, this method has the potential to improve and streamline regulatory assessment of proarrhythmia in preclinical drug development.

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.

High throughput measurement of hERG drug block kinetics using the CiPA dynamic protocol.

The Comprehensive in vitro Proarrhythmic Assay (CiPA) has promoted use of in silico models of drug effects on cardiac repolarization to improve proarrhythmic risk prediction. These models contain a pharmacodynamic component describing drug binding to hERG channels that required in vitro data for kinetics of block, in addition to potency, to constrain them. To date, development and validation has been undertaken using data from manual patch-clamp. The application of this approach at scale requires the development of a high-throughput, automated patch-clamp (APC) implementation. Here, we present a comprehensive analysis of the implementation of the Milnes, or CiPA dynamic protocol, on an APC platform, including quality control and data analysis. Kinetics and potency of block were assessed for bepridil, cisapride, terfenadine and verapamil with data retention/QC pass rate of 21.8% overall, or as high as 50.4% when only appropriate sweep lengths were considered for drugs with faster kinetics. The variability in IC50 and kinetics between manual and APC was comparable to that seen between sites/platforms in previous APC studies of potency. Whilst the experimental success is less than observed in screens of potency alone, it is still significantly greater than manual patch. With the modifications to protocol design, including sweep length, number of repetitions, and leak correction recommended in this study, this protocol can be applied on APC to acquire data comparable to manual patch clamp.

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.

Metabolic and electrolyte abnormalities as risk factors in drug-induced long QT syndrome.

Drug-induced long QT syndrome (diLQTS) is the phenomenon by which the administration of drugs causes prolongation of cardiac repolarisation and leads to an increased risk of the ventricular tachycardia known as torsades de pointes (TdP). In most cases of diLQTS, the primary molecular target is the human ether-à-go-go-related gene protein (hERG) potassium channel, which carries the rapid delayed rectifier current (IKr) in the heart. However, the proarrhythmic risk associated with drugs that block hERG can be modified in patients by a range of environmental- and disease-related factors, such as febrile temperatures, alterations in pH, dyselectrolytaemias such as hypokalaemia and hypomagnesemia and coadministration with other drugs. In this review, we will discuss the clinical occurrence of drug-induced LQTS in the context of these modifying factors as well as the mechanisms by which they contribute to altered hERG potency and proarrhythmic risk.

Multitarget nociceptor sensitization by a promiscuous peptide from the venom of the King Baboon spider.

The King Baboon spider, Pelinobius muticus, is a burrowing African tarantula. Its impressive size and appealing coloration are tempered by reports describing severe localized pain, swelling, itchiness, and muscle cramping after accidental envenomation. Hyperalgesia is the most prominent symptom after bites from P. muticus, but the molecular basis by which the venom induces pain is unknown. Proteotranscriptomic analysis of P. muticus venom uncovered a cysteine-rich peptide, δ/κ-theraphotoxin-Pm1a (δ/κ-TRTX-Pm1a), that elicited nocifensive behavior when injected into mice. In small dorsal root ganglion neurons, synthetic δ/κ-TRTX-Pm1a (sPm1a) induced hyperexcitability by enhancing tetrodotoxin-resistant sodium currents, impairing repolarization and lowering the threshold of action potential firing, consistent with the severe pain associated with envenomation. The molecular mechanism of nociceptor sensitization by sPm1a involves multimodal actions over several ion channel targets, including NaV1.8, KV2.1, and tetrodotoxin-sensitive NaV channels. The promiscuous targeting of peptides like δ/κ-TRTX-Pm1a may be an evolutionary adaptation in pain-inducing defensive venoms.

Pathophysiological metabolic changes associated with disease modify the proarrhythmic risk profile of drugs with potential to prolong repolarisation.

BACKGROUND AND PURPOSE: Hydroxychloroquine, chloroquine and azithromycin are three drugs that were proposed to treat coronavirus disease 2019 (COVID-19). While concern already existed around their proarrhythmic potential, there are little data regarding how altered physiological states encountered in patients such as febrile state, electrolyte imbalances or acidosis might change their risk profiles. EXPERIMENTAL APPROACH: Potency of human ether-à-go-go related gene (hERG) block was measured using high-throughput electrophysiology in the presence of variable environmental factors. These potencies informed simulations to predict population risk profiles. Effects on cardiac repolarisation were verified in human induced pluripotent stem cell-derived cardiomyocytes from multiple individuals. KEY RESULTS: Chloroquine and hydroxychloroquine blocked hERG with IC50 of 1.47 ± 0.07 and 3.78 ± 0.17 µM, respectively, indicating proarrhythmic risk at concentrations effective against severe acute respiratory syndrome-coronovirus-2 (SARS-CoV-2) in vitro. Hypokalaemia and hypermagnesaemia increased potency of chloroquine and hydroxychloroquine, indicating increased proarrhythmic risk. Acidosis significantly reduced potency of all drugs, whereas increased temperature decreased potency of chloroquine and hydroxychloroquine against hERG but increased potency for azithromycin. In silico simulations demonstrated that proarrhythmic risk was increased by female sex, hypokalaemia and heart failure and identified specific genetic backgrounds associated with emergence of arrhythmia. CONCLUSION AND IMPLICATIONS: Our study demonstrates how proarrhythmic risk can be exacerbated by metabolic changes and pre-existing disease. More broadly, the study acts as a blueprint for how high-throughput in vitro screening, combined with in silico simulations, can help guide both preclinical screening and clinical management of patients in relation to drugs with potential to prolong repolarisation.

Heterozygous KCNH2 variant phenotyping using Flp-In HEK293 and high-throughput automated patch clamp electrophysiology.

KCNH2 is one of the 59 medically actionable genes recommended by the American College of Medical Genetics for reporting of incidental findings from clinical genomic sequencing. However, half of the reported KCNH2 variants in the ClinVar database are classified as variants of uncertain significance. In the absence of strong clinical phenotypes, there is a need for functional phenotyping to help decipher the significance of variants identified incidentally. Here, we report detailed methods for assessing the molecular phenotype of any KCNH2 missense variant. The key components of the assay include quick and cost-effective generation of a bi-cistronic vector to co-express Wild-type (WT) and any KCNH2 variant allele, generation of stable Flp-In HEK293 cell lines and high-throughput automated patch clamp electrophysiology analysis of channel function. Stable cell lines take 3-4 weeks to produce and can be generated in bulk, which will then allow up to 30 variants to be phenotyped per week after 48 h of channel expression. This high-throughput functional genomics assay will enable a much more rapid assessment of the extent of loss of function of any KCNH2 variant.

Conserved Role of the Large Conductance Calcium-Activated Potassium Channel, KCa1.1, in Sinus Node Function and Arrhythmia Risk.

BACKGROUND: KCNMA1 encodes the α-subunit of the large-conductance Ca2+-activated K+ channel, KCa1.1, and lies within a linkage interval for atrial fibrillation (AF). Insights into the cardiac functions of KCa1.1 are limited, and KCNMA1 has not been investigated as an AF candidate gene. METHODS: The KCNMA1 gene was sequenced in 118 patients with familial AF. The role of KCa1.1 in normal cardiac structure and function was evaluated in humans, mice, zebrafish, and fly. A novel KCNMA1 variant was functionally characterized. RESULTS: A complex KCNMA1 variant was identified in 1 kindred with AF. To evaluate potential disease mechanisms, we first evaluated the distribution of KCa1.1 in normal hearts using immunostaining and immunogold electron microscopy. KCa1.1 was seen throughout the atria and ventricles in humans and mice, with strong expression in the sinus node. In an ex vivo murine sinoatrial node preparation, addition of the KCa1.1 antagonist, paxilline, blunted the increase in beating rate induced by adrenergic receptor stimulation. Knockdown of the KCa1.1 ortholog, kcnma1b, in zebrafish embryos resulted in sinus bradycardia with dilatation and reduced contraction of the atrium and ventricle. Genetic inactivation of the Drosophila KCa1.1 ortholog, slo, systemically or in adult stages, also slowed the heartbeat and produced fibrillatory cardiac contractions. Electrophysiological characterization of slo-deficient flies revealed bursts of action potentials, reflecting increased events of fibrillatory arrhythmias. Flies with cardiac-specific overexpression of the human KCNMA1 mutant also showed increased heart period and bursts of action potentials, similar to the KCa1.1 loss-of-function models. CONCLUSIONS: Our data point to a highly conserved role of KCa1.1 in sinus node function in humans, mice, zebrafish, and fly and suggest that KCa1.1 loss of function may predispose to AF.

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.

Co-expression of calcium and hERG potassium channels reduces the incidence of proarrhythmic events.

AIMS: Cardiac electrical activity is extraordinarily robust. However, when it goes wrong it can have fatal consequences. Electrical activity in the heart is controlled by the carefully orchestrated activity of more than a dozen different ion conductances. While there is considerable variability in cardiac ion channel expression levels between individuals, studies in rodents have indicated that there are modules of ion channels whose expression co-vary. The aim of this study was to investigate whether meta-analytic co-expression analysis of large-scale gene expression datasets could identify modules of co-expressed cardiac ion channel genes in human hearts that are of functional importance. METHODS AND RESULTS: Meta-analysis of 3653 public human RNA-seq datasets identified a strong correlation between expression of CACNA1C (L-type calcium current, ICaL) and KCNH2 (rapid delayed rectifier K+ current, IKr), which was also observed in human adult heart tissue samples. In silico modelling suggested that co-expression of CACNA1C and KCNH2 would limit the variability in action potential duration seen with variations in expression of ion channel genes and reduce susceptibility to early afterdepolarizations, a surrogate marker for proarrhythmia. We also found that levels of KCNH2 and CACNA1C expression are correlated in human-induced pluripotent stem cell-derived cardiac myocytes and the levels of CACNA1C and KCNH2 expression were inversely correlated with the magnitude of changes in repolarization duration following inhibition of IKr. CONCLUSION: Meta-analytic approaches of multiple independent human gene expression datasets can be used to identify gene modules that are important for regulating heart function. Specifically, we have verified that there is co-expression of CACNA1C and KCNH2 ion channel genes in human heart tissue, and in silico analyses suggest that CACNA1C-KCNH2 co-expression increases the robustness of cardiac electrical activity.