Insights into the dynamics of cyclic diguanosine monophosphate I riboswitch using molecular dynamics simulation

Abstract

Riboswitches are key cis-regulatory elements, present at 5' UTRs of several prokaryotic mRNAs, involved in gene expression regulation by binding selectively to specific ligands followed by conformational changes. However, understanding of ligand-free riboswitch, conformational changes between the ligand-free and ligand-bound riboswitch and binding mechanism of ligand to the aptamer of riboswitch is limited. In the present paper we describe the structural and dynamical properties of ligand-free c-di-GMP I riboswitch aptamer and its possible binding mechanism with c-di-GMP by means of all-atom molecular dynamics (MD) simulations. Various analyses such as principal component analysis, cross correlation dynamics analysis, network analysis and trajectory analyses were carried out. The ligand-free structure shows stable conformation with folded P2, P3 and an unwind P1 helix with open binding pocket at helix-join-helix junction while the ligand-bound structure showed closed binding pocket structure compared to the ligand-free structure. The junction residues significantly showed anti-correlations with P1 helix and weak correlated motions with each other in the open conformation of the ligand-free aptamer of riboswitch. The networks analysis of binding pocket residues suggested interaction among the identified key residues of the binding pocket (G20, A47, C92) and P1 helix illustrating the role of P1 helix in binding of ligand. The structural insights, on c-di-GMP I riboswitch, presented in this paper can be useful for the development of therapeutics against V. cholerae. The understanding of the c-di-GMP I riboswitch at dynamic molecular level provides a potential solution for riboswitch drug design.