Correlation between the dynamics of hydrogen bonds and the local density reorganization in the protein hydration layer

An atomistic molecular dynamics simulation of the protein villin headpiece subdomain or HP-36 has been carried out with explicit water to explore the microscopic inhomogeneity of local density reorganization of the hydration layers of the three alpha-helical segments of the protein. The density reorganization of the hydration layer of helix-3 is found to occur faster than that for the hydration layers of the other two helices. It is noticed that such inhomogeneous density reorganization at the surface of different secondary structures exhibits excellent correlation with the microscopic dynamics of hydrogen bonds between the protein residues and the hydration water. Further, it is observed that the reorientation of water molecules involved in the formation and breaking of protein-water or water-water hydrogen bonds plays an important role in determining the dynamics of local density of the hydration layer. The faster density reorganization of the hydration layer of helix-3 is also consistent with the functionality of HP-36, as helix-3 contains several active site residues.     

Secondary structure sensitivity of hydrogen bond lifetime dynamics in the protein hydration layer

The heterogeneous nature of a protein surface plays an essential role in its biological activity and molecular recognition, and this role is mediated at least partly through the surrounding water molecules. We have performed atomistic molecular dynamics simulations of an aqueous solution of HP-36 to investigate the correlation between the dynamics of the hydration layer water molecules and the lifetimes of protein-water hydrogen bonds. The nonexponential hydrogen bond lifetime correlation functions have been analyzed by using the formalism of Luzar and Chandler, which allowed identification of the quasi-bound states in the surface and quantification of the dynamic equilibrium between quasi-bound and free water molecules in terms of time-dependent rate of interconversion. It is noticed that, irrespective of the structural heterogeneity of different segments of the protein, namely the three alpha-helices, the positively charged amino acid residues form longer-lived hydrogen bonds with water. The overall relaxation behavior of protein-water hydrogen bonds is found to differ significantly among the three helices of the protein. Study of water number density fluctuation reveals that the hydration layer of helix-3 is much less rigid, which can be correlated with faster structural relaxation of the hydrogen bonds between its residues and water. This also agrees excellently with faster translational and rotational motions of water near helix-3, and hence the lower rigidity of its hydration layer. The lower rigidity of the helix-3 hydration layer also correlates well with the biological activity of the protein, as several of the active-site residues of HP-36 are located in helix-3.     

Dynamics of water in the hydration layer of a partially unfolded structure of the protein HP-36

Atomistic molecular dynamics simulations of the folded native structure and a partially unfolded molten globule structure of the protein villin headpiece subdomain or HP-36 have been carried out with explicit solvent to explore the effects of unfolding on the dynamical behavior of water present in the hydration layers of different segments (three alpha-helices) of the protein. The calculations revealed that the unfolding of helix-2 influences the translational and rotational motions of water present in the hydration layers of the three helices in a heterogeneous manner. It is observed that a correlation exists between the unfolding of helix-2 and the microscopic kinetics of protein-water hydrogen bonds formed by its residues. This in turn has an influence on the rigidity of the hydration layers of the helices in the unfolded structure versus that in the folded native structure. These results should provide a microscopic explanation to recent solvation dynamics experiments on folded native and unfolded structures of proteins.     

Aqueous Heterogeneity at the Air/Water Interface Revealed by 2D‐HD‐SFG Spectroscopy

Water molecules interact strongly with each other through hydrogen bonds. This efficient intermolecular coupling causes strong delocalization of molecular vibrations in bulk water. We study intermolecular coupling at the air/water interface and find intermolecular coupling 1) to be significantly reduced and 2) to vary strongly for different water molecules at the interface—whereas in bulk water the coupling is homogeneous. For strongly hydrogen‐bonded OH groups, coupling is roughly half of that of bulk water, due to the lower density in the near‐surface region. For weakly hydrogen‐bonded OH groups that absorb around 3500 cm^?1, which are assigned to the outermost, yet hydrogen‐bonded OH groups pointing towards the liquid, coupling is further reduced by an additional factor of 2. Remarkably, despite the reduced structural constraints imposed by the interfacial hydrogen‐bond environment, the structural relaxation is slow and the intermolecular coupling of these water molecules is weak.     

Local heterogeneous dynamics of water around lysozyme: a computer simulation study

