Molecular view of water dynamics near model peptides

Understanding the intricacy of protein in hydrated deep eutectic solvent: Solvation dynamics, conformational fluctuation dynamics, and stability

Deep eutectic solvents (DESs) are potential biocatalytic media due to their easy preparation, fine-tuneability, biocompatibility, and most importantly, due to their ability to keep protein stable and active. However, there are many unanswered questions and gaps in our knowledge about how proteins behave in these alternate media. Herein, we investigated solvation dynamics, conformational fluctuation dynamics, and stability of human serum albumin (HSA) in 0.5 Acetamide/0.3 Urea/0.2 Sorbitol (0.5Ac/0.3Ur/0.2Sor) DES of varying concentrations to understand the intricacy of protein behaviour in DES. Our result revealed a gradual decrease in the side-chain flexibility and thermal stability of HSA beyond 30 % DES. On the other hand, the associated water dynamics around domain-I of HSA decelerate only marginally with increasing DES content, although viscosity rises considerably. We propose that even though macroscopic solvent properties are altered, a protein feels only an aqueous type of environment in the presence of DES. This is probably the first experimental study to delineate the role of the associated water structure of the enzyme for maintaining its stability inside DES. Although considerable effort is necessary to generalize such claims, it might serve as the basis for understanding why proteins remain stable and active in DES.

The Hydrophobic Effect Studied by Using Interacting Colloidal Suspensions

Interactions between nanoparticles (NPs) determine their self-organization and dynamic processes. In these systems, a quantitative description of the interparticle forces is complicated by the presence of the hydrophobic effect (HE), treatable only qualitatively, and due to the competition between the hydrophobic and hydrophilic forces. Recently, instead, a sort of crossover of HE from hydrophilic to hydrophobic has been experimentally observed on a local scale, by increasing the temperature, in pure confined water and studying the occurrence of this crossover in different water-methanol solutions. Starting from these results, we then considered the idea of studying this process in different nanoparticle solutions. By using photon correlation spectroscopy (PCS) experiments on dendrimer with OH terminal groups (dissolved in water and methanol, respectively), we show the existence of this hydrophobic-hydrophilic crossover with a well defined temperature and nanoparticle volume fraction dependence. In this frame, we have used the mode coupling theory extended model to evaluate the measured time-dependent density correlation functions (ISFs). In this context we will, therefore, show how the measured spectra are strongly dependent on the specificity of the interactions between the particles in solution. The observed transition demonstrates that just the HE, depending sensitively on the system thermodynamics, determines the hydrophobic and hydrophilic interaction properties of the studied nanostructures surface.

Amide Spectral Fingerprints are Hydrogen Bonding-Mediated

The origin of the peculiar amide spectral features of proteins in aqueous solution is investigated, by exploiting a combined theoretical and experimental approach to study UV Resonance Raman (RR) spectra of peptide molecular models, namely N-acetylglycine-N-methylamide (NAGMA) and N-acetylalanine-N-methylamide (NALMA). UVRR spectra are recorded by tuning Synchrotron Radiation at several excitation wavelengths and modeled by using a recently developed multiscale protocol based on a polarizable QM/MM approach. Thanks to the unparalleled agreement between theory and experiment, we demonstrate that specific hydrogen bond interactions, which dominate hydration dynamics around these solutes, play a crucial role in the selective enhancement of amide signals. These results further argue the capability of vibrational spectroscopy methods as valuable tools for refined structural analysis of peptides and proteins in aqueous solution.

Hydration Dynamics of Model Peptides with Different Hydrophobic Character

