The Membrane Dipole Potential and the Roles of Interfacial Water and Lipid Hydrocarbon Chains.

Understanding membrane charge transport processes, including the actions of ion channels, pumps, carriers, and membrane-active peptides, requires a description of the electrostatics of the lipid bilayer. We have simulated a library of different lipid chemistries to reveal the impact of the headgroup, glycerol backbone, and hydrocarbon chains on the membrane dipole potential. We found a strong dependence of the potential on lipid packing, but this was not caused by the packing of lipid polar components, due to cancellation of their electric fields by electrolyte. In contrast, lipid tail contributions were determined by area per lipid, arising from two countering effects. Increased area per lipid leads to chain tilting that increases methylene dipole projections to strengthen the electric field within the bilayer, while at the same time decreasing the electric field from terminal methyl groups. Moreover, electric fields from some nonterminal groups and the terminal methyl group can extend beyond the bilayer center and be canceled by the opposing leaflet. This interleaflet field annulment explains the experimental reduction in dipole potential for unsaturated and branched lipid bilayers, by as much as ∼200 mV, as well as experiments that substitute chain carbons with sulfur. Replacing ester with ether groups (eliminating two carbonyl groups) causes a significant reduction in potential, also by ∼200 mV, in agreement with experiment. We show that the effect can be largely attributed to the loss of aligned water molecules in the glycerol backbone region, lowering the potential inside the bilayer core. When only one of the two carbonyls is removed (using a hybrid ester-ether lipid or a single-chain lipid), most of this reduction in potential was lost, with the single carbonyl group able to maintain full hydration in the interfacial region. While headgroup chemistry can have a major effect (by as much as ±100 mV relative to phosphatidylcholine), anionic headgroups either decrease or increase the dipole potential, with the variation involving perturbation in hydrogen-bonded water molecules and changes in packing of lipid tails. Overall, these results suggest that membrane electrostatics are dominated by aligned water molecules at the polar-hydrocarbon interface and, surprisingly, by the charge distribution of the nonpolar lipid tails, and not the packing of headgroup and glycerol carbonyl dipoles.

Atomistic mechanisms of the regulation of small-conductance Ca2+-activated K+ channel (SK2) by PIP2.

Small-conductance Ca2+-activated K+ channels (SK, KCa2) are gated solely by intracellular microdomain Ca2+. The channel has emerged as a therapeutic target for cardiac arrhythmias. Calmodulin (CaM) interacts with the CaM binding domain (CaMBD) of the SK channels, serving as the obligatory Ca2+ sensor to gate the channels. In heterologous expression systems, phosphatidylinositol 4,5-bisphosphate (PIP2) coordinates with CaM in regulating SK channels. However, the roles and mechanisms of PIP2 in regulating SK channels in cardiomyocytes remain unknown. Here, optogenetics, magnetic nanoparticles, combined with Rosetta structural modeling, and molecular dynamics (MD) simulations revealed the atomistic mechanisms of how PIP2 works in concert with Ca2+-CaM in the SK channel activation. Our computational study affords evidence for the critical role of the amino acid residue R395 in the S6 transmembrane segment, which is localized in propinquity to the intracellular hydrophobic gate. This residue forms a salt bridge with residue E398 in the S6 transmembrane segment from the adjacent subunit. Both R395 and E398 are conserved in all known isoforms of SK channels. Our findings suggest that the binding of PIP2 to R395 residue disrupts the R395:E398 salt bridge, increasing the flexibility of the transmembrane segment S6 and the activation of the channel. Importantly, our findings serve as a platform for testing of structural-based drug designs for therapeutic inhibitors and activators of the SK channel family. The study is timely since inhibitors of SK channels are currently in clinical trials to treat atrial arrhythmias.

Harnessing AlphaFold to reveal hERG channel conformational state secrets.

