Hammerhead Ribozyme Research Articles

An allosteric ribozyme generator and an inverse folding ribozyme generator: Two computer programs for automated computational design of oligonucleotide-sensing allosteric hammerhead ribozymes with YES Boolean logic function based on experimentally validated algorithms

Designing oligonucleotide-sensing ribozymes using computational approaches is advantageous to in vitro selection methods for efficiency and accuracy. Allosteric ribozymes can be computationally designed for various applications in gene therapy, designer gene control systems, biosensors, and molecular computation. Here we present two programs, the allosteric Ribozyme Generator (RG) and the Inverse Folding Ribozyme Generator (IFRG), engineered to generate allosteric ribozymes with YES logic. The RG computes allosteric ribozyme sequences’ secondary structure using the minimal sequence of the hammerhead ribozyme by inserting oligonucleotide binding site (OBS) elements in the second stem. The IFRG program uses inverse folding to generate allosteric ribozyme sequences with OBS bearing distinct sequences and similar folding. For the generation of the OBS sequences, random search algorithms are employed. Allosteric ribozyme sequences generated by the RG can be used as a matrix for the IFRG program. This approach applies RNA-folding algorithms based on applying thermodynamic parameters using the partition function of the RNAfold, and the RNAinverse source codes from the Vienna RNA folding package. The two algorithms apply dynamic programming and random search algorithms to generate in silico allosteric ribozymes with predefined properties within minutes using a personal computer with over 90% accuracy, without high computation power as experimentally validated and published by us previously.

Spectroscopic characterization of Mn2+ and Cd2+ coordination to phosphorothioates in the conserved A9 metal site of the hammerhead ribozyme

Phosphorothioate modifications have widespread use in the field of nucleic acids. As substitution of sulfur for oxygen can alter metal coordination preferences, the phosphorothioate metal-rescue experiment is a powerful method for identifying metal coordination sites that influence specific properties in a large RNAs. The A9/G10.1 metal binding site of the hammerhead ribozyme (HHRz) has previously been shown to be functionally important through phosphorothioate rescue experiments. While an A9-SRp substitution is inhibitory in Mg2+, thiophilic Cd2+ rescues HHRz activity. Mn2+ is also often used in phosphorothioate metal-rescue studies but does not support activity for the A9-SRp HHRz. Here, we use EPR, electron spin-echo envelope modulation (ESEEM), and X-ray absorption spectroscopic methods to directly probe the structural consequences of Mn2+ and Cd2+ coordination to Rp and Sp phosphorothioate modifications at the A9/G10.1 site in the truncated hammerhead ribozyme (tHHRz). The results demonstrate that while Cd2+ does indeed bind to S in the thio-substituted ligand, Mn2+ coordinates to the non‑sulfur oxo group of this phosphorothioate, regardless of isomer. Computational models demonstrate the energetic preference of MnO over MnS coordination in metal-dimethylthiophosphate models. In the case of the tHHRz, the resulting Mn2+ coordination preference of oxygen in either Rp or Sp A9 phosphorothioates differentially tunes catalytic activity, with MnO coordination in the A9-SRp phosphorothioate enzyme being inhibitory.

Rolling circle RNA synthesis catalyzed by RNA.

RNA-catalyzed RNA replication is widely considered a key step in the emergence of life's first genetic system. However, RNA replication can be impeded by the extraordinary stability of duplex RNA products, which must be dissociated for re-initiation of the next replication cycle. Here, we have explored rolling circle synthesis (RCS) as a potential solution to this strand separation problem. We observe sustained RCS by a triplet polymerase ribozyme beyond full-length circle synthesis with strand displacement yielding concatemeric RNA products. Furthermore, we show RCS of a circular Hammerhead ribozyme capable of self-cleavage and re-circularization. Thus, all steps of a viroid-like RNA replication pathway can be catalyzed by RNA alone. Finally, we explore potential RCS mechanisms by molecular dynamics simulations, which indicate a progressive build-up of conformational strain upon RCS with destabilization of nascent strand 5'- and 3'-ends. Our results have implications for the emergence of RNA replication and for understanding the potential of RNA to support complex genetic processes.

