Magic-angle spinning NMR studies of molecular organization in multibilayers formed by 1-octadecanoyl-2-decanoyl-sn-glycero-3-phosphocholine.

Magic-angle spinning 1H and 13C nuclear magnetic resonance (NMR) have been employed to study 50%-by-weight aqueous dispersions of 1-octadecanoyl-2-decanoyl-sn-glycero-3-phosphocholine (C[18]:C[10]PC) and 1-octadecanoyl-2-d19-decanoyl-PC (C[18]:C[10]PC-d19), mixed-chain phospholipids which can form interdigitated multibilayers. The 1H NMR linewidth for methyl protons of the choline headgroup has been used to monitor the liquid crystalline-to-gel (LC-to-G) phase transition and confirm variations between freezing and melting temperatures. Both 1H and 13C spin-lattice relaxation times indicate unusual restrictions on segmental reorientation at megahertz frequencies for C(18):C(10)PC as compared with symmetric-chain species in the LC state; nevertheless each chemical moiety of the mixed-chain phospholipid exhibits motional behavior that may be classified as liquidlike. Two-dimensional nuclear Overhauser spectroscopy (NOESY) on C(18):C(10)PC and C(18):C(10)PC-d19 reveals cross-peaks between the omega-methyl protons of the C18 chain and the N-methyl protons of the phosphocholine headgroup, and several experimental and theoretical considerations argue against an interpretation based on spin diffusion. Using NMR relaxation times and NOESY connectivities along with a computational formalism for four-spin systems (Keepers, J. W., and T. L. James. 1984. J. Magn. Reson. 57:404-426), an estimate of 3.5 A is obtained for the average distance between the omega-methyl protons of the C18 chain and the N-methyl protons of the phosphocholine headgroup. This finding is consistent with a degree of interdigitation similar to that proposed for organized assemblies of gel-state phosphatidylcholine molecules with widely disparate acyl-chain lengths (Hui, S. W., and C.-H. Huang. 1986. Biochemistry. 25:1330-1335); however, acyl-chain bendback or other intermolecular interactions may also contribute to the NOESY results. For multibilayers of C(18):C(10)PC in the gel phase, 13C chemical-shift measurements indicate that trans conformers predominate along both acyl chains. 13C Spin-lattice relaxation times confirm the unusual motional restrictions noted in the LC state; nevertheless, 13C and 1H rotating-frame relaxation times indicate that the interdigitated arrangement enhances chain or bilayer motions which occur at mid-kilohertz frequencies.     

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.     

Repeat motions and backbone flexibility in designed proteins with different numbers of identical consensus tetratricopeptide repeats

The tetratricopeptide repeat (TPR) is a 34-residue helix-turn-helix motif that occurs as three or more tandem repeats in a wide variety of proteins. We have determined the repeat motions and backbone fluctuations of proteins containing two or three consensus TPR repeats (CTPR2 and CPTR3, respectively) using 15N NMR relaxation measurements. Rotational diffusion tensors calculated from these data for each repeat within each TPR protein indicate that there is a high degree of motional correlation between different repeats in the same protein. This is consistent with the prevailing view that repeat proteins, such as CTPR2 and CTPR3, behave as single cooperatively folded domains. The internal motions of backbone NH groups were determined using the Lipari-Szabo model-free formalism. For most residues, there was a clear separation between the influence of internal motion and the influence of global rotational tumbling on the observed magnetic relaxation. The local internal motions are highly restricted in most of the helical elements, with slightly greater flexibility in the linker elements. Comparisons between CTPR2 and CTPR3 indicate that an addition of a TPR repeat to the C-terminus (before the solvation helix) of CTPR2 slightly reduces the flexibility of the preceding helix.     

Tacrine derivatives-acetylcholinesterase interaction: 1H NMR relaxation study

