On the molecular mechanism of drug intercalation into DNA: a simulation study of the intercalation pathway, free energy, and DNA structural changes

Intercalation into DNA (insertion between a pair of base pairs) is a critical step in the function of many anticancer drugs. Despite its importance, a detailed mechanistic understanding of this process at the molecular level is lacking. We have constructed, using extensive atomistic computer simulations and umbrella sampling techniques, a free energy landscape for the intercalation of the anticancer drug daunomycin into a twelve base pair B-DNA. A similar free energy landscape has been constructed for a probable intermediate DNA minor groove-bound state. These allow a molecular level understanding of aspects of the thermodynamics, DNA structural changes, and kinetic pathways of the intercalation process. Key DNA structural changes involve opening the future intercalation site base pairs toward the minor groove (positive roll), followed by an increase in the rise, accompanied by hydrogen bonding changes of the minor groove waters. The calculated intercalation free energy change is -12.3 kcal/mol, in reasonable agreement with the experimental estimate -9.4 kcal/mol. The results point to a mechanism in which the drug first binds to the minor groove and then intercalates into the DNA in an activated process, which is found to be in general agreement with experimental kinetic results.     

A comparative DFT study of some N-based aromatic ligand metal complexes as anticancer agents and analysis of their mode of interaction with DNA base pair

Metal complexed anticancer agents interact with DNA nucleobase pairs (AT and GC) through different types of binding mode such as intercalation, groove binding, covalent binding, etc. Minor and major groove binding mechanism of DNA base pair is the key factor for all kinds of anticancer agent; as metal complexes have a great affinity to bind with DNA nucleobase either through minor or major groove. Ligands in metal complexes also play a vital role during the interaction with DNA base pairs; these ligands directly interact with DNA through different interacting modes. Generally, anticancer agents with less sterically hindered N-based aromatic and planar ligands are the key component for DNA binding; as the structure of such ligands are quite compatible for following intercalation and groove binding mechanism. Since, the experimental investigation for drug-DNA nucleobase complexes are extremely complicated, therefore; quantum mechanical calculations might be very helpful for computing the actual interactions in drug-DNA complexes. Quantum mechanical approaches such as density functional theory (DFT) might be a very important and useful tool to investigate the actual mode of interaction of metal complexed antitumor agents with DNA nucleobase. Herein, we have taken some metal complexes with N-based aromatic ligands as antitumor agents to investigate the proper mode of interaction between drug-DNA complexes.     

Hydroxyl radical reactions with adenine: reactant complexes, transition states, and product complexes

In order to address problems such as aging, cell death, and cancer, it is important to understand the mechanisms behind reactions causing DNA damage. One specific reaction implicated in DNA oxidative damage is hydroxyl free-radical attack on adenine (A) and other nucleic acid bases. The adenine reaction has been studied experimentally, but there are few theoretical results. In the present study, adenine dehydrogenation at various sites, and the potential-energy surfaces for these reactions, are investigated theoretically. Four reactant complexes [A···OH]* have been found, with binding energies relative to A+OH* of 32.8, 11.4, 10.7, and 10.1 kcal mol(-1). These four reactant complexes lead to six transition states, which in turn lie +4.3, -5.4, (-3.7 and +0.8), and (-2.3 and +0.8) kcal mol(-1) below A+OH*, respectively. Thus the lowest lying [A···OH]* complex faces the highest local barrier to formation of the product (A-H)*+H(2)O. Between the transition states and the products lie six product complexes. Adopting the same order as the reactant complexes, the product complexes [(A-H)···H(2)O]* lie at -10.9, -22.4, (-24.2 and -18.7), and (-20.5 and -17.5) kcal mol(-1), respectively, again relative to separated A+OH*. All six A+OH* → (A-H)*+H(2)O pathways are exothermic, by -0.3, -14.7, (-17.4 and -7.8), and (-13.7 and -7.8) kcal mol(-1), respectively. The transition state for dehydrogenation at N(6) lies at the lowest energy (-5.4 kcal mol(-1) relative to A+OH*), and thus reaction is likely to occur at this site. This theoretical prediction dovetails with the observed high reactivity of OH radicals with the NH(2) group of aromatic amines. However, the high barrier (37.1 kcal mol(-1)) for reaction at the C(8) site makes C(8) dehydrogenation unlikely. This last result is consistent with experimental observation of the imidazole ring opening upon OH radical addition to C(8). In addition, TD-DFT computed electronic transitions of the N(6) product around 420 nm confirm that this is the most likely site for hydrogen abstraction by hydroxyl radical.     

