DNA mechanics and its biological impact

Differential Mg2+ deposition on DNA Holliday Junctions dictates the rate and stability of conformational exchange

DNA Holliday junctions (HJs) are crucial intermediates in genetic recombination and genome repair processes, characterized by a dynamic nature and transitioning among multiple conformations on the timescale ranging from sub-milliseconds to seconds. Although the influence of ions on HJ dynamics has been extensively studied, precise quantification of the thermodynamic feasibility of transitions and detailed kinetic cooperativity remain unexplored. Understanding the heterogeneity of stochastic gene recombination using ensemble-averaged experimental techniques is extremely difficult because of its lack of ability to differentiate dynamics and function in a high spatiotemporal resolution. Herein, we developed a new technique that combines single-molecule fluorescence resonance energy transfer (smFRET) experiments and molecular simulation to investigate the kinetic choreography and preferential stability of HJ conformations under ionic conditions that closely mimic the physiological environment relevant to cellular biology. Our findings predict the prevalence of three distinct conformational macrostates in HJ dynamics. At low ion concentrations, HJs transition rapidly among three thermodynamically stable conformational macrostates. However, in a physiological ionic environment, the open conformation becomes predominant. Using a kinetic network model based on the multi-order time correlation function (TCF), we delineated thermodynamic parameters that govern heterogeneous dynamics as a function of divalent ion concentration. Stabilization of conformations due to an ionic environment and activation barriers concertedly affect transition rates between open and closed conformations. Furthermore, we observed a significant enhancement of Mg2+ condensation in the central region of HJs rather than branch ends, leading to a plausible conclusion that the differential stability of conformational states may be governed by the junction region of HJs rather than duplex branches. This study gives a new insight into the complex interplay between the ionic environment and HJ dynamics, offering a comprehensive understanding of their behavior under conditions relevant to cellular biology and roles in key biological processes for creating a heterogeneous nature of life.

Potential of a newly isolated lytic bacteriophage to control Pseudomonas coronafaciens pv. garcae in coffee plants: Molecular characterization with in vitro and ex vivo experiments

Traditionally, control of coffee plant bacterial halo blight (BHB) caused by the phytopathogen Pseudomonas coronafaciens pv. garcae (Pcg) involves frequent spraying of coffee plantations with non-environmentally friendly and potentially bacterial resistance-promoting copper products or with kasugamycin hydrochloride. In this study we report a leap forward in the quest for a new ecofriendly approach, characterizing (both physicochemically and biologically) and testing both in vitro and ex vivo a new lytic phage for Pcg. An in-depth molecular (genomic and DNA structural features) characterization of the phage was also undertaken. Phage PcgS01F belongs to the class Caudoviricetes, Drexlerviridae family and genus Guelphvirus, and presents a siphovirus-like morphotype. Phage PcgS01F showed a latency period of 40 min and a burst size of 46 PFU/host cell, allowing to conclude that it replicates well in Pcg IBSBF-158. At Multiplicity Of Infection (MOI, or the ratio of phage to bacteria) 1000, the performance of phage PcgS01F was much better than at MOI 10, promoting increasing bacterial reductions until the end of the in vitro inactivation assays, stabilizing at a significant 82 % bacterial load reduction. Phage PcgS01F infected and killed Pcg cells ex vivo in coffee plant leaves artificially contaminated, with a maximum of Pcg inactivation of 7.66 log CFU/mL at MOI 1000 after 36 h of incubation. This study provides evidence that the isolated phage is a promising candidate against the causative agent of BHB in coffee plants.

Fuzzy protein-DNA interactions and beyond: A common theme in transcription?

Gene expression regulation requires both diversity and specificity. How can these two contradictory conditions be reconciled? Dynamic DNA recognition mechanisms lead to heterogeneous bound conformations, which can be shifted by the cellular cues. Here we summarise recent experimental evidence on how fuzzy interactions contribute to chromatin remodelling, regulation of DNA replication and repair and transcription factor binding. We describe how the binding mode continuum between DNA and regulatory factors lead to variable, multisite contact patterns; polyelectrolyte competitions; on-the-fly shape readouts; autoinhibition controlled by posttranslational modifications or dynamic oligomerisation mechanisms. Increasing experimental evidence supports the rugged energy landscape of the bound protein-DNA assembly, modulation of which leads to distinct functional outcomes. Recent results suggest the evolutionary conservation of these combinatorial mechanisms with moderate sequence constraints in the malleable transcriptional machinery.

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.

