Crystallographic dissection of the thermal motion of protein-sugar complex

The crystal structure of turkey egg lysozyme (TEL) complexed with di-N-acetylchitobiose (NAG2) was refined at 1.19 A resolution by the full-matrix least-squares method with anisotropic temperature factors, and its thermal motion was evaluated by the TLS method. The average ESDs of atomic parameters of nonhydrogen atoms were 0.030 A for coordinates and 0.025 A(2) for anisotropic temperature factors. The active site cleft of TEL binds the alpha-anomer of NAG2 in a nonproductive binding mode with its pyranose rings parallel to a beta-sheet. The TEL structure was compared with the re-refined 1.12 A structure of native TEL. The RMS difference for equivalent Calpha atoms was 0.103 A and a relatively large difference was observed in the region of residues 104-125 rather than in the beta-sheet region where NAG2 was bound. In contrast, the temperature factor of the beta-sheet region was significantly decreased by the NAG2 binding. The TLS model that describes the rigid body motion in translation, libration, and screw motion was adopted for the evaluation of the molecular motion of TEL and NAG2, and the TLS parameters were determined by the least-squares fit to U(ij). The contribution of the external motion of TEL was estimated to be 55.8% of the observed temperature factor for the native structure and 45.9% for the NAG2 complex. The internal motion of TEL represented with atomic thermal ellipsoids was very similar between the native and complex structures except the NAG2 binding region. In the structure of NAG2, the rigid body motion dominates the thermal motion. The center of rotation of NAG2, 4.45A far from the center of gravity, is on the nitrogen atom of the acetylamino group that is hydrogen bonded to the main-chain peptide groups of Asn49 and Ala107. The rigid body motion of NAG2 indicates that the acetylamino group is most strongly bound to the active site, and the recognition of this group is a crucial step of the substrate binding.     

Full-matrix least-squares refinement of lysozymes and analysis of anisotropic thermal motion

Crystal structures of turkey egg lysozyme (TEL) and human lysozyme (HL) were refined by full-matrix least-squares method using anisotropic temperature factors. The refinement converged at the conventional R-values of 0.104 (TEL) and 0.115 (HL) for reflections with F_o > 0 to the resolution of 1.12 Å and 1.15 Å, respectively. The estimated r.m.s. coordinate errors for protein atoms were 0.031 Å (TEL) and 0.034 Å (HL). The introduction of anisotropic temperature factors markedly reduced the R-value but did not significantly affect the main chain coordinates. The degree of anisotropy of atomic thermal motion has strong positive correlation with the square of distance from the molecular centroid. The ratio of the radial component of thermal ellipsoid to the r.m.s. magnitude of three principal components has negative correlation with the distance from the molecular centroid, suggesting the domination of libration rather than breathing motion. The TLS model was applied to elucidate the characteristics of the rigid-body motion. The TLS tensors were determined by the least-squares fit to observed temperature factors. The profile of the magnitude of reproduced temperature factors by the TLS method well fitted to that of observed B_eqv. However, considerable disagreement was observed in the shape and orientation of thermal ellipsoid for atoms with large temperature factors, indicating the large contribution of local motion. The upper estimate of the external motion, 67% (TEL) and 61% (HL) of B_eqv, was deduced from the plot of the magnitude of TLS tensors determined for main chain atoms which were grouped into shells according to the distance from the center of libration. In the external motion, the translational portion is predominant and the contribution of libration and screw motion is relatively small. The internal motion, estimated by subtracting the upper estimate of the external motion from the observed temperature factor, is very similar between TEL and HL in spite of the difference in 54 of 130 amino acid residues and in crystal packing, being suggested to reflect the intrinsic internal motion of chicken-type lysozymes. Proteins 30:232–243, 1998. © 1998 Wiley-Liss, Inc.     

Atomic resolution structure of prokaryotic phospholipase A2: analysis of internal motion and implication for a catalytic mechanism

