Hydration of krypton and consideration of clathrate models of hydrophobic effects from the perspective of quasi-chemical theory

Ab initio molecular dynamics (AIMD) results on a krypton-water liquid solution are presented and compared to recent XAFS results for the radial hydration structure for a Kr atom in liquid water solution. Though these AIMD calculations have important limitations of scale, the comparisons with the liquid solution results are satisfactory and significantly different from the radial distributions extracted from the data on the solid Kr/H(2)O clathrate hydrate phase. The calculations also produce the coordination number distribution that can be examined for metastable coordination structures suggesting possibilities for clathrate-like organization; none are seen in these results. Clathrate pictures of hydrophobic hydration are discussed, as is the quasi-chemical theory that should provide a basis for clathrate pictures. Outer shell contributions are discussed and estimated; they are positive and larger than the positive experimental hydration free energy of Kr(aq), implying that inner shell contributions must be negative and of comparable size. Clathrate-like inner shell hydration structures on a Kr atom solute are obtained for some, but not all, of the coordination number cases observed in the simulation. The structures found have a delicate stability. Inner shell coordination structures extracted from the simulation of the liquid, and then subjected to quantum chemical optimization, always decomposed. Interactions with the outer shell material are decisive in stabilizing coordination structures observed in liquid solution and in clathrate phases. The primitive quasi-chemical estimate that uses a dielectric model for the influence of the outer shell material on the inner shell equilibria gives a contribution to hydration free energy that is positive and larger than the experimental hydration free energy. The 'what are we to tell students' question about hydrophobic hydration, often answered with structural clathrate pictures, is then considered; we propose an alternative answer that is consistent with successful molecular theories of hydrophobic effects and based upon distinctive observable properties of liquid water. Considerations of parsimony, for instance Ockham's razor, then suggest that additional structural hypotheses in response to 'what are we to tell students' are not required at this stage.     

Modification of enzyme activity in reversed micelles through clathrate hydrate formation.

The physical phenomenon of clathrate hydrate formation in protein-containing reversed micelles is described. Hydrate formation in reversed micelles is a method of adjusting the water to surfactant molar ratio, wo, which influences micellar size. Lipase and alpha-chymotrypsin encapsulated in large reversed micelles of high wo show significant enhancements in activity when the micelle size is reduced through hydrate formation. Alternate methods of micelle size adjustments also show enhancements in activity. The implications for improving the activity of such encapsulated enzymes recovered from fermentation media through phase transfer into reversed micelles are discussed.     

Electrodissociation of clathrate-like structures

In the present work, we develop molecular dynamics (MD) simulations in the NPT (isobaric–isothermic) ensemble to analyse the effect of an external electrostatic field over a cubic methane hydrate crystallite. The amplitude of the field is in the range 0.5–3.0 V/nm. For the simulations, we used the SPC/E rigid water model and a single-site model for methane at a temperature of 248 K and a pressure of 20 bar. When the external electrostatic field is applied, the water dipoles are oriented in such a way that the methane molecules can diffuse far away from the water cages, hence the clathrate dissociation takes place. This last phenomenon was observed for intensities above 1.5 V/nm. Taking the final configuration of each run as input, we develop a new set of MD simulations, and we observe that the stable clathrate is not recovered immediately when the external electrostatic field is turned off due to limitations in the simulation time.     

The thermal properties of binary structure sI clathrate hydrate from molecular dynamics simulation

In this work, we present temperature dependence of lattice parameter and normalised lattice parameter in the atmospheric pressure and 120 bar and also pressure dependence of unit cell volume and normalised unit cell volume at 150 and 250?K for variety guests with different size, polarity and guest–host hydrogen bonding capability such as trimethylene oxide (TMO), ethylene oxide (EO), formaldehyde (FA), cyclobutane (CB), cyclopropane (CP) and ethane (Et) in the large cages with CH₄ in small cages of sI clathrate hydrates by molecular dynamics simulations. The obtained values of lattice parameters for the guest species are compatible with the experimental values. These clathrate hydrates are simulated with TIP4P/ice four-site water potential. Herein, isobaric thermal expansivity and isothermal compressibility are calculated at a temperature range of 50–250?K and a wide pressure range. These structural properties have been compared for guests which they are isoelectronic and have similar masses but with different size and polarity. We use molecular dynamics simulations to relate microscopic guest properties, like guest–host hydrogen bonding to macroscopic sI clathrate hydrate properties. The temperature dependence of thermodynamic properties such as constant-volume and constant-pressure heat capacity is presented in the atmospheric pressure for these guest species.     

