skip to main content


Title: Challenges with Simulating Modified RNA: Insights into Role and Reciprocity of Experimental and Computational Approaches
RNA is critical to a broad spectrum of biological and viral processes. This functional diversity is a result of their dynamic nature; the variety of three-dimensional structures that they can fold into; and a host of post-transcriptional chemical modifications. While there are many experimental techniques to study the structural dynamics of biomolecules, molecular dynamics simulations (MDS) play a significant role in complementing experimental data and providing mechanistic insights. The accuracy of the results obtained from MDS is determined by the underlying physical models i.e., the force-fields, that steer the simulations. Though RNA force-fields have received a lot of attention in the last decade, they still lag compared to their protein counterparts. The chemical diversity imparted by the RNA modifications adds another layer of complexity to an already challenging problem. Insight into the effect of RNA modifications upon RNA folding and dynamics is lacking due to the insufficiency or absence of relevant experimental data. This review provides an overview of the state of MDS of modified RNA, focusing on the challenges in parameterization of RNA modifications as well as insights into relevant reference experiments necessary for their calibration.  more » « less
Award ID(s):
1651877
NSF-PAR ID:
10328456
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Genes
Volume:
13
Issue:
3
ISSN:
2073-4425
Page Range / eLocation ID:
540
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Wei, Guanghong (Ed.)

    Biomolecular condensates are important structures in various cellular processes but are challenging to study using traditional experimental techniques. In silico simulations with residue-level coarse-grained models strike a balance between computational efficiency and chemical accuracy. They could offer valuable insights by connecting the emergent properties of these complex systems with molecular sequences. However, existing coarse-grained models often lack easy-to-follow tutorials and are implemented in software that is not optimal for condensate simulations. To address these issues, we introduce OpenABC, a software package that greatly simplifies the setup and execution of coarse-grained condensate simulations with multiple force fields using Python scripting. OpenABC seamlessly integrates with the OpenMM molecular dynamics engine, enabling efficient simulations with performance on a single GPU that rivals the speed achieved by hundreds of CPUs. We also provide tools that convert coarse-grained configurations to all-atom structures for atomistic simulations. We anticipate that OpenABC will significantly facilitate the adoption of in silico simulations by a broader community to investigate the structural and dynamical properties of condensates.

     
    more » « less
  2. Peptide nucleic acids (PNAs) are a promising group of synthetic analogues of DNA and RNA that offer several distinct advantages over the naturally occurring nucleic acids for applications in biosensing, drug delivery, and nanoelectronics. Because of its structural differences from DNA/RNA, methods to analyze and assess the structure, conformations, and dynamics are needed. In this work, we develop synergistic techniques for the study of the PNA conformation. We use CuQ2, a Cu(II) complex with 8-hydroxyquinoline (HQ), as an alternative base pair and as a spin label in electron paramagnetic resonance (EPR) distance methods. We use molecular dynamics (MD) simulations with newly developed force field parameters for the spin labels to interpret the distance constraints determined by EPR. We complement these methods by UV–vis and circular dichroism measurements and assess the efficacy of the Cu(II) label on a PNA duplex whose backbone is based on aminoethylglycine and a duplex with a hydroxymethyl backbone modification. We show that the Cu(II) label functions efficiently within the standard PNA and the hydroxymethyl-modified PNA and that the MD parameters may be used to accurately reproduce our EPR findings. Through the combination of EPR and MD, we gain new insights into the PNA structure and conformations as well as into the mechanism of orientational selectivity in Cu(II) EPR at X-band. These results present for the first time a rigid Cu(II) spin label used for EPR distance measurements in PNA and the accompanying MD force fields for the spin label. Our studies also reveal that the spin labels have a low impact on the structure of the PNA duplexes. The combined MD and EPR approach represents an important new tool for the characterization of the PNA duplex structure and provides valuable information to aid in the rational application of PNA at large. 
    more » « less
  3. Organic chromophores initiate much of daytime aqueous phase chemistry in the environment. Thus, studying the absorption spectra of commonly used organic photosensitizers is paramount to fully understand their relevance in environmental processes. In this work, we combined UV-Vis spectroscopy, 1 H-NMR spectroscopy, quantum chemical calculations, and molecular dynamics simulations to investigate the absorption spectra of 4-benzoyl benzoic acid (4BBA), a widely used photosensitizer and a common proxy of environmentally relevant chromophores. Solutions of 4BBA at different pH values show that protonated and deprotonated species have an effect on its absorbance spectra. Theoretical calculations of these species in water clusters provide physical and chemical insights into the spectra. Quantum chemical calculations were conducted to analyze the UV-Vis absorbance spectra of 4BBA species using various cluster sizes, such as C 6 H 5 COC 6 H 4 COOH·(H 2 O) n , where n = 8 for relatively small clusters and n = 30 for larger clusters. While relatively small clusters have been successfully used for smaller chromophores, our results indicate that simulations of protonated species of 4BBA require relatively larger clusters of n = 30. A comparison between the experimental and theoretical results shows good agreement in the pH-dependent spectral shift between the hydrated cluster model and the experimental data. Overall, the theoretical and empirical results indicate that the experimental optical spectra of aqueous phase 4BBA can be represented by the acid–base equilibrium of the keto-forms, with a spectroscopically measured p K a of 3.41 ± 0.04. The results summarized here contribute to a molecular-level understanding of solvated organic molecules through calculations restricted to cluster models, and thereby, broader insight into environmentally relevant chromophores. 
    more » « less
  4. de Groot, Bert L. (Ed.)
    Intrinsically disordered proteins (IDPs) are highly dynamic systems that play an important role in cell signaling processes and their misfunction often causes human disease. Proper understanding of IDP function not only requires the realistic characterization of their three-dimensional conformational ensembles at atomic-level resolution but also of the time scales of interconversion between their conformational substates. Large sets of experimental data are often used in combination with molecular modeling to restrain or bias models to improve agreement with experiment. It is shown here for the N-terminal transactivation domain of p53 (p53TAD) and Pup, which are two IDPs that fold upon binding to their targets, how the latest advancements in molecular dynamics (MD) simulations methodology produces native conformational ensembles by combining replica exchange with series of microsecond MD simulations. They closely reproduce experimental data at the global conformational ensemble level, in terms of the distribution properties of the radius of gyration tensor, and at the local level, in terms of NMR properties including 15 N spin relaxation, without the need for reweighting. Further inspection revealed that 10–20% of the individual MD trajectories display the formation of secondary structures not observed in the experimental NMR data. The IDP ensembles were analyzed by graph theory to identify dominant inter-residue contact clusters and characteristic amino-acid contact propensities. These findings indicate that modern MD force fields with residue-specific backbone potentials can produce highly realistic IDP ensembles sampling a hierarchy of nano- and picosecond time scales providing new insights into their biological function. 
    more » « less
  5. Abstract

    Base stacking interactions between adjacent bases in DNA and RNA are important for many biological processes and in biotechnology applications. Previous work has estimated stacking energies between pairs of bases, but contributions of individual bases has remained unknown. Here, we use a Centrifuge Force Microscope for high-throughput single molecule experiments to measure stacking energies between adjacent bases. We found stacking energies strongest between purines (G|A at −2.3 ± 0.2 kcal/mol) and weakest between pyrimidines (C|T at −0.5 ± 0.1 kcal/mol). Hybrid stacking with phosphorylated, methylated, and RNA nucleotides had no measurable effect, but a fluorophore modification reduced stacking energy. We experimentally show that base stacking can influence stability of a DNA nanostructure, modulate kinetics of enzymatic ligation, and assess accuracy of force fields in molecular dynamics simulations. Our results provide insights into fundamental DNA interactions that are critical in biology and can inform design in biotechnology applications.

     
    more » « less