A scenario for the emergence of protoviroids in the RNA world and for their further evolution into viroids and viroid-like RNAs by modular recombinations and mutations.

Viroids are tiny, circular, and noncoding RNAs that are able to replicate and systemically infect plants. The smallest known pathogens, viroids have been proposed to represent survivors from the RNA world that likely preceded the cellular world currently dominating life on the earth. Although the small, circular, and compact nature of viroid genomes, some of which are also endowed with catalytic activity mediated by hammerhead ribozymes, support this proposal, the lack of feasible evolutionary routes and the identification of hammerhead ribozymes in a large number of DNA genomes of organisms along the tree of life have led some to question such a proposal. Here, we reassess the origin and subsequent evolution of viroids by complementing phylogenetic reconstructions with molecular data, including the primary and higher-order structure of the genomic RNAs, their replication, and recombination mechanisms and selected biological information. Features of some viroid-like RNAs found in plants, animals, and possibly fungi are also considered. The resulting evolutionary scenario supports the emergence of protoviroids in the RNA world, mainly as replicative modules, followed by a further increase in genome complexity based on module/domain shuffling and combination and mutation. Such a modular evolutionary scenario would have facilitated the inclusion in the protoviroid genomes of complex RNA structures (or coding sequences, as in the case of hepatitis delta virus and delta-like agents), likely needed for their adaptation from the RNA world to a life based on cells, thus generating the ancestors of current infectious viroids and viroid-like RNAs. Other noninfectious viroid-like RNAs, such as retroviroid-like RNA elements and retrozymes, could also be derived from protoviroids if their reverse transcription and integration into viral or eukaryotic DNA, respectively, are considered as a possible key step in their evolution. Comparison of evidence supporting a general and modular evolutionary model for viroids and viroid-like RNAs with that favoring alternative scenarios provides reasonable reasons to keep alive the hypothesis that these small RNA pathogens may be relics of a precellular world.

Toehold-mediated strand displacement to measure released product from self-cleaving ribozymes.

This paper presents a probe comprising a fluorophore and a quencher, enabling measurement of released product from self-cleaving hammerhead ribozyme, without labeled RNA molecules, regular sampling or use of polyacrylamide gels. The probe is made of two DNA strands; one strand is labeled with a fluorophore at its 5'-end, while the other strand is labeled with a quencher at its 3'-end. These two DNA strands are perfectly complementary, but with a 3'-overhang of the fluorophore strand. These unpaired nucleotides act as a toehold, which is utilized by a detached cleaved fragment (coming from a self-cleaving hammerhead ribozyme) as the starting point for a strand displacement reaction. This reaction causes the separation of the fluorophore strand from the quencher strand, culminating in fluorescence, detectable in a plate reader. Notably, the emitted fluorescence is proportional to the amount of detached cleaved-off RNAs, displacing the DNA quencher strand. This method can replace or complement radio-hazardous unstable 32P as a method of measurement of the product release from ribozyme cleavage reactions; it also eliminates the need for polyacrylamide gels, for the same purpose. Critically, this method allows to distinguish between the total amount of cleaved ribozymes and the amount of detached fragments, resulting from that cleavage reaction.

Production of Recombinant RNA in Escherichia coli Using Eggplant Latent Viroid as a Scaffold.

Similar to viruses, viroids can also be engineered and transformed into useful biotechnological tools. We describe here a viroid-based system to produce large amounts of recombinant RNA in Escherichia coli. A precursor of eggplant latent viroid (ELVd), with the RNA of interest inserted between positions U245 and U246, is co-expressed in E. coli along the chloroplastic isoform of the eggplant tRNA ligase, the enzyme that mediates the circularization of this viroid in the infected plants. In the bacterial cells, the chimeric ELVd-RNA-of-interest precursor self-cleaves through the embedded hammerhead ribozymes, and the monomer is recognized and circularized by the co-expressed tRNA ligase. The resulting circular RNA, likely bound to the tRNA ligase, accumulates to a high concentration in the bacterial cells.