To design safe, selective, and effective new therapies, there must be a deep understanding of the structure and function of the drug target. One of the most difficult problems to solve has been resolution of discrete conformational states of transmembrane ion channel proteins. An example is K v 11.1 (hERG), comprising the primary cardiac repolarizing current, I kr . hERG is a notorious drug anti-target against which all promising drugs are screened to determine potential for arrhythmia. Drug interactions with the hERG inactivated state are linked to elevated arrhythmia risk, and drugs may become trapped during channel closure. However, the structural details of multiple conformational states have remained elusive. Here, we guided AlphaFold2 to predict plausible hERG inactivated and closed conformations, obtaining results consistent with multiple available experimental data. Drug docking simulations demonstrated hERG state-specific drug interactions in good agreement with experimental results, revealing that most drugs bind more effectively in the inactivated state and are trapped in the closed state. Molecular dynamics simulations demonstrated ion conduction for an open but not AlphaFold2 predicted inactivated state that aligned with earlier studies. Finally, we identified key molecular determinants of state transitions by analyzing interaction networks across closed, open, and inactivated states in agreement with earlier mutagenesis studies. Here, we demonstrate a readily generalizable application of AlphaFold2 as an effective and robust method to predict discrete protein conformations, reconcile seemingly disparate data and identify novel linkages from structure to function.

Toward Digital Twin Technology for Precision Pharmacology.

The authors demonstrate the feasibility of technological innovation for personalized medicine in the context of drug-induced arrhythmia. The authors use atomistic-scale structural models to predict rates of drug interaction with ion channels and make predictions of their effects in digital twins of induced pluripotent stem cell-derived cardiac myocytes. The authors construct a simplified multilayer, 1-dimensional ring model with sufficient path length to enable the prediction of arrhythmogenic dispersion of repolarization. Finally, the authors validate the computational pipeline prediction of drug effects with data and quantify drug-induced propensity to repolarization abnormalities in cardiac tissue. The technology is high throughput, computationally efficient, and low cost toward personalized pharmacologic prediction.

Elucidating molecular mechanisms of protoxin-II state-specific binding to the human NaV1.7 channel.

Human voltage-gated sodium (hNaV) channels are responsible for initiating and propagating action potentials in excitable cells, and mutations have been associated with numerous cardiac and neurological disorders. hNaV1.7 channels are expressed in peripheral neurons and are promising targets for pain therapy. The tarantula venom peptide protoxin-II (PTx2) has high selectivity for hNaV1.7 and is a valuable scaffold for designing novel therapeutics to treat pain. Here, we used computational modeling to study the molecular mechanisms of the state-dependent binding of PTx2 to hNaV1.7 voltage-sensing domains (VSDs). Using Rosetta structural modeling methods, we constructed atomistic models of the hNaV1.7 VSD II and IV in the activated and deactivated states with docked PTx2. We then performed microsecond-long all-atom molecular dynamics (MD) simulations of the systems in hydrated lipid bilayers. Our simulations revealed that PTx2 binds most favorably to the deactivated VSD II and activated VSD IV. These state-specific interactions are mediated primarily by PTx2's residues R22, K26, K27, K28, and W30 with VSD and the surrounding membrane lipids. Our work revealed important protein-protein and protein-lipid contacts that contribute to high-affinity state-dependent toxin interaction with the channel. The workflow presented will prove useful for designing novel peptides with improved selectivity and potency for more effective and safe treatment of pain.

Structural modeling of hERG channel-drug interactions using Rosetta.

The human ether-a-go-go-related gene (hERG) not only encodes a potassium-selective voltage-gated ion channel essential for normal electrical activity in the heart but is also a major drug anti-target. Genetic hERG mutations and blockage of the channel pore by drugs can cause long QT syndrome, which predisposes individuals to potentially deadly arrhythmias. However, not all hERG-blocking drugs are proarrhythmic, and their differential affinities to discrete channel conformational states have been suggested to contribute to arrhythmogenicity. We used Rosetta electron density refinement and homology modeling to build structural models of open-state hERG channel wild-type and mutant variants (Y652A, F656A, and Y652A/F656 A) and a closed-state wild-type channel based on cryo-electron microscopy structures of hERG and EAG1 channels. These models were used as protein targets for molecular docking of charged and neutral forms of amiodarone, nifekalant, dofetilide, d/l-sotalol, flecainide, and moxifloxacin. We selected these drugs based on their different arrhythmogenic potentials and abilities to facilitate hERG current. Our docking studies and clustering provided atomistic structural insights into state-dependent drug-channel interactions that play a key role in differentiating safe and harmful hERG blockers and can explain hERG channel facilitation through drug interactions with its open-state hydrophobic pockets.

