Molecular dynamics simulation on the Thermosinus carboxydivorans pfl ZTP riboswitch by ligand binding

Riboswitches are RNA molecules that can regulate gene expression which is affected by ligand-binding during cotranscriptional folding process. However, the role of ligand during the folding is still unclear. In this study, the pfl domain of Thermosinus carboxydivorans ZTP riboswitch was discussed. The ligand is molecule ZMP. We mainly analyzed the change of ZMP-free and ZMP-bound aptamer domain by the dynamics simulation method. Structural features by calculating their RMSD, RMSF, etc. are analyzed. The results demonstrate that the binding domain require the presence of ZMP to maintain a stable fold. It also suggested that ZMP specificly binding to ZTP can generate more hydrogen bonds in the binding domain. Through the calculation of binding free energy decomposition of each nucleotide, molecule ZMP was found to promote the recognition and binding process of ligands by controlling some special nucleotides in the process of ligand binding. At last, the dynamical correlation and components of conformational motions were both applied to explore the effect of molecule ZMP to ZTP riboswitch. In general, ZMP can effectively affect the motions of the pfl riboswitch and facilitate the folding process of the ZTP riboswitch.These results may provide some new ideas for structural changes in riboswitches and their cotranscriptional folding process.