Extreme bendability of DNA less than 100 base pairs long revealed by single-molecule cyclization

Kinetics and dynamics of oligonucleotide hybridization

The hybridization of short nucleic acid strands is a remarkable spontaneous process that is foundational to biotechnology and nanotechnology and plays a crucial role in gene expression, editing and DNA repair. Decades of research into the mechanism of hybridization have resulted in a deep understanding of its thermodynamics, but many questions remain regarding its kinetics and dynamics. Recent advances in experiments and molecular dynamics simulations of nucleic acids are enabling more direct insight into the structural dynamics of hybridization, which can test long-standing assumptions regarding its mechanism. In this Review, we summarize the current state of knowledge of hybridization kinetics, discuss the barriers to a molecular description of hybridization dynamics, and highlight the new approaches that have begun uncovering the dynamics of hybridization and the duplex ensemble. The kinetics and dynamics of hybridization are highly sensitive to the composition of nucleic acids, and we emphasize recent discoveries and open questions on the role of nucleobase sequence and chemical modifications.

Dynamics of nucleosomes and chromatin fibers revealed by single-molecule measurements

The nucleosome is the fundamental structural unit of chromosome fibers. DNA wraps around a histone octamer to form a nucleosome while neighboring nucleosomes interact to form higher-order structures and fit gigabase-long DNAs into a small volume of the nucleus. Nucleosomes interrupt the access of transcription factors to a genomic region and provide regulatory controls of gene expression. Biochemical and physical cues stimulate wrapping-unwrapping and condensation-decondensation dynamics of nucleosomes and nucleosome arrays. Nucleosome dynamics and chromatin fiber organization are influenced by changes in the ionic background within the nucleus, post-translational modifications of histone proteins, and DNA sequence characteristics, such as histone-binding motifs and nucleosome spacing. Biochemical and biophysical measurements, along with in silico simulations, have been extensively used to study the regulatory effects on chromatin dynamics. In particular, single-molecule measurements have revealed novel mechanistic details of nucleosome and chromatin dynamics. This minireview elucidates recent findings on chromatin dynamics from these approaches. [BMB Reports 2025; 58(1): 24-32].

Probing mechanical selection in diverse eukaryotic genomes through accurate prediction of 3D DNA mechanics

Connections between the mechanical properties of DNA and biological functions have been speculative due to the lack of methods to measure or predict DNA mechanics at scale. Recently, a proxy for DNA mechanics, cyclizability, was measured by loop-seq and enabled genome-scale investigation of DNA mechanics. Here, we use this dataset to build a computational model predicting bias-corrected intrinsic cyclizability, with near-perfect accuracy, solely based on DNA sequence. Further, the model predicts intrinsic bending direction in 3D space. Using this tool, we aimed to probe mechanical selection - that is, the evolutionary selection of DNA sequence based on its mechanical properties - in diverse circumstances. First, we found that the intrinsic bend direction of DNA sequences correlated with the observed bending in known protein-DNA complex structures, suggesting that many proteins co-evolved with their DNA partners to capture DNA in its intrinsically preferred bent conformation. We then applied our model to large-scale yeast population genetics data and showed that centromere DNA element II, whose consensus sequence is unknown, leaving its sequence-specific role unclear, is under mechanical selection to increase the stability of inner-kinetochore structure and to facilitate centromeric histone recruitment. Finally, in silico evolution under strong mechanical selection discovered hallucinated sequences with cyclizability values so extreme that they required experimental validation, yet, found in nature in the densely packed mitochondrial(mt) DNA of Namystynia karyoxenos, an ocean-dwelling protist with extreme mitochondrial gene fragmentation. The need to transmit an extraordinarily large amount of mtDNA, estimated to be > 600 Mb, in combination with the absence of mtDNA compaction proteins may have pushed mechanical selection to the extreme. Similarly extreme DNA mechanics are observed in bird microchromosomes, although the functional consequence is not yet clear. The discovery of eccentric DNA mechanics in unrelated unicellular and multicellular eukaryotes suggests that we can predict extreme natural biology which can arise through strong selection. Our methods offer a way to study the biological functions of DNA mechanics in any genome and to engineer DNA sequences with desired mechanical properties.

Bending DNA increases its helical repeat

In all biological systems, DNA is under high mechanical stress from bending and twisting. For example, DNA is tightly bent in nucleosome complexes, virus capsids, bacterial chromosomes, or complexes with transcription factors that regulate gene expression. A structurally and mechanically accurate model of DNA is therefore necessary to understand some of the most fundamental molecular mechanisms in biology including DNA packaging, replication, transcription and gene regulation. An iconic feature of DNA is its double helical nature with an average repeath 0 of ~10.45 base pairs per turn, which is commonly believed to be independent of curvature. We developed a ligation assay on nicked DNA circles of variable curvature that reveals a strong unwinding of DNA to over 11 bp/turn for radii around 3-4 nm. Our work constitutes a major modification of the standard mechanical model of DNA and requires reassessing the molecular mechanisms and energetics of all processes involving tightly bent DNA.

Converting Multiple- to Single-DNA-Tethered Beads and Removing Only-One-End-Tethered DNA in High-Throughput Stretching

S-DNA is a double-stranded DNA that forms under tensions of >65 pN. Here, we report that S-DNA resists the cleavage of Cas12a and the restriction endonuclease SmaI. Taking advantage of this resistance, in magnetic tweezer experiments, we developed an assay to convert multiple-DNA-tethered beads into single-DNA-tethered beads and remove the only-one-end-tethered DNA molecule by cleaving the DNA that does not transition to S-DNA at about 80 pN. When multiple DNA molecules are tethered to a single bead, they share the tension, exist in the B-form, and allow the cleavage. Only-one-end-tethered DNA molecules, free of tension, are also cleaved. In versatile types of experiments, we proved the broad applications of this assay: measuring the correct DNA elasticity and DNA condensation dynamics by avoiding the false results due to interference of only-one-end-tethered DNA molecules and quantifying the accurate cleavage rates of Cas12a and the restriction endonucleases by eliminating the error caused by multiple-DNA-tethered beads. This convenient assay ensures correct and accurate results in high-throughput DNA stretching experiments.

DNA packaging by molecular motors: from bacteriophage to human chromosomes

Dense packaging of genomic DNA is crucial for organismal survival, as DNA length always far exceeds the dimensions of the cells that contain it. Organisms, therefore, use sophisticated machineries to package their genomes. These systems range across kingdoms from a single ultra-powerful rotary motor that spools the DNA into a bacteriophage head, to hundreds of thousands of relatively weak molecular motors that coordinate the compaction of mitotic chromosomes in eukaryotic cells. Recent technological advances, such as DNA proximity-based sequencing approaches, polymer modelling and in vitro reconstitution of DNA loop extrusion, have shed light on the biological mechanisms driving DNA organization in different systems. Here, we discuss DNA packaging in bacteriophage, bacteria and eukaryotic cells, which, despite their extreme variation in size, structure and genomic content, all rely on the action of molecular motors to package their genomes.

