Computational Investigation of RNA CUG Repeats Responsible for Myotonic Dystrophy 1

Myotonic Dystrophy 1 (DM1) is a genetic disease caused by expansion of CTG repeats in DNA. Once transcribed, these repeats form RNA hairpins with repeating 1×1 nucleotide UU internal loop motifs, r(CUG)n, which attract muscleblind-like 1 (MBNL1) protein leading to the disease. In DM1 CUG can be repeated thousands of times, so these structures are intractable to characterization using structural biology. However, inhibition of MBNL1-r(CUG)n binding requires a detailed analysis of the 1×1 UU internal loops. In this contribution we employ regular and umbrella sampling molecular dynamics (MD) simulations to describe the structural and thermodynamic properties of 1×1 UU internal loops. Calculations were run on a reported crystal structure and a designed system, which mimics an infinitely long RNA molecule with continuous CUG repeats. Two-dimensional (2D) potential of mean force (PMF) surfaces were created by umbrella sampling, and the discrete path sampling (DPS) method was utilized to investigate the energy landscape of 1×1 UU RNA internal loops, revealing that 1×1 UU base pairs are dynamic and strongly prefer the anti–anti conformation. Two 2D PMF surfaces were calculated for the 1×1 UU base pairs, revealing several local minima and three syn–anti ↔ anti–anti transformation pathways. Although at room temperature the syn–anti ↔ anti–anti transformation is not observed on the MD time scale, one of these pathways dominates the dynamics of the 1×1 UU base pairs in temperature jump MD simulations. This mechanism has now been treated successfully using the DPS approach. Our results suggest that local minima predicted by umbrella sampling calculations could be stabilized by small molecules, which is of great interest for future drug design. Furthermore, distorted GC/CG conformations may be important in understanding how MBNL1 binds to RNA CUG repeats. Hence we provide new insight into the dynamic roles of RNA loops and their contributions to presently incurable diseases.

[1]  David J Wales,et al.  Thermodynamics and kinetics of aggregation for the GNNQQNY peptide. , 2007, Journal of the American Chemical Society.

[2]  N. Bonini,et al.  Roles of trinucleotide-repeat RNA in neurological disease and degeneration , 2010, Trends in Neurosciences.

[3]  J. Doye,et al.  THE DOUBLE-FUNNEL ENERGY LANDSCAPE OF THE 38-ATOM LENNARD-JONES CLUSTER , 1998, cond-mat/9808265.

[4]  Jessica L. Childs-Disney,et al.  Rational design of ligands targeting triplet repeating transcripts that cause RNA dominant disease: application to myotonic muscular dystrophy type 1 and spinocerebellar ataxia type 3. , 2009, Journal of the American Chemical Society.

[5]  A. Young,et al.  Milestones in Huntington disease , 2011, Movement disorders : official journal of the Movement Disorder Society.

[6]  Y. Kino,et al.  Muscleblind protein, MBNL1/EXP, binds specifically to CHHG repeats. , 2004, Human molecular genetics.

[7]  Thomas Klockgether,et al.  Milestones in ataxia , 2011, Movement disorders : official journal of the Movement Disorder Society.

[8]  M. Swanson,et al.  Failure of MBNL1-dependent post-natal splicing transitions in myotonic dystrophy. , 2006, Human molecular genetics.

[9]  D. Turner,et al.  Optimization of an AMBER Force Field for the Artificial Nucleic Acid, LNA, and Benchmarking with NMR of L(CAAU) , 2013, The journal of physical chemistry. B.

[10]  Abbreviations , 1971 .

[11]  G. Ciccotti,et al.  Numerical Integration of the Cartesian Equations of Motion of a System with Constraints: Molecular Dynamics of n-Alkanes , 1977 .

[12]  Pamela L. Vanegas,et al.  RNA CoSSMos: Characterization of Secondary Structure Motifs—a searchable database of secondary structure motifs in RNA three-dimensional structures , 2011, Nucleic Acids Res..