Two acetylcholinesterase (AChE) inhibitors structurally related to Tacrine, 6-methoxytacrine (1a) and 9-heptylamino-6-methoxytacrine (1b), and their interaction with Electrophorus Electricus AChE were investigated. The complete assignment of the 1H and 13C NMR spectra of 1a and 1b was performed by mono-dimensional and homo- and hetero-correlated two-dimensional NMR experiments. This study was undertaken to elucidate the interaction modes between AChE and 1a and 1b in solution, using NMR. The interaction between the two inhibitors and AChE was studied by the analysis of the motional parameters non-selective and selective spin-lattice relaxation times, thereby allowing the motional state of 1a and 1b, both free and bound with AChE, to be defined. The relaxation data pointed out the ligands molecular moiety most involved in the binding with AChE. The relevant ligand/enzyme interaction constants were also evaluated for both compounds and resulted to be 859 and 5412M(-1) for 1a and1b, respectively.     

Probing protein side chain dynamics via 13C NMR relaxation

Protein side chain dynamics is associated with protein stability, folding, and intermolecular interactions. Detailed dynamics information is crucial for the understanding of protein function and biochemical and biophysical properties, which can be obtained using NMR relaxation techniques. In this review, (13)C relaxation of methine, methylene and methyl groups with and without (1)H decoupling are described briefly for a better understanding of how spin relaxation is associated with motional (dynamics) parameters. Developments in the measurement and interpretation of (13)C auto-relaxation and cross-correlated relaxation data are presented too. Finally, recent progress in the use of (13)C relaxation to probe the dynamics of protein side chains is detailed mainly for the dynamics of non-deuterated proteins on picoseconds-nanosecond timescales.     

13C NMR relaxation studies of RNA base and ribose nuclei reveal a complex pattern of motions in the RNA binding site for human U1A protein

The widespread importance of induced fit and order-disorder transition in RNA recognition by proteins and small molecules makes it imperative that RNA motional properties are characterized quantitatively. Until now, however, very few studies have been dedicated to the systematic characterization of RNA motion and to their changes upon protein or small-molecule binding. The U1A protein-RNA complexes provide some of the best-studied examples of the role of RNA motional changes upon protein binding. Here, we report (13)C NMR relaxation studies of base and ribose dynamics for the RNA internal loop target of human U1A protein located within the 3'-untranslated region (3'-UTR) of the mRNA coding for U1A itself. We also report the semi-quantitative analysis of both fast (nano- to picosecond) and intermediate (micro- to millisecond) motions for this paradigmatic RNA system. We measure (13)C T(1), T(1rho) and heteronuclear nuclear Overhauser effects (NOEs) for sugar and base nuclei, as well as the power dependence of T(1rho) at 500 MHz and 750 MHz, and analyze these results using the model-free formalism. The results provide a much clearer picture of the type of motions experienced by this RNA in the absence of the protein than was provided by the analysis of the structure based solely on NOEs and scalar couplings. They define a model where the RNA internal loop region "breathes" on a micro- to millisecond timescale with respect to the double-helical regions. Superimposed on this slower motion, the residues at the very tip of the loop undergo faster (nano- to picosecond) motions. We hypothesize that these motions allow the RNA to sample multiple conformations so that the protein can select a structure within the ensemble that optimizes intermolecular contacts.     

An NMR spectroscopy study of bendaline-albumin interactions

The complete assignment of the 1H and 13C NMR spectra of bendaline (BNDL) was performed by mono-dimensional and homo- and hetero-correlated two-dimensional NMR experiments. The interaction between bendaline and albumin was also studied by the analysis of the motional parameters spin-lattice relaxation times, allowing the motional state of the BNDL free and bound with albumin to be defined. In absence of albumin the indazolacetic and benzylic moieties are characterized by roughly the same mobility and by positive sigma (cross-relaxation rates) values. In the presence of the macromolecule, the indazolacetic and benzylic moieties and the lysine change their motional behaviour to different extents, as indicated by correlation times. Data obtained in absence and in presence of the protein show that the molecular moiety of the bendaline most involved in the binding with albumin is the fragment H-4 H-5. The binding constant was evaluated at 2.4x10(3)M(-1).     

Dynamics of glycolipids in the liquid-crystalline state. 2H NMR study.