Protein recovery from reversed micellar solutions through contact with a pressurized gas phase.

We describe a new process for the recovery of encapsulated protein from reversed micellar solution in concentrated form. The method involves desolubilization of the protein by decreasing solvent density through gas dissolution. Under appropriate thermodynamic conditions, the micellar water pool can be converted to clathrate hydrates. Protein recovery is facilitated by clathrate hydrate formation, which causes the desolubilized protein to exist in a solid phase, distinct from the micellar supernatant. The process is carried out without any ionic strength or pH modification.     

Dissociation mechanism of gas hydrates (I,II, H) of alkane molecules: a comparative molecular dynamics simulation

Employing NPT molecular dynamics method with consistent valence force field, the dissociation processes of sI, sII and sH gas hydrates are simulated at different temperatures and at a constant pressure of 100 MPa. The dissociation mechanisms of gas hydrates are revealed by analysing the structural snapshots, radial distribution functions and diffusion coefficients at different temperatures. As temperature increases, the diffusion rates of water molecules and guest molecules increase; thus the clathrate skeleton formed by water molecules with hydrogen bonds distorts and breaks down; meanwhile the guest molecules encapsulated in the water cavities are released. The size of guest molecules affects the dissociation behaviour of gas hydrate. In addition, the dissociation behaviour also relies on the structural phase of gas hydrates.     

Collective properties of hydration: long range and specificity of hydrophobic interactions.

We report results of molecular dynamics (MD) simulations of composite model solutes in explicit molecular water solvent, eliciting novel aspects of the recently demonstrated, strong many-body character of hydration. Our solutes consist of identical apolar (hydrophobic) elements in fixed configurations. Results show that the many-body character of PMF is sufficiently strong to cause 1) a remarkable extension of the range of hydrophobic interactions between pairs of solute elements, up to distances large enough to rule out pairwise interactions of any type, and 2) a SIF that drives one of the hydrophobic solute elements toward the solvent rather than away from it. These findings complement recent data concerning SIFs on a protein at single-residue resolution and on model systems. They illustrate new important consequences of the collective character of hydration and of PMF and reveal new aspects of hydrophobic interactions and, in general, of SIFs. Their relevance to protein recognition, conformation, function, and folding and to the observed slight yet significant nonadditivity of functional effects of distant point mutations in proteins is discussed. These results point out the functional role of the configurational and dynamical states (and related statistical weights) corresponding to the complex configurational energy landscape of the two interacting systems: biomolecule + water.     

On the Stability of Clathrate Hydrates Encaging Polar Guest Molecules: Contrast in the Hydrogen Bonds of Methylamine and Methanol Hydrates

Abstract

The stability of clathrate hydrates encaging highly polar guests has been investigated in order to explain the experimental observation that some amines form clathrate hydrates but alcohols act as inhibitor to hydrate formation. We choose methylamine and methanol as guest species and examine the stable structure, at which the total potential energy has a minimum value. At the local minima of those two hydrates, the potential energies of water-water and guest-water, and their hydrogen bonded networks are compared. It is found that methanol does not retain the host lattice structure, while the host-network structure is kept in the presence of methylamine. It is shown that the difference in the magnitude of the partial charge on the hydrogen atom between the hydroxyl and amino groups plays a much more significant role on the stability of both clathrate hydrates than the difference in molecular geometry. This is supported from the result of a methylamine-like model that has the same partial charges on the atoms in the hydrophilic site as methanol.     

Dynamics of tRNA at Different Levels of Hydration