Pioneer factors: Emerging rules of engagement for transcription factors on chromatinized DNA

Pioneering transcription factors (TFs) can drive cell fate changes by binding their DNA motifs in a repressive chromatin environment. Recent structures illustrate emerging rules for nucleosome engagement: TFs distort the nucleosomal DNA to gain access or employ alternative DNA-binding modes with smaller footprints, they preferentially access solvent-exposed motifs near the entry/exit sites, and frequently interact with histones. The extent of TF-histone interactions, in turn, depends on the motif location on the nucleosome, the type of DNA-binding fold, and adjacent domains present. TF-histone interactions can phase TF motifs relative to nucleosomes, and we discuss how these complex and surprisingly diverse interactions between nucleosomes and TFs contribute to function.

Mechanical response of double-stranded DNA: Bend, twist, and overwind

We employed all-atom molecular dynamics simulations to explore the mechanical response of bending, twisting, and overwinding for double-stranded DNA (dsDNA). We analyzed the bending and twisting deformations, as well as their stiffnesses, using the tilt, roll, and twist modes under stretching force. Findings indicate that the roll and twist angles vary linearly with the stretching force but show opposite trends. The tilt, roll, and twist elastic moduli are considered constants, while the coupling between roll and twist modes slightly decreases under stretching force. The effect of the stretching force on the roll and twist modes, including both their deformations and elasticities, exhibits sequence-dependence, with symmetry around the base pair step. Furthermore, we examined the overwinding path and mechanism of dsDNA from the perspective of the stiffness matrix, based on the tilt, roll, and twist modes. The correlations among tilt, roll, and twist angles imply an alternative overwinding pathway via twist-roll coupling when dsDNA is stretched, wherein entropic contribution prevails.

Taming AID mutator activity in somatic hypermutation

Activation-induced cytidine deaminase (AID) initiates somatic hypermutation (SHM) by introducing base substitutions into antibody genes, a process enabling antibody affinity maturation in immune response. How a mutator is tamed to precisely and safely generate programmed DNA lesions in a physiological process remains unsettled, as its dysregulation drives lymphomagenesis. Recent research has revealed several hidden features of AID-initiated mutagenesis: preferential activity on flexible DNA substrates, restrained activity within chromatin loop domains, unique DNA repair factors to differentially decode AID-caused lesions, and diverse consequences of aberrant deamination. Here, we depict the multifaceted regulation of AID activity with a focus on emerging concepts/factors and discuss their implications for the design of base editors (BEs) that install somatic mutations to correct deleterious genomic variants.

DNASCANNER v2: A Web-Based Tool to Analyze the Characteristic Properties of Nucleotide Sequences

Throughout the process of evolution, DNA undergoes the accumulation of distinct mutations, which can often result in highly organized patterns that serve various essential biological functions. These patterns encompass various genomic elements and provide valuable insights into the regulatory and functional aspects of DNA. The physicochemical, mechanical, thermodynamic, and structural properties of DNA sequences play a crucial role in the formation of specific patterns. These properties contribute to the three-dimensional structure of DNA and influence their interactions with proteins, regulatory elements, and other molecules. In this study, we introduce DNASCANNER v2, an advanced version of our previously published algorithm DNASCANNER for analyzing DNA properties. The current tool is built using the FLASK framework in Python language. Featuring a user-friendly interface tailored for nonspecialized researchers, it offers an extensive analysis of 158 DNA properties, including mono/di/trinucleotide frequencies, structural, physicochemical, thermodynamics, and mechanical properties of DNA sequences. The tool provides downloadable results and offers interactive plots for easy interpretation and comparison between different features. We also demonstrate the utility of DNASCANNER v2 in analyzing splice-site junctions, casposon insertion sequences, and transposon insertion sites (TIS) within the bacterial and human genomes, respectively. We also developed a deep learning module for the prediction of potential TIS in a given nucleotide sequence. In the future, we aim to optimize the performance of this prediction model through extensive training on larger data sets.

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.

Effect of temperature on anisotropic bending elasticity of dsRNA: an all-atom molecular dynamics simulation