Water present near the surface of a protein exhibits dynamic properties different from that of water in the pure bulk state. In this work, we have carried out atomistic molecular dynamics simulation of an aqueous solution of hen egg-white lysozyme. Attempts have been made to explore the correlation between the local heterogeneous mobility of water around the protein segments and the rigidity of the hydration layers with the microscopic dynamics of hydrogen bonds formed by water molecules with the protein residues. The kinetics of breaking and reformation of hydrogen bonds involving the surface water molecules have been calculated. It is found that the reformations of broken hydrogen bonds are more frequent for the hydration layers of those segments of the protein that are more rigid. The calculation of the low-frequency vibrational modes of hydration layer water molecules reveals that the protein influences the transverse and longitudinal degrees of freedom of water around it in a differential manner. These findings can be verified by appropriate experimental studies.     

Structure and dynamics of methanol in water: a quantum mechanical charge field molecular dynamics study

An ab initio quantum mechanical charge field molecular dynamics simulation was carried out for one methanol molecule in water to analyze the structure and dynamics of hydrophobic and hydrophilic groups. It is found that water molecules around the methyl group form a cage-like structure whereas the hydroxyl group acts as both hydrogen bond donor and acceptor, thus forming several hydrogen bonds with water molecules. The dynamic analyses correlate well with the structural data, evaluated by means of radial distribution functions, angular distribution functions, and coordination number distributions. The overall ligand mean residence time, τ identifies the methanol molecule as structure maker. The relative dynamics data of hydrogen bonds between hydroxyl of methanol and water molecules prove the existence of both strong and weak hydrogen bonds. The results obtained from the simulation are in excellent agreement with the experimental results for dilute solution of CH(3)OH in water. The overall hydration shell of methanol consists in average of 18 water molecules out of which three are hydrogen bonded.     

Hydrogen bond lifetime dynamics at the interface of a surfactant monolayer

The dynamics of water near the polar headgroups of surfactants in a monolayer adsorbed at the air/water interface is likely to play a decisive role in determining the physical behavior of such organized assemblies. We have carried out an atomistic molecular dynamics (MD) simulation of a monolayer of the anionic surfactant sodium bis(2-ethyl-1-hexyl) sulfosuccinate (aerosol-OT or AOT) adsorbed at the air/water interface. The simulation is performed at room temperature with a surface coverage of that at the critical micelle concentration (78 Angstrom(2)/molecule). Detailed analyses of the lifetime dynamics of surfactant-water (SW) and water-water (WW) hydrogen bonds at the interface have been carried out. The nonexponential hydrogen bond lifetime correlation functions have been analyzed by using the formalism of Luzar and Chandler, which allowed identification of the bound states at the interface and quantification of the dynamic equilibrium between bound and quasi-free water molecules, in terms of time-dependent relaxation rates. It is observed that the water molecules present in the first hydration layer form strong hydrogen bonds with the surfactant headgroups and hence have longer lifetimes. Importantly, it is found that the overall relaxation of the SW hydrogen bonds is faster for those water molecules which form two hydrogen bonds with the surfactant headgroups than those forming one such hydrogen bond. Equally interestingly, it is further noticed that water molecules beyond the first hydration layer form weaker hydrogen bonds than pure bulk water.     

Reconsideration on hydrogen bond strengthening or cleavage of photoexcited coumarin 102 in aqueous solvent: a DFT/TDDFT study

In this work, the excited-state hydrogen bonding dynamics of photoexcited coumarin 102 in aqueous solvent is reconsidered. The electronically excited states of the hydrogen bonded complexes formed by coumarin 102 (C102) chromophore and the hydrogen donating water solvent have been investigated using the time-dependent density functional theory method. Two intermolecular hydrogen bonds between C102 and water molecules are considered. The previous works (Wells et al., J Phys Chem A 2008, 112, 2511) have proposed that one intermolecular hydrogen bond would be strengthened and the other one would be cleaved upon photoexcitation to the electronically excited states. However, our theoretical calculations have demonstrated that both the two intermolecular hydrogen bonds between C102 solute and H(2)O solvent molecules are significantly strengthened in electronically excited states by comparison with those in ground state. Hence, we have confirmed again that intermolecular hydrogen bonds between C102 chromophore and aqueous solvents are strengthened not cleaved upon electronic excitation, which is in accordance with Zhao's works.     

Hydration-dependent protein dynamics revealed by molecular dynamics simulation of crystalline staphylococcal nuclease

