Bullied no more: when and how DNA shoves proteins around

Melting and Bubble Formation in a Double-Stranded DNA: Microscopic Aspects of Early Base-Pair Opening Events and the Role of Water

Despite its rigid structure, DNA is a remarkably flexible molecule. Flexibility is essential for biological functions (such as transcription and gene repair), which require large-amplitude structural changes such as bubble formation. The bubbles thus formed are required to have a certain stability of their own and survive long on the time scale of molecular motions. A molecular understanding of fluctuations leading to quasi-stable structures is not available. Through extensive atomistic molecular dynamics simulations, we identify a sequence of microscopic events that culminate in local bubble formation, which is initiated by base-pair (BP) opening, resulting from the cleavage of native BP hydrogen bonds (HBs). This is followed by the formation of mismatched BPs with non-native contacts. These metastable structures can either revert to their original forms or undergo a flipping transition to form a local bubble that can span across 3-4 BPs. A substantial distortion of the DNA backbone and a disruption of BP stacking are observed because of the structural changes induced by these local perturbations. We also explored how water helps in the entire process. A small number of water molecules undergo rearrangement to stabilize the intermediate states by forming HBs with DNA bases. Water thus acts as a lubricant that counteracts the enthalpic penalty suffered from the loss of native BP contacts. Although the process of bubble formation is reversible, the sequence of steps involved poses an entropic barrier, preventing it from easily retracing the path to the native state.

Non-canonical DNA structures in the human ribosomal DNA

Non-canonical structures (NCS) refer to the various forms of DNA that differ from the B-conformation described by Watson and Crick. It has been found that these structures are usual components of the genome, actively participating in its essential functions. The present review is focused on the nine kinds of NCS appearing or likely to appear in human ribosomal DNA (rDNA): supercoiling structures, R-loops, G-quadruplexes, i-motifs, DNA triplexes, cruciform structures, DNA bubbles, and A and Z DNA conformations. We discuss the conditions of their generation, including their sequence specificity, distribution within the locus, dynamics, and beneficial and detrimental role in the cell.

Extensive Bioinformatics Analyses Reveal a Phylogenetically Conserved Winged Helix (WH) Domain (Zτ) of Topoisomerase IIα, Elucidating Its Very High Affinity for Left-Handed Z-DNA and Suggesting Novel Putative Functions

The dynamic processes operating on genomic DNA, such as gene expression and cellular division, lead inexorably to topological challenges in the form of entanglements, catenanes, knots, "bubbles", R-loops, and other outcomes of supercoiling and helical disruption. The resolution of toxic topological stress is the function attributed to DNA topoisomerases. A prominent example is the negative supercoiling (nsc) trailing processive enzymes such as DNA and RNA polymerases. The multiple equilibrium states that nscDNA can adopt by redistribution of helical twist and writhe include the left-handed double-helical conformation known as Z-DNA. Thirty years ago, one of our labs isolated a protein from Drosophila cells and embryos with a 100-fold greater affinity for Z-DNA than for B-DNA, and identified it as topoisomerase II (gene Top2, orthologous to the human UniProt proteins TOP2A and TOP2B). GTP increased the affinity and selectivity for Z-DNA even further and also led to inhibition of the isomerase enzymatic activity. An allosteric mechanism was proposed, in which topoII acts as a Z-DNA-binding protein (ZBP) to stabilize given states of topological (sub)domains and associated multiprotein complexes. We have now explored this possibility by comprehensive bioinformatic analyses of the available protein sequences of topoII representing organisms covering the whole tree of life. Multiple alignment of these sequences revealed an extremely high level of evolutionary conservation, including a winged-helix protein segment, here denoted as Zτ, constituting the putative structural homolog of Zα, the canonical Z-DNA/Z-RNA binding domain previously identified in the interferon-inducible RNA Adenosine-to-Inosine-editing deaminase, ADAR1p150. In contrast to Zα, which is separate from the protein segment responsible for catalysis, Zτ encompasses the active site tyrosine of topoII; a GTP-binding site and a GxxG sequence motif are in close proximity. Quantitative Zτ-Zα similarity comparisons and molecular docking with interaction scoring further supported the "B-Z-topoII hypothesis" and has led to an expanded mechanism for topoII function incorporating the recognition of Z-DNA segments ("Z-flipons") as an inherent and essential element. We further propose that the two Zτ domains of the topoII homodimer exhibit a single-turnover "conformase" activity on given G(ate) B-DNA segments ("Z-flipins"), inducing their transition to the left-handed Z-conformation. Inasmuch as the topoII-Z-DNA complexes are isomerase inactive, we infer that they fulfill important structural roles in key processes such as mitosis. Topoisomerases are preeminent targets of anti-cancer drug discovery, and we anticipate that detailed elucidation of their structural-functional interactions with Z-DNA and GTP will facilitate the design of novel, more potent and selective anti-cancer chemotherapeutic agents.

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.

