Puzzle of protein dynamical transition

An experimentally representative in-silico protocol for dynamical studies of lyophilised and weakly hydrated amorphous proteins

Characterization of biopolymers in both dry and weakly hydrated amorphous states has implications for the pharmaceutical industry since it provides understanding of the effect of lyophilisation on stability and biological activity. Atomistic Molecular Dynamics (MD) simulations probe structural and dynamical features related to system functionality. However, while simulations in homogenous aqueous environments are routine, dehydrated model assemblies are a challenge with systems investigated in-silico needing careful consideration; simulated systems potentially differing markedly despite seemingly negligible changes in procedure. Here we propose an in-silico protocol to model proteins in lyophilised and weakly hydrated amorphous states that is both more experimentally representative and routinely applicable. Since the outputs from MD align directly with those accessed by neutron scattering, the efficacy of the simulation protocol proposed is shown by validating against experimental neutron data for apoferritin and insulin. This work also highlights that without cooperative experimental and simulative data, development of simulative procedures using MD alone would prove most challenging.

Structure and dynamics of supercooled water in the hydration layer of poly(ethylene glycol)

The statics and dynamics of supercooled water in the hydration layer of poly(ethylene glycol) (PEG) were studied by a combination of quasi-elastic neutron scattering (QENS) and molecular dynamics (MD) simulations. Two samples, that is, hydrogenated PEG/deuterated water (h-PEG/D₂O) and fully deuterated PEG/hydrogenated water (d-PEG/H₂O) with the same molar ratio of ethylene glycol (EG) monomer to water, 1:1, are compared. The QENS data of h-PEG/D₂O show the dynamics of PEG, and that of d-PEG/H₂O reveals the motion of water. The temperature-dependent elastic scattering intensity of both samples has shown transitions at supercooled temperature, and these transition temperatures depend on the energy resolution of the instruments. Therefore, neither one is a phase transition, but undergoes dynamic process. The dynamic of water can be described as an Arrhenius to super-Arrhenius transition, and it reveals the hydrogen bonding network relaxation of hydration water around PEG at supercooled temperature. Since the PEG-water hydrogen bond structural relaxation time from MD is in good agreement with the average relaxation time from QENS (d-PEG/H₂O), MD may further reveal the atomic pictures of the supercooled hydration water. It shows that hydration water molecules form a series of pools around the hydrophilic oxygen atom of PEG. At supercooled temperature, they have a more bond ordered structure than bulk water, proceed a trapping sites diffusion on the PEG surface, and facilitate the structural relaxation of PEG backbone.

Temperature-Dependent Fluorescence of mPlum Fluorescent Protein from 295 to 20 K

The development of bright fluorescent proteins (FPs) emitting beyond 600 nm continues to be of interest both from a fundamental perspective in understanding protein-chromophore interactions and from a practical perspective as these FPs would be valuable for cellular imaging. We previously reported ultrafast spectral observations of the excited-state dynamics in mPlum resulting from interconversion between direct hydrogen bonding and water-mediated hydrogen bonding between the chromophore acylimine carbonyl and the Glu16 side chain. Here, we report temperature-dependent steady-state and time-resolved fluorescence measurements of mPlum and its E16H variant, which does not contain a side-chain permitting hydrogen bonding with the acylimine carbonyl. Lowering the temperature of the system freezes interconversion between the hydrogen-bonding states, thus revealing the spectral signatures of the two states. Analysis of the temperature-dependent spectra assuming Boltzmann populations of the two states yields a 205 cm^-1 energy difference. This value agrees with the predictions from a quantum mechanics/molecular mechanics study of mPlum (198 cm^-1). This study demonstrates the first use of cryogenic spectroscopy to quantify the energetics and timescales of FP chromophore structural states that were only previously obtained from computational methods and further confirms the importance of acylimine hydrogen-bonding dynamics to the fluorescence spectral shifts of red FPs.

Wide-angle X-ray scattering and molecular dynamics simulations of supercooled protein hydration water

Understanding the mechanism responsible for the protein low-temperature crossover observed at T≈ 220 K can help us improve current cryopreservation technologies. This crossover is associated with changes in the dynamics of the system, such as in the mean-squared displacement, whereas experimental evidence of structural changes is sparse. Here we investigate hydrated lysozyme proteins by using a combination of wide-angle X-ray scattering and molecular dynamics (MD) simulations. Experimentally we suppress crystallization by accurate control of the protein hydration level, which allows access to temperatures down to T = 175 K. The experimental data indicate that the scattering intensity peak at Q = 1.54 Å-1, attributed to interatomic distances, exhibits temperature-dependent changes upon cooling. In the MD simulations it is possible to decompose the water and protein contributions and we observe that, while the protein component is nearly temperature independent, the hydration water peak shifts in a fashion similar to that of bulk water. The observed trends are analysed by using the water-water and water-protein radial distribution functions, which indicate changes in the local probability density of hydration water.

Experimental demonstration of the novel "van-Hove integral method (vHI)" for measuring diffusive dynamics by elastic neutron scattering

Quasi-elastic neutron scattering (QENS)-based on the seminal work of Nobel Laureate Brockhouse-has been one of the major methods for studying pico-second to nano-second diffusive dynamics over the past 70 years. This is regarded as an "inelastic" method for dynamics. In contrast, we recently proposed a new neutron-scattering method for dynamics, which uses the elastic line of the scattering to access system dynamics directly in the time domain (Benedetto and Kearley in Sci Rep 9:11284, 2019). This new method has been denoted "vHI" that stands for "van Hove Integral". The reason is that, under certain conditions, the measured elastic intensity corresponds to the running-time integral of the intermediate scattering function, [Formula: see text], up to a time that is inversely proportional to the energy band-width incident on the sample. As a result, [Formula: see text] is accessed from the time derivative of the measured vHI profile. vHI has been supported by numerical and Monte-Carlo simulations, but has been difficult to validate experimentally due to the lack of a suitable instrument. Here we show that vHI works in practice, which we achieved by using a simple modification to the standard QENS backscattering spectrometer methodology. Basically, we varied the neutron-energy band-widths incident at the sample via a step-wise variation of the frequency of the monochromator Doppler-drive. This provides a measurement of the vHI profile at the detectors. The same instrument and sample were also used in standard QENS mode for comparison. The intermediate scattering functions, [Formula: see text], obtained by the two methods-vHI and QENS-are strikingly similar providing a direct experimental validation of the vHI method. Perhaps surprisingly, the counting statistics of the two methods are comparable even though the instrument used was expressly designed for QENS. This shows that the methodology modification adopted here can be used in practice to access vHI profiles at many of the backscattering spectrometers worldwide. We also show that partial integrations of the measured QENS spectrum cannot provide the vHI profile, which clarifies a common misconception. At the same time, we show a novel approach which does access [Formula: see text] from QENS spectra.

