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.


I. Introduction
Biological macromolecules play a crucial role in all processes of life, from catalysis of biochemical reactions, to transport and to genetic codes.Proteins, RNAs and DNA work as small machines executing their functions through stochastic thermal fluctuations.It is well accepted that dynamics of these macromolecules are the key to their biochemical activity and function.Biological macromolecules in their normal state (e.g.under physiological conditions) do not possess a fixed structure.][3][4][5][6][7][8][9][10][11] Currently we have reasonable knowledge of the structure of many proteins, and a clear understanding of the importance of protein dynamics and conformational fluctuations to their function.Yet our knowledge and understanding of the internal protein dynamics is very limited.The same is true for dynamics of RNA and DNA.Even a qualitative atomistic picture of dynamics of biological macromolecules is still absent.This is in contrast to the dynamics of other soft materials (e.g.synthetic polymers), where general classification of various relaxation processes exists and their characteristic behavior is known (although in many cases detailed microscopic mechanisms remain to be explained). 124][15][16][17] They include traditional polymer models, such as worm-like and Rouse/Zimm models with internal friction, 13,15 and colloidal models. 16,17However, as it has been discussed in several papers, these models either fail to describe all the data, or have to involve additional free parameters describing e.g.internal friction that is directly related to conformational fluctuations in biomacromolecules. 11,17any experimental and computational methods have been applied to studies of dynamics of biological macromolecules.Among them, NMR, molecular dynamics (MD) simulations and neutron scattering present the methods that analyze atomic motions directly.9][20] NMR presents very local information on motions of specific groups in protein structures.A very broad array of NMR techniques was able to analyze dynamics of specific groups and residues in the time range from tens of ps up to ms. [21][22][23][24][25][26] Neutron scattering is a unique technique that provides analysis of dynamics not only in a broad time range, but also reveals the geometry of the underlying motions, can probe large scale (e.g.][29][30][31][32][33][34][35][36][37][38][39] The high contrast between neutron scattering of hydrogen and deuterium atoms provides additional advantage for this technique in studies of biological molecules.For example, deuteration of solvent or a part of the molecule helps to drastically reduce their contribution to the neutron scattering spectra and separate the dynamics of the targeted part (containing H-atoms).Moreover, the energy of neutrons used in these experiments is in the meV range and they are not destructive for biological systems (in contrast to X-ray).1][42][43] The results, however, might depend strongly on force-fields and approximations used.The best approach in this case is a combination of MD-simulations with neutron scattering, because the measured experimental parameters can be directly calculated from MD simulated atomic trajectories.This also provides calibration and validation of simulations.

S. Khodadadi
A significant drawback of neutron scattering is the limited frequency (time) range and low accuracy of the measured spectra.In contrast, broadband dielectric spectroscopy provides dynamic measurements in an extremely broad frequency range with very high accuracy.The measured spectra reflect reorientational motions of dipoles and translational motion of charges. 44This technique, however, does not provide microscopic information on the molecular motions involved.][47][48][49][50] In recent years significant progress has been achieved in single molecular studies of biomolecules.2][53][54][55] These studies usually provide dynamics on a time scale longer than milliseconds and length scales of a few nm, missing atomistic details.Combination of the results obtained with these techniques with atomistic details learned from NMR, neutrons and MD-simulations might result in a significant breakthrough in our understanding of biomolecular dynamics.This topic, however, is out of scope of the current review.
Here we present an overview of studies of intra-molecular protein dynamics based mostly on neutron scattering, broadband dielectric spectroscopy and MD simulation results, although results from other techniques are also included.We focus on equilibrium fluctuations and conformational changes in the folded state on time scales from sub-picosecond to microseconds, and discuss the influence of temperature and hydration level on protein dynamics.We intentionally exclude the discussion of the folding/unfolding process, because it presents a separate important topic.In this review we aim to formulate a general atomistic picture of protein dynamics with classification of major relaxation processes that control conformational fluctuations.This classification is based on analogy with dynamics of soft materials, and aims at providing a framework to describe complex dynamics of biological macromolecules.We emphasize specificity of the relaxation processes in proteins that differentiate them from the classical dynamic processes of other soft materials.

