Picosecond dynamical changes on denaturation of yeast phosphoglycerate kinase revealed by quasielastic neutron scattering

Quasielastic neutron scattering experiments performed on yeast phosphoglycerate kinase in the native form and denatured in 1.5 M guanidinium chloride reveal a change in the fast (picosecond time scale) diffusive internal dynamics of the protein. The momentum and energy transfer dependences of the scattering for both states are fitted by an analytical model in which, on the experimentally accessible picosecond time scale and angstrom length scale, the dynamics of a fraction of the nonexchangeable hydrogens in the protein is described as a superposition of vibrations with uniform diffusion in a sphere, the rest of the hydrogens undergoing only vibrational motion. The fraction diffusing changes, from ≈60% in the native protein to ≈82% in the denatured protein. The radius of the sphere also changes slightly, from ≈1.8 Å in the native protein to ≈2.2 Å in the denatured protein. Possible implications of these results for the general protein folding problem are discussed. Proteins 28:380–387, 1997 © 1997 Wiley-Liss, Inc.     

A view of dynamics changes in the molten globule-native folding step by quasielastic neutron scattering

In order to understand the changes in protein dynamics that occur in the final stages of protein folding, we have used neutron scattering to probe the differences between a protein in its folded state and the molten globule states. The internal dynamics of bovine alpha-lactalbumin (BLA) and its molten globules (MBLA) have been examined using incoherent, quasielastic neutron scattering (IQNS). The IQNS results show length scale dependent, pico-second dynamics changes on length scales from 3.3 to 60 A studied. On shorter-length scales, the non-exchangeable protons undergo jump motions over potential barriers, as those involved in side-chain rotamer changes. The mean potential barrier to local jump motions is higher in BLA than in MBLA, as might be expected. On longer length scales, the protons undergo spatially restricted diffusive motions with the diffusive motions being more restricted in BLA than in MBLA. Both BLA and MBLA have similar mean square amplitudes of high frequency motions comparable to the chemical bond vibrational motions. Bond vibrational motions thus do not change significantly upon folding. Interestingly, the quasielastic scattering intensities show pronounced maxima for both BLA and MBLA, suggesting that "clusters" of atoms are moving collectively within the proteins on picosecond time scales. The correlation length, or "the cluster size", of such atom clusters moving collectively is dramatically reduced in the molten globules with the correlation length being 6.9 A in MBLA shorter than that of 18 A in BLA. Such collective motions may be important for the stability of the folded state, and may influence the protein folding pathways from the molten globules.     

Picosecond dynamics in haemoglobin from different species: A quasielastic neutron scattering study

Background

Dynamics in haemoglobin from platypus (Ornithorhynchus anatinus), chicken (Gallus gallus domesticus) and saltwater crocodile (Crocodylus porosus) were measured to investigate response of conformational motions on the picosecond time scale to naturally occurring variations in the amino acid sequence of structurally identical proteins.

Methods

Protein dynamics was measured using incoherent quasielastic neutron scattering. The quasielastic broadening was interpreted first with a simple single Lorentzian approach and then by using the Kneller–Volino Brownian dynamics model.

Results

Mean square displacements of conformational motions, diffusion coefficients of internal dynamics and residence times for jump-diffusion between sites and corresponding effective force constants (resilience) and activation energies were determined from the data.

Conclusions

Modifications of the physicochemical properties caused by mutations of the amino acids were found to have a significant impact on protein dynamics. Activation energies of local side chain dynamics were found to be similar between the different proteins being close to the energy, which is required for the rupture of single hydrogen bond in a protein.

General significance

The measured dynamic quantities showed significant and systematic variations between the investigated species, suggesting that they are the signature of an evolutionary adaptation process stimulated by the different physiological environments of the respective protein.     

Neutron spin echo studies on ferritin: free-particle diffusion and interacting solutions

