Solution structure of the EcoRI DNA octamer containing 5-fluorouracil via restrained molecular dynamics using distance and torsion angle constraints extracted from NMR spectral simulations.

The self-complementary DNA octamer [d(GGAATUFCC)]2, containing the EcoRI recognition sequence with one of the thymines replaced by 5-fluorouracil (UF), was synthesized. Proton homonuclear two-dimensional nuclear Overhauser effect (2D NOE) and double-quantum-filtered correlation (2QF-COSY) spectra, as well as one-dimensional spectra at different temperatures, were recorded for the octamer. Consequently, all proton resonances were assigned. The thermally induced transition from the duplex to single strands has been followed, demonstrating the stability of the duplex containing 5-fluorouracil. Simulations of the 2QF-COSY cross-peaks by means of the programs SPHINX and LINSHA were compared with experimental data, establishing scalar coupling constants for the sugar ring protons and hence sugar pucker parameters. The deoxyribose rings exhibit a dynamic equilibrium of N- and S-type conformers with 75-95% populations of the latter. Two programs used for complete relaxation matrix analysis 2D NOE spectra, CORMA and MARDIGRAS, were modified to account for the influence of the fluorines on dipolar interactions in the proton system. Quantitative assessment of the 2D NOE cross-peak intensities for different mixing times, in conjunction with the program MARDIGRAS, gave a set of interproton distances for each mixing time. The largest and smallest values of each of the interproton distances were chosen as the upper and lower bounds for each distance constraint. The distance bounds define the size of a flat-well potential function term, incorporated into the AMBER force field, which was employed for restrained molecular dynamics calculations. Torsion angle constraints in the form of a flat-well potential were also constructed from the analysis of the sugar pucker data. Several restrained molecular dynamics runs of 35 ps were performed, utilizing 284 experimental distance and torsion angle constraints and two different starting structures, energy-minimized A- and B-DNA. Convergence to similar structures with a root-mean-square deviation of 1.2 A was achieved for the central hexamer of the octamer, starting from A- and B-DNA. The average structure from six different molecular dynamics runs was subjected to final restrained energy minimization. The resulting final structure was in good agreement with the structures derived from different molecular dynamics runs and showed a substantial improvement of the 2D NOE sixth-root residual index in comparison with classical and energy-minimized B-DNA. A detailed analysis of the conformation of the final structure and comparison with structures of similar sequences, obtained by different methods, were performed.     

Solution structure of [d(GTATATAC)]2 via restrained molecular dynamics simulations with nuclear magnetic resonance constraints derived from relaxation matrix analysis of two-dimensional nuclear Overhauser effect experiments

Two-dimensional nuclear Overhauser effect (2D NOE) spectra have been used as the experimental basis for determining the solution structure of the duplex [d(GTATATAC)]2 employing restrained molecular dynamics (rMD) simulations. The MARDIGRAS algorithm has been employed to construct a set of 233 interproton distance constraints via iterative complete relaxation matrix analysis utilizing the peak intensities from the 2D NOE spectra obtained for different mixing times and model structures. The upper and lower bounds for each of the constraints, defining size of a flat-well potential function term used in the rMD simulations, were conservatively chosen as the largest or smallest value calculated by MARDIGRAS. Three different starting models were utilized in several rMD calculations: energy-minimized A-DNA, B-DNA, and a structure containing wrinkled D-DNA in the interior. Considerable effort was made to define the appropriate force constants to be employed with the NOE terms in the AMBER force field, using as criteria the average constraints deviation, the constraints violation energy and the total energy. Of the 233 constraints, one was generated indirectly, but proved to be crucial in defining the structure: the cross-strand A5-H2 A5-H2 distance. As those two protons resonate isochronously for the self-complementary duplex, the distance cannot be determined directly. However, the general pattern of 2D NOE peak intensities, spin-lattice relaxation time (T1) values, and 31P nuclear magnetic resonance spectra lead to use of the A3-H2 A7-H2 distance for A5-H2 A5-H2 as well. Five rMD runs, with different random number seeds, were made for each of the three starting structures with the full distance constraint set. The average structure from all 15 runs and the five-structure averages from each starting structure were all quite similar. Two rMD runs for each starting structure were made with the A5-H2 A5-H2 constraint missing. The average of these six rMD runs revealed differences in structure, compared to that with the full set of constraints, primarily for the middle two base-pairs involving the missing cross-strand constraint but global deviations also were found. Conformational analysis of the resulting structures revealed that the inner four to six base-pairs differed in structure from the termini. Furthermore, an alternating structure was suggested with features alternating for the A-T and T-A steps.     

