Interfacing medicinal chemistry with structural bioinformatics: implications for T box riboswitch RNA drug discovery
BMC Bioinformatics volume 13, Article number: S5 (2012)
The T box riboswitch controls bacterial transcription by structurally responding to tRNA aminoacylation charging ratios. Knowledge of the thermodynamic stability difference between two competing structural elements within the riboswitch, the terminator and the antiterminator, is critical for effective T box-targeted drug discovery.
The ΔG of aminoacyl tRNA synthetase (aaRS) T box riboswitch terminators and antiterminators was predicted using DINAMelt and the resulting ΔΔG (ΔGTerminator - ΔGAntiterminator) values were compared.
Average ΔΔG values did not differ significantly between the bacterial species analyzed, but there were significant differences based on the type of aaRS.
The data indicate that, of the bacteria studied, there is little potential for drug targeting based on overall bacteria-specific thermodynamic differences of the T box antiterminator vs. terminator stability, but that aaRS-specific thermodynamic differences could possibly be exploited for designing drug specificity.
The T box riboswitch is an important regulatory mechanism found in Gram-positive bacteria including many pathogens [1–3]. The riboswitch responds to high levels of uncharged (non-aminoacylated) tRNA to control the transcription of cognate genes (e.g., aminoacyl tRNA synthetase, aaRS, genes) . Cognate, uncharged tRNA binds the 5'-untranslated region of T box mRNA during transcription and, when present in sufficient quantities, results in antitermination (Figure 1a) . The tRNA anticodon binds a specifier sequence in Stem 1, thus providing the cognate specificity, while the uncharged tRNA acceptor end nucleotides bind the first four bases in a seven nucleotide bulge of the highly conserved antiterminator structural element . Aminoacylated-tRNA also binds the specifier sequence, but does not bind the antiterminator element . In the absence of uncharged tRNA bound to the antiterminator element, a more thermodynamically stable stem-loop structure forms (the terminator) followed by factor-independent transcription termination a few nucleotides later . The antiterminator and the terminator are mutually exclusive structural elements due to sharing common nucleotides (Figure 1b). Ligands that target and disrupt the T box riboswitch function could be potential antibacterial agents in light of the critical genes regulated by the T box mechanism [1, 2].
We have been investigating the structure-function relationship of the T box antiterminator element and the key recognition features necessary for ligands to specifically bind the antiterminator and disrupt its function. There are few detailed medicinal chemistry studies of ligands targeting RNA  and this project has required an extensive multidisciplinary approach. The solution structure of antiterminator model RNA AM1A  was determined using molecular modelling with NMR-derived constraints . The structure indicated that the seven-nucleotide bulge of the antiterminator was not pre-ordered for tRNA binding, but rather, that binding of the tRNA acceptor end must require a certain extent of tertiary-structure capture and/or an induced fit. Fluorescence life-time studies confirmed that modest antiterminator structural reorganization occurs upon tRNA binding in a magnesium-dependent manner . While some RNAs have specific divalent metal ion binding sites, for the antiterminator RNA, the Mg2+ binds via a diffuse, outer-sphere interaction . In vitro selection studies of both the antiterminator  and tRNA  indicate that there are likely no direct interactions between the tRNA and the antiterminator other than the known base pairing between the acceptor end nucleotides and the first four nucleotides at the 5'-end of the seven-nucleotide bulge. Given this antiterminator structure-function information, ligands could potentially disrupt tRNA binding simply by competing with the base pairing between the tRNA acceptor end and the antiterminator bulge nucleotides.
Aminoglycosides bind AM1A in a structure-specific manner, most likely via electrostatics [13, 14]. In contrast, two novel classes of heterocyclic compounds, triazoles and oxazolidinones, bind AM1A with structure-specificity and high affinity, but without reliance on electrostatics [15–18] and can alter T box transcription antitermination . A computational, quantitative structure activity relationship analysis has shown that hydrophobic interactions play a significant role in the binding of these compounds to AM1A . Furthermore, the small molecule ligands disrupt tRNA binding to the antiterminator in a structure-specific manner .
From a drug discovery perspective, a key factor to determine is the range of ligand-induced stabilization that can be accommodated without overly stabilizing the antiterminator element and precluding terminator formation. The goal of our T box drug discovery project is to determine the key ligand features that lead to specific antiterminator binding, but that do not result in excessive stabilization of the antiterminator secondary structure. These ligands could then potentially compete with tRNA for binding to the antiterminator, but still allow terminator formation such that transcription of a T box gene critical for bacterial survival would be terminated (Figure 2). The predicted thermodynamic stability (ΔG) of the terminator and antiterminator structural elements have been compared for the B. subtilis tyrS T box . However, there has been no systematic comparison of predicted thermodynamic stability differences for a larger set of T box genes. Using a structural bioinformatics approach, we have analyzed the differences in predicted free energy (ΔΔG) between antiterminators and terminators in a set of aaRS T box genes in order to predict an upper limit of ligand-induced stabilization that can potentially be accommodated.
The thermodynamic stability of T box antiterminator and terminator structural elements was calculated using the DINAMelt server . The DINAMelt server computes the secondary structure and free energy of the folded RNA using a secondary structure folding algorithm and empirically-derived nearest neighbour coefficients . The folding algorithm predicts the minimum energy RNA secondary structure using the available thermodynamic data for base pairing, base stacking and destabilizing energies [22, 23]. The sequences analyzed were predicted to be involved in aaRS T box antiterminator and terminator structural elements from Bacillus cereus (BC), Bacillus subtilis (BS), Clostridium botulinum (CB), Clostridium difficile (CDF), Clostridium perfringens (CPE), Staphylococcus aureus Mu50 (SA), Streptococcus agalactiae (SAG), and, Streptococcus pyogenes (SPY) . The ΔG for each aaRS T box antiterminator and terminator sequence was determined using the Quickfold (RNA 3.0) option on the DINAMelt server and the predicted thermodynamic stability difference calculated from ΔΔG = ΔGTerminator - ΔGAntiterminator. The % suboptimal setting was adjusted as necessary to obtain the lowest free energy antiterminator fold that had the consensus secondary structure of the core seven-nucleotide bulge containing the 5' -UGGN-3' nucleotides that are complementary to the tRNA acceptor end nucleotides .
Results and discussion
The free energy values for the antiterminator and terminator structural elements of the T box genes analyzed were calculated using the DINAMelt webserver . The predicted ΔG values are listed in Additional File 1 and the relative comparisons of terminator vs.antiterminator stability (ΔΔG) are summarized in Figure 3. The overall average ΔΔG for all aaRS studied was -12.8 kcal/mol. The average ΔΔG values did not differ significantly between bacteria when comparing between the pathogenic bacteria nor between pathogenic vs. the non-pathogenic bacteria studied (BS) (Figure 3a). In contrast, there were significant differences in ΔΔG averages between specific aaRS with alanyl aaRS having the smallest average ΔG (-7.8 ± 3.5 kcal/mol) and glycyl aaRS having the largest (-20.1 ± 4.6 kcal/mol) (Figure 3b). Based on these results, the glycyl aaRS, on average, may be best able to accommodate ligand-induced stabilization of the antiterminator and still allow formation of the terminator. An important factor to consider, however, is that the free energy calculations are based on empirically-derived parameters for known RNA structural motifs . Structural motifs, especially in loops and bulges, that have not been previously characterized might contribute to the stability of the RNA elements and not be accounted for in the DINAMelt ΔG calculations. Since the loop of the antiterminator is not highly conserved , most likely there is no structural motif within it that might contribute to the antiterminator stability, however, the possibility cannot be excluded. Further investigation of experimentally-derived free energy values of individual antiterminators and terminators is needed.
The free energy of T box riboswitch antiterminator and terminator elements was predicted and compared for a series of aaRS T box genes. The observed aaRS-specific stability differences between these key riboswitch structural elements could potentially be targeted to effect ligand-specificity in future drug discovery efforts.
Aminoacyl tRNA synthetase
- Mu50 (SA):
Streptococcus agalactiae and, Streptococcus pyogenes (SPY).
Gutierrez-Preciado A, Henkin TM, Grundy FJ, Yanofsky C, Merino E: Biochemical features and functional implications of the RNA-based T box regulatory mechanism. Microbiol Mol Biol Rev. 2009, 73: 36-61. 10.1128/MMBR.00026-08.
Vitreschak AG, Mironov AA, Lyubetsky VA, Gelfand MS: Comparative genomic analysis of T-box regulatory systems in bacteria. RNA. 2008, 14: 717-735. 10.1261/rna.819308.
Green NJ, Grundy FJ, Henkin TM: The T box mechanism: tRNA as a regulatory molecule. FEBS Lett. 2010, 584 (2): 318-324. 10.1016/j.febslet.2009.11.056.
Grundy FJ, Moir TR, Haldeman MT, Henkin TM: Sequence requirements for terminators and antiterminators in the T box transcription antitermination system: disparity between conservation and functional requirements. Nucl Acids Res. 2002, 30 (7): 1646-1655. 10.1093/nar/30.7.1646.
Grundy FJ, Yousef MR, Henkin TM: Monitoring uncharged tRNA during transcription of the Bacillus subtilis glyQS gene. J Mol Biol. 2005, 346: 73-81. 10.1016/j.jmb.2004.11.051.
Thomas JR, Hergenrother PJ: Targeting RNA with Small Molecules. Chem Rev. 2008, 108 (4): 1171-1224. 10.1021/cr0681546.
Gerdeman MS, Henkin TM, Hines JV: In vitro structure-function studies of the Bacillus subtilis tyrS mRNA antiterminator: Evidence for factor independent tRNA acceptor stem binding specificity. Nucleic Acids Res. 2002, 30 (4): 1065-1072. 10.1093/nar/30.4.1065.
Gerdeman MS, Henkin TM, Hines JV: Solution structure of the B. subtilis T box antiterminator RNA: Seven-nucleotide bulge characterized by stacking and flexibility. J Mol Biol. 2003, 326: 189-201. 10.1016/S0022-2836(02)01339-6.
Means JA, Simson CM, Zhou S, Rachford AA, Rack J, Hines JV: Fluorescence probing of T box antiterminator RNA: Insights into riboswitch discernment of the tRNA discriminator base. Biochem Biophys Res Commun. 2009, 389: 616-621. 10.1016/j.bbrc.2009.09.037.
Jack KD, Means JA, Hines JV: Characterizing riboswitch function: Identification of Mg2+ binding site in T box antiterminator RNA. Biochem Biophys Res Commun. 2008, 370: 306-310. 10.1016/j.bbrc.2008.03.079.
Fauzi H, Agyeman A, Hines JV: T box transcription antitermination riboswitch: Influence of nucleotide sequence and orientation on tRNA binding by the antiterminator element. Biochim Biophys Acta. 2009, 1789: 185-191.
Fauzi H, Jack KD, Hines JV: In vitro selection to identify determinants in tRNA for Bacillus subtilis tyrS T box antiterminator mRNA binding. Nucleic Acids Res. 2005, 33: 2595-2602. 10.1093/nar/gki546.
Means JA, Hines JV: Fluorescence resonance energy transfer studies of aminoglycoside binding to a T box antiterminator RNA. Bioorg Med Chem Lett. 2005, 15: 2169-2172. 10.1016/j.bmcl.2005.02.007.
Anupam R, Denapoli L, Muchenditsi AM, Hines JV: Identification of neomycin B binding site in T box antiterminator model RNA. Bioorg Med Chem. 2008, 16: 4466-4470. 10.1016/j.bmc.2008.02.056.
Means JA, Katz SJ, Nayek A, Anupam R, Hines JV, Bergmeier SC: Structure activity studies of oxazolidinone analogs as RNA-binding agents. Bioorg Med Chem Lett. 2006, 16 (13): 3600-3604. 10.1016/j.bmcl.2006.03.068.
Anupam R, Bergmeier SC, Green NJ, Grundy FJ, Henkin TM, Means JA, Nayek A, Hines JV: 4,5-Disubstituted Oxazolidinones: High affinity molecular effectors of RNA function. Bioorg Med Chem Lett. 2008, 18: 3541-3544. 10.1016/j.bmcl.2008.05.015.
Acquaah-Harrison G, Zhou S, Hines JV, Bergmeier SC: Library of 1,4-disubstituted 1,2,3-triazole analogs of oxazolidinone RNA-binding agents. J Comb Chem. 2010, 12: 491-496. 10.1021/cc100029y.
Orac CM, Zhou S, Means JA, Boehme D, Bergmeier SC, Hines JV: Synthesis and stereospecificity of 4,5-disubstituted oxazolidinone ligands binding to T-box riboswitch RNA. J Med Chem. 2011, DOI: 10.1021/jm2006904,
Maciagiewicz I, Zhou S, Bergmeier SC, Hines JV: Structure activity studies of RNA-binding oxazolidinone derivatives. Bioorg Med Chem Lett. 2011, 21: 4524-4527. 10.1016/j.bmcl.2011.05.130.
Zhou S, Acquaah-Harrison G, Bergmeier SC, Hines JV: Anisotropy studies of tRNA - T box antiterminator RNA complex in the presence of 1,4-disubstituted 1,2,3-triazoles. Bioorg Med Chem Lett. 2011, DOI: 10.1016/j.bmcl.2011.1009.1095,
Markham NR, Zuker M: DINAMelt web server for nucleic acid melting prediction. Nucleic Acids Res. 2005, 33 (33 Web Server): W577-W581.
Zuker M, Stiegler P: Optimal computer folding of larger RNA sequences using thermodynamic and auxiliary information. Nucleic Acids Res. 1981, 9 (1): 133-148. 10.1093/nar/9.1.133.
Zuker M, Mathew DH, Turner DH: Algorithms and thermodynamics for RNA secondary structure prediction: a practical guide. RNA Biochemistry and Biotechnology. Edited by: Barciszewski J, Clark BFC. 1999, Kluwer Academic Publishers, 11-43.
We wish to thank the National Institutes of Health (GM073188) for support of this work and the office of the VP for Research, Ohio University and the DAAD for support of FJ on an Internationale Studien- und Ausbildungspartnerschaften (ISAP) grant.
This article has been published as part of BMC Bioinformatics Volume 13 Supplement 2, 2012: Proceedings from the Great Lakes Bioinformatics Conference 2011. The full contents of the supplement are available online at http://www.biomedcentral.com/bmcbioinformatics/supplements/13/S2
JH is a co-inventor on U.S. Patent number 7,005,441 for T box binding assays and ligands.
FJ carried out the sequence folding and free energy predictions. JH conceived of the study, designed and coordinated the study and prepared the manuscript.
Electronic supplementary material
About this article
Cite this article
Jentzsch, F., Hines, J.V. Interfacing medicinal chemistry with structural bioinformatics: implications for T box riboswitch RNA drug discovery. BMC Bioinformatics 13 (Suppl 2), S5 (2012). https://doi.org/10.1186/1471-2105-13-S2-S5