Functional Annotation of Ion Channel Structures by Molecular Simulation

Molecular mechanism of anion permeation through aquaporin 6

Aquaporins (AQPs) are recognized as transmembrane water channels that facilitate selective water permeation through their monomeric pores. Among the AQP family, AQP6 has an intriguing characteristic as an anion channel, which is allosterically controlled by pH conditions and is eliminated by a single amino acid mutation. However, the molecular mechanism of anion permeation through AQP6 remains unclear. Using molecular dynamics simulations in the presence of a transmembrane voltage utilizing an ion concentration gradient, we show that chloride ions permeate through the pore corresponding to the central axis of the AQP6 homotetramer. Under low pH conditions, a subtle opening of the hydrophobic selectivity filter (SF), located near the extracellular part of the central pore, becomes wetted and enables anion permeation. Our simulations also indicate that a single mutation (N63G) in human AQP6, located at the central pore, significantly reduces anion conduction, consistent with experimental data. Moreover, we demonstrate that the pH-sensing mechanism in which the protonation of H184 and H189 under low pH conditions allosterically triggers the gating of the SF region. These results suggest a unique pH-dependent allosteric anion permeation mechanism in AQP6 and could clarify the role of the central pore in some of the AQP tetramers.

Lipid and cholesterols modulate the dynamics of SARS-CoV-2 viral ion channel ORF3a and its pathogenic variants

SARS-CoV-2 accessory protein, ORF3a is a putative ion channel which immensely contributes to viral pathogenicity by modulating host immune responses and virus-host interactions. Relatively high expression of ORF3a in diseased individuals and implication with inflammasome activation, apoptosis and autophagy inhibition, ratifies as an effective target for developing vaccines and therapeutics. Herein, we present the elusive dynamics of ORF3a-dimeric state using all-atoms molecular dynamics (MD) simulations at μ-seconds scale in a heterogeneous lipid-mimetic system in multiple replicates. Additionally, we also explore the effect of non-synonymous pathogenic mutations on ORF3a ion channel activity and viral pathogenicity in different SARS-CoV-2 variants using various structure-based protein stability (ΔΔG) tools and computational saturation mutagenesis. Our study ascertains the role of phosphatidylcholines and cholesterol in modulating the structure of ORF3a, which perturbs the size and flexibility of the polar cavity that allows permeation of large cations. Discrete trend in ion channel pore radius and area per lipid arises the premise that presence of lipids might also affect the overall conformation of ORF3a. MD structural-ensembles, in some replicates rationalize the crucial role of TM2 in maintaining the native structure of ORF3a. We also infer that loss of structural stability primarily grounds for pathogenicity in more than half of the pathogenic variants of ORF3a. Overall, the effect of mutation on alteration of ion permeability of ORF3a, proposed in this study brings mechanistic insights into variant consequences on viral membrane proteins of SARS-CoV-2, which can be utilized for the development of novel therapeutics to treat COVID-19 and other coronavirus diseases.

Inner pore hydration free energy controls the activation of big potassium channels

Hydrophobic gating is an emerging mechanism in regulation of protein ion channels where the pore remains physically open but becomes dewetted to block ion permeation. Atomistic molecular dynamics simulations have played a crucial role in understanding hydrophobic gating by providing the molecular details to complement mutagenesis and structural studies. However, existing studies rely on direct simulations and do not quantitatively describe how the sequence and structural changes may control the delicate liquid-vapor equilibrium of confined water in the pore of the channel protein. To address this limitation, we explore two enhanced sampling methods, namely metadynamics and umbrella sampling, to derive free-energy profiles of pore hydration in both the closed and open states of big potassium (BK) channels, which are important in cardiovascular and neural systems. It was found that metadynamics required substantially longer sampling times and struggled to generate stably converged free-energy profiles due to the slow dynamics of cooperative pore water diffusion even in the barrierless limit. Using umbrella sampling, well-converged free-energy profiles can be readily generated for the wild-type BK channels as well as three mutants with pore-lining mutations experimentally known to dramatically perturb the channel gating voltage. The results show that the free energy of pore hydration faithfully reports the gating voltage of the channel, providing further support for hydrophobic gating in BK channels. Free-energy analysis of pore hydration should provide a powerful approach for quantitative studies of how protein sequence, structure, solution conditions, and/or drug binding may modulate hydrophobic gating in ion channels.

Molecular Dynamics Refinement of Open State Serotonin 5-HT3A Receptor Structures

Pentameric ligand-gated ion channels play an important role in mediating fast neurotransmissions. As a member of this receptor family, cation-selective 5-HT₃ receptors are a clinical target for treating nausea and vomiting associated with chemotherapy and radiation therapy (Thompson and Lummis, 2006). Multiple cryo-electron microscopy (cryo-EM) structures of 5-HT₃ receptors have been determined in distinct functional states (e.g., open, closed, etc.) (Basak et al., 2018; Basak et al., 2018; Polovinkin et al., 2018; Zhang et al., 2015). However, recent work has shown that the transmembrane pores of the open 5-HT₃ receptor structures rapidly collapse and become artificially asymmetric in molecular dynamics (MD) simulations. To avoid this hydrophobic collapse, Dämgen and Biggin developed an equilibration protocol that led to a stable open state structure of the glycine receptor in MD simulations (Dämgen and Biggin, 2020). However, the protocol failed to yield open-like structures of the 5-HT₃ receptor in our simulations. Here, we present a refined equilibration protocol that involves the rearrangement of the transmembrane helices to achieve stable open state structures of the 5-HT₃ receptor that allow both water and ion permeation through the channel. Notably, channel gating is mediated through collective movement of the transmembrane helices, involving not only pore lining M2 helices but also their cross-talk with the adjacent M1 and M3 helices. Thus, the successful application of our refined equilibration protocol underscores the importance of the conformational coupling between the transmembrane helices in stabilizing open-like structures of the 5-HT₃ receptor.

Arginine transportation mechanism through cationic amino acid transporter 1: insights from molecular dynamics studies

Metabolic and signaling mechanisms in mammalian cells are facilitated by the transportation of L-arginine (Arg) across the plasma membrane through cationic amino acid transporter (CAT) proteins. Due to a lack of argininosuccinate synthase (ASS) activity in various tumor cells such as acute myeloid leukemia, acute lymphocytic leukemia, and chronic lymphocytic leukemia, these tumor entities are arginine-auxotrophic and therefore depend on the uptake of the amino acid arginine. Cationic amino acid transporter-1 (CAT-1) is the leading arginine importer expressed in the aforementioned tumor entities. Hence, in the present study, to investigate the transportation mechanism of arginine in CAT-1, we performed molecular dynamics (MD) simulation methods on the modeled human CAT-1. The MM-PBSA approach was conducted to determine the critical residues interacting with arginine within the corresponding binding site of CAT-1. In addition, we found out that the water molecules have the leading role in forming the transportation channel within CAT-1. The conductive structure of CAT-1 was formed only when the water molecules were continuously distributed across the channel. Steered molecular dynamics (SMD) simulation approach showed various energy barriers against arginine transportation through CAT-1, especially while crossing the bottlenecks of the related channel. These findings at the molecular level might shed light on identifying the crucial amino acids in the binding of arginine to eukaryotic CATs and also provide fundamental insights into the arginine transportation mechanisms through CAT-1. Understanding the transportation mechanism of arginine is essential to developing CAT-1 blockers, which can be potential medications for some types of cancers.Communicated by Ramaswamy H. Sarma.

