Molecular Dynamics of Water at the Protein−Solvent Interface

Evaluation of the Coupled Two-Dimensional Main Chain Torsional Potential in Modeling Intrinsically Disordered Proteins

Intrinsically disordered proteins (IDPs) carry out crucial biological functions in essential biological processes of life. Because of the highly dynamic and conformationally heterogeneous nature of the disordered states of IDPs, molecular dynamics simulations are becoming an indispensable tool for the investigation of the conformational ensembles and dynamic properties of IDPs. Nevertheless, there is still no consensus on the most reliable force field in molecular dynamics simulations for IDPs hitherto. In this work, the recently proposed AMBER99SB^2D force field is evaluated in modeling some disordered polypeptides and proteins by checking its ability to reproduce experimental NMR data. The results highlight that when the ildn side-chain corrections are included, AMBER99SB^2D-ildn exhibits reliable results that agree with experiments compared with its predecessors, the AMBER14SB, AMBER99SB, AMBER99SB-ildn, and AMBER99SB^2D force fields, and that decreasing the overall magnitude of protein-protein interactions in favor of protein-water interactions is a key ingredient behind the improvement.

How water mediates the long-range interactions between remote protein molecules

A mechanism of the influence of the presence of one protein molecule on the internal dynamics of another is proposed.

Structure and dynamics of high- and low-density water molecules in the liquid and supercooled regimes

By combining the local structure index with potential energy minimisations we study the local environment of the water molecules for a couple of water models, TIP5P-Ew and SPC/E, in order to characterise low- and high-density "species". Both models show a similar behaviour within the supercooled regime, with two clearly distinguishable populations of unstructured and structured molecules, the fraction of the latter increasing with supercooling. Additionally, for TIP5P-Ew, we find that the structured component vanishes quickly at the normal liquid regime (above the melting temperature). Thus, while SPC/E provides a fraction of structured molecules similar to that found in X-ray experiments, we show that TIP5P-Ew underestimates such value. Moreover, unlike SPC/E, we demonstrate that TIP5P-Ew does not follow the linear dependence of the logarithm of the structured fraction with inverse temperature, as predicted by the two-order parameter model. Finally, we link structure to dynamics by showing that there exists a strong correlation between structural fluctuation and dynamics in the supercooled state with spatial correlations in both static and dynamic quantities.

Dispersed dynamics of solvation in G-quadruplex DNA: comparison of dynamic Stokes shifts of probes in parallel and antiparallel quadruplex structures

G-quadruplex DNA (GqDNA) structures play an important role in many specific cellular functions and are promising anti-tumor targets for small molecules (ligands). Here, we measured the dynamic Stokes shift of a ligand (Hoechst) bound to parallel c-Myc (mPu22) GqDNA over five decades of time from 100 fs to 10 ns, and compared it with the previously reported dynamics of DAPI bound to antiparallel human telomeric (hTelo22) GqDNA (Pal et al 2015 J. Phys. Chem. Lett. 6 1754). Stokes shift data from fluorescence up-conversion and time-correlated single photon counting experiments was combined to cover the broad dynamic range. The results show that the solvation dynamics of Hoechst in parallel mPu22 GqDNA follow a power law relaxation, added to fast 2 ps exponential relaxation, from 100 fs to 10 ns, with only a subtle difference of power law exponents in the two ligand-GqDNA systems (0.06 in Hoechst-mPu22 compared to 0.16 in DAPI-hTelo22). We measured steady-state fluorescence spectra and time-resolved anisotropy decays which confirm the tight binding of Hoechst to parallel mPu22 with a binding constant of ~1  ×  10⁵ M^-1. The molecular docking of Hoechst in parallel GqDNA followed by a 50 ns molecular dynamics (MD) simulation on a Hoechst-GqDNA complex reveals that Hoechst binds to one of the outer G-tetrads by end-stacking near G13 and G4, which is different from the binding site of DAPI inside a groove of antiparallel hTelo22 GqDNA. Reconciling previous experimental and simulation results, we assign the 2 ps component to the hydration dynamics of only weakly perturbed water near mPu22 and the power law relaxation to the coupled motion of water and DNA (i.e. DNA backbone, unpaired bases and loops connecting G-tetrads) which come near the Hoechst inside parallel GqDNA.

