A benchmark for protein dynamics: Ribonuclease A measured by neutron scattering in a large wavevector-energy transfer range

Determination of in cellulo proteome molecular dynamics in different halophilic Archaea

While biophysical studies have unravelled properties of specific proteins in vitro, characterizing globally their native state within the cell remains a challenge. In particular, protein adaptation to harsh intracellular physical and chemical conditions is poorly understood. Extremophiles, which thrive in severe environments, are good models for the study of such adaptation. Five haloarchaeal species, isolated from hypersaline environments, were used to assess correlations between intracellular salt concentrations and molecular dynamics properties. In cellulo protein stability was measured using nano differential scanning fluorimetry, and neutron spectrometry was used to determine molecular dynamics resilience and global flexibility. It was found that high intracellular accumulation of Mg2+ and low intracellular accumulation of K+ were correlated with higher stability and resilience. Sequence traits associated with mean proteome halophilicity, such as decreased hydrophobicity and increased acidity, weighted by the relative abundance of each protein, were also correlated with stability and resilience. Haloferax mediterranei, however, was found to be an exception as its proteome showed the highest in cellulo molecular stability and resilience associated with fewest sequence traits related to halophilicity, highlighting the significance of the intracellular salt environment in determining proteome biophysical properties.

Hydration-induced dynamical changes in lyophilised and weakly hydrated apoferritin: insights from molecular dynamics simulation

The dynamics and functionality of proteins are significantly influenced by their interaction with water. For lyophilised (i.e. h ≤ 0.05 where h = g of H2O per g of protein) and weakly hydrated systems (i.e. h ≤ 0.38) hydration generally enhances protein mobility above the so-called 'dynamical transition' temperature (Td > 220 K). However, water-induced mobility hindrance at low temperatures (T < 175 K) has been reported in various proteins of varying secondary structure; namely green fluorescent protein (GFP), pig liver esterase, lysozyme, ribonuclease A (RNAse A) and apoferritin. By focussing on the dynamic behaviour of the apoferritin molecule, this study proposes mechanisms driving these hydration-induced mobility changes, particularly the less understood hindrance at low temperatures. Using atomistic molecular dynamics (MD) simulations of horse spleen apoferritin in the lyophilised (h = 0.05) and weakly hydrated (h = 0.31) states, we report here the impact of water on protein dynamics as a function of temperature. Through residue-specific mean squared displacement (MSD), radial distribution function (RDF), solvent accessible surface area (SASA), local hydration degree and hydrogen bonding analyses, we demonstrate that while water proximity directly correlates with mobility enhancement at high temperatures, the hydration-induced mobility reduction observed at temperatures below 175 K is primarily propagated through the protein backbone.

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

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

Inelastic neutron scattering and spectroscopy methods to characterize dynamics in colloidal and soft matter systems

Colloidal systems and soft materials are well suited to neutron scattering, and the community has readily adopted elastic scattering techniques to investigate their structure. Due to their unique properties, neutrons may also be used to characterize the dynamics of soft materials over a wide range of length and time scales in situ. Both static structures and an understanding of how molecules move about their equilibrium positions is essential if we are to deliver on the promise of rationally designing soft materials. In this review we introduce the basics of neutron spectroscopy and explore the ways in which inelastic neutron scattering can be used to study colloidal and soft materials. Illustrative examples are chosen that highlight the phenomena suitable for investigation using this suite of techniques.

Functional-State Dependence of Picosecond Protein Dynamics

We examine temperature-dependent picosecond dynamics of two benchmarking proteins lysozyme and cytochrome c using temperature-dependent terahertz permittivity measurements. We find that a double Arrhenius temperature dependence with activation energies E₁ ∼ 0.1 kJ/mol and E₂ ∼ 10 kJ/mol fits the folded and ligand-free state response. The higher activation energy is consistent with the so-called protein dynamical transition associated with beta relaxations at the solvent-protein interface. The lower activation energy is consistent with correlated structural motions. When the structure is removed by denaturing, the lower-activation-energy process is no longer present. Additionally, the lower-activation-energy process is diminished with ligand binding but not for changes in the internal oxidation state. We suggest that the lower-energy activation process is associated with collective structural motions that are no longer accessible with denaturing or binding.

