Physical origin of anharmonic dynamics in proteins: new insights from resolution-dependent neutron scattering on homomeric polypeptides

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.

A Long Journey into the Investigation of the Structure–Dynamics–Function Paradigm in Proteins through the Activities of the Palermo Biophysics Group

An overview of the biophysics activity at the Department of Physics and Chemistry Emilio Segrè of the University of Palermo is given. For forty years, the focus of the research has been on the protein structure–dynamics–function paradigm, with the aim of understanding the molecular basis of the relevant mechanisms and the key role of solvent. At least three research lines are identified; the main results obtained in collaboration with other groups in Italy and abroad are presented. This review is dedicated to the memory of Professors Massimo Ugo Palma, Maria Beatrice Palma Vittorelli, and Lorenzo Cordone, which were the founders of the Palermo School of Biophysics. We all have been, directly or indirectly, their pupils; we miss their enthusiasm for scientific research, their deep physical insights, their suggestions, their strict but always constructive criticisms, and, most of all, their friendship. This paper is dedicated also to the memory of Prof. Hans Frauenfelder, whose pioneering works on nonexponential rebinding kinetics, protein substates, and energy landscape have inspired a large part of our work in the field of protein dynamics.

Temperature-Dependent Dynamics at Protein-Solvent Interfaces

We perform differential scanning calorimetry, broadband dielectric spectroscopy (BDS), and nuclear magnetic resonance (NMR) studies to ascertain the molecular dynamics in mixtures of ethylene glycol with elastin or lysozyme over broad temperature ranges. To focus on the protein-solvent interface, we use mixtures with about equal numbers of amino acids and solvent molecules. The elastin and lysozyme mixtures show similar glass transition steps, which extend over a broad temperature range of 157-185K. The BDS and NMR studies yield fully consistent results for the fastest process P1, which is caused by the structural relaxation of ethylene glycol between the protein molecules and follows an Arrhenius law with an activation energy of Ea=0.63eV. It involves quasi-isotropic reorientation and is very similar in the elastin and lysozyme matrices but different from the alpha and beta relaxations of bulk ethylene glycol. Two slower BDS processes P2 and P3 have protein-dependent time scales, but exhibit a similar Arrhenius-like temperature dependence with an activation energy of Ea~0.81eV. However, P2 and P3 do not have a clear NMR signature. In particular, the NMR results for the lysozyme mixture reveal that the protein backbone does not show isotropic alpha-like motion on the P2 and P3 time scales but only restricted beta-like reorientation. The different activation energies of the P1 and P2/P3 processes do not support an intimate coupling of protein and ethylene glycol dynamics. The present results are compared with previous findings for mixtures of proteins with water or glycerol, implying qualitatively different dynamical couplings at various protein-solvent interfaces.

Universal dynamical onset in water at distinct material interfaces

Interfacial water remains liquid and mobile much below 0 °C, imparting flexibility to the encapsulated materials to ensure their diverse functions at subzero temperatures. However, a united picture that can describe the dynamical differences of interfacial water on different materials and its role in imparting system-specific flexibility to distinct materials is lacking. By combining neutron spectroscopy and isotope labeling, we explored the dynamics of water and the underlying substrates independently below 0 °C across a broad range of materials. Surprisingly, while the function-related anharmonic dynamical onset in the materials exhibits diverse activation temperatures, the surface water presents a universal onset at a common temperature. Further analysis of the neutron experiment and simulation results revealed that the universal onset of water results from an intrinsic surface-independent relaxation: switching of hydrogen bonds between neighboring water molecules with a common energy barrier of ∼35 kJ mol⁻¹.

We demonstrated that the dynamical onset of interfacial water is an intrinsic property of water itself, resulting from a surface independent relaxation process in water with an approximately universal energy barrier of ∼35 kJ mol⁻¹.

Diffusivelike Motions in a Solvent-Free Protein-Polymer Hybrid

The interaction between proteins and hydration water stabilizes protein structure and promotes functional dynamics, with water translational motions enabling protein flexibility. Engineered solvent-free protein-polymer hybrids have been shown to preserve protein structure, function, and dynamics. Here, we used neutron scattering, protein and polymer perdeuteration, and molecular dynamics simulations to explore how a polymer dynamically replaces water. Even though relaxation rates and vibrational properties are strongly modified in polymer coated compared to hydrated proteins, liquidlike polymer dynamics appear to plasticize the conjugated protein in a qualitatively similar way as do hydration-water translational motions.

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.

