Spectroscopic, viscositic and molecular modeling studies on the interaction of 3'-azido-daunorubicin thiosemicarbazone with DNA.

A new daunorubicin has been synthesized and structurally characterized. The interaction of native calf thymus DNA (ctDNA) with 3'-azido-daunorubicin thiosemicarbazone (ADNRT) was investigated under simulated physiological conditions by multi-spectroscopic techniques, viscometric measurements and molecular modeling study. It concluded that ADNRT could intercalate into the base pairs of ctDNA, and the fluorescence quenching by ctDNA was static quenching type. Thermodynamic parameters calculated suggested that the binding of ADNRT to ctDNA was mainly driven by hydrophobic interactions. The relative viscosity of ctDNA increased with the addition of ADNRT, which confirmed the intercalation mode. Furthermore, molecular modeling studies corroborate the above experimental results.

Introduction

As we all know, nucleic acid is the primary intracellular target of antitumor drugs, owing to small molecules interaction with DNA and then damaging the formation of DNA in cancer cells [1]. Thus, the investigation of small molecules binding to DNA is of current general interest and importance in rational design of more powerful and selective anticancer pharmaceuticals and their effectiveness depend on the mode and affinity of the binding [2].

Daunorubicin, one of the most widely used anthracyclines, is used for the clinical treatment of a broad range of human malignancies [3]. Like other intercalating agents, daunorubicin is considered to interfere with the function of topoisomerase II, a nuclear enzyme which can regulate the DNA topology [4]. In spite that a great deal of work has been done [5], the mechanism by which these molecules mediate their relatively selective antitumor effects is not yet clear. What’s more, severe side effects limited the clinical use of daunorubicin, such as cardiotoxicity, myelosuppression and drug resistance [6]. Many efforts have been made to modify the structure of anthracyclines to generate analogs for the sake of reducing the side effects and reversing multidrug resistance. In our previous works, a lead compound derived from daunorubicin, 3′-azido-3′-deamino daunorubicin (ADNR), was synthesized and its bioactivities and the binding to human serum albumin or DNA were also investigated [7-9]. This lead compound exhibited potent anticancer activity in both drug-sensitive (K562) and drug-resistant leukemia cells (K562/Dox), with a 25-fold lower drug resistance index than daunorubicin.

Based on the special bioactivities of ADNR and thiosemicarbazones [10-13], in this paper, we designed and synthesized a hybrid of ADNR with thiosemicabzone, 3′-azido-daunorubicin thiosemicarbazone (ADNRT, inset of Fig. 1), and systematically investigated its interaction with calf thymus DNA (ctDNA) in combination with multi-spectroscopic, viscometric techniques and molecular modeling under simulated physiological conditions (pH 7.4). The characteristics in spectroscopy measurements revealed that ADNRT could bind to ctDNA through intercalation binding modes. The results obtained from the viscosity experiments and molecular modeling validated those conclusions. This work provides valuable information on the binding behavior of daunorubicin with DNA and should be helpful for designing the alternative or even more active daunorubicin.

Fig. 1

Fluorescence spectra of ADNRT in the presence of ctDNA (inset: structure of ADNRT). C ctDNA = 0, 4.6, 9.2, 13.8, 18.4, 23.0 and 27.6 × 10−6 mol · L−1 for curves 1–7; C ADNRT = 3.3 × 10−6 mol · L−1

Materials and Methods

Materials

