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: Coarse-grained modelling of DNA plectoneme pinning in the presence of base-pair mismatches
Abstract Damaged or mismatched DNA bases result in the formation of physical defects in double-stranded DNA. In vivo, defects in DNA must be rapidly and efficiently repaired to maintain cellular function and integrity. Defects can also alter the mechanical response of DNA to bending and twisting constraints, both of which are important in defining the mechanics of DNA supercoiling. Here, we use coarse-grained molecular dynamics (MD) simulation and supporting statistical-mechanical theory to study the effect of mismatched base pairs on DNA supercoiling. Our simulations show that plectoneme pinning at the mismatch site is deterministic under conditions of relatively high force (>2 pN) and high salt concentration (>0.5 M NaCl). Under physiologically relevant conditions of lower force (0.3 pN) and lower salt concentration (0.2 M NaCl), we find that plectoneme pinning becomes probabilistic and the pinning probability increases with the mismatch size. These findings are in line with experimental observations. The simulation framework, validated with experimental results and supported by the theoretical predictions, provides a way to study the effect of defects on DNA supercoiling and the dynamics of supercoiling in molecular detail.  more » « less
Award ID(s):
2019745
PAR ID:
10233853
Author(s) / Creator(s):
; ; ; ;
Date Published:
Journal Name:
Nucleic Acids Research
Volume:
48
Issue:
19
ISSN:
0305-1048
Page Range / eLocation ID:
10713 to 10725
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Schlick, Tamar (Ed.)
    A model for DNA and nucleosomes is introduced with the goal of studying chromosomes from a single base level all the way to higher-order chromatin structures. This model, dubbed the Widely Editable Chromatin Model (WEChroM), reproduces the complex mechanics of the double helix including its bending persistence length and twisting persistence length, and the temperature dependence of the former. The WEChroM Hamiltonian is composed of chain connectivity, steric interactions, and associative memory terms representing all remaining interactions leading to the structure, dynamics, and mechanical characteristics of the B-DNA. Several applications of this model are discussed to demonstrate its applicability. WEChroM is used to investigate the behavior of circular DNA in the presence of positive and negative supercoiling. We show that it recapitulates the formation of plectonemes and of structural defects that relax mechanical stress. The model spontaneously manifests an asymmetric behavior with respect to positive or negative supercoiling, similar to what was previously observed in experiments. Additionally, we show that the associative memory Hamiltonian is also capable of reproducing the free energy of partial DNA unwrapping from nucleosomes. WEChroM is designed to emulate the continuously variable mechanical properties of the 10nm fiber and, by virtue of its simplicity, is ready to be scaled up to molecular systems large enough to investigate the structural ensembles of genes. WEChroM is implemented in the OpenMM simulation toolkits and is freely available for public use. 
    more » « less
  2. Force fields were developed for metal-mediated DNA (mmDNA) structures, using ab-initio methods to parameterize metal coordination. Two mmDNA were considered, comprising of a cytosine/thymine mismatch with coordinated Ag/Hg metal atoms. These basepairs were parameterized with the proposed computational framework and subjected to multiple validation steps. The generated force fields result in enhanced structural stability, with metallated basepairs rotating into the major groove. Our findings show a higher propeller angle associated with metalated base pair, which agrees with previously reported experimental data. Molecular dynamics (MD) simulations showed that the metallated basepairs stabilized the DNA structure, with the mismatch bases locking together via metal coordination. We anticipate the developed force fields can help in unveiling the structural dynamics of long metallo-DNA nanowires. 
    more » « less
  3. Transcription has a mechanical component, as the translocation of the transcription machinery or RNA polymerase (RNAP) on DNA or chromatin is dynamically coupled to the chromatin torsion. This posits chromatin mechanics as a possible regulator of eukaryotic transcription, however, the modes and mechanisms of this regulation are elusive. Here, we first take a statistical mechanics approach to model the torsional response of topology-constrained chromatin. Our model recapitulates the experimentally observed weaker torsional stiffness of chromatin compared to bare DNA and proposes structural transitions of nucleosomes into chirally distinct states as the driver of the contrasting torsional mechanics. Coupling chromatin mechanics with RNAP translocation in stochastic simulations, we reveal a complex interplay of DNA supercoiling and nucleosome dynamics in governing RNAP velocity. Nucleosomes play a dual role in controlling the transcription dynamics. The steric barrier aspect of nucleosomes in the gene body counteracts transcription via hindering RNAP motion, whereas the chiral transitions facilitate RNAP motion via driving a low restoring torque upon twisting the DNA. While nucleosomes with low dissociation rates are typically transcriptionally repressive, highly dynamic nucleosomes offer less of a steric barrier and enhance the transcription elongation dynamics of weakly transcribed genes via buffering DNA twist. We use the model to predict transcription-dependent levels of DNA supercoiling in segments of the budding yeast genome that are in accord with available experimental data. The model unveils a paradigm of DNA supercoiling-mediated interaction between genes and makes testable predictions that will guide experimental design. 
    more » « less
  4. ABSTRACT Formation of the feldspathoid sodalite (Na6Al6Si6O24·2NaCl) by reaction of nepheline (NaAlSiO4) with NaCl-bearing brines was investigated at 3 and 6 kbar and at a constant temperature of 750 °C to determine the brine concentration at which sodalite forms with variation in pressure. The reaction boundary was located by reaction-reversal experiments in the system NaAlSiO4–NaCl–H2O at a brine concentration of 0.16 ± 0.08 XNaCl [= molar ratio NaCl/(NaCl + H2O)] at 3 kbar and at a brine concentration of 0.35 ± 0.03 XNaCl at 6 kbar. Characterization of the sodalite using both X-ray diffraction and infrared spectroscopy after treatment in these brines indicated no obvious evidence of water or hydroxyl incorporation into the cage structure of sodalite. The data from this study were combined with earlier results by Wellman (1970) and Sharp et al. (1989) at lower (1–1.5 kbar) and higher (7–8 kbar) pressures, respectively, on sodalite formation from nepheline and NaCl which models as a concave-down curve in XNaCl – P space. In general, sodalite buffers the concentration of neutral aqueous NaCl° in the brine to relatively low values at P < 4 kbar, but NaCl° increases rapidly at higher pressures. Thermochemical modeling of these data was done to determine the activity of the aqueous NaCl° relative to a 1 molal (m) standard state, demonstrating very low activities (<0.2 m, or 1.2 wt.%) of NaCl° at 3 kbar and lower, but rising to relatively high activities (>20 m, or 54 wt.%) of NaCl° at 6 kbar or higher. The results from this study place constraints on the concentration of NaCl° in brines coexisting with nepheline and sodalite and, because of the relative insensitivity of this reaction to temperature, can provide a convenient geobarometer for those localities where the fluid compositions that formed nepheline and sodalite can be determined independently. 
    more » « less
  5. Peptide nucleic acid (PNA) clamps modified with 2-aminopyridine (M) nucleobase invaded double-stranded DNA under physiological salt conditions. In contrast, PNAs carrying common nucleobases could not fully invade DNA under these conditions. M-modified PNAs may overcome the problematic requirement for low salt concentration, a long-standing DNA invasion problem. 
    more » « less