Structure of DNA sequence d-TGATCA by two-dimensional nuclear magnetic resonance spectroscopy and restrained molecular dynamics

The 5' d-TpG 3' element is a part of DNA sequences involved in regulation of gene expression and is also a site for intercalation of several anticancer drugs. Solution conformation of DNA duplex d-TGATCA containing this element has been investigated by two-dimensional NMR spectroscopy. Using a total of 12 torsional angles and 121 distance constraints, structural refinement has been carried out by restrained molecular dynamics (rMDs) in vacuum up to 100 ps. The structure is characterized by a large positive roll at TpG/CpA base pair step and large negative propeller twist for AT and TA base pairs. The backbone torsional angle, gamma(O5'-C5'-C4'-C3'), of T1 residue adopts a trans-conformation which is corroborated by short intra nucleotide T1H6-T1H5' (3.7A) distance in nuclear overhauser effect spectroscopy (NOESY) spectra while the backbone torsional angle, beta(P-O5'-C5'-C4'), exists in trans as well as gauche state for T1 and C5 residues. There is evidence of significant flexibility of the sugar-phosphate backbone with rapid inter-conversion between two different conformers at TpG/CpA base pair step. The base sequence dependent variations and local structural heterogeneity have important implications in specific recognition of DNA by ligands.     

Structure of DNA hexamer sequence d-CGATCG by two-dimensional nuclear magnetic resonance spectroscopy and restrained molecular dynamics

