Abstract:Mixed double helices formed by RNA and DNA strands, commonly referred to as hybrid duplexes or hybrids, are essential in biological processes like transcription and reverse transcription. They are also important for their applications in CRISPR gene editing and nanotechnology. Yet, despite their significance, the hybrids have been seldom modeled by atomistic molecular dynamics methodology, and there is no benchmark study systematically assessing the force-field performance. Here, we present an extensive benchmark study of the hybrids using contemporary and commonly utilized pairwise additive and polarizable nucleic acid force fields. Our findings indicate that none of the available force-field choices accurately reproduces all the characteristic structural details of the hybrids. The AMBER force fields are unable to populate the C3′-endo (north) pucker of the DNA strand and underestimate inclination. CHARMM force field accurately describes the C3′-endo pucker and inclination but shows base pair instability. The polarizable force fields struggle with accurately reproducing the helical parameters. Some force-field combinations even demonstrate a discernible conflict between the RNA and DNA parameters. In this work, we offer a candid assessment of the force-field performance for mixed DNA/RNA duplexes. We provide guidance on selecting utilizable force-field combinations, as well as highlight potential pitfalls and best practices for obtaining optimal performance.

What problem does this paper attempt to address?

This paper aims to solve the problem of insufficient evaluation of force - field performance in molecular dynamics simulations of DNA/RNA hybrid duplexes. Specifically, the authors hope to systematically evaluate the performance of different force fields in simulating these hybrid duplexes, to determine which combinations of force fields can most accurately reproduce the experimentally observed structural features, and to provide guidance for future research. ### Main problems 1. **Lack of systematic force - field evaluation**: Although DNA/RNA hybrid duplexes are of great significance in biological processes (such as transcription and reverse transcription) and application fields (such as CRISPR gene editing and nanotechnology), there is currently a lack of comprehensive force - field performance evaluation for them. 2. **Limitations of existing force fields**: Existing force fields have many problems when simulating DNA/RNA hybrid duplexes, for example: - The AMBER force field cannot correctly simulate the C3′ - endo conformation (pucker) of the DNA strand, resulting in simulation results biased towards the B - form structure. - Although the CHARMM force field can better describe the C3′ - endo conformation and inclination, it has problems in base - pair stability. - Polarized force fields also face challenges in the accuracy of helical parameters. ### Research objectives - **Evaluate the performance of different force fields**: Conduct extensive benchmark tests on DNA/RNA hybrid duplexes by using currently commonly used pairwise additive and polarized nucleic acid force fields. - **Provide guidance for force - field selection**: Based on the evaluation results, provide researchers with suggestions for selecting appropriate force - field combinations, and at the same time point out potential pitfalls and best practice methods. - **Improve force - field parameters**: Through detailed analysis of the performance of different force fields, provide references for future force - field adjustments to better simulate the structural features of DNA/RNA hybrid duplexes. ### Method overview - **Initial structure selection**: Use the X - ray structures of the Dickerson - Drew B - DNA dodecamer (DD) and Polypurine tract (PPT) hybrid duplexes as starting structures. - **Simulation setup**: Use a variety of non - polarized and polarized force fields for molecular dynamics simulations, including force fields such as AMBER, CHARMM, Drude, and AMOEBA. - **Data analysis**: Calculate and compare helical parameters (such as buckle, propeller, inclination, x - displacement, roll, etc.) under different force fields, and compare them with experimental data. Through these studies, the authors hope to fill the gap in the evaluation of DNA/RNA hybrid duplex force fields and provide valuable references for research in related fields.

Comprehensive Assessment of Force-Field Performance in Molecular Dynamics Simulations of DNA/RNA Hybrid Duplexes

Sensitivity of the RNA Structure to Ion Conditions as Probed by Molecular Dynamics Simulations of Common Canonical RNA Duplexes

Can We Ever Develop an Ideal RNA Force Field? Lessons Learned from Simulations of UUCG RNA Tetraloop and Other Systems.

Multipolar Description of Atom–Atom Electrostatic Interaction Energies in Single/Double-Stranded DNAs

Development of Force Field Parameters for the Simulation of Single- and Double-Stranded DNA Molecules and DNA–Protein Complexes

The development of nucleic acids force fields: From an unchallenged past to a competitive future

Extensive Evaluation of Force Fields for G-Quadruplexes

Multipolar Electrostatics for Hairpin and Pseudoknots in RNA: Improving the Accuracy of Force Field Potential Energy Function

Performance of Molecular Mechanics Force Fields for RNA Simulations: Stability of UUCG and GNRA Hairpins

Benchmarking RNA all-atom force fields using hairpin motifs

C-A test of DNA force fields

Simulations Meet Experiment to Reveal New Insights into DNA Intrinsic Mechanics

RNA force field with accuracy comparable to state-of-the-art protein force fields

Simulations of the molecular dynamics of nucleic acids

Structures and Dynamics of DNA Mini-Dumbbells Are Force Field Dependent

Systematic Comparison of Atomistic Force Fields for the Mechanical Properties of Double-Stranded DNA

Assessing Force Field Potential Energy Function Accuracy Via a Multipolar Description of Atomic Electrostatic Interactions in RNA.

Tumuc1: A New Accurate DNA Force Field Consistent with High-Level Quantum Chemistry

The static and dynamic structural heterogeneities of B-DNA: extending Calladine-Dickerson rules

Self-Consistent Parameterization of DNA Residues for the Non-Polarizable AMBER Force Fields

Coarse-grained modelling of DNA-RNA hybrids