Microsecond timescale backbone conformational dynamics in ubiquitin studied with NMR R1rho relaxation experiments

NMR spin relaxation experiments are used to characterize the dynamics of the backbone of ubiquitin. Chemical exchange processes affecting residues Ile 23, Asn 25, Thr 55, and Val 70 are characterized using on- and off-resonance rotating-frame 15N R1rho relaxation experiments to have a kinetic exchange rate constant of 25,000 sec(-1) at 280 K. The exchange process affecting residues 23, 25, and 55 appears to result from disruption of N-cap hydrogen bonds of the alpha-helix and possibly from repacking of the side chain of Ile 23. Chemical exchange processes affecting other residues on the surface of ubiquitin are identified using 1H-15N multiple quantum relaxation experiments. These residues are located near or at the regions known to interact with various enzymes of the ubiquitin-dependent protein degradation pathway.     

NMR studies of structure, hydrogen exchange, and main-chain dynamics in a disrupted-core mutant of thioredoxin.

Core-packing mutants of proteins often approach molten globule states, and hence may have attributes of folding intermediates. We have studied a core-packing mutant of thioredoxin, L78K, in which a leucine residue is substituted by lysine, using 15N heteronuclear two- and three-dimensional NMR. Chemical shift differences between the mutant and wild-type main-chain resonances reveal that structural changes caused by the mutation are localized within 12 A of the altered side chain. The majority of resonances are unchanged, as are many 1H-1H NOEs indicative of the main-chain fold, suggesting that the structure of L78K is largely similar to wild type. Hydrogen exchange studies reveal that residues comprising the central beta-sheet of both mutant and wild-type proteins constitute a local unfolding unit, but with the unfolding/folding equilibrium approximately 12 times larger in L78K. The dynamics of main-chain NH bonds in L78K were studied by 15N spin relaxation and compared with a previous study of wild type. Order parameters for angular motion of NH bonds in the mutant are on average lower than in wild type, suggesting greater spatial freedom on a rapid time scale, but may also be related to different rotational correlation times in the two proteins. There is also evidence of greater conformational exchange in the mutant. Differences between mutant and wild type in hydrogen exchange and main-chain dynamics are not confined to the vicinity of the mutation. We infer that mispacking of the protein core in one location affects local dynamics and stability throughout.     

Automated NMR determination of protein backbone dihedral angles from cross-correlated spin relaxation

The simultaneous interpretation of a suite of dipole-dipole and dipole-CSA cross-correlation rates involving the backbone nuclei ¹³C, ¹H,¹³CO, ¹⁵N and ¹H^N can be used to resolve the ambiguities associated with each individual cross-correlation rate. The method is based on the transformation of experimental cross-correlation rates via calculated values based on standard peptide plane geometry and solid-state ¹³CO CSA parameters into a dihedral angle probability surface. Triple resonance NMR experiments with improved sensitivity have been devised for the quantification of relaxation interference between ¹H(i)-¹³C(i)/¹⁵N(i)-¹H^N(i) and ¹H(i–1)-¹³C(i–1)/¹⁵N(i)-¹H^N(i) dipole-dipole mechanisms in ¹⁵N,¹³C-labeled proteins. The approach is illustrated with an application to ¹³C,¹⁵N-labeled ubiquitin.     

NMR investigation of the dynamics of tryptophan side-chains in hemoglobins

NMR relaxation measurements of 15N spin-lattice relaxation rate (R(1)), spin-spin relaxation rate (R(2)), and heteronuclear nuclear Overhauser effect (NOE) have been carried out at 11.7T and 14.1T as a function of temperature for the side-chains of the tryptophan residues of 15N-labeled and/or (2H,15N)-labeled recombinant human normal adult hemoglobin (Hb A) and three recombinant mutant hemoglobins, rHb Kempsey (betaD99N), rHb (alphaY42D/betaD99N), and rHb (alphaV96W), in the carbonmonoxy and the deoxy forms as well as in the presence and in the absence of an allosteric effector, inositol hexaphosphate (IHP). There are three Trp residues (alpha14, beta15, and beta37) in Hb A for each alphabeta dimer. These Trp residues are located in important regions of the Hb molecule, i.e. alpha14Trp and beta15Trp are located in the alpha(1)beta(1) subunit interface and beta37Trp is located in the alpha(1)beta(2) subunit interface. The relaxation experiments show that amino acid substitutions in the alpha(1)beta(2) subunit interface can alter the dynamics of beta37Trp. The transverse relaxation rate (R(2)) for beta37Trp can serve as a marker for the dynamics of the alpha(1)beta(2) subunit interface. The relaxation parameters of deoxy-rHb Kemspey (betaD99N), which is a naturally occurring abnormal human hemoglobin with high oxygen affinity and very low cooperativity, are quite different from those of deoxy-Hb A, even in the presence of IHP. The relaxation parameters for rHb (alphaY42D/betaD99N), which is a compensatory mutant of rHb Kempsey, are more similar to those of Hb A. In addition, TROSY-CPMG experiments have been used to investigate conformational exchange in the Trp residues of Hb A and the three mutant rHbs. Experimental results indicate that the side-chain of beta37Trp is involved in a relatively slow conformational exchange on the micro- to millisecond time-scale under certain experimental conditions. The present results provide new dynamic insights into the structure-function relationship in hemoglobin.     

