Glasslike dynamical behavior of the plastocyanin hydration water

Mode coupling behavior and fragile to strong transition of trehalose in a binary mixture with water upon supercooling

We perform molecular dynamics simulations of a binary mixture of water and trehalose with the TIP4P/Ice water model. We analyze the slow dynamics of trehalose molecules in the mildly supercooled region for concentrations of 3.66 and 18.57 wt. %. We previously studied the dynamics of water in the same mixtures. Supercooled TIP4P/Ice water solvating trehalose molecules was found to follow the Mode Coupling Theory (MCT) and to undergo a transition from a fragile to a strong behavior for both concentrations. Here, we show that also the dynamics of trehalose molecules follows the MCT and displays a fragile to strong crossover (FSC). The results show that trehalose in binary mixtures with water shares with it the dynamical behavior typical of glass forming liquids. Moreover, the FSC for trehalose structural relaxation times is found to occur at temperatures close to those previously obtained for water in the same solutions, showing that the dynamics of the solute is strongly coupled to that of the solvent. We also perform a MCT test showing that the trehalose dynamics obeys the MCT time-temperature superposition principle and that the exponents derived from the theory and the ones obtained from fitting procedure of the relaxation times are comparable, confirming that trehalose molecules in supercooled water solutions follow the MCT of glassy dynamics. Moreover, as predicted by the theory, trehalose particles have MCT parameters comparable to those of water in the same mixtures. This is an important result, given that MCT was originally formulated for monoatomic particles.

Structure and slow dynamics of protein hydration water with cryopreserving DMSO and trehalose upon cooling

We study, through molecular dynamics simulations, three aqueous solutions with one lysozyme protein and three different concentrations of trehalose and dimethyl sulfoxide (DMSO). We analyze the structural and dynamical properties of the protein hydration water upon cooling. We find that trehalose plays a major role in modifying the structure of the network of HBs between water molecules in the hydration layer of the protein. The dynamics of hydration water presents, in addition to the α-relaxation, typical of glass formers, a slower long-time relaxation process, which greatly slows down the dynamics of water, particularly in the systems with trehalose, where it becomes dominant at low temperatures. In all the solutions, we observe, from the behavior of the α-relaxation times, a shift of the Mode Coupling Theory crossover temperature and the fragile-to-strong crossover temperature toward higher values with respect to bulk water. We also observe a strong-to-strong crossover from the temperature behavior of the long-relaxation times. In the aqueous solution with only DMSO, the transition shifts to a lower temperature than in the case with only lysozyme reported in the literature. We observe that the addition of trehalose to the mixture has the opposite effect of restoring the original location of the strong-to-strong crossover. In all the solutions analyzed in this work, the observed temperature of the protein dynamical transition is slightly shifted at lower temperatures than that of the strong-to-strong crossover, but their relative order is the same, showing a correlation between the motion of the protein and that of the hydration water.

Time-Resolved Fluorescence and Essential Dynamics Study on the Structural Heterogeneity of p53DBD Bound to the Anticancer p28 Peptide

Time-resolved fluorescence emission was combined with molecular dynamics (MD) simulations to investigate the DNA-binding domain (DBD) of the tumor suppressor p53 alone and its complex with the anticancer peptide p28 (DBD/p28). The fluorescence emission decay of the lone Trp residue, from both DBD and DBD/p28, was well-described by a stretched exponential function. Such a behavior was ascribed to heterogeneity in the Trp relaxation behavior, likely due to the coexistence of different conformational states. The increase of the stretching parameter, on passing from DBD to DBD/p28, indicates a reduced heterogeneity in the Trp146 environment for DBD/p28. Moreover, the effects of p28 on the global dynamics of DBD were analyzed by the essential dynamics method on 30 ns long MD trajectories of both DBD and DBD/p28. We found the establishment of wide-amplitude anharmonic modes throughout the DBD molecule, with a particularly high amplitude being detected in the DNA-binding region. These modes are significantly reduced when DBD is bound to p28, consistently with a structure stabilization. In summary, the results indicate that p28 binding has a strong effect on both the local and global heterogeneity of DBD, thus providing some hints to the understanding of its anticancer activity.