Freezable and Unfreezable Hydration Water: Distinct Contributions to Protein Dynamics Revealed by Neutron Scattering

Hydration water plays a crucial role for activating the protein dynamics required for functional expression. Yet, the details are not understood about how hydration water couples with protein dynamics. A temperature hysteresis of the ice formation of hydration water is a key phenomenon to understand which type of hydration water, unfreezable or freezable hydration water, is crucial for the activation of protein dynamics. Using neutron scattering, we observed a temperature-hysteresis phenomenon in the diffraction peaks of the ice of freezable hydration water, whereas protein dynamics did not show any temperature hysteresis. These results show that the protein dynamics is not coupled with freezable hydration water dynamics, and unfreezable hydration water is essential for the activation of protein dynamics. Decoupling of the dynamics between unfreezable and freezable hydration water could be the cause of the distinct contributions of hydration water to protein dynamics.

From protein and its hydration water dynamics to controlling mechano-elasticity of cellular lipid membranes and cell migration via ionic liquids

In this invited Commentary, as requested, I will walk the reader through my research path starting from my first works on proteins and their hydration water dynamics to my most recent activity on the use of ionic liquids (ILs) as molecular handles to control and manipulate cell membrane mechano-elasticity and cell migration. In doing so I will comment on my research activity on polymers, proteins, natural bioprotectants, phospholipid bilayers, amyloids and cells, which I have carried out by combining several different experimental and computational approaches including neutron scattering, atomic force microscopy, classical molecular dynamics and ab initio calculations, used in tandem with several biological assays and a palette of complementary techniques ranging from calorimetry to static and dynamic light scattering. In parallel to this biophysical/chemical-physical core activity, a smaller portion of my interest and effort has been-I may say always-dedicated to the development of a new neutron scattering method/spectroscopy for dynamics based on "elastic" scattering. I will comment on this instrumental side of my research as well and show its relationship with the biophysical core of my activity. The overall picture that emerges is, from my personal prospective, of a coherent 13-year research path based on curiosity and a problem-solving approach, in which the fundamental importance of inter- and trans-disciplinary approaches and collaborations is emerging on the way, forecasting a prosper and intriguing future for physics in biology and in nanomedicine and bionanotechnology applications.

Advancing predictions of protein stability in the solid state

The β-relaxation associated with the sub-glass transition temperature (T_g,β) is attributed to fast, localised molecular motions which can occur below the primary glass transition temperature (T_g,α). Consistent with T_g,β being observed well-below storage temperatures, the β-relaxation associated motions have been hypothesised to influence protein stability in the solid state and could thus impact the quality of e.g. protein powders for inhalation or reconstitution and injection. Why then do distinct solid state protein formulations with similar aggregation profiles after drying and immediate reconstitution, display different profiles when reconstituted following prolonged storage? Is the value of T_g,β, associated with the β-relaxation process of the system, a reliable parameter for characterising the behaviour of proteins in the solid state? Bearing this in mind, in this work we further explore the different relaxation dynamics of glassy solid state monoclonal antibody formulations using terahertz time-domain spectroscopy and dynamical mechanical analysis. By conducting a 52 week stability study on a series of multi-component spray-dried formulations, an approach for characterising and analysing the solid state dynamics and how these relate to protein stability is outlined.

Quasielastic neutron scattering studies on couplings of protein and water dynamics in hydrated elastin

Performing quasielastic neutron scattering measurements and analyzing both elastic and quasielasic contributions, we study protein and water dynamics of hydrated elastin. At low temperatures, hydration-independent methyl group rotation dominates the findings. It is characterized by a Gaussian distribution of activation energies centered at about E_m = 0.17 eV. At ∼195 K, coupled protein-water motion sets in. The hydration water shows diffusive motion, which is described by a Gaussian distribution of activation energies with E_m = 0.57 eV. This Arrhenius behavior of water diffusion is consistent with previous results for water reorientation, but at variance with a fragile-to-strong crossover at ∼225 K. The hydration-related elastin backbone motion is localized and can be attributed to the cage rattling motion. We speculate that its onset at ∼195 K is related to a secondary glass transition, which occurs when a β relaxation of the protein has a correlation time of τ_β ∼ 100 s. Moreover, we show that its temperature-dependent amplitude has a crossover at the regular glass transition T_g = 320 K of hydrated elastin, where the α relaxation of the protein obeys τ_α ∼ 100 s. By contrast, we do not observe a protein dynamical transition when water dynamics enters the experimental time window at ∼240 K.

A Quantitative Comparison of the Counting Significance of van Hove Integral Spectroscopy and Quasielastic Neutron Scattering

We have recently proposed a new method to access system dynamics via neutron scattering based on measuring the elastic scattered intensity: By varying the energy band-width that impinges on the sample (also known as instrumental energy resolution), the purely elastic-scattering from this variation is the running time-integral of the intermediate scattering function (I(t)) [Benedetto and Kearley, Sci. Rep. 9, 11284, 2019]. In this correspondence we denote our method "vHI", which stands for "van Hove Integral". The method is now widely accepted as "valid" and here we focus on the efficiency of the vHI method compared with the standard quasi-elastic neutron scattering (QENS) method. We use a numerical Monte-Carlo simulation of an instrument that is equally capable of measuring QENS and vHI under identical conditions. For an "experiment" in which the same number of neutrons enter the instrument, we present comparisons between QENS and vHI at three levels of data-reduction. Firstly, at the raw-data level vHI achieves 100 times more neutrons at the detector than QENS. Secondly, vHI has a factor of 2 less statistical error, which would translate to an overall gain of 4 for vHI in counting-time. Lastly, we compare the distortions caused in obtaining the final I(t) via time-Fourier transform (QENS) and polynomial time-derivative (vHI). Here, the statistical error is 10 times smaller for vHI. This last comparison is the most important result where the 10 times smaller residual for vHI gives a net gain in counting time of 100 better than QENS to obtain the same underlying dynamics of the system under study.

Role of hydration water in the onset of protein structural dynamics