The preparation procedure for ADNRT: A drop of acetic acid was added to the mixture of ADNR (27.7 mg, 0.05 mmol) and thiosemicarbazide monohydrochloride (19.1 mg,0.15 mmol) in methanol (3 mL) and the mixture was stirred overnight. TLC showed the reaction was completed. The solvent was removed, then the residue was diluted with saturated aqueous Na2CO3 solution and extracted with CH2Cl2 (5 mL) for three times. The combined extractions were dried by Na2SO4 and concentrated. The crude product was purified on a column of silica gel to get the pure product ADNRT as a yellow solid (29.1 mg, 93 %). HRMS (M + H)+ (ESI) calcd for C28H31N6O9S 627.1867, found 627.1860. 1H NMR (400 MHz, DMSO-d ):13.99 (1H, s), 13.21 (1H, s), 10.08 (1H, s), 8.22 (1H, s), 7.87–7.80 (2H, m), 7.57 (1H, d, J = 8.0 Hz), 7.24 (1H, s), 5.31 (1H, d, J = 4.0 Hz), 5.17 (2H, m), 4.84 (1H, s), 4.10–4.05 (1H, m), 3.94 (1H, s), 3.60 (1H, d, J = 4.0 Hz), 3.49–3.46 (1H, m), 3.13 (1H, d, J = 16.0 Hz), 2.83 (1H, d, J = 20.0 Hz), 2.31 (1H, dd, J = 4.0 Hz, J = 12.0 Hz), 2.15 (1H, dd, J = 4.0 Hz, J = 12.0 Hz), 2.00 (3H, s), 1.69 (1H, dd, J = 4.0 Hz, J = 12.0 Hz), 1.14 (3H, d, J = 8.0 Hz); 13C NMR (100 MHz, DMSO-d ):186.3, 186.1, 179.1, 160.7, 156.2, 155.5, 154.2, 136.1, 135.7, 135.5, 134.5, 119.8, 119.5, 118.9, 110.5, 99.6, 72.2, 70.7, 69.0, 67.0, 56.5, 55.6, 38.1, 33.4, 28.1, 16.9, 11.9.

The stock solution (8.3 × 10−5 mol L−1) of ADNRT (synthesized) was prepared. An appropriate amount of solid ctDNA (Sigma-Aldrich Chemicals Company, USA) without further purification was dissolved into the doubly distilled water and stored at 4 °C. The ratio of the absorbance at 260 nm and 280 nm > 1.8, indicated that ctDNA was sufficiently free from protein [14]. The concentration of DNA (4.6 × 10−4 mol L−1) in stock solution used for binding experiments was determined by the absorption intensity at 260 nm using a molar absorption coefficient value of 6600 L mol−1 cm−1 [15]. All regents were of analytical grade, and doubly distilled water was used throughout.

Methods

Absorption Measurements

The UV absorption spectra were recorded on a TU-1810 spectrophotometer (Puxi Analytic Instrument Ltd., Beijing, China) by adding the increments of ctDNA stock solution into a fixed concentration of ADNRT. The solutions were allowed to incubate for 10 min before the absorption spectra were recorded. The ctDNA solutions of corresponding concentrations were measured as reference solutions. All UV measurements were carried out in Tris–HCl buffer (pH 7.4) at room temperature.

Fluorescence Measurements

All fluorescence measurements were performed on a Varian Cary Eclipse spectrofluorimeter (VARIAN, USA) equipped with a SHP DC-0515 thermostatic bath (Shanghai Hengping Sciences Instrument Co. Ltd., Shanghai, China), and a quartz cuvette of 1 cm path length. 2.0 mL pH 7.4 of Tris–HCl buffer solution, certain volume of drug and vary volume of ctDNA solution were transferred to a 10 mL volumetric flask, and the mixed solution was diluted to the final volume with doubly distilled water. These solutions were allowed to stand for 8 min to equilibrate. The fluorescence emission spectra were measured at 295, 303 and 310 K in the wavelength range of 520–750 nm with an excitation wavelength at 254 nm.

Viscosity Measurements

Viscosity measurements were carried out by Ubbelodhe viscometer, which was immersed in a thermostat water-bath at a constant temperature at 25 ± 0.1 °C. Various concentrations of ADNRT were then added into the viscometer to give a certain r (r = [drug]/[DNA]) value while the ctDNA concentration was constant. The flow time of the samples were repeatedly measured by a digital stopwatch with an accuracy of ±0.20 s after thermal equilibrium was achieved (15 min). The result of each point was the average of at least three readings. The data were presented as η/η 0 versus r, where η and η 0 are the viscosity of ctDNA in the presence and absence of ADNRT. Viscosity values were calculated from observing the flow time of ctDNA containing solutions (t > 180 s) and corrected by buffer solution (t 0), using η = (t- t 0)/t 0 [16].

