The low-temperature inflection observed in neutron scattering measurements of proteins is due to methyl rotation: direct evidence using isotope labeling and molecular dynamics simulations

Molecular Dynamics of Lysozyme Amyloid Polymorphs Studied by Incoherent Neutron Scattering

Lysozyme amyloidosis is a hereditary disease, which is characterized by the deposition of lysozyme amyloid fibrils in various internal organs. It is known that lysozyme fibrils show polymorphism and that polymorphs formed at near-neutral pH have the ability to promote more monomer binding than those formed at acidic pH, indicating that only specific polymorphs become dominant species in a given environment. This is likely due to the polymorph-specific configurational diffusion. Understanding the possible differences in dynamical behavior between the polymorphs is thus crucial to deepen our knowledge of amyloid polymorphism and eventually elucidate the molecular mechanism of lysozyme amyloidosis. In this study, molecular dynamics at sub-nanosecond timescale of two kinds of polymorphic fibrils of hen egg white lysozyme, which has long been used as a model of human lysozyme, formed at pH 2.7 (LP27) and pH 6.0 (LP60) was investigated using elastic incoherent neutron scattering (EINS) and quasi-elastic neutron scattering (QENS). Analysis of the EINS data showed that whereas the mean square displacement of atomic motions is similar for both LP27 and LP60, LP60 contains a larger fraction of atoms moving with larger amplitudes than LP27, indicating that the dynamical difference between the two polymorphs lies not in the averaged amplitude, but in the distribution of the amplitudes. Furthermore, analysis of the QENS data showed that the jump diffusion coefficient of atoms is larger for LP60, suggesting that the atoms of LP60 undergo faster diffusive motions than those of LP27. This study thus characterizes the dynamics of the two lysozyme polymorphs and reveals that the molecular dynamics of LP60 is enhanced compared with that of LP27. The higher molecular flexibility of the polymorph would permit to adjust its conformation more quickly than its counterpart, facilitating monomer binding.

Dynamical Manifestations of Supercooling in Amyloid Hydration

The effect of extreme temperature on amyloidogenic species remains sparsely explored. In a recent study (J. Phys. Chem. Lett., 2019, 10, (10)), we employed exhaustive molecular dynamics simulations to explore the cold thermal response of a putative small amyloid oligomer and to elicit the role of solvent modulation. Herein, we investigate the dynamical response of the hydration waters of the oligomer within the supercooled states. Using NMR-based formalism, we delineate the entropic response in terms of the side-chain conformational entropy that corroborates the weakening of the hydrophobic core with lowering of temperature. The translational dynamics of the protein and hydration waters reveal the coupling of protein dynamical fluctuations with solvent dynamics under supercooled conditions. Probing the translational motion as a space-time correlation indicates glassy dynamics exhibited by hydration waters in the supercooled regime. Caging of the water molecules with lowering of temperature and the resultant hopping dynamics are reflected in the longer β-relaxation timescales of translational motion. Furthermore, we utilized mode-coupling theory (MCT) and derived the ideal glass transition temperature from translational and rotational dynamics, around ∼196 and 209 K, respectively. Interestingly, rotational motion in the supercooled regime deviates from the MCT law, exhibits Arrhenius motion, and marks a fragile-to-strong crossover at 227 K. The low-frequency vibrational modes also coincide with the dynamical transition. This exposition lends dynamical insights into the hydration coupling of an amyloid aggregate under cryogenic conditions.

Zinc determines dynamical properties and aggregation kinetics of human insulin

Protein aggregation is a widespread process leading to deleterious consequences in the organism, with amyloid aggregates being important not only in biology but also for drug design and biomaterial production. Insulin is a protein largely used in diabetes treatment, and its amyloid aggregation is at the basis of the so-called insulin-derived amyloidosis. Here, we uncover the major role of zinc in both insulin dynamics and aggregation kinetics at low pH, in which the formation of different amyloid superstructures (fibrils and spherulites) can be thermally induced. Amyloid aggregation is accompanied by zinc release and the suppression of water-sustained insulin dynamics, as shown by particle-induced x-ray emission and x-ray absorption spectroscopy and by neutron spectroscopy, respectively. Our study shows that zinc binding stabilizes the native form of insulin by facilitating hydration of this hydrophobic protein and suggests that introducing new binding sites for zinc can improve insulin stability and tune its aggregation propensity.

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.

