Structural diversity of supercoiled DNA

Preparation of topoisomers of short circular dsDNA with defined linking number by accurate topological control

Short DNA catenanes [circular double-stranded DNA (dsDNA)] have attracted considerable interest for constructing nanostructures and nanomachines, as well as understanding DNA topology. The study of topoisomers of a circular dsDNA with a definite linking number (Lk) is essential but very difficult for simplifying the complex problems about DNA topology. The topoisomers are difficult to prepare, especially in the case that two strands are completely complementary. In this study, using a model system, we prepared all eight topoisomers (Lk0-Lk7) of a 79-bp-long circular dsDNA (8-14 nm in size) by utilizing aid-DNA to prevent undesired hybridization. By rapid ligation before strand displacement, high selectivity (>75%) for most topoisomers (31% for Lk1) was achieved under the strict topological control. All eight topoisomers with high purity were obtained after purification. Using a gel shift assay with Z-DNA-specific binding proteins, as well as by circular dichroism chromatography and enzymatic digestion, it was found that Z-DNA forms for topoisomers Lk0-Lk6, and Lk0-Lk5 can be converted to Lk6 by topoisomerase I. The approach developed in this study can significantly contribute to DNA or RNA topology, particularly the effect of topological constraints on DNA structures and functions.

RNA sample optimization for cryo-EM analysis

RNAs play critical roles in most biological processes. Although the three-dimensional (3D) structures of RNAs primarily determine their functions, it remains challenging to experimentally determine these 3D structures due to their conformational heterogeneity and intrinsic dynamics. Cryogenic electron microscopy (cryo-EM) has recently played an emerging role in resolving dynamic conformational changes and understanding structure-function relationships of RNAs including ribozymes, riboswitches and bacterial and viral noncoding RNAs. A variety of methods and pipelines have been developed to facilitate cryo-EM structure determination of challenging RNA targets with small molecular weights at subnanometer to near-atomic resolutions. While a wide range of conditions have been used to prepare RNAs for cryo-EM analysis, correlations between the variables in these conditions and cryo-EM visualizations and reconstructions remain underexplored, which continue to hinder optimizations of RNA samples for high-resolution cryo-EM structure determination. Here we present a protocol that describes rigorous screenings and iterative optimizations of RNA preparation conditions that facilitate cryo-EM structure determination, supplemented by cryo-EM data processing pipelines that resolve RNA dynamics and conformational changes and RNA modeling algorithms that generate atomic coordinates based on moderate- to high-resolution cryo-EM density maps. The current protocol is designed for users with basic skills and experience in RNA biochemistry, cryo-EM and RNA modeling. The expected time to carry out this protocol may range from 3 days to more than 3 weeks, depending on the many variables described in the protocol. For particularly challenging RNA targets, this protocol could also serve as a starting point for further optimizations.

Atomic Description of the Reciprocal Action between Supercoils and Melting Bubbles on Linear DNA

Although the mechanical response of DNA to physiological torsion and tension is well characterized, the detailed structures are not yet known. By using molecular dynamics simulations on linear DNA with 300 base-pairs, we provide, for the first time, the conformational phase diagram at atomic resolution. Our simulations also reveal the dynamics and diffusion of supercoils. We observe a new state in negative supercoiling, where denaturation bubbles form in adenine/thymine-rich regions independently of the underlying DNA topology. We thus propose sequence-dependent bubbles could position plectonemes in longer DNA.

Structural basis of Cas9 DNA interrogation with a 5' truncated sgRNA

The efficiency and accuracy of CRISPR-Cas9 targeting varies considerably across genomic targets and remains a persistent issue for using this system in cells. Studies have shown that the use of 5' truncated single guide RNAs (sgRNAs) can reduce the rate of unwanted off-target recognition while still maintaining on-target specificity. However, it is not well-understood how reducing target complementarity enhances specificity or how truncation past 15 nucleotides (nts) prevents full Cas9 activation without compromising on-target binding. Here, we use biochemistry and cryogenic electron microscopy to investigate Cas9 structure and activity when bound to a 14-nt sgRNA. Our structures reveal that the shortened path of the displaced non-target strand (NTS) sterically occludes docking of the HNH L1 linker and prevents proper positioning of the nuclease domains. We show that cleavage inhibition can be alleviated by either artificially melting the protospacer adjacent motif (PAM)-distal duplex or providing a supercoiled substrate. Even though Cas9 forms a stable complex with its target, we find that plasmid cleavage is ∼1000-fold slower with a 14-nt sgRNA than with a full-length 20-nt sgRNA. Our results provide a structural basis for Cas9 target binding with 5' truncated sgRNAs and underline the importance of PAM-distal NTS availability in promoting Cas9 activation.

Spatio-temporal organization of the E. coli chromosome from base to cellular length scales

Escherichia coli has been a vital model organism for studying chromosomal structure, thanks, in part, to its small and circular genome (4.6 million base pairs) and well-characterized biochemical pathways. Over the last several decades, we have made considerable progress in understanding the intricacies of the structure and subsequent function of the E. coli nucleoid. At the smallest scale, DNA, with no physical constraints, takes on a shape reminiscent of a randomly twisted cable, forming mostly random coils but partly affected by its stiffness. This ball-of-spaghetti-like shape forms a structure several times too large to fit into the cell. Once the physiological constraints of the cell are added, the DNA takes on overtwisted (negatively supercoiled) structures, which are shaped by an intricate interplay of many proteins carrying out essential biological processes. At shorter length scales (up to about 1 kb), nucleoid-associated proteins organize and condense the chromosome by inducing loops, bends, and forming bridges. Zooming out further and including cellular processes, topological domains are formed, which are flanked by supercoiling barriers. At the megabase-scale both large, highly self-interacting regions (macrodomains) and strong contacts between distant but co-regulated genes have been observed. At the largest scale, the nucleoid forms a helical ellipsoid. In this review, we will explore the history and recent advances that pave the way for a better understanding of E. coli chromosome organization and structure, discussing the cellular processes that drive changes in DNA shape, and what contributes to compaction and formation of dynamic structures, and in turn how bacterial chromatin affects key processes such as transcription and replication.

