Variation of the NH-C alpha-H coupling constant with dihedral angle in the NMR spectra of peptides

Proton magnetic resonance data and conformational calculations of a series of model compounds containing a NH-C^αH group substituted as in peptides have been used to generate a proton–proton coupling constant–dihedral angle relation for the peptide unit. For those substances used in which the dihedral angle about the N-C^α bond is not fixed, the angle distribution was calculated from conformational theory. Using eight examples in which the number of theoretical assumptions were least, the best values of the coefficients A, B, and C in the expression J(θ) = Acos²θ + B cosθ + Csin²θ were found by a least-squares procedure to be 7.9, ?1.55, and 1.35, respectively. This relation gives reasonable values for the dihedral angles ? in cyclic oligopeptide structures for which the availability of both NMR data and other structural information allow comparison. When applied to N-acetylamino acid N-methylamides having side chains extending beyond C^β, however, agreement with the calculated conformational distribution was found for Leu, Met, and Trp, but observed values of J were larger than expected for Val, He, Phe, and Tyr, These disagreements are considered to be the result of interactions not yet taken into account in the usual conformational calculations.     

Correlation of 2J couplings with protein secondary structure

Geminal two‐bond couplings (²J) in proteins were analyzed in terms of correlation with protein secondary structure. NMR coupling constants measured and evaluated for a total six proteins comprise 3999 values of ²J_CαN′, ²J_C′HN, ²J_HNCα, ²J_C′Cα, ²J_HαC′, ²J_HαCα, ²J_CβC′, ²J_N′Hα, ²J_N′Cβ, and ²J_N′C′, encompassing an aggregate 969 amino‐acid residues. A seamless chain of pattern comparisons across the spectrum datasets recorded allowed the absolute signs of all ²J coupling constants studied to be retrieved. Grouped by their mediating nucleus, C′, N′ or C^α, ²J couplings related to C′ and N′ depend significantly on ?,ψ torsion‐angle combinations. β turn types I, I′, II and II′, especially, can be distinguished on the basis of relative‐value patterns of ²J_CαN′, ²J_HNCα, ²J_C′HN, and ²J_HαC′. These coupling types also depend on planar or tetrahedral bond angles, whereas such dependences seem insignificant for other types. ²J_HαCβ appears to depend on amino‐acid type only, showing negligible correlation with torsion‐angle geometry. Owing to its unusual properties, ²J_CαN′ can be considered a “one‐bond” rather than two‐bond interaction, the allylic analog of ¹J_N′Cα, as it were. Of all protein J coupling types, ²J_CαN′ exhibits the strongest dependence on molecular conformation, and among the ²J types, ²J_HNCα comes second in terms of significance, yet was hitherto barely attended to in protein structure work. Proteins 2010. © 2009 Wiley‐Liss, Inc.     

Conformational dependence of the vicinal proton coupling constant for the Cα?Cβ bond in peptides

We use the nmr data concerning the C^αH? C^βH fragment in eight peptides with rigid side chains to parametrize a Karplus correlation between the vicinal proton J_αβ coupling constant and the dihedral angle θ. When considering molecules containing the fragment C^αH^α? C^βH^βH^β′, the three-dimensional structure of the model peptides does not need to be known with accurate precision, since each set of J_αβ and J_αβ′ coupling constants is then related to the coefficients of the Karplus equation. A good correlation is observed with the Karplus equation, which is in substantial agreement with the J_αβ coupling constants reported for rigid as well as rotating C^α? C^β bonds in peptides.     

One-bond and two-bond J couplings help annotate protein secondary-structure motifs: J-coupling indexing applied to human endoplasmic reticulum protein ERp18

NMR coupling constants, both direct one‐bond (¹J) and geminal two‐bond (²J), are employed to analyze the protein secondary structure of human oxidized ERp18. Coupling constants collected and evaluated for the 18 kDa protein comprise 1268 values of ¹J_CαHα, ¹J_CαCβ, ¹J_CαC′, ¹J_C′N′, ¹J_N′Cα, ¹J_N′HN, ²J_CαN′, ²J_HNCα, ²J_C′HN, and ²J_HαC′. Comparison with ¹J and ²J data from reference proteins and pattern analysis on a per‐residue basis permitted main‐chain φ,ψ torsion‐angle combinations of many of the 149 amino‐acid residues in ERp18 to be narrowed to particular secondary‐structure motifs. J‐coupling indexing is here being developed on statistical criteria and used to devise a ternary grid for interpreting patterns of relative values of J. To account for the influence of the varying substituent pattern in different amino‐acid sidechains, a table of residue‐type specific threshold values was compiled for discriminating small, medium, and large categories of J. For the 15‐residue insertion that distinguishes the ERp18 fold from that of thioredoxin, the J‐coupling data hint at a succession of five isolated Type‐I β turns at progressively shorter sequence intervals, in agreement with the crystal structure. Proteins 2011. © 2010 Wiley‐Liss, Inc.     