Substrate Locking Promotes Dimer-Dimer Docking of an Enzyme Antibiotic Target

Protein dynamics manifested through structural flexibility play a central role in the function of biological molecules. Here we explore the substrate-mediated change in protein flexibility of an antibiotic target enzyme, Clostridium botulinum dihydrodipicolinate synthase. We demonstrate that the substrate, pyruvate, stabilizes the more active dimer-of-dimers or tetrameric form. Surprisingly, there is little difference between the crystal structures of apo and substrate-bound enzyme, suggesting protein dynamics may be important. Neutron and small-angle X-ray scattering experiments were used to probe substrate-induced dynamics on the sub-second timescale, but no significant changes were observed. We therefore developed a simple technique, coined protein dynamics-mass spectrometry (ProD-MS), which enables measurement of time-dependent alkylation of cysteine residues. ProD-MS together with X-ray crystallography and analytical ultracentrifugation analyses indicates that pyruvate locks the conformation of the dimer that promotes docking to the more active tetrameric form, offering insight into ligand-mediated stabilization of multimeric enzymes.

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.

Microbial Rhodopsins

Microbial rhodopsins (MRs) are a large family of photoactive membrane proteins, found in microorganisms belonging to all kingdoms of life, with new members being constantly discovered. Among the MRs are light-driven proton, cation and anion pumps, light-gated cation and anion channels, and various photoreceptors. Due to their abundance and amenability to studies, MRs served as model systems for a great variety of biophysical techniques, and recently found a great application as optogenetic tools. While the basic aspects of microbial rhodopsins functioning have been known for some time, there is still a plenty of unanswered questions. This chapter presents and summarizes the available knowledge, focusing on the functional and structural studies.

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.

Dynamical Transition of Collective Motions in Dry Proteins

Water is widely assumed to be essential for protein dynamics and function. In particular, the well-documented "dynamical" transition at ∼200  K, at which the protein changes from a rigid, nonfunctional form to a flexible, functional state, as detected in hydrogenated protein by incoherent neutron scattering, requires hydration. Here, we report on coherent neutron scattering experiments on perdeuterated proteins and reveal that a transition occurs in dry proteins at the same temperature resulting primarily from the collective heavy-atom motions. The dynamical transition discovered is intrinsic to the energy landscape of dry proteins.

Mobility of a Mononucleotide within a Lipid Matrix: A Neutron Scattering Study

An essential question in studies on the origins of life is how nucleic acids were first synthesized and then incorporated into compartments about 4 billion years ago. A recent discovery is that guided polymerization within organizing matrices could promote a non-enzymatic condensation reaction allowing the formation of RNA-like polymers, followed by encapsulation in lipid membranes. Here, we used neutron scattering and deuterium labelling to investigate 5′-adenosine monophosphate (AMP) molecules captured in a multilamellar phospholipid matrix. The aim of the research was to determine and compare how mononucleotides are captured and differently organized within matrices and multilamellar phospholipid structures and to explore the role of water in organizing the system to determine at which level the system becomes sufficiently anhydrous to lock the AMP molecules into an organized structure and initiate ester bond synthesis. Elastic incoherent neutron scattering experiments were thus employed to investigate the changes of the dynamic properties of AMP induced by embedding the molecules within the lipid matrix. The influence of AMP addition to the lipid membrane organization was determined through diffraction measurement, which also helped us to define the best working Q range for dynamical data analysis with respect to specific hydration. The use of different complementary instruments allowed coverage of a wide time-scale domain, from ns to ps, of atomic mean square fluctuations, providing evidence of a well-defined dependence of the AMP dynamics on the hydration level.

