Ultrafast 2D IR anisotropy of water reveals reorientation during hydrogen-bond switching

Rearrangements of the hydrogen bond network of liquid water are believed to involve rapid and concerted hydrogen bond switching events, during which a hydrogen bond donor molecule undergoes large angle molecular reorientation as it exchanges hydrogen bonding partners. To test this picture of hydrogen bond dynamics, we have performed ultrafast 2D IR spectral anisotropy measurements on the OH stretching vibration of HOD in D(2)O to directly track the reorientation of water molecules as they change hydrogen bonding environments. Interpretation of the experimental data is assisted by modeling drawn from molecular dynamics simulations, and we quantify the degree of molecular rotation on changing local hydrogen bonding environment using restricted rotation models. From the inertial 2D anisotropy decay, we find that water molecules initiating from a strained configuration and relaxing to a stable configuration are characterized by a distribution of angles, with an average reorientation half-angle of 10°, implying an average reorientation for a full switch of ≥20°. These results provide evidence that water hydrogen bond network connectivity switches through concerted motions involving large angle molecular reorientation.     

Structure and hydrogen bonding in neat N-methylacetamide: classical molecular dynamics and Raman spectroscopy studies of a liquid of peptidic fragments

The results of classical molecular dynamics (MD) simulations and Raman spectroscopy studies of neat liquid N-methylacetamide (NMA), the simplest model system relevant to the peptides, are reported as a function of temperature and pressure. The MD simulations predict that near ambient conditions, the molecules form a hydrogen bond network consisting primarily of linear chains. Both the links between molecules within the hydrogen-bonded chains and the associations between chains are stabilized by weak methyl-donated "improper" hydrogen bonds. The three-dimensional structural motifs observed in the liquid show some similarity to protein beta-sheets. The temperature and pressure dependence of the hydrogen bond network, as probed by the mode frequency of the experimentally determined amide-I Raman band, blue shifts on heating and red shifts under compression, respectively, suggesting weakened and enhanced hydrogen bonding in response to temperature and pressure increases. Disruption of the hydrogen-bonding network is clearly observed in the simulation data as temperature is increased, whereas the improper hydrogen bonding is enhanced under compression to reduce the energetic cost of increasing the packing fraction. Because of the neglect of polarizability in the molecular model, the computed dielectric constant is underestimated compared to the experimental value, indicating that the simulation may underestimate dipolar coupling in the liquid.     

Fast centroid molecular dynamics: a force-matching approach for the predetermination of the effective centroid forces

A fast centroid molecular dynamics (CMD) methodology is proposed in which the effective centroid forces are predetermined through a force-matching algorithm applied to a standard path integral molecular dynamics simulation. The resulting method greatly reduces the computational cost of generating centroid trajectories, thus extending the applicability of CMD. The method is applied to the study of liquid para-hydrogen at two state points and liquid ortho-deuterium at one state point. The static and dynamical results are compared to those obtained from full adiabatic CMD simulations and found to be in excellent agreement for all three systems; the transport properties are also compared to experiment and found to have a similar level of agreement.     

Insights into Tris‐(2‐Hydroxylethyl)methylammonium Methylsulfate Aqueous Solutions

We report herein a combined experimental–computational study on tris‐(2‐hydroxylethyl)methylammonium methylsulfate in water solutions, as a representative ionic liquid of the aqueous‐solution behavior of hydroxylammonium‐based ionic liquids. Relevant thermophysical properties were measured as a function of mixture composition and temperature. Classical molecular dynamics simulations were performed to infer microscopic structural features. The reported results for ionic liquid in water‐rich solutions show that it behaves as isolated non‐interacting ions solvated by water molecules, through well‐defined solvation shells, exerting a disrupting effect on the water hydrogen bonding network. Nevertheless, as ionic liquid concentration increase, interionic association increases, even for diluted water solutions, evolving from the typical behavior of strong electrolytes in solution toward large interacting structures. For ionic‐liquid‐rich mixtures, water exerts a minor disrupting effect on the fluid’s structuring because it occupies regions around each ion (developing water–ion hydrogen bonds) but without significantly weakening anion–cation interactions.     

Accelerating the convergence of path integral dynamics with a generalized Langevin equation

The quantum nature of nuclei plays an important role in the accurate modelling of light atoms such as hydrogen, but it is often neglected in simulations due to the high computational overhead involved. It has recently been shown that zero-point energy effects can be included comparatively cheaply in simulations of harmonic and quasiharmonic systems by augmenting classical molecular dynamics with a generalized Langevin equation (GLE). Here we describe how a similar approach can be used to accelerate the convergence of path integral (PI) molecular dynamics to the exact quantum mechanical result in more strongly anharmonic systems exhibiting both zero point energy and tunnelling effects. The resulting PI-GLE method is illustrated with applications to a double-well tunnelling problem and to liquid water.     

