title: Isosteric and nonisosteric base pairs in RNA motifs: molecular dynamics and bioinformatics study of the sarcin-ricin internal loop.
abstract: The sarcin-ricin RNA motif (SR motif) is one of the most prominent recurrent RNA building blocks that occurs in many different RNA contexts and folds autonomously, that is, in a context-independent manner. In this study, we combined bioinformatics analysis with explicit-solvent molecular dynamics (MD) simulations to better understand the relation between the RNA sequence and the evolutionary patterns of the SR motif. A SHAPE probing experiment was also performed to confirm the fidelity of the MD simulations. We identified 57 instances of the SR motif in a nonredundant subset of the RNA X-ray structure database and analyzed their base pairing, base-phosphate, and backbone-backbone interactions. We extracted sequences aligned to these instances from large rRNA alignments to determine the frequency of occurrence for different sequence variants. We then used a simple scoring scheme based on isostericity to suggest 10 sequence variants with a highly variable expected degree of compatibility with the SR motif 3D structure. We carried out MD simulations of SR motifs with these base substitutions. Nonisosteric base substitutions led to unstable structures, but so did isosteric substitutions which were unable to make key base-phosphate interactions. The MD technique explains why some potentially isosteric SR motifs are not realized during evolution. We also found that the inability to form stable cWW geometry is an important factor in the case of the first base pair of the flexible region of the SR motif. A comparison of structural, bioinformatics, SHAPE probing, and MD simulation data reveals that explicit solvent MD simulations neatly reflect the viability of different sequence variants of the SR motif. Thus, MD simulations can efficiently complement bioinformatics tools in studies of conservation patterns of RNA motifs and provide atomistic insight into the role of their different signature interactions.
authors.affiliation: Institute of Biophysics, Academy of Sciences of the Czech Republic , Královopolská 135, 612 65 Brno, Czech Republic.
authors.forename: Marek
authors.initials: M
authors.lastname: Havrila
authors.forename: Kamila
authors.initials: K
authors.lastname: Réblová
authors.forename: Craig L
authors.initials: CL
authors.lastname: Zirbel
authors.forename: Neocles B
authors.initials: NB
authors.lastname: Leontis
authors.forename: Jiří
authors.initials: J
authors.lastname: Šponer
mesh_terms: D020029:Base Pairing; D019295:Computational Biology; D006860:Hydrogen Bonding; D056004:Molecular Dynamics Simulation; D009690:Nucleic Acid Conformation; D059372:Nucleotide Motifs; D012313:RNA / Q000737:chemistry* / Q000378:metabolism; D012335:RNA, Ribosomal / Q000737:chemistry / Q000378:metabolism; D012997:Solvents / Q000737:chemistry
chemical_list: D012335:RNA, Ribosomal; D012997:Solvents; D012313:RNA