The dynamics of proteins are often studied by means of quasielastic neutron scattering (QENS), for example by time-of-flight methods. The spatial dimensions (10-20 nm) present in protein solutions are accessible by neutron scattering. In this article, a systematic study of diffusive dynamics of ferritin and apoferritin (=ferritin without iron core) is presented. Apoferritin consists of a spherical shell built of 24 protein units and carries net negative charge at pH 5. We have studied diffusive dynamics of ferritin solutions by neutron spin echo (NSE). We pay attention to an important feature of this technique compared to other QENS methods, which being the usage of a broad wavelength band. Using a more sophisticated fit function than usually used in NSE, we find as expected in low concentrated systems that the diffusion coefficient approaches the free-particle value of apoferritin and coincides with the diameter of the apoferritin shell (12.2 nm). In interacting solutions, the NSE results reveal that the dynamic picture of this complex liquid is dominated by slowing down of the dynamics. In low-salt solutions, a structure factor peak appears due to ordering of the ferritin molecules on the length scale of several intermolecular distances. We discuss the usage of different NSE fit functions for interacting solutions near the structure factor peak. Comparison of the dependence of elastic and dynamic data on the scattering vector value shows the influence of indirect interactions on the dynamic picture, irrespective of the way of data analysis, which being necessary due to the broad wavelength spectrum.     

Incoherent neutron scattering of copper azurin: a comparison with molecular dynamics simulation results

The low-frequency dynamics of copper azurin has been studied at different temperatures for a dry and deuterium hydrated sample by incoherent neutron scattering and the experimental results have been compared with molecular dynamics (MD) simulations carried out in the same temperature range. Experimental Debye-Waller factors are consistent with a dynamical transition at approximately 200 K which appears partially suppressed in the dry sample. Inelastic and quasielastic scattering indicate that hydration water modulates both vibrational and diffusive motions. The low-temperature experimental dynamical structure factor of the hydrated protein shows an excess of inelastic scattering peaking at about 3 meV and whose position is slightly shifted downwards in the dry sample. Such an excess is reminiscent of the “boson peak” observed in glass-like materials. This vibrational peak is quite well reproduced by MD simulations, although at a lower energy. The experimental quasielastic scattering of the two samples at 300 K shows a two-step relaxation behaviour with similar characteristic times, while the corresponding intensities differ only by a scale factor. Also, MD simulations confirm the two-step diffusive trend, but the slow process seems to be characterized by a decay faster than the experimental one. Comparison with incoherent neutron scattering studies carried out on proteins having different structure indicates that globular proteins display common elastic, quasielastic and inelastic features, with an almost similar hydration dependence, irrespective of their secondary and tertiary structure. Received: 12 October 1998 / Revised version: 19 February 1999 / Accepted: 1 March 1999     

The protein-solvent glass transition