CyPhyRNA-seq: a genome-scale RNA-seq method to detect active self-cleaving ribozymes by capturing RNAs with 2ʹ,3ʹ cyclic phosphates and 5ʹ hydroxyl ends

ABSTRACT Self-cleaving ribozymes are catalytically active RNAs that cleave themselves into a 5′-fragment with a 2′,3′-cyclic phosphate and a 3′-fragment with a 5′-hydroxyl. They are widely applied for the construction of synthetic RNA devices and RNA-based therapeutics. However, the targeted discovery of self-cleaving ribozymes remains a major challenge. We developed a transcriptome-wide method, called cyPhyRNA-seq, to screen for ribozyme cleavage fragments in total RNA extract. This approach employs the specific ligation-based capture of ribozyme 5′-fragments using a variant of the Arabidopsis thaliana tRNA ligase we engineered. To capture ribozyme 3′-fragments, they are enriched from total RNA by enzymatic treatments. We optimized and enhanced the individual steps of cyPhyRNA-seq in vitro and in spike-in experiments. Then, we applied cyPhyRNA-seq to total RNA isolated from the bacterium Desulfovibrio vulgaris and detected self-cleavage of the three predicted type II hammerhead ribozymes, whose activity had not been examined to date. cyPhyRNA-seq can be used for the global analysis of active self-cleaving ribozymes with the advantage to capture both ribozyme cleavage fragments from total RNA. Especially in organisms harbouring many self-cleaving RNAs, cyPhyRNA-seq facilitates the investigation of cleavage activity. Moreover, this method has the potential to be used to discover novel self-cleaving ribozymes in different organisms.

Cross-Chiral, RNA-Catalyzed Exponential Amplification of RNA.

Informational macromolecules in biology are composed of subunits of a single handedness, d-nucleotides in nucleic acids and l-amino acids in proteins. Although this chiral uniformity may be expedient, it is not a chemical necessity, as demonstrated by the recent example of an RNA enzyme that catalyzes the RNA-templated polymerization of RNA molecules of the opposite handedness. This reaction, when carried out iteratively, can provide the basis for exponential amplification of RNA molecules and the information they contain. By carrying out thermal cycling, analogous to the polymerase chain reaction, and supplying oligonucleotide building blocks that comprise both the functional strand of RNA and its complement, cross-chiral exponential amplification was achieved. This process was used to amplify the l-RNA form of the hammerhead ribozyme, catalyzed by the d-RNA form of the polymerase. The resulting l-hammerhead exhibits the expected activity in cleaving a corresponding l-RNA substrate. Exponential amplification was also carried out within individual droplets of a water-in-oil emulsion. The ability to amplify enantio-RNAs, both in bulk solution and within compartments, provides a means to evolve cross-chiral RNA polymerases based on the function of the RNAs they produce.

Satellite-Like W-Elements: Repetitive, Transcribed, and Putative Mobile Genetic Factors with Potential Roles for Biology and Evolution of Schistosoma mansoni.

A large portion of animal and plant genomes consists of noncoding DNA. This part includes tandemly repeated sequences and gained attention because it offers exciting insights into genome biology. We investigated satellite-DNA elements of the platyhelminth Schistosoma mansoni, a parasite with remarkable biological features. Schistosoma mansoni lives in the vasculature of humans causing schistosomiasis, a disease of worldwide importance. Schistosomes are the only trematodes that have evolved separate sexes, and the sexual maturation of the female depends on constant pairing with the male. The schistosome karyotype comprises eight chromosome pairs, males are homogametic (ZZ) and females are heterogametic (ZW). Part of the repetitive DNA of S. mansoni are W-elements (WEs), originally discovered as female-specific satellite DNAs in the heterochromatic block of the W-chromosome. Based on new genome and transcriptome data, we performed a reanalysis of the W-element families (WEFs). Besides a new classification of 19 WEFs, we provide first evidence for stage-, sex-, pairing-, gonad-, and strain-specific/preferential transcription of WEs as well as their mobile nature, deduced from autosomal copies of full-length and partial WEs. Structural analyses suggested roles as sources of noncoding RNA-like hammerhead ribozymes, for which we obtained functional evidence. Finally, the variable WEF occurrence in different schistosome species revealed remarkable divergence. From these results, we propose that WEs potentially exert enduring influence on the biology of S. mansoni. Their variable occurrence in different strains, isolates, and species suggests that schistosome WEs may represent genetic factors taking effect on variability and evolution of the family Schistosomatidae.

