Molecular dynamics (MD) simulations are an important and well-established tool for investigating RNA structural dynamics, but their accuracy relies heavily on the quality of the employed force field (ff). In this work, we present a comprehensive evaluation of widely used pair-additive and polarizable RNA ffs using the challenging UUCG tetraloop (TL) benchmark system. Extensive standard MD simulations, initiated from the NMR structure of the 14-mer UUCG TL, revealed that most ffs did not maintain the native state, instead favoring alternative loop conformations. Notably, three very recent variants of pair-additive ffs, OL3CP-gHBfix21, DES-Amber, and OL3R2.7, successfully preserved the native structure over a 10 × 20 μs time scale. To further assess these ffs, we performed enhanced sampling folding simulations of the shorter 8-mer UUCG TL, starting from the single-stranded conformation. Estimated folding free energies (ΔG°fold) varied significantly among these three ffs, with values of 0.0 ± 0.6, 2.4 ± 0.8, and 7.4 ± 0.2 kcal/mol for OL3CP-gHBfix21, DES-Amber, and OL3R2.7, respectively. The ΔG°fold value predicted by the OL3CP-gHBfix21 ff was closest to experimental estimates, ranging from -1.6 to -0.7 kcal/mol. In contrast, the higher ΔG°fold values obtained using DES-Amber and OL3R2.7 were unexpected, suggesting that key interactions are inaccurately described in the folded, unfolded, or misfolded ensembles. These discrepancies led us to further test DES-Amber and OL3R2.7 ffs on additional RNA and DNA systems, where further performance issues were observed. Our results emphasize the complexity of accurately modeling RNA dynamics and suggest that creating an RNA ff capable of reliably performing across a wide range of RNA systems remains extremely challenging. In conclusion, our study provides valuable insights into the capabilities of current RNA ffs and highlights key areas for future ff development.

• MeSH
• Nucleic Acid Conformation MeSH
• RNA * chemistry   MeSH
• Molecular Dynamics Simulation * MeSH
• Thermodynamics MeSH
• Publication type
• Journal Article MeSH
• Names of Substances
• RNA * MeSH

RNA recognition motifs (RRMs) are a key class of proteins that primarily bind single-stranded RNAs. In this study, we applied standard atomistic molecular dynamics simulations to obtain insights into the intricate binding dynamics between uridine-rich RNAs and TbRGG2 RRM using the recently developed OL3-Stafix AMBER force field, which improves the description of single-stranded RNA molecules. Complementing structural experiments that unveil a primary binding mode with a single uridine bound, our simulations uncover two supplementary binding modes in which adjacent nucleotides encroach upon the binding pocket. This leads to a unique molecular mechanism through which the TbRGG2 RRM is capable of rapidly transitioning the U-rich sequence. In contrast, the presence of non-native cytidines induces stalling and destabilization of the complex. By leveraging extensive equilibrium dynamics and a large variety of binding states, TbRGG2 RRM effectively expedites diffusion along the RNA substrate while ensuring robust selectivity for U-rich sequences despite featuring a solitary binding pocket. We further substantiate our description of the complex dynamics by simulating the fully spontaneous association process of U-rich sequences to the TbRGG2 RRM. Our study highlights the critical role of dynamics and auxiliary binding states in interface dynamics employed by RNA-binding proteins, which is not readily apparent in traditional structural studies but could represent a general type of binding strategy employed by many RNA-binding proteins.

• MeSH
• Nucleic Acid Conformation MeSH
• RNA Recognition Motif * MeSH
• RNA-Binding Proteins * chemistry   metabolism   MeSH
• RNA * metabolism   chemistry   MeSH
• Molecular Dynamics Simulation * MeSH
• Protein Binding MeSH
• Binding Sites MeSH
• Publication type
• Journal Article MeSH
• Names of Substances
• RNA-Binding Proteins * MeSH
• RNA * MeSH