DNA as a Double-Coding Device for Information Conversion and Organization of a Self-Referential Unity

Living systems are capable on the one hand of eliciting a coordinated response to changing environments (also known as adaptation), and on the other hand, they are capable of reproducing themselves. Notably, adaptation to environmental change requires the monitoring of the surroundings, while reproduction requires monitoring oneself. These two tasks appear separate and make use of different sources of information. Yet, both the process of adaptation as well as that of reproduction are inextricably coupled to alterations in genomic DNA expression, while a cell behaves as an indivisible unity in which apparently independent processes and mechanisms are both integrated and coordinated. We argue that at the most basic level, this integration is enabled by the unique property of the DNA to act as a double coding device harboring two logically distinct types of information. We review biological systems of different complexities and infer that the inter-conversion of these two distinct types of DNA information represents a fundamental self-referential device underlying both systemic integration and coordinated adaptive responses.

A high throughput single molecule platform to study DNA supercoiling effect on protein-DNA interactions

DNA supercoiling significantly influences DNA metabolic pathways. To examine its impact on DNA-protein interactions at the single-molecule level, we developed a highly efficient and reliable protocol to modify plasmid DNA at specific sites, allowing us to label plasmids with fluorophores and biotin. We then induced negative and positive supercoiling in these plasmids using gyrase and reverse gyrase, respectively. Comparing supercoiled DNA with relaxed circular DNA, we assessed the effects of supercoiling on CRISPR-Cas9 and mismatch repair protein MutS. We found that negative DNA supercoiling exacerbates off-target effects in DNA unwinding by Cas9. For MutS, we observed both negative and positive DNA supercoiling enhances the binding interaction between MutS and a mismatched base pair but does not affect the rate of ATP-induced sliding clamp formation. These findings not only underscore the versatility of our protocol but also opens new avenues for exploring the intricate dynamics of protein-DNA interactions under the influences of supercoiling.

Structuring lipid nanoparticles, DNA, and protein corona into stealth bionanoarchitectures for in vivo gene delivery

Lipid nanoparticles (LNPs) play a crucial role in addressing genetic disorders, and cancer, and combating pandemics such as COVID-19 and its variants. Yet, the ability of LNPs to effectively encapsulate large-size DNA molecules remains elusive. This is a significant limitation, as the successful delivery of large-size DNA holds immense potential for gene therapy. To address this gap, the present study focuses on the design of PEGylated LNPs, incorporating large-sized DNA, departing from traditional RNA and ionizable lipids. The resultant LNPs demonstrate a unique particle morphology. These particles were further engineered with a DNA coating and plasma proteins. This multicomponent bionanoconstruct exhibits enhanced transfection efficiency and safety in controlled laboratory settings and improved immune system evasion in in vivo tests. These findings provide valuable insights for the design and development of bionanoarchitectures for large-size DNA delivery, opening new avenues for transformative gene therapies.

Nick Induced Dynamics in Supercoiled DNA Facilitates the Protein Target Search Process

A DNA nick, defined as a discontinuity in a double-stranded DNA molecule where the phosphodiester bond between adjacent nucleotides of one strand is absent due to enzyme action, serves as an effective mechanism to alleviate stress in supercoiled DNA. This stress release is essential for the smooth operation of transcriptional machinery. However, the underlying mechanisms and their impact on protein search dynamics, which are crucial for initiating transcription, remain unclear. Through extensive computer simulations, we unravel the molecular picture, demonstrating that intramolecular stress release due to a DNA nick is driven by a combination of writhing and twisting motions, depending on the nick's position. This stress release is quantitatively manifested as a step-like increase in the linking number. Furthermore, we elucidate that the nicked supercoiled minicircles exhibit enhanced torsional dynamics, promoting rapid conformational changes and frequent shifts in the identities of juxtaposed DNA sites on the plectoneme. The dynamics of the juxtaposition sites facilitates communication between protein and DNA, resulting in faster protein diffusion compared with native DNA with the same topology. Our findings highlight the mechanistic intricacies and underscore the importance of DNA nicks in facilitating transcription elongation by actively managing torsional stress during DNA unwinding by the RNA polymerase.

DdrC, a unique DNA repair factor from D. radiodurans, senses and stabilizes DNA breaks through a novel lesion-recognition mechanism

The bacterium Deinococcus radiodurans is known to survive high doses of DNA damaging agents. This resistance is the result of robust antioxidant systems which protect efficient DNA repair mechanisms that are unique to Deinococcus species. The protein DdrC has been identified as an important component of this repair machinery. DdrC is known to bind to DNA in vitro and has been shown to circularize and compact DNA fragments. The mechanism and biological relevance of this activity is poorly understood. Here, we show that the DdrC homodimer is a lesion-sensing protein that binds to two single-strand (ss) or double-strand (ds) breaks. The immobilization of DNA breaks in pairs consequently leads to the circularization of linear DNA and the compaction of nicked DNA. The degree of compaction is directly proportional with the number of available nicks. Previously, the structure of the DdrC homodimer was solved in an unusual asymmetric conformation. Here, we solve the structure of DdrC under different crystallographic environments and confirm that the asymmetry is an endogenous feature of DdrC. We propose a dynamic structural mechanism where the asymmetry is necessary to trap a pair of lesions. We support this model with mutant disruption and computational modeling experiments.

Probing of nucleic acid compaction using low-frequency Raman spectroscopy