Recent Development in Biomedical Applications of Oligonucleotides with Triplex-Forming Ability

A DNA structure, known as triple-stranded DNA, is made up of three oligonucleotide chains that wind around one another to form a triple helix (TFO). Hoogsteen base pairing describes how triple-stranded DNA may be built at certain conditions by the attachment of the third strand to an RNA, PNA, or DNA, which might all be employed as oligonucleotide chains. In each of these situations, the oligonucleotides can be employed as an anchor, in conjunction with a specific bioactive chemical, or as a messenger that enables switching between transcription and replication through the triplex-forming zone. These data are also considered since various illnesses have been linked to the expansion of triplex-prone sequences. In light of metabolic acidosis and associated symptoms, some consideration is given to the impact of several low-molecular-weight compounds, including pH on triplex production in vivo. The review is focused on the development of biomedical oligonucleotides with triplexes.

A 2.8 Å Structure of Zoliflodacin in a DNA Cleavage Complex with Staphylococcus aureus DNA Gyrase

Since 2000, some thirteen quinolones and fluoroquinolones have been developed and have come to market. The quinolones, one of the most successful classes of antibacterial drugs, stabilize DNA cleavage complexes with DNA gyrase and topoisomerase IV (topo IV), the two bacterial type IIA topoisomerases. The dual targeting of gyrase and topo IV helps decrease the likelihood of resistance developing. Here, we report on a 2.8 Å X-ray crystal structure, which shows that zoliflodacin, a spiropyrimidinetrione antibiotic, binds in the same DNA cleavage site(s) as quinolones, sterically blocking DNA religation. The structure shows that zoliflodacin interacts with highly conserved residues on GyrB (and does not use the quinolone water-metal ion bridge to GyrA), suggesting it may be more difficult for bacteria to develop target mediated resistance. We show that zoliflodacin has an MIC of 4 µg/mL against Acinetobacter baumannii (A. baumannii), an improvement of four-fold over its progenitor QPT-1. The current phase III clinical trial of zoliflodacin for gonorrhea is due to be read out in 2023. Zoliflodacin, together with the unrelated novel bacterial topoisomerase inhibitor gepotidacin, is likely to become the first entirely novel chemical entities approved against Gram-negative bacteria in the 21st century. Zoliflodacin may also become the progenitor of a new safer class of antibacterial drugs against other problematic Gram-negative bacteria.

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.

The interplay of supercoiling and thymine dimers in DNA

Thymine dimers are a major mutagenic photoproduct induced by UV radiation. While they have been the subject of extensive theoretical and experimental investigations, questions of how DNA supercoiling affects local defect properties, or, conversely, how the presence of such defects changes global supercoiled structure, are largely unexplored. Here, we introduce a model of thymine dimers in the oxDNA forcefield, parametrized by comparison to melting experiments and structural measurements of the thymine dimer induced bend angle. We performed extensive molecular dynamics simulations of double-stranded DNA as a function of external twist and force. Compared to undamaged DNA, the presence of a thymine dimer lowers the supercoiling densities at which plectonemes and bubbles occur. For biologically relevant supercoiling densities and forces, thymine dimers can preferentially segregate to the tips of the plectonemes, where they enhance the probability of a localized tip-bubble. This mechanism increases the probability of highly bent and denatured states at the thymine dimer site, which may facilitate repair enzyme binding. Thymine dimer-induced tip-bubbles also pin plectonemes, which may help repair enzymes to locate damage. We hypothesize that the interplay of supercoiling and local defects plays an important role for a wider set of DNA damage repair systems.

Searching for New Z-DNA/Z-RNA Binding Proteins Based on Structural Similarity to Experimentally Validated Zα Domain

Z-DNA and Z-RNA are functionally important left-handed structures of nucleic acids, which play a significant role in several molecular and biological processes including DNA replication, gene expression regulation and viral nucleic acid sensing. Most proteins that have been proven to interact with Z-DNA/Z-RNA contain the so-called Zα domain, which is structurally well conserved. To date, only eight proteins with Zα domain have been described within a few organisms (including human, mouse, Danio rerio, Trypanosoma brucei and some viruses). Therefore, this paper aimed to search for new Z-DNA/Z-RNA binding proteins in the complete PDB structures database and from the AlphaFold2 protein models. A structure-based similarity search found 14 proteins with highly similar Zα domain structure in experimentally-defined proteins and 185 proteins with a putative Zα domain using the AlphaFold2 models. Structure-based alignment and molecular docking confirmed high functional conservation of amino acids involved in Z-DNA/Z-RNA, suggesting that Z-DNA/Z-RNA recognition may play an important role in a variety of cellular processes.

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.

Sequence-dependent twist-bend coupling in DNA minicircles