Selective distortion of the planar group CαC'(O)O to a chiral flat tetrahedron in the amino acid alanine

Based on a Cambridge Structural Database (CSD) search, a meta‐analysis of 116 structures of alanine H₃NC_αH(CH₃)C′(O)O and its derivatives H₃NC_αH(CH₃)C′(O)O(H/R/M), protonated, esterified, or coordinated at the carboxylic group, shows that in the first step of a chirality chain, the L configuration at C_α induces (M) and (P) conformations with respect to rotation around the central C′─C_α bond. In the second step, the (M) and (P) conformations selectively distort the planar carboxylic group C_αC'(O_cis)O_trans to asymmetric flat (R) and (S) tetrahedra. High diastereoselectivities are caused by the two players attraction N…O_cis and repulsion O_trans…C_Me, which work together in (L,M,R) configurations but against each other in (L,P,S) configurations.     

A dihedral angle-vicinal proton coupling constant correlation for the α–β bond of amino acid residues

A coupling constant-dihedral angle correlation for the H? Cα? Cβ? H system of amino acid residues in peptides has been derived from a set of model compounds covering the full range of dihedral angles. The expression obtained, J = 11.0 cos₂ θ ?1.4 cos θ + 1.6 sin₂θ, is close to those already used in pmr studies of peptide conformation, and provides a firmer foundation for them. A factor limiting the precision of this and other “Karplus relations” is illustrated.     

1H- and 13C-nmr investigations on cis–trans isomerization of proline peptide bonds and conformation of aromatic side chains in H-Trp-(Pro)n-Tyr-OH peptides

¹H and ¹³C high-resolution nmr spectra of cationic, zwitterionic, and anionic forms of the peptides: H-Trp-(Pro)_n-Tyr-OH, n = 0-5, and H-Trp-Pro-OCH₃ were obtained in D₂O solution. Analysis of H^α(Pro₁), H^α(Trp), C^γ(Pro), H^ε(Tyr), and H^δ(Trp) resonances provided evidence for the presence of two predominant backbone isomers: the all-trans one and another with the Trp-Pro peptide bond in cis conformation; the latter constituted about 0.8 molar fraction of the total peptide (n > 1) concentration. Relative content of these isomers varied in a characteristic way with the number of Pro residues and the ionization state of the peptides. The highest content of the cis (Trp-Pro) isomer, 0.74, was found in the anionic form of H-Trp-Pro-Tyr-OH; it decreased in the order of: anion ? zwitterion ≈ cation, and with the number of Pro residues to reach the value of 0.42 in the cationic form of H-Trp- (Pro)₅-Tyr-OH. Isomerization equilibria about Pro-Pro bond(s) were found to be shifted far (?0.9) in favor of the trans conformation. Interpretation of the measured vicinal coupling constants J_α?β′ and J_α?β″ for C^αH-C^βH₂ proton systems of Trp and Tyr side chains in terms of relative populations of g⁺, g^?, and t staggered rotamers around the χ₁ dihedral angle indicated that in all the peptides studied (a) rotation of Trp indole ring in cis (Trp-Pro) isomers is strongly restricted, and (b) rotation of Tyr phenol ring is relatively free. The most preferred χ₁ rotamer of Trp (0.8-0.9 molar fraction) was assigned as the t one on the basis of a large value of the vicinal coupling constant between the high-field H^β and carbonyl carbon atoms of Trp, estimated for the cis (Pro₁) form of H-Trp-Pro-Tyr-OH from a ¹H, ¹³C correlated spectroscopy ¹H detected multiple quantum experiment. This indicates that cis ? trans equilibrium in the Trp-Pro fragment is governed by nonbonding interactions between the pyrrolidine (Pro) and indole (Trp) rings. A molecular model of the terminal cis Trp-Pro dipeptide fragment is proposed, based on the presented nmr data and the results of our molecular mechanics modeling of low-energy conformers of the peptides, reported elsewhere. © 1993 John Wiley & Sons, Inc.     