Dual Promoters Improve the Rescue of Recombinant Measles Virus in Human Cells.

Reverse genetics is a technology that allows the production of a virus from its complementary DNA (cDNA). It is a powerful tool for analyzing viral genes, the development of novel vaccines, and gene delivery vectors. The standard reverse genetics protocols are laborious, time-consuming, and inefficient for negative-strand RNA viruses. A new reverse genetics platform was established, which increases the recovery efficiency of the measles virus (MV) in human 293-3-46 cells. The novel features compared with the standard system involving 293-3-46 cells comprise (a) dual promoters containing the RNA polymerase II promoter (CMV) and the bacteriophage T7 promoter placed in uni-direction on the same plasmid to enhance RNA transcription; (b) three G nucleotides added just after the T7 promoter to increase the T7 RNA polymerase activity; and (c) two ribozymes, the hairpin hammerhead ribozyme (HHRz), and the hepatitis delta virus ribozyme (HDVrz), were used to cleavage the exact termini of the antigenome RNA. Full-length antigenome cDNA of MV of the wild type IC323 strain or the vaccine AIK-C strain was inserted into the plasmid backbone. Both virus strains were easily rescued from their respective cloned cDNA. The rescue efficiency increased up to 80% compared with the use of the standard T7 rescue system. We assume that this system might be helpful in the rescue of other human mononegavirales.

An Ancient Clade of Penelope-Like Retroelements with Permuted Domains Is Present in the Green Lineage and Protists, and Dominates Many Invertebrate Genomes.

Penelope-like elements (PLEs) are an enigmatic clade of retrotransposons whose reverse transcriptases (RTs) share a most recent common ancestor with telomerase RTs. The single ORF of canonical endonuclease (EN)+ PLEs encodes RT and a C-terminal GIY–YIG EN that enables intrachromosomal integration, whereas EN− PLEs lack EN and are generally restricted to chromosome termini. EN+ PLEs have only been found in animals, except for one case of horizontal transfer to conifers, whereas EN− PLEs occur in several kingdoms. Here, we report a new, deep-branching PLE clade with a permuted domain order, whereby an N-terminal GIY–YIG EN is linked to a C-terminal RT by a short domain with a characteristic CxC motif. These N-terminal EN+ PLEs share a structural organization, including pseudo-LTRs and complex tandem/inverted insertions, with canonical EN+ PLEs from Penelope/Poseidon, Neptune, and Nematis clades, and show insertion bias for microsatellites, but lack canonical hammerhead ribozyme motifs. However, their phylogenetic distribution is much broader. The Naiads, found in numerous invertebrate phyla, can reach tens of thousands of copies per genome. In spiders and clams, Naiads independently evolved to encode selenoproteins containing multiple selenocysteines. Chlamys, which lack the CCHH motif universal to PLE ENs, occur in green algae, spike mosses (targeting ribosomal DNA), and slime molds. Unlike canonical PLEs, RTs of N-terminal EN+ PLEs contain the insertion-in-fingers domain (IFD), strengthening the link between PLEs and telomerases. Additionally, we describe Hydra, a novel metazoan C-terminal EN+ clade. Overall, we conclude that PLE diversity, taxonomic distribution, and abundance are comparable with non-LTR and LTR-retrotransposons.

