The base sequence dependent flexibility of linear single-stranded oligoribonucleotides correlates with the reactivity of the phosphodiester bond

Engineering supramolecular dynamics of self-assembly and turnover of oncogenic microRNAs to drive their synergistic destruction in tumor models

Rationally-engineered functional biomaterials offer the opportunity to interface with complex biology in a predictive, precise, yet dynamic way to reprogram their behaviour and correct shortcomings. Success here may lead to a desired therapeutic effect against life-threatening diseases, such as cancer. Here, we engineered "Crab"-like artificial ribonucleases through coupling of peptide and nucleic acid building blocks, capable of operating alongside and synergistically with intracellular enzymes (RNase H and AGO2) for potent destruction of oncogenic microRNAs. "Crab"-like configuration of two catalytic peptides ("pincers") flanking the recognition oligonucleotide was instrumental here in providing increased catalytic turnover, leading to ≈30-fold decrease in miRNA half-life as compared with that for "single-pincer" conjugates. Dynamic modeling of miRNA cleavage illustrated how such design enabled "Crabs" to drive catalytic turnover through simultaneous attacks at different locations of the RNA-DNA heteroduplex, presumably by producing smaller cleavage products and by providing toeholds for competitive displacement by intact miRNA strands. miRNA cleavage at the 5'-site, spreading further into double-stranded region, likely provided a synergy for RNase H1 through demolition of its loading region, thus facilitating enzyme turnover. Such synergy was critical for sustaining persistent disposal of continually-emerging oncogenic miRNAs. A single exposure to the best structural variant (Crab-p-21) prior to transplantation into mice suppressed their malignant properties and reduced primary tumor volume (by 85 %) in MCF-7 murine xenograft models.

Led-Seq: ligation-enhanced double-end sequence-based structure analysis of RNA

Structural analysis of RNA is an important and versatile tool to investigate the function of this type of molecules in the cell as well as in vitro. Several robust and reliable procedures are available, relying on chemical modification inducing RT stops or nucleotide misincorporations during reverse transcription. Others are based on cleavage reactions and RT stop signals. However, these methods address only one side of the RT stop or misincorporation position. Here, we describe Led-Seq, a new approach based on lead-induced cleavage of unpaired RNA positions, where both resulting cleavage products are investigated. The RNA fragments carrying 2', 3'-cyclic phosphate or 5'-OH ends are selectively ligated to oligonucleotide adapters by specific RNA ligases. In a deep sequencing analysis, the cleavage sites are identified as ligation positions, avoiding possible false positive signals based on premature RT stops. With a benchmark set of transcripts in Escherichia coli, we show that Led-Seq is an improved and reliable approach based on metal ion-induced phosphodiester hydrolysis to investigate RNA structures in vivo.

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.

Photocycle of the LOV-STAS protein from the pathogen Listeria monocytogenes

Listeria monocytogenes, a food-borne bacterial pathogen causing significant human mortality, propagates by expressing genes in response to environmental signals, such as temperature and pH. Listeria gene (lmo0799) encodes a protein homologous to the Bacillus subtilis YtvA, which has a flavin-light, oxygen or voltage (LOV) domain and a Sulfate Transporters Anti-Sigma factor antagonist (STAS) output domain that regulates transcription-initiation factor Sigma B in the bacterial stress response upon exposure to light. This could be significant for the pathogenesis of listeriosis because Sigma B has been linked to virulence of Listeria, and the Listeria Lmo0799 protein has recently been identified as a virulence factor activated by blue light. We have cloned, expressed heterologously in Escherichia coli and purified the full-length LM-LOV-STAS protein. Although it exhibits photochemical activity similar to that of YtvA, LM-LOV-STAS lacks an almost universally conserved arginine in the flavin-binding site, as well as another positively charged residue, a lysine in YtvA. The absence of these positive charges was found to destabilize retention of the flavin mononucleotide (FMN) chromophore in the LM-LOV-STAS protein, particularly at higher temperatures. The unusual sequence of the LM-LOV-STAS protein alters both spectral features and activation/deactivation kinetics, potentially expanding the sensory capacity of this LOV domain, e.g. to detect light plus cold.

Unraveling the structural complexity in a single-stranded RNA tail: implications for efficient ligand binding in the prequeuosine riboswitch

