Slow motions in chicken villin headpiece subdomain probed by cross-correlated NMR relaxation of amide NH bonds in successive residues

The villin headpiece subdomain (HP36) is a widely used system for protein-folding studies. Nuclear magnetic resonance cross-correlated relaxation rates arising from correlated fluctuations of two N-H^N dipole-dipole interactions involving successive residues were measured at two temperatures at which HP36 is at least 99% folded. The experiment revealed the presence of motions slower than overall tumbling of the molecule. Based on the theoretical analysis of the spectral densities we show that the structural and dynamic contributions to the experimental cross-correlated relaxation rate can be separated under certain conditions. As a result, dynamic cross-correlated order parameters describing slow microsecond-to-millisecond motions of N-H bonds in neighboring residues can be introduced for any extent of correlations in the fluctuations of the two bond vectors. These dynamic cross-correlated order parameters have been extracted for HP36. The comparison of their values at two different temperatures indicates that when the temperature is raised, slow motions increase in amplitude. The increased amplitude of these fluctuations may reflect the presence of processes directly preceding the unfolding of the protein.     

Investigation of the backbone dynamics of the IgG-binding domain of streptococcal protein G by heteronuclear two-dimensional 1H-15N nuclear magnetic resonance spectroscopy.

The backbone dynamics of the immunoglobulin-binding domain (B1) of streptococcal protein G, uniformly labeled with 15N, have been investigated by two-dimensional inverse detected heteronuclear 1H-15N NMR spectroscopy at 500 and 600 MHz. 15N T1, T2, and nuclear Overhauser enhancement data were obtained for all 55 backbone NH vectors of the B1 domain at both field strengths. The overall correlation time obtained from an analysis of the T1/T2 ratios was 3.3 ns at 26 degrees C. Overall, the B1 domain is a relatively rigid protein, consistent with the fact that over 95% of the residues participate in secondary structure, comprising a four-stranded sheet arranged in a -1, +3x, -1 topology, on top of which lies a single helix. Residues in the turns and loops connecting the elements of secondary structure tend to exhibit a higher degree of mobility on the picosecond time scale, as manifested by lower values of the overall order parameter. A number of residues at the ends of the secondary structure elements display two distinct internal motions that are faster than the overall rotational correlation time: one is fast (< 20 ps) and lies in the extreme narrowing limit, whereas the other is one to two orders of magnitude slower (1-3 ns) and lies outside the extreme narrowing limit. The slower motion can be explained by large-amplitude (20-40 degrees) jumps in the N-H vectors between states with well-defined orientations that are stabilized by hydrogen bonds.(ABSTRACT TRUNCATED AT 250 WORDS)     

Backbone motions in a crystalline protein from field-dependent 2H-NMR relaxation and line-shape analysis

We have used 2H-nmr to study backbone dynamics of the 2H-labeled, slowly exchanging amide sites of fully hydrated, crystalline hen egg white lysozyme. Order parameters are determined from the residual quadrupole coupling and values increase from S2 = 0.85 at 290 K to S2 = 0.94 at 200 K. Dynamical rates are determined from spin-lattice relaxation at three nmr frequencies (38.8, 61.5, and 76.7 MHz). The approach used here is thus distinct from solution nmr studies where dynamical amplitudes and rates are both determined from relaxation measurements. At temperatures below 250 K, relaxation is independent of the nmr frequency indicating that backbone motions are fast compared to the nmr frequencies. However, as the temperature is increased above 250 K, relaxation is significantly more efficient at the lowest frequency, which shows, in addition, the presence of motions that are slow compared to the nmr frequencies. Using the values of S2 determined from the residual quadrupole coupling and a model-free relaxation formalism that allows for fast and slow internal motions, we conclude that these slow motions have correlation times in the range of 0.1 to 1.0 microsecond and are effectively frozen out at 250 K where fast motions of the amide planes with approximately 15 ps effective correlation times and 9 degrees rms amplitudes dominate relaxation. The fast internal motions increase slightly in amplitude as the temperature rises toward 290 K, but the correlation time, as is also observed in solution nmr studies of RNase H, is approximately constant. These findings are consistent with hypotheses of dynamic glass transitions in hydrated proteins arising from temperature-dependent damping of harmonic modes of motion above the transition point.     

13C Nuclear magnetic relaxation study of segmental dynamics of hyaluronan in aqueous solutions