[13]  M. Hoare,et al.  Statistical mechanics and morphology of very small atomic clusters , 1976 .

[14]  P. Hagerman,et al.  The (CGG)n repeat element within the 5' untranslated region of the FMR1 message provides both positive and negative cis effects on in vivo translation of a downstream reporter. , 2003, Human molecular genetics.

[15]  P. Jin,et al.  Pur α Binds to rCGG Repeats and Modulates Repeat-Mediated Neurodegeneration in a Drosophila Model of Fragile X Tremor/Ataxia Syndrome , 2007, Neuron.

[16]  M. Nadeau,et al.  Proteins : Structure , Function , and Bioinformatics , 2022 .

[17]  David J Wales,et al.  Folding pathways and rates for the three-stranded beta-sheet peptide Beta3s using discrete path sampling. , 2008, The journal of physical chemistry. B.

[18]  T. Cooper,et al.  Pathogenic mechanisms of myotonic dystrophy. , 2009, Biochemical Society transactions.

[19]  George C Schatz,et al.  A dynamic structural model of expanded RNA CAG repeats: a refined X-ray structure and computational investigations using molecular dynamics and umbrella sampling simulations. , 2013, Journal of the American Chemical Society.

[20]  D. Patel,et al.  Structural insights into RNA recognition by the alternative-splicing regulator muscleblind-like MBNL1 , 2008, Nature Structural &Molecular Biology.

[21]  Harry T Orr,et al.  Trinucleotide repeat disorders. , 2007, Annual review of neuroscience.

[22]  M. Karplus,et al.  Hidden complexity of free energy surfaces for peptide (protein) folding. , 2004, Proceedings of the National Academy of Sciences of the United States of America.

[23]  Mark A. Miller,et al.  Energy landscapes, structural topologies and rearrangement mechanisms in clusters of dipolar particles , 2013, 1304.8008.

[24]  K. Usdin,et al.  Fragile X syndrome and Friedreich’s ataxia: two different paradigms for repeat induced transcript insufficiency , 2001, Brain Research Bulletin.

[25]  M. Disney,et al.  NMR spectroscopy and molecular dynamics simulation of r(CCGCUGCGG)₂ reveal a dynamic UU internal loop found in myotonic dystrophy type 1. , 2011, Biochemistry.

[26]  Shankar Kumar,et al.  Multidimensional free‐energy calculations using the weighted histogram analysis method , 1995, J. Comput. Chem..

[27]  M. Carmo-Fonseca,et al.  Triplet repeats, RNA secondary structure and toxic gain-of-function models for pathogenesis , 2001, Brain Research Bulletin.

[28]  David J. Wales,et al.  Free energy surfaces from an extended harmonic superposition approach and kinetics for alanine dipeptide , 2008 .

[29]  Warf Mb,et al.  MBNL binds similar RNA structures in the CUG repeats of myotonic dystrophy and its pre-mRNA substrate cardiac troponin T , 2007 .

[30]  Harry A. Stern,et al.  Revision of AMBER Torsional Parameters for RNA Improves Free Energy Predictions for Tetramer Duplexes with GC and iGiC Base Pairs , 2011, Journal of chemical theory and computation.

[31]  B. Oostra,et al.  Characterization and localization of the FMR-1 gene product associated with fragile X syndrome , 1993, Nature.

[32]  J. Berglund,et al.  Utilizing the GAAA tetraloop/receptor to facilitate crystal packing and determination of the structure of a CUG RNA helix. , 2012, Biochemistry.

[33]  J. Sutcliffe,et al.  Identification of a gene (FMR-1) containing a CGG repeat coincident with a breakpoint cluster region exhibiting length variation in fragile X syndrome , 1991, Cell.

[34]  J. Šponer,et al.  Refinement of the AMBER Force Field for Nucleic Acids: Improving the Description of α/γ Conformers , 2007 .

[35]  I. Kanazawa,et al.  Abnormal gene product identified in hereditary dentatorubral–pallidoluysian atrophy (DRPLA) brain , 1995, Nature Genetics.