Structure and Computational Electrophysiology of Ac-LS3, a Synthetic Ion Channel

Computer simulations are reported on Ac-LS3, a synthetic ion channel, containing 21 residues with a Leu-Ser-Ser-Leu-Leu-Ser-Leu heptad repeat, which forms ions channels upon application of voltage. A hexameric, coiled-coil bundle initially positioned perpendicular to the membrane settled into a stable, tilted structure after 1.5 μs, most likely to improve contacts between the non-polar exterior of the channel and the hydrophobic core of the membrane. Once tilted, the bundle remained in this state during subsequent simulations of nearly 10 μs at voltages ranging from 200 to -100 mV. In contrast, attempts to identify a stable pentameric structure failed, thus supporting the hypothesis that the channel is a hexamer. Results at 100 mV were used to reconstruct the free energy profiles for K⁺ and Cl^- in the channel. This was done by way of several methods in which results of molecular dynamics (MD) simulations were combined with the electrodiffusion model. Two of them developed recently do not require knowledge of the diffusivity. Instead, they utilize one-sided density profiles and committor probabilities. The consistency between different methods is very good, supporting the utility of the newly developed methods for reconstructing free energies of ions in channels. The flux of K⁺, which accounts for most of the current through the channel, calculated directly from MD matches well the total measured current. However, the current of Cl^- is somewhat overestimated, possibly due to a slightly unbalanced force field involving chloride. The current-voltage dependence was also reconstructed by way of a recently developed, efficient method that requires simulations only at a single voltage, yielding good agreement with the experiment. Taken together, the results demonstrate that computational electrophysiology has become a reliable tool for studying how channels mediate ion transport through membranes.

When is a hydrophobic gate not a hydrophobic gate?

The flux of ions through a channel is most commonly regulated by changes that result in steric occlusion of its pore. However, ion permeation can also be prevented by formation of a desolvation barrier created by hydrophobic residues that line the pore. As a result of relatively minor structural changes, confined hydrophobic regions in channels may undergo transitions between wet and dry states to gate the pore closed without physical constriction of the permeation pathway. This concept is referred to as hydrophobic gating, and many examples of this process have been demonstrated. However, the term is also now being used in a much broader context that often deviates from its original meaning. In this Viewpoint, we explore the formal definition of a hydrophobic gate, discuss examples of this process compared with other gating mechanisms that simply exploit hydrophobic residues and/or lipids in steric closure of the pore, and describe the best practice for identification of a hydrophobic gate.

General Pharmacological Activation Mechanism of K+ Channels Bypassing Channel Gates

Under the known pharmacological activation mechanisms, activators allosterically or directly open potassium channel gates. However, herein, molecular dynamics simulations on TREK-1, a member of the channel class gated at the filter, suggested that negatively charged activators act with a gate-independent mechanism where compounds increase currents by promoting ions passing through the central cavity. Then, based on studies of KCNQ2, we uncovered that this noncanonical activation mechanism is shared by the other channel class gated at the helix-bundle crossing. Rational drug design found a novel KCNQ2 agonist, CLE030, which stably binds to the central cavity. Functional analysis, molecular dynamics simulations, and calculations of the potential of mean force revealed that the carbonyl oxygen of CLE030 influences permeant ions in the central cavity to contribute to its activation effects. Together, this study discovered a ligand-to-ion activation mechanism for channels that bypasses their gates and thus is conserved across subfamilies with different gates.

Decisive structural elements in water and ion permeation through mechanosensitive channels of large conductance: insights from molecular dynamics simulation

In this paper, a series of equilibrium molecular dynamics simulations (EMD), steered molecular dynamics (SMD), and computational electrophysiology methods are carried out to explore water and ion permeation through mechanosensitive channels of large conductance (MscL). This research aims to identify the pore-lining side chains of the channel in different conformations of MscL homologs by analyzing the pore size. The distribution of permeating water dipole angles through the pore domains enclosed by VAL21 and GLU104 demonstrated that water molecules are oriented toward the charged oxygen headgroups of GLU104 from their hydrogen atoms to retain this interaction in a stabilized fashion. Although, this behavior was not perceived for VAL21. Numerical assessments of the secondary structure clarified that, during the ion permeation, in addition to the secondary structure alterations, the structure of Tb-MscL would also undergo significant conformational changes. It was elucidated that VAL21, GLU104, and water molecules accomplish a fundamental task in ion permeation. The mentioned residues hinder ion permeation so that the pulling SMD force is increased remarkably when the ions permeate through the domains enclosed by VAL21 and GLU102. The hydration level and potassium diffusivity in the hydrophobic gate of the transmembrane domain were promoted by applying the external electric field. Furthermore, the implementation of an external electric field altered the distribution pattern for potassium ions in the system while intensifying the accumulation of Cl⁻ in the vicinity of ARG11 and ARG98.

Graphical representation of the most determinant pore-lining side chains of Tb-MscL along with the solid surfaces depicting the spatial shape of the interior pore.

Conformational transitions and ligand-binding to a muscle-type nicotinic acetylcholine receptor

Fast synaptic communication requires receptors that respond to the presence of neurotransmitter by opening an ion channel across the post-synaptic membrane. The muscle-type nicotinic acetylcholine receptor from the electric fish, Torpedo, is the prototypic ligand-gated ion channel, yet the structural changes underlying channel activation remain undefined. Here we use cryo-EM to solve apo and agonist-bound structures of the Torpedo nicotinic receptor embedded in a lipid nanodisc. Using both a direct biochemical assay to define the conformational landscape and molecular dynamics simulations to assay flux through the pore, we correlate structures with functional states and elucidate the motions that lead to pore activation of a heteromeric nicotinic receptor. We highlight an underappreciated role for the complementary subunit in channel gating, establish the structural basis for the differential agonist affinities of α/δ versus α /γ sites, and explain why nicotine is less potent at muscle nicotinic receptors compared to neuronal ones.

Dynorphin A induces membrane permeabilization by formation of proteolipidic pores. Insights from electrophysiology and computational simulations

Dynorphins are endogenous neuropeptides that function as ligands for the κ-opioid receptor. In addition to opioid activity, dynorphins can induce several pathological effects such as neurological dysfunctions and cell death. Previous studies have suggested that Dynorphin A (DynA) mediates some pathogenic actions through formation of transient pores in lipid domains of the plasma membrane. Here, we use planar bilayer electrophysiology to show that DynA induces pore formation in negatively charged membranes. We find a large variability in pore conformations showing equilibrium conductance fluctuations, what disregards electroporation as the dominant mechanism of pore formation. Ion selectivity measurements showing cationic selectivity indicate that positive protein charges of DynA are stabilized by phosphatidyl serine negative charges in the formation of combined structures. We complement our study with computational simulations that assess the stability of diverse peptide arrangements in the hydrophobic core of the bilayer. We show that DynA is capable of assembling in charged membranes to form water-filled pores that conduct ions.