The multi-scale dynamics of aqueous solutions of the hydrophilic peptide N-acetyl-glycine-methylamide (NAGMA) have been investigated through extended frequency-range depolarized light scattering (EDLS), which enables the broad-band detection of collective polarizability anisotropy fluctuations. The results have been compared to those obtained for N-acetyl-leucinemethylamide (NALMA), an amphiphilic peptide which shares with NAGMA the same polar backbone, but also contains an apolar group. Our study indicates that the two model peptides induce similar effects on the fast translational dynamics of surrounding water. Both systems slow down the mobility of solvating water molecules by a factor 6–8, with respect to the bulk. Moreover, the two peptides cause a comparable far-reaching spatial perturbation extending to more than two hydration layers in diluted conditions. The observed concentration dependence of the hydration number is explained considering the random superposition of different hydration shells, while no indication of solute aggregation phenomena has been found. The results indicate that the effect on the dynamics of water solvating the amphiphilic peptide is dominated by the hydrophilic backbone. The minor impact of the hydrophobic moiety on hydration features is consistent with structural findings derived by Fourier transform infrared (FTIR) measurements, performed in attenuated total reflectance (ATR) configuration. Additionally, we give evidence that, for both systems, the relaxation mode in the GHz frequency range probed by EDLS is related to solute rotational dynamics. The rotation of NALMA occurs at higher timescales, with respect to the rotation of NAGMA; both processes are significantly slower than the structural dynamics of hydration water, suggesting that solute and solvent motions are uncoupled. Finally, our results do not indicate the presence of super-slow water (relaxation times in the order of tens of picoseconds) around the peptides investigated.

Sequence Modulates Polypeptoid Hydration Water Structure and Dynamics

We use molecular dynamics simulations to investigate the effect of polypeptoid sequence on the structure and dynamics of its hydration waters. Polypeptoids provide an excellent platform to study small-molecule hydration in disordered polymers, as they can be precisely synthesized with a variety of sidechain chemistries. We examine water behavior near a set of peptoid oligomers in which the number and placement of nonpolar versus polar sidechains are systematically varied. To do this, we leverage a new computational workflow enabling accurate sampling of polypeptoid conformations. We find that the hydration waters are less dense, are more tetrahedral, and have slower dynamics compared to bulk water. The magnitude of these shifts increases with the number of nonpolar groups. We also find that shifts in the water structure and dynamics are strongly correlated, suggesting that experimental insight into the dynamics of hydration water obtained by Overhauser dynamic nuclear polarization (ODNP) also contains information about water structural properties. We then demonstrate the ability of ODNP to probe site-specific dynamics of hydration water near these model peptoid systems.

Hydration of Simple Model Peptides in Aqueous Osmolyte Solutions

The biology and chemistry of proteins and peptides are inextricably linked with water as the solvent. The reason for the high stability of some proteins or uncontrolled aggregation of others may be hidden in the properties of their hydration water. In this study, we investigated the effect of stabilizing osmolyte–TMAO (trimethylamine N-oxide) and destabilizing osmolyte–urea on hydration shells of two short peptides, NAGMA (N-acetyl-glycine-methylamide) and diglycine, by means of FTIR spectroscopy and molecular dynamics simulations. We isolated the spectroscopic share of water molecules that are simultaneously under the influence of peptide and osmolyte and determined the structural and energetic properties of these water molecules. Our experimental and computational results revealed that the changes in the structure of water around peptides, caused by the presence of stabilizing or destabilizing osmolyte, are significantly different for both NAGMA and diglycine. The main factor determining the influence of osmolytes on peptides is the structural-energetic similarity of their hydration spheres. We showed that the chosen peptides can serve as models for various fragments of the protein surface: NAGMA for the protein backbone and diglycine for the protein surface with polar side chains.

Dynamics in supramolecular nanomaterials

Self-assembly of amphiphilic small molecules in water leads to nanostructures with customizable structure-property relationships arising from their tunable chemistries. Characterization of these assemblies is generally limited to their static structures -e.g. their geometries and dimensions - but the implementation of tools that provide a deeper understanding of molecular motions has recently emerged. Here, we summarize recent reports showcasing dynamics characterization tools and their application to small molecule assemblies, and we go on to highlight supramolecular systems whose properties are substantially affected by their conformational, exchange, and water dynamics. This review illustrates the importance of considering dynamics in rational amphiphile design.

Behavior of Water Near Multimodal Chromatography Ligands and Its Consequences for Modulating Protein-Ligand Interactions

