Computational design of allosteric ribozymes as molecular biosensors

Nucleic acids have proven to be a very suitable medium for engineering various nanostructures and devices. While synthetic DNAs are commonly used for self-assembly of nanostructures and devices in vitro, functional RNAs, such as ribozymes, are employed both in vitro and in vivo. Allosteric ribozymes have applications in molecular computing, biosensoring, high-throughput screening arrays, exogenous control of gene expression, and others. They switch on and off their catalytic function as a result of a conformational change induced by ligand binding. Designer ribozymes are engineered to respond to different effectors by in vitro selection, rational and computational design methods. Here, I present diverse computational methods for designing allosteric ribozymes with various logic functions that sense oligonucleotides or small molecules. These methods yield the desired ribozyme sequences within minutes in contrast to the in vitro selection methods, which require weeks. Methods for synthesis and biochemical testing of ribozymes are also discussed.     

In vitro selection of allosteric ribozymes: theory and experimental validation

In vitro selection techniques offer powerful and versatile methods to isolate nucleic acid sequences with specific activities from huge libraries. We describe an in vitro selection strategy for the de novo selection of allosteric self-cleaving ribozymes responding to pefloxacin and other quinolone derivatives. Within 16 selection cycles, highly sensitive clones responding to drug levels in the sub-micromolar range were obtained. The morpholine moiety of the quinolone derivatives was required for inhibition of the self-cleavage of the selected ribozymes: modifications of the aromatic system were tolerated better than modifications of the morpholine ring. We also present a theoretical model that analyzes the predicted fraction of ribozymes with a given binding constant and cleavage rate recovered after each selection cycle. This model precisely predicts the actual experimental values obtained with the selection procedure. It can thus be used to determine the optimal conditions for an in vitro selection of an allosteric ribozyme with a desired dissociation constant and cleavage rate for a given application.     

Engineering high-speed allosteric hammerhead ribozymes

Full-length hammerhead ribozymes were subjected to in vitro selection to identify variants that are allosterically regulated by theophylline in the presence of a physiologically relevant concentration of Mg(2+). The population of allosteric ribozymes resulting from 15 rounds of in vitro selection yielded variants with observed rate constants (k (obs)) as high as 8 min(-1) in the presence of theophylline and maximal k (obs) increases of up to 285-fold compared to rate constants measured in the absence of effector. The selected ribozymes have kinetic characteristics that are predicted to be sufficient for cellular gene control applications, but do not exhibit any activity in reporter gene assays. The inability of the engineered RNAs to control gene expression suggests that the in vitro and in vivo folding pathways of the RNAs are different. These results provide several key pieces of information that will aid in future efforts to engineer allosteric ribozymes for gene control applications.     

Monitoring protein modification with allosteric ribozymes

An allosteric ribozyme is an RNA-based enzyme (ribozyme) whose catalytic activity is modulated by molecular recognition of a protein. The direct coupling of a detectable catalytic event to molecular recognition by an allosteric ribozyme enables simple assays for quantitative protein detection. Most significantly, the mode of development and molecular recognition characteristics of allosteric ribozymes are fundamentally different from antibodies, providing them with functional characteristics that complement those of antibodies. Allosteric ribozymes can be developed using native proteins and, therefore, are often sensitive to protein conformation. In contrast, antibodies tend to recognize a series of adjacent amino acids as a consequence of antigen presentation and typically are not sensitive to protein conformation. Unlike antibody development, the development of allosteric ribozymes is a completely in vitro process that allows the specificity of an allosteric ribozyme to be tightly controlled. These significant differences from antibodies allow the pre-programmed development of conformation-state-specific protein detection reagents that can be used to investigate the activation-state of signal transduction components.     

Engineered allosteric ribozymes as biosensor components

RNA and DNA molecules can be engineered to function as molecular switches that trigger catalytic events when a specific target molecule becomes bound. Recent studies on the underlying biochemical properties of these constructs indicate that a significant untapped potential exists for the practical application of allosteric nucleic acids. Engineered molecular switches can be used to report the presence of specific analytes in complex mixtures, making possible the creation of new types of biosensor devices and genetic control elements.     

Controlling the rate of organic reactions: rational design of allosteric Diels-Alderase ribozymes

Allosteric mechanisms are widely used in nature to control the rates of enzymatic reactions, but little is known about RNA catalysts controlled by these principles. The only natural allosteric ribozyme reported to date catalyzes an RNA cleavage reaction, and so do almost all artificial systems. RNA has, however, been shown to accelerate a much wider range of chemical reactions. Here we report that RNA catalysts for organic reactions can be put under the stringent control of effector molecules by straight-forward rational design. This approach uses known RNA sequences with catalytic and ligand-binding properties, and exploits weakly conserved sequence elements and available structural information to induce the formation of alternative, catalytically inactive structures. The potential and general applicability is demonstrated by the design of three different systems in which the rate of a catalytic carbon–carbon bond forming reaction is positively regulated up to 2100-fold by theophylline, tobramycin and a specific mRNA sequence, respectively. Although smaller in size than a tRNA, all three ribozymes show typical features of allosteric metabolic enzymes, namely high rate acceleration and tight allosteric regulation. Not only do these findings demonstrate RNA's power as a catalyst, but also highlight on RNA's capabilities as signaling components in regulatory networks.     

