1
|
Rissone P, Severino A, Pastor I, Ritort F. Universal cold RNA phase transitions. Proc Natl Acad Sci U S A 2024; 121:e2408313121. [PMID: 39150781 DOI: 10.1073/pnas.2408313121] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/26/2024] [Accepted: 07/12/2024] [Indexed: 08/18/2024] Open
Abstract
RNA's diversity of structures and functions impacts all life forms since primordia. We use calorimetric force spectroscopy to investigate RNA folding landscapes in previously unexplored low-temperature conditions. We find that Watson-Crick RNA hairpins, the most basic secondary structure elements, undergo a glass-like transition below [Formula: see text]C where the heat capacity abruptly changes and the RNA folds into a diversity of misfolded structures. We hypothesize that an altered RNA biochemistry, determined by sequence-independent ribose-water interactions, outweighs sequence-dependent base pairing. The ubiquitous ribose-water interactions lead to universal RNA phase transitions below TG, such as maximum stability at [Formula: see text]C where water density is maximum, and cold denaturation at [Formula: see text]C. RNA cold biochemistry may have a profound impact on RNA function and evolution.
Collapse
Affiliation(s)
- Paolo Rissone
- Small Biosystems Lab, Condensed Matter Physics Department, Universitat de Barcelona, Barcelona 08028, Spain
| | - Aurélien Severino
- Small Biosystems Lab, Condensed Matter Physics Department, Universitat de Barcelona, Barcelona 08028, Spain
| | - Isabel Pastor
- Small Biosystems Lab, Condensed Matter Physics Department, Universitat de Barcelona, Barcelona 08028, Spain
| | - Felix Ritort
- Small Biosystems Lab, Condensed Matter Physics Department, Universitat de Barcelona, Barcelona 08028, Spain
- Institut de Nanociència i Nanotecnologia, Universitat de Barcelona, Barcelona 08028, Spain
| |
Collapse
|
2
|
Kim J, Seo M, Lim Y, Kim J. START: A Versatile Platform for Bacterial Ligand Sensing with Programmable Performances. ADVANCED SCIENCE (WEINHEIM, BADEN-WURTTEMBERG, GERMANY) 2024:e2402029. [PMID: 39075726 DOI: 10.1002/advs.202402029] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/26/2024] [Revised: 05/31/2024] [Indexed: 07/31/2024]
Abstract
Recognition of signaling molecules for coordinated regulation of target genes is a fundamental process for biological systems. Cells often rely on transcription factors to accomplish these intricate tasks, yet the subtle conformational changes of protein structures, coupled with the complexity of intertwined protein interaction networks, pose challenges for repurposing these for bioengineering applications. This study introduces a novel platform for ligand-responsive gene regulation, termed START (Synthetic Trans-Acting Riboswitch with Triggering RNA). Inspired by the bacterial ligand sensing system, riboswitch, and the synthetic gene regulator, toehold switch, the START platform enables the implementation of synthetic biosensors for various ligands. Rational sequence design with targeted domain optimization yields high-performance STARTs with a dynamic range up to 67.29-fold and a tunable ligand sensitivity, providing a simple and intuitive strategy for sensor engineering. The START platform also exhibits modularity and composability to allow flexible genetic circuit construction, enabling seamless implementation of OR, AND, and NOT Boolean logic gates for multiple ligand inputs. The START design principle is capable of broadening the suite of synthetic biosensors for diverse chemical and protein ligands, providing a novel riboregulator chassis for synthetic biology and bioengineering applications.
Collapse
Affiliation(s)
- Jeongwon Kim
- Department of Life Sciences, Pohang University of Science and Technology, Pohang, 37673, South Korea
| | - Minchae Seo
- Department of Life Sciences, Pohang University of Science and Technology, Pohang, 37673, South Korea
| | - Yelin Lim
- Department of Life Sciences, Pohang University of Science and Technology, Pohang, 37673, South Korea
| | - Jongmin Kim
- Department of Life Sciences, Pohang University of Science and Technology, Pohang, 37673, South Korea
| |
Collapse
|
3
|
Liang G, Zhao J, Gao Y, Xie T, Zhen J, Pan L, Gong W. Application and evaluation of molecular docking for aptamer and small molecular interaction - A case study with tetracycline antibiotics. Talanta 2024; 266:124942. [PMID: 37459789 DOI: 10.1016/j.talanta.2023.124942] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/20/2023] [Revised: 06/28/2023] [Accepted: 07/12/2023] [Indexed: 09/20/2023]
Abstract
Molecular docking (MD) analysis is currently the most commonly used theoretical simulation method to investigate the interaction of aptamers (receptors) and small molecules (ligands) and understand the recognition mechanism between them at a molecular level. Using the specific aptamers of tetracycline antibiotics (tetracycline (TET), oxytetracycline (OTC), doxycycline (DOC)) as the docking models, three steady-state aptamers of tertiary structures (SATS) were established for each aptamer with the UNAFold and RNAComposer tools. The binding free energy (BFE), docking score (DS), and binding site (base) of the specific ligands (TET, OTC, and DOC) with their respective SATS were obtained by molecular docking. The results revealed one or more binding sites in the established SATS of the aptamers. The BFE and DS of different binding sites of one specific SATS varied significantly. The results also revealed that the site with the highest BFE represented the most dominant binding site, even if it was not the SATS with minimum energy. The BFE values could also be used to evaluate the affinity and specificity of the aptamer to its target. For the first time, this study proposes a method for MD analysis of the aptamer and its target based on different SATS, clarification of the binding mode, and prediction of the binding sites (bases). This study provides a theoretical basis for tailoring; structural optimization; and base modification of aptamers; identifying aptamers with high affinity and specificity.
Collapse
Affiliation(s)
- Gang Liang
- Institute of Quality Standard and Testing Technology, BAAFS (Beijing Academy of Agriculture and Forestry Sciences), Beijing, 100097, China; Risk Assessment Lab for Agro-products (Beijing), Ministry of Agriculture, Beijing, 100097, China.
| | - Jie Zhao
- Institute of Quality Standard and Testing Technology, BAAFS (Beijing Academy of Agriculture and Forestry Sciences), Beijing, 100097, China; Risk Assessment Lab for Agro-products (Beijing), Ministry of Agriculture, Beijing, 100097, China
| | - Yufei Gao
- College of Environmental Science and Engineering, Hebei University of Science and Technology, Shijiazhuang, 050024, China
| | - Tao Xie
- Chengdu University, 610000, Chengdu, China
| | - Jianhui Zhen
- Shijiazhuang Customs Technology Center PR China, Shijiazhuang, 050051, China.
| | - Ligang Pan
- Institute of Quality Standard and Testing Technology, BAAFS (Beijing Academy of Agriculture and Forestry Sciences), Beijing, 100097, China; Risk Assessment Lab for Agro-products (Beijing), Ministry of Agriculture, Beijing, 100097, China
| | - Wenwen Gong
- Institute of Quality Standard and Testing Technology, BAAFS (Beijing Academy of Agriculture and Forestry Sciences), Beijing, 100097, China; Risk Assessment Lab for Agro-products (Beijing), Ministry of Agriculture, Beijing, 100097, China.
| |
Collapse
|
4
|
Wang X, Tan YL, Yu S, Shi YZ, Tan ZJ. Predicting 3D structures and stabilities for complex RNA pseudoknots in ion solutions. Biophys J 2023; 122:1503-1516. [PMID: 36924021 PMCID: PMC10147842 DOI: 10.1016/j.bpj.2023.03.017] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/24/2022] [Revised: 03/05/2023] [Accepted: 03/10/2023] [Indexed: 03/17/2023] Open
Abstract
RNA pseudoknots are a kind of important tertiary motif, and the structures and stabilities of pseudoknots are generally critical to the biological functions of RNAs with the motifs. In this work, we have carefully refined our previously developed coarse-grained model with salt effect through involving a new coarse-grained force field and a replica-exchange Monte Carlo algorithm, and employed the model to predict structures and stabilities of complex RNA pseudoknots in ion solutions beyond minimal H-type pseudoknots. Compared with available experimental data, the newly refined model can successfully predict 3D structures from sequences for the complex RNA pseudoknots including SARS-CoV-2 programming-1 ribosomal frameshifting element and Zika virus xrRNA, and can reliably predict the thermal stabilities of RNA pseudoknots with various sequences and lengths over broad ranges of monovalent/divalent salts. In addition, for complex pseudoknots including SARS-CoV-2 frameshifting element, our analyses show that their thermally unfolding pathways are mainly dependent on the relative stabilities of unfolded intermediate states, in analogy to those of minimal H-type pseudoknots.
Collapse
Affiliation(s)
- Xunxun Wang
- Department of Physics and Key Laboratory of Artificial Micro & Nano-structures of Ministry of Education, School of Physics and Technology, Wuhan University, Wuhan, China
| | - Ya-Lan Tan
- Research Center of Nonlinear Science and School of Mathematical and Physical Sciences, Wuhan Textile University, Wuhan, China
| | - Shixiong Yu
- Department of Physics and Key Laboratory of Artificial Micro & Nano-structures of Ministry of Education, School of Physics and Technology, Wuhan University, Wuhan, China
| | - Ya-Zhou Shi
- Research Center of Nonlinear Science and School of Mathematical and Physical Sciences, Wuhan Textile University, Wuhan, China
| | - Zhi-Jie Tan
- Department of Physics and Key Laboratory of Artificial Micro & Nano-structures of Ministry of Education, School of Physics and Technology, Wuhan University, Wuhan, China.
| |
Collapse
|
5
|
Maksudov F, Kliuchnikov E, Pierson D, Ujwal M, Marx KA, Chanda A, Barsegov V. Therapeutic phosphorodiamidate morpholino oligonucleotides: Physical properties, solution structures, and folding thermodynamics. MOLECULAR THERAPY. NUCLEIC ACIDS 2023; 31:631-647. [PMID: 36910708 PMCID: PMC9996446 DOI: 10.1016/j.omtn.2023.02.007] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/27/2022] [Accepted: 02/10/2023] [Indexed: 02/16/2023]
Abstract
Elucidating the structure-function relationships for therapeutic RNA mimicking phosphorodiamidate morpholino oligonucleotides (PMOs) is challenging due to the lack of information about their structures. While PMOs have been approved by the US Food and Drug Administration for treatment of Duchenne muscular dystrophy, no structural information on these unique, charge-neutral, and stable molecules is available. We performed circular dichroism and solution viscosity measurements combined with molecular dynamics simulations and machine learning to resolve solution structures of 22-mer, 25-mer, and 30-mer length PMOs. The PMO conformational dynamics are defined by the competition between non-polar nucleobases and uncharged phosphorodiamidate groups for shielding from solvent exposure. PMO molecules form non-canonical, partially helical, stable folded structures with a small 1.4- to 1.7-nm radius of gyration, low count of three to six base pairs and six to nine base stacks, characterized by -34 to -51 kcal/mol free energy, -57 to -103 kcal/mol enthalpy, and -23 to -53 kcal/mol entropy for folding. The 4.5- to 6.2-cm3/g intrinsic viscosity and Huggins constant of 4.5-9.9 are indicative of extended and aggregating systems. The results obtained highlight the importance of the conformational ensemble view of PMO solution structures, thermodynamic stability of their non-canonical structures, and concentration-dependent viscosity properties. These principles form a paradigm to understand the structure-properties-function relationship for therapeutic PMOs to advance the design of new RNA-mimic-based drugs.
Collapse
Affiliation(s)
- Farkhad Maksudov
- Department of Chemistry, University of Massachusetts, Lowell, MA 01854, USA
| | | | - Daniel Pierson
- Technical Operations, Sarepta Therapeutics, Cambridge, MA 02142, USA
| | | | - Kenneth A. Marx
- Department of Chemistry, University of Massachusetts, Lowell, MA 01854, USA
- Inciton, Inc., Andover, MA 01854, USA
| | - Arani Chanda
- Technical Operations, Sarepta Therapeutics, Cambridge, MA 02142, USA
- Corresponding author: Arani Chanda, Technical Operations, Sarepta Therapeutics, Cambridge, MA 02142, USA.
| | - Valeri Barsegov
- Department of Chemistry, University of Massachusetts, Lowell, MA 01854, USA
- Inciton, Inc., Andover, MA 01854, USA
- Corresponding author: Valeri Barsegov, Department of Chemistry, University of Massachusetts, Lowell, MA 01854, USA.
| |
Collapse
|
6
|
Watkins AM, Das R. RNA 3D Modeling with FARFAR2, Online. Methods Mol Biol 2023; 2586:233-249. [PMID: 36705908 DOI: 10.1007/978-1-0716-2768-6_14] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/28/2023]
Abstract
Understanding the three-dimensional structure of an RNA molecule is often essential to understanding its function. Sampling algorithms and energy functions for RNA structure prediction are improving, due to the increasing diversity of structural data available for training statistical potentials and testing structural data, along with a steady supply of blind challenges through the RNA-Puzzles initiative. The recent FARFAR2 algorithm enables near-native structure predictions on fairly complex RNA structures, including automated selection of final candidate models and estimation of model accuracy. Here, we describe the use of a publicly available webserver for RNA modeling for realistic scenarios using FARFAR2, available at https://rosie.rosettacommons.org/farfar2 . We walk through two cases in some detail: a simple model pseudoknot from the frameshifting element of beet western yellows virus modeled using the "basic interface" to the webserver and a replication of RNA-Puzzle 20, a metagenomic twister sister ribozyme, using the "advanced interface." We also describe example runs of FARFAR2 modeling including two kinds of experimental data: a c-di-GMP riboswitch modeled with low-resolution restraints from MOHCA-seq experiments and a tandem GA motif modeled with 1H NMR chemical shifts.
Collapse
Affiliation(s)
- Andrew M Watkins
- Department of Biochemistry, Stanford University School of Medicine, Stanford, CA, USA
- Prescient Design, Genentech, South San Francisco, CA, USA
| | - Rhiju Das
- Department of Biochemistry, Stanford University School of Medicine, Stanford, CA, USA.
- Biophysics Program, Stanford University, Stanford, CA, USA.
| |
Collapse
|
7
|
Wang J, Sha CM, Dokholyan NV. Combining Experimental Restraints and RNA 3D Structure Prediction in RNA Nanotechnology. Methods Mol Biol 2023; 2709:51-64. [PMID: 37572272 PMCID: PMC10680996 DOI: 10.1007/978-1-0716-3417-2_3] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/14/2023]
Abstract
Precise RNA tertiary structure prediction can aid in the design of RNA nanoparticles. However, most existing RNA tertiary structure prediction methods are limited to small RNAs with relatively simple secondary structures. Large RNA molecules usually have complex secondary structures, including multibranched loops and pseudoknots, allowing for highly flexible RNA geometries and multiple stable states. Various experiments and bioinformatics analyses can often provide information about the distance between atoms (or residues) in RNA, which can be used to guide the prediction of RNA tertiary structure. In this chapter, we will introduce a platform, iFoldNMR, that can incorporate non-exchangeable imino protons resonance data from NMR as restraints for RNA 3D structure prediction. We also introduce an algorithm, DVASS, which optimizes distance restraints for better RNA 3D structure prediction.
Collapse
Affiliation(s)
- Jian Wang
- Department of Pharmacology, Penn State College of Medicine, Hershey, PA, USA
| | - Congzhou M Sha
- Department of Pharmacology, Penn State College of Medicine, Hershey, PA, USA
- Department of Engineering Science and Mechanics, Penn State University, State College, PA, USA
| | - Nikolay V Dokholyan
- Department of Pharmacology, Penn State College of Medicine, Hershey, PA, USA.
- Department of Engineering Science and Mechanics, Penn State University, State College, PA, USA.
- Department of Biochemistry and Molecular Biology, Penn State College of Medicine, Hershey, PA, USA.
- Department of Chemistry, Penn State University, State College, PA, USA.
- Department of Biomedical Engineering, Penn State University, State College, PA, USA.
| |
Collapse
|
8
|
Yunussova N, Sypabekova M, Zhumabekova Z, Matkarimov B, Kanayeva D. A Novel ssDNA Aptamer Targeting Carcinoembryonic Antigen: Selection and Characterization. BIOLOGY 2022; 11:biology11101540. [PMID: 36290442 PMCID: PMC9598387 DOI: 10.3390/biology11101540] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 09/18/2022] [Revised: 10/15/2022] [Accepted: 10/17/2022] [Indexed: 11/20/2022]
Abstract
One of the major causes of a drastically shorter life expectancy and one of the most prevalent diseases in the world today is cancer. Given the data on the rise in cancer cases throughout the world, it is obvious that, despite the diagnostic techniques currently being used, there is a pressing need to develop precise and sensitive techniques for early diagnosis of the disease. A high degree of affinity and specificity towards particular targets is maintained by the short nucleic acid molecules known as aptamers. Aptamers outperform antibodies due to their unique benefits, such as their simplicity in synthesis and modification, lack of toxicity, and long-term stability. Utilizing an accurate recognition element and a robust signal transduction mechanism, molecular diagnostics can be extremely sensitive and specific. In this study, development of new single-stranded DNA aptamers against CEA for use in cancer diagnostics was accomplished using SELEX and NGS methods. As a result of 12 iterative SELEX rounds, nine aptamer candidates against CEA were developed. NGS comparative analysis revealed that round twelve had an enriched number of aptamers that were specifically bound, as opposed to round eight. Among the selected nine sequences characterized by bioinformatics analysis and ELONA, an aptamer sequence with the highest specificity and affinity for the target protein was identified and further examined. Aptamer sequence (6) was screened in a concentration-dependent assay, specificity analysis was performed, and its potential secondary and tertiary structures were predicted, which enabled us to test one of the possible putative interactions with CEA. Finally, aptamer sequence (6) labelled with a Cy5 fluorescent tag was used in confocal microscopy to observe its binding towards the CEA expressed in HT-29 human colon adenocarcinoma cell line.
Collapse
Affiliation(s)
- Nigara Yunussova
- Ph.D. Program in Life Sciences, Nazarbayev University, 53 Kabanbay Batyr Ave., Astana 010000, Kazakhstan
| | - Marzhan Sypabekova
- National Laboratory Astana, Nazarbayev University, 53 Kabanbay Batyr Ave., Astana 010000, Kazakhstan
| | - Zhazira Zhumabekova
- M.Sc. Program in Biological Sciences, Nazarbayev University, 53 Kabanbay Batyr Ave., Astana 010000, Kazakhstan
| | - Bakhyt Matkarimov
- National Laboratory Astana, Nazarbayev University, 53 Kabanbay Batyr Ave., Astana 010000, Kazakhstan
| | - Damira Kanayeva
- Department of Biology, School of Sciences and Humanities, Nazarbayev University, 53 Kabanbay Batyr Ave., Astana 010000, Kazakhstan
- Correspondence:
| |
Collapse
|
9
|
Zhou L, Wang X, Yu S, Tan YL, Tan ZJ. FebRNA: An automated fragment-ensemble-based model for building RNA 3D structures. Biophys J 2022; 121:3381-3392. [PMID: 35978551 PMCID: PMC9515226 DOI: 10.1016/j.bpj.2022.08.017] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/30/2022] [Revised: 07/19/2022] [Accepted: 08/15/2022] [Indexed: 11/23/2022] Open
Abstract
Knowledge of RNA three-dimensional (3D) structures is critical to understanding the important biological functions of RNAs. Although various structure prediction models have been developed, the high-accuracy predictions of RNA 3D structures are still limited to the RNAs with short lengths or with simple topology. In this work, we proposed a new model, namely FebRNA, for building RNA 3D structures through fragment assembly based on coarse-grained (CG) fragment ensembles. Specifically, FebRNA is composed of four processes: establishing the library of different types of non-redundant CG fragment ensembles regardless of the sequences, building CG 3D structure ensemble through fragment assembly, identifying top-scored CG structures through a specific CG scoring function, and rebuilding the all-atom structures from the top-scored CG ones. Extensive examination against different types of RNA structures indicates that FebRNA consistently gives the reliable predictions on RNA 3D structures, including pseudoknots, three-way junctions, four-way and five-way junctions, and RNAs in the RNA-Puzzles. FebRNA is available on the Web site: https://github.com/Tan-group/FebRNA.
Collapse
Affiliation(s)
- Li Zhou
- Department of Physics and Key Laboratory of Artificial Micro & Nano-structures of Education, School of Physics and Technology, Wuhan University, Wuhan 430072, China
| | - Xunxun Wang
- Department of Physics and Key Laboratory of Artificial Micro & Nano-structures of Education, School of Physics and Technology, Wuhan University, Wuhan 430072, China
| | - Shixiong Yu
- Department of Physics and Key Laboratory of Artificial Micro & Nano-structures of Education, School of Physics and Technology, Wuhan University, Wuhan 430072, China
| | - Ya-Lan Tan
- Research Center of Nonlinear Science, School of Mathematical and Physical Sciences, Wuhan Textile University, Wuhan 430073, China.
| | - Zhi-Jie Tan
- Department of Physics and Key Laboratory of Artificial Micro & Nano-structures of Education, School of Physics and Technology, Wuhan University, Wuhan 430072, China.
| |
Collapse
|
10
|
Matarrese MAG, Loppini A, Nicoletti M, Filippi S, Chiodo L. Assessment of tools for RNA secondary structure prediction and extraction: a final-user perspective. J Biomol Struct Dyn 2022:1-20. [DOI: 10.1080/07391102.2022.2116110] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/14/2022]
Affiliation(s)
- Margherita A. G. Matarrese
- Engineering Department, Campus Bio-Medico University of Rome, Rome, Italy
- Jane and John Justin Neurosciences Center, Cook Children’s Health Care System, TX, USA
- Department of Bioengineering, The University of Texas at Arlington, Arlington, TX, USA
| | - Alessandro Loppini
- Engineering Department, Campus Bio-Medico University of Rome, Rome, Italy
- Center for Life Nano & Neuroscience, Italian Institute of Technology, Rome, Italy
| | - Martina Nicoletti
- Engineering Department, Campus Bio-Medico University of Rome, Rome, Italy
- Center for Life Nano & Neuroscience, Italian Institute of Technology, Rome, Italy
| | - Simonetta Filippi
- Engineering Department, Campus Bio-Medico University of Rome, Rome, Italy
| | - Letizia Chiodo
- Engineering Department, Campus Bio-Medico University of Rome, Rome, Italy
| |
Collapse
|
11
|
Wang T, Simmel FC. Riboswitch-inspired toehold riboregulators for gene regulation in Escherichia coli. Nucleic Acids Res 2022; 50:4784-4798. [PMID: 35446427 PMCID: PMC9071393 DOI: 10.1093/nar/gkac275] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/05/2022] [Accepted: 04/08/2022] [Indexed: 12/24/2022] Open
Abstract
Regulatory RNA molecules have been widely investigated as components for synthetic gene circuits, complementing the use of protein-based transcription factors. Among the potential advantages of RNA-based gene regulators are their comparatively simple design, sequence-programmability, orthogonality, and their relatively low metabolic burden. In this work, we developed a set of riboswitch-inspired riboregulators in Escherichia coli that combine the concept of toehold-mediated strand displacement (TMSD) with the switching principles of naturally occurring transcriptional and translational riboswitches. Specifically, for translational activation and repression, we sequestered anti-anti-RBS or anti-RBS sequences, respectively, inside the loop of a stable hairpin domain, which is equipped with a single-stranded toehold region at its 5' end and is followed by regulated sequences on its 3' side. A trigger RNA binding to the toehold region can invade the hairpin, inducing a structural rearrangement that results in translational activation or deactivation. We also demonstrate that TMSD can be applied in the context of transcriptional regulation by switching RNA secondary structure involved in Rho-dependent termination. Our designs expand the repertoire of available synthetic riboregulators by a set of RNA switches with no sequence limitation, which should prove useful for the development of robust genetic sensors and circuits.
Collapse
Affiliation(s)
- Tianhe Wang
- Physics of Synthetic Biological Systems – E14, Physics Department and ZNN, Technische Universität München, Am Coulombwall 4a, 85748 Garching, Germany
| | | |
Collapse
|
12
|
Sousa DA, Carneiro M, Ferreira D, Moreira FTC, Sales MGFV, Rodrigues LR. Recent advances in the selection of cancer-specific aptamers for the development of biosensors. Curr Med Chem 2022; 29:5850-5880. [PMID: 35209816 DOI: 10.2174/0929867329666220224155037] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/29/2021] [Revised: 11/30/2021] [Accepted: 12/19/2021] [Indexed: 11/22/2022]
Abstract
An early diagnosis has the potential to greatly decrease cancer mortality. For that purpose, specific cancer biomarkers have been molecularly targeted by aptamer sequences to enable an accurate and rapid detection. Aptamer-based biosensors for cancer diagnostics are a promising alternative to those using antibodies, due to their high affinity and specificity to the target molecules and advantageous production. Synthetic nucleic acid aptamers are generated by in vitro Systematic Evolution of Ligands by Exponential enrichment (SELEX) methodologies that have been improved over the years to enhance the efficacy and to shorten the selection process. Aptamers have been successfully applied in electrochemical, optical, photoelectrochemical and piezoelectrical-based detection strategies. These aptasensors comprise a sensitive, accurate and inexpensive option for cancer detection being used as point-of-care devices. This review highlights the recent advances in cancer biomarkers, achievements and optimizations made in aptamer selection, as well as the different aptasensors developed for the detection of several cancer biomarkers.
Collapse
Affiliation(s)
- Diana A Sousa
- CEB- Centre of Biological Engineering, University of Minho, Campus de Gualtar, Braga, Portugal
- MIT-Portugal Program, Lisbon, Portugal
| | - Mariana Carneiro
- CEB- Centre of Biological Engineering, University of Minho, Campus de Gualtar, Braga, Portugal
- BioMark@ISEP, School of Engineering, Polytechnic of Porto, Porto, Portugal
| | - Débora Ferreira
- CEB- Centre of Biological Engineering, University of Minho, Campus de Gualtar, Braga, Portugal
- MIT-Portugal Program, Lisbon, Portugal
| | - Felismina T C Moreira
- CEB- Centre of Biological Engineering, University of Minho, Campus de Gualtar, Braga, Portugal
- BioMark@ISEP, School of Engineering, Polytechnic of Porto, Porto, Portugal
| | - Maria Goreti F V Sales
- CEB- Centre of Biological Engineering, University of Minho, Campus de Gualtar, Braga, Portugal
- MIT-Portugal Program, Lisbon, Portugal
- BioMark@UC, Faculty of Sciences and Technology, University of Coimbra, Coimbra, Portugal
| | - Lígia R Rodrigues
- CEB- Centre of Biological Engineering, University of Minho, Campus de Gualtar, Braga, Portugal
| |
Collapse
|
13
|
Predicting RNA Scaffolds with a Hybrid Method of Vfold3D and VfoldLA. Methods Mol Biol 2021. [PMID: 34086269 DOI: 10.1007/978-1-0716-1499-0_1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/30/2023]
Abstract
The ever-increasing discoveries of noncoding RNA functions draw a strong demand for RNA structure determination from the sequence. In recently years, computational studies for RNA structures, at both the two-dimensional and the three-dimensional levels, led to several highly promising new developments. In this chapter, we describe a hybrid method, which combines the motif template-based Vfold3D model and the loop template-based VfoldLA model, to predict RNA 3D structures. The main emphasis is placed on the definition of motifs and loops, the treatment of no-template motifs, and the 3D structure assembly from templates of motifs and loops. For illustration, we use the ZIKV xrRNA1 as an example to show the template-based prediction of RNA 3D structures from the 2D structure. The web server for the hybrid model is freely accessible at http://rna.physics.missouri.edu/vfold3D2 .
Collapse
|
14
|
Mak CH, Phan ENH. Diagrammatic approaches to RNA structures with trinucleotide repeats. Biophys J 2021; 120:2343-2354. [PMID: 33887227 PMCID: PMC8390803 DOI: 10.1016/j.bpj.2021.04.010] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/03/2020] [Revised: 04/07/2021] [Accepted: 04/09/2021] [Indexed: 11/30/2022] Open
Abstract
Trinucleotide repeat expansion disorders are associated with the overexpansion of (CNG) repeats on the genome. Messenger RNA transcripts of sequences with greater than 60–100 (CNG) tandem units have been implicated in trinucleotide repeat expansion disorder pathogenesis. In this work, we develop a diagrammatic theory to study the structural diversity of these (CNG)n RNA sequences. Representing structural elements on the chain’s conformation by a set of graphs and employing elementary diagrammatic methods, we have formulated a renormalization procedure to re-sum these graphs and arrive at a closed-form expression for the ensemble partition function. With a simple approximation for the renormalization and applied to extended (CNG)n sequences, this theory can comprehensively capture an infinite set of conformations with any number and any combination of duplexes, hairpins, multiway junctions, and quadruplexes. To quantify the diversity of different (CNG)n ensembles, the analytical equations derived from the diagrammatic theory were solved numerically to derive equilibrium estimates for the secondary structural contents of the chains. The results suggest that the structural ensembles of (CNG)n repeat sequence with n ∼60 are surprisingly diverse, and the distribution is sensitive to the ability of the N nucleotide to make noncanonical pairs and whether the (CNG)n sequence can sustain stable quadruplexes. The results show how perturbations in the form of biases on the stabilities of the various structural motifs, duplexes, junctions, helices, and quadruplexes could affect the secondary structures of the chains and how these structures may switch when they are perturbed.
Collapse
Affiliation(s)
- Chi H Mak
- Department of Chemistry, Center of Applied Mathematical Sciences and Department of Quantitative and Computational Biology, University of Southern California, Los Angeles, California.
| | - Ethan N H Phan
- Department of Chemistry, University of Southern California, Los Angeles, California
| |
Collapse
|
15
|
Nanosized Particles Assembled by a Recombinant Virus Protein Are Able to Encapsulate Negatively Charged Molecules and Structured RNA. Polymers (Basel) 2021; 13:polym13060858. [PMID: 33799623 PMCID: PMC7998283 DOI: 10.3390/polym13060858] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/24/2021] [Revised: 03/07/2021] [Accepted: 03/09/2021] [Indexed: 11/19/2022] Open
Abstract
RNA-based molecules have recently become hot candidates to be developed into therapeutic agents. However, successful applications of RNA-based therapeutics might require suitable carriers to protect the RNA from enzymatic degradation by ubiquitous RNases in vivo. Because of their better biocompatibility and biodegradability, protein-based nanoparticles are considered to be alternatives to their synthetic polymer-based counterparts for drug delivery. Hepatitis C virus (HCV) core protein has been suggested to be able to self-assemble into nucleocapsid-like particles in vitro. In this study, the genomic RNA-binding domain of HCV core protein consisting of 116 amino acids (p116) was overexpressed with E. coli for investigation. The recombinant p116 was able to assemble into particles with an average diameter of approximately 27 nm, as visualized by electron microscopy and atomic force microscopy. Measurements with fluorescence spectroscopy, flow cytometry, and fluorescence quenching indicated that the p116-assembled nanoparticles were able to encapsulate small anionic molecules and structured RNA. This study demonstrates methods that exploit the self-assembly nature of a virus-derived protein for nanoparticle production. This study also suggests that the virus-derived protein-assembled particles could possibly be developed into potential carriers for anionic molecular drugs and structured RNA-based therapeutics.
Collapse
|
16
|
Deep Learning Method for RNA Secondary Structure Prediction with Pseudoknots Based on Large-Scale Data. JOURNAL OF HEALTHCARE ENGINEERING 2021. [DOI: 10.1155/2021/6699996] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
Abstract
Traditional machine learning methods are widely used in the field of RNA secondary structure prediction and have achieved good results. However, with the emergence of large-scale data, deep learning methods have more advantages than traditional machine learning methods. As the number of network layers increases in deep learning, there will often be problems such as increased parameters and overfitting. We used two deep learning models, GoogLeNet and TCN, to predict RNA secondary results. And from the perspective of the depth and width of the network, improvements are made based on the neural network model, which can effectively improve the computational efficiency while extracting more feature information. We process the existing real RNA data through experiments, use deep learning models to extract useful features from a large amount of RNA sequence data and structure data, and then predict the extracted features to obtain each base’s pairing probability. The characteristics of RNA secondary structure and dynamic programming methods are used to process the base prediction results, and the structure with the largest sum of the probability of each base pairing is obtained, and this structure will be used as the optimal RNA secondary structure. We, respectively, evaluated GoogLeNet and TCN models based on 5sRNA, tRNA data, and tmRNA data, and compared them with other standard prediction algorithms. The sensitivity and specificity of the GoogLeNet model on the 5sRNA and tRNA data sets are about 16% higher than the best prediction results in other algorithms. The sensitivity and specificity of the GoogLeNet model on the tmRNA dataset are about 9% higher than the best prediction results in other algorithms. As deep learning algorithms’ performance is related to the size of the data set, as the scale of RNA data continues to expand, the prediction accuracy of deep learning methods for RNA secondary structure will continue to improve.
Collapse
|
17
|
Abstract
Novel RNA motif design is of great practical importance for technology and medicine. Increasingly, computational design plays an important role in such efforts. Our coarse-grained RAG (RNA-As-Graphs) framework offers strategies for enumerating the universe of RNA 2D folds, selecting "RNA-like" candidates for design, and determining sequences that fold onto these candidates. In RAG, RNA secondary structures are represented as tree or dual graphs. Graphs with known RNA structures are called "existing", and the others are labeled "hypothetical". By using simplified features for RNA graphs, we have clustered the hypothetical graphs into "RNA-like" and "non-RNA-like" groups and proposed RNA-like graphs as candidates for design. Here, we propose a new way of designing graph features by using Fiedler vectors. The new features reflect graph shapes better, and they lead to a more clustered organization of existing graphs. We show significant increases in K-means clustering accuracy by using the new features (e.g., up to 95% and 98% accuracy for tree and dual graphs, respectively). In addition, we propose a scoring model for top graph candidate selection. This scoring model allows users to set a threshold for candidates, and it incorporates weighing of existing graphs based on their corresponding number of known RNAs. We include a list of top scored RNA-like candidates, which we hope will stimulate future novel RNA design.
Collapse
Affiliation(s)
- Qiyao Zhu
- Courant Institute of Mathematical Sciences, New York University, New York, New York 10012, United States
| | - Tamar Schlick
- Courant Institute of Mathematical Sciences, New York University, New York, New York 10012, United States
- Department of Chemistry, New York University, New York, New York 10003, United States
- NYU-ECNU Center for Computational Chemistry, NYU Shanghai, Shanghai 200062, P. R. China
| |
Collapse
|
18
|
Schmidt M, Hamacher K, Reinhardt F, Lotz TS, Groher F, Suess B, Jager S. SICOR: Subgraph Isomorphism Comparison of RNA Secondary Structures. IEEE/ACM TRANSACTIONS ON COMPUTATIONAL BIOLOGY AND BIOINFORMATICS 2020; 17:2189-2195. [PMID: 31295116 DOI: 10.1109/tcbb.2019.2926711] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/09/2023]
Abstract
RNA aptamer selection during SELEX experiments builds on secondary structural diversity. Advanced structural comparison methods can focus this diversity. We develop SICOR, which uses probabilistic subgraph isomorphisms for graph distances between RNA secondary structure graphs. SICOR outperforms other comparison methods and is applicable to many structural comparisons in experimental design.
Collapse
|
19
|
Navien TN, Thevendran R, Hamdani HY, Tang TH, Citartan M. In silico molecular docking in DNA aptamer development. Biochimie 2020; 180:54-67. [PMID: 33086095 DOI: 10.1016/j.biochi.2020.10.005] [Citation(s) in RCA: 39] [Impact Index Per Article: 9.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/24/2020] [Revised: 09/23/2020] [Accepted: 10/14/2020] [Indexed: 12/21/2022]
Abstract
Aptamers are single-stranded DNA or RNA oligonucleotides generated by SELEX that exhibit binding affinity and specificity against a wide variety of target molecules. Compared to RNA aptamers, DNA aptamers are much more stable and therefore are widely adopted in a number of applications especially in diagnostics. The tediousness and rigor associated with certain steps of the SELEX intensify the efforts to adopt in silico molecular docking approaches together with in vitro SELEX procedures in developing DNA aptamers. Inspired by these endeavors, we carry out an overview of the in silico molecular docking approaches in DNA aptamer generation, by detailing the stepwise procedures as well as shedding some light on the various softwares used. The in silico maturation strategy and the limitations of the in silico approaches are also underscored.
Collapse
Affiliation(s)
- Tholasi Nadhan Navien
- Advanced Medical & Dental Institute (AMDI), Universiti Sains Malaysia, Bertam, 13200, Kepala Batas, Penang, Malaysia
| | - Ramesh Thevendran
- Advanced Medical & Dental Institute (AMDI), Universiti Sains Malaysia, Bertam, 13200, Kepala Batas, Penang, Malaysia
| | - Hazrina Yusof Hamdani
- Advanced Medical & Dental Institute (AMDI), Universiti Sains Malaysia, Bertam, 13200, Kepala Batas, Penang, Malaysia
| | - Thean-Hock Tang
- Advanced Medical & Dental Institute (AMDI), Universiti Sains Malaysia, Bertam, 13200, Kepala Batas, Penang, Malaysia.
| | - Marimuthu Citartan
- Advanced Medical & Dental Institute (AMDI), Universiti Sains Malaysia, Bertam, 13200, Kepala Batas, Penang, Malaysia.
| |
Collapse
|
20
|
3dRNA: Building RNA 3D structure with improved template library. Comput Struct Biotechnol J 2020; 18:2416-2423. [PMID: 33005304 PMCID: PMC7508704 DOI: 10.1016/j.csbj.2020.08.017] [Citation(s) in RCA: 23] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/06/2020] [Revised: 08/18/2020] [Accepted: 08/21/2020] [Indexed: 11/22/2022] Open
Abstract
Most of computational methods of building RNA tertiary structure are template-based. The template-based methods usually can give more accurate 3D structures due to the use of native 3D templates, but they cannot work if the 3D templates are not available. So, a more complete library of the native 3D templates is very important for this type of methods. 3dRNA is a template-based method for building RNA tertiary structure previously proposed by us. In this paper we report improved 3D template libraries of 3dRNA by using two different schemes that give two libraries 3dRNA_Lib1 and 3dRNA_Lib2. These libraries expand the original one by nearly ten times. Benchmark shows that they can significantly increase the accuracy of 3dRNA, especially in building complex and large RNA 3D structures.
Collapse
|
21
|
Emami N, Pakchin PS, Ferdousi R. Computational predictive approaches for interaction and structure of aptamers. J Theor Biol 2020; 497:110268. [PMID: 32311376 DOI: 10.1016/j.jtbi.2020.110268] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/10/2020] [Revised: 03/27/2020] [Accepted: 04/02/2020] [Indexed: 02/07/2023]
Abstract
Aptamers are short single-strand sequences that can bind to their specific targets with high affinity and specificity. Usually, aptamers are selected experimentally via systematic evolution of ligands by exponential enrichment (SELEX), an evolutionary process that consists of multiple cycles of selection and amplification. The SELEX process is expensive, time-consuming, and its success rates are relatively low. To overcome these difficulties, in recent years, several computational techniques have been developed in aptamer sciences that bring together different disciplines and branches of technologies. In this paper, a complementary review on computational predictive approaches of the aptamer has been organized. Generally, the computational prediction approaches of aptamer have been proposed to carry out in two main categories: interaction-based prediction and structure-based predictions. Furthermore, the available software packages and toolkits in this scope were reviewed. The aim of describing computational methods and tools in aptamer science is that aptamer scientists might take advantage of these computational techniques to develop more accurate and more sensitive aptamers.
Collapse
Affiliation(s)
- Neda Emami
- Department of Health Information Technology, School of Management and Medical Informatics, Tabriz University of Medical Sciences, Tabriz, Iran
| | - Parvin Samadi Pakchin
- Research Center for Pharmaceutical Nanotechnology, Biomedicine Institute, Tabriz University of Medical Sciences, Tabriz, Iran
| | - Reza Ferdousi
- Department of Health Information Technology, School of Management and Medical Informatics, Tabriz University of Medical Sciences, Tabriz, Iran; Research Center for Pharmaceutical Nanotechnology, Biomedicine Institute, Tabriz University of Medical Sciences, Tabriz, Iran.
| |
Collapse
|
22
|
In silico design and validation of high-affinity RNA aptamers targeting epithelial cellular adhesion molecule dimers. Proc Natl Acad Sci U S A 2020; 117:8486-8493. [PMID: 32234785 DOI: 10.1073/pnas.1913242117] [Citation(s) in RCA: 39] [Impact Index Per Article: 9.8] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022] Open
Abstract
Nucleic acid aptamers hold great promise for therapeutic applications due to their favorable intrinsic properties, as well as high-throughput experimental selection techniques. Despite the utility of the systematic evolution of ligands by the exponential enrichment (SELEX) method for aptamer determination, complementary in silico aptamer design is highly sought after to facilitate virtual screening and increased understanding of important nucleic acid-protein interactions. Here, with a combined experimental and theoretical approach, we have developed two optimal epithelial cellular adhesion molecule (EpCAM) aptamers. Our structure-based in silico method first predicts their binding modes and then optimizes them for EpCAM with molecular dynamics simulations, docking, and free energy calculations. Our isothermal titration calorimetry experiments further confirm that the EpCAM aptamers indeed exhibit enhanced affinity over a previously patented nanomolar aptamer, EP23. Moreover, our study suggests that EP23 and the de novo designed aptamers primarily bind to EpCAM dimers (and not monomers, as hypothesized in previous published works), suggesting a paradigm for developing EpCAM-targeted therapies.
Collapse
|
23
|
Perry SL, Sing CE. 100th Anniversary of Macromolecular Science Viewpoint: Opportunities in the Physics of Sequence-Defined Polymers. ACS Macro Lett 2020; 9:216-225. [PMID: 35638672 DOI: 10.1021/acsmacrolett.0c00002] [Citation(s) in RCA: 76] [Impact Index Per Article: 19.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
Polymer science has been driven by ever-increasing molecular complexity, as polymer synthesis expands an already-vast palette of chemical and architectural parameter space. Copolymers represent a key example, where simple homopolymers have given rise to random, alternating, gradient, and block copolymers. Polymer physics has provided the insight needed to explore this monomer sequence parameter space. The future of polymer science, however, must contend with further increases in monomer precision, as this class of macromolecules moves ever closer to the sequence-monodisperse polymers that are the workhorses of biology. The advent of sequence-defined polymers gives rise to opportunities for material design, with increasing levels of chemical information being incorporated into long-chain molecules; however, this also raises questions that polymer physics must address. What properties uniquely emerge from sequence-definition? Is this circumstance-dependent? How do we define and think about sequence dispersity? How do we think about a hierarchy of sequence effects? Are more sophisticated characterization methods, as well as theoretical and computational tools, needed to understand this class of macromolecules? The answers to these questions touch on many difficult scientific challenges, setting the stage for a rich future for sequence-defined polymers in polymer physics.
Collapse
Affiliation(s)
- Sarah L. Perry
- Department of Chemical Engineering, University of Massachusetts−Amherst, 686 North Pleasant Street, Amherst, Massachusetts 01003, United States
| | - Charles E. Sing
- Department of Chemical and Biomolecular Engineering, University of Illinois at Urbana−Champaign, 600 South Mathews Avenue Urbana, Illinois 61801, United States
| |
Collapse
|
24
|
Miyamoto S, Noda T. In vitro vRNA-vRNA interactions in the H1N1 influenza A virus genome. Microbiol Immunol 2020; 64:202-209. [PMID: 31840833 DOI: 10.1111/1348-0421.12766] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/30/2019] [Revised: 12/09/2019] [Accepted: 12/12/2019] [Indexed: 01/15/2023]
Abstract
The genome of influenza A virus consists of eight-segmented, single-stranded, negative-sense viral RNAs (vRNAs). Each vRNA contains a central coding region that is flanked by noncoding regions. It has been shown that upon virion formation, the eight vRNAs are selectively packaged into progeny virions through segment-specific packaging signals that are located in both the terminal coding regions and adjacent noncoding regions of each vRNA. Although recent studies using next-generation sequencing suggest that multiple intersegment interactions are involved in genome packaging, contributions of the packaging signals to the intersegment interactions are not fully understood. Herein, using synthesized full-length vRNAs of H1N1 WSN (A/WSN/33 [H1N1]) virus and short vRNAs containing the packaging signal sequences, we performed in vitro RNA binding assays and identified 15 intersegment interactions among eight vRNAs, most of which were mediated by the 3'- and 5'-terminal regions. Interestingly, all eight vRNAs interacted with multiple other vRNAs, in that some bound to different vRNAs through their respective 3'- and 5'-terminal regions. These in vitro findings would be of use in future studies of in vivo vRNA-vRNA interactions during selective genome packaging.
Collapse
Affiliation(s)
- Sho Miyamoto
- Laboratory of Ultrastructural Virology, Institute for Frontier Life and Medical Sciences, Kyoto University, Kyoto, Japan.,Department of Molecular Virology, Graduate School of Medicine, Kyoto University, Kyoto, Japan
| | - Takeshi Noda
- Laboratory of Ultrastructural Virology, Institute for Frontier Life and Medical Sciences, Kyoto University, Kyoto, Japan.,Laboratory of Ultrastructural Virology, Graduate School of Biostudies, Kyoto University, Kyoto, Japan
| |
Collapse
|
25
|
Inverse folding with RNA-As-Graphs produces a large pool of candidate sequences with target topologies. J Struct Biol 2019; 209:107438. [PMID: 31874236 DOI: 10.1016/j.jsb.2019.107438] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/04/2019] [Revised: 12/18/2019] [Accepted: 12/19/2019] [Indexed: 02/07/2023]
Abstract
We present an RNA-As-Graphs (RAG) based inverse folding algorithm, RAG-IF, to design novel RNA sequences that fold onto target tree graph topologies. The algorithm can be used to enhance our recently reported computational design pipeline (Jain et al., NAR 2018). The RAG approach represents RNA secondary structures as tree and dual graphs, where RNA loops and helices are coarse-grained as vertices and edges, opening the usage of graph theory methods to study, predict, and design RNA structures. Our recently developed computational pipeline for design utilizes graph partitioning (RAG-3D) and atomic fragment assembly (F-RAG) to design sequences to fold onto RNA-like tree graph topologies; the atomic fragments are taken from existing RNA structures that correspond to tree subgraphs. Because F-RAG may not produce the target folds for all designs, automated mutations by RAG-IF algorithm enhance the candidate pool markedly. The crucial residues for mutation are identified by differences between the predicted and the target topology. A genetic algorithm then mutates the selected residues, and the successful sequences are optimized to retain only the minimal or essential mutations. Here we evaluate RAG-IF for 6 RNA-like topologies and generate a large pool of successful candidate sequences with a variety of minimal mutations. We find that RAG-IF adds robustness and efficiency to our RNA design pipeline, making inverse folding motivated by graph topology rather than secondary structure more productive.
Collapse
|
26
|
Kimchi O, Cragnolini T, Brenner MP, Colwell LJ. A Polymer Physics Framework for the Entropy of Arbitrary Pseudoknots. Biophys J 2019; 117:520-532. [PMID: 31353036 PMCID: PMC6697467 DOI: 10.1016/j.bpj.2019.06.037] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/08/2018] [Revised: 06/21/2019] [Accepted: 06/27/2019] [Indexed: 11/18/2022] Open
Abstract
The accurate prediction of RNA secondary structure from primary sequence has had enormous impact on research from the past 40 years. Although many algorithms are available to make these predictions, the inclusion of non-nested loops, termed pseudoknots, still poses challenges arising from two main factors: 1) no physical model exists to estimate the loop entropies of complex intramolecular pseudoknots, and 2) their NP-complete enumeration has impeded their study. Here, we address both challenges. First, we develop a polymer physics model that can address arbitrarily complex pseudoknots using only two parameters corresponding to concrete physical quantities-over an order of magnitude fewer than the sparsest state-of-the-art phenomenological methods. Second, by coupling this model to exhaustive enumeration of the set of possible structures, we compute the entire free energy landscape of secondary structures resulting from a primary RNA sequence. We demonstrate that for RNA structures of ∼80 nucleotides, with minimal heuristics, the complete enumeration of possible secondary structures can be accomplished quickly despite the NP-complete nature of the problem. We further show that despite our loop entropy model's parametric sparsity, it performs better than or on par with previously published methods in predicting both pseudoknotted and non-pseudoknotted structures on a benchmark data set of RNA structures of ≤80 nucleotides. We suggest ways in which the accuracy of the model can be further improved.
Collapse
Affiliation(s)
- Ofer Kimchi
- Harvard Graduate Program in Biophysics, Harvard University, Cambridge, Massachusetts.
| | - Tristan Cragnolini
- Department of Chemistry, University of Cambridge, Cambridge, United Kingdom
| | - Michael P Brenner
- School of Engineering and Applied Sciences, Cambridge, Massachusetts; Kavli Institute for Bionano Science and Technology, Harvard University, Cambridge, Massachusetts
| | - Lucy J Colwell
- Department of Chemistry, University of Cambridge, Cambridge, United Kingdom.
| |
Collapse
|
27
|
Abstract
The three-dimensional structures of RNA molecules provide rich and often critical information for understanding their functions, including how they recognize small molecule and protein partners. Computational modeling of RNA 3D structure is becoming increasingly accurate, particularly with the availability of growing numbers of template structures already solved experimentally and the development of sequence alignment and 3D modeling tools to take advantage of this database. For several recent "RNA puzzle" blind modeling challenges, we have successfully identified useful template structures and achieved accurate structure predictions through homology modeling tools developed in the Rosetta software suite. We describe our semi-automated methodology here and walk through two illustrative examples: an adenine riboswitch aptamer, modeled from a template guanine riboswitch structure, and a SAM I/IV riboswitch aptamer, modeled from a template SAM I riboswitch structure.
Collapse
Affiliation(s)
- Andrew M Watkins
- Department of Biochemistry, Stanford University School of Medicine, Stanford, CA, United States
| | - Ramya Rangan
- Biophysics Program, Stanford University, Stanford, CA, United States
| | - Rhiju Das
- Department of Biochemistry, Stanford University School of Medicine, Stanford, CA, United States; Biophysics Program, Stanford University, Stanford, CA, United States.
| |
Collapse
|
28
|
Mak CH, Phan ENH. Topological Constraints and Their Conformational Entropic Penalties on RNA Folds. Biophys J 2019; 114:2059-2071. [PMID: 29742400 PMCID: PMC5961522 DOI: 10.1016/j.bpj.2018.03.035] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/03/2017] [Revised: 03/08/2018] [Accepted: 03/22/2018] [Indexed: 02/01/2023] Open
Abstract
Functional RNAs can fold into intricate structures using a number of different secondary and tertiary structural motifs. Many factors contribute to the overall free energy of the target fold. This study aims at quantifying the entropic costs coming from the loss of conformational freedom when the sugar-phosphate backbone is subjected to constraints imposed by secondary and tertiary contacts. Motivated by insights from topology theory, we design a diagrammatic scheme to represent different types of RNA structures so that constraints associated with a folded structure may be segregated into mutually independent subsets, enabling the total conformational entropy loss to be easily calculated as a sum of independent terms. We used high-throughput Monte Carlo simulations to simulate large ensembles of single-stranded RNA sequences in solution to validate the assumptions behind our diagrammatic scheme, examining the entropic costs for hairpin initiation and formation of many multiway junctions. Our diagrammatic scheme aids in the factorization of secondary/tertiary constraints into distinct topological classes and facilitates the discovery of interrelationships among multiple constraints on RNA folds. This perspective, which to our knowledge is novel, leads to useful insights into the inner workings of some functional RNA sequences, demonstrating how they might operate by transforming their structures among different topological classes.
Collapse
Affiliation(s)
- Chi H Mak
- Department of Chemistry, University of Southern California, Los Angeles, California; Center of Applied Mathematical Sciences, University of Southern California, Los Angeles, California; Department of Biological Sciences, University of Southern California, Los Angeles, California.
| | - Ethan N H Phan
- Department of Chemistry, University of Southern California, Los Angeles, California
| |
Collapse
|
29
|
Pal A, Levy Y. Structure, stability and specificity of the binding of ssDNA and ssRNA with proteins. PLoS Comput Biol 2019; 15:e1006768. [PMID: 30933978 PMCID: PMC6467422 DOI: 10.1371/journal.pcbi.1006768] [Citation(s) in RCA: 36] [Impact Index Per Article: 7.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/13/2018] [Revised: 04/16/2019] [Accepted: 01/01/2019] [Indexed: 02/06/2023] Open
Abstract
Recognition of single-stranded DNA (ssDNA) or single-stranded RNA (ssRNA) is important for many fundamental cellular functions. A variety of single-stranded DNA-binding proteins (ssDBPs) and single-stranded RNA-binding proteins (ssRBPs) have evolved that bind ssDNA and ssRNA, respectively, with varying degree of affinities and specificities to form complexes. Structural studies of these complexes provide key insights into their recognition mechanism. However, computational modeling of the specific recognition process and to predict the structure of the complex is challenging, primarily due to the heterogeneity of their binding energy landscape and the greater flexibility of ssDNA or ssRNA compared with double-stranded nucleic acids. Consequently, considerably fewer computational studies have explored interactions between proteins and single-stranded nucleic acids compared with protein interactions with double-stranded nucleic acids. Here, we report a newly developed energy-based coarse-grained model to predict the structure of ssDNA–ssDBP and ssRNA–ssRBP complexes and to assess their sequence-specific interactions and stabilities. We tuned two factors that can modulate specific recognition: base–aromatic stacking strength and the flexibility of the single-stranded nucleic acid. The model was successfully applied to predict the binding conformations of 12 distinct ssDBP and ssRBP structures with their cognate ssDNA and ssRNA partners having various sequences. Estimated binding energies agreed well with the corresponding experimental binding affinities. Bound conformations from the simulation showed a funnel-shaped binding energy distribution where the native-like conformations corresponded to the energy minima. The various ssDNA–protein and ssRNA–protein complexes differed in the balance of electrostatic and aromatic energies. The lower affinity of the ssRNA–ssRBP complexes compared with the ssDNA–ssDBP complexes stems from lower flexibility of ssRNA compared to ssDNA, which results in higher rate constants for the dissociation of the complex (koff) for complexes involving the former. Quantifying bimolecular self-assembly is pivotal to understanding cellular function. In recent years, a large progress has been made in understanding the structure and biophysics of protein-protein interactions. Particularly, various computational tools are available for predicting these structures and to estimate their stability and the driving forces of their formation. The understating of the interactions between proteins and nucleic acids, however, is still limited, presumably due to the involvement of non-specific interactions as well as the high conformational plasticity that may demand an induced-fit mechanism. In particular, the interactions between proteins and single-stranded nucleic acids (i.e., single-stranded DNA and RNA) is very challenging due to their high flexibility. Furthermore, the interface between proteins and single-stranded nucleic acids is often chemically more heterogeneous than the interface between proteins and double-stranded DNA. In this study, we developed a coarse-grained computational model to predict the structure of complexes between proteins and single-stranded nucleic acids. The model was applied to estimate binding affinities and the estimated binding energies agreed well with the corresponding experimental binding affinities. The kinetics of association as well as the specificity of the complexes between proteins and ssDNA are different than those with ssRNA, mostly due to differences in their conformational flexibility.
Collapse
Affiliation(s)
- Arumay Pal
- Department of Structural Biology, Weizmann Institute of Science, Rehovot, Israel
| | - Yaakov Levy
- Department of Structural Biology, Weizmann Institute of Science, Rehovot, Israel
- * E-mail:
| |
Collapse
|
30
|
Jain S, Saju S, Petingi L, Schlick T. An extended dual graph library and partitioning algorithm applicable to pseudoknotted RNA structures. Methods 2019; 162-163:74-84. [PMID: 30928508 DOI: 10.1016/j.ymeth.2019.03.022] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/14/2019] [Revised: 02/28/2019] [Accepted: 03/22/2019] [Indexed: 12/18/2022] Open
Abstract
Exploring novel RNA topologies is imperative for understanding RNA structure and pursuing its design. Our RNA-As-Graphs (RAG) approach exploits graph theory tools and uses coarse-grained tree and dual graphs to represent RNA helices and loops by vertices and edges. Only dual graphs represent pseudoknotted RNAs fully. Here we develop a dual graph enumeration algorithm to generate an expanded library of dual graph topologies for 2-9 vertices, and extend our dual graph partitioning algorithm to identify all possible RNA subgraphs. Our enumeration algorithm connects smaller-vertex graphs, using all possible edge combinations, to build larger-vertex graphs and retain all non-isomorphic graph topologies, thereby more than doubling the size of our prior library to a total of 110,667 dual graph topologies. We apply our dual graph partitioning algorithm, which keeps pseudoknots and junctions intact, to all existing RNA structures to identify all possible substructures up to 9 vertices. In addition, our expanded dual graph library assigns graph topologies to all RNA graphs and subgraphs, rectifying prior inconsistencies. We update our RAG-3Dual database of RNA atomic fragments with all newly identified substructures and their graph IDs, increasing its size by more than 50 times. The enlarged dual graph library and RAG-3Dual database provide a comprehensive repertoire of graph topologies and atomic fragments to study yet undiscovered RNA molecules and design RNA sequences with novel topologies, including a variety of pseudoknotted RNAs.
Collapse
Affiliation(s)
- Swati Jain
- Department of Chemistry, New York University, 1021 Silver, 100 Washington Square East, New York, NY 10003, USA
| | - Sera Saju
- Department of Chemistry, New York University, 1021 Silver, 100 Washington Square East, New York, NY 10003, USA
| | - Louis Petingi
- Computer Science Department, College of Staten Island, City University of New York, Staten Island, New York, NY 10314, USA
| | - Tamar Schlick
- Department of Chemistry, New York University, 1021 Silver, 100 Washington Square East, New York, NY 10003, USA; Courant Institute of Mathematical Sciences, New York University, 251 Mercer Street, New York, NY 10012, USA; NYU-East China Normal University Center for Computational Chemistry at New York University Shanghai, Room 340, Geography Building, North Zhongshan Road, 3663 Shanghai, China.
| |
Collapse
|
31
|
Jain S, Laederach A, Ramos SBV, Schlick T. A pipeline for computational design of novel RNA-like topologies. Nucleic Acids Res 2018; 46:7040-7051. [PMID: 30137633 PMCID: PMC6101589 DOI: 10.1093/nar/gky524] [Citation(s) in RCA: 20] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/24/2017] [Revised: 05/22/2018] [Accepted: 05/24/2018] [Indexed: 12/11/2022] Open
Abstract
Designing novel RNA topologies is a challenge, with important therapeutic and industrial applications. We describe a computational pipeline for design of novel RNA topologies based on our coarse-grained RNA-As-Graphs (RAG) framework. RAG represents RNA structures as tree graphs and describes RNA secondary (2D) structure topologies (currently up to 13 vertices, ≈260 nucleotides). We have previously identified novel graph topologies that are RNA-like among these. Here we describe a systematic design pipeline and illustrate design for six broad design problems using recently developed tools for graph-partitioning and fragment assembly (F-RAG). Following partitioning of the target graph, corresponding atomic fragments from our RAG-3D database are combined using F-RAG, and the candidate atomic models are scored using a knowledge-based potential developed for 3D structure prediction. The sequences of the top scoring models are screened further using available tools for 2D structure prediction. The results indicate that our modular approach based on RNA-like topologies rather than specific 2D structures allows for greater flexibility in the design process, and generates a large number of candidate sequences quickly. Experimental structure probing using SHAPE-MaP for two sequences agree with our predictions and suggest that our combined tools yield excellent candidates for further sequence and experimental screening.
Collapse
Affiliation(s)
- Swati Jain
- Department of Chemistry, New York University, 1001 Silver, 100 Washington Square East, New York, NY 10003, USA
| | - Alain Laederach
- Department of Biology, University of North Carolina at Chapel Hill, Chapel Hill, NC 27599, USA
| | - Silvia B V Ramos
- Department of Biochemistry and Biophysics, University of North Carolina at Chapel Hill, Chapel Hill, NC 27599, USA
| | - Tamar Schlick
- Department of Chemistry, New York University, 1001 Silver, 100 Washington Square East, New York, NY 10003, USA
- Courant Institute of Mathematical Sciences, New York University, 251 Mercer Street, New York, NY 10012, USA
- NYU-ECNU Center for Computational Chemistry at New York University Shanghai, Room 340, Geography Building, North Zhongshan Road, 3663 Shanghai, China
| |
Collapse
|
32
|
Jain S, Bayrak CS, Petingi L, Schlick T. Dual Graph Partitioning Highlights a Small Group of Pseudoknot-Containing RNA Submotifs. Genes (Basel) 2018; 9:E371. [PMID: 30044451 PMCID: PMC6115904 DOI: 10.3390/genes9080371] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/07/2018] [Revised: 06/26/2018] [Accepted: 06/26/2018] [Indexed: 12/31/2022] Open
Abstract
RNA molecules are composed of modular architectural units that define their unique structural and functional properties. Characterization of these building blocks can help interpret RNA structure/function relationships. We present an RNA secondary structure motif and submotif library using dual graph representation and partitioning. Dual graphs represent RNA helices as vertices and loops as edges. Unlike tree graphs, dual graphs can represent RNA pseudoknots (intertwined base pairs). For a representative set of RNA structures, we construct dual graphs from their secondary structures, and apply our partitioning algorithm to identify non-separable subgraphs (or blocks) without breaking pseudoknots. We report 56 subgraph blocks up to nine vertices; among them, 22 are frequently occurring, 15 of which contain pseudoknots. We then catalog atomic fragments corresponding to the subgraph blocks to define a library of building blocks that can be used for RNA design, which we call RAG-3Dual, as we have done for tree graphs. As an application, we analyze the distribution of these subgraph blocks within ribosomal RNAs of various prokaryotic and eukaryotic species to identify common subgraphs and possible ancestry relationships. Other applications of dual graph partitioning and motif library can be envisioned for RNA structure analysis and design.
Collapse
Affiliation(s)
- Swati Jain
- Department of Chemistry, New York University, New York, NY 10003, USA.
| | - Cigdem S Bayrak
- Department of Chemistry, New York University, New York, NY 10003, USA.
| | - Louis Petingi
- Computer Science Department, College of Staten Island, City University of New York, Staten Island, New York, NY 10314, USA.
| | - Tamar Schlick
- Department of Chemistry, New York University, New York, NY 10003, USA.
- Courant Institute of Mathematical Sciences, New York University, New York, NY 10012, USA.
- NYU-East China Normal University Center for Computational Chemistry, New York University Shanghai, Shanghai 3663, China.
| |
Collapse
|
33
|
Abstract
The structure of RNA has been a natural subject for mathematical modeling, inviting many innovative computational frameworks. This single-stranded polynucleotide chain can fold upon itself in numerous ways to form hydrogen-bonded segments, imperfect with single-stranded loops. Illustrating these paired and non-paired interaction networks, known as RNA's secondary (2D) structure, using mathematical graph objects has been illuminating for RNA structure analysis. Building upon such seminal work from the 1970s and 1980s, graph models are now used to study not only RNA structure but also describe RNA's recurring modular units, sample the conformational space accessible to RNAs, predict RNA's three-dimensional folds, and apply the combined aspects to novel RNA design. In this article, we outline the development of the RNA-As-Graphs (or RAG) approach and highlight current applications to RNA structure prediction and design.
Collapse
Affiliation(s)
- Tamar Schlick
- Department of Chemistry, 100 Washington Square East, Silver Building, New York University, New York, NY 10003, USA; Courant Institute of Mathematical Sciences, New York University, 251 Mercer St., New York, NY 10012, USA; New York University ECNU - Center for Computational Chemistry at NYU Shanghai, 3663 North Zhongshan Road, Shanghai, 200062, China.
| |
Collapse
|
34
|
Lošdorfer Božič A, Micheletti C, Podgornik R, Tubiana L. Compactness of viral genomes: effect of disperse and localized random mutations. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2018; 30:084006. [PMID: 29334364 PMCID: PMC7104904 DOI: 10.1088/1361-648x/aaa7b0] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 09/28/2017] [Revised: 12/29/2017] [Accepted: 01/15/2018] [Indexed: 06/07/2023]
Abstract
Genomes of single-stranded RNA viruses have evolved to optimize several concurrent properties. One of them is the architecture of their genomic folds, which must not only feature precise structural elements at specific positions, but also allow for overall spatial compactness. The latter was shown to be disrupted by random synonymous mutations, a disruption which can consequently negatively affect genome encapsidation. In this study, we use three mutation schemes with different degrees of locality to mutate the genomes of phage MS2 and Brome Mosaic virus in order to understand the observed sensitivity of the global compactness of their folds. We find that mutating local stretches of their genomes' sequence or structure is less disruptive to their compactness compared to inducing randomly-distributed mutations. Our findings are indicative of a mechanism for the conservation of compactness acting on a global scale of the genomes, and have several implications for understanding the interplay between local and global architecture of viral RNA genomes.
Collapse
Affiliation(s)
- Anže Lošdorfer Božič
- Department of Theoretical Physics, Jožef Stefan Institute, SI-1000 Ljubljana, Slovenia
| | | | - Rudolf Podgornik
- Department of Theoretical Physics, Jožef Stefan Institute, SI-1000 Ljubljana, Slovenia
- Department of Physics, Faculty of Mathematics and Physics, University of Ljubljana, SI-1000 Ljubljana, Slovenia
| | - Luca Tubiana
- Faculty of Physics, University of Vienna, A-1090 Vienna, Austria
| |
Collapse
|
35
|
Williams B, Zhao B, Tandon A, Ding F, Weeks KM, Zhang Q, Dokholyan NV. Structure modeling of RNA using sparse NMR constraints. Nucleic Acids Res 2018; 45:12638-12647. [PMID: 29165648 PMCID: PMC5728392 DOI: 10.1093/nar/gkx1058] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/12/2017] [Accepted: 10/18/2017] [Indexed: 01/04/2023] Open
Abstract
RNAs fold into distinct molecular conformations that are often essential for their functions. Accurate structure modeling of complex RNA motifs, including ubiquitous non-canonical base pairs and pseudoknots, remains a challenge. Here, we present an NMR-guided all-atom discrete molecular dynamics (DMD) platform, iFoldNMR, for rapid and accurate structure modeling of complex RNAs. We show that sparse distance constraints from imino resonances, which can be readily obtained from routine NMR experiments and easier to compile than laborious assignments of non-solvent-exchangeable protons, are sufficient to direct a DMD search for low-energy RNA conformers. Benchmarking on a set of RNAs with complex folds spanning up to 56 nucleotides in length yields structural models that recapitulate experimentally determined structures with all-heavy-atom RMSDs ranging from 2.4 to 6.5 Å. This platform represents an efficient approach for high-throughput RNA structure modeling and will facilitate analysis of diverse, newly discovered functional RNAs.
Collapse
Affiliation(s)
- Benfeard Williams
- Department of Biochemistry and Biophysics, University of North Carolina at Chapel Hill, Chapel Hill, NC 27599, USA.,Molecular and Cellular Biophysics Program, University of North Carolina at Chapel Hill, Chapel Hill, NC 27599, USA
| | - Bo Zhao
- Molecular and Cellular Biophysics Program, University of North Carolina at Chapel Hill, Chapel Hill, NC 27599, USA.,Department of Chemistry, University of North Carolina at Chapel Hill, Chapel Hill, NC 27599, USA
| | - Arpit Tandon
- Department of Biochemistry and Biophysics, University of North Carolina at Chapel Hill, Chapel Hill, NC 27599, USA.,Molecular and Cellular Biophysics Program, University of North Carolina at Chapel Hill, Chapel Hill, NC 27599, USA
| | - Feng Ding
- Department of Physics and Astronomy, Clemson University, Clemson, SC 29634, USA
| | - Kevin M Weeks
- Department of Chemistry, University of North Carolina at Chapel Hill, Chapel Hill, NC 27599, USA
| | - Qi Zhang
- Department of Biochemistry and Biophysics, University of North Carolina at Chapel Hill, Chapel Hill, NC 27599, USA.,Molecular and Cellular Biophysics Program, University of North Carolina at Chapel Hill, Chapel Hill, NC 27599, USA
| | - Nikolay V Dokholyan
- Department of Biochemistry and Biophysics, University of North Carolina at Chapel Hill, Chapel Hill, NC 27599, USA.,Molecular and Cellular Biophysics Program, University of North Carolina at Chapel Hill, Chapel Hill, NC 27599, USA
| |
Collapse
|
36
|
Xu X, Chen SJ. Hierarchical Assembly of RNA Three-Dimensional Structures Based on Loop Templates. J Phys Chem B 2018; 122:5327-5335. [PMID: 29258305 DOI: 10.1021/acs.jpcb.7b10102] [Citation(s) in RCA: 33] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The current RNA structure prediction methods cannot keep up the pace of the rapidly increasing number of sequences and the emerging new functions of RNAs. Template-based RNA three-dimensional structure prediction methods are restricted by the limited number of known RNA structures, and traditional motif-based search for the templates does not always lead to successful results. Here we report a new template search and assembly algorithm, the hierarchical loop template-assembly method (VfoldLA). The method searches for templates for single strand loop/junctions instead of the whole motifs, which often renders no available templates, or short fragments (several nucleotides), which requires a long computational time to assemble and refine. The VfoldLA method has the advantage of accounting for local and nonlocal interloop interactions. Benchmark tests indicate that this new method can provide low-resolution predictions for RNA conformations at different levels of structural complexities. Furthermore, the VfoldLA-predicted conformations may also serve as reliable putative models for further structure prediction and refinements. VfoldLA is accessible at http://rna.physics.missouri.edu/vfoldLA .
Collapse
Affiliation(s)
- Xiaojun Xu
- Institute of Bioinformatics and Medical Engineering , Jiangsu University of Technology , Changzhou , Jiangsu 213001 , China.,Department of Physics, Department of Biochemistry, and Informatics Institute , University of Missouri , Columbia , Missouri 65211 , United States
| | - Shi-Jie Chen
- Department of Physics, Department of Biochemistry, and Informatics Institute , University of Missouri , Columbia , Missouri 65211 , United States
| |
Collapse
|
37
|
Takahashi M. Aptamers targeting cell surface proteins. Biochimie 2017; 145:63-72. [PMID: 29198584 DOI: 10.1016/j.biochi.2017.11.019] [Citation(s) in RCA: 35] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/01/2017] [Accepted: 11/28/2017] [Indexed: 02/07/2023]
Abstract
High affinity binders targeting specific cell surface proteins are vital for development of basic and applied biosciences. However, despite sustained efforts to generate such binders by chemicals and antibodies, there are still many cell surface proteins that lack high affinity binders. Nucleic acid aptamers have potential as binding molecules for cell surface proteins, because they form distinct structures that have high affinity and specificity for a wide range of targets. Aptamers are isolated from large combinatorial libraries using a unique iterative selection-amplification process known as systematic evolution of ligands by exponential enrichment (SELEX). Among advantages of this method, purified and complex heterogeneous targets, such as bacteria, viruses, and whole-living cells, can be used for selection of aptamers. Moreover, SELEX allows generation of cell-surface-specific aptamers without prior knowledge of expression profiles in target cells. Therefore, the technology has been widely used as a valid and feasible method to generate aptamers binding to cell surface proteins with intact structure. Herein, this review summarizes and updates iconic SELEX technologies that target membrane proteins.
Collapse
Affiliation(s)
- Masaki Takahashi
- Division of RNA Medical Science, The Institute of Medical Science, The University of Tokyo, 4-6-1, Shirokanedai, Minato-ku, Tokyo, Japan.
| |
Collapse
|
38
|
Jain S, Schlick T. F-RAG: Generating Atomic Coordinates from RNA Graphs by Fragment Assembly. J Mol Biol 2017; 429:3587-3605. [PMID: 28988954 PMCID: PMC5693719 DOI: 10.1016/j.jmb.2017.09.017] [Citation(s) in RCA: 25] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/09/2017] [Revised: 09/12/2017] [Accepted: 09/22/2017] [Indexed: 10/18/2022]
Abstract
Coarse-grained models represent attractive approaches to analyze and simulate ribonucleic acid (RNA) molecules, for example, for structure prediction and design, as they simplify the RNA structure to reduce the conformational search space. Our structure prediction protocol RAGTOP (RNA-As-Graphs Topology Prediction) represents RNA structures as tree graphs and samples graph topologies to produce candidate graphs. However, for a more detailed study and analysis, construction of atomic from coarse-grained models is required. Here we present our graph-based fragment assembly algorithm (F-RAG) to convert candidate three-dimensional (3D) tree graph models, produced by RAGTOP into atomic structures. We use our related RAG-3D utilities to partition graphs into subgraphs and search for structurally similar atomic fragments in a data set of RNA 3D structures. The fragments are edited and superimposed using common residues, full atomic models are scored using RAGTOP's knowledge-based potential, and geometries of top scoring models is optimized. To evaluate our models, we assess all-atom RMSDs and Interaction Network Fidelity (a measure of residue interactions) with respect to experimentally solved structures and compare our results to other fragment assembly programs. For a set of 50 RNA structures, we obtain atomic models with reasonable geometries and interactions, particularly good for RNAs containing junctions. Additional improvements to our protocol and databases are outlined. These results provide a good foundation for further work on RNA structure prediction and design applications.
Collapse
Affiliation(s)
- Swati Jain
- Department of Chemistry, New York University, 1001 Silver, 100 Washington Square East, New York, NY 10003, USA
| | - Tamar Schlick
- Department of Chemistry, New York University, 1001 Silver, 100 Washington Square East, New York, NY 10003, USA; Courant Institute of Mathematical Sciences, New York University, 251 Mercer Street, New York, NY 10012, USA; New York University-East China Normal University Center for Computational Chemistry at New York University Shanghai, Room 340, Geography Building, North Zhongshan Road, 3663 Shanghai, China.
| |
Collapse
|
39
|
Wang J, Mao K, Zhao Y, Zeng C, Xiang J, Zhang Y, Xiao Y. Optimization of RNA 3D structure prediction using evolutionary restraints of nucleotide-nucleotide interactions from direct coupling analysis. Nucleic Acids Res 2017; 45:6299-6309. [PMID: 28482022 PMCID: PMC5499770 DOI: 10.1093/nar/gkx386] [Citation(s) in RCA: 71] [Impact Index Per Article: 10.1] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/27/2016] [Accepted: 04/27/2017] [Indexed: 01/01/2023] Open
Abstract
Direct coupling analysis of nucleotide coevolution provides a novel approach to identify which nucleotides in an RNA molecule are likely in direct contact, and this information obtained from sequence only can be used to predict RNA 3D structures with much improved accuracy. Here we present an efficient method that incorporates this information into current RNA 3D structure prediction methods, specifically 3dRNA. Our method makes much more accurate RNA 3D structure prediction than the original 3dRNA as well as other existing prediction methods that used the direct coupling analysis. In particular our method demonstrates a significant improvement in predicting multi-branch junction conformations, a major bottleneck for RNA 3D structure prediction. We also show that our method can be used to optimize the predictions by other methods. These results indicate that optimization of RNA 3D structure prediction using evolutionary restraints of nucleotide-nucleotide interactions from direct coupling analysis offers an efficient way for accurate RNA tertiary structure predictions.
Collapse
Affiliation(s)
- Jian Wang
- Institute of Biophysics, School of Physics and Key Laboratory of Molecular Biophysics of the Ministry of Education, Huazhong University of Science and Technology, Wuhan 430074, Hubei, China
| | - Kangkun Mao
- Institute of Biophysics, School of Physics and Key Laboratory of Molecular Biophysics of the Ministry of Education, Huazhong University of Science and Technology, Wuhan 430074, Hubei, China
| | - Yunjie Zhao
- Institute of Biophysics and Department of Physics, Central China Normal University, Wuhan 430079, China
| | - Chen Zeng
- Department of Physics, The George Washington University, Washington, DC 20052, USA.,School of Life Sciences, Jianghan University, Wuhan 430056, China
| | - Jianjin Xiang
- Institute of Biophysics, School of Physics and Key Laboratory of Molecular Biophysics of the Ministry of Education, Huazhong University of Science and Technology, Wuhan 430074, Hubei, China
| | - Yi Zhang
- Institute of Biophysics, School of Physics and Key Laboratory of Molecular Biophysics of the Ministry of Education, Huazhong University of Science and Technology, Wuhan 430074, Hubei, China
| | - Yi Xiao
- Institute of Biophysics, School of Physics and Key Laboratory of Molecular Biophysics of the Ministry of Education, Huazhong University of Science and Technology, Wuhan 430074, Hubei, China
| |
Collapse
|
40
|
Borodavka A, Dykeman EC, Schrimpf W, Lamb DC. Protein-mediated RNA folding governs sequence-specific interactions between rotavirus genome segments. eLife 2017; 6:27453. [PMID: 28922109 PMCID: PMC5621836 DOI: 10.7554/elife.27453] [Citation(s) in RCA: 51] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/04/2017] [Accepted: 08/03/2017] [Indexed: 01/06/2023] Open
Abstract
Segmented RNA viruses are ubiquitous pathogens, which include influenza viruses and rotaviruses. A major challenge in understanding their assembly is the combinatorial problem of a non-random selection of a full genomic set of distinct RNAs. This process involves complex RNA-RNA and protein-RNA interactions, which are often obscured by non-specific binding at concentrations approaching in vivo assembly conditions. Here, we present direct experimental evidence of sequence-specific inter-segment interactions between rotavirus RNAs, taking place in a complex RNA- and protein-rich milieu. We show that binding of the rotavirus-encoded non-structural protein NSP2 to viral ssRNAs results in the remodeling of RNA, which is conducive to formation of stable inter-segment contacts. To identify the sites of these interactions, we have developed an RNA-RNA SELEX approach for mapping the sequences involved in inter-segment base-pairing. Our findings elucidate the molecular basis underlying inter-segment interactions in rotaviruses, paving the way for delineating similar RNA-RNA interactions that govern assembly of other segmented RNA viruses. Rotavirus is a highly infectious virus that affects children worldwide, causing severe diarrhoea. Despite the introduction of several highly effective vaccines, more than 200,000 children still die from rotavirus each year. There are currently no drugs that can combat this disease once a child has been infected. Viruses carry the instructions that determine their properties and behavior in molecules of DNA or RNA. Unlike many other viruses, which typically have a single molecule of DNA or RNA, rotavirus has 11 distinct “RNA segments”. After invading a cell the virus begins to replicate itself. During replication, the RNA segments (which consist of two strands of RNA paired together) are copied many times. It is not clear how rotaviruses ‘count’ up to 11 so that each new virus acquires a single copy of each segment. Previous biochemical and structural studies of rotavirus replication suggest that selecting 11 distinct RNA segments must involve the RNAs forming complex interactions with proteins and other RNA molecules. Using a highly sensitive fluorescence-based approach, termed fluorescence cross-correlation spectroscopy, Borodavka et al. now present direct experimental evidence of interactions between the RNA segments that occur via single strands of the rotavirus RNA. These RNA-RNA interactions require the binding of a rotavirus protein NSP2 to the RNA strands, which results in the remodeling of the RNA; this remodeling is required to form stable contacts between different RNA segments. Furthermore, a new experimental approach (called RNA-RNA SELEX) developed by Borodavka et al. identified the parts of the RNA segments that may take part in these interactions. The results presented by Borodavka et al. pave the way for identifying the RNA-RNA interactions that govern how other segmented RNA viruses can package their genetic material. Further work to uncover the entire RNA interaction network in rotaviruses would also accelerate the design of new vaccines and may help us to develop antiviral drugs to treat infections.
Collapse
Affiliation(s)
- Alexander Borodavka
- Astbury Centre for Structural Molecular Biology, School of Molecular and Cellular Biology, University of Leeds, Leeds, United Kingdom.,Department of Chemistry, Center for NanoScience, Nanosystems Initiative Munich (NIM) and Center for Integrated Protein Science Munich (CiPSM), Ludwig-Maximilian University of Munich, Munich, Germany
| | - Eric C Dykeman
- York Centre for Complex Systems Analysis, University of York, York, United Kingdom.,Department of Mathematics, University of York, York, United Kingdom.,Department of Biology, University of York, York, United Kingdom
| | - Waldemar Schrimpf
- Department of Chemistry, Center for NanoScience, Nanosystems Initiative Munich (NIM) and Center for Integrated Protein Science Munich (CiPSM), Ludwig-Maximilian University of Munich, Munich, Germany
| | - Don C Lamb
- Department of Chemistry, Center for NanoScience, Nanosystems Initiative Munich (NIM) and Center for Integrated Protein Science Munich (CiPSM), Ludwig-Maximilian University of Munich, Munich, Germany
| |
Collapse
|
41
|
Jager S, Schiller B, Babel P, Blumenroth M, Strufe T, Hamacher K. StreAM-[Formula: see text]: algorithms for analyzing coarse grained RNA dynamics based on Markov models of connectivity-graphs. Algorithms Mol Biol 2017; 12:15. [PMID: 28572834 PMCID: PMC5450175 DOI: 10.1186/s13015-017-0105-0] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/31/2016] [Accepted: 05/16/2017] [Indexed: 12/05/2022] Open
Abstract
Background In this work, we present a new coarse grained representation of RNA dynamics. It is based on adjacency matrices and their interactions patterns obtained from molecular dynamics simulations. RNA molecules are well-suited for this representation due to their composition which is mainly modular and assessable by the secondary structure alone. These interactions can be represented as adjacency matrices of k nucleotides. Based on those, we define transitions between states as changes in the adjacency matrices which form Markovian dynamics. The intense computational demand for deriving the transition probability matrices prompted us to develop StreAM-\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$T_g$$\end{document}Tg, a stream-based algorithm for generating such Markov models of k-vertex adjacency matrices representing the RNA. Results We benchmark StreAM-\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$T_g$$\end{document}Tg (a) for random and RNA unit sphere dynamic graphs (b) for the robustness of our method against different parameters. Moreover, we address a riboswitch design problem by applying StreAM-\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$T_g$$\end{document}Tg on six long term molecular dynamics simulation of a synthetic tetracycline dependent riboswitch (500 ns) in combination with five different antibiotics. Conclusions The proposed algorithm performs well on large simulated as well as real world dynamic graphs. Additionally, StreAM-\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$T_g$$\end{document}Tg provides insights into nucleotide based RNA dynamics in comparison to conventional metrics like the root-mean square fluctuation. In the light of experimental data our results show important design opportunities for the riboswitch.
Collapse
|
42
|
Wang J, Xiao Y. Using 3dRNA for RNA 3-D Structure Prediction and Evaluation. ACTA ACUST UNITED AC 2017; 57:5.9.1-5.9.12. [PMID: 28463400 DOI: 10.1002/cpbi.21] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
This unit describes how to use 3dRNA to predict RNA 3-D structures from their sequences and secondary (2-D) structures, and how to use 3dRNAscore to evaluate the predicted structures. The predicted RNA 3-D structures can be used to predict or understand their functions and can also be used to find the interactions between the RNA and other molecules. © 2017 by John Wiley & Sons, Inc.
Collapse
Affiliation(s)
- Jian Wang
- Biomolecular Physics and Modeling Group, School of Physics, Huazhong University of Science and Technology, Wuhan, Hubei, China
| | - Yi Xiao
- Biomolecular Physics and Modeling Group, School of Physics, Huazhong University of Science and Technology, Wuhan, Hubei, China
| |
Collapse
|
43
|
Jeddi I, Saiz L. Three-dimensional modeling of single stranded DNA hairpins for aptamer-based biosensors. Sci Rep 2017; 7:1178. [PMID: 28446765 PMCID: PMC5430850 DOI: 10.1038/s41598-017-01348-5] [Citation(s) in RCA: 96] [Impact Index Per Article: 13.7] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/20/2017] [Accepted: 03/24/2017] [Indexed: 12/11/2022] Open
Abstract
Aptamers consist of short oligonucleotides that bind specific targets. They provide advantages over antibodies, including robustness, low cost, and reusability. Their chemical structure allows the insertion of reporter molecules and surface-binding agents in specific locations, which have been recently exploited for the development of aptamer-based biosensors and direct detection strategies. Mainstream use of these devices, however, still requires significant improvements in optimization for consistency and reproducibility. DNA aptamers are more stable than their RNA counterparts for biomedical applications but have the disadvantage of lacking the wide array of computational tools for RNA structural prediction. Here, we present the first approach to predict from sequence the three-dimensional structures of single stranded (ss) DNA required for aptamer applications, focusing explicitly on ssDNA hairpins. The approach consists of a pipeline that integrates sequentially building ssDNA secondary structure from sequence, constructing equivalent 3D ssRNA models, transforming the 3D ssRNA models into ssDNA 3D structures, and refining the resulting ssDNA 3D structures. Through this pipeline, our approach faithfully predicts the representative structures available in the Nucleic Acid Database and Protein Data Bank databases. Our results, thus, open up a much-needed avenue for integrating DNA in the computational analysis and design of aptamer-based biosensors.
Collapse
Affiliation(s)
- Iman Jeddi
- Modeling of Biological Networks and Systems Therapeutics Laboratory, Department of Biomedical Engineering, University of California, 451 East Health Sciences Drive, Davis, CA, 95616, USA
| | - Leonor Saiz
- Modeling of Biological Networks and Systems Therapeutics Laboratory, Department of Biomedical Engineering, University of California, 451 East Health Sciences Drive, Davis, CA, 95616, USA.
| |
Collapse
|
44
|
Bell DR, Cheng SY, Salazar H, Ren P. Capturing RNA Folding Free Energy with Coarse-Grained Molecular Dynamics Simulations. Sci Rep 2017; 7:45812. [PMID: 28393861 PMCID: PMC5385882 DOI: 10.1038/srep45812] [Citation(s) in RCA: 37] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/15/2016] [Accepted: 03/06/2017] [Indexed: 01/25/2023] Open
Abstract
We introduce a coarse-grained RNA model for molecular dynamics simulations, RACER (RnA CoarsE-gRained). RACER achieves accurate native structure prediction for a number of RNAs (average RMSD of 2.93 Å) and the sequence-specific variation of free energy is in excellent agreement with experimentally measured stabilities (R2 = 0.93). Using RACER, we identified hydrogen-bonding (or base pairing), base stacking, and electrostatic interactions as essential driving forces for RNA folding. Also, we found that separating pairing vs. stacking interactions allowed RACER to distinguish folded vs. unfolded states. In RACER, base pairing and stacking interactions each provide an approximate stability of 3-4 kcal/mol for an A-form helix. RACER was developed based on PDB structural statistics and experimental thermodynamic data. In contrast with previous work, RACER implements a novel effective vdW potential energy function, which led us to re-parameterize hydrogen bond and electrostatic potential energy functions. Further, RACER is validated and optimized using a simulated annealing protocol to generate potential energy vs. RMSD landscapes. Finally, RACER is tested using extensive equilibrium pulling simulations (0.86 ms total) on eleven RNA sequences (hairpins and duplexes).
Collapse
Affiliation(s)
- David R. Bell
- Department of Biomedical Engineering, University of Texas at Austin, Austin, Texas 78712, United States
| | - Sara Y. Cheng
- Department of Physics, University of Texas at Austin, Austin, Texas 78712, United States
| | - Heber Salazar
- Department of Biomedical Engineering, University of Texas at Austin, Austin, Texas 78712, United States
| | - Pengyu Ren
- Department of Biomedical Engineering, University of Texas at Austin, Austin, Texas 78712, United States
| |
Collapse
|
45
|
Jedrzejczyk D, Gendaszewska-Darmach E, Pawlowska R, Chworos A. Designing synthetic RNA for delivery by nanoparticles. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2017; 29:123001. [PMID: 28004640 DOI: 10.1088/1361-648x/aa5561] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/06/2023]
Abstract
The rapid development of synthetic biology and nanobiotechnology has led to the construction of various synthetic RNA nanoparticles of different functionalities and potential applications. As they occur naturally, nucleic acids are an attractive construction material for biocompatible nanoscaffold and nanomachine design. In this review, we provide an overview of the types of RNA and nucleic acid's nanoparticle design, with the focus on relevant nanostructures utilized for gene-expression regulation in cellular models. Structural analysis and modeling is addressed along with the tools available for RNA structural prediction. The functionalization of RNA-based nanoparticles leading to prospective applications of such constructs in potential therapies is shown. The route from the nanoparticle design and modeling through synthesis and functionalization to cellular application is also described. For a better understanding of the fate of targeted RNA after delivery, an overview of RNA processing inside the cell is also provided.
Collapse
Affiliation(s)
- Dominika Jedrzejczyk
- Centre of Molecular and Macromolecular Studies, Polish Academy of Sciences, Sienkiewicza 112, 90-363 Lodz, Poland
| | | | | | | |
Collapse
|
46
|
Accurate Classification of RNA Structures Using Topological Fingerprints. PLoS One 2016; 11:e0164726. [PMID: 27755571 PMCID: PMC5068708 DOI: 10.1371/journal.pone.0164726] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/01/2016] [Accepted: 09/29/2016] [Indexed: 12/26/2022] Open
Abstract
While RNAs are well known to possess complex structures, functionally similar RNAs often have little sequence similarity. While the exact size and spacing of base-paired regions vary, functionally similar RNAs have pronounced similarity in the arrangement, or topology, of base-paired stems. Furthermore, predicted RNA structures often lack pseudoknots (a crucial aspect of biological activity), and are only partially correct, or incomplete. A topological approach addresses all of these difficulties. In this work we describe each RNA structure as a graph that can be converted to a topological spectrum (RNA fingerprint). The set of subgraphs in an RNA structure, its RNA fingerprint, can be compared with the fingerprints of other RNA structures to identify and correctly classify functionally related RNAs. Topologically similar RNAs can be identified even when a large fraction, up to 30%, of the stems are omitted, indicating that highly accurate structures are not necessary. We investigate the performance of the RNA fingerprint approach on a set of eight highly curated RNA families, with diverse sizes and functions, containing pseudoknots, and with little sequence similarity-an especially difficult test set. In spite of the difficult test set, the RNA fingerprint approach is very successful (ROC AUC > 0.95). Due to the inclusion of pseudoknots, the RNA fingerprint approach both covers a wider range of possible structures than methods based only on secondary structure, and its tolerance for incomplete structures suggests that it can be applied even to predicted structures. Source code is freely available at https://github.rcac.purdue.edu/mgribsko/XIOS_RNA_fingerprint.
Collapse
|
47
|
Shi YZ, Jin L, Wang FH, Zhu XL, Tan ZJ. Predicting 3D Structure, Flexibility, and Stability of RNA Hairpins in Monovalent and Divalent Ion Solutions. Biophys J 2016; 109:2654-2665. [PMID: 26682822 DOI: 10.1016/j.bpj.2015.11.006] [Citation(s) in RCA: 28] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/03/2015] [Revised: 10/09/2015] [Accepted: 11/06/2015] [Indexed: 10/24/2022] Open
Abstract
A full understanding of RNA-mediated biology would require the knowledge of three-dimensional (3D) structures, structural flexibility, and stability of RNAs. To predict RNA 3D structures and stability, we have previously proposed a three-bead coarse-grained predictive model with implicit salt/solvent potentials. In this study, we further develop the model by improving the implicit-salt electrostatic potential and including a sequence-dependent coaxial stacking potential to enable the model to simulate RNA 3D structure folding in divalent/monovalent ion solutions. The model presented here can predict 3D structures of RNA hairpins with bulges/internal loops (<77 nucleotides) from their sequences at the corresponding experimental ion conditions with an overall improved accuracy compared to the experimental data; the model also makes reliable predictions for the flexibility of RNA hairpins with bulge loops of different lengths at several divalent/monovalent ion conditions. In addition, the model successfully predicts the stability of RNA hairpins with various loops/stems in divalent/monovalent ion solutions.
Collapse
Affiliation(s)
- Ya-Zhou Shi
- Department of Physics and Key Laboratory of Artificial Micro- and Nano-structures of the Ministry of Education, School of Physics and Technology, Wuhan University, Wuhan, China
| | - Lei Jin
- Department of Physics and Key Laboratory of Artificial Micro- and Nano-structures of the Ministry of Education, School of Physics and Technology, Wuhan University, Wuhan, China
| | - Feng-Hua Wang
- Engineering Training Center, Jianghan University, Wuhan, China
| | - Xiao-Long Zhu
- Department of Physics, School of Physics and Information Engineering, Jianghan University, Wuhan, China
| | - Zhi-Jie Tan
- Department of Physics and Key Laboratory of Artificial Micro- and Nano-structures of the Ministry of Education, School of Physics and Technology, Wuhan University, Wuhan, China.
| |
Collapse
|
48
|
Zandi K, Butler G, Kharma N. An Adaptive Defect Weighted Sampling Algorithm to Design Pseudoknotted RNA Secondary Structures. Front Genet 2016; 7:129. [PMID: 27499762 PMCID: PMC4956659 DOI: 10.3389/fgene.2016.00129] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/06/2016] [Accepted: 07/06/2016] [Indexed: 01/18/2023] Open
Abstract
Computational design of RNA sequences that fold into targeted secondary structures has many applications in biomedicine, nanotechnology and synthetic biology. An RNA molecule is made of different types of secondary structure elements and an important RNA element named pseudoknot plays a key role in stabilizing the functional form of the molecule. However, due to the computational complexities associated with characterizing pseudoknotted RNA structures, most of the existing RNA sequence designer algorithms generally ignore this important structural element and therefore limit their applications. In this paper we present a new algorithm to design RNA sequences for pseudoknotted secondary structures. We use NUPACK as the folding algorithm to compute the equilibrium characteristics of the pseudoknotted RNAs, and describe a new adaptive defect weighted sampling algorithm named Enzymer to design low ensemble defect RNA sequences for targeted secondary structures including pseudoknots. We used a biological data set of 201 pseudoknotted structures from the Pseudobase library to benchmark the performance of our algorithm. We compared the quality characteristics of the RNA sequences we designed by Enzymer with the results obtained from the state of the art MODENA and antaRNA. Our results show our method succeeds more frequently than MODENA and antaRNA do, and generates sequences that have lower ensemble defect, lower probability defect and higher thermostability. Finally by using Enzymer and by constraining the design to a naturally occurring and highly conserved Hammerhead motif, we designed 8 sequences for a pseudoknotted cis-acting Hammerhead ribozyme. Enzymer is available for download at https://bitbucket.org/casraz/enzymer.
Collapse
Affiliation(s)
- Kasra Zandi
- Computer Science Department, Concordia UniversityMontreal, QC, Canada
| | - Gregory Butler
- Computer Science Department, Concordia UniversityMontreal, QC, Canada
- Centre for Structural and Functional Genomics, Concordia UniversityMontreal, QC, Canada
| | - Nawwaf Kharma
- Centre for Structural and Functional Genomics, Concordia UniversityMontreal, QC, Canada
- Electrical and Computer Engineering Department, Concordia UniversityMontreal, QC, Canada
| |
Collapse
|
49
|
Bottaro S, Gil-Ley A, Bussi G. RNA folding pathways in stop motion. Nucleic Acids Res 2016; 44:5883-91. [PMID: 27091499 PMCID: PMC4937309 DOI: 10.1093/nar/gkw239] [Citation(s) in RCA: 34] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/02/2016] [Revised: 03/29/2016] [Accepted: 03/29/2016] [Indexed: 11/19/2022] Open
Abstract
We introduce a method for predicting RNA folding pathways, with an application to the most important RNA tetraloops. The method is based on the idea that ensembles of three-dimensional fragments extracted from high-resolution crystal structures are heterogeneous enough to describe metastable as well as intermediate states. These ensembles are first validated by performing a quantitative comparison against available solution nuclear magnetic resonance (NMR) data of a set of RNA tetranucleotides. Notably, the agreement is better with respect to the one obtained by comparing NMR with extensive all-atom molecular dynamics simulations. We then propose a procedure based on diffusion maps and Markov models that makes it possible to obtain reaction pathways and their relative probabilities from fragment ensembles. This approach is applied to study the helix-to-loop folding pathway of all the tetraloops from the GNRA and UNCG families. The results give detailed insights into the folding mechanism that are compatible with available experimental data and clarify the role of intermediate states observed in previous simulation studies. The method is computationally inexpensive and can be used to study arbitrary conformational transitions.
Collapse
Affiliation(s)
- Sandro Bottaro
- Scuola Internazionale Superiore di Studi Avanzati, International School for Advanced Studies, 265, Via Bonomea I-34136 Trieste, Italy
| | - Alejandro Gil-Ley
- Scuola Internazionale Superiore di Studi Avanzati, International School for Advanced Studies, 265, Via Bonomea I-34136 Trieste, Italy
| | - Giovanni Bussi
- Scuola Internazionale Superiore di Studi Avanzati, International School for Advanced Studies, 265, Via Bonomea I-34136 Trieste, Italy
| |
Collapse
|
50
|
Baulin E, Yacovlev V, Khachko D, Spirin S, Roytberg M. URS DataBase: universe of RNA structures and their motifs. DATABASE-THE JOURNAL OF BIOLOGICAL DATABASES AND CURATION 2016; 2016:baw085. [PMID: 27242032 PMCID: PMC4885603 DOI: 10.1093/database/baw085] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 12/08/2015] [Accepted: 05/02/2016] [Indexed: 12/17/2022]
Abstract
The Universe of RNA Structures DataBase (URSDB) stores information obtained from all RNA-containing PDB entries (2935 entries in October 2015). The content of the database is updated regularly. The database consists of 51 tables containing indexed data on various elements of the RNA structures. The database provides a web interface allowing user to select a subset of structures with desired features and to obtain various statistical data for a selected subset of structures or for all structures. In particular, one can easily obtain statistics on geometric parameters of base pairs, on structural motifs (stems, loops, etc.) or on different types of pseudoknots. The user can also view and get information on an individual structure or its selected parts, e.g. RNA–protein hydrogen bonds. URSDB employs a new original definition of loops in RNA structures. That definition fits both pseudoknot-free and pseudoknotted secondary structures and coincides with the classical definition in case of pseudoknot-free structures. To our knowledge, URSDB is the first database supporting searches based on topological classification of pseudoknots and on extended loop classification. Database URL: http://server3.lpm.org.ru/urs/
Collapse
Affiliation(s)
- Eugene Baulin
- Laboratory of Applied Mathematics, Institute of Mathematical Problems of Biology, Russian Academy of Sciences, Pushchino, Moscow Region 142290, Russia Department of Algorithms and Technology of Programming, Faculty of Innovations and High Technology, Moscow Institute of Physics and Technology (State University), Dolgoprudny, Moscow Region 141700, Russia
| | - Victor Yacovlev
- Laboratory of Applied Mathematics, Institute of Mathematical Problems of Biology, Russian Academy of Sciences, Pushchino, Moscow Region 142290, Russia Department of Big Data and Information Retrieval, Faculty of Computer Science, National Research University Higher School of Economics, Moscow 101000, Russia
| | - Denis Khachko
- Laboratory of Applied Mathematics, Institute of Mathematical Problems of Biology, Russian Academy of Sciences, Pushchino, Moscow Region 142290, Russia
| | - Sergei Spirin
- Department of Mathematical Methods in Biology, Belozersky Institute of Physico-Chemical Biology, Lomonosov Moscow State University, Moscow 119992, Russia
| | - Mikhail Roytberg
- Laboratory of Applied Mathematics, Institute of Mathematical Problems of Biology, Russian Academy of Sciences, Pushchino, Moscow Region 142290, Russia Department of Algorithms and Technology of Programming, Faculty of Innovations and High Technology, Moscow Institute of Physics and Technology (State University), Dolgoprudny, Moscow Region 141700, Russia Department of Big Data and Information Retrieval, Faculty of Computer Science, National Research University Higher School of Economics, Moscow 101000, Russia
| |
Collapse
|