Characterization of hydration water in supercooled water-trehalose solutions: The role of the hydrogen bonds network

The structural and dynamical properties of hydration water in aqueous solutions of trehalose are studied with molecular dynamics simulation. We simulate the systems in the supercooled region to investigate how the interaction with the trehalose molecules modifies the hydrogen bond network, the structural relaxation, and the diffusion properties of hydration water. The analysis is performed by considering the radial distribution functions, the residence time of water molecules in the hydration shell, the two body excess entropy, and the hydrogen bond water-water and water-trehalose correlations of the hydration water. The study of the two body excess entropy shows the presence of a fragile to strong crossover in supercooled hydration water also found in the relaxation time of the water-water hydrogen bond correlation function, and this is in agreement with predictions of the mode coupling theory and of previous studies of the oxygen-oxygen density correlators [A. Iorio et al., J. Mol. Liq. 282, 617 (2019); Sci. China: Phys., Mech. Astron. 62, 107011 (2019)]. The water-trehalose hydrogen bond correlation function instead evidences a strong to strong crossover in the relaxation time, and this crossover is related to a trehalose dynamical transition. This signals the role that the strong interplay between the soluted molecules and the surrounding solvent has in determining the dynamical transition common to both components of the system that happens upon cooling and that is similar to the well known protein dynamical transition. We connect our results with the cryoprotecting role of trehalose molecules.

Protein packing defects "heat up" interfacial water

Ligands must displace water molecules from their corresponding protein surface binding site during association. Thus, protein binding sites are expected to be surrounded by non-tightly-bound, easily removable water molecules. In turn, the existence of packing defects at protein binding sites has been also established. At such structural motifs, named dehydrons, the protein backbone is exposed to the solvent since the intramolecular interactions are incompletely wrapped by non-polar groups. Hence, dehydrons are sticky since they depend on additional intermolecular wrapping in order to properly protect the structure from water attack. Thus, a picture of protein binding is emerging wherein binding sites should be both dehydrons rich and surrounded by easily removable water. In this work we shall indeed confirm such a link between structure and dynamics by showing the existence of a firm correlation between the degree of underwrapping of the protein chain and the mobility of the corresponding hydration water molecules. In other words, we shall show that protein packing defects promote their local dehydration, thus producing a region of "hot" interfacial water which might be easily removed by a ligand upon association.

Anomalous dynamics of water confined in MCM-41 at different hydrations

We present the results of molecular dynamics simulations on water confined in a silica pore of 15 Å modelled to mimic the MCM-41 structure. We focus on the dynamical properties of water for different hydration levels of the pore. Density profiles show a well-defined double-layer structure close to the surface. From the analysis of the layers we find that water molecules close to the substrate show an anomalous diffusion which is connected to the behaviour of the residence time distribution. The interaction with the substrate induces temporal disorder. As a consequence strong deviations are found from the Markovian processes that usually determine the long-time diffusion properties. The residence time distribution of the water molecules in the inner part of the pore, far from the hydrophilic surface, is compatible with a Markovian process.

Coupled relaxations at the protein-water interface in the picosecond time scale

The spectral behaviour of a protein and its hydration water has been investigated through neutron scattering. The availability of both hydrogenated and perdeuterated samples of maltose-binding protein (MBP) allowed us to directly measure with great accuracy the signal from the protein and the hydration water alone. Both the spectra of the MBP and its hydration water show two distinct relaxations, a behaviour that is reminiscent of glassy systems. The two components have been described using a phenomenological model that includes two Cole-Davidson functions. In MBP and its hydration water, the two relaxations take place with similar average characteristic times of approximately 10 and 0.2 ps. The common time scales of these relaxations suggest that they may be a preferential route to couple the dynamics of the water hydrogen-bond network around the protein surface with that of protein fluctuations.

Collective dynamics of protein hydration water by brillouin neutron spectroscopy