Computational investigation of dynamical transitions in Trp-cage miniprotein powders

We investigate computationally the dynamical transitions in Trp-cage miniprotein powders, at three levels of hydration: 0.04, 0.26 and 0.4 g water/g protein. We identify two distinct temperatures where transitions in protein dynamics occur. Thermal motions are harmonic and independent of hydration level below T_low ≈ 160 K, above which all powders exhibit harmonic behavior but with a different and enhanced temperature dependence. The second onset, which is often referred to as the protein dynamical transition, occurs at a higher temperature T_D that decreases as the hydration level increases, and at the lowest hydration level investigated here (0.04 g/g) is absent in the temperature range we studied in this work (T ≤ 300 K). Protein motions become anharmonic at T_D, and their amplitude increases with hydration level. Upon heating above T_D, hydrophilic residues experience a pronounced enhancement in the amplitude of their characteristic motions in hydrated powders, whereas it is the hydrophobic residues that experience the more pronounced enhancement in the least hydrated system. The dynamical transition in Trp-cage is a collective phenomenon, with every residue experiencing a transition to anharmonic behavior at the same temperature.

Influence of myelin proteins on the structure and dynamics of a model membrane with emphasis on the low temperature regime

Myelin is an insulating, multi-lamellar membrane structure wrapped around selected nerve axons. Increasing the speed of nerve impulses, it is crucial for the proper functioning of the vertebrate nervous system. Human neurodegenerative diseases, such as multiple sclerosis, are linked to damage to the myelin sheath through demyelination. Myelin exhibits a well defined subset of myelin-specific proteins, whose influence on membrane dynamics, i.e., myelin flexibility and stability, has not yet been explored in detail. In a first paper [W. Knoll, J. Peters, P. Kursula, Y. Gerelli, J. Ollivier, B. Demé, M. Telling, E. Kemner, and F. Natali, Soft Matter 10, 519 (2014)] we were able to spotlight, through neutron scattering experiments, the role of peripheral nervous system myelin proteins on membrane stability at room temperature. In particular, the myelin basic protein and peripheral myelin protein 2 were found to synergistically influence the membrane structure while keeping almost unchanged the membrane mobility. Further insight is provided by this work, in which we particularly address the investigation of the membrane flexibility in the low temperature regime. We evidence a different behavior suggesting that the proton dynamics is reduced by the addition of the myelin basic protein accompanied by negligible membrane structural changes. Moreover, we address the importance of correct sample preparation and characterization for the success of the experiment and for the reliability of the obtained results.

Elasticity and Inverse Temperature Transition in Elastin

Elastin is a structural protein and biomaterial that provides elasticity and resilience to a range of tissues. This work provides insights into the elastic properties of elastin and its peculiar inverse temperature transition (ITT). These features are dependent on hydration of elastin and are driven by a similar mechanism of hydrophobic collapse to an entropically favorable state. Using neutron scattering, we quantify the changes in the geometry of molecular motions above and below the transition temperature, showing a reduction in the displacement of water-induced motions upon hydrophobic collapse at the ITT. We also measured the collective vibrations of elastin gels as a function of elongation, revealing no changes in the spectral features associated with local rigidity and secondary structure, in agreement with the entropic origin of elasticity.

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.

Protein dynamics: from rattling in a cage to structural relaxation

We present an overview of protein dynamics based mostly on results of neutron scattering, dielectric relaxation spectroscopy and molecular dynamics simulations. We identify several major classes of protein motions on the time scale from faster than picoseconds to several microseconds, and discuss the coupling of these processes to solvent dynamics. Our analysis suggests that the microsecond backbone relaxation process might be the main structural relaxation of the protein that defines its glass transition temperature, while faster processes present some localized secondary relaxations. Based on the overview, we formulate a general picture of protein dynamics and discuss the challenges in this field.

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.

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.

