Sliding of alkylating anticancer drugs along the minor groove of DNA: new insights on sequence selectivity

Currently, little is known about the molecular recognition pathways between DNA-alkylating anticancer drugs and their targets despite their pharmacological relevance. In the framework of classical molecular dynamics simulations, here we use umbrella sampling to map the potential of mean force (PMF) associated with sliding along the DNA minor groove of two of these compounds. These are an indole derivative of duocarmycin (DSI) and the putative reactive form of anthramycin (anhydro-anthramycin, IMI). Twenty-three configurations were considered for each drug/DNA complex, corresponding to a movement along ∼3 basepairs. The alkylation site turns out to be the most favorable for DSI, while a barrier of ∼6 kcal/mol separates the reactive configuration of IMI·DNA from the absolute minimum. An analysis of various contributions to the PMF reveals that solvent effects play an important role for the largest and more flexible drug DSI. Instead, the PMF of IMI·DNA overall correlates with changes in the binding enthalpy. Implications of these results on the sequence selectivity of the two drugs are discussed.     

基于亲和探针的药物靶点鉴定技术研究进展

药物靶点的鉴定和相关研究在药学研究领域具有重要的理论指导意义和实用价值。利用亲和探针偶联靶分子的方法是目前发现药 物靶点的主要手段之一。该方法可从分子水平发现药物的作用靶点,从而对药物的分子作用机制提供细胞水平的直接证据。从 DNA 和小 分子药物探针的构筑和应用入手,对近些年鉴定 DNA 损伤识别蛋白的研究进展进行了较为详尽的讨论,并简要介绍目前探索小分子药物 作用靶点的主流技术。作为亲和偶联鉴定药物作用靶点方法的重要组成部分,亲和探针设计的合理性关系到方法本身的可操作性以及鉴 定结果的可靠性。从多个角度对 DNA 探针和小分子药物探针的设计经验进行了较为系统的总结,例如经典的亲和纯化分离方法,以及更 为高效的光激发共价偶联技术等。这些方法和思路为探索 DNA 损伤相关蛋白质的功能以及小分子药物的细胞作用机制提供了丰富的研究 工具,有助于从分子水平理解药物的作用机制。     

Binding of an Indenoisoquinoline to the Topoisomerase-DNA Complex Induces Reduction of Linker Mobility and Strengthening of Protein-DNA Interaction

Long-duration comparative molecular dynamics simulations of the DNA-topoisomerase binary and DNA-topoisomerase-indenoisoquinoline ternary complexes have been carried out. The analyses demonstrated the role of the drug in conformationally stabilizing the protein-DNA interaction. In detail, the protein lips, clamping the DNA substrate, interact more tightly in the ternary complex than in the binary one. The drug also reduces the conformational space sampled by the protein linker domain through an increased interaction with the helix bundle proximal to the active site. A similar alteration of linker domain dynamics has been observed in a precedent work for topotecan but the molecular mechanisms were different if compared to those described in this work. Finally, the indenoisoquinoline keeps Lys532 far from the DNA, making it unable to participate in the religation reaction, indicating that both short- and long-range interactions contribute to the drug poisoning effect.     

Structural and energetic insights into sequence-specific interaction in DNA–drug recognition: development of affinity predictor and analysis of binding selectivity

Although the molecular mechanism and thermodynamic profile of a wide variety of chemical agents have been examined intensively in the past decades in terms of specific recognition of their protein receptors, to date the physicochemical nature of DNA–drug recognition and association still remains largely unexplored. The present study focused on understanding the structural basis, energetic landscape, and biological implications underlying the binding of small-molecule ligands to their cognate or non-cognate DNA receptors. First, a new method to capture the structural features of DNA–drug complex architecture was proposed and then used to correlate the extracted features with binding affinity of the complexes. By employing this method, a statistical regression-based predictor was developed to quantitatively evaluate the interaction potency of drug compounds with DNA in a fast and reliable manner. Subsequently, we use the predictor to examine the binding behavior of a number of structure-available, affinity-known DNA–drug complexes as well as a large pool of randomly generated DNA decoys in complex with the same drugs. It was found that (1) as compared with protein–DNA recognition, small-molecule agents have relatively low specificity in selecting their cognate DNA targets from the background of numerous random decoys; (2) the abundance of A–T base pairs in the DNA core motif exhibits a significant positive correlation with the affinity of drug ligand binding to the DNA receptor; and (3) high affinity seems not to be closely related to high selectivity for a DNA-targeting drug, although high-affinity drug entities have a greater possibility of being ranked computationally as top binders. We hope that this work will provide a preliminary insight into the molecular origin of sequence-specific interactions in DNA–drug recognition.

