Electrophysiological characterization of the hERG R56Q LQTS variant and targeted rescue by the activator RPR260243

Targeted activation of human ether-à-go-go-related gene channels rescues electrical instability induced by the R56Q+/- long QT syndrome variant

AIMS

Long QT syndrome type 2 (LQTS2) is associated with inherited variants in the cardiac human ether-à-go-go-related gene (hERG) K+ channel. However, the pathogenicity of hERG channel gene variants is often uncertain. Using CRISPR-Cas9 gene-edited hiPSC-derived cardiomyocytes (hiPSC-CMs), we investigated the pathogenic mechanism underlying the LQTS-associated hERG R56Q variant and its phenotypic rescue by using the Type 1 hERG activator, RPR260243.

METHODS AND RESULTS

The above approaches enable characterization of the unclear causative mechanism of arrhythmia in the R56Q variant (an N-terminal PAS domain mutation that primarily accelerates channel deactivation) and translational investigation of the potential for targeted pharmacologic manipulation of hERG deactivation. Using perforated patch clamp electrophysiology of single hiPSC-CMs, programmed electrical stimulation showed that the hERG R56Q variant does not significantly alter the mean action potential duration (APD90). However, the R56Q variant increases the beat-to-beat variability in APD90 during pacing at constant cycle lengths, enhances the variance of APD90 during rate transitions, and increases the incidence of 2:1 block. During paired S1-S2 stimulations measuring electrical restitution properties, the R56Q variant was also found to increase the variability in rise time and duration of the response to premature stimulations. Application of the hERG channel activator, RPR260243, reduces the APD variance in hERG R56Q hiPSC-CMs, reduces the variability in responses to premature stimulations, and increases the post-repolarization refractoriness.

CONCLUSION

Based on our findings, we propose that the hERG R56Q variant leads to heterogeneous APD dynamics, which could result in spatial dispersion of repolarization and increased risk for re-entry without significantly affecting the average APD90. Furthermore, our data highlight the antiarrhythmic potential of targeted slowing of hERG deactivation gating, which we demonstrate increases protection against premature action potentials and reduces electrical heterogeneity in hiPSC-CMs.

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.

hERG stereoselective modulation by mexiletine-derived ureas: Molecular docking study, synthesis, and biological evaluation

Long QT syndrome (LQTS) is a disorder of cardiac electrophysiology resulting in life-threatening arrhythmias; nowadays, only a few drugs are available for the management of LQTS. Focusing our attention on LQT2, one of the most common subtypes of LQTS caused by mutations in the human ether-à-go-go-related gene (hERG), in the present work, the stereoselectivity of the recently discovered mexiletine-derived urea 8 was investigated on the hERG potassium channel. According to preliminary in silico predictions, in vitro studies revealed a stereoselective behavior, with the meso form showing the greatest hERG opening activity. In addition, functional studies on guinea pig isolated left atria, aorta, and ileum demonstrated that 8 does not present any cardiac or intestinal liability in our ex vivo studies. Due to its overall profile, (R,S)-8 paves the way for the design and development of a new series of compounds potentially useful in the treatment of both congenital and drug-induced forms of LQTS.

Binding of RPR260243 at the intracellular side of the hERG1 channel pore domain slows closure of the helix bundle crossing gate

The opening and closing of voltage-dependent potassium channels is dependent on a tight coupling between movement of the voltage sensing S4 segments and the activation gate. A specific interaction between intracellular amino- and carboxyl-termini is required for the characteristically slow rate of channel closure (deactivation) of hERG1 channels. Compounds that increase hERG1 channel currents represent a novel approach for prevention of arrhythmia associated with prolonged ventricular repolarization. RPR260243 (RPR), a quinoline oxo-propyl piperidine derivative, inhibits inactivation and dramatically slows the rate of hERG1 channel deactivation. Here we report that similar to its effect on wild-type channels, RPR greatly slows the deactivation rate of hERG1 channels missing their amino-termini, or of split channels lacking a covalent link between the voltage sensor domain and the pore domain. By contrast, RPR did not slow deactivation of C-terminal truncated hERG1 channels or D540K hERG1 mutant channels activated by hyperpolarization. Together, these findings indicate that ability of RPR to slow deactivation requires an intact C-terminus, does not slow deactivation by stabilizing an interaction involving the amino-terminus or require a covalent link between the voltage sensor and pore domains. All-atom molecular dynamics simulations using the cryo-EM structure of the hERG1 channel revealed that RPR binds to a pocket located at the intracellular ends of helices S5 and S6 of a single subunit. The slowing of channel deactivation by RPR may be mediated by disruption of normal S5-S6 interactions.

Ion channel model reduction using manifold boundaries

Mathematical models of voltage-gated ion channels are used in basic research, industrial and clinical settings. These models range in complexity, but typically contain numerous variables representing the proportion of channels in a given state, and parameters describing the voltage-dependent rates of transition between states. An open problem is selecting the appropriate degree of complexity and structure for an ion channel model given data availability. Here, we simplify a model of the cardiac human Ether-à-go-go related gene (hERG) potassium ion channel, which carries cardiac I_Kr, using the manifold boundary approximation method (MBAM). The MBAM approximates high-dimensional model-output manifolds by reduced models describing their boundaries, resulting in models with fewer parameters (and often variables). We produced a series of models of reducing complexity starting from an established five-state hERG model with 15 parameters. Models with up to three fewer states and eight fewer parameters were shown to retain much of the predictive capability of the full model and were validated using experimental hERG1a data collected in HEK293 cells at 37°C. The method provides a way to simplify complex models of ion channels that improves parameter identifiability and will aid in future model development.

Pharmacological activation of the hERG K + channel for the management of the long QT syndrome: A review

In the human heart, the rapid delayed rectifier K⁺ current (I_Kr) contributes significantly to ventricular action potential (AP) repolarization and to set the duration of the QT interval of the surface electrocardiogram (ECG). The pore‐forming (α) subunit of the I_Kr channel is encoded by KCNH2 or human ether‐à‐go‐go‐related gene 1 (hERG1). Impairment of hERG function through either gene mutation (congenital) or pharmacological blockade by diverse drugs in clinical use (acquired) can cause a prolongation of the AP duration (APD) reflected onto the surface ECG as a prolonged QT interval or Long QT Syndrome (LQTS). LQTS can increase the risk of triggered activity of ventricular cardiomyocytes and associated life‐threatening arrhythmia. Current treatments all focus on reducing the incidence of arrhythmia or terminating it after its onset but there is to date no prophylactic treatment for the pharmacological management of LQTS. A new class of hERG modulators (agonists) have been suggested through direct interaction with the hERG channel to shorten the action potential duration (APD) and/or increase the postrepolarisation refractoriness period (PRRP) of ventricular cardiomyocytes protecting thereby against triggered activity and associated arrhythmia. Although promising drug candidates, there remain major obstacles to their clinical development. The aim of this review is to summarize the latest advances as well as the limitations of this proposed pharmacotherapy.