Single-Macromolecule Studies of Eukaryotic Genomic Maintenance

Genomes are self-organized and self-maintained as long, complex macromolecules of chromatin. The inherent heterogeneity, stochasticity, phase separation, and chromatin dynamics of genome operation make it challenging to study genomes using ensemble methods. Various single-molecule force-, fluorescent-, and sequencing-based techniques rooted in different disciplines have been developed to fill critical gaps in the capabilities of bulk measurements, each providing unique, otherwise inaccessible, insights into the structure and maintenance of the genome. Capable of capturing molecular-level details about the organization, conformational changes, and packaging of genetic material, as well as processive and stochastic movements of maintenance factors, a single-molecule toolbox provides an excellent opportunity for collaborative research to understand how genetic material functions in health and malfunctions in disease. In this review, we discuss novel insights brought to genomic sciences by single-molecule techniques and their potential to continue to revolutionize the field-one molecule at a time.

Identification of Eleutherococcus senticosus NAC transcription factors and their mechanisms in mediating DNA methylation of EsFPS, EsSS, and EsSE promoters to regulate saponin synthesis

BACKGROUND

The formation of pharmacologically active components in medicinal plants is significantly impacted by DNA methylation. However, the exact mechanisms through which DNA methylation regulates secondary metabolism remain incompletely understood. Research in model species has demonstrated that DNA methylation at the transcription factor binding site within functional gene promoters can impact the binding of transcription factors to target DNA, subsequently influencing gene expression. These findings suggest that the interaction between transcription factors and target DNA could be a significant mechanism through which DNA methylation regulates secondary metabolism in medicinal plants.

RESULTS

This research conducted a comprehensive analysis of the NAC family in E. senticosus, encompassing genome-wide characterization and functional analysis. A total of 117 EsNAC genes were identified and phylogenetically divided into 15 subfamilies. Tandem duplications and chromosome segment duplications were found to be the primary replication modes of these genes. Motif 2 was identified as the core conserved motif of the genes, and the cis-acting elements, gene structures, and expression patterns of each EsNAC gene were different. EsJUB1, EsNAC047, EsNAC098, and EsNAC005 were significantly associated with the DNA methylation ratio in E. senticosus. These four genes were located in the nucleus or cytoplasm and exhibited transcriptional self-activation activity. DNA methylation in EsFPS, EsSS, and EsSE promoters significantly reduced their activity. The methyl groups added to cytosine directly hindered the binding of the promoters to EsJUB1, EsNAC047, EsNAC098, and EsNAC005 and altered the expression of EsFPS, EsSS, and EsSE genes, eventually leading to changes in saponin synthesis in E. senticosus.

CONCLUSIONS

NAC transcription factors that are hindered from binding by methylated DNA are found in E. senticosus. The incapacity of these NACs to bind to the promoter of the methylated saponin synthase gene leads to subsequent alterations in gene expression and saponin synthesis. This research is the initial evidence showcasing the involvement of EsNAC in governing the impact of DNA methylation on saponin production in E. senticosus.

G-quadruplex modulation by E. coli SSB: A comprehensive study on binding affinities and modes using single-molecule FRET

G-quadruplexes (GQs) are essential guanine-rich secondary structures found in DNA and RNA, playing crucial roles in genomic maintenance and stability. Recent studies have unveiled GQs in the intergenic regions of the E. coli genome, suggesting their biological significance and potential as anti-microbial targets. Here, we investigated the interaction between homo-tetrameric E. coli SSB and GQ-forming single-stranded DNA (ssDNA) sequence with varying lengths. Combining Microscale Thermophoresis (MST) and conventional spectroscopic techniques, we explored E. coli SSB binding to ssDNA and the structural changes of these secondary DNA structures upon protein binding. Subsequently, we have utilized smFRET to probe the conformational changes of GQ-ssDNA structures upon SSB binding. Our results provide detailed insights into SSB's access to various GQ-ssDNA sequencies and the wrapping of this homo-tetrameric protein around GQ-ssDNA in multiple distinct binding modalities. This study sheds light on the intricate details of E. coli SSB's interaction with ssDNA and the resulting widespread conformational changes within these oligonucleotide structures after protein binding. It offers a thorough insight into SSB's accesses to various GQ-ssDNA architectures. The finding demonstrates the multifaceted binding methods through which this homo-tetrameric protein envelops GQ-ssDNA and could prove valuable in deciphering biological processes that involve DNA G-quadruplexes.

Reversible strain-promoted DNA polymerization

Molecular strain can be introduced to influence the outcome of chemical reactions. Once a thermodynamic product is formed, however, reversing the course of a strain-promoted reaction is challenging. Here, a reversible, strain-promoted polymerization in cyclic DNA is reported. The use of nonhybridizing, single-stranded spacers as short as a single nucleotide in length can promote DNA cyclization. Molecular strain is generated by duplexing the spacers, leading to ring opening and subsequent polymerization. Then, removal of the strain-generating duplexers triggers depolymerization and cyclic dimer recovery via enthalpy-driven cyclization and entropy-mediated ring contraction. This reversibility is retained even when a protein is conjugated to the DNA strands, and the architecture of the protein assemblies can be modulated between bivalent and polyvalent states. This work underscores the utility of using DNA not only as a programmable ligand for assembly but also as a route to access restorable bonds, thus providing a molecular basis for DNA-based materials with shape-memory, self-healing, and stimuli-responsive properties.

Dependence of nucleosome mechanical stability on DNA mismatches

The organization of nucleosomes into chromatin and their accessibility are shaped by local DNA mechanics. Conversely, nucleosome positions shape genetic variations, which may originate from mismatches during replication and chemical modification of DNA. To investigate how DNA mismatches affect the mechanical stability and the exposure of nucleosomal DNA, we used an optical trap combined with single-molecule FRET and a single-molecule FRET cyclization assay. We found that a single base-pair C-C mismatch enhances DNA bendability and nucleosome mechanical stability for the 601-nucleosome positioning sequence. An increase in force required for DNA unwrapping from the histone core is observed for single base-pair C-C mismatches placed at three tested positions: at the inner turn, at the outer turn, or at the junction of the inner and outer turn of the nucleosome. The results support a model where nucleosomal DNA accessibility is reduced by mismatches, potentially explaining the preferred accumulation of single-nucleotide substitutions in the nucleosome core and serving as the source of genetic variation during evolution and cancer progression. Mechanical stability of an intact nucleosome, that is mismatch-free, is also dependent on the species as we find that yeast nucleosomes are mechanically less stable and more symmetrical in the outer turn unwrapping compared to Xenopus nucleosomes.