Simultaneous measurement of protein one-bond and two-bond nitrogen-carbon coupling constants using an internally referenced quantitative J-correlated [(15)N,(1)H]-TROSY-HNC experiment

A quantitative J-correlation pulse sequence is described that allows simultaneous determination of one-bond and two-bond nitrogen-carbon coupling constants for protonated or deuterated proteins. Coupling constants are calculated from volume ratios between cross peaks and reference axial peaks observed in a single 3D spectrum. Accurate backbone ¹ J _NC, ¹ J _NC, and ² J _NC coupling constants are obtained for the two [¹⁵N;¹³C]-labeled, medium-sized proteins flavodoxin and xylanase and for the [²H;¹⁵N;¹³C]-labeled, large protein DFPase. A dependence of one-bond and two-bond J _NC values on protein backbone torsion angles is readily apparent, in agreement with previously found correlations. In addition, the experiment is performed on isotropic as well as aligned protein to measure associated ¹⁵N-¹³C residual dipolar couplings.     

Conformational study of ribonucleotides, 2′-deoxyribonucleotides,and arabinonucleotides by carbon-13 nuclear magnetic resonance

The geminal and vicinal ¹³C-³¹P coupling constants have been monitored, as a function of pH, for a series of uracil and cytosine 3′- and 5′-nucleotides with a ribose, arabinose, or 2′-deoxyribose sugar. Data were also obtained for two 3′,5′-diphosphates in the ribose and arabinose series. The geminal J(C5′-P5′) and J(C3′-P3′) couplings show only a small dependence on the ionization state of the phosphate, decreasing by < 0.5 Hz in the pH 5–7 range. For the ribose and arabinose 3′-nucleotides, the vicinal J(C4′-P3′) increase (up to 1.5 Hz) on secondary phosphate ionization in the pH 5–7 range, whereas their J(C2′-P3′) couplings decrease (up to 1.5 Hz) over the same pH range. In contrast for the 2′-deoxyribose molecules, both couplings decrease (～0.5 Hz) on phosphate ionization. The titration curves provide information about the influence of the sugar on the conformation about the C3′? O3′ bond. Some conformational trends could be rationalized by consideration of the sugar-puckerdependent contact interactions between the 3′-phosphate and the substituents on the furanose ring.     

Quantitative f torsion angle analysis in Desulfovibrio vulgaris flavodoxin based on six f related 3 J couplings

Quantitative φ-dihedral angle determinations of non-glycine and non-proline residues in Desulfovibrio vulgaris flavodoxin are carried out on the exclusive basis of ³ J coupling constants. In total 124 ³ JHNH _α , 123 ³ JHNC _′i , 118 ³ JHNC _β , 117 ³ JC′ _i–1Hα , 109 ³ JC′ _i–1C′i , and 103 ³ JC′ _i–1Hβ values form the experimental basis for translating J coupling data into geometry information using various combinations of Karplus parameters given in the literature. In addition, each backbone torsional angle φ is adjusted assuming different models of local geometry, either a rigid torsion, a Gaussian distribution centered at a distinct angle, or a two-site jump model. Numerical optimization is followed by a statistical significance evaluation to assess the results. It is found that experimental coupling constants of most of the residues involved in secondary structure elements agree best with those predicted from rigid local conformations. For dihedral angles in loop regions, mobility effects are not negligible, and a single torsion (Glu 42) is likely to adopt two distinct adjustments. However, α-helix conformations with –60° < φ < –45° give rise to an alternate solution with φ≈+170° with similar statistical significance when using the four traditionally determined proton-involved ³ J couplings. This ambiguity is efficiently avoided only when taking advantage of the complete data set comprising six available experimental ³ J coupling constants and of the degeneracy intrinsic to the Karplus relation. The optimized φ conformations are compared with reference values from the crystal structure of flavodoxin.     

Spatial configuration of single-stranded polynucleotides. Calculations of average dimensions and Nmr coupling constants

