Enhanced Sampling Molecular Dynamics Simulations Research Articles

Gas adsorption and framework flexibility of CALF-20 explored via experiments and simulations

In 2021, Svante, in collaboration with BASF, reported successful scale up of CALF-20 production, a stable MOF with high capacity for post-combustion CO2 capture which exhibits remarkable stability towards water. CALF-20’s success story in the MOF commercialisation space provides new thinking about appropriate structural and adsorptive metrics important for CO2 capture. Here, we combine atomistic-level simulations with experiments to study adsorptive properties of CALF-20 and shed light on its flexible crystal structure. We compare measured and predicted CO2 and water adsorption isotherms and explain the role of water-framework interactions and hydrogen bonding networks in CALF-20’s hydrophobic behaviour. Furthermore, regular and enhanced sampling molecular dynamics simulations are performed with both density-functional theory (DFT) and machine learning potentials (MLPs) trained to DFT energies and forces. From these simulations, the effects of adsorption-induced flexibility in CALF-20 are uncovered. We envisage this work would encourage development of other MOF materials useful for CO2 capture applications in humid conditions.

Base pair dynamics, electrostatics, and thermodynamics at the LTR-III quadruplex:duplex junction

G-quadruplexes (GQs) play key regulatory roles within the human genome and have also been identified to play similar roles in other eukaryotes, bacteria, archaea, and viruses. Human immunodeficiency virus 1, the etiological agent of acquired immunodeficiency syndrome, can form two GQs in its long terminal repeat (LTR) promoter region, each of which act to regulate viral gene expression in opposing manners. The major LTR GQ, called LTR-III, is a distinct hybrid GQ containing a 12-nucleotide duplex loop attached to the quadruplex motif. The resulting quadruplex:duplex junction (QDJ) has been hypothesized to serve as a selective drug targeting site. To better understand the dynamics of this QDJ, we performed conventional and enhanced-sampling molecular dynamics simulations using the Drude-2017 force field. We observed unbiased and reversible formation of additional base pairs in the QDJ, between Ade4:Thy14 and Gua3:Thy14. Both base pairs were electrostatically favored, but geometric constraints within the junction may drive the formation of, and preference for, the Ade4:Thy14 base pair. Finally, we demonstrated that the base pairs are separated only by small energy barriers that may enable transitions between both base-paired states. Together, these simulations provide new insights into the dynamics, electrostatics, and thermodynamics of the LTR-III QDJ.

Interaction of the Immune System TIM-3 Protein with a Model Cellular Membrane Containing Phosphatidyl-Serine Lipids.

T cell transmembrane, Immunoglobulin, and Mucin (TIM) are important immune system proteins which are especially present in T-cells and regulated the immune system by sensing cell engulfment and apoptotic processes. Their role is exerted by the capacity to detect the presence of phosphatidyl-serine lipid polar head in the outer leaflet of cellular membranes (correlated with apoptosis). In this contribution by using equilibrium and enhanced sampling molecular dynamics simulation we unravel the molecular bases and the thermodynamics of TIM, and in particular TIM-3, interaction with phosphatidyl serine in a lipid bilayer. Since TIM-3 deregulation is an important factor of pro-oncogenic tumor micro-environment understanding its functioning at a molecular level may pave the way to the development of original immunotherapeutic approaches.

Unraveling the Interplay of Extracellular Domain Conformational Changes and Parathyroid Hormone Type 1 Receptor Activation in Class B1 G Protein-Coupled Receptors: Integrating Enhanced Sampling Molecular Dynamics Simulations and Markov State Models.

Parathyroid hormone (PTH) type 1 receptor (PTH1R), as a typical class B1 G protein-coupled receptor (GPCR), is responsible for regulating bone turnover and maintaining calcium homeostasis, and its dysregulation has been implicated in the development of several diseases. The extracellular domain (ECD) of PTH1R is crucial for the recognition and binding of ligands, and the receptor may exhibit an autoinhibited state with the closure of the ECD in the absence of ligands. However, the correlation between ECD conformations and PTH1R activation remains unclear. Thus, this study combines enhanced sampling molecular dynamics (MD) simulations and Markov state models (MSMs) to reveal the possible relevance between the ECD conformations and the activation of PTH1R. First, 22 intermediate structures are generated from the autoinhibited state to the active state and conducted for 10 independent 200 ns simulations each. Then, the MSM is constructed based on the cumulative 44 μs simulations with six identified microstates. Finally, the potential interplay between ECD conformational changes and PTH1R activation as well as cryptic allosteric pockets in the intermediate states during receptor activation is revealed. Overall, our findings reveal that the activation of PTH1R has a specific correlation with ECD conformational changes and provide essential insights for GPCR biology and developing novel allosteric modulators targeting cryptic sites.