More than a Confinement: “Soft” and “Hard” Enzyme Entrapment Modulates Biological Catalyst Function

Catalysis makes chemical and biochemical reactions kinetically accessible. From a technological point of view, organic, inorganic, and biochemical catalysis is relevant for several applications, from industrial synthesis to biomedical, material, and food sciences. A heterogeneous catalyst, i.e., a catalyst confined in a different phase with respect to the reagents’ phase, requires either its physical confinement in an immobilization matrix or its physical adsorption on a surface. In this review, we will focus on the immobilization of biological catalysts, i.e., enzymes, by comparing hard and soft immobilization matrices and their effect on the modulation of the catalysts’ function. Indeed, unlike smaller molecules, the catalytic activity of protein catalysts depends on their structure, conformation, local environment, and dynamics, properties that can be strongly affected by the immobilization matrices, which, therefore, not only provide physical confinement, but also modulate catalysis.

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.

A Quasielastic Neutron Scattering Investigation on the Molecular Self-Dynamics of Human Myelin Protein P2

The human myelin protein P2 is a membrane binding protein believed to maintain correct lipid composition and organization in peripheral nerve myelin. Its function is related to its ability to stack membranes, and this function can be enhanced by the P38G mutation, whereby the overall protein structure does not change but the molecular dynamics increase. Mutations in P2 are linked to human peripheral neuropathy. Here, the dynamics of wild-type P2 and the P38G variant were studied using quasielastic neutron scattering on time scales from 10 ps to 1 ns at 300 K. The results suggest that the mutant protein dynamics are increased on both the fastest and the slowest measured time scales, by increasing the dynamics amplitude and/or the portion of atoms participating in the movement.

Disentangling time-focusing from beam divergence: A novel approach for high-flux thermal neutron spectroscopy at continuous and long-pulse sources

We present the concept of a novel time-focusing technique for neutron spectrometers, which allows us to disentangle time-focusing from beam divergence. The core of this approach is a double rotating-crystal monochromator that can be used to extract a larger wavelength band from a white beam, thus providing a higher flux at the sample compared to standard time-of-flight instruments, yet preserving energy resolution and beam collimation. The performances of a spectrometer based on this approach are quantitatively discussed in terms of possible incident wavelengths, flux at the sample, and (Q, E)-resolution. Analytical estimates suggest flux gains of about one order of magnitude at comparable resolutions in comparison to conventional time-of-flight spectrometers. Moreover, the double monochromator configuration natively shifts the sample away from the source line-of-sight, thus significantly improving the signal-to-noise ratio. The latter, in combination with a system that does not increase the beam divergence, brings the further advantage of a cleaner access to the low-Q region, which is recognized to be of fundamental interest for magnetism and for disordered materials, from glasses to biological systems.

Unsaturated lipid bilayers at cryogenic temperature: librational dynamics of chain-labeled lipids from pulsed and CW-EPR

Fully hydrated bilayers of monounsaturated palmitoyloleoylphosphatidylcholine (POPC) and diunsaturated dioleoylphosphatidylcholine (DOPC) lipids have low main phase transition temperatures (271 K for POPC and 253 K for DOPC). Two-pulse echo detected spectra, combined with continuous wave electron paramagnetic resonance spectroscopy, are employed to study the low-temperature lamellar phases of the POPC and DOPC unsaturated bilayers that are usually studied in the fluid state. Phosphatidylcholine spin-labeled at C-5 and C-16 carbon atom positions along the acyl chain were used and the temperature varied over the range 77-270 K. Segmental chain librational oscillations of small amplitude and with correlation time in the subnanosecond to nanosecond range are found in both membranes. The mean-square angular amplitude, α2, of librations increases with temperature, is larger close to the bilayer midplane than close to the first acyl chain segments, and is larger in diunsaturated than in monounsaturated bilayers. In the inner hydrocarbon region of both lipid matrices, α2 increases first slowly and linearly with temperature and then more rapidly, and a dynamical transition is detected in the range 190-210 K. Compared to dipalmitoylphosphatidylcholine bilayers of fully saturated symmetric chain lipids, the presence of double bonds in the acyl chain enhances the intensity of librational motion which is characterized by larger angular variations at the terminal methyl ends. These findings highlight biophysical properties of unsaturated bilayers in the frozen state, including a detailed characterization of segmental chain dynamics and the evidence of a dynamical transition that appears to be a generic feature in hydrated macromolecular systems. These results can also be relevant in regulating membrane physical properties and function at higher physiological temperatures.