The probability distributions of the torsional angles (Φ′, ω′, ω, Φ, and ψ), which fix the structure of nucleotide backbone, have been calculated using the results of energy calculations based on extended Huckel theory (EHT), complete neglect of differential overlap (CNDO), perturbative configuration interaction using localized orbitals (PCILO), and classical potential functions (CPF) methods. Statistical average values of the vicinal ¹H? ¹H, ¹H? ³¹P, and ¹³C? ³¹P nmr coupling constants 〈J〉 have been calculated from the generalized Karplus relations using the probability distribution in the Φ′, Φ, and ψ space. Experimental 〈J〉 values for polyribouridylic acid (polyU) support the theoretical predictions for these torsional angles. Using Monte Carlo technique, random coils of single-stranded polynucleotides have been simulated and the mean-square end-to-end distance 〈r²〉 has been calculated. Molecular orbital methods (EHT, CNDO, and PCILO) suggest considerable flexibility around O? P bonds, leading to fairly small values for the characteristic ratio (C_∞ ～ 4). Observed values of the unperturbed characteristic ratio for polynucleotides are quite large (C_∞ ～ 18) suggesting a relatively rigid nucleotide backbone. The results based on molecular orbital calculations can be reconciled with the experimental values by introducing an additional stabilization of ～2 kcal mol^?1 for the predicted minimum energy ragion (Φ′ ～ 240°, ω′ ～ 290°, ω 290°, Φ 180°, and ψ 60°). Such a stabilization may arise from the association of water molecules and metal ions with the phosphate group and (or) Coulomb interaction between neighboring phosphate groups. The calculations provide a semiquantitative estimate of torsional rigidity in the nucleotide backbone.     

Conformational dependence of <Superscript>13</Superscript>C shielding and coupling constants for methionine methyl groups

Methionine residues fulfill a broad range of roles in protein function related to conformational plasticity, ligand binding, and sensing/mediating the effects of oxidative stress. A high degree of internal mobility, intrinsic detection sensitivity of the methyl group, and low copy number have made methionine labeling a popular approach for NMR investigation of selectively labeled protein macromolecules. However, selective labeling approaches are subject to more limited information content. In order to optimize the information available from such studies, we have performed DFT calculations on model systems to evaluate the conformational dependence of ³ J _CSCC, ³ J _CSCH, and the isotropic shielding, σ_iso. Results have been compared with experimental data reported in the literature, as well as data obtained on [methyl-¹³C]methionine and on model compounds. These studies indicate that relative to oxygen, the presence of the sulfur atom in the coupling pathway results in a significantly smaller coupling constant, ³ J _CSCC/³ J _COCC ~ 0.7. It is further demonstrated that the ³ J _CSCH coupling constant depends primarily on the subtended CSCH dihedral angle, and secondarily on the CSCC dihedral angle. Comparison of theoretical shielding calculations with the experimental shift range of the methyl group for methionine residues in proteins supports the conclusion that the intra-residue conformationally-dependent shift perturbation is the dominant determinant of δ¹³Cε. Analysis of calmodulin data based on these calculations indicates that several residues adopt non-standard rotamers characterized by very large ~100° χ³ values. The utility of the δ¹³Cε as a basis for estimating the gauche/trans ratio for χ³ is evaluated, and physical and technical factors that limit the accuracy of both the NMR and crystallographic analyses are discussed.     

Chirality in amino acids beyond the Cα configuration

Based on a CSD search, a meta‐analysis of 1179 structures of 19 natural amino acids H₃NC_αH(R)C′(O)O and their derivatives H₃NC_αH(R)C′(O)O(H/R/M), protonated, esterified, or coordinated at the carboxylic group, shows that the chirality chain with its two steps, established in the preceding paper for alanine, can be extended to natural amino acids. High diastereoselectivities are observed in the induction from the L configuration at C_α to the ?ψ and +ψ conformations, which in turn distort the planar carboxylic group C_αC′(O_cis)O_trans to asymmetric flat tetrahedra, showing that the chirality chain is an integral part of natural amino acids.     

NMR study of the structure of short-chain peptides in solution.

