Flexibility of short-strand RNA in aqueous solution as revealed by molecular dynamics simulation:are A′-RNA and A-RNA distinct conformational structures?

We use molecular dynamics simulations to compare the conformational structure and dynamics of a 21-base pair RNA sequence initially constructed according to the canonical A-RNA and A'-RNA forms in the presence of counterions and explicit water. Our study aims to add a dynamical perspective to the solid-state structural information that has been derived from X-ray data for these two characteristic forms of RNA. Analysis of the three main structural descriptors commonly used to differentiate between the two forms of RNA namely major groove width, inclination and the number of base pairs in a helical twist over a 30 ns simulation period reveals a flexible structure in aqueous solution with fluctuations in the values of these structural parameters encompassing the range between the two crystal forms and more. This provides evidence to suggest that the identification of distinct A-RNA and A'-RNA structures, while relevant in the crystalline form, may not be generally relevant in the context of RNA in the aqueous phase. The apparent structural flexibility observed in our simulations is likely to bear ramifications for the interactions of RNA with biological molecules (e.g. proteins) and non-biological molecules (e.g. non-viral gene delivery vectors). © CSIRO 2009.

Phenotypic effect of mutations in evolving populations of RNA molecules

Background. The secondary structure of folded RNA sequences is a good model to map phenotype onto genotype, as represented by the RNA sequence. Computational studies of the evolution of ensembles of RNA molecules towards target secondary structures yield valuable clues to the mechanisms behind adaptation of complex populations. The relationship between the space of sequences and structures, the organization of RNA ensembles at mutation-selection equilibrium, the time of adaptation as a function of the population parameters, the presence of collective effects in quasispecies, or the optimal mutation rates to promote adaptation all are issues that can be explored within this framework. Results. We investigate the effect of microscopic mutations on the phenotype of RNA molecules during their in silico evolution and adaptation. We calculate the distribution of the effects of mutations on fitness, the relative fractions of beneficial and deleterious mutations and the corresponding selection coefficients for populations evolving under different mutation rates. Three different situations are explored: the mutation-selection equilibrium (optimized population) in three different fitness landscapes, the dynamics during adaptation towards a goal structure (adapting population), and the behavior under periodic population bottlenecks (perturbed population). Conclusions. The ratio between the number of beneficial and deleterious mutations experienced by a population of RNA sequences increases with the value of the mutation rate µ at which evolution proceeds. In contrast, the selective value of mutations remains almost constant, independent of µ, indicating that adaptation occurs through an increase in the amount of beneficial mutations, with little variations in the average effect they have on fitness. Statistical analyses of the distribution of fitness effects reveal that small effects, either beneficial or deleterious, are well described by a Pareto distribution. These results are robust under changes in the fitness landscape, remarkably when, in addition to selecting a target secondary structure, specific subsequences or low-energy folds are required. A population perturbed by bottlenecks behaves similarly to an adapting population, struggling to return to the optimized state. Whether it can survive in the long run or whether it goes extinct depends critically on the length of the time interval between bottlenecks. © 2010 Stich et al; licensee BioMed Central Ltd.

Forced retroevolution of an RNA bacteriophage

The operator hairpin ahead of the replicase gene in RNA bacteriophage MS2 contains overlapping signals for binding the coat protein and ribosomes. Coat protein binding inhibits further translation of the gene and forms the first step in capsid formation. The hairpin sequence was partially randomized to assess the importance of this structure element for the bacteriophage and to monitor alternative solutions that would evolve on the passaging of mutant phages. The evolutionary reconstruction of the operator failed in the majority of mutants. Instead, a poor imitation developed containing only some of the recognition signals for the coat protein. Three mutants were of particular interest in that they contained double nonsense codons in the lysis reading frame that runs through the operator hairpin. The simultaneous reversion of two stop codons into sense codons has a very low probability of occurring. Therefore the phage solved the problem by deleting the nonsense signals and, in fact, the complete operator, except for the initiation codon of the replicase gene. Several revertants were isolated with activities ranging from 1% to 20% of wild type. The operator, long thought to be a critical regulator, now appears to be a dispensable element. In addition, the results indicate how RNA viruses can be forced to step back to an attenuated form.

