Long-range correlations in the mechanics of small DNA circles under topological stress revealed by multi-scale simulation

Effective Interactions between Double-Stranded DNA Molecules in Aqueous Electrolyte Solutions: Effects of Molecular Architecture and Counterion Valency

A computational investigation of the effects of molecular topology, namely, linear and circular, as well as counterion valency, on the ensuing pairwise effective interactions between DNA molecules in an unlinked state is presented. Umbrella sampling simulations have been performed through the introduction of bias potential along a reaction coordinate defined as the distance between the centers-of-mass of pairs of DNA molecules, and effective pair interaction potentials have been computed by employing the weighted histogram analysis method. An interesting comparison can be drawn between the different DNA topologies studied here, especially with regard to the contrasting effects of divalent counterions on the effective pair potentials: while DNA-DNA repulsion in short center-of-mass distances decreases significantly in the presence of divalent counterion-ions (as compared to monovalent ions) for linear DNA, the opposite effect occurs for the DNA minicircles. This can be attributed to the fact that linear DNA fragments can easily adopt relative orientations that minimize electrostatic and steric repulsions by rotating relative to one another and by exhibiting more pronounced bending due to the presence of free ends.

Tuning Protein Hydrogel Mechanics through Modulation of Nanoscale Unfolding and Entanglement in Postgelation Relaxation

Globular folded proteins are versatile nanoscale building blocks to create biomaterials with mechanical robustness and inherent biological functionality due to their specific and well-defined folded structures. Modulating the nanoscale unfolding of protein building blocks during network formation (in situ protein unfolding) provides potent opportunities to control the protein network structure and mechanics. Here, we control protein unfolding during the formation of hydrogels constructed from chemically cross-linked maltose binding protein using ligand binding and the addition of cosolutes to modulate protein kinetic and thermodynamic stability. Bulk shear rheology characterizes the storage moduli of the bound and unbound protein hydrogels and reveals a correlation between network rigidity, characterized as an increase in the storage modulus, and protein thermodynamic stability. Furthermore, analysis of the network relaxation behavior identifies a crossover from an unfolding dominated regime to an entanglement dominated regime. Control of in situ protein unfolding and entanglement provides an important route to finely tune the architecture, mechanics, and dynamic relaxation of protein hydrogels. Such predictive control will be advantageous for future smart biomaterials for applications which require responsive and dynamic modulation of mechanical properties and biological function.

The dynamic, motile and deformative properties of RNA nanoparticles facilitate the third milestone of drug development

Besides mRNA, rRNA, and tRNA, cells contain many other noncoding RNA that display critical roles in the regulation of cellular functions. Human genome sequencing revealed that the majority of non-protein-coding DNA actually codes for non-coding RNAs. The dynamic nature of RNA results in its motile and deformative behavior. These conformational transitions such as the change of base-pairing, breathing within complemented strands, and pseudoknot formation at the 2D level as well as the induced-fit and conformational capture at the 3D level are important for their biological functions including regulation, translation, and catalysis. The dynamic, motile and catalytic activity has led to a belief that RNA is the origin of life. We have recently reported that the deformative property of RNA nanoparticles enhances their penetration through the leaky blood vessel of cancers which leads to highly efficient tumor accumulation. This special deformative property also enables RNA nanoparticles to pass the glomerulus, overcoming the filtration size limit, resulting in fast renal excretion and rapid body clearance, thus low or no toxicity. The biodistribution of RNA nanoparticles can be further improved by the incorporation of ligands for cancer targeting. In addition to the favorable biodistribution profiles, RNA nanoparticles possess other properties including self-assembly, negative charge, programmability, and multivalency; making it a great material for pharmaceutical applications. The intrinsic negative charge of RNA nanoparticles decreases the toxicity of drugs by preventing nonspecific binding to the negative charged cell membrane and enhancing the solubility of hydrophobic drugs. The polyvalent property of RNA nanoparticles allows the multi-functionalization which can apply to overcome drug resistance. This review focuses on the summary of these unique properties of RNA nanoparticles, which describes the mechanism of RNA dynamic, motile and deformative properties, and elucidates and prepares to welcome the RNA therapeutics as the third milestone in pharmaceutical drug development.

