Simulations of DNA topoisomerase 1B bound to supercoiled DNA reveal changes in the flexibility pattern of the enzyme and a secondary protein-DNA binding site

G-quadruplexes on chromosomal DNA negatively regulates topoisomerase 1 activity

Human DNA topoisomerase 1 (Top1) is a crucial enzyme responsible for alleviating torsional stress on DNA during transcription and replication, thereby maintaining genome stability. Previous researches had found that non-working Top1 interacted extensively with chromosomal DNA in human cells. However, the reason for its retention on chromosomal DNA remained unclear. In this study, we discovered a close association between Top1 and chromosomal DNA, specifically linked to the presence of G-quadruplex (G4) structures. G4 structures, formed during transcription, trap Top1 and hinder its ability to relax neighboring DNAs. Disruption of the Top1-G4 interaction using G4 ligand relieved the inhibitory effect of G4 on Top1 activity, resulting in a further reduction of R-loop levels in cells. Additionally, the activation of Top1 through the use of a G4 ligand enhanced the toxicity of Top1 inhibitors towards cancer cells. Our study uncovers a negative regulation mechanism of human Top1 and highlights a novel pathway for activating Top1.

Transcriptional repression by a secondary DNA binding surface of DNA topoisomerase I safeguards against hypertranscription

Regulation of global transcription output is important for normal development and disease, but little is known about the mechanisms involved. DNA topoisomerase I (TOP1) is an enzyme well-known for its role in relieving DNA supercoils for enabling transcription. Here, we report a non-enzymatic function of TOP1 that downregulates RNA synthesis. This function is dependent on specific DNA-interacting residues located on a conserved protein surface. A loss-of-function knock-in mutation on this surface, R548Q, is sufficient to cause hypertranscription and alter differentiation outcomes in mouse embryonic stem cells (mESCs). Hypertranscription in mESCs is accompanied by reduced TOP1 chromatin binding and change in genomic supercoiling. Notably, the mutation does not impact TOP1 enzymatic activity; rather, it diminishes TOP1-DNA binding and formation of compact protein-DNA structures. Thus, TOP1 exhibits opposing influences on transcription through distinct activities which are likely to be coordinated. This highlights TOP1 as a safeguard of appropriate total transcription levels in cells.

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.

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

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

Genomic Signature of Shifts in Selection in a Subalpine Ant and Its Physiological Adaptations

Understanding how organisms adapt to extreme environments is fundamental and can provide insightful case studies for both evolutionary biology and climate-change biology. Here, we take advantage of the vast diversity of lifestyles in ants to identify genomic signatures of adaptation to extreme habitats such as high altitude. We hypothesized two parallel patterns would occur in a genome adapting to an extreme habitat: 1) strong positive selection on genes related to adaptation and 2) a relaxation of previous purifying selection. We tested this hypothesis by sequencing the high-elevation specialist Tetramorium alpestre and four other phylogenetically related species. In support of our hypothesis, we recorded a strong shift of selective forces in T. alpestre, in particular a stronger magnitude of diversifying and relaxed selection when compared with all other ants. We further disentangled candidate molecular adaptations in both gene expression and protein-coding sequence that were identified by our genome-wide analyses. In particular, we demonstrate that T. alpestre has 1) a higher level of expression for stv and other heat-shock proteins in chill-shock tests and 2) enzymatic enhancement of Hex-T1, a rate-limiting regulatory enzyme that controls the entry of glucose into the glycolytic pathway. Together, our analyses highlight the adaptive molecular changes that support colonization of high-altitude environments.

Direct observation and analysis of TET-mediated oxidation processes in a DNA origami nanochip

DNA methylation and demethylation play a key role in the epigenetic regulation of gene expression; however, a series of oxidation reactions of 5-methyl cytosine (5mC) mediated by ten-eleven translocation (TET) enzymes driving demethylation process are yet to be uncovered. To elucidate the relationship between the oxidative processes and structural factors of DNA, we analysed the behavior of TET-mediated 5mC-oxidation by incorporating structural stress onto a substrate double-stranded DNA (dsDNA) using a DNA origami nanochip. The reactions and behaviors of TET enzymes were systematically monitored by biochemical analysis and single-molecule observation using atomic force microscopy (AFM). A reformative frame-like DNA origami was established to allow the incorporation of dsDNAs as 5mC-containing substrates in parallel orientations. We tested the potential effect of dsDNAs present in the tense and relaxed states within a DNA nanochip on TET oxidation. Based on enzyme binding and the detection of oxidation reactions within the DNA nanochip, it was revealed that TET preferred a relaxed substrate regardless of the modification types of 5-oxidated-methyl cytosine. Strikingly, when a multi-5mCG sites model was deployed to further characterize substrate preferences of TET, TET preferred the fully methylated site over the hemi-methylated site. This analytical modality also permits the direct observations of dynamic movements of TET such as sliding and interstrand transfer by high-speed AFM. In addition, the thymine DNA glycosylase-mediated base excision repair process was characterized in the DNA nanochip. Thus, we have convincingly established the system's ability to physically regulate enzymatic reactions, which could prove useful for the observation and characterization of coordinated DNA demethylation processes at the nanoscale.