Looping of double-stranded DNA molecules with 100-200 base pairs into minicircles, catenanes, and rotaxanes has been suggested as a potential tool for DNA nanotechnologies. However, sharp DNA bending into a minicircle with a diameter of several to ten nanometers occurs with alterations in the DNA helical structure and may lead to defective kink formation that hampers the use of DNA minicircles, catenanes, and rotaxanes in nanoscale DNA applications. Here, we investigated local variations of a helical twist in sharply bent DNA using microsecond-long all-atom molecular dynamics simulations of six different DNA minicircles, focusing on the sequence dependence of the coupling between DNA bending and its helical twist. Twist angles between consecutive base pairs were analyzed at different locations relative to the direction of DNA bending and, among 10 unique dinucleotide steps, we identified four dinucleotide steps with strong twist-bend coupling, the pyrimidine-purine dinucleotide steps of TA/TA, CG/CG, and CA/TG and the purine-purine dinucleotide step of GA/TC. This work suggests the sequence-dependent structural responses of DNA to strong mechanical deformation, providing new molecular-level insights into the structure and stability of sharply bent DNA minicircles for nanoscale applications.

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 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.

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.

Traveling Salesman Finds Random Walk: A Curve Reconstruction Algorithm for Supercoiled DNA

DNA supercoiling plays an important role in a variety of cellular processes, including transcription, replication, and DNA compaction. To fully understand these processes, we must uncover and characterize the dynamics of supercoiled DNA. However, supercoil dynamics are difficult to access because of the wide range of relevant length and timescales. In this work, we present an algorithm to reconstruct the arrangement of identical fluorescent particles distributed around a circular DNA molecule, given their three-dimensional trajectories through time. We find that this curve reconstruction problem is analogous to solving the traveling salesman problem. We demonstrate that our approach converges to the correct arrangement with a sufficiently long observation time. In addition, we show that the time required to accurately reconstruct the fluorophore arrangement is reduced by increasing the fluorophore density or reducing the level of supercoiling. This curve reconstruction algorithm, when paired with next-generation super-resolution imaging systems, could be used to access and thereby advance our understanding of supercoil dynamics.

Dynamics of Ionic Interactions at Protein-Nucleic Acid Interfaces

Molecular association of proteins with nucleic acids is required for many biological processes essential to life. Electrostatic interactions via ion pairs (salt bridges) of nucleic acid phosphates and protein side chains are crucial for proteins to bind to DNA or RNA. Counterions around the macromolecules are also key constituents for the thermodynamics of protein–nucleic acid association. Until recently, there had been only a limited amount of experiment-based information about how ions and ionic moieties behave in biological macromolecular processes. In the past decade, there has been significant progress in quantitative experimental research on ionic interactions with nucleic acids and their complexes with proteins. The highly negatively charged surfaces of DNA and RNA electrostatically attract and condense cations, creating a zone called the ion atmosphere. Recent experimental studies were able to examine and validate theoretical models on ions and their mobility and interactions with macromolecules. The ionic interactions are highly dynamic. The counterions rapidly diffuse within the ion atmosphere. Some of the ions are released from the ion atmosphere when proteins bind to nucleic acids, balancing the charge via intermolecular ion pairs of positively charged side chains and negatively charged backbone phosphates. Previously, the release of counterions had been implicated indirectly by the salt-concentration dependence of the association constant.

Recently, direct detection of counterion release by NMR spectroscopy has become possible and enabled more accurate and quantitative analysis of the counterion release and its entropic impact on the thermodynamics of protein–nucleic acid association. Recent studies also revealed the dynamic nature of ion pairs of protein side chains and nucleic acid phosphates. These ion pairs undergo transitions between two major states. In one of the major states, the cation and the anion are in direct contact and form hydrogen bonds. In the other major state, the cation and the anion are separated by water. Transitions between these states rapidly occur on a picosecond to nanosecond time scale. When proteins interact with nucleic acids, interfacial arginine (Arg) and lysine (Lys) side chains exhibit considerably different behaviors. Arg side chains show a higher propensity to form rigid contacts with nucleotide bases, whereas Lys side chains tend to be more mobile at the molecular interfaces. The dynamic ionic interactions may facilitate adaptive molecular recognition and play both thermodynamic and kinetic roles in protein–nucleic acid interactions.

Chromosomal Organization and Regulation of Genetic Function in Escherichia coli Integrates the DNA Analog and Digital Information

In this article, we summarize our current understanding of the bacterial genetic regulation brought about by decades of studies using the Escherichia coli model. It became increasingly evident that the cellular genetic regulation system is organizationally closed, and a major challenge is to describe its circular operation in quantitative terms. We argue that integration of the DNA analog information (i.e., the probability distribution of the thermodynamic stability of base steps) and digital information (i.e., the probability distribution of unique triplets) in the genome provides a key to understanding the organizational logic of genetic control. During bacterial growth and adaptation, this integration is mediated by changes of DNA supercoiling contingent on environmentally induced shifts in intracellular ionic strength and energy charge. More specifically, coupling of dynamic alterations of the local intrinsic helical repeat in the structurally heterogeneous DNA polymer with structural-compositional changes of RNA polymerase holoenzyme emerges as a fundamental organizational principle of the genetic regulation system. We present a model of genetic regulation integrating the genomic pattern of DNA thermodynamic stability with the gene order and function along the chromosomal OriC-Ter axis, which acts as a principal coordinate system organizing the regulatory interactions in the genome.