Alternative heterocycles for DNA recognition: the benzimidazole/imidazole pair

Boc-protected benzimidazole-pyrrole, benzimidazole-imidazole, and benzimidazole-methoxypyrrole amino acids were synthesized and incorporated into DNA binding polyamides, comprised of N-methyl pyrrole and N-methyl imidazole amino acids, by means of solid-phase synthesis on an oxime resin. These hairpin polyamides were designed to determine the DNA recognition profile of a side-by-side benzimidazole/imidazole pair for the designated six base pair recognition sequence. Equilibrium association constants of the polyamide-DNA complexes were determined at two of the six base pair positions of the recognition sequence by quantitative DNase I footprinting titrations on DNA fragments each containing matched and single base pair mismatched binding sites. The results indicate that the benzimidazole-heterocycle building blocks can replace pyrrole-pyrrole, pyrrole-imidazole, and pyrrole-hydroxypyrrole constructs while retaining relative site specifities and subnanomolar match site affinities. The benzimidazole-containing hairpin polyamides represent a novel class of DNA binding ligands featuring tunable target recognition sequences combined with the favorable properties of the benzimidazole type DNA minor groove binders.     

Interaction of adenine adducts with thymine: a computational study

The existence of DNA adducts bring the danger of carcinogenesis because of mispairing with normal DNA bases. 1,N6-ethenoadenine adducts (epsilonA) and 1,N6-ethanoadenine adducts (EA) have been considered as DNA adducts to study the interaction with thymine, as DNA base. Several different stable conformers for each type of adenine adduct with thymine, [epsilonA(1)-T(I), epsilonA(2)-T(I), epsilonA(3)-T(I) and EA(1)-T(I), EA(2)-T(I), EA(3)-T(I)] and [epsilonA(1)-T(II), epsilonA(2)-T(II), epsilonA(3)-T(II) and EA(1)-T(II), EA(2)-T(II), EA(3)-T(II)], have been considered with regard to their interactions. The differences in their geometrical structures, energetic properties, and hydrogen-bonding strengths have also been compared with Watson-Crick adenine-thymine base pair (A-T). Single-point energy calculations at MP2/6-311++G** levels on B3LYP/6-31+G* optimized geometries have also been carried out to better estimate the hydrogen-bonding strengths. The basis set superposition error corrected hydrogen-bonding strength sequence at MP2/6-311++G**//B3LYP/6-31+G* for the most stable complexes is found to be EA(2)-T(I) (15.30 kcal/mol) > EA(1)-T(II) (14.98 kcal/mol) > EA(3)-T(II) (14.68 kcal/mol) > epsilonA(2)-T(I) (14.54 kcal/mol) > epsilonA(3)-T(II) (14.22 kcal/mol) > epsilonA(3)-T(II) (13.64 kcal/mol) > A-T (13.62 kcal/mol). The calculated reaction enthalpy value for epsilonA(2)-T(I) is 10.05 kcal/mol, which is the highest among the etheno adduct-thymine complexes and about 1.55 kcal/mol more than those obtained for Watson-Crick A-T base pair and the reaction enthalpy value for EA(1)- T(II) is 10.22 kcal/mol, which is highest among the ethano addcut-thymine complexes and about 1.72 kcal/mol more than those obtained for Watson-Crick A-T base pair. The aim of this research is to provide fundamental understanding of adenine adduct and thymine interaction at the molecular level and to aid in future experimental studies toward finding the possible cause of DNA damage.     