An open state of a voltage-gated sodium channel involving a π-helix and conserved pore-facing asparagine

Voltage-gated sodium (Nav) channels play critical roles in propagating action potentials and otherwise manipulating ionic gradients in excitable cells. These channels open in response to membrane depolarization, selectively permeating sodium ions until rapidly inactivating. Structural characterization of the gating cycle in this channel family has proved challenging, particularly due to the transient nature of the open state. A structure from the bacterium Magnetococcus marinus Nav (NavMs) was initially proposed to be open, based on its pore diameter and voltage-sensor conformation. However, the functional annotation of this model, and the structural details of the open state, remain disputed. In this work, we used molecular modeling and simulations to test possible open-state models of NavMs. The full-length experimental structure, termed here the α-model, was consistently dehydrated at the activation gate, indicating an inability to conduct ions. Based on a spontaneous transition observed in extended simulations, and sequence/structure comparison to other Nav channels, we built an alternative π-model featuring a helix transition and the rotation of a conserved asparagine residue into the activation gate. Pore hydration, ion permeation, and state-dependent drug binding in this model were consistent with an open functional state. This work thus offers both a functional annotation of the full-length NavMs structure and a detailed model for a stable Nav open state, with potential conservation in diverse ion-channel families.

Water Nanoconfined in a Hydrophobic Pore: Molecular Dynamics Simulations of Transmembrane Protein 175 and the Influence of Water Models

Water molecules within biological ion channels are in a nanoconfined environment and therefore exhibit behaviors which differ from that of bulk water. Here, we investigate the phenomenon of hydrophobic gating, the process by which a nanopore may spontaneously dewet to form a "vapor lock" if the pore is sufficiently hydrophobic and/or narrow. This occurs without steric occlusion of the pore. Using molecular dynamics simulations with both rigid fixed-charge and polarizable (AMOEBA) force fields, we investigate this wetting/dewetting behavior in the transmembrane protein 175 ion channel. We examine how a range of rigid fixed-charge and polarizable water models affect wetting/dewetting in both the wild-type structure and in mutants chosen to cover a range of nanopore radii and pore-lining hydrophobicities. Crucially, we find that the rigid fixed-charge water models lead to similar wetting/dewetting behaviors, but that the polarizable water model resulted in an increased wettability of the hydrophobic gating region of the pore. This has significant implications for molecular simulations of nanoconfined water, as it implies that polarizability may need to be included if we are to gain detailed mechanistic insights into wetting/dewetting processes. These findings are of importance for the design of functionalized biomimetic nanopores (e.g., sensing or desalination) as well as for furthering our understanding of the mechanistic processes underlying biological ion channel function.

Ionic Transport Triggered by Asymmetric Illumination on 2D Nano-Membrane

Ionic transport and ion sieving are important in the field of separation science and engineering. Based on the rapid development of nanomaterials and nano-devices, more and more phenomena occur on the nanoscale devices in the field of thermology, optics, mechanics, etc. Recently, we experimentally observed a novel ion transport phenomenon in nanostructured graphene oxide membrane (GOM) under asymmetric illumination. We first build a light-induced carriers’ diffusion model based on our previous experimental results. This model can reveal the light-induced ion transport mechanism and predict the carriers’ diffusion behavior under different operational situations and material characters. The voltage difference increases with the rise of illuminate asymmetry, photoresponsivity, recombination coefficient, and carriers’ diffusion coefficient ratio. Finally, we discuss the ion transport behavior with different surface charge densities using MD simulation. Moderate surface charge decreases the ion transport with the same type of charge due to the electrostatic repulsion; however, excess surface charge blocks both cation and anion because a thicker electrical double layer decreases effective channel height. Research here provides referenced operational and material conditions to obtain a greater voltage difference between the membrane sides. Also, the mechanism of ion transport and ion sieving can guide us to modify membrane material according to different aims.

High-Performance Biomimetic Water Channel: The Constructive Interplay of Interaction Parameters and Hydrophilic Doping Levels

In this study, we introduce a superfast biomimetic water channel mimicking the hydrophobicity scales of the Aquaporin (AQP) pore lining. Molecular dynamics simulation is used to scrutinize the impact of hydrophilic doping level in the nanotube and the water-wall interaction strength on water permeability. In the designed biomimetic channel, the constructive interplay of Lennard-Jones (LJ) ε parameters and hydrophilic doping levels increased the possibility of ultrafast water transport. Moreover, a unique set of LJ parameters is discovered for each biomimetic channel with different hydrophilic doping levels, enhancing water permeation. Inside high-performance biomimetic channels, water distribution surprisingly implies a varying pore geometry that narrows down in the middle, mimicking the pattern obtained from GplF pore analysis, evoking the narrow pore induced by the aromatic/arginine selectivity filter. This exciting accordance occurred as a result of tailoring specific hydrophilic arrays within the hydrophobic channel backbone by mimicking the AQP pore interior. The main takeaway of hydrophilic doping arrays implanted within the hydrophobic nanotube is to break the large barrier in the water-wall vdW energy profile into multiple reduced ones to increase water conduction. Consequently, the "water jumping" phenomenon in the middle of the biomimetic channel occurs under specific circumstances. The biomimetic channel with the highest value of water permeability of about 13.67 ± 0.66 × 10^-13 cm³·s^-1 exhibits the best mechanism for artificial water channels (AWCs), serving superfast water transport considering the low entrance barrier and weak water-wall interaction.

Hydrophobic dewetting in gating and regulation of transmembrane protein ion channels

Water is at the heart of almost all biological phenomena, without which no life that we know of would have been possible. It is a misleadingly complex liquid that exists in near coexistence with the vapor phase under ambient conditions. Confinement within a hydrophobic cavity can tip this balance enough to drive a cooperative dewetting transition. For a nanometer-scale pore, the dewetting transition leads to a stable dry state that is physically open but impermeable to ions. This phenomenon is often referred to as hydrophobic gating. Numerous transmembrane protein ion channels have now been observed to utilize hydrophobic gating in their activation and regulation. Here, we review recent theoretical, simulation, and experimental studies that together have started to establish the principles of hydrophobic gating and discuss how channels of various sizes, topologies, and biological functions can utilize these principles to control the thermodynamic properties of water within their interior pores for gating and regulation. Exciting opportunities remain in multiple areas, particularly on direct experimental detection of hydrophobic dewetting in biological channels and on understanding how the cell may control the hydrophobic gating in regulation of ion channels.

Electrophysiological Properties from Computations at a Single Voltage: Testing Theory with Stochastic Simulations