Proteins are the molecular workhorses in a living organism. Their 3D structures are animated by a multitude of equilibrium fluctuations and specific out-of-equilibrium motions that are required for proteins to be biologically active. When studied as a function of temperature, functionally relevant dynamics are observed at and above the so-called protein dynamical transition (~240 K) in hydrated, but not in dry proteins. In this review we present and discuss the main experimental and computational results that provided evidence for the dynamical transition, with a focus on the role of hydration water dynamics in sustaining functional protein dynamics. The coupling and mutual influence of hydration water dynamics and protein dynamics are discussed and the hypotheses illustrated that have been put forward to explain the physical origin of their onsets.

Dynamics from elastic neutron-scattering via direct measurement of the running time-integral of the van Hove distribution function

We present a new neutron-scattering approach to access the van Hove distribution function directly in the time domain, I(t), which reflects the system dynamics. Currently, I(t) is essentially determined from neutron energy-exchange. Our method consists of the straightforward measurement of the running time-integral of I(t), by computing the portion of scattered neutrons corresponding to species at rest within a time t, (conceptually elastic scattering). Previous attempts failed to recognise this connection. Starting from a theoretical standpoint, a practical realisation is assessed via numerical methods and an instrument simulation.

Dynamical transition in molecular glasses and proteins observed by spin relaxation of nitroxide spin probes and labels

In glassy substances and biological media, dynamical transitions are observed in neutron scattering that manifests itself as deviations of the translational mean-squared displacement, 〈x²〉, of hydrogen atoms from harmonic dynamics. In biological media, the deviation occurs at two temperature intervals, at ∼100-150 K and at ∼170-230 K, and it is attributed to the motion of methyl groups in the former case and to the transition from harmonic to anharmonic or diffusive motions in the latter case. In this work, electron spin echo (ESE) spectroscopy-a pulsed version of electron paramagnetic resonance-is applied to study the spin relaxation of nitroxide spin probes and labels introduced in molecular glass former o-terphenyl and in protein lysozyme. The anisotropic contribution to the rate of the two-pulse ESE decay, ΔW, is induced by spin relaxation appearing because of restricted orientational stochastic molecular motion; it is proportional to 〈α²〉τ_c, where 〈α²〉 is the mean-squared angle of reorientation of the nitroxide molecule around the equilibrium position and τ_c is the correlation time of reorientation. The ESE time window allows us to study motions with τ_c < 10^-7 s. For glassy o-terphenyl, the 〈α²〉τ_c temperature dependence shows a transition near 240 K, which is in agreement with the literature data on 〈x²〉. For spin probes of essentially different size, the obtained data were found to be close, which evidences that motion is cooperative, involving a nanocluster of several neighboring molecules. For the dry lysozyme, the 〈α²〉τ_c values below 260 K were found to linearly depend on the temperature in the same way as it was observed in neutron scattering for 〈x²〉. As spin relaxation is influenced only by stochastic motion, the harmonic motions seen in ESE must be overdamped. In the hydrated lysozyme, ESE data show transitions near 130 K for all nitroxides, near 160 K for the probe located in the hydration layer, and near 180 K for the label in the protein interior. For this system, the two latter transitions are not observed in neutron scattering. The ESE-detected transitions are suggested to be related with water dynamics in the nearest hydration shell: with water glass transition near 130 K and with the onset of overall water molecular reorientations near 180 K; the disagreement with neutron scattering is ascribed to the larger time window for ESE-detected motions.

Slow Dynamics around a Protein and Its Coupling to Solvent

Solvent is essential for protein dynamics and function, but its role in regulating the dynamics remains debated. Here, we employ saturation transfer electron spin resonance (ST-ESR) to explore the issue and characterize the dynamics on a longer (from μs to s) time scale than has been extensively studied. We first demonstrate the reliability of ST-ESR by showing that the dynamical changeovers revealed in the spectra agree to liquid-liquid transition (LLT) in the state diagram of the glycerol/water system. Then, we utilize ST-ESR with four different probes to systematically map out the variation in local (site-specific) dynamics around a protein surface at subfreezing temperatures (180-240 K) in 10 mol % glycerol/water mixtures. At highly exposed sites, protein and solvent dynamics are coupled, whereas they deviate from each other when temperature is greater than LLT temperature (∼190 K) of the solvent. At less exposed sites, protein however exhibits a dynamic, which is distinct from the bulk solvent, throughout the temperature range studied. Dominant dynamic components are thus revealed, showing that (from low to high temperatures) the overall structural fluctuation, rotamer dynamics, and internal side-chain dynamics, in turn, dominate the temperature dependence of spin-label motions. The structural fluctuation component is relatively slow, collective, and independent of protein structural segments, which is thus inferred to a fundamental dynamic component intrinsic to protein. This study corroborates that bulk solvent plasticizes protein and facilitates rather than slaves protein dynamics.

Dynamic Neutron Scattering by Biological Systems

Dynamic neutron scattering directly probes motions in biological systems on femtosecond to microsecond timescales. When combined with molecular dynamics simulation and normal mode analysis, detailed descriptions of the forms and frequencies of motions can be derived. We examine vibrations in proteins, the temperature dependence of protein motions, and concepts describing the rich variety of motions detectable using neutrons in biological systems at physiological temperatures. New techniques for deriving information on collective motions using coherent scattering are also reviewed.

Hydrogen-bond dynamics at the bio-water interface in hydrated proteins: a molecular-dynamics study

Water is fundamental to the biochemistry of enzymes. It is well known that without a minimum amount of water, enzymes are not biologically active. Bare minimal solvation for biological function corresponds to about a single layer of water covering enzymes' surfaces. Many contradictory studies on protein-hydration-water-coupled dynamics have been published in recent decades. Following prevailing wisdom, a dynamical crossover in hydration water (at around 220 K for hydrated lysozymes) can trigger larger-amplitude motions of the protein, activating, in turn, biological functions. Here, we present a molecular-dynamics-simulation study on a solvated model protein (hen egg-white lysozyme), in which we determine, inter alia, the relaxation dynamics of the hydrogen-bond network between the protein and its hydration water molecules on a residue-per-residue basis. Hydrogen-bond breakage/formation kinetics is rather heterogeneous in temperature dependence (due to the heterogeneity of the free-energy surface), and is driven by the magnitude of thermal motions of various different protein residues which provide enough thermal energy to overcome energy barriers to rupture their respective hydrogen bonds with water. In particular, arginine residues exhibit the highest number of such hydrogen bonds at low temperatures, losing almost completely such bonding above 230 K. This suggests that hydration water's dynamical crossover, observed experimentally for hydrated lysozymes at ∼220 K, lies not at the origin of the protein residues' larger-amplitude motions, but rather arises as a consequence thereof. This highlights the need for new experimental investigations, and new interpretations to link protein dynamics to functions, in the context of key interrelationships with the solvation layer.