Methyl rotors in flavoproteins

ENDOR evidence shows that methyl groups in flavin behave as quantum locked rotors.

Dynamics measured by neutron scattering correlates with the organization of bioenergetics complexes in natural membranes from hyperthermophile and mesophile bacteria

Various models on membrane structure and organization of proteins and complexes in natural membranes emerged during the last years. However, the lack of systematic dynamical studies to complement structural investigations hindered the establishment of a more complete picture of these systems. Elastic incoherent neutron scattering gives access to the dynamics on a molecular level and was applied to natural membranes extracted from the hyperthermophile Aquifex aeolicus and the mesophile Wolinella succinogenes bacteria. The results permitted to extract a hierarchy of dynamic flexibility and atomic resilience within the samples, which correlated with the organization of proteins in bioenergetics complexes and the functionality of the membranes.

Origin of abrupt rise in deuteron NMR longitudinal relaxation times of protein methyl groups below 90 K

In order to examine the origin of the abrupt change in the temperature dependence of (2)H NMR longitudinal relaxation times observed previously for methyl groups of L69 in the hydrophobic core of villin headpiece protein at around 90 K (Vugmeyster et al. J. Am. Chem. Soc. 2010, 132, 4038-4039), we extended the measurements to several other methyl groups in the hydrophobic core. We show that, for all methyl groups, relaxation times experience a dramatic jump several orders of magnitude around this temperature. Theoretical modeling supports the conclusion that the origin of the apparent transition in the relaxation times is due to the existence of the distribution of conformers distinguished by their activation energy for methyl three-site hops. It is also crucial to take into account the differential contribution of individual conformers into overall signal intensity. When a particular conformer approaches the regime at which its three-site hop rate constant is on the order of the quadrupolar coupling interaction constant, the intensity of the signal due to this conformer experiences a sharp drop, thus changing the balance of the contributions of different conformers into the overall signal. As a result, the observed apparent transition in the relaxation rates can be explained without the assumption of an underlying transition in the rate constants. This work in combination with earlier results also shows that the model based on the distribution of conformers explains the relaxation behavior in the entire temperature range between 300 and 70 K.

Protein surface and core dynamics show concerted hydration-dependent activation

By specifically labeling leucine/valine methyl groups and lysine side chains "inside" and "outside" dynamics of proteins on the nanosecond timescale are compared using neutron scattering. Surprisingly, both groups display similar dynamics as a function of temperature, and the buried hydrophobic core is sensitive to hydration and undergoes a dynamical transition.

Dynamics of protein and its hydration water: neutron scattering studies on fully deuterated GFP

We present a detailed analysis of the picosecond-to-nanosecond motions of green fluorescent protein (GFP) and its hydration water using neutron scattering spectroscopy and hydrogen/deuterium contrast. The analysis reveals that hydration water suppresses protein motions at lower temperatures (<~ 200 K), and facilitates protein dynamics at high temperatures. Experimental data demonstrate that the hydration water is harmonic at temperatures <~ 180-190 K and is not affected by the proteins' methyl group rotations. The dynamics of the hydration water exhibits changes at ~ 180-190 K that we ascribe to the glass transition in the hydrated protein. Our results confirm significant differences in the dynamics of protein and its hydration water at high temperatures: on the picosecond-to-nanosecond timescale, the hydration water exhibits diffusive dynamics, while the protein motions are localized to <~3 Å. The diffusion of the GFP hydration water is similar to the behavior of hydration water previously observed for other proteins. Comparison with other globular proteins (e.g., lysozyme) reveals that on the timescale of 1 ns and at equivalent hydration level, GFP dynamics (mean-square displacements and quasielastic intensity) are of much smaller amplitude. Moreover, the suppression of the protein dynamics by the hydration water at low temperatures appears to be stronger in GFP than in other globular proteins. We ascribe this observation to the barrellike structure of GFP.

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