The electron-mediated spin–spin coupling constant J between the amide NH and the α-CH protons in the dipeptide fragment C^α? CO(NH? C^αH)R? C′ONH? C^α is dependent on the dihedral angle of rotation (Φ) around the N? C bond. Measurement of J in a series of zwitterionic dipeptides H₃N⁺? CHR₁? CONH? CHR₂? CO₂^? (which is conformationally similar to the dipeptide fragment) in TFA solution shows that J is independent of R₁, but dependent on the steric bulk of R₂. The data are interpreted in terms of a model that assumes that what we measure is an average value of J? a thermal average over all the possible rotamers. The groups R₁ and R₂ are, in most cases, sterically kept apart by the trans and planar amide bonds, and hence the independence of J of R₁. This model is consistent with the theoretical calculations done on the dipeptide fragment. The effect of the structural characteristics of the side chains (e.g., the effect of lengthening and branching the side chains) on the J values in dipeptides is discussed in the light of the existing results of theoretical calculations. Study of 〈J〉 values in tripeptides (C₆H₅CH₂OCONH? CHR₁? CONH? CHR₂? CO₂CH₃, essentially three linked peptide units) shows that electrostatic interaction between the two amide bonds modifies the potential energy surface and the 〈J〉 value of a dipeptide subunit in the tripeptides. Also in some cases, direct steric interaction between the two side chains in the two adjacent dipeptide subunits in the tripeptide affects the potential energy surfaces of the individual dipeptide subunits and hence the 〈J〉 values. The influence of the structural characteristics of the side chains of individual amino acids on structure formation at or beyond the dipeptide level is discussed at various points. The J(NH? ^αCH) values of CH₃CONH? CHR? CONH₂ and CH₃CONH? CHR? CO₂CH₃ with the same R are quite different for R = valine, leucine, phenylalanine, methionine, but equal for R = glycine. This, coupled with the fact that one of the carboxamide NH resonances has a chemical shift different from its counterpart in simple amides like CH₃CONH₂ and the other carboxamide NH has the same chemical shift as its counterpart in CH₃CONH₂, suggest the presence of a hydrogen bond in dipeptide CH₃CONH? CHR? CONH₂ with carboxamide NH as the donor. Theoretical evidence for two seven-membered hydrogen-bonded rings with the carboxamide NH as donor and the acetyl oxygen as acceptor is summarized. Our data cannot suggest the number of such hydrogen-bonded rings, nor can they conclude the relative proportion of these rings in a particular dipeptide. A discussion of the difficulty of interpretation is presented and the data are discussed under certain simplifying assumptions.     

Intramolecular steric effects and hydrogen bonding in regular conformations of polyamino acids

A survey has been made, by using computer methods, of the types of helices which polypeptide chains can form, taking into account steric requirements and intramolecular hydrogen-bonding interactions. The influence on these two requirements, of small variations in the bond angles of the peptide residues, or of small changes in the overall dimensions of the helix (pitch and residues per turn), have been assessed for the special case of the α-helix. Criteria for the formation of acceptable hydrogen bonds have also been applied to helices of other types, viz., the 3, γ^?, ω^?, and π-helices. It was shown that the N? H … O and H … O? C angles in hydrogen bonds are sensitive to changes in either the NC^αC′ bond angle or in the rotational angles about the N? C^α and C^α? C′ bonds. However, the variants of the α-helix observed experimentally in myoglobin can all be constructed without distortion of the hydrogen bonds. For α-helices, the steric and hydrogen bonding requirements are more easily fulfilled with an NC^αC′ bond angle of 111°, rather than 109.5°. The decreased stability observed for the left-handed α-helix relative to the right-handed one for L -amino acids is due essentially only to interactions of the C^β atom of the side chains with atoms in adjacent peptide units in the backbone, and interactions with atoms in adjacent turns of the helical backbone are not significantly different in the two helices. Restrictions in the freedom of rotation of bulky side chains may have significant kinetic effects during the formation of the α-helix from the “random coil” state.     

Pressure alters electronic orbital overlap in hydrogen bonds

Pressure-induced changes in ^3h J _NC scalar couplings through hydrogen bonds were investigated in the immunoglobulin binding domain of streptococcal protein G. ¹H, ¹⁵N and ¹³C triple-resonance NMR spectroscopy coupled with the on-line high pressure cell technique was used to monitor ^3h J _NC scalar couplings at 30 and 2000 bar in uniformly labeled ¹⁵N and ¹³C protein isotopes. Both increased and decreased ^3h J _NC scalar couplings were observed at high pressure. No correlation with secondary structure was apparent. The difference in coupling constants as well as pressure-induced chemical shift data suggests a compaction of the helix ends and an increase of the helix pitch at its center in response to pressure. Our data provides the first direct evidence that the electronic orbital overlap in protein backbone hydrogen bonds is altered by pressure.     

HNCA-TOCSY-CANH experiments with alternate <Superscript>13</Superscript>C-<Superscript>12</Superscript>C labeling: a set of 3D experiment with unique supra-sequential information for mainchain resonance assignment

