Role of the linker domain and the 203-214 N-terminal residues in the human topoisomerase I DNA complex dynamics

Molecular Docking and Dynamics Simulation Revealed the Potential Inhibitory Activity of New Drugs against Human Topoisomerase I Receptor

Human Topoisomerase I (hTop1p) is a ubiquitous enzyme that relaxes supercoiled DNA through a conserved mechanism involving transient breakage, rotation, and binding. Htop1p is the molecular target of the chemotherapeutic drug camptothecin (CPT). It causes the hTop1p-DNA complex to slow down the binding process and clash with the replicative machinery during the S phase of the cell cycle, forcing cells to activate the apoptotic response. This gives hTop1p a central role in cancer therapy. Recently, two artesunic acid derivatives (compounds c6 and c7) have been proposed as promising inhibitors of hTop1p with possible antitumor activity. We used several computational approaches to obtain in silico confirmations of the experimental data and to form a comprehensive dynamic description of the ligand-receptor system. We performed molecular docking analyses to verify the ability of the two new derivatives to access the enzyme-DNA interface, and a classical molecular dynamics simulation was performed to assess the capacity of the two compounds to maintain a stable binding pose over time. Finally, we calculated the noncovalent interactions between the two new derivatives and the hTop1p receptor in order to propose a possible inhibitory mechanism like that adopted by CPT.

Swapping of The N-Terminal Domain of Human Topoisomerase 1B with the Corresponding Plasmodium Falciparum Counterpart Strongly Impairs Enzyme Activity

BACKGROUND

DNA topoisomerases 1B are a class of ubiquitous enzyme that solves the topological problems associated with biological processes such as replication, transcription and recombination. Numerous sequence alignment of topoisomerase 1B from different species shows that the lengths of different domains as well as their amino acids sequences are quite different. In the present study a hybrid enzyme, generated by swapping the N-terminal of Plasmodium falciparum into the corresponding domain of the human, has been characterized.

METHODS

The chimeric enzyme was generated using different sets of PCR. The in vitro characterization was carried out using different DNA substrate including radio-labelled oligonucleotides.

RESULTS

The chimeric enzyme displayed slower relaxation activity, cleavage and re-ligation kinetics strongly perturbed when compared to the human enzyme.

CONCLUSION

These results indicate that the N-terminal domain has a crucial role in modulating topoisomerase activity in different species.

The human DNA topoisomerase I mutant Gly717Asp: Higher religation rate is not always associated with camptothecin resistance

DNA topoisomerases are key enzyme responsible for modulating the topological state of the DNA by breaking and rejoining of DNA strand. Characterization of a Gly717Asp mutation in the human topoisomerase was performed using several catalytic assays. The mutant enzyme was shown to have comparable cleavage and fast religation rate as compared to the wild-type protein. Addition of the anticancer drug camptothecin significantly reduced the religation step. The simulative approaches and analysis of the cleavage/religation equilibrium indicate that the mutation is able to modify the architecture of the drug binding site, increasing the persistence of the drug for the enzyme-DNA covalent complex. Taken together these results indicate that the structure modification of the drug binding site is the key reason for the increasing CPT persistence and furthermore provide the possibility for new anti-cancer drug discovery.

Topoisomerase 1B poisons: Over a half-century of drug leads, clinical candidates, and serendipitous discoveries

Topoisomerases are DNA processing enzymes that relieve supercoiling (torsional strain) in DNA, are necessary for normal cellular division, and act by nicking (and then religating) DNA strands. Type 1B topoisomerase (Top1) is overexpressed in certain tumors, and the enzyme has been extensively investigated as a target for cancer chemotherapy. Various chemical agents can act as "poisons" of the enzyme's religation step, leading to Top1-DNA lesions, DNA breakage, and eventual cellular death. In this review, agents that poison Top1 (and have thus been investigated for their anticancer properties) are surveyed, including natural products (such as camptothecins and indolocarbazoles), semisynthetic camptothecin and luotonin derivatives, and synthetic compounds (such as benzonaphthyridines, aromathecins, and indenoisoquinolines), as well as targeted therapies and conjugates. Top1 has also been investigated as a therapeutic target in certain viral and parasitic infections, as well as autoimmune, inflammatory, and neurological disorders, and a summary of literature describing alternative indications is also provided. This review should provide both a reference for the medicinal chemist and potentially offer clues to aid in the development of new Top1 poisons.

