Sequence-Dependent Persistence Lengths of DNA

Constructing Stiffness Tunable DNA Hydrogels Based on DNA Modules with Adjustable Rigidity

DNA hydrogel represents a potent material for crafting biological scaffolds, but the toolbox to systematically regulate the mechanical property is still limited. Herein, we have provided a strategy to tune the stiffness of DNA hydrogel through manipulating the rigidity of DNA modules. By introducing building blocks with higher molecular rigidity and proper connecting fashion, DNA hydrogel stiffness could be systematically elevated. These hydrogels showed excellent dynamic properties and biocompatibility, thus exhibiting great potential in three-dimensional (3D) cell culture. This study has offered a systematic method to explore the structure-property relationship, which may contribute to the development of more intelligent and personalized biomedical platforms.

Insights into elastic properties of coarse-grained DNA models: q-stiffness of cgDNA vs cgDNA

Coarse-grained models have emerged as valuable tools to simulate long DNA molecules while maintaining computational efficiency. These models aim at preserving interactions among coarse-grained variables in a manner that mirrors the underlying atomistic description. We explore here a method for testing coarse-grained vs all-atom models using stiffness matrices in Fourier space (q-stiffnesses), which are particularly suited to probe DNA elasticity at different length scales. We focus on a class of coarse-grained rigid base DNA models known as cgDNA and its most recent version, cgDNA+. Our analysis shows that while cgDNA+ closely follows the q-stiffnesses of the all-atom model, the original cgDNA shows some deviations for twist and bending variables, which are rather strong in the q → 0 (long length scale) limit. The consequence is that while both cgDNA and cgDNA+ give a suitable description of local elastic behavior, the former misses some effects that manifest themselves at longer length scales. In particular, cgDNA performs poorly on twist stiffness, with a value much lower than expected for long DNA molecules. Conversely, the all-atom and cgDNA+ twist are strongly length scale dependent: DNA is torsionally soft at a few base pair distances but becomes more rigid at distances of a few dozen base pairs. Our analysis shows that the bending persistence length in all-atom and cgDNA+ is somewhat overestimated.

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.

Temperature-dependent elasticity of DNA, RNA, and hybrid double helices

Nucleic acid double helices in their DNA, RNA, and DNA-RNA hybrid form play a fundamental role in biology and are main building blocks of artificial nanostructures, but how their properties depend on temperature remains poorly understood. Here, we report thermal dependence of dynamic bending persistence length, twist rigidity, stretch modulus, and twist-stretch coupling for DNA, RNA, and hybrid duplexes between 7°C and 47°C. The results are based on all-atom molecular dynamics simulations using different force field parameterizations. We first demonstrate that unrestrained molecular dynamics can reproduce experimentally known mechanical properties of the duplexes at room temperature. Beyond experimentally known features, we also infer the twist rigidity and twist-stretch coupling of the hybrid duplex. As for the temperature dependence, we found that increasing temperature softens all the duplexes with respect to bending, twisting, and stretching. The relative decrease of the stretch moduli is 0.003-0.004/°C, similar for all the duplex variants despite their very different stretching stiffness, whereas RNA twist stiffness decreases by 0.003/°C, and smaller values are found for the other elastic moduli. The twist-stretch couplings are nearly unaffected by temperature. The stretching, bending, and twisting stiffness all include an important entropic component. Relation of our results to the two-state model of DNA flexibility is discussed. Our work provides temperature-dependent elasticity of nucleic acid duplexes at the microsecond scale relevant for initial stages of protein binding.

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.

Programmable all-DNA hydrogels based on rolling circle and multiprimed chain amplification products

Soft, biocompatible, and tunable materials offer biomedical engineers and material scientists programmable matrices for a variety of biomedical applications. In this regard, DNA hydrogels have emerged as highly promising biomaterials that offer programmable self-assembly, superior biocompatibility, and the presence of specific molecular identifiable structures. Many types of DNA hydrogels have been developed, yet the programmability of the DNA building blocks has not been fully exploited, and further efforts must be directed toward understanding how to finely tune their properties in a predictable manner. Herein, we develop physically crosslinked all-DNA hydrogels with tunable morphology and controllable biodegradation, based on rolling circle amplification and multiprimed chain amplification products. Through molecular engineering of the DNA sequences and their nano-/microscale architectures, the precursors self-assemble in a controlled manner to produce soft hydrogels in an efficient, cost-effective, and highly tunable manner. Notably, we develop a novel DNA microladder architecture that serves as a framework for modulating the hydrogel properties, including over an order of magnitude change in pore size and up to 50% change in biodegradation rate. Overall, we demonstrate how the properties of this DNA-based biomaterial can be tuned by modulating the amounts of rigid double-stranded DNA chains compared to flexible single-stranded DNA chains, as well as through the precursor architecture. Ultimately, this work opens new avenues for the development of programmable and biodegradable soft materials in which DNA functions not only as a store of genetic information but also as a versatile polymeric biomaterial and molecularly engineered macroscale scaffold.

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.

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.

Coarse Graining DNA: Symmetry, Nonlocal Elasticity, and Persistence Length

