Hydrophobicity profiles for protein sequence analysis.

Hydrophobic interactions are a major force in protein folding and numerous hydropathy scales have been developed to quantify the relative hydrophobicity of the amino acids. Hydropathy profiles can be used to examine the surface features of proteins in order to generate hypotheses that can be confirmed experimentally. This unit describes the application of hydrophobicity plots to typical problems and provides suggested uses for a few selected scales.

Protein secondary structure prediction.

This unit describes procedures developed for predicting protein structure from the amino acid sequence. The first of the four sections is an overview and brief history of structure prediction schemes. The second section describes four distinct prediction schemes, with emphasis on their differences. In the third part each prediction scheme is used to evaluate three proteins that have different folding patterns. The final section is a comparison of the prediction results and suggestions for secondary structure prediction.

Analysis of the three-dimensional antigenic structure of giant ragweed allergen, Amb t 5.

The ragweed allergens Amb t 5 and Amb a 5 are among the smallest inhaled protein allergens known, containing a single, immunodominant T-cell epitope. In this study we analyzed the B-cell epitope structure of Amb t 5. The three-dimensional structures of Amb t 5 and Amb a 5 have been determined by NMR spectroscopy, providing a rare opportunity to analyze three-dimensional antigenic sites. Amb t 5 residues likely to be important for antigenicity were identified by examining the surface area of Amb t 5 accessible to a probe of the size of an antibody molecule. After changing these residues to the corresponding Amb a 5 residues, recombinant proteins were purified and tested for loss of antigenic activity. Inhibition radio-immunoassays, using sera from 8 individuals who had received immunotherapy with giant ragweed extract, allowed the mutations to be divided into three groups: (1) mutations that had little or no effect on antibody binding, (2) mutations that caused a loss of antigenic activity to a different degree in different sera and (3) mutations that drastically reduced antigenic activity in all sera tested. This last set of mutations clustered in the third loop of Amb t 5, suggesting that antibody recognition of Amb t 5, like T-cell recognition, is primarily directed towards a single, immunodominant site.

Solution structure of human CTLA-4 and delineation of a CD80/CD86 binding site conserved in CD28.

The structure of human CTLA-4 reveals that residues Met 99, Tyr 100 and Tyr 104 of the M99YPPPY104 motif are adjacent to a patch of charged surface residues on the A'GFCC' face of the protein. Mutation of these residues, which are conserved in the CTLA-4/CD28 family, significantly reduces binding to CD80 and/or CD86, implicating this patch as a ligand binding site.

Molecular modeling of CD28 and three-dimensional analysis of residue conservation in the CD28/CD152 family.

CD28/CD152-CD80/CD86 receptor-ligand interactions result in costimulatory signals critical for optimal T cell activation. CD28/CD152 and CD80/CD86 are members of the immunoglobulin superfamily (IgSF). Despite common receptor-ligand interactions, both receptor and ligand pairs share only limited sequence identity. A detailed molecular model of the extracellular Ig-like domain of human CD28 was constructed using a combination of different modeling methods. The model was based on the solution structure of CD152 and sequence comparison of the CD28/CD152 family. Assessment of the model revealed good stereochemical quality and sequence-structure compatibility. The CD28 model was used to map surface residues, N-linked glycosylation sites, and to compare residue conservation in CD28 and CD152. The location of N-linked glycosylation sites in CD28/CD152 restricts the surface area available for binding. Rigorous sequence conservation in CD28 and CD152 is limited to core IgSF consensus positions and surface residues implicated in ligand binding. Other surface residues vary greatly in CD28/CD152. Residues critical for ligand binding are surrounded by surface patches conserved only in either CD28 or CD152.

Characterization of NADP+ binding to perdeuterated MurB: backbone atom NMR assignments and chemical-shift changes.

