Effect of conformational states on protein dynamical transition

Incoherent Neutron Scattering and Terahertz Time-Domain Spectroscopy on Protein and Hydration Water

Incoherent inelastic and quasi-elastic neutron scattering (INS) and terahertz time-domain spectroscopy (THz-TDS) are spectroscopy methods that directly detect molecular dynamics, with an overlap in the measured energy regions of each method. Due to the different characteristics of their probes (i.e., neutron and light), the information obtained and the sample conditions suitable for each method differ. In this review, we introduce the differences in the quantum beam properties of the two methods and their associated advantages and disadvantages in molecular spectroscopy. Neutrons are scattered via interaction with nuclei; one characteristic of neutron scattering is a large incoherent scattering cross-section of a hydrogen atom. INS records the auto-correlation functions of atomic positions. By using the difference in neutron scattering cross-sections of isotopes in multi-component systems, some molecules can be selectively observed. In contrast, THz-TDS observes the cross-correlation function of dipole moments. In water-containing biomolecular samples, the absorption of water molecules is particularly large. While INS requires large-scale experimental facilities, such as accelerators and nuclear reactors, THz-TDS can be performed at the laboratory level. In the analysis of water molecule dynamics, INS is primarily sensitive to translational diffusion motion, while THz-TDS observes rotational motion in the spectrum. The two techniques are complementary in many respects, and a combination of the two is very useful in analyzing the dynamics of biomolecules and hydration water.

Universality and Structural Implications of the Boson Peak in Proteins

The softness and rigidity of proteins are reflected in the structural dynamics, which are in turn affected by the environment. The characteristic low-frequency vibrational spectrum of a protein, known as boson peak, is an indication of the structural rigidity of the protein at a cryogenic temperature or dehydrated conditions. In this article, the effect of hydration, temperature, and pressure on the boson peak and volumetric properties of a globular protein are evaluated by using inelastic neutron scattering and molecular dynamics simulation. Hydration, pressurization, and cooling shift the boson peak position to higher energy and depress the peak intensity and decreases the protein and cavity volumes. We found the correlation between the boson peak and cavity volume in a protein. A decrease of cavity volume means the increase of rigidity, which is the origin of the boson peak shift. Boson peak is the universal property of a protein, which is rationalized by the correlation.

Determination of Dynamical Heterogeneity from Dynamic Neutron Scattering of Proteins

Motional displacements of hydrogen (H) in proteins can be measured using incoherent neutron-scattering methods. These displacements can also be calculated numerically using data from molecular dynamics simulations. An enormous amount of data on the average mean-square motional displacement (MSD) of H as a function of protein temperature, hydration, and other conditions has been collected. H resides in a wide spectrum of sites in a protein. Some H are tightly bound to molecular chains, and the H motion is dictated by that of the chain. Other H are quite independent. As a result, there is a distribution of motions and MSDs of H within a protein that is denoted dynamical heterogeneity. The goal of this paper is to incorporate a distribution of MSDs into models of the H incoherent intermediate scattering function, I(Q,t), that is calculated and observed. The aim is to contribute information on the distribution as well as on the average MSD from comparison of the models with simulations and experiment. For example, we find that simulations of I(Q,t) in lysozyme are well reproduced if the distribution of MSDs is bimodal with two broad peaks rather than a single broad peak.

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.

Effect of Temperature and Hydration Level on Purple Membrane Dynamics Studied Using Broadband Dielectric Spectroscopy from Sub-GHz to THz Regions

To investigate the effects of temperature and hydration on the dynamics of purple membrane (PM), we measured the broadband complex dielectric spectra from 0.5 GHz to 2.3 THz using a vector network analyzer and terahertz time-domain spectroscopy from 233 to 293 K. In the lower temperature region down to 83 K, the complex dielectric spectra in the THz region were also obtained. The complex dielectric spectra were analyzed through curve fitting using several model functions. We found that the hydrated states of one relaxational mode, which was assigned as the coupled motion of water molecules with the PM surface, began to overlap with the THz region at approximately 230 K. On the other hand, the relaxational mode was not observed for the dehydrated state. On the basis of this result, we conclude that the protein-dynamical-transition-like behavior in the THz region is due to the onset of the overlap of the relaxational mode with the THz region. Temperature hysteresis was observed in the dielectric spectrum at 263 K when the hydration level was high. It is suggested that the hydration water behaves similarly to supercooled liquid at that temperature. The third hydration layer may be partly formed to observe such a phenomenon. We also found that the relaxation time is slower than that of a globular protein, lysozyme, and the microscopic environment in the vicinity of the PM surface is suggested to be more heterogeneous than lysozyme. It is proposed that the spectral overlap of the relaxational mode and the low-frequency vibrational mode is necessary for the large conformational change of protein.