II. Dynamics of soft materials
Soft Matter includes a wide range of materials from polymers and glass-forming liquids, to colloids and biological systems.The particular properties that differentiate soft materials from other materials include the existence of a large number of metastable states with comparable potential energies that are separated by relatively small energy barriers.So, there are constant fluctuations between different conformational states that result in a delicate balance between the entropic and enthalpic contributions to the free energy, and large thermodynamically equilibrium fluctuations, significant dynamic heterogeneity and strong sensitivity to small external forces and perturbations. 12,56These properties emphasize the distinct characteristics of the soft materials: permanent transitions between multiple meta-stable conformational states under relatively small perturbations or due to equilibrium fluctuations.So, the dynamics, i.e. motions of molecular (or other structural) units are the key to the main macroscopic properties of Soft Matter.This is directly applicable to biological macromolecules.They exhibit all the properties listed above and are very sensitive to slight variations in temperature, pH and other external perturbations.
The dynamical processes that take place in soft materials are very complex and the microscopic mechanisms of many of them remain poorly understood.7][58][59][60][61]  of the dynamic processes in soft materials with the example of polymer -polyisoprene (PI).One of them is fast dynamics that includes low-frequency collective vibrations, the so-called boson peak, and fast picosecond fluctuations. 62,63The former present excess vibrational modes relative to what is expected in a simple Debye model for acoustic modes. 64,65Debye-like behavior in density of vibrational states is usually observed in crystals, but all disordered systems show excess vibrational density of states.The nature of these excess vibrations remains a topic of active discussion, but most authors ascribe them to heterogeneity in elastic constants 66,67 or to soft modes. 68,69The fast picosecond fluctuations are usually ascribed to rattling of a structural unit in a cage formed by neighboring units and is considered to be a precursor of structural relaxation. 70,71ts characteristic time scale has weak temperature variations reflecting very low energy barriers. 72econdary relaxations are usually chemically specific and present mainly local conformational changes. 12The characteristic relaxation time has Arrhenius temperature dependence, t = t 0 exp(E/k B T), where t is the relaxation time, E a is the activation energy, k B is the Boltzman constant and T is temperature (Fig. 1).The term secondary relaxation includes a broad class of relaxation processes from some intra-molecular and side-group motions (e.g.methyl group rotation), to localized processes that are reminiscent of structural relaxation, the so-called Johari-Goldstein b-relaxation. 12Most of the secondary relaxations separate clearly from the main structural relaxation only at lower temperatures close to T g , although some motions, such as methyl group rotation, remain faster than structural dynamics up to very high temperatures (Fig. 1).
The main structural relaxation (often called a-relaxation) is the main relaxation process that defines transition from a solid ( jammed) to a liquid (flow) state.This relaxation process leads to the glass transition and is a collective motion that involves many structural units.It controls viscosity and diffusion in nonpolymeric systems.It is called segmental relaxation in polymers where viscosity is defined by chain relaxation (see below).The main structural relaxation in most of the systems exhibits strongly non-Arrhenius temperature dependence that is traditionally described by the Vogel-Fulcher-Tammann (VFT) equation: t = t 0 exp[B/(T À T 0 )], where B and T 0 are material dependent constants (Fig. 1).While the mechanism of this sharp slowing down of structural relaxation time with temperature decrease still remains a subject of active discussions, 12,56,73,74 the role of dynamic cooperativity and heterogeneity becomes the center of attention in current research. 56,57,75,76he main structural relaxation is the longest relaxation process in non-polymeric systems.However, chain connectivity in polymers leads to an even slower relaxation process in macromolecules.Chain relaxation (Fig. 1) is the slowest relaxation process in polymers that defines their viscosity and macromolecular diffusion.Chain relaxation depends on the length of the chain (e.g.molecular weight) and is usually described by Rouse and reptation models. 77

III. Dynamics of biomolecules
There is an analogy between dynamics of biological macromolecules (proteins, DNA and RNA) and other soft materials with many similarities, while there are also clear differences due to the globular nature of the proteins.Unlike in many other soft materials, a generally accepted picture for the protein dynamics even on a qualitative level is still missing.][5]78 There is a clear analogy between this approach and energy landscape models used for other soft materials. 58In this review, however, we will follow the same classification for dynamics of the soft materials as has been described in the previous section on example of a polymer.Studies discussed in this review cover a broad temperature range, far outside of the usual ambient temperature range.Although such a broad temperature range is not directly relevant to most of the biological functions and processes (except bio-and cryo-preservation), it is a traditional approach in physics that helps to disentangle various dynamic processes, analyze potential energy landscape, separate enthalpic and entropic contributions.

III.1 Fast dynamics
0][81][82] Neutron scattering spectra are usually presented as a dynamic structure factor S(Q,n), which is the measured intensity normalized to the sample and scattering geometry; here Q -is the scattering wave vector and hn is the change in energy between the incoming and scattered neutrons (h is the Planck constant and n is the frequency).The boson peak appears as an inelastic peak in S(Q,n) at frequencies B1 THz, while fast relaxation appears as a quasielastic broadening at lower Fig. 1 Traditional relaxation map of a soft material is illustrated with the example of the polymer polyisoprene. 27It includes (i) chain and (ii) segmental relaxation processes with strongly non-Arrhenius temperature variations; (iii) a secondary relaxation; (iv) methyl group rotation with Arrhenius temperature dependence; and (v) a fast picosecond relaxation with a rather weak temperature dependence of the characteristic relaxation time. 27The characteristic relaxation times are usually estimated from the maximum of the loss spectra (e.g.mechanical, dielectric, etc.), width of quasielastic neutron or light scattering spectra, or time decay of correlation functions of the corresponding relaxation process.
frequencies (Fig. 2).Analysis of neutron and light scattering spectra demonstrates that the fast dynamics depend strongly on hydration level (Fig. 2) and temperature. 33,39,83,84According to simulations, 85 the boson peak vibrations present collective motions of amino acids that involve the whole protein and even water of hydration.Moreover, collective protein modes seem to propagate to the hydration water. 86It has been discovered recently that the boson peak in many proteins has a universal spectral shape that is similar to the spectral shape observed in glassforming systems. 87The properties of the fast picosecond relaxation are similar to those of the fast relaxation known in all soft materials.It depends strongly on temperature and solvent, and can be ascribed to a ''rattling'' of amino acids in a cage formed by neighbor residues and solvent molecules. 39It can also be described as small-scale conformational fluctuations between different conformational states separated by a very low-energy barrier.The amplitude of these fluctuations is relatively small with mean-squared displacements (MSD) hr 2 i B 0.15-0.4Å 2 under ambient conditions. 37,39,88[91]