While the behavior of double-stranded DNA at mesoscopic scales is fairly well understood, less is known about its relation to the rich mechanical properties in the base-pair scale, which is crucial, for instance, to understand DNA-protein interactions and the nucleosome diffusion mechanism. Here, by employing the rigid base-pair model, we connect its microscopic parameters to the persistence length. Combined with all-atom molecular dynamic simulations, our scheme identifies relevant couplings between different degrees of freedom at each coarse-graining step. This allows us to clarify how the scale dependence of the elastic moduli is determined in a systematic way encompassing the role of previously unnoticed off-site couplings between deformations with different parity.

A single strand: A simplified approach to DNA origami

Just as a single polypeptide strand can self-fold into a complex 3D structure, a single strand of DNA can self-fold into DNA origami. Most DNA origami structures (i.e., the scaffold-staple and DNA tiling systems) utilize hundreds of short single-stranded DNA. As such, these structures come with challenges inherent to intermolecular construction. Many assembly challenges involving intermolecular interactions can be resolved if the origami structure is constructed from one DNA strand, where folding is not concentration dependent, the folded structure is more resistant to nuclease degradation, and the synthesis can be achieved at an industrial scale at a thousandth of the cost. This review discusses the design principles and considerations employed in single-stranded DNA origami and its potential benefits and drawbacks.

DNA-Mediated Control of Protein Function in Semi-Synthetic Systems

The development of strategies for controlling protein function in a precise and predictable manner has the potential to revolutionize catalysis, diagnostics, and medicine. In this regard, the use of DNA has emerged as a powerful approach for modulating protein activity. The programmable nature of DNA allows for constructing sophisticated architectures wherein proteins can be placed with control over position, orientation, and stoichiometry. This ability is especially useful considering that the properties of proteins can be influenced by their local environment or their proximity to other functional molecules. Here, we chronicle the different strategies that have been developed to interface DNA with proteins in semi-synthetic systems. We further delineate the unique applications unlocked by the unprecedented level of structural control that DNA affords. We end by outlining outstanding challenges in the area and discuss future research directions towards potential solutions.

Conformational Change of Transcription Factors from Search to Specific Binding: A lac Repressor Case Study

In a process known as facilitated diffusion, DNA-binding proteins find their target sites by combining three-dimensional diffusion and one-dimensional scanning of the DNA. Following the trade-off between speed and stability, agile exploration of DNA requires loose binding, whereas, at the DNA target site, the searching protein needs to establish tight interactions with the DNA. To enable both efficient search and stable binding, DNA-binding proteins and DNA often switch conformations upon recognition. Here, we study the one-dimensional diffusion and DNA binding of the dimeric lac repressor (LacI), which was reported to adopt two different conformations when binding different conformations of DNA. Using coarse-grained molecular dynamic simulations, we studied the diffusion and the sequence-specific binding of these conformations of LacI, as well as their truncated or monomeric variants, with two DNA conformations: straight and bent. The simulations were compared to experimental observables. This study supports that linear diffusion along DNA combines tight rotation-coupled groove tracking and rotation-decoupled hopping, where the protein briefly dissociates and reassociates just a few base pairs away. Tight groove tracking is crucial for target-site recognition, while hopping speeds up the overall search process. We investigated the diffusion of different LacI conformations on DNA and show how the flexibility of LacI's hinge regions ensures agility on DNA as well as faithful groove tracking. If the hinge regions instead form α-helices at the protein-DNA interface, tight groove tracking is not possible. On the contrary, the helical hinge region is essential for tight binding to bent, specific DNA, for the formation of the specific complex. Based on our study of different encounter complexes, we argue that the conformational change in LacI and DNA bending are somewhat coupled. Our findings underline the importance of two distinct protein conformations for facilitated diffusion and specific binding, respectively.

Brownian ratchet for directional nanoparticle transport by repetitive stretch-relaxation of DNA

Brownian motion subject to a periodic and asymmetric potential can be biased by external, nonequilibrium fluctuations, leading to directional movement of Brownian particles. Sequence-dependent flexibility variation along double-stranded DNA has been proposed as a tool to develop periodic and asymmetric potentials for DNA binding of cationic nanoparticles with sizes below tens of nanometers. Here, we propose that repetitive stretching and relaxation of a long, double-stranded DNA molecule with periodic flexibility gradient can induce nonequilibrium fluctuations that tune the amplitude of asymmetric potentials for DNA-nanoparticle binding to result in directional transport of nanometer-sized particles along DNA. Realization of the proposed Brownian ratchet was proven by Brownian dynamics simulations of coarse-grained models of a single, long DNA molecule with flexibility variation and a cationic nanoparticle.

Modulating membrane fusion through the design of fusogenic DNA circuits and bilayer composition

Membrane fusion is a ubiquitous phenomenon linked to many biological processes, and represents a crucial step in liposome-based drug delivery strategies. The ability to control, ever more precisely, membrane fusion pathways would thus be highly valuable for next generation nano-medical solutions and, more generally, the design of advanced biomimetic systems such as synthetic cells. In this article, we present fusogenic nanostructures constructed from synthetic DNA which, different from previous solutions, unlock routes for modulating the rate of fusion and making it conditional to the presence of soluble DNA molecules, thus demonstrating how membrane fusion can be controlled through simple DNA-based molecular circuits. We then systematically explore the relationship between lipid-membrane composition, its biophysical properties, and measured fusion efficiency, linking our observations to the stability of transition states in the fusion pathway. Finally, we observe that specific lipid compositions lead to the emergence of complex bilayer architectures in the fusion products, such as nested morphologies, which are accompanied by alterations in biophysical behaviour. Our findings provide multiple, orthogonal strategies to program lipid-membrane fusion, which leverage the design of either the fusogenic DNA constructs or the physico/chemical properties of the membranes, and could thus be valuable in applications where some design parameters are constrained by other factors such as material cost and biocompatibility, as it is often the case in biotechnological applications.