Kinetic Study of the Avocado Sunblotch Viroid Self-Cleavage Reaction Reveals Compensatory Effects between High-Pressure and High-Temperature: Implications for Origins of Life on Earth

Simple SummaryViroids remain the smallest infectious agents ever discovered. They are found in plants and consist of single-stranded non-coding circular RNA. Due to their simplicity, viroids are considered relics of an ancient RNA World that may have originated in the deep seas near hydrothermal vents where temperature and pressure are both elevated. To test this hypothesis, a synthetic avocado sunblotch viroid, whose structure contain an autocatalytic hammerhead ribozyme, was subjected to increased pressure (from atmospheric pressure to 300 MPa) at different temperatures (0–65 °C) and the reaction rate constant of the catalytic activity was calculated for each condition. The results obtained allowed calculation of the positive activation volume of this viroid and revealed a compensatory effect between pressure and temperature. In conclusion, these results not only exemplify the plasticity of RNA and support the RNA World hypothesis, but also highlight the usefulness of the hydrostatic pressure in understanding the structure–function relationships of biomacromolecules.A high pressure apparatus allowing one to study enzyme kinetics under pressure was used to study the self-cleavage activity of the avocado sunblotch viroid. The kinetics of this reaction were determined under pressure over a range up to 300 MPa (1–3000 bar). It appears that the initial rate of this reaction decreases when pressure increases, revealing a positive ΔV≠ of activation, which correlates with the domain closure accompanying the reaction and the decrease of the surface of the viroid exposed to the solvent. Although, as expected, temperature increases the rate of the reaction whose energy of activation was determined, it appeared that it does not significantly influence the ΔV≠ of activation and that pressure does not influence the energy of activation. These results provide information about the structural aspects or this self-cleavage reaction, which is involved in the process of maturation of this viroid. The behavior of ASBVd results from the involvement of the hammerhead ribozyme present at its catalytic domain, indeed a structural motif is very widespread in the ancient and current RNA world.

Enhancement of the Catalytic Activity of Hammerhead Ribozymes by Organic Cations.

Catalytic turnover is important for the application of ribozymes to biotechnology. However, the turnover is often impaired because of the intrinsic high stability of base pairs with cleaved RNA products. Here, organic cations were used as additives to improve the catalytic performance of hammerhead ribozyme constructs that exhibit different kinetic behaviors. Kinetic analysis of substrate cleavage demonstrated that bulky cations, specifically tetra-substituted ammonium ions containing pentyl groups or a benzyl group, have the ability to greatly increase the turnover rate of the ribozymes. Thermal stability analysis of RNA structures revealed that the bulky cations promote the dissociation of cleaved products and refolding of incorrectly folded structures with small disruption of the catalytic structure. The use of bulky cations is a convenient method for enhancing the catalytic activity of hammerhead ribozymes, and the approach may be useful for advancing ribozyme technologies.

Ribozyme-processed guide RNA enhances virus-mediated plant genome editing.

In virus-induced gene-editing system, subgenomic promoters have been used to express guide RNAs (gRNAs). However, the transcription initiation site of the subgenomic promoters remains elusive. Here, we examined the sequence of gRNAs expressed by subgenomic promoters and found the variable length of overhangs at 5'-end of gRNAs. The overhangs at 5'-end of gRNA decrease the cleavage activity of SpCas9. To overcome this problem, we inserted hammerhead ribozyme between the subgenomic promoter and gRNA and confirmed that gRNAs with a precise 5'-end increase the editing efficacy in wild tobacco. This system will be widely used for editing target genes in plants with high efficiency.

A single transcript CRISPR/Cas9 mediated mutagenesis of CaERF28 confers anthracnose resistance in chilli pepper (Capsicum annuum L.).

