Monte Carlo simulations on short single-stranded oligonucleotides. I. Application to RNA trimers

Conformational flexibility of two RNA trimers explored by computational tools and database search

Two RNA sequences, AAA and AUG, were studied by the conformational search program CICADA and by molecular dynamics (MD) in the framework of the AMBER force field, and also via thorough PDB database search. CICADA was used to provide detailed information about conformers and conformational interconversions on the energy surfaces of the above molecules. Several conformational families were found for both sequences. Analysis of the results shows differences, especially between the energy of the single families, and also in flexibility and concerted conformational movement. Therefore, several MD trajectories (altogether 16 ns) were run to obtain more details about both the stability of conformers belonging to different conformational families and about the dynamics of the two systems. Results show that the trajectories strongly depend on the starting structure. When the MD start from the global minimum found by CICADA, they provide a stable run, while MD starting from another conformational family generates a trajectory where several different conformational families are visited. The results obtained by theoretical methods are compared with the thorough database search data. It is concluded that all except for the highest energy conformational families found in theoretical result also appear in experimental data. Registry numbers: adenylyl-(3' --> 5')-adenylyl-(3' --> 5')-adenosine [917-44-2] adenylyl-(3' --> 5')-uridylyl-(3' --> 5')-guanosine [3494-35-7].

Common structural features of UUCG and UACG tetraloops in very short hairpins determined by UV absorption, Raman, IR and NMR spectroscopies

Thermodynamic and structural properties of two UNCG tetraloops in very short hairpin octamers, 5'-r(GCUUCGGC)-3' and 5'-r(GCUACGGC)-3', have been studied by means of various physical techniques. Melting profiles of both octamers, obtained from UV absorption spectra taken as a function of temperature, are consistent with a monophasic, progressive and completely reversible order-to-disorder transition and confirm their unusual structural stability (Tm > 51 degrees C). The 1H, 13C and 31P NMR chemical shifts and coupling constants of the UACG loop nucleotides are comparable with those reported previously for UUCG loops, i.e. 2'-endo/anti conformation of the second and third nucleotide of the loop as well as the syn orientation of the ultimate guanine base and the A-type double helical conformation of the hairpin stem. Simulation of quantitative NOESY volumes shows that the UACG octamer adopts a very rigid compact structure which is well represented by an average order parameter of 0.9. Three base-pairs and four additional strong hydrogen bonds are undoubtedly responsible for such limited flexibility. Raman and infrared spectra as a function of temperature reflect the order-to-disorder transition, as well. Vibrational conformational markers in low temperature spectra of both octamers indicate the hairpin structure as the major conformer in aqueous phase. These spectra further support the structural features of most of the nucleotides involved in the tetraloops and clearly demonstrate the structural similarities of the phosphodiester backbone in both hairpins. Consequently, on the basis of all present results, one can deduce that the conformational features of the UUCG and UACG tetraloops seem to be inherent to the UNCG type tetraloops, regardless of either the nature of the tetraloop second base or the stem length.

Thermodynamic and structural properties of r(ACC) as revealed by ultraviolet electronic absorption, circular dichroism, 1H-NMR spectroscopy and Monte Carlo simulations

UV absorption, circular dichroism (CD) and 1H NMR, associated with Monte Carlo (MC) molecular structure simulations have been applied to the study of the trinucleoside diphosphate: r(ACC). The MC study which has been conducted as a function of temperature, is based on random variations of the nucleotide conformational angles, i.e. phosphodiester chain torsional angles and sugar pucker pseudorotational angles. All of the chemical bond lengths and valence angles remained fixed during the structural simulation, except those of the sugar pucker. Six different initial structures have been selected in order to explore the molecular conformational space as completely as possible. This simulation procedure led to distinct families of equilibrium conformations at 283, 298 and 318 K. The thermodynamical parameters such as variations in entropy, enthalpy and also melting temperature (delta SX0, delta HX0 and Tm) of the stacking (X) equilibrium were obtained from UV absorption and circular dichroism (CD) spectra recorded over a 80K temperature range. Chemical shifts (delta), vicinal coupling constants (3Jk,l), and cross-relaxation rate (sigma k,l) of trimers were measured at 400.13 MHz over a range of concentrations (2-13 mM) and temperatures (283-333K). Least-squares fitting of the experimental chemical shifts to simple models of association (A) and stacking equilibria allowed separation of the variations in the delta values (delta delta X and delta delta A) due to either phenomenon. The three NMR data sets (delta delta X, 3Jk,l, and sigma k,l) were then evaluated for the minima conformers obtained with the MC stimulations. Theoretical values of delta delta X were estimated using the results of an ab initio study while the coupling constant data were simulated with Karplus-type equations. Finally, the relaxation data were simulated from the distance matrices using treatment for cases of both slow conformational exchange accompanied by rapid small-amplitude fluctuations about the minima structures. A consistent picture of the large amplitude deformations (torsional angle variation) of these trimers has emerged from the present study. Optimized conformational blends at 283,296 and 318K were obtained by least-squares fitting of the experimental data to the theoretical ones, while considering the populations as adjustable parameters. As it would be expected, the right-handed helical conformation (A-RNA type) is found to be the major stacked species, in the temperature range of 283 to 318K. Limited evidence for bulged structures has been obtained, whereas novel reverse-stacked and half-stacked conformers also presented theoretical data compatible with the NMR observables of aqueous r(ACC).