Free-energy landscape of glycerol permeation through aquaglyceroporin GlpF determined from steered molecular dynamics simulations

Molecular dynamics study of water transport through AQP5-R188C mutant causing palmoplantar keratoderma (PPK) using the gating mechanism concept

It is widely known that any disruption to the water regulation in aquaporins (AQPs) leads to numerous important diseases. However, studies of dynamics and energetics of disease-causing mutations in the aquaporins on the molecular level are still limited. In the present work, the effects of a skin disease-causing mutant, R188C, on the structure of AQP5 and water transport mechanism within this mutated aquaporin are investigated using the concept of gating mechanism. Our results have revealed that the R188C mutation causes a remarkable increase in the pore radius inside the selectivity filter (SF) region facilitating the passage of water molecules. This observation is supported by plotting the free energy profiles of water molecules transport and calculating permeability values through AQP5-R188C, such that the energy barrier in the SF region of the pores was substantially reduced by this mutation, and therefore, the translocation of water molecules was improved. The total averaged osmotic permeability for R188C has been computed as about 11-fold of the wild-type permeability. However, a comparison between the osmotic permeability values related to the open conformation of CE revealed that this coefficient for AQP5-R188C is about 6.5 times larger than that of wt-AQP5, which can be a more accurate value according to the gating mechanism associated with the constriction region of the aquaporin.

Computational Studies of Molecular Permeation through Connexin26 Channels

A signal property of connexin channels is the ability to mediate selective diffusive movement of molecules through plasma membrane(s), but the energetics and determinants of molecular movement through these channels have yet to be understood. Different connexin channels have distinct molecular selectivities that cannot be explained simply on the basis of size or charge of the permeants. To gain insight into the forces and interactions that underlie selective molecular permeation, we investigated the energetics of two uncharged derivatized sugars, one permeable and one impermeable, through a validated connexin26 (Cx26) channel structural model, using molecular dynamics and associated analytic tools. The system is a Cx26 channel equilibrated in explicit membrane/solvent, shown by Brownian dynamics to reproduce key conductance characteristics of the native channel. The results are consistent with the known difference in permeability to each molecule. The energetic barriers extend through most of the pore length, rather than being highly localized as in ion-specific channels. There is little evidence for binding within the pore. Force decomposition reveals how, for each tested molecule, interactions with water and the Cx26 protein vary over the length of the pore and reveals a significant contribution from hydrogen bonding and interaction with K(+). The flexibility of the pore width varies along its length, and the tested molecules have differential effects on pore width as they pass through. Potential sites of interaction within the pore are defined for each molecule. The results suggest that for the tested molecules, differences in hydrogen bonding and entropic factors arising from permeant flexibility substantially contribute to the energetics of permeation. This work highlights factors involved in selective molecular permeation that differ from those that define selectivity among atomic ions.

Polarizable empirical force field for acyclic polyalcohols based on the classical Drude oscillator

A polarizable empirical force field for acyclic polyalcohols based on the classical Drude oscillator is presented. The model is optimized with an emphasis on the transferability of the developed parameters among molecules of different sizes in this series and on the condensed-phase properties validated against experimental data. The importance of the explicit treatment of electronic polarizability in empirical force fields is demonstrated in the cases of this series of molecules with vicinal hydroxyl groups that can form cooperative intra- and intermolecular hydrogen bonds. Compared to the CHARMM additive force field, improved treatment of the electrostatic interactions avoids overestimation of the gas-phase dipole moments resulting in significant improvement in the treatment of the conformational energies and leads to the correct balance of intra- and intermolecular hydrogen bonding of glycerol as evidenced by calculated heat of vaporization being in excellent agreement with experiment. Computed condensed phase data, including crystal lattice parameters and volumes and densities of aqueous solutions are in better agreement with experimental data as compared to the corresponding additive model. Such improvements are anticipated to significantly improve the treatment of polymers in general, including biological macromolecules.

Glycerol modulates water permeation through Escherichia coli aquaglyceroporin GlpF

Among aquaglyceroporins that transport both water and glycerol across the cell membrane, Escherichia coli glycerol uptake facilitator (GlpF) is the most thoroughly studied. However, one question remains: Does glycerol modulate water permeation? This study answers this fundamental question by determining the three-dimensional potential of mean force of glycerol along the permeation path through GlpF's conducting pore. There is a deep well near the Asn-Pro-Ala (NPA) motifs (6.5kcal/mol below the bulk level) and a barrier near the selectivity filter (10.1kcal/mol above the well bottom). This profile owes its existence to GlpF's perfect steric arrangement: The glycerol-protein van der Waals interactions are attractive near the NPA but repulsive elsewhere in the conducting pore. In light of the single-file nature of waters and glycerols lining up in GlpF's amphipathic pore, it leads to the following conclusion: Glycerol modulates water permeation in the μM range. At mM concentrations, GlpF is glycerol-saturated and a glycerol residing in the well occludes the conducting pore. Therefore, water permeation is fully correlated to glycerol dissociation that has an Arrhenius activation barrier of 6.5kcal/mol. Validation of this theory is based on the existent in vitro data, some of which have not been given the proper attention they deserved: The Arrhenius activation barriers were found to be 7kcal/mol for water permeation and 9.6kcal/mol for glycerol permeation; The presence of up to 100mM glycerol did not affect the kinetics of water transport with very low permeability, in apparent contradiction with the existent theories that predicted high permeability (0M glycerol).

In silico study of Aquaporin V: Effects and affinity of the central pore-occluding lipid