We have found a secreted phospholipase A(2) (PLA(2), EC 3.1.1.4) from Streptomyces violaceoruber A-2688, which is the first PLA(2) identified in prokaryote, and determined its tertiary structure by NMR and X-ray analyses. In this study, we collected the X-ray diffraction data of the bacterial PLA(2) at room temperature (297 K) using conventional MoK(alpha) radiation and refined the structure at a 1.05 A resolution. The atomic resolution analysis led us to introduce disordered conformations and hydrogen atoms into a full anisotropic model. The molecular motion, which is expressed as the sum of rigid-body motion and internal motion of protein, is roughly estimated as the thermal motion when the X-ray diffraction data are collected at room temperature. In this study, we applied a TLS (rigid-body motion in terms of translation, libration, and screw motions) model to analyze the rigid-body motion of the bacterial PLA(2) and calculated the internal motion by subtracting the estimate of the rigid-body motion from the observed anisotropic temperature factor. We also subjected the TLS model to estimate the internal motion of the bovine pancreatic PLA(2) using the anisotropic temperature factor deposited in the Protein Data Bank. Both results indicate that the localization of regions exhibiting larger internal motion in the bacterial PLA(2) is almost the same as that in the bovine pancreatic PLA(2), suggesting that although the tertiary structure of the bacterial PLA(2) is strikingly different from that of the bovine pancreatic PLA(2), the internal motion, which is associated with the calcium(II) ion-binding, phospholipid-binding, and allosteric interfacial activation, is commonly observed in both PLA(2)s.     

Critical evaluation of simple network models of protein dynamics and their comparison with crystallographic B-factors

In recent years, elastic network models (ENM) have been widely used to describe low-frequency collective motions in proteins. These models are often validated and calibrated by fitting mean-square atomic displacements estimated from x-ray crystallography (B-factors). We show that a proper calibration procedure must account for the rigid-body motion and constraints imposed by the crystalline environment on the protein. These fundamental aspects of protein dynamics in crystals are often ignored in currently used ENMs, leading to potentially erroneous network parameters. Here we develop an ENM that properly takes the rigid-body motion and crystalline constraints into account. Its application to the crystallographic B-factors reveals that they are dominated by rigid-body motion and thus are poorly suited for the calibration of models for internal protein dynamics. Furthermore, the translation libration screw (TLS) model that treats proteins as rigid bodies is considerably more successful in interpreting the experimental B-factors than ENMs. This conclusion is reached on the basis of a comparative study of various models of protein dynamics. To evaluate their performance, we used a data set of 330 protein structures that combined the sets previously used in the literature to test and validate different models. We further propose an extended TLS model that treats the bulk of the protein as a rigid body while allowing for flexibility of chain ends. This model outperforms other simple models of protein dynamics in interpreting the crystallographic B-factors.     

Crystallographic analysis of the thermal motion of the inclusion complex of cyclomaltoheptaose (beta-cyclodextrin) with hexamethylenetetramine

The crystal structure of the inclusion complex of cyclomaltoheptaose (beta-cyclodextrin) with hexamethylenetetramine was determined at temperatures of 123, 173, 223, and 293 K. The rigid-body motion of the host and guest molecules was evaluated by means of the TLS method that represents the molecular motion in terms of translation, libration, and screw motion. In increasing the temperature from 123 to 293 K, the amplitude of the rigid body vibration of the host molecule was increased from 1.0 to 1.3 degrees in the rotational motion and from 0.16 to 0.17 A in the translational motion. The cyclomaltoheptaose molecule has the flexibility in seven alpha-(1-->4)-linkages, and each glucose unit was in the rotational vibration around an axis through two glycosidic oxygen atoms. As a result, the rigid-body parameters of cyclomaltoheptaose were considered to be overestimated because of including the contribution from the local motion of glucose units. In contrast, for the guest molecule having no structural flexibility, the TLS analysis demonstrated that the atomic thermal vibration was mostly derived from the rigid body motion. The rotational amplitude of hexamethylenetetramine was changed from 5.2 to 6.6 degrees in increasing the temperature from 123 to 293 K, while the change of the translational amplitude was from 0.20 to 0.23 A. The translational motion of the guest molecule was hindered by the inside wall of the host cavity. The molecular motion was characterized by the rotational vibration around the axis through two nitrogen atoms that were involved in the hydrogen-bond formation.     