Sequential requirements for distinct Polθ domains during theta-mediated end joining

DNA polymerase θ (Polθ) plays a central role in a DNA double-strand break repair pathway termed theta-mediated end joining (TMEJ). TMEJ functions by pairing short-sequence "microhomologies" (MHs) in single-stranded DNA at each end of a break and subsequently initiating DNA synthesis. It is not known how the Polθ helicase domain (HD) and polymerase domain (PD) operate to bring together MHs and facilitate repair. To resolve these transient processes in real time, we utilized in vitro single-molecule FRET approaches and biochemical analyses. We find that the Polθ-HD mediates the initial capture of two ssDNA strands, bringing them in close proximity. The Polθ-PD binds and stabilizes pre-annealed MHs to form a synaptic complex (SC) and initiate repair synthesis. Individual synthesis reactions show that Polθ is inherently non-processive, accounting for complex mutational patterns during TMEJ. Binding of Polθ-PD to stem-loop-forming sequences can substantially limit synapsis, depending on the available dNTPs and sequence context.

Systematic Comparison of Atomistic Force Fields for the Mechanical Properties of Double-Stranded DNA

The response of double-stranded DNA to external mechanical stress plays a central role in its interactions with the protein machinery in the cell. Modern atomistic force fields have been shown to provide highly accurate predictions for the fine structural features of the duplex. In contrast, and despite their pivotal function, less attention has been devoted to the accuracy of the prediction of the elastic parameters. Several reports have addressed the flexibility of double-stranded DNA via all-atom molecular dynamics, yet the collected information is insufficient to have a clear understanding of the relative performance of the various force fields. In this work, we fill this gap by performing a systematic study in which several systems, characterized by different sequence contexts, are simulated with the most popular force fields within the AMBER family, bcs1 and OL15, as well as with CHARMM36. Analysis of our results, together with their comparison with previous work focused on bsc0, allows us to unveil the differences in the predicted rigidity between the newest force fields and suggests a roadmap to test their performance against experiments. In the case of the stretch modulus, we reconcile these differences, showing that a single mapping between sequence-dependent conformation and elasticity via the crookedness parameter captures simultaneously the results of all force fields, supporting the key role of crookedness in the mechanical response of double-stranded DNA.

Use of DNA forceps to measure receptor-ligand dissociation equilibrium constants in a single-molecule competition assay

The ability of biophysicists to decipher the behavior of individual biomolecules has steadily improved over the past thirty years. However, it still remains unclear how an ensemble of data acquired at the single-molecule level compares with the data acquired on an ensemble of the same molecules. We here propose an assay to tackle this question in the context of dissociation equilibrium constant measurements. A sensor is built by engrafting a receptor and a ligand onto a flexible dsDNA scaffold and mounting this assembly on magnetic tweezers. This way, looking at the position of the magnetic bead enables one to determine in real-time if the two molecular partners are associated or not. Next, to quantify the affinity of the scrutinized single-receptor for a given competitor, various amounts of the latter molecule are introduced in solution and the equilibrium response of the sensor is monitored throughout the titration protocol. Proofs of concept are established for the binding of three rapamycin analogs to the FKBP12 cis-trans prolyl isomerase. For each of these drugs the mean affinity constant obtained on a ten of individual receptors agrees with the one previously determined in a bulk assay. Furthermore, experimental contingencies are sufficient to explain the dispersion observed over the single-molecule values.

Molecular Characterization and Genome Mechanical Features of Two Newly Isolated Polyvalent Bacteriophages Infecting Pseudomonas syringae pv. garcae

Coffee plants have been targeted by a devastating bacterial disease, a condition known as bacterial blight, caused by the phytopathogen Pseudomonas syringae pv. garcae (Psg). Conventional treatments of coffee plantations affected by the disease involve frequent spraying with copper- and kasugamycin-derived compounds, but they are both highly toxic to the environment and stimulate the appearance of bacterial resistance. Herein, we report the molecular characterization and mechanical features of the genome of two newly isolated (putative polyvalent) lytic phages for Psg. The isolated phages belong to class Caudoviricetes and present a myovirus-like morphotype belonging to the genuses Tequatrovirus (PsgM02F) and Phapecoctavirus (PsgM04F) of the subfamilies Straboviridae (PsgM02F) and Stephanstirmvirinae (PsgM04F), according to recent bacterial viruses' taxonomy, based on their complete genome sequences. The 165,282 bp (PsgM02F) and 151,205 bp (PsgM04F) genomes do not feature any lysogenic-related (integrase) genes and, hence, can safely be assumed to follow a lytic lifestyle. While phage PsgM02F produced a morphogenesis yield of 124 virions per host cell, phage PsgM04F produced only 12 virions per host cell, indicating that they replicate well in Psg with a 50 min latency period. Genome mechanical analyses established a relationship between genome bendability and virion morphogenesis yield within infected host cells.

The impact of the sequence-dependent physical properties of DNA on chromatin dynamics

The local mechanical properties of DNA depend on local sequence. Here we review recent genomic, structural, and computational efforts at deciphering the "mechanical code", i.e., the mapping between sequence and mechanics. We then discuss works that suggest how evolution has exploited the mechanical code to control the energetics of DNA-deforming biological processes such as nucleosome organization, transcription factor binding, DNA supercoiling, gene regulation, and 3D chromatin organization. As a whole, these recent works suggest that DNA sequence in diverse organisms can encode regulatory information governing diverse processes via the mechanical code.

Predictions of DNA mechanical properties at a genomic scale reveal potentially new functional roles of DNA flexibility

Mechanical properties of DNA have been implied to influence many of its biological functions. Recently, a new high-throughput method, called loop-seq, which allows measuring the intrinsic bendability of DNA fragments, has been developed. Using loop-seq data, we created a deep learning model to explore the biological significance of local DNA flexibility in a range of different species from different kingdoms. Consistently, we observed a characteristic and largely dinucleotide-composition-driven change of local flexibility near transcription start sites. In the presence of a TATA-box, a pronounced peak of high flexibility can be observed. Furthermore, depending on the transcription factor investigated, flanking-sequence-dependent DNA flexibility was identified as a potential factor influencing DNA binding. Compared to randomized genomic sequences, depending on species and taxa, actual genomic sequences were observed both with increased and lowered flexibility. Furthermore, in Arabidopsis thaliana, mutation rates, both de novo and fixed, were found to be associated with relatively rigid sequence regions. Our study presents a range of significant correlations between characteristic DNA mechanical properties and genomic features, the significance of which with regard to detailed molecular relevance awaits further theoretical and experimental exploration.