The effects of pressure on the energy landscape of proteins

Protein dynamics is characterized by fluctuations among different conformational substates, i.e. the different minima of their energy landscape. At temperatures above ~200 K, these fluctuations lead to a steep increase in the thermal dependence of all dynamical properties, phenomenon known as Protein Dynamical Transition. In spite of the intense studies, little is known about the effects of pressure on these processes, investigated mostly near room temperature. We studied by neutron scattering the dynamics of myoglobin in a wide temperature and pressure range. Our results show that high pressure reduces protein motions, but does not affect the onset temperature for the Protein Dynamical Transition, indicating that the energy differences and barriers among conformational substates do not change with pressure. Instead, high pressure values strongly reduce the average structural differences between the accessible conformational substates, thus increasing the roughness of the free energy landscape of the system.

Ergodicity breaking of iron displacement in heme proteins

We present a model of the dynamical transition of atomic displacements in proteins. Increased mean-square displacement at higher temperatures is caused by the softening of the force constant for atomic/molecular displacements by electrostatic and van der Waals forces from the protein-water thermal bath. Displacement softening passes through a nonergodic dynamical transition when the relaxation time of the force-force correlation function enters, with increasing temperature, the instrumental observation window. Two crossover temperatures are identified. The lower crossover, presently connected to the glass transition, is related to the dynamical unfreezing of rotations of water molecules within nanodomains polarized by charged surface residues of the protein. The higher crossover temperature, usually assigned to the dynamical transition, marks the onset of water translations. All crossovers are ergodicity breaking transitions depending on the corresponding observation windows. Allowing stretched exponential relaxation of the protein-water thermal bath significantly improves the theory-experiment agreement when applied to solid protein samples studied by Mössbauer spectroscopy.

Low-Temperature Decoupling of Water and Protein Dynamics Measured by Neutron Scattering

Water plays a major role in biosystems, greatly contributing to determine their structure, stability, and function. It is well known, for instance, that proteins require a minimum amount of water to be fully functional. Despite many years of intensive research, however, the detailed nature of protein-hydration water interactions is still partly unknown. The widely accepted "protein dynamical transition" scenario is based on perfect coupling between the dynamics of proteins and that of their hydration water, which has never been probed in depth experimentally. I present here high-resolution elastic neutron scattering measurements of the atomistic dynamics of lysozyme in water. The results show for the first time that the dynamics of proteins and of their hydration water are actually decoupled at low temperatures. This important result challenges the "protein dynamical transition" scenario and requires a new model to link protein dynamics to the dynamics of its hydration water.

Solvent-Driven Dynamical Crossover in the Phenylalanine Side-Chain from the Hydrophobic Core of Amyloid Fibrils Detected by 2H NMR Relaxation

Aromatic residues are important markers of dynamical changes in proteins' hydrophobic cores. In this work we investigated the dynamics of the F19 side-chain in the core of amyloid fibrils across a wide temperature range of 300 to 140 K. We utilized solid-state 2H NMR relaxation to demonstrate the presence of a solvent-driven dynamical crossover between different motional regimes, often also referred to as the dynamical transition. In particular, the dynamics are dominated by small-angle fluctuations at low temperatures and by π-flips of the aromatic ring at high temperatures. The crossover temperature is more than 43 degrees lower for the hydrated state of the fibrils compared to the dry state, indicating that interactions with water facilitate π-flips. Further, crossover temperatures are shown to be very sensitive to polymorphic states of the fibrils, such as the 2-fold and 3-fold symmetric morphologies of the wild-type protein as well as D23N mutant protofibrils. We speculate that these differences can be attributed, at least partially, to enhanced interactions with water in the 3-fold polymorph, which has been shown to have a water-accessible cavity. Combined with previous studies of methyl group dynamics, the results highlight the presence of multiple dynamics modes in the core of the fibrils, which was originally believed to be quite rigid.

Solvent-Slaved Dynamic Processes Observed by Tryptophan Phosphorescence of Human Serum Albumin

Despite extensive experimental and computational efforts to understand the nature of the hierarchy of protein fluctuations and the modulating role of the protein hydration shell, a detailed microscopic description of the dynamics of the protein-solvent system has yet to be achieved. By using single tryptophan protein phosphorescence, we follow site-specific internal protein dynamics over a broad temperature range and demonstrate three independent dynamic processes. Process I is seen at temperatures below the bulk solvent T_g, has low activation energy, and is likely due to fast vibrations that may be enabled by water mobility on the protein surface. Process II is observed above 170 K, with activation energy typical of β relaxations in a glass; it has the same temperature dependence as fluctuations of hydration shell waters. Process III is observed at T > 200 K; it has super-Arrhenius temperature dependence and closely follows the primary relaxation of the bulk. The fluorescence of pyranine bound to the protein reports on the mobility of water in the hydration shell; it reveals a shift in emission spectra with increasing temperature, indicative of a changing H-bond network at the surface of the protein. These results support a model of solvent-slaved protein dynamics.

Perturbation of hydration layer in solvated proteins by external electric and electromagnetic fields: Insights from non-equilibrium molecular dynamics

Given the fundamental role of water in governing the biochemistry of enzymes, and in regulating their wider biological activity (e.g., by local water concentration surrounding biomolecules), the influence of extraneous electric and electromagnetic (e/m) fields thereon is of central relevance to biophysics and, more widely, biology. With the increase in levels of local and atmospheric microwave-frequency radiation present in modern life, as well as other electric-field exposure, the impact upon hydration-water layers surrounding proteins, and biomolecules generally, becomes a particularly pertinent issue. Here, we present a (non-equilibrium) molecular-dynamics-simulation study on a model protein (hen egg-white lysozyme) hydrated in water, in which we determine, inter alia, translational self-diffusivities for both hen egg-white lysozyme and its hydration layer together with relaxation dynamics of the hydrogen-bond network between the protein and its hydration-layer water molecules on a residue-per-residue basis. Crucially, we perform this analysis both above and below the dynamical-transition temperature (at ∼220 K), at 300 and 200 K, respectively, and we compare the effects of external static-electric and e/m fields with linear-response-régime (r.m.s.) intensities of 0.02 V/Å. It was found that the translational self-diffusivity of hen egg-white lysozyme and its hydration-water layer are increased substantially in static fields, primarily due to the induced electrophoretic motion, whilst the water-protein hydrogen-bond-network-rearrangement kinetics can also undergo rather striking accelerations, primarily due to the enhancement of a larger-amplitude local translational and rotational motion by charged and dipolar residues, which serves to promote hydrogen-bond breakage and re-formation kinetics. These external-field effects are particularly evident at 200 K, where they serve to induce the protein- and solvation-layer-response effects redolent of dynamical transition at a lower temperature (∼200 K) vis-à-vis the zero-field case (∼220 K).