Monitoring post-translational modification of proteins with allosteric ribozymes

An allosteric hammerhead ribozyme activated specifically by the unphosphorylated form of the protein kinase ERK2 was created through a rational design strategy that relies on molecular recognition of ERK2 to decrease the formation of an alternate, inactive ribozyme conformer. Neither closely related mitogen-activated protein kinases (MAPKs) nor the phosphorylated form of ERK2 induced ribozyme activity. The ribozyme quantitatively detected ERK2 added to mammalian cell lysates and also functioned quantitatively in a multiplexed solution-phase assay. This same strategy was used to construct a second ribozyme selectively activated by the phosphorylated (active) form of ERK2. This approach is generally applicable to the development of ribozymes capable of monitoring post-translational modification of specific proteins.     

Computational and experimental validation of morin as adenosine deaminase inhibitor

Adenosine deaminase (ADA) is one of the major enzymes involved in purin metabolism, it has a significant role in cell growth and differentiation. Over-activity of ADA has been noticed in some pathology, like malignancy and inflammation and makes it an attractive target for the development of drugs for such diseases. In the present study, ADA inhibitory activity of morin, a bioactive flavonoid, was assessed through computational and biophysical methods. The enzyme kinetics data showed that morin is a competitive inhibitor of ADA. Binding energy calculated from ITC analysis was ?7.11?kcal/mol. Interaction of morin with ADA was also studied using fluorescence quenching method. Molecular docking studies revealed the structural details of the interaction. Molecular dynamics study in explicit solvent was also conducted to assess the structural stability of protein ligand complex.     

Deprotonation stimulates productive folding in allosteric TRAP hammerhead ribozymes

Hammerhead ribozymes in crystals change conformation in response to deprotonation of the nucleophilic 2' OH, thereby aligning the hydroxyl for in-line displacement at the scissile phosphate. Published data do not address whether deprotonation affects folding in solution. Allosteric hammerhead "TRAPs," when activated by the appropriate oligonucleotide, show the expected log-linear relation between initial cleavage rate and pH. In contrast, attenuated TRAPs shows biphasic kinetics in which a rapid burst is followed by slow cleavage that is nearly independent of pH. Attenuated ribozymes are stimulated by urea at both low and high pH, confirming that rearrangement of secondary structure is rate-limiting for the attenuated ribozymes once they have folded. Plots of burst magnitude versus pH in the absence of urea show a sharp transition around pH 8.3, which is near the kinetic pKa for the cleavage reaction in Mg2+. Raising the pH after folding at pH 7.5 did not activate attenuated ribozymes even when the RNA was incubated at the elevated pH for extended periods prior to addition of Mg2+. In contrast, lowering the pH after folding at pH 9.5 rapidly re-established attenuation. Deprotonation of the ribozyme-substrate complex thus appears to alter the folding landscape such that a metastable "pre-activated" complex forms before the thermodynamically more stable attenuated state can be attained. From the initial partition into active and inactive conformers, we estimate that this deprotonation contributes approximately 1.2 kcal/mol toward stabilization of the active fold at a crucial step during folding of the TRAP. Assuming that the nucleophilic 2' OH is the relevant acid, its deprotonation would thus serve a dual role of favoring productive fold and enhancing the nucleophilicity of this oxygen.     

Engineered allosteric ribozymes that respond to specific divalent metal ions

In vitro selection was used to isolate five classes of allosteric hammerhead ribozymes that are triggered by binding to certain divalent metal ion effectors. Each of these ribozyme classes are similarly activated by Mn2+, Fe2+, Co2+, Ni2+, Zn2+ and Cd2+, but their allosteric binding sites reject other divalent metals such as Mg2+, Ca2+ and Sr2+. Through a more comprehensive survey of cations, it was determined that some metal ions (Be2+, Fe3+, Al3+, Ru2+ and Dy2+) are extraordinarily disruptive to the RNA structure and function. Two classes of RNAs examined in greater detail make use of conserved nucleotides within the large internal bulges to form critical structures for allosteric function. One of these classes exhibits a metal-dependent increase in rate constant that indicates a requirement for the binding of two cation effectors. Additional findings suggest that, although complex allosteric functions can be exhibited by small RNAs, larger RNA molecules will probably be required to form binding pockets that are uniquely selective for individual cation effectors.     