It has been shown previously that two types of motion are adequate to describe the partially relaxed 2H NMR line shapes (inversion recovery experiment) for the backbone portion of the glycolipid 1,2-di-O-tetradecyl-3-O-(beta-D-glucopyranosyl)-sn-glycerol (beta-DTGL) in the highly ordered gel phase (Auger, M.A., D. Carrier, I.C.P. Smith, and H. C. Jarrell. 1990. J. Am. Chem. Soc. 112:1373-1381). This study extends the latter investigation to the more fluid liquid-crystalline phase, where more complex motions are anticipated. Analyses of the powder line shapes and oriented sample relaxation data for both the glycerol backbone and head group regions of this lipid have been performed. The dynamics of glycerol at the C3 position in the gel state have been described by large angle jumps about the C2-C3 bond with a correlation time in the fast-limit motional regime (omega o tau c much less than 1) and site populations 0.46, 0.34, and 0.20. The present data show that in the liquid-crystalline phase the internal jump rate is maintained, and two additional motions are necessary to describe the dependence of the relaxation rate on the orientation of the director with respect to the magnetic field direction. These are rotation about the molecular long axis with a correlation time in the slow-limit motional regime very near to the T1 minimum (omega o tau c approximately 0.65), and molecular fluctuations about the order director (modeled by a Maier-Saupe restoration potential). This treatment was also extended to the glucose head group where additional segmental motion about the glycosidic bond has been reported previously. While the two motions dominating relaxation at the glycerol C3 segment reproduce the general relaxation features of the glucose head group, the results suggest that additional motion about the glycosidic linkage must be present. This study is a stringent test of the motional model chosen earlier because relaxation data were obtained at two 2H NMR frequencies using two relaxation experiments (T1Z and T1Q) and two types of sample preparation (oriented and dispersed multibilayers). The results strongly uphold the choice of model and indicate the utility of both oriented samples and the T1Q experiment.     

The analysis of NMR relaxation data in terms of multiple internal motions

A novel formalism for estimating the complex motions of proteins and other flexible macromolecules from NMR relaxation measurements is applied to ¹³C NMR relaxation data on the Bovine Pancreatic Trypsin Inhibitor (M. W. 6,500). Six experimental parameters measured at two field strengths are accounted for by a minimum of three motions at each carbon group. Low frequency components make small but finite contribution to the relaxation of all resonances, suggesting a general low frequency distortion of the backbone. Rotational diffusion of the protein makes a relatively minor contribution to the relaxation process. For aliphatic groups, rotation of side chains dominates relaxation.     

Backbone dynamics of TEM-1 determined by NMR: evidence for a highly ordered protein

Backbone dynamics of TEM-1 beta-lactamase (263 amino acids, 28.9 kDa) were studied by 15N nuclear magnetic resonance relaxation at 11.7, 14.1, and 18.8 T. The high quality of the spectra allowed us to measure the longitudinal relaxation rate (R1), the transverse relaxation rate (R2), and the {1H}-15N NOE for up to 227 of the 250 potentially observable backbone amide groups. The model-free formalism was used to determine internal motional parameters using an axially anisotropic model. TEM-1 exhibits a small prolate axial anisotropy (D(parallel)/D(perpendicular) = 1.23 +/- 0.01) and a global correlation time (tau(m)) of 12.41 +/- 0.01 ns. The unusually high average generalized order parameter (S2) of 0.90 +/- 0.02 indicates that TEM-1 is one of the most ordered proteins studied by liquid-state NMR to date. Although the omega-loop has a high degree of order in the picosecond-to-nanosecond time scale (mean S2 value of 0.90 +/- 0.02), we observed the presence of microsecond-to-millisecond time scale motions for this loop, as for the vicinity of the active site. These motions could be relevant for the catalytic function of TEM-1. Amide exchange experiments were also performed, and several amide groups were not exchanged after 12 days, an indication that global motions in TEM-1 are also very limited. Although detailed dynamics characterization by NMR cannot be readily applied to TEM-1 in the presence of relevant substrates, the unusual picosecond-to-nanosecond dynamics behavior of TEM-1 presented here will be essential to the validation and improvement of future molecular dynamics simulations of TEM-1 in the presence of functionally relevant substrates.     

