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.     

Differences in internal dynamics of actin under different structural states detected by neutron scattering

F-actin, a helical polymer formed by polymerization of the monomers (G-actin), plays crucial roles in various aspects of cell motility. Flexibility of F-actin has been suggested to be important for such a variety of functions. Understanding the flexibility of F-actin requires characterization of a hierarchy of dynamical properties, from internal dynamics of the actin monomers through domain motions within the monomers and relative motions between the monomers within F-actin to large-scale motions of F-actin as a whole. As a first step toward this ultimate purpose, we carried out elastic incoherent neutron scattering experiments on powders of F-actin and G-actin hydrated with D₂O and characterized the internal dynamics of F-actin and G-actin. Well established techniques and analysis enabled the extraction of mean-square displacements and their temperature dependence in F-actin and in G-actin. An effective force constant analysis with a model consisting of three energy states showed that two dynamical transitions occur at ∼150 K and ∼245 K, the former of which corresponds to the onset of anharmonic motions and the latter of which couples with the transition of hydration water. It is shown that behavior of the mean-square displacements is different between G-actin and F-actin, such that G-actin is “softer” than F-actin. The differences in the internal dynamics are detected for the first time between the different structural states (the monomeric state and the polymerized state). The different behavior observed is ascribed to the differences in dynamical heterogeneity between F-actin and G-actin. Based on structural data, the assignment of the differences observed in the two samples to dynamics of specific loop regions involved in the polymerization of G-actin into F-actin is proposed.     

Equilibrium fluctuations in myoglobin and lysozyme

The angular dependencies of inelastic intensities of Rayleigh scattering of Moessbauer radiation were measured for myoglobin and lysozyme (in the hydration range h = 0.05-0.7). The data were fitted within the framework of model, when two types of intraglobular motions were taken into account: individual motions of small side-chain groups and cooperative motions of segments. The best agreement with the experiment at h > 0.05 was obtained when individual motions of small groups together with the cooperative motions of alpha-helices and beta-sheets for lysozyme, and alpha-helices for myoglobin were considered. At further hydration (h = 0.45), mean-square displacements (x2) of both types of motions strongly increase with the increase in hydration degree, while the motions with a large correlation radius (not less than macromolecule radius) remain nearly the same as for h = 0.05. The results of the study of the radial distribution function deduced by Fourier-transform from the diffuse x-ray measurements together with RSMR data allow one to conclude that the water during protein hydration competes with the intramolecular hydrogen bonds, loosens the protein and increases the internal dynamics. Concurrently, water arranges the ordering of macromolecule, which takes the native structure at h = 0.4-0.7. The analysis of auto and cross-correlation functions of bending fluctuations of alpha-helices in the large domain of lysozyme performed by molecular dynamics allows one to come to the final conclusion that it is the difference in the structural organization of myoglobin and lysozyme and not the presence of SS-bonds in lysozyme macromolecule that is responsible for different structural fluctuations in these proteins.     

The dynamical transition of proteins,concepts and misconceptions

The dynamics of hydrated proteins and of protein crystals can be studied within a wide temperature range, since the water of hydration does not crystallize at low temperature. Instead it turns into an amorphous glassy state below 200 K. Extending the temperature range facilitates the spectral separation of different molecular processes. The conformational motions of proteins show an abrupt enhancement near 180 K, which has been called a "dynamical transition". In this contribution various aspects of the transition are critically reviewed: the role of the instrumental resolution function in extracting displacements from neutron elastic scattering data and the question of the appropriate dynamic model, discrete transitions between states of different energy versus continuous diffusion inside a harmonic well, are discussed. A decomposition of the transition involving two motional components is performed: rotational transitions of methyl groups and small scale librations of side-chains, induced by water at the protein surface. Both processes create an enhancement of the observed amplitude. The onset occurs, when their time scale becomes compatible with the resolution of the spectrometer. The reorientational rate of hydration water follows a super-Arrhenius temperature dependence, a characteristic feature of a dynamical transition. It occurs only with hydrated proteins, while the torsional motion of methyl groups takes place also in the dehydrated or solvent-vitrified system. Finally, the role of fast hydrogen bond fluctuations contributing to the amplitude enhancement is discussed.     

The essential dynamics of Cu, Zn superoxide dismutase: suggestion of intersubunit communication.

