Liquid-vapor oscillations of water in hydrophobic nanopores

Water plays a key role in biological membrane transport. In ion channels and water-conducting pores (aquaporins), one dimensional confinement in conjunction with strong surface effects changes the physical behavior of water. In molecular dynamics simulations of water in short (0.8 nm) hydrophobic pores the water density in the pore fluctuates on a nanosecond time scale. In long simulations (460 ns in total) at pore radii ranging from 0.35 nm to 1.0 nm we quantify the kinetics of oscillations between a liquid-filled and a vapor-filled pore. This behavior can be explained as capillary evaporation alternating with capillary condensation, driven by pressure fluctuations in the water outside the pore. The free energy difference between the two states depends linearly on the radius. The free energy landscape shows how a metastable liquid state gradually develops with increasing radius. For radii larger than ca. 0.55 nm it becomes the globally stable state and the vapor state vanishes. One dimensional confinement affects the dynamic behavior of the water molecules and increases the self diffusion by a factor of two to three compared to bulk water. Permeabilities for the narrow pores are of the same order of magnitude as for biological water pores. Water flow is not continuous but occurs in bursts. Our results suggest that simulations aimed at collective phenomena such as hydrophobic effects may require simulation times longer than 50 ns. For water in confined geometries, it is not possible to extrapolate from bulk or short time behavior to longer time scales.Comment: 20 pages, 4 figures, 3 tables; to be published in Proc. Natl. Acad. Sci. US

Intermittent permeation of cylindrical nanopores by water

Molecular Dynamics simulations of water molecules in nanometre sized cylindrical channels connecting two reservoirs show that the permeation of water is very sensitive to the channel radius and to electric polarization of the embedding material. At threshold, the permeation is {\emph{intermittent}} on a nanosecond timescale, and strongly enhanced by the presence of an ion inside the channel, providing a possible mechanism for gating. Confined water remains surprisingly fluid and bulk-like. Its behaviour differs strikingly from that of a reference Lennard-Jones fluid, which tends to contract into a highly layered structure inside the channel.Comment: 4 pages, 4 figure

A Hydrophobic Gate in an Ion Channel: The Closed State of the Nicotinic Acetylcholine Receptor