The influence of hydration on the nanosecond timescale dynamics of tRNA is investigated using neutron scattering spectroscopy. Unlike protein dynamics, the dynamics of tRNA is not affected by methyl group rotation. This allows for a simpler analysis of the influence of hydration on the conformational motions in RNA. We find that hydration affects the dynamics of tRNA significantly more than that of lysozyme. Both the characteristic length scale and the timescale of the conformational motions in tRNA depend strongly on hydration. Even the characteristic temperature of the so-called “dynamical transition” appears to be hydration-dependent in tRNA. The amplitude of the conformational motions in fully hydrated tRNA is almost twice as large as in hydrated lysozyme. We ascribe these differences to a more open and flexible structure of hydrated RNA, and to a larger fraction and different nature of hydrophilic sites. The latter leads to a higher density of water that makes the biomolecule more flexible. All-atom molecular-dynamics simulations are used to show that the extent of hydration is greater in tRNA than in lysozyme. We propose that water acts as a “lubricant” in facilitating enhanced motion in solvated RNA molecules.     

NMR relaxation studies on the hydrate layer of intrinsically unstructured proteins

Intrinsically unstructured/disordered proteins (IUPs) exist in a disordered and largely solvent-exposed, still functional, structural state under physiological conditions. As their function is often directly linked with structural disorder, understanding their structure-function relationship in detail is a great challenge to structural biology. In particular, their hydration and residual structure, both closely linked with their mechanism of action, require close attention. Here we demonstrate that the hydration of IUPs can be adequately approached by a technique so far unexplored with respect to IUPs, solid-state NMR relaxation measurements. This technique provides quantitative information on various features of hydrate water bound to these proteins. By freezing nonhydrate (bulk) water out, we have been able to measure free induction decays pertaining to protons of bound water from which the amount of hydrate water, its activation energy, and correlation times could be calculated. Thus, for three IUPs, the first inhibitory domain of calpastatin, microtubule-associated protein 2c, and plant dehydrin early responsive to dehydration 10, we demonstrate that they bind a significantly larger amount of water than globular proteins, whereas their suboptimal hydration and relaxation parameters are correlated with their differing modes of function. The theoretical treatment and experimental approach presented in this article may have general utility in characterizing proteins that belong to this novel structural class.     

Myosin-induced volume increase of the hyper-mobile water surrounding actin filaments

Microwave dielectric spectroscopy can measure the rotational mobility of water molecules that hydrate proteins and the hydration-shell volume. Using this technique, we have recently shown that apart from typical hydrating water molecules with lowered mobility there are other water molecules around the actin filaments (F-actin) which have a much higher mobility than that of bulk water [Biophys. J. 85 (2003) 3154]. We report here that the volume of this water component (hyper-mobile water) markedly increases without significant change of the volume of the ordinary hydration shell when the myosin motor-domain (S1, myosin subfragment-1) binds to F-actin. No hyper-mobile component was found in the hydration shell of S1 itself. The present results strongly suggest that the solvent space around S1 bound to F-actin is diffusionally asymmetric, which supports our model of force generation by actomyosin proposed previously [op. cit.].     

Interfacial Water at Protein Surfaces: Wide-Line NMR and DSC Characterization of Hydration in Ubiquitin Solutions

Wide-line ¹H-NMR and differential scanning calorimetry measurements were done in aqueous solutions and on lyophilized samples of human ubiquitin between −70°C and +45°C. The measured properties (size, thermal evolution, and wide-line NMR spectra) of the protein-water interfacial region are substantially different in the double-distilled and buffered-water solutions of ubiquitin. The characteristic transition in water mobility is identified as the melting of the nonfreezing/hydrate water. The amount of water in the low-temperature mobile fraction is 0.4 g/g protein for the pure water solution. The amount of mobile water is higher and its temperature dependence more pronounced for the buffered solution. The specific heat of the nonfreezing/hydrate water was evaluated using combined differential scanning calorimetry and NMR data. Considering the interfacial region as an independent phase, the values obtained are 5.0-5.8 J·g⁻¹·K⁻¹, and the magnitudes are higher than that of pure/bulk water (4.2 J·g⁻¹·K⁻¹). This unexpected discrepancy can only be resolved in principle by assuming that hydrate water is in tight H-bond coupling with the protein matrix. The specific heat for the system composed of the protein molecule and its hydration water is 2.3 J·g⁻¹·K⁻¹. It could be concluded that the protein ubiquitin and its hydrate layer behave as a highly interconnected single phase in a thermodynamic sense.     