Investigating the allosteric activation mechanism of protein phosphatase 2A using enhanced sampling molecular dynamics simulations

Investigating the allosteric activation mechanism of protein phosphatase 2A using enhanced sampling molecular dynamics simulations

Computational Probing of the Folding Mechanism of Human Telomeric G-Quadruplex DNA.

The human telomeric (htel) sequences in the terminal regions of human telomeres form diverse G-quadruplex (GQ) structures. Despite much experimental efforts to elucidate the folding pathways of htel GQ, no comprehensive model of htel GQ folding has been presented. Here, we describe folding pathways of the htel GQ determined by state-of-the-art enhanced sampling molecular dynamics simulation at the all-atom level. Briefly, GQ folding is initiated by the formation of a single-hairpin and then followed by the formation of double-hairpins, which then branch via distinct folding pathways to produce different GQ topologies (antiparallel chair, antiparallel basket, hybrids 1 and 2, and parallel propeller). In addition to these double-hairpin states, three-triad and two-tetrad structures in antiparallel backbone alignment serve as key intermediates that connect the GQ folding and transition between two different GQs.

Involvement of the alpha-subunit N-terminus in the mechanism of the Na+,K+-ATPase.

Previous studies have shown that cytoplasmic K+ release and the associated E2→E1 conformational change of the Na+,K+-ATPase is a major rate-determining step of the enzyme's ion pumping cycle and hence a prime site of acute regulatory intervention. From the ionic strength dependence of the enzyme's distribution between the E2 and E1 states, it has also been found that E2 is stabilized by an electrostatic attraction. Any disruption of this electrostatic attraction would, thus, have profound effects on the rate of ion pumping. The aim of this paper is to identify the location of this interaction. Using enhanced-sampling molecular dynamics simulations with a predicted N-terminal structure added to the X-ray crystal structure of the Na+,K+-ATPase, a previously postulated salt bridge between Lys32 and Glu233 (rat sequence numbering) of the enzyme's α-subunit can be excluded. The residues never approach closely enough to form a salt bridge. In contrast, strong interactions with anionic lipid head groups were seen. To investigate the possibility of a protein-lipid interaction experimentally, the surface charge density of Na+,K+-ATPase-containing membrane fragments was estimated from zeta potential measurements to be 0.019 (± 0.001) C m-2. This is in good agreement with the charge density previously determined to be responsible for stabilization of the E2 state of 0.023 (± 0.009) C m-2 and the membrane charge density estimated here from published electron-microscopic images of 0.018C m-2. The results are, therefore, consistent with an interaction of the Na+,K+-ATPase α-subunit N-terminus with negatively-charged lipid head groups of the neighbouring cytoplasmic membrane surface as the origin of the electrostatic interaction stabilising the E2 state.

Determinants of sugar-induced influx in the mammalian fructose transporter GLUT5

In mammals, glucose transporters (GLUT) control organism-wide blood-glucose homeostasis. In human, this is accomplished by 14 different GLUT isoforms, that transport glucose and other monosaccharides with varying substrate preferences and kinetics. Nevertheless, there is little difference between the sugar-coordinating residues in the GLUT proteins and even the malarial Plasmodium falciparum transporter PfHT1, which is uniquely able to transport a wide range of different sugars. PfHT1 was captured in an intermediate ‘occluded’ state, revealing how the extracellular gating helix TM7b has moved to break and occlude the sugar-binding site. Sequence difference and kinetics indicated that the TM7b gating helix dynamics and interactions likely evolved to enable substrate promiscuity in PfHT1, rather than the sugar-binding site itself. It was unclear, however, if the TM7b structural transitions observed in PfHT1 would be similar in the other GLUT proteins. Here, using enhanced sampling molecular dynamics simulations, we show that the fructose transporter GLUT5 spontaneously transitions through an occluded state that closely resembles PfHT1. The coordination of D-fructose lowers the energetic barriers between the outward- and inward-facing states, and the observed binding mode for D-fructose is consistent with biochemical analysis. Rather than a substrate-binding site that achieves strict specificity by having a high affinity for the substrate, we conclude GLUT proteins have allosterically coupled sugar binding with an extracellular gate that forms the high-affinity transition-state instead. This substrate-coupling pathway presumably enables the catalysis of fast sugar flux at physiological relevant blood-glucose concentrations.