Sequence specific fluorescence detection of double strand DNA

Methods for the fluorescent detection of specific sequences of double strand DNA in homogeneous solution may be useful in the field of human genetics. A series of hairpin polyamides with tetramethyl rhodamine (TMR) attached to an internal pyrrole ring were synthesized, and the fluorescence properties of the polyamide-fluorophore conjugates in the presence and absence of duplex DNA were examined. We observe weak TMR fluorescence in the absence of DNA. Addition of >/=1:1 match DNA affords a significant fluorescence increase over equimolar mismatch DNA for each polyamide-TMR conjugate. Polyamide-fluorophore conjugates offer a new class of sensors for the detection of specific DNA sequences without the need for denaturation. The polyamide-dye fluorescence-based method can be used to screen in parallel the interactions between aromatic ring pairs and the minor groove of DNA even when the binding site contains a non-Watson-Crick DNA base pair. A ranking of the specificity of three polyamide ring pairs-Py/Py, Im/Py, and Im/Im-was established for all 16 possible base pairs of A, T, G, and C in the minor groove. We find that Im/Im is an energetically favorable ring pair for minor groove recognition of the T.G base pair.     

Azole-bridged diplatinum anticancer compounds. Modulating DNA flexibility to escape repair mechanism and avoid cross resistance

Dinuclear azole-bridged Pt compounds bind to DNA helices, forming intrastrand crosslinks between adjacent guanines in a similar way to cisplatin. Their cytotoxic profile is, however, different from that of first and second generation Pt drugs in that they lack cross resistance in cisplatin-resistant cell lines. In contrast to cisplatin, which induces a large kink in DNA duplex, structural NMR studies and molecular dynamics simulations have shown that azole-bridged diplatinum compounds induce only small structural changes in double-stranded DNA. These structural differences have been invoked to explain the different cytotoxic profile of these compounds. Here, we show that in addition to the small structural changes in DNA, dinuclear Pt compounds also affect DNA minor groove flexibility in a different way than cisplatin. Free-energy calculations on azole-bridged diplatinum DNA adducts reveal that opening of the minor groove requires a higher free-energy cost (DeltaG ~ 7-15 kcal/mol) than in the corresponding cisplatin-DNA adduct (DeltaG ~ 0 kcal/mol). This could prevent minor groove binding proteins from binding to diplatinum-DNA adducts thus leading to a different cellular response than cisplatin and possibly decreasing the activity of excision repair enzymes. Although the development of drug resistance is a highly complex mechanism, our findings provide an additional rationale for the improved cytotoxic activity of these compounds in cell lines resistant to cisplatin.     

Computation of relative binding free energy for an inhibitor and its analogs binding with Erk kinase using thermodynamic integration MD simulation

In the present study, we carried out thermodynamic integration molecular dynamics simulation for a pair of analogous inhibitors binding with Erk kinase to investigate how computation performs in reproducing the relative binding free energy. The computation with BCC-AM1 charges for ligands gave ?1.1?kcal/mol, deviated from experimental value of ?2.3?kcal/mol by 1.2?kcal/mol, in good agreement with experimental result. The error of computed value was estimated to be 0.5?kcal/mol. To obtain convergence, switching vdw interaction on and off required approximately 10 times more CPU time than switching charges. Residue-based contributions and hydrogen bonding were analyzed and discussed. Furthermore, subsequent simulation using RESP charge for ligand gave ΔΔG of ?1.6?kcal/mol. The computed results are better than the result of ?5.6?kcal/mol estimated using PBSA method in a previous study. Based on these results, we further carried out computations to predict ΔΔG for five new analogs, focusing on placing polar and nonpolar functional groups at the meta site of benzene ring shown in the Fig.?1, to see if these ligands have better binding affinity than the above ligands. The computations resulted that a ligand with polar –OH group has better binding affinity than the previous examined ligand by ~2.0?kcal/mol and two other ligands have better affinity by ~1.0?kcal/mol. The predicted better inhibitors of this kind should be of interest to experimentalist for future experimental enzyme and/or cell assays.     

