The stability of a model substrate for topoisomerase 1-mediated DNA religation depends on the presence of mismatched base pairs

Entrapment of DNA topoisomerase-DNA complexes by nucleotide/nucleoside analogs

Topoisomerases are well-validated targets for cancer chemotherapy and DNA topoisomerase 1 (Top1) is the sole target of the camptothecin (CPT) class of anticancer drugs. Over the last 20 years, multiple studies have shown Top1 activity is modulated by non-native DNA structures and this can lead to trapping of Top1 cleavage complexes (Top1cc) and conversion to DNA double strand breaks. Among the perturbations to DNA structure that generate Top1cc are nucleoside analogs that are incorporated into genomic DNA during replication including cytarabine, gemcitabine, and 5-fluoro-2'-deoxyuridine (FdU). We review the literature summarizing the role of Top1cc in mediating the DNA damaging and cytotoxic activities of nucleoside analogs. We also summarize studies demonstrating distinct differences between Top1cc induced by nucleoside analogs and CPTs, particularly with regard to DNA repair. Collectively, these studies demonstrate that, while Top1 is a common target for both Top1 poisons such as CPT and nucleoside analogs such as FdU, these agents are not redundant. In recent years, studies have shown that Top1 poisons and nucleoside analogs together with other anti-cancer drugs such as cisplatin cause replication stress and the DNA repair pathways that modulate the cytotoxic activities of these compounds are being elucidated. We present an overview of this evolving literature, which has implications for how targeting of Top1 with nucleoside analogs can be used more effectively for cancer treatment.

Using correlated motions to determine sufficient sampling times for molecular dynamics

Here we present a time-dependent correlation method that provides insight into how long a system takes to grow into its equal-time (Pearson) correlation. We also show a usage of an extant time-lagged correlation method that indicates the time for parts of a system to become decorrelated, relative to equal-time correlation. Given a completed simulation (or set of simulations), these tools estimate (i) how long of a simulation of the same system would be sufficient to observe the same correlated motions, (ii) if patterns of observed correlated motions indicate events beyond the timescale of the simulation, and (iii) how long of a simulation is needed to observe these longer timescale events. We view this method as a decision-support tool that will aid researchers in determining necessary sampling times. In principle, this tool is extendable to any multidimensional time series data with a notion of correlated fluctuations; however, here we limit our discussion to data from molecular-dynamics simulations.

Visualizing correlated motion with HDBSCAN clustering

Correlated motion analysis provides a method for understanding communication between and dynamic similarities of biopolymer residues and domains. The typical equal-time correlation matrices-frequently visualized with pseudo-colorings or heat maps-quickly convey large regions of highly correlated motion but hide more subtle similarities of motion. Here we propose a complementary method for visualizing correlations within proteins (or general biopolymers) that quickly conveys intuition about which residues have a similar dynamic behavior. For grouping residues, we use the recently developed non-parametric clustering algorithm HDBSCAN. Although the method we propose here can be used to group residues using correlation as a similarity matrix-the most straightforward and intuitive method-it can also be used to more generally determine groups of residues which have similar dynamic properties. We term these latter groups "Dynamic Domains", as they are based not on spatial closeness but rather closeness in the column space of a correlation matrix. We provide examples of this method across three human proteins of varying size and function-the Nf-Kappa-Beta essential modulator, the clotting promoter Thrombin and the mismatch repair protein (dimer) complex MutS-alpha. Although the examples presented here are from all-atom molecular dynamics simulations, this visualization technique can also be used on correlations matrices built from any ensembles of conformations from experiment or computation.

All-atom MD indicates ion-dependent behavior of therapeutic DNA polymer

Understanding the efficacy of and creating delivery mechanisms for therapeutic nucleic acids requires understanding structural and kinetic properties which allow these polymers to promote the death of cancerous cells. One molecule of interest is a 10 mer of FdUMP (5-fluoro-2'-deoxyuridine-5'-O-monophosphate) - also called F10. Here we investigate the structural and kinetic behavior of F10 in intracellular and extracellular solvent conditions along with non-biological conditions that may be efficacious in in vitro preparations of F10 delivery systems. From our all-atom molecular dynamics simulations totaling 80 microseconds, we predict that F10's phosphate groups form close-range interactions with calcium and zinc ions, with calcium having the highest affinity of the five ions investigated. We also predict that F10's interactions with magnesium, potassium and sodium are almost exclusively long-range interactions. In terms of intramolecular interactions, we find that F10 is least structured (in terms of hydrogen bonds among bases) in the 150 mM NaCl (extracellular-like solvent conditions) and most structured in 150 mM ZnCl2. Kinetically, we see that F10 is unstable in the presence of magnesium, sodium or potassium, finding stable kinetic traps in the presence of calcium or zinc.