A 300-ps molecular dynamics simulation of the whole Cu, Zn superoxide dismutase dimer has been carried out in water, and the trajectory has been analyzed by the essential dynamics method. The results indicate that the motion is defined by few preferred directions identified by the first four to six eigenvectors and that the motion of the two monomers at each instant is not symmetrical. The vectors symmetrical to the eigenvectors are significantly sampled, suggesting that, on average, the motions of the two subunits will exchange. Large intra- and intersubunit motions involving different subdomains of the protein are observed. A mechanical coupling between the two subunits is also suggested, because displacements of the loops surrounding the active site in one monomer are correlated with the motion of parts of the second toward the intersubunit interface.     

Dynamics of C-phycocyanin in various deuterated trehalose/water environments measured by quasielastic and elastic neutron scattering

The molecular understanding of protein stabilization by the disaccharide trehalose in extreme temperature or hydration conditions is still debated. In the present study, we investigated the role of trehalose on the dynamics of the protein C-phycocyanin (C-PC) by neutron scattering. To single out the motions of C-PC hydrogen (H) atoms in various trehalose/water environments, measurements were performed in deuterated trehalose and heavy water (D(2)O). We report that trehalose decreases the internal C-PC dynamics, as shown by a reduced diffusion coefficient of protein H atoms. By fitting the Elastic Incoherent Structure Factor-which gives access to the "geometry" of the internal proton motions-with the model of diffusion inside a sphere, we found that the presence of trehalose induces a significantly higher proportion of immobile C-PC hydrogens. We investigated, by elastic neutron scattering, the mean square displacements (MSDs) of deuterated trehalose/D(2)O-embedded C-PC as a function of temperature in the range of 40-318 K. Between 40 and approximately 225 K, harmonic MSDs of C-PC are slightly smaller in samples containing trehalose. Above a transition temperature of approximately 225 K, we observed anharmonic motions in all trehalose/water-coated C-PC samples. In the hydrated samples, MSDs are not significantly changed by addition of 15% trehalose but are slightly reduced by 30% trehalose. In opposition, no dynamical transition was detected in dry trehalose-embedded C-PC, whose hydrogen motions remain harmonic up to 318 K. These results suggest that a role of trehalose would be to stabilize proteins by inhibiting some fluctuations at the origin of protein unfolding and denaturation.     

Evolution of the internal dynamics of two globular proteins from dry powder to solution.

Myoglobin and lysozyme picosecond internal dynamics in solution is compared to that in hydrated powders by quasielastic incoherent neutron scattering. This technique is sensitive to the motions of the nonexchangeable hydrogen atoms in a sample. Because these are homogeneously distributed throughout the protein structure, the average dynamics of the protein is described. We first propose an original data treatment to deal with the protein global motions in the case of solution samples. The validity of this treatment is checked by comparison with classical measurements of the diffusion constants. The evolution with the scattering vector of the width and relative contribution of the quasielastic component was then used to derive information on the amount of local diffusive motions and their characteristic average relaxation time. From dry powder to coverage by one water layer, the surface side chains progressively acquire the possibility to diffuse locally. On subsequent hydration, the main effect of water is to improve the rate of these diffusive motions. Motions with higher average amplitude occur in solution, about three times more than for a hydrated powder at complete coverage, with a shorter average relaxation time, approximately 4.5 ps compared to 9.4 ps for one water monolayer.     

Correlations of atomic movements in lysozyme crystals.

Diffuse scattering data have been collected on two crystal forms of lysozyme, tetragonal and triclinic, using synchrotron radiation. The observed diffraction patterns were simulated using an exact theory for simple model crystals which relates the diffuse scattering intensity distribution to the amplitudes and correlations of atomic movements. Although the mean square displacements in the tetragonal form are twice that in the triclinic crystal, the predominant component of atomic movement in both crystals is accounted for by short-range coupled motions where displacement correlations decay exponentially as a function of atomic separation, with a relaxation distance of approximately 6 A. Lattice coupled movements with a correlation distance approximately 50 A account for only about 5-10% of the total atomic mean square displacements in the protein crystals. The results contradict various presumptions that the disorder in protein crystals can be modeled predominantly by elastic vibrations or rigid body movements.     

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.     

Enzyme dynamics and activity: time-scale dependence of dynamical transitions in glutamate dehydrogenase solution.

We have examined the temperature dependence of motions in a cryosolution of the enzyme glutamate dehydrogenase (GDH) and compared these with activity. Dynamic neutron scattering was performed with two instruments of different energy resolution, permitting the separate determination of the average dynamical mean square displacements on the sub-approximately 100 ps and sub-approximately 5 ns time scales. The results demonstrate a marked dependence on the time scale of the temperature profile of the mean square displacement. The lowest temperature at which anharmonic motion is observed is heavily dependent on the time window of the instrument used to observe the dynamics. Several dynamical transitions (inflexions of the mean squared displacement) are observed in the slower dynamics. Comparison with the temperature profile of the activity of the enzyme in the same solvent reveals dynamical transitions that have no effect on GDH function.     