We use stochastic simulations to investigate the performance of two recently developed methods for calculating the free energy profiles of ion channels and their electrophysiological properties, such as current–voltage dependence and reversal potential, from molecular dynamics simulations at a single applied voltage. These methods require neither knowledge of the diffusivity nor simulations at multiple voltages, which greatly reduces the computational effort required to probe the electrophysiological properties of ion channels. They can be used to determine the free energy profiles from either forward or backward one-sided properties of ions in the channel, such as ion fluxes, density profiles, committor probabilities, or from their two-sided combination. By generating large sets of stochastic trajectories, which are individually designed to mimic the molecular dynamics crossing statistics of models of channels of trichotoxin, p7 from hepatitis C and a bacterial homolog of the pentameric ligand-gated ion channel, GLIC, we find that the free energy profiles obtained from stochastic simulations corresponding to molecular dynamics simulations of even a modest length are burdened with statistical errors of only 0.3 kcal/mol. Even with many crossing events, applying two-sided formulas substantially reduces statistical errors compared to one-sided formulas. With a properly chosen reference voltage, the current–voltage curves can be reproduced with good accuracy from simulations at a single voltage in a range extending for over 200 mV. If possible, the reference voltages should be chosen not simply to drive a large current in one direction, but to observe crossing events in both directions.

On the Functional Annotation of Open-Channel Structures in the Glycine Receptor

The glycine receptor (GlyR) is by far the best-characterized pentameric ligand-gated ion channel, with several high-resolution structures from X-ray crystallography, cryoelectron microscopy (cryo-EM), and modeling. Nonetheless, the significance of the currently available open-pore conformations is debated due to their diversity in the pore geometry. Here, we discuss the physiological significance of existing models of the GlyR active state based on conductance and selectivity measurements by computational electrophysiology. The results support the conclusion that the original cryo-EM reconstruction of the active state obtained in detergents as well as its subsequent refinement by molecular dynamics simulations are likely to be non-physiological as they feature artificially dilated ion pores. In addition, the calculations indicate that a physiologically relevant open pore should be constricted within a radius of 2.5 and 2.8 Å, which is consistent with previous modeling, electrophysiology measurements, and the most recent cryo-EM structures obtained in a native lipid membrane environment.

Molecular Simulations of Hydrophobic Gating of Pentameric Ligand Gated Ion Channels: Insights into Water and Ions

Ion channels are proteins which form gated nanopores in biological membranes. Many channels exhibit hydrophobic gating, whereby functional closure of a pore occurs by local dewetting. The pentameric ligand gated ion channels (pLGICs) provide a biologically important example of hydrophobic gating. Molecular simulation studies comparing additive vs polarizable models indicate predictions of hydrophobic gating are robust to the model employed. However, polarizable models suggest favorable interactions of hydrophobic pore-lining regions with chloride ions, of relevance to both synthetic carriers and channel proteins. Electrowetting of a closed pLGIC hydrophobic gate requires too high a voltage to occur physiologically but may inform designs for switchable nanopores. Global analysis of ∼200 channels yields a simple heuristic for structure-based prediction of (closed) hydrophobic gates. Simulation-based analysis is shown to provide an aid to interpretation of functional states of new channel structures. These studies indicate the importance of understanding the behavior of water and ions within the nanoconfined environment presented by ion channels.

Effects of amino acid modifications on the permeability of the pentameric sarcolipin channel

Sarcolipin (SLN) is an important transmembrane (TM) protein encoded by long noncoding RNA. SLN is expressed in the sarcoplasmic reticulum and regulates cardiac and skeletal muscle contractions. SLN forms a pentameric hydrophobic ligand-gated ion channel. The protonation of Glu7 (protonated SLN, pSLN) and mutation of Thr18 to Ala18 (T18A) have been reported to exert a significant influence on the permeability of the channel. In this study, the altered permeability of both the pSLN and T18A pentameric channels was simulated. Combined with molecular dynamics simulation, the free-energy landscape for single ions, computational electrophysiology, diffusion coefficient, and pore geometrical characteristic analyses were performed to further understand the properties of amino acid modifications in the SLN pentameric channel. The results suggest that both the pSLN and T18A pentameric channels form stable hydrophobic ligand-gated channels. The TM voltage has a positive effect on the permeability of water molecules and ions. By using pSLN and T18A, our study provides helpful information on the pore-forming mechanism of SLN and furthers our understanding of the regulatory mechanisms underlying the permeation of ions and water molecules in the pentameric SLN channel.

Computational Study of the Ion and Water Permeation and Transport Mechanisms of the SARS-CoV-2 Pentameric E Protein Channel

Coronavirus disease 2019 (COVID-19) is caused by a novel coronavirus (SARS-CoV-2) and represents the causative agent of a potentially fatal disease that is a public health emergency of international concern. Coronaviruses, including SARS-CoV-2, encode an envelope (E) protein, which is a small, hydrophobic membrane protein; the E protein of SARS-CoV-2 shares a high level of homology with severe acute respiratory syndrome coronavirus (SARS-CoV). In this study, we provide insights into the function of the SARS-CoV-2 E protein channel and the ion and water permeation mechanisms using a combination of in silico methods. Based on our results, the pentameric E protein promotes the penetration of cation ions through the channel. An analysis of the potential mean force (PMF), pore radius and diffusion coefficient reveals that Leu10 and Phe19 are the hydrophobic gates of the channel. In addition, the pore exhibits a clear wetting/dewetting transition with cation selectivity under transmembrane voltage, indicating that it is a hydrophobic voltage-dependent channel. Overall, these results provide structure-based insights and molecular dynamic information that are needed to understand the regulatory mechanisms of ion permeability in the pentameric SARS-CoV-2 E protein channel.

Water in Nanopores and Biological Channels: A Molecular Simulation Perspective

This Review explores the dynamic behavior of water within nanopores and biological channels in lipid bilayer membranes. We focus on molecular simulation studies, alongside selected structural and other experimental investigations. Structures of biological nanopores and channels are reviewed, emphasizing those high-resolution crystal structures, which reveal water molecules within the transmembrane pores, which can be used to aid the interpretation of simulation studies. Different levels of molecular simulations of water within nanopores are described, with a focus on molecular dynamics (MD). In particular, models of water for MD simulations are discussed in detail to provide an evaluation of their use in simulations of water in nanopores. Simulation studies of the behavior of water in idealized models of nanopores have revealed aspects of the organization and dynamics of nanoconfined water, including wetting/dewetting in narrow hydrophobic nanopores. A survey of simulation studies in a range of nonbiological nanopores is presented, including carbon nanotubes, synthetic nanopores, model peptide nanopores, track-etched nanopores in polymer membranes, and hydroxylated and functionalized nanoporous silica. These reveal a complex relationship between pore size/geometry, the nature of the pore lining, and rates of water transport. Wider nanopores with hydrophobic linings favor water flow whereas narrower hydrophobic pores may show dewetting. Simulation studies over the past decade of the behavior of water in a range of biological nanopores are described, including porins and β-barrel protein nanopores, aquaporins and related polar solute pores, and a number of different classes of ion channels. Water is shown to play a key role in proton transport in biological channels and in hydrophobic gating of ion channels. An overall picture emerges, whereby the behavior of water in a nanopore may be predicted as a function of its hydrophobicity and radius. This informs our understanding of the functions of diverse channel structures and will aid the design of novel nanopores. Thus, our current level of understanding allows for the design of a nanopore which promotes wetting over dewetting or vice versa. However, to design a novel nanopore, which enables fast, selective, and gated flow of water de novo would remain challenging, suggesting a need for further detailed simulations alongside experimental evaluation of more complex nanopore systems.