Hydrogen Mobility and Protein-Water Interactions in Proteins in the Solid State

In this work the groundwork is laid for characterizing the mobility of hydrogen-hydrogen pairs (proton-proton radial vectors) in proteins in the solid state that contain only residual water. In this novel approach, we introduce new ways of analyzing and interpreting data: 1) by representing hydrogen mobility (HM) and melting diagram (MD) data recorded by wide-line ¹ H NMR spectroscopic analysis as a function of fundamental temperature (thermal excitation energy); 2) by suggesting a novel mode of interpretation of these parameters that sheds light on details of protein-water interactions, such as the exact amount of water molecules and the distribution of barrier potentials pertaining to their rotational and surface translational mobility; 3) by relying on directly determined physical observables. We illustrate the power of this approach by studying the behavior of two proteins, the structured enzyme lysozyme and the intrinsically disordered ERD14.

Investigating Structure and Dynamics of Proteins in Amorphous Phases Using Neutron Scattering

In order to increase shelf life and minimize aggregation during storage, many biotherapeutic drugs are formulated and stored as either frozen solutions or lyophilized powders. However, characterizing amorphous solids can be challenging with the commonly available set of biophysical measurements used for proteins in liquid solutions. Therefore, some questions remain regarding the structure of the active pharmaceutical ingredient during freezing and drying of the drug product and the molecular role of excipients. Neutron scattering is a powerful technique to study structure and dynamics of a variety of systems in both solid and liquid phases. Moreover, neutron scattering experiments can generally be correlated with theory and molecular simulations to analyze experimental data. In this article, we focus on the use of neutron techniques to address problems of biotechnological interest. We describe the use of small-angle neutron scattering to study the solution structure of biological molecules and the packing arrangement in amorphous phases, that is, frozen glasses and freeze-dried protein powders. In addition, we discuss the use of neutron spectroscopy to measure the dynamics of glassy systems at different time and length scales. Overall, we expect that the present article will guide and prompt the use of neutron scattering to provide unique insights on many of the outstanding questions in biotechnology.

Fast Motions of Key Methyl Groups in Amyloid-β Fibrils

Amyloid-β (Aβ) peptide is the major component of plaques found in Alzheimer's disease patients. Using solid-state 2H NMR relaxation performed on selectively deuterated methyl groups, we probed the dynamics in the threefold symmetric and twofold symmetric polymorphs of native Aβ as well as the protofibrils of the D23N mutant. Specifically, we investigated the methyl groups of two leucine residues that belong to the hydrophobic core (L17 and L34) as well as M35 residues belonging to the hydrophobic interface between the cross-β subunits, which has been previously found to be water-accessible. Relaxation measurements performed over 310-140 K and two magnetic field strengths provide insights into conformational variability within and between polymorphs. Core packing variations within a single polymorph are similar to what is observed for globular proteins for the core residues, whereas M35 exhibits a larger degree of variability. M35 site is also shown to undergo a solvent-dependent dynamical transition in which slower amplitude motions of methyl axes are activated at high temperature. The motions, modeled as a diffusion of methyl axis, have activation energy by a factor of 2.7 larger in the twofold compared with the threefold polymorph, whereas D23N protofibrils display a value similar to the threefold polymorph. This suggests enhanced flexibility of the hydrophobic interface in the threefold polymorph. This difference is only observed in the hydrated state and is absent in the dry fibrils, highlighting the role of solvent at the cavity. In contrast, the dynamic behavior of the core is hydration-independent.

Elastic Scattering Spectroscopy (ESS): an Instrument-Concept for Dynamics of Complex (Bio-) Systems From Elastic Neutron Scattering

A new type of neutron-scattering spectroscopy is presented that is designed specifically to measure dynamics in bio-systems that are difficult to obtain in any other way. The temporal information is largely model-free and is analogous to relaxation processes measured with dielectric spectroscopy, but provides additional spacial and geometric aspects of the underlying dynamics. Numerical simulations of the basic instrument design show the neutron beam can be highly focussed, giving efficiency gains that enable the use of small samples. Although we concentrate on continuous neutron sources, the extension to pulsed neutron sources is proposed, both requiring minimal data-treatment and being broadly analogous with dielectric spectroscopy, they will open the study of dynamics to new areas of biophysics.

Protein dynamics as seen by (quasi) elastic neutron scattering

BACKGROUND

Elastic and quasielastic neutron scattering studies proved to be efficient probes of the atomic mean square displacement (MSD), a fundamental parameter for the characterization of the motion of individual atoms in proteins and its evolution with temperature and compositional environment.

SCOPE OF REVIEW

We present a technical overview of the different types of experimental situations and the information quasi-elastic neutron scattering approaches can make available. In particular, MSD can crucially depend on the time scale over which the averaging (building of the "mean") takes place, being defined by the instrumental resolution. Due to their high neutron scattering cross section, hydrogen atoms can be particularly sensitively observed with little interference by the other atoms in the sample. A few examples, including new data, are presented for illustration.

MAJOR CONCLUSIONS

The incoherent character of neutron scattering on hydrogen atoms restricts the information obtained to the self-correlations in the motion of individual atoms, simplifying at the same time the data analysis. On the other hand, the (often overlooked) exploration of the averaging time dependent character of MSD is crucial for unambiguous interpretation and can provide a wealth of information on micro- and nanoscale atomic motion in proteins.

GENERAL SIGNIFICANCE

By properly exploiting the broad range capabilities of (quasi)elastic neutron scattering techniques to deliver time dependent characterization of atomic displacements, they offer a sensitive, direct and simple to interpret approach to exploration of the functional activity of hydrogen atoms in proteins. Partial deuteration can add most valuable selectivity by groups of hydrogen atoms. "This article is part of a Special Issue entitled "Science for Life" Guest Editor: Dr. Austen Angell, Dr. Salvatore Magazù and Dr. Federica Migliardo".

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

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