(13)C spin-lattice relaxation times (T(1)) and nuclear Overhauser enhancements (NOE) were measured as a function of temperature and magnetic field strength for the hetero-polysaccharide hyaluronan in water solutions. The relaxation data of the endocyclic ring carbons were successfully interpreted in terms of chain segmental motions by using the bimodal time-correlation function of Dejean de la Batie, Laupretre and Monnerie. On the basis of the calculated correlation times for segmental motion and amplitudes of librational motions of the C-H vectors at the various carbon sites of the HA repeating unit, we concluded that intramolecular hydrogen bonding of the secondary structure of HA plays a major role in the conformational flexibility of this carbohydrate molecule. The internal rotation of the free hydroxymethyl groups about the exocyclic C-5-C-6 bonds superimposed on segmental motion has been described as a diffusion process of restricted amplitude. The rate and amplitude of the internal rotation indicate that the hydroxymethyl groups are not involved in intramolecular hydrogen bonding. Finally, the motional parameters describing the local dynamics of the HA chain were correlated with the secondary structure of HA in aqueous solutions.     

Structure and dynamics of an amphiphilic peptide in a lipid bilayer: a molecular dynamics study

A molecular dynamics simulation of a simple model membrane system composed of a single amphiphilic helical peptide (ace-K2GL16K2A-amide) in a fully hydrated 1,2-dimyristoyl-sn-glycero-3-phosphocholine bilayer was performed for a total of 1060 ps. The secondary structure of the peptide and its stability were described in terms of average dihedral angles, phi and psi, and the C alpha torsion angles formed by backbone atoms; by the average translation per residue along the helix axis; and by the intramolecular peptide hydrogen bonds. The results indicated that residues 6 through 15 remain in a stable right-handed alpha-helical conformation, whereas both termini exhibit substantial fluctuations. A change in the backbone dihedral angles for residues 16 and 17 is accompanied by the loss of two intramolecular hydrogen bonds, leading to a local but long-lived disruption of the helix. The dynamics of the peptide was characterized in terms of local and global helix motions. The local motions of the N-H bond angles were described in terms of the autocorrelation functions of P2[cos thetaNH(t, t + tau)] and reflected the different degrees of local peptide order as well as a variation in time scale for local motions. The chi1 and chi2 dihedral angles of the leucine side chains underwent frequent transitions between potential minima. No connection between the side-chain positions and their mobility was observed, however. In contrast, the lysine side chains displayed little mobility during the simulation. The global peptide motions were characterized by the tilting and bending motions of the helix. Although the peptide was initially aligned parallel to the bilayer normal, during the simulation it was observed to tilt away from the normal, reaching an angle of approximately 25 degrees by the end of the simulation. In addition, a slight bend of the helix was detected. Finally, the solvation of the peptide backbone and side-chain atoms was also investigated.     

Backbone dynamics and thermodynamics of Borrelia outer surface protein A

Nuclear spin relaxation experiments performed at 298K, 308K and 318K are used to characterize the intramolecular dynamics and thermodynamics of outer surface protein A (OspA), a key protein in the life-cycle of Borrelia burgdorferi, the causative agent of Lyme disease. It has recently been demonstrated that OspA specifically binds to the gut of the intermediate tick host (Ixodes scapularis), and that this interaction is mediated, at least in part, by residues in the C-terminal domain of OspA that are largely inaccessible to solvent in all X-ray structures of this protein. Our analysis of 15N relaxation parameters in OspA shows that the putative-binding region contains and is surrounded by flexible residues, which could facilitate accessibility to solvent and ligands. In addition, residues with similar activation energies are clustered in a manner that suggests locally collective motions. We have used molecular modeling to show that these collective motions are consistent with a hinge-bending mechanism that exposes residues implicated in binding. Characteristic temperatures describing the energy landscape of the OspA backbone are derived from the temperature dependence of the N-H bond vector order parameters, and a comparison is made between the N and C-terminal globular domains and the unusual single-layer beta-sheet connecting them. The average characteristic temperatures in the three regions indicate that, with an increase in temperature, a larger increase in accessible conformational states occurs for N-H bond vectors in the single-layer central beta-sheet than for bond vectors in the globular N and C-terminal domains. These conformational states are accessible without disruption of hydrogen bonds, providing a conformational entropic gain, upon increase in temperature, without a significant enthalpic penalty. This increase in heat capacity may help to explain the unexpected thermal stability of the unusual single-layer beta-sheet.     

Refined solution structure and backbone dynamics of the archaeal MC1 protein