Solution conformation and dynamics of a fungal cell wall polysaccharide isolated from Microsporum gypseum

The conformational and dynamical features of a branched mannan isolated from a fungal cell wall have been analysed by homo and heteronuclear NMR methods, employing different magnetic fields. 1HNMR cross relaxation times have been obtained for this polysaccharide and have been interpreted qualitatively using different motional models. 13C NMR relaxation parameters (T1, T2, NOE) have also been measured and interpreted using different approximations based on the Lipari and Szabo model free approach. The analysis of the data indicate the existence of important flexibility for the different linkages of the polysaccharide. Motions in the range of 4–6 ns contribute to the relaxation of the macromolecule, although faster internal motions in the 500 ps and 100 ps timescales are also present. These time scales indicate that segmental motions as well as internal motions around the glycosidic linkages are the major sources of relaxation for this molecule at 318 K. Molecular dynamics simulations have also been performed. The obtained results also indicate that the polysaccharide possess a substantial amount of conformational freedom.     

Identification of Slow Motions in the Reduced Recombinant High-potential Iron Sulfur Protein I (HiPIP I) from Ectothiorhodospira Halophila via 15N Rotating-frame NMR Relaxation Measurements

Rotating-frame 15N relaxation rate (R1) NMR experiments have been performed in order to study the dynamic behavior of the reduced recombinant high-potential iron-sulfur protein iso I (HiPIP I) from Ectothiorhodospira halophila, in the s to ms time range. Measurements of R1 were performed as a function of the effective spin-lock magnetic field amplitude by using both on and off-resonance radio frequency irradiation. The two data sets provided consistent results and were fit globally in order to identify possible exchange processes in an external loop of the reduced HiPIP I. The loop consists of residues 43-45 and the correlation time of the exchange process was determined to be 50 ± 8 s for the backbone nitrogen of Gln 44.     

Estimating the time scale of chemical exchange of proteins from measurements of transverse relaxation rates in solution

Chemical (conformational) exchange on the ms-s time scale is reliably identified by the observation of transverse relaxation rates, R_ex, that depend upon the strength of the effective field (_1eff=B_1eff) used in spin lock or CPMG experiments. In order to determine if the exchange correlation time, _ex, is the fast or slow limit, measurements of (i) signal line shape and (ii) temperature dependence of R_ex have been commonly used in studies of stable, small molecules. However, these approaches are often not applicable to proteins, because sample stability and solubility, respectively, limit the temperature range and signal sensitivity of experiments. Herein we use a complex, but general, two-site exchange equation to show when the simple fast exchange equations for R_ex are good approximations, in the case of proteins. We then present a simple empirical equation that approximately predicts R_ex in all exchange regimes, and explains these results in a clear, straightforward manner. Finally we show how one can reliably determine whether _ex is in the fast or slow exchange limit.     

Comparison of the dynamics of the membrane-bound form of fd coat protein in micelles and in bilayers by solution and solid-state nitrogen-15 nuclear magnetic resonance spectroscopy

Solid-state and solution 15N nuclear magnetic resonance experiments on uniformly and specifically 15N labeled coat protein in phospholipid bilayers and in detergent micelles are used to describe the dynamics of the membrane-bound form of the protein. The residues in the N- and C-terminal portions of the coat protein in both phospholipid bilayers and in detergent micelles are mobile, while those in the hydrophobic midsection are immobile. There is evidence for a gradient of mobility in the C-terminal region of the coat protein in micelles; at 25 degrees C only the last two residues are mobile on the 10(9)-Hz timescale, while the last six to eight residues appear to be mobile on slower timescales and highly mobile at higher temperatures. Since all of the C-terminal residues are immobile in the virus particles, the mobility of these residues in the membrane-bound form of the protein may be important for the formation of protein-DNA interactions in the assembly process.     