Spatiotemporal Coupling of DNA Supercoiling and Genomic Sequence Organization-A Timing Chain for the Bacterial Growth Cycle?

In this article we describe the bacterial growth cycle as a closed, self-reproducing, or autopoietic circuit, reestablishing the physiological state of stationary cells initially inoculated in the growth medium. In batch culture, this process of self-reproduction is associated with the gradual decline in available metabolic energy and corresponding change in the physiological state of the population as a function of "travelled distance" along the autopoietic path. We argue that this directional alteration of cell physiology is both reflected in and supported by sequential gene expression along the chromosomal OriC-Ter axis. We propose that during the E. coli growth cycle, the spatiotemporal order of gene expression is established by coupling the temporal gradient of supercoiling energy to the spatial gradient of DNA thermodynamic stability along the chromosomal OriC-Ter axis.

Atomistic Molecular Dynamics Simulations of DNA in Complex 3D Arrangements for Comparison with Lower Resolution Structural Experiments

Atomic-level computer simulations are a very useful tool for describing the structure and dynamics of complex biomolecules such as DNA and for providing detail at a resolution where experimental techniques cannot arrive. Molecular dynamics (MD) simulations of mechanically distorted DNA caused by agents like supercoiling and protein binding are computationally challenging due to the large size of the associated systems and timescales. However, nowadays they are achievable thanks to the efficient usage of GPU and to the improvements of continuum solvation models. This together with the concurrent improvements in the resolution of single-molecule experiments, such as atomic force microscopy (AFM), makes possible the convergence between the two. Here we present detailed protocols for doing so: for performing molecular dynamics (MD) simulations of DNA adopting complex three-dimensional arrangements and for comparing the outcome of the calculations with single-molecule experimental data with a lower resolution than atomic.

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.

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.

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

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

Topological tuning of DNA mobility in entangled solutions of supercoiled plasmids

Ring polymers in dense solutions are among the most intriguing problems in polymer physics. Because of its natural occurrence in circular form, DNA has been extensively used as a proxy to study the fundamental physics of ring polymers in different topological states. Yet, torsionally constrained-such as supercoiled-topologies have been largely neglected so far. The applicability of existing theoretical models to dense supercoiled DNA is thus unknown. Here, we address this gap by coupling large-scale molecular dynamics simulations with differential dynamic microscopy of entangled supercoiled DNA plasmids. We find that, unexpectedly, larger supercoiling increases the size of entangled plasmids and concomitantly induces an enhancement in DNA mobility. These findings are reconciled as due to supercoiling-driven asymmetric and double-folded plasmid conformations that reduce interplasmid entanglements and threadings. Our results suggest a way to topologically tune DNA mobility via supercoiling, thus enabling topological control over the (micro)rheology of DNA-based complex fluids.

Base-pair resolution analysis of the effect of supercoiling on DNA flexibility and major groove recognition by triplex-forming oligonucleotides

In the cell, DNA is arranged into highly-organised and topologically-constrained (supercoiled) structures. It remains unclear how this supercoiling affects the detailed double-helical structure of DNA, largely because of limitations in spatial resolution of the available biophysical tools. Here, we overcome these limitations, by a combination of atomic force microscopy (AFM) and atomistic molecular dynamics (MD) simulations, to resolve structures of negatively-supercoiled DNA minicircles at base-pair resolution. We observe that negative superhelical stress induces local variation in the canonical B-form DNA structure by introducing kinks and defects that affect global minicircle structure and flexibility. We probe how these local and global conformational changes affect DNA interactions through the binding of triplex-forming oligonucleotides to DNA minicircles. We show that the energetics of triplex formation is governed by a delicate balance between electrostatics and bonding interactions. Our results provide mechanistic insight into how DNA supercoiling can affect molecular recognition, that may have broader implications for DNA interactions with other molecular species.

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