III.2 Methyl dynamics
Various side group motions (e.g.methyl group rotations, phenyl ring flips) exist in proteins.These relaxations usually have a broad distribution of relaxation times with Arrhenius-like temperature dependence.3][94][95] All proteins contain a significant number of methyl groups.They provide proteins hydrophobicity, and they also play an important role in facilitating protein dynamics.Combined analysis of NMR and MD-simulations data suggests that methyl group dynamics can be used as an entropy meter for entire protein, 96 because their NMR order parameter correlates well with conformational motions of side groups.][94][95] These observations suggest that methyl groups might play the role of internal plasticizers by facilitating protein dynamics even under extreme conditions of low hydration or temperature.The methyl groups show broad distribution of relaxation times that depends on the chemical structure of the residues and their position in the protein.It has been noticed in lysozyme and in myoglobin that the fastest methyl groups are placed around the active sites. 19,39,42,94,97,98Apparently there might be a particular role that nature assigns to methyl groups in proteins.It is possible that methyl groups and their position in proteins play the key role in providing conformational flexibility and additional entropy for active sites.It is important to notice that RNA has no methyl groups and only one nucleic acid of DNA has a single methyl group (Fig. 3C).

III.3 Sub-nanosecond process: coupled protein-solvent relaxation
3][104][105] The time scale of the protein process and its temperature dependence appear to be very similar to those of the hydration water (Fig. 5).The latter is B2-4 times slower than bulk water. 102,103However, analysis of the Q dependence revealed that hydration water shows diffusive LG -lysozyme in glycerol; LT -lysozyme in trehalose). 39,90(A) As the hydration level increases the Boson peak shifts to higher frequencies and the quasielastic contribution coming mainly from the fast relaxation process strongly increases.(B) At low T = 150 K the fast relaxation process is strongly suppressed in LG.While as the temperature is increased to ambient T B 295 K. (C) WL has the most contribution from the fast relaxation process.This behavior of the protein's fast relaxation follows the behavior of the solvent. 79,89his journal is © The Royal Society of Chemistry 2015 translational motion, while this protein process is a localized motion, usually described as a ''localized diffusion''. 6,102nalysis of MSD (Fig. 3D) demonstrates this point very clearly: MSD slowly increases due to this process, but saturates at time scales longer than B100-200 ps.Analysis of the atomic trajectories in MD-simulations 6 revealed Brownian-like motions in a confined space, with radius of confinement B1-1.5 Å.This leads to the name of the process -'localized diffusion'.
Since the pioneering work of Doster and co-workers, 32 the protein relaxation spectra have also been analyzed in susceptibility presentation: where n(n,T) = [exp (hn/kT) À 1] À1 is the Bose population factor.This presentation is similar to dielectric or mechanical loss spectra and has several advantages in comparison to the traditional  94 Proteins have a significant number of methyl groups in comparison to DNA (one out of four nucleic acids has a single methyl group) and RNA (no methyl in primary structure of RNA).(D) An example of MSD of hydrogen atoms in lysozyme involved in the localized diffusion motion as a function of time obtained from MD-simulations. 6It demonstrates fast dynamics and caging (plateau) at times shorter than B3 ps, followed by the localized diffusion with an almost logarithmic increase in MSD up to B200 ps.MSD saturates at longer times indicating localized motion.dynamic structure factor.The susceptibility spectra of the protein process have symmetric spectral shape that is usually well described by a Cole-Cole distribution function, where o is angular frequency, with the stretching parameter a B 0.25-0.4(Fig. 4). 32,39,106Such a strongly stretched process suggests significant dynamic heterogeneity in this motion.The appearance of this process and its relaxation time (defined from the frequency of the maximum t = 1/o max ) depends strongly on protein hydration, 39,48 and solvent in general. 91,107112 MD-simulations demonstrate that not only surface, but also core residues are involved in this motion.However, according to various simulations this process might involve mostly side group motions with much smaller involvement of the backbone. 9,40,41According to analysis of MD-simulations of lysozyme and Green Fluorescence Protein (GFP), this process can be described as an over-damped diffusive motion in a harmonic potential well. 6,113So, these motions are not jumps between different well-defined conformational states, but rather present a diffusion in a restricted space with a particular friction.5][116] However, NMR, neutron scattering and simulations all reveal the existence of protein dynamics on the same time scale B10-50 ps. 6,18,38,39,48,84,101,102,106,117Coincidence of the characteristic time scales and similar temperature dependences of hydration water and protein dynamics (Fig. 5B and C) suggest strong coupling of those two processes.Apparently the translational motion of the hydration water molecules is strongly Fig. 5 (A) Temperature dependence of characteristic relaxation times of hydrated lysozyme, 48,122 myoglobin 123 and tRNA, 122 as estimated from broadband dielectric (closed symbols) and neutron scattering (open symbols) spectroscopies.In all cases the relaxation time was estimated from the susceptibility maximum (Fig. 4) or dielectric loss maximum, t = (2pn max ) À1 .(B) Temperature dependence of characteristic relaxation times estimated from quasielastic neutron scattering spectra of tRNA, 122 lysozyme, 48,122 and DNA 124 hydrated with D 2 O; (C) characteristic relaxation times of hydration water of RNA, 121 lysozyme, 103 and DNA, 120 estimated using quasielastic neutron scattering spectra.
coupled to this localized motion of proteins.This has been clearly demonstrated from MD-simulations where the translational motion of water molecules has been artificially restricted. 118This restriction of water motions strongly suppressed the protein relaxation on the sub-nanosecond time scale.
We emphasize that a similar relaxation process has also been observed in RNA and DNA dynamics (Fig. 4B), but in t-RNA it has a significantly larger localization length B7 Å. 119 In the case of nucleic acids it exhibits stronger temperature dependence, with t-RNA showing slower dynamics than proteins like lysozyme and myoglobin, and DNA showing even slower dynamics than t-RNA (Fig. 5B).In all these cases it has been found that dynamics of hydration water shows the same difference (Fig. 5C). 103,120,121This result suggests that dynamics of biomolecules not only strongly coupled to hydration water dynamics, but also strongly influence the latter.There is no simple 'slaving' of biomolecular dynamics by its hydration water, rather there is a coupled motion where the chemical and 3-D structure of biomolecules play important roles.