Molecular dynamics simulations of crystalline Staphylococcal nuclease in full and minimal hydration states were performed to study hydration effects on protein dynamics at temperatures ranging from 100 to 300 K. In a full hydration state (hydration ratio in weight, h=0.49), gaps are fully filled with water molecules, whereas only crystal waters are included in a minimal hydration state (h=0.09). The inflection of the atomic mean-square fluctuation of protein as a function of temperature, known as the glass-like transition, is observed at approximately 220 K in both cases, which is more significant in the full hydration state. By examining the temperature dependence of residual fluctuation, we found that the increase of fluctuations in the loop and terminal regions, which are exposed to water, is much greater than that in other regions in the full hydration state, but the mobilities of the corresponding regions are relatively restricted in the minimal hydration state by intermolecular contact. The atomic mean-square fluctuation of water molecules in the full hydration state at 300 K is 1 order of magnitude greater than that in the minimal hydration state. Above the transition temperature, most water molecules in the full hydration state behave like bulk water and act as a lubricant for protein dynamics. In contrast, water molecules in the minimal hydration state tend to form more hydrogen bonds with the protein, restricting the fluctuation of these water molecules to the level of the protein. Thus, intermolecular interaction and solvent mobility are important to understand the glass-like transition in proteins.     

Sensitivity of hydrogen bond lifetime dynamics to the presence of ethanol at the interface of a phospholipid bilayer

Atomistic molecular dynamics simulations of a fully hydrated liquid crystalline lamellar phase of a dimyrystoylphosphatidylcholine lipid bilayer containing ethanol at 1:1 composition as well as of the pure lamellar phase of the bilayer have been performed. Detailed analyses have been carried out to investigate the effects of ethanol, if any, on the lifetime dynamics of lipid-water and water-water hydrogen bonds in the hydration layer of the lipid headgroups. The nonexponential hydrogen bond lifetime correlation functions have been analyzed by using the formalism of Luzar and Chandler, which allowed the identification of the bound states at the bilayer interface and the quantification of the dynamic equilibrium between the bound and the free water molecules, in terms of time-dependent relaxation rates. The calculations show that the overall relaxation of phosphate-water hydrogen bonds is faster in the presence of ethanol. Studies of the residence time and the number fluctuation of the hydration layer water molecules reveal that the presence of ethanol molecules decreases the rigidity of the lipid hydration layer.     

Dynamics of protein and peptide hydration

Biological processes often involve the surfaces of proteins, where the structural and dynamic properties of the aqueous solvent are modified. Information about the dynamics of protein hydration can be obtained by measuring the magnetic relaxation dispersion (MRD) of the water (2)H and (17)O nuclei or by recording the nuclear Overhauser effect (NOE) between water and protein protons. Here, we use the MRD method to study the hydration of the cyclic peptide oxytocin and the globular protein BPTI in deeply supercooled solutions. The results provide a detailed characterization of water dynamics in the hydration layer at the surface of these biomolecules. More than 95% of the water molecules in contact with the biomolecular surface are found to be no more than two-fold motionally retarded as compared to bulk water. In contrast to small nonpolar molecules, the retardation factor for BPTI showed little or no temperature dependence, suggesting that the exposed nonpolar residues do not induce clathrate-like hydrophobic hydration structures. New NOE data for oxytocin and published NOE data for BPTI were analyzed, and a mutually consistent interpretation of MRD and NOE results was achieved with the aid of a new theory of intermolecular dipolar relaxation that accounts explicitly for the dynamic perturbation at the biomolecular surface. The analysis indicates that water-protein NOEs are dominated by long-range dipolar couplings to bulk water, unless the monitored protein proton is near a partly or fully buried hydration site where the water molecule has a long residence time.     

Monolayer of aerosol-OT surfactants adsorbed at the air/water interface: an atomistic computer simulation study

An atomistic molecular dynamics (MD) simulation has been carried out to investigate the structural and dynamical properties of a monolayer of the anionic surfactant sodium bis(2-ethyl-1-hexyl) sulfosuccinate (aerosol-OT or AOT) adsorbed at the air/water interface. The simulation is performed at room temperature and at a surface coverage corresponding to that at its critical micelle concentration (78 A(2)/molecule). The estimated thickness of the adsorbed layer is in good agreement with neutron reflection data. The study shows that the surfactants exhibit diffusive motion in the plane of the interface. It is observed that the surfactant monolayer has a strong influence in restricting both the translational and reorientational motions of the water molecules close to the interface. A drastic difference in the dipolar reorientational motion of water molecules in the aqueous layer is observed with a small variation of the distance from the surfactant headgroups. It has been observed that the water molecules in the first hydration layer (region 1) form strong hydrogen bonds with surfactant headgoups. This results in the slower structural relaxation of water-water hydrogen bonds in the first hydration layer compared to that in the pure bulk water. Most interestingly, we notice that the water molecules present in the layer immediately after the first hydration layer form weaker hydrogen bonds and thus relax faster than even pure bulk water.     