Single-stranded RNAs (ssRNAs) are ubiquitous RNA elements that serve diverse functional roles. Much of our understanding of ssRNA conformational behavior is limited to structures in which ssRNA directly engages in tertiary interactions or is recognized by proteins. Little is known about the structural and dynamic behavior of free ssRNAs at atomic resolution. Here, we report the collaborative application of nuclear magnetic resonance (NMR) and replica exchange molecular dynamics (REMD) simulations to characterize the 12 nt ssRNA tail derived from the prequeuosine riboswitch. NMR carbon spin relaxation data and residual dipolar coupling measurements reveal a flexible yet stacked core adopting an A-form-like conformation, with the level of order decreasing toward the terminal ends. An A-to-C mutation within the polyadenine tract alters the observed dynamics consistent with the introduction of a dynamic kink. Pre-ordering of the tail may increase the efficacy of ligand binding above that achieved by a random-coil ssRNA. The REMD simulations recapitulate important trends in the NMR data, but suggest more internal motions than inferred from the NMR analysis. Our study unmasks a previously unappreciated level of complexity in ssRNA, which we believe will also serve as an excellent model system for testing and developing computational force fields.

PNAzymes that are artificial RNA restriction enzymes

DNA-cleaving restriction enzymes are well-known tools in biomedical and biotechnological research. There are, however, no corresponding enzymes known for RNA cleavage. There has been an ongoing development of artificial ribonucleases, including some attempts at sequence selectivity. However, so far these systems have displayed modest rates of cleavage, and in most cases, the cleaver has been used in excess or in stoichiometric amounts. In the current work, we present PNA-based systems (PNAzymes) that carry a Cu(II)-2,9-dimethylphenanthroline group and that act as site and sequence specific RNases. The general basis for the systems is that the target is cleaved at a nonbase paired region (RNA bulge) which is formed in the substrate upon binding of the PNAzyme. With this copper based system, cleavage takes place at virtually only one site and with a half-life of down to 30 min under stoichiometric conditions. Efficient turnover of RNA-substrate is shown with a 100-fold excess of substrate, thus, demonstrating true enzyme behavior. In addition, alteration of the sequence in the RNA bulge or a mismatch in the base-pairing region leads to substantial decreases in rate showing both kinetic resolution and binding discrimination in the substrate selectivity. The selectivity is further demonstrated by the substrates, with two potential cleavage sites differing in only one base, are cleaved only at the site that either does not have a mismatch or is kinetically preferred. We suggest that these systems can serve as a basis for construction of RNA restriction enzymes for in vitro manipulations.

[An RNA sequence determines the speed of its splitting by artificial ribonucleases]

Phosphodiester bonds in RNA situated between similar nucleotides but in different sequences (context) were split under the action of artificial and natural ribonucleases with different speeds, and the reason for this phenomenon has not yet been fully revealed. In this study, the influence of one-nucleotide substitution on the sensitivity to splitting of the phosphodiester bonds in linear and structured RNA with homologous sequences is studied for the first time. It is indicated that the introduction of one-nucleotide substitution in the RNA sequence significantly (up to 10 times) changes the speed of the splitting of the bonds that are separated from the substitution point not only by 1-3, but also 6-8 nucleotides, by artificial ribonucleases. The observed regularities may be explained by the fact that the introduction of a one-nucleotide substitution significantly changes the stacking interactions and the net of hydrogen bonds in the RNA molecule. The applied value of this study consists of the ability of using low-molecular artificial ribonucleases with the aim of choosing the region of the binding of the oligonucleotide in the construction of a conjugate for the site-directed cutting of RNA, because the choice of a phosphodiester bond (motif) easily subjected to splitting significantly determines the effectiveness of artificial ribonucleases of directed action.

Intra- and intermolecular interactions influence the reactivity of RNA oligonucleotides

The transesterification of RNA oligonucleotides was studied over a wide pH range. The rate constants obtained indicate that, under neutral conditions, oligonucleotides with an adenosine moiety as the 5'-linked nucleoside can be up to 1000-fold more reactive than the reference oligonucleotide, a 13-mer oligo-U (1). Experiments with the modified oligonucleotide with N6,N6-dimethyladenosine (9) as the 5'-linked nucleoside moiety suggest that the exocyclic amino group is involved in the reaction, influencing the reactivity of the neighboring phosphodiester bond. In addition to such intramolecular interactions, weak intermolecular interactions most probably contribute to the reactivity.

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.

Molecular dynamics simulations of RNA: an in silico single molecule approach