The 3D structure of methanogen chromosomal protein 1 (MC1), determined with heteronuclear NMR methods, agrees with its function in terms of the shape and nature of the binding surface, whereas the 3D structure determined with homonuclear NMR does not. The structure features five loops, which show a large distribution in the ensemble of 3D structures. Evidence for the fact that this distribution signifies internal mobility on the nanosecond time scale was provided by using (15)N-relaxation and molecular dynamics simulations. Structural variations of the arm (11 residues) induced large shape anisotropy variations on the nanosecond time scale that ruled out the use of the model-free formalism to analyze the relaxation data. The backbone dynamics analysis of MC1 was achieved by comparison with 20 ns molecular dynamics trajectories. Two β-bulges showed that hydrogen bond formation correlated with ? and ψ dihedral angle transitions. These jumps were observed on the nanosecond time scale, in agreement with a large decrease in (15)N-NOE for Gly17 and Ile89. One water molecule bridging NH(Glu87) and CO(Val57) through hydrogen bonding contributed to these dynamics. Nanosecond slow motions observed in loops LP3 (35-42) and LP5 (67-77) reflected the lack of stable hydrogen bonds, whereas the other loops, LP1 (10-14), LP2 (22-24), and LP4 (50-53), were stabilized by several hydrogen bonds. Dynamics are often directly related to function. Our data strongly suggest that residues belonging to the flexible regions of MC1 could be involved in the interaction with DNA.     

Fast internal main-chain dynamics of human ubiquitin.

The fast internal dynamics of human ubiquitin have been studied by the analysis of 15N relaxation of backbone amide nitrogens. The amide 15N resonances have been assigned by use of heteronuclear multiple-quantum spectroscopy. Spin lattice relaxation times at 60.8 and 30.4 MHz and the steady-state nuclear Overhauser effect at 60.8 MHz have been determined for 67 amide 15N sites in the protein using two-dimensional spectroscopy. These data have been analyzed in terms of the model free treatment of Lipari and Szabo [Lipari, G., & Szabo, A. (1982) J. Am. Chem. Soc. 104, 4546-4559]. The global motion of the protein is shown to be isotropic and is characterized by a correlation time of 4.1 ns rad-1. The generalized order parameters (S2) of backbone amide N-H vectors in the globular region of the protein range from 0.5 to 0.95. No apparent correlation between secondary structure and generalized order parameters is observed. There is, however, a strong correlation between the magnitude of the generalized order parameters of a given N-H vector and the presence of hydrogen bonding of the amide hydrogen or its peptide bond associated carbonyl. Using a chemical shift tensor breadth of 160 ppm, the N-H vectors of peptide linkages participating in one or more hydrogen bonds to the main chain show an average generalized order parameter of 0.80 (SD 0.06), while those amide NH of peptide linkages free of hydrogen-bonding interactions with the main chain show an average order parameter of 0.69 (SD 0.06).(ABSTRACT TRUNCATED AT 250 WORDS)     

Accuracy and precision of NMR relaxation experiments and MD simulations for characterizing protein dynamics

Model-free parameters obtained from nuclear magnetic resonance (NMR) relaxation experiments and molecular dynamics (MD) simulations commonly are used to describe the intramolecular dynamical properties of proteins. To assess the relative accuracy and precision of experimental and simulated model-free parameters, three independent data sets derived from backbone ¹⁵N NMR relaxation experiments and two independent data sets derived from MD simulations of Escherichia coli ribonuclease HI are compared. The widths of the distributions of the differences between the order parameters for pairs of NMR data sets are congruent with the uncertainties derived from statistical analyses of individual data sets; thus, current protocols for analyzing NMR data encapsulate random uncertainties appropriately. Large differences in order parameters for certain residues are attributed to systematic differences between samples for intralaboratory comparisons and unknown, possibly magnetic field-dependent, experimental effects for interlaboratory comparisons. The widths of distributions of the differences between the order parameters for two NMR sets are similar to widths of distributions for an NMR and an MD set or for two MD sets. The linear correlations between the order parameters for an MD set and an NMR set are within the range of correlations observed between pairs of NMR sets. These comparisons suggest that the NMR and MD generalized order parameters for the backbone amide N—H bond vectors are of comparable accuracy for residues exhibiting motions on a fast time scale (<100 ps). Large discrepancies between NMR and MD order parameters for certain residues are attributed to the occurrence of “rare” motional events over the simulation trajectories, the disruption of an element of secondary structure in one of the simulations, and lack of consensus among the experimental data sets. Consequently, (easily detectable) severe distortions of local protein structure and infrequent motional events in MD simulations appear to be the most serious artifacts affecting the accuracy and precision, respectively, of MD order parameters relative to NMR values. In addition, MD order parameters for motions on a fast (<100 ps) timescale are more precisely determined than their NMR counterparts, thereby permitting more detailed dynamic characterization of biologically important residues by MD simulation than is sometimes possible by experimental methods. Proteins 28:481–493, 1997. © 1997 Wiley-Liss, Inc.     