Selective Integrin α5β1 Targeting through Spatially Constrained Multivalent DNA-Based Nanoparticles

Targeting cells specifically based on receptor expression levels remains an area of active research to date. Selective binding of receptors cannot be achieved by increasing the individual binding strength, as this does not account for differing distributions of receptor density across healthy and diseased cells. Engaging receptors above a threshold concentration would be desirable in devising selective diagnostics. Integrins are prime target candidates as they are readily available on the cell surface and have been reported to be overexpressed in diseases. Insights into their spatial organization would therefore be advantageous to design selective targeting agents. Here, we investigated the effect of activation method on integrin α₅β₁ clustering by immunofluorescence and modeled the global neighbor distances with input from an immuno-staining assay and image processing of microscopy images. This data was used to engineer spatially-controlled DNA-scaffolded bivalent ligands, which we used to compare trends in spatial-selective binding observed across HUVEC, CHO and HeLa in resting versus activated conditions in confocal microscopy images. For HUVEC and CHO, the data demonstrated an improved selectivity and localisation of binding for smaller spacings ~7 nm and ~24 nm, in good agreement with the model. A deviation from the mode predictions for HeLa was observed, indicative of a clustered, instead of homogeneous, integrin organization. Our findings demonstrate how low-technology imaging methods can guide the design of spatially controlled ligands to selectively differentiate between cell type and integrin activation state.

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.

Nanoscale structures and mechanics of peptide nucleic acids

Peptide nucleic acids (PNAs) are charge-neutral polyamide oligomers having extremely favorable thermal stability and high affinity to cell membranes when coupled with cationic cell-penetrating peptides (CPPs), as well as the encouraging antisense and antigene activity in cell-free systems. The study of the mechanical properties of short PNA molecules is rare both in experiments and theoretical calculations. Here, we studied the microscopic structures and elastic properties; namely, persistence length, stretch modulus, twist-stretch coupling, and structural crookedness of double-stranded PNA (dsPNA) and their hybrid derivatives using all-atom MD simulation and compared them with those of double-stranded DNA (dsDNA) and double-stranded RNA (dsRNA). The stretch modulus of the dsPNA is found to be ∼160 pN, an order of magnitude lower than that of dsDNA and smaller than dsRNA, respectively. Similarly, the persistence length of dsPNA is found to be ∼35 nm, significantly smaller than those of dsDNA and dsRNA. The PNA-DNA and PNA-RNA hybrid duplexes have elastic properties lying between that of dsPNA and dsDNA/dsRNA. We argue that the neutral backbones of the PNA make it less stiff than dsDNA and dsRNA molecules. Measurement of structural crookedness and principal component analysis additionally support the bending flexibility of dsPNA. Detailed analysis of the helical-rise coupled to helical-twist indicates that the PNA-DNA hybrid over-winds like dsDNA, while PNA-PNA and PNA-RNA unwind like dsRNA upon stretching. Because of the highly flexible nature of PNA, it can bind other biomolecules by adopting a wide range of conformations and is believed to be crucial for future nanobiotechnology research studies.

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.

Predicting the configuration and energy of DNA in a nucleosome by coarse-grain modelling

We present a novel algorithm which uses a coarse-grained model and an energy minimisation procedure to predict the sequence-dependent DNA configuration in a nucleosome together with its energetic cost.

Potential of Mean Force for DNA Wrapping Around a Cationic Nanoparticle

Sharp bending and wrapping of DNA around proteins and nanoparticles (NPs) has been of extensive research interest. Here, we present the potential of mean force (PMF) for wrapping a DNA double helix around a cationic NP using coarse-grained models of a double-stranded DNA and a cationic NP. Starting from a NP wrapped around by DNA, the PMF was calculated along the distance between the center of the NP and one end of the DNA molecule. A relationship between the distance and the extent of DNA wrapping is used to calculate the PMF as a function of DNA wrapping around a NP. In particular, the PMF was compared for two DNA sequences of (AT)₂₅/(AT)₂₅ and (AC)₂₅/(GT)₂₅, for which the persistence lengths are different by ∼10 nm. The simulation results provide solid evidence of the thermodynamic preference for complex formation of a cationic NP with more flexible DNA over the less flexible DNA. Furthermore, we estimated the elastic energy of DNA bending, which was in good order-of-magnitude agreement with the theoretical prediction of elastic rods. This work suggests that the variation of sequence-dependent DNA flexibility can be utilized in DNA nanotechnologies, in which the position and dynamics of NPs are regulated on large-scale DNA structures, or the structural transformation of DNA is triggered by the sequence-dependent binding of NPs.

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.