Backbone-atom resonances have been assigned for both the substrate-free and the NADP+-complexed forms of UDP-N-acetylenolpyruvylglucosamine reductase (MurB), a monomeric, 347-residue (38.5 kDa) flavoenzyme essential for bacterial cell-wall biosynthesis. NMR studies were performed using perdeuterated, uniformly 13C/15N-labeled samples of MurB. In the case of substrate-free MurB, one or more backbone atoms have been assigned for 334 residues (96%). The assigned backbone atoms include 309 1HN and 15N atoms (94%), 315 13CO atoms (91%), 331 13C(alpha) atoms (95%), and 297 13C(beta) atoms (93%). For NADP+-complexed MurB, one or more backbone atoms have been assigned for 313 residues (90%); these include 283 1HN and 15N atoms (86%), 305 13CO atoms (88%), 310 13C(alpha) atoms (89%), and 269 13C(beta) atoms (84%). The strategies used for obtaining resonance assignments are described in detail. Information on the secondary structure in solution for both the substrate-free and NADP+-complexed forms of the enzyme has been derived both from 13C(alpha) and 13C(beta) chemical-shift deviations from random-coil values and from 1HN-1HN NOEs. These data are compared to X-ray crystallographic structures of substrate-free MurB and MurB complexed with the UDP-N-acetylglucosamine enolpyruvate (UNAGEP) substrate. NADP+ binding induces significant chemical-shift changes in residues both within the known UNAGEP and FAD binding pockets and within regions known to undergo conformational changes upon UNAGEP binding. The NMR data indicate that NADP+ and UNAGEP utilize the same binding pocket and, furthermore, that the binding of NADP+ induces structural changes in MurB. Finally, many of the residues within the UNAGEP/NADP+ binding pocket were difficult to assign due to dynamic processes which weaken and/or broaden the respective resonances. Overall, our results are consistent with MurB having a flexible active site.

The three-dimensional solution structure of the SH2 domain from p55blk kinase.

Signal transduction in B cells is mediated, in part, by the interaction of the cytoplasmic components of the antigen receptor complex and various members of the src family tyrosine kinases. Key to this process appears to be the interaction of the tyrosine kinase SH2 domains with the tyrosine-phosphorylated cytoplasmic domain of Ig-alpha, a disulfide-bonded heterodimeric (with Ig-beta or Ig-gamma) transmembrane protein that noncovalently associates with the antigen receptor immunoglobin chains. In addition to binding to the phosphorylated cytoplasmic domains of Ig-alpha and Ig-beta, blk and fyn(T), two members of the src family kinases, have been shown to bind overlapping but distinct sets of phosphoproteins [Malek & Desiderio (1993) J. Biol. Chem. 268. 22557-22565]. A comparison of their three-dimensional structures may elucidate the apparently subtle differences required for phosphoprotein discrimination. To begin characterizing the blk/fyn/phosphosphoprotein interactions, we have determined the three-dimensional solution structure of the SH2 domain of blk kinase by nuclear magnetic resonance (NMR) spectroscopy. 1H, 13C, and 15N resonances of the SH2 domain of blk kinase were assigned by analysis of multidimensional, double- and triple-resonance NMR experiments. Twenty structures of the blk SH2 domain were refined with the program X-PLOR using a total of 2080 experimentally derived conformational restraints. The structures converged to a root-mean-squared (rms) distance deviation of 0.51 and 0.95 A for the backbone atoms and for the non-hydrogen atoms, respectively. The blk SH2 domain adopts the prototypical SH2 fold. Structurally, blk SH2 is most similar to the crystal structure of the v-src SH2 domain [Waksman et al. (1993) Nature 358.646-653] and superimposes on the crystal structure with an rmsd of 1.52 A for the backbone atoms. The largest deviations occur in the four loops interconnecting beta-strands A-E, which are the least well-defined regions in the NMR structure. Exclusion of these loops lowers this rmsd to 0.82 A. The conformation of the BC loop in the blk SH2 domain is similar to the open conformation in the apo lck SH2 domain, suggesting that, like the lck SH2 domain, the blk SH2 domain may have a gated phosphopeptide binding site. Finally, it is proposed that the amino acid substitution of Lys 88 (blk) for Glu [fyn(T)] is important for the observed differences in specificity between blk and fyn(T) SH2 domains.

Simulated annealing with restrained molecular dynamics using a flexible restraint potential: theory and evaluation with simulated NMR constraints.

A new functional representation of NMR-derived distance constraints, the flexible restraint potential, has been implemented in the program CONGEN (Bruccoleri RE, Karplus M, 1987, Biopolymers 26:137-168) for molecular structure generation. In addition, flat-bottomed restraint potentials for representing dihedral angle and vicinal scalar coupling constraints have been introduced into CONGEN. An effective simulated annealing (SA) protocol that combines both weight annealing and temperature annealing is described. Calculations have been performed using ideal simulated NMR constraints, in order to evaluate the use of restrained molecular dynamics (MD) with these target functions as implemented in CONGEN. In this benchmark study, internuclear distance, dihedral angle, and vicinal coupling constant constraints were calculated from the energy-minimized X-ray crystal structure of the 46-amino acid polypeptide crambin (ICRN). Three-dimensional structures of crambin that satisfy these simulated NMR constraints were generated using restrained MD and SA. Polypeptide structures with extended backbone and side-chain conformations were used as starting conformations. Dynamical annealing calculations using extended starting conformations and assignments of initial velocities taken randomly from a Maxwellian distribution were found to adequately sample the conformational space consistent with the constraints. These calculations also show that loosened internuclear constraints can allow molecules to overcome local minima in the search for a global minimum with respect to both the NMR-derived constraints and conformational energy. This protocol and the modified version of the CONGEN program described here are shown to be reliable and robust, and are applicable generally for protein structure determination by dynamical simulated annealing using NMR data.