Picosecond timescale rigid-helix and side-chain motions in deoxymyoglobin

The contribution of rigidbody motions to the atomic trajectories in a 100 ps molecular dynamics simulation of deoxymyoglobin is examined. Two typesof rigid-body motions are considered: one in which the helices are rigid units and one in which the side-chains are rigid units. Using a quaternionbased algorithm, fits of the rigid reference structures are made to each time frame of the simulation to derive trajectories of the rigid-body motions. The fitted trajectories are analysed in terms of atomic position fluctuations, mean-square displacements as a function of time, velocity autocorrelation functions and densities of states. The results are compared with the corresponding quantities calculated from the full trajectory. The relative contribution of the rigid helix motions to the helix atom dynamics depends on which quantity is examined and on which subset of atoms is chosen: rigid-helix motions contribute 86% of the rms helix backbone atomic position fluctuations, but 30% of the helix,: atom (backbone and side-chain) mean square displacements and only 1.1% of total kinetic energy. Only very low-frequency motions contribute to the rigid-helix dynamics; the rigid-body analysis allows characteristic rigid-helix vibrations to be identified and described. Treating the side-chains as rigid bodies is foundto be an excellent approximation to both their diffusive and vibrationalmean-square displacements: 96% of side-chain atom mean-square displacements originate from rigid side-Chain motions. However, the errors in theside-chain atomic positional fits are not always small. An analysis is madeof factors contributing to the positional error for different types of side-chain. © Wiley-Liss, Inc.     

Projection Methods for the Analysis of Complex Motions in Macromolecules

Abstract

In studies of macromolecular dynamics it is often desirable to analyze complex motions in terms of a small number of coordinates. Only for simple types of motion, e.g., rigid-body motions, these coordinates can be easily constructed from the Cartesian atomic coordinates. This article presents an approach that is applicable to infinitesimal or approximately infinitesimal motions, e.g., Cartesian velocities, normal modes, or atomic fluctuations. The basic idea is to characterize the subspace of interesting motions by a set of (possibly linearly dependent) vectors describing elementary displacements, and then project the dynamics onto this subspace. Often the elementary displacements can be found by physical intuition. The restriction to small displacements facilitates the study of complicated coupled motions and permits the construction of collective-motion subspaces that do not correspond to any set of generalized coordinates.

As an example for this technique, we analyze the low-frequency normal modes of proteins up to ≈ 20 THz (600 cm^?1) in order to see what kinds of motions occupy which frequency range. This kind of analysis is useful for the interpretation of spectroscopic measurements on proteins, e.g., inelastic neutron scattering experiments.     

Refinement of the influenza virus hemagglutinin by simulated annealing