Multivalent Molecular Tweezers Disrupt the Essential NDC80 Interaction with Microtubules.

Binding of microtubule filaments by the conserved Ndc80 protein is required for kinetochore-microtubule attachments in cells and the successful distribution of the genetic material during cell division. The reversible inhibition of microtubule binding is an important aspect of the physiological error correction process. Small molecule inhibitors of protein-protein interactions involving Ndc80 are therefore highly desirable, both for mechanistic studies of chromosome segregation and also for their potential therapeutic value. Here, we report on a novel strategy to develop rationally designed inhibitors of the Ndc80 Calponin-homology domain using Supramolecular Chemistry. With a multiple-click approach, lysine-specific molecular tweezers were assembled to form covalently fused dimers to pentamers with a different overall size and preorganization/stiffness. We identified two dimers and a trimer as efficient Ndc80 CH-domain binders and have shown that they disrupt the interaction between Ndc80 and microtubules at low micromolar concentrations without affecting microtubule dynamics. NMR spectroscopy allowed us to identify the biologically important lysine residues 160 and 204 as preferred tweezer interaction sites. Enhanced sampling molecular dynamics simulations provided a rationale for the binding mode of multivalent tweezers and the role of pre-organization and secondary interactions in targeting multiple lysine residues across a protein surface.

Interplay between VSD, pore and membrane lipids in electromechanical coupling in HCN channels.

Hyperpolarized-activated and Cyclic Nucleotide-gated (HCN) channels are the only members of the voltage-gated ion channel superfamily in mammals that open upon hyperpolarization, conferring them pacemaker properties that are instrumental for rhythmic firing of cardiac and neuronal cells. Activation of their voltage-sensor domains (VSD) upon hyperpolarization occurs through a downward movement of the S4 helix bearing the gating charges, which triggers a break in the alpha-helical hydrogen bonding pattern at the level of a conserved Serine residue. Previous structural and molecular simulation studies had however failed to capture pore opening that should be triggered by VSD activation, presumably because of a low VSD/pore electromechanical coupling efficiency and the limited timescales accessible to such techniques. Here, we have used advanced modeling strategies, including enhanced sampling molecular dynamics simulations exploiting comparisons between non-domain swapped voltage-gated ion channel structures trapped in closed and open states to trigger pore gating and characterize electromechanical coupling in HCN1. We propose that the coupling mechanism involves the reorganization of the interfaces between the VSD helices, in particular S4, and the pore-forming helices S5 and S6, subtly shifting the balance between hydrophobic and hydrophilic interactions in a 'domino effect' during activation and gating in this region. Remarkably, our simulations reveal state-dependent occupancy of lipid molecules at this emergent coupling interface, suggesting a key role of lipids in hyperpolarization-dependent gating. Our model provides a rationale for previous observations and a possible mechanism for regulation of HCN channels by the lipidic components of the membrane.

Toward Overcoming Pyrethroid Resistance in Mosquito Control: The Role of Sodium Channel Blocker Insecticides.

Diseases spread by mosquitoes lead to the death of 700,000 people each year. The main way to reduce transmission is vector control by biting prevention with chemicals. However, the most commonly used insecticides lose efficacy due to the growing resistance. Voltage-gated sodium channels (VGSCs), membrane proteins responsible for the depolarizing phase of an action potential, are targeted by a broad range of neurotoxins, including pyrethroids and sodium channel blocker insecticides (SCBIs). Reduced sensitivity of the target protein due to the point mutations threatened malaria control with pyrethroids. Although SCBIs-indoxacarb (a pre-insecticide bioactivated to DCJW in insects) and metaflumizone-are used in agriculture only, they emerge as promising candidates in mosquito control. Therefore, a thorough understanding of molecular mechanisms of SCBIs action is urgently needed to break the resistance and stop disease transmission. In this study, by performing an extensive combination of equilibrium and enhanced sampling molecular dynamics simulations (3.2 μs in total), we found the DIII-DIV fenestration to be the most probable entry route of DCJW to the central cavity of mosquito VGSC. Our study revealed that F1852 is crucial in limiting SCBI access to their binding site. Our results explain the role of the F1852T mutation found in resistant insects and the increased toxicity of DCJW compared to its bulkier parent compound, indoxacarb. We also delineated residues that contribute to both SCBIs and non-ester pyrethroid etofenprox binding and thus could be involved in the target site cross-resistance.