Substitutions in Spodoptera exigua topoisomerase I modulate its relaxation activity and camptothecin sensitivity

BACKGROUND

Topoisomerase I (Top I) is referred as the cellular target of the camptothecins (CPTs) which are now being explored as potential pesticides for insect control. Three amino acid substitutions, including L530P, A653T and S729T, in Top Is of insects were found in our previous studies. In order to investigate the effect of these three substitutions, a comparative analysis was conducted between the wild-type and mutant Top Is in Spodoptera exigua Hübner.

RESULTS

The optimal salt concentration of A653T and S729T was 150 mm, which is consistent with that of the wild-type Top I, while the mutant L530P showed maximum relaxation activity at a lower KCl concentration (100 mm). The mutated L530P and A653T Top Is showed a higher relaxation efficiency owing to an increased relaxation velocity toward the negatively supercoiled plasmid pBR322 DNA, which rendered L530P and A653T resistant to CPTs, whereas mutant S729T exhibited sensitivity to CPTs as a result of a decreased relaxation activity toward plasmid pBR322 DNA.

CONCLUSIONS

These results suggested that the polymorphism in Top I of insects was related to the biological activity of CPTs, which provided the basic information for reasonable usage of CPTs to control insect pests. © 2016 Society of Chemical Industry.

An In Silico Approach of Coumarin-Derived Inhibitors for Human DNA Topoisomerase I

Human topoisomerase I (Htopo I) is a vital target for anti-cancer agents; however, available anti-cancer agents are linked with several limitations. Therefore, designing novel inhibitors for Htopo I is significant. The rationale behind the current study is to identify novel coumarin inhibitors for Htopo I using in silico approaches and predict drug leads for in vitro studies. Using molecular docking and molecular dynamics, the binding affinities of 75 coumarins were compared with a known Htopo I inhibitor, topotecan. Docking studies predict three coumarins T1L25, T2L25, and T3L25 as most potent inhibitors for Htopo I. T2L25 gives the best grid score (–295 kJ mol–1), which is very comparable with that of topotecan (–302 kJ mol–1). The binding of these coumarins occurs preferentially via a planar geometry, and ligands bind at the binding site parallel to the axis of base pairing. NHCOCH3-substituted ligands are more favourable for binding when compared with the other substitute groups considered. The binding free energies calculated from molecular mechanics Poisson–Boltzmann surface area (MM-PBSA) method imply that T3L25 possesses the highest binding affinity when compared with the other two ligands. However, T1L25 and T2L25 have comparable binding free energies according to MM-PBSA calculations. Additionally, other calculated properties also support the suitability of these three derivatives as inhibitors for Htopo I. Therefore, the current study theoretically predicts three coumarin derivatives T1L25, T2L25, and T3L25 as potent inhibitors for Htopo I. These findings could lead to exploring novel non-camptothecin inhibitors for Htopo I.

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.

Sequence selectivity of the cleavage sites induced by topoisomerase I inhibitors: a molecular dynamics study

Topoisomerase IB (Top1) inhibitors, such as camptothecin (CPT), stabilize the Top1-DNA cleavage complex in a DNA sequence-dependent manner. The sequence selectivity of Top1 inhibitors is important for targeting specific genomic sequences of therapeutic value. However, the molecular mechanisms underlying this selectivity remain largely unknown. We performed molecular dynamics simulations to delineate structural, dynamic and energetic features that contribute to the differential sequence selectivity of the Top1 inhibitors. We found the sequence selectivity of CPT to be highly correlated with the drug binding energies, dynamic and structural properties of the linker domain. Chemical insights, gained by per-residue binding energy analysis revealed that the non-polar interaction between CPT and nucleotide at the +1 position of the cleavage site was the major (favorable) contributor to the total binding energy. Mechanistic insights gained by a potential of mean force analysis implicated that the drug dissociation step was associated with the sequence selectivity. Pharmaceutical insights gained by our molecular dynamics analyses explained why LMP-776, an indenoisoquinoline derivative under clinical development at the National Institutes of Health, displays different sequence selectivity when compared with camptothecin and its clinical derivatives.

