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
• konformace nukleové kyseliny MeSH
• RNA * chemie   MeSH
• simulace molekulární dynamiky * MeSH
• termodynamika MeSH
• Publikační typ
• časopisecké články MeSH
• Názvy látek
• 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.

• Klíčová slova
• ionizable lipid, lipid nanocarrier, lipid nanoparticle, liposome, molecular simulation, vesicle,
• MeSH
• léčivé přípravky chemie   MeSH
• lékové transportní systémy metody   MeSH
• lidé MeSH
• lipidy * chemie   MeSH
• nanočástice chemie   MeSH
• nosiče léků chemie   MeSH
• počítačová simulace MeSH
• simulace molekulární dynamiky * MeSH
• strojové učení MeSH
• Check Tag
• lidé MeSH
• Publikační typ
• časopisecké články MeSH
• přehledy MeSH
• Názvy látek
• léčivé přípravky MeSH
• lipidy * MeSH
• nosiče léků MeSH

Guanine quadruplexes (GQs) play crucial roles in various biological processes, and understanding their folding pathways provides insight into their stability, dynamics, and functions. This knowledge aids in designing therapeutic strategies, as GQs are potential targets for anticancer drugs and other therapeutics. Although experimental and theoretical techniques have provided valuable insights into different stages of the GQ folding, the structural complexity of GQs poses significant challenges, and our understanding remains incomplete. This study introduces a novel computational protocol for folding an entire GQ from single-strand conformation to its native state. By combining two complementary enhanced sampling techniques, we were able to model folding pathways, encompassing a diverse range of intermediates. Although our investigation of the GQ free energy surface (FES) is focused solely on the folding of the all-anti parallel GQ topology, this protocol has the potential to be adapted for the folding of systems with more complex folding landscapes.

• Klíčová slova
• DNA quadruplex, computational folding, enhanced sampling, kinetic partitioning mechanism, metadynamics, molecular dynamics, nudged elastic band, pathCV, transition path sampling,
• MeSH
• DNA chemie   MeSH
• G-kvadruplexy * MeSH
• konformace nukleové kyseliny MeSH
• simulace molekulární dynamiky MeSH
• termodynamika MeSH
• Publikační typ
• časopisecké články MeSH
• Názvy látek
• DNA 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
• fosfáty * MeSH
• molekulární konformace MeSH
• nukleotidové motivy MeSH
• RNA * chemie   MeSH
• simulace molekulární dynamiky MeSH
• Publikační typ
• časopisecké články MeSH
• Názvy látek
• fosfáty * MeSH
• RNA * MeSH

Recognition of single-stranded RNA (ssRNA) by RNA recognition motif (RRM) domains is an important class of protein-RNA interactions. Many such complexes were characterized using nuclear magnetic resonance (NMR) and/or X-ray crystallography techniques, revealing ensemble-averaged pictures of the bound states. However, it is becoming widely accepted that better understanding of protein-RNA interactions would be obtained from ensemble descriptions. Indeed, earlier molecular dynamics simulations of bound states indicated visible dynamics at the RNA-RRM interfaces. Here, we report the first atomistic simulation study of spontaneous binding of short RNA sequences to RRM domains of HuR and SRSF1 proteins. Using a millisecond-scale aggregate ensemble of unbiased simulations, we were able to observe a few dozen binding events. HuR RRM3 utilizes a pre-binding state to navigate the RNA sequence to its partially disordered bound state and then to dynamically scan its different binding registers. SRSF1 RRM2 binding is more straightforward but still multiple-pathway. The present study necessitated development of a goal-specific force field modification, scaling down the intramolecular van der Waals interactions of the RNA which also improves description of the RNA-RRM bound state. Our study opens up a new avenue for large-scale atomistic investigations of binding landscapes of protein-RNA complexes, and future perspectives of such research are discussed.

• MeSH
• HuR protein metabolismus   MeSH
• motiv rozpoznávající RNA genetika   MeSH
• proteiny vázající RNA * metabolismus   MeSH
• RNA * chemie   MeSH
• RRM proteiny metabolismus   MeSH
• vazba proteinů MeSH
• vazebná místa MeSH
• Publikační typ
• časopisecké články MeSH
• práce podpořená grantem MeSH
• Názvy látek
• HuR protein MeSH
• proteiny vázající RNA * MeSH
• RNA * MeSH
• RRM proteiny MeSH