The Interaction of Factor Xa and IXa with Non-Activated Antithrombin in Michaelis Complex: Insights from Enhanced-Sampling Molecular Dynamics Simulations.

The interaction between coagulation factors Xa and IXa and the activated state of their inhibitor, antithrombin (AT),have been investigated using X-ray diffraction studies. However, only mutagenesis data are available for non-activated AT. Our aim was to propose a model based on docking and advanced-sampling molecular dynamics simulations that can reveal the conformational behavior of the systems when AT is not binding a pentasaccharide. We built the initial structure for non-activated AT-FXa and AT-FIXa complexes using HADDOCK 2.4. The conformational behavior was studied using Gaussian accelerated molecular dynamics simulations. In addition to the docked complexes, two systems based on the X-ray structures were also simulated, with and without the ligand. The simulations revealed large variability in conformation for both factors. In the docking-based complex of AT-FIXa, conformations with stable Arg150-AT interactions can exist for longer time periods but the system also has a higher tendency for reaching states with very limited interaction with the "exosite" of AT. By comparing simulations with or without the pentasaccharide, we were able to gain insights into the effects of conformational activation on the Michaelis complexes. RMSF analysis and correlation calculations for the alpha-carbon atoms revealed important details of the allosteric mechanisms. Our simulations provide atomistic models for better understanding the conformational activation mechanism of AT against its target factors.

On the Binding Mode and Molecular Mechanism of Enzymatic Polyethylene Terephthalate Degradation

Enzymatic degradation of polyethylene terephthlate (PET) by polyester hydrolases is currently subject to intensive research, as it is considered as a potential eco-friendly recycling method for plastic waste. However, the substrate-binding mode and the molecular mechanism of enzymatic PET hydrolysis are still under intense investigation, and controversial hypotheses have been presented. To help unravel the inherent mechanism of biocatalytic PET degradation at the atomic level, we performed solid-state NMR measurements of a cutinase from Thermobifida fusca (TfCut2) embedded in trehalose glasses together with chemically synthesized, amorphous 13C(═O)-labeled oligomeric PET. The resulting ternary enzyme-PET-trehalose glassy system enabled advanced solid-state NMR methods for real-time tracking of the enzymatic PET degradation and the investigation of PET chain dynamics. Combined with enhanced-sampling molecular dynamics simulations, specific enzyme–substrate interactions during the degradation process could also be monitored. Our results demonstrate that the PET chain is first cleaved by TfCut2 in blocks of at least one repeat unit and further to terephthalic acid and ethylene glycol. Moreover, the second step (formation of final hydrolysis products) appears to be rate-limiting in such reactions. The observed dynamic changes and interfacial protein contacts of 13C-labeled PET carbonyl groups suggest that only one PET repeat unit is bound to the enzyme during the degradation process while the rest of the PET chain is only loosely confined to the active site. These results, not accessible by using conventional solution enzyme samples and small nonhydrolyzable substrates, provide a better understanding of the biocatalytic PET degradation mechanism of polyester hydrolases.

Self-assembling peptide biomaterials: Insights from spontaneous and enhanced sampling molecular dynamics simulations

Peptide self-assembly is the process by which peptide molecules aggregate into low dimensional (1D, 2D) or 3D ordered materials with potential applications ranging from drug delivery to electronics. Short peptides are particularly good candidates for forming supramolecular assemblies due to the relatively simple structure and ease of modulating their self-assembly process to achieve required material properties. The experimental resolution of fibrous peptide-based nanomaterials as 3D atomic coordinates remains challenging. For surface-mediated peptide assembly in particular, it is typically not feasible to resolve multiple conformationally distinct surface bound peptide structures by experiment. The mechanisms of peptide self-assembly also remain elusive due to the interchange of complex interactions and multiple time and length scales involved in the self-assembly process. Peptide self-assembly in solution, or mediated by surfaces, is driven by specific interactions between the peptides and water, competing interactions within the peptide and/or between peptide aggregate units and, in the latter case, an interplay of the interactions between peptides and solvent molecules for adsorption onto a proximal surface. Computational methodologies have proven beneficial in elucidating the structures formed during peptide self-assembly and the molecular mechanisms driving it, and hence have scope in facilitating the development of functional peptide-based nanomaterials for medical or biotechnological applications. In this perspective, computational methods that have provided molecular insights into the mechanisms of formation of peptide biomaterials, and the all-atom-resolved structures of peptide assemblies are presented. Established and recently emerged molecular simulation approaches are reviewed with a focus on applications relevant to peptide assembly, including all-atom and coarse-grained “brute force” molecular dynamics methods as well as the enhanced sampling methodologies: umbrella sampling, steered and replica exchange molecular dynamics, and variants of metadynamics. These approaches have been shown to contribute all-atom details not yet available experimentally, to advance our understanding of peptide self-assembly processes and biomaterial formation. The scope of this review includes a summary of the current state of the computational methods, in terms of their strengths and limitations for application to self-assembling peptide biomaterials.