Highly Permeable DNA Supramolecular Hydrogel Promotes Neurogenesis and Functional Recovery after Completely Transected Spinal Cord Injury

Regeneration after severe spinal cord injury cannot occur naturally in mammals. Transplanting stem cells to the injury site is a highly promising method, but it faces many challenges because it relies heavily on the microenvironment provided by both the lesion site and delivery material. Although mechanical properties, biocompatibility, and biodegradability of delivery materials have been extensively explored, their permeability has rarely been recognized. Here, a DNA hydrogel is designed with extremely high permeability to repair a 2 mm spinal cord gap in Sprague-Dawley rats. The rats recover basic hindlimb function with detectable motor-evoked potentials, and a renascent neural network is formed via the proliferation and differentiation of both implanted and endogenous stem cells. The signal at the lesion area is conveyed by, on average, 15 newly formed synapses. This hydrogel system offers great potential in clinical trials. Further, it should be easily adaptable to other tissue regeneration applications.

Integration host factor bends and bridges DNA in a multiplicity of binding modes with varying specificity

Nucleoid-associated proteins (NAPs) are crucial in organizing prokaryotic DNA and regulating genes. Vital to these activities are complex nucleoprotein structures, however, how these form remains unclear. Integration host factor (IHF) is an Escherichia coli NAP that creates very sharp bends in DNA at sequences relevant to several functions including transcription and recombination, and is also responsible for general DNA compaction when bound non-specifically. We show that IHF–DNA structural multimodality is more elaborate than previously thought, and provide insights into how this drives mechanical switching towards strongly bent DNA. Using single-molecule atomic force microscopy and atomic molecular dynamics simulations we find three binding modes in roughly equal proportions: ‘associated’ (73° of DNA bend), ‘half-wrapped’ (107°) and ‘fully-wrapped’ (147°), only the latter occurring with sequence specificity. We show IHF bridges two DNA double helices through non-specific recognition that gives IHF a stoichiometry greater than one and enables DNA mesh assembly. We observe that IHF-DNA structural multiplicity is driven through non-specific electrostatic interactions that we anticipate to be a general NAP feature for physical organization of chromosomes.

Effect of DNA Flexibility on Complex Formation of a Cationic Nanoparticle with Double-Stranded DNA

We present extensive molecular dynamics simulations of a cationic nanoparticle and a double-stranded DNA molecule to discuss the effect of DNA flexibility on the complex formation of a cationic nanoparticle with double-stranded DNA. Martini coarse-grained models were employed to describe double-stranded DNA molecules with two different flexibilities and cationic nanoparticles with three different electric charges. As the electric charge of a cationic nanoparticle increases, the degree of DNA bending increases, eventually leading to the wrapping of DNA around the nanoparticle at high electric charges. However, a small increase in the persistence length of DNA by 10 nm requires a cationic nanoparticle with a markedly increased electric charge to bend and wrap DNA around. Thus, a more flexible DNA molecule bends and wraps around a cationic nanoparticle with an intermediate electric charge, whereas a less flexible DNA molecule binds to a nanoparticle with the same electric charge without notable bending. This work provides solid evidence that a small difference in DNA flexibility (as small as 10 nm in persistence length) has a substantial influence on the complex formation of DNA with proteins from a biological perspective and suggests that the variation of sequence-dependent DNA flexibility can be utilized in DNA nanotechnology as a new tool to manipulate the structure of DNA molecules mediated by nanoparticle binding.

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.

Smart and Functionalized Development of Nucleic Acid-Based Hydrogels: Assembly Strategies, Recent Advances, and Challenges

Nucleic acid-based hydrogels that integrate intrinsic biological properties of nucleic acids and mechanical behavior of their advanced assemblies are appealing bioanalysis and biomedical studies for the development of new-generation smart biomaterials. It is inseparable from development and incorporation of novel structural and functional units. This review highlights different functional units of nucleic acids, polymers, and novel nanomaterials in the order of structures, properties, and functions, and their assembly strategies for the fabrication of nucleic acid-based hydrogels. Also, recent advances in the design of multifunctional and stimuli-responsive nucleic acid-based hydrogels in bioanalysis and biomedical science are discussed, focusing on the applications of customized hydrogels for emerging directions, including 3D cell cultivation and 3D bioprinting. Finally, the key challenge and future perspectives are outlined.

CRISPR-Cas System: The Powerful Modulator of Accessory Genomes in Prokaryotes

Being frequently exposed to foreign nucleic acids, bacteria and archaea have developed an ingenious adaptive defense system, called CRISPR-Cas. The system is composed of the Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR) array, together with CRISPR (cas)-associated genes. This system consists of a complex machinery that integrates fragments of foreign nucleic acids from viruses and mobile genetic elements (MGEs), into CRISPR arrays. The inserted segments (spacers) are transcribed and then used by cas proteins as guide RNAs for recognition and inactivation of the targets. Different types and families of CRISPR-Cas systems consist of distinct adaptation and effector modules with evolutionary trajectories, partially independent. The origin of the effector modules and the mechanism of spacer integration/deletion is far less clear. A review of the most recent data regarding the structure, ecology, and evolution of CRISPR-Cas systems and their role in the modulation of accessory genomes in prokaryotes is proposed in this article. The CRISPR-Cas system's impact on the physiology and ecology of prokaryotes, modulation of horizontal gene transfer events, is also discussed here. This system gained popularity after it was proposed as a tool for plant and animal embryo editing, in cancer therapy, as antimicrobial against pathogenic bacteria, and even for combating the novel coronavirus - SARS-CoV-2; thus, the newest and promising applications are reviewed as well.

