RNA folding landscapes from explicit solvent all-atom simulations

Atomically detailed simulations of RNA folding have proven very challenging in view of the difficulties of developing realistic force fields and the intrinsic computational complexity of sampling rare conformational transitions. To tackle both these issues, we extend to RNA an enhanced path sampling method previously successfully applied to proteins. In this scheme, the information about the RNA's native structure is harnessed by a soft history-dependent biasing force, which is added to the atomistic force field, thus promoting the generation of productive folding trajectories. Here, we report on the results of simulations in explicit solvent of RNA molecules from 20 to 47 nucleotides long and increasing topological complexity. From a statistical analysis of the folding pathways we infer that, differently from proteins, the underlying free energy landscape is significantly frustrated, even for relatively small chains with a simple topology. The folding mechanism and the thermodynamics are in agreement with the available experiments and some of the existing coarse-grained models. This scheme provides a fully microscopic characterization of RNA folding, relating the kinetics and dynamics of the transition to the chemistry of the chain and its solvent. Therefore, it provides a transferable framework that sets the stage for future translational applications