We have applied the method of simulated annealing to the refinement of the 3 A resolution crystal structure of the influenza virus hemagglutinin glycoprotein, using the program X-PLOR. Two different methods were introduced into X-PLOR to treat the non-crystallographic symmetry present in this and in other crystal structures. In the first, only the unique protomer atoms are refined; by application of the non-crystallographic symmetry operators to the protomer atoms, the X-ray structure factor derivatives are effectively averaged, and a non-bonded energy term models the interactions of the protomer with its neighbors in the oligomer without explicit refinement of the other protomers in the crystallographic asymmetric unit. In the second method, the entire asymmetric unit is refined, but an effective energy term is added to the empirical energy that restrains symmetry-related atomic positions to their average values after least-squares superposition. Several other modifications and additions were made to previously published X-PLOR protocols, including weighting of the X-ray terms, maintenance of the temperature of the molecular dynamics simulation, treatment of charged groups, changes in the values of certain empirical energy parameters, and the use of N-linked carbohydrate empirical energy parameters. The hemagglutinin refinement proceeded in several stages. An initial round of simulated annealing of the monomer was followed by rigid-body refinement of the 3-fold non-crystallographic symmetry axis position and a second round of monomer refinement. A third round was performed on the trimer using non-crystallographic symmetry restraints in all regions except those in lattice contacts showing obvious derivations from 3-fold symmetry. The refinement was completed with several rounds of conventional positional and isotropic temperature factor refinement needed to correct bad model geometry introduced by high-temperature molecular dynamics in regions of weak electron density. This structure was then used as the basis for refinement of three crystallographically isomorphous hemagglutinin structures, including complexes with the influenza virus receptor, sialic acid. Model geometry comparable to well-refined high-resolution structures was obtained with relatively little manual intervention, demonstrating the ability of simulated annealing refinement to produce highly idealized structures at moderate resolution.     

Effect of anisotropy and anharmonicity on protein crystallographic refinement. An evaluation by molecular dynamics

Molecular dynamics simulations are employed to determine the errors introduced by anharmonicity and anisotropy in the structure and temperature factors obtained for proteins by refinement of X-ray diffraction data. Simulations (25 ps and 300 ps) of metmyoglobin are used to generate time-averaged diffraction data at 1.5 A resolution. The crystallographic restrained-parameter least-squares refinement program PROLSQ is used to refine models against these simulated data. The resulting atomic positions and isotropic temperature factors are compared with the average structure and fluctuations calculated directly from the simulations. It is found that significant errors in the atomic positions and fluctuations are introduced by the refinement, and that the errors increase with the magnitude of the atomic fluctuations. Of particular interest is the fact that the refinement generally underestimates the atomic motions. Moreover, while the actual fluctuations go up to a mean-square value of about 5 A2, the X-ray results never go above approximately 2 A2. This systematic deviation in the motional parameters appears to be due to the use of a single-site isotropic model for the atomic fluctuations. Many atoms have multiple peaks in their probability distribution functions. For some atoms, the multiple peaks are seen in difference electron density maps and it is possible to include these in the refinement as disordered residues. However, for most atoms the refinement fits only one peak and neglects the rest, leading to the observed errors in position and temperature factor. The use of strict stereochemical restraints is inconsistent with the average dynamical structure; nevertheless, refinement with tight restraints results in structures that are comparable to those obtained with loose restraints and better than those obtained with no restraints. The results support the use of tight stereochemical restraints, but indicate that restraints on the variation of temperature factors are too restrictive.     

An atomic resolution structure for human fibroblast growth factor 1

A 1.10-A atomic resolution X-ray structure of human fibroblast growth factor 1 (FGF-1), a member of the beta-trefoil superfold, has been determined. The beta-trefoil is one of 10 fundamental protein superfolds and is the only superfold to exhibit 3-fold structural symmetry (comprising 3 "trefoil" units). The quality of the diffraction data permits unambiguous assignment of Asn, Gln, and His rotamers, Pro ring pucker, as well as refinement of atomic anisotropic displacement parameters (ADPs). The FGF-1 structure exhibits numerous core-packing defects, detectable using a 1.0-A probe radius. In addition to contributing to the relatively low thermal stability of FGF-1, these defects may also permit domain motions within the structure. The availability of refined ADPs allows a translation/libration/screw (TLS) analysis of putative rigid body domains. The TLS analysis shows that beta-strands 6-12 together form a rigid body, and there is a clear demarcation in TLS motions between the adjacent carboxyl- and amino-termini. Although separate from beta-strands 6-12, the individual beta-strands 1-5 do not exhibit correlated motions; thus, this region appears to be comparatively flexible. The heparin-binding contacts of FGF-1 are located within beta-strands 6-12; conversely, a significant portion of the receptor-binding contacts are located within beta-strands 1-5. Thus, the observed rigid body motion in FGF-1 appears related to the ligand-binding functionalities.     