Multiscale QM/MM molecular dynamics study on the first steps of guanine damage by free hydroxyl radicals in solution

Understanding the damage of DNA bases from hydrogen abstraction by free OH radicals is of particular importance to understanding the indirect effect of ionizing radiation. Previous studies address the problem with truncated DNA bases as ab initio quantum simulations required to study such electronic-spin-dependent processes are computationally expensive. Here, for the first time, we employ a multiscale and hybrid quantum mechanical-molecular mechanical simulation to study the interaction of OH radicals with a guanine-deoxyribose-phosphate DNA molecular unit in the presence of water, where all of the water molecules and the deoxyribose-phosphate fragment are treated with the simplistic classical molecular mechanical scheme. Our result illustrates that the presence of water strongly alters the hydrogen-abstraction reaction as the hydrogen bonding of OH radicals with water restricts the relative orientation of the OH radicals with respect to the DNA base (here, guanine). This results in an angular anisotropy in the chemical pathway and a lower efficiency in the hydrogen-abstraction mechanisms than previously anticipated for identical systems in vacuum. The method can easily be extended to single- and double-stranded DNA without any appreciable computational cost as these molecular units can be treated in the classical subsystem, as has been demonstrated here.     

On the origin of the stabilization of the zwitterionic resting state of cysteine proteases: a theoretical study

Papain-like cysteine proteases are ubiquitous proteolytic enzymes. The protonated His199/deprotonated Cys29 ion pair (cathepsin B numbering) in the active site is essential for their proper functioning. The presence of this ion pair stands in contrast to the corresponding intrinsic residue p K a values, indicating a strong influence of the enzyme environment. In the present work we show by molecular dynamics simulations on quantum mechanical/molecular mechanical (QM/MM) potentials that the ion pair is stabilized by a complex hydrogen bond network which comprises several amino acids situated in the active site of the enzyme and 2-4 water molecules. QM/MM reaction path computations for the proton transfer from His199 to the thiolate of the Cys29 moiety indicate that the ion pair is about 32-36 kJ mol (-1) more stable than the neutral form if the whole hydrogen bonding network is active. Without any hydrogen bonding network the ion pair is predicted to be significantly less stable than the neutral form. QM/MM charge deletion analysis and QM model calculations are used to quantify the stabilizing effect of the active-site residues and the L1 helix in favor of the zwitterionic form. The active-site water molecules contribute about 30 kJ mol (-1) to the overall stabilization. Disruption of the hydrogen bonding network upon substrate binding is expected to enhance the nucleophilic reactivity of the thiolate.     

Hydrogen-Bond Structure and Low-Frequency Dynamics of Electrolyte Solutions: Hydration Numbers from ab Initio Water Reorientation Dynamics and Dielectric Relaxation Spectroscopy

We present an atomistic simulation scheme for the determination of the hydration number (h) of aqueous electrolyte solutions based on the calculation of the water dipole reorientation dynamics. In this methodology, the time evolution of an aqueous electrolyte solution generated from ab initio molecular dynamics simulations is used to compute the reorientation time of different water subpopulations. The value of h is determined by considering whether the reorientation time of the water subpopulations is retarded with respect to bulk-like behavior. The application of this computational protocol to magnesium chloride (MgCl₂) solutions at different concentrations (0.6–2.8 mol kg⁻¹) gives h values in excellent agreement with experimental hydration numbers obtained using GHz-to-THz dielectric relaxation spectroscopy. This methodology is attractive because it is based on a well-defined criterion for the definition of hydration number and provides a link with the molecular-level processes responsible for affecting bulk solution behavior. Analysis of the ab initio molecular dynamics trajectories using radial distribution functions, hydrogen bonding statistics, vibrational density of states, water-water hydrogen bonding lifetimes, and water dipole reorientation reveals that MgCl₂ has a considerable influence on the hydrogen bond network compared with bulk water. These effects have been assigned to the specific strong Mg-water interaction rather than the Cl-water interaction.     

Circular hydrogen bond networks on the surface of beta-ribofuranose in aqueous solution

This paper examines the hydration structure on the surface of beta-ribofuranose in aqueous solution, using the ab initio molecular dynamics method. In particular, we focus on circular hydrogen bond networks involving two ribofuranose oxygens and three water molecules. In our simulations, the circular hydrogen bond networks near the ring oxygen of beta-ribofuranose are found to be significantly influenced by the orientation of the hydroxymethyl group. The arrangements of hydrogen bonds observed in the circular hydrogen bond networks are both homodromic and antidromic. To explain these observations, we analyze the electronic properties of the first-hydration-shell water molecules and the OH groups of beta-ribofuranose, using the centers of their maximally localized Wannier functions. The dipole moments of the proton-accepting first-hydration-shell water molecules in our well-defined circular hydrogen bond networks are found to increase by about 0.3 D compared with that of liquid water, indicating the relatively strong polarization effects created by the interactions between the OH groups of the solute and the surrounding water molecules. Our analysis also implies that circular H-bond networks cannot be fully explained from a simple geometrical point of view.     