RNA molecules are now known to be involved in the processing of genetic information at all levels, taking on a wide variety of central roles in the cell. Understanding how RNA molecules carry out their biological functions will require an understanding of structure and dynamics at the atomistic level, which can be significantly improved by combining computational simulation with experiment. This review provides a critical survey of the state of molecular dynamics (MD) simulations of RNA, including a discussion of important current limitations of the technique and examples of its successful application. Several types of simulations are discussed in detail, including those of structured RNA molecules and their interactions with the surrounding solvent and ions, catalytic RNAs, and RNA-small molecule and RNA-protein complexes. Increased cooperation between theorists and experimentalists will allow expanded judicious use of MD simulations to complement conceptually related single molecule experiments. Such cooperation will open the door to a fundamental understanding of the structure-function relationships in diverse and complex RNA molecules. .

Non-identical electronic characters of the internucleotidic phosphates in RNA modulate the chemical reactivity of the phosphodiester bonds

We here show that the electronic properties and the chemical reactivities of the internucleotidic phosphates in the heptameric ssRNAs are dissimilar in a sequence-specific manner because of their non-identical microenvironments, in contrast with the corresponding isosequential ssDNAs. This has been evidenced by monitoring the delta H8(G) shifts upon pH-dependent ionization (pK(a1)) of the central 9-guaninyl (G) to the 9-guanylate ion (G-), and its electrostatic effect on each of the internucleotidic phosphate anions, as measured from the resultant delta 31P shifts (pKa2) in the isosequential heptameric ssRNAs vis-à-vis ssDNAs: [d/r(5'-Cp1Ap2Q1p3Gp4Q2p5Ap6C-3'): Q1 = Q2 = A (5a/5b) or C (8a/8b), Q(1) = A, Q(2) = C (6a/6b), Q1 = C, Q2 = A (7a/7b)]. These oligos with single ionizable G in the centre are chosen because of the fact that the pseudoaromatic character of G can be easily modulated in a pH-dependent manner by its transformation to G- (the 2'-OH to 2-O- ionization effect is not detectable below pH 11.6 as evident from the N(1-Me)-G analog), thereby modulating/titrating the nature of the electrostatic interactions of G to G- with the phosphates, which therefore constitute simple models to interrogate how the variable pseudoaromatic characters of nucleobases under different sequence context (J. Am. Chem. Soc., 2004, 126, 8674-8681) can actually influence the reactivity of the internucleotide phosphates as a result of modulation of sequence context-specific electrostatic interactions. In order to better understand the impact of the electrostatic effect of the G to G- on the tunability of the electronic character of internucleotidic phosphates in the heptameric ssRNAs 5b, 6b, 7b and 8b, we have also performed their alkaline hydrolysis at pH 12.5 at 20 degrees C, and have identified the preferences of the cleavage sites at various phosphates, which are p2, p3 and p4 (Fig.3). The results of these alkaline hydrolysis studies have been compared with the hydrolysis of analogous N(1-Me)-G heptameric ssRNA sequences 5c, 7c and 8c under identical conditions in order to establish the role of the electrostatic effect of the 9-guanylate ion (and the 2'-OH to 2-O- ionization) on the internucleotidic phosphate. It turned out that the relative alkaline hydrolysis rate at those particular phosphates (p2, p3 and p4) in the N(1-Me)-G heptamers was reduced from 16-78% compared to those in the native counterparts [Fig. 4, and ESI 2 (Fig. S11)]. Thus, these physico-chemical studies have shown that those p2, p3 and p4 phosphates in the native heptameric RNAs, which show pKa2 as well as more deshielding (owing to weaker 31P screening) in the alkaline pH compared to those at the neutral pH, are more prone to the alkaline hydrolysis because of their relatively enhanced electrophilic character resulting from weaker 31P screening. This screening effect originates as a result of the systematic charge repulsion effect between the electron cloud in the outermost orbitals of phosphorus and the central guanylate ion, leading to delocalization of the phosphorus p(pi) charge into its dpi orbitals. It is thus likely that, just as in the non-enzymatic hydrolysis, the enzymatic hydrolysis of a specific phosphate in RNA by general base-catalysis in RNA-cleaving proteins (RNase A, RNA phosphodiesterase or nuclease) can potentially be electrostatically influenced by tuning the transient charge on the nucleobase in the steric proximity or as a result of specific sequence context owing to nearest-neighbor interactions.

Chemical models for ribozyme action

