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: Self-Consistent Parameterization of DNA Residues for the Non-Polarizable AMBER Force Fields
Fixed-charge (non-polarizable) forcefields are accurate and computationally efficient tools for modeling the molecular dynamics of nucleic acid polymers, particularly DNA, well into the µs timescale. The continued utility of these forcefields depends in part on expanding the residue set in step with advancing nucleic acid chemistry and biology. A key step in parameterizing new residues is charge derivation which is self-consistent with the existing residues. As atomic charges are derived by fitting against molecular electrostatic potentials, appropriate structural models are critical. Benchmarking against the existing charge set used in current AMBER nucleic acid forcefields, we report that quantum mechanical models of deoxynucleosides, even at a high level of theory, are not optimal structures for charge derivation. Instead, structures from molecular mechanics minimization yield charges with up to 6-fold lower RMS deviation from the published values, due to the choice of such an approach in the derivation of the original charge set. We present a contemporary protocol for rendering self-consistent charges as well as optimized charges for a panel of nine non-canonical residues that will permit comparison with literature as well as studying the dynamics of novel DNA polymers.  more » « less
Award ID(s):
2028902
PAR ID:
10350325
Author(s) / Creator(s):
; ; ; ;
Date Published:
Journal Name:
Life
Volume:
12
Issue:
5
ISSN:
2075-1729
Page Range / eLocation ID:
666
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Xeno-nucleic acids (XNAs) are synthetic genetic polymers with backbone structures composed of non-ribose or non-deoxyribose sugars. Phosphonomethylthreosyl nucleic acid (pTNA), a type of XNA that does not base pair with DNA or RNA, has been suggested as a possible genetic material for storing synthetic biology information in cells. A critical step in this process is the synthesis of XNA episomes using laboratory-evolved polymerases to copy DNA information into XNA. Here, we investigate the polymerase recognition of pTNA nucleotides using X-ray crystallography to capture the post-catalytic complex of engineered polymerases following the sequential addition of two pTNA nucleotides onto the 3′-end of a DNA primer. High-resolution crystal structures reveal that the polymerase mediates Watson–Crick base pairing between the extended pTNA adducts and the DNA template. Comparative analysis studies demonstrate that the sugar conformation and backbone position of pTNA are structurally more similar to threose nucleic acid than DNA even though pTNA and DNA share the same six-atom backbone repeat length. Collectively, these findings provide new insight into the structural determinants that guide the enzymatic synthesis of an orthogonal genetic polymer, and may lead to the discovery of new variants that function with enhanced activity. 
    more » « less
  2. Peptide amphiphiles are a class of molecules that can self-assemble into a variety of supramolecular structures, including high-aspect-ratio nanofibers. It is challenging to model and predict the charges in these supramolecular nanofibers because the ionization state of the peptides are not fixed but liable to change due to the acid-base equilibrium that is coupled to the structural organization of the peptide amphiphile molecules. Here, we have developed a theoretical model to describe and predict the amount of charge found on self-assembled peptide amphiphiles as a function of pH and ion concentration. In particular, we computed the amount of charge of peptide amphiphiles nanofibers with the sequence C 16 − V 2 A 2 E 2 . In our theoretical formulation, we consider charge regulation of the carboxylic acid groups, which involves the acid-base chemical equilibrium of the glutamic acid residues and the possibility of ion condensation. The charge regulation is coupled with the local dielectric environment by allowing for a varying dielectric constant that also includes a position-dependent electrostatic solvation energy for the charged species. We find that the charges on the glutamic acid residues of the peptide amphiphile nanofiber are much lower than the same functional group in aqueous solution. There is a strong coupling between the charging via the acid-base equilibrium and the local dielectric environment. Our model predicts a much lower degree of deprotonation for a position-dependent relative dielectric constant compared to a constant dielectric background. Furthermore, the shape and size of the electrostatic potential as well as the counterion distribution are quantitatively and qualitatively different. These results indicate that an accurate model of peptide amphiphile self-assembly must take into account charge regulation of acidic groups through acid–base equilibria and ion condensation, as well as coupling to the local dielectric environment. 
    more » « less
  3. 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
  4. Living cells regulate the dynamics of developmental events through interconnected signaling systems that activate and deactivate inert precursors. This suggests that similarly, synthetic biomaterials could be designed to develop over time by using chemical reaction networks to regulate the availability of assembling components. Here we demonstrate how the sequential activation or deactivation of distinct DNA building blocks can be modularly coordinated to form distinct populations of self-assembling polymers using a transcriptional signaling cascade of synthetic genes. Our building blocks are DNA tiles that polymerize into nanotubes, and whose assembly can be controlled by RNA molecules produced by synthetic genes that target the tile interaction domains. To achieve different RNA production rates, we use a strategy based on promoter “nicking” and strand displacement. By changing the way the genes are cascaded and the RNA levels, we demonstrate that we can obtain spatially and temporally different outcomes in nanotube assembly, including random DNA polymers, block polymers, and as well as distinct autonomous formation and dissolution of distinct polymer populations. Our work demonstrates a way to construct autonomous supramolecular materials whose properties depend on the timing of molecular instructions for self-assembly, and can be immediately extended to a variety of other nucleic acid circuits and assemblies. 
    more » « less
  5. Abstract The RNA programmed non-specific (trans) nuclease activity of CRISPR-Cas Type V and VI systems has opened a new era in the field of nucleic acid-based detection. Here, we report on the enhancement of trans-cleavage activity of Cas12a enzymes using hairpin DNA sequences as FRET-based reporters. We discover faster rate of trans-cleavage activity of Cas12a due to its improved affinity (Km) for hairpin DNA structures, and provide mechanistic insights of our findings through Molecular Dynamics simulations. Using hairpin DNA probes we significantly enhance FRET-based signal transduction compared to the widely used linear single stranded DNA reporters. Our signal transduction enables faster detection of clinically relevant double stranded DNA targets with improved sensitivity and specificity either in the presence or in the absence of an upstream pre-amplification step. 
    more » « less