Water network dynamics at the critical moment of a peptide's beta-turn formation: a molecular dynamics study

All-atom molecular dynamics simulations for a single molecule of Leu-Enkephalin in aqueous solution have been used to study the role of the water network during the formation of beta-turns. We give a detailed account of the intramolecular hydrogen bonding, the water-peptide hydrogen bonding, and the orientation and residence times of water molecules focusing on the short critical periods of transition to the stable beta-turns. These studies suggest that, when intramolecular hydrogen bonding between the first and fourth residue of the beta-turn is not present, the disruption of the water network and the establishment of water bridges constitute decisive factors in the formation and stability of the beta-turn. Finally, we provide possible explanations and mechanisms for the formations of different kinds of beta-turns.     

Influence of temperature on the structure of liquid water

Molecular dynamics simulations have been carried out for liquid water at 7 different temperatures to understand the nature of hydrogen bonding at molecular level through the investigation of the effects of temperature on the geometry of water molecules. The changes in bond length and bond angle of water molecules from gaseous state to liquid state have been observed, and the change in the bond angle of water molecules in liquid against temperature has been revealed, which has not been seen in literature so far. The analysis of the radial distribution functions and the coordinate numbers shows that, on an average, each water molecule in liquid acts as both receptor and donor, and forms at least two hydrogen bonds with its neigbors. The analysis of the results also indicates that the water molecules form clusters in liquid.     

Water order profiles on phospholipid/cholesterol membrane bilayer surfaces

Water is pivotal in the stabilization of macromolecular biological structures, although the dynamic ensemble structure of water near to molecular surfaces has yet to be fully understood. We show, through molecular simulation and fluorescence measurements, that water at the membrane surface is substantially more ordered than bulk water, due to a loss of hydrogen bonding between water molecules, coupled with an alignment of lipid and water dipole moments. Ordering of the water leads to a gradient in the effective dielectric permittivity, which is evident in both the molecular simulations and the fluorescence measurements. A lower effective dielectric permittivity was correlated with a decreasing degree of hydrogen bonding over the same spatial range. The water molecules closest to the lipid headgroup oxygen atoms form hydrogen bonds which exhibit a mean lifetime of 6.3 ps, compared with a mean lifetime of water-water hydrogen bonds of less than 2 ps. Membranes made up purely of phosphatidylcholine (PC) were compared with those made with a PC/cholesterol ratio relevant to cell membranes. Clear differences were found between these membrane configurations. These observations point to molecular structural differences in the surface environments of membranes and may underlie regional differences in the surface biophysical properties of membrane microdomains.     

Hydrogen bonding in liquid alcohols: a computer simulation study

A series of molecular dynamics simulations has been performed to investigate hydrogen bonding in liquid alcohols. The systems considered have been methanol, ethanol, ethylene glycol and glycerol at 298 K. The hydrogen bonding statistics as well as the mean lifetime of the hydrogen bonds are analyzed. The results are compared with those corresponding to liquid water.     

Investigation of the crystal structure of two amphiphilic diols

The crystal and molecular structures of (R)-4'-(2,3-dihydroxypropyloxy)-4-cyanobiphenyl, 1 and 4'-(cis-cis-3,5-dihydroxycyclohexyloxy)-4-cyanobiphenyl, 5 have been determined. The packing of compound 1 occurs in sheets, with a dense hydrogen bonding network within the layer of one sheet and to the next sheet, whereby the cyano groups are incorporated into the hydrogen bonding network. The hydrogen bonding scheme of compound 5 involves a water molecule; there is some disorder with regard to the hydrogen atoms bound to the water molecules.     

Intermolecular polarizability dynamics of aqueous formamide liquid mixtures studied by molecular dynamics simulations

A molecular dynamics simulation study is presented for the relaxation of the polarizability anisotropy in liquid mixtures of formamide and water, using a dipolar induction scheme that involves the intrinsic polarizability and first hyperpolarizability tensors of the molecules, and the dipole-quadrupole polarizability of water species. The long time diffusive decay of the collective polarizability anisotropy correlations exhibits a substantial slowing down as the formamide mole fraction increases in the mixture. The diffusive times for the polarizability relaxation obtained from the authors' simulations are in good agreement with optical Kerr effect experimental data, and they are found to correlate nearly linearly with the estimated mean lifetimes of the hydrogen bonds within the mixture, suggesting that the relaxation of the hydrogen bond network is responsible to some extent for the collective relaxation of the polarizability anisotropy of the mixture. The short time behavior of the polarizability anisotropy relaxation was investigated by computing the nuclear response function, R(t), which is very rapidly dominated by the formamide contribution as it is added to water, due to the much larger polarizability anisotropy of formamide molecules compared to that of water. Several contributions to the Raman spectrum were also analyzed as a function of composition, and the dynamical origin of the different bands was determined.     

