Structural determinants of RNA function

Theoretical studies of an exceptionally stable RNA tetraloop: observation of convergence from an incorrect NMR structure to the correct one using unrestrained molecular dynamics

We report on the results of five independent and unrestrained molecular dynamics simulations of an RNA tetraloop, r(GGACUUCGGUCC), and its related structures with the loop UUCG sugars changed to deoxyribose. Two separate NMR structures have been reported for the loop portion of this molecule, with the second refinement resulting in a slightly different and more accurate conformation for the loop. The root-mean-square deviation (RMSd) between the two NMR structures, for the loop portions only, is 2.5 A. Our simulations, starting from the two NMR structures, demonstrate that this tetraloop is a very stable and rigid structure with both nanosecond length simulations staying very close to the initial structures. Additionally, both simulations preserved most, if not all, of the NMR-derived interactions and violated very few of the nuclear Overhauser effect (NOE)-derived distances used in the structure refinements. However, when the two NMR structures were simulated with deoxyriboses in the loops instead of the native riboses, the flexibility of the systems increased and we observed a conversion from the incorrect to the correct loop conformation in the simulation which started in the incorrect loop conformation. When the riboses were subsequently re-introduced back into the structure which underwent the conversion, the agreement between this simulation and the one starting from the correct NMR structure was a remarkably low 0.5 A, demonstrating an almost complete convergence from the incorrect to the correct structure using unrestrained molecular dynamics.

Determining the conformation of RNAs in solution. Application to a retroviral system: structure of the HIV-1 primer binding site region and effect of tRNA(3Lys) binding

RNAs play a crucial and central role in a large variety of biological functions obviously linked to the wide variety of structures that they can adopt. Understanding the function of RNAs thus requires the knowledge of their two- and three-dimensional structures. We describe in detail the way to access the secondary structure of RNAs, by combining sequence comparison, secondary structure prediction by computer and, mainly, experimental data obtained by probing with chemicals and ribonucleases. These approaches were used to investigate secondary structure of the region containing the primer binding site of HIV-1 genomic RNA either free or involved in the binary complex with the replication primer tRNA(3Lys).

A graph-topological approach to recognition of pattern and similarity in RNA secondary structures

Secondary and tertiary RNA structures play an important role in many biological processes. Therefore the necessity arises to find similar higher-order structures for different but functionally homologous RNA sequences. We propose here a graph-topological approach to the problem, which shows two main features: simplified graph representation which allows the recognition of similarity of RNA secondary structures with the same branching look despite minor differences. This allows comparison among foldings from different sequences, and "pruning" of the secondary structures not shared by all the sequences since the early stages of the search. (b) The graph representation is encoded by the Randić topological index, and the search for the folding similarity is reduced to checking the identity of single numbers. These characteristics make this approach significantly different, less depending on empirical criteria, and less computationally heavy then previous methods, where the folding consensus has been measured by an alignment procedure or correlation of strings representing the secondary structures. Some U2 snRNA and viroid sequences are studied by this approach, which is imbedded in our previous search method based on genetic algorithms.

Molecular recognition in the bovine immunodeficiency virus Tat peptide-TAR RNA complex

BACKGROUND

In lentiviruses such as human immunodeficiency virus (HIV) and bovine immunodeficiency virus (BIV), the Tat (trans-activating) protein enhances transcription of the viral RNA by complexing to the 5'-end of the transcribed mRNA, at a region known as TAR (the trans-activation response element). Identification of the determinants that account for specific molecular recognition requires a high resolution structure of the Tat peptide-TAR RNA complex.

RESULTS

We report here on the structural characterization of a complex of the recognition domains of BIV Tat and TAR in aqueous solution using a combination of NMR and molecular dynamics. The 17-mer Tat peptide recognition domain folds into a beta-hairpin and penetrates in an edge-on orientation deep into a widened major groove of the 28-mer TAR RNA recognition domain in the complex. The RNA fold is defined, in part, by two uracil bulged bases; U12 has a looped-out conformation that widens the major groove and U10 forms a U.AU base triple that buttresses the RNA helix. Together, these bulged bases induce a approximately 40 degree bend between the two helical stems of the TAR RNA in the complex. A set of specific intermolecular hydrogen bonds between arginine side chains and the major-groove edge of guanine residues contributes to sequence specificity. These peptide-RNA contacts are complemented by other intermolecular hydrogen bonds and intermolecular hydrophobic packing contacts involving glycine and isoleucine side chains.

CONCLUSIONS

We have identified a new structural motif for protein-RNA recognition, a beta-hairpin peptide that interacts with the RNA major groove. Specificity is associated with formation of a novel RNA structural motif, a U.AU base triple, which facilitates hydrogen bonding of an arginine residue to a guanine and to a backbone phosphate. These results should facilitate the design of inhibitors that can disrupt HIV Tat-TAR association.

A genetic algorithm to search for optimal and suboptimal RNA secondary structures

Genetic algorithms are a search method used in solving problems by selection, recombination and mutation of tentative solutions, until the better ones are achieved. They are very efficient when the 'building block' hypothesis is effective for the solutions, which means that a better solution can be obtained by assembling short 'motifs' or 'schemata' that can be retrieved in some other worse solutions. The additive nature of the secondary structure free energy rules suggests the validity of this hypothesis, and therefore the likely power of a genetic algorithm approach to search for RNA secondary structures. We describe in detail an original genetic algorithm specific for this problem. The sharing function used to obtain differentiated solutions is also described. It results in a greater effectiveness of the algorithm in retrieving a large number of suboptimal RNA foldings besides the optimal one. RNA sequences of different length are used to test the method. The PSTV viroid sequence has been studied.

Stabilizing effects of the RNA 2'-substituent: crystal structure of an oligodeoxynucleotide duplex containing 2'-O-methylated adenosines

BACKGROUND

The stability of hybrids of 2'-O-methyl-ribonucleotides with complementary RNA is considerably higher than that of the corresponding DNA.RNA duplexes. The 2'-O-modified ribonucleotides are thus an attractive class of compounds for antisense applications. Understanding how these substituents stabilize the structure of the hybrid duplex may be important in the design of ribonucleotides with novel properties.

RESULTS

The crystal structure of a dimer of the self-complementary DNA strand d(GCGT)O2'mer(A)d(TACGC), which has a 2'-O-methylated ribonucleotide incorporated at position 5, was determined at 2.1 A resolution. This strand forms a duplex with an overall A-type conformation; the methyl groups of the two modified adenosines point into the relatively wide minor groove. Both 2'-methoxy groups are hydrogen-bonded to solvent molecules. These results allowed us to build a model of a fully 2'-O-methylated RNA double helix.

CONCLUSIONS

Insertion of 2'-O-modified RNA residues into a stretch of DNA can nucleate a local A-type conformation, in part because modification with a bulky residue at this position stabilizes a C3'-endo type sugar pucker. The increased stability of fully 2'-O-methylated RNA may result from hydrophobic interactions between substituents in the minor groove. As the 2'-O-methyl groups are directed into the minor groove, it may be worthwhile to introduce tailor-made 2'-O-substituents into RNA; it might be possible to design groups that both stabilize the hybrid duplexes and carry a nuclease function, further improving the efficacy of these modified RNAs in antisense applications.