Molecular Dynamics of Water at the Protein−Solvent Interface

Dimensionality of diffusive exploration at the protein interface in solution

The dynamics of water are critically important to the energies of interaction between proteins and substrates and determine the efficiency of transport at the interface. The magnetic field dependence of the nuclear spin-lattice relaxation rate constant 1/T(1) of water protons provides a direct characterization of water diffusional dynamics at the protein interface. We find that the surface-average translational correlation time is 30-40 ps and the magnetic field dependence of the water proton 1/T(1) is characteristic of two-dimensional diffusion of water in the protein interfacial region. The reduced dimensionality substantially increases the intermolecular re-encounter probability and the efficiency of the surface exploration by the small molecule, water in this case. We propose a comprehensive theory of the translational effects of a small diffusing particle confined in the vicinity of a spherical macromolecule as a function of the relative size of the two particles. We show that the change in the apparent dimensionality of the diffusive exploration is a general result of the small diffusing particle encountering a much larger particle that presents a diffusion barrier. Examination of the effects of the size of the confinement relative to the macromolecule size reveals that the reduced dimensionality characterizing the small-molecule diffusion persists to remarkably small radius ratios. The experimental results on several different proteins in solution support the proposed theoretical model, which may be generalized to other small-particle-large-body systems like vesicles and micelles.

Site-dependent excited-state dynamics of a fluorescent probe bound to avidin and streptavidin

The excited-state dynamics of biotin-spacer-Lucifer-Yellow (LY) constructs bound to avidin (Avi) and streptavidin (Sav) was investigated using femtosecond spectroscopy. Two different locations in the proteins, identified by molecular dynamics simulations of Sav, namely the entrance of the binding pocket and the protein surface, were probed by varying the length of the spacer. A reduction of the excited-state lifetime, stronger in Sav than in Avi, was observed with the long spacer construct. Transient absorption measurements show that this effect originates from an electron transfer quenching of LY, most probably by a nearby tryptophan residue. The local environment of the LY chromophore could be probed by measuring the time-dependent polarisation anisotropy and Stokes shift of the fluorescence. Substantial differences in both dynamics were observed. The fluorescence anisotropy decays analysed by using the wobbling-in-a-cone model reveal a much more constrained environment of the chromophore with the short spacer. Moreover, the dynamic Stokes shift is multiphasic in all cases, with a approximately 1 ps component that can be ascribed to diffusive motion of bulk-like water molecules, and with slower components with time constants varying not only with the spacer, but with the protein as well. These slow components, which depend strongly on the local environment of the probe, are ascribed to the motion of the hydration layer coupled to the conformational dynamics of the protein.

Viscoelastic subdiffusion: from anomalous to normal

We study viscoelastic subdiffusion in bistable and periodic potentials within the generalized Langevin equation approach. Our results justify the (ultra)slow fluctuating rate view of the corresponding bistable non-Markovian dynamics which displays bursting and anticorrelation of the residence times in two potential wells. The transition kinetics is asymptotically stretched exponential when the potential barrier V0 several times exceeds thermal energy k(B)T [V(0) approximately (2-10)k(B)T] and it cannot be described by the non-Markovian rate theory (NMRT). The well-known NMRT result approximates, however, ever better with the increasing barrier height, the most probable logarithm of the residence times. Moreover, the rate description is gradually restored when the barrier height exceeds a fuzzy borderline which depends on the power-law exponent of free subdiffusion alpha . Such a potential-free subdiffusion is ergodic. Surprisingly, in periodic potentials it is not sensitive to the barrier height in the long time asymptotic limit. However, the transient to this asymptotic regime is extremally slow and it does profoundly depend on the barrier height. The time scale of such subdiffusion can exceed the mean residence time in a potential well or in a finite spatial domain by many orders of magnitude. All these features are in sharp contrast with an alternative subdiffusion mechanism involving jumps among traps with the divergent mean residence time in these traps.

The effect of nanometre-scale structure on interfacial energy