Hydration of Cytidine, 2′-Deoxycytidine and their Phosphate Salts in the Aqueous Solutions. A Molecular Dynamics Computer Simulation Study

Abstract

To compare the hydration pattern of the cytidine (Cyd) and 2′-deoxycytidine (dCyd) in the aqueous solutions at the level of microscopic interactions, Molecular Dynamics (MD) computer simulations have been undertaken. The results indicate that the hydration of the heterocyclic base moiety in cytidine and 2′-deoxycytidine has a hydrophobic character. None of the three potential Watson—;Crick base pair centres hydrogen bonds with the water molecules and the formation of something akin to a clathrate cage structure of water around base moieties of nucleosides in the aqueous solution is suggested. In contrast, the hydration of Cyd and dCyd sugar moieties shows a hydrophilic character and the three-dimensional networks of H-bonds involving all hydrophilic centres are formed differently around the ribose and 2′-deoxyribose. The sugar hydroxyl groups participate in the hydrogen bonding with water both as H-donor and as H-acceptor. Their donor-acceptor abilities have been evaluated and compared. The coordination numbers, the geometrical data of the first hydration shell, and the number of hydrogen bonds have been calculated. The changes in the pattern of hydration with the increased concentration of nucleosides and upon nucleoside protonation are discussed. The analysis of the pairwise interaction energies are also presented.     

State of water and its diffusion across lipid bilayers: effect of hydration degree

The state of adsorbed water (estimated from the dependence of the shape of the 1H NMR spectrum on the angle between the normal to the bilayers and the direction of the magnetic field) and the diffusion of water molecules in the direction of the normal to the bilayers (estimated by 1H NMR spectroscopy with the impulse gradient of magnetic field) in microscopically oriented dioleoylphosphatidylcholine bilayers have been studied depending on hydration. The dependences of the shape of the NMR spectrum on angle differ qualitatively only at concentrations of water greater and less than the concentration that is achieved upon hydration from saturated vapors chi(eq) (about 23 weight %). At concentrations below chi(eq), all water present in samples enters the hydrate shells of polar "heads" of lipids or is in the state of "rapid exchange" with the water of hydrate shells, with the result that the signal of spin echo for water is observed only in a narrow range of angles close to the "magic angle", 54 degrees C. At concentrations above xhi(eq), the signal of spin echo for water is retained at all orientations, indicating probably that part of water between the bilayers ("quasi-free water") is in the state of a "slow exchange" with water "bound" to polar "heads". It was found that the coefficient of self-diffusion of water across the system of bilayers inversely depends on the degree of hydration, which is described in the Tanner model with consideration of the self-diffusion of water molecules in the hydrophobic moiety of the bilayer. The permeability of the bilayer, the coefficient of distribution of molecules between the water and lipid phases, and the coefficient of self-diffusion of water in the hydrophobic moiety of the bilayer were estimated.     

Hydrophobic tendency of polar group hydration as a major force in type I antifreeze protein recognition

The random network model of water quantitatively describes the different hydration heat capacities of polar and apolar solutes in terms of distortions of the water-water hydrogen bonding angle in the first hydration shell (Gallagher and Sharp, JACS 2003;125:9853). The distribution of this angle in pure water is bimodal, with a low-angle population and high-angle population. Polar solutes increase the high-angle population while apolar solutes increase the low-angle population. The ratio of the two populations quantifies the hydrophobicity of the solute and provides a sensitive measure of water structural distortions. This method of analysis is applied to study hydration of type I thermal hysteresis protein (THP) from winter flounder and three quadruple mutants of four threonine residues at positions 2, 13, 24, and 35. Wild-type and two mutants (VVVV and AAAA) have antifreeze (thermal hysteresis) activity, while the other mutant (SSSS) has no activity. The analysis reveals significant differences in the hydration structure of the ice-binding site. For the SSSS mutant, polar groups have a typical polar-like hydration, that is, more high-angle H-bonds than bulk water. For the wild-type and active mutants, polar groups have unusual, very apolar-like hydration, that is, more low-angle H-bonds than bulk water. This pattern of hydration was seen previously in the structurally distinct type III THPs (Yang & Sharp Biophys Chem 2004;109:137), suggesting for the first time a general mechanism for different THP classes. The specific shape, residue size, and clustering of both polar and apoler groups are essential for an active ice binding surface.     