Thickness of the hydration layer of a protein from molecular dynamics simulation

Water molecules around a protein exhibit slow dynamics with respect to that of pure bulk water. One important issue in protein hydration is the thickness of the hydration layer (i.e., the distance from the protein surface up to which the water dynamics is influenced by the protein). Estimation of thickness is crucial to understand better the properties of "biological water" and the role that it plays in guiding the protein's function. We have performed an atomistic molecular dynamics simulation of an aqueous solution of the protein villin headpiece subdomain or HP-36 to estimate the thickness of its hydration water. In particular, several dynamical properties of water around different segments (three alpha-helices) of the protein have been calculated by varying the thickness of the hydration layers. It is found that in general the influence of the helices on water properties extends beyond the first hydration layer. However, the heterogeneous nature of water among the first hydration layers of the three helices diminishes as the thickness is increased. It indicates that, for a small protein such as HP-36, the thickness of "biological water" is uniform for different segments of the protein.     

Intermolecular hydrogen bonds in hetero-complexes of biologically active aromatic ligands: Monte Carlo simulations results

We report results of the Monte Carlo simulations of systems containing heterodimers of biological active ligands and water molecules. The study was designed to identify the possible formation of intermolecular hydrogen bonds in such systems in order to investigate the molecular mechanisms of hetero-association of aromatic ligands in aqueous solution. The geometry optimization and the calculation of the atomic charges of free ligands were carried out at DFT/B3LYP level of theory. Monte Carlo simulations with Metropolis algorithm were used to determine the low energy conformations of heterodimers in water clusters. The analysis of the Monte Carlo simulation results allows us to describe in detail the hydration properties of all investigated heterodimers and to determine the intermolecular hydrogen bonds between the functional donor–acceptor groups for some of hetero-associates under investigation. In the case of heterodimers without intermolecular hydrogen bonds, the additional stabilization of these hetero-complexes can be explained by the formation the water bridges between donor and acceptor groups of the ligands.     

Dynamic properties of water around a protein-DNA complex from molecular dynamics simulations

Formation of protein-DNA complex is an important step in regulation of genes in living organisms. One important issue in this problem is the role played by water in mediating the protein-DNA interactions. In this work, we have carried out atomistic molecular dynamics simulations to explore the heterogeneous dynamics of water molecules present in different regions around a complex formed between the DNA binding domain of human TRF1 protein and a telomeric DNA. It is demonstrated that such heterogeneous water motions around the complex are correlated with the relaxation time scales of hydrogen bonds formed by those water molecules with the protein and DNA. The calculations reveal the existence of a fraction of extraordinarily restricted water molecules forming a highly rigid thin layer in between the binding motifs of the protein and DNA. It is further proved that higher rigidity of water layers around the complex originates from more frequent reformations of broken water-water hydrogen bonds. Importantly, it is found that the formation of the complex affects the transverse and longitudinal degrees of freedom of surrounding water molecules in a nonuniform manner.     

Hydration state of nonionic surfactant monolayers at the liquid/vapor interface: structure determination by vibrational sum frequency spectroscopy

The OH stretching region of water molecules in the vicinity of nonionic surfactant monolayers has been investigated using vibrational sum frequency spectroscopy (VSFS) under the polarization combinations ssp, ppp, and sps. The surface sensitivity of the VSFS technique has allowed targeting the few water molecules present at the surface with a net orientation and, in particular, the hydration shell around alcohol, sugar, and poly(ethylene oxide) headgroups. Dramatic differences in the hydration shell of the uncharged headgroups were observed, both in comparison to each another and in comparison to the pure water surface. The water molecules around the rigid glucoside and maltoside sugar rings were found to form strong hydrogen bonds, similar to those observed in tetrahedrally coordinated water in ice. In the case of the poly(ethylene oxide) surfactant monolayer a significant ordering of both strongly and weakly hydrogen bonded water was observed. Moreover, a band common to all the surfactants studied, clearly detected at relatively high frequencies in the polarization combinations ppp and sps, was assigned to water species located in proximity to the surfactant hydrocarbon tail phase, with both hydrogen atoms free from hydrogen bonds. An orientational analysis provided additional information on the water species responsible for this band.     