Generating new ligand-binding RNAs by affinity maturation and disintegration of allosteric ribozymes

Allosteric ribozymes are engineered RNAs that operate as molecular switches whose rates of catalytic activity are modulated by the binding of specific effector molecules. New RNA molecular switches can be created by using "allosteric selection," a molecular engineering process that combines modular rational design and in vitro evolution strategies. In this report, we describe the characterization of 3',5'-cyclic nucleotide monophosphate (cNMP)-dependent hammerhead ribozymes that were created using allosteric selection (Koizumi et al., Nat Struct Biol, 1999, 6:1062-1071). Artificial phylogeny data generated by random mutagenesis and reselection of existing cGMP-, cCMP-, and cAMP-dependent ribozymes indicate that each is comprised of distinct effector-binding and catalytic domains. In addition, patterns of nucleotide covariation and direct mutational analysis both support distinct secondary-structure organizations for the effector-binding domains. Guided by these structural models, we were able to disintegrate each allosteric ribozyme into separate ligand-binding and catalytic modules. Examinations of the independent effector-binding domains reveal that each retains its corresponding cNMP-binding function. These results validate the use of allosteric selection and modular engineering as a means of simultaneously generating new nucleic acid structures that selectively bind ligands. Furthermore, we demonstrate that the binding affinity of an allosteric ribozyme can be improved through random mutagenesis and allosteric selection under conditions that favor tighter binding. This "affinity maturation" effect is expected to be a valuable attribute of allosteric selection as future endeavors seek to apply engineered allosteric ribozymes as biosensor components and as controllable genetic switches.     

Design of allosteric hammerhead ribozymes activated by ligand-induced structure stabilization.

BACKGROUND: Ribozymes can function as allosteric enzymes that undergo a conformational change upon ligand binding to a site other than the active site. Although allosteric ribozymes are not known to exist in nature, nucleic acids appear to be well suited to display such advanced forms of kinetic control. Current research explores the mechanisms of allosteric ribozymes as well as the strategies and methods that can be used to create new controllable enzymes. RESULTS: In this study, we exploit the modular nature of certain functional RNAs to engineer allosteric ribozymes that are activated by flavin mononucleotide (FMN) or theophylline. By joining an FMN- or theophylline-binding domain to a hammerhead ribozyme by different stem II elements, we have identified a minimal connective bridge comprised of a G.U wobble pair that is responsive to ligand binding. Binding of FMN or theophylline to its allosteric site induces a conformational change in the RNA that stabilizes the wobble pair and ultimately favors the active form of the catalytic core. These ligand-sensitive ribozymes exhibit rate enhancements of more than 100-fold in the presence of FMN and of approximately 40-fold in the presence of theophylline. CONCLUSIONS: An adaptive strategy for modular rational design has proven to be an effective approach to the engineering of novel allosteric ribozymes. This strategy was used to create allosteric ribozymes that function by a mechanism involving ligand-induced structure stabilization. Conceivably, similar engineering strategies and allosteric mechanisms could be used to create a variety of novel allosteric ribozymes that function with other effector molecules.     

Working at the cutting edge: the creation of allosteric ribozymes

The appropriate folding of catalytic RNA is a prerequisite for effective catalysis. A novel ribozyme, the maxizyme, has been generated and its activity can be controlled allosterically. The maxizymes work both in vitro and in vivo indicating the potential utility of this novel class of ribozyme as a gene-inactivating agent with a biosensor function.     

Computational prediction and experimental validation of microRNAs in the brown alga Ectocarpus siliculosus

We used an in silico approach to predict microRNAs (miRNAs) genome-wide in the brown alga Ectocarpus siliculosus. As brown algae are phylogenetically distant from both animals and land plants, our approach relied on features shared by all known organisms, excluding sequence conservation, genome localization and pattern of base-pairing with the target. We predicted between 500 and 1500 miRNAs candidates, depending on the values of the energetic parameters used to filter the potential precursors. Using quantitative polymerase chain reaction assays, we confirmed the existence of 22 miRNAs among 72 candidates tested, and of 8 predicted precursors. In addition, we compared the expression of miRNAs and their precursors in two life cycle states (sporophyte, gametophyte) and under salt stress. Several miRNA precursors, Argonaute and DICER messenger RNAs were differentially expressed in these conditions. Finally, we analyzed the gene organization and the target functions of the predicted candidates. This showed that E. siliculosus miRNA genes are, like plant miRNA genes, rarely clustered and, like animal miRNA genes, often located in introns. Among the predicted targets, several widely conserved functional domains are significantly overrepresented, like kinesin, nucleotide-binding/APAF-1, R proteins and CED-4 (NB-ARC) and tetratricopeptide repeats. The combination of computational and experimental approaches thus emphasizes the originality of molecular and cellular processes in brown algae.     