Employing all-atom molecular dynamics simulations, we examined the temperature-dependent behavior of bending elasticity in double-stranded RNA (dsRNA). Specifically, we focused on the bending persistence length and its constituent components, namely, the tilt and roll stiffness. Our results revealed a near-linear decrease in these stiffness components as a function of temperature, thereby highlighting the increased flexibility of dsRNA at elevated temperatures. Furthermore, our data revealed a significant anisotropy in dsRNA bending elasticity, which diminished with increasing temperature, attributable to marked disparities in tilt and roll stiffness components. We delineated the underlying biophysical mechanisms and corroborated our findings with extant literature. These observations offer salient implications for advancing our understanding of nucleic acid elasticity, and are pertinent to potential medical applications.

Rapid and ultra-sensitive lateral flow assay for pathogens based on multivalent aptamer and magnetic nanozyme

Ultra-sensitive LFA methods for pathogen detection commonly depended on tedious and time-consuming nucleic acid amplification. Here, a high affinity multivalent aptamer (multi-Apt) for S. aureus was obtained through exquisite engineering design. The scaffold and conformation of the multi-Apt were found to be key factors in the detection signal of aptsensors. After optimization, the binding affinity of the multi-Apt to S. aureus was improved by more than 8-fold from 135.9 nM to 16.77 nM. By the joint use of the multi-Apt and a multifunctional nanozyme Fe3O4@MOF@PtPd, a fast and ultra-sensitive LFA for S. aureus was developed (termed MA-MN LFA). In this method, a Fe3O4@MOF@PtPd nanozyme was modified with vancomycin and could efficiently capture and separate S. aureus. Moreover, the multi-Apt worked together with the nanozyme to bind with S. aureus to form a ternary complex at the same time, which simply the fabrication of LFA strip. The developed MA-MN LFA could detect S. aureus as low as 2 CFU/mL within 30 min and a wide linear range of 10-1 × 108 CFU/mL was obtained. The detection is easily operated, fast (can be completed within 30 min) and versatile for Gram-positive pathogens, thus has great potential as a powerful tool in pathogen detection.

Employing bimodal representations to predict DNA bendability within a self-supervised pre-trained framework

The bendability of genomic DNA, which measures the DNA looping rate, is crucial for numerous biological processes of DNA. Recently, an advanced high-throughput technique known as 'loop-seq' has made it possible to measure the inherent cyclizability of DNA fragments. However, quantifying the bendability of large-scale DNA is costly, laborious, and time-consuming. To close the gap between rapidly evolving large language models and expanding genomic sequence information, and to elucidate the DNA bendability's impact on critical regulatory sequence motifs such as super-enhancers in the human genome, we introduce an innovative computational model, named MIXBend, to forecast the DNA bendability utilizing both nucleotide sequences and physicochemical properties. In MIXBend, a pre-trained language model DNABERT and convolutional neural network with attention mechanism are utilized to construct both sequence- and physicochemical-based extractors for the sophisticated refinement of DNA sequence representations. These bimodal DNA representations are then fed to a k-mer sequence-physicochemistry matching module to minimize the semantic gap between each modality. Lastly, a self-attention fusion layer is employed for the prediction of DNA bendability. In conclusion, the experimental results validate MIXBend's superior performance relative to other state-of-the-art methods. Additionally, MIXBend reveals both novel and known motifs from the yeast. Moreover, MIXBend discovers significant bendability fluctuations within super-enhancer regions and transcription factors binding sites in the human genome.

Influence of temperature on bend, twist and twist-bend coupling of dsDNA

The temperature-dependent bend and twist elasticities of dsDNA, as well as their couplings, were explored through all-atom molecular dynamics simulations. Three rotational parameters, tilt, roll, and twist, were employed to assess the bend and twist elasticities through their stiffness matrix. Our analysis indicates that the bend and twist stiffnesses decrease as the temperature rises, primarily owing to entropic influences stemming from thermodynamic fluctuations. Furthermore, the couplings between these rotational parameters also exhibit a decline with increasing temperature, although the roll-twist coupling displays greater strength than the tilt-roll and tilt-twist couplings, attributed to its more robust correction component. We elucidated the influence of temperature on bend and twist elasticities based on the comparisons between various models and existing data.

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.

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.

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.

Interplay of the mechanical and structural properties of DNA nanostructures determines their electrostatic interactions with lipid membranes