Kink formation is essential in highly bent DNA complexed with gene regulatory proteins such as histones to release the bending stress stored within the DNA duplex. Local opening of the double-stranded DNA creates a sharp turn along the specific sequence, which leads to the global bending of the DNA strand. Despite the critical role of kink formation, it is still challenging to predict the position of kink formation for a given DNA sequence. In this study, we propose a theoretical model and perform molecular dynamics simulations to quantify the sequence-dependent kink probability of a strongly bent DNA. By incorporating the elastic bending energy and the sequence-specific thermodynamic parameters, we investigate the importance of the DNA sequence on kink formation. We find that the sequence with TA dinucleotide repeats flanked by GC steps increases the kink propensity by more than an order of magnitude under the same bending stress. The number of base pairs involved in the local opening is found to be coupled with the sequence-specific bubble formation free energy. Our study elucidates the molecular origin of the sequence heterogeneity on kink formation, which is fundamental to understanding protein-DNA recognition.

Measuring Internal Forces in Single-Stranded DNA: Application to a DNA Force Clamp

We present a new method for calculating internal forces in DNA structures using coarse-grained models and demonstrate its utility with the oxDNA model. The instantaneous forces on individual nucleotides are explored and related to model potentials, and using our framework, internal forces are calculated for two simple DNA systems and for a recently published nanoscopic force clamp. Our results highlight some pitfalls associated with conventional methods for estimating internal forces, which are based on elastic polymer models, and emphasize the importance of carefully considering secondary structure and ionic conditions when modeling the elastic behavior of single-stranded DNA. Beyond its relevance to the DNA nanotechnological community, we expect our approach to be broadly applicable to calculations of internal force in a variety of structures-from DNA to protein-and across other coarse-grained simulation models.

A review on nonlinear DNA physics

The study and the investigation of structural and dynamical properties of complex systems have attracted considerable interest among scientists in general and physicists and biologists in particular. The present review paper represents a broad overview of the research performed over the nonlinear dynamics of DNA, devoted to some different aspects of DNA physics and including analytical, quantum and computational tools to understand nonlinear DNA physics. We review in detail the semi-discrete approximation within helicoidal Peyrard-Bishop model and show that localized modulated solitary waves, usually called breathers, can emerge and move along the DNA. Since living processes occur at submolecular level, we then discuss a quantum treatment to address the problem of how charge and energy are transported on DNA and how they may play an important role for the functioning of living cells. While this problem has attracted the attention of researchers for a long time, it is still poorly understood how charge and energy transport can occur at distances comparable to the size of macromolecules. Here, we review a theory based on the mechanism of 'self-trapping' of electrons due to their interaction with mechanical (thermal) oscillation of the DNA structure. We also describe recent computational models that have been developed to capture nonlinear mechanics of DNA in vitro and in vivo, possibly under topological constraints. Finally, we provide some conjectures on potential future directions for this field.

How global DNA unwinding causes non-uniform stress distribution and melting of DNA

DNA unwinding is an important process that controls binding of proteins, gene expression and melting of double-stranded DNA. In a series of all-atom MD simulations on two DNA molecules containing a transcription start TATA-box sequence we demonstrate that application of a global restraint on the DNA twisting dramatically changes the coupling between helical parameters and the distribution of deformation energy along the sequence. Whereas only short range nearest-neighbor coupling is observed in the relaxed case, long-range coupling is induced in the globally restrained case. With increased overall unwinding the elastic deformation energy is strongly non-uniformly distributed resulting ultimately in a local melting transition of only the TATA box segment during the simulations. The deformation energy tends to be stored more in cytidine/guanine rich regions associated with a change in conformational substate distribution. Upon TATA box melting the deformation energy is largely absorbed by the melting bubble with the rest of the sequences relaxing back to near B-form. The simulations allow us to characterize the structural changes and the propagation of the elastic energy but also to calculate the associated free energy change upon DNA unwinding up to DNA melting. Finally, we design an Ising model for predicting the local melting transition based on empirical parameters. The direct comparison with the atomistic MD simulations indicates a remarkably good agreement for the predicted necessary torsional stress to induce a melting transition, for the position and length of the melted region and for the calculated associated free energy change between both approaches.

The emergence of sequence-dependent structural motifs in stretched, torsionally constrained DNA

