Conformational studies of nucleic acids. V. Sequence specificities in the conformational energetics of oligonucleotides: the homo-tetramers

Duplex DNA Retains the Conformational Features of Single Strands: Perspectives from MD Simulations and Quantum Chemical Computations

Molecular dynamics simulations and geometry optimizations carried out at the quantum level as well as by quantum mechanical/molecular mechanics methods predict that short, single-stranded DNA oligonucleotides adopt conformations very similar to those observed in crystallographic double-stranded B-DNA, with rise coordinates close to ≈3.3 Å. In agreement with the experimental evidence, the computational results show that DNA single strands rich in adjacent purine nucleobases assume more regular arrangements than poly-thymine. The preliminary results suggest that single-stranded poly-cytosine DNA should also retain a substantial helical order in solution. A comparison of the structures of single and double helices confirms that the B-DNA motif is a favorable arrangement also for single strands. Indeed, the optimal geometry of the complementary single helices is changed to a very small extent in the formation of the duplex.

Single-Stranded DNA Oligonucleotides Retain Rise Coordinates Characteristic of Double Helices

The structures of single-stranded DNA oligonucleotides from dimeric to hexameric sequences have been thoroughly investigated. Computations performed at the density functional level of theory including dispersion forces and solvation show that single-stranded helices adopt conformations very close to crystallographic B-DNA, with rise coordinates amounting up to 3.3 Å. Previous results, suggesting that single strands should be shorter than double helices, largely originated from the incompleteness of the adopted basis set. Although sensible deviations with respect to standard B-DNA are predicted, computations indicate that sequences rich in stacked adenines are the most ordered ones, favoring the B-DNA pattern and inducing regular arrangements also on flanking nucleobases. Several structural properties of double helices rich in adenine are indeed already reflected by the corresponding single strands.

Conformations of the sugar-phosphate backbone in helical DNA crystal structures

The DNA backbone geometry was analyzed for 96 crystal structures of oligodeoxynucleotides. The ranges and mean values of the torsion angles for the observed subclasses of the A-, B-, and Z-DNA conformational types were determined by analyzing distributions of the torsion angles and scattergrams relating pairs of angles.

A complete conformational map for RNA

A simple stereochemical framework for understanding RNA structure has remained elusive to date. We present a comprehensive conformational map for two nucleoside-5',3'-diphosphates and for a truncated dinucleotide derived from a grid search of all potential conformers using hard sphere steric exclusion criteria to define allowed conformers. The eight-dimensional conformational space is presented as a series of two-dimensional projections. These projections reveal several well-defined allowed and disallowed regions which correlate well with data obtained from X-ray crystallography of both large and small RNA molecules. Furthermore, the two-dimensional projections show that consecutive and ribose ring-proximal torsion angles are interdependent, while more distant torsion angles are not. Remarkably, using steric criteria alone, it is possible to generate a predictive conformational map for RNA.

Creating RNA bulges: cleavage of RNA in RNA/DNA duplexes by metal ion catalysis

The manipulation of a single-stranded RNA target by forming different RNA/antisense hybrids demonstrates the possibility of cleaving the RNA strand within duplexes. This was achieved using the sequence composition of the antisense oligonucleotide, an approach that results in various bulges [unpaired base(s)] in the RNA target, which is then cleavable at these specific bulge sites under free metal ion or metal complex catalysis. RNA cleavages promoted by metal ions were performed under mild conditions and characterized by separating the RNA fragments carrying end label. The observed products result from intramolecular transesterification causing RNA strand scission. No detectable cleavage of the RNA was observed with either a fully complementary RNA/antisense hybrid or a bulged base in the antisense strand. A molecular modeling study of the RNA backbone suggests that the local conformation of the RNA backbone at a bulge in such hybrid duplexes greatly facilitates the metal-assisted catalytic cleavage. Endonucleolytic RNA cleavage within an RNA/antisense hybrid by metal complexes attached to the antisense oligonucleotide might lead to a new approach in antisense technology with artificial ribonucleases which operate with catalytic turnover.

A Monte Carlo method for generating structures of short single-stranded DNA sequences

A Monte Carlo method has been developed for generating the conformations of short single-stranded DNAs from arbitrary starting states. The chain conformers are constructed from energetically favorable arrangements of the constituent mononucleotides. Minimum energy states of individual dinucleotide monophosphate molecules are identified using a torsion angle minimizer. The glycosyl and acyclic backbone torsions of the dimers are allowed to vary, while the sugar rings are held fixed in one of the two preferred puckered forms. A total of 108 conformationally distinct states per dimer are considered in this first stage of minimization. The torsion angles within 5 kcal/mole of the global minimum in the resulting optimized states are then allowed to vary by +/- 10 degrees in an effort to estimate the breadth of the different local minima. The energies of a total of 2187 (3(7)) angle combinations are examined per local conformational minimum. Finally, the energies of all dinucleotide conformers are scaled so that the populations of differently puckered sugar rings in the theoretical sample match those found in nmr solution studies. This last step is necessitated by limitations in the theoretical methods to predict DNA sugar puckering accurately. The conformer populations of the individual acyclic torsion angles in the composite dimer ensembles are found to be in good agreement with the distributions of backbone conformations deduced from nmr coupling constants and the frequencies of glycosyl conformations in x-ray crystal structures, suggesting that the low energy states are reasonable. The low energy dimer forms (consisting of 150-325 conformational states per dimer step) are next used as variables in a Monte Carlo algorithm, which generates the conformations of single-stranded d(CXnG) chains, where X = A, T and n = 3, 4, 5. The oligonucleotides are built sequentially from the 5' end of the chain using random numbers to select the conformations of overlapping dimer units. The simulations are very fast, involving a total of 10(6) conformations per chain sequence. The potential errors in the buildup procedure are minimized by taking advantage of known rotational interdependences in the sugar-phosphate backbone. The distributions of oligonucleotide conformations are examined in terms of the magnitudes, positions, and orientations of the end-to-end vectors of the chains. The differences in overall flexibility and extension of the oligomers are discussed in terms of the conformations of the constituent dinucleotide steps, while the general methodology is discussed and compared with other nucleic acid model building techniques.

An ellipticine derivative (oxazolopyridocarbazolium) 3' linked to tetrathymidylate stacks intramolecularly with the nearest thymine at low concentration and head-to-tail intermolecularly at high concentration

The solution conformation of a tetrathymidylate linked through an ester bond to an ellipticine derivative oxazolopyridocarbazolium (OPC) at the 3' position was investigated using one- and two-dimensional nmr experiments. Since the total electric charge of the OPC ring may influence self-association, we first determined the pKa of the oxazole cyclic acidic function. Nuclear Overhauser effect spectroscopy experiments showed that, at low concentration, the OPC stacks intramolecularly with the nearest thymine at the 3' end. At highest concentration, however, the OPC rings are self-associated. The stacking constant was calculated using 1H chemical shift dilution experiment. The conformational model suggested by P-nmr was tested by molecular mechanics computations.