Neutron scattering reveals a complex dynamics in polypeptide chains, with two main onsets of anharmonicity whose physical origin and biological role are still debated. In this study the dynamics of strategically selected homomeric polypeptides is investigated with elastic neutron scattering using different energy resolutions and compared with that of a real protein. Our data spotlight the dependence of anharmonic transition temperatures and fluctuation amplitudes on energy resolution, which we quantitatively explain in terms of a two-site model for the protein-hydration water energy landscape. Experimental data strongly suggest that the protein dynamical transition is not a mere resolution effect but is due to a real physical effect. Activation barriers and free energy values obtained for the protein dynamical transition allow us to make a connection with the two-well interaction potential of supercooled-confined water proposed to explain a low-density→high-density liquid-liquid transition.

Thermal fluctuations of haemoglobin from different species: adaptation to temperature via conformational dynamics

Thermodynamic stability, configurational motions and internal forces of haemoglobin (Hb) of three endotherms (platypus, Ornithorhynchus anatinus; domestic chicken, Gallus gallus domesticus and human, Homo sapiens) and an ectotherm (salt water crocodile, Crocodylus porosus) were investigated using circular dichroism, incoherent elastic neutron scattering and coarse-grained Brownian dynamics simulations. The experimental results from Hb solutions revealed a direct correlation between protein resilience, melting temperature and average body temperature of the different species on the 0.1 ns time scale. Molecular forces appeared to be adapted to permit conformational fluctuations with a root mean square displacement close to 1.2 Å at the corresponding average body temperature of the endotherms. Strong forces within crocodile Hb maintain the amplitudes of motion within a narrow limit over the entire temperature range in which the animal lives. In fully hydrated powder samples of human and chicken, Hb mean square displacements and effective force constants on the 1 ns time scale showed no differences over the whole temperature range from 10 to 300 K, in contrast to the solution case. A complementary result of the study, therefore, is that one hydration layer is not sufficient to activate all conformational fluctuations of Hb in the pico- to nanosecond time scale which might be relevant for biological function. Coarse-grained Brownian dynamics simulations permitted to explore residue-specific effects. They indicated that temperature sensing of human and chicken Hb occurs mainly at residues lining internal cavities in the β-subunits.

Anharmonic transitions in nearly dry L-cysteine I

Two special dynamical transitions of universal character have recently been observed in macromolecules (lysozyme, myoglobin, bacteriorhodopsin, DNA and RNA) at T* ~100-150 K and T(D) ~180-220 K. The underlying mechanisms governing these transitions have been the subject of debate. In the present work, a survey is reported on the temperature dependence of structural, vibrational and thermodynamical properties of a nearly anhydrous amino acid (orthorhombic polymorph of the amino acid l-cysteine at a hydration level of 3.5%). The temperature dependence of x-ray powder diffraction patterns, Raman spectra and specific heat revealed these two transitions at T* = 70 K and T(D) = 230 K for this sample. The data were analyzed considering amino acid-amino acid, amino acid-water, water-water phonon-phonon interactions and molecular rotor activation. Our results indicated that the two referred temperatures define the triggering of very simple and particular events that govern all the interactions of the biomolecular: activation of CH(2) rigid rotors (T < T* ), phonon-phonon interactions between specific amino acid and water dimer vibrational modes (T* < T < T(D)), and water rotational barriers surpassing (T > T(D)).

Role of methyl groups in dynamics and evolution of biomolecules

Recent studies have discovered strong differences between the dynamics of nucleic acids (RNA and DNA) and proteins, especially at low hydration and low temperatures. This difference is caused primarily by dynamics of methyl groups that are abundant in proteins, but are absent or very rare in RNA and DNA. In this paper, we present a hypothesis regarding the role of methyl groups as intrinsic plasticizers in proteins and their evolutionary selection to facilitate protein dynamics and activity. We demonstrate the profound effect methyl groups have on protein dynamics relative to nucleic acid dynamics, and note the apparent correlation of methyl group content in protein classes and their need for molecular flexibility. Moreover, we note the fastest methyl groups of some enzymes appear around dynamical centers such as hinges or active sites. Methyl groups are also of tremendous importance from a hydrophobicity/folding/entropy perspective. These significant roles, however, complement our hypothesis rather than preclude the recognition of methyl groups in the dynamics and evolution of biomolecules.Electronic supplementary material The online version of this article (doi:10.1007/s10867-012-9268-6) contains supplementary material, which is available to authorized users.