Natural surfaces are often structured with nanometre-scale domains, yet a framework providing a quantitative understanding of how nanostructure affects interfacial energy, gamma(SL), is lacking. Conventional continuum thermodynamics treats gamma(SL) solely as a function of average composition, ignoring structure. Here we show that, when a surface has domains commensurate in size with solvent molecules, gamma(SL) is determined not only by its average composition but also by a structural component that causes gamma(SL) to deviate from the continuum prediction by a substantial amount, as much as 20% in our system. By contrasting surfaces coated with either molecular- (<2 nm) or larger-scale domains (>5 nm), we find that whereas the latter surfaces have the expected linear dependence of gamma(SL) on surface composition, the former show a markedly different non-monotonic trend. Molecular dynamics simulations show how the organization of the solvent molecules at the interface is controlled by the nanostructured surface, which in turn appreciably modifies gamma(SL).

Understanding dynamics of myoglobin in heterogeneous aqueous environments using coupled water fractions

This work presents an analysis of near environment of myoglobin (Mb) in different aqueous solutions (in the presence of NaCl, sucrose, trehalose, urea, and glycerol) using the coupled water fractions measured using a quartz crystal microbalance (QCM). The secondary structural features of the protein from circular dichroic (CD) spectroscopy and the coupled water fractions give important clues to the overall dynamics of the protein. Using time resolved fluorescence, these leads have been applied to understand the observed lifetime relaxations of Mb. Though the time scales of observation of coupled water and the lifetimes are very different, our study suggests that the trends in coupled water fraction seem to be good indicators for regulation of the relaxation dynamics of the protein. The relaxations generally show a triphasic distribution of time scales. The initial relaxation in the picoseconds time scale represents the local motions of coupled water followed by a slightly slower decay in hundreds of picoseconds attributable to coupled water-'quasi free' water interactions. The third nanosecond lifetime is due to changes in transitions in isomers of hydrated protein. The dynamics of coupled water in Mb with NaCl is the fastest (around 21 ps) and is slowest in glycerol (250 ps). The results strongly indicate that it is the resident times of water molecules that play a dominant role in the overall stability of protein in a particular hydrated isomer and not just always the number of such water molecules in the hydrated protein.

Transport properties of water at functionalized molecular interfaces

Understanding transport properties of solvent such as diffusion and viscosity at interfaces with biomacromolecules and hard materials is of fundamental importance to both biology and biotechnology. Our study utilizes equilibrium molecular dynamics simulations to calculate solvent transport properties at a model peptide and microarray surface. Both diffusion and selected components of viscosity are considered. Solvent diffusion is found to be affected near the peptide and surface. The stress-stress correlation function of solvent near the hard surface exhibits long time memory. Both diffusion and viscosity are shown to be closely correlated with the density distribution function of water along the microarray surface.

Relaxation dynamics of nucleosomal DNA

Recent experimental and theoretical evidence demonstrates that proteins and water in the hydration layer can follow complex stretched exponential or power law relaxation dynamics. Here, we report on a 50 ns all atom molecular dynamics (MD) simulation of the yeast nucleosome, where the interactions between DNA, histones, surrounding water and ions are explicitly included. DNA interacts with the histone core in 14 locations, approximately every 10.4 base pairs. We demonstrate that all sites of interaction exhibit anomalously slow power law relaxation, extending up to 10 ns, while fast exponential relaxation dynamics of hundreds of picoseconds applies to DNA regions outside these locations. The appearance of 1/f(alpha) noise or pink noise in DNA dynamics is ubiquitous. For histone-bound nucleotide dynamics alpha --> 1 and is a signature of complexity of the protein-DNA interactions. For control purposes two additional DNA simulations free of protein are conducted. Both utilize the same sequence of DNA, as found the in the nucleosome. In one simulation the initial conformation of the double helix is a straight B-form. In the other, the initial conformation is super helical. Neither of these simulations exhibits the variation of alpha as a function of position, the measure of power law for dynamical behavior, which we observe in the nucleosome simulation. The unique correspondence (high alpha to DNA-histone interaction sites, low alpha to free DNA sites), suggests that alpha may be an important and new quantification of protein-DNA interactions for future experiments.

Detection of nanosecond time scale side-chain jumps in a protein dissolved in water/glycerol solvent