Temperature and salt controlled tuning of protein clusters

The formation of molecular assemblies in protein solutions is of strong interest both from a fundamental viewpoint and for biomedical applications. While ordered and desired protein assemblies are indispensable for some biological functions, undesired protein condensation can induce serious diseases. As a common cofactor, the presence of salt ions is essential for some biological processes involving proteins, and in aqueous suspensions of proteins can also give rise to complex phase diagrams including homogeneous solutions, large aggregates, and dissolution regimes. Here, we systematically study the cluster formation approaching the phase separation in aqueous solutions of the globular protein BSA as a function of temperature (T), the protein concentration (c_p) and the concentrations of the trivalent salts YCl₃ and LaCl₃ (c_s). As an important complement to structural, i.e. time-averaged, techniques we employ a dynamical technique that can detect clusters even when they are transient on the order of a few nanoseconds. By employing incoherent neutron spectroscopy, we unambiguously determine the short-time self-diffusion of the protein clusters depending on c_p, c_s and T. We determine the cluster size in terms of effective hydrodynamic radii as manifested by the cluster center-of-mass diffusion coefficients D. For both salts, we find a simple functional form D(c_p, c_s, T) in the parameter range explored. The calculated inter-particle attraction strength, determined from the microscopic and short-time diffusive properties of the samples, increases with salt concentration and temperature in the regime investigated and can be linked to the macroscopic behavior of the samples.

The Dynamics of Hydrated Proteins Are the Same as Those of Highly Asymmetric Mixtures of Two Glass-Formers

Customarily, the studies of dynamics of hydrated proteins are focused on the fast hydration water ν-relaxation, the slow structural α-relaxation responsible for a single glass transition, and the protein dynamic transition (PDT). Guided by the analogy with the dynamics of highly asymmetric mixtures of molecular glass-formers, we explore the possibility that the dynamics of hydrated proteins are richer than presently known. By providing neutron scattering, dielectric relaxation, calorimetry, and deuteron NMR data in two hydrated globular proteins, myoglobin and BSA, and the fibrous elastin, we show the presence of two structural α-relaxations, α1 and α2, and the hydration water ν-relaxation, all coupled together with interconnecting properties. There are two glass transition temperatures T _g ^α1and T _g ^α2 corresponding to vitrification of the α1 and α2 processes. Relaxation time τ_α2(T) of the α2-relaxation changes its Arrhenius temperature dependence to super-Arrhenius on crossing T _g ^α1 from below. The ν-relaxation responds to the two vitrifications by changing the T-dependence of its relaxation time τ_ν(T) on crossing consecutively T _g ^α2 and T _g ^α1. It generates the PDT at T _d where τ_ν(T _d) matches about five times the experimental instrument timescale τ_exp, provided that T _d > T _g ^α1. This condition is satisfied by the hydrated globular proteins considered in this paper, and the ν-relaxation is in the liquid state with τ_ν(T) having the super-Arrhenius temperature dependence. However, if T _d < T _g ^α1, the ν-relaxation fails to generate the PDT because it is in the glassy state and τ_ν(T) has Arrhenius T-dependence, as in the case of hydrated elastin. Overall, the dynamics of hydrated proteins are the same as the dynamics of highly asymmetric mixtures of glass-formers. The results from this study have expanded the knowledge of the dynamic processes and their properties in hydrated proteins, and impact on research in this area is expected.

Dynamics of proteins in solution

The dynamics of proteins in solution includes a variety of processes, such as backbone and side-chain fluctuations, interdomain motions, as well as global rotational and translational (i.e. center of mass) diffusion. Since protein dynamics is related to protein function and essential transport processes, a detailed mechanistic understanding and monitoring of protein dynamics in solution is highly desirable. The hierarchical character of protein dynamics requires experimental tools addressing a broad range of time- and length scales. We discuss how different techniques contribute to a comprehensive picture of protein dynamics, and focus in particular on results from neutron spectroscopy. We outline the underlying principles and review available instrumentation as well as related analysis frameworks.