[36]  David J Wales,et al.  Calculating rate constants and committor probabilities for transition networks by graph transformation. , 2009, The Journal of chemical physics.

[37]  Angel E García,et al.  High-resolution reversible folding of hyperstable RNA tetraloops using molecular dynamics simulations , 2013, Proceedings of the National Academy of Sciences.

[38]  Harry A. Stern,et al.  Reparameterization of RNA χ Torsion Parameters for the AMBER Force Field and Comparison to NMR Spectra for Cytidine and Uridine , 2010, Journal of chemical theory and computation.

[39]  L. Malinina,et al.  Structural dynamics of double-helical RNAs composed of CUG/CUG- and CUG/CGG-repeats , 2012, Journal of biomolecular structure & dynamics.

[40]  T. Cheatham,et al.  Determination of Alkali and Halide Monovalent Ion Parameters for Use in Explicitly Solvated Biomolecular Simulations , 2008, The journal of physical chemistry. B.

[41]  M. Searle,et al.  Structure and folding dynamics of a DNA hairpin with a stabilising d(GNA) trinucleotide loop: influence of base pair mis-matches and point mutations on conformational equilibria. , 2007, Organic & biomolecular chemistry.

[42]  S. Richard,et al.  Sam68 sequestration and partial loss of function are associated with splicing alterations in FXTAS patients , 2010, The EMBO journal.

[43]  S. Sugano,et al.  Structural basis for the sequence-specific RNA-recognition mechanism of human CUG-BP1 RRM3 , 2009, Nucleic acids research.

[44]  M. Karplus,et al.  The topology of multidimensional potential energy surfaces: Theory and application to peptide structure and kinetics , 1997 .

[45]  Peter Sanders Algorithm Engineering , 2010, Informatik-Spektrum.

[46]  Peter Salamon,et al.  Observation time scale, free-energy landscapes, and molecular symmetry , 2013, Proceedings of the National Academy of Sciences.

[47]  Brian B. Gibbens,et al.  Non-ATG–initiated translation directed by microsatellite expansions , 2010, Proceedings of the National Academy of Sciences.

[48]  Eric T. Wang,et al.  Design of a bioactive small molecule that targets the myotonic dystrophy type 1 RNA via an RNA motif-ligand database and chemical similarity searching. , 2012, Journal of the American Chemical Society.

[49]  David J. Wales,et al.  Some further applications of discrete path sampling to cluster isomerization , 2004 .

[50]  Mark A. Miller,et al.  Archetypal energy landscapes , 1998, Nature.

[51]  S. Warren,et al.  Trinucleotide repeat expansion and human disease. , 1995, Annual review of genetics.

[52]  宁北芳,et al.  疟原虫var基因转换速率变化导致抗原变异[英]/Paul H, Robert P, Christodoulou Z, et al//Proc Natl Acad Sci U S A , 2005 .

[53]  Jessica L. Childs-Disney,et al.  Rationally designed small molecules targeting the RNA that causes myotonic dystrophy type 1 are potently bioactive. , 2012, ACS chemical biology.

[54]  Benjamin L Miller,et al.  Dynamic combinatorial selection of molecules capable of inhibiting the (CUG) repeat RNA-MBNL1 interaction in vitro: discovery of lead compounds targeting myotonic dystrophy (DM1). , 2008, Journal of the American Chemical Society.

[55]  Juan J. Marugan,et al.  Induction and Reversal of Myotonic Dystrophy Type 1 Pre-mRNA Splicing Defects by Small Molecules , 2013, Nature Communications.

[56]  Jacek Krol,et al.  Structural Diversity of Triplet Repeat RNAs*♦ , 2010, The Journal of Biological Chemistry.

[57]  S. Mirkin,et al.  Expansion of the (CTG)(n) repeat in the 5'-UTR of a reporter gene impedes translation. , 2000, Nucleic acids research.

[58]  B. Roux The calculation of the potential of mean force using computer simulations , 1995 .