Temperature-induced DNA density transition in phage λ capsid revealed with contrast-matching SANS

Structural details of a genome packaged in a viral capsid are essential for understanding how the structural arrangement of a viral genome in a capsid controls its release dynamics during infection, which critically affects viral replication. We previously found a temperature-induced, solid-like to fluid-like mechanical transition of packaged λ-genome that leads to rapid DNA ejection. However, an understanding of the structural origin of this transition was lacking. Here, we use small-angle neutron scattering (SANS) to reveal the scattering form factor of dsDNA packaged in phage λ capsid by contrast matching the scattering signal from the viral capsid with deuterated buffer. We used small-angle X-ray scattering and cryoelectron microscopy reconstructions to determine the initial structural input parameters for intracapsid DNA, which allows accurate modeling of our SANS data. As result, we show a temperature-dependent density transition of intracapsid DNA occurring between two coexisting phases-a hexagonally ordered high-density DNA phase in the capsid periphery and a low-density, less-ordered DNA phase in the core. As the temperature is increased from 20 °C to 40 °C, we found that the core-DNA phase undergoes a density and volume transition close to the physiological temperature of infection (~37 °C). The transition yields a lower energy state of DNA in the capsid core due to lower density and reduced packing defects. This increases DNA mobility, which is required to initiate rapid genome ejection from the virus capsid into a host cell, causing infection. These data reconcile our earlier findings of mechanical DNA transition in phage.

Twist-stretch relations in nucleic acids

Nucleic acids are highly deformable helical molecules constantly stretched, twisted and bent in their biological functioning. Single molecule experiments have shown that double stranded (ds)-RNA and standard ds-DNA have opposite twist-stretch patterns and stretching properties when overwound under a constant applied load. The key structural features of the A-form RNA and B-form DNA helices are here incorporated in a three-dimensional mesoscopic Hamiltonian model which accounts for the radial, bending and twisting fluctuations of the base pairs. Using path integral techniques which sum over the ensemble of the base pair fluctuations, I compute the average helical repeat of the molecules as a function of the load. The obtained twist-stretch relations and stretching properties, for short A- and B-helical fragments, are consistent with the opposite behaviors observed in kilo-base long molecules.

A method for assaying DNA flexibility

The transcription, replication, packaging, and repair of genetic information ubiquitously involves DNA:protein interactions and other biological processes that require local mechanical distortions of DNA. The energetics of such DNA-deforming processes are thus dependent on the local mechanical properties of DNA such as bendability or torsional rigidity. Such properties, in turn, depend on sequence, making it possible for sequence to regulate diverse biological processes by controlling the local mechanical properties of DNA. A deeper understanding of how such a "mechanical code" can encode broad regulatory information has historically been hampered by the absence of technology to measure in high throughput how local DNA mechanics varies with sequence along large regions of the genome. This was overcome in a recently developed technique called loop-seq. Here we describe a variant of the loop-seq protocol, that permits making rapid flexibility measurements in low-throughput, without the need for next-generation sequencing. We use our method to validate a previous prediction about how the binding site for the bacterial transcription factor Integration Host Factor (IHF) might serve as a rigid roadblock, preventing efficient enhancer-promoter contacts in IHF site containing promoters in E. coli, which can be relieved by IHF binding.

Examining DNA breathing with pyDNA-EPBD

MOTIVATION

The two strands of the DNA double helix locally and spontaneously separate and recombine in living cells due to the inherent thermal DNA motion. This dynamics results in transient openings in the double helix and is referred to as "DNA breathing" or "DNA bubbles." The propensity to form local transient openings is important in a wide range of biological processes, such as transcription, replication, and transcription factors binding. However, the modeling and computer simulation of these phenomena, have remained a challenge due to the complex interplay of numerous factors, such as, temperature, salt content, DNA sequence, hydrogen bonding, base stacking, and others.

RESULTS

We present pyDNA-EPBD, a parallel software implementation of the Extended Peyrard-Bishop-Dauxois (EPBD) nonlinear DNA model that allows us to describe some features of DNA dynamics in detail. The pyDNA-EPBD generates genomic scale profiles of average base-pair openings, base flipping probability, DNA bubble probability, and calculations of the characteristically dynamic length indicating the number of base pairs statistically significantly affected by a single point mutation using the Markov Chain Monte Carlo algorithm.

AVAILABILITY AND IMPLEMENTATION

pyDNA-EPBD is supported across most operating systems and is freely available at https://github.com/lanl/pyDNA_EPBD. Extensive documentation can be found at https://lanl.github.io/pyDNA_EPBD/.

Assessing base-resolution DNA mechanics on the genome scale

Intrinsic DNA properties including bending play a crucial role in diverse biological systems. A recent advance in a high-throughput technology called loop-seq makes it possible to determine the bendability of hundred thousand 50-bp DNA duplexes in one experiment. However, it's still challenging to assess base-resolution sequence bendability in large genomes such as human, which requires thousands of such experiments. Here, we introduce 'BendNet'-a deep neural network to predict the intrinsic DNA bending at base-resolution by using loop-seq results in yeast as training data. BendNet can predict the DNA bendability of any given sequence from different species with high accuracy. To explore the utility of BendNet, we applied it to the human genome and observed DNA bendability is associated with chromatin features and disease risk regions involving transcription/enhancer regulation, DNA replication, transcription factor binding and extrachromosomal circular DNA generation. These findings expand our understanding on DNA mechanics and its association with transcription regulation in mammals. Lastly, we built a comprehensive resource of genomic DNA bendability profiles for 307 species by applying BendNet, and provided an online tool to assess the bendability of user-specified DNA sequences (http://www.dnabendnet.com/).

Mechanical codes of chemical-scale specificity in DNA motifs

In gene transcription, certain sequences of double-stranded (ds)DNA play a vital role in nucleosome positioning and expression initiation. That dsDNA is deformed to various extents in these processes leads us to ask: Could the genomic DNA also have sequence specificity in its chemical-scale mechanical properties? We approach this question using statistical machine learning to determine the rigidity between DNA chemical moieties. What emerges for the polyA, polyG, TpA, and CpG sequences studied here is a unique trigram that contains the quantitative mechanical strengths between bases and along the backbone. In a way, such a sequence-dependent trigram could be viewed as a DNA mechanical code. Interestingly, we discover a compensatory competition between the axial base-stacking interaction and the transverse base-pairing interaction, and such a reciprocal relationship constitutes the most discriminating feature of the mechanical code. Our results also provide chemical-scale understanding for experimental observables. For example, the long polyA persistence length is shown to have strong base stacking while its complement (polyAc) exhibits high backbone rigidity. The mechanical code concept enables a direct reading of the physical interactions encoded in the sequence which, with further development, is expected to shed new light on DNA allostery and DNA-binding drugs.