Enthalpic and Entropic Contributions to Hydrophobicity

Hydrophobic hydration plays a key role in a vast variety of biological processes, ranging from the formation of cells to protein folding and ligand binding. Hydrophobicity scales simplify the complex process of hydration by assigning a value describing the averaged hydrophobic character to each amino acid. Previously published scales were not able to calculate the enthalpic and entropic contributions to the hydrophobicity directly. We present a new method, based on Molecular Dynamics simulations and Grid Inhomogeneous Solvation Theory, that calculates hydrophobicity from enthalpic and entropic contributions. Instead of deriving these quantities from the temperature dependence of the free energy of hydration or as residual of the free energy and the enthalpy, we directly obtain these values from the phase space occupied by water molecules. Additionally, our method is able to identify regions with specific enthalpic and entropic properties, allowing to identify so-called "unhappy water" molecules, which are characterized by weak enthalpic interactions and unfavorable entropic constraints.

Atomically resolved three-dimensional structures of electrolyte aqueous solutions near a solid surface

Interfacial liquid layers play a central role in a variety of phenomena ranging from friction to molecular recognition. Liquids near a solid surface form an interfacial layer where the molecular structure is different from that of the bulk. Here we report atomic resolution three-dimensional images of electrolyte solutions near a mica surface that demonstrate the existence of three types of interfacial structures. At low concentrations (0.01-1 M), cations are adsorbed onto the mica. The cation layer is topped by a few hydration layers. At higher concentrations, the interfacial layer extends several nanometres into the liquid. It involves the alternation of cation and anion planes. Fluid Density Functional calculations show that water molecules are a critical factor for stabilizing the structure of the interfacial layer. The interfacial layer stabilizes a crystal-like structure compatible with liquid-like ion and solvent mobilities. At saturation, some ions precipitate and small crystals are formed on the mica.

THz frequency spectrum of protein-solvent interaction energy using a recurrence plot-based Wiener-Khinchin method

The dynamics of a protein and the water surrounding it are coupled via nonbonded energy interactions. This coupling can exhibit a complex, nonlinear, and nonstationary nature. The THz frequency spectrum for this interaction energy characterizes both the vibration spectrum of the water hydrogen bond network, and the frequency range of large amplitude modes of proteins. We use a Recurrence Plot based Wiener-Khinchin method RPWK to calculate this spectrum, and the results are compared to those determined using the classical auto-covariance-based Wiener-Khinchin method WK. The frequency spectra for the total nonbonded interaction energy extracted from molecular dynamics simulations between the β-Lactamase Inhibitory Protein BLIP, and water molecules within a 10 Å distance from the protein surface, are calculated at 150, 200, 250, and 310 K, respectively. Similar calculations are also performed for the nonbonded interaction energy between the residues 49ASP, 53TYR, and 142PHE in BLIP, with water molecules within 10 Å from each residue respectively at 150, 200, 250, and 310 K. A comparison of the results shows that RPWK performs better than WK, and is able to detect some frequency data points that WK fails to detect. This points to the importance of using methods capable of taking the complex nature of the protein-solvent energy landscape into consideration, and not to rely on standard linear methods. In general, RPWK can be a valuable addition to the analysis tools for protein molecular dynamics simulations. Proteins 2016; 84:1549-1557. © 2016 Wiley Periodicals, Inc.

Hydration of proteins and nucleic acids: Advances in experiment and theory. A review

BACKGROUND

Most biological processes involve water, and the interactions of biomolecules with water affect their structure, function and dynamics.

SCOPE OF REVIEW

This review summarizes the current knowledge of protein and nucleic acid interactions with water, with a special focus on the biomolecular hydration layer. Recent developments in both experimental and computational methods that can be applied to the study of hydration structure and dynamics are reviewed, including software tools for the prediction and characterization of hydration layer properties.

MAJOR CONCLUSIONS

In the last decade, important advances have been made in our understanding of the factors that determine how biomolecules and their aqueous environment influence each other. Both experimental and computational methods contributed to the gradually emerging consensus picture of biomolecular hydration.