Transmission of dynamic supercoiling in linear and multi-way branched DNAs and its regulation revealed by a fluorescent G-quadruplex torsion sensor

DNA supercoiling is an important regulator of gene activity. The transmission of transcription-generated supercoiling wave along a DNA helix provides a way for a gene being transcribed to communicate with and regulate its neighboring genes. Currently, the dynamic behavior of supercoiling transmission remains unclear owing to the lack of a suitable tool for detecting the dynamics of supercoiling transmission. In this work, we established a torsion sensor that quantitatively monitors supercoiling transmission in real time in DNA. Using this sensor, we studied the transmission of transcriptionally generated negative supercoiling in linear and multi-way DNA duplexes. We found that transcription-generated dynamic supercoiling not only transmits along linear DNA duplex but also equally diverges at and proceeds through multi-way DNA junctions. We also show that such a process is regulated by DNA–protein interactions and non-canonical DNA structures in the path of supercoiling transmission. These results imply a transcription-coupled mechanism of dynamic supercoiling-mediated intra- and inter-chromosomal signal transduction pathway and their regulation in DNA.

BIOPHYSICS MEETS GENE THERAPY: HOW EXPLORING SUPERCOILING-DEPENDENT STRUCTURAL CHANGES IN DNA LED TO THE DEVELOPMENT OF MINIVECTOR DNA

Supercoiling affects every aspect of DNA function (replication, transcription, repair, recombination, etc.), yet the vast majority of studies on DNA and crystal structures of the molecule utilize short linear duplex DNA, which cannot be supercoiled. To study how supercoiling drives DNA biology, we developed and patented methods to make milligram quantities of tiny supercoiled circles of DNA called minicircles. We used a collaborative and multidisciplinary approach, including computational simulations (both atomistic and coarse-grained), biochemical experimentation, and biophysical methods to study these minicircles. By determining the three-dimensional conformations of individual supercoiled DNA minicircles, we revealed the structural diversity of supercoiled DNA and its highly dynamic nature. We uncovered profound structural changes, including sequence-specific base-flipping (where the DNA base flips out into the solvent), bending, and denaturing in negatively supercoiled minicircles. Counterintuitively, exposed DNA bases emerged in the positively supercoiled minicircles, which may result from inside-out DNA (Pauling-like, or "P-DNA"). These structural changes strongly influence how enzymes interact with or act on DNA. We hypothesized that, because of their small size and lack of bacterial sequences, these small supercoiled DNA circles may be efficient at delivering DNA into cells for gene therapy applications. "Minivectors," as we named them for this application, have proven to have therapeutic potential. We discovered that minivectors efficiently transfect a wide range of cell types, including many clinically important cell lines that are refractory to transfection with conventional plasmid vectors. Minivectors can be aerosolized for delivery to lungs and transfect human cells in culture to express RNA or genes. Importantly, minivectors demonstrate no obvious vector-associated toxicity. Minivectors can be repeatedly delivered and are long-lasting without integrating into the genome. Requests from colleagues around the world for minicircle and minivector DNA revealed a demand for our invention. We successfully obtained start-up funding for Twister Biotech, Inc. to help fulfill this demand, providing DNA for those who needed it, with a long-term goal of developing human therapeutics. In summary, what started as a tool for studying DNA structure has taken us in new and unanticipated directions.

Investigating DNA supercoiling in eukaryotic genomes

Supercoiling is a fundamental property of DNA, generated by polymerases and other DNA-binding proteins as a consequence of separating/bending the DNA double helix. DNA supercoiling plays a key role in gene expression and genome organization, but has proved difficult to study in eukaryotes because of the large, complex and chromatinized genomes. Key approaches to study DNA supercoiling in eukaryotes are (1) centrifugation-based or electrophoresis-based techniques in which supercoiled plasmids extracted from eukaryotic cells form a compacted writhed structure that migrates at a rate proportional to the level of DNA supercoiling; (2) in vivo approaches based on the preferential intercalation of psoralen molecules into under-wound DNA. Here, we outline the principles behind these techniques and discuss key discoveries, which have confirmed the presence and functional potential of unconstrained DNA supercoiling in eukaryotic genomes.

Atomistic Molecular Dynamics Simulations of DNA Minicircle Topoisomers: A Practical Guide to Setup, Performance, and Analysis