Intrinsic dynamics of the partly unstructured PX domain from the Sendai virus RNA polymerase cofactor P

Despite their evident importance for function, dynamics of intrinsically unstructured proteins are poorly understood. Sendai virus phosphoprotein, cofactor of the RNA polymerase, contains a partly unstructured protein domain. The phosphoprotein X domain (PX) is responsible for binding the polymerase to the nucleocapsid assembling the viral RNA. For RNA synthesis, the interplay of the dynamics of the unstructured and structured PX subdomains is thought to drive progression of the RNA polymerase along the nucleocapsid. Here we present a detailed study of the dynamics of PX using hydrogen/deuterium exchange and different NMR relaxation measurements. In the unstructured subdomain, large amplitude fast motions were found to be fine-tuned by the presence of residues with short side chains. In the structured subdomain, where fast motions of both backbone and side chains are fairly restricted, the first helix undergoes slow conformational exchange corresponding to a local unfolding event. The other two helices, which represent the nucleocapsid binding site, were found to be more stable and to reorient with respect to each other, as probed by slow conformational exchange identified for residues on the third helix. The study illustrates the intrinsically differential dynamics of this partly unstructured protein and proposes the relation between these dynamics and its function.     

Characterization of NMR relaxation-active motions of a partially folded A-state analogue of ubiquitin

The dominant dynamics of a partially folded A-state analogue of ubiquitin that give rise to NMR 15N spin relaxation have been investigated using molecular dynamics (MD) computer simulations and reorientational quasiharmonic analysis. Starting from the X-ray structure of native ubiquitin with a protonation state corresponding to a low pH, the A-state analogue was generated by a MD simulation of a total length of 33 ns in a 60%/40% methanol/water mixture using a variable temperature scheme to control and speed up the structural transformation. The N-terminal half of the A-state analogue consists of loosely coupled native-like secondary structural elements, while the C-terminal half is mostly irregular in structure. Analysis of dipolar N-H backbone correlation functions reveals reorientational amplitudes and time-scale distributions that are comparable to those observed experimentally. Thus, the trajectory provides a realistic picture of a partially folded protein that can be used for gaining a better understanding of the various types of reorientational motions that are manifested in spin-relaxation parameters of partially folded systems. For this purpose, a reorientational quasiharmonic reorientational analysis was performed on the final 5 ns of the trajectory of the A-state analogue, and for comparison on a 5 ns trajectory of native ubiquitin. The largest amplitude reorientational modes show a markedly distinct behavior for the two states. While for native ubiquitin, such motions have a more local character involving loops and the C-terminal end of the polypeptide chain, the A-state analogue shows highly collective motions in the nanosecond time-scale range corresponding to larger-scale movements between different segments. Changes in reorientational backbone entropy between the A-state analogue and the native state of ubiquitin, which were computed from the reorientational quasiharmonic analyses, are found to depend significantly on motional correlation effects.     

Backbone dynamics of Escherichia coli thioesterase/protease I: evidence of a flexible active-site environment for a serine protease

Escherichia coli thioesterase/protease I (TEP-I) is a member of a novel subclass of the lipolytic enzymes with a distinctive GDSLS motif. In addition to possessing thioesterase and protease activities, TEP-I also exhibits arylesterase activity. We have determined the (15)N nuclear magnetic spin relaxation rates, R(1) and R(2), and the steady state (1)H-(15)N heteronuclear Overhauser effect, measured at both 11.74 T and 14.09 T, of (u-(15)N) TEP-I. These data were analyzed using model-free formalism (with axially symmetric rotational diffusion anisotropy) to extract the backbone dynamics of TEP-I. The results reveal that the core structure of the central beta-sheet and the long alpha-helices are rigid, while the binding pocket appears to be rather flexible. The rigid core serves as a scaffold to anchor the essential loops, which form the binding pocket. The most flexible residues display large amplitude fast (ps/ns time-scale) motion and lie on one stripe whose orientation is presumed to be the ligand-binding orientation. We also detected the presence of several residues displaying slow (microseconds/ms time-scale) conformational exchanging processes. These residues lie around the binding pocket and are oriented perpendicularly to the orientation of the flexible stripe. Two of the putative catalytic triads, Ser10 and His157, and their neighbors show motion on the microseconds/ms time-scale, suggesting that their slow motion may have a role in catalysis, in addition to their possible roles in ligand binding. The presence of a flexible substrate-binding pocket may also facilitate binding to a wide range of substrates and confer the versatile functional property of this protein.     

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.     