Sequence-specific trapping of topoisomerase I by DNA binding polyamide-camptothecin conjugates

Hairpin pyrrole-imidazole polyamides are synthetic ligands that bind in the minor groove of DNA with affinities and specificities comparable to those of DNA binding proteins. Three polyamide-camptothecin conjugates 1-3 with linkers varying in length between 7, 13, and 18 atoms were synthesized to trap the enzyme Topoisomerase I and induce cleavage at predetermined DNA sites. One of these, polyamide-camptothecin conjugate 3 at nanomolar concentration (50 nM) in the presence of Topo I (37 degrees C), induces DNA cleavage between three and four base pairs from the polyamide binding site in high yield (77%).     

Association and insertion complexes of nickel with water and methanol studied using Kohn-Sham theory

Summary The complexes formed by simple association of the nickel atom with water and methanol, and those formed by insertion of nickel into O-H and O-C bonds, have been studied by the Linear Combination of Gaussian-type orbitalsdensity functional (LCGTO-DF) method. The binding energies calculated for theassociation complexes are 7 kcal/mol for both Ni(H₂O) and Ni(CH₃OH). These association complexes have equilibrium geometries characterized by a tetrahedral arrangement of bonds and lone pair around the oxygen atom. Theinsertion complexes HNiOH and HNiOCH₃ are more stable than the association complexes by 14 and 18 kcal/mol, respectively. The H₃CNiOH insertion complex is still more stable, by an additional 14 kcal/mol. In all five complexes, the lowest singlet and triplet states are nearly degenerate. Comparison of calculated frequencies with infrared spectra of matrix isolated species [MA Park (1988) PhD thesis, Rice University] indicates that only triplet states are seen in the experiment.     

寡聚酰胺为探针的电化学DNA生物传感器

DNA分子中的碱基对可以长程传递电荷, DNA分子中的碱基π堆积结构为电荷的长程传递提供了良好的通道. 电荷在DNA分子中的传递受碱基序列的影响, 利用这种性质可以构建DNA碱基错配检测的电化学传感器. 寡聚酰胺能和DNA以小沟绑定方式高亲和力地结合, 并且具有序列识别功能, 本文以带有硝基官能团的寡聚酰胺分子为电化学探针, 设计了电化学DNA生物传感器. 结果显示, 寡聚酰胺与DNA修饰电极作用后, 电化学响应显著增强, 并且可以作为检测DNA碱基错配的电化学探针分子.     

Toward the exact solution of the electronic Schrödinger equation for noncovalent molecular interactions: worldwide distributed quantum monte carlo calculations

Quantum Monte Carlo (QMC) calculations on the stacked (st) and Watson/Crick (wc) bound adenine/thymine (A/T) and cytosine/guanine (C/G) DNA base pair complexes were made possible with the first large scale distributed computing project in ab initio quantum chemistry, Quantum Monte Carlo at Home (QMC@HOME). The results for the interaction energies (wc-A/T = 15.7 kcal/mol, wc-C/G = 30.2 kcal/mol, st-A/T = 13.1 kcal/mol, st-C/G = 19.6 kcal/mol) are in very good agreement with the best known coupled-cluster based estimates. The accuracy of these values is further supported by calculations on the S22 benchmark set of noncovalently bound systems, for which we obtain a small mean absolute deviation of 0.68 kcal/mol. Our results support previous claims that the stacking energies are of comparable magnitude to the interactions of the commonly discussed hydrogen-bonded motif. Furthermore, we show that QMC can serve as an advantageous alternative to conventional wave function methods for large noncovalently bound systems. We also investigated in detail all technical parameters of the QMC simulations and recommend a careful optimization procedure of the Jastrow correlation factors in order to obtain numerically stable and reliable results.     