The Role of Hydrogen Bonding in the Folding/Unfolding Process of Hydrated Lysozyme: A Review of Recent NMR and FTIR Results

The biological activity of proteins depends on their three-dimensional structure, known as the native state. The main force driving the correct folding mechanism is the hydrophobic effect and when this folding kinetics is altered, aggregation phenomena intervene causing the occurrence of illnesses such as Alzheimer and Parkinson’s diseases. The other important effect is performed by water molecules and by their ability to form a complex network of hydrogen bonds whose dynamics influence the mobility of protein amino acids. In this work, we review the recent results obtained by means of spectroscopic techniques, such as Fourier Transform Infrared (FTIR) and Nuclear Magnetic Resonance (NMR) spectroscopies, on hydrated lysozyme. In particular, we explore the Energy Landscape from the thermal region of configurational stability up to that of the irreversible denaturation. The importance of the coupling between the solute and the solvent will be highlighted as well as the different behaviors of hydrophilic and hydrophobic moieties of protein amino acid residues.

Direct Experimental Characterization of Contributions from Self-Motion of Hydrogen and from Interatomic Motion of Heavy Atoms to Protein Anharmonicity

One fundamental challenge in biophysics is to understand the connection between protein dynamics and its function. Part of the difficulty arises from the fact that proteins often present local atomic motions and collective dynamics on the same time scales, and challenge the experimental identification and quantification of different dynamic modes. Here, by taking lyophilized proteins as the example, we combined deuteration technique and neutron scattering to separate and characterize the self-motion of hydrogen and the collective interatomic motion of heavy atoms (C, O, N) in proteins on the pico-to-nanosecond time scales. We found that hydrogen atoms present an instrument-resolution-dependent onset for anharmonic motions, which can be ascribed to the thermal activation of local side-group motions. However, the protein heavy atoms exhibit an instrument-resolution-independent anharmonicity around 200 K, which results from unfreezing of the relaxation of the protein structures on the laboratory equilibrium time (100-1000 s), softening of the entire bio-macromolecules.

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.

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.

Low-Temperature Dynamics of Chain-Labeled Lipids in Ester- and Ether-Linked Phosphatidylcholine Membranes

Continuous wave electron paramagnetic resonance spectroscopy and two-pulse echo detected spectra of chain-labeled lipids are used to study the dynamics of frozen lipid membranes over the temperature range 77-260 K. Bilayers of ester-linked dihexadecanoylphosphatidylcholine (DPPC) with noninterdigitated chains and ether-linked dihexadecyl phosphatidylcholine (DHPC) with interdigitated chains are considered. Rapid stochastic librations of small angular amplitude are found in both lipid matrices. In noninterdigitated DPPC bilayers, the mean-square angular amplitude, [Formula: see text], of the motion increases with temperature and it is larger close to the chain termini than close to the polar/apolar interface. In contrast, in interdigitated DHPC lamellae, [Formula: see text] is small and temperature and label-position independent at low temperature and increases steeply at high temperature. The rotational correlation time, τ_c, of librations lies in the subnanosecond range for DPPC and in the nanosecond range for DHPC. In all membrane samples, the temperature dependence of [Formula: see text] resembles that of the mean-square atomic displacement revealed by neutron scattering and a dynamical transition is detected in the range 210-240 K. The results highlight the librational oscillations and the glass-like behavior in bilayer and interdigitated lipid membranes.

Critical structural fluctuations of proteins upon thermal unfolding challenge the Lindemann criterion

Internal subnanosecond timescale motions are key for the function of proteins, and are coupled to the surrounding solvent environment. These fast fluctuations guide protein conformational changes, yet their role for protein stability, and for unfolding, remains elusive. Here, in analogy with the Lindemann criterion for the melting of solids, we demonstrate a common scaling of structural fluctuations of lysozyme protein embedded in different environments as the thermal unfolding transition is approached. By combining elastic incoherent neutron scattering and advanced molecular simulations, we show that, although different solvents modify the protein melting temperature, a unique dynamical regime is attained in proximity of thermal unfolding in all solvents that we tested. This solvation shell-independent dynamical regime arises from an equivalent sampling of the energy landscape at the respective melting temperatures. Thus, we propose that a threshold for the conformational entropy provided by structural fluctuations of proteins exists, beyond which thermal unfolding is triggered.

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.

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".

Molecular Dynamics Simulations of a Powder Model of the Intrinsically Disordered Protein Tau