Multimodal chromatography is a powerful approach for purifying proteins that uses ligands containing multiple modes of interaction. Recent studies have shown that selectivity in multimodal chromatographic separations is a function of the ligand structure and geometry. Here, we performed molecular dynamics simulations to explore how the ligand structure and geometry affect ligand-water interactions and how these differences in solution affect the nature of protein-ligand interactions. Our investigation focused on three chromatography ligands: Capto MMC, Nuvia cPrime, and Prototype 4, a structural variant of Nuvia cPrime. First, the solvation characteristics of each ligand were quantified via three metrics: average water density, fluctuations, and residence time. We then explored how solvation was perturbed when the ligand was bound to the protein surface and found that the probability of the phenyl ring dewetting followed the order: Capto MMC > Prototype 4 > Nuvia cPrime. To explore how these differences in dewetting affect protein-ligand interactions, we calculated the probability of each ligand binding to different types of residues on the protein surface and found that the probability of binding to a hydrophobic residue followed the same order as the dewetting behavior. This study illustrates the role that wetting and dewetting play in modulating protein-ligand interactions.

First-passage fingerprints of water diffusion near glutamine surfaces

The extent to which biological interfaces affect the dynamics of water plays a key role in the exchange of matter and chemical interactions that are essential for life. The density and the mobility of water molecules depend on their proximity to biological interfaces and can play an important role in processes such as protein folding and aggregation. In this work, we study the dynamics of water near glutamine surfaces-a system of interest in studies of neurodegenerative diseases. Combining molecular-dynamics simulations and stochastic modelling, we study how the mean first-passage time and related statistics of water molecules escaping subnanometer-sized regions vary from the interface to the bulk. Our analysis reveals a dynamical complexity that reflects underlying chemical and geometrical properties of the glutamine surfaces. From the first-passage time statistics of water molecules, we infer their space-dependent diffusion coefficient in directions normal to the surfaces. Interestingly, our results suggest that the mobility of water varies over a longer length scale than the chemical potential associated with the water-protein interactions. The synergy of molecular dynamics and first-passage techniques opens the possibility for extracting space-dependent diffusion coefficients in more complex, inhomogeneous environments that are commonplace in living matter.

Structural and dynamical heterogeneities at glutamine-water interfaces

The behavior of water at the surfaces of solid amino acid crystals has received little attention despite its importance in nucleation processes. In this work, we take a first step to fill this gap by using molecular dynamics simulations to study the structural and dynamical properties of water near the (100), (010) and (001) surfaces of l-glutamine crystals. These highly hydrophilic surfaces serve as excellent model systems for interrogating the behavior of water. Despite having the same molecular composition, water at each surface displays characteristic structural, orientational and dynamical correlations. This behavior is tuned by how the different chemical groups of amino acids make contact with the liquid phase. All three surfaces yield a glassy layer of interfacial water which is reflected in different ways such as the presence of a rotationally arrested layer of water molecules and substantial slow down of the diffusion of water near the interface. By increasing the concentration of molecules in solution, we show that the binding of glutamine molecules to the crystal surface creates a crowded environment involving pockets of trapped water molecules altering the water dynamics in a highly non-trivial manner suggesting that the solvent dynamics may have important implications on crystal nucleation.

Dynamics Properties of Photosynthetic Microorganisms Probed by Incoherent Neutron Scattering

Studies on the dynamical properties of photosynthetic membranes of land plants and purple bacteria have been previously performed by neutron spectroscopy, revealing a tight coupling between specific photochemical reactions and macromolecular dynamics. Here, we probed the intrinsic dynamics of biotechnologically useful mutants of the green alga Chlamydomonas reinhardtii by incoherent neutron scattering coupled with prompt chlorophyll fluorescence experiments. We brought to light that single amino acid replacements in the plastoquinone (PQ)-binding niche of the photosystem II D1 protein impair electron transport (ET) efficiency between quinones and confer increased flexibility to the host membranes, expanding to the entire cells. Hence, a more flexible environment in the PQ-binding niche has been associated to a less efficient ET. A similar function/dynamics relationship was also demonstrated in Rhodobacter sphaeroides reaction centers having inhibited ET, indicating that flexibility at the quinones region plays a crucial role in evolutionarily distant organisms. Instead, a different functional/dynamical correlation was observed in algal mutants hosting a single amino acid replacement residing in a D1 domain far from the PQ-binding niche. Noteworthy, this mutant displayed the highest degree of flexibility, and besides having a nativelike ET efficiency in physiological conditions, it acquired novel, to our knowledge, phenotypic traits enabling it to preserve a high maximal quantum yield of photosystem II photochemistry in extreme habitats. Overall, in the nanosecond timescale, the degree of the observed flexibility is related to the mutation site; in the picosecond timescale, we highlighted the presence of a more pronounced dynamic heterogeneity in all mutants compared to the native cells, which could be related to a marked chemically heterogeneous environment.