DNA looping by protamine follows a nonuniform spatial distribution

DNA looping plays an important role in cells in both regulating and protecting the genome. Often, studies of looping focus on looping by prokaryotic transcription factors like lac repressor or by structural maintenance of chromosomes proteins such as condensin. Here, however, we are interested in a different looping method whereby condensing agents (charge ≥+3) such as protamine proteins neutralize the DNA, causing it to form loops and toroids. We considered two previously proposed mechanisms for DNA looping by protamine. In the first mechanism, protamine stabilizes spontaneous DNA fluctuations, forming randomly distributed loops along the DNA. In the second mechanism, protamine binds and bends the DNA to form a loop, creating a distribution of loops that is biased by protamine binding. To differentiate between these mechanisms, we imaged both spontaneous and protamine-induced loops on short-length (≤1 μm) DNA fragments using atomic force microscopy. We then compared the spatial distribution of the loops to several model distributions. A random looping model, which describes the mechanism of spontaneous DNA folding, fit the distribution of spontaneous loops, but it did not fit the distribution of protamine-induced loops. Specifically, it failed to predict a peak in the spatial distribution of loops at an intermediate location along the DNA. An electrostatic multibinding model, which was created to mimic the bind-and-bend mechanism of protamine, was a better fit of the distribution of protamine-induced loops. In this model, multiple protamines bind to the DNA electrostatically within a particular region along the DNA to coordinate the formation of a loop. We speculate that these findings will impact our understanding of protamine’s in vivo role for looping DNA into toroids and the mechanism of DNA condensation by condensing agents more broadly.

Accurate modeling of DNA conformational flexibility by a multivariate Ising model

The sequence-dependent structure and deformability of DNA play a major role for binding of proteins and regulation of gene expression. So far, most efforts to model DNA flexibility are based on unimodal harmonic stiffness models at base-pair resolution. However, multimodal behavior due to distinct conformational substates also contributes significantly to the conformational flexibility of DNA. Moreover, these local substates are correlated to their nearest-neighbor substates. A description for DNA elasticity which includes both multimodality and nearest-neighbor coupling has remained a challenge, which we solve by combining our multivariate harmonic approximation with an Ising model for the substates. In a series of applications to DNA fluctuations and protein-DNA complexes, we demonstrate substantial improvements over the unimodal stiffness model. Furthermore, our multivariate Ising model reveals a mechanical destabilization for adenine (A)-tracts to undergo nucleosome formation. Our approach offers a wide range of applications to determine sequence-dependent deformation energies of DNA and to investigate indirect readout contributions to protein-DNA recognition.

Length-scale-dependent elasticity in DNA from coarse-grained and all-atom models

We investigate the influence of nonlocal couplings on the torsional and bending elasticities of DNA. Such couplings have been observed in the past by several simulation studies. Here, we use a description of DNA conformations based on the variables tilt, roll, and twist. Our analysis of both coarse-grained (oxDNA) and all-atom models indicates that these share strikingly similar features: there are strong off-site couplings for tilt-tilt and twist-twist, while they are much weaker in the roll-roll case. By developing an analytical framework to estimate bending and torsional persistence lengths in nonlocal DNA models, we show how off-site interactions generate a length-scale-dependent elasticity. Based on the simulation-generated elasticity data, the theory predicts a significant length-scale-dependent effect on torsional fluctuations but only a modest effect on bending fluctuations. These results are in agreement with experiments probing DNA mechanics from single base pair to kilobase pair scales.

Bottom-Up Coarse-Grained Modeling of DNA

Recent advances in methodology enable effective coarse-grained modeling of deoxyribonucleic acid (DNA) based on underlying atomistic force field simulations. The so-called bottom-up coarse-graining practice separates fast and slow dynamic processes in molecular systems by averaging out fast degrees of freedom represented by the underlying fine-grained model. The resulting effective potential of interaction includes the contribution from fast degrees of freedom effectively in the form of potential of mean force. The pair-wise additive potential is usually adopted to construct the coarse-grained Hamiltonian for its efficiency in a computer simulation. In this review, we present a few well-developed bottom-up coarse-graining methods, discussing their application in modeling DNA properties such as DNA flexibility (persistence length), conformation, "melting," and DNA condensation.

A perspective on the molecular simulation of DNA from structural and functional aspects

As genetic material, DNA not only carries genetic information by sequence, but also affects biological functions ranging from base modification to replication, transcription and gene regulation through its structural and dynamic properties and variations. The motion and structural properties of DNA involved in related biological processes are also multi-scale, ranging from single base flipping to local DNA deformation, TF binding, G-quadruplex and i-motif formation, TAD establishment, compartmentalization and even chromosome territory formation, just to name a few. The sequence-dependent physical properties of DNA play vital role in all these events, and thus it is interesting to examine how simple sequence information affects DNA and the formation of the chromatin structure in these different hierarchical orders. Accordingly, molecular simulations can provide atomistic details of interactions and conformational dynamics involved in different biological processes of DNA, including those inaccessible by current experimental methods. In this perspective, which is mainly based on our recent studies, we provide a brief overview of the atomistic simulations on how the hierarchical structure and dynamics of DNA can be influenced by its sequences, base modifications, environmental factors and protein binding in the context of the protein-DNA interactions, gene regulation and structural organization of chromatin. We try to connect the DNA sequence, the hierarchical structures of DNA and gene regulation.