While DNA supercoiling is ubiquitous in vivo, the structure of supercoiled DNA is more challenging to study experimentally than simple linear sequences because the DNA must have a closed topology in order to sustain superhelical stress. DNA minicircles, which are closed circular double-stranded DNA sequences typically containing between 60 and 500 base pairs, have proven to be useful biochemical tools for the study of supercoiled DNA mechanics. We present detailed protocols for constructing models of DNA minicircles in silico, for performing atomistic molecular dynamics (MD) simulations of supercoiled minicircle DNA, and for analyzing the results of the calculations. These simulations are computationally challenging due to the large system sizes. However, improvements in parallel computing software and hardware promise access to improve conformational sampling and simulation timescales. Given the concurrent improvements in the resolution of experimental techniques such as atomic force microscopy (AFM) and cryo-electron microscopy, the study of DNA minicircles will provide a more complete understanding of both the structure and the mechanics of supercoiled DNA.

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.

The Bacillus subtilis Conjugative Plasmid pLS20 Encodes Two Ribbon-Helix-Helix Type Auxiliary Relaxosome Proteins That Are Essential for Conjugation

Bacterial conjugation is the process by which a conjugative element (CE) is transferred horizontally from a donor to a recipient cell via a connecting pore. One of the first steps in the conjugation process is the formation of a nucleoprotein complex at the origin of transfer (oriT), where one of the components of the nucleoprotein complex, the relaxase, introduces a site- and strand specific nick to initiate the transfer of a single DNA strand into the recipient cell. In most cases, the nucleoprotein complex involves, besides the relaxase, one or more additional proteins, named auxiliary proteins, which are encoded by the CE and/or the host. The conjugative plasmid pLS20 replicates in the Gram-positive Firmicute bacterium Bacillus subtilis. We have recently identified the relaxase gene and the oriT of pLS20, which are separated by a region of almost 1 kb. Here we show that this region contains two auxiliary genes that we name aux1LS20 and aux2LS20 , and which we show are essential for conjugation. Both Aux1LS20 and Aux2LS20 are predicted to contain a Ribbon-Helix-Helix DNA binding motif near their N-terminus. Analyses of the purified proteins show that Aux1LS20 and Aux2LS20 form tetramers and hexamers in solution, respectively, and that they both bind preferentially to oriTLS20 , although with different characteristics and specificities. In silico analyses revealed that genes encoding homologs of Aux1LS20 and/or Aux2LS20 are located upstream of almost 400 relaxase genes of the RelLS20 family (MOBL) of relaxases. Thus, Aux1LS20 and Aux2LS20 of pLS20 constitute the founding member of the first two families of auxiliary proteins described for CEs of Gram-positive origin.

Conformational rearrangement of 1,2-d(GG) intrastrand cis-diammineplatinum crosslinked DNA is driven by counter-ion penetration within the minor groove of the modified site

The major structural aberrations of DNA induced by a cis-diammineplatinum (II) 1,2-d(GG) intrastrand cross-link (CPT) have been known for decades. To gain deeper insights into the structural dynamics of the sequence-dependent DNA distortions adjacent to the CPT adduct, we employed molecular modeling and molecular dynamics (MD) simulations. The structural dynamics of native (N-DNA) and cisPt 1,2-d(GG) crosslinked (CPT-DNA) in the form of symmetric 36 nt d(G₂T₁₅G*G*T₁₅G₂)●C₂A₁₅CCA₁₅C₂) oligonucleotide duplexes is compared. The selected sequence context enabled tracking of the origin of the DNA axis curvature at the YpR flexible points (N-DNA), the enhancement of axis bending, and further distortions due to steric/electrostatic perturbations arising from the CPT-crosslink. In addition to the known structural distortions of CPT-DNA: helix bend towards the major groove; local helix unwinding; high roll angle between cross-linked guanine bases; and adoption of A-form DNA on the 5'-side of the CPT-crosslink (TpG junction); our results show the existence of a singular irreversible and reproducible conformational rearrangement, not previously observed, resulting in two stable CPT-DNA1 and CPT-DNA2 conformers. The CPT-DNA2 conformation presents an enhanced DNA axis bend and a wider and shallower minor grove with increased solvent accessibility within the modified site. It is concluded that the polymorphous (unstable) DNA environment near the cisPt 1,2-d(GG) unit in synergy with specific dynamic events, such as prolonged minor groove retention of particular Na⁺ ions and water redistribution within the d(TG*G*T) site, together with the formation of extra and more stable H-bonds between Pt(NH₃)₂ amines and neighboring nucleotides, are cooperatively responsible for the initiation of the conformational rearrangement leading to the CPT-DNA2 conformer, which, surprisingly, closely resembles the HMGB1-bound CPT-DNA structure. Graphical abstract Superimposed averaged structures of normal (N-DNA, green) and cisplatin intrastrand cross-linked (CPT-DNA, orange). Global DNA axes: N-DNA (blue beads); CPT-DNA (red beads); PT (yellow sphere); crosslinked dGs viewed from the minor groove (bold).

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.

Predicting conformational ensembles and genome-wide transcription factor binding sites from DNA sequences