The double-helical structure of DNA results from canonical base pairing and stacking interactions. However, variations from steady-state conformations resulting from mechanical perturbations in cells have physiological relevance but their dependence on sequence remains unclear. Here, we use molecular dynamics simulations showing sequence differences result in markedly different structural motifs upon physiological twisting and stretching. We simulate overextension on different sequences of DNA ((AA)12, (AT)12, (CC)12 and (CG)12) with supercoiling densities at 200 and 50 mM salt concentrations. We find that DNA denatures in the majority of stretching simulations, surprisingly including those with over-twisted DNA. GC-rich sequences are observed to be more stable than AT-rich ones, with the specific response dependent on the base pair order. Furthermore, we find that (AT)12 forms stable periodic structures with non-canonical hydrogen bonds in some regions and non-canonical stacking in others, whereas (CG)12 forms a stacking motif of four base pairs independent of supercoiling density. Our results demonstrate that 20-30% DNA extension is sufficient for breaking B-DNA around and significantly above cellular supercoiling, and that the DNA sequence is crucial for understanding structural changes under mechanical stress. Our findings have important implications for the activities of protein machinery interacting with DNA in all cells.

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.

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.

Base-Pair Mismatch Can Destabilize Small DNA Loops through Cooperative Kinking

Base-pair mismatch can relieve mechanical stress in highly strained DNA molecules, but how it affects their kinetic stability is not known. Using single-molecule fluorescence resonance energy transfer, we measured the lifetimes of tightly bent DNA loops with and without base-pair mismatch. Surprisingly, for loops captured by stackable sticky ends which leave single-stranded DNA breaks (or nicks) upon annealing, the mismatch decreased the loop lifetime despite reducing the overall bending stress, and the decrease was largest when the mismatch was placed at the DNA midpoint. These findings suggest that base-pair mismatch increases bending stress at the opposite side of the loop through an allosteric mechanism known as cooperative kinking. Based on this mechanism, we present a three-state model that explains the apparent dichotomy between thermodynamic and kinetic stability.

Long-range temporal correlation in Auditory Brainstem Responses to Spoken Syllable/da/

The speech auditory brainstem response (sABR) is an objective clinical tool to diagnose particular impairments along the auditory brainstem pathways. We explore the scaling behavior of the brainstem in response to synthetic /da/ stimuli using a proposed pipeline including Multifractal Detrended Moving Average Analysis (MFDMA) modified by Singular Value Decomposition. The scaling exponent confirms that all normal sABR are classified into the non-stationary process. The average Hurst exponent is H = 0:77 ± 0:12 at 68% confidence interval indicating long-range correlation which shows the first universality behavior of sABR. Our findings exhibit that fluctuations in the sABR series are dictated by a mechanism associated with long-term memory of the dynamic of the auditory system in the brainstem level. The q-dependency of h(q) demonstrates that underlying data sets have multifractal nature revealing the second universality behavior of the normal sABR samples. Comparing Hurst exponent of original sABR with the results of the corresponding shuffled and surrogate series, we conclude that its multifractality is almost due to the long-range temporal correlations which are devoted to the third universality. Finally, the presence of long-range correlation which is related to the slow timescales in the subcortical level and integration of information in the brainstem network is confirmed.

Dynamics of supercoiled DNA with complex knots: large-scale rearrangements and persistent multi-strand interlocking

Knots and supercoiling are both introduced in bacterial plasmids by catalytic processes involving DNA strand passages. While the effects on plasmid organization has been extensively studied for knotting and supercoiling taken separately, much less is known about their concurrent action. Here, we use molecular dynamics simulations and oxDNA, an accurate mesoscopic DNA model, to study the kinetic and metric changes introduced by complex (five-crossing) knots and supercoiling in 2 kbp-long DNA rings. We find several unexpected results. First, the conformational ensemble is dominated by two distinct states, differing in branchedness and knot size. Secondly, fluctuations between these states are as fast as the metric relaxation of unknotted rings. In spite of this, certain boundaries of knotted and plectonemically-wound regions can persist over much longer timescales. These pinned regions involve multiple strands that are interlocked by the cooperative action of topological and supercoiling constraints. Their long-lived character may be relevant for the simplifying action of topoisomerases.

The temperature dependence of the helical twist of DNA