Dynamic character of the complex of human blood coagulation factor VIIa with the extracellular domain of human tissue factor: a normal mode analysis

As an attempt to investigate the dynamic interactions between plasma serine protease, coagulation factor VIIa (VIIa) and its cofactor, tissue factor (TF), we performed normal mode analysis (NMA) of the complex of VIIa with soluble TF (the extracellular part of TF; sTF). We compared fluctuations of Calpha atoms of VIIa or sTF derived from NMA in the VIIa-sTF complex with those of VIIa or sTF in an uncomplexed condition. The atomic fluctuations of the Calpha atoms of sTF complexed with VIIa did not significantly differ from those of sTF without VIIa. In contrast, the atomic fluctuations of VIIa complexed with sTF were much smaller than those of VIIa without sTF. These results suggest that domain motions of VIIa molecule alone are markedly dampened in the VIIa-sTF complex and that the sTF molecule is relatively more rigid than the VIIa molecule. This may indicate functions of TF as a cofactor.     

Structure and internal mobility of proteins: a molecular dynamics study of hen egg white lysozyme

The structure and internal motions of the protein hen egg white lysozyme are studied by analysis of simulation and experimental data. A molecular dynamics simulation and an energy minimization of the protein in vacuum have been made and the results compared with high-resolution structures and temperature factors of hen egg white lysozyme in two different crystal forms and of the homologous protein human lysozyme. The structures obtained from molecular dynamics and energy minimization have root-mean-square deviations for backbone atoms of 2.3 Å and 1.1–1.3 Å, respectively, relative to the crystal structures; the different crystal structures have root-mean-square deviations of 0.73–0.81 Å for the backbone atoms. In comparing the backbone dihedral angles, the difference between the dynamics and the crystal structure on which it is based is the same as that between any two crystal structures. The internal fluctuations of atomic positions calculated from the molecular dynamics trajectory agree well with the temperature factors from the three structures. Simulation and crystal results both show that there are large motions for residues involved in exposed turns of the backbone chain, relatively smaller motions for residues involved in the middle of helices or β-sheet structures, and relatively small motions of residues near disulfide bridges. Also, both the simulation and crystal data show that side-chain atoms have larger fluctuations than main-chain atoms. Moreover, the regions that have large deviations among the x-ray crystal structures, which indicates flexibility, are found to have large fluctuations in the simulation.     

Exploring the structural dynamics of the E.coli chaperonin GroEL using translation-libration-screw crystallographic refinement of intermediate states

Large rigid-body domain movements are critical to GroEL-mediated protein folding, especially apical domain elevation and twist associated with the formation of a folding chamber upon binding ATP and co-chaperonin GroES. Here, we have modeled the anisotropic displacements of GroEL domains from various crystallized states, unliganded GroEL, ATPgammaS-bound, ADP-AlFx/GroES-bound, and ADP/GroES bound, using translation-libration-screw (TLS) analysis. Remarkably, the TLS results show that the inherent motions of unliganded GroEL, a polypeptide-accepting state, are biased along the transition pathway that leads to the folding-active state. In the ADP-AlFx/GroES-bound folding-active state the dynamic modes of the apical domains become reoriented and coupled to the motions of bound GroES. The ADP/GroES complex exhibits these same motions, but they are increased in magnitude, potentially reflecting the decreased stability of the complex after nucleotide hydrolysis. Our results have allowed the visualization of the anisotropic molecular motions that link the static conformations previously observed by X-ray crystallography. Application of the same analyses to other macromolecules where rigid body motions occur may give insight into the large scale dynamics critical for function and thus has the potential to extend our fundamental understanding of molecular machines.     

