Pentameric ligand-gated ion channels control synaptic neurotransmission by converting chemical signals into electrical signals. Agonist binding leads to rapid signal transduction via an allosteric mechanism, where global protein conformational changes open a pore across the nerve cell membrane. We use all-atom molecular dynamics with a swarm-based string method to solve for the minimum free-energy gating pathways of the proton-activated bacterial GLIC channel. We describe stable wetted/open and dewetted/closed states, and uncover conformational changes in the agonist-binding extracellular domain, ion-conducting transmembrane domain, and gating interface that control communication between these domains. Transition analysis is used to compute free-energy surfaces that suggest allosteric pathways; stabilization with pH; and intermediates, including states that facilitate channel closing in the presence of an agonist. We describe a switching mechanism that senses proton binding by marked reorganization of subunit interface, altering the packing of β-sheets to induce changes that lead to asynchronous pore-lining M2 helix movements. These results provide molecular details of GLIC gating and insight into the allosteric mechanisms for the superfamily of pentameric ligand-gated channels.

A multitude of biological processes requires the participation of specific cations, such as H + , Na + , K + , Ca 2+ , and Mg 2+ . Many of these processes can take place only when proteins have the ability to discriminate between different ions with a very high fidelity. How this is possible is a fundamental question that has fascinated scientists for a long time. At the most fundamental level, it is anticipated that ion selectivity must result from a delicate balance of strong interactions. Yet, identifying and quantifying the key microscopic factors is difficult, as many of these cannot be directly measured by experiments. Theory and computations can contribute by providing a virtual route to complement the missing information. Because ion selectivity is often dominated by thermodynamic factors, detailed molecular dynamics (MD) free energy simulations become important tools. This was vividly illustrated early on with studies of ion solvation (Straatsma and Berendsen, 1988;Åqvist, 1990) and ion-selective systems (Lybrand et al., 1986;Grootenhuis and Kollman, 1989;Åqvist, 1992). These pioneering studies inspired our own efforts.In this Perspective, we aim to present our understanding of ion selectivity as it has evolved over approximately 15 years from studies based on various specific structures: gramicidin A channels (e.g.

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

Excitatory amino acid transporters (EAATs) remove glutamate from synapses. They maintain an efficient synaptic transmission and prevent glutamate from reaching neurotoxic levels. Glutamate transporters couple the uptake of one glutamate to the cotransport of three sodium ions and one proton and the countertransport of one potassium ion. The molecular mechanism for this coupled uptake of glutamate and its co-and counter-transported ions is not known. In a crystal structure of the bacterial glutamate transporter homolog, Glt Ph , only two cations are bound to the transporter, and there is no indication of the location of the third sodium site. In experiments using voltage clamp fluorometry and simulations based on molecular dynamics combined with grand canonical Monte Carlo and free energy simulations performed on different isoforms of Glt Ph as well on a homology model of EAAT3, we sought to locate the third sodium-binding site in EAAT3. Both experiments and computer simulations suggest that T370 and N451 (T314 and N401 in Glt Ph ) form part of the third sodium-binding site. Interestingly, the sodium bound at T370 forms part of the binding site for the amino acid substrate, perhaps explaining both the strict coupling of sodium transport to uptake of glutamate and the ion selectivity of the affinity for the transported amino acid in EAATs.excitatory amino acid transporters | fluorescence | the sodium/aspartate symporter from Pyrococcus horikoshii (Glt Ph ) | simulations G lutamate, the main excitatory neurotransmitter in the central nervous system, is removed from the extracellular synaptic space by the glutamate excitatory amino acid transporters EAAT1-5 (1, 2). These transporters thereby maintain an efficient synaptic communication between neurons and prevent extracellular glutamate from reaching neurotoxic levels (1, 2). EAATs are trimeric proteins in which each subunit functions as an independent transporter (3, 4). Each subunit has eight transmembrane domains and two membrane inserted hair-pin (HP) loops (5-7). EAATs use the Na + and K + gradients in driving the uptake of glutamate against a concentration gradient (8). The uptake of one glutamate is coupled to the cotransport of three Na + ions and one H + ion and the countertransport of one K + ion (9). How the thermodynamically coupled transport of glutamate, H + , Na + , and K + ions is accomplished by glutamate transporters is not known. Here we present evidence for a binding site for Na + ions that suggests a mechanism for the coupling of sodium and glutamate transport.At least one extracellular Na + ion appears to bind before glutamate can bind, and at least one extracellular Na + ion appears to bind after glutamate has bound to the transporter (10-12). For example, in the absence of glutamate, a fluorophore attached to a cysteine at position A430 on HP2 in EAAT3 reported voltageand Na + -dependent fluorescence changes, consistent with Na + binding to the glutamate-free transporter and inducing a conformational change in HP2. Li + also supports glutamat...

Methods directly evaluating the hydration structure and thermodynamics of physiologically relevant cations (Na(+), K(+), Cl(-), etc.) have wide ranging applications in the fields of inorganic, physical, and biological chemistry. All-atom simulations based on accurate potential energy surfaces appear to offer a viable option for assessing the chemistry of ion solvation. Although MD and free energy simulations of ion solvation with classical force fields have proven their usefulness, a number of challenges still remain. One of them is the difficulty of force field benchmarking and validation against structural and thermodynamic data obtained for a condensed phase. Hybrid quantum mechanical/molecular mechanical (QM/MM) models combined with sampling algorithms have the potential to provide an accurate solvation model and to incorporate the effects from the surrounding, which is often missing in gas-phase ab initio computations. Herein, we report the results from QM/MM free energy simulations of Na(+)/K(+) and Cl(-)/Br(-) hydration where we simultaneously characterized the relative thermodynamics of ion solvation and changes in the solvation structure. The Flexible Inner Region Ensemble Separator (FIRES) method was used to impose a spatial separation between QM region and the outer sphere of solvent molecules treated with the CHARMM27 force field. FEP calculations based on QM/MM simulations utilizing the CHARMM/deMon2k interface were performed with different basis set combinations for K(+)/Na(+) and Cl(-)/Br(-) perturbations to establish the dependence of the computed free energies on the basis set level. The dependence of the computed relative free energies on the size of the QM and MM regions is discussed. The current methodology offers an accurate description of structural and thermodynamic aspects of the hydration of alkali and halide ions in neat solvents and can be used to obtain thermodynamic data on ion solvation in condensed phase along with underlying structural properties of the ion-solvent system.

Many enzymes require a specific monovalent cation (M(+)), that is either Na(+) or K(+), for optimal activity. While high selectivity M(+) sites in transport proteins have been extensively studied, enzyme M(+) binding sites generally have lower selectivity and are less characterized. Here we study the M(+) binding site of the model enzyme E. coli β-galactosidase, which is about 10 fold selective for Na(+) over K(+). Combining data from X-ray crystallography and computational models, we find the electrostatic environment predominates in defining the Na(+) selectivity. In this lower selectivity site rather subtle influences on the electrostatic environment become significant, including the induced polarization effects of the M(+) on the coordinating ligands and the effect of second coordination shell residues on the charge distribution of the primary ligands. This work expands the knowledge of ion selectivity in proteins to denote novel mechanisms important for the selectivity of M(+) sites in enzymes.

No abstract

We present a new QM/MM interface for fast and efficient simulations of organic and biological molecules. The CHARMM/deMon interface has been developed and tested to perform minimization and atomistic simulations for multi-particle systems. The current features of this QM/MM interface include readability for molecular dynamics, tested compatibility with Free Energy Perturbation simulations (FEP) using the dual topology/single coordinate method. The current coupling scheme uses link atoms, but further extensions of the code to incorporate other available schemes are planned. We report the performance of different levels of theory for the treatment of the QM region, while the MM region was represented by a classical force-field (CHARMM27) or a polarizable force-field based on a simple Drude model. The current QM/MM implementation can be coupled to the dual-thermostat method and the VV2 integrator to run molecular dynamics simulations.