Influence of adjacent bases on the stacking-unstacking process of single-stranded oligonucleotides

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.

Structural modifications as tools in mechanistic studies of the cleavage of RNA phosphodiester linkages

The cleavage of RNA phosphodiester bonds by RNase A and hammerhead ribozyme at neutral pH fundamentally differs from the spontaneous reactions of these bonds under the same conditions. While the predominant spontaneous reaction is isomerization of the 3',5'-phosphodiester linkages to their 2',5'-counterparts, this reaction has never been reported to compete with the enzymatic cleavage reaction, not even as a minor side reaction. Comparative kinetic measurements with structurally modified di-nucleoside monophosphates and oligomeric phosphodiesters have played an important role in clarification of mechanistic details of the buffer-independent and buffer-catalyzed reactions. More recently, heavy atom isotope effects and theoretical calculations have refined the picture. The primary aim of all these studies has been to form a solid basis for mechanistic analyses of the action of more complicated catalytic machineries. In other words, to contribute to conception of a plausible unified picture of RNA cleavage by biocatalysts, such as RNAse A, hammerhead ribozyme and DNAzymes. In addition, structurally modified trinucleoside monophosphates as transition state models for Group I and II introns have clarified some features of the action of large ribozymes.

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.

Intramolecular base stacking of dinucleoside monophosphate anions in aqueous solution

Time-dependent motions of 32 deoxyribodinucleoside and ribodinucleoside monophosphate anions in aqueous solution at 310 K were monitored during 40 ns using classical molecular dynamics (MD). In all studied molecules, spontaneous stacking/unstacking transitions occurred on a time-scale of 10 ns. To facilitate the structural analysis of the sampled configurations we defined a reaction coordinate for the nucleobase stacking that considers both the angle between the planes of the two nucleobases and the distance between their mass-centers. Additionally, we proposed a physically meaningful transient point on this coordinate that separates the stacked and unstacked states. We applied this definition to calculate free energies for stacking of all pairwise combinations of adenine, thymine (uracil), cytosine and guanine moieties embedded in studied dinucleosides monophosphate anions. The stacking equilibrium constants decreased in the order 5'-AG-3' > GA ~ GG ~ AA > GT ~ TG ~ AT ~ GC ~ AC > CG ~ TA > CA ~ TC ~ TT ~ CT ~ CC. The stacked conformations of AG occurred 10 times more frequently than its unstacked conformations. On the other hand, the last five base combinations showed a greater preference for the unstacked than the stacked state. The presence of an additional 2'-OH group in the RNA-based dinucleoside monophosphates increased the fraction of stacked complexes but decreased the compactness of the stacked state. The calculated MD trajectories were also used to reveal prevailing mutual orientation of the nucleobase dipoles in the stacked state.

Nonenzymatic RNA ligation in water

We describe the nonenzymatic ligation of RNA oligomers in water. Dimers and tetramers are formed in a time-, pH-, and temperature-dependent reaction. Ligation efficiency depends on oligonucleotide length and sequence and is strongly enhanced by adenine-based nucleotide cofactors. Ligation of short RNA fragments could have liberated the prebiotic polymerization systems from the thermodynamically demanding task of reaching a (pre)genetically meaningful size by stepwise addition of one precursor monomer at the time.

Molecular complexity favors the evolution of ribopolymers

We have explored the stability of selected ribo oligomers in water and have determined the physical-chemical conditions in which the key 3'-phosphoester bond is more stable when embedded in the polymer than when present in the monomer. In these conditions, the spontaneous formation and the survival of ribo polymers are potentially favored. A narrow pH range was identified in which complex sequences resist degradation markedly more than monotonous ones, thus potentially favoring the evolution of sequence-based genetic information. Given that the founding property of a polymer is to maintain its polymeric form and its sequence information, these findings support the view that the evolution of pregenetic molecular information occurred based on intrinsic properties of nucleic polymers.

Predicting helical coaxial stacking in RNA multibranch loops

The hypothesis that RNA coaxial stacking can be predicted by free energy minimization using nearest-neighbor parameters is tested. The results show 58.2% positive predictive value (PPV) and 65.7% sensitivity for accuracy of the lowest free energy configuration compared with crystal structures. The probability of each stacking configuration can be predicted using a partition function calculation. Based on the dependence of accuracy on the calculated probability of the stacks, a probability threshold of 0.7 was chosen for predicting coaxial stacks. When scoring these likely stacks, the PPV was 66.7% at a sensitivity of 51.9%. It is observed that the coaxial stacks of helices that are not separated by unpaired nucleotides can be predicted with a significantly higher accuracy (74.0% PPV, 66.1% sensitivity) than the coaxial stacks mediated by noncanonical base pairs (55.9% PPV, 36.5% sensitivity). It is also shown that the prediction accuracy does not show any obvious trend with multibranch loop complexity as measured by three different parameters.

Molecular dynamics applied to nucleic acids

In this Account, we focus on molecular dynamics (MD) simulations involving fully solvated nucleic acids. Historically, MD simulations were first applied to proteins and several years later to nucleic acids. The first MD simulations of DNA were carried out in vacuo, but nowadays fully solvated systems are common practice. Recently, technical improvements have made it possible to conduct accurate MD simulations of highly charged nucleic acids. The state-of-the-art of MD simulations and a number of applications on various nucleic acid systems are discussed.

Influence of sodium ions on the dynamics and structure of single-stranded DNA oligomers: a molecular dynamics study

The effects of sodium counterion presence and chain length on the structure and dynamics of single DNA strands of polythymidylate were studied by means of molecular dynamics simulations. The importance of the base-base stacking phenomenon increases with the chain length and partially reduces the flexibility of the strand. Sodium ions directly interact with the phosphate groups and keto oxygens of the thymine bases, complexes showing lifetimes below 400 ps. Simultaneous phosphate and keto complexes were observed for one of the sodium ions with lifetimes around 1 ns. The implications of such complexes in the folding process experienced by the strand are considered. Structurally, cation inner- and outer-sphere complexes were observed in the coordination of phosphate groups. For the inner-sphere complexes, the structural information retrieved from the simulations is in very good agreement with experimental data. The diffusion properties of the sodium ions also reflect both types of coordination modes.

Thermodynamic parameters for DNA sequences with dangling ends

The thermodynamic contributions to duplex formation of all 32 possible single-nucleotide dangling ends on a Watson-Crick pair are reported. In most instances, dangling ends are stabilizing with free energy contributions ranging from +0.48 (GT(A)) to-0.96 kcal/mol (). In comparison, Watson-Crick nearest-neighbor increments range from -0. 58 (TA/AT) to -2.24 (GC/CG) kcal/mol. Hence, in some cases, a dangling end contributes as much to duplex stability as a Watson-Crick A-T base pair. The implications of these results for DNA probe design are discussed. Analysis of the sequence dependence of dangling-end stabilities show that the nature of the closing base pair largely determines the stabilization. For a given closing base pair, however, adenine dangling ends are always more or equally as stable as the other dangling nucleotides. Moreover, 5' dangling ends are more or equally as stabilizing as their 3' counterparts. Comparison of DNA with RNA dangling-end motifs shows that DNA motifs with 5' dangling ends contribute to stability equally or more than their RNA counterparts. Conversely, RNA 3' dangling ends contribute to stability equally or more than their DNA counterparts. This data set has been incorporated into a DNA secondary structure prediction algorithm (DNA MFOLD) (http://mfold2.wustl.edu/mfold/dna/for m1.cgi) as well as a DNA hybridization prediction algorithm (HYTHERtrade mark) (http://jsl1.chem.wayne.edu/Hyther/hythermenu .html).