skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Nucleic acid folding simulations using a physics-based atomistic free energy model
Performing full-resolution atomistic simulations of nucleic acid folding has remained a challenge for biomolecular modeling. Understanding how nucleic acids fold and how they transition between different folded structures as they unfold and refold has important implications for biology. This paper reports a theoretical model and computer simulation of the ab initio folding of DNA inverted repeat sequences. The formulation is based on an all-atom conformational model of the sugar-phosphate backbone via chain closure, and it incorporates three major molecular-level driving forces—base stacking, counterion-induced backbone self-interactions, and base pairing—via separate analytical theories designed to capture and reproduce the effects of the solvent without requiring explicit water and ions in the simulation. To accelerate computational throughput, a mixed numerical/analytical algorithm for the calculation of the backbone conformational volume is incorporated into the Monte Carlo simulation, and special stochastic sampling techniques were employed to achieve the computational efficiency needed to fold nucleic acids from scratch. This paper describes implementation details, benchmark results, and the advantages and technical challenges with this approach.  more » « less
Award ID(s):
1664801
PAR ID:
10593488
Author(s) / Creator(s):
Publisher / Repository:
American Institute of Physics
Date Published:
Journal Name:
The Journal of Chemical Physics
Volume:
156
Issue:
17
ISSN:
0021-9606
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. RNA macromolecules, like proteins, fold to assume shapes that are intimately connected to their broadly recognized biological functions; however, because of their high charge and dynamic nature, RNA structures are far more challenging to determine. We introduce an approach that exploits the high brilliance of x-ray free-electron laser sources to reveal the formation and ready identification of angstrom-scale features in structured and unstructured RNAs. Previously unrecognized structural signatures of RNA secondary and tertiary structures are identified through wide-angle solution scattering experiments. With millisecond time resolution, we observe an RNA fold from a dynamically varying single strand through a base-paired intermediate to assume a triple-helix conformation. While the backbone orchestrates the folding, the final structure is locked in by base stacking. This method may help to rapidly characterize and identify structural elements in nucleic acids in both equilibrium and time-resolved experiments. 
    more » « less
  2. Many proteins have slow folding times in vitro that are physiologically untenable. To combat this challenge, ATP-dependent chaperonins are thought to possess the unique ability to catalyze protein folding. Performing quantitative model selection using protein folding and unfolding data, we here show that short nucleic acids containing G-quadruplex (G4) structure can also catalyze protein folding. Performing the experiments as a function of temperature demonstrates that the G4 reshapes the underlying driving forces of protein folding. As short nucleic acids can catalyze protein folding without the input of ATP, the ability of the cell to fold proteins is far higher than previously anticipated. 
    more » « less
  3. Abstract Electron paramagnetic resonance (EPR) has become an important tool to probe conformational changes in nucleic acids. An array of EPR labels for nucleic acids are available, but they often come at the cost of long tethers, are dependent on the presence of a particular nucleotide or can be placed only at the termini. Site directed incorporation of Cu2+-chelated to a ligand, 2,2′dipicolylamine (DPA) is potentially an attractive strategy for site-specific, nucleotide independent Cu2+-labelling in DNA. To fully understand the potential of this label, we undertook a systematic and detailed analysis of the Cu2+-DPA motif using EPR and molecular dynamics (MD) simulations. We used continuous wave EPR experiments to characterize Cu2+ binding to DPA as well as optimize Cu2+ loading conditions. We performed double electron-electron resonance (DEER) experiments at two frequencies to elucidate orientational selectivity effects. Furthermore, comparison of DEER and MD simulated distance distributions reveal a remarkable agreement in the most probable distances. The results illustrate the efficacy of the Cu2+-DPA in reporting on DNA backbone conformations for sufficiently long base pair separations. This labelling strategy can serve as an important tool for probing conformational changes in DNA upon interaction with other macromolecules. 
    more » « less
  4. Abstract Proteins and nucleic acids are key components in many processes in living cells, and interactions between proteins and nucleic acids are often crucial pathway components. In many cases, large flexibility of proteins as they interact with nucleic acids is key to their function. To understand the mechanisms of these processes, it is necessary to consider the 3D atomic structures of such protein–nucleic acid complexes. When such structures are not yet experimentally determined, protein docking can be used to computationally generate useful structure models. However, such docking has long had the limitation that the consideration of flexibility is usually limited to small movements or to small structures. We previously developed a method of flexible protein docking which could model ordered proteins which undergo large‐scale conformational changes, which we also showed was compatible with nucleic acids. Here, we elaborate on the ability of that pipeline, Flex‐LZerD, to model specifically interactions between proteins and nucleic acids, and demonstrate that Flex‐LZerD can model more interactions and types of conformational change than previously shown. 
    more » « less
  5. null (Ed.)
    Abstract The helical structures of DNA and RNA were originally revealed by experimental data. Likewise, the development of programs for modeling these natural polymers was guided by known structures. These nucleic acid polymers represent only two members of a potentially vast class of polymers with similar structural features, but that differ from DNA and RNA in the backbone or nucleobases. Xeno nucleic acids (XNAs) incorporate alternative backbones that affect the conformational, chemical, and thermodynamic properties of XNAs. Given the vast chemical space of possible XNAs, computational modeling of alternative nucleic acids can accelerate the search for plausible nucleic acid analogs and guide their rational design. Additionally, a tool for the modeling of nucleic acids could help reveal what nucleic acid polymers may have existed before RNA in the early evolution of life. To aid the development of novel XNA polymers and the search for possible pre-RNA candidates, this article presents the proto-Nucleic Acid Builder (https://github.com/GT-NucleicAcids/pnab), an open-source program for modeling nucleic acid analogs with alternative backbones and nucleobases. The torsion-driven conformation search procedure implemented here predicts structures with good accuracy compared to experimental structures, and correctly demonstrates the correlation between the helical structure and the backbone conformation in DNA and RNA. 
    more » « less