Derivation of mean-square displacements for protein dynamics from elastic incoherent neutron scattering

The derivation of mean-square displacements from elastic incoherent neutron scattering (EINS) of proteins is examined, with the aid of experiments on camphor-bound cytochrome P450cam and complementary molecular dynamics simulations. It is shown that a q(4) correction to the elastic incoherent structure factor (where q is the scattering vector) can be simply used to reliably estimate from the experiment both the average mean-square atomic displacement, <Δr(2)> of the nonexchanged hydrogen atoms in the protein and its variance, σ(2). The molecular dynamics simulation results are in broad agreement with the experimentally derived <Δr(2)> and σ(2) derived from EINS on instruments at two different energy resolutions, corresponding to dynamics on the ∼100 ps and ∼1 ns time scales. Significant dynamical heterogeneity is found to arise from methyl-group rotations. The easy-to-apply q(4) correction extends the information extracted from elastic incoherent neutron scattering experiments and should be of wide applicability.

Further evidence that interfacial water is the main "driving force" of protein dynamics: a neutron scattering study on perdeuterated C-phycocyanin

The fundamental role of hydration water (also called interfacial water) is widely recognized in protein flexibility, especially in the existence of the so-called protein "dynamical transition" at around 220 K. In the present study, we take advantage of perdeuterated C-phycocyanin (CPC) and elastic incoherent neutron scattering (EINS) to distinguish between protein dynamics and interfacial water dynamics. Powders of hydrogenated (hCPC) and perdeuterated (dCPC) CPC protein have been hydrated, respectively, with D(2)O or H(2)O and measured by EINS to separately probe protein dynamics (hCPC/D(2)O) and water dynamics (dCPC/H(2)O) at different time- and length-scales. We find that "fast" (<20 ps) local mean-square displacements (MSD) of both protein and interfacial water coincide all along the temperature range, with the same dynamical transition temperature at ~220 K. On higher resolution (<400 ps), two different types of motions can be separated: (i) localized motions with the same amplitude for CPC and hydration water and two transitions at ~170 and ~240 K for both; (ii) large scale fluctuations exhibiting for both water molecules and CPC protein a single transition at ~240 K, with a significantly higher amplitude for the interfacial water than for CPC. Moreover, by comparing these motions with bulk water MSD measured under the same conditions, we show no coupling between bulk water dynamics and protein dynamics all along the temperature range. These results show that interfacial water is the main "driving force" governing both local and large scale motions in proteins.

Low-temperature (293-60 K) Raman spectra of the crystalline tert-butylisocyanide complexes, Cr(CO)5(CNBu-t) and cis-Cr(CO)4(CNBu-t)2

Low-temperature (293–60 K) Raman spectra have been recorded in the 3500–50 cm–1 region for the crystalline tert-butylisocyanide (CNBu-t) complexes, Cr(CO)5(CNBu-t) (I) and cis-Cr(CO)4(CNBu-t)2 (II). Several of the factor group splittings that are predicted for the νCO, νCN, νCrC, νCrCN and δCrCO modes on the basis of the P21/a (C2h5) (I) and Pccn (D2h10) (II) crystallographic symmetries can be detected at very low temperatures. No phase transitions are observed for either complex throughout the temperature range investigated.

Crowding induces differences in the diffusion of thermophilic and mesophilic proteins: a new look at neutron scattering results