The capability of current force fields to reproduce RNA structural dynamics is limited. Several methods have been developed to take advantage of experimental data in order to enforce agreement with experiments. Here, we extend an existing framework which allows arbitrarily chosen force-field correction terms to be fitted by quantification of the discrepancy between observables back-calculated from simulation and corresponding experiments. We apply a robust regularization protocol to avoid overfitting and additionally introduce and compare a number of different regularization strategies, namely, L1, L2, Kish size, relative Kish size, and relative entropy penalties. The training set includes a GACC tetramer as well as more challenging systems, namely, gcGAGAgc and gcUUCGgc RNA tetraloops. Specific intramolecular hydrogen bonds in the AMBER RNA force field are corrected with automatically determined parameters that we call gHBfixopt. A validation involving a separate simulation of a system present in the training set (gcUUCGgc) and new systems not seen during training (CAAU and UUUU tetramers) displays improvements regarding the native population of the tetraloop as well as good agreement with NMR experiments for tetramers when using the new parameters. Then, we simulate folded RNAs (a kink-turn and L1 stalk rRNA) including hydrogen bond types not sufficiently present in the training set. This allows a final modification of the parameter set which is named gHBfix21 and is suggested to be applicable to a wider range of RNA systems.

• MeSH
• RNA ribozomální MeSH
• RNA * chemie   MeSH
• simulace molekulární dynamiky * MeSH
• vodík MeSH
• vodíková vazba MeSH
• Publikační typ
• časopisecké články MeSH
• Názvy látek
• RNA ribozomální MeSH
• RNA * MeSH
• vodík MeSH

High-yielding and selective prebiotic syntheses of RNA and DNA nucleotides involve UV irradiation to promote the key reaction steps and eradicate biologically irrelevant isomers. While these syntheses were likely enabled by UV-rich prebiotic environment, UV-induced formation of photodamages in polymeric nucleic acids, such as cyclobutane pyrimidine dimers (CPDs), remains the key unresolved issue for the origins of RNA and DNA on Earth. Here, we demonstrate that substitution of adenine with 2,6-diaminopurine enables repair of CPDs with yields reaching 92%. This substantial self-repairing activity originates from excellent electron donating properties of 2,6-diaminopurine in nucleic acid strands. We also show that the deoxyribonucleosides of 2,6-diaminopurine and adenine can be formed under the same prebiotic conditions. Considering that 2,6-diaminopurine was previously shown to increase the rate of nonenzymatic RNA replication, this nucleobase could have played critical roles in the formation of functional and photostable RNA/DNA oligomers in UV-rich prebiotic environments.

• MeSH
• 2-aminopurin analogy a deriváty   farmakologie   MeSH
• adenin MeSH
• DNA účinky léků   účinky záření   MeSH
• nukleotidy MeSH
• nukleové kyseliny MeSH
• oprava DNA účinky léků   MeSH
• pyrimidinové dimery MeSH
• RNA chemie   MeSH
• simulace molekulární dynamiky MeSH
• ultrafialové záření škodlivé účinky   MeSH
• Publikační typ
• časopisecké články MeSH
• práce podpořená grantem MeSH
• Názvy látek
• 2-aminopurin MeSH
• 2,6-diaminopurine MeSH Prohlížeč
• adenin MeSH
• DNA MeSH
• nukleotidy MeSH
• nukleové kyseliny MeSH
• pyrimidinové dimery MeSH
• RNA MeSH

The widespread Mn2+-sensing yybP-ykoY riboswitch controls the expression of bacterial Mn2+ homeostasis genes. Here, we first determine the crystal structure of the ligand-bound yybP-ykoY riboswitch aptamer from Xanthomonas oryzae at 2.96 Å resolution, revealing two conformations with docked four-way junction (4WJ) and incompletely coordinated metal ions. In >100 µs of MD simulations, we observe that loss of divalents from the core triggers local structural perturbations in the adjacent docking interface, laying the foundation for signal transduction to the regulatory switch helix. Using single-molecule FRET, we unveil a previously unobserved extended 4WJ conformation that samples transient docked states in the presence of Mg2+. Only upon adding sub-millimolar Mn2+, however, can the 4WJ dock stably, a feature lost upon mutation of an adenosine contacting Mn2+ in the core. These observations illuminate how subtly differing ligand preferences of competing metal ions become amplified by the coupling of local with global RNA dynamics.