The effects of ligand binding on the backbone dynamics of the kringle 1 domain of human plasminogen

The internal motions of the backbone nitrogen atoms of the kringle 1 domain of human plasminogen (K1(Pg)) were examined in the absence and presence of the ligand, epsilon-aminocaproic acid. These dynamic properties were determined from (15)N NMR relaxation data in terms of the extended model-free parameters. The model of isotropic reorientation was found sufficient to account for overall molecular tumbling for both apo and EACA-bound K1(Pg). The global rotational correlation time (tau(m)) for apo-K1(Pg) was 5.87(+/-0.01) ns, while the tau(m) for ligand-bound K1(Pg) was 5.20(+/-0.01) ns, suggesting that perhaps some small degree of aggregation occurred in the apo form of the kringle module. Complexation of K1(Pg) with ligand mainly reduced those internal motions that occurred on a 100 ps to 5 ns time-scale. The magnitude of the chemical exchange was also attenuated upon ligand binding. These data are consistent with studies employing other approaches that suggest that the binding pocket is preformed in K1(Pg).     

Crystal structure of human serum albumin at 2.5 A resolution.

A new triclinic crystal form of human serum albumin (HSA), derived either from pool plasma (pHSA) or from a Pichia pastoris expression system (rHSA), was obtained from polyethylene glycol 4000 solution. Three-dimensional structures of pHSA and rHSA were determined at 2.5 A resolution from the new triclinic crystal form by molecular replacement, using atomic coordinates derived from a multiple isomorphous replacement work with a known tetragonal crystal form. The structures of pHSA and rHSA are virtually identical, with an r.m. s. deviation of 0.24 A for all Calpha atoms. The two HSA molecules involved in the asymmetric unit are related by a strict local twofold symmetry such that the Calpha atoms of the two molecules can be superimposed with an r.m.s. deviation of 0.28 A in pHSA. Cys34 is the only cysteine with a free sulfhydryl group which does not participate in a disulfide linkage with any external ligand. Domains II and III both have a pocket formed mostly of hydrophobic and positively charged residues and in which a very wide range of compounds may be accommodated. Three tentative binding sites for long-chain fatty acids, each with different surroundings, are located at the surface of each domain.     

Nucleic-acid structural deformability deduced from anisotropic displacement parameters

The growing numbers of very well resolved nucleic-acid crystal structures with anisotropic displacement parameters provide an unprecedented opportunity to learn about the natural motions of DNA and RNA. Here we report a new Monte-Carlo approach that takes direct account of this information to extract the distortions of covalent structure, base pairing, and dinucleotide geometry intrinsic to regularly organized double-helical molecules. We present new methods to test the validity of the anisotropic parameters and examine the apparent deformability of a variety of structures, including several A, B, and Z DNA duplexes, an AB helical intermediate, an RNA, a ligand-DNA complex, and an enzyme-bound DNA. The rigid-body parameters characterizing the positions of the bases in the structures mirror the mean parameters found when atomic motion is taken into account. The base-pair fluctuations intrinsic to a single structure, however, differ from those extracted from collections of nucleic-acid structures, although selected base-pair steps undergo conformational excursions along routes suggested by the ensembles. The computations reveal surprising new molecular insights, such as the stiffening of DNA and concomitant separation of motions of contacted nucleotides on opposite strands by the binding of Escherichia coli endonuclease VIII, which suggest how the protein may direct enzymatic action.     

Domain motions of glucosamine-6P synthase: comparison of the anisotropic displacements in the crystals and the catalytic hinge-bending rotation