The ligand binding domain of the human retinoic acid receptor gamma is predominantly alpha-helical with a Trp residue in the ligand binding site.

Retinoic acid exerts its many biological effects by interaction with a nuclear protein, the retinoic acid receptor (RAR). The details of this interaction are unknown due mainly to the lack of sufficient quantities of pure functional receptor protein for biochemical and structural studies. We have recently subcloned the D and E domains of human RAR gamma for expression in Escherichia coli. Using nickel-chelation affinity chromatography with a polyhistidine amino-terminal tail, purification of the DE peptide with a pI of 5.18 was accomplished to greater than 98% purity. Scatchard analysis and fluorescence quenching techniques using the purified protein indicate a very high percentage of functional molecules ( > 95%) with a Kd for retinoic acid (t-RA) of 0.6 +/- 0.1 nM. Circular dichroism spectra of the purified domains predict a predominantly alpha-helical structure (approximately 56%) with little beta sheet present. No significant changes in these structural characteristics were observed upon binding of t-RA. Inspection of the amino acid sequence within these domains identified a single tryptophan residue at position 227. Modeling the amino acid sequence in this region as an alpha-helical structure indicates that this tryptophan is adjacent to alanine 234, which corresponds to alanine 225 in RAR beta that has previously been linked to the ligand binding site. Fluorescence of this tryptophan was quenched in a dose-dependent manner on the addition of t-RA, confirming that Trp-227 is within the ligand binding site. Tryptophan fluorescence quenching analysis also demonstrates that a single retinoic acid molecule is bound per receptor and suggests that receptor-ligand interactions occur within the amino-terminal portion of the predominantly alpha-helical ligand binding domain.

Refined solution structure of human profilin I.

Profilin is a ubiquitous eukaryotic protein that binds to both cytosolic actin and the phospholipid phosphatidylinositol-4,5-bisphosphate. These dual competitive binding capabilities of profilin suggest that profilin serves as a link between the phosphatidyl inositol cycle and actin polymerization, and thus profilin may be an essential component in the signaling pathway leading to cytoskeletal rearrangement. The refined three-dimensional solution structure of human profilin I has been determined using multidimensional heteronuclear NMR spectroscopy. Twenty structures were selected to represent the solution conformational ensemble. This ensemble of structures has root-mean-square distance deviations from the mean structure of 0.58 A for the backbone atoms and 0.98 A for all non-hydrogen atoms. Comparison of the solution structure of human profilin to the crystal structure of bovine profilin reveals that, although profilin adopts essentially identical conformations in both states, the solution structure is more compact than the crystal structure. Interestingly, the regions that show the most structural diversity are located at or near the actin-binding site of profilin. We suggest that structural differences are reflective of dynamical properties of profilin that facilitate favorable interactions with actin. The global folding pattern of human profilin also closely resembles that of Acanthamoeba profilin I, reflective of the 22% sequence identity and approximately 45% sequence similarity between these two proteins.

Immunologic and molecular characterization of Amb p V allergens from Ambrosia psilostachya (western Ragweed) pollen.

We have purified and characterized the Amb p V allergen (A1 variant) from western ragweed (Ambrosia psilostachya) pollen. This allergen was found to be highly cross-reactive with the Amb a VA1 allergen from short ragweed (A. artemisiifolia) pollen in a competitive double-Ab radioimmunoassay (DARIA) and the two allergens showed concordant allergenic potency in histamine-release experiments. We cloned and sequenced several Amb p V genes from western ragweed pollen and flowers by direct PCR of genomic DNA. The amino acid sequences deduced from the nucleotide sequences indicated the presence of multiple forms of Amb p V that could be broadly classified into two groups: Amb p VA and Amb p VB variants. The sequences of the Amb p VA variants are highly homologous to Amb a V (about 90% identity) and very similar to the protein sequence that we obtained. The Amb p VB variants share approximately 65% amino acid homology with Amb a V and have five to seven cysteine residues as compared with the eight found in Amb a V and Amb t V. Two cysteine residues that form disulfide bonds in other Amb Vs (positions 19 and 43 in Amb a V) are replaced by serine and alanine in the Amb p VB1 and Amb p VB2 variants. We have generated model structures of Amb p VA1, VA2, VA3, and VB1 variants from the nuclear magnetic resonance-derived structure of Amb a VA1 by homology modeling. Comparison of antigenic epitopes predicted for the structures of Amb p V variants and Amb a VA1 explains the observed cross-reactivity of the two ragweed proteins and suggests the epitopes likely to be involved in Ab recognition.