Substituted Oligosaccharides as Protein Mimics: Deep Learning Free Energy Landscapes.

Protein-protein complexes power the majority of cellular processes. Interfering with the formation of such complexes using well-designed mimics is a difficult, yet actively pursued, research endeavor. Due to the limited availability of results on the conformational preferences of oligosaccharides compared to polypeptides, the former have been much less explored than the latter as protein mimics, despite interesting ADMET characteristics. In this work, the conformational landscapes of a series of 956 substituted glucopyranose oligomers of lengths 3 to 12 designed as protein interface mimics are revealed using microsecond-time-scale, enhanced-sampling molecular dynamics simulations. Deep convolutional networks are trained on these large conformational ensembles, to predict the stability of longer oligosaccharide structures from those of their constituent trimer motifs. Deep generative adversarial networks are then designed to suggest plausible conformations for oligosaccharide mimics of arbitrary length and substituent sequences that can subsequently be used as input to docking simulations. Analyzing the performance of the neural networks also yields insights into the intricate collective effects that dominate oligosaccharide conformational dynamics.

Amide-to-ester substitution as a stable alternative to N-methylation for increasing membrane permeability in cyclic peptides

Naturally occurring peptides with high membrane permeability often have ester bonds on their backbones. However, the impact of amide-to-ester substitutions on the membrane permeability of peptides has not been directly evaluated. Here we report the effect of amide-to-ester substitutions on the membrane permeability and conformational ensemble of cyclic peptides related to membrane permeation. Amide-to-ester substitutions are shown to improve the membrane permeability of dipeptides and a model cyclic hexapeptide. NMR-based conformational analysis and enhanced sampling molecular dynamics simulations suggest that the conformational transition of the cyclic hexapeptide upon membrane permeation is differently influenced by an amide-to-ester substitution and an amide N-methylation. The effect of amide-to-ester substitution on membrane permeability of other cyclic hexapeptides, cyclic octapeptides, and a cyclic nonapeptide is also investigated to examine the scope of the substitution. Appropriate utilization of amide-to-ester substitution based on our results will facilitate the development of membrane-permeable peptides.

Gating interactions steer loop conformational changes in the active site of the L1 metallo-β-lactamase.

β-Lactam antibiotics are the most important and widely used antibacterial agents across the world. However, the widespread dissemination of β-lactamases among pathogenic bacteria limits the efficacy of β-lactam antibiotics. This has created a major public health crisis. The use of β-lactamase inhibitors has proven useful in restoring the activity of β-lactam antibiotics, yet, effective clinically approved inhibitors against class B metallo-β-lactamases are not available. L1, a class B3 enzyme expressed by Stenotrophomonas maltophilia, is a significant contributor to the β-lactam resistance displayed by this opportunistic pathogen. Structurally, L1 is a tetramer with two elongated loops, α3-β7 and β12-α5, present around the active site of each monomer. Residues in these two loops influence substrate/inhibitor binding. To study how the conformational changes of the elongated loops affect the active site in each monomer, enhanced sampling molecular dynamics simulations were performed, Markov State Models were built, and convolutional variational autoencoder-based deep learning was applied. The key identified residues (D150a, H151, P225, Y227, and R236) were mutated and the activity of the generated L1 variants was evaluated in cell-based experiments. The results demonstrate that there are extremely significant gating interactions between α3-β7 and β12-α5 loops. Taken together, the gating interactions with the conformational changes of the key residues play an important role in the structural remodeling of the active site. These observations offer insights into the potential for novel drug development exploiting these gating interactions.

