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

Lipid-mediated delivery of active pharmaceutical ingredients (API) opened new possibilities in advanced therapies. By encapsulating an API into a lipid nanocarrier (LNC), one can safely deliver APIs not soluble in water, those with otherwise strong adverse effects, or very fragile ones such as nucleic acids. However, for the rational design of LNCs, a detailed understanding of the composition-structure-function relationships is missing. This review presents currently available computational methods for LNC investigation, screening, and design. The state-of-the-art physics-based approaches are described, with the focus on molecular dynamics simulations in all-atom and coarse-grained resolution. Their strengths and weaknesses are discussed, highlighting the aspects necessary for obtaining reliable results in the simulations. Furthermore, a machine learning, i.e., data-based learning, approach to the design of lipid-mediated API delivery is introduced. The data produced by the experimental and theoretical approaches provide valuable insights. Processing these data can help optimize the design of LNCs for better performance. In the final section of this Review, state-of-the-art of computer simulations of LNCs are reviewed, specifically addressing the compatibility of experimental and computational insights.

• Keywords
• ionizable lipid, lipid nanocarrier, lipid nanoparticle, liposome, molecular simulation, vesicle,
• MeSH
• Pharmaceutical Preparations chemistry   administration & dosage   MeSH
• Drug Delivery Systems * methods   MeSH
• Humans MeSH
• Lipids * chemistry   MeSH
• Nanoparticles chemistry   MeSH
• Drug Carriers * chemistry   MeSH
• Computer Simulation MeSH
• Molecular Dynamics Simulation MeSH
• Machine Learning MeSH
• Check Tag
• Humans MeSH
• Publication type
• Journal Article MeSH
• Review MeSH
• Names of Substances
• Pharmaceutical Preparations MeSH
• Lipids * MeSH
• Drug Carriers * MeSH

Doxorubicin (DOX) is a widely used chemotherapeutic agent known for intercalating into DNA. However, the exact modes of DOX interactions with various DNA structures remain unclear. Using molecular dynamics (MD) simulations, we explored DOX interactions with DNA duplexes (dsDNA), G-quadruplex, and nucleosome. DOX predominantly stacks on terminal bases of dsDNA and occasionally binds into its minor groove. In the G-quadruplex, DOX stacks on planar tetrads but does not spontaneously intercalate into these structures. Potential of mean force calculations indicate that while intercalation is the most energetically favorable interaction mode for DOX in dsDNA, the process requires overcoming a significant energy barrier. In contrast, DOX spontaneously intercalates into bent nucleosomal DNA, due to the increased torsional stress. This preferential intercalation of DOX into regions with higher torsional stress provides new insights into its mechanism of action and underscores the importance of DNA tertiary and quaternary structures in therapies utilizing DNA intercalation.

• Keywords
• DNA, G‐quadruplex, MD simulations, doxorubicin intercalation, intercalation, nucleosome,
• MeSH
• DNA * chemistry   MeSH
• Doxorubicin * chemistry   MeSH
• G-Quadruplexes MeSH
• Intercalating Agents chemistry   MeSH
• Nucleic Acid Conformation MeSH
• Nucleosomes * chemistry   MeSH
• Antibiotics, Antineoplastic * chemistry   MeSH
• Molecular Dynamics Simulation * MeSH
• Publication type
• Journal Article MeSH
• Names of Substances
• DNA * MeSH
• Doxorubicin * MeSH
• Intercalating Agents MeSH
• Nucleosomes * MeSH
• Antibiotics, Antineoplastic * MeSH

The transition from B-DNA to A-DNA occurs in many protein-DNA interactions or in DNA/RNA hybrid duplexes, and thus plays a role in many important biomolecular processes that convey the biological function of DNA. However, the stability of A-DNA is severely underestimated in current AMBER force fields such as OL15, OL21 or bsc1, potentially leading to unstable or deformed protein-DNA complexes. In this study, we refine the deoxyribose dihedral potential to increase the stability of the north (N) puckering present in A-DNA. The new parameters, termed OL24, model A/B equilibrium in B-DNA duplexes in water in good agreement with nuclear magnetic resonance (NMR) experiment. They also improve the description of DNA/RNA hybrids and the transition of the DNA duplex to the A-form in concentrated ethanol solutions. These refinements significantly improve the modeling of protein-DNA complexes, increasing their structural stability and A-form population, while maintaining accurate representation of canonical B-DNA duplexes. Overall, the new parameters should allow more reliable modeling of the thermodynamic equilibrium between A- and B-DNA forms and the interactions of DNA with proteins.