Covalent Complex of DNA and Bacterial Topoisomerase: Implications in Antibacterial Drug Development

A topoisomerase-DNA transient covalent complex can be a druggable target for novel topoisomerase poison inhibitors that represent a new class of antibacterial or anticancer drugs. Herein, we have investigated molecular features of the functionally important Escherichia coli topoisomerase I (EctopoI)-DNA covalent complex (EctopoIcc) for molecular simulations, which is very useful in the development of new antibacterial drugs. To demonstrate the usefulness of our approach, we used a model small molecule (SM), NSC76027, obtained from virtual screening. We examined the direct binding of NSC76027 to EctopoI as well as inhibition of EctopoI relaxation activity of this SM via experimental techniques. We then performed molecular dynamics (MD) simulations to investigate the dynamics and stability of EctopoIcc and EctopoI-NSC76027-DNA ternary complex. Our simulation results show that NSC76027 forms a stable ternary complex with EctopoIcc. EctopoI investigated here also serves as a model system for investigating a complex of topoisomerase and DNA in which DNA is covalently attached to the protein.

Topoisomerase IB: a relaxing enzyme for stressed DNA

DNA topoisomerase I enzymes relieve the torsional strain in DNA; they are essential for fundamental molecular processes such as DNA replication, transcription, recombination, and chromosome condensation; and act by cleaving and then religating DNA strands. Over the past few decades, scientists have focused on the DNA topoisomerases biological functions and established a unique role of Type I DNA topoisomerases in regulating gene expression and DNA chromosome condensation. Moreover, the human enzyme is being investigated as a target for cancer chemotherapy. The active site tyrosine is responsible for initiating two transesterification reactions to cleave and then religate the DNA backbone, allowing the release of superhelical tension. The different steps of the catalytic mechanism are affected by various inhibitors; some of them prevent the interaction between the enzyme and the DNA while others act as poisons, leading to TopI-DNA lesions, breakage of DNA, and eventually cellular death. In this review, our goal is to provide an overview of mechanism of human topoisomerase IB action together with the different types of inhibitors and their effect on the enzyme functionality.

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

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

Advantages of an optical nanosensor system for the mechanistic analysis of a novel topoisomerase I targeting drug: a case study

The continuous need for the development of new small molecule anti-cancer drugs calls for easily accessible sensor systems for measuring the effect of vast numbers of new drugs on their potential cellular targets. Here we demonstrate the use of an optical DNA biosensor to unravel the inhibitory mechanism of a member of a new family of small molecule human topoisomerase I inhibitors, the so-called indeno-1,5-naphthyridines. By analysing human topoisomerase I catalysis on the biosensor in the absence or presence of added drug complemented with a few traditional assays, we demonstrate that the investigated member of the indeno-1,5-naphthyridine family inhibited human topoisomerase I activity by blocking enzyme-DNA dissociation. To our knowledge, this represents the first characterized example of a small molecule drug that inhibits a post-ligation step of catalysis. The elucidation of a completely new and rather surprising drug mechanism-of-action using an optical real time sensor highlights the value of this assay system in the search for new topoisomerase I targeting small molecule drugs.

A single nucleotide resolution model for large-scale simulations of double stranded DNA

The computational modelling of DNA is becoming crucial in light of new advances in DNA nano-technology, single-molecule experiments and in vivo DNA tampering. Here we present a mesoscopic model for double stranded DNA (dsDNA) at the single nucleotide level which retains the characteristic helical structure, while being able to simulate large molecules - up to a million base pairs - for time-scales which are relevant to physiological processes. This is made possible by an efficient and highly-parallelised implementation of the model which we discuss here. The model captures the main characteristics of DNA, such as the different persistence lengths for double and single strands, pitch, torsional rigidity and the presence of major and minor grooves. The model constitutes a starting point for the future implementation of further features, such as sequence specificity and electrostatic repulsion. We show that the behaviour of the presented model compares favourably with single molecule experiments where dsDNA is manipulated by external forces or torques. We finally present some results on the kinetics of denaturation of linear DNA and supercoiling of closed dsDNA molecules.

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.

Molecular basis of the recognition of FMN by a HAD phosphatase TON_0338

The haloalkaloic acid dehalogenase (HAD) phosphatase from Thermococcus onnurineus NA1 (TON_0338), has phosphatase activity the flavin mono-nucleotide (FMN). The molecular origin and structural motifs for the activity deficiency of double-tryptophan mutant have not been rationalized at atomic resolution. Molecular dynamics (MD) simulations and the molecular mechanics/Generalized-Born surface area (MM/GBSA) free energy calculations were used to explore the effects of mutations on the changes in both structural flexibility and conformational dynamics. The non-polar solvation energy plays an indispensable role in the binding process of TON_0338 and FMN. The tryptophan sandwich structure provides a primary function to anchor FMN and keeps FMN well bound to TON_0338. The double-tryptophan mutation has influences on the secondary structures of TON_0338 and changes the conformation, which would lead to reduced activity of W58A/W61A-FMN binding. The present study provides important insights into the structure-function relationships of TON_0338 protein, which could contribute to further understanding about the HAD phosphatases.

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.