Compaction of nucleic acids, namely DNA and RNA, determines their functions and involvement in vital cell processes including transcription, replication, DNA repair and translation. However, experimental probing of the compaction of nucleic acids is not straightforward. In this study, we suggest an approach for this probing using low-frequency Raman spectroscopy. Specifically, we show theoretically, computationally and experimentally the quantifiable correlation between the low-frequency Raman intensity from nucleic acids, magnitude of thermal fluctuations of atomic positions, and the compaction state of biomolecules. Noteworthily, we highlight that the LF Raman intensity differs by an order of magnitude for different samples of DNA, and even for the same sample in the course of long-term storage. The feasibility of the approach is further shown by assessment of the DNA compaction in the nuclei of plant cells. We anticipate that the suggested approach will enlighten compaction of nucleic acids and their dynamics during the key processes of the cell life cycle and under various factors, facilitating advancement of molecular biology and medicine.

Angle between DNA linker and nucleosome core particle regulates array compaction revealed by individual-particle cryo-electron tomography

The conformational dynamics of nucleosome arrays generate a diverse spectrum of microscopic states, posing challenges to their structural determination. Leveraging cryogenic electron tomography (cryo-ET), we determine the three-dimensional (3D) structures of individual mononucleosomes and arrays comprising di-, tri-, and tetranucleosomes. By slowing the rate of condensation through a reduction in ionic strength, we probe the intra-array structural transitions that precede inter-array interactions and liquid droplet formation. Under these conditions, the arrays exhibite irregular zig-zag conformations with loose packing. Increasing the ionic strength promoted intra-array compaction, yet we do not observe the previously reported regular 30-nanometer fibers. Interestingly, the presence of H1 do not induce array compaction; instead, one-third of the arrays display nucleosomes invaded by foreign DNA, suggesting an alternative role for H1 in chromatin network construction. We also find that the crucial parameter determining the structure adopted by chromatin arrays is the angle between the entry and exit of the DNA and the corresponding tangents to the nucleosomal disc. Our results provide insights into the initial stages of intra-array compaction, a critical precursor to condensation in the regulation of chromatin organization.

Reactions of deoxyribonucleotide bases with sulfooxymethyl or halomethyl polycyclic aromatic hydrocarbons induce unwinding of DNA supercoils

Torsional stress in double-stranded DNA enables and regulates facets of chromosomal metabolism, replication, and transcription and requires regulatory enzymatic systems including topoisomerases and histone methyltransferases. As such, this machinery may be subject to deleterious effects from reactive mutagens, including ones from carcinogenic polycyclic aromatic hydrocarbon (PAH) adduct formation with DNA. Supercoiled plasmid DNA was investigated for its torsional responses to adducts formed in vitro from PAH benzylic carbocation reactive intermediates created spontaneously by release of leaving groups. PAH sulfate esters were found to (1) unwind DNA in a concentration dependent manner, and (2) provide maximum unwinding in a pattern consistent with known carcinogenicities of the parent PAHs, that is, 6-methylbenzo[a]pyrene > 7,12-methylbenz[a]anthracene > 3-methylcholanthrene > 9-methylanthracene > 7-methylbenz[a]anthracene > 1-methylpyrene. Supercoil unwinding was demonstrated to be dependent on the presence of sulfate or chloride leaving groups such that reactive carbocations were generated in situ by hydrolysis. In silico modeling of intercalative complex topology showed PAH benzylic carbocation reactive functional groups in alignment with target nucleophiles on guanine bases in a 5'-dCdG-3' pocket in agreement with known formation of nucleotide adducts. Inhibitory or modulatory effects on PAH-induced supercoil unwinding were seen with ascorbic acid and an experimental antineoplastic agent Antineoplaston A10 in agreement with their known anticarcinogenic properties. In summary, the reactive PAH intermediates studied here undoubtedly participate in well-known mutational mechanisms such as frameshifts and apurinic site generation. However, they are also capable of random disruption of chromosomal supercoiling in a manner consistent with the known carcinogenicities of the parent compounds, and this mechanism may represent an additional detrimental motif worthy of further study for a more complete understanding of chemical carcinogenicity.

Supercoiled ring polymers under shear flow

We apply monomer-resolved computer simulations of supercoiled ring polymers under shear, taking full account of the hydrodynamic interactions, accompanied, in parallel, by simulations in which these are switched off. The combination of bending and torsional rigidities inherent in these polymers, in conjunction with hydrodynamics, has a profound impact on their flow properties. In contrast to their flexible counterparts, which dramatically deform and inflate under shear [Liebetreu et al., Commun. Mater. 2020, 1, 4], supercoiled rings undergo only weak changes in their overall shape and they display both a reduced propensity to tumbling (at fixed Weissenberg number) and a much stronger orientational resistance with respect to their flexible counterparts. In the presence of hydrodynamic interactions, the coupling of the polymer to solvent flow is capable of bringing about a topological transformation of writhe to twist at strong shear upon conservation of the overall linking number.

Structural basis of DNA crossover capture by Escherichia coli DNA gyrase

DNA supercoiling must be precisely regulated by topoisomerases to prevent DNA entanglement. The interaction of type IIA DNA topoisomerases with two DNA molecules, enabling the transport of one duplex through the transient double-stranded break of the other, remains elusive owing to structures derived solely from single linear duplex DNAs lacking topological constraints. Using cryo-electron microscopy, we solved the structure of Escherichia coli DNA gyrase bound to a negatively supercoiled minicircle DNA. We show how DNA gyrase captures a DNA crossover, revealing both conserved molecular grooves that accommodate the DNA helices. Together with molecular tweezer experiments, the structure shows that the DNA crossover is of positive chirality, reconciling the binding step of gyrase-mediated DNA relaxation and supercoiling in a single structure.

Unveiling Sticholysin II and plasmid DNA interaction: Implications for developing non-viral vectors

