The flexibility of locally melted DNA

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.

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.

Influence of a Single Deuterium Substitution for Protium on the Frequency Generation of Different-Size Bubbles in IFNA17

The influence of a single 2H/1H replacement on the frequency generation of different-size bubbles in the human interferon alpha-17 gene (IFNA17) under various energies was studied by a developed algorithm and mathematical modeling without simplifications or averaging. This new approach showed the efficacy of researching DNA bubbles and open states both when all hydrogen bonds in nitrogenous base pairs are protium and after an 2H-substitution. After a single deuterium substitution under specific energies, it was demonstrated that the non-coding region of IFNA17 had a more significant regulatory role in bubble generation in the whole gene than the promoter had. It was revealed that a single deuterium substitution for protium has an influence on the frequency generation of DNA bubbles, which also depends on their size and is always higher for the smaller bubbles under the largest number of the studied energies. Wherein, compared to the natural condition under the same critical value of energy, the bigger raises of the bubble frequency occurrence (maximums) were found for 11-30 base pair (bp) bubbles (higher by 319%), 2-4 bp bubbles (higher by 300%), and 31 bp and over ones (higher by 220%); whereas the most significant reductions of the indicators (minimums) were observed for 11-30 bp bubbles (lower by 43%) and bubbles size over 30 bp (lower by 82%). In this study, we also analyzed the impact of several circumstances on the AT/GC ratio in the formation of DNA bubbles, both under natural conditions and after a single hydrogen isotope exchange. Moreover, based on the obtained data, substantial positive and inverse correlations were revealed between the AT/GC ratio and some factors (energy values, size of DNA bubbles). So, this modeling and variant of the modified algorithm, adapted for researching DNA bubbles, can be useful to study the regulation of replication and transcription in the genes under different isotopic substitutions in the nucleobases.

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.

Kineto-Mechanistic Investigation of Effect of Macromolecular Crowding on the Breathing of DNA Bubble

Macromolecular crowding along with hydrogen bonding or stacking interactions and hydration reportedly has enormous repercussions on elementary biochemical processes, such as the folding of proteins or nucleic acids involving the stability of DNA base pairing. By using the mismatch-induced DNA bubble as a mesoscopic model, the complex interplay of macromolecular crowding on the dynamical fluctuations at the bubble region within the thermodynamic limit has been monitored using single-molecule fluorescence resonance energy transfer (sm-FRET). These single-molecule experimental results have been further corroborated using physical models such as "scaled particle theory" (SPT) and "Gaussian cloud model" (GCM), to predict the biological activity of DNA. The two-state fluctuation of the DNA bubble has been visualized as a function of the nature, size, and concentration of the crowder. The influence of crowders on the DNA conformation has been investigated with the help of the m-factor, the eccentricity, and the kinetic and thermodynamic parameters without any prior assumption. The clear effect of crowding on the dynamics of such a simple biomolecular system emphasizes the power of single-molecule methods and the dependency of the radius of gyration of the co-solute as well as the preferential interaction with the crowder on the distinct conformational states adopted by the bubble. This study provides an idea and hypothesizes the preferential propensity of the DNA bubble to adopt a conformation with the single-stranded domains being far apart, independent of the crowder size, that may be beneficial for efficient recognition by proteins for an uninterrupted procession of the biological process of the central dogma.

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.

The interaction of proanthocyanidins with DNA molecules studied by atomic force microscopy and spectroscopic method

The interactions of the naturally available botanicals with DNA molecules have received considerable attention owing to the potential to develop for medicinal agents. In this study, the interaction of proanthocyanidins with DNA molecules was studied by atomic force microscopy (AFM) and spectroscopic method. The AFM observation indicated that proanthocyanidins induced DNA molecules from double helix chains to the thick rope and toroids. The heights of the formed DNA structures are more than eight times than that of DNA double helix. Spectroscopic measurement results revealed that proanthocyanidins intercalated between the base pairs of DNA in the intercalative binding mode, which resulted in unwinding the DNA helix, twisting the DNA strands and finally transforming into the supercoiled DNA structures. All these results implied that DNA molecule is an important interaction target of proanthocyanidins, and the formed compact DNA structures have biological significance on the gene expression and regulation.

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

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

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.

Improving therapeutic potential of non-viral minimized DNA vectors