In solution, the correlation time of the overall protein tumbling, tau(R), plays a role of a natural dynamics cutoff-internal motions with correlation times on the order of tau ( R ) or longer cannot be reliably identified on the basis of spin relaxation data. It has been proposed some time ago that the 'observation window' of solution experiments can be expanded by changing the viscosity of solvent to raise the value of tau(R). To further explore this concept, we prepared a series of samples of alpha-spectrin SH3 domain in solvent with increasing concentration of glycerol. In addition to the conventional (15)N labeling, the protein was labeled in the Val, Leu methyl positions ((13)CHD(2) on a deuterated background). The collected relaxation data were used in asymmetric fashion: backbone (15)N relaxation rates were used to determine tau(R) across the series of samples, while methyl (13)C data were used to probe local dynamics (side-chain motions). In interpreting the results, it has been initially suggested that addition of glycerol leads only to increases in tau(R), whereas local motional parameters remain unchanged. Thus the data from multiple samples can be analyzed jointly, with tau(R) playing the role of experimentally controlled variable. Based on this concept, the extended model-free model was constructed with the intent to capture the effect of ns time-scale rotameric jumps in valine and leucine side chains. Using this model, we made a positive identification of nanosecond dynamics in Val-23 where ns motions were already observed earlier. In several other cases, however, only tentative identification was possible. The lack of definitive results was due to the approximate character of the model-contrary to what has been assumed, addition of glycerol led to a gradual 'stiffening' of the protein. This and other observations also shed light on the interaction of the protein with glycerol, which is one of the naturally occurring osmoprotectants. In particular, it has been found that the overall protein tumbling is controlled by the bulk solvent, and not by a thin solvation layer which contains a higher proportion of water.

Interfacial water at protein surfaces: wide-line NMR and DSC characterization of hydration in ubiquitin solutions

Wide-line 1H-NMR and differential scanning calorimetry measurements were done in aqueous solutions and on lyophilized samples of human ubiquitin between -70 degrees C and +45 degrees 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 x g(-1) x K(-1), and the magnitudes are higher than that of pure/bulk water (4.2 J x g(-1) x K(-1)). 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 x g(-1) x K(-1). It could be concluded that the protein ubiquitin and its hydrate layer behave as a highly interconnected single phase in a thermodynamic sense.

Structural analysis and the effect of cyclo(His-Pro) dipeptide on neurotoxins--a dynamics and density functional theory study

The switching propensity and maximum probability of occurrence of the side chain imidazole group in the dipeptide cyclo(His-Pro) (CHP) were studied by applying molecular dynamics simulations and density functional theory. The atomistic behaviour of CHP with the neurotoxins glutamate (E) and paraquat (Pq) were also explored; E and Pq engage in hydrogen bond formation with the diketopiperazine (DKP) ring of the dipeptide, with which E shows a profound interaction, as confirmed further by NH and CO stretching vibrational frequencies. The effect of CHP was found to be greater on E than on Pq neurotoxin. A ring puckering study indicated a twist boat conformation for the six-membered DKP ring. Molecular electrostatic potential (MESP) mapping was also used to explore the hydrogen bond interactions prevailing between the neurotoxins and the DKP ring. The results of this study reveal that the DKP ring of the dipeptide CHP can be expected to play a significant role in reducing effects such as oxidative stress and cell death caused by neurotoxins.

Hydration water dynamics near biological interfaces

We performed classical molecular dynamics simulations using both fixed-charge and polarizable water and protein force fields to contrast the hydration dynamics near hydrophilic and amphiphilic peptides as a function of temperature. The high peptide concentrations we use serve as a model for the surface of folded proteins where hydration layers around each residue overlap significantly. Through simulation we determine that there are notable differences in the water dynamics analyzed from the outer and inner hydration layer regions of the amphiphilic peptide solution that explains the experimentally observed presence of two translational relaxations, while the hydrophilic peptide solution shows only a single non-Arrhenius translational process with no distinction between hydration layers. Given that water dynamics for the amphiphilic peptide system reproduces all known rotational and translational hydration dynamical anomalies exhibited by hydration water near protein surfaces, our analysis provides strong evidence that dynamical signatures near biological interfaces arises because of frustration in the hydration dynamics induced by chemical heterogeneity, as opposed to just topological roughness, of the protein surface.

Fractional Fokker-Planck subdiffusion in alternating force fields

The fractional Fokker-Planck equation for subdiffusion in time-dependent force fields is derived from the underlying continuous time random walk. Its limitations are discussed and it is then applied to the study of subdiffusion under the influence of a time-periodic rectangular force. As a main result, we show that such a force does not affect the universal scaling relation between the anomalous current and diffusion when applied to the biased dynamics: in the long-time limit, subdiffusion current and anomalous diffusion are immune to the driving. This is in sharp contrast with the unbiased case when the subdiffusion coefficient can be strongly enhanced, i.e., a zero-frequency response to a periodic driving is present.