Propagation of experimental uncertainties using the Lipari-Szabo model-free analysis of protein dynamics

In this paper we make use of the graphical procedure previously described [Jin, D. et al. (1997) J. Am. Chem. Soc., 119, 6923–6924] to analyze NMR relaxation data using the Lipari-Szabo model-free formalism. The graphical approach is advantageous in that it allows the direct visualization of the experimental uncertainties in the motional parameter space. Some general rules describing the relationship between the precision of the relaxation measurements and the precision of the model-free parameters and how this relationship changes with the overall tumbling time (m) are summarized. The effect of the precision in the relaxation measurements on the detection of internal motions not close to the extreme narrowing limit is analyzed. We also show that multiple timescale internal motions may be obscured by experimental uncertainty, and that the collection of relaxation data at very high field strength can improve the ability to detect such deviations from the simple Lipari-Szabo model.     

A theoretical study of rotational diffusion models for rod-shaped viruses. The influence of motion on 31P nuclear magnetic resonance lineshapes and transversal relaxation.

Information about the interaction between nucleic acids and coat proteins in intact virus particles may be obtained by studying the restricted backbone dynamics of the incapsulated nucleic acids using 31P nuclear magnetic resonance (NMR) spectroscopy. In this article, simulations are carried out to investigate how reorientation of a rod-shaped virus particle as a whole and isolated nucleic acid motions within the virion influence the 31P NMR lineshape and transversal relaxation dominated by the phosphorus chemical shift anisotropy. Two opposite cases are considered on a theoretical level. First, isotropic rotational diffusion is used as a model for mobile nucleic acids that are loosely or partially bound to the protein coat. The effect of this type of diffusion on lineshape and transversal relaxation is calculated by solving the stochastic Liouville equation by an expansion in spherical functions. Next, uniaxial rotational diffusion is assumed to represent the mobility of phosphorus in a virion that rotates as a rigid rod about its length axis. This type of diffusion is approximated by an exchange process among discrete sites. As turns out from these simulations, the amplitude and the frequency of the motion can only be unequivocally determined from experimental data by a combined analysis of the lineshape and the transversal relaxation. In the fast motional region both the isotropic and the uniaxial diffusion model predict the same transversal relaxation as the Redfield theory. For very slow motion, transversal relaxation resembles the nonexponential relaxation as observed for water molecules undergoing translational diffusion in a magnetic field gradient. In this frequency region T2e is inversely proportional to the cube root of the diffusion coefficient. In addition to the isotropic and uniaxial diffusion models, a third model is presented, in which fast restricted nucleic acid backbone motions dominating the lineshape are superimposed on a slow rotation of the virion about its length axis, dominating transversal relaxation. In an accompanying article the models are applied to the 31P NMR results obtained for bacteriophage M13 and tobacco mosaic virus.     

Dynamics in a protein-lipid complex: nuclear magnetic resonance measurements on the headgroup of cardiolipin when bound to cytochrome c.

Deuterium and phosphorus nuclear magnetic resonance (NMR) has been used to investigate the dynamics of slow motional processes induced in bilayer cardiolipin upon binding with cytochrome c. 31P NMR line shapes suggest that protein binding induces less restricted, isotropic-like motions in the lipid phosphates within the ms time scale of this measurement. However, these motions impart rapid transverse relaxation to methylene deuterons adjacent to the phosphate in the lipid headgroup and so did not feature strongly in the NMR line shapes recorded from these nuclei by using the quadrupolar echo. Nonetheless, motional characteristics of the headgroup deuterons were accessible to a dynamic NMR approach using the Carr-Purcell-Meiboom-Gill multiple-pulse experiment. Compared to the well-studied case of deuterons in fatty acyl chains of bilayer phosphatidylcholine, the motions determining the 2H spin transverse relaxation in the headgroup of bilayer cardiolipin were much faster, having a lower limit in the 5-10 kHz range. On binding with cytochrome c, the T2 effecting motions in the cardiolipin headgroup became faster still, with rates comparable to the residual quadrupolar coupling frequency of the headgroup deuterons (approximately 25 kHz) and so coincided with the time scale for recording the quadrupolar echo (approximately 40 microseconds). It is concluded that the headgroup of cardiolipin does not exclusively report localized dynamic information but is particularly sensitive to collective motions occurring throughout the bilayer molecules. Although the rates of collective modes of motion may be dependent on the lipid type in pure lipid bilayers, these low-frequency fluctuations appear to occupy a similar dynamic range in a variety of lipid-protein systems, including the natural membranes.     