Molecular Docking

The rigid molecular docking study was performed by Sybyl 6.9 software. The crystal structure of the B-DNA was extracted from the structure having Protein Data Bank identifier 425D. The structure of ADNRT was sketched and minimized with the help of the Trops force field using Gasteiger–Marsili charges. The macromolecule file was prepared for docking by adding polar hydrogen atoms and Gasteiger charges. After DNA was enclosed in the grid defined with 0.375 Å spacing, the grid map was calculated using the FlexX program. Visualization of the docked pose was done by using PyMol.

Results and Discussion

Interaction Between ADNRT and ctDNA

Spectral Characteristics of ADNRT Binding to ctDNA

It is known that the intrinsic fluorescence of DNA is of little practical use [17], whereas the titled molecule displays luminescent property (λ ex/λ em = 254/590 nm). Hence, the fluorescence emission spectra of ADNRT in the absence and presence of ctDNA were studied. As shown in Fig. 1, with the increasing amounts of ctDNA, the fluorescence intensity of ADNRT decreased without apparent shift of maximum fluorescence emission wavelength, which indicated the binding of ADNRT to ctDNA indeed existed.

The formation of ADNRT-ctDNA complex was further confirmed by UV absorption spectra (Fig. 2). The UV absorption spectra of ADNRT showed an intense absorption band, assigned to the π-π* transition of the anthraquinone chromophore [18]. It was apparent that as the concentration of ctDNA increased, the absorption peak at 230 nm decreased significantly with slight red shift (2 nm). These spectral changes suggested that ADNRT was likely to interact with ctDNA through some binding modes. Moreover, it is necessary to judge the binding mode of ADNRT with ctDNA by other methods.

Fig. 2

Absorption spectra of ADNRT in presence of ctDNA at different concentrations (versus corresponding concentration of ctDNA) at pH = 7.4; Conditions: C ADNRT = 8.3 × 10−7 mol · L−1, from 1 to 7: C ctDNA = 0, 4.6, 9.2, 13.8, 18.4, 23.0 and 27.6 × 10−6 mol · L−1

The Quenching Mechanism of Fluorescence of ADNRT by ctDNA

The mechanisms of fluorescence quenching are usually classified as dynamic quenching and static quenching, which can be distinguished by examination of the temperature on the Stern-Volmer equation (Eq. (1)):where F 0 and F are the fluorescence intensities in the absence and presence of ctDNA, respectively, K q is the quenching rate constant of biomolecular, τ0 is the average lifetime of the fluorophore and K SV is the Stern-Volmer quenching constant which can measure the fluorescence quenching efficiency.

It was assumed that the interaction of ADNRT with ctDNA proceeded via a dynamic quenching. The K SV value was obtained from the slope of the F 0/F versus [DNA] linear plot (Fig. 3a), and the values of K SV and K q at the three temperatures (295, 303 and 310 K) are listed in Table 1. As can be seen, the values of K SV decreased with the increasing temperature, indicating that the fluorescence quenching of ADNRT by ctDNA was static [19].

Fig. 3

a Stern-Volmer curves of ctDNA quenching the fluorescence of ADNRT. b Double logarithmic curves of ctDNA quenching the fluorescence of ADNRT

Table 1

Stern-Volmer quenching constants for the interaction of ADNRT with ctDNA at different temperatures

T(K)	Stern-Volmer equation	K SV(L · mol−1)	K q(L · mol−1 · s−1)	R
298	Y = 1.0002 + 2.209 × 104[DNA]	2.209 × 104	2.209 × 1012	0.9982
304	Y = 0.9915 + 1.809 × 104[DNA]	1.809 × 104	1.809 × 1012	0.9980
310	Y = 0.9899 + 1.389 × 104[DNA]	1.389 × 104	1.389 × 1012	0.9979

Binding Constant and Binding Sites

