Entropy and dynamics of water in hydration layers of a bilayer.
07 Nov 2010-Journal of Chemical Physics (American Institute of Physics)-Vol. 133, Iss: 17, pp 174704-174704
TL;DR: The translational diffusion of water in the vicinity of the head groups is found to be in a subdiffusive regime and the rotational diffusion constant increases going away from the interface, supported by the slower reorientational relaxation of the dipole vector and OH bond vector of interfacial water.
Abstract: We compute the entropy and transport properties of water in the hydration layer of dipalmitoylphosphatidylcholine bilayer by using a recently developed theoretical scheme [two-phase thermodynamic model, termed as 2PT method; S.-T. Lin et al., J. Chem. Phys. 119, 11792 (2003)] based on the translational and rotational velocity autocorrelation functions and their power spectra. The weights of translational and rotational power spectra shift from higher to lower frequency as one goes from the bilayer interface to the bulk. Water molecules near the bilayer head groups have substantially lower entropy (48.36 J/mol/K) than water molecules in the intermediate region (51.36 J/mol/K), which have again lower entropy than the molecules (60.52 J/mol/K) in bulk. Thus, the entropic contribution to the free energy change (TΔS) of transferring an interface water molecule to the bulk is 3.65 kJ/mol and of transferring intermediate water to the bulk is 2.75 kJ/mol at 300 K, which is to be compared with 6.03 kJ/mol for melting of ice at 273 K. The translational diffusion of water in the vicinity of the head groups is found to be in a subdiffusive regime and the rotational diffusion constant increases going away from the interface. This behavior is supported by the slower reorientational relaxation of the dipole vector and OH bond vector of interfacial water. The ratio of reorientational relaxation time for Legendre polynomials of order 1 and 2 is approximately 2 for interface, intermediate, and bulk water, indicating the presence of jump dynamics in these water molecules.
Citations
More filters
TL;DR: The calculations reveal that the inhomogeneous water ordering and density distributions around the helices are correlated with their relative hydrophobicity, and identifies the existence of a narrow relatively dehydrated region containing randomly organized "quasi-free" water molecules beyond the first layer of "bound" waters at the protein surface.
Abstract: Atomistic molecular dynamics simulation of an aqueous solution of the small protein HP-36 has been carried out with explicit solvent at room temperature. Efforts have been made to explore the influence of the protein on the relative packing and ordering of water molecules around its secondary structures, namely, three α-helices. The calculations reveal that the inhomogeneous water ordering and density distributions around the helices are correlated with their relative hydrophobicity. Importantly, we have identified the existence of a narrow relatively dehydrated region containing randomly organized “quasi-free” water molecules beyond the first layer of “bound” waters at the protein surface. These water molecules with relatively weaker binding energies form the transition state separating the “bound” and “free” water molecules at the interface. Further, increased contribution of solid-like caging motions of water molecules around the protein is found to be responsible for reduced fluidity of the hydration layer. Interestingly, we notice that the hydration layer of helix-3 is more fluidic with relatively higher entropy as compared to the hydration layers of the other two helical segments. Such characteristics of helix-3 hydration layer correlate well with the activity of HP-36, as helix-3 contains the active site of the protein.
21 citations
TL;DR: Investigation of the effects of DMSO on the H-bond dynamics at the lipid-water interface using a combination of ultrafast two-dimensional infrared spectroscopy and molecular dynamics simulations takes an important step towards characterizing the environments around lipid membranes, which are essential to numerous biological processes.
Abstract: Dimethyl sulfoxide (DMSO) is a common cosolvent and cryopreservation agent used to freeze cells and tissues. DMSO alters the H-bond structure of water, but its interactions with biomolecules, and s...
21 citations
TL;DR: In this paper, the authors studied the thermodynamics of filling of a (6,6) carbon nanotube in a temperature range from 273 to 353K and with different strengths of the nanotubes-water interaction.
Abstract: Spontaneous entry of water molecules inside single-wall carbon nanotubes (SWCNTs) has been confirmed by both simulations and experiments. Using molecular dynamics simulations, we have studied the thermodynamics of filling of a (6,6) carbon nanotube in a temperature range from 273 to 353K and with different strengths of the nanotube-water interaction. From explicit energy and entropy calculations using the two-phase thermodynamics method, we have presented a thermodynamic understanding of the filling behaviour of a nanotube. We show that both the energy and the entropy of transfer decrease with increasing temperature. On the other hand, scaling down the attractive part of the carbon-oxygen interaction results in increased energy of transfer while the entropy of transfer increases slowly with decreasing the interaction strength. Our results indicate that both energy and entropy favour water entry into (6,6) SWCNTs. Our results are compared with those of several recent studies of water entry into carbon nanotubes.
20 citations
TL;DR: In this paper, the authors present thermodynamic integration simulations for the binding of mannose and methanoic acid onto the calcite surface producing free energy of binding values of −2.89 and −1.64 kJ mol−1, respectively.
Abstract: We present thermodynamic integration simulations for the binding of mannose and methanoic acid onto the {10.4} calcite surface producing free energy of binding values of −2.89 and −1.64 kJ mol–1, respectively. We extract the entropy of binding from vacuum-based simulations and use these values to determine the entropy of binding for surface water molecules which is ∼6 J mol–1 K–1.
19 citations
Cites methods from "Entropy and dynamics of water in hy..."
...This has been successfully implemented by Lin et al(16,17) to calculate the entropy of water in complicated systems such as a bilayer.(18) Other methods in the literature include the construction of the non-ideal component of the entropy via configurational sampling(19) where an approximation is made that of “hypothetical scanning”....
[...]
TL;DR: The entropy and potential energy of interfacial water show distinct changes at the bilayers melting transition, indicating a strong correlation in the thermodynamic state of water and the accompanying first-order phase transition of the bilayer membrane.
Abstract: We investigate the effect of bilayer melting transition on thermodynamics and dynamics of interfacial water using molecular dynamics simulation with the two-phase thermodynamic model. We show that the diffusivity of interface water depicts a dynamic crossover at the chain melting transition following an Arrhenius behavior until the transition temperature. The corresponding change in the diffusion coefficient from the bulk to the interface water is comparable with experimental observations found recently for water near 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) vesicles Phys. Chem. Chem. Phys. 13, 7732 (2011)]. The entropy and potential energy of interfacial water show distinct changes at the bilayer melting transition, indicating a strong correlation in the thermodynamic state of water and the accompanying first-order phase transition of the bilayer membrane. DOI: 10.1103/PhysRevLett.110.018303
18 citations
References
More filters
TL;DR: In this paper, a method is described to realize coupling to an external bath with constant temperature or pressure with adjustable time constants for the coupling, which can be easily extendable to other variables and to gradients, and can be applied also to polyatomic molecules involving internal constraints.
Abstract: In molecular dynamics (MD) simulations the need often arises to maintain such parameters as temperature or pressure rather than energy and volume, or to impose gradients for studying transport properties in nonequilibrium MD A method is described to realize coupling to an external bath with constant temperature or pressure with adjustable time constants for the coupling The method is easily extendable to other variables and to gradients, and can be applied also to polyatomic molecules involving internal constraints The influence of coupling time constants on dynamical variables is evaluated A leap‐frog algorithm is presented for the general case involving constraints with coupling to both a constant temperature and a constant pressure bath
25,256 citations
TL;DR: A parallel message-passing implementation of a molecular dynamics program that is useful for bio(macro)molecules in aqueous environment is described and can handle rectangular periodic boundary conditions with temperature and pressure scaling.
8,195 citations
TL;DR: The design includes an extraction of virial and periodic boundary conditions from the loops over pairwise interactions, and special software routines to enable rapid calculation of x–1/2.
Abstract: GROMACS 3.0 is the latest release of a versatile and very well optimized package for molecular simulation. Much effort has been devoted to achieving extremely high performance on both workstations and parallel computers. The design includes an extraction of virial and periodic boundary conditions from the loops over pairwise interactions, and special software routines to enable rapid calculation of x–1/2. Inner loops are generated automatically in C or Fortran at compile time, with optimizations adapted to each architecture. Assembly loops using SSE and 3DNow! Multimedia instructions are provided for x86 processors, resulting in exceptional performance on inexpensive PC workstations. The interface is simple and easy to use (no scripting language), based on standard command line arguments with self-explanatory functionality and integrated documentation. All binary files are independent of hardware endian and can be read by versions of GROMACS compiled using different floating-point precision. A large collection of flexible tools for trajectory analysis is included, with output in the form of finished Xmgr/Grace graphs. A basic trajectory viewer is included, and several external visualization tools can read the GROMACS trajectory format. Starting with version 3.0, GROMACS is available under the GNU General Public License from http://www.gromacs.org.
6,375 citations
01 Jan 1981
TL;DR: In this article, a three-point charge model (on hydrogen and oxygen positions) with a Lennard-Jones 6-12 potential on the oxygen positions only was developed, and parameters for the model were determined from 12 molecular dynamics runs covering the two-dimensional parameter space of charge and oxygen repulsion.
Abstract: For molecular dynamics simulations of hydrated proteins a simple yet reliable model for the intermolecular potential for water is required. Such a model must be an effective pair potential valid for liquid densities that takes average many-body interactions into account. We have developed a three-point charge model (on hydrogen and oxygen positions) with a Lennard-Jones 6–12 potential on the oxygen positions only. Parameters for the model were determined from 12 molecular dynamics runs covering the two-dimensional parameter space of charge and oxygen repulsion. Both potential energy and pressure were required to coincide with experimental values. The model has very satisfactory properties, is easily incorporated into protein-water potentials, and requires only 0.25 sec computertime per dynamics step (for 216 molecules) on a CRAY-1 computer.
5,336 citations
TL;DR: The first and second papers in this series, which make it possible to interpret entropy data in terms of a physical picture, are applied to binary solutions, and equations are derived relating energy and volume changes when a solution is formed to the entropy change for the process as discussed by the authors.
Abstract: The ideas of the first and second papers in this series, which make it possible to interpret entropy data in terms of a physical picture, are applied to binary solutions, and equations are derived relating energy and volume changes when a solution is formed to the entropy change for the process. These equations are tested against data obtained by various authors on mixtures of normal liquids, and on solutions of non‐polar gases in normal solvents. Good general agreement is found, and it is concluded that in such solutions the physical picture of molecules moving in a ``normal'' manner in each others' force fields is adequate. As would be expected, permanent gases, when dissolved in normal liquids, loosen the forces on neighboring solvent molecules producing a solvent reaction which increases the partial molal entropy of the solute.Entropies of vaporization from aqueous solutions diverge strikingly from the normal behavior established for non‐aqueous solutions. The nature of the deviations found for non‐polar solutes in water, together with the large effect of temperature upon them, leads to the idea that the water forms frozen patches or microscopic icebergs around such solute molecules, the extent of the iceberg increasing with the size of the solute molecule. Such icebergs are apparently formed also about the non‐polar parts of the molecules of polar substances such as alcohols and amines dissolved in water, in agreement with Butler's observation that the increasing insolubility of large non‐polar molecules is an entropy effect. The entropies of hydration of ions are discussed from the same point of view, and the conclusion is reached that ions, to an extent which depends on their sizes and charges, may cause a breaking down of water structure as well as a freezing or saturation of the water nearest them. Various phenomena recorded in the literature are interpreted in these terms. The influence of temperature on certain salting‐out coefficients is interpreted in terms of entropy changes. It appears that the salting‐out phenomenon is at least partly a structural effect. It is suggested that structural influences modify the distribution of ions in an electrolytesolution, and reasons are given for postulating the existence of a super‐lattice structure in solutions of LaCl3 and of EuCl3. An example is given of a possible additional influence of structural factors upon reacting tendencies in aqueous solutions.
2,572 citations