By a detailed experimental study of THz dynamics in the ribonuclease protein, we could detect the propagation of coherent collective density fluctuations within the protein hydration shell. The emerging picture indicates the presence of both a dispersing mode, traveling with a speed greater than 3000 m/s, and a nondispersing one, characterized by an almost constant energy of 6-7 meV. In agreement with molecular dynamics simulations [Phys. Rev. Lett. 2002, 89, 275501], the features of the dispersion curves closely resemble those observed in pure liquid water [Phys. Rev. E: Stat. Phys., Plasmas, Fluids, Relat. Interdiscip. Top. 2004, 69, 061203]. On the contrary, the observed damping factors are much larger than in bulk water, with the dispersing mode becoming overdamped at Q = 0.6 A(-1) already. Such novel experimental findings are discussed as a dynamic signature of the disordering effect induced by the protein surface on the local structure of water.

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.

Instantaneous normal modes and the protein glass transition

In the instantaneous normal mode method, normal mode analysis is performed at instantaneous configurations of a condensed-phase system, leading to modes with negative eigenvalues. These negative modes provide a means of characterizing local anharmonicities of the potential energy surface. Here, we apply instantaneous normal mode to analyze temperature-dependent diffusive dynamics in molecular dynamics simulations of a small protein (a scorpion toxin). Those characteristics of the negative modes are determined that correlate with the dynamical (or glass) transition behavior of the protein, as manifested as an increase in the gradient with T of the average atomic mean-square displacement at approximately 220 K. The number of negative eigenvalues shows no transition with temperature. Further, although filtering the negative modes to retain only those with eigenvectors corresponding to double-well potentials does reveal a transition in the hydration water, again, no transition in the protein is seen. However, additional filtering of the protein double-well modes, so as to retain only those that, on energy minimization, escape to different regions of configurational space, finally leads to clear protein dynamical transition behavior. Partial minimization of instantaneous configurations is also found to remove nondiffusive imaginary modes. In summary, examination of the form of negative instantaneous normal modes is shown to furnish a physical picture of local diffusive dynamics accompanying the protein glass transition.

How protein surfaces induce anomalous dynamics of hydration water

Water around biomolecules slows down with respect to pure water, and both rotation and translation exhibit anomalous time dependence in the hydration shell. The origin of such behavior remains elusive. We use molecular dynamics simulations of water dynamics around several designed protein models to establish the connection between the appearance of the anomalous dynamics and water-protein interactions. For the first time we quantify the separate effect of protein topological and energetic disorder on the hydration water dynamics. When a static protein structure is simulated, we show that both types of disorder contribute to slow down water diffusion, and that allowing for protein motion, increasing the spatial dimensionality of the interface, reduces the anomalous character of hydration water. The rotation of water is, instead, altered by the energetic disorder only; indeed, when electrostatic interactions between the protein and water are switched off, water reorients even faster than in the bulk. The dynamics of water is also related to the collective structure--à voir the hydrogen bond (H-bond) network--formed by the solvent enclosing the protein surface. We show that, as expected for a full hydrated protein, when the protein surface offers pinning sites (charged or polar sites), the superficial water-water H-bond network percolates throughout the whole surface, hindering the water diffusion, whereas it does not when the protein surface lacks electrostatic interactions with water and the water diffusion is enhanced.

Molecular view of water dynamics near model peptides

Incoherent quasi-elastic neutron scattering (QENS) has been used to measure the dynamics of water molecules in solutions of a model protein backbone, N-acetyl-glycine-methylamide (NAGMA), as a function of concentration, for comparison with results for water dynamics in aqueous solutions of the N-acetyl-leucine-methylamide (NALMA) hydrophobic peptide at comparable concentrations. From the analysis of the elastic incoherent structure factor, we find significant fractions of elastic intensity at high and low concentrations for both solutes, which corresponds to a greater population of protons with rotational time scales outside the experimental resolution (>13 ps). The higher-concentration solutions show a component of the elastic fraction that we propose is due to water motions that are strongly coupled to the solute motions, while for low-concentration solutions an additional component is activated due to dynamic coupling between inner and outer hydration layers. An important difference between the solute types at the highest concentration studied is found from stretched exponential fits to their experimental intermediate scattering functions, showing more pronounced anomalous diffusion signatures for NALMA, including a smaller stretched exponent beta and a longer structural relaxation time tau than those found for NAGMA. The more normal water diffusion exhibited near the hydrophilic NAGMA provides experimental support for an explanation of the origin of the anomalous diffusion behavior of NALMA as arising from frustrated interactions between water molecules when a chemical interface is formed upon addition of a hydrophobic side chain, inducing spatial heterogeneity in the hydration dynamics in the two types of regions of the NALMA peptide. We place our QENS measurements on model biological solutes in the context of other spectroscopic techniques and provide both confirming as well as complementary dynamic information that attempts to give a unifying molecular view of hydration dynamics signatures near peptides and proteins.