Bis(dl‐cysteinium) oxalate

In the title compound, 2C₃H₈NO₂S⁺·C₂O₄²⁻, the oxalate anion occupies an inversion centre and is coordinated to cysteine molecules of different chirality (l and d ) via O—H...O and N—H...O hydrogen bonds, the resulting cysteine–oxalate stoichiometry in the crystal structure being 2:1. The oxalate anion is completely deprotonated, whereas cysteine has a positively charged –NH₃⁺ group and a neutral protonated carboxyl group. The structure is built from infinite hydrogen‐bonded triple layers, consisting of an oxalate layer in the middle with layers of l ‐ and d ‐cysteine molecules on either side. The thiol groups are at the external sides of the layers and form S—H...O hydrogen bonds with the carboxyl groups of neighbouring cysteine molecules. An interesting feature of the structure is the occurrence of short S...S contacts between SH groups of molecules in neighbouring layers, which form not S—H...S but S—H...O intermolecular hydrogen bonds. Due to the effects of crystal packing and intermolecular hydrogen‐bond formation, the conformation of the cysteine cation in the title structure is different from that calculated theoretically for an individual cation, as well as from those of cysteine zwitterions in crystals of pure cysteine.     

Hydration of poly(ethylene oxide)s in aqueous solution as studied by dielectric relaxation measurements

The hydration state of poly(ethylene oxide)s (PEOs) in aqueous solutions was investigated using dielectric relaxation measurements at 25 degrees C over a frequency range up to 20 GHz, which is the relaxation frequency of water molecules in a bulk state. The dielectric relaxation spectra obtained indicated decomposition into two major and one minor relaxation modes with relaxation times of 8.3, 22, and 250 ps, respectively. The two major modes were attributed to rotational relaxation of water molecules belonging to the bulk state and water molecules hydrogen bonded to ethylene oxide (EO) monomer units. The number of hydration water molecules per EO unit depended on the molar mass of PEO (M) and reached a constant value of 3.7 at M > 1500, which agrees with the value obtained by other experiments.     

Time‐dependent density functional theory study on the electronic excited‐state geometric structure,infrared spectra,and hydrogen bonding of a doubly hydrogen‐bonded complex

The geometric structures and infrared (IR) spectra in the electronically excited state of a novel doubly hydrogen‐bonded complex formed by fluorenone and alcohols, which has been observed by IR spectra in experimental study, are investigated by the time‐dependent density functional theory (TDDFT) method. The geometric structures and IR spectra in both ground state and the S₁ state of this doubly hydrogen‐bonded FN‐2MeOH complex are calculated using the DFT and TDDFT methods, respectively. Two intermolecular hydrogen bonds are formed between FN and methanol molecules in the doubly hydrogen‐bonded FN‐2MeOH complex. Moreover, the formation of the second intermolecular hydrogen bond can make the first intermolecular hydrogen bond become slightly weak. Furthermore, it is confirmed that the spectral shoulder at around 1700 cm^?1 observed in the IR spectra should be assigned as the doubly hydrogen‐bonded FN‐2MeOH complex from our calculated results. The electronic excited‐state hydrogen bonding dynamics is also studied by monitoring some vibraitonal modes related to the formation of hydrogen bonds in different electronic states. As a result, both the two intermolecular hydrogen bonds are significantly strengthened in the S₁ state of the doubly hydrogen‐bonded FN‐2MeOH complex. The hydrogen bond strengthening in the electronically excited state is similar to the previous study on the singly hydrogen‐bonded FN‐MeOH complex and play important role on the photophysics of fluorenone in solutions. © 2009 Wiley Periodicals, Inc. J Comput Chem 2009     

Conformational changes of a calix[8]arene derivative at the air-water interface

The particular behavior of a p-tert-butyl calix[8]arene derivative (C8A) has been studied at the air-water interface using surface pressure-area isotherms, surface potential-area isotherms, film relaxation measurements, Brewster angle microscopy (BAM), and infrared spectroscopy for Langmuir-Blodgett films. Thus, it is observed that the properties of the film, for example, isotherms, domain formation, and FTIR spectra, recorded during the first compression cycle differ appreciably from those during the second compression and following cycles. The results obtained are interpreted on the basis of the conformational changes of the C8A molecules by surface pressure, allowing us to inquire into the inter- and intramolecular interactions (hydrogen bonds) of those molecules. Thus, the compression induces changes in the kind of hydrogen bonds from intra- and intermolecular with other C8A molecules to hydrogen bonds with water molecules.