III.4 Nanosecond dynamics: conformational jumps
Detailed analysis of MD-simulations suggested the existence of jump-like motions in proteins on the nanosecond time scale. 6,9eutron scattering on protein solutions also detected relaxation processes on the time scale B10-30 ns. 10,125,126They present rather large atomic displacements B8 Å, and were ascribed to domain-like motions. 10,28,29,1255,125 A similar process has recently been reported for an intrinsically disordered protein Myelin basic protein. 11Apparently this nanosecond process does not necessarily present domain-like motion, but can be ascribed to more general conformational changes in proteins.In contrast to above-described ''localized diffusion'' process on the sub-nanosecond time scale, the nanosecond process presents jumps between well-defined conformational states. 6,10,28,29,125Protein relaxations in the ns time range has also been detected in many NMR studies. 18,19Recent detailed MD-simulation studies suggested broad distribution of energy barriers for conformational motions of side groups. 127,128They span from B10 kJ mol À1 up to B35 kJ mol À1 and higher. 127ethyl groups dominate the lower energy barriers motions, while aromatic side groups usually have higher energy barriers. 127,128here is also a heterogeneity in the behavior of hydrophobic and hydrophilic groups. 128][131][132][133][134] While this dielectric process is often ascribed to hydration water strongly bound to proteins, 115,123,129,132,134,135 its large dielectric amplitude would require hundreds of tightly bound hydration water molecules per protein with dynamics of about B10 2 -10 3 times slower than bulk water molecules. 41Such a behavior of hydration water contradicts current neutron scattering, NMR and simulations studies. 6,10,18,19,28,29,38,41,101,125,136hus the nanosecond relaxation process should be attributed to the internal protein dynamics.The latter cause water continuum motions that might induce additional large dielectric response.
Temperature dependence of the nanosecond process has been analyzed mostly by dielectric relaxation spectroscopy. 48,130,131,133t shows slightly non-Arrhenius behavior similar to the subnanosecond relaxation process, and is B10 2 -10 3 times slower than the latter (Fig. 6).This process in dielectric spectra slows down strongly with a decrease in the hydration level. 48,130,131,133he process is also strongly stretched and its characteristic relaxation time, according to dielectric spectroscopy, is comparable in different globular proteins, regardless of their size. 137This is consistent with a rather universal distribution of energy barriers for side group motions of different proteins suggested from MD-simulation analysis in ref. 127.This process definitely presents conformational jumps between well-defined states, and analysis of MD-simulation results suggests that the nanosecond process involves more side group motions than the motion of the backbone. 9,41,128It is most probably dominated by their conformational jumps with relatively large amplitude, and energy barriers controlled by the chemical structure, interactions and steric constraints.It involves the protein surface and core. 41,128ccording to neutron scattering, motions of secondary structures, e.g.domain-like motions in the case of multi-domain proteins, also appear in the ns time range.Its temperature dependence is most probably controlled by the behavior of the sub-nanosecond localized diffusion process that effectively defines a local friction for the nanosecond process.More detailed studies using MD-simulations and neutron scattering would be very helpful in unravelling microscopic details of the proteins nanosecond relaxation.

III.5 Microsecond dynamics: main structural relaxation?
Very recent dielectric spectroscopy studies of several hydrated protein powders revealed a relaxation process in the microsecond time scale. 130,131,133This process has an anomalously large dielectric amplitude, De B 10 3 -10 4 . 130,131,133This is B10-100 times larger than the dielectric amplitude of relaxation in water, one of the most polar solvent.It is not possible to analyze this process in solution using dielectric spectroscopy because its characteristic time scale (B1-20 ms) is much longer than the protein tumbling time (B10-100 ns). 116,129,132,134,137his also creates a problem for studying this process using traditional NMR techniques.However there are NMR studies confirming the existence of backbone conformational dynamics at the microsecond timescale. 138,139Current neutron scattering spectroscopy is limited to hundreds of ns and cannot reach this time scale, and current MD-simulations are usually limited to hundreds of ns.Nevertheless, a few simulations that reach millisecond times indeed identified protein relaxations in microsecond time scales (Fig. 7). 9,26In the case of bovine pancreatic trypsin inhibitor (BPTI), it has been shown that rearrangement of the sulfide bridge accompanied by proximal loop motions appear in the microsecond time range. 9However for ubiquitin, it has been shown that the jumps in the microsecond time range are between multiple conformational states, not just between two favorable states, and present the main component of the backbone dynamics. 26hus analyses of MD simulations and NMR studies suggest that the microsecond process presents the main backbone motion, its rotation, and jumps between well-defined basins in potential energy landscape of the protein structures. 9,26,138,139owever the general microscopic picture of this relaxation remains unknown.Its anomalously large dielectric amplitude provides some hints to its nature.Analysis of temperature dependence of the microsecond process in dielectric relaxation spectra suggests direct connection of this process to the glass transition in the hydrated protein system. 130,131,1331][142] In that case the microsecond relaxation presents the main structural relaxation of hydrated proteins.This connection makes it especially intriguing to understand the microscopic nature of this process.Moreover, this process appears on time scales usual for many biochemical processes and might be directly connected to the proteins function.Advanced MD-simulations, NMR and complementary scattering and spectroscopic studies might provide detailed microscopic information on the mechanism of the microsecond relaxation, and its dependence on protein structure, temperature and hydration.