[59]  F. Walker Huntington's disease , 2007, The Lancet.

[60]  D. Turner,et al.  Benchmarking AMBER Force Fields for RNA: Comparisons to NMR Spectra for Single-Stranded r(GACC) Are Improved by Revised χ Torsions , 2011, The journal of physical chemistry. B.

[61]  T. Cooper,et al.  Aberrant regulation of insulin receptor alternative splicing is associated with insulin resistance in myotonic dystrophy , 2001, Nature Genetics.

[62]  Jacek Krol,et al.  RNA structure of trinucleotide repeats associated with human neurological diseases. , 2003, Nucleic acids research.

[63]  Alexander D. MacKerell,et al.  Free energy and structural pathways of base flipping in a DNA GCGC containing sequence. , 2002, Journal of molecular biology.

[64]  J Grigsby,et al.  Neuropathology of fragile X-associated tremor/ataxia syndrome (FXTAS). , 2006, Brain : a journal of neurology.

[65]  Jorge Nocedal,et al.  On the limited memory BFGS method for large scale optimization , 1989, Math. Program..

[66]  M. Nakamori,et al.  From dynamic combinatorial ‘hit’ to lead: in vitro and in vivo activity of compounds targeting the pathogenic RNAs that cause myotonic dystrophy , 2012, Nucleic acids research.

[67]  D. Turner,et al.  The Nuclear Magnetic Resonance of CCCC RNA Reveals a Right-Handed Helix, and Revised Parameters for AMBER Force Field Torsions Improve Structural Predictions from Molecular Dynamics , 2013, Biochemistry.

[68]  George C Schatz,et al.  Structure of the myotonic dystrophy type 2 RNA and designed small molecules that reduce toxicity. , 2014, ACS chemical biology.

[69]  Sergei V. Krivov,et al.  Free energy disconnectivity graphs: Application to peptide models , 2002 .

[70]  Xiang-Jun Lu,et al.  3DNA: a versatile, integrated software system for the analysis, rebuilding and visualization of three-dimensional nucleic-acid structures , 2008, Nature Protocols.

[71]  A. Baranger,et al.  MBNL1–RNA Recognition: Contributions of MBNL1 Sequence and RNA Conformation , 2012, Chembiochem : a European journal of chemical biology.

[72]  Sergei V Krivov,et al.  One-dimensional free-energy profiles of complex systems: progress variables that preserve the barriers. , 2006, The journal of physical chemistry. B.

[73]  P. Kollman,et al.  A Second Generation Force Field for the Simulation of Proteins, Nucleic Acids, and Organic Molecules , 1995 .

[74]  R. Kierzek,et al.  Structural insights into CUG repeats containing the ‘stretched U–U wobble’: implications for myotonic dystrophy , 2009, Nucleic acids research.

[75]  Michal Otyepka,et al.  Performance of Molecular Mechanics Force Fields for RNA Simulations: Stability of UUCG and GNRA Hairpins. , 2010, Journal of chemical theory and computation.

[76]  D. Turner,et al.  Expanded CUG repeat RNAs form hairpins that activate the double-stranded RNA-dependent protein kinase PKR. , 2000, RNA.

[77]  David J Wales,et al.  Folding of the GB1 hairpin peptide from discrete path sampling. , 2004, The Journal of chemical physics.

[78]  R. Rosenfeld Nature , 2009, Otolaryngology--head and neck surgery : official journal of American Academy of Otolaryngology-Head and Neck Surgery.

[79]  B. Byrne,et al.  Recruitment of human muscleblind proteins to (CUG)n expansions associated with myotonic dystrophy , 2000, The EMBO journal.

[80]  D. Wales,et al.  Energy landscapes and global thermodynamics for alanine peptides. , 2013, The Journal of chemical physics.

[81]  A. Emery,et al.  The muscular dystrophies , 2002, The Lancet.

[82]  Lennart Nilsson,et al.  Conformational Preferences of Modified Uridines: Comparison of AMBER Derived Force Fields , 2014, J. Chem. Inf. Model..