The tragic deaths of three patients in a recent AAV-based X-linked myotubular myopathy clinical trial highlight once again the pressing need for safe and reliable gene delivery vectors. Non-viral minimized DNA vectors offer one possible way to meet this need. Recent pre-clinical results with minimized DNA vectors have yielded promising outcomes in cancer therapy, stem cell therapy, stem cell reprograming, and other uses. Broad clinical use of these vectors, however, remains to be realized. Further advances in vector design and production are ongoing. An intriguing and promising potential development results from manipulation of the specific shape of non-viral minimized DNA vectors. By improving cellular uptake and biodistribution specificity, this approach could impact gene therapy, DNA nanotechnology, and personalized medicine.

Softening of DNA near melting as disappearance of an emergent property

Near the melting transition the bending elastic constant κ, an emergent property of double-stranded DNA (dsDNA), is shown not to follow the rodlike scaling for small-length N. The reduction in κ with temperature is determined by the denatured bubbles for a continuous transition, e.g., when the two strands are Gaussian, but by the broken bonds near the open end in a Y-like configuration for a first-order transition as for strands with excluded volume interactions. In the latter case, a lever rule is operational, implying a phase coexistence although dsDNA is known to be a single phase.

Strongly Bent Double-Stranded DNA: Reconciling Theory and Experiment

The strong bending of polymers is poorly understood. We propose a general quantitative framework of polymer bending that includes both the weak and strong bending regimes on the same footing, based on a single general physical principle. As the bending deformation increases beyond a certain (polymer-specific) point, the change in the convexity properties of the effective bending energy of the polymer makes the harmonic deformation energetically unfavorable: in this strong bending regime the energy of the polymer varies linearly with the average bending angle as the system follows the convex hull of the deformation energy function. For double-stranded DNA, the effective bending deformation energy becomes non-convex for bends greater than ~ 2° per base-pair, equivalent to the curvature of a closed circular loop of ~ 160 base pairs. A simple equation is derived for the polymer loop energy that covers both the weak and strong bending regimes. The theory shows quantitative agreement with recent DNA cyclization experiments on short DNA fragments, while maintaining the expected agreement with experiment in the weak bending regime. Counter-intuitively, cyclization probability (j-factor) of very short DNA loops is predicted to increase with decreasing loop length; the j-factor reaches its minimum for loops of ≃ 45 base pairs. Atomistic simulations reveal that the attractive component of the short-range Lennard-Jones interaction between the backbone atoms can explain the underlying non-convexity of the DNA effective bending energy, leading to the linear bending regime. Applicability of the theory to protein-DNA complexes, including the nucleosome, is discussed.

CRISPR-Cas12a Nucleases Bind Flexible DNA Duplexes without RNA/DNA Complementarity

Cas12a (also known as "Cpf1") is a class 2 type V-A CRISPR-associated nuclease that can cleave double-stranded DNA at specific sites. The Cas12a effector enzyme comprises a single protein and a CRISPR-encoded small RNA (crRNA) and has been used for genome editing and manipulation. Work reported here examined in vitro interactions between the Cas12a effector enzyme and DNA duplexes with varying states of base-pairing between the two strands. The data revealed that in the absence of complementarity between the crRNA guide and the DNA target-strand, Cas12a binds duplexes with unpaired segments. These off-target duplexes were bound at the Cas12a site responsible for RNA-guided double-stranded DNA binding but were not cleaved due to the lack of RNA/DNA hybrid formation. Such promiscuous binding was attributed to increased DNA flexibility induced by the unpaired segment present next to the protospacer-adjacent-motif. The results suggest that target discrimination of Cas12a can be influenced by flexibility of the DNA. As such, in addition to the linear sequence, flexibility and other physical properties of the DNA should be considered in Cas12a-based genome engineering applications.

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.

How does temperature impact the conformation of single DNA molecules below melting temperature?

The double stranded DNA molecule undergoes drastic structural changes during biological processes such as transcription during which it opens locally under the action of RNA polymerases. Local spontaneous denaturation could contribute to this mechanism by promoting it. Supporting this idea, different biophysical studies have found an unexpected increase in the flexibility of DNA molecules with various sequences as a function of the temperature, which would be consistent with the formation of a growing number of locally denatured sequences. Here, we take advantage of our capacity to detect subtle changes occurring on DNA by using high throughput tethered particle motion to question the existence of bubbles in double stranded DNA under physiological salt conditions through their conformational impact on DNA molecules ranging from several hundreds to thousands of base pairs. Our results strikingly differ from previously published ones, as we do not detect any unexpected change in DNA flexibility below melting temperature. Instead, we measure a bending modulus that remains stable with temperature as expected for intact double stranded DNA.