III.6 Dynamics at longer times
51][52][53][54][55] These dynamics are usually ascribed to chain-like motions using worm-like, or Rouse/Zimm models with internal friction.The latter is obviously controlled by the internal dynamics discussed above.However, the connection between atomistic details of dynamics on ps-ms time scales to these longer time dynamics remains to be studied, and is out of scope of the present review.

IV. General picture of protein dynamics
Based on the presented overview of the experimental and computational results we propose a general picture of internal protein dynamics (Fig. 8) and analyze its similarity and differences with dynamics of other soft materials.First of all, there is clearly a broad spectrum of conformational fluctuations and transitions in all hydrated biomolecules.The internal dynamics under ambient conditions span over a wide time range from faster than picosecond to longer than microsecond, with atomic displacements ranging from smaller than 0.5 Å to B10 Å.As already emphasized by many researchers, 9,143,144 protein samples contain a large amount of conformational states and resemble much closer a liquid than a solid state.Therefore the use of a protein crystallographic structure as a state of active protein is misleading.
On the fast time scale, B1 ps and faster, there are small conformational fluctuations reflecting 'rattling in a cage' with usual amplitude MSD B0.2-0.4Å 2 .These motions have extremely low energy barriers of B2-5 kJ mol À1 , and are characteristic for dynamics of all soft materials. 69,71,72][147] Several relaxation processes are observed on the time scale B10-50 ps at ambient temperature (Fig. 8).They include methyl group rotations that apparently play an important role in facilitating protein dynamics and remain fast even at low hydration levels and lower temperatures. 18,19,42,84,94,95,102,148e speculate that the absence (or very low number) of methyl This journal is © The Royal Society of Chemistry 2015 groups in nucleic acids (RNA and DNA) might be one of the reasons for much stronger slowing down of their dynamics with a decrease in temperature (Fig. 5).Methyl group dynamics in proteins is similar to their dynamics in synthetic polymers and can be considered as a usual secondary relaxation.3,25,138 Translational and rotational dynamics of hydration water at ambient temperature also appears on the same time scale (Fig. 5 and 6).NMR, neutron scattering and MD-simulations all suggest that the hydration water is slowed down on average B2-4 times relative to the dynamics of bulk water. 6,102,116,1260][151][152] Moreover, most of the hydration water molecules exchange with the bulk water on a time scale much faster than one nanosecond. 41The characteristic activation energy barrier for hydration water motion is B15-18 kJ mol À1 , but it increases with a decrease in temperature, leading to the VFT-like temperature dependence of the characteristic relaxation time (Fig. 5, 6 and 9).The VFT-like behavior is typical for Soft Matter and is usually ascribed to the cooperative nature of the relaxation process that involves several molecular units in a single relaxation event.We emphasize that dynamics of hydration water and its temperature dependence are different for proteins, t-RNA and DNA (Fig. 5C), reflecting the importance of the biomolecular structure in the behavior of hydration water. 122oteins 'localized diffusion' also appears on the same time scale and presents conformational fluctuations of residues, side groups and backbone with an amplitude of B1-3 Å.Its time scale and temperature dependence (Fig. 5, 6, 8 and 9) suggest strong coupling of the localized diffusion to the dynamics of hydration water.Essentially this relaxation process presents a coupled hydration water -biomolecule relaxation, and its VFT-like behavior reflects the cooperative nature of this process.However, in contrast to translational motions of water, this process in proteins presents localized fluctuations (see, e.g.Fig. 3D).These conformational fluctuations involve the entire protein, its surface and core. 41,112,153The extremely strong stretching of its relaxation spectrum reflects significant dynamic heterogeneity apparently caused by the difference in the chemical structure and position of residues in proteins.7,130 Its characteristic time scale and temperature dependence are dictated by the behavior of hydration water that together with internal rigidity/flexibility of the biomacromolecules control the friction of the sub-nanosecond process.
The sub-nanosecond coupled protein-solvent process is associated with the dynamic transition observed as a sharp increase in the MSD of hydrated proteins at T B 220 K (Fig. 3B).Moreover, strong correlations of behavior of this process with biochemical or enzymatic activity of proteins 31,39,108,110,[154][155][156][157] Fig. 8 Proposed qualitative picture of hydrated protein dynamics at ambient temperature.Fast dynamics present small scale (B0.1-0.5 Å) conformational fluctuations on a ps time scale.Intra-basin transitions present kinds of secondary relaxations.They include a coupled protein-solvent process on the sub-nanosecond time scale (localized diffusion) with an amplitude of 1-3 Å, methyl and other side group motions on time scales from B10 ps to many ns, domain motions and other conformational jumps with amplitude extending to several Å.Inter-basin transitions that can be considered as segmental relaxation in proteins appear on a microsecond time scale and should have an amplitude of several Å.The right axis shows the corresponding energy barriers that vary from a few kJ mol À1 for fast dynamics to more than 50 kJ mol À1 for segmental relaxation.
suggest its importance for protein functions.Although this process is much faster than the characteristic time of biochemical reactions, it might present necessary precursors for slower protein dynamics, in particular, for motions directly involved in protein functions.The strength of this process reflects conformational flexibility of the protein on a ns time scale that might be critical for protein activity.In particular, the temperature dependence of the slower nanosecond process follows well the temperature dependence of the localized diffusion 48,130 (Fig. 6 and 9).
The symmetrically stretched spectral shape and localized nature of this process are similar to the behavior of a secondary relaxation in other Soft Matter.However, its non-Arrhenius temperature dependence differs from temperature variations usually observed for secondary relaxation.Moreover, usual secondary relaxation would merge with structural relaxation at time scales B10 À6 -10 À8 s, while the localized diffusion remains as a well separated process in biomacromolecules even at t B 10 À11 s.The reason for the VFT-like behavior of this process might be its coupling to the structural relaxation of hydration water.This process is certainly not the main structural relaxation responsible for the glass transition of hydrated proteins at T B 190-170 K.In fully hydrated lysozyme, this process reaches a relaxation time of B10 À4 s at T B 190-170 K, 131,133,142 instead of the usual structural relaxation of B10 2 -10 3 s at T g .9][160] Thus we will classify this subnanosecond process as a specific secondary relaxation in biomacromolecules (Fig. 8).This assignment is supported by the dominating contribution of side group motions to this process (Fig. 7).
Larger scale conformational jumps with displacements on scale B1-10 Å appear on the nanosecond time scale 6,8,10,11,28,125 (Fig. 6-9).The characteristic activation energy of this process should be B20-40 kJ mol À1 at ambient temperature, and its characteristic time scale shows VFT behavior that is similar to the VFT dependence of the sub-nanosecond process, but is slower than the latter by B10 2 -10 3 times.It depends strongly on protein hydration, even stronger than the localized diffusion.Detailed analysis of neutron scattering suggests that in the case of proteins phosphoglycerate kinase 125 and alcohol dehydrogenase 28 the nanosecond relaxation is well described by domain motions.Detection of a similar process in intrinsically disordered proteins 11,161 suggests that it can have a more general nature.It would be of great importance to understand whether this process can be described as rigid like motions and twist/bending of secondary structures, or more complicated motions have to be involved to describe the nanosecond process in various proteins.At the same time, NMR and MD-simulations studies reveal strong contribution of side group motions in this time range. 9,18,19,40,128They also present conformational jumps between different conformations of the side groups and show broad distribution of characteristic time scales, depending on the chemical structure and position of the residue in the protein.The nanosecond process is also analogous to secondary relaxation in other soft materials, and its temperature dependence most probably is controlled by the temperature dependence of the coupled solvent-protein process.
Recent dielectric relaxation studies on hydrated protein powders and MD-simulations extended to the millisecond time range both detected protein relaxations in the microsecond time range (Fig. 6-8).This process has an activation energy of the order of B50 kJ mol À1 at ambient temperature, exhibits slightly non-Arrhenius temperature dependence and has a strong dependence on hydration. 130,131,133Its possible relation to the glass transition temperature of hydrated proteins makes especially intriguing to unravel the microscopic nature of this process.It might be the main structural (segmental) relaxation of this system, analogous to structural relaxation in usual soft materials, and might present the motions directly involved in protein functions.This process can be considered as interbasin jumps (Fig. 8).Unfortunately, very little is known currently about this process.Studies on this process are complicated due to several reasons: (i) its time scale is longer than the characteristic protein tumbling time in solutions; (ii) modern neutron scattering spectroscopy cannot reach this time; and (iii) only a few simulations have been performed at long enough time to equilibrate the system and detect the process.Dielectric spectroscopy reveals an anomalously large amplitude of this process, 130,131,133 which remains a puzzle.At the same time, MD-simulations suggest that this process presents the main relaxation of the backbone. 9,26However even the length scale of this relaxation is not known, but obviously it should be larger than the length scale of the nanosecond process.3][164] Thus the microsecond process indeed can be segmental relaxation -the main structural relaxation in proteins.
The structural relaxation in folded proteins will have similarity and difference with segmental relaxation in synthetic polymers.The latter usually form random coils and on a long time scale segments can have a MSD comparable to the size of the coil, i.e.MSD B R g 2 .This is not possible for motions of residues in the protein as long as it stays folded.So, dynamics of residues will be always more localized/restricted than segmental dynamics in synthetic polymers.However, on the time scale of structural relaxation segments in synthetic polymers move only B1-2 Å, and this is possible even for residues in folded proteins.Backbone motions are the key for the segmental dynamics in polymers, and the microsecond process in proteins presents the backbone motions. 9,26,138,139Thus the proposed assignment of this process to structural relaxation of proteins seems justified.More studies are required to shed light on this probably most biologically relevant process.It is possible that this process reflects large scale changes in secondary structures and partial protein unfolding.We want to emphasize the strong similarity of the proposed picture to the analysis of atomistic millisecond simulations of BPTI (Fig. 7): one can identify sub-nanosecond and nanosecond processes as large amplitude motions of side groups and the microsecond relaxation process dominated by the backbone motions.Based on this analysis we can propose the following hierarchy of protein dynamics (Fig. 8): at a time scale of about picosecond there are small conformational fluctuations with a characteristic length scale of B0.3-0.5 Å, and low energy barriers of B2-5 kJ mol À1 .This is characteristic for the fast dynamics of all soft materials and presents some rattling in a cage formed by neighboring structural units.Methyl group dynamics appear on time scales of several ps and seem to facilitate protein dynamics even at very low temperatures or hydration levels.Methyl rotation has low energy barriers for rotation (B10-20 kJ mol À1 ) and remains fast even at lower temperatures or low hydrations (Fig. 9).Dynamics of hydration water also appear on a time scale of B10-50 ps, but it has much stronger temperature dependence than methyl groups (Fig. 9).Dynamics of hydration water is coupled to protein's subnanosecond process, and the latter presents a kind of diffusion localized in a harmonic potential with the length scale limited to 1-3 Å (Fig. 8).The friction coefficient of this motion is defined by hydration water and internal protein friction (flexibility).As a result, the characteristic time scale of the localized diffusion depends strongly on the hydration level and follows temperature dependence of hydration water dynamics (Fig. 9).This motion might be a precursor for conformational jumps that occur on a nanosecond time scale and has an amplitude of B2-8 Å (Fig. 8).It has characteristic energy barriers B20-40 kJ mol À1 and presents a kind of intrabasin conformational jumps.Its temperature variations seem to be defined by the temperature dependence of the sub-nanosecond process (Fig. 9).But it depends on hydration even stronger than the latter. 48,130,131Both these processes, localized diffusion and nanosecond relaxation, involve mostly side group motions in the protein core and surface (Fig. 7).Thus they might be considered as secondary relaxations in biomolecules.Finally, the inter-basin relaxation appears at the microsecond time scale (Fig. 8).This microsecond process involves mostly the backbone rotation 9,26 (Fig. 7) and apparently can be considered as segmental (structural) relaxation of the protein.It has a characteristic energy barrier 450 kJ mol À1 and depends strongly on hydration. 130,131This process might be the most directly related to protein functions, folding and other biologically relevant processes, including glass transition.
16][17][51][52][53][54][55] These studies include different optical and pulling experiments with characteristic time scales usually longer than ms and displacements in the range 10-100 Å. Analysis of these experiments usually can be done using coarse-grained models where atomistic details of studied processes are hidden in friction terms and other phenomenological parameters.For example, authors of ref. 15 analyzed retraction of an extended DNA chain on the millisecond time scale using a classical worm-like model.This process reflects the chain-like motion of DNA and should be significantly slower than segmental (structural) relaxation.The relaxation processes discussed in this review at ps-ms time scales define the internal friction coefficient for such worm-like motions.The importance of internal conformational dynamics for coarse-grained models of proteins (e.g.colloidal-like approach, Zimm-like models) has also been emphasized in many papers. 11,17At present, these two worlds of protein studies, atomistically detailed analysis presented here and coarse-grained single molecular studies, do not overlap much.It would be extremely important to bridge these two approaches, and in this way broaden and deepen our understanding of the biomolecular dynamics.