Local segmental dynamics of cis-1,4-polybutadiene,polypropylene and polyethylene terephthalate via molecular dynamics simulations

The local segmental dynamics of cis-1,4-polybutadiene, polypropylene and polyethylene terephthalate have been investigated via isothermal-isobaric molecular dynamics simulations. The simulation pressure was 1 atm for all systems, with all simulation temperatures being at least 150 K above the polymer's glass transition temperature. The trajectories have been analysed via autocorrelation functions (ACFs) of chord vectors spanning different number of chain backbone bonds. Inverse Laplace transformations of these ACFs using the CONTIN algorithm afforded the corresponding distribution of relaxation times (DRTs) for the simulated dynamics. All DRTs illustrated a peak on fast timescales corresponding to short length scale segmental motion and a peak at longer timescales corresponding to longer length scale relaxations. A third peak, intermediate between the fast and slow processes, appears as the relaxation of chord vectors spanning increasing number of backbone bonds is considered. The temperature dependence of the relaxation dynamics is also investigated.     

Relationships between protein structure and dynamics from a database of NMR-derived backbone order parameters

The amplitude of protein backbone NH group motions on a time-scale faster than molecular tumbling may be determined by analysis of (15)N NMR relaxation data according to the Lipari-Szabo model free formalism. An internet-accessible database has been compiled containing 1855 order parameters from 20 independent NMR relaxation studies on proteins whose three-dimensional structures are known. A series of statistical analyses has been performed to identify relationships between the structural features and backbone dynamics of these proteins. Comparison of average order parameters for different amino acid types indicates that amino acids with small side-chains tend to have greater backbone flexibility than those with large side-chains. In addition, the motions of a given NH group are also related to the sizes of the neighboring amino acids in the primary sequence. The secondary structural environment appears to influence backbone dynamics relatively weakly, with only subtle differences between the order parameter distributions of loop structures and regular hydrogen bonded secondary structure elements. However, NH groups near helix termini are more mobile on average than those in the central regions of helices. Tertiary structure influences are also relatively weak but in the expected direction, with more exposed residues being more flexible on average than residues that are relatively inaccessible to solvent.     

Main chain and side chain dynamics of oxidized flavodoxin from Cyanobacterium anabaena.

Oxidized flavodoxin from Cyanobacterium anabaena PCC 7119 is used as a model system to investigate the fast internal dynamics of a flavin-bearing protein. Virtually complete backbone and side chain resonance NMR assignments of an oxidized flavodoxin point mutant (C55A) have been determined. Backbone and side chain dynamics in flavodoxin (C55A) were investigated using (15)N amide and deuterium methyl NMR relaxation methods. The squared generalized order parameters (S(NH)(2)) for backbone amide N-H bonds are found to be uniformly high ( approximately 0.923 over 109 residues in regular secondary structure), indicating considerable restriction of motion in the backbone of the protein. In contrast, methyl-bearing side chains are considerably heterogeneous in their amplitude of motion, as indicated by obtained symmetry axis squared generalized order parameters (S(axis)(2)). However, in comparison to nonprosthetic group-bearing proteins studied with these NMR relaxation methods, the side chains of oxidized flavodoxin are unusually rigid.     

Concordance of residual dipolar couplings, backbone order parameters and crystallographic B-factors for a small alpha/beta protein: a unified picture of high probability, fast atomic motions in proteins

Using ensemble refinement of the third immunoglobulin binding domain (GB3) of streptococcal protein G (a small alpha/beta protein of 56 residues), we demonstrate that backbone (N-H, N-C', Calpha-Halpha, Calpha-C') residual dipolar coupling data in five independent alignment media, generalized order parameters from 15N relaxation data, and B-factors from a high-resolution (1.1A), room temperature crystal structure are entirely consistent with one another within experimental error. The optimal ensemble size representation is between four and eight, as assessed by complete cross-validation of the residual dipolar couplings. Thus, in the case of GB3, all three observables reflect the same low-amplitude anisotropic motions arising from fluctuations in backbone phi/psi torsion angles in the picosecond to nanosecond regime in both solution and crystalline environments, yielding a unified picture of fast, high-probability atomic motions in proteins. An understanding of these motions is crucial for understanding the impact of protein dynamics on protein function, since they provide part of the driving force for triggered conformational changes that occur, for example, upon ligand binding, signal transduction and enzyme catalysis.