The general anesthetic etomidate and fenamate mefenamic acid oppositely affect GABAAR and GlyR: a structural explanation

GABA and glycine act as inhibitory neurotransmitters in the CNS. Inhibitory neurotransmission is mediated via activation of ionotropic GABA_A and glycine receptors. We used a modeling approach to explain the opposite effects of the general anesthetic etomidate (ETM) and fenamate mefenamic acid (MFA) on GABA- and glycine-activated currents recorded in isolated cerebellar Purkinje cells and hippocampal pyramidal neurons, respectively. These drugs potentiated GABA_ARs but blocked GlyRs. We built a homology model of α1β GlyR based on the cryo-EM structure of open α1 GlyR, used the α1β3γ2 GABA_AR structure from the PDB, and applied Monte-Carlo energy minimization to optimize models of receptors and ligand-receptor complexes. In silico docking suggests that ETM/MFA bind at the transmembrane β( +)/α( -) intersubunit interface in GABA_AR. Our models predict that the bulky side chain of the highly conserved Arg^19' residue at the plus interface side wedges the interface and maintains the conducting receptor state. We hypothesized that MFA/ETM binding at the β( +)/α( -) interface leads to prolongation of receptor life-time in the open state. Having analyzed different GABA_AR and GlyR structures available in the PDB, we found that mutual arrangement of the Arg^19' and Gln^-26' side chains at the plus and minus interface sides, respectively, plays an important role when the receptor switches from the open to closed state. We show that this process is accompanied by narrowing of the intersubunit interfaces, leading to extrusion of the Arg^19' side chain from the interface. Our models allow us to explain the lack of GlyR potentiation in our electrophysiological experiments.

Thermodynamics and Kinetics of Ion Permeation in Wild-Type and Mutated Open Active Conformation of the Human α7 Nicotinic Receptor

Molecular studies of human pentameric ligand-gated ion channels (LGICs) expressed in neurons and at neuromuscular junctions are of utmost importance in the development of therapeutic strategies for neurological disorders. We focus here on the nicotinic acetylcholine receptor nAChR-α7, a homopentameric channel widely expressed in the human brain, with a proven role in a wide spectrum of disorders including schizophrenia and Alzheimer’s disease. By exploiting an all-atom structural model of the full (transmembrane and extracellular) protein in the open, agonist-bound conformation we recently developed, we evaluate the free energy and the mean first passage time of single-ion permeation using molecular dynamics simulations and the milestoning method with Voronoi tessellation. The results for the wild-type channel provide the first available mapping of the potential of mean force in the full-length α7 nAChR, reveal its expected cationic nature, and are in good agreement with simulation data for other channels of the LGIC family and with experimental data on nAChRs. We then investigate the role of a specific mutation directly related to ion selectivity in LGICs, the E-1′ → A-1′ substitution at the cytoplasmatic selectivity filter. We find that the mutation strongly affects sodium and chloride permeation in opposite directions, leading to a complete inversion of selectivity, at variance with the limited experimental results available that classify this mutant as cationic. We thus provide structural determinants for the observed cationic-to-anionic inversion, revealing a key role of the protonation state of residue rings far from the mutation, in the proximity of the hydrophobic channel gate.

Electric Field Induced Wetting of a Hydrophobic Gate in a Model Nanopore Based on the 5-HT3 Receptor Channel

In this study we examined the influence of a transmembrane voltage on the hydrophobic gating of nanopores using molecular dynamics simulations. We observed electric field induced wetting of a hydrophobic gate in a biologically inspired model nanopore based on the 5-HT3 receptor in its closed state, with a field of at least ∼100 mV nm-1 (corresponding to a supra-physiological potential difference of ∼0.85 V across the membrane) required to hydrate the pore. We also found an unequal distribution of charged residues can generate an electric field intrinsic to the nanopore which, depending on its orientation, can alter the effect of the external field, thus making the wetting response asymmetric. This wetting response could be described by a simple model based on water surface tension, the volumetric energy contribution of the electric field, and the influence of charged amino acids lining the pore. Finally, the electric field response was used to determine time constants characterizing the phase transitions of water confined within the nanopore, revealing liquid-vapor oscillations on a time scale of ∼5 ns. This time scale was largely independent of the water model employed and was similar for different sized pores representative of the open and closed states of the pore. Furthermore, our finding that the threshold voltage required for hydrating a hydrophobic gate depends on the orientation of the electric field provides an attractive perspective for the design of rectifying artificial nanopores.

Mechanism of ligand activation of a eukaryotic cyclic nucleotide-gated channel

Cyclic nucleotide-gated (CNG) channels convert cyclic nucleotide (CN) binding and unbinding into electrical signals in sensory receptors and neurons. The molecular conformational changes underpinning ligand activation are largely undefined. We report both closed- and open-state atomic cryo-EM structures of a full-length Caenorhabditis elegans cyclic GMP-activated channel TAX-4, reconstituted in lipid nanodiscs. These structures, together with computational and functional analyses and a mutant channel structure, reveal a double-barrier hydrophobic gate formed by two S6 amino acids in the central cavity. cGMP binding produces global conformational changes that open the cavity gate located ~52 Å away but do not alter the structure of the selectivity filter-the commonly presumed activation gate. Our work provides mechanistic insights into the allosteric gating and regulation of CN-gated and nucleotide-modulated channels and CNG channel-related channelopathies.

Induced Polarization in Molecular Dynamics Simulations of the 5-HT3 Receptor Channel

Ion channel proteins form water-filled nanoscale pores within lipid bilayers, and their properties are dependent on the complex behavior of water in a nanoconfined environment. Using a simplified model of the pore of the 5-HT3 receptor (5HT3R) which restrains the backbone structure to that of the parent channel protein from which it is derived, we compare additive with polarizable models in describing the behavior of water in nanopores. Molecular dynamics simulations were performed with four conformations of the channel: two closed state structures, an intermediate state, and an open state, each embedded in a phosphatidylcholine bilayer. Water density profiles revealed that for all water models, the closed and intermediate states exhibited strong dewetting within the central hydrophobic gate region of the pore. However, the open state conformation exhibited varying degrees of hydration, ranging from partial wetting for the TIP4P/2005 water model to complete wetting for the polarizable AMOEBA14 model. Water dipole moments calculated using polarizable force fields also revealed that water molecules remaining within dewetted sections of the pore resemble gas phase water. Free energy profiles for Na+ and for Cl- ions within the open state pore revealed more rugged energy landscapes using polarizable force fields, and the hydration number profiles of these ions were also sensitive to induced polarization resulting in a substantive reduction of the number of waters within the first hydration shell of Cl- while it permeates the pore. These results demonstrate that induced polarization can influence the complex behavior of water and ions within nanoscale pores and provides important new insights into their chemical properties.

A Refined Open State of the Glycine Receptor Obtained via Molecular Dynamics Simulations