Combinatorial determination of sequence specificity for nanomolar DNA-binding hairpin polyamides

Development of sequence-specific DNA-binding drugs is an important pharmacological goal, given the fact that numerous existing DNA-directed chemotherapeutic drugs rely on the strength and selectivity of their DNA interactions for therapeutic activity. Among the DNA-binding antibiotics, hairpin polyamides represent the only class of small molecules that can practically bind any predetermined DNA sequence. DNA recognition by these ligands depends on their side-by-side amino acid pairings in the DNA minor groove. Extensive studies have revealed that these molecules show extremely high affinity for sequence-directed, minor groove interaction. However, the specificity of such interactions in the presence of a large selection of sequences such as the human genome is not known. We used the combinatorial selection method restriction endonuclease protection, selection, and amplification (REPSA) to determine the DNA binding specificity of two hairpin polyamides, ImPyPyPy-gamma-PyPyPyPy-beta-Dp and ImPyPyPy-gamma-ImPyPyPy-beta-Dp, in the presence of more than 134 million different sequences. These were verified by restriction endonuclease protection assays and DNase I footprinting analysis. Our data showed that both hairpin polyamides preferentially selected DNA sequences having consensus recognition sites as defined by the Dervan pairing rules. These consensus sequences were rather degenerate, as expected, given that the stacked pyrrole-pyrrole amino acid pairs present in both polyamides are unable to discriminate between A.T and T.A base pairs. However, no individual sequence within these degenerate consensus sequences was preferentially selected by REPSA, indicating that these hairpin polyamides are truly consensus-specific DNA-binding ligands. We also discovered a preference for overlapping consensus binding sites among the sequences selected by the hairpin polyamide ImPyPyPy-gamma-PyPyPyPy-beta-Dp, and confirmed by DNase I footprinting that these complex sites provide higher binding affinity. These data suggest that multiple hairpin polyamides can cooperatively bind to their highest-affinity sites.     

Base mutation analysis of topoisomerase II-idarubicin-DNA ternary complex formation. Evidence for enzyme subunit cooperativity in DNA cleavage.

Antitumor drugs, such as anthracyclines, interfere with mammalian DNA topoisomerase II by forming a ternary complex, DNA-drug-enzyme, in which DNA strands are cleaved and covalently linked to the enzyme. In this work, a synthetic 36-bp DNA oligomer derived from SV40 and mutated variants were used to determine the effects of base mutations on DNA cleavage levels produced by murine topoisomerase II with and without idarubicin. Although site competition could affect cleavage levels, mutation effects were rather similar among several cleavage sites. The major sequence determinants of topoisomerase II DNA cleavage without drugs are up to five base pairs apart from the strand cut, suggesting that DNA protein contacts involving these bases are particularly critical for DNA site recognition. Cleavage sites with adenines at positions -1 were detected without idarubicin only under conditions favouring enzyme binding to DNA, showing that these sites are low affinity sites for topoisomerase II DNA cleavage and/or binding. Moreover, the results indicated that the sequence 5'-(A)TA/(A)-3' (the slash indicates the cleaved bond, parenthesis indicate conditioned preference) from -3 to +1 positions constitutes the complete base sequence preferred by anthracyclines. An important finding was that mutations that improve the fit to the above consensus on one strand can also increase cleavage on the opposite strand, suggesting that a drug molecule may effectively interact with one enzyme subunit only and trap the whole dimeric enzyme. These findings documented that DNA recognition by topoisomerase II may occur at one or the other strand, and not necessarily at both of them, and that the two subunits can act cooperatively to cleave a double helix.     

Triple helix-forming oligonucleotides conjugated to indolocarbazole poisons direct topoisomerase I-mediated DNA cleavage to a specific site.