T-DNA-free homozygous mutant lines developed through a single transcript CRISPR/Cas9 system harboring the desired modification in the CaERF28 locus exhibited significantly enhanced resistance to the anthracnose pathogen Colletotrichum truncatum coupled with the improved expression of defense responsive genes. Anthracnose, caused by Colletotrichum species, is a major disease of chilli (Capsicum annuum) accounting for significant pre- and post-harvest yield losses across the tropical and subtropical regions of the world. Management of chilli anthracnose using traditional methods have not met with noticeable success. In the present study, we have demonstrated an enhanced anthracnose resistance through a single transcript unit CRISPR/Cas9 mediated alteration of the susceptibility gene CaERF28 in C. annuum. A construct with a single Pol II promoter-driven expression of Cas9, sgRNA and a hammerhead ribozyme (RZ) was designed to modify the CaERF28 gene in the susceptible chilli genotype Arka Lohit. Fourty-five C-ERF28-induced mutant lines (72.5%) were identified from 62 T0 transgenic plants. Further, simultaneously targeted multiple sites within CaERF28 showed increased mutation (85.7%) efficiency. DNA sequence analysis showed that these plants harboured multiple InDels at the target site. The allelic mutants of C-ERF28 were transferred to the following generations by simple Mendelian inheritance. Segregation in the T1 and T2 generations resulted in the identification of T-DNA free and marker-free C-ERF28 mutant lines. Five homozygous mutants demonstrated enhanced resistance to anthracnose compared to wild type as demonstrated by reduced spore count and fungal growth as well as induced expression of defense-related genes. Our results demonstrated that the STU-CRISPR/Cas9 mediated editing of the CaERF28 gene is a rapid, safe and versatile approach for enhancing anthracnose resistance in chilli pepper and pave way for its utilization in the improvement of other solanaceous crops.

The hammerhead self-cleaving motif as a precursor to complex endonucleolytic ribozymes.

Connections between distinct catalytic RNA motifs through networks of mutations that retain catalytic function (neutral networks) were likely central to the evolution of biocatalysis. Despite suggestions that functional RNAs collectively form an interconnected web of neutral networks, little evidence has emerged to demonstrate the existence of such intersecting networks in naturally occurring RNAs. Here we show that neutral networks of two naturally occurring, seemingly unrelated endonucleolytic ribozymes, the hammerhead (HH) and hairpin (HP), intersect. Sequences at the intersection of these networks exhibit catalytic functions corresponding to both ribozymes by potentially populating both catalytic folds and enable a smooth crossover between the two. Small and structurally simple endonucleolytic motifs like the HH ribozyme could, through mutational walks along their neutral networks, encounter novel catalytic phenotypes, and structurally flexible, bifunctional sequences at the intersection of these networks could have acted as nodes for evolutionary diversification in an RNA world. Considering the simplicity and small size of the HH ribozyme, we propose that this self-cleaving motif could have been a precursor to other more complex endonucleolytic ribozymes. More generally, our results suggest that RNAs that possess distinct sequences, structures, and catalytic functions, can potentially share evolutionary history through mutational connections in sequence space.

A Threose Nucleic Acid Enzyme with RNA Ligase Activity.

Threose nucleic acid (TNA) has been considered a potential RNA progenitor in evolution due to its chemical simplicity and base pairing property. Catalytic TNA sequences with RNA ligase activities might have facilitated the transition to the RNA world. Here we report the isolation of RNA ligase TNA enzymes by in vitro selection. The identified TNA enzyme T8-6 catalyzes the formation of a 2'-5' phosphoester bond between a 2',3'-diol and a 5'-triphosphate group, with a kobs of 1.1 × 10-2 min-1 (40 mM Mg2+, pH 9.0). For efficient reaction, T8-6 requires UA|GA at the ligation junction and tolerates variations at other substrate positions. Functional RNAs such as hammerhead ribozyme can be prepared by T8-6-catalyzed ligation, with site-specific introduction of a 2'-5' linkage. Together, this work provides experimental support for TNA as a plausible pre-RNA genetic polymer and also offers an alternative molecular tool for biotechnology.

A synthetic RNA-based biosensor for fructose-1,6-bisphosphate that reports glycolytic flux

