With both catalytic and genetic functions, ribonucleic acid (RNA) is perhaps the most pluripotent chemical species in molecular biology, and its functions are intimately linked to its structure and dynamics. Computer simulations, and in particular atomistic molecular dynamics (MD), allow structural dynamics of biomolecular systems to be investigated with unprecedented temporal and spatial resolution. We here provide a comprehensive overview of the fast-developing field of MD simulations of RNA molecules. We begin with an in-depth, evaluatory coverage of the most fundamental methodological challenges that set the basis for the future development of the field, in particular, the current developments and inherent physical limitations of the atomistic force fields and the recent advances in a broad spectrum of enhanced sampling methods. We also survey the closely related field of coarse-grained modeling of RNA systems. After dealing with the methodological aspects, we provide an exhaustive overview of the available RNA simulation literature, ranging from studies of the smallest RNA oligonucleotides to investigations of the entire ribosome. Our review encompasses tetranucleotides, tetraloops, a number of small RNA motifs, A-helix RNA, kissing-loop complexes, the TAR RNA element, the decoding center and other important regions of the ribosome, as well as assorted others systems. Extended sections are devoted to RNA–ion interactions, ribozymes, riboswitches, and protein/RNA complexes. Our overview is written for as broad of an audience as possible, aiming to provide a much-needed interdisciplinary bridge between computation and experiment, together with a perspective on the future of the field.

Recent computational efforts have shown that the current potential energy models used in molecular dynamics are not accurate enough to describe the conformational ensemble of RNA oligomers and suggest that molecular dynamics should be complemented with experimental data. We here propose a scheme based on the maximum entropy principle to combine simulations with bulk experiments. In the proposed scheme the noise arising from both the measurements and the forward models used to back calculate the experimental observables is explicitly taken into account. The method is tested on RNA nucleosides and is then used to construct chemically consistent corrections to the Amber RNA force field that allow a large set of experimental data on nucleosides and dinucleosides to be correctly reproduced. The transferability of these corrections is assessed against independent data on tetranucleotides and displays a previously unreported agreement with experiments. This procedure can be applied to enforce multiple experimental data on multiple systems in a self-consistent framework thus suggesting a new paradigm for force field refinement.

The computational study of conformational transitions in RNA and proteins with atomistic molecular dynamics often requires suitable enhanced sampling techniques. We here introduce a novel method where concurrent metadynamics are integrated in a Hamiltonian replica-exchange scheme. The ladder of replicas is built with different strengths of the bias potential exploiting the tunability of well-tempered metadynamics. Using this method, free-energy barriers of individual collective variables are significantly reduced compared with simple force-field scaling. The introduced methodology is flexible and allows adaptive bias potentials to be self-consistently constructed for a large number of simple collective variables, such as distances and dihedral angles. The method is tested on alanine dipeptide and applied to the difficult problem of conformational sampling in a tetranucleotide.

The computational study of conformational transitions in nucleic acids still faces many challenges. For example, in the case of single stranded RNA tetranucleotides, agreement between simulations and experiments is not satisfactory due to inaccuracies in the force fields commonly used in molecular dynamics simulations. We here use experimental data collected from high-resolution X-ray structures to attempt an improvement of the latest version of the AMBER force field. A modified metadynamics algorithm is used to calculate correcting potentials designed to enforce experimental distributions of backbone torsion angles. Replica-exchange simulations of tetranucleotides including these correcting potentials show significantly better agreement with independent solution experiments for the oligonucleotides containing pyrimidine bases. Although the proposed corrections do not seem to be portable to generic RNA systems, the simulations revealed the importance of the α and ζ backbone angles for the modulation of the RNA conformational ensemble. The correction protocol presented here suggests a systematic procedure for force-field refinement.

We introduce a method for predicting RNA folding pathways, with an application to the most important RNA tetraloops. The method is based on the idea that ensembles of three-dimensional fragments extracted from high-resolution crystal structures are heterogeneous enough to describe metastable as well as intermediate states. These ensembles are first validated by performing a quantitative comparison against available solution nuclear magnetic resonance (NMR) data of a set of RNA tetranucleotides. Notably, the agreement is better with respect to the one obtained by comparing NMR with extensive all-atom molecular dynamics simulations. We then propose a procedure based on diffusion maps and Markov models that makes it possible to obtain reaction pathways and their relative probabilities from fragment ensembles. This approach is applied to study the helix-to-loop folding pathway of all the tetraloops from the GNRA and UNCG families. The results give detailed insights into the folding mechanism that are compatible with available experimental data and clarify the role of intermediate states observed in previous simulation studies. The method is computationally inexpensive and can be used to study arbitrary conformational transitions.