A multiscale predictive digital twin for neurocardiac modulation.

Cardiac function is tightly regulated by the autonomic nervous system (ANS). Activation of the sympathetic nervous system increases cardiac output by increasing heart rate and stroke volume, while parasympathetic nerve stimulation instantly slows heart rate. Importantly, imbalance in autonomic control of the heart has been implicated in the development of arrhythmias and heart failure. Understanding of the mechanisms and effects of autonomic stimulation is a major challenge because synapses in different regions of the heart result in multiple changes to heart function. For example, nerve synapses on the sinoatrial node (SAN) impact pacemaking, while synapses on contractile cells alter contraction and arrhythmia vulnerability. Here, we present a multiscale neurocardiac modelling and simulator tool that predicts the effect of efferent stimulation of the sympathetic and parasympathetic branches of the ANS on the cardiac SAN and ventricular myocardium. The model includes a layered representation of the ANS and reproduces firing properties measured experimentally. Model parameters are derived from experiments and atomistic simulations. The model is a first prototype of a digital twin that is applied to make predictions across all system scales, from subcellular signalling to pacemaker frequency to tissue level responses. We predict conditions under which autonomic imbalance induces proarrhythmia and can be modified to prevent or inhibit arrhythmia. In summary, the multiscale model constitutes a predictive digital twin framework to test and guide high-throughput prediction of novel neuromodulatory therapy. KEY POINTS: A multi-layered model representation of the autonomic nervous system that includes sympathetic and parasympathetic branches, each with sparse random intralayer connectivity, synaptic dynamics and conductance based integrate-and-fire neurons generates firing patterns in close agreement with experiment. A key feature of the neurocardiac computational model is the connection between the autonomic nervous system and both pacemaker and contractile cells, where modification to pacemaker frequency drives initiation of electrical signals in the contractile cells. We utilized atomic-scale molecular dynamics simulations to predict the association and dissociation rates of noradrenaline with the ß-adrenergic receptor. Multiscale predictions demonstrate how autonomic imbalance may increase proclivity to arrhythmias or be used to terminate arrhythmias. The model serves as a first step towards a digital twin for predicting neuromodulation to prevent or reduce disease.

Elucidating Molecular Mechanisms of Protoxin-2 State-specific Binding to the Human NaV1.7 Channel.

Human voltage-gated sodium (hNaV) channels are responsible for initiating and propagating action potentials in excitable cells and mutations have been associated with numerous cardiac and neurological disorders. hNaV1.7 channels are expressed in peripheral neurons and are promising targets for pain therapy. The tarantula venom peptide protoxin-2 (PTx2) has high selectivity for hNaV1.7 and serves as a valuable scaffold to design novel therapeutics to treat pain. Here, we used computational modeling to study the molecular mechanisms of the state-dependent binding of PTx2 to hNaV1.7 voltage-sensing domains (VSDs). Using Rosetta structural modeling methods, we constructed atomistic models of the hNaV1.7 VSD II and IV in the activated and deactivated states with docked PTx2. We then performed microsecond-long all-atom molecular dynamics (MD) simulations of the systems in hydrated lipid bilayers. Our simulations revealed that PTx2 binds most favorably to the deactivated VSD II and activated VSD IV. These state-specific interactions are mediated primarily by PTx2's residues R22, K26, K27, K28, and W30 with VSD as well as the surrounding membrane lipids. Our work revealed important protein-protein and protein-lipid contacts that contribute to high-affinity state-dependent toxin interaction with the channel. The workflow presented will prove useful for designing novel peptides with improved selectivity and potency for more effective and safe treatment of pain.

Elucidation of a dynamic interplay between a beta-2 adrenergic receptor, its agonist, and stimulatory G protein.