Non-viral gene delivery systems offer significant potential for gene therapy due to their versatility, safety, and cost advantages over viral vectors. However, their effectiveness can be hindered by the challenge of efficiently releasing the genetic cargo from endosomes to prevent degradation in lysosomes. To overcome this obstacle, functional components can be incorporated into these systems. Sticholysin II (StII) is one of the pore-forming proteins derived from the sea anemone Stichodactyla helianthus, known for its high ability to permeabilize cellular and model membranes. In this study, we aimed to investigate the interaction between StII, and a model plasmid (pDNA) as an initial step towards designing an improved vector with enhanced endosomal escape capability. The electrophoretic mobility shift assay (EMSA) confirmed the formation of complexes between StII and pDNA. Computational predictions identified specific residues involved in the StII-DNA interaction interface, highlighting the importance of electrostatic interactions and hydrogen bonds in mediating the binding. Atomic force microscopy (AFM) of StII-pDNA complexes revealed the presence of nodular fiber and toroid shapes. These complexes were found to have a predominantly micrometer size, as confirmed by dynamic light scattering (DLS) measurements. Despite increase in the overall charge, the complexes formed at the evaluated nitrogen-to-phosphorus (N/P) ratios still maintained a negative charge. Moreover, StII retained its pore-forming capacity regardless of its binding to the complexes. These findings suggest that the potential ability of StII to permeabilize endosomal membranes could be largely maintained when combined with nucleic acid delivery systems. Additionally, the still remaining negative charge of the complexes would enable the association of another positively charged component to compact pDNA. However, to minimize non-specific cytotoxic effects, it is advisable to explore methods to regulate the protein's activity in response to the microenvironment.

Nonmonotonic electrophoretic mobility of rodlike polyelectrolytes by multivalent coions in added salt

It is well established that when multivalent counterions or salts are added to a solution of highly charged polyelectrolytes (PEs), correlation effects can cause charge inversion of the PE, leading to electrophoretic mobility (EM) reversal. In this work, we use coarse-grained molecular-dynamics simulations to unravel the less understood effect of coion valency on EM reversal for rigid DNA-like PEs. We find that EM reversal induced by multivalent counterions is suppressed with increasing coion valency in the salt added and eventually vanishes. Further, we find that EM is enhanced at fixed low salt concentrations for salts with monovalent counterions when multivalent coions with increasing valency are introduced. However, increasing the salt concentration causes a crossover that leads to EM reversal which is enhanced by increasing coion valency at high salt concentration. Remarkably, this multivalent coion-induced EM reversal persists even for low values of PE linear charge densities where multivalent counterions alone cannot induce EM reversal. These results facilitate tuning PE-PE interactions and self-assembly with both coion and counterion valencies.

DNA supercoiling in bacteria: state of play and challenges from a viewpoint of physics based modeling

DNA supercoiling is central to many fundamental processes of living organisms. Its average level along the chromosome and over time reflects the dynamic equilibrium of opposite activities of topoisomerases, which are required to relax mechanical stresses that are inevitably produced during DNA replication and gene transcription. Supercoiling affects all scales of the spatio-temporal organization of bacterial DNA, from the base pair to the large scale chromosome conformation. Highlighted in vitro and in vivo in the 1960s and 1970s, respectively, the first physical models were proposed concomitantly in order to predict the deformation properties of the double helix. About fifteen years later, polymer physics models demonstrated on larger scales the plectonemic nature and the tree-like organization of supercoiled DNA. Since then, many works have tried to establish a better understanding of the multiple structuring and physiological properties of bacterial DNA in thermodynamic equilibrium and far from equilibrium. The purpose of this essay is to address upcoming challenges by thoroughly exploring the relevance, predictive capacity, and limitations of current physical models, with a specific focus on structural properties beyond the scale of the double helix. We discuss more particularly the problem of DNA conformations, the interplay between DNA supercoiling with gene transcription and DNA replication, its role on nucleoid formation and, finally, the problem of scaling up models. Our primary objective is to foster increased collaboration between physicists and biologists. To achieve this, we have reduced the respective jargon to a minimum and we provide some explanatory background material for the two communities.

Suitability of double-stranded DNA as a molecular standard for the validation of analytical ultracentrifugation instruments

To address the current lack of validated molecular standards for analytical ultracentrifugation (AUC), we investigated the suitability of double-stranded DNA molecules. We compared the hydrodynamic properties of linear and circular DNA as a function of temperature. Negatively supercoiled, nicked, and linearized 333 and 339 bp minicircles were studied. We quantified the hydrodynamic properties of these DNAs at five different temperatures, ranging from 4 to 37 °C. To enhance the precision of our measurements, each sample was globally fitted over triplicates and five rotor speeds. The exceptional stability of DNA allowed each sample to be sedimented repeatedly over the course of several months without aggregation or degradation, and with excellent reproducibility. The sedimentation and diffusion coefficients of linearized and nicked minicircle DNA demonstrated a highly homogeneous sample, and increased with temperature, indicating a decrease in friction. The sedimentation of linearized DNA was the slowest; supercoiled DNA sedimented the fastest. With increasing temperature, the supercoiled samples shifted to slower sedimentation, but sedimented faster than nicked minicircles. These results suggest that negatively supercoiled DNA becomes less compact at higher temperatures. The supercoiled minicircles, as purified from bacteria, displayed heterogeneity. Therefore, supercoiled DNA isolated from bacteria is unsuitable as a molecular standard. Linear and nicked samples are well suited as a molecular standard for AUC and have exceptional colloidal stability in an AUC cell. Even after sixty experiments at different speeds and temperatures, measured over the course of 4 months, all topological states of DNA remained colloidal, and their concentrations remained essentially unchanged.

Two Different Ways of Stress Release in Supercoiled DNA Minicircles under DNA Nick