NMR solution structure and backbone dynamics of domain III of the E protein of tick-borne Langat flavivirus suggests a potential site for molecular recognition

Flaviviruses cause many human diseases, including dengue fever, yellow fever, West Nile viral encephalitis, and hemorrhagic fevers, and are transmitted to their vertebrate hosts by infected mosquitoes and ticks. Domain III of the envelope protein (E-D3) is considered to be the primary viral determinant involved in the virus-host-cell receptor interaction, and thus represents an excellent target for antiviral drug development. Langat (LGT) virus is a naturally attenuated BSL-2 TBE virus and is a model for the pathogenic BSL-3 and BSL-4 viruses in the serogroup. We have determined the solution structure of LGT-E-D3 using heteronuclear NMR spectroscopy. The backbone dynamics of LGT-E-D3 have been investigated using 15N relaxation measurements. A detailed analysis of the solution structure and dynamics of LGT-E-D3 suggests potential residues that could form a surface for molecular recognition, and thereby represent a target site for antiviral therapeutics design.     

Characterizing and controlling the inherent dynamics of cyclophilin-A

With the recent advances in NMR relaxation techniques, protein motions on functionally important timescales can be studied at atomic resolution. Here, we have used NMR-based relaxation experiments at several temperatures and both 600 and 900 MHz to characterize the inherent dynamics of the enzyme cyclophilin-A (CypA). We have discovered multiple chemical exchange processes within the enzyme that form a “dynamic continuum” that spans 20–30 Å comprising active site residues and residues proximal to the active site. By combining mutagenesis with these NMR relaxation techniques, a simple method of counting the dynamically sampled conformations has been developed. Surprisingly, a combination of point mutations has allowed for the specific regulation of many of the exchange processes that occur within CypA, suggesting that the dynamics of an enzyme may be engineered.     

NMR solution structure and dynamics of motilin in isotropic phospholipid bicellar solution

The structure and dynamics of the gastrointestinal peptide hormone motilin, consisting of 22 amino acid residues, have been studied in the presence of isotropic q=0.5 phospholipid bicelles. The NMR solution structure of the peptide in acidic bicelle solution was determined from 203 NOE-derived distance constraints and six backbone torsion angle constraints. Dynamic properties for the ¹³C-¹H vector in Leu10 were determined for motilin specifically labeled with ¹³C at this position by analysis of multiple-field relaxation data. The structure reveals an ordered -helical conformation between Glu9 and Lys20. The N-terminus is also well structured with a turn resembling that of a classical -turn. The ¹³C dynamics clearly show that motilin tumbles slowly in solution, with a correlation time characteristic of a large object. It was also found that motilin has a large degree of local flexibility as compared with what has previously been reported in SDS micelles. The results show that motilin interacts with the bicelle, displaying motional properties of a peptide bound to a membrane. In comparison, motilin in neutral bicelles seems less structured and more flexible. This study shows that the small isotropic bicelles are well suited for use as membrane-mimetic for structural as well as dynamical investigations of membrane-bound peptides by high-resolution NMR.     

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.     

Comparison of backbone dynamics of monomeric and domain-swapped stefin A

Three-dimensional domain swapping has been observed in increasing number of proteins and has been implicated in the initial stages of protein aggregation, including that of the cystatins. Stefin A folds as a monomer under native conditions, while under some denaturing conditions domain-swapped dimer is formed. We have determined the backbone dynamics of the monomeric and domain-swapped dimeric forms of stefin A by (15)N relaxation using a model-free approach. The overall correlation times of the molecules were determined to be 4.6 +/- 0.1 ns and 9.2 +/- 0.2 ns for the monomer and the dimer, respectively. In the monomer, decreased order parameters indicate an increased mobility for the N-terminal trunk, the first and the second binding loops. At the opposite side of the molecule, the loop connecting the alpha-helix with strand B, the beginning of strand B and the loop connecting strands C and D show increased localized mobility. In the domain-swapped dimer, a distinctive feature of the structure is the concatenation of strands B and C into a single long beta-strand. The newly formed linker region between strands B and C, which substitutes for the first binding loop in the monomer, has order parameters typical for the remainder of the beta-strands. Thus, the interaction between subunits that occurs on domain-swapping has consequences for the dynamics of the protein at long-range from the site of conformational change, where an increased rigidity in the newly formed linker region is accompanied by an increased mobility of loops remote from that site.     

Protein dynamics studied by rotating frame 15N spin relaxation times