Dynamical and orientational structural crossovers in low-temperature glycerol

Mean-square displacements of hydrogen atoms in glass-forming materials and proteins, as reported by incoherent elastic neutron scattering, show kinks in their temperature dependence. This crossover, known as the dynamical transition, connects two approximately linear regimes. It is often assigned to the dynamical freezing of subsets of molecular modes at the point of equality between their corresponding relaxation times and the instrumental observation window. The origin of the dynamical transition in glass-forming glycerol is studied here by extensive molecular dynamics simulations. We find the dynamical transition to occur for both the center-of-mass translations and the molecular rotations at the same temperature, insensitive to changes of the observation window. Both the translational and rotational dynamics of glycerol show a dynamic crossover from the structural to a secondary relaxation at the temperature of the dynamical transition. A significant and discontinuous increase in the orientational Kirkwood factor and in the dielectric constant is observed in the same range of temperatures. No indication is found of a true thermodynamic transition to an ordered low-temperature phase. We therefore suggest that all observed crossovers are dynamic in character. The increase in the dielectric constant is related to the dynamic freezing of dipolar domains on the time scale of simulations.

Quasielastic neutron scattering in biology: Theory and applications

Neutrons scatter quasielastically from stochastic, diffusive processes, such as overdamped vibrations, localized diffusion and transitions between energy minima. In biological systems, such as proteins and membranes, these relaxation processes are of considerable physical interest. We review here recent methodological advances and applications of quasielastic neutron scattering (QENS) in biology, concentrating on the role of molecular dynamics simulation in generating data with which neutron profiles can be unambiguously interpreted. We examine the use of massively-parallel computers in calculating scattering functions, and the application of Markov state modeling. The decomposition of MD-derived neutron dynamic susceptibilities is described, and the use of this in combination with NMR spectroscopy. We discuss dynamics at very long times, including approximations to the infinite time mean-square displacement and nonequilibrium aspects of single-protein dynamics. Finally, we examine how neutron scattering and MD can be combined to provide information on lipid nanodomains. This article is part of a Special Issue entitled "Science for Life" Guest Editor: Dr. Austen Angell, Dr. Salvatore Magazù and Dr. Federica Migliardo.

Contrasting dynamics of fragile and non-fragile polyalcohols through the glass, and dynamical, transitions: A comparison of neutron scattering and dielectric relaxation data for sorbitol and glycerol

BACKGROUND

Glycerol and sorbitol are glass-forming hydrogen-bonded systems characterized by intriguing properties which make these systems very interesting also from the applications point of view. The goal of this work is to relate the hydrogen-bonded features, relaxation dynamics, glass transition properties and fragility of these systems, in particular to seek insight into their very different liquid fragilities.

METHODS

The comparison between glycerol and sorbitol is carried out by collecting the elastic incoherent neutron scattering (EINS) intensity as a function of temperature and of the instrumental energy resolution.

RESULTS

Intensity data vs temperature and resolution are analyzed in terms of thermal restraint and Resolution Elastic Neutron Scattering (RENS) approaches.

CONCLUSIONS

The number of OH groups, which are related to the connecting sites, is a significant parameter both in the glass transition and in the dynamical transition. On the other hand, the disordered nature of sorbitol is confirmed by the existence of different relaxation processes.

GENERAL SIGNIFICANCE

From the applications point of view, glycerol and sorbitol have remarkable bioprotectant properties which make these systems useful in different technological and industrial fields. Furthermore, polyols are rich in glassforming liquid phenomenology and highly deserving of study in their own right. The comparison of EINS and calorimetric data on glycerol and sorbitol helps provide a connection between structural relaxation, dynamical transition, glass transition, and fragility. The evaluation of the inflection point in the elastic intensity behavior as a function of temperature and instrumental energy resolution provides a confirmation of the validity of the RENS approach. This article is part of a Special Issue entitled "Science for Life" Guest Editor: Dr. Austen Angell, Dr. Salvatore Magazù and Dr. Federica Migliardo.

Broadband Dielectric Spectroscopy on Lysozyme in the Sub-Gigahertz to Terahertz Frequency Regions: Effects of Hydration and Thermal Excitation

We have performed dielectric spectral measurements of lysozyme in a solid state to understand the effects of hydration and thermal excitation on the low-frequency dynamics of protein. Dielectric measurements were performed under changing hydration conditions at room temperature in the frequency region of 0.5 GHz to 1.8 THz. We also studied the temperature dependence (83 to 293 K) of the complex dielectric spectra in the THz frequency region (0.3 THz to 1.8 THz). Spectral analyses were performed using model functions for the complex dielectric constant. To reproduce the spectra, we found that two relaxational modes and two underdamped modes are necessary together with an ionic conductivity term in the model function. At room temperature, the two relaxational modes have relaxation times of ∼20 ps and ∼100 ps. The faster component has a major spectral intensity and is suggested to be due to coupled water-protein motion. The two underdamped modes are necessary to reproduce the temperature dependence of the spectra in the THz region satisfactorily. The protein dynamical transition is a well-known behavior in the neutron-scattering experiment for proteins, where the atomic mean-square displacement shows a sudden change in the temperature dependence at approximately 200 K, when the samples are hydrated. A similar behavior has also been observed in the temperature dependence of the absorption spectra of protein in the THz frequency region. From our broadband dielectric spectroscopic measurements, we conclude that the increase in the spectral intensities in the THz region at approximately 200 K is due to a spectral blue-shift of the fast relaxational mode.

Dynamics of hydrated proteins and bio-protectants: Caged dynamics, β-relaxation, and α-relaxation

BACKGROUND

The properties of the three dynamic processes, α-relaxation, ν-relaxation, and caged dynamics in aqueous mixtures and hydrated proteins are analogous to corresponding processes found in van der Waals and polymeric glass-formers apart from minor differences.

METHODS

Collection of various experimental data enables us to characterize the structural α-relaxation of the protein coupled to hydration water (HW), the secondary or ν-relaxation of HW, and the caged HW process.

RESULTS

From the T-dependence of the ν-relaxation time of hydrated myoglobin, lysozyme, and bovine serum albumin, we obtain Ton at which it enters the experimental time windows of Mössbauer and neutron scattering spectroscopies, coinciding with protein dynamical transition (PDT) temperature Td. However, for all systems considered, the α-relaxation time at Ton or Td is many orders of magnitude longer. The other step change of the mean-square-displacement (MSD) at Tg_alpha originates from the coupling of the nearly constant loss (NCL) of caged HW to density. The coupling of the NCL to density is further demonstrated by another step change at the secondary glass temperature Tg_beta in two bio-protectants, trehalose and sucrose.