It is generally believed that DNA nick is an effective way to release stress in supercoiled DNA, resulting from the twisting motion that individual strands rotate around the axis of the DNA helix. Here, we use MD simulations based on the oxDNA model to investigate the relaxation of 336 bp supercoiled minicircular DNA under DNA nick. Our simulations show that stress release, characterized by the abrupt decrease in linking number, may be induced by two types of DNA motion depending on the nick position. Except for the twisting motion, there is a writhing motion, that is, double strands collectively rotating with one plectoneme removal, which may occur in the process of DNA relaxation with the nick position in the loop region. Moreover, the writhing motion is more likely to occur in the DNA with relatively high hardness, such as C-G pairs. Our simulation results uncover the relationship between structural transformation, stress release, and DNA motion during the dynamic process under DNA nick, indicating the influence of nick position on the relaxation of the supercoiled DNA.

An associative memory Hamiltonian model for DNA and nucleosomes

A model for DNA and nucleosomes is introduced with the goal of studying chromosomes from a single base level all the way to higher-order chromatin structures. This model, dubbed the Widely Editable Chromatin Model (WEChroM), reproduces the complex mechanics of the double helix including its bending persistence length and twisting persistence length, and their respective temperature dependence. The WEChroM Hamiltonian is composed of chain connectivity, steric interactions, and associative memory terms representing all remaining interactions leading to the structure, dynamics, and mechanical characteristics of the B-DNA. Several applications of this model are discussed to demonstrate its applicability. WEChroM is used to investigate the behavior of circular DNA in the presence of positive and negative supercoiling. We show that it recapitulates the formation of plectonemes and of structural defects that relax mechanical stress. The model spontaneously manifests an asymmetric behavior with respect to positive or negative supercoiling, similar to what was previously observed in experiments. Additionally, we show that the associative memory Hamiltonian is also capable of reproducing the free energy of partial DNA unwrapping from nucleosomes. WEChroM is designed to emulate the continuously variable mechanical properties of the 10nm fiber and, by virtue of its simplicity, is ready to be scaled up to molecular systems large enough to investigate the structural ensembles of genes. WEChroM is implemented in the OpenMM simulation toolkits and is freely available for public use.

Conformational Change of Nucleosome Arrays prior to Phase Separation

Chromatin phase transition serves as a regulatory mechanism for eukaryotic transcription. Understanding this process requires the characterization of the nucleosome array structure in response to external stimuli prior to phase separation. However, the intrinsic flexibility and heterogeneity hinders the arrays' structure determination. Here we exploit advances in cryogenic electron tomography (cryo-ET) to determine the three-dimensional (3D) structure of each individual particle of mono-, di-, tri-, and tetranucleosome arrays. Statistical analysis reveals the ionic strength changes the angle between the DNA linker and nucleosome core particle (NCP), which regulate the overall morphology of nucleosome arrays. The finding that one-third of the arrays in the presence of H1 contain an NCP invaded by foreign DNA suggests an alternative function of H1 in constructing nucleosomal networks. The new insights into the nucleosome conformational changes prior to the intermolecular interaction stage extends our understanding of chromatin phase separation regulation.

DNA supercoiling-induced shapes alter minicircle hydrodynamic properties

DNA in cells is organized in negatively supercoiled loops. The resulting torsional and bending strain allows DNA to adopt a surprisingly wide variety of 3-D shapes. This interplay between negative supercoiling, looping, and shape influences how DNA is stored, replicated, transcribed, repaired, and likely every other aspect of DNA activity. To understand the consequences of negative supercoiling and curvature on the hydrodynamic properties of DNA, we submitted 336 bp and 672 bp DNA minicircles to analytical ultracentrifugation (AUC). We found that the diffusion coefficient, sedimentation coefficient, and the DNA hydrodynamic radius strongly depended on circularity, loop length, and degree of negative supercoiling. Because AUC cannot ascertain shape beyond degree of non-globularity, we applied linear elasticity theory to predict DNA shapes, and combined these with hydrodynamic calculations to interpret the AUC data, with reasonable agreement between theory and experiment. These complementary approaches, together with earlier electron cryotomography data, provide a framework for understanding and predicting the effects of supercoiling on the shape and hydrodynamic properties of DNA.

Recognition of DNA Supercoil Handedness during Catenation Catalyzed by Type II Topoisomerases

Although the presence of catenanes (i.e., intermolecular tangles) in chromosomal DNA stabilizes interactions between daughter chromosomes, a lack of resolution can have serious consequences for genomic stability. In all species, from bacteria to humans, type II topoisomerases are the enzymes primarily responsible for catenating/decatenating DNA. DNA topology has a profound influence on the rate at which these enzymes alter the superhelical state of the double helix. Therefore, the effect of supercoil handedness on the ability of human topoisomerase IIα and topoisomerase IIβ and bacterial topoisomerase IV to catenate DNA was examined. Topoisomerase IIα preferentially catenated negatively supercoiled over positively supercoiled substrates. This is opposite to its preference for relaxing (i.e., removing supercoils from) DNA and may prevent the enzyme from tangling the double helix ahead of replication forks and transcription complexes. The ability of topoisomerase IIα to recognize DNA supercoil handedness during catenation resides in its C-terminal domain. In contrast to topoisomerase IIα, topoisomerase IIβ displayed little ability to distinguish DNA geometry during catenation. Topoisomerase IV from three bacterial species preferentially catenated positively supercoiled substrates. This may not be an issue, as these enzymes work primarily behind replication forks. Finally, topoisomerase IIα and topoisomerase IV maintain lower levels of covalent enzyme-cleaved DNA intermediates with catenated over monomeric DNA. This allows these enzymes to perform their cellular functions in a safer manner, as catenated daughter chromosomes may be subject to stress generated by the mitotic spindle that could lead to irreversible DNA cleavage.

Multiscale Unfolding: Illustratively Visualizing the Whole Genome at a Glance