Summary Conformational rate processes in aqueous solutions of uniformly ¹⁵N-labeled pancreatic trypsin inhibitor (BPTI) at 36°C were investigated by measuring the rotating frame relaxation times of the backbone ¹⁵N spins as a function of the spin-lock power. Two different intramolecular exchange processes were identified. A first local rate process involved the residues Cys³⁸ and Arg³⁹, had a correlation time of about 1.3 ms, and was related to isomerization of the chirality of the disulfide bond Cys¹⁴-Cys³⁸. A second, faster motional mode was superimposed on the disulfide bond isomerization and was tentatively attributed to local segmental motions in the polypeptide sequence-Cys¹⁴-Ala¹⁵-Lys¹⁶-. The correlation time for the overall rotational tumbling of the protein was found to be 2 ns, using the assumption that relaxation is dominated by dipolar coupling and chemical shift anistropy modulated by isotropic molecular reorientation.Abbreviations BPTI basic pancreatic trypsin inhibitor - 2D two-dimensional - COSY 2D correlation spectroscopy - TOCSY 2D total correlation spectroscopy - RF radio frequency - CW continuous wave - TPPI time-proportional phase incrementation - CSA chemical shift anisotropy - T₁ longitudinal relaxation time - T₂ transverse relaxation time - T₁ relaxation time in the rotating frame , correlation time for overall rotational reorientation of the protein - _ex ^s , _ex ^f , correlation times for two conformational exchange processes (slow and fast).     

Refined solution structure and backbone dynamics of 15N-labeled C12A-p8MTCP studied by NMR relaxation

MTCP1 (for Mature-T-Cell Proliferation) was the first gene unequivocally identified in the group of uncommon leukemias with a mature phenotype. The three-dimensional solution structure of the human p8^MTCP protein encoded by the MTCP1 oncogene has been previously determined by homonuclear proton two-dimensional NMR methods at 600 MHz: it consists of an original scaffold comprising three -helices, associated with a new cysteine motif. Two of the helices are covalently paired by two disulfide bridges, forming an -hairpin which resembles an antiparallel coiled-coil. The third helix is orientated roughly parallel to the plane defined by the -antiparallel motif and appears less well defined. In order to gain more insight into the details of this new scaffold, we uniformly labeled with nitrogen-15 a mutant of this protein (C12A-p8^MTCP1) in which the unbound cysteine at position 12 has been replaced by an alanine residue, thus allowing reproducibly high yields of recombinant protein. The refined structure benefits from 211 additional NOEs, extracted from ¹⁵N-edited 3D experiments, and from a nearly complete set of angular restraints allowing the estimation of the helical content of the structured part of the protein. Moreover, measurements of ¹⁵ N spin relaxation times and heteronuclear ¹⁵ N¹HNOEs provided additional insights into the dynamics of the protein backbone. The analysis of the linear correlation between J(0) and J() was used to interpret relaxation parameters. It appears that the apparent relative disorder seen in helix III is not simply due to a lack of experimental constraints, but associated with substantial contributions of sub-nanosecond motions in this segment.     

Insights into the local residual entropy of proteins provided by NMR relaxation.

A simple model is used to illustrate the relationship between the dynamics measured by NMR relaxation methods and the local residual entropy of proteins. The expected local dynamic behavior of well-packed extended amino acid side chains are described by employing a one-dimensional vibrator that encapsulates both the spatial and temporal character of the motion. This model is then related to entropy and to the generalized order parameter of the popular "model-free" treatment often used in the analysis of NMR relaxation data. Simulations indicate that order parameters observed for the methyl symmetry axes in, for example, human ubiquitin correspond to significant local entropies. These observations have obvious significance for the issue of the physical basis of protein structure, dynamics, and stability.     

Structural and dynamics characterization of norovirus protease

Norovirus protease is an essential enzyme for proteolytic maturation of norovirus nonstructural proteins and has been implicated as a potential target for antiviral drug development. Although X‐ray structural studies of the protease give us wealth of structural information including interactions of the protease with its substrate and dimeric overall structure, the role of protein dynamics in the substrate recognition and the biological relevance of the protease dimer remain unclear. Here we determined the solution NMR structure of the 3C‐like protease from Norwalk virus (NV 3CLpro), a prototype strain of norovirus, and analyzed its backbone dynamics and hydrodynamic behavior in solution. ¹⁵N spin relaxation and analytical ultracentrifugation analyses demonstrate that NV 3CLpro is predominantly a monomer in solution. Solution structure of NV 3CLpro shows significant structural variation in C‐terminal domain compared with crystal structures and among lower energy structure ensembles. Also, ¹⁵N spin relaxation and Carr–Purcell–Meiboom–Gill (CPMG)‐based relaxation dispersion analyses reveal the dynamic properties of residues in the C‐terminal domain over a wide range of timescales. In particular, the long loop spanning residues T123–G133 show fast motion (ps‐ns), and the residues in the bII–cII region forming the large hydrophobic pocket (S2 site) undergo conformational exchanges on slower timescales (μs–ms), suggesting their important role in substrate recognition.     