Replacement of the human topoisomerase linker domain with the plasmodial counterpart renders the enzyme camptothecin resistant

A human/plasmodial hybrid enzyme, generated by swapping the human topoisomerase IB linker domain with the corresponding domain of the Plasmodium falciparum enzyme, has been produced and characterized. The hybrid enzyme displays a relaxation activity comparable to the human enzyme, but it is characterized by a much faster religation rate. The hybrid enzyme is also camptothecin resistant. A 3D structure of the hybrid enzyme has been built and its structural-dynamical properties have been analyzed by molecular dynamics simulation. The analysis indicates that the swapped plasmodial linker samples a conformational space much larger than the corresponding domain in the human enzyme. The large linker conformational variability is then linked to important functional properties such as an increased religation rate and a low drug reactivity, demonstrating that the linker domain has a crucial role in the modulation of the topoisomerase IB activity.

A kinetic clutch governs religation by type IB topoisomerases and determines camptothecin sensitivity

Type IB topoisomerases (Top1Bs) relax excessive DNA supercoiling associated with replication and transcription by catalyzing a transient nick in one strand to permit controlled rotation of the DNA about the intact strand. The natural compound camptothecin (CPT) and the cancer chemotherapeutics derived from it, irinotecan and topotecan, are highly specific inhibitors of human nuclear Top1B (nTop1). Previous work on vaccinia Top1B led to an elegant model that describes a straightforward dependence of rotation and religation on the torque caused by supercoiling. Here, we used a single-molecule DNA supercoil relaxation assay to measure the torque dependence of nTop1 and its inhibition by CPT. For comparison, we also examined mitochondrial Top1B and an N-terminal deletion mutant of nTop1. Despite substantial sequence homology in their core domains, nTop1 and mitochondrial Top1B exhibit dramatic differences in sensitivity to torque and CPT, with the N-terminal deletion mutant of nTop1 showing intermediate characteristics. In particular, nTop1 displays nearly torque-independent religation probability, distinguishing it from other Top1B enzymes studied to date. Kinetic modeling reveals a hitherto unobserved torque-independent transition linking the DNA rotation and religation phases of the enzymatic cycle. The parameters of this transition determine the torque sensitivity of religation and the efficiency of CPT binding. This "kinetic clutch" mechanism explains the molecular basis of CPT sensitivity and more generally provides a framework with which to interpret Top1B activity and inhibition.

Activities of topoisomerase I in its complex with SRSF1

Human DNA topoisomerase I (topo I) catalyzes DNA relaxation and phosphorylates SRSF1. Whereas the structure of topo I complexed with DNA has been resolved, the structure of topo I in the complex with SRSF1 and structural determinants of topo I activities in this complex are not known. The main obstacle to resolving the structure is a contribution of unfolded domains of topo I and SRSF1 in formation of the complex. To overcome this difficulty, we employed a three-step strategy: identifying the interaction regions, modeling the complex, and validating the model with biochemical methods. The binding sites in both topo I and SRSF1 are localized in the structured regions as well as in the unfolded domains. One observes cooperation between the binding sites in topo I but not in SRSF1. Our results indicate two features of the unfolded RS domain of SRSF1 containing phosphorylated residues that are critical for the kinase activity of topo I: its spatial arrangement relative to topo I and the organization of its sequence. The efficiency of phosphorylation of SRSF1 depends on the length and flexibility of the spacer between the two RRM domains that uniquely determine an arrangement of the RS domain relative to topo I. The spacer also influences inhibition of DNA nicking, a prerequisite for DNA relaxation. To be phosphorylated, the RS domain has to include a short sequence recognized by topo I. A lack of this sequence in the mutants of SRSF1 or its spatial inaccessibility in SRSF9 makes them inadequate as topo I/kinase substrates.