Structure and Dynamics of the Hydration Shell: Spatially Decomposed Time Correlation Approach

Molecular simulations provide insight into solvation structures and dynamics with unparalleled spatial and temporal resolution. Here, we take advantage of this fact and develop a set of generally applicable computational tools for a detailed analysis of the hydration shell around an ionic or molecular solute. These tools allow us to quantify and visualize orientationally resolved radial distribution functions as well as distance-resolved orientational time-correlation functions of water molecules surrounding the solute. Such a detailed view of the hydration shells allows us to unravel important structural and dynamical features, which are not accessible when employing standard analysis techniques.

Ion-Specific Confined Water Dynamics in Convex Nanopores of Gemini Surfactant Lyotropic Liquid Crystals

The impact of pore geometry and functionality on the dynamics of water nanoconfined in porous media are the subject of some debate. We report the synthesis and small-angle X-ray scattering (SAXS) characterization of a series of perdeuterated gemini surfactant lyotropic liquid crystals (LLCs), in which convex, water-filled nanopores of well-defined dimensions are lined with carboxylate functionalities. Quasielastic neutron scattering (QENS) measurements of the translational water dynamics in these dicarboxylate LLC nanopores as functions of the surfactant hydration state and the charge compensating counterion (Na⁺, K⁺, NMe₄⁺) reveal that the measured dynamics depend primarily on surfactant hydration, with an unexpected counterion dependence that varies with hydration number. We rationalize these trends in terms of a balance between counterion-water attractions and the nanopore volume excluded by the counterions. On the basis of electron density maps derived from SAXS analyses of these LLCs, we directly show that the volume excluded by the counterions depends on both their size and spatial distribution in the water-filled channels. The translational water dynamics in the convex pores of these LLCs are also slower than those reported in the concave pores of AOT reverse micelles, implying that water dynamics also depend on the nanopore curvature.

Computational discovery of chemically patterned surfaces that effect unique hydration water dynamics

The interactions of water with solid surfaces govern their apparent hydrophobicity/hydrophilicity, influenced at the molecular scale by surface coverage of chemical groups of varied nonpolar/polar character. Recently, it has become clear that the precise patterning of surface groups, and not simply average surface coverage, has a significant impact on the structure and thermodynamics of hydration layer water, and, in turn, on macroscopic interfacial properties. Here we show that patterning also controls the dynamics of hydration water, a behavior frequently thought to be leveraged by biomolecules to influence functional dynamics, but yet to be generalized. To uncover the role of surface heterogeneities, we couple a genetic algorithm to iterative molecular dynamics simulations to design the patterning of surface functional groups, at fixed coverage, to either minimize or maximize proximal water diffusivity. Optimized surface configurations reveal that clustering of hydrophilic groups increases hydration water mobility, while dispersing them decreases it, but only if hydrophilic moieties interact with water through directional, hydrogen-bonding interactions. Remarkably, we find that, across different surfaces, coverages, and patterns, both the chemical potential for inserting a methane-sized hydrophobe near the interface and, in particular, the hydration water orientational entropy serve as strong predictors for hydration water diffusivity, suggesting that these simple thermodynamic quantities encode the way surfaces control water dynamics. These results suggest a deep and intriguing connection between hydration water thermodynamics and dynamics, demonstrating that subnanometer chemical surface patterning is an important design parameter for engineering solid-water interfaces with applications spanning separations to catalysis.

Hydration water dynamics around a protein surface: a first passage time approach

