TPP riboswitch aptamer: Role of Mg2+ ions, ligand unbinding, and allostery

Abstract

Riboswitches are non-coding RNAs that regulate gene expression in response to the binding of metabolites. Their abundance in bacteria makes them ideal drug targets. The prokaryotic thiamine pyrophosphate (TPP) riboswitch regulates gene expression in a wide range of bacteria by undergoing conformational changes in response to the binding of TPP. Although an experimental structure for the aptamer domain of the riboswitch is now available, details of the conformational changes that occur during the binding of the ligand, and the factors that govern these conformational changes, are still not clear. This study employs microsecond-scale molecular dynamics simulations to provide insights into the functioning of the riboswitch aptamer in atomistic detail. A mechanism for the transmission of conformational changes from the ligand-binding site to the P1 switch helix is proposed. Mg²⁺ ions in the binding site play a critical role in anchoring the ligand to the riboswitch. Finally, modeling the egress of TPP from the binding site reveals a two-step mechanism for TPP unbinding. Findings from this study can motivate the design of future studies aimed at modulating the activity of this drug target.

Keywords: Cation-RNA interactions; Ligand unbinding; Molecular dynamics; Ribonucleic acid (RNA); Thiamine pyrophosphate riboswitch.