Effect of ligand binding on riboswitch folding: Theory and simulations

The effect of ligand binding on the conformational transitions of the add A-riboswitch in cellular environments is investigated theoretically within the framework of the generalized Langevin equation combined with steered molecular dynamics simulations. Results for the transition path time distribution provide an estimate of the transit times, which are difficult to determine experimentally. The time for the conformational transitions of the riboswitch aptamer is longer for the ligand bound state as compared to that of the unbound one. The transition path time of the riboswitch follows a counterintuitive trend as it decreases with an increase in the barrier height. The mean transition path time of either transitions of the riboswitch in the ligand bound/unbound state increases with an increase in the complexity of the surrounding environment due to the caging effect. The results of the probability density function, transition path time distribution, and mean transition path time obtained from the theory qualitatively agree with those obtained from the simulations and with earlier experimental and theoretical studies.