Motional displacements in proteins: The origin of wave-vector-dependent values

The average mean-square displacement, 〈r(2)〉, of H atoms in a protein is frequently determined using incoherent neutron-scattering experiments. 〈r(2)〉 is obtained from the observed elastic incoherent dynamic structure factor, S(i)(Q,ω=0), assuming the form S(i)(Q,ω=0) =exp(-Q(2)〈r(2)〉/3). This is often referred to as the Gaussian approximation (GA) to S(i)(Q,ω=0). 〈r(2)〉 obtained in this way depends on the value of the wave vector, Q considered. Equivalently, the observed S(i)(Q,ω=0) deviates from the GA. We investigate the origin of the Q dependence of 〈r(2)〉 by evaluating the scattering functions in different approximations using molecular dynamics (MD) simulation of the protein lysozyme. We find that keeping only the Gaussian term in a cumulant expansion of S(Q,ω) is an accurate approximation and is not the origin of the Q dependence of 〈r(2)〉. This is demonstrated by showing that the term beyond the Gaussian is negligible and that the GA is valid for an individual atom in the protein. Rather, the Q dependence (deviation from the GA) arises from the dynamical heterogeneity of the H in the protein. Specifically it arises from representing, in the analysis of data, this diverse dynamics by a single average scattering center that has a single, average 〈r(2)〉. The observed Q dependence of 〈r(2)〉 can be used to provide information on the dynamical heterogeneity in proteins.

Elastic and conformational softness of a globular protein

Flexibility, or softness, is crucial for protein function and consists of a conformational component, involving jumps between potential wells, and an elastic component, involving fluctuations within the wells. Combining molecular dynamics simulation with incoherent neutron scattering and light scattering measurements on green fluorescent protein, we reveal a relationship between the intrawell fluctuations and elastic moduli of the protein. This finding leads to a simple means of experimentally separating the conformational from the elastic atomic displacements.

Dynamical coupling of intrinsically disordered proteins and their hydration water: comparison with folded soluble and membrane proteins

Hydration water is vital for various macromolecular biological activities, such as specific ligand recognition, enzyme activity, response to receptor binding, and energy transduction. Without hydration water, proteins would not fold correctly and would lack the conformational flexibility that animates their three-dimensional structures. Motions in globular, soluble proteins are thought to be governed to a certain extent by hydration-water dynamics, yet it is not known whether this relationship holds true for other protein classes in general and whether, in turn, the structural nature of a protein also influences water motions. Here, we provide insight into the coupling between hydration-water dynamics and atomic motions in intrinsically disordered proteins (IDP), a largely unexplored class of proteins that, in contrast to folded proteins, lack a well-defined three-dimensional structure. We investigated the human IDP tau, which is involved in the pathogenic processes accompanying Alzheimer disease. Combining neutron scattering and protein perdeuteration, we found similar atomic mean-square displacements over a large temperature range for the tau protein and its hydration water, indicating intimate coupling between them. This is in contrast to the behavior of folded proteins of similar molecular weight, such as the globular, soluble maltose-binding protein and the membrane protein bacteriorhodopsin, which display moderate to weak coupling, respectively. The extracted mean square displacements also reveal a greater motional flexibility of IDP compared with globular, folded proteins and more restricted water motions on the IDP surface. The results provide evidence that protein and hydration-water motions mutually affect and shape each other, and that there is a gradient of coupling across different protein classes that may play a functional role in macromolecular activity in a cellular context.

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.

Dynamic behavior of oligomeric inorganic pyrophosphatase explored by quasielastic neutron scattering

The purpose of this investigation is to determine whether a large oligomeric protein, inorganic pyrophosphatase (IPPase) from Thermococcus thioreducens with quaternary structural complexity, would have distinguishable dynamic characteristics compared to those of the small simple monomeric model protein, lysozyme. In this study, the β-relaxational dynamics of the two proteins, IPPase and lysozyme, are compared in the 10 ps to 0.5 ns time interval using quasi-elastic neutron scattering (QENS). Both of the protein dynamics show a characteristic logarithmic-like decay in the intermediate scattering function (ISF) of the hydrogen atoms. Distinguishable dynamical behavior found between two proteins reveals local flexibility and conformational substates unique to oligomeric structures. Moreover, the temperature dependence of the mean square displacement (MSD) of the hydrogen atoms in protein molecules, which is a traditional way to determine the "softness" of the protein molecule, is measured and shows no difference for the two proteins.