Pentameric ligand-gated ion channels are key players in mediating fast neurotransmission. Glycine receptors are chloride-selective members of this receptor family that mediate inhibitory synaptic transmission and are implicated in neurological disorders including autism and hyperekplexia. They have been structurally characterized by both X-ray crystallography and cryoelectron microscopy (cryo-EM) studies, with the latter giving rise to what was proposed as a possible open state. However, recent work has questioned the physiological relevance of this open state structure, since it rapidly collapses in molecular dynamics simulations. Here, we show that the collapse can be avoided by a careful equilibration protocol that reconciles the more problematic regions of the original density map and gives a stable open state that shows frequent selective chloride permeation. The protocol developed in this work provides a means to refine open-like structures of the whole pentameric ligand-gated ion channel superfamily and reconciles the previous issues with the cryo-EM structure.

Enhanced wettability of long narrow carbon nanotubes in a double-walled hetero-structure: unraveling the effects of a boron nitride nanotube as the exterior

Studying the structure and dynamics of nano-confined water inside carbon nanotubes has consistently attracted the wide-spread interest of researchers. In the present work, molecular dynamics simulations indicated internal nonwetting behavior for the central region of the long and narrow single-wall carbon nanotube (5,5) (SWNT) and showed that continuous single-file water molecules are not formed through it. Unlike the SWNT, by adding boron nitride nanotubes (6,6) as an outer wall to the SWNT, a continuously long single-file water chain is formed through the double-walled carbon and boron nitride hetero-nanotube (DWHNT) and thorough internal wetting of the DWHNT is observed. The position and the number of water molecules, electrostatic potential heatmap of the nanotube's wall, free energy profile of nano-confined water, and number of hydrogen bonds between them confirmed the aforementioned results and complete internal wetting of the DWHNT. After using the boron nitride nanotube (6,6) as the outer wall, an homogeneous electrostatic potential distribution in the DWHNT and increase in the hydrophilic characteristics of the nano-channel wall are observed, bringing about gradual trapping of more water molecules through it. Finally, water molecules occupied the central region of the DWHNT and a thorough single-file water chain is formed inside the nano-channel. Water dipole orientation inside the DWHNT and their radial distribution function asserted the occurrence of the liquid-solid quasi-phase transition of single-file water molecules confined inside the long and narrow carbon nanotube (5,5) under ambient conditions.

Unique structure and function of viral rhodopsins

Recently, two groups of rhodopsin genes were identified in large double-stranded DNA viruses. The structure and function of viral rhodopsins are unknown. We present functional characterization and high-resolution structure of an Organic Lake Phycodnavirus rhodopsin II (OLPVRII) of group 2. It forms a pentamer, with a symmetrical, bottle-like central channel with the narrow vestibule in the cytoplasmic part covered by a ring of 5 arginines, whereas 5 phenylalanines form a hydrophobic barrier in its exit. The proton donor E42 is placed in the helix B. The structure is unique among the known rhodopsins. Structural and functional data and molecular dynamics suggest that OLPVRII might be a light-gated pentameric ion channel analogous to pentameric ligand-gated ion channels, however, future patch clamp experiments should prove this directly. The data shed light on a fundamentally distinct branch of rhodopsins and may contribute to the understanding of virus-host interactions in ecologically important marine protists.

Conduction through a narrow inward-rectifier K+ channel pore

G-protein–gated inwardly rectifying potassium channels are important mediators of inhibitory neurotransmission. Based on microsecond-scale molecular dynamics simulations, Bernsteiner et al. propose novel gating details that may enable K⁺ flux via a direct knock-on mechanism.

Inwardly rectifying potassium (Kir) channels play a key role in controlling membrane potentials in excitable and unexcitable cells, thereby regulating a plethora of physiological processes. G-protein–gated Kir channels control heart rate and neuronal excitability via small hyperpolarizing outward K⁺ currents near the resting membrane potential. Despite recent breakthroughs in x-ray crystallography and cryo-EM, the gating and conduction mechanisms of these channels are poorly understood. MD simulations have provided unprecedented details concerning the gating and conduction mechanisms of voltage-gated K⁺ and Na⁺ channels. Here, we use multi-microsecond–timescale MD simulations based on the crystal structures of GIRK2 (Kir3.2) bound to phosphatidylinositol-4,5-bisphosphate to provide detailed insights into the channel’s gating dynamics, including insights into the behavior of the G-loop gate. The simulations also elucidate the elementary steps that underlie the movement of K⁺ ions through an inward-rectifier K⁺ channel under an applied electric field. Our simulations suggest that K⁺ permeation might occur via direct knock-on, similar to the mechanism recently shown for K_v channels.

Water and hydrophobic gates in ion channels and nanopores

Simulations of water behaviour have been used to probe hydrophobic gates in BEST1 and TMEM175, which can reveal important design principles for the engineering of gates in novel biomimetic nanopores.

Ion channel proteins form nanopores in biological membranes which allow the passage of ions and water molecules. Hydrophobic constrictions in such pores can form gates, i.e. energetic barriers to water and ion permeation. Molecular dynamics simulations of water in ion channels may be used to assess whether a hydrophobic gate is closed (i.e. impermeable to ions) or open. If there is an energetic barrier to water permeation then it is likely that a gate will also be impermeable to ions. Simulations of water behaviour have been used to probe hydrophobic gates in two recently reported ion channel structures: BEST1 and TMEM175. In each of these channels a narrow region is formed by three consecutive rings of hydrophobic sidechains and in both cases such analysis demonstrates that the crystal structures correspond to a closed state of the channel. In silico mutations of BEST1 have also been used to explore the effect of changes in the hydrophobicity of the gating constriction, demonstrating that substitution of hydrophobic sidechains with more polar sidechains results in an open gate which allows water permeation. A possible open state of the TMEM175 channel was modelled by the in silico expansion of the hydrophobic gate resulting in the wetting of the pore and free permeation of potassium ions through the channel. Finally, a preliminary study suggests that a hydrophobic gate motif can be transplanted in silico from the BEST1 channel into a simple β-barrel pore template. Overall, these results suggest that simulations of the behaviour of water in hydrophobic gates can reveal important design principles for the engineering of gates in novel biomimetic nanopores.

A heuristic derived from analysis of the ion channel structural proteome permits the rapid identification of hydrophobic gates

Ion channels are nanoscale protein pores in cell membranes. An exponentially increasing number of structures for channels means that computational methods for predicting their functional state are needed. Hydrophobic gates in ion channels result in local dewetting of pores, which functionally closes them to water and ion permeation. We use simulations of water behavior within nearly 200 different ion channel structures to explore how the radius and hydrophobicity of pores determine their hydration vs. dewetting behavior. Machine learning-assisted analysis of these simulations allowed us to propose a simple model for this relationship and present an easy method for rapidly predicting the functional state of new channel structures as they emerge.