A stochastic noise-driven dynamic model is proposed to study the diffusion of water molecules around a protein surface, under the effect of thermal fluctuations that arise due to the collision of water molecules with the surrounding environment. The underlying dynamics of such a system may be described in the framework of the generalized Langevin equation, where the thermal fluctuations are assumed to be algebraically correlated in time, which governs the non-Markovian behavior of the system. Results of the calculations of mean-square displacement and the velocity autocorrelation function reveal that the hydration water around the protein surface follows subdiffusive dynamics at long times. Analytical expressions for the first passage time distribution, survival probability, mean residence time and mean first passage time of water molecules are derived for different boundary conditions, to analyze hydration water dynamics under the effect of thermally correlated noise. The results depict a unimodal distribution of the first passage time unlike Brownian motion. The survival probability of hydration water follows a stretched exponential decay for both boundary conditions. The mean residence time of the hydration water molecule for different initial positions increases with increase in the complexity/heterogeneity of the surrounding environment for both boundary conditions. The mean first passage time of the water molecule to reach the absorbing/reflecting boundary follows an asymptotic power law with respect to the thickness of the hydration layer, and increases with increase in the complexity/heterogeneity of the environment.

Anomalous water dynamics at surfaces and interfaces: synergistic effects of confinement and surface interactions

In nature, water is often found in contact with surfaces that are extended on the scale of molecule size but small on a macroscopic scale. Examples include lipid bilayers and reverse micelles as well as biomolecules like proteins, DNA and zeolites, to name a few. While the presence of surfaces and interfaces interrupts the continuous hydrogen bond network of liquid water, confinement on a mesoscopic scale introduces new features. Even when extended on a molecular scale, natural and biological surfaces often have features (like charge, hydrophobicity) that vary on the scale of the molecular diameter of water. As a result, many new and exotic features, which are not seen in the bulk, appear in the dynamics of water close to the surface. These different behaviors bear the signature of both water-surface interactions and of confinement. In other words, the altered properties are the result of the synergistic effects of surface-water interactions and confinement. Ultrafast spectroscopy, theoretical modeling and computer simulations together form powerful synergistic approaches towards an understanding of the properties of confined water in such systems as nanocavities, reverse micelles (RMs), water inside and outside biomolecules like proteins and DNA, and also between two hydrophobic walls. We shall review the experimental results and place them in the context of theory and simulations. For water confined within RMs, we discuss the possible interference effects propagating from opposite surfaces. Similar interference is found to give rise to an effective attractive force between two hydrophobic surfaces immersed and kept fixed at a separation of d, with the force showing an exponential dependence on this distance. For protein and DNA hydration, we shall examine a multitude of timescales that arise from frustration effects due to the inherent heterogeneity of these surfaces. We pay particular attention to the role of orientational correlations and modification of the same due to interaction with the surfaces.

Spatially Heterogeneous Surface Water Diffusivity around Structured Protein Surfaces at Equilibrium

Hydration water on the surface of a protein is thought to mediate the thermodynamics of protein-ligand interactions. For hydration water to play a role beyond modulating global protein solubility or stability, the thermodynamic properties of hydration water must reflect on the properties of the heterogeneous protein surface and thus spatially vary over the protein surface. A potent read-out of local variations in thermodynamic properties of hydration water is its equilibrium dynamics spanning picosecond to nanosecond time scales. In this study, we employ Overhauser dynamic nuclear polarization (ODNP) to probe the equilibrium hydration water dynamics at select sites on the surface of Chemotaxis Y (CheY) in dilute solution. ODNP reports on site-specific hydration water dynamics within 5-10 Å of a label tethered to the biomolecular surface on two separate time scales of motion, corresponding to diffusive water (DW) and protein-water coupled motions, referred to as bound water (BW). We find DW dynamics to be highly heterogeneous across the surface of CheY. We identify a significant correlation between DW dynamics and the local hydropathy of the CheY protein surface, empirically determined by molecular dynamics (MD) simulations, and find the more hydrophobic sites to be hydrated with slower diffusing water. Furthermore, we compare the hydration water dynamics on different polypeptides and liposome surfaces and find the DW dynamics on globular proteins to be significantly more heterogeneous than on intrinsically disordered proteins (IDPs), peptides, and liposomes. The heterogeneity in the hydration water dynamics suggests that structured proteins have the capacity to encode information into the surrounding hydration shell.

Water Dynamics in the Hydration Shells of Biomolecules