Liquid structures of water, methanol, and hydrogen fluoride at ambient conditions from first principles molecular dynamics simulations with a dispersion corrected density functional

Using first principles molecular dynamics simulations in the isobaric-isothermal ensemble (T = 300 K, p = 1 atm) with the Becke-Lee-Yang-Parr exchange/correlation functional and a dispersion correction due to Grimme, the hydrogen bonding networks of pure liquid water, methanol, and hydrogen fluoride are probed. Although an accurate density is found for water with this level of electronic structure theory, the average liquid densities for both hydrogen fluoride and methanol are overpredicted by 50 and 25%, respectively. The radial distribution functions indicate somewhat overstructured liquid phases for all three compounds. The number of hydrogen bonds per molecule in water is about twice as high as for methanol and hydrogen fluoride, though the ratio of cohesive energy over number of hydrogen bonds is lower for water. An analysis of the hydrogen-bonded aggregates revealed the presence of mostly linear chains in both hydrogen fluoride and methanol, with a few stable rings and chains spanning the simulation box in the case of hydrogen fluoride. Only an extremely small fraction of smaller clusters was found for water, indicating that its hydrogen bond network is significantly more extensive. A special form of water with on average about two hydrogen bonds per molecule yields a hydrogen-bonding environment significantly different from the other two compounds.     

Direct observation of the substitution effects on the hydrogen bridge dynamics in selected Schiff bases--a comparative molecular dynamics study

We have studied substituent effects on the properties of the intramolecular hydrogen bond of some ortho-hydroxy Schiff bases using density functional theory (DFT) based first-principle molecular dynamics (FPMD) and path integral molecular dynamics. The studied compounds possess a strong intramolecular hydrogen bond (r((O???N)) ≤ 2.6 A?), which can be tuned by substitution to either (i) enhance the basicity of the acceptor moiety by induction effects or (ii) decrease the hydrogen bond length through steric repulsion. DFT calculations and FPMD were employed to investigate structural and dynamical properties of the selected molecules, while quantum effects on the structural properties were assessed using path integral FPMD. The simulations were performed in vacuo and in the solid state to study the influence of the environment on the hydrogen bond and spectroscopic properties. We give computational support to the suggestion that induction effects are less effective to tune the intramolecular hydrogen bond properties of the discussed ortho-hydroxy Schiff bases than the steric or the environmental effects.     

All‐atom Molecular Dynamic Simulations Combined with the Chemical Shifts Study on the Weak Interactions of Ethanol‐water System

All‐atom molecular dynamics (MD) simulation combined with chemical shifts was performed to investigate the interactions over the entire concentration range of the ethanol (EtOH)‐water system. The results of the simulation were adopted to explain the NMR experiments by hydrogen bonding analysis. The strong hydrogen bonds and weak C–H···O contacts coexist in the mixtures through the analysis of the radial distribution functions. And the liquid structures in the whole concentration of EtOH‐water mixtures can be classified into three regions by the statistic analysis of the hydrogen‐bonding network in the MD simulations. Moreover, the chemical shifts of the hydrogen atom are in agreement with the statistical results of the average number hydrogen bonds in the MD simulations. Interestingly, the excess relative extent of η_rel^E calculated by the MD simulations and chemical shifts in the EtOH aqueous solutions shows the largest deviation at x_EtOH≈0.18. The excess properties present good agreement with the excess enthalpy in the concentration dependence.     

Rigid quantum Monte Carlo simulations of condensed molecular matter: water clusters in the n=2-->8 range

The numerical advantage of quantum Monte Carlo simulations of rigid bodies relative to the flexible simulations is investigated for some simple systems. The results show that if high frequency modes in molecular condensed matter are predominantly in the ground state, the convergence of path integral simulations becomes nonuniform. Rigid body quantum parallel tempering simulations are necessary to accurately capture thermodynamic phenomena in the temperature range where the dynamics are influenced by intermolecular degrees of freedom; the stereographic projection path integral adapted for quantum simulations of asymmetric tops is a significantly more efficient strategy compared with Cartesian coordinate simulations for molecular condensed matter under these conditions. The reweighted random series approach for stereographic path integral Monte Carlo is refined and implemented for the quantum simulation of water clusters treated as an assembly of rigid asymmetric tops.