On the properties of aqueous amide solutions through classical molecular dynamics simulations

The structure and dynamics of infinitely diluted aqueous amide solutions is studied for 13 compounds in the NVT ensemble using classical molecular dynamics simulations. The aim of this work is to provide valuable insights into the effect of amides on liquid water properties in order to understand the amides role in the kinetic inhibition of clathrate hydrate formation in natural gas mixtures. The OPLS-AA forcefield is used to describe the amides, with parameters obtained through fitting of computed B3LYP/6-311++g^{* *} data when not available in the literature, and the SPC-E model is applied for water molecules. Structural properties of the solutions are analyzed via calculated radial distribution functions and dynamic properties are studied with the computed mean square displacements and velocity autocorrelation functions. Most of the studied compounds show a remarkable structuring effect on the surrounding water with strong interactions resulting from hydrogen bonding between solute and solvent molecules. Hydrophobic and hydrophilic synergistic effects influence the amide–water interaction and the properties of the water solvation shells around amides.     

Hydration of RNA base pairs

The hydration patterns around the RNA Watson-Crick and non-Watson-Crick base pairs in crystals are analyzed and described. The results indicate that (i) the base pair hydration is mostly "in-plane"; (ii) eight hydration sites surround the Watson-Crick G-C and A-U base pairs, with five in the deep and three in the shallow groove, an observation which extends the characteristic isostericity of Watson-Crick pairs; (iii) while the hydration around G-C base pairs is well defined, the hydration around A-U base pairs is more diffuse; (iv) the hydration sites close to the phosphate groups are the best defined and the most recurrent ones; (v) a string of water molecules links the two shallow groove 2'-hydroxyl groups, and (vi) the water molecules fit into notches, the size and accessibility of which are almost as important as the number and strength of the hydrophilic groups lining the cavity. Residence times of water molecules at specific hydration sites, inferred from molecular dynamics simulations, are discussed in the light of present data.     

The Stability of Clathrate Hydrates: Temperature Dependence of Dissociation Pressure in Xe and Ar Hydrate

Abstract

The thermodynamic stability of clathrate hydrates I and II encaging xenon or argon has been investigated by examining the temperature dependence of the dissociation pressure. The evaluation of the stability is made based on the generalized van der Waals and Platteeuw theory developed by Tanaka and Kiyohara [J. Chem. Phys. 98, 4098 (1993)]. In the new treatment, the free energy of formation of hydrates in equilibrium with ice is calculated by taking the coupling of the host lattice vibrations with guests into consideration. The predicted dissociation pressures of Xe and Ar hydrates agree well with experiments in higher temperature range. A poor agreement between experiment and calculation for Ar clathrate hydrate at low temperature is improved by the use of a quantum mechanical partition function for a harmonic oscillator in evaluating the free energy difference between ice and empty hydrate.     

Vibrational spectra of deuterated methane and water molecules in structure I clathrate hydrate from ab initio MD simulation

The deuteration of the lattice molecules in clathrate hydrates is a widely used experimental technique to clearly separate the vibrational modes. However, the effect of the deuteration on the vibrational spectra and molecular motions is not fully understood. Since the guest–host coupling may change the vibrational spectra, a detailed analysis of the vibrational spectra of deuterated clathrate hydrate is significant in the understanding of the mechanism of the vibrational shift. In this study, the vibrational spectra of the deuterated methane hydrates were calculated by ab initio molecular dynamics simulation. The intramolecular vibrational frequency of the methane in D₂O lattice and deuterated methane in H₂O lattice was calculated and compared with the pure methane hydrate. The bending, rocking and overtone of the bending mode was also reported. The effect of coupling of the rattling motions of guest and host molecules on the vibrational spectra was revealed.