Surface hydration amplifies single-well protein atom diffusion propagating into the macromolecular core

The effect of surface hydration water on internal protein motion is of fundamental interest in molecular biophysics. Here, by decomposing the picosecond to nanosecond atomic motion in molecular dynamics simulations of lysozyme at different hydration levels into three components--localized single-well diffusion, methyl group rotation, and nonmethyl jumps--we show that the effect of surface hydration is mainly to increase the volume of the localized single-well diffusion. These diffusive motions are coupled in such a way that the hydration effect propagates from the protein surface into the dry core.

Use of neutrons reveals the dynamics of cell surface glycosaminoglycans

Cell surface glycosaminoglycans (GAG), such as heparan sulfate (HS) and heparin, are key multifunctional cell regulators, which are involved in numerous molecular events associated with tumor growth, metastasis, pathogen attachment, and immune response. GAG dynamically bind and regulate the activities of many signaling proteins such as growth factors, chemokines, and cytokines. GAG-binding interactions with proteins rely on the coupling between the geometry, flexibility, and rigidity of the polysaccharide chain. Understanding GAG dynamics at the molecular level can therefore provide fundamental insights into GAG function in a cellular context. Elastic incoherent neutron scattering is a powerful tool for the exploration of fast molecular motions in biological macromolecules. Recently, the technique was used to evaluate HS flexibility and rigidity on different timescales between the picosecond (ps) and the nanosecond (ns). Here, neutron spectroscopy experimental procedures are presented, with emphasis on the practical details necessary to prepare samples, run neutron scattering experiments, and extract the dynamics parameters from the data.

Study of Protein Dynamics vs. Amyloid Formation

Protein fibril formation has been often associated to manifestation of serious and devastating amyloyd diseases, including Alzheimer and BSA. The main mechanism for the formation of amyloid fibrils is the accumulation of protein aggregates in body´s organs.In this paper, we try to compare the dynamical behaviour of two amyloidogenic proteins, the Insulin and the Myoglobin. Insulin has been chosen for its pharmacological extensive use in diabete´s therapy, while Myoglobin is used as control, since its dynamics is now largely known. The investigation has been performed through incoherent elastic neutron scattering over a wide temperature range. Our results suggest an enhanced stiffness of Insulin with respect to Myoglobin.

Collective dynamics of protein hydration water by brillouin neutron spectroscopy

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

Dynamics of hydration water in deuterated purple membranes explored by neutron scattering

The function and dynamics of proteins depend on their direct environment, and much evidence has pointed to a strong coupling between water and protein motions. Recently however, neutron scattering measurements on deuterated and natural-abundance purple membrane (PM), hydrated in H2O and D2O, respectively, revealed that membrane and water motions on the ns-ps time scale are not directly coupled below 260 K (Wood et al. in Proc Natl Acad Sci USA 104:18049-18054, 2007). In the initial study, samples with a high level of hydration were measured. Here, we have measured the dynamics of PM and water separately, at a low-hydration level corresponding to the first layer of hydration water only. As in the case of the higher hydration samples previously studied, the dynamics of PM and water display different temperature dependencies, with a transition in the hydration water at 200 K not triggering a transition in the membrane at the same temperature. Furthermore, neutron diffraction experiments were carried out to monitor the lamellar spacing of a flash-cooled deuterated PM stack hydrated in H2O as a function of temperature. At 200 K, a sudden decrease in lamellar spacing indicated the onset of long-range translational water diffusion in the second hydration layer as has already been observed on flash-cooled natural-abundance PM stacks hydrated in D2O (Weik et al. in J Mol Biol 275:632-634, 2005), excluding thus a notable isotope effect. Our results reinforce the notion that membrane-protein dynamics may be less strongly coupled to hydration water motions than the dynamics of soluble proteins.