The tau protein, whose aggregates are involved in Alzheimer's disease, is an intrinsically disordered protein (IDP) that regulates microtubule activity in neurons. An IDP lacks a single, well-defined structure and, rather, constantly exchanges among multiple conformations. In order to study IDP dynamics, the combination of experimental techniques, such as neutron scattering, and computational techniques, such as molecular dynamics (MD) simulations, is a powerful approach. Amorphous hydrated powder samples have been very useful for studying protein internal dynamics experimentally, e.g., using neutron scattering. Thus, there is demand for realistic in silico models of hydrated protein powders. Here we present an MD simulation analysis of a powder hydrated at 0.4 g water/g protein of the IDP tau in the temperature range 20-300 K. By comparing with neutron scattering data, we identify the protein-water interface as the predominant feature determining IDP dynamics. The so-called protein dynamical transition is shown to be attenuated, but not suppressed, in the parts of the protein that are not exposed to the solvent. In addition, we find similarities in the mean-squared displacements of the core of a globular protein and "dry" clusters formed by the IDP in hydrated powders. Thus, the ps to ns dynamics of proteins in hydrated powders originate mainly from those residues in contact with solvent. We propose that by measuring the dynamics of protein assemblies, such as aggregates, one might assess qualitatively their state of hydration.

Dynamics of the Peripheral Membrane Protein P2 from Human Myelin Measured by Neutron Scattering--A Comparison between Wild-Type Protein and a Hinge Mutant

Myelin protein P2 is a fatty acid-binding structural component of the myelin sheath in the peripheral nervous system, and its function is related to its membrane binding capacity. Here, the link between P2 protein dynamics and structure and function was studied using elastic incoherent neutron scattering (EINS). The P38G mutation, at the hinge between the β barrel and the α-helical lid, increased the lipid stacking capacity of human P2 in vitro, and the mutated protein was also functional in cultured cells. The P38G mutation did not change the overall structure of the protein. For a deeper insight into P2 structure-function relationships, information on protein dynamics in the 10 ps to 1 ns time scale was obtained using EINS. Values of mean square displacements mainly from protein H atoms were extracted for wild-type P2 and the P38G mutant and compared. Our results show that at physiological temperatures, the P38G mutant is more dynamic than the wild-type P2 protein, especially on a slow 1-ns time scale. Molecular dynamics simulations confirmed the enhanced dynamics of the mutant variant, especially within the portal region in the presence of bound fatty acid. The increased softness of the hinge mutant of human myelin P2 protein is likely related to an enhanced flexibility of the portal region of this fatty acid-binding protein, as well as to its interactions with the lipid bilayer surface requiring conformational adaptations.

Thermodynamic properties of bulk and confined water

The thermodynamic response functions of water display anomalous behaviors. We study these anomalous behaviors in bulk and confined water. We use nuclear magnetic resonance (NMR) to examine the configurational specific heat and the transport parameters in both the thermal stable and the metastable supercooled phases. The data we obtain suggest that there is a behavior common to both phases: that the dynamics of water exhibit two singular temperatures belonging to the supercooled and the stable phase, respectively. One is the dynamic fragile-to-strong crossover temperature (T(L) ≃ 225 K). The second, T* ∼ 315 ± 5 K, is a special locus of the isothermal compressibility K(T)(T, P) and the thermal expansion coefficient α(P)(T, P) in the P-T plane. In the case of water confined inside a protein, we observe that these two temperatures mark, respectively, the onset of protein flexibility from its low temperature glass state (T(L)) and the onset of the unfolding process (T*).

Coupling between overall rotational diffusion and domain motions in proteins and its effect on dielectric spectra

In this work, we formulate a closed-form solution of the model of a semirigid molecule for the case of fluctuating and reorienting molecular electric dipole moment. We illustrate with numeric calculations the impact of protein domain motions on dielectric spectra using the example of the 128 kDa protein dimer of Enzyme I. We demonstrate that the most drastic effect occurs for situations when the characteristic time of protein domain dynamics is comparable to the time of overall molecular rotational diffusion. We suggest that protein domain motions could be a possible explanation for the high-frequency contribution that accompanies the major relaxation dispersion peak in the dielectric spectra of protein aqueous solutions. We propose that the presented computational methodology could be used for the simultaneous analysis of dielectric spectroscopy and nuclear magnetic resonance data. Proteins 2015; 83:1571-1581. © 2015 Wiley Periodicals, Inc.

Intramolecular vibrations in low-frequency normal modes of amino acids: L-alanine in the neat solid state