V. Current challenges and future perspectives
Despite significant efforts towards studies on biomolecular dynamics, we are still missing a general atomistic picture of internal protein dynamics, classification of different relaxation processes and their connection to protein function and activity.There is no accepted theory or model of the protein dynamics even on a qualitative phenomenological level.In this review we attempt to formulate a general atomistic picture of biomolecular dynamics with clear separation of various relaxation processes.The presented picture (Fig. 8) is consistent with an earlier review by Henzler-Wildman and Kern, 4 but provides more specific classification and more microscopic details based on neutron scattering, MD-simulations and dielectric spectroscopy studies.
The presented review did not touch the important question of heterogeneity in protein dynamics caused by difference in rigidity of secondary structures, flexibility of loops and disordered regions of the protein structure.Based on hydrogen exchange rate studies, ideas of knots in protein structures with strongly suppressed dynamics have been proposed. 165MD simulations also revealed a smaller amplitude of motions in the protein core relative to the motions on the protein surface. 41,112The intrinsic dynamic heterogeneity might be a key to protein stability and unfolding, enzymatic activity and other protein functions.Unfortunately, this topic of protein research remains in its early stage (see e.g., ref. 128), and no systematic picture can be drawn at this time.Creating a general picture of protein dynamics and finding details of intrinsic heterogeneities is one of the challenges in future research on protein dynamics.
In this review we did not discuss the role of secondary and tertiary structure in protein dynamics, flexibility and stability.7][168] Analysis of their dynamics demonstrated that they have higher flexibility and large amplitude of motions on a picosecond-nanosecond time scale than regular globular proteins. 36,161,169,170Whether intrinsically disordered proteins exhibit dynamics significantly different from the dynamics of globular proteins considered here remains to be studied.
We also did not discuss the role of solvents in biomolecular dynamics.This is another important topic that is directly related to the field of bio-preservation.46][147]171 But detailed understanding of the role of the solvents in protein dynamics, activity and stability remains a great challenge that affects many life science and biotechnological applications.Developing a general concept of protein dynamics is one of the keys to address this challenge.
One of the important challenges is to reach longer time scales (milliseconds) with MD-simulations and with experimental techniques that provide microscopic information, such as neutron or X-ray scattering.While the former is feasible with developments of more powerful computers, faster algorithms and proper coarsegraining, the latter does not look promising in the near future.However, development of intense X-ray sources with short pulses opens a new possibility for analysis of kinetics of conformational changes in proteins. 172This might help in analysis of dynamics of biological macromolecules on time scales longer than microsecond.Also detailed dielectric and NMR studies combined with MD-simulations might provide the required information on protein dynamics on longer time scales, once the millisecond time scale becomes more accessible for atomistic MD.In addition, combining these atomistic studies with optical and single molecular studies might provide a breakthrough in understanding the connection between conformational dynamics and enzymatic activity of proteins.
The main challenge, in our view, remains the development of analytical models describing protein dynamics.They most probably will be based on traditional Langevin dynamics with frictions estimated from experimental data or atomistic MD-simulations.These models should provide a description of atomic motions, motions of secondary structures and loops, and even partial folding/unfolding of secondary structures due to thermodynamically equilibrium fluctuations.We hope that the presented review and the proposed general picture will help to develop such a model of the internal protein dynamics that will describe the sub-nanosecond localized diffusion, nanosecond conformational jumps and microsecond inter-basin transitions in a unified way.