Relaxation kinetics and the glassiness of native proteins: coupling of timescales

We provide evidence that the onset of functional dynamics of folded proteins with elevated temperatures is associated with the effective sampling of its energy landscape under physiological conditions. The analysis is based on data describing the relaxation phenomena governing the backbone dynamics of bovine pancreatic trypsin inhibitor derived from molecular dynamics simulations, previously reported by us. By representing the backbone dynamics of the folded protein by three distinct regimes, it is possible to decompose its seemingly complex dynamics, described by a stretch exponential decay of the backbone motions. Of these three regimes, one is associated with the slow timescales due to the activity along the envelope of the energy surface defining the folded protein. Another, with fast timescales, is due to the activity along the pockets decorating the folded-state envelope. The intermediate regime emerges at temperatures where jumps between the pockets become possible. It is at the temperature window where motions corresponding to all three timescales become operative that the protein becomes active.

The influence of solvent composition on global dynamics of human butyrylcholinesterase powders: a neutron-scattering study

A major result of incoherent elastic neutron-scattering experiments on protein powders is the strong dependence of the intramolecular dynamics on the sample environment. We performed a series of incoherent elastic neutron-scattering experiments on lyophilized human butyrylcholinesterase (HuBChE) powders under different conditions (solvent composition and hydration degree) in the temperature range from 20 to 285 K to elucidate the effect of the environment on the enzyme atomic mean-square displacements. Comparing D(2)O- with H(2)O-hydrated samples, we were able to investigate protein as well as hydration water molecular dynamics. HuBChE lyophilized from three distinct buffers showed completely different atomic mean-square displacements at temperatures above approximately 200 K: a salt-free sample and a sample containing Tris-HCl showed identical small-amplitude motions. A third sample, containing sodium phosphate, displayed highly reduced mean-square displacements at ambient temperature with respect to the other two samples. Below 200 K, all samples displayed similar mean-square displacements. We draw the conclusion that the reduction of intramolecular protein mean-square displacements on an Angstrom-nanosecond scale by the solvent depends not only on the presence of salt ions but also on their type.

Structure, dynamics and reactions of protein hydration water

The apparent simplicity of the water molecule belies the wide range of fascinating protein phenomena in which it participates. We review recent computer simulation work on buried, internal water molecules, discussing the thermodynamics of water molecule binding and the participation of water in proton transfer reactions. Surface water molecules are also considered, with emphasis on the modification of average solvent structure on a protein surface, the role of water in the protein dynamical 'glass' transition and a simplified description of the protein motions thereby activated.

Principal components of the protein dynamical transition

Proteins exhibit a solvent-driven dynamical transition at 180-220 K, manifested by nonlinearity in the temperature dependence of the average mean-square displacement. Here, molecular dynamics simulations of hydrated myoglobin show that the onset of the transition at approximately 180 K is characterized by the appearance of a single double-well principal component mode involving a global motion of two groups of helices. As the temperature is raised a few more quasiharmonic and multiminimum components successively appear. The results indicate an underlying simplicity in the protein dynamical transition.

The role of dynamics in enzyme activity

Although protein function is thought to depend on flexibility, precisely how the dynamics of the molecule and its environment contribute to catalytic mechanisms is unclear. We review experimental and computational work relating to enzyme dynamics and function, including the role of solvent. The evidence suggests that fast motions on the 100 ps timescale, and any motions coupled to these, are not required for enzyme function. Proteins where the function is electron transfer, proton tunneling, or ligand binding may have different dynamical dependencies from those for enzymes, and enzymes with large turnover numbers may have different dynamical dependencies from those that turn over more slowly. The timescale differences between the fastest anharmonic fluctuations and the barrier-crossing rate point to the need to develop methods to resolve the range of motions present in enzymes on different time- and lengthscales.