We present Multiscale Unfolding, an interactive technique for illustratively visualizing multiple hierarchical scales of DNA in a single view, showing the genome at different scales and demonstrating how one scale spatially folds into the next. The DNA's extremely long sequential structure-arranged differently on several distinct scale levels-is often lost in traditional 3D depictions, mainly due to its multiple levels of dense spatial packing and the resulting occlusion. Furthermore, interactive exploration of this complex structure is cumbersome, requiring visibility management like cut-aways. In contrast to existing temporally controlled multiscale data exploration, we allow viewers to always see and interact with any of the involved scales. For this purpose we separate the depiction into constant-scale and scale transition zones. Constant-scale zones maintain a single-scale representation, while still linearly unfolding the DNA. Inspired by illustration, scale transition zones connect adjacent constant-scale zones via level unfolding, scaling, and transparency. We thus represent the spatial structure of the whole DNA macro-molecule, maintain its local organizational characteristics, linearize its higher-level organization, and use spatially controlled, understandable interpolation between neighboring scales. We also contribute interaction techniques that provide viewers with a coarse-to-fine control for navigating within our all-scales-in-one-view representations and visual aids to illustrate the size differences. Overall, Multiscale Unfolding allows viewers to grasp the DNA's structural composition from chromosomes to the atoms, with increasing levels of "unfoldedness," and can be applied in data-driven illustration and communication.

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.

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

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

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.

Twisting DNA by salt

The structure and properties of DNA depend on the environment, in particular the ion atmosphere. Here, we investigate how DNA twist -one of the central properties of DNA- changes with concentration and identity of the surrounding ions. To resolve how cations influence the twist, we combine single-molecule magnetic tweezer experiments and extensive all-atom molecular dynamics simulations. Two interconnected trends are observed for monovalent alkali and divalent alkaline earth cations. First, DNA twist increases monotonously with increasing concentration for all ions investigated. Second, for a given salt concentration, DNA twist strongly depends on cation identity. At 100 mM concentration, DNA twist increases as Na+ < K+ < Rb+ < Ba2+ < Li+ ≈ Cs+ < Sr2+ < Mg2+ < Ca2+. Our molecular dynamics simulations reveal that preferential binding of the cations to the DNA backbone or the nucleobases has opposing effects on DNA twist and provides the microscopic explanation of the observed ion specificity. However, the simulations also reveal shortcomings of existing force field parameters for Cs+ and Sr2+. The comprehensive view gained from our combined approach provides a foundation for understanding and predicting cation-induced structural changes both in nature and in DNA nanotechnology.

Revisiting DNA Sequence-Dependent Deformability in High-Resolution Structures: Effects of Flanking Base Pairs on Dinucleotide Morphology and Global Chain Configuration

DNA carries more than the list of biochemical ingredients that drive the basic functions of living systems. The sequence of base pairs includes a multitude of structural and energetic signals, which determine the degree to which the long, threadlike molecule moves and how it responds to proteins and other molecules that control its processing and govern its packaging. The chemical composition of base pairs directs the spatial disposition and fluctuations of successive residues. The observed arrangements of these moieties in high-resolution protein–DNA crystal structures provide one of the best available estimates of the natural, sequence-dependent structure and deformability of the double-helical molecule. Here, we update the set of knowledge-based elastic potentials designed to describe the observed equilibrium structures and configurational fluctuations of the ten unique base-pair steps. The large number of currently available structures makes it possible to characterize the configurational preferences of the DNA base-pair steps within the context of their immediate neighbors, i.e., tetrameric context. Use of these knowledge-based potentials shows promise in accounting for known effects of sequence in long chain molecules, e.g., the degree of curvature reported in classic gel mobility studies and the recently reported sequence-dependent responses of supercoiled minicircles to nuclease cleavage.

emDNA - A Tool for Modeling Protein-decorated DNA Loops and Minicircles at the Base-pair Step Level