Sheila
Khodadadi received her PhD in Polymer Science in 2009 from the University of Akron, USA.In 2010, she was awarded two years NRC Postdoctoral Research Associateship to conduct her research at NIST Center for Neutron Research, National Institute of Standards and Technology, USA.Currently, she is a researcher at Delft University of Technology in the Netherlands.Her main research interest is understanding protein-solvent dynamics and interactions in pharmaceutically relevant conditions.A. P. Sokolov Alexei Sokolov received an MS in Physics in 1981 from the Novosibirsk State University, and a PhD in Physics in 1986 from the Russian Academy of Sciences.He worked several years in Germany before joining faculty at the University of Akron, USA, in 1998.In 2009 he accepted Governor's Chair position at the University of Tennessee and Oak Ridge National Laboratory, where he leads the Soft Matter group.His current research interest focuses primarily on dynamics of soft materials, including topics of glass transition, polymer dynamics and dynamics of biological macromolecules, composite materials and materials for energy applications.
Fig. 1 presents an overview of the general picture This journal is © The Royal Society of Chemistry 2015

Fig. 2
Fig. 2 Dynamic structure factor S(Q,n) of lysozyme at (A) different hydration levels 39 (h refers to grams of water, D 2 O, per one gram of protein), and (B, C) different temperatures and solvents (DL -dry lysozyme; WL -hydrated lysozyme;LG -lysozyme in glycerol; LT -lysozyme in trehalose).39,90(A) As the hydration level increases the Boson peak shifts to higher frequencies and the quasielastic contribution coming mainly from the fast relaxation process strongly increases.(B) At low T = 150 K the fast relaxation process is strongly suppressed in LG.While as the temperature is increased to ambient T B 295 K. (C) WL has the most contribution from the fast relaxation process.This behavior of the protein's fast relaxation follows the behavior of the solvent.79,89