The role of momentum transfer during incoherent neutron scattering is explained by the energy landscape model

We recently introduced a model of incoherent quasielastic neutron scattering (QENS) that treats the neutrons as wave packets of finite length and the protein as a random walker in the free energy landscape. We call the model ELM for "energy landscape model." In ELM, the interaction of the wave packet with a proton in a protein provides the dynamic information. During the scattering event, the momentum [Formula: see text] is transferred by the wave packet to the struck proton and its moiety, exerting the force [Formula: see text] The resultant energy [Formula: see text] is stored elastically and returned to the neutron as it exits. The energy is given by [Formula: see text], where [Formula: see text] is the ambient temperature and [Formula: see text] ([Formula: see text] 91 K Å) is a new elastobaric coefficient. Experiments yield the scattering intensity (dynamic structure factor) [Formula: see text] as a function of [Formula: see text] and [Formula: see text] To test our model, we use published data on proteins where only thermal vibrations are active. ELM competes with the currently accepted theory, here called the spatial motion model (SMM), which explains [Formula: see text] by motions in real space. ELM is superior to SMM: It can explain the experimental angular and temperature dependence, whereas SMM cannot do so.

Quasielastic neutron scattering in biology: Theory and applications

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

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

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

Dynamics of 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.

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.

Dynamics and rigidity in an intrinsically disordered protein, β-casein

The emergence of intrinsically disordered proteins (IDPs) as a recognized structural class has forced the community to confront a new paradigm of structure, dynamics, and mechanical properties for proteins. We present novel data on the similarities and differences in the dynamics and nanomechanical properties of IDPs and other biomacromolecules on the picosecond time scale. An IDP, β-casein (CAS), has been studied in a calcium bound and unbound state using neutron and light scattering techniques. We show that CAS partially folds and stiffens upon calcium binding, but in the unfolded state, it is softer than folded proteins such as green fluorescence protein (GFP). We also see that some localized diffusive motions in CAS have a larger amplitude than in GFP at this time scale but are still smaller than those observed in tRNA. In spite of these differences, CAS dynamics are consistent with the classes of motions seen in folded protein on this time scale.

Long-time mean-square displacements in proteins

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

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.

Intrinsic mean-square displacements in proteins

The thermal mean-square displacement (MSD) of hydrogen in proteins and its associated hydration water is measured by neutron scattering experiments and used an indicator of protein function. The observed MSD as currently determined depends on the energy resolution width of the neutron scattering instrument employed. We propose a method for obtaining the intrinsic MSD of H in the proteins, one that is independent of the instrument resolution width. The intrinsic MSD is defined as the infinite time value of (r(2)) that appears in the Debye-Waller factor. The method consists of fitting a model to the resolution broadened elastic incoherent structure factor or to the resolution dependent MSD. The model contains the intrinsic MSD, the instrument resolution width, and a rate constant characterizing the motions of H in the protein. The method is illustrated by obtaining the intrinsic MSD (r(2)) of heparan sulphate (HS-0.4), ribonuclease A, and staphysloccal nuclase (SNase) from data in the literature.

Dynamical transition in proteins and non-Gaussian behavior of low-frequency modes in self-consistent normal mode analysis

Self-consistent normal mode analysis (SCNMA) is applied to heme c type cytochrome f to study temperature-dependent protein motion. Classical normal mode analysis assumes harmonic behavior and the protein mean-square displacement has a linear dependence on temperature. This is only consistent with low-temperature experimental results. To connect the protein vibrational motions between low and physiological temperatures, we have incorporated a fitted set of anharmonic potentials into SCNMA. In addition, quantum harmonic-oscillator theory has been used to calculate the displacement distribution for individual vibrational modes. We find that the modes involving soft bonds exhibit significant non-Gaussian dynamics at physiological temperature, which suggests that it may be the cause of the non-Gaussian behavior of the protein motions probed by elastic incoherent neutron scattering. The combined theory displays a dynamical transition caused by the softening of few "torsional" modes in the low-frequency regime ( <50 cm(-1) or <6 meV or >0.6 ps). These modes change from Gaussian to a classical distribution upon heating. Our theory provides an alternative way to understand the microscopic origin of the protein dynamical transition.