RNA-based sensors for intracellular metabolites are a promising solution to the emerging issue of metabolic heterogeneity. However, their development, i.e., the conversion of an aptamer into an invivo-functional intracellular metabolite sensor, still harbors challenges. Here, we accomplished this for the glycolytic flux-signaling metabolite, fructose-1,6-bisphosphate (FBP). Starting from invitro selection of an aptamer, we constructed device libraries with a hammerhead ribozyme as actuator. Using high-throughput screening in yeast with fluorescence-activated cell sorting (FACS), next-generation sequencing, and genetic-environmental perturbations to modulate the intracellular FBP levels, we identified a sensor that generates ratiometric fluorescent readout. An abrogated response in sensor mutants and occurrence of two sensor conformations-revealed by RNA structural probing-indicated invivo riboswitching activity. Microscopy showed that the sensor can differentiate cells with different glycolytic fluxes within yeast populations, opening research avenues into metabolic heterogeneity. We demonstrate the possibility to generate RNA-based sensors for intracellular metabolites for which no natural metabolite-binding RNA element exits.

Extension of the Variational Free Energy Profile and Multistate Bennett Acceptance Ratio Methods for High-Dimensional Potential of Mean Force Profile Analysis.

We redevelop the variational free energy profile (vFEP) method using a cardinal B-spline basis to extend the method for analyzing free energy surfaces (FESs) involving three or more reaction coordinates. We also implemented software for evaluating high-dimensional profiles based on the multistate Bennett acceptance ratio (MBAR) method which constructs an unbiased probability density from global reweighting of the observed samples. The MBAR method takes advantage of a fast algorithm for solving the unbinned weighted histogram (UWHAM)/MBAR equations which replaces the solution of simultaneous equations with a nonlinear optimization of a convex function. We make use of cardinal B-splines and multiquadric radial basis functions to obtain smooth, differentiable MBAR profiles in arbitrary high dimensions. The cardinal B-spline vFEP and MBAR methods are compared using three example systems that examine 1D, 2D, and 3D profiles. Both methods are found to be useful and produce nearly indistinguishable results. The vFEP method is found to be 150 times faster than MBAR when applied to periodic 2D profiles, but the MBAR method is 4.5 times faster than vFEP when evaluating unbounded 3D profiles. In agreement with previous comparisons, we find the vFEP method produces superior FESs when the overlap between umbrella window simulations decreases. Finally, the associative reaction mechanism of hammerhead ribozyme is characterized using 3D, 4D, and 6D profiles, and the higher-dimensional profiles are found to have smaller reaction barriers by as much as 1.5 kcal/mol. The methods presented here have been implemented into the FE-ToolKit software package along with new methods for network-wide free energy analysis in drug discovery.

Engineering Allosteric Ribozymes to Detect Thiamine Pyrophosphate in Whole Blood.

Thiamine deficiency contributes to several human diseases including Alzheimer's. As its biologically active form, thiamine pyrophosphate (TPP) has been considered as a potential biomarker for Alzheimer's disease (AD) based on several clinical reports that apparently lower blood TPP levels were found in patients with mild cognitive impairment to AD. However, highly sensitive and high-throughput detection of TPP in biological fluids remains an analytical challenge. Here, we report engineering RNA-based sensors to quantitatively measure TPP concentrations in whole blood samples with a detection limit down to a few nM. By fusing a TPP-specific aptamer with the hammerhead ribozyme for in vitro selection, we isolated an allosteric ribozyme with an EC50 value (68 nM) similar to the aptamer's KD value (50 nM) for TPP, which for the first time demonstrates the possibility to maintain the effector binding affinity of the aptamer in such engineered allosteric RNA constructs. Meanwhile, we developed a new blood sample preparation protocol to be compatible with RNA. By coupling the TPP-induced ribozyme cleavage event with isothermal amplification, we achieved fluorescence monitoring of whole blood TPP levels through the "mix-and-read" operation with high-throughput potential. We expect that the engineered TPP-sensing RNAs will facilitate clinical research on AD as well as other thiamine-related diseases.