1
|
Jawad B, Poudel L, Podgornik R, Steinmetz NF, Ching WY. Molecular mechanism and binding free energy of doxorubicin intercalation in DNA. Phys Chem Chem Phys 2019; 21:3877-3893. [PMID: 30702122 DOI: 10.1039/c8cp06776g] [Citation(s) in RCA: 54] [Impact Index Per Article: 10.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023]
Abstract
The intercalation process of binding doxorubicin (DOX) in DNA is studied by extensive MD simulations. Many molecular factors that control the binding affinity of DOX to DNA to form a stable complex are inspected and quantified by employing continuum solvation models for estimating the binding free energy. The modified MM-PB(GB)SA methodology provides a complete energetic profile of ΔGele, ΔGvDW, ΔGpolar, ΔGnon-polar, TΔStotal, ΔGdeform, ΔGcon, and ΔGion. To identify the sequence specificity of DOX, two different DNA sequences, d(CGATCG) or DNA1 and d(CGTACG) or DNA2, with one molecule (1 : 1 complex) or two molecule (2 : 1 complex) configurations of DOX were selected in this study. Our results show that the DNA deformation energy (ΔGdeform), the energy cost from translational and rotational entropic contributions (TΔStran+rot), the total electrostatic interactions (ΔGpolar-PB/GB + ΔGele) of incorporation, the intramolecular electrostatic interactions (ΔGele) and electrostatic polar solvation interactions (ΔGpolar-PB/GB) are all unfavorable to the binding of DOX to DNA. However, they are overcome by at least five favorable interactions: the van der Waals interactions (ΔGvDW), the non-polar solvation interaction (ΔGnon-polar), the vibrational entropic contribution (TΔSvib), and the standard concentration dependent free energies of DOX (ΔGcon) and the ionic solution (ΔGion). Specifically, the van der Waals interaction appears to be the major driving force to form a stable DOX-DNA complex. We also predict that DOX has stronger binding to DNA1 than DNA2. The DNA deformation penalty and entropy cost in the 2 : 1 complex are less than those in the 1 : 1 complex, thus they indicate that the 2 : 1 complex is more stable than the 1 : 1 complex. We have calculated the total binding free energy (BFE) (ΔGt-sim) using both MM-PBSA and MM-GBSA methods, which suggests a more stable DOX-DNA complex at lower ionic concentration. The calculated BFE from the modified MM-GBSA method for DOX-DNA1 and DOX-DNA2 in the 1 : 1 complex is -9.1 and -5.1 kcal mol-1 respectively. The same quantities from the modified MM-PBSA method are -12.74 and -8.35 kcal mol-1 respectively. The value of the total BFE ΔGt-sim in the 1 : 1 complex is in reasonable agreement with the experimental value of -7.7 ± 0.3 kcal mol-1.
Collapse
Affiliation(s)
- Bahaa Jawad
- Department of Physics and Astronomy, University of Missouri-Kansas City, Kansas City, Missouri 64110, USA.
| | | | | | | | | |
Collapse
|
2
|
Scaglioni L, Mondelli R, Artali R, Sirtori FR, Mazzini S. Nemorubicin and doxorubicin bind the G-quadruplex sequences of the human telomeres and of the c-MYC promoter element Pu22. Biochim Biophys Acta Gen Subj 2016; 1860:1129-38. [PMID: 26922833 DOI: 10.1016/j.bbagen.2016.02.011] [Citation(s) in RCA: 34] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/09/2015] [Revised: 02/12/2016] [Accepted: 02/21/2016] [Indexed: 12/22/2022]
Abstract
BACKGROUND Intra-molecular G-quadruplex structures are present in the guanine rich regions of human telomeres and were found to be prevalent in gene promoters. More recently, the targeting of c-MYC transcriptional control has been suggested, because the over expression of the c-MYC oncogene is one of the most common aberration found in a wide range of human tumors. METHODS The interaction of nemorubicin and doxorubicin with DNA G-quadruplex structures has been studied by NMR, ESI-MS and molecular modelling, in order to obtain further information about the complex and the multiple mechanisms of action of these drugs. RESULTS AND CONCLUSIONS Nemorubicin intercalates between A3 and G4 of d(TTAGGGT)4 and form cap-complex at the G6pT7 site. The presence of the adenine in this sequence is important for the stabilization of the complex, as was shown by the interaction with d(TTGGGTT)4 and d(TTTGGGT)4, which form only a 1:1 complex. The interaction of doxorubicin with d(TTAGGGT)4 is similar, but the complex appears less stable. Nemorubicin also binds with high efficiency the c-MYC G-quadruplex sequence Pu22, to form a very well defined complex. Two nemorubicin molecules bind to the 3'-end and to the 5'-end, forming an additional plane of stacking over each external G-tetrad. The wild type c-MYCPu22 sequence forms with nemorubicin the same complex. GENERAL SIGNIFICANCE Nemorubicin and doxorubicin, not only intercalate into the duplex DNA, but also result in significant ligands for G-quadruplex DNA segments, stabilizing their structure; this may in part explain the multiple mechanisms of action of their antitumor activity.
Collapse
Affiliation(s)
- Leonardo Scaglioni
- Department of Food, Environmental and Nutritional Sciences, Division of Chemistry and Molecular Biology, University of Milan, Via Celoria 2, 20133 Milano, Italy
| | - Rosanna Mondelli
- Department of Food, Environmental and Nutritional Sciences, Division of Chemistry and Molecular Biology, University of Milan, Via Celoria 2, 20133 Milano, Italy
| | | | - Federico Riccardi Sirtori
- Nerviano, Medical Sciences, Oncology-Chemical Core, Technologies Department, viale Pasteur, 10, 20014 Nerviano, Milano, Italy
| | - Stefania Mazzini
- Department of Food, Environmental and Nutritional Sciences, Division of Chemistry and Molecular Biology, University of Milan, Via Celoria 2, 20133 Milano, Italy.
| |
Collapse
|
3
|
Mazzini S, Scaglioni L, Mondelli R, Caruso M, Sirtori FR. The interaction of nemorubicin metabolite PNU-159682 with DNA fragments d(CGTACG)(2), d(CGATCG)(2) and d(CGCGCG)(2) shows a strong but reversible binding to G:C base pairs. Bioorg Med Chem 2012; 20:6979-88. [PMID: 23154079 DOI: 10.1016/j.bmc.2012.10.033] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/20/2012] [Revised: 10/11/2012] [Accepted: 10/14/2012] [Indexed: 11/20/2022]
Abstract
The antitumor anthracycline nemorubicin is converted by human liver microsomes to a major metabolite, PNU-159682 (PNU), which was found to be much more potent than its parent drug toward cultured tumor cells and in vivo tumor models. The mechanism of action of nemorubicin appears different from other anthracyclines and until now is the object of studies. In fact PNU is deemed to play a dominant, but still unclear, role in the in vivo antitumor activity of nemorubicin. The interaction of PNU with the oligonucleotides d(CGTACG)(2), d(CGATCG)(2) and d(CGCGCG)(2) was studied with a combined use of (1)H and (31)P NMR spectroscopy and by ESI-mass experiments. The NMR studies allowed to establish that the intercalation between the base pairs of the duplex leads to very stable complexes and at the same time to exclude the formation of covalent bonds. Melting experiments monitored by NMR, allowed to observe with high accuracy the behaviour of the imine protons with temperature, and the results showed that the re-annealing occurs after melting. The formation of reversible complexes was confirmed by HPLC-tandem mass spectra, also combined with endonuclease P1digestion. The MS/MS spectra showed the loss of neutral PNU before breaking the double helix, a behaviour typical of intercalators. After digestion with the enzyme, the spectra did not show any compound with PNU bound to the bases. The evidence of a reversible process appears from both proton and phosphorus NOESY spectra of PNU bound to d(CGTACG)(2) and to d(CGATCG)(2). The dissociation rate constants (k(off)) of the slow step of the intercalation process, measured by (31)P NMR NOE-exchange experiments, showed that the kinetics of the process is slower for PNU than for doxorubicin and nemorubicin, leading to a 10- to 20-fold increase of the residence time of PNU into the intercalation sites, with respect to doxorubicin. A relevant number of NOE interactions allowed to derive a model of the complexes in solution from restrained MD calculations. The conformation of PNU bound to the oligonucleotides was also derived from the coupling constant values.
Collapse
Affiliation(s)
- Stefania Mazzini
- DeFENS-Department of Food, Environmental and Nutritional Sciences, via Celoria 2, 20133 Milano, Italy.
| | | | | | | | | |
Collapse
|
4
|
Interaction between double helix DNA fragments and the new antitumor agent sabarubicin, Men10755. Bioorg Med Chem 2010; 18:1497-506. [DOI: 10.1016/j.bmc.2010.01.012] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/31/2009] [Revised: 01/05/2010] [Accepted: 01/06/2010] [Indexed: 11/19/2022]
|
5
|
Agrawal P, Govil G, Barthwal R. Studies on drug-DNA complexes, adriamycin-d-(TGATCA)(2) and 4'-epiadriamycin-d-(CGATCG)(2), by phosphorus-31 nuclear magnetic resonance spectroscopy. MAGNETIC RESONANCE IN CHEMISTRY : MRC 2009; 47:390-397. [PMID: 19170249 DOI: 10.1002/mrc.2398] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/27/2023]
Abstract
The complexes of adriamycin-d-(TGATCA)(2) and 4'-epiadriamycin-d-(CGATCG)(2) are studied by one- and two-dimensional (31)P nuclear magnetic resonance spectroscopy (NMR) at 500 MHz in the temperature range 275-328 K and as a function of drug to DNA ratio (0.0-2.0). The binding of drug to DNA is clearly evident in (31)P-(31)P exchange NOESY spectra that shows two sets of resonances in slow chemical exchange. The phosphate resonances at the intercalating steps, T1pG2/C1pG2 and C5pA6/C5pG6, shift downfield up to 1.7 ppm and that at the adjacent step shift downfield up to 0.7 ppm, whereas the central phosphate A3pT4 is relatively unaffected. The variations of chemical shift with drug to DNA ratio and temperature as well as linewidths are different in each of the two complexes. These observations reflect change in population of B(I)/B(II) conformation, stretching of backbone torsional angle zeta, and distortions in O-P-O bond angles that occur on binding of drug to DNA. To the best of our knowledge, there are no solution studies on 4'-epiadriamycin, a better tolerated drug, and binding of daunomycin or its analogue to d-(TGATCA)(2) hexamer sequence. The studies report the use of (31)P NMR as a tool to differentiate various complexes. The specific differences may well be the reasons that are responsible for different antitumor action of these drugs due to different binding ability and distortions in DNA.
Collapse
Affiliation(s)
- Prashansa Agrawal
- Department of Biotechnology, Indian Institute of Technology Roorkee, Roorkee 247667, India
| | | | | |
Collapse
|
6
|
Scaglioni L, Mazzini S, Mondelli R, Dallavalle S, Gattinoni S, Tinelli S, Beretta GL, Zunino F, Ragg E. Interaction between double helix DNA fragments and a new topopyrone acting as human topoisomerase I poison. Bioorg Med Chem 2008; 17:484-91. [PMID: 19097800 DOI: 10.1016/j.bmc.2008.12.005] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/03/2008] [Revised: 11/28/2008] [Accepted: 12/02/2008] [Indexed: 11/16/2022]
Abstract
A water soluble derivative (2) of topopyrones was selected for NMR studies directed to elucidate the mode of binding with specific oligonucleotides. Topopyrone 2 can intercalate into the CG base pairs, but the residence time into the double helix is very short and a fast chemical exchange averaging occurs at room temperature between the free and bound species. The equilibria involved become slow below room temperature, thus allowing to measure a mean lifetime of the complex of ca. 7 ms at 15 degrees C. Structural models of the complex with d(CGTACG)(2) were developed on the basis of DOSY, 2D NOESY and (31)P NMR experiments. Topopyrone 2 presents a strong tendency to self-associate. In the presence of oligonucleotide a certain number of ligand molecules are found to externally stack to the double-helix, in addition to a small fraction of the same ligand intercalated. The external binding to the ionic surface of the phosphoribose chains may thus represents the first step of the intercalation process.
Collapse
Affiliation(s)
- Leonardo Scaglioni
- Dipartimento di Scienze Molecolari Agroalimentari, Università di Milano, via Celoria 2, 20133 Milano, Italy
| | | | | | | | | | | | | | | | | |
Collapse
|
7
|
Hertweck C, Luzhetskyy A, Rebets Y, Bechthold A. Type II polyketide synthases: gaining a deeper insight into enzymatic teamwork. Nat Prod Rep 2007; 24:162-90. [PMID: 17268612 DOI: 10.1039/b507395m] [Citation(s) in RCA: 386] [Impact Index Per Article: 22.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
This review covers advances in understanding of the biosynthesis of polyketides produced by type II PKS systems at the genetic, biochemical and structural levels.
Collapse
Affiliation(s)
- Christian Hertweck
- Leibniz Institute for Natural Product Research and Infection Biology, HKI, Beutenbergstrasse 11a, 07745 Jena, Germany
| | | | | | | |
Collapse
|
8
|
Kolakowski B, Battaglini F, Lee YS, Klironomos G, Mikkelsen SR. Comparison of an intercalating dye and an intercalant-enzyme conjugate for DNA detection in a microtiter-based assay. Anal Chem 1996; 68:1197-200. [PMID: 8651491 DOI: 10.1021/ac950878m] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/01/2023]
Abstract
Two methods have been developed for the detection of DNA immobilized on the surface of microtiter wells. An intercalating dye, 3,6-diaminoacridine, is used in stain and rinse solutions, so that measured absorbance values (450 nm) reflect the sum of DNA-bound and free dye. With diaminoacridine, signal increases of 0.056 +/- 0.010 were achieved on immobilizing double-stranded calf thymus DNA. An intercalant-enzyme conjugate, consisting of an average of four daunomycin moieties covalently bound to each glucose oxidase, was shown to provide a 10-fold signal enhancement (optimum 0.25 microM, with rinsing and peroxidase-o-dianisidine detection) compared to diaminoacridine, due to catalytic amplification; signals of 0.50 +/- 0.05 were obtained. This conjugate possesses 56% of the activity of native glucose oxidase and was prepared using water-soluble carbodiimide and N-hydroxysuccinimide reagents. Single-stranded DNA was immobilized onto avidin-coated polystyrene plates and commercially available (Covalink) plates possessing secondary amine groups. Following hybridization with complementary DNA, detection was performed with the daunomycin-glucose oxidase conjugate. Both immobilization methods showed optimum DNA concentrations of 0.10 microgram/mL, and maximum signal intensities were obtained when > 0.5 microgram/mL complementary DNA was present in the hybridization solution. Some nonspecific binding of the intercalant-enzyme conjugate was suggested by results obtained with avidin-coated polystyrene plates, but not with Covalink plates.
Collapse
Affiliation(s)
- B Kolakowski
- Department of Chemistry and Biochemistry, Concordia University, Montreal, Canada
| | | | | | | | | |
Collapse
|
9
|
|
10
|
Rizzo V, Battistini C, Vigevani A, Sacchi N, Razzano G, Arcamone F, Garbesi A, Colonna FP, Capobianco M, Tondelli L. Association of anthracyclines and synthetic hexanucleotides. Structural factors influencing sequence specificity. J Mol Recognit 1989; 2:132-41. [PMID: 2636902 DOI: 10.1002/jmr.300020306] [Citation(s) in RCA: 17] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
The equilibrium and kinetic aspects of the interaction between four anthracyclines and two synthetic self-complementary hexanucleotides was investigated by fluorescence detection. Two of the studied anthracyclines are widely used antitumor drugs: doxorubicin (1, formerly adriamycin) and daunorubicin (2, formerly daunomycin). The other two, 9-deoxydoxorubicin (3) and 3'-deamino-3'-hydroxy-4'-epidoxorubicin (4), are doxorubicin analogues with modifications of the chemical groups that have been proposed as responsible for sequence specificity (Chen, K.-X., Gresh, N. and Pullman, B. (1985). J. Biomol. Struct. Dyn. 3, 445-466). One of the oligonucleotides, d(CGTACG), is identical to that used in the high resolution x-ray structure determination of the daunorubicin intercalative complex (Wang, A. H.-J., Ughetto, G., Quigley, G. J. & Rich, A. (1987). Biochemistry 26, 1152-1163). Binding to this hexanucleotide is compared with intercalation into the d(CGCGCG) duplex, revealing sequence preferences of the four anthracyclines. Taking into account the anthracycline aggregation and the dissociation of the hexanucleotide double standard form, results can be interpreted with a model that assumes complete fluorescence quenching at intercalative sites containing the CG base pair, and a large residual fluorescence after intercalation within the TpA fragment. All four anthracyclines show preferential intercalation at sites near the ends of both hexanucleotide duplexes, partly as a result of positive cooperativity in the formation of di-intercalated species at these sites. Within the limits of experimental error, complete site specificity for the CpG fragment is found in the intercalation of 1 and 2 into d(CGTACG) duplex, whereas analogues 3 and 4 give increasing evidence of intercalation at other sites including the fluorescence-preserving TpA fragment. Site specificity is less pronounced in the association with d(CGCGCG), when cooperativity is taken into account. Kinetic data corroborate the results of equilibrium studies and are interpreted with a mechanism that includes formation of an intermediate bound species followed by drug redistribution to preferential sites. Finally, from a comparison of pertinent site binding constants, approximate free energy contributions to sequence specific DNA interaction, due to C9-OH on the aglycone and -NH3+ on daunosamine, are estimated not to exceed 2 kcal/mol.
Collapse
Affiliation(s)
- V Rizzo
- Farmitalia Carlo Erba/Research & Development-Erbamont Group, Milan, Italy
| | | | | | | | | | | | | | | | | | | |
Collapse
|