Dynamic Manipulation of DNA-Programmed Crystals Embedded in a Polyelectrolyte Hydrogel

DNA is a powerful tool for programming the three-dimensional organization of nanomaterials, where the specificity of nucleotide base-pairing can enable precise, complex, and dynamically addressable structures like colloidal crystals. However, because these DNA-programmed materials are often only stable in solution, their organization can be easily disrupted by changes to its local environment. Methods to stabilize these materials have been developed, but often come at the expense of altering or permanently fixing the materials' structures, removing many of the benefits of using DNA interactions to program assembly. Thus, these methods limit the application of DNA-assembled structures as dynamic and programmable material components. Here, a method is presented to resolve these drawbacks for DNA-grafted nanoparticles, also known as Programmable Atom Equivalents (PAEs), by embedding assembled lattices within a hydrogel matrix. The preformed lattices are exposed to polymerizable residues that electrostatically bind to the charged backbone of the DNA ligands and form a continuous, permeating gel network that stabilizes the colloidal crystals upon introduction of a radical initiator. After embedding PAEs in a hydrogel, deformation of the macroscopic matrix results in concomitant deformation of the PAE lattices, allowing superlattice structural changes to be induced by chemical methods (such as changing solute concentration to alter swelling pressure) or by application of mechanical strain. Changes to the structure of the PAE lattices are reversible and repeatable over multiple cycles and can be either isotropic (such as by swelling) or anisotropic (such as by mechanical deformation). This method of embedding nanoparticle crystals inside of a flexible and environmentally responsive hydrogel is therefore a useful tool in extending the utility of PAEs and other micro- and nanostructures assembled with DNA.

Mechanical Flexibility of DNA: A Quintessential Tool for DNA Nanotechnology

The mechanical properties of DNA have enabled it to be a structural and sensory element in many nanotechnology applications. While specific base-pairing interactions and secondary structure formation have been the most widely utilized mechanism in designing DNA nanodevices and biosensors, the intrinsic mechanical rigidity and flexibility are often overlooked. In this article, we will discuss the biochemical and biophysical origin of double-stranded DNA rigidity and how environmental and intrinsic factors such as salt, temperature, sequence, and small molecules influence it. We will then take a critical look at three areas of applications of DNA bending rigidity. First, we will discuss how DNA's bending rigidity has been utilized to create molecular springs that regulate the activities of biomolecules and cellular processes. Second, we will discuss how the nanomechanical response induced by DNA rigidity has been used to create conformational changes as sensors for molecular force, pH, metal ions, small molecules, and protein interactions. Lastly, we will discuss how DNA's rigidity enabled its application in creating DNA-based nanostructures from DNA origami to nanomachines.

Stretching DNA origami: effect of nicks and Holliday junctions on the axial stiffness

The axial stiffness of DNA origami is determined as a function of key nanostructural characteristics. Different constructs of two-helix nanobeams with specified densities of nicks and Holliday junctions are synthesized and stretched by fluid flow. Implementing single particle tracking to extract force–displacement curves enables the measurement of DNA origami stiffness values at the enthalpic elasticity regime, i.e. for forces larger than 15 pN. Comparisons between ligated and nicked helices show that the latter exhibit nearly a two-fold decrease in axial stiffness. Numerical models that treat the DNA helices as elastic rods are used to evaluate the local loss of stiffness at the locations of nicks and Holliday junctions. It is shown that the models reproduce the experimental data accurately, indicating that both of these design characteristics yield a local stiffness two orders of magnitude smaller than the corresponding value of the intact double-helix. This local degradation in turn leads to a macroscopic loss of stiffness that is evaluated numerically for multi-helix DNA bundles.

Coarse-grained modelling of DNA plectoneme pinning in the presence of base-pair mismatches

Damaged or mismatched DNA bases result in the formation of physical defects in double-stranded DNA. In vivo, defects in DNA must be rapidly and efficiently repaired to maintain cellular function and integrity. Defects can also alter the mechanical response of DNA to bending and twisting constraints, both of which are important in defining the mechanics of DNA supercoiling. Here, we use coarse-grained molecular dynamics (MD) simulation and supporting statistical-mechanical theory to study the effect of mismatched base pairs on DNA supercoiling. Our simulations show that plectoneme pinning at the mismatch site is deterministic under conditions of relatively high force (>2 pN) and high salt concentration (>0.5 M NaCl). Under physiologically relevant conditions of lower force (0.3 pN) and lower salt concentration (0.2 M NaCl), we find that plectoneme pinning becomes probabilistic and the pinning probability increases with the mismatch size. These findings are in line with experimental observations. The simulation framework, validated with experimental results and supported by the theoretical predictions, provides a way to study the effect of defects on DNA supercoiling and the dynamics of supercoiling in molecular detail.