G protein-coupled receptors (GPCRs) represent the largest group of membrane receptors for transmembrane signal transduction. Ligand-induced activation of GPCRs triggers G protein activation followed by various signaling cascades. Understanding the structural and energetic determinants of ligand binding to GPCRs and GPCRs to G proteins is crucial to the design of pharmacological treatments targeting specific conformations of these proteins to precisely control their signaling properties. In this study, we focused on interactions of a prototypical GPCR, beta-2 adrenergic receptor (ß2AR), with its endogenous agonist, norepinephrine (NE), and the stimulatory G protein (Gs). Using molecular dynamics (MD) simulations, we demonstrated the stabilization of cationic NE, NE(+), binding to ß2AR by Gs protein recruitment, in line with experimental observations. We also captured the partial dissociation of the ligand from ß2AR and the conformational interconversions of Gs between closed and open conformations in the NE(+)-ß2AR-Gs ternary complex while it is still bound to the receptor. The variation of NE(+) binding poses was found to alter Gs α subunit (Gsα) conformational transitions. Our simulations showed that the interdomain movement and the stacking of Gsα α1 and α5 helices are significant for increasing the distance between the Gsα and ß2AR, which may indicate a partial dissociation of Gsα The distance increase commences when Gsα is predominantly in an open state and can be triggered by the intracellular loop 3 (ICL3) of ß2AR interacting with Gsα, causing conformational changes of the α5 helix. Our results help explain molecular mechanisms of ligand and GPCR-mediated modulation of G protein activation.

Structural modeling of the hERG potassium channel and associated drug interactions.

The voltage-gated potassium channel, KV11.1, encoded by the human Ether-à-go-go-Related Gene (hERG), is expressed in cardiac myocytes, where it is crucial for the membrane repolarization of the action potential. Gating of the hERG channel is characterized by rapid, voltage-dependent, C-type inactivation, which blocks ion conduction and is suggested to involve constriction of the selectivity filter. Mutations S620T and S641A/T within the selectivity filter region of hERG have been shown to alter the voltage dependence of channel inactivation. Because hERG channel blockade is implicated in drug-induced arrhythmias associated with both the open and inactivated states, we used Rosetta to simulate the effects of hERG S620T and S641A/T mutations to elucidate conformational changes associated with hERG channel inactivation and differences in drug binding between the two states. Rosetta modeling of the S641A fast-inactivating mutation revealed a lateral shift of the F627 side chain in the selectivity filter into the central channel axis along the ion conduction pathway and the formation of four lateral fenestrations in the pore. Rosetta modeling of the non-inactivating mutations S620T and S641T suggested a potential molecular mechanism preventing F627 side chain from shifting into the ion conduction pathway during the proposed inactivation process. Furthermore, we used Rosetta docking to explore the binding mechanism of highly selective and potent hERG blockers - dofetilide, terfenadine, and E4031. Our structural modeling correlates well with much, but not all, existing experimental evidence involving interactions of hERG blockers with key residues in hERG pore and reveals potential molecular mechanisms of ligand interactions with hERG in an inactivated state.

Rearrangement of a unique Kv1.3 selectivity filter conformation upon binding of a drug.

We report two structures of the human voltage-gated potassium channel (Kv) Kv1.3 in immune cells alone (apo-Kv1.3) and bound to an immunomodulatory drug called dalazatide (dalazatide-Kv1.3). Both the apo-Kv1.3 and dalazatide-Kv1.3 structures are in an activated state based on their depolarized voltage sensor and open inner gate. In apo-Kv1.3, the aromatic residue in the signature sequence (Y447) adopts a position that diverges 11 Å from other K+ channels. The outer pore is significantly rearranged, causing widening of the selectivity filter and perturbation of ion binding within the filter. This conformation is stabilized by a network of intrasubunit hydrogen bonds. In dalazatide-Kv1.3, binding of dalazatide to the channel's outer vestibule narrows the selectivity filter, Y447 occupies a position seen in other K+ channels, and this conformation is stabilized by a network of intersubunit hydrogen bonds. These remarkable rearrangements in the selectivity filter underlie Kv1.3's transition into the drug-blocked state.