The structure and function of biomolecules are strongly influenced by their hydration shells. Structural fluctuations and molecular excitations of hydrating water molecules cover a broad range in space and time, from individual water molecules to larger pools and from femtosecond to microsecond time scales. Recent progress in theory and molecular dynamics simulations as well as in ultrafast vibrational spectroscopy has led to new and detailed insight into fluctuations of water structure, elementary water motions, electric fields at hydrated biointerfaces, and processes of vibrational relaxation and energy dissipation. Here, we review recent advances in both theory and experiment, focusing on hydrated DNA, proteins, and phospholipids, and compare dynamics in the hydration shells to bulk water.

On the coupling between the dynamics of protein and water

The solvation entropy of flexible protein regions is higher than that of rigid regions and contributes differently to the overall thermodynamic stability.

Molecular properties of aqueous solutions: a focus on the collective dynamics of hydration water

When a solute is dissolved in water, their mutual interactions determine the molecular properties of the solute on one hand, and the structure and dynamics of the surrounding water particles (the so-called hydration water) on the other. The very existence of soft matter and its peculiar properties are largely due to the wide variety of possible water-solute interactions. In this context, water is not an inert medium but rather an active component, and hydration water plays a crucial role in determining the structure, stability, dynamics, and function of matter. This review focuses on the collective dynamics of hydration water in terms of retardation with respect to the bulk, and of the number of molecules whose dynamics is perturbed. Since water environments are in a dynamic equilibrium, with molecules continuously exchanging from around the solute towards the bulk and vice versa, we examine the ability of different techniques to measure the water dynamics on the basis of the explored time scales and exchange rates. Special emphasis is given to the collective dynamics probed by extended depolarized light scattering and we discuss whether and to what extent the results obtained in aqueous solutions of small molecules can be extrapolated to the case of large biomacromolecules. In fact, recent experiments performed on solutions of increasing complexity clearly indicate that a reductionist approach is not adequate to describe their collective dynamics. We conclude this review by presenting current ideas that are being developed to describe the dynamics of water interacting with macromolecules.

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.

Nuclear quantum effects in water exchange around lithium and fluoride ions

We employ classical and ring polymer molecular dynamics simulations to study the effect of nuclear quantum fluctuations on the structure and the water exchange dynamics of aqueous solutions of lithium and fluoride ions. While we obtain reasonably good agreement with experimental data for solutions of lithium by augmenting the Coulombic interactions between the ion and the water molecules with a standard Lennard-Jones ion-oxygen potential, the same is not true for solutions of fluoride, for which we find that a potential with a softer repulsive wall gives much better agreement. A small degree of destabilization of the first hydration shell is found in quantum simulations of both ions when compared with classical simulations, with the shell becoming less sharply defined and the mean residence time of the water molecules in the shell decreasing. In line with these modest differences, we find that the mechanisms of the exchange processes are unaffected by quantization, so a classical description of these reactions gives qualitatively correct and quantitatively reasonable results. We also find that the quantum effects in solutions of lithium are larger than in solutions of fluoride. This is partly due to the stronger interaction of lithium with water molecules, partly due to the lighter mass of lithium and partly due to competing quantum effects in the hydration of fluoride, which are absent in the hydration of lithium.

Anomalously Rapid Hydration Water Diffusion Dynamics Near DNA Surfaces

The emerging Overhauser effect dynamic nuclear polarization (ODNP) technique measures the translational mobility of water within the vicinity (5-15 Å) of preselected sites. The work presented here expands the capabilities of the ODNP technique and illuminates an important, previously unseen, property of the translational diffusion dynamics of water at the surface of DNA duplexes. We attach nitroxide radicals (i.e., spin labels) to multiple phosphate backbone positions of DNA duplexes, allowing ODNP to measure the hydration dynamics at select positions along the DNA surface. With a novel approach to ODNP analysis, we isolate the contributions of water molecules at these sites that undergo free translational diffusion from water molecules that either loosely bind to or exchange protons with the DNA. The results reveal that a significant population of water in a localized volume adjacent to the DNA surface exhibits fast, bulk-like characteristics and moves unusually rapidly compared to water found in similar probe volumes near protein and membrane surfaces. Control studies show that the observation of these characteristics are upheld even when the DNA duplex is tethered to streptavidin or the mobility of the nitroxides is altered. This implies that, as compared to protein or lipid surfaces, it is an intrinsic feature of the DNA duplex surface that it interacts only weakly with a significant fraction of the surface hydration water network. The displacement of this translationally mobile water is energetically less costly than that of more strongly bound water by up to several kBT and thus can lower the activation barrier for interactions involving the DNA surface.

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.