GENERAL SIGNIFICANCE

An improved knowledge of the structural and thermodynamic properties of the hydration layer will enable a detailed understanding of the various biological processes in which it is involved, with implications for a wide range of applications, including protein-structure prediction and structure-based drug design.

Water Determines the Structure and Dynamics of Proteins

Water is an essential participant in the stability, structure, dynamics, and function of proteins and other biomolecules. Thermodynamically, changes in the aqueous environment affect the stability of biomolecules. Structurally, water participates chemically in the catalytic function of proteins and nucleic acids and physically in the collapse of the protein chain during folding through hydrophobic collapse and mediates binding through the hydrogen bond in complex formation. Water is a partner that slaves the dynamics of proteins, and water interaction with proteins affect their dynamics. Here we provide a review of the experimental and computational advances over the past decade in understanding the role of water in the dynamics, structure, and function of proteins. We focus on the combination of X-ray and neutron crystallography, NMR, terahertz spectroscopy, mass spectroscopy, thermodynamics, and computer simulations to reveal how water assist proteins in their function. The recent advances in computer simulations and the enhanced sensitivity of experimental tools promise major advances in the understanding of protein dynamics, and water surely will be a protagonist.

Effect of temperature on the structure and hydration layer of TATA-box DNA: A molecular dynamics simulation study

DNA within the living cells experiences a diverse range of temperature, ranging from freezing condition to hot spring water. How the structure, the mechanical properties of DNA, and the solvation dynamics around DNA changes with the temperature is important to understand the functionality of DNA under those acute temperature conditions. In that notion, we have carried out molecular dynamics simulations of a DNA oligomer, containing TATA-box sequence for three different temperatures (250K, 300K and 350K). We observed that the structure of the DNA, in terms of backbone torsion angles, sugar pucker, base pair parameters, and base pair step parameters, did not show any unusual properties within the studied range of temperatures, but significant structural alteration was noticed between BI and BII forms at higher temperature. As expected, the flexibility of the DNA, in terms of the torsional rigidity and the bending rigidity is highly temperature dependent, confirming that flexibility increases with increase in temperature. Additionally, the groove widths of the studied DNA showed temperature sensitivity, specifically, the major groove width decreases and the minor groove width increases, respectively, with the increase in temperature. We observed that at higher temperature, water around both the major and the minor groove of the DNA is less structured. However, the water dynamics around the minor groove of the DNA is more restricted as compared to the water around the major groove throughout the studied range of temperatures, without any anomalous behavior.

Effects of protein-DNA complex formation on the intermolecular vibrational density of states of interfacial water

Single-stranded DNAs (ss-DNAs) are formed as intermediates during DNA metabolic processes. ss-DNA binding (SSB) proteins specifically bind to the single-stranded segments of the DNA and protect it from being degraded. We have performed room temperature molecular dynamics simulations of the aqueous solution of two DNA binding K homology (KH) domains (KH3 and KH4) of the far upstream element (FUSE) binding protein (FBP) complexed with two ss-DNA oligomers. Efforts have been made to explore the influence of complex formation on low-frequency vibrational density of states of the surface water molecules. It is revealed that increased back scattering of water confined around the complexed structures leads to significant blue shifts of the band corresponding to the O···O···O bending or restricted transverse motions of water, the effect being more for the bridged water molecules. Importantly, it is demonstrated that the formation of such complexed structures of a similar type may often influence the transverse and longitudinal degrees of freedom of the surrounding water molecules in a nonuniform manner.

Impact of secondary structure and hydration water on the dielectric spectrum of poly-alanine and possible relation to the debate on slaved versus slaving water