DNA is the carrier of all cellular genetic information and increasingly used in nanotechnology. Quantitative understanding and optimization of its functions requires precise experimental characterization and accurate modeling of DNA properties. A defining feature of DNA is its helicity. DNA unwinds with increasing temperature, even for temperatures well below the melting temperature. However, accurate quantitation of DNA unwinding under external forces and a microscopic understanding of the corresponding structural changes are currently lacking. Here we combine single-molecule magnetic tweezers measurements with atomistic molecular dynamics and coarse-grained simulations to obtain a comprehensive view of the temperature dependence of DNA twist. Experimentally, we find that DNA twist changes by ΔTw(T) = (-11.0 ± 1.2)°/(°C·kbp), independent of applied force, in the range of forces where torque-induced melting is negligible. Our atomistic simulations predict ΔTw(T) = (-11.1 ± 0.3)°/(°C·kbp), in quantitative agreement with experiments, and suggest that the untwisting of DNA with temperature is predominantly due to changes in DNA structure for defined backbone substates, while the effects of changes in substate populations are minor. Coarse-grained simulations using the oxDNA framework yield a value of ΔTw(T) = (-6.4 ± 0.2)°/(°C·kbp) in semi-quantitative agreement with experiments.

Bend-Induced Twist Waves and the Structure of Nucleosomal DNA

Recent work indicates that twist-bend coupling plays an important role in DNA micromechanics. Here we investigate its effect on bent DNA. We provide an analytical solution of the minimum-energy shape of circular DNA, showing that twist-bend coupling induces sinusoidal twist waves. This solution is in excellent agreement with both coarse-grained simulations of minicircles and nucleosomal DNA data, which is bent and wrapped around histone proteins in a superhelical conformation. Our analysis shows that the observed twist oscillation in nucleosomal DNA, so far attributed to the interaction with the histone proteins, is an intrinsic feature of free bent DNA, and should be observable in other protein-DNA complexes.

DNA elasticity from coarse-grained simulations: The effect of groove asymmetry

It is well established that many physical properties of DNA at sufficiently long length scales can be understood by means of simple polymer models. One of the most widely used elasticity models for DNA is the twistable worm-like chain (TWLC), which describes the double helix as a continuous elastic rod with bending and torsional stiffness. An extension of the TWLC, which has recently received some attention, is the model by Marko and Siggia, who introduced an additional twist-bend coupling, expected to arise from the groove asymmetry. By performing computer simulations of two available versions of oxDNA, a coarse-grained model of nucleic acids, we investigate the microscopic origin of twist-bend coupling. We show that this interaction is negligible in the oxDNA version with symmetric grooves, while it appears in the oxDNA version with asymmetric grooves. Our analysis is based on the calculation of the covariance matrix of equilibrium deformations, from which the stiffness parameters are obtained. The estimated twist-bend coupling coefficient from oxDNA simulations is G=30±1 nm. The groove asymmetry induces a novel twist length scale and an associated renormalized twist stiffness κ_t≈80 nm, which is different from the intrinsic torsional stiffness C≈110 nm. This naturally explains the large variations on experimental estimates of the intrinsic stiffness performed in the past.

Chromosomal organization of transcription: in a nutshell

Early studies of transcriptional regulation focused on individual gene promoters defined specific transcription factors as central agents of genetic control. However, recent genome-wide data propelled a different view by linking spatially organized gene expression patterns to chromosomal dynamics. Therefore, the major problem in contemporary molecular genetics concerned with transcriptional gene regulation is to establish a unifying model that reconciles these two views. This problem, situated at the interface of polymer physics and network theory, requires development of an integrative methodology. In this review, we discuss recent achievements in classical model organism E. coli and provide some novel insights gained from studies of a bacterial plant pathogen, D. dadantii. We consider DNA topology and the basal transcription machinery as key actors of regulation, in which activation of functionally relevant genes is coupled to and coordinated with the establishment of extended chromosomal domains of coherent transcription. We argue that the spatial organization of genome plays a fundamental role in its own regulation.

Evaluating the role of coherent delocalized phonon-like modes in DNA cyclization