Differential Effects of Strand Asymmetry on the Energetics and Structural Flexibility of DNA Internal Loops

Internal loops within structured nucleic acids disrupt local base stacking and destabilize neighboring helical domains; however, these structural motifs also expand the conformational and functional capabilities of structured nucleic acids. Variations in the size, distribution of loop nucleotides on opposing strands (strand asymmetry), and sequence alter their biophysical properties. Here, the thermodynamics and structural flexibility of oligo-T-rich DNA internal loops were systematically investigated in terms of loop size and strand asymmetry. From optical melting experiments, a thermodynamic prediction model is proposed for the energetic penalty of internal loops that accounts for diminishing enthalpic and increasing entropic contributions due to loop size and strand asymmetry for bulges, asymmetric loops, and symmetric loops. These single-stranded domains become less sequence-dependent and more polymeric as the loop size increases. Single-molecule fluorescence resonance energy transfer studies reveal a gradual transition in conformation and structural flexibility from an elongated domain to an increasingly flexible bend that results from increasing strand asymmetry. The findings provide a framework for understanding the thermodynamic and conformational effects of internal loops for the rational design of functional DNA nanostructures.

Revisiting the Anomalous Bending Elasticity of Sharply Bent DNA

Several recent experiments suggest that sharply bent DNA has a surprisingly high bending flexibility, but the cause of this flexibility is poorly understood. Although excitation of flexible defects can explain these results, whether such excitation can occur with the level of DNA bending in these experiments remains unclear. Intriguingly, the DNA contained preexisting nicks in most of these experiments but whether nicks might play a role in flexibility has never been considered in the interpretation of experimental results. Here, using full-atom molecular dynamics simulations, we show that nicks promote DNA basepair disruption at the nicked sites, which drastically reduces DNA bending energy. In addition, lower temperatures suppress the nick-dependent basepair disruption. In the absence of nicks, basepair disruption can also occur but requires a higher level of DNA bending. Therefore, basepair disruption inside B-form DNA can be suppressed if the DNA contains preexisting nicks. Overall, our results suggest that the reported mechanical anomaly of sharply bent DNA is likely dependent on preexisting nicks, therefore the intrinsic mechanisms of sharply bent nick-free DNA remain an open question.

Rigidity of melting DNA

The temperature dependence of DNA flexibility is studied in the presence of stretching and unzipping forces. Two classes of models are considered. In one case the origin of elasticity is entropic due to the polymeric correlations, and in the other the double-stranded DNA is taken to have an intrinsic rigidity for bending. In both cases single strands are completely flexible. The change in the elastic constant for the flexible case due to thermally generated bubbles is obtained exactly. For the case of intrinsic rigidity, the elastic constant is found to be proportional to the square root of the bubble number fluctuation.

Fast Prediction of DNA Melting Bubbles Using DNA Thermodynamic Stability

DNA melting bubbles are the basis of many DNA-protein interactions, such as those in regulatory DNA regions driving gene expression, DNA replication and bacterial horizontal gene transfer. Bubble formation is affected by DNA duplex stability and thermally induced duplex destabilization (TIDD). Although prediction of duplex stability with the nearest neighbor (NN) method is much faster than prediction of TIDD with the Peyrard-Bishop-Dauxois (PBD) model, PBD predicted TIDD defines regulatory DNA regions with higher accuracy and detail. Here, we considered that PBD predicted TIDD is inherently related to the intrinsic duplex stabilities of destabilization regions. We show by regression modeling that NN duplex stabilities can be used to predict TIDD almost as accurately as is predicted with PBD. Predicted TIDD is in fact ascribed to non-linear transformation of NN duplex stabilities in destabilization regions as well as effects of neighboring regions relative to destabilization size. Since the prediction time of our models is over six orders of magnitude shorter than that of PBD, the models present an accessible tool for researchers. TIDD can be predicted on our webserver at http://tidd.immt.eu.

Allostery through protein-induced DNA bubbles