Topoisomerase I is an ubiquitous DNA-cleaving enzyme and an important therapeutic target in cancer chemotherapy for camptothecins as well as for indolocarbazole antibiotics such as rebeccamycin. To achieve a sequence-specific cleavage of DNA by topoisomerase I, a triple helix-forming oligonucleotide was covalently linked to indolocarbazole-type topoisomerase I poisons. The three indolocarbazole-oligonucleotide conjugates investigated were able to direct topoisomerase I cleavage at a specific site based upon sequence recognition by triplex formation. The efficacy of topoisomerase I-mediated DNA cleavage depends markedly on the intrinsic potency of the drug. We show that DNA cleavage depends also upon the length of the linker arm between the triplex-forming oligonucleotide and the drug. Based on a known structure of the DNA-topoisomerase I complex, a molecular model of the oligonucleotide conjugates bound to the DNA-topoisomerase I complex was elaborated to facilitate the design of a potent topoisomerase I inhibitor-oligonucleotide conjugate with an optimized linker between the two moieties. The resulting oligonucleotide-indolocarbazole conjugate at 10 nM induced cleavage at the triple helix site 2-fold more efficiently than 5 microM of free indolocarbazole, while the other drug-sensitive sites were not cleaved. The rational design of drug-oligonucleotide conjugates carrying a DNA topoisomerase poison may be exploited to improve the efficacy and selectivity of chemotherapeutic cancer treatments by targeting specific genes and reducing drug toxicity.     

DNA sequence recognition by bispyrazinonaphthalimides antitumor agents

Bifunctional DNA intercalating agents have long attracted considerable attention as anticancer agents. One of the lead compounds in this category is the dimeric antitumor drug elinafide, composed of two tricyclic naphthalimide chromophores separated by an aminoalkyl linker chain optimally designed to permit bisintercalation of the drug into DNA. In an effort to optimize the DNA recognition capacity, different series of elinafide analogues have been prepared by extending the surface of the planar drug chromophore which is important for DNA sequence recognition. We report here a detailed investigation of the DNA sequence preference of three tetracyclic monomeric or dimeric pyrazinonaphthalimide derivatives. Melting temperature measurements and surface plasmon resonance (SPR) studies indicate that the dimerization of the tetracyclic planar chromophore considerably augments the affinity of the drug for DNA, polynucleotides, or hairpin oligonucleotides and promotes selective interaction with G.C sites. The (CH(2))(2)NH(CH(2))(3)NH(CH(2))(2) connector stabilizes the drug-DNA complexes. The methylation of the two nitrogen atoms of this linker chain reduces the binding affinity and increases the dissociation rates of the drug-DNA complexes by a factor of 10. DNase I footprinting experiments were used to investigate the sequence selectivity of the drugs, demonstrating highly preferential binding to G.C-rich sequences. It also served to select a high-affinity site encompassing the sequence 5'-GACGGCCAG which was then introduced into a biotin-labeled hairpin oligonucleotide to accurately measure the binding parameters by SPR. The affinity constant of the unmethylated dimer for this sequence is 500 times higher than that of the monomer compound and approximately 10 times higher than that of the methylated dimer. The DNA groove accessibility was also probed with three related oligonucleotides carrying G --> c(7)G, G --> I, and C --> M substitutions. The level of drug binding to the two hairpin oligonucleotides containing 7-deazaguanine (c(7)G) or 5-methylcytosine (M) residues is unchanged or only slightly reduced compared to that of the unmodified target. In contrast, incorporation of inosine (I) residues considerably decreases the extent of drug binding or even abolishes the interaction as is the case with the monomer. The pyrazinonaphthalimide derivatives are thus much more sensitive to the deletion of the exocyclic guanine 2-amino group exposed in the minor groove of the duplex than to the modification of the major groove elements. The complementary SPR footprinting methodology combining site selection and quantitative DNA affinity analysis constitutes a reliable method for dissecting the DNA sequence selectivity profile of reversible DNA binding small molecules.     

Interaction of Antidepressants with the Serotonin and Norepinephrine Transporters: MUTATIONAL STUDIES OF THE S1 SUBSTRATE BINDING POCKET*