Using extensive molecular dynamics simulations of a single eight-residue alanine polypeptide in explicit water, we investigate the influence of α-helix formation on the dielectric spectrum. For this, we project long equilibrium trajectories into folded and unfolded states and thereby obtain dielectric spectra representative for disordered as well α-helical conformations without the need to change any other system parameter such as pH or temperature. The absorption spectrum in the α-helical state exhibits a feature in the sub-GHz range that is significantly stronger than in the unfolded state. As we show by an additional decomposition into peptide and water contributions, this slow dielectric mode, the relaxation time of which matches the independently determined peptide rotational relaxation time, is mostly caused by peptide polarization correlations, but also contains considerable contributions from peptide-water correlations. In contrast, the peptide spectral contribution shows no features in the GHz range where bulk water absorbs, not even in the peptide-water correlation part, we conclude that hydration water around Ala8 is more influenced by peptide polarization relaxation effects than the other way around. A further decomposition into water-self and water-collective polarization correlations shows that the dielectric response of hydration water is, in contrast to electrolyte solutions, retarded and that this retardation is mostly due to collective effects, the self relaxation of hydration water molecules is only slightly slowed down compared to bulk water. We find the dynamic peptide-water polarization cross correlations to be rather long-ranged and to extend more than one nanometer away from the peptide-water interface into the water hydration shell, in qualitative agreement with previous simulation studies and recent THz absorption experiments.

Albumin (BSA) Adsorption over Graphene in Aqueous Environment: Influence of Orientation, Adsorption Protocol, and Solvent Treatment

We report 150 ns explicit solvent MD simulations of the adsorption on graphene of albumin (BSA) in two orientations and using two different adsorption protocols, i.e., free and forced adsorption. Our results show that free adsorption occurs with little structural rearrangements. Even taking adsorption to an extreme, by forcing it with a 5 nN downward force applied during the initial 20 ns, we show that along a particular orientation BSA is able to preserve the structural properties of the majority of its binding sites. Furthermore, in all the cases considered in this work, the ibuprofen binding site has shown a strong resilience to structural changes. Finally, we compare these results with implicit solvent simulations and find that the latter predicts an extreme protein unfolding upon adsorption. The origin of this discrepancy is attributed to a poor description of the water entropic forces at interfaces in the implicit solvent methods.

Microscopic dynamics of water around unfolded structures of barstar at room temperature

The breaking of the native structure of a protein and its influences on the dynamic response of the surrounding solvent is an important issue in protein folding. In this work, we have carried out atomistic molecular dynamics simulations to unfold the protein barstar at two different temperatures (400 K and 450 K). The two unfolded forms obtained at such high temperatures are further studied at room temperature to explore the effects of nonuniform unfolding of the protein secondary structures along two different pathways on the microscopic dynamical properties of the surface water molecules. It is demonstrated that though the structural transition of the protein in general results in less restricted water motions around its segments, but there are evidences of formation of new conformational motifs upon unfolding with increasingly confined environment around them, thereby resulting in further restricted water mobility in their hydration layers. Moreover, it is noticed that the effects of nonuniform unfolding of the protein segments on the relaxation times of the protein-water (PW) and the water-water (WW) hydrogen bonds are correlated with hindered hydration water motions. However, the kinetics of breaking and reformation of such hydrogen bonds are found to be influenced differently at the interface. It is observed that while the effects of unfolding on the PW hydrogen bond kinetics seem to be minimum, but the kinetics involving the WW hydrogen bonds around the protein segments exhibit noticeably heterogeneous characteristics. We believe that this is an important observation, which can provide valuable insights on the origin of heterogeneous influence of unfolding of a protein on the microscopic properties of its hydration water.

Signatures of protein thermal denaturation and local hydrophobicity in domain specific hydration behavior: a comparative molecular dynamics study

We investigate, using atomistic molecular dynamics simulations, the association of surface hydration accompanying local unfolding in the mesophilic protein Yfh1 under a series of thermal conditions spanning its cold and heat denaturation temperatures. The results are benchmarked against the thermally stable protein, Ubq, and behavior at the maximum stability temperature. Local unfolding in Yfh1, predominantly in the beta sheet regions, is in qualitative agreement with recent solution NMR studies; the corresponding Ubq unfolding is not observed. Interestingly, all domains, except for the beta sheet domains of Yfh1, show increased effective surface hydrophobicity with increase in temperature, as reflected by the density fluctuations of the hydration layer. Velocity autocorrelation functions (VACF) of oxygen atoms of water within the hydration layers and the corresponding vibrational density of states (VDOS) are used to characterize alteration in dynamical behavior accompanying the temperature dependent local unfolding. Enhanced caging effects accompanying transverse oscillations of the water molecules are found to occur with the increase in temperature preferentially for the beta sheet domains of Yfh1. Helical domains of both proteins exhibit similar trends in VDOS with changes in temperature. This work demonstrates the existence of key signatures of the local onset of protein thermal denaturation in solvent dynamical behavior.