A deep learning algorithm to translate and classify cardiac electrophysiology.

The development of induced pluripotent stem cell-derived cardiomyocytes (iPSC-CMs) has been a critical in vitro advance in the study of patient-specific physiology, pathophysiology, and pharmacology. We designed a new deep learning multitask network approach intended to address the low throughput, high variability, and immature phenotype of the iPSC-CM platform. The rationale for combining translation and classification tasks is because the most likely application of the deep learning technology we describe here is to translate iPSC-CMs following application of a perturbation. The deep learning network was trained using simulated action potential (AP) data and applied to classify cells into the drug-free and drugged categories and to predict the impact of electrophysiological perturbation across the continuum of aging from the immature iPSC-CMs to the adult ventricular myocytes. The phase of the AP extremely sensitive to perturbation due to a steep rise of the membrane resistance was found to contain the key information required for successful network multitasking. We also demonstrated successful translation of both experimental and simulated iPSC-CM AP data validating our network by prediction of experimental drug-induced effects on adult cardiomyocyte APs by the latter.

Molecular determinants of pro-arrhythmia proclivity of d- and l-sotalol via a multi-scale modeling pipeline.

Drug isomers may differ in their proarrhythmia risk. An interesting example is the drug sotalol, an antiarrhythmic drug comprising d- and l- enantiomers that both block the hERG cardiac potassium channel and confer differing degrees of proarrhythmic risk. We developed a multi-scale in silico pipeline focusing on hERG channel - drug interactions and used it to probe and predict the mechanisms of pro-arrhythmia risks of the two enantiomers of sotalol. Molecular dynamics (MD) simulations predicted comparable hERG channel binding affinities for d- and l-sotalol, which were validated with electrophysiology experiments. MD derived thermodynamic and kinetic parameters were used to build multi-scale functional computational models of cardiac electrophysiology at the cell and tissue scales. Functional models were used to predict inactivated state binding affinities to recapitulate electrocardiogram (ECG) QT interval prolongation observed in clinical data. Our study demonstrates how modeling and simulation can be applied to predict drug effects from the atom to the rhythm for dl-sotalol and also increased proarrhythmia proclivity of d- vs. l-sotalol when accounting for stereospecific beta-adrenergic receptor blocking.

Molecular Dynamics Simulations Based on Polarizable Models Show that Ion Permeation Interconverts between Different Mechanisms as a Function of Membrane Thickness.

Different mechanisms have been proposed to explain the permeation of charged compounds through lipid membranes. Overall, it is expected that an ion-induced defect permeation mechanism, where substantial membrane deformations accompany ion movement, should be dominant in thin membranes but that a solubility-diffusion mechanism, where ions partition into the membrane core with large associated dehydration energy costs, becomes dominant in thicker membranes. However, while this physical picture is intuitively reasonable, capturing the interconversion between these two permeation mechanisms in molecular dynamics (MD) simulations based on atomic models is challenging. In particular, simulations relying on nonpolarizable force fields are artificially unfavorable to the solubility-diffusion mechanism, as induced polarization of the nonpolar hydrocarbon is ignored, causing overestimated free energy costs for charged molecules to enter into this region of the membrane. In this study, all-atom MD simulations based on nonpolarizable and polarizable force fields are used to quantitatively characterize the permeation process for the arginine side chain analog methyl-guanidinium through bilayer membranes of mono-unsaturated phosphatidylcholine lipids with and without cholesterol, resulting in thicknesses spanning from ∼24 to ∼42 Å. With simulations based on a nonpolarizable force field, ion translocation can take place solely through an ion-induced defect mechanism, with free energy barriers increasing linearly from 14 to 40 kcal/mol, depending on the thickness. However, with simulations based on a polarizable force field, ion translocation is predominantly dominated by an ion-induced defect mechanism in thin membranes, which progressively converts to a solubility-diffusion mechanism as the membranes get thicker. The transition between the two mechanisms occurs at a thickness of ∼29 Å, with lipid tails of 22 or more carbon atoms. This situation appears to represent the upper limit for ion-induced defect permeation within the current polarizable models. Beyond this thickness, it becomes energetically preferable for the ion to dehydrate and partition into the membrane core-a phenomenon that cannot be captured using the nonpolarizable models. Induced electronic polarizability therefore leads not just to a shift in permeation energetics but to an interconversion between two strikingly different physical mechanisms. The result highlights the importance of induced polarizability in modeling lipid membranes.