Described here is a set of three-dimensional (3D) NMR experiments that rely on CACA-TOCSY magnetization transfer via the weak ³ \textJ_{\textCa\textCa} ^{ 3} {\text{J}}_{{{\text{C}}\alpha {\text{C}}\alpha }} coupling. These pulse sequences, which resemble recently described ¹³C detected CACA-TOCSY (Takeuchi et al. 2010) experiments, are recorded in ¹H₂O, and use ¹H excitation and detection. These experiments require alternate ¹³C-¹²C labeling together with perdeuteration, which allows utilizing the small ³ \textJ_{\textCa\textCa} ^{ 3} {\text{J}}_{{{\text{C}}\alpha {\text{C}}\alpha }} scalar coupling that is otherwise masked by the stronger ¹J_CC couplings in uniformly ¹³C labeled samples. These new experiments provide a unique assignment ladder-mark that yields bidirectional supra-sequential information and can readily straddle proline residues. Unlike the conventional HNCA experiment, which contains only sequential information to the ^{1 3} \textC^a ^{ 1 3} {\text{C}}^{\alpha } of the preceding residue, the 3D hnCA-TOCSY-caNH experiment can yield sequential correlations to alpha carbons in positions i−1, i + 1 and i−2. Furthermore, the 3D hNca-TOCSY-caNH and Hnca-TOCSY-caNH experiments, which share the same magnetization pathway but use a different chemical shift encoding, directly couple the ¹⁵N-¹H spin pair of residue i to adjacent amide protons and nitrogens at positions i−2, i−1, i + 1 and i + 2, respectively. These new experimental features make protein backbone assignments more robust by reducing the degeneracy problem associated with the conventional 3D NMR experiments.     

An NMR investigation of putative interresidue H-bonding in methyl α-cellobioside in solution

Methyl α-cellobioside (methyl β-d-glucopyranosyl-(1→4)-α-d-glucopyranoside) was labeled with ¹³C at C4′ for use in NMR studies in DMSO-d₆ solvent to attempt the detection of a trans-H-bond J-coupling (^3hJ_CCOH) between C4′ and OH3. Analysis of the OH3 signal at 600 MHz revealed only the presence of two homonuclear J-couplings: ³J_H3,OH3 and a smaller, longer range J_HH. No evidence for ^3hJ_C4′,OH3 was found. The longer range J_HH was traced to ⁴J_H4,OH3 based on 2D ¹H–¹H COSY data and inspection of the H2 and H4 signal lineshapes. A limited set of DFT calculations was performed on a methyl cellobioside mimic to evaluate the structural dependencies of ⁴J_H2,O3H and ⁴J_H4,O3H on the H3–C3–O3–H torsion angle. Computed couplings range from about −0.7 to about +1.1 Hz, with maximal values observed when the C–H and O–H bonds are roughly diaxial.     

Conformational dynamics in mixed α/β-oligonucleotides containing polarity reversals: A molecular dynamics study using time-averaged restraints*

Nucleic acid duplexes featuring a single alpha-anomeric thymidine inserted into each DNA strand via 3′-3′ and 5′-5′ phosphodiester linkages exhibit local conformational dynamics that are not adequately depicted by conventional restrained molecular dynamics (rMD) methods. We have used molecular dynamics with time-averaged NMR restraints (MDtar) to explore its applicability to describing the conformational dynamics of two α-containing duplexes – d(GCGAAT-3′-3′-αT-5′-5′-CGC)₂ and d(ATGG-3′-3′-αT-5′-5′-GCTC)?r(gagcaccau). In contrast to rMD, enforcing NOE-based distance restraints over a period of time in MDtar rather than instantaneously results in better agreement with the experimental NOE and J-data. This conclusion is based on the dramatic decreases in average distance and coupling constant violations (Δd _av, J _rms, and ΔJ _av) and improvements in sixth-root R-factors (R ^x). In both duplexes, the deoxyribose ring puckering behavior predicted independently by pseudorotation analysis is portrayed remarkably well using this approach compared to rMD. This indicates that the local dynamic behavior is encoded within the NOE data, although this is not obvious from the local R ^x values. In both systems, the backbone torsion angles comprising the 3′-3′ linkage as well as the (high S-) sugars of the α-nucleotide and preceding residue (α?1) are relatively static, while the conformations of the 5′-5′ linkage and the sugar in the neighboring β-nucleotide (α+1) show enhanced flexibility. To reduce the large ensembles generated by MDtar to more manageable clusters we utilized the PDQPRO program. The resulting PDQPRO clusters (in both cases, 13 structures and associated probabilities extracted from a pool of 300 structures) adequately represent the structural and dynamic characteristics predicted by the experimental data.