DNA shape is emerging as an important determinant of transcription factor binding beyond just the DNA sequence. The only tool for large scale DNA shape estimates, DNAshape was derived from Monte-Carlo simulations and predicts four broad and static DNA shape features, Propeller twist, Helical twist, Minor groove width and Roll. The contributions of other shape features e.g. Shift, Slide and Opening cannot be evaluated using DNAshape. Here, we report a novel method DynaSeq, which predicts molecular dynamics-derived ensembles of a more exhaustive set of DNA shape features. We compared the DNAshape and DynaSeq predictions for the common features and applied both to predict the genome-wide binding sites of 1312 TFs available from protein interaction quantification (PIQ) data. The results indicate a good agreement between the two methods for the common shape features and point to advantages in using DynaSeq. Predictive models employing ensembles from individual conformational parameters revealed that base-pair opening - known to be important in strand separation - was the best predictor of transcription factor-binding sites (TFBS) followed by features employed by DNAshape. Of note, TFBS could be predicted not only from the features at the target motif sites, but also from those as far as 200 nucleotides away from the motif.

Discovery of a new family of relaxases in Firmicutes bacteria

Antibiotic resistance is a serious global problem. Antibiotic resistance genes (ARG), which are widespread in environmental bacteria, can be transferred to pathogenic bacteria via horizontal gene transfer (HGT). Gut microbiomes are especially apt for the emergence and dissemination of ARG. Conjugation is the HGT route that is predominantly responsible for the spread of ARG. Little is known about conjugative elements of Gram-positive bacteria, including those of the phylum Firmicutes, which are abundantly present in gut microbiomes. A critical step in the conjugation process is the relaxase-mediated site- and strand-specific nick in the oriT region of the conjugative element. This generates a single-stranded DNA molecule that is transferred from the donor to the recipient cell via a connecting channel. Here we identified and characterized the relaxosome components oriT and the relaxase of the conjugative plasmid pLS20 of the Firmicute Bacillus subtilis. We show that the relaxase gene, named relLS20, is essential for conjugation, that it can function in trans and provide evidence that Tyr26 constitutes the active site residue. In vivo and in vitro analyses revealed that the oriT is located far upstream of the relaxase gene and that the nick site within oriT is located on the template strand of the conjugation genes. Surprisingly, the RelLS20 shows very limited similarity to known relaxases. However, more than 800 genes to which no function had been attributed so far are predicted to encode proteins showing significant similarity to RelLS20. Interestingly, these putative relaxases are encoded almost exclusively in Firmicutes bacteria. Thus, RelLS20 constitutes the prototype of a new family of relaxases. The identification of this novel relaxase family will have an important impact in different aspects of future research in the field of HGT in Gram-positive bacteria in general, and specifically in the phylum of Firmicutes, and in gut microbiome research.

DNA Shape versus Sequence Variations in the Protein Binding Process

The binding process of a protein with a DNA involves three stages: approach, encounter, and association. It has been known that the complexation of protein and DNA involves mutual conformational changes, especially for a specific sequence association. However, it is still unclear how the conformation and the information in the DNA sequences affects the binding process. What is the extent to which the DNA structure adopted in the complex is induced by protein binding, or is instead intrinsic to the DNA sequence? In this study, we used the multiscale simulation method to explore the binding process of a protein with DNA in terms of DNA sequence, conformation, and interactions. We found that in the approach stage the protein can bind both the major and minor groove of the DNA, but uses different features to locate the binding site. The intrinsic conformational properties of the DNA play a significant role in this binding stage. By comparing the specific DNA with the nonspecific in unbound, intermediate, and associated states, we found that for a specific DNA sequence, ∼40% of the bending in the association forms is intrinsic and that ∼60% is induced by the protein. The protein does not induce appreciable bending of nonspecific DNA. In addition, we proposed that the DNA shape variations induced by protein binding are required in the early stage of the binding process, so that the protein is able to approach, encounter, and form an intermediate at the correct site on DNA.

The dynamic interplay between DNA topoisomerases and DNA topology

Topological properties of DNA influence its structure and biochemical interactions. Within the cell, DNA topology is constantly in flux. Transcription and other essential processes, including DNA replication and repair, not only alter the topology of the genome but also introduce additional complications associated with DNA knotting and catenation. These topological perturbations are counteracted by the action of topoisomerases, a specialized class of highly conserved and essential enzymes that actively regulate the topological state of the genome. This dynamic interplay among DNA topology, DNA processing enzymes, and DNA topoisomerases is a pervasive factor that influences DNA metabolism in vivo. Building on the extensive structural and biochemical characterization over the past four decades that has established the fundamental mechanistic basis of topoisomerase activity, scientists have begun to explore the unique roles played by DNA topology in modulating and influencing the activity of topoisomerases. In this review we survey established and emerging DNA topology-dependent protein-DNA interactions with a focus on in vitro measurements of the dynamic interplay between DNA topology and topoisomerase activity.

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.