The protein dynamical transition and its connection with the liquid-glass transition (GT) of hydration water and aqueous solvents are reviewed. The protein solvation shell exhibits a regular glass transition, characterized by steps in the specific heat and the thermal expansion coefficient at the calorimetric glass temperature T_G ≈ 170 K. It implies that the time scale of the structural α-relaxation has reached the experimental time window of 1–100 s. The protein dynamical transition, identified from elastic neutron scattering experiments by enhanced amplitudes of molecular motions exceeding the vibrational level [1], probes the α-process on a shorter time scale. The corresponding liquid-glass transition occurs at higher temperatures, typically 240 K. The GT is generally associated with diverging viscosities, the freezing of long-range translational diffusion in the supercooled liquid. Due to mutual hydrogen bonding, both, protein- and solvent relaxational degrees of freedom slow down in paralled near the GT. However, the freezing of protein motions, where surface-coupled rotational and librational degrees of freedom are arrested, is better characterized as a rubber-glass transition. In contrast, internal protein modes such as the rotation of side chains are not affected. Moreover, ligand binding experiments with myoglobin in various glass-forming solvents show, that only ligand entry and exit rates depend on the local viscosity near the protein surface, but protein-internal ligand migration is not coupled to the solvent. The GT leads to structural arrest on a macroscopic scale due to the microscopic cage effect on the scale of the intermolecular distance. Mode coupling theory provides a theoretical framework to understand the microcopic nature of the GT even in complex systems. The role of the α- and β-process in the dynamics of protein hydration water is evaluated. The protein-solvent GT is triggered by hydrogen bond fluctuations, which give rise to fast β-processes. High-frequency neutron scattering spectra indicate increasing hydrogen bond braking above T_G.     

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 H(2)O and D(2)O, 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 H(2)O 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 D(2)O (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.     

Biophysical aspects of neutron scattering from vibrational modes of proteins

This review describes a major portion of the published work on neutron scattering experiments aimed at measuring large scale motions in proteins. The importance of these motions for enzyme function and oxygen transport is indicated. The theory applicable to each type of neutron scattering measurement is given and results are discussed with a view to biological relevance. New experiments are suggested and a comparison of neutron scattering data is made with results from other techniques such as raman scattering, infrared absorption, photolysis and molecular dynamics simulations.     

Temperature dependence of protein dynamics: computer simulation analysis of neutron scattering properties

The temperature dependence of the internal dynamics of an isolated protein, bovine pancreatic trypsin inhibitor, is examined using normal mode analysis and molecular dynamics (MD) simulation. It is found that the protein exhibits marked anharmonic dynamics at temperatures of approximately 100-120 K, as evidenced by departure of the MD-derived average mean square displacement from that of the harmonic model. This activation of anharmonic dynamics is at lower temperatures than previously detected in proteins and is found in the absence of solvent molecules. The simulation data are also used to investigate neutron scattering properties. The effects are determined of instrumental energy resolution and of approximations commonly used to extract mean square displacement data from elastic scattering experiments. Both the presence of a distribution of mean square displacements in the protein and the use of the Gaussian approximation to the dynamic structure factor lead to quantified underestimation of the mean square displacement obtained.     

Picosecond dynamics of T and R forms of aspartate transcarbamylase: a neutron scattering study

E. coli aspartate transcarbamylase (ATCase) is a 310 kDa allosteric enzyme which catalyses the first committed step in pyrimidine biosynthesis. The binding of its substrates, carbamylphosphate and aspartate, induces significant conformational changes. This enzyme shows homotropic cooperative interactions between the catalytic sites for the binding of aspartate. This property is explained by a quaternary structure transition from T state (aspartate low affinity) to R state (aspartate high affinity) accompanied by a 5% increase of radius of gyration of ATCase. The same quaternary structure change is observed upon binding of the bisubstrate analogue PALA (N-(phosphonacetyl)-L-aspartate. Owing to the large incoherent neutron scattering cross-section of the hydrogen atom and the abundance of this element in proteins, inelastic neutron scattering gives a global view of protein dynamics as sensed via the individual motions of its hydrogen atoms. We present neutron scattering results of the local dynamics (few angstroms), at short time (few tens of picoseconds), of ATCase in T and R forms. Compared to the T form, we observe an increased mobility of the protein in the R form that we associate to an increase of accessible surface area to the solvent. Beyond this specific result, this highlights the key role of the accessible surface area (ASA) in dynamic contribution to inelastic neutron data in the picosecond time scale. In particular, we want to stress out (i) that a difference at the picosecond time scale does not allow to conclude to a difference in the dynamics at a longer time scale and to address whether the T state is looser than the R state (ii) how challenging is, any comparison in terms of general dynamics (tense or relaxed) between dynamic values deduced from experimental neutron data on proteins with different sequences and therefore ASA. This caveat holds particularly when comparing dynamics of a mesophile with the corresponding extremophile.     

Dynamics of apomyoglobin in the α-to-β transition and of partially unfolded aggregated protein

Changes of molecular dynamics in the α-to-β transition associated with amyloid fibril formation were explored on apomyoglobin (ApoMb) as a model system. Circular dichroism, neutron and X-ray scattering experiments were performed as a function of temperature on the protein, at different solvent conditions. A significant change in molecular dynamics was observed at the α-to-β transition at about 55°C, indicating a more resilient high temperature β structure phase. A similar effect at approximately the same temperature was observed in holo-myoglobin, associated with partial unfolding and protein aggregation. A study in a wide temperature range between 20 and 360 K revealed that a dynamical transition at about 200 K for motions in the 50 ps time scale exists also for a hydrated powder of heat-denatured aggregated ApoMb.     

The relation of the X-ray B-factor to protein dynamics: insights from recent dynamic solid-state NMR data

In addressing the potential use of B-factors derived from X-ray scattering data of proteins for the understanding the (functional) dynamics of proteins, we present a comparison of B-factors of five different proteins (SH3 domain, Crh, GB1, ubiquitin and thioredoxin) with data from recent solid-state nuclear magnetic resonance experiments reflecting true (rotational) dynamics on well-defined timescales. Apart from trivial correlations involving mobile loop regions and chain termini, we find no significant correlation of B-factors with the dynamic data on any of the investigated timescales, concluding that there is no unique and general correlation of B-factors with the internal reorientational dynamics of proteins.     

Protein interactions and dynamics probed by quantum chemistry, computer simulations and neutron experiments

We review some recent experiments and calculations on aspects of the structure and dynamics of proteins and related systems. The use of quantum chemical techniques to determine geometries and energies of supramolecular complexes of biological interest is illustrated, and the concomitant development of empirical energy functions for use in protein simulations outlined. We describe how simulations of crystalline peptides and amino-acids using an empirical force field can be combined with appropriate coherent and incoherent inelastic neutron scattering experiments to elucidate the characteristics of lattice vibrations and diffusive atomic motions in the crystals. The application of molecular dynamics simulations to the interpretation of incoherent neutron scattering experiments on proteins is examined and the resulting ideas on the general characteristics of protein motion discussed in terms of their functional implications.     

Internal motions of actin characterized by quasielastic neutron scattering

Quasielastic neutron scattering (QENS) experiments were carried out on powders of F-actin and G-actin hydrated with D₂O to characterize the internal dynamics on the picosecond time scale and the Ångstrom length scale. To investigate the effects of hydration, the measurements were done on samples at hydration ratio (h) of 0.4 (mg D₂O/mg protein), containing only the first layer of hydration water, and at h = 1.0, containing more layers of water. The QENS spectra, obtained from the measurements at two energy resolutions of 110 and 15 μeV, indicated that the internal motions of both F-actin and G-actin have distributions of motions with distinct correlation times and amplitudes. Increasing hydration changes relative populations of these distinct motions. The effects of hydration were shown to be different between F-actin and G-actin. Elastic incoherent neutron scattering measurements provided the concerted results. The observed effects were interpreted in terms of the dynamical heterogeneity of the actin molecule: in G-actin, more surface loops become flexible and undergo diffusive motions of large amplitudes, whereas in F-actin the molecular interactions that keep the polymerized state suppress the large motions of the surface loops involved with polymerization so that the population of atoms undergoing large motions can increase only to a lesser degree.     

Coupling of laser excitation and inelastic neutron scattering: attempt to probe the dynamics of light-induced C-phycocyanin dynamics

Excitation energy transfer (EET) in light-harvesting antennae is a highly efficient key event in photosynthesis, where light-induced dynamics of the antenna pigment-protein complexes may play a functional role. So far, however, the relationship between EET and protein dynamics remains unknown. C-phycocyanin (C-PC) is the main pigment/protein complex present in the cyanobacterial antenna, called "phycobilisome". The aim of the present study was to investigate light-induced C-PC internal thermal motions (ps timescale) measured by inelastic neutron scattering. To synchronize the beginning of the laser flash (6 ns duration) with that of the neutron test pulse ( approximately 87 mus duration), we developed a novel type of "time-resolved" experimental setup on MIBEMOL time-of-flight neutron spectrometer (LLB, France). Data acquisition has been modified to get quasi-simultaneously "light" and "dark" measurements (with and without laser, respectively) and eliminate many spurious effects that could occur on the sample during the experiment. The study was carried out on concentrated C-PC ( approximately 135 g/L protein in D(2)O phosphate buffer), contained in an aluminium/sapphire sample holder (almost "transparent" for neutrons) and homogeneously illuminated inside an "integrating sphere". We observed very similar incoherent dynamical structure factors of C-PC with or without light. The vibrational density of states showed two very slightly increased vibrational modes with light, at approximately 30 and approximately 50 meV ( approximately 240 and approximately 400 cm(-1), respectively). These effects have to be verified by further experiments before probing any temporal evolution, by introducing a time delay between the laser flash and the neutron test pulse.     

Nanoscale protein domain motion and long-range allostery in signaling proteins—a view from neutron spin echo spectroscopy

Many cellular proteins are multi-domain proteins. Coupled domain–domain interactions in these multidomain proteins are important for the allosteric relay of signals in the cellular signaling networks. We have initiated the application of neutron spin echo spectroscopy to the study of nanoscale protein domain motions on submicrosecond time scales and on nanometer length scale. Our NSE experiments reveal the activation of protein domain motions over a long distance of over more than 100 Å in a multidomain scaffolding protein NHERF1 upon binding to another protein, Ezrin. Such activation of nanoscale protein domain motions is correlated with the allosteric assembly of multi-protein complexes by NHERF1 and Ezrin. Here, we summarize the theoretical framework that we have developed, which uses simple concepts from nonequilibrium statistical mechanics to interpret the NSE data, and employs a mobility tensor to describe nanoscale protein domain motion. Extracting nanoscale protein domain motion from the NSE does not require elaborate molecular dynamics simulations, nor complex fits to rotational motion, nor elastic network models. The approach is thus more robust than multiparameter techniques that require untestable assumptions. We also demonstrate that an experimental scheme of selective deuteration of a protein subunit in a complex can highlight and amplify specific domain dynamics from the abundant global translational and rotational motions in a protein. We expect NSE to provide a unique tool to determine nanoscale protein dynamics for the understanding of protein functions, such as how signals are propagated in a protein over a long distance to a distal domain.     

Crowding in extremophiles: linkage between solvation and weak protein-protein interactions, stability and dynamics, provides insight into molecular adaptation

The study of the molecular adaptation of microorganisms to extreme environments (solvent, temperature, etc.) has provided tools to investigate the complex relationships between protein-solvent and protein-protein interactions, protein stability and protein dynamics, and how they are modulated by the crowded environment of the cell. We have evaluated protein-solvent and protein-protein interactions by solution experiments (analytical ultracentrifugation, small angle neutron and X-ray scattering, density) and crystallography, and protein dynamics by energy resolved neutron scattering. This review concerns work from our laboratory on (i) proteins from extreme halophilic Archaea, and (ii) psychrophile, mesophile, thermophile and hyperthermophile bacterial cells.     

Neutron Spin-Echo Studies of Hemoglobin and Myoglobin: Multiscale Internal Dynamics

Neutron spin-echo spectroscopy was used to study structural fluctuations that occur in hemoglobin (Hb) and myoglobin (Mb) in solution. Using neutron spin-echo data up to a very high momentum transfer q (∼ 0.62 Å^− 1), we characterized the internal dynamics of these proteins at the levels of dynamic pair correlation function and self-correlation function in the time range of several picoseconds to a few nanoseconds. In the same protein solution, data transition from pair correlation motion to self-correlation motion as the momentum transfer q increases. At low q, coherent scattering dominates; at high q, observations are largely due to incoherent scattering. The low q data were interpreted in terms of an effective diffusion coefficient; on the other hand, the high q data were interpreted in terms of mean square displacements. Comparison of data from the two homologous proteins collected at different temperatures and protein concentrations was used to assess the contributions made by translational and rotational diffusion and internal modes of motion to the data. The temperature dependence of decay times can be attributed to changes in the viscosity and temperature of the solvent, as predicted by the Stokes-Einstein relationship. This is true for contributions from both diffusive and internal modes of motion, indicating an intimate relationship between the internal dynamics of the proteins and the viscosity of the solvent. Viscosity change associated with protein concentration can account for changes in diffusion observed at different concentrations, but is apparently not the only factor involved in the changes in internal dynamics observed with change in protein concentration. Data collected at high q indicate that internal modes in Mb are generally faster than those in Hb, perhaps due to the greater surface-to-volume ratio of Mb and the fact that surface groups tend to exhibit faster motion than buried groups. Comparison of data from Hb and data from Mb at low q indicates an unexpectedly rapid motion of Hb αβ dimers relative to one another. Dynamic motion of subunits is increasingly perceived as important to the allosteric behavior of Hb. Our data demonstrate that this motion is highly sensitive to protein concentration, temperature, and solvent viscosity, indicating that great care needs to be exercised in interpreting its effect on protein function.     

Combination of Neutron Scattering and Molecular Dynamics to Determine Internal Motions in Biomolecules

Abstract

The forms and frequencies of atomic dynamics on the pico- and nanosecond timescales are accessible experimentally using incoherent neutron scattering. Molecular dynamics simulations cover the same space and time domains and neutron scattering intensities can be calculated from the simulations for direct comparison with experiment. To illustrate the complementarity of neutron scattering and molecular dynamics we examine measured and simulation-derived elastic incoherent scattering profiles from myoglobin and from the crystalline alanine dipeptide. Elastic incoherent scattering gives information on the geometry of the volume accessible to the atoms in the samples. The simulation-derived dipeptide elastic scattering profiles are in reasonable accord with experiment, deviations being due to the sampling limitations in the simulations and experimental detector normalisation procedures. The simulated dynamics is decomposed, revealing characteristic profiles due to rotational diffusional and translational vibrational motions of the methyl groups. In myoglobin, for which the timescale of the simulation matches more closely that accessible to the experiment, good agreement is seen for the elastic incoherent structure factor. This indicates that the space sampled by the hydrogen atoms in the protein on the timescale <100 ps is well represented by the simulation. Part of the helix atom fluctuations can be described in terms of rigid helix motions.     

Probing the flexibility of the bacterial reaction center: the wild-type protein is more rigid than two site-specific mutants

Experimental and theoretical studies have stressed the importance of flexibility for protein function. However, more local studies of protein dynamics, using temperature factors from crystallographic data or elastic models of protein mechanics, suggest that active sites are among the most rigid parts of proteins. We have used quasielastic neutron scattering to study the native reaction center protein from the purple bacterium Rhodobacter sphaeroides, over a temperature range of 4-260 K, in parallel with two nonfunctional mutants both carrying the mutations L212Glu/L213Asp --> Ala/Ala (one mutant carrying, in addition, the M249Ala --> Tyr mutation). The so-called dynamical transition temperature, Td, remains the same for the three proteins around 230 K. Below Td the mean square displacement, u2, and the dynamical structure factor, S(Q,omega), as measured respectively by backscattering and time-of-flight techniques are identical. However, we report that above Td, where anharmonicity and diffusive motions take place, the native protein is more rigid than the two nonfunctional mutants. The higher flexibility of both mutant proteins is demonstrated by either their higher u2 values or the notable quasielastic broadening of S(Q,omega) that reveals the diffusive nature of the motions involved. Remarkably, we demonstrate here that in proteins, point genetic mutations may notably affect the overall protein dynamics, and this effect can be quantified by neutron scattering. Our results suggest a new direction of investigation for further understanding of the relationship between fast dynamics and activity in proteins. Brownian dynamics simulations we have carried out are consistent with the neutron experiments, suggesting that a rigid core within the native protein is specifically softened by distant point mutations. L212Glu, which is systematically conserved in all photosynthetic bacteria, seems to be one of the key residues that exerts a distant control over the rigidity of the core of the protein.