[83]  Matthew D Disney,et al.  Chemical correction of pre-mRNA splicing defects associated with sequestration of muscleblind-like 1 protein by expanded r(CAG)-containing transcripts. , 2012, ACS chemical biology.

[84]  Faraday Discuss , 1985 .

[85]  J. Berglund,et al.  The four Zn fingers of MBNL1 provide a flexible platform for recognition of its RNA binding elements , 2011, BMC Molecular Biology.

[86]  W. L. Jorgensen,et al.  Comparison of simple potential functions for simulating liquid water , 1983 .

[87]  N. Bonini,et al.  RNA toxicity is a component of ataxin-3 degeneration in Drosophila , 2008, Nature.

[88]  A. Roses,et al.  Polyglutamine Domains Are Substrates of Tissue Transglutaminase: Does Transglutaminase Play a Role in Expanded CAG/Poly‐Q Neurodegenerative Diseases? , 1997, Journal of neurochemistry.

[89]  F. Young Biochemistry , 1955, The Indian Medical Gazette.

[90]  D. Wales Discrete path sampling , 2002 .

[91]  D. Wales,et al.  Energy landscapes, folding mechanisms, and kinetics of RNA tetraloop hairpins. , 2014, Journal of the American Chemical Society.

[92]  B. Schoser,et al.  Expression of RNA CCUG repeats dysregulates translation and degradation of proteins in myotonic dystrophy 2 patients. , 2009, American Journal of Pathology.

[93]  K Schulten,et al.  VMD: visual molecular dynamics. , 1996, Journal of molecular graphics.

[94]  M. Swanson,et al.  Muscleblind-like 1 interacts with RNA hairpins in splicing target and pathogenic RNAs , 2007, Nucleic acids research.

[95]  Andrés Marzal,et al.  Computing the K Shortest Paths: A New Algorithm and an Experimental Comparison , 1999, WAE.

[96]  S. Rice,et al.  ADVANCES IN CHEMICAL PHYSICS , 2002 .

[97]  David J Wales,et al.  Evolution of the potential energy landscape with static pulling force for two model proteins. , 2012, The journal of physical chemistry. B.

[98]  J. Lueck,et al.  Muscle Chloride Channel Dysfunction in Two Mouse Models of Myotonic Dystrophy , 2007, The Journal of general physiology.

[99]  A. Baranger,et al.  A simple ligand that selectively targets CUG trinucleotide repeats and inhibits MBNL protein binding , 2009, Proceedings of the National Academy of Sciences.

[100]  Michal Otyepka,et al.  Computer Folding of RNA Tetraloops? Are We There Yet? , 2013, Journal of chemical theory and computation.

[101]  J. Lueck,et al.  Correction of ClC-1 splicing eliminates chloride channelopathy and myotonia in mouse models of myotonic dystrophy. , 2007, The Journal of clinical investigation.

[102]  Alexander D. MacKerell,et al.  Computational approaches for investigating base flipping in oligonucleotides. , 2006, Chemical reviews.

[103]  M. Goodman,et al.  Analysis of Strand Slippage in DNA Polymerase Expansions of CAG/CTG Triplet Repeats Associated with Neurodegenerative Disease* , 1998, The Journal of Biological Chemistry.

[104]  K. Sobczak,et al.  CAG repeats mimic CUG repeats in the misregulation of alternative splicing , 2011, Nucleic acids research.

[105]  M. Disney,et al.  Myotonic dystrophy type 1 RNA crystal structures reveal heterogeneous 1 × 1 nucleotide UU internal loop conformations. , 2011, Biochemistry.

[106]  S. Delaney,et al.  Structure of even/odd trinucleotide repeat sequences modulates persistence of non-B conformations and conversion to duplex. , 2011, Biochemistry.

[107]  C. Simmerling,et al.  An Improved Reaction Coordinate for Nucleic Acid Base Flipping Studies. , 2009, Journal of chemical theory and computation.