• MeSH
• bakteriální RNA chemie   genetika   metabolismus   MeSH
• Escherichia coli genetika   MeSH
• hořčík metabolismus   MeSH
• konformace nukleové kyseliny MeSH
• krystalografie rentgenová MeSH
• Lactococcus lactis genetika   metabolismus   MeSH
• ligandy MeSH
• mangan metabolismus   MeSH
• molekulární konformace MeSH
• molekulární modely MeSH
• mutace MeSH
• regulace genové exprese u bakterií MeSH
• riboswitch fyziologie   MeSH
• signální transdukce * MeSH
• simulace molekulární dynamiky MeSH
• vazebná místa MeSH
• Xanthomonas metabolismus   MeSH
• Publikační typ
• časopisecké články MeSH
• práce podpořená grantem MeSH
• Research Support, N.I.H., Extramural MeSH
• Research Support, U.S. Gov't, Non-P.H.S. MeSH
• Názvy látek
• bakteriální RNA MeSH
• hořčík MeSH
• ligandy MeSH
• mangan MeSH
• riboswitch MeSH

Guanine quadruplexes (G4s) are non-canonical nucleic acids structures common in important genomic regions. Parallel-stranded G4 folds are the most abundant, but their folding mechanism is not fully understood. Recent research highlighted that G4 DNA molecules fold via kinetic partitioning mechanism dominated by competition amongst diverse long-living G4 folds. The role of other intermediate species such as parallel G-triplexes and G-hairpins in the folding process has been a matter of debate. Here, we use standard and enhanced-sampling molecular dynamics simulations (total length of ∼0.9 ms) to study these potential folding intermediates. We suggest that parallel G-triplex per se is rather an unstable species that is in local equilibrium with a broad ensemble of triplex-like structures. The equilibrium is shifted to well-structured G-triplex by stacked aromatic ligand and to a lesser extent by flanking duplexes or nucleotides. Next, we study propeller loop formation in GGGAGGGAGGG, GGGAGGG and GGGTTAGGG sequences. We identify multiple folding pathways from different unfolded and misfolded structures leading towards an ensemble of intermediates called cross-like structures (cross-hairpins), thus providing atomistic level of description of the single-molecule folding events. In summary, the parallel G-triplex is a possible, but not mandatory short-living (transitory) intermediate in the folding of parallel-stranded G4.

• MeSH
• DNA chemie   genetika   metabolismus   MeSH
• G-kvadruplexy * MeSH
• guanin chemie   metabolismus   MeSH
• jednovláknová DNA chemie   genetika   metabolismus   MeSH
• kinetika MeSH
• konformace nukleové kyseliny * MeSH
• lidé MeSH
• sekvence nukleotidů MeSH
• simulace molekulární dynamiky * MeSH
• zvířata MeSH
• Check Tag
• lidé MeSH
• zvířata MeSH
• Publikační typ
• časopisecké články MeSH
• práce podpořená grantem MeSH
• Názvy látek
• DNA MeSH
• guanin MeSH
• jednovláknová DNA MeSH

Molecular dynamics (MD) simulations became a leading tool for investigation of structural dynamics of nucleic acids. Despite recent efforts to improve the empirical potentials (force fields, ffs), RNA ffs have persisting deficiencies, which hamper their utilization in quantitatively accurate simulations. Previous studies have shown that at least two salient problems contribute to difficulties in the description of free-energy landscapes of small RNA motifs: (i) excessive stabilization of the unfolded single-stranded RNA ensemble by intramolecular base-phosphate and sugar-phosphate interactions and (ii) destabilization of the native folded state by underestimation of stability of base pairing. Here, we introduce a general ff term (gHBfix) that can selectively fine-tune nonbonding interaction terms in RNA ffs, in particular, the H bonds. The gHBfix potential affects the pairwise interactions between all possible pairs of the specific atom types, while all other interactions remain intact; i.e., it is not a structure-based model. In order to probe the ability of the gHBfix potential to refine the ff nonbonded terms, we performed an extensive set of folding simulations of RNA tetranucleotides and tetraloops. On the basis of these data, we propose particular gHBfix parameters to modify the AMBER RNA ff. The suggested parametrization significantly improves the agreement between experimental data and the simulation conformational ensembles, although our current ff version still remains far from being flawless. While attempts to tune the RNA ffs by conventional reparametrizations of dihedral potentials or nonbonded terms can lead to major undesired side effects, as we demonstrate for some recently published ffs, gHBfix has a clear promising potential to improve the ff performance while avoiding introduction of major new imbalances.

• MeSH
• RNA chemie   MeSH
• simulace molekulární dynamiky * MeSH
• vodíková vazba MeSH
• Publikační typ
• časopisecké články MeSH
• Názvy látek
• RNA MeSH