This paper presents a theoretical analysis of the low-frequency phonons of L-alanine by using the solid-state density functional theory at the Γ point. We are particularly interested in the intramolecular vibrations accessing low-frequency phonons via harmonic coupling with intermolecular vibrations. A new mode-analysis method is introduced to quantify the vibrational characteristics of such intramolecular vibrations. We find that the torsional motions of COO(-) are involved in low-frequency phonons, although COO(-) is conventionally assumed to undergo localized torsion. We also find the broad distributions of intramolecular vibrations relevant to important functional groups of amino acids, e.g., the COO(-) and NH3(+) torsions, in the low-frequency phonons. The latter finding is illustrated by the concept of frequency distribution of vibrations. These findings may lead to immediate implications in other amino acid systems.

Translational diffusion of hydration water correlates with functional motions in folded and intrinsically disordered proteins

Hydration water is the natural matrix of biological macromolecules and is essential for their activity in cells. The coupling between water and protein dynamics has been intensively studied, yet it remains controversial. Here we combine protein perdeuteration, neutron scattering and molecular dynamics simulations to explore the nature of hydration water motions at temperatures between 200 and 300 K, across the so-called protein dynamical transition, in the intrinsically disordered human protein tau and the globular maltose binding protein. Quasi-elastic broadening is fitted with a model of translating, rotating and immobile water molecules. In both experiment and simulation, the translational component markedly increases at the protein dynamical transition (around 240 K), regardless of whether the protein is intrinsically disordered or folded. Thus, we generalize the notion that the translational diffusion of water molecules on a protein surface promotes the large-amplitude motions of proteins that are required for their biological activity.

Hydration water plasticizes protein structures and is essential for their biological functions, such as enzymatic catalysis. Here, the authors use neutron scattering and molecular dynamics simulations to study hydration water at the dynamical transition of folded and disordered proteins.

Effect of temperature on the low-frequency vibrational spectrum and relative structuring of hydration water around a single-stranded DNA

Molecular dynamics simulations of the single-stranded DNA oligomer (5'-CGCGAAT TCGCG-3') in aqueous solution have been carried out at different temperatures between 160 K and 300 K. The effects of temperature on the low-frequency vibrational spectrum and local structural arrangements of water molecules hydrating the DNA strand have been explored in detail. The low-frequency density of states distributions reveal that increasingly trapped transverse water motions play a dominant role in controlling the band corresponding to O⋯O⋯O bending or transverse oscillations of hydration water at supercooled temperatures. In addition, presence of a broad band around 260 (±20) cm(-1) under supercooled conditions indicates transformation from high density liquid-like structuring of hydration water at higher temperatures to that of a low density liquid at lower temperatures. It is found that long-range correlations between the supercooled hydration water molecules arise due to such local structural transition around the DNA oligomer.

Experimental evidence for a liquid-liquid crossover in deeply cooled confined water

In this work we investigate, by means of elastic neutron scattering, the pressure dependence of mean square displacements (MSD) of hydrogen atoms of deeply cooled water confined in the pores of a three-dimensional disordered SiO2 xerogel; experiments have been performed at 250 and 210 K from atmospheric pressure to 1200 bar. The "pressure anomaly" of supercooled water (i.e., a mean square displacement increase with increasing pressure) is observed in our sample at both temperatures; however, contrary to previous simulation results and to the experimental trend observed in bulk water, the pressure effect is smaller at lower (210 K) than at higher (250 K) temperature. Elastic neutron scattering results are complemented by differential scanning calorimetry data that put in evidence, besides the glass transition at about 170 K, a first-order-like endothermic transition occurring at about 230 K that, in view of the neutron scattering results, can be attributed to a liquid-liquid crossover. Our results give experimental evidence for the presence, in deeply cooled confined water, of a crossover occurring at about 230 K (at ambient pressure) from a liquid phase predominant at 210 K to another liquid phase predominant at 250 K; therefore, they are fully consistent with the liquid-liquid transition hypothesis.

The role of water in protein's behavior: The two dynamical crossovers studied by NMR and FTIR techniques

The role the solvent plays in determining the biological activity of proteins is of primary importance. Water is the solvent of life and proteins need at least a water monolayer covering their surface in order to become biologically active. We study how the properties of water and the effect of its coupling with the hydrophilic moieties of proteins govern the regime of protein activity. In particular we follow, by means of Fourier Transform Infrared spectroscopy, the thermal evolution of the amide vibrational modes of hydrated lysozyme in the temperature interval 180 K < T < 350 K. In such a way we are able to observe the thermal limit of biological activity characterizing hydrated lysozyme. Finally we focus on the region of lysozyme thermal denaturation by following the evolution of the proton Nuclear Magnetic Resonance (NMR) spectra for 298 K < T < 366 K with the High-Resolution Magic Angle Spinning probe. Our data suggest that the hydrogen bond coupling between hydration water and protein hydrophilic groups is crucial in triggering the main mechanisms that define the enzymatic activity of proteins.

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.