Translational hydration water dynamics drives the protein glass transition

Experimental and computer simulation studies have revealed the presence of a glass-like transition in the internal dynamics of hydrated proteins at approximately 200 K involving an increase of the amplitude of anharmonic dynamics. This increase in flexibility has been correlated with the onset of protein activity. Here, we determine the driving force behind the protein transition by performing molecular dynamics simulations of myoglobin surrounded by a shell of water. A dual heat bath method is used with which, in any given simulation, the protein and solvent are held at different temperatures, and sets of simulations are performed varying the temperature of the components. The results show that the protein transition is driven by a dynamical transition in the hydration water that induces increased fluctuations primarily in side chains in the external regions of the protein. The water transition involves activation of translational diffusion and occurs even in simulations where the protein atoms are held fixed.

Energy resolution and dynamical heterogeneity effects on elastic incoherent neutron scattering from molecular systems

Incoherent neutron scattering is widely used to probe picosecond-nanosecond time scale dynamics of molecular systems. In systems of spatially confined atoms the relatively high intensity of elastic incoherent neutron scattering is often used to obtain a first estimate of the dynamics present. For many complex systems, however, experimental elastic scattering is difficult to interpret unambiguously using analytical dynamical models that go beyond the determination of an average mean-square displacement. To circumvent this problem a description of the scattering is derived here that encompasses a variety of analytical models in a common framework. The framework describes the time-converged part of the dynamic structure factor [the elastic incoherent scattering function (EISF)] and lends itself to practical use by explicitly incorporating effects due to the finite energy resolution of the instrument used. The dependence of the elastic scattering on wave vector is examined, and it is shown how heterogeneity in the distribution of mean-square displacements can be related to deviations of the scattering from Gaussian behavior. In this case, a correction to fourth order in the scattering vector can be used to extract the variance of the distribution of mean-square displacements. The formalism is used in a discussion of measurements on dynamics accompanying the glass transition in molecular systems. By fitting to experimental data obtained on a protein solution the present methodology is used to show how the existence of a temperature-dependent relaxation frequency can lead to a transition in the measured mean-square displacement in the absence of an EISF change.

Ferricytochrome c encapsulated in silica hydrogels: correlation between active site dynamics and solvent structure

Ferricytochrome c encapsulated in silica hydrogels has been prepared by the sol-gel technique following, with some modifications, the procedure originally developed by Zink et al. A suitable preparation of hydrogels enables to have both 'wet' and 'dry' samples. Wet samples have a high water content: as the temperature is lowered below approximately 260 K water freezes and the samples crack. On the contrary, dry samples have a low water content (hydration h approximately 0.35): in these conditions water does not freeze even at cryogenic temperatures and the samples remain transparent and non-cracking. The dynamics of ferricytochrome c and its dependence on the surrounding medium have been studied by optical absorption spectroscopy in the temperature range 10-300 K. At each temperature, spectra were collected both in the Soret region and in the near infrared at approximately 1.45 microm (the water overtone band); this enables to probe the local dynamics of the protein active site as well as the 'structure' of water molecules present in the sample. The data show that sol-gel encapsulation 'per se' does not alter the protein active site dynamics, but rather introduces an increased local heterogeneity. At difference, we find a correlation between active site dynamics and water structure: in the wet hydrogel, freezing of water quenches the ensemble of soft modes linearly coupled to the Soret transition; while, in the dry hydrogel, water does not freeze, and an active site dynamic behavior-similar to the non-freezing water/glycerol solution-is observed.

Ferricytochrome c encapsulated in silica hydrogels: correlation between active site dynamics and solvent structure