For a static quenching interaction, it is often assumed that the binding capability of DNA at each binding site is equal [20]. The binding constant K and the number of binding site n can be calculated by using the double logarithm regression equation,

The values of K and n were obtained from the intercept and slope of the plot of log[(F 0-F)/F] versus log[DNA] (Fig. 3b). The corresponding results at different temperatures are summarized in Table 2. The values of n approximately equal to 1 indicated that there was just a single binding site between ctDNA and ADNRT. The increasing trend of K with temperature indicated that the capacity of ADNRT binding to DNA was enhanced with the temperature rising [21].

Table 2

The binding constants, number of binding sites and thermodynamic parameters of the interaction of ADNRT with ctDNA at different temperatures

T(K)	K(L · mol−1)	n	R	ΔG(kJ · mol−1)	ΔH(kJ · mol−1)	ΔS(J · mol−1 · K−1)
295	1.319 × 104	0.953	0.9972	−23.616
303	3.980 × 104	1.076	0.9987	−25.911	61.020	286.900
310	4.273 × 104	1.108	0.9972	−27.919

Thermodynamic Parameters and Nature of Binding Forces

The acting forces between a small molecule and macromolecule mainly include hydrogen bonds, van der Waals force, electrostatic force and hydrophobic interaction force. To further characterize the interaction forces, the thermodynamic parameters which are the main evidence to determine the binding mode were analyzed. If the enthalpy change (ΔH) does not vary significantly within the range of temperature, the thermodynamic parameters ΔH and ΔS can be estimated using van’t Hoff plots [22]:

The free energy change (ΔG) was estimated from the Gibbs’ equation [23]:

The values of ΔH and ΔS were obtained from the slope (−ΔH/R, R is the gas constant) and intercept (ΔS/R). The values of ΔH, ΔS and ΔG at different temperature are given in Table 2. The negative values of ΔG revealed that the interaction processes was spontaneous, and the positive ΔH and ΔS values associated indicated that the hydrophobic interaction played a major role [24].

Effect of Ionic Strength on the Spectrum of ADNRT-ctDNA

In order to prove the binding mode of ADNRT with ctDNA, the effect of different concentrations of NaCl (from 0 to 0.25 mol L−1) on the fluorescence of ADNRT-ctDNA system and free compound was investigated. It represented that the fluorescence intensity was not dependent on ionic strength. When NaCl exists in the system, the electrostatic repulsion between the negatively charged phosphate skeletons on adjacent nucleotides is reduced with the increasing concentration of Na+ [25]. Apparently, the groove bound molecules can be released from the helix by increasing the ionic strength, whereas it is difficult for the intercalation bound molecules to be released, owing that a small molecule binding in the groove of DNA duplex exposes much more to the solvent surrounding than it does for the intercalation [26]. As seen from Fig. 4, the results demonstrated that the effect of ionic strength on the ADNRT–ctDNA system was very limited, so the interaction between ADNRT and ctDNA was intercalative binding.

Fig. 4

Effect of ionic strength on the fluorescence intensity of ADNRT and ADNRT-ctDNA system. C ADNRT = 0.1 × 10−6 mol · L−1; C ctDNA = 9.1 × 10−6 mol · L−1; C NaCl = 0.01, 0.02, 0.03, 0.04, 0.05, 0.10, 0.15, 0.20, 0.25 mol · L−1

Iodide Quenching Studies

Fluorescence quenching is a reliable method to study the binding of small molecules to DNA. A highly negatively charged quencher such as I- ion is expected to be repelled by the negatively charged phosphate backbone of DNA [27, 28]. Therefore, the molecules buried by intercalating into the helix will be protected from being quenched by I-, while the free aqueous complex and groove binding molecule should be quenched readily [9]. Figure 5a depicts the Stern-Volmer plots of the free ADNRT and the ADNRT-ctDNA system in the presence of KI. The values of K SV for free ADNRT and the ADNRT-ctDNA system were 419.23 L mol−1 and 303.85 L mol−1, respectively. It was obvious that the quenching effect for free ADNRT was stronger than that of the bound ADNRT, suggesting that an intercalation binding should be the interaction mode of ADNRT with ctDNA [29].