Examining DNA Breathing with pyDNA-EPBD

Motivation

The two strands of the DNA double helix locally and spontaneously separate and recombine in living cells due to the inherent thermal DNA motion.This dynamics results in transient openings in the double helix and is referred to as "DNA breathing" or "DNA bubbles." The propensity to form local transient openings is important in a wide range of biological processes, such as transcription, replication, and transcription factors binding. However, the modeling and computer simulation of these phenomena, have remained a challenge due to the complex interplay of numerous factors, such as, temperature, salt content, DNA sequence, hydrogen bonding, base stacking, and others.

Results

We present pyDNA-EPBD, a parallel software implementation of the Extended Peyrard-Bishop- Dauxois (EPBD) nonlinear DNA model that allows us to describe some features of DNA dynamics in detail. The pyDNA-EPBD generates genomic scale profiles of average base-pair openings, base flipping probability, DNA bubble probability, and calculations of the characteristically dynamic length indicating the number of base pairs statistically significantly affected by a single point mutation using the Markov Chain Monte Carlo (MCMC) algorithm.

Single cohesin molecules generate force by two distinct mechanisms

Spatial organization of DNA is facilitated by cohesin protein complexes that move on DNA and extrude DNA loops. How cohesin works mechanistically as a molecular machine is poorly understood. Here, we measure mechanical forces generated by conformational changes in single cohesin molecules. We show that bending of SMC coiled coils is driven by random thermal fluctuations leading to a ~32 nm head-hinge displacement that resists forces up to 1 pN; ATPase head engagement occurs in a single step of ~10 nm and is driven by an ATP dependent head-head movement, resisting forces up to 15 pN. Our molecular dynamic simulations show that the energy of head engagement can be stored in a mechanically strained conformation of NIPBL and released during disengagement. These findings reveal how single cohesin molecules generate force by two distinct mechanisms. We present a model, which proposes how this ability may power different aspects of cohesin-DNA interaction.

Reversible Regulation of Long-Distance Charge Transport in DNA Nanowires by Dynamically Controlling Steric Conformation

Understanding of DNA-mediated charge transport (CT) is significant for exploring circuits at the molecular scale. However, the fabrication of robust DNA wires remains challenging due to the persistence length and natural flexibility of DNA molecules. Moreover, CT regulation in DNA wires often relies on predesigned sequences, which limit their application and scalability. Here, we addressed these issues by preparing self-assembled DNA nanowires with lengths of 30-120 nm using structural DNA nanotechnology. We employed these nanowires to plug individual gold nanoparticles into a circuit and measured the transport current in nanowires with an optical imaging technique. Contrary to the reported cases with shallow or no length dependence, a fair current attenuation was observed with increasing nanowire length, which experimentally confirmed the prediction of the incoherent hopping model. We also reported a mechanism for the reversible CT regulation in DNA nanowires, which involves dynamic transitions in the steric conformation.

Temperature dependence of DNA elasticity: An all-atom molecular dynamics simulation study

We used all-atom molecular dynamics simulation to investigate the elastic properties of double-stranded DNA (dsDNA). We focused on the influences of temperature on the stretch, bend, and twist elasticities, as well as the twist-stretch coupling, of the dsDNA over a wide range of temperature. The results showed that the bending and twist persistence lengths, together with the stretch and twist moduli, decrease linearly with temperature. However, the twist-stretch coupling behaves in a positive correction and enhances as the temperature increases. The potential mechanisms of how temperature affects dsDNA elasticity and coupling were investigated by using the trajectories from atomistic simulation, in which thermal fluctuations in structural parameters were analyzed in detail. We analyzed the simulation results by comparing them with previous simulation and experimental data, which are in good agreement. The prediction about the temperature dependence of dsDNA elastic properties provides a deeper understanding of DNA elasticities in biological environments and potentially helps in the further development of DNA nanotechnology.

Optical Tweezers to Force Information out of Biological and Synthetic Systems One Molecule at a Time

Over the last few decades, in vitro single-molecule manipulation techniques have enabled the use of force and displacement as controlled variables in biochemistry. Measuring the effect of mechanical force on the real-time kinetics of a biological process gives us access to the rates, equilibrium constants and free-energy landscapes of the mechanical steps of the reaction; this information is not accessible by ensemble assays. Optical tweezers are the current method of choice in single-molecule manipulation due to their versatility, high force and spatial and temporal resolutions. The aim of this review is to describe the contributions of our lab in the single-molecule manipulation field. We present here several optical tweezers assays refined in our laboratory to probe the dynamics and mechano-chemical properties of biological molecular motors and synthetic molecular devices at the single-molecule level.

DNA mechanical flexibility controls DNA potential to activate cGAS-mediated immune surveillance

DNA is well-documented to stimulate immune response. However, the nature of the DNA to activate immune surveillance is less understood. Here, we show that the activation of cyclic GMP-AMP synthase (cGAS) depends on DNA mechanical flexibility, which is controlled by DNA-sequence, -damage and -length. Consistently, DNA-sequence was shown to control cGAS activation. Structural analyses revealed that a conserved cGAS residue (mouse R222 or human R236) contributed to the DNA-flexibility detection. And the residue substitution neutralised the flexibility-controlled DNA-potential to activate cGAS, and relaxed the DNA-length specificity of cGAS. Moreover, low dose radiation was shown to mount cGAS-mediated acute immune surveillance (AIS) via repairable (reusable) DNAs in hrs. Loss of cGAS-mediated AIS decreased the regression of local and abscopal tumours in the context of focal radiation and immune checkpoint blockade. Our results build a direct link between immunosurveillance and DNA mechanical feature.

Mechanical Constraint Effect on DNA Persistence Length

Persistence length is a significant criterion to characterize the semi-flexibility of DNA molecules. The mechanical constraints applied on DNA chains in new single-molecule experiments play a complex role in measuring DNA persistence length; however, there is a difficulty in quantitatively characterizing the mechanical constraint effects due to their complex interactions with electrostatic repulsions and thermal fluctuations. In this work, the classical buckling theory of Euler beam and Manning's statistical theories of electrostatic force and thermal fluctuation force are combined for an isolated DNA fragment to formulate a quantitative model, which interprets the relationship between DNA persistence length and critical buckling length. Moreover, this relationship is further applied to identify the mechanical constraints in different DNA experiments by fitting the effective length factors of buckled fragments. Then, the mechanical constraint effects on DNA persistence lengths are explored. A good agreement among the results by theoretical models, previous experiments, and present molecular dynamics simulations demonstrates that the new superposition relationship including three constraint-dependent terms can effectively characterize changes in DNA persistence lengths with environmental conditions, and the strong constraint-environment coupling term dominates the significant changes of persistence lengths; via fitting effective length factors, the weakest mechanical constraints on DNAs in bulk experiments and stronger constraints on DNAs in single-molecule experiments are identified, respectively. Moreover, the consideration of DNA buckling provides a new perspective to examine the bendability of short-length DNA.