Temperature dependence of dynamics of hydrated myoglobin. Comparison of force field calculations with neutron scattering data

Molecular dynamics is used to probe the atomic motions of the carboxy-myoglobin protein as a function of temperature. Simulations of 150 picoseconds in length are carried out on the protein at 20, 60, 100, 180, 220, 240, 260, 280, 300, 320 and 340 K. The simulations attempt to mimic neutron scattering experiments very closely by including a partial hydration shell around the protein. Theoretical elastic, quasielastic and inelastic neutron scattering data are derived from the trajectories and directly compared with experiment. Compared to experiment, the simulation-derived elastic scattering curves show a decrease in intensity as a function of the scattering wavevector, q2. The inelastic and quasielastic spectra show that the inelastic peak is shifted to lower frequency than the experimental value, while quasielastic behavior is in good agreement with experiment. This suggests that the theoretical model is too flexible in the harmonic limit (low temperature), but accurately reproduces high-temperature behavior. Time correlation functions of the intermediate scattering function are determined. At low temperature there is one fast decay process, and at high temperatures there is an additional slow relaxation process that is due to quasielastic scattering. The average atomic fluctuations show that the protein behaves harmonically at low temperatures. At approximately 210 K, a glass-like transition in atomic fluctuations is seen. Above the transition temperature, the atomic fluctuations exhibit both harmonic and anharmonic behavior. Comparison of protein mobility behavior with experiment indicate the fluctuations derived from simulations are larger in the harmonic region. However, the anharmonic region agrees very well with experiment. The anharmonicity is large at all temperatures, with a gradual monotonic increase from 0.5 at 20 K to greater than 0.7 at 340 K without a noticeable change at the glass transition temperature. Heavy-atom dihedral transitions are monitored as a function of temperature. Trends in the type of dihedral transitions that occur with temperature are clearly visible. Dihedral transitions involving backbone atoms occur only above the glass transition temperature. The overall protein behavior results suggest that at low temperatures there is purely vibrational motion with one fast decay process, and above the glass transition temperature there is more anharmonic motion with a fast and a slower relaxation process occurring simultaneously.(ABSTRACT TRUNCATED AT 400 WORDS)     

Deformation-enhanced fluctuations in the red cell skeleton with theoretical relations to elasticity, connectivity, and spectrin unfolding.

To assess local elasticity in the red cell's spectrin-actin network, nano-particles were tethered to actin nodes and their constrained thermal motions were tracked. Cells were both immobilized and controllably deformed by aspiration into a micropipette. Since the network is well-appreciated as soft, thermal fluctuations even in an unstressed portion of network were expected to be many tens of nanometers based on simple equipartition ideas. Real-time particle tracking indeed reveals such root-mean-squared motions for 40-nm fluorescent beads either tethered to actin directly within a cell ghost or connected to actin from outside a cell via glycophorin. Moreover, the elastically constrained displacements are significant on the scale of the network's internodal distance of approximately 60-80 nm. Surprisingly, along the aspirated projection-where the network is axially extended by as much as twofold or more-fluctuations in the axial direction are increased by almost twofold relative to motions in the unstressed network. The molecular basis for such strain softening is discussed broadly in terms of force-driven transitions. Specific considerations are given to 1) protein dissociations that reduce network connectivity, and 2) unfolding kinetics of a localized few of the red cell's approximately 10(7) spectrin repeats.     

Direct visualization of dynamics and co-operative conformational changes within RecA filaments that appear to be associated with the hydrolysis of adenosine 5'-O-(3-thiotriphosphate).

Highly co-operative structural transitions and conformational changes can be directly observed in bundles of filaments formed by the RecA protein of Escherichia coli. These filaments have been formed with RecA protein, DNA and the ATP analog adenosine 5'-O-(3-thiotriphosphate) (ATP-gamma-S). We show that while ATP-gamma-S has frequently been called non-hydrolyzable in the RecA literature, it is hydrolyzed by RecA with a kcat of about 0.01 to 0.005 min-1. This rate of ATP-gamma-S hydrolysis is significant to structural studies conducted on a time scale of hours. It has been shown that RecA subunits may be seen in different conformations within one particular form of RecA bundle. We now show that additional structural transitions take place within these bundles when they are allowed to incubate at 37 degrees C for several hours. This is the same time scale on which ATP-gamma-S is being hydrolyzed, and the suggestion that the observed structural transitions arise from the hydrolysis of ATP-gamma-S is supported by the fact that when the hydrolysis of ATP-gamma-S is inhibited (at 4 degrees C), the transitions are not observed. The transitions that occur are highly co-operative, with filaments as a whole changing their state over lengths of several thousand Angstroms. This shows that RecA filaments have an internal co-operativity, and we suggest that this is important to their function in vivo. The motions of subunits that we visualize appear to be mainly rotational, and this can be used to infer information about the motions of RecA subunits associated with the RecA ATPase that occurs during the DNA strand exchange reaction.     