All-Atom MD Predicts Magnesium-Induced Hairpin in Chemically Perturbed RNA Analog of F10 Therapeutic

Given their increasingly frequent usage, understanding the chemical and structural properties which allow therapeutic nucleic acids to promote the death of cancer cells is critical for medical advancement. One molecule of interest is a 10-mer of FdUMP (5-fluoro-2'-deoxyuridine-5'-O-monophosphate) also called F10. To investigate causes of structural stability, we have computationally restored the 2' oxygen on each ribose sugar of the phosphodiester backbone, creating FUMP[10]. Microsecond time-scale, all-atom, simulations of FUMP[10] in the presence of 150 mM MgCl₂ predict that the strand has a 45% probability of folding into a stable hairpin-like secondary structure. Analysis of 16 μs of data reveals phosphate interactions as likely contributors to the stability of this folded state. Comparison with polydT and polyU simulations predicts that FUMP[10]'s lowest order structures last for one to 2 orders of magnitude longer than similar nucleic acid strands. Here we provide a brief structural and conformational analysis of the predicted structures of FUMP[10], and suggest insights into its stability via comparison to F10, polydT, and polyU.

MutSα's Multi-Domain Allosteric Response to Three DNA Damage Types Revealed by Machine Learning

MutSα is a key component in the mismatch repair (MMR) pathway. This protein is responsible for initiating the signaling pathways for DNA repair or cell death. Herein we investigate this heterodimer's post-recognition, post-binding response to three types of DNA damage involving cytotoxic, anti-cancer agents-carboplatin, cisplatin, and FdU. Through a combination of supervised and unsupervised machine learning techniques along with more traditional structural and kinetic analysis applied to all-atom molecular dynamics (MD) calculations, we predict that MutSα has a distinct response to each of the three damage types. Via a binary classification tree (a supervised machine learning technique), we identify key hydrogen bond motifs unique to each type of damage and suggest residues for experimental mutation studies. Through a combination of a recently developed clustering (unsupervised learning) algorithm, RMSF calculations, PCA, and correlated motions we predict that each type of damage causes MutSα to explore a specific region of conformation space. Detailed analysis suggests a short range effect for carboplatin-primarily altering the structures and kinetics of residues within 10 angstroms of the damaged DNA-and distinct longer-range effects for cisplatin and FdU. In our simulations, we also observe that a key phenylalanine residue-known to stack with a mismatched or unmatched bases in MMR-stacks with the base complementary to the damaged base in 88.61% of MD frames containing carboplatinated DNA. Similarly, this Phe71 stacks with the base complementary to damage in 91.73% of frames with cisplatinated DNA. This residue, however, stacks with the damaged base itself in 62.18% of trajectory frames with FdU-substituted DNA and has no stacking interaction at all in 30.72% of these frames. Each drug investigated here induces a unique perturbation in the MutSα complex, indicating the possibility of a distinct signaling event and specific repair or death pathway (or set of pathways) for a given type of damage.

Uncovering Large-Scale Conformational Change in Molecular Dynamics without Prior Knowledge

As the length of molecular dynamics (MD) trajectories grows with increasing computational power, so does the importance of clustering methods for partitioning trajectories into conformational bins. Of the methods available, the vast majority require users to either have some a priori knowledge about the system to be clustered or to tune clustering parameters through trial and error. Here we present non-parametric uses of two modern clustering techniques suitable for first-pass investigation of an MD trajectory. Being non-parametric, these methods require neither prior knowledge nor parameter tuning. The first method, HDBSCAN, is fast-relative to other popular clustering methods-and is able to group unstructured or intrinsically disordered systems (such as intrinsically disordered proteins, or IDPs) into bins that represent global conformational shifts. HDBSCAN is also useful for determining the overall stability of a system-as it tends to group stable systems into one or two bins-and identifying transition events between metastable states. The second method, iMWK-Means, with explicit rescaling followed by K-Means, while slower than HDBSCAN, performs well with stable, structured systems such as folded proteins and is able to identify higher resolution details such as changes in relative position of secondary structural elements. Used in conjunction, these clustering methods allow a user to discern quickly and without prior knowledge the stability of a simulated system and identify both local and global conformational changes.

All-Atom Molecular Dynamics Reveals Mechanism of Zinc Complexation with Therapeutic F10