"Flexible hinge" dynamics in mismatched DNA revealed by fluorescence correlation spectroscopy

Altered unwinding/bending fluctuations at DNA lesion sites are implicated as plausible mechanisms for damage sensing by DNA-repair proteins. These dynamics are expected to occur on similar timescales as one-dimensional (1D) diffusion of proteins on DNA if effective in stalling these proteins as they scan DNA. We examined the flexibility and dynamics of DNA oligomers containing 3 base pair (bp) mismatched sites specifically recognized in vitro by nucleotide excision repair protein Rad4 (yeast ortholog of mammalian XPC). A previous Forster resonance energy transfer (FRET) study mapped DNA conformational distributions with cytosine analog FRET pair primarily sensitive to DNA twisting/unwinding deformations (Chakraborty et al. Nucleic Acids Res. 46: 1240-1255 (2018)). These studies revealed B-DNA conformations for nonspecific (matched) constructs but significant unwinding for mismatched constructs specifically recognized by Rad4, even in the absence of Rad4. The timescales of these unwinding fluctuations, however, remained elusive. Here, we labeled DNA with Atto550/Atto647N FRET dyes suitable for fluorescence correlation spectroscopy (FCS). With these probes, we detected higher FRET in specific, mismatched DNA compared with matched DNA, reaffirming unwinding/bending deformations in mismatched DNA. FCS unveiled the dynamics of these spontaneous deformations at ~ 300 µs with no fluctuations detected for matched DNA within the ~ 600 ns-10 ms FCS time window. These studies are the first to visualize anomalous unwinding/bending fluctuations in mismatched DNA on timescales that overlap with the < 500 µs "stepping" times of repair proteins on DNA. Such "flexible hinge" dynamics at lesion sites could arrest a diffusing protein to facilitate damage interrogation and recognition.

Non-linear Hamiltonian models for DNA

Nucleic acids' physical properties have been investigated by theoretical methods based both on fully atomistic representations and on coarse-grained models, e.g., the worm-like-chain, taken from polymer physics. In this review article, I discuss an intermediate (mesoscopic) approach and show how to build a three-dimensional Hamiltonian model which accounts for the main interactions responsible for the stability of the helical molecules. While the 3D mesoscopic model yields a sufficiently detailed description of the helix at the level of the base pair, it also allows one to predict the thermodynamical and structural properties of molecules in solution. Relying on the idea that the base pair fluctuations can be conceived as trajectories, I have built over the past years a computational method based on the time-dependent path integral formalism to derive the partition function. While the main features of the method are presented, I focus here in particular on a newly developed statistical method to set the maximum amplitude of the base pair fluctuations, a key parameter of the theory. Some applications to the calculation of DNA flexibility properties are discussed together with the available experimental data.

Pseudomonas-tailed lytic phages: genome mechanical analysis and putative correlation with virion morphogenesis yield

Aim: To unveil a putative correlation between phage genome flexibility and virion morphogenesis yield. Materials & methods: A deeper analysis of the mechanical properties of three Pseudomonas aeruginosa lytic phage genomes was undertaken, together with full genome cyclizability calculations. Results & conclusion: A putative correlation was established among phage genome flexibility, eclipse timeframe and virion particle morphogenesis yield, with a more flexible phage genome leading to a higher burst size and a more rigid phage genome leading to lower burst sizes. The results obtained are highly relevant to understand the influence of the phage genome plasticity on the virion morphogenesis yield inside the infected bacterial host cells and assumes particular relevance in the actual context of bacterial resistance to antibiotics.

Accurate Sequence-Dependent Coarse-Grained Model for Conformational and Elastic Properties of Double-Stranded DNA

We introduce MADna, a sequence-dependent coarse-grained model of double-stranded DNA (dsDNA), where each nucleotide is described by three beads localized at the sugar, at the base moiety, and at the phosphate group, respectively. The sequence dependence is included by considering a step-dependent parametrization of the bonded interactions, which are tuned in order to reproduce the values of key observables obtained from exhaustive atomistic simulations from the literature. The predictions of the model are benchmarked against an independent set of all-atom simulations, showing that it captures with high fidelity the sequence dependence of conformational and elastic features beyond the single step considered in its formulation. A remarkably good agreement with experiments is found for both sequence-averaged and sequence-dependent conformational and elastic features, including the stretching and torsion moduli, the twist–stretch and twist–bend couplings, the persistence length, and the helical pitch. Overall, for the inspected quantities, the model has a precision comparable to atomistic simulations, hence providing a reliable coarse-grained description for the rationalization of single-molecule experiments and the study of cellular processes involving dsDNA. Owing to the simplicity of its formulation, MADna can be straightforwardly included in common simulation engines. Particularly, an implementation of the model in LAMMPS is made available on an online repository to ease its usage within the DNA research community.

Coordinated DNA and histone dynamics drive accurate histone H2A.Z exchange

Nucleosomal histone H2A is exchanged for its variant H2A.Z by the SWR1 chromatin remodeler, but the mechanism and timing of histone exchange remain unclear. Here, we quantify DNA and histone dynamics during histone exchange in real time using a three-color single-molecule FRET assay. We show that SWR1 operates with timed precision to unwrap DNA with large displacement from one face of the nucleosome, remove H2A-H2B from the same face, and rewrap DNA, all within 2.3 s. This productive DNA unwrapping requires full SWR1 activation and differs from unproductive, smaller-scale DNA unwrapping caused by SWR1 binding alone. On an asymmetrically positioned nucleosome, SWR1 intrinsically senses long-linker DNA to preferentially exchange H2A.Z on the distal face as observed in vivo. The displaced H2A-H2B dimer remains briefly associated with the SWR1-nucleosome complex and is dissociated by histone chaperones. These findings reveal how SWR1 coordinates DNA unwrapping with histone dynamics to rapidly and accurately place H2A.Z at physiological sites on chromatin.

Characterizing the free-energy landscapes of DNA origamis

