Received date: February 25, 2016; Accepted date: March 08, 2016; Published date: March 10, 2016
Citation: Zhang H, Yang F, Wu Q, Si-Yan L, Liu P, et al. (2016) Synthesis, DNA-binding Properties and Quantitative Structure-Activity Relationships on Ruthenium(II) Complexes with Calf-Thymus DNA. Med chem (Los Angeles) 6:137-142. doi:10.4172/2161-0444.1000338
Copyright: © 2016 Zhang H, et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Visit for more related articles at Medicinal Chemistry
A series of novel ruthenium(II) complexes with electron-donor or electron-acceptor groups in intercalative ligands, [Ru(phen)2(o-MOPIP)]2+(1), [Ru(phen)2(o-MPIP)]2+(2), [Ru(phen)2(o-CPIP)]2+(3) and [Ru(phen)2(o-NPIP)]2+(4) have been synthesized and characterized with elementary, ES-MS, 1H NMR, electronic absorption and emission spectra. The binding properties of these complexes to CT-DNA have been investigated by spectroscopy and viscosity experiments. It’s illustrated that these complexes bind to DNA in a non-classical intercalation mode and their intrinsic binding constants (Kb) for 1, 2, 3 and 4 are calculated as 1.1, 0.35, 0.53 and 1.7 × 105 M-1, respectively. The Quantitative Structure-Activity Relationships (QSAR) of these ruthenium complexes, as well as some other ruthenium complexes congers has been investigated, and a linearity equation have been obtained: logKb=0.2429π+0.0429π2+0.2907σ+0.6 389I+4.3491 (n=12; R=0.9338; F=11.9134; p=0.0030). This results show that the electron-acceptor group and a large hydrophobic group will enhance the DNA binding affinity of ruthenium complexes.
Ruthenium; DNA; Quantitative structure-Activity relationships
Phen: 1,10-phenanthroline; o-MOP: o- (2- methoxylphenyl) imidazo [4,5-f][1,10] phenanthroline; o-MP: o- (2-methylphenyl) imidazo [4,5-f][1,10] phenanthroline; o-CP: o- (2-chlorophenyl) imidazo [4,5-f][1,10] phenanthroline; o-NP: 2-(2-nitrophenyl) imidazo [4,5-f][1,10] phenanthroline
For years, many attentions have been focused on the interaction of octahedral Ru (II) complexes with DNA owing to their potential utility as DNA probes, molecular light switches and chemotherapy drugs and photodynamic therapy for tumors [1-13]. For one thing, DNA has long been considered the main target for anticancer drugs. In general, Ru (II) complexes can bind to DNA in three non- covalent modes: intercalation binding, groove binding and electrostatic binding. It’s know that complexes with an enlarged aromatic ligand(intercalating ligand) can bind to DNA with high affinity(104~106), while those complex such as Ru(bpy)32+ can bind to DNA mostly in electrostatic .
For the last decade, a number of ruthenium complex with 2-phenylimidazo [4,5-f][1,10]-phenanthroline(PIP) and its derivates as intercalating ligand have been synthesized and their DNA binding properties have been investigated thoroughly. Ji et al. indicate that the binding affinity of these complexes depended not only on the conformations of DNA, but also on the structures of intercalative ligands [15-18]. The factors included the enlarged aromatic ring, intramolecular hydrogen bond and the planarity properties of intercalative ligand will enhance the binding affinity of these ruthenium to DNA, and the electronic effects (the donor/acceptor electron properties of substituent group on intercalative litgand) is also one factor influencing the bind of complexes to DNA [19,20].
A recently studies try to explain the DNA-binding affinity by computational calculations with density functional theory (DFT). It’s shown that the energy of these complex’ frontier molecular orbit, that is the highest occupied molecular orbit (HOMO) and the lowest unoccupied molecular orbit (LUMO) is varied and when they intercalating in the DNA base pairs, the energy of the transition conformation will different. According to the frontier molecular theory, an electron will transfer more easily from a high HOMO to a lower LUMO, and resulting those complex have the lowest LUMO (in generally, the HOMO of DNA is higher than that of ruthenium complexes) will bind to DNA the strongest [21-23]. But this is still confused since the binding energy of complexes to DNA is not known and the optimal conformation of the supermolecular complex- DNA is also not been illuminated.
In 1930’s, Hammett indicated that the activity of an organic reaction is in relating to the substituent effects [24-26]. Based on these, Fujita and Hansch developed a linear Hansch equation to elucidate the relationship between the bioactivity/physical activity and the structure of organic molecules , and which is so called the linear free energy relationships and thus have been utilized extensively in agrochemistry, pharmaceutical chemistry, toxicology  for its excellent predictable ability. More recently, Prasanna S et al. successful discerned the structural and physicochemical requirements for selective COX-2 over COX-1 inhibition among the fused pyrazole ring systems by Hansch method . However, there are still no reports focused on quantity structure-activity relationship on DNA-binding properties of ruthenium(II) complexes [30,31].
In this paper, a series of Ru (II) complexes with electron-donor or electron- acceptor substituents in the intercalative ligands, [Ru(phen)2(o-MOP)]2+1, [Ru(phen)2(o-MP)]2+2, [Ru(phen)2(o- CP)]2+3 and [Ru(phen)2(o-NP)] 2+4 (Scheme 1) were synthesized and characterized. The DNA-binding properties of these complexes have been investigated by the spectroscopic and viscosity experiments. The quantity structural-activity relationship of these ruthenium complexes, as well as some other analogues has also been investigated.
CT-DNA was purchased from the Sino-American Biotechnology Company. All reagents and solvents were purchased commercially (AR, Acros Inc., and Sigma Inc., etc.) and used without further purification unless otherwise noted. Doubly distilled water was used to prepare buffers. The concentration of calf thymus DNA was determined spectrophotomertrically using the molar absorptivity 6600 M-1·cm-1 (260 nm) (The ratio of UV absorbance at 260 and 280 nm is in the range of 1.8-1.9:1).
Synthesize and characterization
[Ru(phen)2Cl2]·2H2O were prepared following the literature procedure . Ru(II) complexes 1, 2, 3 and 4 were synthesized by refluxing Ru(phen)2Cl2 and o-MOP (o-MP, o-CP or o-NP) in ethylene glycol under an argon atmosphere with high yield. Each complex was obtained as a PF6- salt and purified with column chromatography.
Electrospray mass spectrometry (ES MS) has recently been shown to be a powerful tool for measuring the molecular mass of non-volatile and thermally unstable compounds . The ES-MS for 1, 2, 3 and 4 exhibits a fragment ion peak of (M+1PF6)+at 933.0, 917.0, 936.9 and 947.9 (m/z), respectively (Table 1). The fragment ion peaks of M2+of these complexes appear at 394.3, 386.4, 396.4 and 401.8 for 1, 2, 3 and 4, respectively, and the resolution of this peak for these complexes shows that the species is doubly charged and the isotopic distribution corresponds to the calculated one.
|No.||Substituent Group||Physical Parameter||logKb|
a. σ is the abbreviation of electron parameter. All data are cited from reference .
b. π is the abbreviation of hydrophobic parameter. All data are cited from reference .
c. I is the abbreviation of an indactive parameter. The data is 1 when there is a hydrogen bond in the intercalative ligand, or it will be 0.
d. These data is cited from reference .
e. These data is cited from reference .
f. These data is cited from reference 
g. These data is cited from reference .
Table 1: The intrinsic binding constant and the physical parameters of Ruthenium(II) complex.
The electronic absorption spectra of these Ru (II) complexes in Tris-buffer are characterized by an intense ligand-centered transition in the UV region and a metal-to-ligand charge transfer transition (MLCT) in the visible region. The lowest-energy absorption bands ascribed to the MLCT transitions for 1, 2, 3 and 4 are 455, 452, 453 and 453 nm respectively. The intense and sharp bands at 263, 263, 263 and 261 nm in UV region for 1, 2, 3 and 4, respectively, are attributed to the intraligand π→π* transition via comparison with the spectra of [Ru(bpy)3]2+. Little variation in the energy of the MLCT bands with the electronic effect in intercalative ligand was observed.
Ru (II) complexes 1, 2 and 3 emit fluorescence in Tris-buffer in the range of 500 – 700 nm at room temperature, with the maximum at 589, 588 and 589, respectively, and only a very weak fluorescence was observed for complexes 4 at the same conditions (the maximum is at 588 nm).
2-(2-methoxylphenyl) imidazo [4,5–f][1,10] phenanthroline(o- MOP) (1a): The ligand 2-(2-methoxylphenyl) imidazo [4,5-f][1,10] phenanthroline (o-MOP) was prepared by the method similar to that in reference , and with some modification.
A solution of phenanthraquinone (0.26 g, 1.2 mmol), o-anisaldehyde (0.24 g, 1.8 mmol) and ammonium acetate (1.9 g, 25 mmol) in 10 cm3 glacial acetic acid was refluxed for 2 hour. The cooled deep red solution was diluted with 25 cm3 water, and neutralized with ammonmium hydroxide. Then the mixture was filtered and the precipitates were washed with water and acetone, then dried and purified by chromatography over 60-80 mesh SiO2 using methanol as an eluent, yields: 0.35 g, 84%. Calculated for C20H14N4O·H2O (%): C: 69.7; H: 4.69; N: 16.3; Found(%): C: 69.3; H: 4.66; N: 16.2. ES-MS (in DMSO, m/z): 326.7 (calc. 326.4).
2-(2-methylphenyl) imidazo [4,5–f][1,10] phenanthroline(o- MP) (2a): O-MP was synthesized by the same method as above, but with phenanthraquinone (0.26 g, 1.2 mmol) and o-tolualdehyde (0.22 g, 1.8 mmole), yield: 0.31 g, 78%. Calculated for C20H14N4·H2O(%): C: 73.2; H: 4.91; N: 17.1; Found(%): C: 73.6; H: 4.97; N: 17.4. ES-MS (in DMSO, m/z): 310.6 (calc. 310.4).
2-(2-chlorophenyl) imidazo [4,5–f][1,10] phenanthroline(o-CP) (3a): O-CP was synthesized by the same method as above, but with phenanthraquinone (0.26 g, 1.2 mmol) and 2-chlorobenzaldehyde (0.25 g, 1.8 mmole), yield: 0.32 g, 76%. Calculated for C19ClH11N4·H2O(%): C: 65.4; H: 3.76; N: 16.1; Found(%): C: 65.6; H: 3.82; N: 16.4. ES-MS (in DMSO, m/z): 331.0 (calc. 330.8).
2-(2-nitrophenyl) imidazo [4,5–f][1,10] phenanthroline(o-NP) (4a): O-NP was synthesized by the same method as above, but with phenanthraquinone (0.26 g, 1.2 mmol) and 2-nitrobenzaldehyde (0.27 g, 1.8 mole), yield: 0.38 g, 87%. Calculated for C19H11N5O2·H2O(%): C: 63.5; H: 3.65 N: 19.5; Found(%): C: 63.9; H: 3.62 N: 19.1. ES-MS (in DMSO, m/z): 341.6 (calc. 341.3).
[Ru(phen)2(o-MOP)]2+(1): [Ru(phen)2 (o-MOP)]2+was synthesized by the literature [31,33] and with some modification. [Ru(phen)2Cl2]·2H2O (0.09 g, 0.17 mmol) and 1a (0.058 g, 0.17 mmol) were added to 10 cm3 ethylene glycol. The mixture was refluxed for 2 h under an argon atmosphere. The cooled reaction mixture was diluted with water (20 cm3) and filtered to remove solid impurities. The complex was then separated from soluble impurities by precipitation with NH4PF6. The precipitated complex was dried, dissolved in a small amount of acetonitrile, and purified by chromatography over alumina oxide using MeCN-toluene (2:1, v/v) as an eluent, yield: 0.16 g, 84%. Calculated for C44F12H30N8OP2Ru·2H2O(%): C: 47.4; H: 3.08; N: 10.1; Found(%): C: 47.7; H: 3.11; N: 10.4; 1H NMR (DMSO-d6, δ ppm): 9.32 (1H, d); 9.08 (1H, d); 8.77 (4H, d); 8.39 (4H, s); 8.21 (2H, d); 8.14 (2H, t); 8.12 (2H, d); 8.02(2H, 2d); 7.81 (6H, m); 7.74 (1H, t); 7.35 (1H, d); 7.24(1H, t); 4.04(3H, s); ES MS of the PF6- salt in MeCN: m/z 933.0(M+1PF6)+(calc: 932.8); 394.3 (M) 2+(calc: 393.9). Absorption UV-Vis, in water at pH 7.2 λmax(ε/104 M-1cm-1): 263(8.6), 455(1.7). No corrected emission maximum in water at pH 7.2: 589.4 nm.
[Ru(phen)2(o-NP)]2+(2): [Ru(phen)2(o-MP)]2+was prepared by the above-mentioned method but with 2a (0.056 g, 0.17 mmol); yield: 0.14 g, 77%. Calculated for C44F12H30N8P2Ru·2H2O(%): C: 48.1; H: 3.12; N: 10.2; Found (%): C: 48.4; H: 3.12; N: 10.6; 1H NMR (DMSO-d6, δ ppm): 9.06(2H, d); 8.78 (4H, d); 8.39 (4H, s); 8.17 (2H, d); 8.08 (2H, d), 8.04 (2H, d); 7.85(6H, m); 7.51(3H, m); 7.26(1H, t); 2.70(3H, s); ES MS of the PF6 - salt in MeCN: m/z 917.0(M+1PF6)+(calc: 916.8); 386.4 (M) 2+(calc: 385.9). Absorption UV-Vis, in water at pH 7.2 λmax(ε/104 M-1cm-1): 264(9.9), 452(1.7). No corrected emission maximum in water at pH 7.2: 587.8 nm.
[Ru(phen)2(o-CP)]2+(3): [Ru(phen)2(o-CP)]2+was prepared by the above-mentioned method but with 3a (0.059 g; 0.17 mmol); yield: 0.15 g, 79%. Calculated for C43F12ClH27N8P2Ru·4H2O(%): C: 44.8; H: 3.06; N: 9.71; Found (%): C: 44.5; H: 3.10; N: 9.81; 1H NMR (DMSO-d6, δ ppm): 9.03(2H, d); 8.77 (4H, d); 8.38 (4H, s); 8.13 (2H, d); 8.07 (2H, d), 8.00 (2H, d); 7.93(1H, d); 7.77(6H, m); 7.64(2H, m); 7.26(1H, m); ES MS of the PF6- salt in MeCN: m/z 936.9 (M+1PF6)+(calc: 937.2); 396.4 (M) 2+(calc: 396.2). Absorption UV-Vis, in water at pH 7.2 λmax(ε/104 M-1cm-1): 263(11.6) 454(2.1). No corrected emission maximum in water at pH 7.2: 589.4 nm.
[Ru(phen)2(o-NP)]2+(4): [Ru(phen)2(o-NP)]2+was prepared by the above-mentioned method but with 4a (0.061 g; 0.17 mmol); yield: 0.14 g, 74%. Calculated for C43F12H27N9O2P2Ru·2H2O (%): C: 45.8; H: 2.77; N: 11.2; Found (%): C: 45.5; H: 2.81; N: 11.8; 1H NMR (DMSO-d6, δ ppm): 8.87(2H, d); 8.77 (4H, d); 8.38 (4H, s); 8.21 (1H, d); 8.12 (2H, d), 8.10 (2H, d); 7.85(2H, d); 7.76 (6H, m); 7.68 (3H, m); ES MS of the PF-6 salt in MeCN: m/z 947.9(M+1PF6)+(calc: 947.8); 802.3 (M-H+)+(calc: 801.8); 401.8 (M) 2+(calc: 401.4). Resolution of the peak 401.8 shows that the species is double charged and the isotopic distribution corresponds to the calculated one. Absorption UV-Vis, in water at pH 7.2 λmax (ε/104 M-1cm-1): 262(12.2) 454(2.2). No corrected emission maximum in water at pH 7.2: 588.0 nm.
Microanalyses were carried out on an Elementar Vario EL elemental analyser. Electrospray mass spectra (ESI-MS) were recorded on a LCQ system (Finnigan MAT, USA). The spray voltage, tube lens offset, capillary voltage and capillary temperature were set at 4.50 kV, 30.00 V, 23.00 V and 200 °C, respectively, and the quoted m/z values are for the major peaks in the isotope distribution. Emission spectra were measured on a Shimadzu RF-5000 spectrofluorophotometer and UV- Visble absorption was recorded on a Shimadzu UVPC-3000 spectrophotometer. Viscosity experiments were performed on an Ulbbelodhe viscometer, immersed in a thermostatted water-bath maintained at 30.0 ± 0.1°C. Data were presented as (η/η0)1 / 3 vs. the concentration of [Ru]/[DNA]. Viscosity values were calculated from the observed flow time of DNA-containing solutions (t > 100 s) corrected for the flow time of buffer alone (t0), i.e., η=t-t0.
DNA-binding properties of Ru (II) complexes
Electronic absorption spectra: In general, the complex binding to DNA in an intercalation mode exhibits a red and hypochromism shift in the absorption spectra, and the extents of spectral change are closely correlative to the DNA-binding affinities of these complexes. The spectral shifts in an intercalation mode are usually greater than those in groove binding mode. In the presence of double helix calf thymus DNA (CT- DNA), the electronic absorption spectra for all of these complexes exhibit obviously hypochromism, and the hypochromism values for 1, 2, 3 and 4 at MLCT absorption band (452~455 nm) are 12, 9, 9 and 21%, respectively.
In order to clarify the DNA-binding affinities of these complexes, the intrinsic binding constants were calculated according to equation (1) , through a plot of
[DNA]/ Σa-Σf vs. [DNA]
[DNA]/ Σa-Σf=[DNA]/Σb-Σf+1/Kb( Σa-Σf) (1)
where [DNA] is the concentration of DNA in base pairs, Σa, Σf and Σb are respectively the apparent extinction coefficient (Aobsd/[M]), the extinction coefficient for free metal (M) complex and the extinction coefficient for the metal(M) complex in the fully bound form. In plots of [DNA]/ Σa-Σf versus [DNA], Kb is given by the ratio of slope to intercept. The calculated values for 1, 2, 3 and 4 at MLCT absorption band are 1.1, 0.35, 0.53 and 1.7 × 105 M-1, respectively. These values are smaller than those for [Ru(bpy)2dppz]2+(> 106 M-1)  and [Ru(ip)2dppz]2+(2.1 × 107 M-1) . Such DNA-binding constants suggest that the interaction of these complexes with DNA should be in an intercalation mode.
Emission spectra: The interaction of Ru (II) complexes with double helix CT-DNA was monitored via luminescence. All ruthenium complexes 1-4 emit luminescence in the range 500-700 with the maxium near 600 nm at room temperature. Upon the addition of CT-DNA, the emission spectra of all of these complexes enhanced obviously (Figure 1). The emission of complex 4 exhibits pronounced enhancement (Figure 2), and its emission intensity increases steadily to8.5 times relative to that of the original and reaches saturation at ca. [DNA] / [Ru]=8:1. However, the emission intensities increase by 2.2, 1.7 and 1.5 for complexes 1, 2 and 3, respectively (Figure 2). The enhancement of emission intensities of these complexes can be attributed to the hydrophobic environment inside the DNA helix, which reduces the accessibility of water molecules and makes the mobility of the complexes be restricted at the binding site.
Viscosity behaviors: The viscosity experiments, being sensitive to the change of length of double helix DNA, were considered as one of the most unambiguous methods to determine the binding mode of complex to DNA in absence of crystal data . In general, the relative viscosity of DNA in presence of complex in an intercalation mode will be increased, because the intercalative ligand will separate the base pairs of DNA, and thus lengthen the DNA helix. On the contrary, a partial and/or non-classical intercalation of complex will reduce the relative viscosity of DNA, since the binding ligand may bend (or kink) the DNA helix and reduce its effective length . The experiments on relative viscosity of rod-like CT-DNA in the presence of complexes 1, 2, 3 and 4, as well as [Ru(bpy)3]2+, were carried out and the results were shown in Figure 3.
The viscosity of DNA remains almost unchanged upon addition of [Ru(bpy)3]2+, which is consistent with an electrostatic association. However, in the presence of Ru (II) complexes 1, 2, 3 and 4 respectively, the relative viscosity of rod-like DNA was increased (Figure 3), because the stacking interaction of these complexes with the base pairs of DNA lengthens the DNA helix, indicating these complexes can bind to DNA in intercalation mode.
Quantitative structure-activity relationships on ruthenium(II) complexes
Quantitative Structure-Activity Relationships was carried out on MatLab 6.5 for these newly synthesized ruthenium(II) complexes, as well as some congers from references. The data of intrinsic binding constant and the physical parameters were listed in Table 1.
Firstly, we try to draw a plot with the logKb versus electronic parameter (σ) on Origin 6.0, as shown in Figure 4, and the corresponding equation is as follows:
LogKb=0.6726 (± 0.1751)σ+4.5001 (± 0.09799) (2)
n=7; r=0.8643; SD=0.1683; p=0.0121
Outlinear: X=-OCH3, Y=-H; X=-CH3, Y=-OH; X=-H, Y=-Phen; X=-H, Y=-OH; X=-H, Y=-OCH3
Considering there is hydrogen bond exists in some of these ruthenium complexes, we import indicative variable (IH-bonding) to indicate hydrogen bond exists in the intercalative ligand. The value of IH-bonding is 1 if there is hydrogen bond, regardless it’s intramolecular or intermolecular hydrogen bond, and the value of IH-bonding is 0 if there is not, thus we get model 2:
n=12; R=0.7441; F=5.5824; p=0.0265
It’s obviously the hydrogen bond contribute to the DNA binding of these ruthenium complexes, since the coefficient IH-bonding is positive. Encouraged, we considered that the hydrophobic parameter may also contribute to the DNA- binding properties of these complexes, thus we obtained the model 3:
n=12; R=0.9338; F=11.9134; p=0.0030
In this model 3, it’s obviously the coefficient of electronic parameter(σ) is positive, indicating the electron acceptor group on intercalative ligand of ruthenium complexes will enhance the binding affinity of ruthenium complexes to DNA, while an electron donor group will decrease the binding affinity. The positive coefficient for the hydrophobic parameter(π) indicate a hydrophobic group in the intercalative ligand will increase the DNA-binding affinity of ruthenium complexes, while hydrophilic group will decrease the DNA-g affinity.
A series of ruthenium(II) have been synthesized, and the bindingbehavior of these ruthenium(II) complexes with calf-thymus DNA have been investigated, and the results show that these complexes can bind to DNA in intercalating mode. The further studies on the quantity structure-activity relationship of these ruthenium complexes, as well as some from reference was investigated, and a QSAR equation was obtained: logKb=0.2429π+0.0429π2+ 0.2907σ+0.6389I+4.3491 (n=12; R=0.9338; F = 11.9134; p = 0.0030). It’s shown that the DNA-binding affinity of ruthenium complexes in studied depended on the electronic effect, hydrophobic effect and hydrogen bond, and an electron withdraw group in the intercalative ligand will increase the DNA-binding affinity of ruthenium complexes, while an electron donor group will decrease the DNA-binding affinity. In addition, hydrogen bond is important to obtain a high DNA-binding ruthenium complex.
This work was supported by the National Nature Science Foundation of China (81572926), the Provincial Major Scientific Research Projects in Universities of Guangdong Province (2014KZDXM053), the Science and Technology Project of Guangdong Province(2014A020212312), the Guangzhou City Science and Technology Plan (2013J4100072), the Science and Technology Projects of Yuexiu District(2014-WS-039), and the Joint Natural Sciences Fund of the Department of Science and Technology and the First Affiliated Hospital of Guangdong Pharmaceutical University (GYFYLH201309).