How to Connect Cardiac Excitation to the Atomic Interactions of Ion Channels

Emerging methods to model cardiac ion channel and myocyte electrophysiology

In the field of cardiac electrophysiology, modeling has played a central role for many decades. However, even though the effort is well-established, it has recently seen a rapid and sustained evolution in the complexity and predictive power of the models being created. In particular, new approaches to modeling have allowed the tracking of parallel and interconnected processes that span from the nanometers and femtoseconds that determine ion channel gating to the centimeters and minutes needed to describe an arrhythmia. The connection between scales has brought unprecedented insight into cardiac arrhythmia mechanisms and drug therapies. This review focuses on the generation of these models from first principles, generation of detailed models to describe ion channel kinetics, algorithms to create and numerically solve kinetic models, and new approaches toward data gathering that parameterize these models. While we focus on application of these models for cardiac arrhythmia, these concepts are widely applicable to model the physiology and pathophysiology of any excitable cell.

Models of the cardiac L-type calcium current: A quantitative review

The L-type calcium current (I CaL ) plays a critical role in cardiac electrophysiology, and models ofI CaL are vital tools to predict arrhythmogenicity of drugs and mutations. Five decades of measuring and modelingI CaL have resulted in several competing theories (encoded in mathematical equations). However, the introduction of new models has not typically been accompanied by a data-driven critical comparison with previous work, so that it is unclear which model is best suited for any particular application. In this review, we describe and compare 73 published mammalianI CaL models and use simulated experiments to show that there is a large variability in their predictions, which is not substantially diminished when grouping by species or other categories. We provide model code for 60 models, list major data sources, and discuss experimental and modeling work that will be required to reduce this huge list of competing theories and ultimately develop a community consensus model ofI CaL . This article is categorized under: Cardiovascular Diseases > Computational Models Cardiovascular Diseases > Molecular and Cellular Physiology.

Exploring cardiac form and function: A length-scale computational biology approach

The ability to adequately pump blood throughout the body is the result of tightly regulated feedback mechanisms that exist across many spatial scales in the heart. Diseases which impede the function at any one of the spatial scales can cause detrimental cardiac remodeling and eventual heart failure. An overarching goal of cardiac research is to use engineered heart tissue in vitro to study the physiology of diseased heart tissue, develop cell replacement therapies, and explore drug testing applications. A commonality within the field is to manipulate the flow of mechanical signals across the various spatial scales to direct self-organization and build functional tissue. Doing so requires an understanding of how chemical, electrical, and mechanical cues can be used to alter the cellular microenvironment. We discuss how mathematical models have been used in conjunction with experimental techniques to explore various structure-function relations that exist across numerous spatial scales. We highlight how a systems biology approach can be employed to recapitulate in vivo characteristics in vitro at the tissue, cell, and subcellular scales. Specific focus is placed on the interplay between experimental and theoretical approaches. Various modeling methods are showcased to demonstrate the breadth and power afforded to the systems biology approach. An overview of modeling methodologies exemplifies how the strengths of different scientific disciplines can be used to supplement and/or inspire new avenues of experimental exploration. This article is categorized under: Models of Systems Properties and Processes > Mechanistic Models Models of Systems Properties and Processes > Cellular Models Models of Systems Properties and Processes > Organ, Tissue, and Physiological Models.

A Molecularly Detailed NaV1.5 Model Reveals a New Class I Antiarrhythmic Target

•

Antiarrhythmic therapies remain suboptimal due to our inability to predict how drug interactions with ion channels will affect the ability of the tissue to initiate and sustain an arrhythmia.

•

We built a computational framework that allows for in silico design of precision-targeted therapeutic agents that simultaneously assesses antiarrhythmic markers of success and failure at multiple spatial and time scales. Using this framework, a novel in silico mexiletine “booster” was designed that may dramatically improve the efficacy of mexiletine in suppression of arrhythmia triggers.

•

These results provide a roadmap for the design of novel molecular-based therapy to treat myriad arrhythmia syndromes, including ventricular tachycardia, heart failure arrhythmias, and inherited arrhythmia syndromes.

•

In summary, computational modeling approaches to drug discovery represent a novel tool to design and test precision-targeted therapeutic agents. By exploiting nontraditional ion channel drug targets, an entirely new dimension can be added to the wide parameter space of traditional antiarrhythmic drugs to develop more precision-targeted and potent Class I therapeutic agents.

Antiarrhythmic treatment strategies remain suboptimal due to our inability to predict how drug interactions with ion channels will affect the ability of the tissues to initiate and sustain an arrhythmia. We built a multiscale molecular model of the Na⁺ channel domain III (domain III voltage-sensing domain) to highlight the molecular underpinnings responsible for mexiletine drug efficacy. This model predicts that a hyperpolarizing shift in the domain III voltage-sensing domain is critical for drug efficacy and may be leveraged to design more potent Class I molecules. The model was therefore used to design, in silico, a theoretical mexiletine booster that can dramatically rescue a mutant resistant to the potent antiarrhythmic effects of mexiletine. Our framework provides a strategy for in silico design of precision-targeted therapeutic agents that simultaneously assesses antiarrhythmic markers of success and failure at multiple spatial and time scales. This approach provides a roadmap for the design of novel molecular-based therapy to treat myriad arrhythmia syndromes, including ventricular tachycardia, heart failure arrhythmias, and inherited arrhythmia syndromes.

The Structural Basis of IKs Ion-Channel Activation: Mechanistic Insights from Molecular Simulations

Relating ion channel (iCh) structural dynamics to physiological function remains a challenge. Current experimental and computational techniques have limited ability to explore this relationship in atomistic detail over physiological timescales. A framework associating iCh structure to function is necessary for elucidating normal and disease mechanisms. We formulated a modeling schema that overcomes the limitations of current methods through applications of artificial intelligence machine learning. Using this approach, we studied molecular processes that underlie human IKs voltage-mediated gating. IKs malfunction underlies many debilitating and life-threatening diseases. Molecular components of IKs that underlie its electrophysiological function include KCNQ1 (a pore-forming tetramer) and KCNE1 (an auxiliary subunit). Simulations, using the IKs structure-function model, reproduced experimentally recorded saturation of gating-charge displacement at positive membrane voltages, two-step voltage sensor (VS) movement shown by fluorescence, iCh gating statistics, and current-voltage relationship. Mechanistic insights include the following: 1) pore energy profile determines iCh subconductance; 2) the entire protein structure, not limited to the pore, contributes to pore energy and channel subconductance; 3) interactions with KCNE1 result in two distinct VS movements, causing gating-charge saturation at positive membrane voltages and current activation delay; and 4) flexible coupling between VS and pore permits pore opening at lower VS positions, resulting in sequential gating. The new modeling approach is applicable to atomistic scale studies of other proteins on timescales of physiological function.