High-resolution NMR studies of A- and G-containing oligonucleotides

What is the conformation of physiologically-active dinucleoside polyphosphates in solution? Conformational analysis of free dinucleoside polyphosphates by NMR and molecular dynamics simulations

Dinucleoside polyphosphates, or dinucleotides (Np(n)N'; N, N' = A, U, G, C; n = 2-7), are naturally occurring ubiquitous physiologically active compounds. Despite the interest in dinucleotides, and the relevance of their conformation to their biological function, the conformation of dinucleotides has been insufficiently studied. Therefore, here we performed conformational analysis of a series of Np(n)N' Na(+) salts (N = A, G, U, C; N' = A, G, U, C; n = 2-5) by various NMR techniques. All studied dinucleotides, except for Up(4/5)U, formed intramolecular base stacking interactions in aqueous solutions as indicated by NMR. The conformation around the glycosidic angle in Np(n)N's was found to be anti/high anti and the preferred conformation around the C4'-C5', C5'-O5' bonds was found to be gauche-gauche (gg). The ribose moiety in Np(n)N's showed a small preference for the S conformation, but when attached to cytosine the ribose ring preferred the N conformation. However, no predominant conformation was observed for the ribose moiety in any of the dinucleotides. Molecular dynamics simulations of Ap(2)A and Ap(4)A Na(+) salts supported the experimental results. In addition, three modes of base-stacking were found for Ap(2/4)A: α-α, β-β and α-β, which exist in equilibrium, while none is dominant. We conclude that natural, free Np(n)N's (n = 2-5) at physiological pH exist mostly in a folded (stacked), rather than extended conformation, in several interconverting stacking modes. Intramolecular base stacking of Np(n)N's does not alter the conformation of each of the nucleotide moieties, which remains the same as that of the mononucleotides in solution.

Critical effect of the N2 amino group on structure, dynamics, and elasticity of DNA polypurine tracts

Unrestrained 5-20-ns explicit-solvent molecular dynamics simulations using the Cornell et al. force field have been carried out for d[GCG(N)11GCG]2 (N, purine base) considering guanine*cytosine (G*C), adenine*thymine (A*T), inosine*5-methyl-cytosine (I*mC), and 2-amino-adenine*thymine (D*T) basepairs. The simulations unambiguously show that the structure and elasticity of N-tracts is primarily determined by the presence of the amino group in the minor groove. Simulated A-, I-, and AI-tracts show almost identical structures, with high propeller twist and minor groove narrowing. G- and D-tracts have small propeller twisting and are partly shifted toward the A-form. The elastic properties also differ between the two groups. The sequence-dependent electrostatic component of base stacking seems to play a minor role. Our conclusions are entirely consistent with available experimental data. Nevertheless, the propeller twist and helical twist in the simulated A-tract appear to be underestimated compared to crystallographic studies. To obtain further insight into the possible force field deficiencies, additional multiple simulations have been made for d(A)10, systematically comparing four major force fields currently used in DNA simulations and utilizing B and A-DNA forms as the starting structure. This comparison shows that the conclusions of the present work are not influenced by the force field choice.

An energy model that predicts the correct folding of both the tRNA and the 5S RNA molecules

A new set of energy values to predict the secondary structures in RNA molecules has been derived through a multiple-step refinement procedure. It achieves more than 80% success in predicting the cloverleaf pattern in tRNA (200 sequences tested) and more than 60% success in predicting the consensus folding of 5S RNA (100 sequences). Improvements in our initial program for predicting secondary structures, based on the principle of the "incompatibility islets" made possible the work on 5S RNA. The program was speeded up by introducing a dynamic grouping of the islets into three disjoint blocks. The novel features in the energy model include i) an evaluation of the contribution of odd pairs according to their position within a segment ii) a penalty for internal loops related to their dissymmetry iii) a bonus for bulge loops when the two terminal paired bases at the junction point are both pyrimidines.