The RNA chaperone Hfq plays crucial roles in bacterial gene expression and is a major facilitator of small regulatory RNA (sRNA) action. The toroidal architecture of the Hfq hexamer presents three well-characterized surfaces that allow it to bind sRNAs to stabilize them and engage target transcripts. Hfq-interacting sRNAs are categorized into two classes based on the surfaces they use to bind Hfq. By characterizing a systematic alanine mutant library of Hfq to identify amino acid residues that impact survival of Escherichia coli experiencing nitrogen (N) starvation, we corroborated the important role of the three RNA-binding surfaces for Hfq function. We uncovered two, previously uncharacterized, conserved residues, V22 and G34, in the hydrophobic core of Hfq, to have a profound impact on Hfq's RNA-binding activity in vivo. Transcriptome-scale analysis revealed that V22A and G34A Hfq mutants cause widespread destabilization of both sRNA classes, to the same extent as seen in bacteria devoid of Hfq. However, the alanine substitutions at these residues resulted in only modest alteration in stability and structure of Hfq. We propose that V22 and G34 have impact on Hfq function, especially critical under cellular conditions when there is an increased demand for Hfq, such as N starvation.

• MeSH
• RNA, Bacterial * metabolism   genetics   MeSH
• Nitrogen metabolism   MeSH
• Escherichia coli genetics   metabolism   MeSH
• Hydrophobic and Hydrophilic Interactions MeSH
• Conserved Sequence MeSH
• RNA, Small Untranslated * metabolism   genetics   MeSH
• Molecular Chaperones genetics   metabolism   chemistry   MeSH
• Mutation MeSH
• Host Factor 1 Protein * genetics   chemistry   metabolism   MeSH
• Escherichia coli Proteins * genetics   chemistry   metabolism   MeSH
• Gene Expression Regulation, Bacterial MeSH
• RNA Stability * MeSH
• Gene Expression Profiling MeSH
• Transcriptome MeSH
• Publication type
• Journal Article MeSH
• Names of Substances
• RNA, Bacterial * MeSH
• Nitrogen MeSH
• Hfq protein, E coli MeSH Browser
• RNA, Small Untranslated * MeSH
• Molecular Chaperones MeSH
• Host Factor 1 Protein * MeSH
• Escherichia coli Proteins * MeSH

Molecular dynamics (MD) simulations represent an established tool to study RNA molecules. The outcome of MD studies depends, however, on the quality of the force field (ff). Here we suggest a correction for the widely used AMBER OL3 ff by adding a simple adjustment of the nonbonded parameters. The reparameterization of the Lennard-Jones potential for the -H8···O5'- and -H6···O5'- atom pairs addresses an intranucleotide steric clash occurring in the type 0 base-phosphate interaction (0BPh). The nonbonded fix (NBfix) modification of 0BPh interactions (NBfix0BPh modification) was tuned via a reweighting approach and subsequently tested using an extensive set of standard and enhanced sampling simulations of both unstructured and folded RNA motifs. The modification corrects minor but visible intranucleotide clash for the anti nucleobase conformation. We observed that structural ensembles of small RNA benchmark motifs simulated with the NBfix0BPh modification provide better agreement with experiments. No side effects of the modification were observed in standard simulations of larger structured RNA motifs. We suggest that the combination of OL3 RNA ff and NBfix0BPh modification is a viable option to improve RNA MD simulations.

• MeSH
• Phosphates * MeSH
• Molecular Conformation MeSH
• Nucleotide Motifs MeSH
• RNA * chemistry   MeSH
• Molecular Dynamics Simulation MeSH
• Publication type
• Journal Article MeSH
• Names of Substances
• Phosphates * MeSH
• RNA * MeSH

Holliday junction (HJ) is a noncanonical four-way DNA structure with a prominent role in DNA repair, recombination, and DNA nanotechnology. By rearranging its four arms, HJ can adopt either closed or open state. With enzymes typically recognizing only a single state, acquiring detailed knowledge of the rearrangement process is an important step toward fully understanding the biological function of HJs. Here, we carried out standard all-atom molecular dynamics (MD) simulations of the spontaneous opening-closing transitions, which revealed complex conformational transitions of HJs with an involvement of previously unconsidered "half-closed" intermediates. Detailed free-energy landscapes of the transitions were obtained by sophisticated enhanced sampling simulations. Because the force field overstabilizes the closed conformation of HJs, we developed a system-specific modification which for the first time allows the observation of spontaneous opening-closing HJ transitions in unbiased MD simulations and opens the possibilities for more accurate HJ computational studies of biological processes and nanomaterials.

• MeSH
• DNA * MeSH
• DNA, Cruciform * MeSH
• Molecular Conformation MeSH
• DNA Repair MeSH
• Publication type
• Journal Article MeSH
• Research Support, Non-U.S. Gov't MeSH
• Names of Substances
• DNA * MeSH
• DNA, Cruciform * MeSH