Phage G Structure at 6.1 Å Resolution, Condensed DNA, and Host Identity Revision to a Lysinibacillus

Phage G has the largest capsid and genome of any known propagated phage. Many aspects of its structure, assembly, and replication have not been elucidated. Herein, we present the dsDNA-packed and empty phage G capsid at 6.1 and 9 Å resolution, respectively, using cryo-EM for structure determination and mass spectrometry for protein identification. The major capsid protein, gp27, is identified and found to share the HK97-fold universally conserved in all previously solved dsDNA phages. Trimers of the decoration protein, gp26, sit on the 3-fold axes and are thought to enhance the interactions of the hexameric capsomeres of gp27, for other phages encoding decoration proteins. Phage G's decoration protein is longer than what has been reported in other phages, and we suspect the extra interaction surface area helps stabilize the capsid. We identified several additional capsid proteins, including a candidate for the prohead protease responsible for processing gp27. Furthermore, cryo-EM reveals a range of partially full, condensed DNA densities that appear to have no contact with capsid shell. Three analyses confirm that the phage G host is a Lysinibacillus, and not Bacillus megaterium: identity of host proteins in our mass spectrometry analyses, genome sequence of the phage G host, and host range of phage G.

A multi-modal coarse grained model of DNA flexibility mappable to the atomistic level

We present a new coarse grained method for the simulation of duplex DNA. The algorithm uses a generalized multi-harmonic model that can represent any multi-normal distribution of helical parameters, thus avoiding caveats of current mesoscopic models for DNA simulation and representing a breakthrough in the field. The method has been parameterized from accurate parmbsc1 atomistic molecular dynamics simulations of all unique tetranucleotide sequences of DNA embedded in long duplexes and takes advantage of the correlation between helical states and backbone configurations to derive atomistic representations of DNA. The algorithm, which is implemented in a simple web interface and in a standalone package reproduces with high computational efficiency the structural landscape of long segments of DNA untreatable by atomistic molecular dynamics simulations.

Translational nucleosome positioning: A computational study

About three-quarters of eukaryotic DNA is wrapped into nucleosomes; DNA spools with a protein core. The affinity of a given DNA stretch to be incorporated into a nucleosome is known to depend on the base-pair sequence-dependent geometry and elasticity of the DNA double helix. This causes the rotational and translational positioning of nucleosomes. In this study we ask the question whether the latter can be predicted by a simple coarse-grained DNA model with sequence-dependent elasticity, the rigid base-pair model. Whereas this model is known to be rather robust in predicting rotational nucleosome positioning, we show that the translational positioning is a rather subtle effect that is dominated by the guanine-cytosine content dependence of entropy rather than energy. A correct qualitative prediction within the rigid base-pair framework can only be achieved by assuming that DNA elasticity effectively changes on complexation into the nucleosome complex. With that extra assumption we arrive at a model which gives an excellent quantitative agreement to experimental in vitro nucleosome maps, under the additional assumption that nucleosomes equilibrate their positions only locally.

DNA Nunchucks: Nanoinstrumentation for Single-Molecule Measurement of Stiffness and Bending

Bending of double-stranded DNA (dsDNA) has important applications in biology and engineering, but measurement of DNA bend angles is notoriously difficult and rarely dynamic. Here we introduce a nanoscale instrument that makes dynamic measurement of the bend in short dsDNAs easy enough to be routine. The instrument works by embedding the ends of a dsDNA in stiff, fluorescently labeled DNA nanotubes, thereby mechanically magnifying their orientations. The DNA nanotubes are readily confined to a plane and imaged while freely diffusing. Single-molecule bend angles are rapidly and reliably extracted from the images by a neural network. We find that angular variance across a population increases with dsDNA length, as predicted by the worm-like chain model, although individual distributions can differ significantly from one another. For dsDNAs with phased A₆-tracts, we measure an intrinsic bend of 17 ± 1° per A₆-tract, consistent with other methods, and a length-dependent angular variance that indicates A₆-tracts are (80 ± 30)% stiffer than generic dsDNA.

SerraNA: a program to determine nucleic acids elasticity from simulation data

AT-rich motifs can generate extreme mechanical properties, which are critical for creating strong global bends when phased properly.

Polymer models for the mechanisms of chromatin 3D folding: review and perspective

In this perspective paper, classical physical models for mammalian interphase chromatin folding are reviewed.

The static and dynamic structural heterogeneities of B-DNA: extending Calladine-Dickerson rules

We present a multi-laboratory effort to describe the structural and dynamical properties of duplex B-DNA under physiological conditions. By processing a large amount of atomistic molecular dynamics simulations, we determine the sequence-dependent structural properties of DNA as expressed in the equilibrium distribution of its stochastic dynamics. Our analysis includes a study of first and second moments of the equilibrium distribution, which can be accurately captured by a harmonic model, but with nonlocal sequence-dependence. We characterize the sequence-dependent choreography of backbone and base movements modulating the non-Gaussian or anharmonic effects manifested in the higher moments of the dynamics of the duplex when sampling the equilibrium distribution. Contrary to prior assumptions, such anharmonic deformations are not rare in DNA and can play a significant role in determining DNA conformation within complexes. Polymorphisms in helical geometries are particularly prevalent for certain tetranucleotide sequence contexts and are always coupled to a complex network of coordinated changes in the backbone. The analysis of our simulations, which contain instances of all tetranucleotide sequences, allow us to extend Calladine-Dickerson rules used for decades to interpret the average geometry of DNA, leading to a set of rules with quantitative predictive power that encompass nonlocal sequence-dependence and anharmonic fluctuations.