Solution conformation of self-complementary DNA duplex d-CGATCG, containing 5' d-CpG 3' site for intercalation of anticancer drug, daunomycin and adriamycin, has been investigated by nuclear magnetic resonance (NMR) spectroscopy. Complete resonance assignments of all the protons (except some H5'/H5" protons) have been obtained following standard procedures based on double quantum filtered correlation spectroscopy (dQF COSY) and two-dimensional nuclear Overhauser effect (NOE) spectra. Analysis of sums of coupling constants in one-dimensional NMR spectra, cross peak patterns in dQF COSY spectra and inter proton distances shows that the DNA sequence assumes a conformation close to the B-DNA family. The deoxyribose sugar conformation is in dynamic equilibrium with predominantly S-type conformer and a minor N-type conformer with N<-->S equilibrium varying with temperature. At 325 K, the mole fraction of the N-conformer increases for some of the residues by approximately 9%. Using a total of 10 spin-spin coupling constants and 112 NOE intensities, structural refinement has been carried out using Restrained Molecular Dynamics (rMD) with different starting structures, potential functions and rMD protocols. It is observed that pseudorotation phase angle of deoxyribose sugar for A3 and T4 residues is approximately 180 degrees and approximately 120 degrees, respectively while all other residues are close to C2'endo-conformation. A large propeller twist (approximately -18 degrees) and smallest twist angle (approximately 31 degrees) at A3pT4 step, in the middle of the sequence, a wider (12 A) and shallower (3.0 A) major groove with glycosidic bond rotation as high anti at both the ends of hexanucleotide are observed. The structure shows base-sequence dependent variations and hence strong local structural heterogeneity, which may have implications in ligand binding.     

Solution structure of neurotoxin I from the sea anemone Stichodactyla helianthus. A nuclear magnetic resonance, distance geometry, and restrained molecular dynamics study

The three-dimensional structure of the sea anemone polypeptide Stichodactyla helianthus neurotoxin I in aqueous solution has been determined using distance geometry and restrained molecular dynamics simulations based on NMR data acquired at 500 MHz. A set of 470 nuclear Overhauser enhancement values was measured, of which 216 were used as distance restraints in the structure determination along with 15 dihedral angles derived from coupling constants. After restrained molecular dynamics refinement, the eight structures that best fit the input data form a closely related family. They describe a structure that consists of a core of twisted, four-stranded, antiparallel beta-sheet encompassing residues 1-3, 19-24, 29-34, and 40-47, joined by three loops, two of which are well defined by the NMR data. The third loop, encompassing residues 7-16, is poorly defined by the data and is assumed to undergo conformational averaging in solution. Pairwise root mean square displacement values for the backbone heavy atoms of the eight best structures are 1.3 +/- 0.2A when the poorly defined loop is excluded and 3.6 +/- 1.0A for all backbone atoms. Refinement using restrained molecular dynamics improved the quality of the structures generated by distance geometry calculations with respect to the number of nuclear Overhauser enhancements violated, the size of the total distance violations and the total potential energies of the structures. The family of structures for S. heliathus neurotoxin I is compared with structures of related sea anemone proteins that also bind to the voltage-gated sodium channel.     

Solution structure of human insulin-like growth factor 1: a nuclear magnetic resonance and restrained molecular dynamics study

The solution structure of human insulin-like growth factor 1 has been investigated with a combination of nuclear magnetic resonance and restrained molecular dynamics methods. The results show that the solution structure is similar to that of insulin, but minor differences exist. The regions homologous to insulin are well-defined, while the remainder of the molecule exhibits greater disorder. The resultant structures have been used to visualize the sites for interaction with a number of physiologically important proteins.     

Solution conformation of purine-pyrimidine DNA octamers using nuclear magnetic resonance, restrained molecular dynamics and NOE-based refinement

The solution structures of two alternating purine-pyrimidine octamers, [d(G-T-A-C-G-T-A-C)]2 and the reverse sequence [d(C-A-T-G-C-A-T-G)]2, are investigated by using nuclear magnetic resonance spectroscopy and restrained molecular dynamics calculations. Chemical shift assignments are obtained for non-exchangeable protons by a combination of two-dimensional correlation and nuclear Overhauser enhancement (NOE) spectroscopy experiments. Distances between protons are estimated by extrapolating distances derived from time-dependent NOE measurements to zero mixing time. Approximate dihedral angles are determined within the deoxyribose ring from coupling constants observed in one and two-dimensional spectra. Sets of distance and dihedral determinations for each of the duplexes form the bases for structure determination. Molecular dynamics is then used to generate structures that satisfy the experimental restraints incorporated as effective potentials into the total energy. Separate runs start from classical A and B-form DNA and converge to essentially identical structures. To circumvent the problems of spin diffusion and differential motion associated with distance measurements within molecules, models are improved by NOE-based refinement in which observed NOE intensities are compared to those calculated using a full matrix analysis procedure. The refined structures generally have the global features of B-type DNA. Some, but not all, variations in dihedral angles and in the spatial relationships of adjacent base-pairs are observed to be in synchrony with the alternating purine-pyrimidine sequence.     

Solution conformation of human neuropeptide Y by 1H nuclear magnetic resonance and restrained molecular dynamics.

The solution structure of human neuropeptide Y has been solved by conventional two-dimensional NMR techniques followed by distance-geometry and molecular-dynamics methods. The conformation obtained is composed of two short contiguous alpha-helices comprising residues 15-26 and 28-35, linked by a hinge inducing a 100 degree angle. The first helix (15-26) is connected to a polyproline stretch (residues 1-10) by a tight hairpin (residues 11-14). The helices and the polyproline stretch are packed together by hydrophobic interactions. This structure is related to that of the homologous avian pancreatic polypeptide and bovine pancreatic polypeptide. The C- and N-terminii, known to be involved in the biological activity for respectively the receptor binding and activation, are close together in space. The side chains of residues Arg33, Arg35 and Tyr36 on the one hand, and Tyr1 and Pro2 on the other, form a continuous solvent-exposed surface of 4.9 mm2 which is supposed to interact with the receptor for neuropeptide Y.     

Three-dimensional structure of bradykinin in SDS micelles. Study using nuclear magnetic resonance, distance geometry, and restrained molecular mechanics and dynamics

The conformational properties of bradykinin in five molar excess sodium dodecyl sulfate (SDS) micelles have been examined by two-dimensional nuclear magnetic resonance (NMR) techniques at 500 MHz. Detailed structural information for bradykinin in SDS was obtained from quantitative 2-D nuclear Overhauser enhancement (n.O.e.) analyses, distance geometry, and restrained molecular mechanics and dynamics calculations. The conformation of bradykinin in SDS micelles, as determined by these methods, is characterized by a beta-turn-like structure at residues 6-9. A detailed comparison of the structures derived from distance geometry and restrained molecular mechanics and dynamics is also presented.     

Solution structure studies of d(AC)4.d(GT)4 via restrained molecular dynamics simulations with NMR constraints derived from two-dimensional NOE and double-quantum-filtered COSY experiments

The structure of d(AC)4.d(GT)4 is investigated by constrained molecular dynamics simulations. The constraints include proton pair distances derived from 2D NOE intensities by using the iterative relaxation matrix analysis algorithm MARDIGRAS and sugar pucker phases and amplitudes derived from double-quantum-filtered COSY spectra. Molecular dynamics runs on simulated intensity and distance sets as well as the experimental data were carried out to determine the effects of starting structure, distance constraint derivation, energy functions, and experimental errors on the end result. It was found that structural details could not be elucidated within about 1.5-A overall atomic deviation. This limitation is due in part to the accuracy of the experimental data but, more importantly, is attributable to the quantity of experimental constraints available and to imperfections in the force field utilized in the molecular dynamics calculations. Within the limits of the method, some structural characteristics of d(AC)4.d(GT)4 could be elucidated.     

Three-dimensional structure of the wild-type lac Pribnow promoter DNA in solution. Two-dimensional nuclear magnetic resonance studies and distance geometry calculations

The solution structure of a 12 base-pair DNA duplex containing the wt-lac promoter Pribnow sequence TATGTT has been studied by two-dimensional nuclear magnetic resonance spectroscopy. Proton assignments for the 24 sugar and base residues were obtained from two-dimensional correlated nuclear magnetic resonance and two-dimensional nuclear Overhauser effect spectra in both 2H2O and H2O, and by two-dimensional relayed coherence transfer nuclear magnetic resonance spectroscopy experiments. Time-dependent, two-dimensional nuclear Overhauser effect spectra were used to determine the initial cross-relaxation rates between 212 pairs of assigned protons, leading to 212 interproton distances in the double helix (8 to 9 per nucleotide). These distance constraints, and known bond lengths and angles, were entered into a distance matrix. After smoothing the bounds of the distance matrix, 12 trial matrices within the bounds constraints were independently generated and embedded in three-dimensional space using a distance geometry algorithm, to generate 12 trial structures. These trial structures were then refined until they no longer violated the distance matrix. The resulting structures are very similar at the local base-pair and nearest-neighbor base-pair level, but exhibit increasing variation at more distant and global levels. At the nearest-neighbor level, the A to T step and the G to T step within the Pribnow hexamer, as well as the G to T step preceding the hexamer, all exhibit very low screw pitch, i.e. 5(+/- 6) degrees. Conversely, the T to G step in the center of the promoter has a large screw pitch (47(+/- 2) degrees) and the T to G step at the 3' end of the promoter has a very large screw pitch (60(+/- 3) degrees). The limitations of nuclear magnetic resonance spectroscopy distance determination of structure are discussed in terms of resolution and spectral overlap of two-dimensional nuclear Overhauser effect crosspeaks. In the present duplex, the inability to measure several 1'-2' and 1'-2" distances resulted in underdetermination of the precise local sugar conformation for seven of the 24 residues, although the spatial position of all sugars was well defined.     

Three-dimensional structure of potato carboxypeptidase inhibitor in solution. A study using nuclear magnetic resonance, distance geometry, and restrained molecular dynamics

The solution conformation of potato carboxypeptidase inhibitor (CPI) has been investigated by 1H NMR spectroscopy. The spectrum is assigned in a sequential manner by using two-dimensional NMR techniques to identify through-bond and through-space (less than 5 A) connectivities. A set of 309 approximate interproton distance restraints is derived from the two-dimensional nuclear Overhauser enhancement spectra and used as the basis of a three-dimensional structure determination by a combination of metric matrix distance geometry and restrained molecular dynamics calculations. A total of 11 converged distance geometry structures were computed and refined by using restrained molecular dynamics. The average atomic root mean square (rms) difference between the final 11 structures and the mean structure obtained by averaging their coordinates is 1.4 +/- 0.3 A for residues 2-39 and 0.9 +/- 0.2 A for residues 5-37. The corresponding values for all atoms are 1.9 +/- 0.3 and 1.4 +/- 0.2 A, respectively. The larger values for residues 2-38 relative to those for residues 5-37 arise from the fact that the positions of the N- (residues 1-4) and C- (residues 38-39) terminal tails are rather poorly determined, whereas those of the core of the protein (residues 5-37) are well determined by the experimental interproton distance data. The computed structures are very close to the X-ray structure of CPI in its complex with carboxypeptidase, and the backbone atomic rms difference between the mean of the computed structures and the X-ray structure is only 1.2 A. Nevertheless, there are some real differences present which are evidenced by significant deviations between the experimental upper interproton distance limits and the corresponding interproton distances derived from the X-ray structure. These principally occur in two regions, residues 18-20 and residues 28-30, the latter comprising part of the region of secondary contacts between CPI and carboxypeptidase in the X-ray structure.     

Hairpin structures in DNA containing arabinofuranosylcytosine. A combination of nuclear magnetic resonance and molecular dynamics

Nuclear magnetic resonance (NMR) and model-building studies were carried out on the hairpin form of the octamer d(CGaCTAGCG) (aC = arabinofuranosylcytosine), referred to as the TA compound. The nonexchangeable protons of the TA compound were assigned by means of nuclear Overhauser effect spectroscopy (NOESY) and correlated spectroscopy (COSY). From a detailed analysis of the coupling data and of the NOESY spectra the following conclusions are reached: (i) The hairpin consists of a stem of three Watson-Crick type base pairs, and the two remaining residues, T(4) and dA(5), participate in a loop. (ii) All sugar rings show conformational flexibility although a strong preference for the S-type (C2'-endo) conformer is observed. (iii) The thymine does not stack upon the 3' side of the stem as expected, but swings into the minor groove. (This folding principle of the loop involves an unusual alpha t conformer in residue T(4).) (iv) At the 5'-3' loop-stem junction a stacking discontinuity occurs as a consequence of a sharp turn in that part of the backbone, caused by the unusual beta + and gamma t torsion angles in residue dG(6). (v) The A base slides over the 5' side of the stem to stack upon the aC(3) residue at the 3' side of the stem in an antiparallel fashion. On the basis of J couplings and a set of approximate proton-proton distances from NOE cross peaks, a model for the hairpin was constructed. This model was then refined by using an iterative relaxation matrix approach (IRMA) in combination with restrained molecular dynamics calculations. The resulting final model satisfactorily explains all the distance constraints.     

The three-dimensional structure of alpha1-purothionin in solution: combined use of nuclear magnetic resonance, distance geometry and restrained molecular dynamics

The determination of the three-dimensional solution structure of α1-purothionin using a combination of metric matrix distance geometry and restrained molecular dynamics calculations based on n.m.r. data is presented. The experimental data comprise complete sequence-specific proton resonance assignments, a set of 310 approximate interproton distance restraints derived from nuclear Overhauser effects, 27 Ø backbone torsion angle restraints derived from vicinal coupling constants, 4 distance restraints from hydrogen bonds and 12 distance restraints from disulphide bridges. The average atomic rms difference between the final nine converged structures and the mean structure obtained by averaging their coordinates is 1.5 ± 0.1 å for the backbone atoms and 2.0 ± 0.1 å for all atoms. The overall shape of α1-purothionin is that of the capital letter L, similar to that of crambin, with the longer arm comprising two approximately parallel α-helices and the shorter arm a strand and a mini anti-parallel β sheet.     

The determination of the three-dimensional structure of barley serine proteinase inhibitor 2 by nuclear magnetic resonance, distance geometry and restrained molecular dynamics

The solution structure of the 64 residue structured domain (residues 20-83) of barley serine proteinase inhibitor 2 (BSPI-2) is determined on the basis of 403 interproton distance, 34 phi backbone torsion angle and 26 hydrogen bonding restraints derived from n.m.r. measurements. A total of 11 converged structures were computed using a metric matrix distance geometry algorithm and refined by restrained molecular dynamics. The average rms difference between the final 11 structures and the mean structure obtained by averaging their coordinates is 1.4 +/- 0.2 A for the backbone atoms and 2.1 +/- 0.1 A for all atoms. The overall structure, which is almost identical to that found by X-ray crystallography, is disc shaped and consists of a central four component mixed parallel and antiparallel beta-sheet flanked by a 13 residue alpha-helix on one side and the reactive site loop on the other.     

Solution conformations of a trimannoside from nuclear magnetic resonance and molecular dynamics simulations

N-linked oligosaccharides often act as ligands for receptor proteins in a variety of cell recognition processes. Knowledge of the solution conformations, as well as protein-bound conformations, of these oligosaccharides is required to understand these important interactions. In this paper we present a model for the solution conformations sampled by a simple trimannoside, methyl 3, 6-di-O-(alpha-D-mannopyranosyl)-alpha-D-mannopyranoside, which contains two of the most commonly found glycosidic linkages in N-linked oligosaccharides. This model was derived from simulated annealing protocols incorporating distance restraints extracted from NOESY spectra along with torsional restraints computed from three-bond (1)H-(13)C coupling constants measured across the glycosidic bonds. The model was refined in light of unrestrained molecular dynamics simulations conducted in the presence of solvent water. The resulting model depicts a molecule undergoing conformational averaging in solution, adopting four major and two minor conformations. The four major conformations arise from a pair of two-state transitions, one each at the alpha(1-->3) and alpha(1-->6) linkages, whereas the minor conformations result from an additional transition of the alpha(1-->6) linkage. Our data also suggest that the alpha(1-->3) transition is fast and changes the molecular shape slightly, whereas the alpha(1-->6) is much slower and alters the molecular shape dramatically.     

Determination of protein structures from nuclear magnetic resonance data using a restrained molecular dynamics approach: the lac repressor DNA binding domain

A procedure is described to determine from NMR data the three-dimensional structure of biomolecules. This procedure combines model building with a restrained Molecular Dynamics algorithm, in which distance information from NOEs is incorporated in the form of pseudo potentials. The method has been applied to the N-terminal DNA-binding domain or "headpiece" (amino acids 1-51) of the lac repressor from E. coli, for which no crystal structure is available. The spatial structure of the headpiece is discussed in terms of known physical and biochemical data and of its DNA binding properties.     

The conformations of hirudin in solution: a study using nuclear magnetic resonance, distance geometry and restrained molecular dynamics

The solution conformations of the protein hirudin have been investigated by the combined use of distance geometry and restained molecular dynamics calculations. The basis for the structure determination comprised 359 approximate inter-proton distance restrains and 10 phi backbone torsion angle restrains derived from n.m.r. measurements. It is shown that hirudin is composed of three domains: a central core made up of residues 3-30, 37-46 and 56-57; a protruding 'finger' (residues 31-36) consisting of the tip of an antiparallel beta sheet, and an exposed loop (residues 47-55). The structure of each individual domain is relatively well defined with average backbone atomic r.m.s. differences of <2 A between the final seven converged restrained dynamic structures and the mean structure obtained by averaging their coordinates. The orientation of the two minor domains relative to the central core, however, could not be determined as no long-range (i-h >5) interdomain proton-proton contacts could be observed in the two-dimensional nuclear Overhauser enhancement spectra. From the restrained molecular dynamics calculations it appears that the two minor domains exhibit large rigid-body motions relative to the central core.     

Combined procedure of distance geometry and restrained molecular dynamics techniques for protein structure determination from nuclear magnetic resonance data: application to the DNA binding domain of lac repressor from Escherichia coli

The technique of two-dimensional nuclear magnetic resonance (2D-NMR) has recently assumed an active role in obtaining information on structures of polypeptides, small proteins, sugars, and DNA fragments in solution. In order to generate spatial structures from the atom-atom distance information obtained by the NMR method, different procedures have been developed. Here we introduce a combined procedure of distance geometry (DG) and molecular dynamics (MD) calculations for generating 3D structures that are consistent with the NMR data set and have reasonable internal energies. We report the application of the combined procedure on the lac repressor DNA binding domain (headpiece) using a set of 169 NOE and 17 "hydrogen bond" distance constraints. Eight of ten structures generated by the distance geometry algorithm were refined within 10 ps MD simulation time to structures with low internal energies that satisfied the distance constraints. Although the combination of DG and MD was designed to combine the good sampling properties of the DG algorithm with an efficient method of lowering the internal energy of the molecule, we found that the MD algorithm contributes significantly to the sampling as well.     

The solution structure of a DNA hairpin containing a loop of three thymidines determined by nuclear magnetic resonance and molecular mechanics.

We have determined by two-dimensional nuclear magnetic resonance studies and molecular mechanics calculations the three-dimensional solution structure of a 21 residue oligonucleotide capable of forming a hairpin structure with a loop of three thymidine residues. This structure is in equilibrium with a duplex form. At 33 degrees C, low ionic strength and in the presence of MgCl2 the hairpin form dominates in solution. Six Watson-Crick base pairs are formed topped by the loop structure. The residues 1-3 and 18-21 are not complementary and form dangling ends. Distance constraints have been derived from nuclear Overhauser enhancement measurements. These, together with molecular mechanics calculations, have been used to determine the structure. We do not observe stacking of thymidine residues either over the 3' or the 5' end of the stem.