Advancing the use of therapeutic nucleic acids requires understanding the chemical and structural properties that allow these polymers to promote the death of malignant cells. Here we explore Zn²⁺ complexation by the fluoropyrimidine polymer F10, which has strong activities in multiple preclinical models of cancer. Delivery of fluoropyrimidine FdUMP in the 10-residue polymer F10 rather than the nucleobase (5-fluorouracil) allows consideration of metal ion binding effects on drug delivery. The differences in metal ion interactions with fluoropyrimidine compared to normal DNA results in conformation changes that affect protein binding, cell uptake, and codelivery of metals such as zinc, and the cytoxicity thereof. Microsecond-time-scale, all-atom simulations of F10 predict that zinc selectively stabilizes the polymer via interactions with backbone phosphate groups and suggest a mechanism of complexation for the zinc-base interactions shown in previous experimental work. The positive zinc ions are attracted to the negatively charged phosphate groups. Once the Zn²⁺ ions are near F10, they cause the base's N3 nitrogen to deprotonate. Subsequently, magnesium atoms displace zinc from their interactions with phosphate, freeing the zinc ions to interact with the FdU bases by forming weak interactions with the O4 oxygen and the fluorine attached to C5. These interactions of magnesium with phosphate groups and zinc with nucleobases agree with previous experimental results and are seen in MD simulations only when magnesium is introduced after N3 deprotonation, indicating a specific order of metal binding events. Additionally, we predict interactions between zinc and F10's O2 atoms, which were not previously observed. By comparison to 10mers of polyU and polydT, we also predict that the presence of fluorine increases the binding affinity of zinc to F10 relative to analogous strands of RNA and DNA consisting of only native nucleotides.

Replication-dependent irreversible topoisomerase 1 poisoning is responsible for FdUMP[10] anti-leukemic activity

Previous studies have indicated that 5-Fluoro-2'-deoxyuridine-5'-O-monophosphate 10mer (FdUMP[10]) displays strong antileukemic activity through the dual targeting of thymidylate synthase (TS) and DNA topoisomerase 1 (Top1). The present studies were undertaken to clarify the relationship between the induction of a thymineless state and the formation of Top1 cleavage complexes (Top1CC) for inducing cell death and to clarify the role of DNA replication for induction of lethal DNA double-strand breaks (DSBs) in FdUMP[10]-treated acute myeloid leukemia (AML) cells. Human promyelocytic (HL60) and AML (KG1a, Molm13, THP-1) cells were synchronized by serum starvation and treated with FdUMP[10] with thymidine (Thy) rescue. Cells were assayed for TS inhibition, DNA DSBs, Top1CC, and apoptosis to clarify the interrelationship of TS inhibition and Top1CC for cell death. FdUMP[10] induced a thymineless state in AML cells and exogenous Thy administered within the first 18 hours of treatment rescued FdUMP[10]-induced Top1CC formation, γH2AX phosphorylation, and apoptosis induction. Exogenous Thy was not effective after cells had committed to mitosis and undergone cell division in the presence of FdUMP[10]. FdUMP[10] treatment resulted in Chk1 activation, and Chk1 inhibition enhanced FdUMP[10]-induced DNA damage and apoptosis. Jnk-signaling was required for FdUMP[10]-induced apoptosis in promyelocytic HL60 cells and in THP1 cells, but was antiapoptotic in Molm13 and to a lesser extent KG1a AML cells. The results are consistent with FdUMP[10] inducing a thymineless state, leading to misincorporation of FdU into genomic DNA of proliferating cells. Top1CC form in cells upon re-entry into S-phase, resulting in DNA double-strand breaks, and initiating apoptotic signaling that can be either muted or enhanced by Jnk-signaling depending on cell type.

Cooperative stabilization of Zn(2+):DNA complexes through netropsin binding in the minor groove of FdU-substituted DNA

The simultaneous binding of netropsin in the minor groove and Zn(2+) in the major groove of a DNA hairpin that includes 10 consecutive FdU nucleotides at the 3'-terminus (3'FdU) was demonstrated based upon NMR spectroscopy, circular dichroism (CD), and computational modeling studies. The resulting Zn(2+)/netropsin: 3'FdU complex had very high thermal stability with aspects of the complex intact at 85 °C, conditions that result in complete dissociation of Mg(2+) complexes. CD and (19)F NMR spectroscopy were consistent with Zn(2+) binding in the major groove of the DNA duplex and utilizing F5 and O4 of consecutive FdU nucleotides as ligands with FdU nucleotides hemi-deprotonated in the complex. Netropsin is bound in the minor groove of the DNA duplex based upon 2D NOESY data demonstrating contacts between AH2 (1)H and netropsin (1)H resonances. The Zn(2+)/netropsin: 3'FdU complex displayed increased cytotoxicity towards PC3 prostate cancer (PCa) cells relative to the constituent components or separate complexes (e.g. Zn(2+):3'FdU) indicating that this new structural motif may be therapeutically useful for PCa treatment. An animated interactive 3D complement (I3DC) is available in Proteopedia at http://proteopedia.org/w/Journal:JBSD:32.