Fig. 5

a Fluorescence quenching plots of ADNRT by KI in the absence and presence of ctDNA. C ctDNA = 9.2 × 10−6 mol · L−1, C ADNRT = 3.3 × 10−6 mol · L−1; b Effect of dsDNA and ssDNA on the ADNRT fluorescence intensity, C ADNRT = 3.3 × 10−6 mol · L−1

Comparison of the Interaction of ADNRT with dsDNA and ssDNA

The behaviors of native DNA and denatured DNA were compared. Double-strand DNA was converted into single-strand DNA with the opening of its double helix by incubation at 100 °C for 30 min followed by rapid cooling in ice water [30]. Generally, thermal denaturation involves the rupture of hydrogen bonds and no covalent bonds are expected to be broken. If the interaction between ADNRT and DNA belonged to a groove binding mode, the extent of the fluorescence quenching of ADNRT would be stronger by ssDNA than that by dsDNA [31]. The results of the comparison experiments are given in Fig. 5b. The fluorescence quenching of ADNRT by ssDNA was smaller than that by dsDNA, which also supported that ADNRT intercalated into the helix stack.

Viscosity Measurements

Optical photophysical studies provide necessary but not sufficient clues to explain a binding between DNA and the complex, while hydrodynamic measurements that are sensitive to the length change are regarded as the least ambiguous tests of a binding model in solution [32]. Thus, viscosity measurements were carried out as an effective tool to further clarify the binding mode of ADNRT to ctDNA. An intercalator is generally known to cause a significant increase in the viscosity of a DNA solution due to lengthen the DNA helix as base pairs are separated to accommodate the binding ligand [33]. In contrast, a partial, non-classical ligand intercalation in grooves causes a bend in DNA helix reducing its effective length and thereby its viscosity [34]. As illustrated in Fig. 6, the relative viscosities of the ctDNA increased steadily upon the increasing concentrations of ADNRT. Such behavior further confirmed that ADNRT bound to DNA through an intercalative binding mode.

Fig. 6

Effect of increasing amount of ADNRT on the relatively viscosity of ctDNA

Molecular Docking

Molecular docking technique is an attractive scaffold to understand the drug-DNA interactions which can help to design and discover the rational drug. As shown in Fig. 7, ADNRT was successively docked with DNA duplex of sequence dodecamer (PDB ID: 425D) to search the optimal energy ranked results along with preferred orientation of the ligand inside the DNA groove. It revealed that the chromophore of ADNRT seemed to be mainly intercalated between guanine and cytosine base pairs of DNA. In addition, hydrogen bondings were prone to form between the amino group H of G or, the carbonyl O of C, and the intercalator. In this case, the hydroxy O of ADNRT was at a distance of 2.8 Å from azyl H atom of G10, the hydroxy H atom was 2.3 Å from the O2 of C15, the carbonyl O atom was 1.9 Å from the N-2 H atom of G10. The results of molecular modeling indicated that the interaction between ADNRT and ctDNA was intercalation dominated and hydrogen bonding forces might play an essential role in the binding. This result was consistent with the above spectroscopic and viscosity experiments.

Fig. 7

Computational docking model of the binding of ADNRT to ctDNA

Conclusion

The interaction between ADNRT and ctDNA was studied by spectral analysis, viscosity method and molecular docking. Hypochromic effects in UV and fluorescence spectroscopy showed that ADNRT could interact with ctDNA through static quenching mechanism. The binding constants and the number of binding sites of ADNRT binding to ctDNA were measured at different temperatures and the thermodynamic parameters were calculated as well. It was found that hydrophobic force played a significant role in the binding. The intercalative binding was much more reasonable by taking account of ionic strength effects, tests for the ability of ADNRT binding with dsDNA and ssDNA, KI quenching studies, and the viscosity measurements. The results obtained from molecular docking corroborated the experimental results obtained from spectroscopic investigations. This combination of spectroscopic and molecular modeling methods were expected to provide important insight into designing new antiviral anthracyclines.