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 IKr, 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 5-state hERG model with 15 parameters. Models with up to 3 fewer states and 8 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

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

Human Ether-à-go-go (hERG) channels contribute to cardiac repolarization, and inherited variants or drug block are associated with long QT syndrome type 2 (LQTS2) and arrhythmia. Therefore, hERG activator compounds present a therapeutic opportunity for targeted treatment of LQTS. However, a limiting concern is over-activation of hERG resurgent current during the action potential and abbreviated repolarization. Activators that slow deactivation gating (type I), such as RPR260243, may enhance repolarizing hERG current during the refractory period, thus ameliorating arrhythmogenicity with reduced early repolarization risk. Here, we show that, at physiological temperature, RPR260243 enhances hERG channel repolarizing currents conducted in the refractory period in response to premature depolarizations. This occurs with little effect on the resurgent hERG current during the action potential. The effects of RPR260243 were particularly evident in LQTS2-associated R56Q mutant channels, whereby RPR260243 restored WT-like repolarizing drive in the early refractory period and diastolic interval, combating attenuated protective currents. In silico kinetic modeling of channel gating predicted little effect of the R56Q mutation on hERG current conducted during the action potential and a reduced repolarizing protection against afterdepolarizations in the refractory period and diastolic interval, particularly at higher pacing rates. These simulations predicted partial rescue from the arrhythmic effects of R56Q by RPR260243 without risk of early repolarization. Our findings demonstrate that the pathogenicity of some hERG variants may result from reduced repolarizing protection during the refractory period and diastolic interval with limited effect on action potential duration, and that the hERG channel activator RPR260243 may provide targeted antiarrhythmic potential in these cases