Long-range translational coupling in single-stranded RNA bacteriophages: an evolutionary analysis

In coliphage MS2 RNA a long-distance interaction (LDI) between an internal segment of the upstream coat gene and the start region of the replicase gene prevents initiation of replicase synthesis in the absence of coat gene translation. Elongating ribosomes break up the repressor LDI and thus activate the hidden initiation site. Expression studies on partial MS2 cDNA clones identified base pairing between 1427-1433 and 1738-1744, the so-called Min Jou (MJ) interaction, as the molecular basis for the long-range coupling mechanism. Here, we examine the biological significance of this interaction for the control of replicase gene translation. The LDI was disrupted by mutations in the 3'-side and the evolutionary adaptation was monitored upon phage passaging. Two categories of pseudorevertants emerged. The first type had restored the MJ interaction but not necessarily the native sequence. The pseudorevertants of the second type acquired a compensatory substitution some 80 nt downstream of the MJ interaction that stabilizes an adjacent LDI. In one examined case we confirmed that the second site mutations had restored coat-replicase translational coupling. Our results show the importance of translational control for fitness of the phage. They also reveal that the structure that buries the replicase start extends to structure elements bordering the MJ interaction.

LacI-mediated sequence-specific affinity purification of plasmid DNA for therapeutic applications

Affinity purification of plasmid DNA is an attractive option for the biomanufacture of therapeutic plasmids, which are strictly controlled for levels of host protein, DNA, RNA, and endotoxin. Plasmid vectors are considered to be a safer alternative than viruses for gene therapy, but milligram quantities of DNA are required per dose. Previous affinity approaches have involved triplex DNA formation and a sequence-specific zinc finger protein. We present a more generically applicable protein-based approach, which exploits the lac operator, present in a wide diversity of plasmids, as a target sequence. We used a GFP/His-tagged Lacl protein, which is precomplexed with the plasmid, and the resulting complex was immobilized on a solid support (TALON resin). Ensuing elution gives plasmid DNA, in good yield (>80% based on recovered starting material, 35-50% overall process), free from detectable RNA and protein and with minimal genomic DNA contamination. Such an affinity-based process should enhance plasmid purity and ultimately, after appropriate development, may simplify the biomanufacturing process of therapeutic plasmids.

Evaluation of DNA enzymes targeted against the RNA component of human telomerase