Comparative simulative analysis of single and double stranded truncated octahedral DNA nanocages

Spacefill view of double (DSL) and single (SSL) stranded linkers DNA cages. The blue atoms represent the shared cages scaffold, while the yellow atoms show the single stranded DNA oligonucleotides shaping the double stranded linkers of the DSL cage.

Refinement of the Sugar-Phosphate Backbone Torsion Beta for AMBER Force Fields Improves the Description of Z- and B-DNA

Z-DNA duplexes are a particularly complicated test case for current force fields. We performed a set of explicit solvent molecular dynamics (MD) simulations with various AMBER force field parametrizations including our recent refinements of the ε/ζ and glycosidic torsions. None of these force fields described the ZI/ZII and other backbone substates correctly, and all of them underpredicted the population of the important ZI substate. We show that this underprediction can be attributed to an inaccurate potential for the sugar-phosphate backbone torsion angle β. We suggest a refinement of this potential, β(OL1), which was derived using our recently introduced methodology that includes conformation-dependent solvation effects. The new potential significantly increases the stability of the dominant ZI backbone substate and improves the overall description of the Z-DNA backbone. It also has a positive (albeit small) impact on another important DNA form, the antiparallel guanine quadruplex (G-DNA), and improves the description of the canonical B-DNA backbone by increasing the population of BII backbone substates, providing a better agreement with experiment. We recommend using β(OL1) in combination with our previously introduced corrections, εζ(OL1) and χ(OL4), (the combination being named OL15) as a possible alternative to the current β torsion potential for more accurate modeling of nucleic acids.

Insights into the effects of mutations on Cren7-DNA binding using molecular dynamics simulations and free energy calculations

A novel, highly conserved chromatin protein, Cren7 is involved in regulating essential cellular processes such as transcription, replication and repair. Although mutations in the DNA-binding loop of Cren7 destabilize the structure and reduce DNA-binding activity, the details are not very clear. Focusing on the specific Cren7-dsDNA complex (PDB code ), we applied molecular dynamics (MD) simulations and the molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) free energy calculations to explore the structural and dynamic effects of W26A, L28A, and K53A mutations in comparison to the wild-type protein. The energetic analysis indicated that the intermolecular van der Waals interaction and nonpolar solvation energy play an important role in the binding process of Cren7 and dsDNA. Compared with the wild type Cren7, all the studied mutants W26A, L28A, and K53A have obviously reduced binding free energies with dsDNA in the reduction of the polar and/or nonpolar interactions. These results further elucidated the previous experiments to understand the Cren7-DNA interaction comprehensively. Our work also would provide support for an understanding of the interactions of proteins with nucleic acids.

Multiscale modelling of DNA mechanics

Mechanical properties of DNA are important not only in a wide range of biological processes but also in the emerging field of DNA nanotechnology. We review some of the recent developments in modeling these properties, emphasizing the multiscale nature of the problem. Modern atomic resolution, explicit solvent molecular dynamics simulations have contributed to our understanding of DNA fine structure and conformational polymorphism. These simulations may serve as data sources to parameterize rigid base models which themselves have undergone major development. A consistent buildup of larger entities involving multiple rigid bases enables us to describe DNA at more global scales. Free energy methods to impose large strains on DNA, as well as bead models and other approaches, are also briefly discussed.

Recovery of the poisoned topoisomerase II for DNA religation: coordinated motion of the cleavage core revealed with the microsecond atomistic simulation

Type II topoisomerases resolve topological problems of DNA double helices by passing one duplex through the reversible double-stranded break they generated on another duplex. Despite the wealth of information in the cleaving operation, molecular understanding of the enzymatic DNA ligation remains elusive. Topoisomerase poisons are widely used in anti-cancer and anti-bacterial therapy and have been employed to entrap the intermediates of topoisomerase IIβ with religatable DNA substrate. We removed drug molecules from the structure and conducted molecular dynamics simulations to investigate the enzyme-mediated DNA religation. The drug-unbound intermediate displayed transitions toward the resealing-compliant configuration: closing distance between the cleaved DNA termini, B-to-A transformation of the double helix, and restoration of the metal-binding motif. By mapping the contact configurations and the correlated motions between enzyme and DNA, we identified the indispensable role of the linker preceding winged helix domain (WHD) in coordinating the movements of TOPRIM, the nucleotide-binding motifs, and the bound DNA substrate during gate closure. We observed a nearly vectorial transition in the recovery of the enzyme and identified the previously uncharacterized roles of Asn508 and Arg677 in DNA rejoining. Our findings delineate the dynamic mechanism of the DNA religation conducted by type II topoisomerases.