Control of RNA Function by Conformational Design
Total Page:16
File Type:pdf, Size:1020Kb
Load more
Recommended publications
-
Design of Nucleic Acid Strands with Long Low-Barrier Folding Pathways
Nat Comput DOI 10.1007/s11047-016-9587-9 Design of nucleic acid strands with long low-barrier folding pathways Anne Condon1 · Bonnie Kirkpatrick2 · Ján Maˇnuch1 © The Author(s) 2017. This article is published with open access at Springerlink.com Abstract A major goal of natural computing is to design fully developed, that exploit base pairing interactions of biomolecules, such as nucleic acid sequences, that can be nucleic acids. Prominent examples include DNA strand dis- used to perform computations. We design sequences of placement systems (DSDs) (Seelig et al. 2006) and RNA nucleic acids that are “guaranteed” to have long folding origami systems (Geary and Andersen 2014). Our work here pathways relative to their length. This particular sequences is motivated by the goal of computing with a single RNA with high probability follow low-barrier folding pathways sequence as the nucleic acids of the sequence interact with that visit a large number of distinct structures. Long fold- each other. ing pathways are interesting, because they demonstrate that RNA sequences form folded structures in which pairs of natural computing can potentially support long and com- nucleic acids biochemically bond to each other. These bonds plex computations. Formally, we provide the first scalable change the physical energy of the sequence, and a given designs of molecules whose low-barrier folding pathways, sequence prefers to assume low-energy folded structures. with respect to a simple, stacked pair energy model, grow Folding is a dynamic process, constrained by kinetics, dur- superlinearly with the molecule length, but for which all sig- ing which an RNA sequence will move through a sequence nificantly shorter alternative folding pathways have an energy of structures with each differing from the previous one by barrier that is 2 − times that of the low-barrier pathway for the addition or removal of a single base pair. -
Algorithms for Nucleic Acid Sequence Design
Algorithms for Nucleic Acid Sequence Design Thesis by Joseph N. Zadeh In Partial Fulfillment of the Requirements for the Degree of Doctor of Philosophy California Institute of Technology Pasadena, California 2010 (Defended December 8, 2009) ii © 2010 Joseph N. Zadeh All Rights Reserved iii Acknowledgements First and foremost, I thank Professor Niles Pierce for his mentorship and dedication to this work. He always goes to great lengths to make time for each member of his research group and ensures we have the best resources available. Professor Pierce has fostered a creative environment of learning, discussion, and curiosity with a particular emphasis on quality. I am grateful for the tremendously positive influence he has had on my life. I am fortunate to have had access to Professor Erik Winfree and his group. They have been very helpful in pushing the limits of our software and providing fun test cases. I am also honored to have two other distinguished researchers on my thesis committee: Stephen Mayo and Paul Rothemund. All of the work presented in this thesis is the result of collaboration with extremely talented individuals. Brian Wolfe and I codeveloped the multiobjective design algorithm (Chapter 3). Brian has also been instru- mental in finessing details of the single-complex algorithm (Chapter 2) and contributing to the parallelization of NUPACK’s core routines. I would also like to thank Conrad Steenberg, the NUPACK software engineer (Chapter 4), who has significantly improved the performance of the site and developed robust secondary structure drawing code. Another codeveloper on NUPACK, Justin Bois, has been a good friend, mentor, and reliable coding partner. -
A Global Sampling Approach to Designing and Reengineering RNA
i “rnaensign” — 2017/7/1 — 22:02 — page 1 — #1 i i i Published online ?? Nucleic Acids Research, ??, Vol. ??, No. ?? 1–11 doi:10.1093/nar/gkn000 Aglobalsamplingapproachtodesigningandreengineering RNA secondary structures 1,2, 1, 3 1 1 Alex Levin ⇤,MieszkoLis ⇤,YannPonty ,CharlesW.O’Donnell ,SrinivasDevadas , 1,2, 1,2,4, Bonnie Berger †,J´erˆomeWaldisp¨uhl † 1 Computer Science and Artificial Intelligence Laboratory, MIT, Cambridge, MA 02139, USA. and 2 Department of Mathematics, MIT, Cambridge, MA 02139, USA. and 3 Laboratoire d’Informatique, Ecole´ Polytechnique, Palaiseau, France. and 4 School of Computer Science, McGill University, Montreal, QC H2A 3A7, Canada. Received ??; Revised ??; Accepted ?? ABSTRACT mechanisms (1), reprogramming cellular behavior (2), and designing logic circuits (3). The development of algorithms for designing artificial Here, we aim to design RNA sequences that fold into RNA sequences that fold into specific secondary structures specific secondary structures (a.k.a. inverse folding). Even has many potential biomedical and synthetic biology in this case, efficient computational formulations remain applications. To date, this problem remains computationally difficult, with no exact solutions known. Instead, the solutions difficult, and current strategies to address it resort available today rely on local search strategies and heuristics. to heuristics and stochastic search techniques. The Indeed, the computational difficulty of the RNA design most popular methods consist of two steps: First a problem was proven by M. Schnall-Levin et al. (4). random seed sequence is generated; next, this seed is One of the first and most widely known programs for progressively modified (i.e. mutated) to adopt the desired the RNA inverse folding problem is RNAinverse (5). -
Degeneracy and Genetic Assimilation in RNA Evolution Reza Rezazadegan1* and Christian Reidys1,2
Rezazadegan and Reidys BMC Bioinformatics (2018) 19:543 https://doi.org/10.1186/s12859-018-2497-3 RESEARCH ARTICLE Open Access Degeneracy and genetic assimilation in RNA evolution Reza Rezazadegan1* and Christian Reidys1,2 Abstract Background: The neutral theory of Motoo Kimura stipulates that evolution is mostly driven by neutral mutations. However adaptive pressure eventually leads to changes in phenotype that involve non-neutral mutations. The relation between neutrality and adaptation has been studied in the context of RNA before and here we further study transitional mutations in the context of degenerate (plastic) RNA sequences and genetic assimilation. We propose quasineutral mutations, i.e. mutations which preserve an element of the phenotype set, as minimal mutations and study their properties. We also propose a general probabilistic interpretation of genetic assimilation and specialize it to the Boltzmann ensemble of RNA sequences. Results: We show that degenerate sequences i.e. sequences with more than one structure at the MFE level have the highest evolvability among all sequences and are central to evolutionary innovation. Degenerate sequences also tend to cluster together in the sequence space. The selective pressure in an evolutionary simulation causes the population to move towards regions with more degenerate sequences, i.e. regions at the intersection of different neutral networks, and this causes the number of such sequences to increase well beyond the average percentage of degenerate sequences in the sequence space. We also observe that evolution by quasineutral mutations tends to conserve the number of base pairs in structures and thereby maintains structural integrity even in the presence of pressure to the contrary. -
Folding RNA/DNA Hybrid Duplexes
Vol. 00 no. 00 BIOINFORMATICS Pages 1–2 Folding RNA/DNA hybrid duplexes Ronny Lorenz 1,∗, Ivo L Hofacker 1,2 and Stephan H Bernhart 2,∗ 1Institute for Theoretical Chemistry, University Vienna, Wahringer¨ Straße 17, 1090 Vienna, Austria 2Research group BCB, Faculty of Computer Science, University Vienna 3Department of Bioinformatics, University of Leipzig, Haertelstraße 16-18, 04109 Leipzig, Germany Received on XXXXX; revised on XXXXX; accepted on XXXXX Associate Editor: XXXXXXX ABSTRACT (for a review see e.g. [11]) or even their dimers, there is, to our Motivation: While there are numerous programs that can predict knowledge, no program that also includes the possibility to predict RNA or DNA secondary structures, a program that predicts RNA/DNA the full secondary structure of RNA/DNA hetero-dimers. This may hetero-dimers is still missing. The lack of easy to use tools for be due to the lack of compiled energy parameters, as there are only predicting their structure may be in part responsible for the small stacking energies available today. We introduce here a possibility number of reports of biologically relevant RNA/DNA hetero-dimers. to predict the secondary structure of such hetero-dimers within the Results: We present here an extension to the widely used widely used ViennaRNA Package. ViennaRNA Package [6] for the prediction of the structure of RNA/DNA hetero-dimers. Availability: http://www.tbi.univie.ac.at/˜ronny/RNA/vrna2.html 2 APPROACH Contact: [email protected] 2.1 Adaptation of the ViennaRNA Package For the prediction of RNA/DNA hetero-dimers, three distinct energy 1 INTRODUCTION parameter sets (RNA, DNA and mixed) are necessary. -
Supplemental Material Predicting Transfer RNA Gene Activity From
Supplemental Material Predicting transfer RNA gene activity from sequence and genome context 1 1,2 1 1 Bryan P. Thornlow , Joel Armstrong , Andrew D. Holmes , Jonathan M. Howard , Russell B. 1,2 1,2 Corbett-Detig , Todd M. Lowe 1 2 D epartment of Biomolecular Engineering, University of California, Santa Cruz, CA 95064; G enomics Institute, University of California, Santa Cruz, CA 95064 Supplemental Material Supplemental Text Supplemental Methods Supplemental Figure S1 Supplemental Figure S2 Supplemental Figure S3 Supplemental Figure S4 Supplemental Figure S5 Supplemental Figure S6 Supplemental Figure S7 Supplemental Figure S8 Supplemental Figure S9 Supplemental Figure S10 Supplemental Figure S11 Supplemental References SUPPLEMENTAL TEXT: Possible tRNA gene classification and assembly errors. Ideally, we would observe that all species contained at least one tRNA gene per expected anticodon (excluding anticodons for which no active tRNA gene is observed in human or mouse). We found three exceptions to this rule. C. hircus (goat) is the only species without a predicted active tRNA-Ser-TGA gene. C. hircus is also one of only three species that has any tRNA-Ser-GGA genes, and it has 27 such genes, one of which was predicted active. These 27 tRNA-Ser-GGA genes do exhibit some sequence similarity, as 21 of these share a 10-base 3’ flanking motif (AGAAGGAAAT). They may therefore reflect a very recent proliferation event. However, it is more likely that they represent assembly errors, because we have previously shown that the immediate 3’ flanking regions of tRNA genes experience high mutation rates (Thornlow et al. 2018), these 27 tRNA genes are found on 15 different chromosomes, and Ser-GGA is not a standard eukaryotic tRNA (Grosjean et al. -
(12) United States Patent (10) Patent No.: US 7.058,515 B1 Selifonov Et Al
US007058515B1 (12) United States Patent (10) Patent No.: US 7.058,515 B1 Selifonov et al. (45) Date of Patent: *Jun. 6, 2006 (54) METHODS FOR MAKING CHARACTER (56) References Cited STRINGS, POLYNUCLEOTIDES AND POLYPEPTIDES HAVING DESRED CHARACTERISTICS U.S. PATENT DOCUMENTS 4.959,312 A 9, 1990 Sirokin (75) Inventors: Sergey A. Selifonov, Los Altos, CA 4,994,379 A 2/1991 Chang (US); Willem P. C. Stemmer, Los 5,023,171 A 6, 1991 Ho et al. ....................... 435/6 Gatos, CA (US); Claes Gustafsson, 5,043,272 A 8/1991 Hartly et al. Belmont, CA (US); Matthew Tobin, 5,066,584 A 11/1991 Gyllensten et al. San Jose, CA (US); Stephen del 5,264,563 A 11/1993 Huse Cardayre, Belmont, CA (US); Phillip SE A Rhen A. Patten, Menlo Park, CA (US); 3.) A 96 MTS al. Jeremy Minshull, Menlo Park, CA 5.512.463 A 4/1996 Stemmer ................... 435,912 (US); Lorraine J. Giver, Sunnyvale, 5.514.588 A 5/1996 Varadaraj.................... 435/262 CA (US) 5,519,319 A 5/1996 Smith et al. 5,521,077 A 5/1996 Kholsa et al. (73) Assignee: Maxygen, Inc., Redwood City, CA 5,565,350 A 10, 1996 Kmiec (US) 5,580,730 A 12/1996 Okamoto 5,605,793 A 2/1997 Stemmer ....................... 435/6 (*) Notice: Subject to any disclaimer, the term of this 5,612.205 A 3/1997 Kay et al. patent is extended or adjusted under 35 3. 6. 3. Site U.S.C. 154(b) by 0 days. 5,717,085w - A 2f1998 Little et al.a This patent is Subject to a terminal dis- 29, Stil 435/6 claimer. -
Ribotoolkit: an Integrated Platform for Analysis and Annotation of Ribosome Profiling Data to Decode Mrna Translation at Codon R
W218–W229 Nucleic Acids Research, 2020, Vol. 48, Web Server issue Published online 19 May 2020 doi: 10.1093/nar/gkaa395 RiboToolkit: an integrated platform for analysis and annotation of ribosome profiling data to decode mRNA translation at codon resolution Qi Liu1,2,*, Tanya Shvarts3, Piotr Sliz2,3 and Richard I. Gregory 1,2,4,5,6,* 1Stem Cell Program, Division of Hematology/Oncology, Boston Children’s Hospital, Boston, MA 02115, USA, 2Department of Biological Chemistry and Molecular Pharmacology, Harvard Medical School, Boston, MA 02115, USA, 3Computational Health Informatics Program, Boston Children’s Hospital, Boston, MA 02115, USA, 4Department of Pediatrics, Harvard Medical School, Boston, MA 02115, USA, 5Harvard Initiative for RNA Medicine, Boston, MA 02115, USA and 6Harvard Stem Cell Institute, Cambridge, MA 02138, USA Received March 06, 2020; Revised April 23, 2020; Editorial Decision May 02, 2020; Accepted May 15, 2020 ABSTRACT INTRODUCTION Ribosome profiling (Ribo-seq) is a powerful technol- Ribosome profiling (Ribo-seq), also known as ribosome ogy for globally monitoring RNA translation; ranging footprinting, has revolutionized the ‘translatomics’ field from codon occupancy profiling, identification of ac- by mapping the position of ribosome-protected fragments tively translated open reading frames (ORFs), to the (RPFs), which typically range in length from 26 to 34 nu- quantification of translational efficiency under vari- cleotides (nt), over the entire transcriptome (1,2). The sci- entific community has employed Ribo-seq to answer awide ous physiological or experimental conditions. How- range of questions, ranging from the identification of trans- ever, analyzing and decoding translation information lated open reading frames (ORFs) to the quantification from Ribo-seq data is not trivial. -
Nanoparticle‐Mediated DNA and Mrna Delivery
REVIEW www.advhealthmat.de Next-Generation Vaccines: Nanoparticle-Mediated DNA and mRNA Delivery William Ho, Mingzhu Gao, Fengqiao Li, Zhongyu Li, Xue-Qing Zhang,* and Xiaoyang Xu* organism-based vaccines have wiped out or Nucleic acid vaccines are a method of immunization aiming to elicit immune nearly eradicated many once great killers responses akin to live attenuated vaccines. In this method, DNA or messenger of humanity, including smallpox, polio, RNA (mRNA) sequences are delivered to the body to generate proteins, which measles, mumps, rubella, diphtheria, per- [1–4] mimic disease antigens to stimulate the immune response. Advantages of tussis, and tetanus. However, the quick emergence of diseases such as SARS- nucleic acid vaccines include stimulation of both cell-mediated and humoral CoV-2, H1N1 as well as quickly evolving immunity, ease of design, rapid adaptability to changing pathogen strains, and deadly diseases like Ebola create a challenge customizable multiantigen vaccines. To combat the SARS-CoV-2 pandemic, for conventional vaccines such as live at- and many other diseases, nucleic acid vaccines appear to be a promising tenuated viral vaccines (LAV), and inacti- [5,6] method. However, aid is needed in delivering the fragile DNA/mRNA payload. vated/killed viral vaccines, which with Many delivery strategies have been developed to elicit effective immune the traditional vaccine development path- way may take on average over 10 years to stimulation, yet no nucleic acid vaccine has been FDA-approved for human develop, or with Ebola requiring an acceler- use. Nanoparticles (NPs) are one of the top candidates to mediate successful ated 5-year development,[7] and even more DNA/mRNA vaccine delivery due to their unique properties, including time needed to scale up manufacturing and unlimited possibilities for formulations, protective capacity, simultaneous stockpile for a large country. -
A Rule-Based Modelling and Simulation Tool for DNA Strand Displacement Systems
This is an electronic reprint of the original article. This reprint may differ from the original in pagination and typographic detail. Gautam, Vinay; Long, Shiting; Orponen, Pekka RuleDSD: A rule-based modelling and simulation tool for DNA strand displacement systems Published in: Proceedings of the 13th International Joint Conference on Biomedical Engineering Systems and Technologies DOI: 10.5220/0008979101580167 Published: 01/01/2020 Document Version Peer reviewed version Published under the following license: CC BY-NC-ND Please cite the original version: Gautam, V., Long, S., & Orponen, P. (2020). RuleDSD: A rule-based modelling and simulation tool for DNA strand displacement systems. In E. De Maria, A. Fred, & H. Gamboa (Eds.), Proceedings of the 13th International Joint Conference on Biomedical Engineering Systems and Technologies: Volume 3 BIOINFORMATICS (pp. 158-167) https://doi.org/10.5220/0008979101580167 This material is protected by copyright and other intellectual property rights, and duplication or sale of all or part of any of the repository collections is not permitted, except that material may be duplicated by you for your research use or educational purposes in electronic or print form. You must obtain permission for any other use. Electronic or print copies may not be offered, whether for sale or otherwise to anyone who is not an authorised user. Powered by TCPDF (www.tcpdf.org) RuleDSD: A Rule-based Modelling and Simulation Tool for DNA Strand Displacement Systems Vinay Gautam a, Shiting Long and Pekka Orponen b Department of Computer Science, Aalto University, 00076 Aalto, Finland fvinay.gautam, shiting.long, pekka.orponeng@aalto.fi Keywords: DNA Strand Displacement, DSD Modelling and Simulation, PySB, Rule-based Model, BioNetGen. -
Comparative Ribosome Profiling Uncovers a Dominant Role for Translational Control in Toxoplasma Gondii
Comparative ribosome profiling uncovers a dominant role for translational control in Toxoplasma gondii The MIT Faculty has made this article openly available. Please share how this access benefits you. Your story matters. Citation Hassan, Musa A. et al. "Comparative ribosome profiling uncovers a dominant role for translational control in Toxoplasma gondii." BMC Genomics 18 (December 2017): 961 © 2017 The Author(s) As Published http://dx.doi.org/10.1186/s12864-017-4362-6 Publisher BioMed Central Version Final published version Citable link http://hdl.handle.net/1721.1/114431 Terms of Use Creative Commons Attribution Detailed Terms http://creativecommons.org/licenses/by/4.0/ Hassan et al. BMC Genomics (2017) 18:961 DOI 10.1186/s12864-017-4362-6 RESEARCH ARTICLE Open Access Comparative ribosome profiling uncovers a dominant role for translational control in Toxoplasma gondii Musa A. Hassan1,2*, Juan J. Vasquez3,8, Chew Guo-Liang4, Markus Meissner5,6 and T. Nicolai Siegel3,6,7 Abstract Background: The lytic cycle of the protozoan parasite Toxoplasma gondii, which involves a brief sojourn in the extracellular space, is characterized by defined transcriptional profiles. For an obligate intracellular parasite that is shielded from the cytosolic host immune factors by a parasitophorous vacuole, the brief entry into the extracellular space is likely to exert enormous stress. Due to its role in cellular stress response, we hypothesize that translational control plays an important role in regulating gene expression in Toxoplasma during the lytic cycle. Unlike transcriptional profiles, insights into genome-wide translational profiles of Toxoplasma gondii are lacking. Methods: We have performed genome-wide ribosome profiling, coupled with high throughput RNA sequencing, in intracellular and extracellular Toxoplasma gondii parasites to investigate translational control during the lytic cycle. -
Design of Nucleic Acid Sequences for DNA Computing Based on a Thermodynamic Approach
Title Design of nucleic acid sequences for DNA computing based on a thermodynamic approach Author(s) Tanaka, Fumiaki; Kameda, Atsushi; Yamamoto, Masahito; Ohuchi, Azuma Nucleic Acids Research, 33(3), 903-911 Citation https://doi.org/10.1093/nar/gki235 Issue Date 2005 Doc URL http://hdl.handle.net/2115/64504 Type article File Information gki235.pdf Instructions for use Hokkaido University Collection of Scholarly and Academic Papers : HUSCAP Published online February 8, 2005 Nucleic Acids Research, 2005, Vol. 33, No. 3 903–911 doi:10.1093/nar/gki235 Design of nucleic acid sequences for DNA computing based on a thermodynamic approach Fumiaki Tanaka1,*, Atsushi Kameda2, Masahito Yamamoto1,2 and Azuma Ohuchi1,2 1Graduate School of Engineering, Hokkaido University, North 13, West 8, Kita-ku, Sapporo 060-8628, Japan and 2CREST, Japan Science and Technology Corporation, 4-1-8, Honmachi, Kawaguchi, Saitama, 332-0012, Japan Received October 19, 2004; Revised and Accepted January 20, 2005 ABSTRACT order to satisfy the constraints based on the physicochemical properties of nucleic acids. In particular, it is essential to We have developed an algorithm for designing prevent undesired hybridization. In DNA computing, multiple multiple sequences of nucleic acids that have a sequences need to be designed that do not hybridize non- uniform melting temperature between the sequence specifically with each other (10,11), while in RNA secondary and its complement and that do not hybridize non- structure design, a single sequence needs to be designed that specifically with each other based on the minimum folds into the desired secondary structure (12–15). For free energy (DGmin).