Allostery through DNA is increasingly recognized as an important modulator of DNA functions. Here, we show that the coalescence of protein-induced DNA bubbles can mediate allosteric interactions that drive protein aggregation. We propose that such allostery may regulate DNA's flexibility and the assembly of the transcription machinery. Mitochondrial transcription factor A (TFAM), a dual-function protein involved in mitochondrial DNA (mtDNA) packaging and transcription initiation, is an ideal candidate to test such a hypothesis owing to its ability to locally unwind the double helix. Numerical simulations demonstrate that the coalescence of TFAM-induced bubbles can explain experimentally observed TFAM oligomerization. The resulting melted DNA segment, approximately 10 base pairs long, around the joints of the oligomers act as flexible hinges, which explains the efficiency of TFAM in compacting DNA. Since mitochondrial polymerase (mitoRNAP) is involved in melting the transcription bubble, TFAM may use the same allosteric interaction to both recruit mitoRNAP and initiate transcription.

Atomic force microscopy study of DNA flexibility on short length scales: smooth bending versus kinking

The apparently anomalous flexibility of DNA on short length scales has attracted a lot of attention in recent years. We use atomic force microscopy (AFM) in solution to directly study the DNA bending statistics for small lengths down to one helical turn. The accuracy of experimental estimates could be improved due to a large data volume and a refined algorithm for image processing and measuring bend angles. It is found that, at length scales beyond two helical turns (7 nm), DNA is well described by the harmonic worm-like chain (WLC) model with the bending persistence length of 56 nm. Below this threshold, the AFM data are also described by the WLC model assuming that the accuracy of measured bend angles is limited by the physical width of the double helix. We conclude that the double helical DNA behaves as a uniform elastic rod even at very short length scales. Strong bends due to kinks, melting bubbles and other deviations from the WLC model are statistically negligible.

Effect of temperature on the intrinsic flexibility of DNA and its interaction with architectural proteins

The helical structure of double-stranded DNA is destabilized by increasing temperature. Above a critical temperature (the melting temperature), the two strands in duplex DNA become fully separated. Below this temperature, the structural effects are localized. Using tethered particle motion in a temperature-controlled sample chamber, we systematically investigated the effect of increasing temperature on DNA structure and the interplay between this effect and protein binding. Our measurements revealed that (1) increasing temperature enhances DNA flexibility, effectively leading to more compact folding of the double-stranded DNA chain, and (2) temperature differentially affects different types of DNA-bending chromatin proteins from mesophilic and thermophilic organisms. Thus, our findings aid in understanding genome organization in organisms thriving at moderate as well as extreme temperatures. Moreover, our results underscore the importance of carefully controlling and measuring temperature in single-molecule DNA (micromanipulation) experiments.

Role of microscopic flexibility in tightly curved DNA

The genetic material in living cells is organized into complex structures in which DNA is subjected to substantial contortions. Here we investigate the difference in structure, dynamics, and flexibility between two topological states of a short (107 base pair) DNA sequence in a linear form and a covalently closed, tightly curved circular DNA form. By employing a combination of all-atom molecular dynamics (MD) simulations and elastic rod modeling of DNA, which allows capturing microscopic details while monitoring the global dynamics, we demonstrate that in the highly curved regime the microscopic flexibility of the DNA drastically increases due to the local mobility of the duplex. By analyzing vibrational entropy and Lipari–Szabo NMR order parameters from the simulation data, we propose a novel model for the thermodynamic stability of high-curvature DNA states based on vibrational untightening of the duplex. This novel view of DNA bending provides a fundamental explanation that bridges the gap between classical models of DNA and experimental studies on DNA cyclization, which so far have been in substantial disagreement.

BIOPHYSICAL PROPERTIES OF NUCLEIC ACIDS AT SURFACES RELEVANT TO MICROARRAY PERFORMANCE