The regulatory role of DNA supercoiling in nucleoprotein complex assembly and genetic activity

We argue that dynamic changes in DNA supercoiling in vivo determine both how DNA is packaged and how it is accessed for transcription and for other manipulations such as recombination. In both bacteria and eukaryotes, the principal generators of DNA superhelicity are DNA translocases, supplemented in bacteria by DNA gyrase. By generating gradients of superhelicity upstream and downstream of their site of activity, translocases enable the differential binding of proteins which preferentially interact with respectively more untwisted or more writhed DNA. Such preferences enable, in principle, the sequential binding of different classes of protein and so constitute an essential driver of chromatin organization.

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

It is well established that gene regulation can be achieved through activator and repressor proteins that bind to DNA and switch particular genes on or off, and that complex metabolic networks determine the levels of transcription of a given gene at a given time. Using three complementary computational techniques to study the sequence-dependence of DNA denaturation within DNA minicircles, we have observed that whenever the ends of the DNA are constrained, information can be transferred over long distances directly by the transmission of mechanical stress through the DNA itself, without any requirement for external signalling factors. Our models combine atomistic molecular dynamics (MD) with coarse-grained simulations and statistical mechanical calculations to span three distinct spatial resolutions and timescale regimes. While they give a consensus view of the non-locality of sequence-dependent denaturation in highly bent and supercoiled DNA loops, each also reveals a unique aspect of long-range informational transfer that occurs as a result of restraining the DNA within the closed loop of the minicircles.

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.

The Dynamic Interplay Between DNA Topoisomerases and DNA Topology

Topological properties of DNA influence its structure and biochemical interactions. Within the cell DNA topology is constantly in flux. Transcription and other essential processes including DNA replication and repair, alter the topology of the genome, while introducing additional complications associated with DNA knotting and catenation. These topological perturbations are counteracted by the action of topoisomerases, a specialized class of highly conserved and essential enzymes that actively regulate the topological state of the genome. This dynamic interplay among DNA topology, DNA processing enzymes, and DNA topoisomerases, is a pervasive factor that influences DNA metabolism in vivo. Building on the extensive structural and biochemical characterization over the past four decades that established the fundamental mechanistic basis of topoisomerase activity, the unique roles played by DNA topology in modulating and influencing the activity of topoisomerases have begun to be explored. In this review we survey established and emerging DNA topology dependent protein-DNA interactions with a focus on in vitro measurements of the dynamic interplay between DNA topology and topoisomerase activity.

Topological diversity of chromatin fibers: Interplay between nucleosome repeat length, DNA linking number and the level of transcription

The spatial organization of nucleosomes in 30-nm fibers remains unknown in detail. To tackle this problem, we analyzed all stereochemically possible configurations of two-start chromatin fibers with DNA linkers L = 10–70 bp (nucleosome repeat length NRL = 157–217 bp). In our model, the energy of a fiber is a sum of the elastic energy of the linker DNA, steric repulsion, electrostatics, and the H4 tail-acidic patch interaction between two stacked nucleosomes. We found two families of energetically feasible conformations of the fibers—one observed earlier, and the other novel. The fibers from the two families are characterized by different DNA linking numbers—that is, they are topologically different. Remarkably, the optimal geometry of a fiber and its topology depend on the linker length: the fibers with linkers L = 10n and 10n + 5 bp have DNA linking numbers per nucleosome ΔLk ≈ −1.5 and −1.0, respectively. In other words, the level of DNA supercoiling is directly related to the length of the inter-nucleosome linker in the chromatin fiber (and therefore, to NRL). We hypothesize that this topological polymorphism of chromatin fibers may play a role in the process of transcription, which is known to generate different levels of DNA supercoiling upstream and downstream from RNA polymerase. A genome-wide analysis of the NRL distribution in active and silent yeast genes yielded results consistent with this assumption.

Searching target sites on DNA by proteins: Role of DNA dynamics under confinement

DNA-binding proteins (DBPs) rapidly search and specifically bind to their target sites on genomic DNA in order to trigger many cellular regulatory processes. It has been suggested that the facilitation of search dynamics is achieved by combining 3D diffusion with one-dimensional sliding and hopping dynamics of interacting proteins. Although, recent studies have advanced the knowledge of molecular determinants that affect one-dimensional search efficiency, the role of DNA molecule is poorly understood. In this study, by using coarse-grained simulations, we propose that dynamics of DNA molecule and its degree of confinement due to cellular crowding concertedly regulate its groove geometry and modulate the inter-communication with DBPs. Under weak confinement, DNA dynamics promotes many short, rotation-decoupled sliding events interspersed by hopping dynamics. While this results in faster 1D diffusion, associated probability of missing targets by jumping over them increases. In contrast, strong confinement favours rotation-coupled sliding to locate targets but lacks structural flexibility to achieve desired specificity. By testing under physiological crowding, our study provides a plausible mechanism on how DNA molecule may help in maintaining an optimal balance between fast hopping and rotation-coupled sliding dynamics, to locate target sites rapidly and form specific complexes precisely.