The nicotinic acetylcholine receptor (nAChR) is the prototypic member of the `Cys-loop' superfamily of ligand-gated ion channels which mediate synaptic neurotransmission, and whose other members include receptors for glycine, gamma-aminobutyric acid, and serotonin. Cryo-electron microscopy has yielded a three dimensional structure of the nAChR in its closed state. However, the exact nature and location of the channel gate remains uncertain. Although the transmembrane pore is constricted close to its center, it is not completely occluded. Rather, the pore has a central hydrophobic zone of radius about 3 A. Model calculations suggest that such a constriction may form a hydrophobic gate, preventing movement of ions through a channel. We present a detailed and quantitative simulation study of the hydrophobic gating model of the nicotinic receptor, in order to fully evaluate this hypothesis. We demonstrate that the hydrophobic constriction of the nAChR pore indeed forms a closed gate. Potential of mean force (PMF) calculations reveal that the constriction presents a barrier of height ca. 10 kT to the permeation of sodium ions, placing an upper bound on the closed channel conductance of 0.3 pS. Thus, a 3 A radius hydrophobic pore can form a functional barrier to the permeation of a 1 A radius Na+ ion. Using a united atom force field for the protein instead of an all atom one retains the qualitative features but results in differing conductances, showing that the PMF is sensitive to the detailed molecular interactions.Comment: Accepted by Physical Biology; includes a supplement and a supplementary mpeg movie can be found at http://sbcb.bioch.ox.ac.uk/oliver/download/Movies/watergate.mp

Topological dissection of the membrane transport protein Mhp1 derived from cysteine accessibility and mass spectrometry

Cys accessibility and quantitative intact mass spectrometry (MS) analyses have been devised to study the topological transitions of Mhp1, the membrane protein for sodium-linked transport of hydantoins from Microbacterium liquefaciens. Mhp1 has been crystallised in three forms (outward-facing open, outward-facing occluded with substrate bound, and inward-facing open). We show that one natural cysteine residue, Cys327, out of three, has an enhanced solvent accessibility in the inward-facing (relative to the outward-facing) form. Reaction of the purified protein, in detergent, with the thiol-reactive N-ethylmalemide (NEM), results in modification of Cys327, suggesting that Mhp1 adopts predominantly inward-facing conformations. Addition of either sodium ions or the substrate 5-benzyl-L-hydantoin (L-BH) does not shift this conformational equilibrium, but, systematic co-addition of the two results in an attenuation of labelling, indicating a shift toward outward-facing conformations that can be interpreted using conventional enzyme kinetic analyses. Such measurements can afford the Km for each ligand as well as the stoichiometry of ion-substrate coupled conformational changes. Mutations that perturb the substrate binding site either result in the protein being unable to adopt outward-facing conformations or in a global destabilisation of structure. The methodology combines covalent labeling, mass spectrometry and kinetic analyses in a straightforward workflow applicable to a range of systems, enabling the interrogation of changes in a protein’s conformation required for function at varied concentrations of substrates, and the consequences of mutations on these conformational transitions

Inherent Dynamics of the Acid-Sensing Ion Channel 1 Correlates with the Gating Mechanism

A combination of computational and experimental approaches reveals the dynamics of ASIC1 gating, involving a deformation of the channel that triggers “twist-to-open” motions of the channel pore

InBO3 and ScBO3 at high pressures: an ab initio study of elastic and thermodynamic properties

We have theoretically investigated the elastic properties of calcite-type orthoborates ABO(3) (A= Sc and In) at high pressure by means of ab initio total-energy calculations. From the elastic stiffness coefficients, we have obtained the elastic moduli (B, G and E), Poisson's ratio (nu), B/G ratio, universal elastic anisotropy index (A(U)), Vickers hardness, and sound wave velocities for both orthoborates. Our simulations show that both borates are more resistive to volume compression than to shear deformation (B > G). Both compounds are ductile and become more ductile, with an increasing elastic anisotropy, as pressure increases. We have also calculated some thermodynamic properties, like Debye temperature and minimum thermal conductivity. Finally, we have evaluated the theoretical mechanical stability of both borates at high hydrostatic pressures. It has been found that the calcite-type structure of InBO3 and ScBO3 becomes mechanically unstable at pressures beyond 56.2 and 57.7 GPa, respectively. (C) 2016 Elsevier Ltd. All rights reserved.This study is supported by the Spanish MICINN projects MAT2013-46649-C4-2-P/3-P and MAT2015-71070-REDC. H.M.O., A.M., and P.R-H. acknowledge computing time provided by Red Espanola de Supercomputacion (RES) and MALTA-Cluster. J.A.S. acknowledges Juan de la Cierva fellowship program for financial support.Gomis, O.; Ortiz, HM.; Sans Tresserras, JÁ.; Manjón Herrera, FJ.; Santamaría-Pérez, D.; Rodríguez-Hernández, P.; Muñoz, A. (2016). InBO3 and ScBO3 at high pressures: an ab initio study of elastic and thermodynamic properties. Journal of Physics and Chemistry of Solids. 98:198-208. https://doi.org/10.1016/j.jpcs.2016.07.002S1982089

Structural and elastic properties of defect chalcopyrite HgGa2S4 under high pressure

In this work, we focus on the study of the structural and elastic properties of mercury digallium sulfide (HgGa2S4) at high pressures. This compound belongs to the family of AB(2)X(4) ordered-vacancy compounds and exhibits a tetragonal defect chalcopyrite structure. X-ray diffraction measurements at room temperature have been performed under compression up to 15.1 GPa in a diamond anvil cell. Our measurements have been complemented and compared with ab initio total energy calculations. The axial compressibility and the equation of state of the low-pressure phase of HgGa2S4 have been experimentally and theoretically determined and compared to other related ordered-vacancy compounds. The pressure dependence of the theoretical cation-anion and vacancy-anion distances and compressibilities in HgGa2S4 are reported and discussed in comparison to other related ordered-vacancy compounds. Finally, the pressure dependence of the theoretical elastic constants and elastic moduli of HgGa2S4 has been studied. Our calculations indicate that the low-pressure phase of HgGa2S4 becomes mechanically unstable above 13.8 GPa. (C) 2013 Elsevier B. V. All rights reserved.This study was supported by the Spanish government MEC under Grants No: MAT2010-21270-C04-01/03/04 and CTQ2009-14596-C02-01, by the Comunidad de Madrid and European Social Fund (S2009/PPQ-1551 4161893), by MALTA Consolider Ingenio 2010 Project (CSD2007-00045), by Generalitat Valenciana (GVA-ACOMP-2013-1012), and by the Vicerrectorado de Investigacion y Desarrollo of the Universidad Politecnica de Valencia (UPV2011-0914 PAID-05-11 and UPV2011-0966 PAID-06-11). E.P-G., A. M., and P.R-H. acknowledge computing time provided by Red Espa ola de Supercomputacion (RES) and MALTA-Cluster. J.A.S. acknowledges Juan de la Cierva fellowship program for financial support.Gomis Hilario, O.; Santamaría-Pérez, D.; Vilaplana Cerda, RI.; Luna Molina, R.; Sans, JA.; Manjón Herrera, FJ.; Errandonea, D.... (2014). Structural and elastic properties of defect chalcopyrite HgGa2S4 under high pressure. Journal of Alloys and Compounds. 583:70-78. https://doi.org/10.1016/j.jallcom.2013.08.123S707858

Minimum Free Energy Path of Ligand-Induced Transition in Adenylate Kinase

Large-scale conformational changes in proteins involve barrier-crossing transitions on the complex free energy surfaces of high-dimensional space. Such rare events cannot be efficiently captured by conventional molecular dynamics simulations. Here we show that, by combining the on-the-fly string method and the multi-state Bennett acceptance ratio (MBAR) method, the free energy profile of a conformational transition pathway in Escherichia coli adenylate kinase can be characterized in a high-dimensional space. The minimum free energy paths of the conformational transitions in adenylate kinase were explored by the on-the-fly string method in 20-dimensional space spanned by the 20 largest-amplitude principal modes, and the free energy and various kinds of average physical quantities along the pathways were successfully evaluated by the MBAR method. The influence of ligand binding on the pathways was characterized in terms of rigid-body motions of the lid-shaped ATP-binding domain (LID) and the AMP-binding (AMPbd) domains. It was found that the LID domain was able to partially close without the ligand, while the closure of the AMPbd domain required the ligand binding. The transition state ensemble of the ligand bound form was identified as those structures characterized by highly specific binding of the ligand to the AMPbd domain, and was validated by unrestrained MD simulations. It was also found that complete closure of the LID domain required the dehydration of solvents around the P-loop. These findings suggest that the interplay of the two different types of domain motion is an essential feature in the conformational transition of the enzyme