We show how coarse-grained modelling combined with umbrella sampling using distance-based order parameters can be applied to compute the free-energy landscapes associated with mechanical deformations of large DNA nanostructures. We illustrate this approach for the strong bending of DNA nanotubes and the potentially bistable landscape of twisted DNA origami sheets. The homogeneous bending of the DNA nanotubes is well described by the worm-like chain model; for more extreme bending the nanotubes reversibly buckle with the bending deformations localized at one or two "kinks". For a twisted one-layer DNA origami, the twist is coupled to the bending of the sheet giving rise to a free-energy landscape that has two nearly-degenerate minima that have opposite curvatures. By contrast, for a two-layer origami, the increased stiffness with respect to bending leads to a landscape with a single free-energy minimum that has a saddle-like geometry. The ability to compute such landscapes is likely to be particularly useful for DNA mechanotechnology and for understanding stress accumulation during the self-assembly of origamis into higher-order structures.

Single-Molecule Methods for Investigating the Double-Stranded DNA Bendability

The various DNA-protein interactions associated with the expression of genetic information involve double-stranded DNA (dsDNA) bending. Due to the importance of the formation of the dsDNA bending structure, dsDNA bending properties have long been investigated in the biophysics field. Conventionally, DNA bendability is characterized by innate averaging data from bulk experiments. The advent of single-molecule methods, such as atomic force microscopy, optical and magnetic tweezers, tethered particle motion, and single-molecule fluorescence resonance energy transfer measurement, has provided valuable tools to investigate not only the static structures but also the dynamic properties of bent dsDNA. Here, we reviewed the single-molecule methods that have been used for investigating dsDNA bendability and new findings related to dsDNA bending. Single-molecule approaches are promising tools for revealing the unknown properties of dsDNA related to its bending, particularly in cells.

OUP accepted manuscript

DNA mechanical properties play a critical role in every aspect of DNA-dependent biological processes. Recently a high throughput assay named loop-seq has been developed to quantify the intrinsic bendability of a massive number of DNA fragments simultaneously. Using the loop-seq data, we develop a software tool, DNAcycP, based on a deep-learning approach for intrinsic DNA cyclizability prediction. We demonstrate DNAcycP predicts intrinsic DNA cyclizability with high fidelity compared to the experimental data. Using an independent dataset from in vitro selection for enrichment of loopable sequences, we further verified the predicted cyclizability score, termed C-score, can well distinguish DNA fragments with different loopability. We applied DNAcycP to multiple species and compared the C-scores with available high-resolution chemical nucleosome maps. Our analyses showed that both yeast and mouse genomes share a conserved feature of high DNA bendability spanning nucleosome dyads. Additionally, we extended our analysis to transcription factor binding sites and surprisingly found that the cyclizability is substantially elevated at CTCF binding sites in the mouse genome. We further demonstrate this distinct mechanical property is conserved across mammalian species and is inherent to CTCF binding DNA motif.

Sequence-dependent twist-bend coupling in DNA minicircles

Looping of double-stranded DNA molecules with 100-200 base pairs into minicircles, catenanes, and rotaxanes has been suggested as a potential tool for DNA nanotechnologies. However, sharp DNA bending into a minicircle with a diameter of several to ten nanometers occurs with alterations in the DNA helical structure and may lead to defective kink formation that hampers the use of DNA minicircles, catenanes, and rotaxanes in nanoscale DNA applications. Here, we investigated local variations of a helical twist in sharply bent DNA using microsecond-long all-atom molecular dynamics simulations of six different DNA minicircles, focusing on the sequence dependence of the coupling between DNA bending and its helical twist. Twist angles between consecutive base pairs were analyzed at different locations relative to the direction of DNA bending and, among 10 unique dinucleotide steps, we identified four dinucleotide steps with strong twist-bend coupling, the pyrimidine-purine dinucleotide steps of TA/TA, CG/CG, and CA/TG and the purine-purine dinucleotide step of GA/TC. This work suggests the sequence-dependent structural responses of DNA to strong mechanical deformation, providing new molecular-level insights into the structure and stability of sharply bent DNA minicircles for nanoscale applications.

Kinetics of DNA looping by Anabaena sensory rhodopsin transducer (ASRT) by using DNA cyclization assay

DNA cyclization assay together with single-molecule FRET was employed to monitor protein-mediated bending of a short dsDNA (~ 100 bp). This method provides a simple and easy way to monitor the structural change of DNA in real-time without necessitating prior knowledge of the molecular structures for the optimal dye-labeling. This assay was applied to study how Anabaena sensory rhodopsin transducer (ASRT) facilitates loop formation of DNA as a possible mechanism for gene regulation. The ASRT-induced DNA looping was maximized at 50 mM of Na⁺, while Mg²⁺ also played an essential role in the loop formation.

Compression and Stretching of Confined Linear and Ring Polymers by Applying Force

We use Langevin dynamics to study the deformations of linear and ring polymers in different confinements by applying compression and stretching forces on their two sides. Our results show that the compression deformations are the results of an interplay among of polymer rigidity, degree of confinement, and force applied. When the applied force is beyond the threshold required for the buckling transition, the semiflexible chain under the strong confinement firstly buckles; then comes helical deformation. However, under the same force loading, the semiflexible chain under the weaker confinement exhibits buckling instability and shrinks from the folded ends/sides until it becomes three-folded structures. This happens because the strong confinement not only strongly reduces the buckling wavelength, but also increases the critical buckling force threshold. For the weakly confined polymers, in compression process, the flexible linear polymer collapses into condensed states under a small external force, whereas the ring polymer only shows slight shrinkage, due to the excluded volume interactions of two strands in the crowded states. These results are essential for understanding the deformations of the ring biomacromolecules and polymer chains in mechanical compression or driven transport.

DNA sequence and methylation prescribe the inside-out conformational dynamics and bending energetics of DNA minicircles

Eukaryotic genome and methylome encode DNA fragments' propensity to form nucleosome particles. Although the mechanical properties of DNA possibly orchestrate such encoding, the definite link between 'omics' and DNA energetics has remained elusive. Here, we bridge the divide by examining the sequence-dependent energetics of highly bent DNA. Molecular dynamics simulations of 42 intact DNA minicircles reveal that each DNA minicircle undergoes inside-out conformational transitions with the most likely configuration uniquely prescribed by the nucleotide sequence and methylation of DNA. The minicircles' local geometry consists of straight segments connected by sharp bends compressing the DNA's inward-facing major groove. Such an uneven distribution of the bending stress favors minimum free energy configurations that avoid stiff base pair sequences at inward-facing major grooves. Analysis of the minicircles' inside-out free energy landscapes yields a discrete worm-like chain model of bent DNA energetics that accurately account for its nucleotide sequence and methylation. Experimentally measuring the dependence of the DNA looping time on the DNA sequence validates the model. When applied to a nucleosome-like DNA configuration, the model quantitatively reproduces yeast and human genomes' nucleosome occupancy. Further analyses of the genome-wide chromatin structure data suggest that DNA bending energetics is a fundamental determinant of genome architecture.