Networks for the allosteric control of protein kinases

The allosteric regulation of protein kinases serves as an efficient strategy for molecular communication, event coupling and interconversion between catalytic states. Recent co-crystal structures have revealed novel ways in which kinases control activity and substrate specificity following phosphorylation, dimerization, or binding to regulatory proteins, substrates and scaffolds. In addition, hydrogen exchange coupled with mass spectrometry is emerging as a complementary strategy to probe the solution behavior of kinases; recent results have shown that allosteric regulation may involve transitions in protein motions as well as structural rearrangements.     

Molecular dynamics study of the KcsA channel at 2.0-A resolution: stability and concerted motions within the pore

The stability of the KcsA channel accommodating more than one ion in the pore has been studied with molecular dynamics. We have used the very last X-ray structure of the KcsA channel at 2.0-A resolution determined by Zhou et al. [Nature 414 (2001) 43]. In this channel, six of the seven experimentally evidenced sites have been considered. We show that the protein remains very stable in the presence of four K+ ions (three in the selectivity filter and one in the cavity). The locations and the respective distances of the different K+ ions and water molecules (W), calculated within our KWKWKK sequence, also fits well with the experimental observations. The analysis of the K+ ions and water molecules displacements shows concerted file motions on the simulated time scale (approximately 1 ns), which could act as precursor to the diffusion of K+ ions inside the channel. A simple one-dimensional dynamical model is used to interpret the concerted motions of the ions and water molecules in the pore leading ultimately to ion transfer.     

Classification and annotation of the relationship between protein structural change and ligand binding

The causal relationship between protein structural change and ligand binding was classified and annotated for 839 nonredundant pairs of crystal structures in the Protein Data Bank—one with and the other without a bound low-molecular-weight ligand molecule. Protein structural changes were first classified into either domain or local motions depending on the size of the moving protein segments. Whether the protein motion was coupled with ligand binding was then evaluated based on the location of the ligand binding site and by application of the linear response theory of protein structural change. Protein motions coupled with ligand binding were further classified into either closure or opening motions. This classification revealed the following: (i) domain motions coupled with ligand binding are dominated by closure motions, which can be described by the linear response theory; (ii) local motions frequently accompany order-disorder or α-helix-coil conformational transitions; and (iii) transferase activity (Enzyme Commission   number 2) is the predominant function among coupled domain closure motions. This could be explained by the closure motion acting to insulate the reaction site of these enzymes from environmental water.     

A physical picture of atomic motions within the Dickerson DNA dodecamer in solution derived from joint ensemble refinement against NMR and large-angle X-ray scattering data

The structure and dynamics of the Dickerson DNA dodecamer [5'd(CGCGAATTCGCG)2] in solution have been investigated by joint simulated annealing refinement against NMR and large-angle X-ray scattering data (extending from 0.25 to 3 A-1). The NMR data comprise an extensive set of hetero- and homonuclear residual dipolar coupling and 31P chemical shift anisotropy restraints in two alignment media, supplemented by NOE and 3J coupling data. The NMR and X-ray scattering data cannot be fully ascribed to a single structure representation, indicating the presence of anisotropic motions that impact the experimental observables in different ways. Refinement with ensemble sizes (Ne) of >or=2 to represent the atomic motions reconciles all the experimental data within measurement error. Cross validation against both the dipolar coupling and X-ray scattering data suggests that the optimal ensemble size required to account for the current data is 4. The resulting ensembles permit one to obtain a detailed view of the conformational space sampled by the dodecamer in solution and permit one to analyze fluctuations in helicoidal parameters, sugar puckers, and BI-BII backbone transitions and to obtain quantitative metrics of atomic motion such as generalized order parameters and thermal B factors. The calculated order parameters are in good agreement with experimental order parameters obtained from 13C relaxation measurements. Although DNA behaves as a relatively rigid rod with a persistence length of approximately 150 bp, dynamic conformational heterogeneity at the base pair level is functionally important since it readily permits optimization of intermolecular protein-DNA interactions.     

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.