A Computational Pipeline to Predict Cardiotoxicity: From the Atom to the Rhythm.

RATIONALE: Drug-induced proarrhythmia is so tightly associated with prolongation of the QT interval that QT prolongation is an accepted surrogate marker for arrhythmia. But QT interval is too sensitive a marker and not selective, resulting in many useful drugs eliminated in drug discovery. OBJECTIVE: To predict the impact of a drug from the drug chemistry on the cardiac rhythm. METHODS AND RESULTS: In a new linkage, we connected atomistic scale information to protein, cell, and tissue scales by predicting drug-binding affinities and rates from simulation of ion channel and drug structure interactions and then used these values to model drug effects on the hERG channel. Model components were integrated into predictive models at the cell and tissue scales to expose fundamental arrhythmia vulnerability mechanisms and complex interactions underlying emergent behaviors. Human clinical data were used for model framework validation and showed excellent agreement, demonstrating feasibility of a new approach for cardiotoxicity prediction. CONCLUSIONS: We present a multiscale model framework to predict electrotoxicity in the heart from the atom to the rhythm. Novel mechanistic insights emerged at all scales of the system, from the specific nature of proarrhythmic drug interaction with the hERG channel, to the fundamental cellular and tissue-level arrhythmia mechanisms. Applications of machine learning indicate necessary and sufficient parameters that predict arrhythmia vulnerability. We expect that the model framework may be expanded to make an impact in drug discovery, drug safety screening for a variety of compounds and targets, and in a variety of regulatory processes.

Selectivity filter modalities and rapid inactivation of the hERG1 channel.

The human ether-á-go-go-related gene (hERG1) channel conducts small outward K+ currents that are critical for cardiomyocyte membrane repolarization. The gain-of-function mutation N629D at the outer mouth of the selectivity filter (SF) disrupts inactivation and K+-selective transport in hERG1, leading to arrhythmogenic phenotypes associated with long-QT syndrome. Here, we combined computational electrophysiology with Markov state model analysis to investigate how SF-level gating modalities control selective cation transport in wild-type (WT) and mutant (N629D) hERG1 variants. Starting from the recently reported cryogenic electron microscopy (cryo-EM) open-state channel structure, multiple microseconds-long molecular-dynamics (MD) trajectories were generated using different cation configurations at the filter, voltages, electrolyte concentrations, and force-field parameters. Most of the K+ permeation events observed in hERG1-WT simulations occurred at microsecond timescales, influenced by the spontaneous dehydration/rehydration dynamics at the filter. The SF region displayed conductive, constricted, occluded, and dilated states, in qualitative agreement with the well-documented flickering conductance of hERG1. In line with mutagenesis studies, these gating modalities resulted from dynamic interaction networks involving residues from the SF, outer-mouth vestibule, P-helices, and S5-P segments. We found that N629D mutation significantly stabilizes the SF in a state that is permeable to both K+ and Na+, which is reminiscent of the SF in the nonselective bacterial NaK channel. Increasing the external K+ concentration induced "WT-like" SF dynamics in N629D, in qualitative agreement with the recovery of flickering currents in experiments. Overall, our findings provide an understanding of the molecular mechanisms controlling selective transport in K+ channels with a nonconventional SF sequence.

Atomistic Simulations of Membrane Ion Channel Conduction, Gating, and Modulation.