Sub-diffusion and population dynamics of water confined in soft environments

We have studied by using molecular dynamics computer simulations the dynamics of water confined in ionic surfactant phases, ranging from well ordered lamellar structures to micelles at low and high water loading, respectively. We have analysed in depth the main dynamical features in terms of mean-squared displacements and intermediate scattering functions, and found clear evidence of sub-diffusive behaviour. We have identified water molecules lying at the charged interface with the hydrophobic confining matrix as the main factor responsible for this unusual feature, and given a comprehensive picture of dynamics based on a very precise analysis of lifetimes at the interface. We conclude by providing, for the first time to our knowledge, a unique framework for rationalizing the existence of important dynamical heterogeneities in fluids adsorbed in soft confining environments.

Elucidating the fundamental forces in protein crystal formation: the case of crambin

This study demonstrates the feasibility of periodic all-electron hybrid density functional theory calculations in the description of protein crystals, using crambin as a test case.

Molecular simulations of proteins have been usually accomplished through empirical or semi-empirical potentials, due to the large size and inherent complexity of these biological systems. On the other hand, a theoretical description of proteins based on quantum-mechanical methods would however provide an unbiased characterization of their electronic properties, possibly offering a link between these and the ultimate biological activity. Yet, such approaches have been historically hindered by the large amount of requested computational power. Here we demonstrate the feasibility of periodic all-electron density functional theory calculations in the description of the crystal of the protein crambin (46 aminoacids), which is determined with exceptional structural accuracy. We have employed the hybrid B3LYP functional, coupled to an empirical description of London interactions (D*) to simulate the crambin crystal with an increasing amount of lattice water molecules in the cell (up to 172H₂O per cell). The agreement with the experiment is good for both protein geometry and protein–water interactions. The energetics was computed to predict crystal formation energies, protein–water and protein–protein interaction energies. We studied the role of dispersion interactions which are crucial for holding the crambin crystal in place. B3LYP-D* electrostatic potential and dipole moment of crambin as well as the electronic charge flow from crambin to the solvating water molecules (0.0015e per H₂O) have also been predicted. These results proved that quantum-mechanical simulations of small proteins, both free and in their crystalline state, are now feasible in a reasonable amount of time, by programs capable of exploiting high performance computing architectures, allowing the study of protein properties not easily amenable through classical force fields.

Diffusion of Hydration Water around Intrinsically Disordered Proteins

Hydration water dynamics around globular proteins have attracted considerable attention in the past decades. This work investigates the hydration water dynamics around partially/fully intrinsically disordered proteins and compares it to that of the globular proteins via molecular dynamics simulations. The translational diffusion of the hydration water is examined by evaluating the mean-square displacement and the velocity autocorrelation function, while the rotational diffusion is probed through the dipole-dipole time correlation function. The results reveal that the translational and rotational motions of water molecules at the surface of intrinsically disordered proteins/regions are less restricted as compared to those around globular proteins/ordered regions, which is reflected in their higher diffusion coefficient and lower orientational relaxation time. The restricted mobility of hydration water in the vicinity of the protein leads to a sublinear diffusion in a heterogeneous interface. A positive correlation between the mean number of hydrogen bonds and the diffusion coefficient of hydration water implies higher mobility of water molecules at the surface of disordered proteins, which is due to their higher number of hydrogen bonds. Enhanced hydration water mobility around disordered proteins/regions is also related to their higher hydration capacity, low hydrophobicity, and increased internal protein motions. Thus, we generalize that the intrinsically disordered proteins/regions are associated with higher hydration water mobility as compared to globular protein/ordered regions, which may help to elucidate their varied functional specificity.

Structural and functional stabilization of protein entities: state-of-the-art