Fig. 3 (
Fig. 3 (A) Dynamic structure factor S(Q, E) of lysozyme, tRNA and DNA.Data for lysozyme and tRNA are obtained from ref. 92, for DNA from ref. 99.The relaxation process appears as a quasielastic scattering, i.e. extra broadening in comparison to the resolution function of the spectrometer (line).Lysozyme spectra exhibit relaxation contribution even in the dry state and at T = 200 K, while no quasielastic scattering is observed for dry tRNA and DNA, or for hydrated tRNA and DNA below 210 K.This broadening for lysozyme is related to the methyl group rotation.(B) Mean squared displacements demonstrate the onset of methyl group rotation at T B 100 K for lysozyme as marked by the arrow.This onset is absent in MSD of tRNA.A sharp rise of MSD in hydrated biomacromolecules at T above B220 K is usually called the dynamic transition.(C) Methyl groups in different classes of biological macromolecules are presented in numbers of methyl groups per molecular weight.94Proteins have a significant number of methyl groups in comparison to DNA (one out of four nucleic acids has a single methyl group) and RNA (no methyl in primary structure of RNA).(D) An example of MSD of hydrogen atoms in lysozyme involved in the localized diffusion motion as a function of time obtained from MD-simulations.6It demonstrates fast dynamics and caging (plateau) at times shorter than B3 ps, followed by the localized diffusion with an almost logarithmic increase in MSD up to B200 ps.MSD saturates at longer times indicating localized motion.

Fig. 4 (
Fig. 4 (A) Neutron scattering susceptibility spectra of the localized motion in hydrated lysozyme (h B 0.35) at different temperature after correction for methyl group rotation. 106(B) Neutron scattering susceptibility of the localized motion in hydrated t-RNA (h B 0.65) at different temperatures. 122Symbols are experimental data and lines are the fit to Cole-Cole distribution functions.Names at the top present different neutron scattering spectrometers used in these studies.

Fig. 6
Fig. 6 Temperature dependence of the relaxation times observed by dielectric spectroscopy for hydrated lysozyme (h B 0.4; filled symbols) and for hydrated myoglobin (h B 0.33; half-filled symbols) powders.Neutron scattering data for hydrated lysozyme (h B 0.4) are shown as crossed circles (data taken from ref.48 and 130).Circles are localized diffusion (sub-nanosecond coupled protein-solvent motion), triangles are domain motions (nanosecond process) and squares are segmental relaxation (microsecond dielectric process).For comparison the dielectric relaxation times in protein solutions at 298 K are also included: RNase A (open triangles 116 and star 137 ), lysozyme 137 (open circle).Characteristic t of the nanosecond relaxation estimated from neutron spin echo measurements are also shown: alcohol dehydrogenase 28 (open diamond); PGK bound (open square) and unbound to the substrate (open pentagon) at 283 K (data from ref. 125).

Fig. 7
Fig. 7 Dynamical content of the P 2 internal correlation function and its decomposition into side chain and backbone contribution for bovine pancreatic trypsin inhibitor (BPTI) obtained from millisecond simulations. 9It demonstrates how dynamics of different protein parts spread over a wide range of time scale.Most important, it clearly indicates existence of three relaxation modes: (i) sub-nanosecond; (ii) nanosecond and (iii) microsecond relaxations.The latter is the only one dominated by the backbone motions.Adopted with permission from ref. 9.

Fig. 9
Fig. 9 Schematic of temperature dependence of different internal motions of protein including: (i) main structural (segmental) relaxation reaching t B 100 s at the glass transition temperature T g of hydrated proteins; (ii) domain motions, conformational jumps with nanosecond relaxation times at room temperature; (iii) coupled protein-solvent motions, presents ''localized diffusion'' of residues and side group motions with a relaxation time of tens of picosecond at room temperature; (iv) methyl group rotation with Arrhenius temperature dependence; and (v) fast picosecond relaxation with a rather weak temperature dependence of the characteristic relaxation time.Processes (ii)-(iv) are apparently kinds of secondary relaxation in proteins.