Structural and dynamical effects induced by the anticancer drug topotecan on the human topoisomerase I - DNA complex

BACKGROUND

Human topoisomerase I catalyzes the relaxation of DNA supercoils in fundamental cell processes like transcription, replication and chromosomal segregation. It is the only target of the camptothecin family of anticancer drugs. Among these, topotecan has been used to treat lung and ovarian carcinoma for several years. Camptothecins reversibly binds to the covalent intermediate DNA-enzyme, stabilizing the cleavable complex and reducing the religation rate. The stalled complex then collides with the progression of the replication fork, producing lethal double strand DNA breaks and eventually cell death.

METHODOLOGY/PRINCIPAL FINDINGS

Long lasting molecular dynamics simulations of the DNA-topoisomerase I binary complex and of the DNA-topoisomerase-topotecan ternary complex have been performed and compared. The conformational space sampled by the binary complex is reduced by the presence of the drug, as observed by principal component and cluster analyses. This conformational restraint is mainly due to the reduced flexibility of residues 633-643 (the region connecting the linker to the core domain) that causes an overall mobility loss in the ternary complex linker domain. During the simulation, DNA/drug stacking interactions are fully maintained, and hydrogen bonds are maintained with the enzyme. Topotecan keeps the catalytic residue Lys532 far from the DNA, making it unable to participate to the religation reaction. Arg364 is observed to interact with both the B and E rings of topotecan with two stable direct hydrogen bonds. An interesting constrain exerted by the protein on the geometrical arrangement of topotecan is also observed.

CONCLUSIONS/SIGNIFICANCE

Atomistic-scale understanding of topotecan interactions with the DNA-enzyme complex is fundamental to the explaining of its poisonous effect and of the drug resistance observed in several single residue topoisomerase mutants. We observed significant alterations due to topotecan in both short-range interactions and long-range protein domain communications.

Evidence of the crucial role of the linker domain on the catalytic activity of human topoisomerase I by experimental and simulative characterization of the Lys681Ala mutant

The functional and structural-dynamical properties of the Lys681Ala mutation in the human topoisomerase IB linker domain have been investigated by catalytic assays and molecular dynamics simulation. The mutant is characterized by a comparable cleavage and a strongly reduced religation rate when compared to the wild type protein. The mutant also displays perturbed linker dynamics, as shown by analysis of the principal components of the motion, and a reduced electrostatic interaction with DNA. Inspection of the inter atomic distances in proximity of the active site shows that in the mutant the distance between the amino group of Lys532 side chain and the 5′ OH of the scissile phosphate is longer than the wild type enzyme, providing an atomic explanation for the reduced religation rate of the mutant. Taken together these results indicate the existence of a long range communication between the linker domain and the active site region and points out the crucial role of the linker in the modulation of the catalytic activity.

UV-vis spectra of the anticancer camptothecin family drugs in aqueous solution: specific spectroscopic signatures unraveled by a combined computational and experimental study

The ultraviolet-visible absorption spectrum of camptothecin (CPT) has been been recorded in aqueous solution at pH 5.3, where the equilibrium among the different CPT forms is shifted toward the lactonic one. Time-dependent density functional theory (TD-DFT) computations lead to a remarkable reproduction of the experimental spectrum only upon addition of explicit water molecules in interaction with specific moieties of the camptothecin molecule. Molecular dynamics (MD) simulations enforcing boundary periodic conditions for CPT embedded with 865 water molecules, with a force field derived from DFT computations, show that the experimental spectrum is due to the contributions of CPT molecules with different solvation patterns. A similar solvent effect is observed for several CPT derivatives, including the clinically relevant SN-38 and topotecan drugs. The quantitative agreement between TD-DFT/MD computations and experimental data allow us to identify specific spectroscopic signatures diagnostic of the drug environment and to develop procedures that can be used to monitor the drug-DNA/protein interaction.

DNA topoisomerases: harnessing and constraining energy to govern chromosome topology

