Sequence-Dependent Kink Formation in Short DNA Loops: Theory and Molecular Dynamics Simulations

Cytosine methylation regulates DNA bendability depending on the curvature

Cytosine methylation plays an essential role in many biological processes, such as nucleosome inactivation and regulation of gene expression. The modulation of DNA mechanics may be one of the regulatory mechanisms influenced by cytosine methylation. However, it remains unclear how methylation influences DNA mechanics. Here, we show that methylation has contrasting effects on the bending property of dsDNA depending on DNA curvature. We directly applied bending force on 30 base pairs of dsDNA using a D-shaped DNA nanostructure and measured the degree of bending using single-molecule fluorescence resonance energy transfer without surface immobilization. When dsDNA is weakly bent, methylation increases the stiffness of dsDNA. The stiffness of dsDNA increased by approximately 8% with a single methylation site for 30 bp dsDNA. When dsDNA is highly bent by a strong force, it forms a kink, i.e., a sharp bending of dsDNA. Under strong bending, methylation destabilizes the non-kink form compared with the kink form, which makes dsDNA near the kink region apparently more bendable. However, if the kink region is methylated, the kink form is destabilized, and dsDNA becomes stiffer. As a result, methylation increases the stiffness of weakly bent dsDNA and concurrently can promote kink formation, which may stabilize the nucleosome structure. Our results provide new insight into the effect of methylation, showing that cytosine methylation has opposite effects on DNA mechanics depending on its curvature and methylation location.

Molecular dynamics simulation‐guided toehold mediated strand displacement probe for single‐nucleotide variants detection

Single nucleotide variant (SNV) has become an emerging biomarker for various diseases such as cancers and infectious diseases. Toehold‐mediated strand displacement (TMSD), the core reaction of DNA nanotechnology, has been widely leveraged to identify SNVs. However, inappropriate choice of mismatch location results in poor discrimination ability. Here, we comprehensively investigate the effect of mismatch location on TMSD kinetics by molecular dynamic simulation tool oxDNA through umbrella sampling and forward flux sampling disclosing that mismatches at the border of the toehold and branch migration domain yield the lowest TMSD reaction rate. Nine disease‐related SNVs (SARS‐CoV‐2‐D614G, EGFR‐L858R, EGFR‐T790M, KRAS‐G12R, etc.) were tested experimentally showing a good agreement with simulation. The best choice of mismatch location enables high discrimination factor with a median of 124 for SNV and wild type. Coupling with a probe‐sink system, a low variant allele frequency of 0.1% was detected with 3 S/N. We successfully used the probes to detect SNVs with high confidence in the PCR clones of constructed plasmids. This work provides mechanistic insights into TMSD process at the single‐nucleotide level and can be a guidance for the design of TMSD system with fine‐tuning kinetics for various applications in biosensors and nanotechnology.

A Primer on the oxDNA Model of DNA: When to Use it, How to Simulate it and How to Interpret the Results

The oxDNA model of Deoxyribonucleic acid has been applied widely to systems in biology, biophysics and nanotechnology. It is currently available via two independent open source packages. Here we present a set of clearly documented exemplar simulations that simultaneously provide both an introduction to simulating the model, and a review of the model's fundamental properties. We outline how simulation results can be interpreted in terms of-and feed into our understanding of-less detailed models that operate at larger length scales, and provide guidance on whether simulating a system with oxDNA is worthwhile.