Identification of the poly-L-proline-binding site on human profilin.

Profilin is a ubiquitous protein that has been implicated in the signaling pathway leading to cytoskeletal rearrangement in cells. An unusual property of profilin is its high binding affinity for poly-L-proline (PLP). This binding property is conserved in the profilins from diverse species with little sequence homology. We have monitored the binding of PLP to profilin by fluorescence and nuclear magnetic resonance spectroscopies. NMR spectroscopy has identified several residues whose amide nitrogen and amide hydrogen chemical shifts are significantly perturbed by binding of PLP. The affected residues are located at various locations throughout profilin's primary structure; however, mapping the location of the affected residues onto the recently determined three-dimensional solution structure of human profilin indicates that the effects of PLP binding are highly localized. Poly-L-proline binds profilin at the hydrophobic interface between profilin's NH2- and COOH-terminal helices and the upper face of its antiparallel beta-sheet. In contrast, residues located on the opposite side of the profilin structure are unaffected. The extent of the potential interaction surface of the PLP-profilin complex suggests that as few as 6 contiguous prolines would be sufficient for binding profilin. Examination of sequence data bases indicates that stretches of prolines of this length and longer occur in numerous regulatory proteins, suggesting that the ability of profilin to bind polyproline may be an important component of its signaling capabilities.

Solution structure of an isolated antibody VL domain.

The solution structure of the isolated VL domain of the anti-digoxin antibody 26-10 has been determined using data derived from heteronuclear multi-dimensional nuclear magnetic resonance (n.m.r.) experiments. Analytical ultracentrifugation and n.m.r. data demonstrate that the VL domain is only weakly associating (Kd = 2.5 (+/- 0.7) mM) and that it experiences a rapid monomer/dimer equilibrium under the n.m.r. experimental conditions. Therefore, the results reported here represent the first structure determination of an antibody VL domain in the absence of fixed quaternary interactions. The structure determination is based on 930 proton-proton distance constraints, 113 dihedral angle constraints, and 46 hydrogen bond constraints. Eighty initial structures were calculated with the variable target function program DIANA; of these, 31 were accepted on the basis of satisfaction of constraints (no distance constraint violations > 0.5 A; target function < 3.0 A2). Accepted DIANA structures were refined by restrained energy minimization using the X-PLOR program. The 15 best energy-minimized DIANA structures were chosen as a representative ensemble of solution conformations. The average root-mean-square differences (r.m.s.d.) between the individual structures of this ensemble and the mean coordinates is 0.85 (+/- 0.10) A for all backbone atoms and 1.29 (+/- 0.10) A for all heavy atoms. For beta-strands A, B, C, D, E and F, the average backbone atom r.m.s.d. to the mean structure is 0.46 (+/- 0.06) A. A higher-resolution ensemble, with all backbone atom and all heavy atom r.m.s.d.s. to the mean coordinates of 0.54 (+/- 0.08) A and 0.98 (+/- 0.12) A, respectively, was obtained by X-PLOR simulated annealing refinement of the 15 energy-minimized DIANA structures. A detailed analysis of the original ensemble of 15 energy-minimized DIANA structures is presented, as this ensemble retains a broader, and possibly more realistic, sampling of conformation space. The backbone atom and all heavy atom r.m.s.d.s between the mean energy-minimized DIANA structure and the X-ray derived coordinates of the VL domain within the Fab/digoxin complex are 1.05 A and 1.56 A, respectively. Subtle differences between the solution and X-ray structures occur primarily in CDR2, CDR3, beta-strands A, F and G, and localized regions of hydrophobic packing. Overall, these results demonstrate that the 26-10 VL domain conformation is determined primarily by intradomain interactions, and that quaternary VL-VH association induces relatively minor conformational adjustments.