DNA topoisomerases are a diverse set of essential enzymes responsible for maintaining chromosomes in an appropriate topological state. Although they vary considerably in structure and mechanism, the partnership between topoisomerases and DNA has engendered commonalities in how these enzymes engage nucleic acid substrates and control DNA strand manipulations. All topoisomerases can harness the free energy stored in supercoiled DNA to drive their reactions; some further use the energy of ATP to alter the topology of DNA away from an enzyme-free equilibrium ground state. In the cell, topoisomerases regulate DNA supercoiling and unlink tangled nucleic acid strands to actively maintain chromosomes in a topological state commensurate with particular replicative and transcriptional needs. To carry out these reactions, topoisomerases rely on dynamic macromolecular contacts that alternate between associated and dissociated states throughout the catalytic cycle. In this review, we describe how structural and biochemical studies have furthered our understanding of DNA topoisomerases, with an emphasis on how these complex molecular machines use interfacial interactions to harness and constrain the energy required to manage DNA topology.

Thr729 in human topoisomerase I modulates anti-cancer drug resistance by altering protein domain communications as suggested by molecular dynamics simulations

The role of Thr729 in modulating the enzymatic function of human topoisomerase I has been characterized by molecular dynamics (MD) simulation. In detail, the structural–dynamical behaviour of the Thr729Lys and the Thr729Pro mutants have been characterized because of their in vivo and in vitro functional properties evidenced in the accompanying paper. Both mutants can bind to the DNA substrate and are enzymatically active, but while Thr729Lys is resistant even at high concentration of the camptothecin (CPT) anti-cancer drug, Thr729Pro shows only a mild reduction in drug sensitivity and in DNA binding. MD simulations show that the Thr729Lys mutation provokes a structural perturbation of the CPT-binding pocket. On the other hand, the Thr729Pro mutant maintains the wild-type structural scaffold, only increasing its rigidity. The simulations also show the complete abolishment, in the Thr729Lys mutant, of the protein communications between the C-terminal domain (where the active Tyr723 is located) and the linker domain, that plays an essential role in the control of the DNA rotation, thus explaining the distributive mode of action displayed by this mutant.

The open state of human topoisomerase I as probed by molecular dynamics simulation

The open state of human topoisomerase I has been probed by molecular dynamics simulation, starting from the coordinates of the closed structure of the protein complexed with DNA, after elimination of the 22-bp DNA duplex oligonucleotide. A repulsion force between the two lips of the protein has been introduced for a short time to induce destabilization of the local minimum, after which an unperturbed simulation has been carried out for 10 ns. The simulation shows that the protein undergoes a large conformational change due to rearrangements in the orientation of the protein domains, which however move as a coherent unit, fully maintaining their secondary and tertiary structures. Despite movements between the domains as large as 80-90 A, the catalytic pentad remains preassembled, the largest deviation of the active site backbone atoms from the starting crystallographic structure being only 1.7 A. Electrostatic calculation of the open protein structure shows that the protein displays a vast positive region with the active site residues located nearly at its center, in a conformation perfectly suited to interact with the negatively charged supercoiled DNA substrate.

Alterations in linker flexibility suppress DNA topoisomerase I mutant-induced cell lethality

Eukaryotic DNA topoisomerase I (Top1p) catalyzes changes in DNA topology via the formation of a covalent enzyme-DNA intermediate, which is reversibly stabilized by the anticancer agent camptothecin (CPT). Crystallographic studies of the 70-kDa C terminus of human Top1p bound to duplex DNA describe a monomeric protein clamp circumscribing the DNA helix. The structures, which lack the N-terminal domain, comprise the conserved clamp, an extended linker domain, and the conserved C-terminal active site Tyr domain. CPT bound to the covalent Top1p-DNA complex limits linker flexibility, allowing structural determination of this domain. We previously reported that mutation of Ala(653) to Pro in the linker increases the rate of enzyme-catalyzed DNA religation, thereby rendering Top1A653Pp resistant to CPT (Fiorani, P., Bruselles, A., Falconi, M., Chillemi, G., Desideri, A., and Benedetti P. (2003) J. Biol. Chem. 278, 43268-43275). Molecular dynamics studies suggested mutation-induced increases in linker flexibility alter Top1p catalyzed DNA religation. To address the functional consequences of linker flexibility on enzyme catalysis and drug sensitivity, we investigated the interactions of the A653P linker mutation with a self-poisoning T718A mutation within the active site of Top1p. The A653P mutation suppressed the lethal phenotype of Top1T718Ap in yeast, yet did not restore enzyme sensitivity to CPT. However, the specific activity of the double mutant was decreased in vivo and in vitro, consistent with a decrease in DNA binding. These findings support a model where changes in the flexibility or orientation of the linker alter the geometry of the active site and thereby the kinetics of DNA cleavage/religation catalyzed by Top1p.