Both clinical and analytical metrics produced by microarray-based assay technology have recognized problems in reproducibility, reliability and analytical sensitivity. These issues are often attributed to poor understanding and control of nucleic acid behaviors and properties at solid-liquid interfaces. Nucleic acid hybridization, central to DNA and RNA microarray formats, depends on the properties and behaviors of single strand (ss) nucleic acids (e.g., probe oligomeric DNA) bound to surfaces. ssDNA's persistence length, radius of gyration, electrostatics, conformations on different surfaces and under various assay conditions, its chain flexibility and curvature, charging effects in ionic solutions, and fluorescent labeling all influence its physical chemistry and hybridization under assay conditions. Nucleic acid (e.g., both RNA and DNA) target interactions with immobilized ssDNA strands are highly impacted by these biophysical states. Furthermore, the kinetics, thermodynamics, and enthalpic and entropic contributions to DNA hybridization reflect global probe/target structures and interaction dynamics. Here we review several biophysical issues relevant to oligomeric nucleic acid molecular behaviors at surfaces and their influences on duplex formation that influence microarray assay performance. Correlation of biophysical aspects of single and double-stranded nucleic acids with their complexes in bulk solution is common. Such analysis at surfaces is not commonly reported, despite its importance to microarray assays. We seek to provide further insight into nucleic acid-surface challenges facing microarray diagnostic formats that have hindered their clinical adoption and compromise their research quality and value as genomics tools.

DNA flexibility on short length scales probed by atomic force microscopy

Unusually high bending flexibility has been recently reported for DNA on short length scales. We use atomic force microscopy (AFM) in solution to obtain a direct estimate of DNA bending statistics for scales down to one helical turn. It appears that DNA behaves as a Gaussian chain and is well described by the wormlike chain model at length scales beyond 3 helical turns (10.5 nm). Below this threshold, the AFM data exhibit growing noise because of experimental limitations. This noise may hide small deviations from the Gaussian behavior, but they can hardly be significant.

Poly(dA:dT)-rich DNAs are highly flexible in the context of DNA looping

Large-scale DNA deformation is ubiquitous in transcriptional regulation in prokaryotes and eukaryotes alike. Though much is known about how transcription factors and constellations of binding sites dictate where and how gene regulation will occur, less is known about the role played by the intervening DNA. In this work we explore the effect of sequence flexibility on transcription factor-mediated DNA looping, by drawing on sequences identified in nucleosome formation and ligase-mediated cyclization assays as being especially favorable for or resistant to large deformations. We examine a poly(dA:dT)-rich, nucleosome-repelling sequence that is often thought to belong to a class of highly inflexible DNAs; two strong nucleosome positioning sequences that share a set of particular sequence features common to nucleosome-preferring DNAs; and a CG-rich sequence representative of high G+C-content genomic regions that correlate with high nucleosome occupancy in vivo. To measure the flexibility of these sequences in the context of DNA looping, we combine the in vitro single-molecule tethered particle motion assay, a canonical looping protein, and a statistical mechanical model that allows us to quantitatively relate the looping probability to the looping free energy. We show that, in contrast to the case of nucleosome occupancy, G+C content does not positively correlate with looping probability, and that despite sharing sequence features that are thought to determine nucleosome affinity, the two strong nucleosome positioning sequences behave markedly dissimilarly in the context of looping. Most surprisingly, the poly(dA:dT)-rich DNA that is often characterized as highly inflexible in fact exhibits one of the highest propensities for looping that we have measured. These results argue for a need to revisit our understanding of the mechanical properties of DNA in a way that will provide a basis for understanding DNA deformation over the entire range of biologically relevant scenarios that are impacted by DNA deformability.

Critical bending torque of DNA is a materials parameter independent of local base sequence

Short double-stranded DNA molecules exhibit a softening transition under large bending which is quantitatively described by a critical bending torque τ_{c} at which the molecule develops a kink. Through equilibrium measurements of the elastic energy of short (∼10 nm), highly stressed DNA molecules with a nick at the center we determine τ_{c} for different sequences around the nick. We find that τ_{c} is a robust materials parameter essentially independent of sequence. The measurements also show that, at least for nicked DNA, the local structure at the origin of the softening transition is not a single-stranded "bubble."

Euler buckling and nonlinear kinking of double-stranded DNA

The bending stiffness of double-stranded DNA (dsDNA) at high curvatures is fundamental to its biological activity, yet this regime has been difficult to probe experimentally, and literature results have not been consistent. We created a 'molecular vise' in which base-pairing interactions generated a compressive force on sub-persistence length segments of dsDNA. Short dsDNA strands (<41 base pairs) resisted this force and remained straight; longer strands became bent, a phenomenon called 'Euler buckling'. We monitored the buckling transition via Förster Resonance Energy Transfer (FRET) between appended fluorophores. For low-to-moderate concentrations of monovalent salt (up to ∼150 mM), our results are in quantitative agreement with the worm-like chain (WLC) model of DNA elasticity, without the need to invoke any 'kinked' states. Greater concentrations of monovalent salts or 1 mM Mg(2+) induced an apparent softening of the dsDNA, which was best accounted for by a kink in the region of highest curvature. We tested the effects of all single-nucleotide mismatches on the DNA bending. Remarkably, the propensity to kink correlated with the thermodynamic destabilization of the mismatched DNA relative the perfectly complementary strand, suggesting that the kinked state is locally melted. The molecular vise is exquisitely sensitive to the sequence-dependent linear and nonlinear elastic properties of dsDNA.