We have carried out a series of extended unbiased molecular dynamics (MD) simulations (up to 10 μs long, ∼162 μs in total) complemented by replica-exchange with the collective variable tempering (RECT) approach for several human telomeric DNA G-quadruplex (GQ) topologies with TTA propeller loops. We used different AMBER DNA force-field variants and also processed simulations by Markov State Model (MSM) analysis. The slow conformational transitions in the propeller loops took place on a scale of a few μs, emphasizing the need for long simulations in studies of GQ dynamics. The propeller loops sampled similar ensembles for all GQ topologies and for all force-field dihedral-potential variants. The outcomes of standard and RECT simulations were consistent and captured similar spectrum of loop conformations. However, the most common crystallographic loop conformation was very unstable with all force-field versions. Although the loss of canonical γ-trans state of the first propeller loop nucleotide could be related to the indispensable bsc0 α/γ dihedral potential, even supporting this particular dihedral by a bias was insufficient to populate the experimentally dominant loop conformation. In conclusion, while our simulations were capable of providing a reasonable albeit not converged sampling of the TTA propeller loop conformational space, the force-field description still remained far from satisfactory.

Over two-thirds of integral membrane proteins of known structure assemble into oligomers. Yet, the forces that drive the association of these proteins remain to be delineated, as the lipid bilayer is a solvent environment that is both structurally and chemically complex. In this study we reveal how the lipid solvent defines the dimerization equilibrium of the CLC-ec1 Cl-/H+ antiporter. Integrating experimental and computational approaches, we show that monomers associate to avoid a thinned-membrane defect formed by hydrophobic mismatch at their exposed dimerization interfaces. In this defect, lipids are strongly tilted and less densely packed than in the bulk, with a larger degree of entanglement between opposing leaflets and greater water penetration into the bilayer interior. Dimerization restores the membrane to a near-native state and therefore, appears to be driven by the larger free-energy cost of lipid solvation of the dissociated protomers. Supporting this theory, we demonstrate that addition of short-chain lipids strongly shifts the dimerization equilibrium towards the monomeric state, and show that the cause of this effect is that these lipids preferentially solvate the defect. Importantly, we show that this shift requires only minimal quantities of short-chain lipids, with no measurable impact on either the macroscopic physical state of the membrane or the protein's biological function. Based on these observations, we posit that free-energy differentials for local lipid solvation define membrane-protein association equilibria. With this, we argue that preferential lipid solvation is a plausible cellular mechanism for lipid regulation of oligomerization processes, as it can occur at low concentrations and does not require global changes in membrane properties.

Adsorption of enzymes on solid surfaces may lead to conformational changes that reduce their catalytic conversion activity and are thus detrimental to the efficiency of biotechnology or biosensing applications. This work is a joint theoretical and experimental endeavor in which we identify and quantify the conformational changes that 1 Page 2 of 46 ACS Paragon Plus Environment ACS Biomaterials Science & Engineering chymotrypsin undergoes when in contact with the surface of amorphous silica nanoparticles. For this purpose, we use circular dichroism spectroscopy, standard molecular dynamics and advanced-sampling methods. Only the combination of these techniques allowed us to pinpoint a destabilization effect of silica on specific structural motifs of chymotrypsin. They are linked by the possibility of theoretically predicting CD spectra, allowing us to elucidate the source of the experimentally observed spectral changes.We find that chymotrypsin loses part of its helical content upon adsorption, with minor perturbation of its overall tertiary structure, associated to changes in the aromatic interactions. We demonstrate that the C-terminal helical fragment is unfolded as an isolated oligopeptide in pure water, folded as an α-helix as terminus of chymotrypsin in solution, and again partly disordered when the protein is adsorbed on silica. We believe that the joint methodology introduced in this manuscript has a direct general applicability to investigate any biomolecule -inorganic surface system. Methods to theoretically predict Circular Dichroism spectra from atomistic simulations were compared and improved. The drawbacks of the approaches are discussed; in particular the limited capability of advanced-sampling MD schemes to explore the conformational phase space of large proteins, and the dependency of the predicted ellipticity bands on the choice of calculation parameters.