Characterization of the three-dimensional solution structure of human profilin: 1H, 13C, and 15N NMR assignments and global folding pattern.

Human profilin is a 15-kDa protein that plays a major role in the signaling pathway leading to cytoskeletal rearrangement. Essentially complete assignment of the 1H, 13C, and 15N resonances of human profilin have been made by analysis of multidimensional, double- and triple-resonance nuclear magnetic resonance (NMR) experiments. The deviation of the 13C alpha and 13C beta chemical shifts from their respective random coil values were analyzed and correlate well with the secondary structure determined from the NMR data. Twenty structures of human profilin were refined in the program X-PLOR using a total of 1186 experimentally derived conformational restraints. The structures converged to a root mean squared distance deviation of 1.5 A for the backbone atoms. The resultant conformational ensemble indicates that human profilin is an alpha/beta protein comprised of a seven-stranded, antiparallel beta-sheet and three helices. The secondary structure elements for human profilin are quite similar to those found in Acanthamoeba profilin I [Archer, S. J., Vinson, V. K., Pollard, T. D., & Torchia, D. A. (1993), Biochemistry 32, 6680-6687], suggesting that the three-dimensional structure of Acanthamoeba profilin I should be analogous to that determined here for human profilin. The structure determination of human profilin has facilitated the sequence alignment of lower eukaryotic and human profilins and provides a framework upon which the various functionalities of profilin can be explored. At least one element of the actin-binding region of human profilin is an alpha-helix. Two mechanisms by which phosphatidylinositol 4,5-bisphosphate can interfere with actin-binding by human profilin are proposed.

Relaxation study of the backbone dynamics of human profilin by two-dimensional 1H-15N NMR.

The dynamic properties of 111 backbone HN sites in uncomplexed human profilin, a protein of 139 residues, have been characterized by two-dimensional inverse-detected 1H-15N NMR spectroscopy. Heteronuclear (1H)-15N nuclear Overhauser effects and 15N longitudinal and transverse relaxation rates have been analyzed in terms of model-free spectral density functions and exchange contributions to transverse relaxation rates. Relatively high mobilities on the nanosecond time-scale are observed for Asp26 and Ser27, which form part of a loop connecting beta-strands A and B, and for Thr92 through Ala95, which are in a loop connecting beta-strands E and F. Significant exchange contributions, indicative of motions on the microsecond to millisecond time-scale, have been obtained for 30 residues. These include Leu77, Asp80 and Gly81 of a loop between beta-strands D and E, Ser84 and Met85 of beta-strand E, Gly121 of a loop connecting beta-strand G and the C-terminal helix, and Gln138, which is next to the C-terminal residue Tyr139. Some of the regions showing high flexibility in profilin are known to be involved in poly-L-proline binding.

Solution structure of the calcium channel antagonist omega-conotoxin GVIA.

The three-dimensional solution structure is reported for omega-conotoxin GVIA, which is a potent inhibitor of presynaptic calcium channels in vertebrate neuromuscular junctions. Structures were generated by a hybrid distance geometry and restrained molecular dynamics approach using interproton distance, torsion angle, and hydrogen-bonding constraints derived from 1H NMR data. Conformations of GVIA with low constraint violations converged to a common peptide fold. The secondary structure in the peptide is an antiparallel triple-stranded beta-sheet containing a beta-hairpin and three tight turns. The NMR data are consistent with the region of the peptide from residues S9 to C16 being more dynamic than the rest of the peptide. The peptide has an amphiphilic structure with a positively charged hydrophilic side and an opposite side that contains a small hydrophobic region. Residues that are thought to be important in binding and function are located on the hydrophilic face of the peptide.

Expression and analysis of recombinant Amb a V and Amb t V allergens. Comparison with native proteins by immunological assays and NMR spectroscopy.

The Amb V allergens are small, highly disulfide-bonded ragweed pollen allergens that serve as useful models for understanding the molecular basis of the human immune response. We have produced recombinant Amb a V and Amb t V (from short and giant ragweed pollens, respectively) in Escherichia coli and have compared their structural and functional characteristics to those of the native proteins. Recombinant Amb t V was indistinguishable from native Amb t V as determined by NMR spectroscopy and antibody-binding studies. Whereas inhibition analysis showed that recombinant Amb a V possessed only approximately 50% of the antibody-binding activity of native Amb a V, the two proteins were similarly effective in stimulating Amb a V-specific T-cells. Our results demonstrate that even highly homologous proteins exhibit different abilities to fold into their native three-dimensional conformations and establish the potential and limits of expressing the recombinant Amb V allergens intracellularly in E. coli.

