1
|
Ichijo R, Kamimura T, Kawai G. Interaction between a fluoroquinolone derivative KG022 and RNAs: Effect of base pairs 3′ adjacent to the bulged residues. Front Mol Biosci 2023; 10:1145528. [PMID: 36999159 PMCID: PMC10043337 DOI: 10.3389/fmolb.2023.1145528] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/16/2023] [Accepted: 02/27/2023] [Indexed: 03/18/2023] Open
Abstract
RNA-targeted small molecules are a promising modality in drug discovery. Recently, we found that a fluoroquinolone derivative, KG022, can bind to RNAs with bulged C or G. To clarify the RNA specificity of KG022, we analyzed the effect of the base pair located at the 3′side of the bulged residue. It was found that KG022 prefers G-C and A-U base pairs at the 3′side. Solution structures of the complexes of KG022 with the four RNA molecules with bulged C or G and G-C or A-U base pairs at the 3′side of the bulged residue were determined to find that the fluoroquinolone moiety is located between two purine bases, and this may be the mechanism of the specificity. This work provides an important example of the specificity of RNA-targeted small molecules.
Collapse
Affiliation(s)
- Rika Ichijo
- Graduate School of Engineering, Chiba Institute of Technology, Chiba, Japan
| | | | - Gota Kawai
- Graduate School of Engineering, Chiba Institute of Technology, Chiba, Japan
- *Correspondence: Gota Kawai,
| |
Collapse
|
2
|
Bae JH, Zhang DY. Predicting stability of DNA bulge at mononucleotide microsatellite. Nucleic Acids Res 2021; 49:7901-7908. [PMID: 34308470 PMCID: PMC8373066 DOI: 10.1093/nar/gkab616] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/02/2021] [Revised: 06/28/2021] [Accepted: 07/07/2021] [Indexed: 11/14/2022] Open
Abstract
Mononucleotide microsatellites are clinically and forensically crucial DNA sequences due to their high mutability and abundance in the human genome. As a mutagenic intermediate of an indel in a microsatellite and a consequence of probe hybridization after such mutagenesis, a bulge with structural degeneracy sliding within a microsatellite is formed. Stability of such dynamic bulges, however, is still poorly understood despite their critical role in cancer genomics and neurological disease studies. In this paper, we have built a model that predicts the thermodynamics of a sliding bulge at a microsatellite. We first identified 40 common bulge states that can be assembled into any sliding bulges, and then characterized them with toehold exchange energy measurement and the partition function. Our model, which is the first to predict the free energy of sliding bulges with more than three repeats, can infer the stability penalty of a sliding bulge of any sequence and length with a median prediction error of 0.22 kcal/mol. Patterns from the prediction clearly explain landscapes of microsatellites observed in the literature, such as higher mutation rates of longer microsatellites and C/G microsatellites.
Collapse
Affiliation(s)
- Jin H Bae
- Department of Bioengineering, Rice University, Houston, TX 77005, USA
| | - David Yu Zhang
- Department of Bioengineering, Rice University, Houston, TX 77005, USA.,Systems, Synthetic, and Physical Biology, Rice University, Houston, TX 77005, USA
| |
Collapse
|
3
|
Mathews DH. How to benchmark RNA secondary structure prediction accuracy. Methods 2019; 162-163:60-67. [PMID: 30951834 DOI: 10.1016/j.ymeth.2019.04.003] [Citation(s) in RCA: 21] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/31/2018] [Revised: 03/24/2019] [Accepted: 04/01/2019] [Indexed: 11/18/2022] Open
Abstract
RNA secondary structure prediction is widely used. As new methods are developed, these are often benchmarked for accuracy against existing methods. This review discusses good practices for performing these benchmarks, including the choice of benchmarking structures, metrics to quantify accuracy, the importance of allowing flexibility for pairs in the accepted structure, and the importance of statistical testing for significance.
Collapse
Affiliation(s)
- David H Mathews
- Center for RNA Biology, Department of Biochemistry & Biophysics, and Department of Biostatistics & Computational Biology, University of Rochester Medical Center, 601 Elmwood Avenue, Box 712, Rochester, NY 14642, United States.
| |
Collapse
|
4
|
Zuber J, Cabral BJ, McFadyen I, Mauger DM, Mathews DH. Analysis of RNA nearest neighbor parameters reveals interdependencies and quantifies the uncertainty in RNA secondary structure prediction. RNA (NEW YORK, N.Y.) 2018; 24:1568-1582. [PMID: 30104207 PMCID: PMC6191722 DOI: 10.1261/rna.065102.117] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/28/2017] [Accepted: 08/07/2018] [Indexed: 05/08/2023]
Abstract
RNA secondary structure prediction is often used to develop hypotheses about structure-function relationships for newly discovered RNA sequences, to identify unknown functional RNAs, and to design sequences. Secondary structure prediction methods typically use a thermodynamic model that estimates the free energy change of possible structures based on a set of nearest neighbor parameters. These parameters were derived from optical melting experiments of small model oligonucleotides. This work aims to better understand the precision of structure prediction. Here, the experimental errors in optical melting experiments were propagated to errors in the derived nearest neighbor parameter values and then to errors in RNA secondary structure prediction. To perform this analysis, the optical melting experimental values were systematically perturbed within the estimates of experimental error and alternative sets of nearest neighbor parameters were then derived from these error-bounded values. Secondary structure predictions using either the perturbed or reference parameter sets were then compared. This work demonstrated that the precision of RNA secondary structure prediction is more robust than suggested by previous work based on perturbation of the nearest neighbor parameters. This robustness is due to correlations between parameters. Additionally, this work identified weaknesses in the parameter derivation that makes accurate assessment of parameter uncertainty difficult. Considerations for experimental design are provided to mitigate these weaknesses are provided.
Collapse
Affiliation(s)
- Jeffrey Zuber
- Department of Biochemistry and Biophysics and Center for RNA Biology, University of Rochester Medical Center, Rochester, New York 14642, USA
| | - B Joseph Cabral
- Computational Sciences, Moderna Therapeutics, Cambridge, Massachusetts 02141, USA
| | - Iain McFadyen
- Computational Sciences, Moderna Therapeutics, Cambridge, Massachusetts 02141, USA
| | - David M Mauger
- Computational Sciences, Moderna Therapeutics, Cambridge, Massachusetts 02141, USA
| | - David H Mathews
- Department of Biochemistry and Biophysics and Center for RNA Biology, University of Rochester Medical Center, Rochester, New York 14642, USA
- Department of Biostatistics and Computational Biology, University of Rochester Medical Center, Rochester, New York 14642, USA
| |
Collapse
|
5
|
Patil KM, Toh DFK, Yuan Z, Meng Z, Shu Z, Zhang H, Ong A, Krishna MS, Lu L, Lu Y, Chen G. Incorporating uracil and 5-halouracils into short peptide nucleic acids for enhanced recognition of A-U pairs in dsRNAs. Nucleic Acids Res 2018; 46:7506-7521. [PMID: 30011039 PMCID: PMC6125629 DOI: 10.1093/nar/gky631] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/05/2018] [Revised: 06/09/2018] [Accepted: 07/03/2018] [Indexed: 01/16/2023] Open
Abstract
Double-stranded RNA (dsRNA) structures form triplexes and RNA-protein complexes through binding to single-stranded RNA (ssRNA) regions and proteins, respectively, for diverse biological functions. Hence, targeting dsRNAs through major-groove triplex formation is a promising strategy for the development of chemical probes and potential therapeutics. Short (e.g., 6-10 mer) chemically-modified Peptide Nucleic Acids (PNAs) have been developed that bind to dsRNAs sequence specifically at physiological conditions. For example, a PNA incorporating a modified base thio-pseudoisocytosine (L) has an enhanced recognition of a G-C pair in an RNA duplex through major-groove L·G-C base triple formation at physiological pH, with reduced pH dependence as observed for C+·G-C base triple formation. Currently, an unmodified T base is often incorporated into PNAs to recognize a Watson-Crick A-U pair through major-groove T·A-U base triple formation. A substitution of the 5-methyl group in T by hydrogen and halogen atoms (F, Cl, Br, and I) causes a decrease of the pKa of N3 nitrogen atom, which may result in improved hydrogen bonding in addition to enhanced base stacking interactions. Here, we synthesized a series of PNAs incorporating uracil and halouracils, followed by binding studies by non-denaturing polyacrylamide gel electrophoresis, circular dichroism, and thermal melting. Our results suggest that replacing T with uracil and halouracils may enhance the recognition of an A-U pair by PNA·RNA2 triplex formation in a sequence-dependent manner, underscoring the importance of local stacking interactions. Incorporating bromouracils and chlorouracils into a PNA results in a significantly reduced pH dependence of triplex formation even for PNAs containing C bases, likely due to an upshift of the apparent pKa of N3 atoms of C bases. Thus, halogenation and other chemical modifications may be utilized to enhance hydrogen bonding of the adjacent base triples and thus triplex formation. Furthermore, our experimental and computational modelling data suggest that PNA·RNA2 triplexes may be stabilized by incorporating a BrUL step but not an LBrU step, in dsRNA-binding PNAs.
Collapse
Affiliation(s)
- Kiran M Patil
- Division of Chemistry and Biological Chemistry, School of Physical and Mathematical Sciences, Nanyang Technological University, 21 Nanyang Link, Singapore 637371
| | - Desiree-Faye Kaixin Toh
- Division of Chemistry and Biological Chemistry, School of Physical and Mathematical Sciences, Nanyang Technological University, 21 Nanyang Link, Singapore 637371
| | - Zhen Yuan
- Division of Chemistry and Biological Chemistry, School of Physical and Mathematical Sciences, Nanyang Technological University, 21 Nanyang Link, Singapore 637371
| | - Zhenyu Meng
- Division of Chemistry and Biological Chemistry, School of Physical and Mathematical Sciences, Nanyang Technological University, 21 Nanyang Link, Singapore 637371
| | - Zhiyu Shu
- Division of Chemistry and Biological Chemistry, School of Physical and Mathematical Sciences, Nanyang Technological University, 21 Nanyang Link, Singapore 637371
| | - Haiping Zhang
- School of Biological Sciences, Nanyang Technological University, 60 Nanyang Drive, Singapore 637551
| | - Alan Ann Lerk Ong
- Division of Chemistry and Biological Chemistry, School of Physical and Mathematical Sciences, Nanyang Technological University, 21 Nanyang Link, Singapore 637371
| | - Manchugondanahalli S Krishna
- Division of Chemistry and Biological Chemistry, School of Physical and Mathematical Sciences, Nanyang Technological University, 21 Nanyang Link, Singapore 637371
| | - Lanyuan Lu
- School of Biological Sciences, Nanyang Technological University, 60 Nanyang Drive, Singapore 637551
| | - Yunpeng Lu
- Division of Chemistry and Biological Chemistry, School of Physical and Mathematical Sciences, Nanyang Technological University, 21 Nanyang Link, Singapore 637371
| | - Gang Chen
- Division of Chemistry and Biological Chemistry, School of Physical and Mathematical Sciences, Nanyang Technological University, 21 Nanyang Link, Singapore 637371
| |
Collapse
|
6
|
de Oliveira Martins E, Weber G. An asymmetric mesoscopic model for single bulges in RNA. J Chem Phys 2017; 147:155102. [PMID: 29055303 DOI: 10.1063/1.5006948] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022] Open
Abstract
Simple one-dimensional DNA or RNA mesoscopic models are of interest for their computational efficiency while retaining the key elements of the molecular interactions. However, they only deal with perfectly formed DNA or RNA double helices and consider the intra-strand interactions to be the same on both strands. This makes it difficult to describe highly asymmetric structures such as bulges and loops and, for instance, prevents the application of mesoscopic models to determine RNA secondary structures. Here we derived the conditions for the Peyrard-Bishop mesoscopic model to overcome these limitations and applied it to the calculation of single bulges, the smallest and simplest of these asymmetric structures. We found that these theoretical conditions can indeed be applied to any situation where stacking asymmetry needs to be considered. The full set of parameters for group I RNA bulges was determined from experimental melting temperatures using an optimization procedure, and we also calculated average opening profiles for several RNA sequences. We found that guanosine bulges show the strongest perturbation on their neighboring base pairs, considerably reducing the on-site interactions of their neighboring base pairs.
Collapse
Affiliation(s)
- Erik de Oliveira Martins
- Departamento de Física, Universidade Federal de Minas Gerais, 31270-901 Belo Horizonte, MG, Brazil
| | - Gerald Weber
- Departamento de Física, Universidade Federal de Minas Gerais, 31270-901 Belo Horizonte, MG, Brazil
| |
Collapse
|
7
|
Crowther CV, Jones LE, Morelli JN, Mastrogiacomo EM, Porterfield C, Kent JL, Serra MJ. Influence of two bulge loops on the stability of RNA duplexes. RNA (NEW YORK, N.Y.) 2017; 23:217-228. [PMID: 27872162 PMCID: PMC5238796 DOI: 10.1261/rna.056168.116] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/29/2016] [Accepted: 11/13/2016] [Indexed: 05/24/2023]
Abstract
Fifty-three RNA duplexes containing two single nucleotide bulge loops were optically melted in 1 M NaCl in order to determine the thermodynamic parameters ΔH°, ΔS°, ΔG°37, and TM for each duplex. Because of the large number of possible combinations and lack of sequence effects observed previously, we limited our initial investigation to adenosine bulges, the most common naturally occurring bulge. For example, the following duplexes were investigated: 5'GGCAXYAGGC/3'CCG YX CCG, 5'GGCAXY GCC/3'CCG YXACGG, and 5'GGC XYAGCC/3'CCGAYX CGG. The identity of XY (where XY are Watson-Crick base pairs) and the total number of base pairs in the terminal and central stems were varied. As observed for duplexes with a single bulge loop, the effect of the two bulge loops on duplex stability is primarily influenced by non-nearest neighbor interactions. In particular, the stability of the stems influences the destabilization of the duplex by the inserted bulge loops. The model proposed to predict the influence of multiple bulge loops on duplex stability suggests that the destabilization of each bulge is related to the stability of the adjacent stems. A database of RNA secondary structures was examined to determine the naturally occurring abundance of duplexes containing multiple bulge loops. Of the 2000 examples found in the database, over 65% of the two bulge loops occur within 3 base pairs of each other. A database of RNA three-dimensional structures was examined to determine the structure of duplexes containing two single nucleotide bulge loops. The structures of the bulge loops are described.
Collapse
Affiliation(s)
- Claire V Crowther
- Department of Chemistry, Allegheny College, Meadville, Pennsylvania 16335, USA
| | - Laura E Jones
- Department of Chemistry, Allegheny College, Meadville, Pennsylvania 16335, USA
| | - Jessica N Morelli
- Department of Chemistry, Allegheny College, Meadville, Pennsylvania 16335, USA
| | | | - Claire Porterfield
- Department of Chemistry, Allegheny College, Meadville, Pennsylvania 16335, USA
| | - Jessica L Kent
- Department of Chemistry, Allegheny College, Meadville, Pennsylvania 16335, USA
| | - Martin J Serra
- Department of Chemistry, Allegheny College, Meadville, Pennsylvania 16335, USA
| |
Collapse
|
8
|
Truitt AR, Choi BE, Li J, Soto AM. Effect of Mutations on the Binding of Kanamycin-B to RNA Hairpins Derived from the Mycobacterium tuberculosis Ribosomal A-Site. Biochemistry 2015; 54:7425-37. [PMID: 26560864 DOI: 10.1021/acs.biochem.5b00710] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Kanamycin is an aminoglycoside antibiotic used in the treatment of drug-resistant tuberculosis. Mutations at the rRNA A-site have been associated with kanamycin resistance in Mycobacterium tuberculosis clinical isolates. Understanding the effect of these mutations on the conformation of the M. tuberculosis A-site is critical for understanding the mechanisms of antibiotic resistance in M. tuberculosis. In this work, we have studied RNA hairpins derived from the M. tuberculosis A-site, the wild type and three mutants at the following positions (M. tuberculosis/Escherichia coli numbering): A1400/1408 → G, C1401/1409 → U, and the double mutant G1483/1491 C1401/1409 → UA. Specifically, we used circular dichroism, ultraviolet spectroscopy, and fluorescence spectroscopy to characterize the conformation, stability, and binding affinity of kanamycin-B and other aminoglycoside antibiotics for these RNA hairpins. Our results show that the mutations affect the conformation of the decoding site, with the mutations at position 1401/1409 resulting in significant destabilizations. Interestingly, the mutants bind paromomycin with weaker affinity than the wild type, but they bind kanamycin-B with similar affinity than the wild type. The results suggest that the presence of mutations does not prevent kanamycin-B from binding. Instead, kanamycin may promote different interactions with a third partner in the mutants compared to the wild type. Furthermore, our results with longer and shorter hairpins suggest that the region of the A-site that varies among organisms may have modulating effects on the binding and interactions of the A-site.
Collapse
Affiliation(s)
- Amber R Truitt
- Department of Chemistry and ‡Molecular Biology, Biochemistry and Bioinformatics Program, Towson University , 8000 York Road, Towson, Maryland 21252, United States
| | - Bok-Eum Choi
- Department of Chemistry and ‡Molecular Biology, Biochemistry and Bioinformatics Program, Towson University , 8000 York Road, Towson, Maryland 21252, United States
| | - Jenny Li
- Department of Chemistry and ‡Molecular Biology, Biochemistry and Bioinformatics Program, Towson University , 8000 York Road, Towson, Maryland 21252, United States
| | - Ana Maria Soto
- Department of Chemistry and ‡Molecular Biology, Biochemistry and Bioinformatics Program, Towson University , 8000 York Road, Towson, Maryland 21252, United States
| |
Collapse
|
9
|
Dishler AL, McMichael EL, Serra MJ. Determination of the secondary structure of group II bulge loops using the fluorescent probe 2-aminopurine. RNA (NEW YORK, N.Y.) 2015; 21:975-984. [PMID: 25805856 PMCID: PMC4408803 DOI: 10.1261/rna.048306.114] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 09/29/2014] [Accepted: 01/17/2015] [Indexed: 06/04/2023]
Abstract
Eleven RNA hairpins containing 2-aminopurine (2-AP) in either base-paired or single nucleotide bulge loop positions were optically melted in 1 M NaCl; and, the thermodynamic parameters ΔH°, ΔS°, ΔG°37, and TM for each hairpin were determined. Substitution of 2-AP for an A (adenosine) at a bulge position (where either the 2-AP or A is the bulge) in the stem of a hairpin, does not affect the stability of the hairpin. For group II bulge loops such as AA/U, where there is ambiguity as to which of the A residues is paired with the U, hairpins with 2-AP substituted for either the 5' or 3' position in the hairpin stem have similar stability. Fluorescent melts were performed to monitor the environment of the 2-AP. When the 2-AP was located distal to the hairpin loop on either the 5' or 3' side of the hairpin stem, the change in fluorescent intensity upon heating was indicative of an unpaired nucleotide. A database of phylogenetically determined RNA secondary structures was examined to explore the presence of naturally occurring bulge loops embedded within a hairpin stem. The distribution of bulge loops is discussed and related to the stability of hairpin structures.
Collapse
Affiliation(s)
- Abigael L Dishler
- Department of Chemistry, Allegheny College, Meadville, Pennsylvania 16335, USA
| | | | - Martin J Serra
- Department of Chemistry, Allegheny College, Meadville, Pennsylvania 16335, USA
| |
Collapse
|
10
|
Mouzakis KD, Dethoff EA, Tonelli M, Al-Hashimi H, Butcher SE. Dynamic motions of the HIV-1 frameshift site RNA. Biophys J 2015; 108:644-54. [PMID: 25650931 PMCID: PMC4317556 DOI: 10.1016/j.bpj.2014.12.006] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/07/2014] [Revised: 11/11/2014] [Accepted: 12/05/2014] [Indexed: 12/13/2022] Open
Abstract
The HIV-1 frameshift site (FS) plays a critical role in viral replication. During translation, the HIV-1 FS transitions from a 3-helix to a 2-helix junction RNA secondary structure. The 2-helix junction structure contains a GGA bulge, and purine-rich bulges are common motifs in RNA secondary structure. Here, we investigate the dynamics of the HIV-1 FS 2-helix junction RNA. Interhelical motions were studied under different ionic conditions using NMR order tensor analysis of residual dipolar couplings. In 150 mM potassium, the RNA adopts a 43°(±4°) interhelical bend angle (β) and displays large amplitude, anisotropic interhelical motions characterized by a 0.52(±0.04) internal generalized degree of order (GDOint) and distinct order tensor asymmetries for its two helices (η = 0.26(±0.04) and 0.5(±0.1)). These motions are effectively quenched by addition of 2 mM magnesium (GDOint = 0.87(±0.06)), which promotes a near-coaxial conformation (β = 15°(±6°)) of the two helices. Base stacking in the bulge was investigated using the fluorescent purine analog 2-aminopurine. These results indicate that magnesium stabilizes extrahelical conformations of the bulge nucleotides, thereby promoting coaxial stacking of helices. These results are highly similar to previous studies of the HIV transactivation response RNA, despite a complete lack of sequence similarity between the two RNAs. Thus, the conformational space of these RNAs is largely determined by the topology of their interhelical junctions.
Collapse
Affiliation(s)
- Kathryn D Mouzakis
- Department of Biochemistry, University of Wisconsin-Madison, Madison, Wisconsin
| | - Elizabeth A Dethoff
- Department of Chemistry, University of North Carolina, Chapel Hill, North Carolina
| | - Marco Tonelli
- Department of Biochemistry, University of Wisconsin-Madison, Madison, Wisconsin
| | | | - Samuel E Butcher
- Department of Biochemistry, University of Wisconsin-Madison, Madison, Wisconsin.
| |
Collapse
|
11
|
Kent JL, McCann MD, Phillips D, Panaro BL, Lim GF, Serra MJ. Non-nearest-neighbor dependence of stability for group III RNA single nucleotide bulge loops. RNA (NEW YORK, N.Y.) 2014; 20:825-34. [PMID: 24742935 PMCID: PMC4024637 DOI: 10.1261/rna.043232.113] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/30/2013] [Accepted: 02/21/2014] [Indexed: 05/24/2023]
Abstract
Thirty-five RNA duplexes containing single nucleotide bulge loops were optically melted and the thermodynamic parameters for each duplex determined. The bulge loops were of the group III variety, where the bulged nucleotide is either a AG/U or CU/G, leading to ambiguity to the exact position and identity of the bulge. All possible group III bulge loops with Watson-Crick nearest-neighbors were examined. The data were used to develop a model to predict the free energy of an RNA duplex containing a group III single nucleotide bulge loop. The destabilization of the duplex by the group III bulge could be modeled so that the bulge nucleotide leads to the formation of the Watson-Crick base pair rather than the wobble base pair. The destabilization of an RNA duplex caused by the insertion of a group III bulge is primarily dependent upon non-nearest-neighbor interactions and was shown to be dependent upon the stability of second least stable stem of the duplex. In-line structure probing of group III bulge loops embedded in a hairpin indicated that the bulged nucleotide is the one positioned further from the hairpin loop irrespective of whether the resulting stem formed a Watson-Crick or wobble base pair. Fourteen RNA hairpins containing group III bulge loops, either 3' or 5' of the hairpin loop, were optically melted and the thermodynamic parameters determined. The model developed to predict the influence of group III bulge loops on the stability of duplex formation was extended to predict the influence of bulge loops on hairpin stability.
Collapse
|
12
|
Murray MH, Hard JA, Znosko BM. Improved model to predict the free energy contribution of trinucleotide bulges to RNA duplex stability. Biochemistry 2014; 53:3502-8. [PMID: 24853497 PMCID: PMC4051427 DOI: 10.1021/bi500204e] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
![]()
Trinucleotide
bulges in RNA commonly occur in nature. Yet, little
data exists concerning the thermodynamic parameters of this motif.
Algorithms that predict RNA secondary structure from sequence currently
attribute a constant free energy value of 3.2 kcal/mol to all trinucleotide
bulges, regardless of bulge sequence. To test the accuracy of this
model, RNA duplexes that contain frequent naturally occurring trinucleotide
bulges were optically melted, and their thermodynamic parameters—enthalpy,
entropy, free energy, and melting temperature—were determined.
The thermodynamic data were used to derive a new model to predict
the free energy contribution of trinucleotide bulges to RNA duplex
stability: ΔG°37, trint bulge = ΔG°37, bulge + ΔG°37, AU + ΔG°37, GU. The parameter ΔG°37, bulge is variable depending upon the purine
and pyrimidine composition of the bulge, ΔG°37, AU is a 0.49 kcal/mol penalty for an A-U
closing pair, and ΔG° 37, GU is a −0.56 kcal/mol bonus for a G-U closing pair. With both
closing pair and bulge sequence taken into account, this new model
predicts free energy values within 0.30 kcal/mol of the experimental
value. The new model can be used by algorithms that predict RNA free
energies as well as algorithms that use free energy minimization to
predict RNA secondary structure from sequence.
Collapse
Affiliation(s)
- Meghan H Murray
- Department of Chemistry, Saint Louis University , Saint Louis, Missouri 63103, United States
| | | | | |
Collapse
|
13
|
Lim GF, Merz GE, McCann MD, Gruskiewicz JM, Serra MJ. Stability of single-nucleotide bulge loops embedded in a GAAA RNA hairpin stem. RNA (NEW YORK, N.Y.) 2012; 18:807-14. [PMID: 22345128 PMCID: PMC3312567 DOI: 10.1261/rna.028308.111] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/23/2011] [Accepted: 12/11/2011] [Indexed: 05/24/2023]
Abstract
Forty-six RNA hairpins containing combinations of 3' or 5' bulge loops and a 3' or 5' fluorescein label were optically melted in 1 M NaCl, and the thermodynamic parameters ΔH°, ΔS°, ΔG°(37), and T(M) for each hairpin were determined. The bulge loops were of the group I variety, in which the identity of the bulge is known, and the group II variety, in which the bulged nucleotide is identical to one of its nearest neighbors, leading to ambiguity as to the exact position of the bulge. The fluorescein label at either the 3' end or 5' end of the hairpin did not significantly influence the stability of the hairpin. As observed with bulge loops inserted into a duplex motif, the insertion of a bulge loop into the stem of a hairpin loop was destabilizing. The model developed to predict the influence of bulge loops on the stability of duplex formation was extended to predict the influence of bulge loops on hairpin stability. Specifically, the influence of the bulge is related to the stability of the hairpin stem distal from the hairpin loop.
Collapse
Affiliation(s)
- Geoffrey F.S. Lim
- Department of Chemistry, Allegheny College, Meadville, Pennsylvania 16335, USA
| | - Gregory E. Merz
- Department of Chemistry, Allegheny College, Meadville, Pennsylvania 16335, USA
| | - Michael D. McCann
- Department of Chemistry, Allegheny College, Meadville, Pennsylvania 16335, USA
| | | | - Martin J. Serra
- Department of Chemistry, Allegheny College, Meadville, Pennsylvania 16335, USA
| |
Collapse
|
14
|
Vanegas PL, Horwitz TS, Znosko BM. Effects of non-nearest neighbors on the thermodynamic stability of RNA GNRA hairpin tetraloops. Biochemistry 2012; 51:2192-8. [PMID: 22329761 DOI: 10.1021/bi300008j] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
Abstract
Currently, several models for predicting the secondary structure of RNA exist, one of which is free energy minimization using the Nearest Neighbor Model. This model predicts the lowest-free energy secondary structure from a primary sequence by summing the free energy contributions of the Watson-Crick nearest neighbor base pair combinations and any noncanonical secondary structure motif. The Nearest Neighbor Model also assumes that the free energy of the secondary structure motif is dependent solely on the identities of the nucleotides within the motif and the motif's nearest neighbors. To test the current assumption of the Nearest Neighbor Model that the non-nearest neighbors do not affect the stability of the motif, we optically melted different stem-loop oligonucleotides to experimentally determine their thermodynamic parameters. In each of these oligonucleotides, the hairpin loop sequence and the adjacent base pairs were held constant, while the first or second non-nearest neighbors were varied. The experimental results show that the thermodynamic contributions of the hairpin loop were dependent upon the identity of the first non-nearest neighbor, while the second non-nearest neighbor had a less obvious effect. These results were then used to create an updated model for predicting the thermodynamic contributions of a hairpin loop to the overall stability of the stem-loop structure.
Collapse
Affiliation(s)
- Pamela L Vanegas
- Department of Chemistry, Saint Louis University, Saint Louis, Missouri 63103, United States
| | | | | |
Collapse
|