Tetraplatinated artificial oligopeptides afford high affinity intercalation into dsDNA

This paper reports the binding of an artificial tetrapeptide to which are tethered four Pt(II) complexes (i.e., [Pt(tpy)(py)]48+) with a 12 base pair duplex DNA oligonucleotide. Isothermal Titration Calorimetry reveals that two tetrametallic peptides stoichiometrically bind to each DNA duplex with a binding constant, KB, of 1.7 x 106 M-1, with a change in free energy of -8.5 kcal/mol. This KB represents an affinity 2 orders of magnitude greater than that of the monometallic analogue [Pt(tpy)(pic)]2+ for the same dsDNA sequence. The metalated peptides bind by intercalation into the DNA, partially unwinding the helix while stabilizing the structure, causing an increase in the dsDNA melting temperature of 25 degrees C.     

DNA-binding of semirigid binuclear ruthenium complex delta,delta-[mu-(11,11'-bidppz)(phen)(4)ru(2)](4+): extremely slow intercalation kinetics

We here report a remarkably slow rearrangement of binding modes for a binuclear ruthenium(II) complex upon interaction with DNA. It has been previously shown that Delta,Delta-[mu-(11,11'-bidppz)(phen)4Ru2]4+ binds to DNA in one of the grooves. However, we find that this is only an initial, metastable, binding mode, which is extremely slowly reorganized into an intercalative binding geometry. The slow rearrangement and dissociation, revealed by flow linear dichroism and fluorescence spectroscopy, are concluded to be a result from the complex being threaded through the DNA, with one of the bridging aromatic dppz ligands intercalated between the base pairs of the DNA, placing one metal center in the minor groove and one in the major groove. A negative LD, a high luminescence quantum yield, and long luminescence lifetimes, similar to the intercalating complex Delta-[Ru(phen)2dppz]2+, indicate intercalation of the bidppz moiety. The unique slow dissociation of the complex in its final DNA-binding mode suggests that this class of threading, partially intercalated binuclear complexes may be interesting in the context of cancer therapy. Also, their unique optical and photophysical properties could make such complexes, either alone or scaffolded by DNA structures, of interest for the development of nanometer-sized molecular optoelectronic devices.     

How and how much can Hoechst 33258 cause unwinding in a DNA duplex?

The effect of Hoechst 33258 binding on the geometry of a DNA duplex (plasmid pBR322) has been examined using topoisomerase II relaxation followed by gel electrophoresis. Of this drug-DNA system, fluorescence, optical absorption, and calorimetric measurements were also made at various drug and DNA concentrations and in the same buffer as that for the topoisomerase reaction. It has been confirmed that there are two modes of drug-DNA interaction. When the drug concentration is much lower than the DNA base pair concentration, the Hoechst 33258 molecule binds in the minor groove of the DNA duplex and occupies a site formed of five continuous base pair sequences that contain no G.C pair. Here, the equilibrium constant K1 is 1.8 x 10(7) M-1 (at 37 degrees C), and the enthalpy of binding delta H1 is -865 cal/mol. When the drug concentration is much higher, on the other hand, it shows another binding mode which is much weaker, so that K2 = 2.25 x 10(4) M-1 and delta H2 is -464 cal/mol, which gives fluorescence quenching, which has no base pair preference, and which causes an unwinding of the duplex by 1 degree.     

Use of MM-PBSA in reproducing the binding free energies to HIV-1 RT of TIBO derivatives and predicting the binding mode to HIV-1 RT of efavirenz by docking and MM-PBSA.

In this work, a new ansatz is presented that combines molecular dynamics simulations with MM-PBSA (Molecular Mechanics Poisson-Boltzmann/surface area) to rank the binding affinities of 12 TIBO-like HIV-1 RT inhibitors. Encouraging results have been obtained not only for the relative binding free energies, but also for the absolute ones, which have a root-mean-square deviation of 1.0 kcal/mol (the maximum error is 1.89 kcal/mol). Since the root-mean-square error is rather small, this approach can be reliably applied in ranking the ligands from the databases for this important target. Encouraged by the results, we decided to apply MM-PBSA combined with molecular docking to determine the binding mode of efavirenz SUSTIVA(TM) another promising HIV-1 RT inhibitor for which no ligand-protein crystal structure had been published at the time of this work. To proceed, we define the following ansatz: Five hundred picosecond molecular dynamics simulations were first performed for the five binding modes suggested by DOCK 4.0, and then MM-PBSA was carried out for the collected snapshots. MM-PBSA successfully identified the correct binding mode, which has a binding free energy about 7 kcal/mol more favorable than the second best mode. Moreover, the calculated binding free energy (-13.2 kcal/mol) is in reasonable agreement with experiment (-11.6 kcal/mol). In addition, this procedure was also quite successful in modeling the complex and the structure of the last snapshot was quite close to that of the measured 2,3 A resolution crystal (structure the root-mean-square deviation of the 54 C(alpha) around the binding site and the inhibitor is 1.1 A). We want to point out that this result was achieved without prior knowledge of the structure of the efavirenz/RT complex. Therefore, molecular docking combined with MD simulations followed by MM-PBSA analysis is an attractive approach for modeling protein complexes a priori.     

Free energy simulations of uncatalyzed DNA replication fidelity: structure and stability of T.G and dTTP.G terminal DNA mismatches flanked by a single dangling nucleotide

A reference system for DNA replication fidelity was studied by free energy perturbation (FEP) and linear interaction energy (LIE) methods. The studied system included a hydrated duplex DNA with the 5'-CG dangling end of the templating strand, and dCTP4-.Mg2+ or dTTP4-.Mg2+ inserted opposite the dangling G to form a correct (i.e., Watson-Crick) or incorrect (i.e., wobble) base pair, respectively. The average distance between the 3'-terminal oxygen of the primer strand and the alpha-phosphorus of dNTP was found to be 0.2 A shorter for the correct base pair than for the incorrect base pair. Binding of the incorrect dNTP was found to be disfavored by 0.4 kcal/mol relative to the correct dNTP. We estimated that improved binding and more near-attack configurations sampled by the correct base pair should translate in aqueous solution and in the absence of DNA polymerase into a six times faster rate for the incorporation of the correct dNTP into DNA. The accuracy of the calculated binding free energy difference was verified by examining the relative free energy for melting duplex DNA containing GC and GT terminal base pairs flanked by a 5' dangling C. The calculated LIE and FEP free energies of 1.7 and 1.1 kcal/mol, respectively, compared favorably with the experimental estimate of 1.4 kcal/mol obtained using the nearest neighbor parameters. To decompose the calculated free energies into additive electrostatic and van der Waals contributions and to provide a set of rigorous theoretical data for the parametrization of the LIE method, we suggested a variant of the FEP approach, for which we coined a binding-relevant free energy (BRFE) acronym. BRFE approach is characterized by its unique perturbation pathway and by its exclusion of the intramolecular energy of a rigid part of the ligand from the total potential energy.     

In silico footprinting of ligands binding to the minor groove of DNA

The sequence selectivity of small molecules binding to the minor groove of DNA can be predicted by "in silico footprinting". Any potential ligand can be docked in the minor groove and then moved along it using simple simulation techniques. By applying a simple scoring function to the trajectory after energy minimization, the preferred binding site can be identified. We show application to all known noncovalent binding modes, namely 1:1 ligand:DNA binding (including hairpin ligands) and 2:1 side-by-side binding, with various DNA base pair sequences and show excellent agreement with experimental results from X-ray crystallography, NMR, and gel-based footprinting.