Membrane ion channels are the fundamental electrical components in the nervous system. Recent developments in X-ray crystallography and cryo-EM microscopy have revealed what these proteins look like in atomic detail but do not tell us how they function. Molecular dynamics simulations have progressed to the point that we can now simulate realistic molecular assemblies to produce quantitative calculations of the thermodynamic and kinetic quantities that control function. In this review, we summarize the state of atomistic simulation methods for ion channels to understand their conduction, activation, and drug modulation mechanisms. We are at a crossroads in atomistic simulation, where long time scale observation can provide unbiased exploration of mechanisms, supplemented by biased free energy methodologies. We illustrate the use of these approaches to describe ion conduction and selectivity in voltage-gated sodium and acid-sensing ion channels. Studies of channel gating present a significant challenge, as activation occurs on longer time scales. Enhanced sampling approaches can ensure convergence on minimum free energy pathways for activation, as illustrated here for pentameric ligand-gated ion channels that are principal to nervous system function and the actions of general anesthetics. We also examine recent studies of local anesthetic and antiepileptic drug binding to a sodium channel, revealing sites and pathways that may offer new targets for drug development. Modern simulations thus offer a range of molecular-level insights into ion channel function and modulation as a learning platform for mechanistic discovery and drug development.

A demonstration of modularity, reuse, reproducibility, portability and scalability for modeling and simulation of cardiac electrophysiology using Kepler Workflows.

Multi-scale computational modeling is a major branch of computational biology as evidenced by the US federal interagency Multi-Scale Modeling Consortium and major international projects. It invariably involves specific and detailed sequences of data analysis and simulation, often with multiple tools and datasets, and the community recognizes improved modularity, reuse, reproducibility, portability and scalability as critical unmet needs in this area. Scientific workflows are a well-recognized strategy for addressing these needs in scientific computing. While there are good examples if the use of scientific workflows in bioinformatics, medical informatics, biomedical imaging and data analysis, there are fewer examples in multi-scale computational modeling in general and cardiac electrophysiology in particular. Cardiac electrophysiology simulation is a mature area of multi-scale computational biology that serves as an excellent use case for developing and testing new scientific workflows. In this article, we develop, describe and test a computational workflow that serves as a proof of concept of a platform for the robust integration and implementation of a reusable and reproducible multi-scale cardiac cell and tissue model that is expandable, modular and portable. The workflow described leverages Python and Kepler-Python actor for plotting and pre/post-processing. During all stages of the workflow design, we rely on freely available open-source tools, to make our workflow freely usable by scientists.

Structural basis for antiarrhythmic drug interactions with the human cardiac sodium channel.

The human voltage-gated sodium channel, hNaV1.5, is responsible for the rapid upstroke of the cardiac action potential and is target for antiarrhythmic therapy. Despite the clinical relevance of hNaV1.5-targeting drugs, structure-based molecular mechanisms of promising or problematic drugs have not been investigated at atomic scale to inform drug design. Here, we used Rosetta structural modeling and docking as well as molecular dynamics simulations to study the interactions of antiarrhythmic and local anesthetic drugs with hNaV1.5. These calculations revealed several key drug binding sites formed within the pore lumen that can simultaneously accommodate up to two drug molecules. Molecular dynamics simulations identified a hydrophilic access pathway through the intracellular gate and a hydrophobic access pathway through a fenestration between DIII and DIV. Our results advance the understanding of molecular mechanisms of antiarrhythmic and local anesthetic drug interactions with hNaV1.5 and will be useful for rational design of novel therapeutics.

Challenges and advances in atomistic simulations of potassium and sodium ion channel gating and permeation.

Ion channels are implicated in many essential physiological events such as electrical signal propagation and cellular communication. The advent of K+ and Na+ ion channel structure determination has facilitated numerous investigations of molecular determinants of their behaviour. At the same time, rapid development of computer hardware and molecular simulation methodologies has made computational studies of large biological molecules in all-atom representation tractable. The concurrent evolution of experimental structural biology with biomolecular computer modelling has yielded mechanistic details of fundamental processes unavailable through experiments alone, such as ion conduction and ion channel gating. This review is a short survey of the atomistic computational investigations of K+ and Na+ ion channels, focusing on KcsA and several voltage-gated channels from the KV and NaV families, which have garnered many successes and engendered several long-standing controversies regarding the nature of their structure-function relationship. We review the latest advancements and challenges facing the field of molecular modelling and simulation regarding the structural and energetic determinants of ion channel function and their agreement with experimental observations.