Anisotropic rotational diffusion in model-free analysis for a ternary DHFR complex

Model-free analysis has been extensively used to extract information on motions in proteins over a wide range of timescales from NMR relaxation data. We present a detailed analysis of the effects of rotational anisotropy on the model-free analysis of a ternary complex for dihydrofolate reductase (DHFR). Our findings show that the small degree of anisotropy exhibited by DHFR (D_||/D=1.18) introduces erroneous motional models, mostly exchange terms, to over 50% of the NH spins analyzed when isotropic tumbling is assumed. Moreover, there is a systematic change in S², as large as 0.08 for some residues. The significant effects of anisotropic rotational diffusion on model-free motional parameters are in marked contrast to previous studies and are accentuated by lowering of the effective correlation time using isotropic tumbling methods. This is caused by the preponderance of NH vectors aligned perpendicular to the principal diffusion tensor axis and is readily detected because of the high quality of the relaxation data. A novel procedure, COPED (COmparison of Predicted and Experimental Diffusion tensors) is presented for distinguishing genuine motions from the effects of anisotropy by comparing experimental relaxation data and data predicted from hydrodynamic analyses. The procedure shows excellent agreement with the slow motions detected from the axially symmetric model-free analysis and represents an independent procedure for determining rotational diffusion and slow motions that can confirm or refute established procedures that rely on relaxation data. Our findings show that neglect of even small degrees of rotational diffusion anisotropy can introduce significant errors in model-free analysis when the data is of high quality. These errors can hinder our understanding of the role of internal motions in protein function.     

NMR investigations of Mn(II)-leucine-enkephalin complexes in [2H6]-DMSO solution.

Structural and kinetic features of the Mn(II)-Leu-enkephalin binding equilibria were delineated by measuring 13C and 1H NMR spin-lattice relaxation rates. The temperature dependence of such rates showed that some carbons were experiencing slow exchange regimes such that kinetic parameters at room temperature could be calculated (k(off) = 1400 sec-1, delta H* = 12.0 kcal/mol, delta S* = -9.9 e.u.). The paramagnetic rates of fast exchanging carbons were interpreted by the Solomon-Bloembergen-Morgan theory to provide structural parameters. The terminal carboxyl and amino groups were shown to be the binding sites. The motional correlation time (tau c = 0.6 nsec at 298 K) was calculated by measuring selective and double-selective 1H spin-lattice relaxation rates for the free peptide. The number of coordinated ligands was evaluated by considering the distance of the Leu CO in the complex at 2.54 A, as shown by molecular models. Finally, carbon-Mn(II) distances were calculated and the molecular model of the 1:1 complex was built.     

Angular variances for internal bond rotations of side chains in GXG-based tripeptides derived from (13)C-NMR relaxation measurements: Implications to protein folding