Mechanistic studies of the action of catalytic ribonucleic acids, ribozymes, are highly challenging, because even a slight structural change can dramatically affect the chain folding. This, in turn, alters the binding properties of the catalytic core, making identification of the real origin of the observed influence on rate difficult. Unambiguous structure-reactivity correlations based on studies with structurally simplified chemical models may help to distinguish between alternative mechanistic interpretations. The results of such model studies are reviewed. The topics include intramolecular cleavage of RNA phosphodiester bonds by solvent-derived species, general acids/bases and metal ions, effect of molecular environment on their hydrolytic stability and trinucleoside monophosphates as models for large ribozymes.

Artificial ribonucleases

Mimicking the action of enzymes by simpler and more robust man-made catalysts has long inspired bioorganic chemists. During the past decade, mimics for RNA-cleaving enzymes, ribonucleases, or, more precisely, mimics of ribozymes that cleave RNA in sequence-selective rather than base-selective manner, have received special attention. These artificial ribonucleases are typically oligonucleotides (or their structural analogs) that bear a catalytically active conjugate group and catalyze sequence-selective hydrolysis of RNA phosphodiester bonds.

Enhanced RNA cleavage within bulge-loops by an artificial ribonuclease

Cleavage of phosphodiester bonds by small ribonuclease mimics within different bulge-loops of RNA was investigated. Bulge-loops of different size (1–7 nt) and sequence composition were formed in a 3′ terminal fragment of influenza virus M2 RNA (96 nt) by hybridization of complementary oligodeoxynucleotides. Small bulges (up to 4 nt) were readily formed upon oligonucleotide hybridization, whereas hybridization of the RNA to the oligonucleotides designed to produce larger bulges resulted in formation of several alternative structures. A synthetic ribonuclease mimic displaying Pyr–Pu cleavage specificity cleaved CpA motifs located within bulges faster than similar motifs within the rest of the RNA. In the presence of 10 mM MgCl₂, 75% of the cleavage products resulted from the attack of this motif. Thus, selective RNA cleavage at a single target phosphodiester bond was achieved by using bulge forming oligonucleotides and a small ribonuclease A mimic.

Synthetic metallonucleases for RNA cleavage

Synthetic metallonucleases are versatile metal ion catalysts that use multiple catalytic strategies for the cleavage of RNA. Recent work in the design of more active metallonucleases combines a single metal ion with functional groups that interact with RNA, including amino acid fragments or additional metal ions. Rate enhancements by multifunctional catalysts for cleavage of simple model substrates with good leaving groups are as high as 10(6) but somewhat lower (10(5)) for real RNA. However, cleavage of RNA substrates is complicated by different binding modes and steric interactions that can interfere with catalysis. Antisense oligonucleotides, peptides and small molecules that act as RNA recognition agents increase the strength of substrate binding, but not necessarily the catalytic rate constant. In general, catalytic strategies used by synthetic metallonucleases are probably not optimized. A better grasp of the mechanism of RNA cleavage by metal ions and more effort on positioning the metal ion complex with respect to the cleavage site may lead to improved catalysts.

Stabilisation of the transition state of phosphodiester bond cleavage within linear single-stranded oligoribonucleotides

The effect of base sequence on the stability of the transition state (TS) of phosphodiester bond cleavage within linear single-stranded oligoribonucleotides has been studied in order to better understand why the reactivity of some phosphodiester bonds is enhanced compared to an unconstrained linkage. Molecular dynamics simulations of 3.0 ns were carried out for 14 oligonucleotides that contain in the place of the scissile phosphodiester bond a phosphorane structure mimicking the TS of the bond cleavage. The hydrolytic stability of the same oligonucleotides had previously been reported. Both the non-bridging oxyanions and the leaving 5[prime or minute]-oxygen of the pentacoordinated phosphorane moiety were observed to form hydrogen bonds with solvent water molecules in a similar way with all the compounds studied. In addition, water mediated hydrogen bonds between the phosphorane non-bridging oxyanions and the bases of the 3[prime or minute]-flanking sequence were detected with some of the compounds, but not with the most labile ones. Hence, it seems that the enhanced cleavage of some internucleosidic linkages does not result from the TS stabilisation by hydrogen bonding. With heterooligomers, the stacking of bases next to the cleavage site was observed to be enhanced on going from the initial state to the TS, whereas within uracil homooligomer, having initially negligible stacking, no change in the magnitude of stacking was seen. Accordingly, while strong stacking in the initial state is known to retard the phosphodiester bond cleavage, it may in the TS accelerate the reaction. Therefore, enhanced stacking on going from the initial state to transition state appears to be a factor that markedly contributes to the hydrolytic stability of phosphodiester bonds within oligonucleotides and may, at least partly, explain accelerated cleavage compared to fully unconstrained bonds, such as those in polyuridylic acid.