Pairing interactions between nucleobases and ligands in aptamer:ligand complexes of riboswitches: Crystal structure analysis, classification, optimal structures and accurate interaction energies

Abstract

ABSTRACTIn the present work, sixty-seven crystal structures of the aptamer domains of RNA riboswitches, are chosen for analysis of the structure and strength of hydrogen bonding (pairing) interactions between nucleobases constituting the aptamer binding pockets and the bound ligands. A total of eighty unique base:ligand hydrogen-bonded pairs containing at least two hydrogen bonds were identified through visual inspection. Classification of these contacts in terms of the interacting edge of the aptamer nucleobase revealed that interactions involving the Watson-Crick edge are the most common, followed by the sugar edge of purines and the Hoogsteen edge of uracil. Alternatively, classification in terms of the chemical constitution of the ligand yields five unique classes of base:ligand pairs: base:base, base:amino acid, base:sugar, base:phosphate and base:other. Further, quantum mechanical (QM) geometry optimizations revealed that sixty seven out of eighty pairs exhibit stable geometries and optimal deviations from their macromolecular crystal occurrences. This indicates that these contacts are well-defined RNA aptamer:ligand interaction motifs. QM calculated interaction energies of base:ligand pairs reveal rich hydrogen bonding landscape, ranging from weak interactions (base:other, –3 kcal/mol) to strong (base:phosphate, –48 kcal/mol) contacts. The analysis was further extended to study the biological importance of base:ligand interactions in the binding pocket of the tetrahydrofolate riboswitch and thiamine pyrophosphate riboswitch. Overall, our study helps in understanding the structural and energetic features of base:ligand pairs in riboswitches, which could aid in developing meaningful hypotheses in context of RNA:ligand recognition. This can, in turn contribute towards current efforts to develop antimicrobials that target RNAs.