Computational prediction of efficient splice sites for trans-splicing ribozymes

Group I introns have been engineered into trans-splicing ribozymes capable of replacing the 3'-terminal portion of an external mRNA with their own 3'-exon. Although this design makes trans-splicing ribozymes potentially useful for therapeutic application, their trans-splicing efficiency is usually too low for medical use. One factor that strongly influences trans-splicing efficiency is the position of the target splice site on the mRNA substrate. Viable splice sites are currently determined using a biochemical trans-tagging assay. Here, we propose a rapid and inexpensive alternative approach to identify efficient splice sites. This approach involves the computation of the binding free energies between ribozyme and mRNA substrate. We found that the computed binding free energies correlate well with the trans-splicing efficiency experimentally determined at 18 different splice sites on the mRNA of chloramphenicol acetyl transferase. In contrast, our results from the trans-tagging assay correlate less well with measured trans-splicing efficiency. The computed free energy components suggest that splice site efficiency depends on the following secondary structure rearrangements: hybridization of the ribozyme's internal guide sequence (IGS) with mRNA substrate (most important), unfolding of substrate proximal to the splice site, and release of the IGS from the 3'-exon (least important). The proposed computational approach can also be extended to fulfill additional design requirements of efficient trans-splicing ribozymes, such as the optimization of 3'-exon and extended guide sequences.     

A novel strategy for selection of allosteric ribozymes yields RiboReporter sensors for caffeine and aspartame

We have utilized in vitro selection technology to develop allosteric ribozyme sensors that are specific for the small molecule analytes caffeine or aspartame. Caffeine- or aspartame-responsive ribozymes were converted into fluorescence-based RiboReporter™ sensor systems that were able to detect caffeine or aspartame in solution over a concentration range from 0.5 to 5 mM. With read-times as short as 5 min, these caffeine- or aspartame-dependent ribozymes function as highly specific and facile molecular sensors. Interestingly, successful isolation of allosteric ribozymes for the analytes described here was enabled by a novel selection strategy that incorporated elements of both modular design and activity-based selection methods typically used for generation of catalytic nucleic acids.     

Adaptive combinatorial design to explore large experimental spaces: approach and validation

Systems biology requires mathematical tools not only to analyse large genomic datasets, but also to explore large experimental spaces in a systematic yet economical way. We demonstrate that two-factor combinatorial design (CD), shown to be useful in software testing, can be used to design a small set of experiments that would allow biologists to explore larger experimental spaces. Further, the results of an initial set of experiments can be used to seed further 'Adaptive' CD experimental designs. As a proof of principle, we demonstrate the usefulness of this Adaptive CD approach by analysing data from the effects of six binary inputs on the regulation of genes in the N-assimilation pathway of Arabidopsis. This CD approach identified the more important regulatory signals previously discovered by traditional experiments using far fewer experiments, and also identified examples of input interactions previously unknown. Tests using simulated data show that Adaptive CD suffers from fewer false positives than traditional experimental designs in determining decisive inputs, and succeeds far more often than traditional or random experimental designs in determining when genes are regulated by input interactions. We conclude that Adaptive CD offers an economical framework for discovering dominant inputs and interactions that affect different aspects of genomic outputs and organismal responses.     

Computational modeling and experimental validation of odor detection behaviors of classically conditioned parasitic wasp,Microplitis croceipes

A prototype chemical sensor named Wasp hound® that utilizes five classically conditioned parasitoid wasps, Microplitis croceipes (Cresson) (Hymenoptera: Braconidae), to detect volatile odors was successfully implemented in a previous study. To improve the odor‐detecting ability of Wasp Hound®, searching behaviors of an individual wasp in a confined area are studied and modeled through stochastic differential equations in this paper. The wasps are conditioned to 20 mg of coffee when associated with food and subsequently, tested to 5, 10, 20, and 40 mg of coffee. A stochastic model is developed and validated based on three positive behavioral responses (walking, rotation around odor source, and self‐rotation) from conditioned wasps at four different test dosages. The model is capable to reproducing the behaviors of conditioned wasps, and can be used to improve the ability of Wasp Hound® to assess changes in odor concentration. The model simulation results show the behaviors of conditioned wasps are significantly different when tested at different coffee dosages. We conjecture that the searching behaviors of conditioned wasps are based on the temporal and spatial neuron activity of olfactory receptor neurons and glomeruli, which are strongly correlated to the training dosages. The overall results demonstrate the utility of mathematical models for interpreting experimental observations, gaining novel insights into the dynamic behavior of classically conditioned wasps, as well as broadening the practical uses of Wasp Hound. © 2014 American Institute of Chemical Engineers Biotechnol. Prog., 31:596–606, 2015