Ion channel proteins control ionic flux across biological membranes through conformational changes in their transmembrane pores. An exponentially increasing number of channel structures captured in different conformational states are now being determined; however, these newly resolved structures are commonly classified as either open or closed based solely on the physical dimensions of their pore, and it is now known that more accurate annotation of their conductive state requires additional assessment of the effect of pore hydrophobicity. A narrow hydrophobic gate region may disfavor liquid-phase water, leading to local dewetting, which will form an energetic barrier to water and ion permeation without steric occlusion of the pore. Here we quantify the combined influence of radius and hydrophobicity on pore dewetting by applying molecular dynamics simulations and machine learning to nearly 200 ion channel structures. This allows us to propose a simple simulation-free heuristic model that rapidly and accurately predicts the presence of hydrophobic gates. This not only enables the functional annotation of new channel structures as soon as they are determined, but also may facilitate the design of novel nanopores controlled by hydrophobic gates.

CHAP: A Versatile Tool for the Structural and Functional Annotation of Ion Channel Pores

The control of ion channel permeation requires the modulation of energetic barriers or “gates” within their pores. However, such barriers are often simply identified from the physical dimensions of the pore. Such approaches have worked well in the past, but there is now evidence that the unusual behavior of water within narrow hydrophobic pores can produce an energetic barrier to permeation without requiring steric occlusion of the pathway. Many different ion channels have now been shown to exploit “hydrophobic gating” to regulate ion flow, and it is clear that new tools are required for more accurate functional annotation of the increasing number of ion channel structures becoming available. We have previously shown how molecular dynamics simulations of water can be used as a proxy to predict hydrophobic gates, and we now present a new and highly versatile computational tool, the Channel Annotation Package (CHAP) that implements this methodology.

•

Some ion channels exhibit hydrophobic gating via dewetting of the central pore.

•

This cannot be predicted from their pore radius alone.

•

The hydrophobicity of a pore contributes to this dewetting effect.

•

We have developed a new tool (CHAP) that combines these measurements.

•

CHAP now enables the rapid functional annotation of ion channel structures.

A lipid gating mechanism for the channel-forming O antigen ABC transporter

Extracellular glycan biosynthesis is a widespread microbial protection mechanism. In Gram-negative bacteria, the O antigen polysaccharide represents the variable region of outer membrane lipopolysaccharides. Fully assembled lipid-linked O antigens are translocated across the inner membrane by the WzmWzt ABC transporter for ligation to the lipopolysaccharide core, with the transporter forming a continuous transmembrane channel in a nucleotide-free state. Here, we report its structure in an ATP-bound conformation. Large structural changes within the nucleotide-binding and transmembrane regions push conserved hydrophobic residues at the substrate entry site towards the periplasm and provide a model for polysaccharide translocation. With ATP bound, the transporter forms a large transmembrane channel with openings toward the membrane and periplasm. The channel's periplasmic exit is sealed by detergent molecules that block solvent permeation. Molecular dynamics simulation data suggest that, in a biological membrane, lipid molecules occupy this periplasmic exit and prevent water flux in the transporter's resting state.

Closed-Locked and Apo-Resting State Structures of the Human α7 Nicotinic Receptor: A Computational Study

Nicotinic acetylcholine receptors, belonging to the Cys-loop superfamily of ligand-gated ion channels (LGICs), are membrane proteins present in neurons and at neuromuscular junctions. They are responsible for signal transmission, and their function is regulated by neurotransmitters, agonists, and antagonists drugs. A detailed knowledge of their conformational transition in response to ligand binding is critical to understanding the basis of ligand-receptor interaction, in view of new pharmacological approaches to control receptor activity. However, the scarcity of experimentally derived structures of human channels makes this perspective extremely challenging. To contribute overcoming this issue, we have recently reported structural models for the open and the desensitized states of the human α7 nicotinic receptor. Here, we provide all-atom structural models of the same receptor in two different nonconductive states. The first structure, built via homology modeling and relaxed with extensive Molecular Dynamics simulations, represents the receptor bound to the natural antagonist α-conotoxin ImI. After comparison with available experimental data and computational models of other eukaryotic LGICs, we deem it consistent with the "closed-locked" state. The second model, obtained with simulations from the spontaneous relaxation of the open, agonist-bound α7 structure after ligand removal, recapitulates the characteristics of the apo-resting state of the receptor. These results add to our previous work on the active and desensitized state conformations, contributing to the structural characterization of the conformational landscape of the human α7 receptor and suggesting benchmarks to discriminate among conformations found in experiments or in simulations of LGICs. In particular key interactions at the interface between the extracellular domain and the transmembrane domain are identified, that could be critical to the α7 receptor function.

Cryo-EM reveals two distinct serotonin-bound conformations of full-length 5-HT3A receptor

Serotonin receptor (5-HT_3AR)¹, a cationic pentameric ligand-gated ion channel (pLGIC), is the clinical target for management of nausea and vomiting associated with radiation and chemotherapies². Upon binding, serotonin induces a global conformational change encompassing the ligand-binding extracellular domain (ECD), the transmembrane domain (TMD), and the intracellular domain (ICD), the molecular details of which are unclear. Here, we present two serotonin-bound structures of the full-length 5-HT_3AR in distinct conformations at 3.32 Å and 3.89 Å resolutions that reveal the mechanism underlying channel activation. When compared to Apo-5-HT_3AR, serotonin-bound states underwent a large twisting motion in the ECD and TMD leading to the opening of a 165 Å long permeation pathway. Notably, this motion results in creation of lateral portals for ion permeation at the interface of the TMD and ICD. Combined with molecular dynamics simulations, these structures provide novel insights into conformational coupling across domains and functional modulation.

Conformational transitions of the serotonin 5-HT3 receptor

The serotonin 5-HT₃ receptor is a pentameric ligand-gated ion channel (pLGIC). It belongs to a large family of receptors that function as allosteric signal transducers across the plasma membrane^1,2; upon binding of neurotransmitter molecules to extracellular sites, the receptors undergo complex conformational transitions that result in transient opening of a pore permeable to ions. 5-HT₃ receptors are therapeutic targets for emesis and nausea, irritable bowel syndrome and depression³. In spite of several reported pLGIC structures^4-8, no clear unifying view has emerged on the conformational transitions involved in channel gating. Here we report four cryo-electron microscopy structures of the full-length mouse 5-HT₃ receptor in complex with the anti-emetic drug tropisetron, with serotonin, and with serotonin and a positive allosteric modulator, at resolutions ranging from 3.2 Å to 4.5 Å. The tropisetron-bound structure resembles those obtained with an inhibitory nanobody⁵ or without ligand⁹. The other structures include an 'open' state and two ligand-bound states. We present computational insights into the dynamics of the structures, their pore hydration and free-energy profiles, and characterize movements at the gate level and cation accessibility in the pore. Together, these data deepen our understanding of the gating mechanism of pLGICs and capture ligand binding in unprecedented detail.

Hydrophobic gating in BK channels

The gating mechanism of transmembrane ion channels is crucial for understanding how these proteins control ion flow across membranes in various physiological processes. Big potassium (BK) channels are particularly interesting with large single-channel conductance and dual regulation by membrane voltage and intracellular Ca²⁺. Recent atomistic structures of BK channels failed to identify structural features that could physically block the ion flow in the closed state. Here, we show that gating of BK channels does not seem to require a physical gate. Instead, changes in the pore shape and surface hydrophobicity in the Ca²⁺-free state allow the channel to readily undergo hydrophobic dewetting transitions, giving rise to a large free energy barrier for K⁺ permeation. Importantly, the dry pore remains physically open and is readily accessible to quaternary ammonium channel blockers. The hydrophobic gating mechanism is also consistent with scanning mutagenesis studies showing that modulation of pore hydrophobicity is correlated with activation properties.