Glioblastoma Multiforme (GBM) is a highly malignant form of brain cancer for which there is currently no effective cure. Consequently, developing new therapies and elucidating effective targets is crucial for this fatal disease. In recent years, DNA enzymes, deoxyribonucleic acid molecules with enzymatic activity, have emerged. In the same manner as ribozymes, DNA enzymes are able to effect cleavage of RNA in a sequence-specific manner, and operate with catalytic efficiency. In this study, two DNA enzymes were designed to target the template region of human telomerase RNA (hTR), utilising the 10-23 and 8-17 catalytic motifs elucidated by Santoro and Joyce (1997). Telomerase is an RNA-dependent DNA polymerase, which stabilises telomere lengths by adding hexameric repeats (TTAGGG in humans) to chromosome termini, thus preventing the telomere shortening that usually occurs during mitotic cell division. Telomerase activity, whilst absent in normal somatic tissues, is present in almost 90% of all tumours. Thus, there is speculation that telomerase may be the much sought universal target for therapeutic intervention in cancer. In vitro cleavage assays showed both DNA enzymes to be catalytically competent. Unmodified phosphodiester (PO) backbone DNA enzymes were rapidly degraded in the presence of serum, with a half-life of 10 minutes. The common approach of introducing phosphorothioate (PS) linkages was used in an effort to overcome this instability. As a result of concurrent activity and stability studies on the DNA enzymes with various numbers of PS linkages, the DNA enzymes with a PO core and PS arms were chosen for use in further cell work. The cleavage activity of both was shown to be specific and affected by temperature, pH, MgCI2 concentration and enzyme concentration. Both DNA enzyme motifs reduced telomerase activity in cell lysates, as assessed by the telomerase repeat amplification protocol (TRAP) with an IC50 of 100nM. DNA enzymes being polyanionic molecules do not readily cross biological barriers. Cellular association of naked DNA enzyme was inefficient at less than 2%. Cellular delivery of the DNA enzymes was effectively improved using commercial cationic lipid formulations. However, the lipid-mediated delivery of DNA enzymes to U87-MG cells over a 4-hour period did not significantly inhibit cell proliferation compared to controls. This is possibly due to an expected lag period between the inhibition of telomere maintenance and cell death. Therefore, biodegradable polymer microspheres were investigated as a potential delivery option for prolonged and sustained delivery. In vitro release profiles showed that after an initial burst, sustained release of DNA enzymes was observed over 35 days. Finally, the efficacy and specificity of the DNA enzymes were demonstrated in a luciferase based reporter assay. Specific inhibition of luciferase expression was displayed at 10nM. Thus DNA enzymes have potential against endogenous cellular targets.

On the structural repertoire of pools of short, random RNA sequences

A detailed knowledge of the mapping between sequence and structure spaces in populations of RNA molecules is essential to better understand their present-day functional properties, to envisage a plausible early evolution of RNA in a prebiotic chemical environment and to improve the design of in vitro evolution experiments, among others. Analysis of natural RNAs, as well as in vitro and computational studies, show that certain RNA structural motifs are much more abundant than others, pointing out a complex relation between sequence and structure. Within this framework, we have investigated computationally the structural properties of a large pool (10 molecules) of single-stranded, 35 nt-long, random RNA sequences. The secondary structures obtained are ranked and classified into structure families. The number of structures in main families is analytically calculated and compared with the numerical results. This permits a quantification of the fraction of structure space covered by a large pool of sequences. We further show that the number of structural motifs and their frequency is highly unbalanced with respect to the nucleotide composition: simple structures such as stem-loops and hairpins arise from sequences depleted in G, while more complex structures require an enrichment of G. In general, we observe a strong correlation between subfamilies-characterized by a fixed number of paired nucleotides-and nucleotide composition. Our results are compared to the structural repertoire obtained in a second pool where isolated base pairs are prohibited. © 2008 Elsevier Ltd. All rights reserved.

Efficient HIV-1 inhibition by a 16 nt-long RNA aptamer designed by combining in vitro selection and in silico optimisation strategies

The human immunodeficiency virus type-1 (HIV-1) genome contains multiple, highly conserved structural RNA domains that play key roles in essential viral processes. Interference with the function of these RNA domains either by disrupting their structures or by blocking their interaction with viral or cellular factors may seriously compromise HIV-1 viability. RNA aptamers are amongst the most promising synthetic molecules able to interact with structural domains of viral genomes. However, aptamer shortening up to their minimal active domain is usually necessary for scaling up production, what requires very time-consuming, trial-and-error approaches. Here we report on the in vitro selection of 64 nt-long specific aptamers against the complete 5' -untranslated region of HIV-1 genome, which inhibit more than 75% of HIV-1 production in a human cell line. The analysis of the selected sequences and structures allowed for the identification of a highly conserved 16 nt-long stem-loop motif containing a common 8 nt-long apical loop. Based on this result, an in silico designed 16 nt-long RNA aptamer, termed RNApt16, was synthesized, with sequence 5'-CCCCGGCAAGGAGGGG-3-'. The HIV-1 inhibition efficiency of such an aptamer was close to 85%, thus constituting the shortest RNA molecule so far described that efficiently interferes with HIV-1 replication.