The dynamical basis underlying the increased thermal stability of thermophilic proteins remains uncertain. Here, we challenge the new paradigm established by neutron scattering experiments in solution, in which the adaptation of thermophilic proteins to high temperatures lies in the lower sensitivity of their flexibility to temperature changes. By means of a combination of molecular dynamics and Brownian dynamics simulations, we report a reinterpretation of those experiments and show evidence that under crowding conditions, such as in vivo, thermophilic and homolog mesophilic proteins have diffusional properties with different thermal behavior.

Redox-promoting protein motions in rubredoxin

Proteins are dynamic objects, constantly undergoing conformational fluctuations, yet the linkage between internal protein motion and function is widely debated. This study reports on the characterization of temperature-activated collective and individual atomic motions of oxidized rubredoxin, a small 53 residue protein from thermophilic Pyrococcus furiosus (RdPf). Computational modeling allows detailed investigations of protein motions as a function of temperature, and neutron scattering experiments are used to compare to computational results. Just above the dynamical transition temperature which marks the onset of significant anharmonic motions of the protein, the computational simulations show both a significant reorientation of the average electrostatic force experienced by the coordinated Fe(3+) ion and a dramatic rise in its strength. At higher temperatures, additional anharmonic modes become activated and dominate the electrostatic fluctuations experienced by the ion. At 360 K, close to the optimal growth temperature of P. furiosus, simulations show that three anharmonic modes including motions of two conserved residues located at the protein active site (Ile7 and Ile40) give rise to the majority of the electrostatic fluctuations experienced by the Fe(3+) ion. The motions of these residues undergo displacements which may facilitate solvent access to the ion.

Comparative dynamics of leucine methyl groups in FMOC-leucine and in a protein hydrophobic core probed by solid-state deuteron nuclear magnetic resonance over 7-324 K temperature range

Quantitative dynamics of methyl groups in 9-fluorenylmethyloxycarbonyl-leucine (FMOC-leu) have been analyzed and compared with earlier studies of methyl dynamics in chicken villin headpiece subdomain protein (HP36) labeled at L69, a key hydrophobic core position. A combination of deuteron solid-state nuclear magnetic resonance experiments over the temperature range of 7-324 K and computational modeling indicated that while the two compounds show the same modes of motions, there are marked differences in the best-fit parameters of these motions. One of the main results is that the crossover observed in the dynamics of the methyl groups in the HP36 sample at 170 K is absent in FMOC-leu. A second crossover at around 95-88 K is present in both samples. The differences in the behavior of the two compounds suggest that some of the features of methyl dynamics reflect the complexity of the protein hydrophobic core and are not determined solely by local interactions.

Temperature- and hydration-dependent internal dynamics of stripped human erythrocyte vesicles studied by incoherent neutron scattering

BACKGROUND

We focus on temperature- and hydration-dependence of internal molecular motions in stripped human red blood cell (RBC) vesicles, widely used as a model system for more complex biomembranes.

METHODS

We singled out picosecond local motions of the non-exchangeable hydrogen atoms of RBC vesicles by performing elastic and quasielastic incoherent neutron scattering measurements in dry and heavy water (D₂O)-hydrated RBC powders.

RESULTS

In dry stripped RBCs, hydrogen motions remained harmonic all along the measured temperature range (100-310K) and mean-square displacements (MSDs) exhibited no temperature transition up to 310K. In contrast, MSDs of hydrated stripped RBCs (h ≈ 0.38g D₂O/g dry powder) exhibited a pronounced transition near 260K, with the sharp rise of anharmonic diffusive motions of hydrogen atoms. This transition at ~260K was correlated with both the onset of nonvibrational (harmonic and nonharmonic) motions and the melting of crystallized hydration water.

GENERAL SIGNIFICANCE

In conclusion, we have shown that MSDs in human RBC vesicles are temperature-and hydration-dependent. These results provide insight into biomembrane internal dynamics at picosecond timescale and nanometer length scale. Such motions have been shown to act as the "lubricant" of larger conformational changes on a slower, millisecond timescale that are necessary for important biological processes.