CONCLUSIONS

The structural α-relaxation plays no role in PDT. Since PDT is simply due to the ν-relaxation of HW, the term PDT is a misnomer. NCL of caged dynamics is coupled to density and show transitions at lower temperature, Tg_beta and Tg_alpha.

GENERAL SIGNIFICANCE

The so-called protein dynamical transition (PDT) of hydrated proteins is not caused by the structural α-relaxation of the protein but by the secondary ν-relaxation of hydration water. "This article is part of a Special Issue entitled "Science for Life" Guest Editor: Dr. Austen Angell, Dr. Salvatore Magazù and Dr. Federica Migliardo".

A neutron spectrometer concept implementing RENS for studies in life sciences

BACKGROUND

Resolution Elastic Neutron Scattering (RENS) method involves performing elastic scattering intensity scans as a function of the instrumental energy resolution and as a function of temperature.

METHODS

In the framework of RENS, numerical simulation and experimental data show that in the measured elastic scattering law against the logarithm of the instrumental energy resolution an inflection point occurs when the resolution time intersects the system relaxation time; conversely, in the measured elastic scattering law against temperature an inflection point turns up when the system relaxation time intersects the resolution time.

RESULTS

For practical implementation of the RENS technique, a dedicated neutron spectrometer would be needed. Here we propose a concept of such a spectrometer that utilizes mechanical velocity selection of both incident and scattered neutrons over a wide angular range. The instrument is able to collect intensity scans vs energy resolution where the instrumental resolution time changes crisscrossing the system relaxation time, and intensity scans vs temperature where the system relaxation time changes intersecting the instrumental resolution time.

CONCLUSIONS

We propose a RENS spectrometer concept based on velocity selection of incident neutrons and wide-angle velocity selection of scattered neutrons achieved by the same rotating collimator-type mechanical device with the optimized shape of blades.

GENERAL SIGNIFICANCE

RENS spectrometer is strongly appealing and innovative because of the simultaneous data collection as a function of energy resolution, wide wavevector range and temperature. Such a spectrometer would be the first practical implementation of RENS concept with a broad range of applications in Life Sciences. This article is part of a Special Issue entitled "Science for Life" Guest Editor: Dr. Austen Angell, Dr. Salvatore Magazù and Dr. Federica Migliardo.

Dynamics of Hydrophobic Core Phenylalanine Residues Probed by Solid-State Deuteron NMR

We conducted a detailed investigation of the dynamics of two phenylalanine side chains in the hydrophobic core of the villin headpiece subdomain protein (HP36) in the hydrated powder state over the 298-80 K temperature range. Our main tools were static deuteron NMR measurements of longitudinal relaxation and line shapes supplemented with computational modeling. The temperature dependence of the relaxation times reveals the presence of two main mechanisms that can be attributed to the ring-flips, dominating at high temperatures, and small-angle fluctuations, dominating at low temperatures. The relaxation is nonexponential at all temperatures with the extent of nonexponentiality increasing from higher to lower temperatures. This behavior suggests a distribution of conformers with unique values of activation energies. The central values of the activation energies for the ring-flipping motions are among the smallest reported for aromatic residues in peptides and proteins and point to a very mobile hydrophobic core. The analysis of the widths of the distributions, in combination with the earlier results on the dynamics of flanking methyl groups (Vugmeyster et al. J. Phys. Chem. B 2013, 117, 6129-6137), suggests that the hydrophobic core undergoes slow concerted fluctuations. There is a pronounced effect of dehydration on the ring-flipping motions, which shifts the distribution toward more rigid conformers. The crossover temperature between the regions of dominance of the small-angle fluctuations and ring-flips shifts from 195 K in the hydrated protein to 278 K in the dry one. This result points to the role of solvent in softening the core and highlights aromatic residues as markers of the protein dynamical transitions.

Nanometer-sized dynamic entities in an aqueous system

Using neutron spin-echo and backscattering spectroscopy, we have found that at low temperatures water molecules in an aqueous solution engage in center-of-mass dynamics that are different from both the main structural relaxations and the well-known localized motions in the transient cages of the nearest neighbor molecules. While the latter localized motions are known to take place on the picosecond time scale and Angstrom length scale, the slower motions that we have observed are found on the nanosecond time scale and nanometer length scale. They are associated with the slow secondary relaxations, or excess wing dynamics, in glass-forming liquids. Our approach, therefore, can be applied to probe the characteristic length scale of the dynamic entities associated with slow dynamics in glass-forming liquids, which presently cannot be studied by other experimental techniques.

Hydration-dependent dynamic crossover phenomenon in protein hydration water

The characteristic relaxation time τ of protein hydration water exhibits a strong hydration level h dependence. The dynamic crossover is observed when h is higher than the monolayer hydration level hc=0.2-0.25 and becomes more visible as h increases. When h is lower than hc, τ only exhibits Arrhenius behavior in the measured temperature range. The activation energy of the Arrhenius behavior is insensitive to h, indicating a local-like motion. Moreover, the h dependence of the crossover temperature shows that the protein dynamic transition is not directly or solely induced by the dynamic crossover in the hydration water.

Temperature-dependent dynamics of dry and hydrated β-casein studied by quasielastic neutron scattering

β-Casein is a component of casein micelle with amphillic nature and is recognized as a "natively disordered" protein that lacks secondary structures. In this study, the temperature and hydration effects on the dynamics of β-casein are explored by quasielastic neutron scattering (QENS). An upturn in the mean square displacement (MSD) of hydrated β-casein indicates an increase of protein flexibility at a temperature of ~225 K. Another increase in MSD at ~100 K, observed in both dry and hydrated β-casein, is ascribed to the methyl group rotations, which are not sensitive to hydration. QENS analysis in the energy domain reveals that the fraction of hydrogen atoms participating in motion in a sphere of diffusion is highly hydration dependent and increases with temperature. In the time domain analysis, a logarithmic-like decay is observed in the range of picosecond to nanosecond (β-relaxation time) in the dynamics of hydrated β-casein. This dynamical behavior has been observed in hydrated globular and oligomeric proteins. Our temperature-dependent QENS experiments provide evidence that lack of a secondary structure in β-casein results in higher flexibility in its dynamics and easier reversible thermal unfolding compared to other rigid biomolecules.

The importance of individual protein molecule dynamics in developing and assessing solid state protein preparations