BK channels are regulated by membrane voltage and intracellular Ca2+ but the structural features that block the ion flow in the closed state remain unknown. Here authors use molecular dynamics simulation and show that a physical gate is not required; instead ion flow is regulated by hydrophobic dewetting due to changes in pore shape and surface hydrophobicity.

The dual-gate model for pentameric ligand-gated ion channels activation and desensitization

Pentameric ligand-gated ion channels (pLGICs) mediate fast neurotransmission in the nervous system. Their dysfunction is associated with psychiatric, neurological and neurodegenerative disorders such as schizophrenia, epilepsy and Alzheimer's disease. Understanding their biophysical and pharmacological properties, at both the functional and the structural level, thus holds many therapeutic promises. In addition to their agonist-elicited activation, most pLGICs display another key allosteric property, namely desensitization, in which they enter a shut state refractory to activation upon sustained agonist binding. While the activation mechanisms of several pLGICs have been revealed at near-atomic resolution, the structural foundation of desensitization has long remained elusive. Recent structural and functional data now suggest that the activation and desensitization gates are distinct, and are located at both sides of the ion channel. Such a 'dual gate mechanism' accounts for the marked allosteric effects of channel blockers, a feature illustrated herein by theoretical kinetics simulations. Comparison with other classes of ligand- and voltage-gated ion channels shows that this dual gate mechanism emerges as a common theme for the desensitization and inactivation properties of structurally unrelated ion channels.

Permeating disciplines: Overcoming barriers between molecular simulations and classical structure-function approaches in biological ion transport

Ion translocation across biological barriers is a fundamental requirement for life. In many cases, controlling this process-for example with neuroactive drugs-demands an understanding of rapid and reversible structural changes in membrane-embedded proteins, including ion channels and transporters. Classical approaches to electrophysiology and structural biology have provided valuable insights into several such proteins over macroscopic, often discontinuous scales of space and time. Integrating these observations into meaningful mechanistic models now relies increasingly on computational methods, particularly molecular dynamics simulations, while surfacing important challenges in data management and conceptual alignment. Here, we seek to provide contemporary context, concrete examples, and a look to the future for bridging disciplinary gaps in biological ion transport. This article is part of a Special Issue entitled: Beyond the Structure-Function Horizon of Membrane Proteins edited by Ute Hellmich, Rupak Doshi and Benjamin McIlwain.

Conformational landscapes of membrane proteins delineated by enhanced sampling molecular dynamics simulations

The expansion of computational power, better parameterization of force fields, and the development of novel algorithms to enhance the sampling of the free energy landscapes of proteins have allowed molecular dynamics (MD) simulations to become an indispensable tool to understand the function of biomolecules. The temporal and spatial resolution of MD simulations allows for the study of a vast number of processes of interest. Here, we review the computational efforts to uncover the conformational free energy landscapes of a subset of membrane proteins: ion channels, transporters and G-protein coupled receptors. We focus on the various enhanced sampling techniques used to study these questions, how the conclusions come together to build a coherent picture, and the relationship between simulation outcomes and experimental observables.

Un-gating and allosteric modulation of a pentameric ligand-gated ion channel captured by molecular dynamics

Pentameric ligand-gated ion channels (pLGICs) mediate intercellular communication at synapses through the opening of an ion pore in response to the binding of a neurotransmitter. Despite the increasing availability of high-resolution structures of pLGICs, a detailed understanding of the functional isomerization from closed to open (gating) and back is currently missing. Here, we provide the first atomistic description of the transition from open to closed (un-gating) in the glutamate-gated chloride channel (GluCl) from Caenorhabditis Elegans. Starting with the active-state structure solved in complex with the neurotransmitter L-glutamate and the positive allosteric modulator (PAM) ivermectin, we analyze the spontaneous relaxation of the channel upon removal of ivermectin by explicit solvent/membrane Molecular Dynamics (MD) simulations. The μs-long trajectories support the conclusion that ion-channel deactivation is mediated by two distinct quaternary transitions, i.e. a global receptor twisting followed by the radial expansion (or blooming) of the extracellular domain. At variance with previous models, we show that pore closing is exclusively regulated by the global twisting, which controls the position of the β1-β2 loop relative to the M2-M3 loop at the EC/TM domain interface. Additional simulations with L-glutamate restrained to the crystallographic binding mode and ivermectin removed indicate that the same twisting isomerization is regulated by agonist binding at the orthosteric site. These results provide a structural model for gating in pLGICs and suggest a plausible mechanism for the pharmacological action of PAMs in this neurotransmitter receptor family. The simulated un-gating converges to the X-ray structure of GluCl resting state both globally and locally, demonstrating the predictive character of state-of-art MD simulations.

A BEST example of channel structure annotation by molecular simulation

An increasing number of ion channel structures are being determined. This generates a need for computational tools to enable functional annotation of channel structures. However, several studies of ion channel and model pores have indicated that the physical dimensions of a pore are not always a reliable indicator of its conductive status. This is due to the unusual behavior of water within nano-confined spaces, resulting in a phenomenon referred to as “hydrophobic gating”. We have recently demonstrated how simulating the behavior of water within an ion channel pore can be used to predict its conductive status. In this addendum to our study, we apply this method to compare the recently solved structure of a mutant of the bestrophin chloride channel BEST1 with that of the wild-type channel. Our results support the hypothesis of a hydrophobic gate within the narrow neck of BEST1. This provides further validation that this simulation approach provides the basis for an accurate and computationally efficient tool for the functional annotation of ion channel structures.

Voltage Gating of a Biomimetic Nanopore: Electrowetting of a Hydrophobic Barrier

It is desirable that nanopores that are components of biosensors are gated, i.e., capable of controllable switching between closed (impermeable) and open (permeable) states. A central hydrophobic barrier within a nanopore may act as a voltage-dependent gate via electrowetting, i.e., changes in nanopore surface wettability by application of an electric field. We use "computational electrophysiology" simulations to demonstrate and characterize electrowetting of a biomimetic nanopore containing a hydrophobic gate. We show that a hydrophobic gate in a model β-barrel nanopore can be functionally opened by electrowetting at voltages that do not electroporate lipid bilayers. During the process of electrowetting, voltage-induced alignment of water dipoles occurs within the hydrophobic gate region of the nanopore, with water entry preceding permeation of ions through the opened nanopore. When the ionic imbalance that generates a transbilayer potential is dissipated, water is expelled from the hydrophobic gate and the nanopore recloses. The open nanopore formed by electrowetting of a "featureless" β-barrel is anionic selective due to the transmembrane dipole potential resulting from binding of Na⁺ ions to the headgroup regions of the surrounding lipid bilayer. Thus, hydrophobic barriers can provide voltage-dependent gates in designed biomimetic nanopores. This extends our understanding of hydrophobic gating in synthetic and biological nanopores, providing a framework for the design of functional nanopores with tailored gating functionality.