• MeSH
• DNA, A-Form * chemistry   MeSH
• DNA, B-Form * chemistry   MeSH
• Deoxyribose * chemistry   MeSH
• DNA * chemistry   MeSH
• Nucleic Acid Conformation MeSH
• Molecular Dynamics Simulation MeSH
• Thermodynamics MeSH
• Publication type
• Journal Article MeSH
• Names of Substances
• DNA, A-Form * MeSH
• DNA, B-Form * MeSH
• Deoxyribose * MeSH
• DNA * MeSH

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 hybrid duplexes 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 polypurine tract (PPT) and Dickerson-Drew dodecamer hybrid duplexes 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 hybrid duplexes. The AMBER force fields are unable to populate the C3'-endo (north) pucker of the DNA strand and underestimate inclination. The 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 and also highlight potential pitfalls and best practices for obtaining optimal performance.

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

Guanine quadruplexes (GQs) are non-canonical nucleic acid structures involved in many biological processes. GQs formed in single-stranded regions often need to be unwound by cellular machinery, so their mechanochemical properties are important. Here, we performed steered molecular dynamics simulations of human telomeric GQs to study their unfolding. We examined four pulling regimes, including a very slow setup with pulling velocity and force load accessible to high-speed atomic force microscopy. We identified multiple factors affecting the unfolding mechanism, i.e.,: (i) the more the direction of force was perpendicular to the GQ channel axis (determined by GQ topology), the more the base unzipping mechanism happened, (ii) the more parallel the direction of force was, GQ opening and cross-like GQs were more likely to occur, (iii) strand slippage mechanism was possible for GQs with an all-anti pattern in a strand, and (iv) slower pulling velocity led to richer structural dynamics with sampling of more intermediates and partial refolding events. We also identified that a GQ may eventually unfold after a force drop under forces smaller than those that the GQ withstood before the drop. Finally, we found out that different unfolding intermediates could have very similar chain end-to-end distances, which reveals some limitations of structural interpretations of single-molecule spectroscopic data.

• MeSH
• G-Quadruplexes * MeSH
• Guanine * chemistry   MeSH
• Humans MeSH
• Mechanical Phenomena MeSH
• Molecular Dynamics Simulation MeSH
• Telomere MeSH
• Check Tag
• Humans MeSH
• Publication type
• Journal Article MeSH
• Research Support, Non-U.S. Gov't MeSH
• Names of Substances
• Guanine * MeSH

Advances in molecular dynamics (MD) software alongside enhanced computational power and hardware have allowed for MD simulations to significantly expand our knowledge of biomolecular structure, dynamics, and interactions. Furthermore, it has allowed for the extension of conformational sampling times from nanoseconds to the microsecond level and beyond. This has not only made convergence of conformational ensembles through comprehensive sampling possible but consequently exposed deficiencies and allowed the community to overcome limitations in the available force fields. The reproducibility and accuracy of the force fields are imperative in order to produce biologically relevant data. The Amber nucleic acid force fields have been used widely since the mid-1980s, and improvement of these force fields has been a community effort with several artifacts revealed, corrected, and reevaluated by various research groups. Here, we focus on the Amber force fields for use with double-stranded DNA and present the assessment of two recently developed force field parameter sets (OL21 and Tumuc1). Extensive MD simulations were performed with six test systems and two different water models. We observe the improvement of OL21 and Tumuc1 compared to previous generations of the Amber DNA force. We did not detect any significant improvement in the performance of Tumuc1 compared to OL21 despite the reparameterization of bonded force field terms in the former; however, we did note discrepancies in Tumuc1 when modeling Z-DNA sequences.

• MeSH
• DNA * chemistry   MeSH
• Molecular Conformation MeSH
• Reproducibility of Results MeSH
• Molecular Dynamics Simulation MeSH
• DNA, Z-Form * MeSH
• Publication type
• Journal Article MeSH
• Names of Substances
• DNA * MeSH
• DNA, Z-Form * MeSH