Accessibility of tobacco lipid transfer protein cavity revealed by 15N NMR relaxation studies and molecular dynamics simulations

Plant LTP1 are small helical proteins stabilized by four disulfide bridges and are characterized by the presence of an internal cavity, in which various hydrophobic ligands can be inserted. Recently, we have determined the solution structure of the recombinant tobacco LTP1_1. Unexpectedly, despite a global fold very similar to the structures already known for cereal seed LTP1, its binding properties are different: Tobacco LTP1_1 is able to bind only one monoacylated lipid, whereas cereal LTP1 can bind either one or two. The 3D structure of tobacco LTP1_1 revealed the presence of a hydrophobic cluster, not observed on cereal LTP1 structures, which may hinder one of the two entrances of the cavity defined for wheat LTP1. To better understand the mechanism of lipid entrance for tobacco LTP1_1 and to define the regions of the protein monitoring the accessibility of the cavity, we have complemented our structural data by the study of the internal dynamics of tobacco LTP1_1, using (15)N magnetic relaxation rate data and MD simulations at room and high temperatures. This work allowed us to define two regions of the protein experiencing the largest motions. These two regions delineate a portal that opens up during the simulation constituting a unique entrance of the hydrophobic cavity, in contrast with wheat LTP1 where two routes were detected. The hydrophobic interactions resulting from a few point mutations are strong enough to completely block the second portal so that the accessibility of the cavity is restricted to one entrance, explaining why this particular LTP1 binds only one lipid molecule.     

A solution NMR study of the binding kinetics and the internal dynamics of an HIV-1 protease-substrate complex

NMR studies of the binding of a substrate to an inactive HIV-1 protease construct, containing an active site mutation PR(D25N), are reported. Substrate titration measurements monitored by HSQC spectra and a (15)N-edited NOESY experiment show that the chromogenic substrate analog of the capsid/p2 cleavage site binds to PR(D25N) with an equilibrium dissociation constant, K(D), of 0.27 +/- 0.05 mM, and upper limits of the association and dissociation rate constants, 2 x 10(4) M(-1)s(-1) and 10 s(-1), respectively, at 20 degrees C, pH 5.8. This association rate constant is not in the diffusion limit, suggesting that association is controlled by a rare event, such as opening of the protease flaps. Analysis of (15)N relaxation experiments reveals a slight reduction of S(2) values in the flap region, indicating a small increase in the amplitude of internal motion on the sub-nsec timescale. In addition, several residues in the flap region are mobile on the conformational exchange timescale, msec-microsec. Flap dynamics of the protease-substrate complex are compared with those of protease-inhibitor complexes, and the implications of these results for substrate-binding models are discussed.     

Internal dynamics of d(CGCAAATTTGCG)2: a comparison of NMR relaxation measurements with a molecular dynamics simulation

We report the analysis of a 250 ps molecular dynamics simulation of the dodecamer d(CGCAAATTT-GCG)₂ immersed in a rectangular box of 3469 water molecules with 22 Na⁺ counterions. The internal dynamics of the molecule were investigated by studying the relevant autocorrelation functions related to the ¹³C-NMR relaxation parameters of the C1′-H1′ bonds of the sugar rings. The calculated effective correlation times τ _e (∼13 ps) and the order parameter S² (∼0.82) of the Lipari and Szabo formalism (Lipari and Szabo 1982a, b) are in satisfactory agreement with those determined previously by NMR (Gaudin et al. 1995, 1996). ¹H-¹H NOE buildups have also been measured experimentally and agree with those computed from the simulation. These results validate the simulation, and a more detailed analysis of the internal dynamics of the dodecamer was undertaken. Analysis of the distributions and of the autocorrelation functions of the glycosidic angle flucuations χ shows that the rotational motion of the sugar rings about their glycosidic bond conforms to a restricted diffusion mechanism. The amplitude of the motions and the diffusion constant are 20° and 17.10⁹ rad²s^–1 respectively. These values are in good agreement with ¹³C NMR data. Furthermore the simulation allows us to rule out another model also consistent with the experiment, consisting of a two-state jump between a syn and an anti conformation. Received: 19 November 1996 / Accepted: 17 March 1997