cgDNAweb: a web interface to the cgDNA sequence-dependent coarse-grain model of double-stranded DNA

The sequence-dependent statistical mechanical properties of fragments of double-stranded DNA is believed to be pertinent to its biological function at length scales from a few base pairs (or bp) to a few hundreds of bp, e.g. indirect read-out protein binding sites, nucleosome positioning sequences, phased A-tracts, etc. In turn, the equilibrium statistical mechanics behaviour of DNA depends upon its ground state configuration, or minimum free energy shape, as well as on its fluctuations as governed by its stiffness (in an appropriate sense). We here present cgDNAweb, which provides browser-based interactive visualization of the sequence-dependent ground states of double-stranded DNA molecules, as predicted by the underlying cgDNA coarse-grain rigid-base model of fragments with arbitrary sequence. The cgDNAweb interface is specifically designed to facilitate comparison between ground state shapes of different sequences. The server is freely available at cgDNAweb.epfl.ch with no login requirement.

Identifying Physical Causes of Apparent Enhanced Cyclization of Short DNA Molecules with a Coarse-Grained Model

DNA cyclization is a powerful technique to gain insight into the nature of DNA bending. While the wormlike chain model provides a good description of small to moderate bending fluctuations, it is expected to break down for large bending. Recent cyclization experiments on strongly bent shorter molecules indeed suggest enhanced flexibility over and above that expected from the wormlike chain. Here, we use a coarse-grained model of DNA to investigate the subtle thermodynamics of DNA cyclization for molecules ranging from 30 to 210 base pairs. As the molecules get shorter, we find increasing deviations between our computed equilibrium j-factor and the classic wormlike chain predictions of Shimada and Yamakawa for a torsionally aligned looped molecule. These deviations are due to sharp kinking, first at nicks, and only subsequently in the body of the duplex. At the shortest lengths, substantial fraying at the ends of duplex domains is the dominant method of relaxation. We also estimate the dynamic j-factor measured in recent FRET experiments. We find that the dynamic j-factor is systematically larger than its equilibrium counterpart—with the deviation larger for shorter molecules—because not all the stress present in the fully cyclized state is present in the transition state. These observations are important for the interpretation of recent cyclization experiments, suggesting that measured anomalously high j-factors may not necessarily indicate non-WLC behavior in the body of duplexes.

TacoxDNA: A user-friendly web server for simulations of complex DNA structures, from single strands to origami

Simulations of nucleic acids at different levels of structural details are increasingly used to complement and interpret experiments in different fields, from biophysics to medicine and materials science. However, the various structural models currently available for DNA and RNA and their accompanying suites of computational tools can be very rarely used in a synergistic fashion. The tacoxDNA webserver and standalone software package presented here are a step toward a long-sought interoperability of nucleic acids models. The webserver offers a simple interface for converting various common input formats of DNA structures and setting up molecular dynamics (MD) simulations. Users can, for instance, design DNA rings with different topologies, such as knots, with and without supercoiling, by simply providing an XYZ coordinate file of the DNA centre-line. More complex DNA geometries, as designable in the cadnano, CanDo and Tiamat tools, can also be converted to all-atom or oxDNA representations, which can then be used to run MD simulations. Though the latter are currently geared toward the native and LAMMPS oxDNA representations, the open-source package is designed to be further expandable. TacoxDNA is available at http://tacoxdna.sissa.it. © 2019 Wiley Periodicals, Inc.

Twist-bend coupling and the statistical mechanics of the twistable wormlike-chain model of DNA: Perturbation theory and beyond

The simplest model of DNA mechanics describes the double helix as a continuous rod with twist and bend elasticity. Recent work has discussed the relevance of a little-studied coupling G between twisting and bending, known to arise from the groove asymmetry of the DNA double helix. Here the effect of G on the statistical mechanics of long DNA molecules subject to applied forces and torques is investigated. We present a perturbative calculation of the effective torsional stiffness C_{eff} for small twist-bend coupling. We find that the "bare" G is "screened" by thermal fluctuations, in the sense that the low-force, long-molecule effective free energy is that of a model with G=0 but with long-wavelength bending and twisting rigidities that are shifted by G-dependent amounts. Using results for torsional and bending rigidities for freely fluctuating DNA, we show how our perturbative results can be extended to a nonperturbative regime. These results are in excellent agreement with numerical calculations for Monte Carlo "triad" and molecular dynamics "oxDNA" models, characterized by different degrees of coarse graining, validating the perturbative and nonperturbative analyses. While our theory is in generally good quantitative agreement with experiment, the predicted torsional stiffness does systematically deviate from experimental data, suggesting that there are as-yet-uncharacterized aspects of DNA twisting-stretching mechanics relevant to low-force, long-molecule mechanical response, which are not captured by widely used coarse-grained models.