New Hints on the pH-Driven Tautomeric Equilibria of the Topotecan Anticancer Drug in Aqueous Solutions from an Integrated Spectroscopic and Quantum-Mechanical Approach

The equilibria between the different forms of the topotecan anticancer drug have been studied at moderately acidic and physiological pH by an integrated computational tool rooted in the density functional theory and its time-dependent extension together with the polarizable continuum model. The results allow an unbiased selection between the different possible tautomeric forms and provide invaluable complements to experimental data. The ultraviolet-visible topotecan spectrum, recorded at moderately acidic pH, is accurately reproduced only by TD-DFT computations including solvent effects. Comparison of the experimental and calculated bands of the UV-vis spectrum at physiological pH indicates the presence of an equilibrium among different forms that is tuned by the microenvironment embedding the drug. The quantitative agreement between TD-DFT/PCM computations and experiments allows the identification of unequivocal spectroscopic signatures for different forms of topotecan.

Effect on DNA relaxation of the single Thr718Ala mutation in human topoisomerase I: a functional and molecular dynamics study

The functional and dynamical properties of the human topoisomerase I Thr718Ala mutant have been compared to that of the wild-type enzyme using functional assays and molecular dynamics (MD) simulations. At physiological ionic strength, the cleavage and religation rates, evaluated on oligonucleotides containing the preferred topoisomerase I DNA sequence, are almost identical for the wild-type and the mutated enzymes, as is the cleavage/religation equilibrium. On the other hand, the Thr718Ala mutant shows a decreased efficiency in a DNA plasmid relaxation assay. The MD simulation, carried out on the enzyme complexed with its preferred DNA substrate, indicates that the mutant has a different dynamic behavior compared to the wild-type enzyme. Interestingly, no changes are observed in the proximity of the mutation site, whilst a different flexibility is detected in regions contacting the DNA scissile strand, such as the linker and the V-shaped α helices. Taken together, the functional and simulation results indicate a direct communication between the mutation site and regions located relatively far away, such as the linker domain, that with their altered flexibility confer a reduced DNA relaxation efficiency. These results provide evidence that the comprehension of the topoisomerase I dynamical properties are an important element in the understanding of its complex catalytic cycle.

Human DNA topoisomerase I: relaxation, roles, and damage control

Human DNA topoisomerase I is an essential enzyme involved in resolving the torsional stress associated with DNA replication, transcription, and chromatin condensation. The catalytic cycle of the enzyme consists of DNA cleavage to form a covalent enzyme-DNA intermediate, DNA relaxation, and finally, re-ligation of the phosphate backbone to restore the continuity of the DNA. Structure/function studies have elucidated a flexible enzyme that relaxes DNA through coordinated, controlled movements of distinct enzyme domains. The cellular roles of topoisomerase I are apparent throughout the nucleus, but the concentration of processes acting on ribosomal DNA results in topoisomerase I accumulation in the nucleolus. Although the activity of topoisomerase I is required in these processes, the enzyme can also have a deleterious effect on cells. In the event that the final re-ligation step of the reaction cycle is prevented, the covalent topoisomerase I-DNA intermediate becomes a toxic DNA lesion that must be repaired. The complexities of the relaxation reaction, the cellular roles, and the pathways that must exist to repair topoisomerase I-mediated DNA damage highlight the importance of continued study of this essential enzyme.