The serotonin transporter (SERT) and the norepinephrine transporter (NET) are sodium-dependent neurotransmitter transporters responsible for reuptake of released serotonin and norepinephrine, respectively, into nerve terminals in the brain. A wide range of inhibitors of SERT and NET are used as treatment of depression and anxiety disorders or as psychostimulant drugs of abuse. Despite their clinical importance, the molecular mechanisms by which various types of antidepressant drugs bind and inhibit SERT and NET are still elusive for the majority of the inhibitors, including the molecular basis for SERT/NET selectivity. Mutational analyses have suggested that a central substrate binding site (denoted the S1 pocket) also harbors an inhibitor binding site. In this study, we determine the effect of mutating six key S1 residues in human SERT (hSERT) and NET (hNET) on the potency of 15 prototypical SERT/NET inhibitors belonging to different drug classes. Analysis of the resulting drug sensitivity profiles provides novel information on drug binding modes in hSERT and hNET and identifies specific S1 residues as important molecular determinants for inhibitor potency and hSERT/hNET selectivity.     

Infrared linear dichroism studies of DNA-drug complexes: quantitative determination of the drug-induced restriction of the B-A transition.

The B-A transition of films or fibers of NaDNA occurs at a relative humidity of 75-85%. The fraction of DNA that changed the conformation from B to A form can be determined quantitatively by infrared linear dichroism. DNA-binding drugs can 'freeze' a fraction of DNA in the B form. This fraction of DNA is in the B form and cannot be converted to A-DNA even at a reduced relative humidity of 54%. The 'freezing' potentiality of various drugs can be described by the 'freezing' index, FI, expressed in base pairs per added drug. Drugs with a high value of FI (more than eight base pairs per drug) were observed among both intercalating and groove-binding drugs. High values of FI imply restriction of the conformational flexibility of DNA significantly going beyond the binding site of the drug. This long-range effect of drugs on the conformational flexibility of DNA may be connected with the molecular mechanism of drug action. The freezing index FI is a new quantitative parameter of drug-DNA interaction that should be considered as a valuable tool for drug design.     

Thermodynamics aspects of interactions between small ligands and DNA

DNA is a molecular target for many anticancer and antiviral drugs. Therefore, a clear understanding of the interaction of small molecules with DNA is important in the rational design of ligands that can bind to DNA with high affinity and selectivity. There are several methods to investigate interactions between drug and DNA. Some of them measures changing into DNA structures, such as lengthening and untwisting of helix of DNA. Other techniques measure changing in drug environment. With the increasing availability of sensitive microcalorimeters, particular interest has arisen in the thermodynamics of drug-DNA interaction. Using such methods permit direct determination of enthalpy changes associated with reactions. One experiment permits to obtain also binding constant, hence an almost complete thermodynamic profile can be established. This profile offers key insights into the molecular forces that drive complex formation and permit to estimate which kind of interaction are responsible of forming these complexes.     

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.     

Topoisomerase II poisoning by ICRF-193

Antineoplastic bis(dioxopiperazine)s, such as meso-2,3-bis(2,6-dioxopiperazin-4-yl)butane (ICRF-193), are widely believed to be only catalytic inhibitors of topoisomerase II. However, topoisomerase inhibitors have little or no antineoplastic activity unless they are topoisomerase poisons, a special subclass of topoisomerase-targeting drugs that stabilize topoisomerase-DNA strand passing intermediates and thus cause the topoisomerase to become a cytotoxic DNA-damaging agent. Here we report that ICRF-193 is a very significant topoisomerase II poison. Detection of topoisomerase II poisoning by ICRF-193 required the use of a chaotropic protein denaturant in the topoisomerase poisoning assays. ICRF-193 caused dose-dependent cross-linking of human topoisomerase IIbeta to DNA and stimulated topoisomerase IIbeta-mediated DNA cleavage at specific sites on (32)P-end-labeled DNA. Human topoisomerase IIalpha-mediated DNA cleavage was stimulated to a lesser extent by ICRF-193. In vivo experiments with MCF-7 cells also showed the requirement of a chaotropic protein denaturant in the assays and selectivity for the beta-isozyme of human topoisomerase II. Studies with two topoisomerase IIbeta-negative cell model systems confirmed significant topoisomerase II poisoning by ICRF-193 in the wild type cells and were consistent with beta-isozyme selectivity. Common use of only the detergent, SDS, in assays may have led to failure to detect topoisomerase II poisoning by ICRF-193 in earlier studies.     

The EcoRI restriction endonuclease with bacteriophage lambda DNA. Equilibrium binding studies.