Processing protein solutions into the solid state is a common approach for generating stable amorphous protein mixtures that are suitable for long-term storage. Great care is typically given to protecting the protein native structure during the various drying steps that render it into the amorphous solid state. However, many studies illustrate that chemical and physical degradations still occur in spite of this amorphous material having good glassy properties and it being stored at temperatures below its glass transition temperature (Tg). Because of these persistent issues and recent biophysical studies that have refined the debate ascribing meaning to the molecular dynamical transition temperature and Tg of protein molecules, we provide an updated discussion on the impact of assessing and managing localized, individual protein molecule nondiffusive motions in the context of proteins being prepared into bulk amorphous mixtures. Our aim is to bridge the pharmaceutical studies addressing bulk amorphous preparations and their glassy behavior, with the biophysical studies historically focused on the nondiffusive internal protein dynamics and a protein's activity, along with their combined efforts in assessing the impact of solvent hydrogen-bonding networks on local stability. We also provide recommendations for future research efforts in solid-state formulation approaches.

Communication: Protein dynamical transition vs. liquid-liquid phase transition in protein hydration water

In this work, we compare experimental data on myoglobin hydrated powders from elastic neutron scattering, broadband dielectric spectroscopy, and differential scanning calorimetry. Our aim is to obtain new insights on the connection between the protein dynamical transition, a fundamental phenomenon observed in proteins whose physical origin is highly debated, and the liquid-liquid phase transition (LLPT) possibly occurring in protein hydration water and related to the existence of a low temperature critical point in supercooled water. Our results provide a consistent thermodynamic/dynamic description which gives experimental support to the LLPT hypothesis and further reveals how fundamental properties of water and proteins are tightly related.

Does a dry protein undergo a glass transition?

Bovine serum albumin (BSA) with extremely low hydration level 0.04, which is usually defined as dry, has been investigated in the temperature range between 200 and 340 K by incoherent inelastic neutron scattering using the neutron time-of-flight spectrometer FOCUS (PSI, Switzerland). Anomalous temperature behavior has been revealed for relaxational and low-frequency vibrational dynamics of BSA in the vicinity of 250 K. The mean-square atomic displacement has been shown to exhibit a change in the slope of temperature dependence near the same temperature. The presented results point out that the glass-like transition occurs in the dry protein.

Vanishing amplitude of backbone dynamics causes a true protein dynamical transition: 2H NMR studies on perdeuterated C-phycocyanin

Using a combination of H2 nuclear magnetic resonance (NMR) methods, we study internal rotational dynamics of the perdeuterated protein C-phycocyanin (CPC) in dry and hydrated states over broad temperature and dynamic ranges with high angular resolution. Separating H2 NMR signals from methyl deuterons, we show that basically all backbone deuterons exhibit highly restricted motion occurring on time scales faster than microseconds. The amplitude of this motion increases when a hydration shell exists, while it decreases upon cooling and vanishes near 175 K. We conclude that the vanishing of the highly restricted motion marks a dynamical transition, which is independent of the time window and of a fundamental importance. This conclusion is supported by results from experimental and computational studies of the proteins myoglobin and elastin. In particular, we argue based on findings in molecular dynamics simulations that the behavior of the highly restricted motion of proteins at the dynamical transition resembles that of a characteristic secondary relaxation of liquids at the glass transition, namely the nearly constant loss. Furthermore, H2 NMR studies on perdeuterated CPC reveal that, in addition to highly restricted motion, small fractions of backbone segments exhibit weakly restricted dynamics when temperature and hydration are sufficiently high.

Scaling analysis of bio-molecular dynamics derived from elastic incoherent neutron scattering experiments

Numerous neutron scattering studies of bio-molecular dynamics employ a qualitative analysis of elastic scattering data and atomic mean square displacements. We provide a new quantitative approach showing that the intensity at zero energy exchange can be a rich source of information of bio-structural fluctuations on a pico- to nano-second time scale. Elastic intensity scans performed either as a function of the temperature (back-scattering) and∕or by varying the instrumental resolution (time of flight spectroscopy) yield the activation parameters of molecular motions and the approximate structural correlation function in the time domain. The two methods are unified by a scaling function, which depends on the ratio of correlation time and instrumental resolution time. The elastic scattering concept is illustrated with a dynamic characterization of alanine-dipeptide, protein hydration water, and water-coupled protein motions of lysozyme, per-deuterated c-phycocyanin (CPC) and hydrated myoglobin. The complete elastic scattering function versus temperature, momentum exchange, and instrumental resolution is analyzed instead of focusing on a single cross-over temperature of mean square displacements at the apparent onset temperature of an-harmonic motions. Our method predicts the protein dynamical transition (PDT) at Td from the collective (α) structural relaxation rates of the solvation shell as input. By contrast, the secondary (β) relaxation enhances the amplitude of fast local motions in the vicinity of the glass temperature Tg. The PDT is specified by step function in the elastic intensity leading from elastic to viscoelastic dynamic behavior at a transition temperature Td.

Dynamics of supercooled water in a biological model system of the amino acid l-lysine

Lysine solutions establish a new relaxation behaviour of supercooled interfacial water.

Hydration dependence of myoglobin dynamics studied with elastic neutron scattering, differential scanning calorimetry and broadband dielectric spectroscopy

In this work we present a thorough investigation of the hydration dependence of myoglobin dynamics. The study is performed on D2O-hydrated protein powders in the hydration range 0<h<0.5 (h≡gr[D2O]/gr[protein]) and in the temperature range 20-300K. The protein equilibrium fluctuations are investigated with Elastic Neutron Scattering using the spectrometer IN13 at ILL (Grenoble), while the relaxations of the protein+hydration water system are investigated with Broadband Dielectric Spectroscopy; finally, Differential Scanning Calorimetry is used to obtain a thermodynamic description of the system. The effect of increasing hydration is to speed up the relaxations of the myoglobin+hydration water system and, thermodynamically, to decrease the glass transition temperature; these effects tend to saturate at h values greater than ~0.3. Moreover, the calorimetric scans put in evidence the occurrence of an endothermic peak whose onset temperature is located at ~230K independent of hydration. From the point of view of the protein equilibrium fluctuations, while the amplitude of anharmonic mean square displacements is found to increase with hydration, their onset temperature (i.e. the onset temperature of the well known "protein dynamical transition") is hydration independent. On the basis of the above results, the relevance of protein+hydration water relaxations and of the thermodynamic state of hydration water to the onset of the protein dynamical transition is discussed.