Nucleic acids and lipids function in close proximity in biological processes, as well as in nanoengineered constructs for therapeutic applications. As both molecules carry a rich charge profile, and frequently coexist in complex ionic solutions, the electrostatics surely play a pivotal role in interactions between them. Here we discuss how each component of a DNA/ion/lipid system determines its electrostatic attachment. We examine membrane binding of a library of DNA molecules varying from nanoengineered DNA origami through plasmids to short DNA domains, demonstrating the interplay between the molecular structure of the nucleic acid and the phase of lipid bilayers. Furthermore, the magnitude of DNA/lipid interactions is tuned by varying the concentration of magnesium ions in the physiologically relevant range. Notably, we observe that the structural and mechanical properties of DNA are critical in determining its attachment to lipid bilayers and demonstrate that binding is correlated positively with the size, and negatively with the flexibility of the nucleic acid. The findings are utilized in a proof-of-concept comparison of membrane interactions of two DNA origami designs - potential nanotherapeutic platforms - showing how the results can have a direct impact on the choice of DNA geometry for biotechnological applications.

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.

Modulation of Biofilm Mechanics by DNA Structure and Cell Type

Deoxyribonucleic acid (DNA) evolved as a tool for storing and transmitting genetic information within cells, but outside the cell, DNA can also serve as “construction material” present in microbial biofilms or various body fluids, such as cystic fibrosis, sputum, and pus. In the present work, we investigate the mechanics of biofilms formed from Pseudomonas aeruginosa Xen 5, Staphylococcus aureus Xen 30, and Candida albicans 1408 using oscillatory shear rheometry at different levels of compression and recreate these mechanics in systems of entangled DNA and cells. The results show that the compression-stiffening and shear-softening effects observed in biofilms can be reproduced in DNA networks with the addition of an appropriate number of microbial cells. Additionally, we observe that these effects are cell-type dependent. We also identify other mechanisms that may significantly impact the viscoelastic behavior of biofilms, such as the compression-stiffening effect of DNA cross-linking by bivalent cations (Mg²⁺, Ca²⁺, and Cu²⁺) and the stiffness-increasing interactions of P. aeruginosa Xen 5 biofilm with Pf1 bacteriophage produced by P. aeruginosa. This work extends the knowledge of biofilm mechanobiology and demonstrates the possibility of modifying biopolymers toward obtaining the desired biophysical properties.

Polymer models reveal how chromatin modification can modulate force at the kinetochore

A key feature of chromosome segregation is the ability to sense tension between sister kinetochores. DNA between sister kinetochores must be packaged in a way that sustains tension propagation from one kinetochore to its sister, approximately 1 micron away. A molecular bottlebrush consisting of a primary axis populated with a crowded array of side chains provides a means to build tension over length scales considerably larger than the stiffness of the individual elements, that is, DNA polymer. Evidence for the bottlebrush organization of chromatin between sister kinetochores comes from genetic, cell biological, and polymer modeling of the budding yeast centromere. In this study, we have used polymer dynamic simulations of the bottlebrush to recapitulate experimental observations of kinetochore structure. Several aspects of the spatial distribution of kinetochore proteins and their response to perturbation lack a mechanistic understanding. Changes in physical parameters of bottlebrush, DNA stiffness, and DNA loops directly impact the architecture of the inner kinetochore. This study reveals that the bottlebrush is an active participant in building tension between sister kinetochores and proposes a mechanism for chromatin feedback to the kinetochore.

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.

DNA Bending Force Facilitates Z-DNA Formation under Physiological Salt Conditions

Z-DNA, a noncanonical helical structure of double-stranded DNA (dsDNA), plays pivotal roles in various biological processes, including transcription regulation. Mechanical stresses on dsDNA, such as twisting and stretching, help to form Z-DNA. However, the effect of DNA bending, one of the most common dsDNA deformations, on Z-DNA formation is utterly unknown. Here, we show that DNA bending induces the formation of Z-DNA, that is, more Z-DNA is formed as the bending force becomes stronger. We regulated the bending force on dsDNA by using D-shaped DNA nanostructures. The B-Z transition was observed by single-molecule fluorescence resonance energy transfer. We found that as the bending force became stronger, Z-DNA was formed at lower Mg²⁺ concentrations. When dsDNA contained cytosine methylations, the B-Z transition occurred at 78 mM Mg²⁺ (midpoint) in the absence of the bending force. However, the B-Z transition occurred at a 28-fold lower Mg²⁺ concentration (2.8 mM) in the presence of the bending force. Monte Carlo simulation suggested that the B-Z transition stabilizes the bent form via the formation of the B-Z junction with base extrusion, which effectively releases the bending stress on DNA. Our results clearly show that the bending force facilitates the B-Z transition under physiological salt conditions.

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.