Heterogeneity of protein substates visualized by spin-label EPR

The energy landscape of proteins is characterized by a hierarchy of substates, which give rise to conformational heterogeneity at low temperatures. In multiply spin-labeled membranous Na,K-ATPase, this heterogeneous population of conformations is manifest by strong inhomogeneous broadening of the electron paramagnetic resonance (EPR) line shapes and nonexponential spin-echo decays, which undergo a transition to homogeneous broadening and exponential relaxation at higher temperatures (previous study). In this study, we apply these EPR methods to small water-soluble proteins, of the type for which the existence of conformational substates is well established. Both α-helical and β-sheet aqueous proteins that are spin-labeled on a single cysteine residue display spin-echo decays with a single phase-memory time T2M and conventional EPR line shapes with predominantly homogeneous broadening, over a broad range of temperatures from 77 K to ∼ 250 K or higher. Above ∼ 200 K, the residual inhomogeneous broadening is reduced almost to zero. In contrast, both the proteins and the spin label alone, when in a glycerol-water mixture below the glass transition, display heterogeneity in spin-echo phase-memory time and a stronger inhomogeneous broadening of the conventional line shapes, similar to multiply spin-labeled membranous Na,K-ATPase below 200 K. Above 200 K (or the glass transition), a single phase-memory time and predominantly homogeneous broadening are found in both spin-label systems. The results are discussed in terms of solvent-mediated protein transitions, the ability of single spin-label sites to detect conformational heterogeneity, and the desirability of exploring multiple sites for proteins with the size and complexity of the Na,K-ATPase.

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.

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.

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.

Long-time mean-square displacements in proteins

We propose a method for obtaining the intrinsic, long-time mean square displacement (MSD) of atoms and molecules in proteins from finite-time molecular dynamics (MD) simulations. Typical data from simulations are limited to times of 1 to 10 ns, and over this time period the calculated MSD continues to increase without a clear limiting value. The proposed method consists of fitting a model to MD simulation-derived values of the incoherent intermediate neutron scattering function, I(inc)(Q,t), for finite times. The infinite-time MSD, <r(2)>, appears as a parameter in the model and is determined by fits of the model to the finite-time I(inc)(Q,t). Specifically, the <r(2)> is defined in the usual way in terms of the Debye-Waller factor as I(Q,t=∞)=exp(-Q(2)<r(2)>/3). The method is illustrated by obtaining the intrinsic MSD <r(2)> of hydrated lysozyme powder (h=0.4 g water/g protein) over a wide temperature range. The intrinsic <r(2)> obtained from data out to 1 and to 10 ns is found to be the same. The intrinsic <r(2)> is approximately twice the value of the MSD that is reached in simulations after times of 1 ns which correspond to those observed using neutron instruments that have an energy resolution width of 1 μeV.

Protein dynamics by neutron scattering

Neutron scattering techniques represent a powerful tool for characterizing both the structure and dynamical properties of bio-systems, for example, proteins and membranes interacting with their solvents. In this paper, Elastic Neutron Scattering (ENS) data collected at the Institut Laue-Langevin (Grenoble, France) on dry and D2O hydrated lysozyme by varying hydration level are presented, and compared with previously published data on the same protein system, also with the addition of bio-protectants. The data have been collected with three different spectrometers, i.e. IN13, IN10 and IN4. This set of ENS data gives direct access to the temperature behavior of both (i) the Mean Square Displacement (MSD) and (ii) the characteristic system relaxation time. As a result, an explicative hypothesis on the relationship between the so-called "protein dynamical transition" (PDT) and the "fragile-to-strong dynamical crossover" (FSC) is formulated. Furthermore, by taking into proper account the effect of the finite instrumental energy resolution of the used spectrometers, the vibrational MSD of dry and hydrated lysozyme is calculated. The vibrational MSD of the lysozyme in the dry state resulted to be higher than the one in the hydrated state; the latter reaches the former at a temperature value of T=220K that corresponds to the temperature at which the FSC occurs. As a result, a cage effect resulting from the hydration water on the protein surface is hypothesized and subsequently linked to the FSC.