The EcoRI restriction endonuclease was found by the filter binding technique to form stable complexes, in the absence of Mg2+, with the DNA from derivatives of bacteriophage lambda that either contain or lack EcoRI recognition sites. The amount of complex formed at different enzyme concentrations followed a hyperbolic equilibrium-binding curve with DNA molecules containing EcoRI recognition sites, but a sigmoidal equilibrium-binding curve was obtained with a DNA molecule lacking EcoRI recognition sites. The EcoRI enzyme displayed the same affinity for individual recognition sites on lambda DNA, even under conditions where it cleaves these sites at different rates. The binding of the enzyme to a DNA molecule lacking EcoRI sites was decreased by Mg2+. These observations indicate that (a) the EcoRI restriction enzyme binds preferentially to its recognition site on DNA, and that different reaction rates at different recognition sites are due to the rate of breakdown of this complex; (b) the enzyme also binds to other DNA sequences, but that two molecules of enzyme, in a different protein conformation, are involved in the formation of the complex at non-specific consequences; (c) the different affinities of the enzyme for the recognition site and for other sequences on DNA, coupled with the different protein conformations, account for the specificity of this enzyme for the cleavage of DNA at this recognition site; (d) the decrease in the affinity of the enzyme for DNA, caused by Mg2+, liberates binding energy from the DNA-protein complex that can be used in the catalytic reaction.     

Dna binding-independent anti-proliferative action of benzazolo[3,2-α]quinolinium DNA intercalators

The proposed mechanism of action of the antineoplastic drug 3-nitrobenzothiazolo[3,2-& agr;]quinolinium chloride (NBQ-2) involves its interaction with DNA by intercalation and inhibition of topoisomerase II activity by arresting the enzyme in a covalent cleavage complex. In an attempt to identify some structural determinants for activity and develop a molecular structure/cytotoxicity correlation, four new structural analogs of the antitumor NBQ-2 were prepared and their cytotoxic activity and DNA binding properties were investigated. The cytotoxic activity was evaluated against six different human tumor cell lines: U937, K-562, HL-60, HT-29, HeLa, and A431. The results showed that these new drugs elicit pronounced cytotoxic effects against U937, K-562, HL-60 and A431 while HeLa and HT-29 were less sensitive to the new drugs. This apparent selectivity was different to that of m-AMSA, a drug currently used for cancer treatment. Since the interaction of NBQ-2 to DNA by intercalation has been proposed as the initial step leading to its antineoplastic activity, DNA binding and changes in DNA contour length induced by the new NBQ-2 structural analogs were also investigated using calf thymus and human DNA. The drug, 7-(1-propenyl)-3-nitrobenzimidazolo[3,2-& agr;]quinolinium chloride (NBQ-59) was the most cytotoxic agent of the analog series (IC₅₀ = 16 & mgr;M for HL-60 cells), however, it demonstrated the weakest binding to DNA (K_int = 0.9 × 10₅ M^-1 for calf thymus DNA). NBQ-59 was also found to be a poor intercalator into the DNA double helix. Therefore, our results suggest that DNA binding is not the primary mechanism of drug action for this family of compounds. In addition structural determinants important for cytotoxicity of the benzazolo quinolinium chlorides were suggested by our results. In particular, the nitro group in the 3 position does not seem to be necessary for bioactivity, while substitutions in the benzazolo moiety have striking effects on the biological activity of the drugs.     

Ligand‐based pharmacophore modelling and screening of DNA minor groove binders targeting Staphylococcus aureus

The recognition of DNA by small molecules is of special importance in the design of new drugs. Many natural and synthetic compounds have the ability to interact with the minor groove of DNA. In the present study, identification of minor groove binding compounds was attained by the combined approach of pharmacophore modelling, virtual screening and molecular dynamics approach. Experimentally reported 32 minor groove binding compounds were used to develop the pharmacophore model. Based on the fitness score, best three pharmacophore hypotheses were selected and used as template for screening the compounds from drug bank database. This pharmacophore‐based screening provides many compounds with the same pharmacological properties. All these compounds were subjected to four phases of docking protocols with combined Glide‐quantum‐polarized ligand docking approach. Molecular dynamics results indicated that selected compounds are more active and showed good interaction in the binding site of DNA. Based on the scoring parameters and energy values, the best compounds were selected, and antibacterial activity of these compounds was identified using in vitro antimicrobial techniques. Copyright © 2014 John Wiley & Sons, Ltd.