Stretching and bending fluctuations of short DNA molecules

Recent measurements of the distribution of end-to-end distance in short DNA molecules infer cooperative stretching fluctuations. The assumptions underlying the analysis can be questioned if transient, thermally induced defects producing a localized decrease in bending stiffness are present in thermal equilibrium, such as regions in which DNA melts locally (bubbles), sustains large-angle bends (kinks), or can locally transform into an alternative (S-DNA) state. We study a generalized discrete worm-like chain model for DNA, capable of describing these experiments, showing that the model yields accurate fits to available experimental data. Our results indicate that DNA bending arising from such localized defects, rather than solely stretching, can be an equal contributor to end-to-end distance fluctuations for 35-bp DNA and contributes nontrivially to such fluctuations at all scales below the persistence length. The analysis suggests that such fluctuations should exhibit a scale-dependent cooperativity, specifically relevant in determining the behavior of short chains, but which saturates rapidly to a length-independent value for longer DNA, to ensure a consistent physical description of DNA across multiple scales. Our approach provides a minimal, yet accurate, coarse-grained description of DNA at the subpersistence length scales of current experimental interest.

Strong bending of the DNA double helix

During the past decade, the issue of strong bending of the double helix has attracted a lot of attention. Here, we overview the major experimental and theoretical developments in the field sorting out reliably established facts from speculations and unsubstantiated claims. Theoretical analysis shows that sharp bends or kinks have to facilitate strong bending of the double helix. It remains to be determined what is the critical curvature of DNA that prompts the appearance of the kinks. Different experimental and computational approaches to the problem are analyzed. We conclude that there is no reliable evidence that any anomalous behavior of the double helix happens when DNA fragments in the range of 100 bp are circularized without torsional stress. The anomaly starts at the fragment length of about 70 bp when sharp bends or kinks emerge in essentially every molecule. Experimental data and theoretical analysis suggest that kinks may represent openings of isolated base pairs, which had been experimentally detected in linear DNA molecules. The calculation suggests that although the probability of these openings in unstressed DNA is close to 10⁻⁵, it increases sharply in small DNA circles reaching 1 open bp per circle of 70 bp.

Polymorphism of DNA conformation inside the bacteriophage capsid

Double-stranded DNA bacteriophage genomes are packaged into their icosahedral capsids at the highest densities known so far (about 50 % w:v). How the molecule is folded at such density and how its conformation changes upon ejection or packaging are fascinating questions still largely open. We review cryo-TEM analyses of DNA conformation inside partially filled capsids as a function of the physico-chemical environment (ions, osmotic pressure, temperature). We show that there exists a wide variety of DNA conformations. Strikingly, the different observed structures can be described by some of the different models proposed over the years for DNA organisation inside bacteriophage capsids: either spool-like structures with axial or concentric symmetries, or liquid crystalline structures characterised by a DNA homogeneous density. The relevance of these conformations for the understanding of DNA folding and unfolding upon ejection and packaging in vivo is discussed.

DNA kinks and bubbles: temperature dependence of the elastic energy of sharply bent 10-nm-size DNA molecules

A 10-nm-long DNA molecule can bend through large angles reversibly. Past the linear regime, its equilibrium nonlinear bending elasticity is governed by a critical bending torque τ(c)≈30pN×nm at which the molecule develops a kink. This nonlinearity has long been attributed to the nucleation of a bubble or melted region in the molecule. Here we measure the temperature dependence of the critical bending torque for nicked DNA, and determine that the entropy associated with the kink in the nonlinear regime is negligible. Thus in the case of nicked DNA the kink is not a bubble, but a compact region deformed beyond a yield strain. We further argue that, with our boundary conditions, the same is likely true for intact DNA. The present measurements confirm that the critical bending torque τ(c) is a materials parameter of DNA mechanics analogous to the bending modulus B≈200pN×nm.