The innate flexibility of a DNA sequence is quantified by the Jacobson-Stockmayer's J-factor, which measures the propensity for DNA loop formation. Recent studies of ultra-short DNA sequences revealed a discrepancy of up to six orders of magnitude between experimentally measured and theoretically predicted J-factors. These large differences suggest that, in addition to the elastic moduli of the double helix, other factors contribute to loop formation. Here, we develop a new theoretical model that explores how coherent delocalized phonon-like modes in DNA provide single-stranded "flexible hinges" to assist in loop formation. We combine the Czapla-Swigon-Olson structural model of DNA with our extended Peyrard-Bishop-Dauxois model and, without changing any of the parameters of the two models, apply this new computational framework to 86 experimentally characterized DNA sequences. Our results demonstrate that the new computational framework can predict J-factors within an order of magnitude of experimental measurements for most ultra-short DNA sequences, while continuing to accurately describe the J-factors of longer sequences. Further, we demonstrate that our computational framework can be used to describe the cyclization of DNA sequences that contain a base pair mismatch. Overall, our results support the conclusion that coherent delocalized phonon-like modes play an important role in DNA cyclization.

Analyzing DNA curvature and its impact on the ionic environment: application to molecular dynamics simulations of minicircles

We propose a method for analyzing the magnitude and direction of curvature within nucleic acids, based on the curvilinear helical axis calculated by Curves+. The method is applied to analyzing curvature within minicircles constructed with varying degrees of over- or under-twisting. Using the molecular dynamics trajectories of three different minicircles, we are able to quantify how curvature varies locally both in space and in time. We also analyze how curvature influences the local environment of the minicircles, notably via increased heterogeneity in the ionic distributions surrounding the double helix. The approach we propose has been integrated into Curves+ and the utilities Canal (time trajectory analysis) and Canion (environmental analysis) and can be used to study a wide variety of static and dynamic structural data on nucleic acids.

Influence of DNA sequence on the structure of minicircles under torsional stress

The sequence dependence of the conformational distribution of DNA under various levels of torsional stress is an important unsolved problem. Combining theory and coarse-grained simulations shows that the DNA sequence and a structural correlation due to topology constraints of a circle are the main factors that dictate the 3D structure of a 336 bp DNA minicircle under torsional stress. We found that DNA minicircle topoisomers can have multiple bend locations under high torsional stress and that the positions of these sharp bends are determined by the sequence, and by a positive mechanical correlation along the sequence. We showed that simulations and theory are able to provide sequence-specific information about individual DNA minicircles observed by cryo-electron tomography (cryo-ET). We provided a sequence-specific cryo-ET tomogram fitting of DNA minicircles, registering the sequence within the geometric features. Our results indicate that the conformational distribution of minicircles under torsional stress can be designed, which has important implications for using minicircle DNA for gene therapy.

Interference between Triplex and Protein Binding to Distal Sites on Supercoiled DNA

We have explored the interdependence of the binding of a DNA triplex and a repressor protein to distal recognition sites on supercoiled DNA minicircles using MD simulations. We observe that the interaction between the two ligands through their influence on their DNA template is determined by a subtle interplay of DNA mechanics and electrostatics, that the changes in flexibility induced by ligand binding play an important role and that supercoiling can instigate additional ligand-DNA contacts that would not be possible in simple linear DNA sequences.

Protein/DNA interactions in complex DNA topologies: expect the unexpected

DNA supercoiling results in compacted DNA structures that can bring distal sites into close proximity. It also changes the local structure of the DNA, which can in turn influence the way it is recognised by drugs, other nucleic acids and proteins. Here, we discuss how DNA supercoiling and the formation of complex DNA topologies can affect the thermodynamics of DNA recognition. We then speculate on the implications for transcriptional control and the three-dimensional organisation of the genetic material, using examples from our own simulations and from the literature. We introduce and discuss the concept of coupling between the multiple length-scales associated with hierarchical nuclear structural organisation through DNA supercoiling and topology.

Protein/DNA interactions in complex DNA topologies: expect the unexpected

DNA supercoiling results in compacted DNA structures that can bring distal sites into close proximity. It also changes the local structure of the DNA, which can in turn influence the way it is recognised by drugs, other nucleic acids and proteins. Here, we discuss how DNA supercoiling and the formation of complex DNA topologies can affect the thermodynamics of DNA recognition. We then speculate on the implications for transcriptional control and the three-dimensional organisation of the genetic material, using examples from our own simulations and from the literature. We introduce and discuss the concept of coupling between the multiple length-scales associated with hierarchical nuclear structural organisation through DNA supercoiling and topology.