Stretching of Long Double-Stranded DNA and RNA Described by the Same Approach

We propose an approach to help interpret polymer force-extension curves that exhibit plateau regimes. When coupled to a bead-spring dynamic model, the approach accurately reproduces a variety of experimental force-extension curves of long double-stranded DNA and RNA, including torsionally constrained and unconstrained DNA and negatively supercoiled DNA. A key feature of the model is a specific nonconvex energy function of the spring. We provide an algorithm to obtain the five required parameters of the model from experimental force-extension curves. The applicability of the approach to the force-extension curves of double-stranded (ds) DNA of variable GC content as well as to a DNA/RNA hybrid structure is explored and confirmed. We use the approach to explain counterintuitive sequence-dependent trends and make predictions. In the plateau region of the force-extension curves, our molecular dynamics simulations show that the polymer separates into a mix of weakly and strongly stretched states without forming macroscopically distinct phases. The distribution of these states is predicted to depend on the sequence.

Genesis of Open States Zones in a DNA Molecule Depends on the Localization and Value of the Torque

The formation and dynamics of the open states in a double-stranded DNA molecule are largely determined by its mechanical parameters. The main one is the torque. However, the experimental study of DNA dynamics and the occurrence of open states is limited by the spatial resolution of available biophysical instruments. Therefore, in this work, on the basis of a mechanical mathematical model of DNA, calculations of the torque effect on the process of occurrence and dynamics of open states were carried out for the interferon alpha 17 gene. It was shown that torsion action leads to the occurrence of rotational movements of nitrogenous bases. This influence is nonlinear, and an increase in the amplitude of the torsion action does not lead to an automatic increase in the amplitude of rotational movements and an increase in the zones’ open states. Calculations with a constant torsion moment demonstrate that open states zones are more often formed at the boundaries of the gen and in regions with a predominance of A–T pairs. It is shown, that for the occurrence of open states in the part of the gene that contains a small number of A–T pairs, a large amount of torque is required. When the torque is applied to a certain region of the gene, the probability of the formation of the open state depends on the content of A–T pairs in this region, the size of this region, and on the exposure time. For this mathematical model, open states zones can be closed when the torsion action stops. The simulation results showed that the values of the torsion moment required for the appearance of open states zones, in some cases, are close to experimentally measured (13–15 pN·nm).

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.

Influence of Environmental Parameters on the Stability of the DNA Molecule

Fluctuations in viscosity within the cell nucleus have wide limits. When a DNA molecule passes from the region of high viscosity values to the region of low values, open states, denaturation bubbles, and unweaving of DNA strands can occur. Stabilization of the molecule is provided by energy dissipation-dissipation due to interaction with the environment. Separate sections of a DNA molecule in a twisted state can experience supercoiling stress, which, among other things, is due to complex entropic effects caused by interaction with a solvent. In this work, based on the numerical solution of a mechanical mathematical model for the interferon alpha 17 gene and a fragment of the Drosophila gene, an analysis of the external environment viscosity influence on the dynamics of the DNA molecule and its stability was carried out. It has been shown that an increase in viscosity leads to a rapid stabilization of the angular vibrations of nitrogenous bases, while a decrease in viscosity changes the dynamics of DNA: the rate of change in the angular deviations of nitrogenous bases increases and the angular deformations of the DNA strands increase at each moment of time. These processes lead to DNA instability, which increases with time. Thus, the paper considers the influence of the external environment viscosity on the dissipation of the DNA nitrogenous bases' vibrational motion energy. Additionally, the study on the basis of the described model of the molecular dynamics of physiological processes at different indicators of the rheological behavior of nucleoplasm will allow a deeper understanding of the processes of nonequilibrium physics of an active substance in a living cell to be obtained.

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.

More Than Just Letters and Chemistry: Genomics Goes Mechanics

Although ubiquitously thought of as a simple string of letters, DNA exhibits complex physicochemical properties. As a result, DNA can store information beyond the extensively studied explicit genetic message. The mechanical code of DNA has not been studied systematically in a genome-wide context until recent groundbreaking work by Basu et al.