Within the context of biomedicine and pharmaceutical sciences, the issue of (therapeutic) protein stabilization assumes particular relevance. Stabilization of protein and protein-like molecules translates into preservation of both structure and functionality during storage and/or targeting, and such stabilization is mostly attained through establishment of a thermodynamic equilibrium with the (micro)environment. The basic thermodynamic principles that govern protein structural transitions and the interactions of the protein molecule with its (micro)environment are, therefore, tackled in a systematic fashion. Highlights are given to the major classes of (bio)therapeutic molecules, viz. enzymes, recombinant proteins, (macro)peptides, (monoclonal) antibodies and bacteriophages. Modification of the microenvironment of the biomolecule via multipoint covalent attachment onto a solid surface followed by hydrophilic polymer co-immobilization, or physical containment within nanocarriers, are some of the (latest) strategies discussed aiming at full structural and functional stabilization of said biomolecules.

Comparative study of hydration shell dynamics around a hyperactive antifreeze protein and around ubiquitin

The hydration layer surrounding a protein plays an essential role in its biochemical function and consists of a heterogeneous ensemble of water molecules with different local environments and different dynamics. What determines the degree of dynamical heterogeneity within the hydration shell and how this changes with temperature remains unclear. Here, we combine molecular dynamics simulations and analytic modeling to study the hydration shell structure and dynamics of a typical globular protein, ubiquitin, and of the spruce budworm hyperactive antifreeze protein over the 230-300 K temperature range. Our results show that the average perturbation induced by both proteins on the reorientation dynamics of water remains moderate and changes weakly with temperature. The dynamical heterogeneity arises mostly from the distribution of protein surface topographies and is little affected by temperature. The ice-binding face of the antifreeze protein induces a short-ranged enhancement of water structure and a greater slowdown of water reorientation dynamics than the non-ice-binding faces whose effect is similar to that of ubiquitin. However, the hydration shell of the ice-binding face remains less tetrahedral than the bulk and is not "ice-like". We finally show that the hydrogen bonds between water and the ice-binding threonine residues are particularly strong due to a steric confinement effect, thereby contributing to the strong binding of the antifreeze protein on ice crystals.

Weighted protein residue networks based on joint recurrences between residues

BACKGROUND

Weighted and un-weighted protein residue networks can predict key functional residues in proteins based on the closeness centrality C and betweenness centrality B values for each residue. A static snapshot of the protein structure, and a cutoff distance, are used to define edges between the network nodes. In this work we apply the weighted network approach to study the β-Lactamase Inhibitory Protein (BLIP). Joint recurrences extracted from molecular dynamics MD trajectory positions of the protein residue carbon alpha atoms are used to define edge weights between nodes, and no cutoff distance is used. The results for B and C from our approach are compared with those extracted from an un-weighted network, and a weighted network that uses interatomic contacts to define edge weights between nodes, respectively.

RESULTS

The joint recurrence weighted network approach performs well in pointing out key protein residues. Furthermore, it seems to emphasize residues with medium to high relative solvent accessibility that lie in loop regions between secondary structure elements of the protein.

CONCLUSIONS

Protein residue networks that use joint recurrences extracted from molecular dynamics simulations of a solvated protein perform well in pointing to hotspot residues and hotspot clusters. This approach uses no distance cutoff threshold, and does not exclude any interactions between the residues, including water-mediated interactions.

Power-Law Solvation Dynamics in G-Quadruplex DNA: Role of Hydration Dynamics on Ligand Solvation inside DNA