Structural diversity of supercoiled DNA

By regulating access to the genetic code, DNA supercoiling strongly affects DNA metabolism. Despite its importance, however, much about supercoiled DNA (positively supercoiled DNA, in particular) remains unknown. Here we use electron cryo-tomography together with biochemical analyses to investigate structures of individual purified DNA minicircle topoisomers with defined degrees of supercoiling. Our results reveal that each topoisomer, negative or positive, adopts a unique and surprisingly wide distribution of three-dimensional conformations. Moreover, we uncover striking differences in how the topoisomers handle torsional stress. As negative supercoiling increases, bases are increasingly exposed. Beyond a sharp supercoiling threshold, we also detect exposed bases in positively supercoiled DNA. Molecular dynamics simulations independently confirm the conformational heterogeneity and provide atomistic insight into the flexibility of supercoiled DNA. Our integrated approach reveals the three-dimensional structures of DNA that are essential for its function.

DNA supercoiling strongly affects its metabolism. By electron cryo-tomography, biochemical assays and molecular dynamics simulations, here the authors show that supercoiled DNA minicircles adopt unique and wide distributions of three-dimensional conformations, many with disrupted base pairs.

The Potato Nucleotide-binding Leucine-rich Repeat (NLR) Immune Receptor Rx1 Is a Pathogen-dependent DNA-deforming Protein

Plant nucleotide-binding leucine-rich repeat (NLR) proteins enable cells to respond to pathogen attack. Several NLRs act in the nucleus; however, conserved nuclear targets that support their role in immunity are unknown. Previously, we noted a structural homology between the nucleotide-binding domain of NLRs and DNA replication origin-binding Cdc6/Orc1 proteins. Here we show that the NB-ARC (nucleotide-binding, Apaf-1, R-proteins, and CED-4) domain of the Rx1 NLR of potato binds nucleic acids. Rx1 induces ATP-dependent bending and melting of DNA in vitro, dependent upon a functional P-loop. In situ full-length Rx1 binds nuclear DNA following activation by its cognate pathogen-derived effector protein, the coat protein of potato virus X. In line with its obligatory nucleocytoplasmic distribution, DNA binding was only observed when Rx1 was allowed to freely translocate between both compartments and was activated in the cytoplasm. Immune activation induced by an unrelated NLR-effector pair did not trigger an Rx1-DNA interaction. DNA binding is therefore not merely a consequence of immune activation. These data establish a role for DNA distortion in Rx1 immune signaling and define DNA as a molecular target of an activated NLR.

Consistent rationalization of type-2 topoisomerases' unknotting, decatenating, supercoil-relaxing actions and their scaling relation

How type-2 topoisomerases discern global topology from local properties of DNA is not known precisely but the hypothesis that the enzymes selectively pass double-helix strands at hook-like juxtapositions is promising. Building upon an investigation of unknotting and decatenating using an improved wormlike DNA model, here we focus primarily on the enzymes' action in narrowing the distribution of linking number (Lk) in supercoiled DNA. Consistent with experiments, with selective passage at a hooked juxtaposition, the simulated narrowing factor RLk diminishes with decreasing DNA circle size but approaches an asymptotic RLk ≈ 1.7-1.8 for circle size ≳3.5 kb. For the larger DNA circles, we found that (RLk - 1) ≈ 0.42log10RK ≈ 0.68log10RL and thus RK ≈ (RL)(1.6) holds for the computed RLk and knot and catenane reduction factors RK and RL attained by selective passage at different juxtaposition geometries. Remarkably, this general scaling relation is essentially identical to that observed experimentally for several type-2 topoisomerases from a variety of organisms, indicating that the different disentangling powers of the topoisomerases likely arise from variations in the hooked geometries they select. Taken together, our results suggest strongly that type-2 topoisomerases recognize not only the curvature of the G-segment but also that of the T-segment.

Baculoviruses and nucleosome management

Negatively-supercoiled-ds DNA molecules, including the genomes of baculoviruses, spontaneously wrap around cores of histones to form nucleosomes when present within eukaryotic nuclei. Hence, nucleosome management should be essential for baculovirus genome replication and temporal regulation of transcription, but this has not been documented. Nucleosome mobilization is the dominion of ATP-dependent chromatin-remodeling complexes. SWI/SNF and INO80, two of the best-studied complexes, as well as chromatin modifier TIP60, all contain actin as a subunit. Retrospective analysis of results of AcMNPV time course experiments wherein actin polymerization was blocked by cytochalasin D drug treatment implicate actin-containing chromatin modifying complexes in decatenating baculovirus genomes, shutting down host transcription, and regulating late and very late phases of viral transcription. Moreover, virus-mediated nuclear localization of actin early during infection may contribute to nucleosome management.