Because of its roles in human physiology, Aquaporin V (AQP5), a major intrinsic protein, has been a subject of many in vitro studies. In particular, a 2008 experiment produced its crystal structure at 2.0Å resolution, which is in a tetrameric conformation consisting of four protomers. Each protomer forms an amphipathic pore that is fit for water permeation. The tetramer has a pore along its quasi-symmetry axis formed by quadruplets of hydrophobic residues (every protomer contributes equally to the quadruplets). A lipid, phosphatidylserine (PS6), is bound to AQP5 in the central pore, totally occluding it. A 2009 experiment showed that AQP5 facilitates not only permeation of water but also permeation of hydrophobic gas molecules across the cell membrane. In this article, we present an in silico study of AQP5 to elucidate the effects of PS6's binding to and dissociating from AQP5's central pore. Computing the lipid's chemical-potential along its dissociation path, we find that PS6 inhibits the function of the central pore with an IC(50) in the micromolar range. Examining the central pore and the interstices between two adjacent protomers, we propose that nonpolar gas molecules (O(2)) permeate through AQP5's hydrophobic central pore when un-occluded.

Molecular dynamics simulation and bioinformatics study on yeast aquaporin Aqy1 from Pichia pastoris

In the present study, an equilibrated system for the Aqy1 tetramer was developed, and molecular biophysics modeling showed that the Aqy1 channel was blocked by Tyr-31 in the N-terminus, which was also supported by the free energy profiles. However, bioinformatics analysis of the amino acid sequence of Aqy1 indicated this Tyr-31 is not conserved across all fungi. Analysis of the equilibrated structure showed that the central pore along the four-fold axis of the tetramers is formed with hydrophobic amino acid residues. In particular, Phe-90, Trp-198, and Phe-202 form the narrowest part of the pore. Therefore, water molecules are not expected to translocate through the central pore, a hypothesis that we confirmed by molecular dynamics simulations. Each monomer of the Aqy1 tetramers forms a channel whose walls consist mostly of hydrophilic residues, transporting through the selectivity filter containing Arg-227, His-212, Phe-92, and Ala-221, and the two conserved Asn-Pro-Ala (NPA) motifs containing asparagines 224 and 112. In summary, not all fungal aquaporins share the same gating mechanism by a tyrosine residue in the N-terminus, and the structural analysis in the present study should aid our understanding of aquaporin structure and its functional implications.

Water transport in human aquaporin-4: molecular dynamics (MD) simulations

Aquaporin-4 (AQP4) is the predominant water channel in the central nervous system, where it has been reported to be involved in many pathophysiological roles including water transport. In this paper, the AQP4 tetramer was modeled from its PDB structure file, embedded in a palmitoyl-oleoyl-phosphatidyl-choline (POPC) lipid bilayer, solvated in water, then minimized and equilibrated by means of molecular dynamics simulations. Analysis of the equilibrated structure showed that the central pore along the fourfold axis of the tetramers is formed with hydrophobic amino acid residues. In particular, Phe-195, Leu-191 and Leu-75, form the narrowest part of the pore. Therefore water molecules are not expected to transport through the central pore, which was confirmed by MD simulations. Each monomer of the AQP4 tetramers forms a channel whose walls consist mostly of hydrophilic residues. There are eight water molecules in single file observed in each of the four channels, transporting through the selectivity filter containing Arg-216, His-201, Phe-77, Ala-210, and the two conserved Asn-Pro-Ala (NPA) motifs containing Asn-213 and Asn-97. By using Brownian dynamics fluctuation-dissipation-theorem (BD-FDT), the overall free-energy profile was obtained for water transporting through AQP4 for the first time, which gives a complete map of the entire channel of water permeation.

Mercury inhibits the L170C mutant of aquaporin Z by making waters clog the water channel

We conduct in silico experiments of the L170C mutant of the Escherichia coli aquaporin Z (AQPZ) with and without mercury bonded to residue Cys 170. We find that bonding mercury to Cys 170 does not induce consequential structural changes to the protein. We further find that mercury does not stick in the middle of the water channel to simply occlude water permeation, but resides on the wall of the water pore. However, we observe that the water permeation coefficient of L170C-Hg(+) (with one mercury ion bonded to Cys 170) is approximately half of that of the mercury-free L170C. We examine the interactions between the mercury ion and the waters in its vicinity and find that five to six waters are strongly attracted by the mercury ion, occluding the space of the water channel. Therefore we conclude that mercury, at low concentration, inhibits AQPZ-L170C mutant by making water molecules clog the water channel.

Exploring the free-energy landscapes of biological systems with steered molecular dynamics

We perform steered molecular dynamics (SMD) simulations and use the Brownian dynamics fluctuation-dissipation-theorem (BD-FDT) to accurately compute the free-energy profiles for several biophysical processes of fundamental importance: hydration of methane and cations, binding of benzene to T4-lysozyme L99A mutant, and permeation of water through aquaglyceroporin. For each system, the center-of-mass of the small molecule (methane, ion, benzene, and water, respectively) is steered (pulled) at a given speed over a period of time, during which the system transitions from one macroscopic state/conformation (State A) to another one (State B). The mechanical work of pulling the system is measured during the process, sampling a forward pulling path. Then the reverse pulling is conducted to sample a reverse path from B back to A. Sampling a small number of forward and reverse paths, we are able to accurately compute the free-energy profiles for all the afore-listed systems that represent various important aspects of biological physics. The numerical results are in excellent agreement with the experimental data and/or other computational studies available in the literature.