Advantage of being a dimer for Serratia marcescens endonuclease

The monomer and dimer of the bacterium Serratia marcescens endonuclease (SMnase) are each catalytically active, and the two subunits of the dimer function independently of each other. Nature, however, chooses the dimer form instead of the monomer. In order to explain this, we performed molecular dynamics (MD) simulations of both model-built complexes of a subunit of SMnase and the dimer with DNA in aqueous solution. We estimated the electrostatic binding energy, analyzed the distribution and dynamics of water around the complexes, identified water clusters in the protein, and related the dynamics of water to the protein's function. We find that the dimer form has an electrostatic advantage over the monomer to associate with DNA. Although Mg(2+) remains hexa-coordinated during the simulation, the binding pathway of DNA to Mg(2+) changes from inner-sphere binding in the monomer to outer-sphere in the dimer, which may be more energetically favorable. In addition, two water clusters in the active site of each monomer and in the dimer complex were identified and localized in two regions, named the "stabilizing" and "working" regions. Water in the "working" region in the dimer complex has larger fluctuations than that in the monomer.

Dinâmica molecular: teoria e aplicações em planejamento de fármacos

Dinâmica Molecular (DM) é uma ferramenta computacional poderosa usada em Química Medicinal para o planejamento racional de fármacos. DM é uma extensão da Mecânica Molecular, onde o comportamento dinâmico de um sistema molecular é simulado através da integração numérica das equações de movimento. Esta técnica tem sido usada extensivamente para auxiliar e complementar o planejamento de novos ligantes de um alvo terapêutico, bem como estimar a sua potência. Este artigo enfoca a teoria básica da DM clássica e suas importantes aplicações no planejamento racional de potenciais compostos bioativos, particularmente compostos com atividade anti-HIV.

Exploring the conserved water site and hydration of a coiled-coil trimerisation motif: a MD simulation study

The solvent structure and dynamics around ccbeta-p, a 17-residue peptide that forms a parallel three-stranded alpha-helical coiled coil in solution, was analysed through 10 ns explicit solvent molecular dynamics (MD) simulations at 278 and 330 K. Comparison with two corresponding simulations of the monomeric form of ccbeta-p was used to investigate the changes of hydration upon coiled-coil formation. Pronounced peaks in the solvent density distribution between residues Arg8 and Glu13 of neighbouring helices show the presence of water bridges between the helices of the ccbeta-p trimer; this is in agreement with the water sites observed in X-ray crystallography experiments. Interestingly, this water site is structurally conserved in many three-stranded coiled coils and, together with the Arg and Glu residues, forms part of a motif that determines three-stranded coiled-coil formation. Our findings show that little direct correlation exists between the solvent density distribution and the temporal ordering of water around the trimeric coiled coil. The MD-calculated effective residence times of up to 40 ps show rapid exchange of surface water molecules with the bulk phase, and indicate that the solvent distribution around biomolecules requires interpretation in terms of continuous density distributions rather than in terms of discrete molecules of water. Together, our study contributes to understanding the principles of three-stranded coiled-coil formation.

Comment on 'Diffusion of water and sodium counter-ions in nanopores of a β-lactoglobulin crystal: a molecular dynamics study'

Malek et al (2005 Nanotechnology 16 S522) studied the diffusion of water and sodium counter-ions in a β-lactoglobulin (βLG) crystal from a 5 ns molecular dynamics simulation. They presented water density in the main pore, and water mobility in the core and surface zones. Unfortunately, their results were incorrect and misleading. In this comment, we provide the results from our analysis for the same system.

Reply to Comment on 'Diffusion of water and sodium counter-ions in nanopores of β-lactoglobulin crystal: a molecular simulation study'

The analysis in Hu and Jiang's Comment to our paper cannot reveal long-time diffusion, and incorrectly led the authors to conclude that the diffusion in beta-lactoglobuline is anomalous. In this context, the limitations of applying a mean-square displacement analysis to short, heterogeneous pore channels are discussed. A more appropriate approach based on first-passage time analysis is illustrated by a detailed analysis of water motion in a natural membrane protein channel. The partitioning and the motion of water molecules between core and surface hydration layers is discussed. Finally, the calculation of the water density profile is commented upon.

