1
|
Krepl M, Vögele J, Kruse H, Duchardt-Ferner E, Wöhnert J, Sponer J. An intricate balance of hydrogen bonding, ion atmosphere and dynamics facilitates a seamless uracil to cytosine substitution in the U-turn of the neomycin-sensing riboswitch. Nucleic Acids Res 2018; 46:6528-6543. [PMID: 29893898 PMCID: PMC6061696 DOI: 10.1093/nar/gky490] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/20/2018] [Accepted: 05/23/2018] [Indexed: 01/04/2023] Open
Abstract
The neomycin sensing riboswitch is the smallest biologically functional RNA riboswitch, forming a hairpin capped with a U-turn loop-a well-known RNA motif containing a conserved uracil. It was shown previously that a U→C substitution of the eponymous conserved uracil does not alter the riboswitch structure due to C protonation at N3. Furthermore, cytosine is evolutionary permitted to replace uracil in other U-turns. Here, we use molecular dynamics simulations to study the molecular basis of this substitution in the neomycin sensing riboswitch and show that a structure-stabilizing monovalent cation-binding site in the wild-type RNA is the main reason for its negligible structural effect. We then use NMR spectroscopy to confirm the existence of this cation-binding site and to demonstrate its effects on RNA stability. Lastly, using quantum chemical calculations, we show that the cation-binding site is altering the electronic environment of the wild-type U-turn so that it is more similar to the cytosine mutant. The study reveals an amazingly complex and delicate interplay between various energy contributions shaping up the 3D structure and evolution of nucleic acids.
Collapse
Affiliation(s)
- Miroslav Krepl
- Institute of Biophysics of the Czech Academy of Sciences, Kralovopolska 135, 612 65 Brno, Czech Republic
| | - Jennifer Vögele
- Institute of Molecular Biosciences and Center for Biomolecular Magnetic Resonance (BMRZ), Goethe-University Frankfurt, Max-von-Laue-Str. 9, 60438 Frankfurt, Germany
| | - Holger Kruse
- Institute of Biophysics of the Czech Academy of Sciences, Kralovopolska 135, 612 65 Brno, Czech Republic
| | - Elke Duchardt-Ferner
- Institute of Molecular Biosciences and Center for Biomolecular Magnetic Resonance (BMRZ), Goethe-University Frankfurt, Max-von-Laue-Str. 9, 60438 Frankfurt, Germany
| | - Jens Wöhnert
- Institute of Molecular Biosciences and Center for Biomolecular Magnetic Resonance (BMRZ), Goethe-University Frankfurt, Max-von-Laue-Str. 9, 60438 Frankfurt, Germany
| | - Jiri Sponer
- Institute of Biophysics of the Czech Academy of Sciences, Kralovopolska 135, 612 65 Brno, Czech Republic
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University Olomouc, 17. listopadu 12, 771 46 Olomouc, Czech Republic
| |
Collapse
|
2
|
Chattopadhyay A, Zheng M, Waller MP, Priyakumar UD. A Probabilistic Framework for Constructing Temporal Relations in Replica Exchange Molecular Trajectories. J Chem Theory Comput 2018; 14:3365-3380. [PMID: 29791153 DOI: 10.1021/acs.jctc.7b01245] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Knowledge of the structure and dynamics of biomolecules is essential for elucidating the underlying mechanisms of biological processes. Given the stochastic nature of many biological processes, like protein unfolding, it is almost impossible that two independent simulations will generate the exact same sequence of events, which makes direct analysis of simulations difficult. Statistical models like Markov chains, transition networks, etc. help in shedding some light on the mechanistic nature of such processes by predicting long-time dynamics of these systems from short simulations. However, such methods fall short in analyzing trajectories with partial or no temporal information, for example, replica exchange molecular dynamics or Monte Carlo simulations. In this work, we propose a probabilistic algorithm, borrowing concepts from graph theory and machine learning, to extract reactive pathways from molecular trajectories in the absence of temporal data. A suitable vector representation was chosen to represent each frame in the macromolecular trajectory (as a series of interaction and conformational energies), and dimensionality reduction was performed using principal component analysis (PCA). The trajectory was then clustered using a density-based clustering algorithm, where each cluster represents a metastable state on the potential energy surface (PES) of the biomolecule under study. A graph was created with these clusters as nodes with the edges learned using an iterative expectation maximization algorithm. The most reactive path is conceived as the widest path along this graph. We have tested our method on RNA hairpin unfolding trajectory in aqueous urea solution. Our method makes the understanding of the mechanism of unfolding in the RNA hairpin molecule more tractable. As this method does not rely on temporal data, it can be used to analyze trajectories from Monte Carlo sampling techniques and replica exchange molecular dynamics (REMD).
Collapse
Affiliation(s)
- Aditya Chattopadhyay
- Centre for Computational Natural Sciences and Bioinformatics , International Institute of Information Technology , Hyderabad 500032 , India
| | - Min Zheng
- Centre for Multiscale Theory and Computation , Westfälische Wilhelms-Universität Münster , Münster , Germany
| | - Mark P Waller
- Department of Physics and International Centre for Quantum and Molecular Structures , Shanghai University , Shanghai , 200444 , People's Republic of China
| | - U Deva Priyakumar
- Centre for Computational Natural Sciences and Bioinformatics , International Institute of Information Technology , Hyderabad 500032 , India
| |
Collapse
|
3
|
Šponer J, Bussi G, Krepl M, Banáš P, Bottaro S, Cunha RA, Gil-Ley A, Pinamonti G, Poblete S, Jurečka P, Walter NG, Otyepka M. RNA Structural Dynamics As Captured by Molecular Simulations: A Comprehensive Overview. Chem Rev 2018; 118:4177-4338. [PMID: 29297679 PMCID: PMC5920944 DOI: 10.1021/acs.chemrev.7b00427] [Citation(s) in RCA: 366] [Impact Index Per Article: 52.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/14/2017] [Indexed: 12/14/2022]
Abstract
With both catalytic and genetic functions, ribonucleic acid (RNA) is perhaps the most pluripotent chemical species in molecular biology, and its functions are intimately linked to its structure and dynamics. Computer simulations, and in particular atomistic molecular dynamics (MD), allow structural dynamics of biomolecular systems to be investigated with unprecedented temporal and spatial resolution. We here provide a comprehensive overview of the fast-developing field of MD simulations of RNA molecules. We begin with an in-depth, evaluatory coverage of the most fundamental methodological challenges that set the basis for the future development of the field, in particular, the current developments and inherent physical limitations of the atomistic force fields and the recent advances in a broad spectrum of enhanced sampling methods. We also survey the closely related field of coarse-grained modeling of RNA systems. After dealing with the methodological aspects, we provide an exhaustive overview of the available RNA simulation literature, ranging from studies of the smallest RNA oligonucleotides to investigations of the entire ribosome. Our review encompasses tetranucleotides, tetraloops, a number of small RNA motifs, A-helix RNA, kissing-loop complexes, the TAR RNA element, the decoding center and other important regions of the ribosome, as well as assorted others systems. Extended sections are devoted to RNA-ion interactions, ribozymes, riboswitches, and protein/RNA complexes. Our overview is written for as broad of an audience as possible, aiming to provide a much-needed interdisciplinary bridge between computation and experiment, together with a perspective on the future of the field.
Collapse
Affiliation(s)
- Jiří Šponer
- Institute of Biophysics of the Czech Academy of Sciences , Kralovopolska 135 , Brno 612 65 , Czech Republic
| | - Giovanni Bussi
- Scuola Internazionale Superiore di Studi Avanzati , Via Bonomea 265 , Trieste 34136 , Italy
| | - Miroslav Krepl
- Institute of Biophysics of the Czech Academy of Sciences , Kralovopolska 135 , Brno 612 65 , Czech Republic
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science , Palacky University Olomouc , 17. listopadu 12 , Olomouc 771 46 , Czech Republic
| | - Pavel Banáš
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science , Palacky University Olomouc , 17. listopadu 12 , Olomouc 771 46 , Czech Republic
| | - Sandro Bottaro
- Structural Biology and NMR Laboratory, Department of Biology , University of Copenhagen , Copenhagen 2200 , Denmark
| | - Richard A Cunha
- Scuola Internazionale Superiore di Studi Avanzati , Via Bonomea 265 , Trieste 34136 , Italy
| | - Alejandro Gil-Ley
- Scuola Internazionale Superiore di Studi Avanzati , Via Bonomea 265 , Trieste 34136 , Italy
| | - Giovanni Pinamonti
- Scuola Internazionale Superiore di Studi Avanzati , Via Bonomea 265 , Trieste 34136 , Italy
| | - Simón Poblete
- Scuola Internazionale Superiore di Studi Avanzati , Via Bonomea 265 , Trieste 34136 , Italy
| | - Petr Jurečka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science , Palacky University Olomouc , 17. listopadu 12 , Olomouc 771 46 , Czech Republic
| | - Nils G Walter
- Single Molecule Analysis Group and Center for RNA Biomedicine, Department of Chemistry , University of Michigan , Ann Arbor , Michigan 48109 , United States
| | - Michal Otyepka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science , Palacky University Olomouc , 17. listopadu 12 , Olomouc 771 46 , Czech Republic
| |
Collapse
|
4
|
Sripathi KN, Banáš P, Réblová K, Šponer J, Otyepka M, Walter NG. Wobble pairs of the HDV ribozyme play specific roles in stabilization of active site dynamics. Phys Chem Chem Phys 2015; 17:5887-900. [PMID: 25631765 DOI: 10.1039/c4cp05083e] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
The hepatitis delta virus (HDV) is the only known human pathogen whose genome contains a catalytic RNA motif (ribozyme). The overall architecture of the HDV ribozyme is that of a double-nested pseudoknot, with two GU pairs flanking the active site. Although extensive studies have shown that mutation of either wobble results in decreased catalytic activity, little work has focused on linking these mutations to specific structural effects on catalytic fitness. Here we use molecular dynamics simulations based on an activated structure to probe the active site dynamics as a result of wobble pair mutations. In both wild-type and mutant ribozymes, the in-line fitness of the active site (as a measure of catalytic proficiency) strongly depends on the presence of a C75(N3H3+)N1(O5') hydrogen bond, which positions C75 as the general acid for the reaction. Our mutational analyses show that each GU wobble supports catalytically fit conformations in distinct ways; the reverse G25U20 wobble promotes high in-line fitness, high occupancy of the C75(N3H3+)G1(O5') general-acid hydrogen bond and stabilization of the G1U37 wobble, while the G1U37 wobble acts more locally by stabilizing high in-line fitness and the C75(N3H3+)G1(O5') hydrogen bond. We also find that stable type I A-minor and P1.1 hydrogen bonding above and below the active site, respectively, prevent local structural disorder from spreading and disrupting global conformation. Taken together, our results define specific, often redundant architectural roles for several structural motifs of the HDV ribozyme active site, expanding the known roles of these motifs within all HDV-like ribozymes and other structured RNAs.
Collapse
Affiliation(s)
- Kamali N Sripathi
- Department of Medicinal Chemistry, University of Michigan, Ann Arbor, Michigan 48109-1065, USA
| | | | | | | | | | | |
Collapse
|
5
|
Kellerman DL, Simmons KS, Pedraza M, Piccirilli JA, York DM, Harris ME. Determination of hepatitis delta virus ribozyme N(-1) nucleobase and functional group specificity using internal competition kinetics. Anal Biochem 2015; 483:12-20. [PMID: 25937290 DOI: 10.1016/j.ab.2015.04.024] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/18/2015] [Revised: 04/17/2015] [Accepted: 04/23/2015] [Indexed: 12/16/2022]
Abstract
Biological catalysis involves interactions distant from the site of chemistry that can position the substrate for reaction. Catalysis of RNA 2'-O-transphosphorylation by the hepatitis delta virus (HDV) ribozyme is sensitive to the identity of the N(-1) nucleotide flanking the reactive phosphoryl group. However, the interactions that affect the conformation of this position, and in turn the 2'O nucleophile, are unclear. Here, we describe the application of multiple substrate internal competition kinetic analyses to understand how the N(-1) nucleobase contributes to HDV catalysis and test the utility of this approach for RNA structure-function studies. Internal competition reactions containing all four substrate sequence variants at the N(-1) position in reactions using ribozyme active site mutations at A77 and A78 were used to test a proposed base-pairing interaction. Mutants A78U, A78G, and A79G retain significant catalytic activity but do not alter the specificity for the N(-1) nucleobase. Effects of nucleobase analog substitutions at N(-1) indicate that U is preferred due to the ability to donate an H-bond in the Watson-Crick face and avoid minor groove steric clash. The results provide information essential for evaluating models of the HDV active site and illustrate multiple substrate kinetic analyses as a practical approach for characterizing structure-function relationships in RNA reactions.
Collapse
Affiliation(s)
- Daniel L Kellerman
- Department of Biochemistry, Case Western Reserve University School of Medicine, Cleveland, OH 44106, USA
| | - Kandice S Simmons
- Department of Biochemistry, Case Western Reserve University School of Medicine, Cleveland, OH 44106, USA
| | - Mayra Pedraza
- Department of Biochemistry, Case Western Reserve University School of Medicine, Cleveland, OH 44106, USA
| | - Joseph A Piccirilli
- Department of Chemistry and Department of Biochemistry and Molecular Biology, University of Chicago, Chicago, IL 60637, USA
| | - Darrin M York
- Center for Integrative Proteomics Research, BioMaPS Institute for Quantitative Biology and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, NJ 08854, USA
| | - Michael E Harris
- Department of Biochemistry, Case Western Reserve University School of Medicine, Cleveland, OH 44106, USA.
| |
Collapse
|
6
|
Sripathi KN, Tay WW, Banáš P, Otyepka M, Šponer J, Walter NG. Disparate HDV ribozyme crystal structures represent intermediates on a rugged free-energy landscape. RNA (NEW YORK, N.Y.) 2014; 20:1112-28. [PMID: 24854621 PMCID: PMC4114689 DOI: 10.1261/rna.044982.114] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 02/23/2014] [Accepted: 04/16/2014] [Indexed: 06/03/2023]
Abstract
The hepatitis delta virus (HDV) ribozyme is a member of the class of small, self-cleaving catalytic RNAs found in a wide range of genomes from HDV to human. Both pre- and post-catalysis (precursor and product) crystal structures of the cis-acting genomic HDV ribozyme have been determined. These structures, together with extensive solution probing, have suggested that a significant conformational change accompanies catalysis. A recent crystal structure of a trans-acting precursor, obtained at low pH and by molecular replacement from the previous product conformation, conforms to the product, raising the possibility that it represents an activated conformer past the conformational change. Here, using fluorescence resonance energy transfer (FRET), we discovered that cleavage of this ribozyme at physiological pH is accompanied by a structural lengthening in magnitude comparable to previous trans-acting HDV ribozymes. Conformational heterogeneity observed by FRET in solution appears to have been removed upon crystallization. Analysis of a total of 1.8 µsec of molecular dynamics (MD) simulations showed that the crystallographically unresolved cleavage site conformation is likely correctly modeled after the hammerhead ribozyme, but that crystal contacts and the removal of several 2'-oxygens near the scissile phosphate compromise catalytic in-line fitness. A cis-acting version of the ribozyme exhibits a more dynamic active site, while a G-1 residue upstream of the scissile phosphate favors poor fitness, allowing us to rationalize corresponding changes in catalytic activity. Based on these data, we propose that the available crystal structures of the HDV ribozyme represent intermediates on an overall rugged RNA folding free-energy landscape.
Collapse
Affiliation(s)
- Kamali N. Sripathi
- Department of Medicinal Chemistry, University of Michigan, Ann Arbor, Michigan 48109-1065, USA
| | - Wendy W. Tay
- Program in Chemical Biology, University of Michigan, Ann Arbor, Michigan 48109-1055, USA
| | - Pavel Banáš
- Regional Centre of Advance Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University, 771 46 Olomouc, Czech Republic
| | - Michal Otyepka
- Regional Centre of Advance Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University, 771 46 Olomouc, Czech Republic
| | - Jiří Šponer
- Institute of Biophysics, Academy of Sciences of the Czech Republic, 612 65 Brno, Czech Republic
- Masaryk University, Campus Bohunice, 625 00 Brno, Czech Republic
| | - Nils G. Walter
- Department of Chemistry, Single Molecule Analysis Group, University of Michigan, Ann Arbor, Michigan 48109-1055, USA
| |
Collapse
|
7
|
Weigand JE, Gottstein-Schmidtke SR, Demolli S, Groher F, Duchardt-Ferner E, Wöhnert J, Suess B. Sequence elements distal to the ligand binding pocket modulate the efficiency of a synthetic riboswitch. Chembiochem 2014; 15:1627-37. [PMID: 24954073 DOI: 10.1002/cbic.201402067] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/01/2014] [Indexed: 01/16/2023]
Abstract
Synthetic riboswitches can serve as sophisticated genetic control devices in synthetic biology, regulating gene expression through direct RNA-ligand interactions. We analyzed a synthetic neomycin riboswitch, which folds into a stem loop structure with an internal loop important for ligand binding and regulation. It is closed by a terminal hexaloop containing a U-turn and a looped-out adenine. We investigated the relationship between sequence, structure, and biological activity in the terminal loop by saturating mutagenesis, ITC, and NMR. Mutants corresponding to the canonical U-turn fold retained biological activity. An improvement of stacking interactions in the U-turn led to an RNA element with slightly enhanced regulatory activity. For the first position of the U-turn motif and the looped out base, sequence-activity relationships that could not initially be explained on the basis of the structure of the aptamer-ligand complex were observed. However, NMR studies of these mutants revealed subtle relationships between structure and dynamics of the aptamer in its free or bound state and biological activity.
Collapse
Affiliation(s)
- Julia E Weigand
- Department of Biology, Technical University Darmstadt, Schnittspahnstrasse 10, 64287 Darmstadt (Germany)
| | | | | | | | | | | | | |
Collapse
|
8
|
Riccitelli N, Lupták A. HDV family of self-cleaving ribozymes. PROGRESS IN MOLECULAR BIOLOGY AND TRANSLATIONAL SCIENCE 2014; 120:123-71. [PMID: 24156943 DOI: 10.1016/b978-0-12-381286-5.00004-4] [Citation(s) in RCA: 31] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
The hepatitis delta virus (HDV) ribozymes are catalytic RNAs capable of cleaving their own sugar-phosphate backbone. The HDV virus possesses the ribozymes in both sense and antisense genomic transcripts, where they are essential for processing during replication. These ribozymes have been the subject of intense biochemical scrutiny and have yielded a wealth of mechanistic insights. In recent years, many HDV-like ribozymes have been identified in nearly all branches of life. The ribozymes are implicated in a variety of biological events, including episodic memory in mammals and retrotransposition in many eukaryotes. Detailed analysis of additional HDV-like ribozyme isolates will likely reveal many more biological functions and provide information about the evolution of this unique RNA.
Collapse
Affiliation(s)
- Nathan Riccitelli
- Department of Chemistry, University of California, Irvine, California, USA
| | | |
Collapse
|
9
|
Riccitelli NJ, Delwart E, Lupták A. Identification of minimal HDV-like ribozymes with unique divalent metal ion dependence in the human microbiome. Biochemistry 2014; 53:1616-26. [PMID: 24555915 DOI: 10.1021/bi401717w] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023]
Abstract
HDV-like self-cleaving ribozymes have been found in a wide variety of organisms, implicated in diverse biological processes, and their activity typically shows a strong divalent metal dependence, but little metal specificity. Recent studies suggested that very short variants of these ribozymes exist in nature, but their distribution and biochemical properties have not been established. To map out the distribution of small HDV-like ribozymes, the drz-Spur-3 sequence was minimized to yield a core construct for structure-based bioinformatic searches. These searches revealed several microbial ribozymes, particularly in the human microbiome. Kinetic profile of the smallest ribozyme revealed two distinct metal binding sites, only one of which promotes fast catalysis. Furthermore, this ribozyme showed markedly reduced activity in Ca(2+), even in the presence of physiological Mg(2+) concentrations. Our study substantially expands the number of microbial HDV-like ribozymes and provides an example of cleavage regulation by divalent metals.
Collapse
Affiliation(s)
- Nathan J Riccitelli
- Department of Chemistry, ∥Department of Pharmaceutical Sciences, and ⊥Department of Molecular Biology and Biochemistry, University of California-Irvine , Irvine, California 92697, United States
| | | | | |
Collapse
|
10
|
Abstract
This chapter gives an overview over the current methods for automated modeling of RNA structures, with emphasis on template-based methods. The currently used approaches to RNA modeling are presented with a side view on the protein world, where many similar ideas have been used. Two main programs for automated template-based modeling are presented: ModeRNA assembling structures from fragments and MacroMoleculeBuilder performing a simulation to satisfy spatial restraints. Both approaches have in common that they require an alignment of the target sequence to a known RNA structure that is used as a modeling template. As a way to find promising template structures and to align the target and template sequences, we propose a pipeline combining the ParAlign and Infernal programs on RNA family data from Rfam. We also briefly summarize template-free methods for RNA 3D structure prediction. Typically, RNA structures generated by automated modeling methods require local or global optimization. Thus, we also discuss methods that can be used for local or global refinement of RNA structures.
Collapse
Affiliation(s)
- Kristian Rother
- Laboratory of Bioinformatics and Protein Engineering, International Institute of Molecular and Cell Biology, Warsaw, Poland,
| | | | | | | |
Collapse
|
11
|
Bouchard P, Legault P. Structural insights into substrate recognition by the Neurospora Varkud satellite ribozyme: importance of U-turns at the kissing-loop junction. Biochemistry 2013; 53:258-69. [PMID: 24325625 PMCID: PMC3893828 DOI: 10.1021/bi401491g] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023]
Abstract
![]()
Substrate
recognition by the Neurospora Varkud
satellite ribozyme depends on the formation of a magnesium-dependent
kissing-loop interaction between the stem-loop I (SLI) substrate and
stem-loop V (SLV) of the catalytic domain. From mutagenesis studies,
it has been established that this I/V kissing-loop interaction involves
three Watson–Crick base pairs and is associated with a structural
rearrangement of the SLI substrate that facilitates catalysis. Here,
we report the NMR structural characterization of this I/V kissing-loop
using isolated stem-loops. NMR studies were performed on different
SLI/SLV complexes containing a common SLV and shiftable, preshifted,
or double-stranded SLI variants. These studies confirm the presence
of three Watson–Crick base pairs at the kissing-loop junction
and provide evidence for the structural rearrangement of shiftable
SLI variants upon SLV binding. NMR structure determination of an SLI/SLV
complex demonstrates that both the SLI and SLV loops adopt U-turn
structures, which facilitates intermolecular Watson–Crick base
pairing. Several other interactions at the I/V interface, including
base triples and base stacking, help create a continuously stacked
structure. These NMR studies provide a structural basis to understand
the stability of the I/V kissing-loop interaction and lead us to propose
a kinetic model for substrate activation in the VS ribozyme.
Collapse
Affiliation(s)
- Patricia Bouchard
- Département de Biochimie et Médecine Moléculaire, Université de Montréal , C.P. 6128, Succursale Centre-Ville, Montréal, Quebec H3C 3J7, Canada
| | | |
Collapse
|
12
|
Todd GC, Walter NG. Secondary structure of bacteriophage T4 gene 60 mRNA: implications for translational bypassing. RNA (NEW YORK, N.Y.) 2013; 19:685-700. [PMID: 23492219 PMCID: PMC3677283 DOI: 10.1261/rna.037291.112] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/13/2023]
Abstract
Translational bypassing is a unique phenomenon of bacteriophage T4 gene 60 mRNA wherein the bacterial ribosome produces a single polypeptide chain from a discontinuous open reading frame (ORF). Upon reaching the 50-nucleotide untranslated region, or coding gap, the ribosome either dissociates or bypasses the interruption to continue translating the remainder of the ORF, generating a subunit of a type II DNA topoisomerase. Mutational and computational analyses have suggested that a compact structure, including a stable hairpin, forms in the coding gap to induce bypassing, yet direct evidence is lacking. Here we have probed the secondary structure of gene 60 mRNA with both Tb³⁺ ions and the selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) reagent 1M7 under conditions where bypassing is observed. The resulting experimentally informed secondary structure models strongly support the presence of the predicted coding gap hairpin and highlight the benefits of using Tb³⁺ as a second, complementary probing reagent. Contrary to several previously proposed models, however, the rest of the coding gap is highly reactive with both probing reagents, suggesting that it forms only a short stem-loop. Mutational analyses coupled with functional assays reveal that two possible base-pairings of the coding gap with other regions of the mRNA are not required for bypassing. Such structural autonomy of the coding gap is consistent with its recently discovered role as a mobile genetic element inserted into gene 60 mRNA to inhibit cleavage by homing endonuclease MobA.
Collapse
Affiliation(s)
- Gabrielle C. Todd
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109-1055, USA
- Program in Chemical Biology, University of Michigan, Ann Arbor, Michigan 48109-1055, USA
| | - Nils G. Walter
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109-1055, USA
- Corresponding authorE-mail
| |
Collapse
|
13
|
Lee TS, Giambaşu G, Harris ME, York DM. Characterization of the Structure and Dynamics of the HDV Ribozyme at Different Stages Along the Reaction Path. J Phys Chem Lett 2011; 2:2538-2543. [PMID: 22200005 PMCID: PMC3244300 DOI: 10.1021/jz201106y] [Citation(s) in RCA: 25] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/28/2023]
Abstract
The structure and dynamics of the hepatitis delta virus ribozyme (HDVr) are studies using molecular dynamics simulations at several stages along its catalytic reaction path, including reactant, activated precursor, transition state mimic and product states, departing from an initial structure based on the C75U mutant crystal structure (PDB: 1VC7). Results of five 350 ns molecular dynamics simulations reveal a spontaneous rotation of U-1 that leads to an in-line conformation and support the role of protonated C75 as the general acid in the transition state. Our results provide rationale for the interpretation of several important experimental results, and make experimentally testable predictions regarding the roles of key active site residues that are not obvious from any available crystal structures.
Collapse
|
14
|
McDowell SE, Jun JM, Walter NG. Long-range tertiary interactions in single hammerhead ribozymes bias motional sampling toward catalytically active conformations. RNA (NEW YORK, N.Y.) 2010; 16:2414-2426. [PMID: 20921269 PMCID: PMC2995402 DOI: 10.1261/rna.1829110] [Citation(s) in RCA: 24] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/16/2009] [Accepted: 08/30/2010] [Indexed: 05/29/2023]
Abstract
Enzymes generally are thought to derive their functional activity from conformational motions. The limited chemical variation in RNA suggests that such structural dynamics may play a particularly important role in RNA function. Minimal hammerhead ribozymes are known to cleave efficiently only in ∼ 10-fold higher than physiologic concentrations of Mg(2+) ions. Extended versions containing native loop-loop interactions, however, show greatly enhanced catalytic activity at physiologically relevant Mg(2+) concentrations, for reasons that are still ill-understood. Here, we use Mg(2+) titrations, activity assays, ensemble, and single molecule fluorescence resonance energy transfer (FRET) approaches, combined with molecular dynamics (MD) simulations, to ask what influence the spatially distant tertiary loop-loop interactions of an extended hammerhead ribozyme have on its structural dynamics. By comparing hammerhead variants with wild-type, partially disrupted, and fully disrupted loop-loop interaction sequences we find that the tertiary interactions lead to a dynamic motional sampling that increasingly populates catalytically active conformations. At the global level the wild-type tertiary interactions lead to more frequent, if transient, encounters of the loop-carrying stems, whereas at the local level they lead to an enrichment in favorable in-line attack angles at the cleavage site. These results invoke a linkage between RNA structural dynamics and function and suggest that loop-loop interactions in extended hammerhead ribozymes-and Mg(2+) ions that bind to minimal ribozymes-may generally allow more frequent access to a catalytically relevant conformation(s), rather than simply locking the ribozyme into a single active state.
Collapse
|
15
|
Efficient inhibition of hepatitis B virus replication by hepatitis delta virus ribozymes delivered by targeting retrovirus. Virol J 2010; 7:61. [PMID: 20236514 PMCID: PMC2850903 DOI: 10.1186/1743-422x-7-61] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/20/2009] [Accepted: 03/17/2010] [Indexed: 12/14/2022] Open
Abstract
Background Hepatitis delta virus (HDV) ribozyme is an attractive molecular tool that can specifically recognize and catalyze the self-cleavage of the viral RNA phosphodiester backbone. However, a major obstacle in the medical application of the HDV ribozyme is the lack of specificity in the delivery of the ribozyme to defined target cells. Results The objective of this study was to determine whether retroviral vectors can deliver the HDV ribozyme into the target cells and to elucidate whether HDV ribozyme plays a role in hepatitis B virus (HBV) replication. In our study, the transduction of helper-free pseudotyped retrovirus, which showed a broad host range, in human hepatoma cells was performed under 2 conditions, that is, in the presence of polymerized human serum albumin (pHSA) and in the absence of pHSA. The transduction ability in the presence of pHSA was higher than in the absence of pHSA. Moreover, HBsAg and HBeAg levels after transductions with pHSA were significantly lower than those in the absence of pHSA, thus indicating that the recombinant retrovirus had HBV-specific cleavage activity and targeted HepG2215 cells. Conclusions These data suggest that this system provides a new approach for targeting hepatocytes and has a great potential in gene therapy for HBV infection.
Collapse
|
16
|
Ditzler MA, Otyepka M, Šponer J, Walter NG. Molecular dynamics and quantum mechanics of RNA: conformational and chemical change we can believe in. Acc Chem Res 2010; 43:40-7. [PMID: 19754142 PMCID: PMC2808146 DOI: 10.1021/ar900093g] [Citation(s) in RCA: 142] [Impact Index Per Article: 9.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
Abstract
Structure and dynamics are both critical to RNA’s vital functions in biology. Numerous techniques can elucidate the structural dynamics of RNA, but computational approaches based on experimental data arguably hold the promise of providing the most detail. In this Account, we highlight areas wherein molecular dynamics (MD) and quantum mechanical (QM) techniques are applied to RNA, particularly in relation to complementary experimental studies.
We have expanded on atomic-resolution crystal structures of RNAs in functionally relevant states by applying explicit solvent MD simulations to explore their dynamics and conformational changes on the submicrosecond time scale. MD relies on simplified atomistic, pairwise additive interaction potentials (force fields). Because of limited sampling, due to the finite accessible simulation time scale and the approximated force field, high-quality starting structures are required. Despite their imperfection, we find that currently available force fields empower MD to provide meaningful and predictive information on RNA dynamics around a crystallographically defined energy minimum. The performance of force fields can be estimated by precise QM calculations on small model systems. Such calculations agree reasonably well with the Cornell et al. AMBER force field, particularly for stacking and hydrogen-bonding interactions. A final verification of any force field is accomplished by simulations of complex nucleic acid structures. The performance of the Cornell et al. AMBER force field generally corresponds well with and augments experimental data, but one notable exception could be the capping loops of double-helical stems. In addition, the performance of pairwise additive force fields is obviously unsatisfactory for inclusion of divalent cations, because their interactions lead to major polarization and charge-transfer effects neglected by the force field. Neglect of polarization also limits, albeit to a lesser extent, the description accuracy of other contributions, such as interactions with monovalent ions, conformational flexibility of the anionic sugar−phosphate backbone, hydrogen bonding, and solute polarization by solvent. Still, despite limitations, MD simulations are a valid tool for analyzing the structural dynamics of existing experimental structures. Careful analysis of MD simulations can identify problematic aspects of an experimental RNA structure, unveil structural characteristics masked by experimental constraints, reveal functionally significant stochastic fluctuations, evaluate the structural role of base ionization, and predict structurally and potentially functionally important details of the solvent behavior, including the presence of tightly bound water molecules. Moreover, combining classical MD simulations with QM calculations in hybrid QM/MM approaches helps in the assessment of the plausibility of chemical mechanisms of catalytic RNAs (ribozymes). In contrast, the reliable prediction of structure from sequence information is beyond the applicability of MD tools. The ultimate utility of computational studies in understanding RNA function thus requires that the results are neither blindly accepted nor flatly rejected, but rather considered in the context of all available experimental data, with great care given to assessing limitations through the available starting structures, force field approximations, and sampling limitations. The examples given in this Account showcase how the judicious use of basic MD simulations has already served as a powerful tool to help evaluate the role of structural dynamics in biological function of RNA.
Collapse
Affiliation(s)
- Mark A. Ditzler
- Biophysics, University of Michigan, 930 North University Avenue, Ann Arbor, Michigan 48109-1055
- Department of Chemistry, Single Molecule Analysis Group, University of Michigan, 930 North University Avenue, Ann Arbor, Michigan 48109-1055
| | - Michal Otyepka
- Department of Physical Chemistry, Faculty of Science, Palacky University Olomouc, tr. Svobody 26, 771 46 Olomouc, Czech Republic
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Kralovopolska 135, 612 65 Brno, Czech Republic
| | - Jiřì Šponer
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Kralovopolska 135, 612 65 Brno, Czech Republic
| | - Nils G. Walter
- Department of Chemistry, Single Molecule Analysis Group, University of Michigan, 930 North University Avenue, Ann Arbor, Michigan 48109-1055
| |
Collapse
|
17
|
Ditzler MA, Sponer J, Walter NG. Molecular dynamics suggest multifunctionality of an adenine imino group in acid-base catalysis of the hairpin ribozyme. RNA (NEW YORK, N.Y.) 2009; 15:560-75. [PMID: 19223444 PMCID: PMC2661834 DOI: 10.1261/rna.1416709] [Citation(s) in RCA: 32] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/08/2023]
Abstract
Despite numerous structural and biochemical investigations, the catalytic mechanism of hairpin ribozyme self-cleavage remains elusive. To gain insight into the coupling of active site dynamics with activity of this small catalytic RNA, we analyzed a total of approximately 300 ns of molecular dynamics (MD) simulations. Our simulations predict improved global stability for an in vitro selected "gain of function" mutation, which is validated by native gel electrophoretic mobility shift assay. We observe that active site nucleobases and water molecules stabilize a geometry favorable to catalysis through a dynamic hydrogen bonding network. Simulations in which A38 is unprotonated show its N1 move into close proximity of the active site 2'-OH, indicating that A38 may act as a general base during cleavage, a role that has generally been discounted due to the longer distances observed in crystal structures involving inactivating substrate analogs. By contrast, simulations in which N1 of A38 is protonated place N1 in close proximity to the 5'-oxygen leaving group, which supports the proposal that A38 serves as a general acid. In analogy to protein enzymes, we discuss a plausible mechanism in which A38 acts bifunctionally and shuttles a proton directly from the 2'-OH to the 5'-oxygen. Furthermore, our simulations suggest an important role for protonation of N1 of A38 in promoting a favorable geometry similar to that observed in transition-state analog crystal structures, and support previously proposed roles of A38, G8, and long residency water molecules in transition-state stabilization.
Collapse
|
18
|
Smedile A, Ciancio A, Rizzetto M. Hepatitis D Virus. CLINICAL VIROLOGY 2009:1291-1306. [DOI: 10.1128/9781555815981.ch56] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/02/2025]
|
19
|
Walter NG, Perumal S. The Small Ribozymes: Common and Diverse Features Observed through the FRET Lens. SPRINGER SERIES IN BIOPHYSICS 2009; 13:103-127. [PMID: 21796234 DOI: 10.1007/978-3-540-70840-7_5] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
The hammerhead, hairpin, HDV, VS and glmS ribozymes are the five known, naturally occurring catalytic RNAs classified as the "small ribozymes". They share common reaction chemistry in cleaving their own backbone by phosphodiester transfer, but are diverse in their secondary and tertiary structures, indicating that Nature has found at least five independent solutions to a common chemical task. Fluorescence resonance energy transfer (FRET) has been extensively used to detect conformational changes in these ribozymes and dissect their reaction pathways. Common and diverse features are beginning to emerge that, by extension, highlight general biophysical properties of non-protein coding RNAs.
Collapse
Affiliation(s)
- Nils G Walter
- Department of Chemistry, Single Molecule Analysis Group, University of Michigan, Ann Arbor, MI 48109
| | | |
Collapse
|
20
|
Abstract
Many recently discovered RNA functions rely on highly complex multistep conformational transitions that occur in response to an array of cellular signals. These dynamics accompany and guide, for example, RNA cotranscriptional folding, ligand sensing and signaling, site-specific catalysis in ribozymes, and the hierarchically ordered assembly of ribonucleoproteins. RNA dynamics are encoded by both the inherent properties of RNA structure, spanning many motional modes with a large range of amplitudes and timescales, and external trigger factors, ranging from proteins, nucleic acids, metal ions, metabolites, and vitamins to temperature and even directional RNA biosynthesis itself. Here, we review recent advances in our understanding of RNA dynamics as highlighted by biophysical tools.
Collapse
|
21
|
Shao J, Tanner SW, Thompson N, Cheatham TE. Clustering Molecular Dynamics Trajectories: 1. Characterizing the Performance of Different Clustering Algorithms. J Chem Theory Comput 2007; 3:2312-34. [DOI: 10.1021/ct700119m] [Citation(s) in RCA: 614] [Impact Index Per Article: 34.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Jianyin Shao
- Departments of Medicinal Chemistry, Pharmaceutics and Pharmaceutical Chemistry, and Bioengineering, College of Pharmacy, University of Utah, 2000 East 30 South, Skaggs Hall 201, Salt Lake City, Utah 84112
| | - Stephen W. Tanner
- Departments of Medicinal Chemistry, Pharmaceutics and Pharmaceutical Chemistry, and Bioengineering, College of Pharmacy, University of Utah, 2000 East 30 South, Skaggs Hall 201, Salt Lake City, Utah 84112
| | - Nephi Thompson
- Departments of Medicinal Chemistry, Pharmaceutics and Pharmaceutical Chemistry, and Bioengineering, College of Pharmacy, University of Utah, 2000 East 30 South, Skaggs Hall 201, Salt Lake City, Utah 84112
| | - Thomas E. Cheatham
- Departments of Medicinal Chemistry, Pharmaceutics and Pharmaceutical Chemistry, and Bioengineering, College of Pharmacy, University of Utah, 2000 East 30 South, Skaggs Hall 201, Salt Lake City, Utah 84112
| |
Collapse
|
22
|
McDowell SE, Špačková N, Šponer J, Walter NG. Molecular dynamics simulations of RNA: an in silico single molecule approach. Biopolymers 2007; 85:169-84. [PMID: 17080418 PMCID: PMC2018183 DOI: 10.1002/bip.20620] [Citation(s) in RCA: 114] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
RNA molecules are now known to be involved in the processing of genetic information at all levels, taking on a wide variety of central roles in the cell. Understanding how RNA molecules carry out their biological functions will require an understanding of structure and dynamics at the atomistic level, which can be significantly improved by combining computational simulation with experiment. This review provides a critical survey of the state of molecular dynamics (MD) simulations of RNA, including a discussion of important current limitations of the technique and examples of its successful application. Several types of simulations are discussed in detail, including those of structured RNA molecules and their interactions with the surrounding solvent and ions, catalytic RNAs, and RNA-small molecule and RNA-protein complexes. Increased cooperation between theorists and experimentalists will allow expanded judicious use of MD simulations to complement conceptually related single molecule experiments. Such cooperation will open the door to a fundamental understanding of the structure-function relationships in diverse and complex RNA molecules. .
Collapse
Affiliation(s)
- S. Elizabeth McDowell
- Biophysics Research Division, Single Molecule Analysis Group, University of Michigan, 930 N. University Avenue, Ann Arbor, MI 48109-1055
| | - Nad'a Špačková
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Královopolská 135, 612 65 Brno, Czech Republic
| | - Jiří Šponer
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Královopolská 135, 612 65 Brno, Czech Republic
- Institute of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic, Prague
| | - Nils G. Walter
- Department of Chemistry, Single Molecule Analysis Group, University of Michigan, 930 N. University Avenue, Ann Arbor, MI 48109-1055
| |
Collapse
|
23
|
Tinsley RA, Walter NG. Long-range impact of peripheral joining elements on structure and function of the hepatitis delta virus ribozyme. Biol Chem 2007; 388:705-15. [PMID: 17570823 DOI: 10.1515/bc.2007.088] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/22/2023]
Abstract
The HDV ribozyme is an RNA enzyme from the human pathogenic hepatitis delta virus (HDV) that has recently also been identified in the human genome. It folds into a compact, nested double-pseudoknot. We examined here the functional relevance of the capping loop L4 and the helical crossover J1/2, which tightly interlace the two helical stacks of the ribozyme. Peripheral structural elements such as these are present in cis-acting, but not trans-acting ribozymes, which may explain the order-of-magnitude decrease in cleavage activity observed in trans-acting ribozymes with promise in gene therapy applications. Comparison of a systematic set of cis- and trans-acting HDV ribozymes shows that the absence of either L4 or J1/2 significantly and independently impacts catalytic activity. Using terbium(III) footprinting and affinity studies, as well as distance measurements based on time-resolved fluorescence resonance energy transfer, we find that J1/2 is most important for conferring structural properties similar to those of the cis-acting ribozyme. Our results are consistent with a model in which removal of either a helical crossover or surprisingly a capping loop induces greater dynamics and expansion of the catalytic core at long range, impacting local and global folding, as well as catalytic function.
Collapse
Affiliation(s)
- Rebecca A Tinsley
- Department of Chemistry, University of Michigan, Ann Arbor, MI 48109-1055, USA
| | | |
Collapse
|