Computational modeling of nucleic acids plays an important role in molecular biology, enhancing our general understanding of the relationship between structure and function. Biophysical studies have provided a wealth of information on how double-helical DNA responds to proteins and other molecules in its local environment but far less understanding of the larger scale structural responses found in protein-decorated loops and minicircles. Current computational models of DNA range from detailed all-atom molecular dynamics studies, which produce rich time and spatially dependent depictions of small DNA fragments, to coarse-grained simulations, which sacrifice detailed physical and chemical information to treat larger scale systems. The treatment of DNA used here, at the base-pair step level with rigid-body parameters, allows one to develop quality models hundreds of base pairs long from local, sequence-specific features found from experiment. The emDNA software takes advantage of this framework, producing elastically optimized models of DNA at thermal equilibrium with built-in or user-generated elastic models. This versatile program, in combination with case studies included in this article, allows users of any skill level to develop and investigate mesoscale models of their own design. The functionality of emDNA includes a tool to incorporate experiment-specific configurations, e.g. protein-bound and/or melted DNA from known high-resolution structures, within higher-order 3D models by fixing the orientation and position of user-specified base pairs. The software provides a new avenue into multiscale genetic modeling, giving a wide range of users a deeper understanding of DNA mesoscale organization and the opportunity to pose new questions in genetic research. The publicly available emDNA software, including build instructions and usage information, is available on GitHub (https://nicocvn.github.io/emDNA/).

DNA microcircles - The promising tool for in vivo studies of the behavior of non-canonical DNA

The paper discusses the reasons for the resurrection of the term DNA microcircles through the change of its definition to "topologically closed DNA circles with the length less than 1 Kbp" from the entire population of circular DNA that holds the name of minicircles. The possible applications of such tool for in vivo studies of non-canonical DNA are also discussed. Prospective for in vivo and in vitro studies of non-canonical DNA cloned into microcircles are demonstrated. A method of stepwise elongation or shortening of plasmids is discussed.

Non-equilibrium structural dynamics of supercoiled DNA plasmids exhibits asymmetrical relaxation

Many cellular processes occur out of equilibrium. This includes site-specific unwinding in supercoiled DNA, which may play an important role in gene regulation. Here, we use the Convex Lens-induced Confinement (CLiC) single-molecule microscopy platform to study these processes with high-throughput and without artificial constraints on molecular structures or interactions. We use two model DNA plasmid systems, pFLIP-FUSE and pUC19, to study the dynamics of supercoiling-induced secondary structural transitions after perturbations away from equilibrium. We find that structural transitions can be slow, leading to long-lived structural states whose kinetics depend on the duration and direction of perturbation. Our findings highlight the importance of out-of-equilibrium studies when characterizing the complex structural dynamics of DNA and understanding the mechanisms of gene regulation.

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.

Ordered Packing of β-Sheet Nanofibrils into Nanotubes: Multi-hierarchical Assembly of Designed Short Peptides

Although it is well-known proteins and their complexes are hierarchically organized and highly ordered structures, it remains a major challenge to replicate their hierarchical self-assembly process and to fabricate multihierarchical architectures with well-defined shapes and monodisperse characteristic sizes via peptide self-assembly. Here we describe an amphiphilic short peptide Ac-I₃GGHK-NH₂ that first preassembles into thin, left-handed β-sheet nanofibrils, followed by their ordered packing into right-handed nanotubes. The key intermediate morphology and structures featuring the hierarchical process are simultaneously demonstrated. Further mechanistic exploration with the variants Ac-I₃GGGK-NH₂, Ac-I₃GGFK-NH₂, and Ac-I₃GG^DH^DK-NH₂ reveals the vital role of multiple His-His side chain interactions between nanofibrils in mediating higher-order assembly and architectures. Altogether, our findings not only advance current understanding of hierarchical assembly of peptides and proteins but also afford a paradigm of how to take advantage of side chain interactions to construct higher-order assemblies with enhanced complexities.

The biological applications of DNA nanomaterials: current challenges and future directions

DNA, a genetic material, has been employed in different scientific directions for various biological applications as driven by DNA nanotechnology in the past decades, including tissue regeneration, disease prevention, inflammation inhibition, bioimaging, biosensing, diagnosis, antitumor drug delivery, and therapeutics. With the rapid progress in DNA nanotechnology, multitudinous DNA nanomaterials have been designed with different shape and size based on the classic Watson-Crick base-pairing for molecular self-assembly. Some DNA materials could functionally change cell biological behaviors, such as cell migration, cell proliferation, cell differentiation, autophagy, and anti-inflammatory effects. Some single-stranded DNAs (ssDNAs) or RNAs with secondary structures via self-pairing, named aptamer, possess the ability of targeting, which are selected by systematic evolution of ligands by exponential enrichment (SELEX) and applied for tumor targeted diagnosis and treatment. Some DNA nanomaterials with three-dimensional (3D) nanostructures and stable structures are investigated as drug carrier systems to delivery multiple antitumor medicine or gene therapeutic agents. While the functional DNA nanostructures have promoted the development of the DNA nanotechnology with innovative designs and preparation strategies, and also proved with great potential in the biological and medical use, there is still a long way to go for the eventual application of DNA materials in real life. Here in this review, we conducted a comprehensive survey of the structural development history of various DNA nanomaterials, introduced the principles of different DNA nanomaterials, summarized their biological applications in different fields, and discussed the current challenges and further directions that could help to achieve their applications in the future.

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.

Atomic force microscopy-A tool for structural and translational DNA research

Atomic force microscopy (AFM) is a powerful imaging technique that allows for structural characterization of single biomolecules with nanoscale resolution. AFM has a unique capability to image biological molecules in their native states under physiological conditions without the need for labeling or averaging. DNA has been extensively imaged with AFM from early single-molecule studies of conformational diversity in plasmids, to recent examinations of intramolecular variation between groove depths within an individual DNA molecule. The ability to image dynamic biological interactions in situ has also allowed for the interaction of various proteins and therapeutic ligands with DNA to be evaluated-providing insights into structural assembly, flexibility, and movement. This review provides an overview of how innovation and optimization in AFM imaging have advanced our understanding of DNA structure, mechanics, and interactions. These include studies of the secondary and tertiary structure of DNA, including how these are affected by its interactions with proteins. The broader role of AFM as a tool in translational cancer research is also explored through its use in imaging DNA with key chemotherapeutic ligands, including those currently employed in clinical practice.

Investigating site-selection mechanisms of retroviral integration in supercoiled DNA braids

We theoretically study the integration of short viral DNA in a DNA braid made up by two entwined double-stranded DNA molecules. We show that the statistics of single integration events substantially differ in the straight and buckled, or plectonemic, phase of the braid and are more likely in the latter. We further discover that integration is most likely close to plectoneme tips, where the larger bending energy helps overcome the associated energy barrier and that successive integration events are spatio-temporally correlated, suggesting a potential mechanistic explanation of clustered integration sites in host genomes. The braid geometry we consider provides a novel experimental set-up to quantify integration in a supercoiled substrate in vitro, and to better understand the role of double-stranded DNA topology during this process.

A Close Look at Molecular Self-Assembly with the Transmission Electron Microscope

Molecular self-assembly is pervasive in the formation of living and synthetic materials. Knowledge gained from research into the principles of molecular self-assembly drives innovation in the biological, chemical, and materials sciences. Self-assembly processes span a wide range of temporal and spatial domains and are often unintuitive and complex. Studying such complex processes requires an arsenal of analytical and computational tools. Within this arsenal, the transmission electron microscope stands out for its unique ability to visualize and quantify self-assembly structures and processes. This review describes the contribution that the transmission electron microscope has made to the field of molecular self-assembly. An emphasis is placed on which TEM methods are applicable to different structures and processes and how TEM can be used in combination with other experimental or computational methods. Finally, we provide an outlook on the current challenges to, and opportunities for, increasing the impact that the transmission electron microscope can have on molecular self-assembly.

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.

Topoisomerase II Poisons: Converting Essential Enzymes into Molecular Scissors

The extensive length, compaction, and interwound nature of DNA, together with its controlled and restricted movement in eukaryotic cells, create a number of topological issues that profoundly affect all of the functions of the genetic material. Topoisomerases are essential enzymes that modulate the topological structure of the double helix, including the regulation of DNA under- and overwinding and the removal of tangles and knots from the genome. Type II topoisomerases alter DNA topology by generating a transient double-stranded break in one DNA segment and allowing another segment to pass through the DNA gate. These enzymes are involved in a number of critical nuclear processes in eukaryotic cells, such as DNA replication, transcription, and recombination, and are required for proper chromosome structure and segregation. However, because type II topoisomerases generate double-stranded breaks in the genetic material, they also are intrinsically dangerous enzymes that have the capacity to fragment the genome. As a result of this dualistic nature, type II topoisomerases are the targets for a number of widely prescribed anticancer drugs. This article will describe the structure and catalytic mechanism of eukaryotic type II topoisomerases and will go on to discuss the actions of topoisomerase II poisons, which are compounds that stabilize DNA breaks generated by the type II enzyme and convert these essential enzymes into "molecular scissors." Topoisomerase II poisons represent a broad range of structural classes and include anticancer drugs, dietary components, and environmental chemicals.

On valorization of solvent extracts of Terminalia arjuna (arjuna) upon DNA scission and free radical scavenging improves coupling responses and cognitive functions under in vitro conditions

Chronic diseases have been treated using the phytochemical concepts of ethnomedicinal plant-derived herbal products. Terminalia arjuna, a significant ethnomedicinal plant, was revisited and reconnoitred for antioxidant, free radical scavenging, and DNA nicking inhibiting activity under H_2O2 conditions using 21 solvent extracts. Ferric reducing antioxidant power (FRAP), 2,2-diphenyl-1-picrylhydrazyl (DPPH), 2,2′-azino-bis-3-ethylbenzothiazoline-6-sulphonic acid, and nitrous oxide scavenging (%) were found to have a strong positive association and interaction (PCA 1 explains 84.54% variation) with ethanol bark (Etoh-AB) (Meoh-AF). TPC (144.67–1794 µg/mL GAE) and TFC (2.5–34 µM Fe(II)/g were highest in Etoh-AB. In a pattern of combined solvent extracts, Etoh-AB had the highest antioxidant capacity, accompanied by Etoh-AL ≥ Meoh-AB ≥ Dw-AF. With R² = 0.94, the DNA nicking inhibition behaviour parameters relative front, relative quantity, band (%), and lane (%) formed a positive significant (p < 0.01) connection. For the first time, we show that Etoh-AB nicks supercoiled, circular plasmid DNA in a way that is comparable to normal antioxidants. Normal antioxidants with the ability to prevent DNA nicking include Butylated hydroxy anisole < Butylated hydroxy toluene < ascorbic acid < and Gallic acid. Gallic acid (m/z 170.0208 g/mol) and Ellagic acid (m/z 302.0063 g/mol were present in high concentrations in solvent extracts. 0.48 mg was found to be the effective concentration for inhibiting relative DNA nicking. The current study is the first of its kind to show that steroid concentrations are higher in bark fractions of acetone, ethanol, and methanol. Furthermore, T. arjuna solvent extracts provide a wealth of information on phytochemical profiling, antioxidant ability, and DNA nicking inhibition, which may be useful for exploring the natural way and further research to develop a remedy against geriatric chronic disease. Despite the fact that ethanol is very close to methanol in terms of solvent toxicity, the current study identified it as the preferred solvent. Thus, the current research revisits previous studies and explores the potentiality of non-polar and polar aprotic and polar protic solvent systems, which lend credence to bioactive compounds that may be useful in isolating and formulating safe and cost effective herbal medicament for livestocks and aquaculture, and drugs for deoxygenerative human diseases, and can also be investigated further to instil environmental frugality.

Biomolecular Modeling and Simulation: A Prospering Multidisciplinary Field

We reassess progress in the field of biomolecular modeling and simulation, following up on our perspective published in 2011. By reviewing metrics for the field's productivity and providing examples of success, we underscore the productive phase of the field, whose short-term expectations were overestimated and long-term effects underestimated. Such successes include prediction of structures and mechanisms; generation of new insights into biomolecular activity; and thriving collaborations between modeling and experimentation, including experiments driven by modeling. We also discuss the impact of field exercises and web games on the field's progress. Overall, we note tremendous success by the biomolecular modeling community in utilization of computer power; improvement in force fields; and development and application of new algorithms, notably machine learning and artificial intelligence. The combined advances are enhancing the accuracy andscope of modeling and simulation, establishing an exemplary discipline where experiment and theory or simulations are full partners.

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.

WASP: a software package for correctly characterizing the topological development of ribbon structures

We introduce the Writhe Application Software Package (WASP) which can be used to characterisze the topology of ribbon structures, the underlying mathematical model of DNA, Biopolymers, superfluid vorticies, elastic ropes and magnetic flux ropes. This characterization is achieved by the general twist–writhe decomposition of both open and closed ribbons, in particular through a quantity termed the polar writhe. We demonstrate how this decomposition is far more natural and straightforward than artificial closure methods commonly utilized in DNA modelling. In particular, we demonstrate how the decomposition of the polar writhe into local and non-local components distinctly characterizes the local helical structure and knotting/linking of the ribbon. This decomposition provides additional information not given by alternative approaches. As example applications, the WASP routines are used to characterise the evolving topology (writhe) of DNA minicircle and open ended plectoneme formation magnetic/optical tweezer simulations, and it is shown that the decomponsition into local and non-local components is particularly important for the detection of plectonemes. Finally it is demonstrated that a number of well known alternative writhe expressions are actually simplifications of the polar writhe measure.