Isn't self-assembly a misnomer? Multi-disciplinary arguments in favor of co-assembly

Self-assembly is a ubiquitous physicochemical phenomenon. It is inherent to molecular recognition effects in the biological domain, and thus presents a basis for understanding the constitution and dynamic organization of living beings. However, it is argued in this paper that the very notion of self-assembly presents an incorrect and misleading one for the physical effects that it stands for. Self-assembly presents a misnomer as each self-assembly process not only depends on, but also includes a certain degree of ordering in the immediate surrounding of the 'self-assembling' entities. Claims are provided to support the idea that every self-assembly process is, in fact, a co-assembly event, as it implies mutual changes and structuring of both the 'self-assembled' system and its nearest surrounding. This point of view is related to the systemic observation that there could be no physical qualities independent on the context of their existence. In order to illustrate that the concept of co-assembly presents only an instance of a more profound and systemic philosophy of life, the examples of the interplay between internal and external factors of development are provided from the fields of quantum and classical physics, biology, evolutionary theory and psychology.

Hierarchical organization of eglin c native state dynamics is shaped by competing direct and water-mediated interactions

The native state dynamics of the small globular serine protease inhibitor eglin c has been studied in a long 336 ns computer simulation in explicit solvent. We have elucidated the energy landscape explored during the course of the simulation by using Principal Component Analysis. We observe several basins in the energy landscape in which the system lingers for extended periods. Through an iterative process we have generated a tree-like hierarchy of states describing the observed dynamics. We observe a range of divergent contact types including salt bridges, hydrogen bonds, hydrophilic interactions, and hydrophobic interactions, pointing to the frustration between competing interactions. Additionally, we find evidence of competing water-mediated interactions. Divergence in water-mediated interactions may be found to supplement existing direct contacts, but they are also found to be independent of such changes. Water-mediated contacts facilitate interactions between residues of like charge as observed in the simulation. Our results provide insight into the complexity of the dynamic native state of a globular protein and directly probe the residual frustration in the native state.

Dynamics at the protein-water interface from 17O spin relaxation in deeply supercooled solutions

Most of the decisive molecular events in biology take place at the protein-water interface. The dynamical properties of the hydration layer are therefore of fundamental importance. To characterize the dynamical heterogeneity and rotational activation energy in the hydration layer, we measured the (17)O spin relaxation rate in dilute solutions of three proteins in a wide temperature range extending down to 238 K. We find that the rotational correlation time can be described by a power-law distribution with exponent 2.1-2.3. Except for a small fraction of secluded hydration sites, the dynamic perturbation in the hydration layer is the same for all proteins and does not differ in any essential way from the hydration shell of small organic solutes. In both cases, the dynamic perturbation factor is <2 at room temperature and exhibits a maximum near 262 K. This maximum implies that, at low temperatures, the rate of water molecule rotation has a weaker temperature dependence in the hydration layer than in bulk water. We attribute this difference to the temperature-independent constraints that the protein surface imposes on the water H-bond network. The free hydration layer studied here differs qualitatively from confined water in solid protein powder samples.

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

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

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

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

Long-range influence of carbohydrates on the solvation dynamics of water--answers from terahertz absorption measurements and molecular modeling simulations

We present new terahertz (THz) spectroscopic measurements of solvated sugars and compare the effect of two disaccharides (trehalose and lactose) and one monosaccharide (glucose) with respect to the solute-induced changes in the sub-picosecond network dynamics of the hydration water. We found that the solute affects the fast collective network motions of the solvent, even beyond the first solvation layer. For all three carbohydrates, we find an increase of 2-4% in the THz absorption coefficient of the hydration water in comparison to bulk water. Concentration-dependent changes in the THz absorption between 2.1 and 2.8 THz of the solute-water mixture were measured with a precision better than 1% and were used to deduce a dynamical hydration shell, which extends from the surface up to 5.7 +/- 0.4 and 6.5 +/- 0.9 A for the disaccharides lactose and trehalose, respectively, and 3.7 +/- 0.9 A for the glucose. This exceeds the values for the static hydration shell as determined, for example, by scattering, where the long-range structure was found to be not significantly affected by the solute beyond the first hydration shell. When comparing all three carbohydrates, we found that the solute-induced change in the THz absorption depends on the product of molar concentration of the solute and the number of hydrogen bonds between the carbohydrate and water molecules. We can conclude that the long-range influence on the sub-picosecond collective water network motions of the hydration water is directly correlated with the average number of hydrogen bonds between the molecule and adjacent water molecules for carbohydrates. This implies that monosaccharides have a smaller influence on the surrounding water molecules than disaccharides. This could explain the bioprotection mechanism of sugar-water mixtures, which has been found to be more effective for disaccharides than for monosaccharides.