Integrating comparative modeling and accelerated simulations reveals conformational and energetic basis of actomyosin force generation

Muscle contraction is performed by arrays of contractile proteins in the sarcomere. Serious heart diseases, such as cardiomyopathy, can often be results of mutations in myosin and actin. Direct characterization of how small changes in the myosin-actin complex impact its force production remains challenging. Molecular dynamics (MD) simulations, although capable of studying protein structure-function relationships, are limited owing to the slow timescale of the myosin cycle as well as a lack of various intermediate structures for the actomyosin complex. Here, employing comparative modeling and enhanced sampling MD simulations, we show how the human cardiac myosin generates force during the mechanochemical cycle. Initial conformational ensembles for different myosin-actin states are learned from multiple structural templates with Rosetta. This enables us to efficiently sample the energy landscape of the system using Gaussian accelerated MD. Key myosin loop residues, whose substitutions are related to cardiomyopathy, are identified to form stable or metastable interactions with the actin surface. We find that the actin-binding cleft closure is allosterically coupled to the myosin motor core transitions and ATP-hydrolysis product release from the active site. Furthermore, a gate between switch I and switch II is suggested to control phosphate release at the prepowerstroke state. Our approach demonstrates the ability to link sequence and structural information to motor functions.

Stabilized SARS-CoV-2 spike S2 immunogens designed using insights from enhanced sampling molecular dynamics simulations.

The SARS-CoV-2 spike protein is an important molecular machine that facilitates membrane fusion to host cells and is the major antigenic component of vaccines. Notably, full-length prefusion-stabilized spike proteins produce a strong neutralizing antibody response against S1, which contains the N-terminal domain (NTD) and receptor binding domain (RBD). Unfortunately, emerging variants of concern contain numerous mutations, insertions, and deletions within S1 that impact the binding of neutralizing antibodies which may diminish the efficacy of current vaccines. In contrast, S2 is the most conserved region of the spike and is capable of eliciting antibodies that bind broadly to diverse β-coronavirus spikes. Although S2-based immunogens may have potential benefits, a challenge to its usage as an alternative vaccine strategy is the general instability of S2 in the absence of S1. Using enhanced sampling molecular dynamics simulations, we observed that the S2 head adopts an open conformation where individual protomers separate from one another at the apex. Extensive contact map analysis reveals that the formation and breakage of contacts between three helices mediates these opening motions. Through introducing mutations in these helices, we designed a stabilized S2 construct that maintains a closed conformation and has enhanced thermal stability and protein expression. These studies represent the first simulations characterizing S2 conformational dynamics and demonstrate the potential of using enhanced sampling molecular dynamics simulations in the design of stabilized SARS-CoV-2 S2 immunogens.

Voltage-dependent transitions between resting and activated states of the four Cav1.1 voltage sensors observed in MD simulation.

The sliding helix hypothesis describes how, upon depolarization and repolarization of the membrane, the positive gating charges of the S4 helices of voltage-sensing domains (VSDs) one by one move across the membrane electrical field and thereby effect the opening and closing of ion channels. Voltage-gated calcium channels contain four VSDs, which together regulate gating of a common pore. Although these four VSDs are highly homologous, they differ from each other, e.g. in the number of gating charges and in the number and position of their countercharges. Recent atomic-resolution structures of voltage-gated calcium and sodium channels delineate the endpoints of the VSD movement upon depolarization, showing the S4 helix in the activated up-state. While the first structures of VSDs in resting states have been solved for homotetrameric (prokaryotic) channels, resting state structures of the (pseudo-)heterotetrameric eukaryotic sodium and calcium channels are still lacking. Accordingly, structural models showing the transitions from the resting to the activated state and back are still elusive. Here we addressed this problem using enhanced sampling molecular dynamics simulations of the CaV1.1 calcium channel exposed to a strong membrane potential. Indeed, starting from the known up-state, the S4 helices of the four VSDs moved down to positions closely matching the experimentally determined resting state of the bacterial NaVAb channel. Importantly, the trajectories, the transiently formed ionic interactions between gating- and countercharges, and the speed of the S4 transitions differed substantially between the four VSDs. This is consistent with data from recent mutagenesis and voltage-clamp fluorometry studies on CaV1.1 and allows relating the individual VSD movements with their putative roles in the regulation of distinct channel gating properties, as well as their role in the activation of excitation-contraction coupling in skeletal muscle.