Ferricytochrome c encapsulated in silica hydrogels has been prepared by the sol-gel technique following, with some modifications, the procedure originally developed by Ellerby et al. (Science 255 1113 (1992)). A suitable preparation of hydrogels enables having both 'wet' and 'dry' samples. Wet samples have a high water content: as the temperature is lowered below approximately 260 K, water freezes and the samples crack. On the contrary, dry samples have a low water content (hydration h approximately equal 0.35): in these conditions water does not freeze even at cryogenic temperatures and the samples remain transparent and non-cracking. The dynamics of ferricytochrome c and its dependence on the surrounding medium have been studied by optical absorption spectroscopy in the temperature range 10-300 K. At each temperature, spectra were collected both in the Soret region and in the near infrared at approximately 1.45 microm (the water overtone band); this enables probing the local dynamics of the protein active site as well as the 'structure' of water molecules present in the sample. The data show that sol-gel encapsulation 'per se' does not alter the protein active site dynamics, but rather introduces an increased local heterogeneity. We find a correlation between active site dynamics and water structure: in the wet hydrogel, freezing of water quenches the ensemble of soft modes linearly coupled to the Soret transition; while, in the dry hydrogel, water does not freeze and an active site dynamic behavior--similar to the non-freezing water/glycerol solution--is observed.

Relaxation kinetics and the glassiness of proteins: the case of bovine pancreatic trypsin inhibitor

Folded proteins may be regarded as soft active matter under physiological conditions. The densely packed hydrophobic interior, the relatively molten hydrophilic exterior, and the spacer connecting these put together a large number of locally homogeneous regions. For the case of the bovine pancreatic trypsin inhibitor, with the aid of molecular dynamics simulations, we have demonstrated that the kinetics of the relaxation of the internal motions is highly concerted, manifesting the protein's heterogeneity, which may arise from variations in density, local packing, or the local energy landscape. This behavior is characterized in a stretched exponential decay described by an exponent of approximately 0.4 at physiological temperatures. Due to the trapped conformations, configurational entropy becomes smaller, and the associated stretch exponent drops to half of its value below the glass transition range. The temperature dependence of the inverse relaxation time closely follows the Vogel-Tamman-Fulcher expression when the protein is biologically active.

Vibrational dynamics of myoglobin determined by the phonon-assisted Mössbauer effect

The phonon-assisted Mössbauer effect is used to determine the partial phonon density of states of the iron within the active center of deoxymyoglobin, carboxymyoglobin, and dry and wet metmyoglobin between 40 and 300 K. Between 0 and 1 meV the iron density of states increases quadratically with the energy, as in a Debye solid. Mean sound velocities are extracted from this slope. Between 1 and 3 meV a nearly quadratic "Debye-like" increase follows due to the similar strength of intermolecular and intramolecular forces. Above 3 meV, optical vibrations are characteristic for the iron-ligand conformation. The overall mean square displacements of the heme iron atom obtained from the density of states agree well with the values of Mössbauer absorption experiments below 180 K. In the physiological temperature regime the data confirm the existence of harmonic vibrations in addition to the protein specific dynamics measured by Mössbauer absorption. In the Debye energy regime the mean square displacement of the iron is in agreement with that of the hydrogens measured by incoherent neutron scattering demonstrating the global character of these modes. At higher energies the vibration of the heavy iron atom at 33 meV in metmyoglobin is as large as that of the lightweight hydrogens at that energy. A freeze dried, rehydrated (h=0.38 g H2O/g protein) metmyoglobin sample shows an excess of states above the Debye law between 1 and 3 meV, similar to neutron scattering experiments. The room temperature density of states below 3 meV exhibit an increase of the density compared to the low temperature data, which can be interpreted as mode softening.

Glassy dynamics and enzymatic activity of lysozyme

There has been much interest in the analogies between dynamic processes in proteins and in other complex systems such as viscous liquids and glasses. We have investigated the dynamics of protons along chains of hydrogen-bonded water molecules adsorbed on the surface of the globular protein lysozyme. The hydration dependence of the dielectric relaxation time is fitted by a modified Vogel-Fulcher-Tamman equation, in which the variable temperature has been replaced with hydration. We find that the relaxation time diverges at a singular hydration that coincides with the critical water content required to trigger lysozyme enzymatic activity. This surprising correlation suggests a direct coupling between protein function and glasslike behavior, with possible implications for the turnover number of the enzyme.