Water dynamics in protein hydration shells: the molecular origins of the dynamical perturbation

Protein hydration shell dynamics play an important role in biochemical processes including protein folding, enzyme function, and molecular recognition. We present here a comparison of the reorientation dynamics of individual water molecules within the hydration shell of a series of globular proteins: acetylcholinesterase, subtilisin Carlsberg, lysozyme, and ubiquitin. Molecular dynamics simulations and analytical models are used to access site-resolved information on hydration shell dynamics and to elucidate the molecular origins of the dynamical perturbation of hydration shell water relative to bulk water. We show that all four proteins have very similar hydration shell dynamics, despite their wide range of sizes and functions, and differing secondary structures. We demonstrate that this arises from the similar local surface topology and surface chemical composition of the four proteins, and that such local factors alone are sufficient to rationalize the hydration shell dynamics. We propose that these conclusions can be generalized to a wide range of globular proteins. We also show that protein conformational fluctuations induce a dynamical heterogeneity within the hydration layer. We finally address the effect of confinement on hydration shell dynamics via a site-resolved analysis and connect our results to experiments via the calculation of two-dimensional infrared spectra.

Nontrivial behavior of water in the vicinity and inside lipid bilayers as probed by molecular dynamics simulations

The atomic-scale diffusion of water in the presence of several lipid bilayers mimicking biomembranes is characterized via unconstrained molecular dynamics (MD) simulations. Although the overall water dynamics corresponds well to literature data, namely, the efficient braking near polar head groups of lipids, a number of interesting and biologically relevant details observed in this work have not been sufficiently discussed so far; for instance, the fact that waters "sense" the membrane unexpectedly early, before water density begins to decrease. In this "transitional zone" the velocity distributions of water and their H-bonding patterns deviate from those in the bulk solution. The boundaries of this zone are well preserved even despite the local (<1 nm size) perturbation of the lipid bilayer, thus indicating a decoupling of the surface and bulk dynamics of water. This is in excellent agreement with recent experimental data. Near the membrane surface, water movement becomes anisotropic, that is, solvent molecules preferentially move outward the bilayer. Deep in the membrane interior, the velocities can even exceed those in the bulk solvent and undergo large-scale fluctuations. The analysis of MD trajectories of individual waters in the middle part of the acyl chain region of lipids reveals a number of interesting rare phenomena, such as the fast (ca. 50 ps) breakthrough across the membrane or long-time (up to 750 ps) "roaming" between lipid leaflets. The analysis of these events was accomplished to delineate the mechanisms of spontaneous water permeation inside the hydrophobic membrane core. It was shown that such nontrivial dynamics of water in an "alien" environment is driven by the dynamic heterogeneities of the local bilayer structure and the formation of transient atomic-scale "defects" in it. The picture observed in lipid bilayers is drastically different from that in a primitive membrane mimic, a hydrated cyclohexane slab. The possible biological impact of such phenomena in equilibrated lipid bilayers is discussed.

Theoretical study of HOCl-catalyzed keto-enol tautomerization of β-cyclopentanedione in an explicit water environment