Consequences of local inter-strand dehybridization for large-amplitude bending fluctuations of double-stranded DNA

The wormlike chain model of DNA bending accurately reproduces single-molecule force-extension profiles of long (kilobase) chains. These bending statistics over large scales do not, however, establish a unique microscopic model for elasticity at the 1-10 basepair (bp) scale, which holds particular interest in biological contexts. Here, we examine a class of microscopic models which allow for disruption of base pairing (i.e., a "melt" or "kink", generically an "excitation") and consequently enhanced local flexibility. We first analyze the effect on the excitation free energy of integrating out the spatial degrees of freedom in a wormlike chain. Based on this analysis, we present a formulation of these models that ensures consistency with the well-established thermodynamics of melting in long chains. Using a new method to calculate cyclization statistics of short chains from enhanced-sampling Monte Carlo simulations, we compute J-factors of a meltable wormlike chain over a broad range of chain lengths, including very short molecules (30 bp) that have not yet been explored experimentally. For chains longer than about 120 bp, including most molecules studied to date in the laboratory, we find that melting excitations have little impact on cyclization kinetics. Strong signatures of melting, which might be resolved within typical experimental scatter, emerge only for shorter chains.

DNA Bending through Large Angles Is Aided by Ionic Screening

We used adaptive umbrella sampling on a modified version of the roll angle to simulate the bending of DNA dodecamers. Simulations were carried out with the AMBER and CHARMM force fields for 10 sequences in which the central base pair step was varied. On long length scales, the DNA behavior was found to be consistent with the worm-like chain model. Persistence lengths calculated directly from the simulated structures and indirectly through the use of sequence-independent coarse-grained models based on simulation data were similar to literature values. On short length scales, the free energy cost of bending DNA was found to be consistent with the worm-like chain model for small and intermediate bending angles. At large angles, the bending free energy as a function of the roll angle became linear, suggesting a relative increase in flexibility at larger roll angles. Counterions congregated on the concave side of the highly bent DNA and screened the repulsion of the phosphate groups, facilitating the bending.

Macroscopic loop formation in circular DNA denaturation

The statistical mechanics of DNA denaturation under fixed linking number is qualitatively different from that of unconstrained DNA. Quantitatively different melting scenarios are reached from two alternative assumptions, namely, that the denatured loops are formed at the expense of (i) overtwist or (ii) supercoils. Recent work has shown that the supercoiling mechanism results in a picture similar to Bose-Einstein condensation where a macroscopic loop appears at T{c} and grows steadily with temperature, while the nature of the denatured phase for the overtwisting case has not been studied. By extending an earlier result, we show here that a macroscopic loop appears in the overtwisting scenario as well. We calculate its size as a function of temperature and show that the fraction of the total sum of microscopic loops decreases above T{c}, with a cusp at the critical point.

Enhanced bubble formation in looped short double-stranded DNA

Recent experiments have shown the double-stranded (ds) DNAs readily bend and loop over the scale much shorter than their persistence length (50 nm). In an effort to unveil this seemingly surprising phenomenon, we study the emergence of bubbles in short ds DNA loops by simulating the breathing DNA model. We analyze the bubble size distributions and the melting curves for varying contour lengths, which are critically compared with those of linear DNA of the same lengths. We analytically evaluate the free energies associated with double-strand bending and single-strand bubble formation to explain the simulation data. It is found that in shorter looped DNA the bubbles are more easily initiated and formed to release the large bending energy, giving rise to melting at a lower temperature and a lower contour length.

Temperature dependence of DNA persistence length

We have determined the temperature dependence of DNA persistence length, a, using two different methods. The first approach was based on measuring the j-factors of short DNA fragments at various temperatures. Fitting the measured j-factors by the theoretical equation allowed us to obtain the values of a for temperatures between 5°C and 42°C. The second approach was based on measuring the equilibrium distribution of the linking number between the strands of circular DNA at different temperatures. The major contribution into the distribution variance comes from the fluctuations of DNA writhe in the nicked circular molecules which are specified by the value of a. The computation-based analysis of the measured variances was used to obtain the values of a for temperatures up to 60°C. We found a good agreement between the results obtained by these two methods. Our data show that DNA persistence length strongly depends on temperature and accounting for this dependence is important in quantitative comparison between experimental results obtained at different temperatures.