The study of backbone and side-chain internal motions in proteins and peptides is crucial to having a better understanding of protein/peptide "structure" and to characterizing unfolded and partially folded states of proteins and peptides. To achieve this, however, requires establishing a baseline for internal motions and motional restrictions for all residues in the fully, solvent-exposed "unfolded state." GXG-based tripeptides are the simpliest peptides where residue X is fully solvent exposed in the context of an actual peptide. In this study, a series of GXG-based tripeptides has been synthesized with X being varied to include all twenty common amino acid residues. Proton-coupled and -decoupled (13)C-nmr relaxation measurements have been performed on these twenty tripeptides and various motional models (Lipari-Szabo model free approach, rotational anisotropic diffusion, rotational fluctuations within a potential well, rotational jump model) have been used to analyze relaxation data for derivation of angular variances and motional correlation times for backbone and side-chain chi(1) and chi(2) bonds and methyl group rotations. At 298 K, backbone motional correlation times range from about 50 to 85 ps, whereas side-chain motional correlation times show a much broader spread from about 18 to 80 ps. Angular variances for backbone phi,psi bond rotations range from 11 degrees to 23 degrees and those for side chains vary from 5 degrees to 24 degrees for chi(1) bond rotations and from 5 degrees to 27 degrees for chi(2) bond rotations. Even in these peptide models of the "unfolded state," side-chain angular variances can be as restricted as those for backbone and beta-branched (valine, threonine, and isoleucine) and aromatic side chains display the most restricted motions probably due to steric hinderence with backbone atoms. Comparison with motional data on residues in partially folded, beta-sheet-forming peptides indicates that side-chain motions of at least hydrophobic residues are less restricted in the partially folded state, suggesting that an increase in side-chain conformational entropy may help drive early-stage protein folding. Copyright 1999 John Wiley & Sons, Inc.     

Trifluoroethanol increases the stability of Delta(5)-3-ketosteroid isomerase. 15N NMR relaxation studies

In the equilibrium unfolding process of Delta(5)-3-ketosteroid isomerase from Pseudomonas testosteroni by urea, it was observed that the enzyme stability increases by 2.5 kcal/mol in the presence of 5% trifluoroethanol (TFE). To elucidate the increased enzyme stability by TFE, the backbone dynamics of Delta(5)-3-ketosteroid isomerase were studied in the presence and absence of 5% TFE by (15)N NMR relaxation measurements, and the motional parameters (S(2), tau(e), and R(ex)) were extracted from the relaxation data using the model-free formalism. The presence of 5% TFE causes little change or a slight increase in the order parameters (S(2)) for a number of residues, which are located mainly in the dimer interface region. However, the majority of the residues exhibit reduced order parameters in the presence of 5% TFE, indicating that high frequency (pico- to nanosecond) motions are generally enhanced by TFE. The results suggest that the entropy can be an important factor for the enzyme stability, and the increase in entropy by TFE is partially responsible for the increased stability of Delta(5)-3-ketosteroid isomerase.     

An NMR database for simulations of membrane dynamics

Computational methods are powerful in capturing the results of experimental studies in terms of force fields that both explain and predict biological structures. Validation of molecular simulations requires comparison with experimental data to test and confirm computational predictions. Here we report a comprehensive database of NMR results for membrane phospholipids with interpretations intended to be accessible by non-NMR specialists. Experimental 13C-1H and 2H NMR segmental order parameters (S(CH) or S(CD)) and spin-lattice (Zeeman) relaxation times (T(1Z)) are summarized in convenient tabular form for various saturated, unsaturated, and biological membrane phospholipids. Segmental order parameters give direct information about bilayer structural properties, including the area per lipid and volumetric hydrocarbon thickness. In addition, relaxation rates provide complementary information about molecular dynamics. Particular attention is paid to the magnetic field dependence (frequency dispersion) of the NMR relaxation rates in terms of various simplified power laws. Model-free reduction of the T(1Z) studies in terms of a power-law formalism shows that the relaxation rates for saturated phosphatidylcholines follow a single frequency-dispersive trend within the MHz regime. We show how analytical models can guide the continued development of atomistic and coarse-grained force fields. Our interpretation suggests that lipid diffusion and collective order fluctuations are implicitly governed by the viscoelastic nature of the liquid-crystalline ensemble. Collective bilayer excitations are emergent over mesoscopic length scales that fall between the molecular and bilayer dimensions, and are important for lipid organization and lipid-protein interactions. Future conceptual advances and theoretical reductions will foster understanding of biomembrane structural dynamics through a synergy of NMR measurements and molecular simulations.     

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.