G-quadruplex DNA (GqDNA) structures act as promising anticancer targets for small-molecules (ligands). Solvation dynamics of a ligand (DAPI: 4',6-diamidino-2-phenylindole) inside antiparallel-GqDNA is studied through direct comparison of time-resolved experiments to molecular dynamics (MD) simulation. Dynamic Stokes shifts of DAPI in GqDNA prepared in H2O buffer and D2O are compared to find the effect of water on ligand solvation. Experimental dynamics (in H2O) is then directly compared with the dynamics computed from 65 ns simulation on the same DAPI-GqDNA complex. Ligand solvation follows power-law relaxation (summed with fast exponential relaxation) from ~100 fs to 10 ns. Simulation results show relaxation below ~5 ps is dominated by water motion, while both water and DNA contribute comparably to dictate long-time power-law dynamics. Ion contribution is, however, found to be negligible. Simulation results also suggest that anomalous solvation dynamics may have origin in subdiffusive motion of perturbed water near GqDNA.

Concentration dependence of hydration water in a model peptide

The molecular dynamics of aqueous solutions of a model amphiphilic peptide is studied as a function of concentration by broad-band light scattering experiments. Similarly to protein aqueous solutions, a considerable retardation, of about a factor 6-8, of hydration water dynamics with respect to bulk water is found, showing a slight dependence on solute concentration. Conversely, the average number of water molecules perturbed by the presence of peptide, i.e. the hydration number, appears to be strongly modified by adding solute. Its behaviour, decreasing upon increasing concentration, can be interpreted considering the random close-to-contact condition experienced by solute particles. Overall, the present findings support the view of a "long range" effect of peptides on the surrounding water, extending beyond the first two hydration shells.

Changes of water hydrogen bond network with different externalities

It is crucial to uncover the mystery of water cluster and structural motif to have an insight into the abundant anomalies bound to water. In this context, the analysis of influence factors is an alternative way to shed light on the nature of water clusters. Water structure has been tentatively explained within different frameworks of structural models. Based on comprehensive analysis and summary of the studies on the response of water to four externalities (i.e., temperature, pressure, solutes and external fields), the changing trends of water structure and a deduced intrinsic structural motif are put forward in this work. The variations in physicochemical and biological effects of water induced by each externality are also discussed to emphasize the role of water in our daily life. On this basis, the underlying problems that need to be further studied are formulated by pointing out the limitations attached to current study techniques and to outline prominent studies that have come up recently.

Water confined in self-assembled ionic surfactant nano-structures

We present a coarse-grained model for ionic surfactants in explicit aqueous solutions, and study by computer simulation both the impact of water content on the morphology of the system, and the consequent effect of the formed interfaces on the structural features of the absorbed fluid. On increasing the hydration level under ambient conditions, the model exhibits a series of three distinct phases: lamellar, cylindrical and micellar. We characterize the different structures in terms of diffraction patterns and neutron scattering static structure factors. We demonstrate that the rate of variation of the nano-metric sizes of the self-assembled water domains shows peculiar changes in the different phases. We also analyse in depth the structure of the water/confining matrix interfaces, the implications of their tunable degree of curvature, and the properties of water molecules in different restricted environments. Finally, we compare our results with experimental data and their impact on a wide range of important scientific and technological domains, where the behavior of water at the interface with soft materials is crucial.

Boundary condition effects on the dynamic and electric properties of hydration layers

Water solvation has a central role in several biochemical processes ranging from protein folding to biomolecular recognition and enzyme catalysis. Because of its importance, the structure and dynamics of hydration layers around biological macromolecules have been the targets of a great number of experimental and computational studies. In the present contribution, we have investigated the effects of periodic boundary conditions (PBCs), as used in conjunction with molecular dynamics (MD) simulations, on the dynamic and electric properties of water layers. In particular, we have systematically performed MD simulations of neat water and biomolecules in aqueous solutions by imposing a different external dielectric constant, a generally overlooked parameter in PBC simulations. The effect of the system size has also been addressed. Overall, our results consistently indicate that the dipole moment properties of water layers, and specifically the dipole moment fluctuations and the reorientational correlation functions, can be sensitive to the choice of the external boundary conditions, whereas other molecular properties, such as the self-diffusion coefficient and the reorientational relaxation times, are not affected. We think that our investigation may help to assess appropriate simulation conditions for modeling the aqueous environment of relevant biochemical systems and processes.

Wetting behavior of nonpolar nanotubes in simple dipolar liquids for varying nanotube diameter and solute-solvent interactions

Atomistic simulations of model nonpolar nanotubes in a Stockmayer liquid are carried out for varying nanotube diameter and nanotube-solvent interactions to investigate solvophobic interactions in generic dipolar solvents. We have considered model armchair type single-walled nonpolar nanotubes with increasing radii from (5,5) to (12,12). The interactions between solute and solvent molecules are modeled by the well-known Lennard-Jones and repulsive Weeks-Chandler-Andersen potentials. We have investigated the density profiles and microscopic arrangement of Stockmayer molecules, orientational profiles of their dipole vectors, time dependence of their occupation, and also the translational and rotational motion of solvent molecules in confined environments of the cylindrical nanopores and also in their external peripheral regions. The present results of structural and dynamical properties of Stockmayer molecules inside and near atomistically rough nonpolar surfaces including their wetting and dewetting behavior for varying interactions provide a more generic picture of solvophobic effects experienced by simple dipolar liquids without any specific interactions such as hydrogen bonds.

Systematic evaluation of bundled SPC water for biomolecular simulations

How accurate is bundled SPC water as inner shell solvent for hybrid all-atom/coarse-grained simulations?

Hydration and Nanoconfined Water: Insights from Computer Simulations

The comprehension of the structure and behavior of water at interfaces and under nanoconfinement represents an issue of major concern in several central research areas like hydration, reaction dynamics and biology. From one side, water is known to play a dominant role in the structuring, the dynamics and the functionality of biological molecules, governing main processes like protein folding, protein binding and biological function. In turn, the same principles that rule biological organization at the molecular level are also operative for materials science processes that take place within a water environment, being responsible for the self-assembly of molecular structures to create synthetic supramolecular nanometrically-sized materials. Thus, the understanding of the principles of water hydration, including the development of a theory of hydrophobicity at the nanoscale, is imperative both from a fundamental and an applied standpoint. In this work we present some molecular dynamics studies of the structure and dynamics of water at different interfaces or confinement conditions, ranging from simple model hydrophobic interfaces with different geometrical constraints (in order to single out curvature effects), to self-assembled monolayers, proteins and phospholipid membranes. The tendency of the water molecules to sacrifice the lowest hydrogen bond (HB) coordination as possible at extended interfaces is revealed. This fact makes the first hydration layers to be highly oriented, in some situations even resembling the structure of hexagonal ice. A similar trend to maximize the number of HBs is shown to hold in cavity filling, with small subnanometric hydrophobic cavities remaining empty while larger cavities display an alternation of filled and dry states with a significant inner HB network. We also study interfaces with complex chemical and geometrical nature in order to determine how different conditions affect the local hydration properties. Thus, we show some results for protein hydration and, particularly, some preliminary studies on membrane hydration. Finally, calculations of a local hydrophobicity measure of relevance for binding and self-assembly are also presented. We then conclude with a few words of further emphasis on the relevance of this kind of knowledge to biology and to the design of new materials by highlighting the context-dependent and non-additive nature of different non-covalent interactions in an aqueous nanoenvironment, an issue that is usually greatly overlooked.

The non-uniform early structural response of globular proteins to cold denaturing conditions: a case study with Yfh1

The mechanism of cold denaturation in proteins is often incompletely understood due to limitations in accessing the denatured states at extremely low temperatures. Using atomistic molecular dynamics simulations, we have compared early (nanosecond timescale) structural and solvation properties of yeast frataxin (Yfh1) at its temperature of maximum stability, 292 K (Ts), and the experimentally observed temperature of complete unfolding, 268 K (Tc). Within the simulated timescales, discernible "global" level structural loss at Tc is correlated with a distinct increase in surface hydration. However, the hydration and the unfolding events do not occur uniformly over the entire protein surface, but are sensitive to local structural propensity and hydrophobicity. Calculated infrared absorption spectra in the amide-I region of the whole protein show a distinct red shift at Tc in comparison to Ts. Domain specific calculations of IR spectra indicate that the red shift primarily arises from the beta strands. This is commensurate with a marked increase in solvent accessible surface area per residue for the beta-sheets at Tc. Detailed analyses of structure and dynamics of hydration water around the hydrophobic residues of the beta-sheets show a more bulk water like behavior at Tc due to preferential disruption of the hydrophobic effects around these domains. Our results indicate that in this protein, the surface exposed beta-sheet domains are more susceptible to cold denaturing conditions, in qualitative agreement with solution NMR experimental results.