Glucosamine-6-phosphate synthase channels ammonia over 18 A from glutamine at the glutaminase site to fructose-6P at the synthase site. We have modeled the anisotropic displacements of the glutaminase and synthase domains from the two crystallized states, the enzyme in complex with fructose-6P or in complex with glucose-6P and a glutamine affinity analog, using TLS (rigid-body motion in terms of translation, libration, and screw motions) refinement implemented in REFMAC. The domains displacements in the crystal lattices are compared to the movement of the glutaminase domain relative to the synthase domain that occurs during the catalytic cycle upon glutamine binding, which was visualized by comparing the two structures. This movement was analyzed by the program DYNDOM as a 22.8 degrees rotation around an effective hinge axis running approximately parallel to helix 300-317 of the synthase domain, the glutaminase loop that covers the glutaminase site upon glutamine binding acting as the mechanical hinge.     

Dynamic transition associated with the thermal denaturation of a small Beta protein

We studied the temperature dependence of the picosecond internal dynamics of an all-beta protein, neocarzinostatin, by incoherent quasielastic neutron scattering. Measurements were made between 20 degrees C and 71 degrees C in heavy water solution. At 20 degrees C, only 33% of the nonexchanged hydrogen atoms show detectable dynamics, a number very close to the fraction of protons involved in the side chains of random coil structures, therefore suggesting a rigid structure in which the only detectable diffusive movements are those involving the side chains of random coil structures. At 61.8 degrees C, although the protein structure is still native, slight dynamic changes are detected that could reflect enhanced backbone and beta-sheet side-chain motions at this higher temperature. Conversely, all internal dynamics parameters (amplitude of diffusive motions, fraction of immobile scatterers, mean-squared vibration amplitude) rapidly change during heat-induced unfolding, indicating a major loss of rigidity of the beta-sandwich structure. The number of protons with diffusive motion increases markedly, whereas the volume occupied by the diffusive motion of protons is reduced. At the half-transition temperature (T = 71 degrees C) most of backbone and beta-sheet side-chain hydrogen atoms are involved in picosecond dynamics.     

Relationship of DNA structure to internal dynamics: correlation of helical parameters from NOE-based NMR solution structures of d(GCGTACGC)(2) and d(CGCTAGCG)(2) with (13)C order parameters implies conformational coupling in dinucleotide units

The coupling between the conformational properties of double-stranded DNA and its internal dynamics has been examined. The solution structures of the isomeric DNA oligomers d(GCGTACGC)(2) (UM) and d(CGCTAGCG)(2) (CTSYM) were determined with (1)H NMR spectroscopy by utilizing distance restraints from total relaxation matrix analysis of NOESY cross-peak intensities in restrained molecular dynamics calculations. The root-mean-square deviation of the coordinates for the ensemble of structures was 0.13 A for UM and 0.49 A for CTSYM, with crystallographic equivalent R(c)=0.41 and 0.39 and sixth-root residual R(x)=0.11 and 0.10 for UM and CTSYM, respectively. Both UM and CTSYM are B-form with straight helical axes and show sequence-dependent variations in conformation. The internal dynamics of UM and CTSYM were previously determined by analysis of (13)C relaxation parameters in the context of the Lipari & Szabo model-free formalism. Helical parameters for the two DNA oligomers were examined for linear correlations with the order parameters (S(2)) of groups of (13)C spins in base-pairs and dinucleotide units of UM and CTSYM. Correlations were found for six interstrand base-pair parameters tip, y-displacement, inclination, buckle and stretch with various combinations of S(2) for atoms in Watson-Crick base-pairs and for two inter-base-pair parameters, rise and roll with various combinations of S(2) for atoms in dinucleotides. The correlations for the interstrand base-pair helical parameters indicate that the conformations of the deoxyribose residues of each strand are dynamically coupled. Also, the inter-base-pair separation has a profound effect on the local internal motions available to the DNA, supporting the idea that rise is a principal degree of freedom for DNA conformational variability. The correlations indicate collective atomic motions of spins that may represent specific motional modes in DNA, and that base sequence has a predictable effect on the relative order of groups of spins both in the bases and in the deoxyribose ring of the DNA backbone. These observations suggest that an important functional outcome of DNA base sequence is the modulation of both the conformation and dynamic behavior of the DNA backbone.