Loop-seq: A high-throughput technique to measure the mesoscale mechanical properties of DNA

The local mechanical properties of the DNA polymer influence molecular processes in biology that require mechanical deformations of DNA. Lack of suitable high-throughput experimental techniques had precluded measuring how these properties might vary with sequence along the vast lengths of genomes. Here, we present a detailed protocol for a recently developed experimental technique called loop-seq, which measures at least one local mechanical property of DNA-its propensity to cyclize-in genome-scale throughput. Loop-seq has been used to obtain experimentally derived genome-wide maps of a physical property of DNA. Such measurements have revealed that diverse DNA-deforming processes involved in chromatin organization at various genomic loci are regulated by the genetically encoded, sequence-dependent variations in the mechanical properties of DNA.

Supercoiling and looping promote DNA base accessibility and coordination among distant sites

DNA in cells is supercoiled and constrained into loops and this supercoiling and looping influence every aspect of DNA activity. We show here that negative supercoiling transmits mechanical stress along the DNA backbone to disrupt base pairing at specific distant sites. Cooperativity among distant sites localizes certain sequences to superhelical apices. Base pair disruption allows sharp bending at superhelical apices, which facilitates DNA writhing to relieve torsional strain. The coupling of these processes may help prevent extensive denaturation associated with genomic instability. Our results provide a model for how DNA can form short loops, which are required for many essential processes, and how cells may use DNA loops to position nicks to facilitate repair. Furthermore, our results reveal a complex interplay between site-specific disruptions to base pairing and the 3-D conformation of DNA, which influences how genomes are stored, replicated, transcribed, repaired, and many other aspects of DNA activity.

Base pair fluctuations in helical models for nucleic acids

A statistical method is developed to estimate the maximum amplitude of the base pair fluctuations in a three dimensional mesoscopic model for nucleic acids. The base pair thermal vibrations around the helix diameter are viewed as a Brownian motion for a particle embedded in a stable helical structure. The probability to return to the initial position is computed, as a function of time, by integrating over the particle paths consistent with the physical properties of the model potential. The zero time condition for the first-passage probability defines the constraint to select the integral cutoff for various macroscopic helical conformations, obtained by tuning the twist, bending, and slide motion between adjacent base pairs along the molecule stack. Applying the method to a short homogeneous chain at room temperature, we obtain meaningful estimates for the maximum fluctuations in the twist conformation with ∼10.5 base pairs per helix turn, typical of double stranded DNA helices. Untwisting the double helix, the base pair fluctuations broaden and the integral cutoff increases. The cutoff is found to increase also in the presence of a sliding motion, which shortens the helix contour length, a situation peculiar of dsRNA molecules.

A Brownian ratchet model for DNA loop extrusion by the cohesin complex

The cohesin complex topologically encircles DNA to promote sister chromatid cohesion. Alternatively, cohesin extrudes DNA loops, thought to reflect chromatin domain formation. Here, we propose a structure-based model explaining both activities. ATP and DNA binding promote cohesin conformational changes that guide DNA through a kleisin N-gate into a DNA gripping state. Two HEAT-repeat DNA binding modules, associated with cohesin's heads and hinge, are now juxtaposed. Gripping state disassembly, following ATP hydrolysis, triggers unidirectional hinge module movement, which completes topological DNA entry by directing DNA through the ATPase head gate. If head gate passage fails, hinge module motion creates a Brownian ratchet that, instead, drives loop extrusion. Molecular-mechanical simulations of gripping state formation and resolution cycles recapitulate experimentally observed DNA loop extrusion characteristics. Our model extends to asymmetric and symmetric loop extrusion, as well as z-loop formation. Loop extrusion by biased Brownian motion has important implications for chromosomal cohesin function.

A molecular view of DNA flexibility

DNA dynamics can only be understood by taking into account its complex mechanical behavior at different length scales. At the micrometer level, the mechanical properties of single DNA molecules have been well-characterized by polymer models and are commonly quantified by a persistence length of 50 nm (~150 bp). However, at the base pair level (~3.4 Å), the dynamics of DNA involves complex molecular mechanisms that are still being deciphered. Here, we review recent single-molecule experiments and molecular dynamics simulations that are providing novel insights into DNA mechanics from such a molecular perspective. We first discuss recent findings on sequence-dependent DNA mechanical properties, including sequences that resist mechanical stress and sequences that can accommodate strong deformations. We then comment on the intricate effects of cytosine methylation and DNA mismatches on DNA mechanics. Finally, we review recently reported differences in the mechanical properties of DNA and double-stranded RNA, the other double-helical carrier of genetic information. A thorough examination of the recent single-molecule literature permits establishing a set of general 'rules' that reasonably explain the mechanics of nucleic acids at the base pair level. These simple rules offer an improved description of certain biological systems and might serve as valuable guidelines for future design of DNA and RNA nanostructures.

Mechanical properties of DNA and DNA nanostructures: comparison of atomistic, Martini and oxDNA models

The flexibility and stiffness of small DNA molecules play a fundamental role ranging from several biophysical processes to nano-technological applications. Here, we estimate the mechanical properties of short double-stranded DNA (dsDNA) with lengths ranging from 12 base-pairs (bp) to 56 bp, paranemic crossover (PX) DNA and hexagonal DNA nanotubes (DNTs) using two widely used coarse-grained models - Martini and oxDNA. To calculate the persistence length (Lp) and the stretch modulus (γ) of the dsDNA, we incorporate the worm-like chain and elastic rod model, while for the DNTs, we implement our previously developed theoretical framework. We compare and contrast all of the results with previously reported all-atom molecular dynamics (MD) simulations and experimental results. The mechanical properties of dsDNA (Lp ∼ 50 nm, γ ∼ 800-1500 pN), PX DNA (γ ∼ 1600-2000 pN) and DNTs (Lp ∼ 1-10 μm, γ ∼ 6000-8000 pN) estimated using the Martini soft elastic network and oxDNA are in very good agreement with the all-atom MD and experimental values, while the stiff elastic network Martini reproduces values of Lp and γ which are an order of magnitude higher. The high flexibility of small dsDNA is also depicted in our calculations. However, Martini models proved inadequate to capture the salt concentration effects on the mechanical properties with increasing salt molarity. oxDNA captures the salt concentration effect on the small dsDNA mechanics. But it is found to be ineffective for reproducing the salt-dependent mechanical properties of DNTs. Also, unlike Martini, the time evolved PX DNA and DNT structures from the oxDNA models are comparable to the all-atom MD simulated structures. Our findings provide a route to study the mechanical properties of DNA and DNA based nanostructures with increased time and length scales and has a remarkable implication in the context of DNA 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.