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 \u27domino effect\u27 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

Calmodulin acts as a state-dependent switch to control a cardiac potassium channel opening

Calmodulin (CaM) and phosphatidylinositol 4,5-bisphosphate (PI

Does Proton Conduction in the Voltage-Gated H+ Channel hHv1 Involve Grotthuss-Like Hopping via Acidic Residues?

Hv1s are ubiquitous highly selective voltage-gated proton channels involved in male fertility, immunology, and the invasiveness of certain forms of breast cancer. The mechanism of proton extrusion in Hv1 is not yet understood, while it constitutes the first step toward the design of high-affinity drugs aimed at this important pharmacological target. In this contribution, we explore the details of the mechanism via an integrative approach, using classical and QM/MM molecular dynamics simulations of a monomeric hHv1 model. We propose that protons localize in three binding sites along the channel lumen, formed by three pairs of conserved negatively charged residues lining the pore: D174/E153, D112/D185, and E119/D123. Local rearrangements, involving notably a dihedral transition of F150, a conserved phenylalanine lining the permeation pathway, appear to allow protons to hop from one acidic residue to the next through a bridging water molecule. These results constitute a first attempt at rationalizing hHv1 selectivity for H+ and the role played by D112 in this process. They pave the way for further quantitative characterization of H+ transport in hHv1

Canine CNGA3 Gene Mutations Provide Novel Insights Into Human Achromatopsia-Associated Channelopathies and Treatment

Cyclic nucleotide-gated (CNG) ion channels are key mediators underlying signal transduction in retinal and olfactory receptors. Genetic defects in CNGA3 and CNGB3, encoding two structurally related subunits of cone CNG channels, lead to achromatopsia (ACHM). ACHM is a congenital, autosomal recessive retinal disorder that manifests by cone photoreceptor dysfunction, severely reduced visual acuity, impaired or complete color blindness and photophobia. Here, we report the first canine models for CNGA3-associated channelopathy caused by R424W or V644del mutations in the canine CNGA3 ortholog that accurately mimic the clinical and molecular features of human CNGA3-associated ACHM. These two spontaneous mutations exposed CNGA3 residues essential for the preservation of channel function and biogenesis. The CNGA3-R424W results in complete loss of cone function in vivoand channel activity confirmed by in vitro electrophysiology. Structural modeling and molecular dynamics (MD) simulations revealed R424-E306 salt bridge formation and its disruption with the R424W mutant. Reversal of charges in a CNGA3-R424E-E306R double mutant channel rescued cGMP-activated currents uncovering new insights into channel gating. The CNGA3-V644del affects the C-terminal leucine zipper (CLZ) domain destabilizing intersubunit interactions of the coiled-coil complex in the MD simulations; the in vitro experiments showed incompetent trimeric CNGA3 subunit assembly consistent with abnormal biogenesis of in vivochannels. These newly characterized large animal models not only provide a valuable system for studying cone-specific CNG channel function in health and disease, but also represent prime candidates for proof-of-concept studies of CNGA3 gene replacement therapy for ACHM patients

A Structurally Precise Mechanism Links an Epilepsy-Associated KCNC2 Potassium Channel Mutation to Interneuron Dysfunction

De novo heterozygous variants in KCNC2 encoding the voltage-gated potassium (K+) channel subunit Kv3.2 are a recently described cause of developmental and epileptic encephalopathy (DEE). A de novo variant in KCNC2 c.374G \u3e A (p.Cys125Tyr) was identified via exome sequencing in a patient with DEE. Relative to wild-type Kv3.2, Kv3.2-p.Cys125Tyr induces K+ currents exhibiting a large hyperpolarizing shift in the voltage dependence of activation, accelerated activation, and delayed deactivation consistent with a relative stabilization of the open conformation, along with increased current density. Leveraging the cryogenic electron microscopy (cryo-EM) structure of Kv3.1, molecular dynamic simulations suggest that a strong π-π stacking interaction between the variant Tyr125 and Tyr156 in the α-6 helix of the T1 domain promotes a relative stabilization of the open conformation of the channel, which underlies the observed gain of function. A multicompartment computational model of a Kv3-expressing parvalbumin-positive cerebral cortex fast-spiking γ-aminobutyric acidergic (GABAergic) interneuron (PV-IN) demonstrates how the Kv3.2-Cys125Tyr variant impairs neuronal excitability and dysregulates inhibition in cerebral cortex circuits to explain the resulting epilepsy

Cryo-EM structure of the human Kv3.1 channel reveals gating control by the cytoplasmic T1 domain

Kv3 channels have distinctive gating kinetics tailored for rapid repolarization in fast-spiking neurons. Malfunction of this process due to genetic variants in the KCNC1 gene causes severe epileptic disorders, yet the structural determinants for the unusual gating properties remain elusive. Here, we present cryo-electron microscopy structures of the human Kv3.1a channel, revealing a unique arrangement of the cytoplasmic tetramerization domain T1 which facilitates interactions with C-terminal axonal targeting motif and key components of the gating machinery. Additional interactions between S1/S2 linker and turret domain strengthen the interface between voltage sensor and pore domain. Supported by molecular dynamics simulations, electrophysiological and mutational analyses, we identify several residues in the S4/S5 linker which influence the gating kinetics and an electrostatic interaction between acidic residues in α6 of T1 and R449 in the pore-flanking S6T helices. These findings provide insights into gating control and disease mechanisms and may guide strategies for the design of pharmaceutical drugs targeting Kv3 channels

Sharing data from molecular simulations

Given the need for modern researchers to produce open, reproducible scientific output, the lack of standards and best practices for sharing data and workflows used to produce and analyze molecular dynamics (MD) simulations has become an important issue in the field. There are now multiple well-established packages to perform molecular dynamics simulations, often highly tuned for exploiting specific classes of hardware, each with strong communities surrounding them, but with very limited interoperability/transferability options. Thus, the choice of the software package often dictates the workflow for both simulation production and analysis. The level of detail in documenting the workflows and analysis code varies greatly in published work, hindering reproducibility of the reported results and the ability for other researchers to build on these studies. An increasing number of researchers are motivated to make their data available, but many challenges remain in order to effectively share and reuse simulation data. To discuss these and other issues related to best practices in the field in general, we organized a workshop in November 2018 (https://bioexcel.eu/events/workshop-on-sharing-data-from-molecular-simulations/). Here, we present a brief overview of this workshop and topics discussed. We hope this effort will spark further conversation in the MD community to pave the way toward more open, interoperable, and reproducible outputs coming from research studies using MD simulations