Water/hydrocarbon interfaces: effect of hydrocarbon branching on single-molecule relaxation

Water/hydrocarbon interfaces are studied using molecular dynamics simulations in order to understand the effect of hydrocarbon branching on the dynamics of the system at and away from the interface. A recently proposed procedure for studying the intrinsic structure of the interface in such systems is utilized, and dynamics are probed in the usual laboratory frame as well as the intrinsic frame. The use of these two frames of reference leads to insight into the effect of capillary waves at the interface on dynamics. The systems were partitioned into zones with a width of 5 A, and a number of quantities of dynamical relevance, namely, the residence times, mean squared displacements, the velocity auto correlation functions, and orientational time correlations for molecules of both phases, were calculated in the laboratory and intrinsic frames at and away from the interface. For the aqueous phase, translational motion is found to be (a) diffusive at long times and not anomalous as in proteins or micelles, (b) faster at the interface than in the bulk, and (c) faster upon reduction of the effect of capillary waves. The rotational motion of water is (a) more anisotropic at the interface than in the bulk and (b) dependent on the orientation of the covalent O-H bond with respect to the plane of the interface. The effect of hydrocarbon branching on aqueous dynamics was found to be small, a result similar to the effect on the interfacial water structure. The hydrocarbon phase shows a larger variation for all dynamical probes, a trend consistent with their interfacial structure.

Hydration effect on low-frequency protein dynamics observed in simulated neutron scattering spectra

Hydration effects on protein dynamics were investigated by comparing the frequency dependence of the calculated neutron scattering spectra between full and minimal hydration states at temperatures between 100 and 300 K. The protein boson peak is observed in the frequency range 1–4 meV at 100 K in both states. The peak frequency in the minimal hydration state shifts to lower than that in the full hydration state. Protein motions with a frequency higher than 4 meV were shown to undergo almost harmonic motion in both states at all temperatures simulated, whereas those with a frequency lower than 1 meV dominate the total fluctuations above 220 K and contribute to the origin of the glass-like transition. At 300 K, the boson peak becomes buried in the quasielastic contributions in the full hydration state but is still observed in the minimal hydration state. The boson peak is observed when protein dynamics are trapped within a local minimum of its energy surface. Protein motions, which contribute to the boson peak, are distributed throughout the whole protein. The fine structure of the dynamics structure factor is expected to be detected by the experiment if a high resolution instrument (<∼20 μeV) is developed in the near future.

Anomalous escape governed by thermal 1/f noise

We present an analytic study for subdiffusive escape of overdamped particles out of a cusp-shaped parabolic potential well which are driven by thermal, fractional Gaussian noise with a 1/omega 1-alpha power spectrum. This long-standing challenge becomes mathematically tractable by use of a generalized Langevin dynamics via its corresponding non-Markovian, time-convolutionless master equation: We find that the escape is governed asymptotically by a power-law whose exponent depends exponentially on the ratio of barrier height and temperature. This result is in distinct contrast to a description with a corresponding subdiffusive fractional Fokker-Planck approach, thus providing experimentalists an amenable testbed to differentiate between the two escape scenarios.

Mapping hydration dynamics around a protein surface

Protein surface hydration is fundamental to its structure and activity. We report here the direct mapping of global hydration dynamics around a protein in its native and molten globular states, using a tryptophan scan by site-specific mutations. With 16 tryptophan mutants and in 29 different positions and states, we observed two robust, distinct water dynamics in the hydration layer on a few ( approximately 1-8 ps) and tens to hundreds of picoseconds ( approximately 20-200 ps), representing the initial local relaxation and subsequent collective network restructuring, respectively. Both time scales are strongly correlated with protein's structural and chemical properties. These results reveal the intimate relationship between hydration dynamics and protein fluctuations and such biologically relevant water-protein interactions fluctuate on picosecond time scales.