Proton resonance assignments and three-dimensional solution structure of the ragweed allergen Amb a V by nuclear magnetic resonance spectroscopy.

Essentially complete assignment of the proton resonances in the allergenic protein Amb a V has been made by analysis of two-dimensional NMR experiments. Conformational constraints were obtained in three forms: interproton distances derived from NOE cross-peak intensities of NOESY spectra, torsion angle constraints derived from J-coupling constants of COSY and PE-COSY spectra, and hydrogen bond constraints derived from hydrogen-exchange experiments. Conformations of Amb a V with low constraint violations were generated using dynamic simulated annealing in the program XPLOR. The refined structures are comprised of a C-terminal alpha-helix, a small segment of antiparallel beta-sheet, and several loops. A hydrophobic core exists at the interface of the alpha-helix and beta-sheet. The derived structure accounts for the several anomalous proton chemical shifts that are observed. The structure determined here for Amb a V is topologically similar to the structure determined previously for the homologous allergenic protein Amb t V [Metzler, W. J., Valentine, K., Roebber, M., Friedrichs, M. S., Marsh, D., & Mueller, L. (1992) Biochemistry 31, 5117-5127]; however, significant differences exist in the packing of side chains in the hydrophobic core of the molecules. Comparison of the detailed structural features of these two proteins will allow us to suggest surface substructures for the Amb V allergens that are likely to participate in B cell epitopes.

Determination of the three-dimensional solution structure of ragweed allergen Amb t V by nuclear magnetic resonance spectroscopy.

Analysis of two-dimensional NMR experiments has afforded essentially complete assignment of all proton resonances in the allergenic protein Amb t V. Conformational constraints were obtained from the NMR data in three forms: interproton distances derived from NOE cross-peak intensities of NOESY spectra, torsion angle constraints derived from J-coupling constants of COSY and PE-COSY spectra, and hydrogen bond constraints derived from hydrogen-exchange experiments. Conformations of Amb t V with low constraint violations were generated using dynamic simulated annealing in the program XPLOR. The refined structures are comprised of a C-terminal alpha-helix, a short stretch of triple-stranded antiparallel beta-sheet, and several loops. In addition, the cystine partners of the four disulfide linkages (for which there are no biochemical data) have been assigned. The refined structures of Amb t V will allow us to suggest surface substructures for the Amb V allergens that are likely to participate in B cell epitopes and will assist us in defining the Ia/T cell epitopes that interact with the MHC class II (or Ia) molecule and the T cell receptor leading to the induction of the immune response to Amb t V.

Determining local conformational variations in DNA. Nuclear magnetic resonance structures of the DNA duplexes d(CGCCTAATCG) and d(CGTCACGCGC) generated using back-calculation of the nuclear Overhauser effect spectra, a distance geometry algorithm and constrained molecular dynamics.

Two-dimensional nuclear magnetic resonance (n.m.r.) spectroscopy and a variety of computational techniques have been used to generate three-dimensional structures of the two DNA duplexes d(CGCCTAATCG) and d(CGTCACGCGC). The central six base-pairs in these two decamers contain all ten dinucleotide pairs in DNA and thus, represent a model system for investigating how the local structure of DNA varies with base sequence. Resonance assignments were made for the non-exchangeable base protons and most of the C-1'-C-4' sugar protons in both decamers. Three-dimensional structures were generated using a distance geometry algorithm and these initial structures were refined by optimizing the fit of back-calculated spectra against the experimental two-dimensional nuclear Overhauser effect (NOE) spectra. This back-calculation procedure consists of calculating NOE cross relaxation rates for a given structure by solution of the Bloch equations, and directly accounts for spin diffusion effects. Use of this refinement procedure eliminates some assumptions that have been invoked when generating structures of DNA oligomers from n.m.r. data. Constrained energy minimization and constrained quenched molecular dynamics calculation were also performed on both decamers to help generate energetically favorable structures consistent with the experimental data. Analysis of the local conformational parameters of helical twist, helical rise, propeller twist, displacement and the alpha, beta, gamma, epison and zeta backbone torsion angles in these structures shows that these parameters span a large range of values relative to the X-ray data of nucleic acids. However, the glycosidic and pseudorotation angles are quite well defined in these structures. The implications that these results have for determination of local structural variations of DNA in solution, such as those predicted by Callidine's rules, are discussed. Our results differ significantly from some previous studies on determining local conformations of nucleic acids and comparisons with these studies are made.