The mechanism of acid-catalyzed keto-enol tautomerization of β-cyclopentanedione (CPD) in solution is studied computationally. Reaction profiles are first calculated for a limited solvation environment using ab initio and density functional methods. Barrier heights for systems including up to one hydration shell of explicit water molecules depend strongly on the number of waters involved in proton transfer and to a lesser but significant extent on the number of waters forming hydrogen bonds with waters in the proton-transfer chain (each such water reduces the barrier by 4.4 kcal/mol on average). Barriers of 8-13 kcal/mol were obtained when a full or nearly full hydration shell was present, consistent with calculations for nonacid-catalyzed keto-enol tautomerization of related molecules. The presence of HOCl reduced the barrier by 4.5 kcal/mol in relation to the gas phase, consistent with the well-known principle that keto-enol tautomerization can be acid- or base-catalyzed. The reaction was also modeled beginning with snapshots of reactant conformations taken from a 300 K molecular dynamics simulation of CPD, HOCl, and 324 explicit waters. Reaction profiles were calculated at a QM/MM level with waters in the first hydration shell either fixed or energy-minimized at each step along the reaction coordinate. A substantial variation in barrier height was observed in both cases, depending primarily on electrostatic interactions (hydrogen bonding) with first-hydration-shell waters and, to a lesser extent, on electrostatic interactions with more distant waters and geometric distortion effects. For the lowest barriers, the extent of barrier reduction by waters involved in proton transfer is consistent with the limited solvation results, but further barrier reduction due to hydrogen bonding to waters involved in proton transfer is not observed. It is postulated that this is because highly flexible structures such as extensive hydrogen bonding networks optimal for reaction are entropically disfavored and so may not contribute significantly to the observed reaction rate.

Nonlinear scaling of surface water diffusion with bulk water viscosity of crowded solutions

The translational hydration dynamics within 0.5-1.5 nm of the surface of a DPPC liposome, a model biomacromolecular surface, is analyzed by the recently developed Overhauser dynamic nuclear polarization (ODNP) technique. We find that dramatic changes to the bulk solvent cause only weak changes in the surface hydration dynamics. Specifically, both a >10-fold increase in bulk viscosity and the restriction of diffusion by confinement on a multiple nm length-scale change the local translational diffusion coefficient of the surface water surrounding the lipid bilayer by <2.5-fold. By contrast, previous ODNP studies have shown that changes to the biomacromolecular surface induced by folding, binding, or aggregation can cause local hydration dynamics to vary by factors of up to 30. We suggest that the surface topology and chemistry at the ≤1.5 nm scale, rather than the characteristics of the solvent, nearly exclusively determine the macromolecule's surface hydration dynamics.

Comparison of the structural stability and dynamic properties of recombinant anthrax protective antigen and its 2-fluorohistidine-labeled analogue

Protective antigen (PA) is the primary protein antigenic component of both the currently used anthrax vaccine and related recombinant vaccines under development. An analogue of recombinant PA (2-FHis rPA) has been recently shown to block the key steps of pore formation in the process of inducing cytotoxicity in cells, and thus can potentially be used as an antitoxin or a vaccine. This rPA analogue was produced by fermentation to incorporate the unnatural amino acid 2-fluorohistidine (2-FHis). In this study, the effects of 2-FHis labeling on rPA antigen's conformational stability and dynamic properties were investigated by various biophysical techniques. Temperature/pH stability profiles of rPA and 2-FHis rPA were analyzed by the empirical phase diagram (EPD) approach, and physical stability differences between them were identified. Results showed that rPA and 2-FHis rPA had similar stability at pH 7-8. With decreasing solution pH, however, 2-FHis rPA was found to be more stable. Dynamic sensitive measurements of the two proteins at pH 5 found that 2-FHis rPA was more dynamic and/or differentially hydrated under acidic pH conditions. The biophysical characterization and stability data provide information useful for the potential development of 2-FHis rPA as a more stable rPA vaccine candidate.

Neutron spin echo measurements of monolayer and capillary condensed water in MCM-41 at low temperatures

Neutron spin echo measurements of monolayer and capillary condensed heavy water (D(2)O) confined in MCM-41 C10 (pore diameter 2.10 nm) were performed in a temperature range of 190-298 K. The intermediate scattering functions were analyzed by the Kohlrausch-Williams-Watts stretched exponential function. The relaxation times of confined D(2)O in the capillary condensed state follow remarkably well the Vogel-Fulcher-Tammann equation between 298 and 220 K, whereas below 220 K they show an Arrhenius type behavior. That is, the fragile-to-strong (FTS) dynamic crossover occurs, which has never been seen in experiments on bulk water. On the other hand, for monolayer D(2)O, the FTS dynamic crossover was not observed in the temperature range measured. The FTS dynamic crossover observed in capillary condensed water would take place in the central region of the pore, not near the pore surface. Because the tetrahedral-like water structure in the central region of the pore is more preserved than that near the pore surface, the FTS dynamic crossover would be concerned with the tetrahedral-like water structure.