Jit Karmakar1, Promita Nandy1, Saurabh Das1, Debalina Bhattacharya2,3, Parimal Karmakar2, Samaresh Bhattacharya1. 1. Department of Chemistry, Inorganic Chemistry Section, Jadavpur University, Kolkata 700 032, India. 2. Department of Life Science and Biotechnology, Jadavpur University, Kolkata 700 032, India. 3. Department of Microbiology, Maulana Azad College, Kolkata 700 013, India.
Abstract
A family of three water-soluble half-sandwich arene-ruthenium complexes, depicted as C 1 -C 3 , having the general formula [Ru(p-cymene)(L)Cl]Cl has been synthesized, where L represents (1H-benzo[d]imidazol-2-yl)guanidine (L 1 ) or (benzo[d]oxazol-2-yl)guanidine (L 2 ) or (benzo[d]thiazol-2-yl)guanidine (L 3 ). The crystal structure of complex C 3 has been determined. The complexes show several absorption bands in the visible and ultraviolet regions, and they also show prominent emission in the visible region while excited near 400 nm. Studies on the interaction of ligands L 1 -L 3 and complexes C 1 -C 3 with calf thymus DNA reveal that the complexes are better DNA binders than the ligands, which is attributable to the imposed planarity of the ruthenium-bound guanidine-based ligand, enabling it to serve as a better intercalator. Molecular docking studies show that the complexes effectively bind with DNA through electrostatic and H-bonding interactions and partial intercalation of the guanidine-based ligands. Cytotoxicity studies carried out on two carcinoma cell lines (PC3 and A549) and on two non-cancer cell lines (BPH1 and WI-38) show a marked improvement in antitumor activity owing to complex formation, which is attributed to improvement in cellular uptake on complex formation. The C 1 complex is found to exhibit the most prominent activity against the PC3 cell line. Inclusion of the guanidine-based ligands in the half-sandwich ruthenium-arene complexes is found to be effective for displaying selective cytotoxicity to cancer cells and also for convenient tracing of the complexes in cells due to their prominent emissive nature.
A family of three water-soluble half-sandwich arene-ruthenium complexes, depicted as C 1 -C 3 , having the general formula [Ru(p-cymene)(L)Cl]Cl has been synthesized, where L represents (1H-benzo[d]imidazol-2-yl)guanidine (L 1 ) or (benzo[d]oxazol-2-yl)guanidine (L 2 ) or (benzo[d]thiazol-2-yl)guanidine (L 3 ). The crystal structure of complex C 3 has been determined. The complexes show several absorption bands in the visible and ultraviolet regions, and they also show prominent emission in the visible region while excited near 400 nm. Studies on the interaction of ligands L 1 -L 3 and complexesC 1 -C 3 with calf thymus DNA reveal that the complexesare better DNA binders than the ligands, which is attributable to the imposed planarity of the ruthenium-bound guanidine-based ligand, enabling it to serve as a better intercalator. Molecular docking studies show that the complexes effectively bind with DNA through electrostatic and H-bonding interactions and partial intercalation of the guanidine-based ligands. Cytotoxicity studiescarried out on two carcinomacell lines (PC3 and A549) and on two non-cancercell lines (BPH1 and WI-38) show a marked improvement in antitumor activity owing to complex formation, which is attributed to improvement in cellular uptake on complex formation. The C 1 complex is found to exhibit the most prominent activity against the PC3cell line. Inclusion of the guanidine-based ligands in the half-sandwich ruthenium-arenecomplexes is found to be effective for displaying selective cytotoxicity to cancercells and also for convenient tracing of the complexes in cells due to their prominent emissive nature.
Development
of a new series of therapeutic agents and modification
of any existing seriesare an essential and challenging aspect of
research related to the treatment of cancer. Platinum metal-based
species, especially coordination and organometalliccomplexes, are
widely employed as chemotherapeutics in combating this dreadful disease.[1] However, relatively low selectivity and adverse
side effects of majority of these species have led to new initiatives
toward development of better chemotherapeutic agents, particularly
of new platinum metal-based anticancercomplexes with minimal side
effects and high selectivity and cytotoxicity toward cancercells.Among the platinum metal-based species, rutheniumcompounds have
found a very special position owing to their prominent anticancer
activities.[2] Ruthenium-based molecular
speciesare found to be promising candidates for the development of
novel anticancer agents, primarily as they can bind DNA in several
possible modes, a property usually not found in platinum-based species.
Hence, the ruthenium-based species also have the potential to treat
platinum-resistant cancers. Among the different oxidation states of
ruthenium, the +2 state is most preferred for antitumor activity due
to stability of the ruthenium(II)complexes in vitro. Proper choice
of a ligand scaffold is crucial for inducing the desired DNA binding
and antitumor activity in the ruthenium(II)complexes. In this context,
half-sandwich ruthenium–arenecomplexesare particularly noteworthy.[3,4] The presence of the planararene moiety primarily blocks one face
of the complex and thus directs most of the reactivity toward the
other side. Besides, para-cymene, a heavily used
arene moiety, is known to cause distortion in DNA conformation that
eventually leads to DNA damage.[5] In half-sandwich
ruthenium–arenecomplexes, there is ample scope of varying
the ligand/ligand combination to occupy the remaining three coordination
sites on the metalcenter. In the present work, where our main objective
was to develop a new family of half-sandwich ruthenium–arenecomplexes for efficient DNA binding and anticancer activity, a group
of three guanidine-based compounds was selected (Chart ) as ancillary ligands. These
ligands have several important features, which are favorable for developing
efficient DNA binding and anticancer agents. They have multiple N–H
bonds, which are expected to favor solubility in water and binding
with pyrimidine or purine bases via H-bonding and thus mislead the
transcription process, resulting in DNA damage. In addition, the near-planar
geometry of these ligands may enable them to function as intercalating
moieties in the derived ruthenium complexes.[3c,3h,6] The skeleton of the selected ligands has
close resemblance with that of 4′,6-diamidino-2-phenylindole
(DAPI), a conjugated system that is known to efficiently stain DNA
present in a cell,[7] which further encouraged
us to use this group of ligands, with the hope that the mixed-ligand
half-sandwich complexes derived from them may exhibit luminescence
property and thus, they will be useful for measuring the capacity
for cellular uptake and track a molecule within a cell. Among the
three chosen ligands, synthesis of L and L and complexes
of L with few metals are known
in the literature,[8,9] while L is, to the best of our knowledge, new and coordination
chemistry of L and L appears to remain unexplored. The L ligand is known to coordinate metalcenters as a bidentate N,N-donor, forming a stable six-membered chelate
ring (I, X = NH),[9] and the L and L ligands are likely to display a similar mode of binding (I, X = O and S). As a source of ruthenium(II) and the arene
moiety the dimeric[{Ru(p-cymene)Cl2}2]compound was utilized. Reaction of the selected guanidine-based
ligands (L–L) with [{Ru(p-cymene)Cl2}2] indeed afforded half-sandwich ruthenium–arenecomplexescontaining the ligands L–L. Herein, we
describe the formation and characterization of these complexes, their
DNA binding properties, and their cytotoxicity toward selected cancercell lines.
Chart 1
Experimental Section
Materials
Ruthenium trichloride was
purchased from Arora Matthey, Kolkata, India. α-Phellandrene,
[Ru(bpy)3](ClO4)2, and dicyandiamide
were purchased from Sigma-Aldrich, USA. o-Phenylenediamine,
2-aminophenol, and 2-aminothiophenol were procured from Spectrochem,
India. The guanidine-based ligands (L–L) were synthesized
by reaction between o-phenylenediamine (or 2-aminophenol
or 2-aminothiophenol) and dicyandiamide following a reported protocol.[8] [{Ru(p-cymene)Cl2}2] was synthesized by following a published procedure.[10] All other chemicals and solvents were reagent-grade
commercial materials and were used as received.
Physical Measurements
Microanalyses
(C, H, and N) were performed on a Heraeus Carlo Erba 1108 elemental
analyzer. Magnetic susceptibilities were measured using a Sherwood
MK-1 balance. NMR spectra were recorded in CDCl3 solution
on Bruker Avance DPX 300 and 400 NMR spectrometers. IR spectra were
obtained on a PerkinElmer Spectrum Two spectrometer with samples prepared
as KBr pellets. Mass spectra were recorded with a Micromass LCT electrospray
(Qtof Micro YA263) mass spectrometer. Electronic spectra were recorded
on a PerkinElmer LAMBDA 25 spectrophotometer. Steady-state emission
spectra were collected on a PerkinElmer LS 55 fluorescence spectrometer,
and the quantum yields were determined by a relative method using
[Ru(bpy)3]2+ as the standard. Solution electrical
conductivities were measured using an Elico CM 180 conductivity meter
with a solute concentration of ca. 10–3 M. Geometry optimization by the density functional theory (DFT)
method and electronic spectral analysis by TDDFT calculation were
performed using the Gaussian 09 (B3LYP/SDD-6-31G) package.[11]
X-ray Crystallography
Single crystals
of [Ru(p-cymene)(L)Cl]Cl (C) were grown
by diffusion of diethyl ether vapor into a solution of the complex
in acetonitrile. Selected crystal data and data collection parameters
are given in Table . Data were collected on a Bruker SMART CCD diffractometer using
graphite monochromated Mo Kα radiation (λ = 0.71073 Å)
at 296 K. X-ray structure solution and refinement were done using
the SHELX-97 package.[12] H atoms were added
at calculated positions.
Table 1
Crystal Data and
Details of the Structure
Determination for Complex C
empirical formula
C18H24Cl2N4ORuS
formula mass
516.45
crystal system
orthorhombic
space group
P212121
a (Å)
7.5864(4)
b (Å)
12.9878(7)
c (Å)
21.6163(12)
V (Å3)
2129.9(2)
Z
4
Dcalcd (g/cm3)
1.611
F(000)
1048
crystal size (mm)
0.16 × 0.18 × 0.24
T (K)
296
μ (mm–1)
1.101
R1a
0.0452
wR2b
0.1175
GOFc
1.03
R1 =
Σ||Fo| – |Fc||/Σ|Fo|.
wR2 =
[Σ[w(Fo2 – Fc2)2]/Σ[w(Fo2)2]]1/2.
GOF = [Σ[w(Fo2 – Fc2)2]/(M – N)]1/2, where M is the number of reflections
and N is the number of parameters refined.
R1 =
Σ||Fo| – |Fc||/Σ|Fo|.wR2 =
[Σ[w(Fo2 – Fc2)2]/Σ[w(Fo2)2]]1/2.GOF = [Σ[w(Fo2 – Fc2)2]/(M – N)]1/2, where M is the number of reflections
and N is the number of parameters refined.
Synthesis of Complexes
The [Ru(p-cymene)(L)Cl]Clcomplexes (C–C) were synthesized
by following a general procedure as described below.The guanidine-based
ligand (L–L) (0.2 mmol) was dissolved in hot methanol
(50 mL). To it was added a solution of [{Ru(p-cymene)Cl2}2] (50 mg, 0.08 mmol) in dichloromethane (10 mL).
The resulting solution was heated at reflux for 4 h to generate a
yellowish-orange solution. The solvent was evaporated to almost one-fourth
of its initial volume, diethyl ether (50 mL) was added to it, and
the mixture was kept in a freezer for 12 h. Orange crystalline solid
was separated, which was collected by filtration, washed with dichloromethane
followed with diethyl ether, and dried in air.
Phosphate
buffer was prepared using NaH2PO4 and Na2HPO4 using triple-distilled water. Sodium nitrate
(AR) was used to maintain the ionic strength of the medium. Calf thymus
DNA, purchased from Sisco Research Laboratories, India, was dissolved
in triple-distilled watercontaining 120 mM NaCl, 35 mM KCl, and 5
mM CaCl2. Absorbance was recorded at 260 and 280 nm; A260/A280 was determined.
The ratio being between 1.8 and 1.9 suggests that the DNA was sufficiently
free from protein. It was characterized by measuring its CD spectrum
at 260 nm using a CD spectropolarimeter (J815, JASCO). Concentration
was determined in terms of nucleotide, considering the molar extinction
coefficient at 260 nm to be 6600 M–1 cm–1.50 μM L (ligand)
and 50 μM C (complex)
were titrated separately with calf thymus DNA at constant pH and ionic
strength of the medium. For the interaction of L followed by fluorescence spectroscopy, excitation
was done at 295 nm and emission was recorded at 331 nm. For the complex,
the excitation was done at 425 nm and emission was recorded at 496
nm. Ionic strength was maintained using NaCl and NaNO3.
The interaction of the compounds with DNA during titration led to
a decrease in fluorescence in the case of the ligand and to an increase
in fluorescence for the complex at the respective wavelengths where
they were followed. The interaction of compounds with DNA could be
realized with eq .[13−18]L represents compounds and Kd is the
dissociation constant for the interaction whose
reciprocal provides the apparent binding constant (Kapp). Equation is obtained from eq where the reciprocal of the change in absorbance was plotted
against the reciprocal of (CD – C0).[13−18]CD refers to concentration of calf thymus
DNA and C0 refers to concentration of
compounds. Using eq , ΔFmax could be determined along
with Kapp (1/Kd) from the intercept and slope.[13−18]ΔF represents the change
in fluorescence of the compounds interacting with the calf thymus
DNA during titration, while ΔFmax indicates the maximum possible change in fluorescence.ΔF/ΔFmax was plotted against CD. Equations and 4 were used to fit the data to a non-linear square
fit that provides another set of values for the apparent binding constant.[13−18] Titrations were also analyzed using a modified form of the Scatchard
equation [eq ].[19] The overall binding constant (K′) and site size (n) were determined.r denotes the ratio of the
concentration of the compound bound to DNA to the total concentration
of DNA present in the reaction mixture at any point of the titration
(Cb/CD); Cb is the concentration of the bound compound,
while Cf is that of the free compound.
“n” provides binding stoichiometry
in terms of the bound compound per nucleotide, while “nb” being the reciprocal of “n” denotes binding site size in terms of the number
of nucleotides bound to a compound. “nb” was obtained by plotting ΔF/ΔFmax against CD/[compound]. K′ may also be obtained
by multiplying Kapp with “nb” and is compared with values obtained
from a modified form of the Scatchard equation.
Molecular Docking Studies
Molecular
docking studies on complexesC1–C were performed using HEX 6.3 software to
identify possible binding sites in biomolecules. The three guanidine-based
ligands (L–L) were also docked using AutoDockTools 1.5.6
software. The coordinates of each rutheniumcomplex were taken from
its optimized structure as a .mol file and converted to a .pdb format
using PyMOL 2.4 software. The crystal structure of B-DNA (PDB ID: 1BNA) was retrieved from
the Protein Data Bank (http://www.rcsb.org./pdb). Visualization of the docked systems was performed using BIOVIA
Discovery Studio Visualizer (DSV) 2020 software. Default parameters
were used for docking calculations with the correlation type shape
only, FFT mode at the 3D level, and grid dimension of 0.6 with receptor
range 180 and ligand range 180 with twist range 360 and distance range
40.
Cell Culture
PC3, BPH1, A549, and
WI-38 lung fibroblast cells were cultured in RPMI or DMEM medium (GIBCO,
Invitrogen, Carlsbad, CA, US) supplemented with 10% fetal bovine serum
(GIBCO), 100 IU/mL penicillin, and 100 μg/mL streptomycin at
37 °C in a humidified atmosphere containing 5% CO2 (Heraeus, Thermo Scientific, MA, USA). All cell lines were procured
from the National Centre for Cell Science in Pune, India. Cells were
seeded in 96 well plates for 24 h prior to treatment with compounds.
MTT Assay
The antiproliferative
effect of three complexes and the guanidine-based ligands on four
cell lines, PC3, BPH1, A549, and WI-38 was determined by the MTT assay.
Cells were seeded at a density of 1 × 104 cells per
well in a 24-well plate. Next, the cells were exposed to the complex
and its ligand at different concentrations for another 24 h. After
incubation, cells were washed with 1× PBS twice. Thereafter,
they were treated with 0.5 mg mL–1 MTT solution
(SRL) and incubated for 3–4 h at 37 °C until an insoluble
purple-colored formazan product developed. The resulting product was
dissolved in MTT extraction buffer and the OD was measured at 570
nm using a microplate reader (Epoch). The percentage survival was
calculated considering the untreated cells as 100%.
Single-Cell Gel Electrophoresis/Comet Assay
A comet
assay was performed after treating the cells (PC3 and BPH1)
with complex C for 24 h at the
lower (20 μM) and higher (80 μM) concentrations of the
IC50 dose of 39.5 μM for PC3. Briefly, 1 × 105 cells mL–1 were mixed with 0.7% LMPA and
embedded onto frosted slides. The slides were then dipped in a lysis
solution [2.5 M NaCl, 100 mM EDTA, 10 mM Tris–HCl (pH 10)]
that contains freshly added 1% Triton-X 100 and 10% DMSO and incubated
for 1 h at 4 °C and placed into a horizontal electrophoresis
tank filled with freshly prepared buffer (1 mM EDTA, 300 mM NaOH).
Next, electrophoresis was performed for 20 min at a fixed voltage
of 25 V and 300 mA. After that, slides were washed with a neutralization
buffer (0.4 M Tris–HCl, pH 7.5) followed by staining with 20
mg mL–1 ethidium bromide (SRL, India) for 15 min.
The slides were then washed three times with 1× PBS and observed
under a fluorescence microscope (model: Leica, Germany). Around 50
comets per slide were counted for both the cell lines. An extension
of each comet was analyzed using a computerized image analysis system
(Kometsoftware 5.5) that gave % of tail DNA.[20]
DAPI Staining
After exposure with
complex C at lower and higher
concentrations of the IC50 dose of PC3 (39.5 μM)
for 24 h, both the cells (PC3 and BPH1) were washed several times
with 1× PBS and stained with 0.2 mg mL–1 DAPI
in Vecta shield (Vector Laboratories Inc.). The percentage of cells
with ruptured and decondensed nuclei was counted under a fluorescence
microscope (Leica) and photographs were taken at 40× magnification.
Measurement of Intracellular Reactive Oxygen
Species
The production of intracellularreactive oxygen species
(ROS) was estimated using a fluorescent dye, DCFDA. Approximately,
3 × 105 cells per well were seeded in 35 mm plates,
and after 24 h of seeding, cells were incubated with 20 μM DCFDA
(Sigma) dye for 1 h at 37 °C under dark conditions followed by
the treatment of complex C for
24 h, at the lower (20 μM) and higher (80 μM) concentrations
of the IC50 doses of PC3 (39.5 μM). Cells without
the complex were used as control. Fluorescence intensity was measured
in a fluorescence spectrophotometer (model Hitachi, USA) at excitation
and emission wavelengths of 504 and 529 nm, respectively. To nullify
the autofluorescence of the complex which may interfere with the DCFDA
dye, a set of experiments without cells were performed simultaneously.
Cellular Imaging Study
Both the
cell linesPC3 and BPH1 were seeded in a cover slip for overnight.
Next day, cells were incubated with 3 mM complex C for 1 h in 37 °C in a CO2 incubator.
After incubation, cells were washed several times with 1× PBS
under dark conditions. Cells were then stained with DAPI in Vecta
shield and observed under a fluorescence microscope.
Results and Discussion
Synthesis and Characterization
As
delineated in the Introduction, the first
task of the present study was to synthesize a group of arene–rutheniumcomplexes using the chosen guanidine-based ancillary ligands (L–L). Accordingly, reactions of these ligands (L–L) with [{Ru(p-cymene)Cl2}2] were carried out in 5:1 methanol/dichloromethane, which furnished
the desired complexes of type [Ru(p-cymene)(L)Cl]Cl
in decent yields. The three complexes obtained with ligands L, L, and L are depicted,
respectively, as C, C, and C. Preliminary characterization (microanalysis, mass, IR, and NMR)
data of these complexes agreed well with their compositions. In order
to ascertain the coordination mode of the guanidine-based ligands
in these complexes, the crystal structure of C was determined by X-ray crystallography.[21] The structure is presented in Figure , and some selected bond distances
and anglesare provided in Table . The structure reveals that the guanidine-based ligand
(L) is coordinated to ruthenium
as a neutral N,N-donor, forming a six-membered chelate ring (I, X = S and M = Ru) with a bite angle of ∼83°.
The Ru(L) fragment of the complex
is found to be nearly planar, as envisaged. The p-cymene moiety is bound to ruthenium in the usual π-fashion,
and a chloride ion has taken up the sixth coordination site on the
metalcenter. Another isolated chloride ion was located outside the
coordination sphere. The bond parameters around ruthenium and within
the coordinated organic ligands are all found to be quite usual.[3,9] Based on the similarity of the synthetic method and properties (vide
infra), the other two complexes (C and C) are assumed to
have similar structures as C.
Figure 1
Crystal structure of the [Ru(p-cymene)(L)Cl]+ complex.
Table 2
Selected Bond Distances and Bond Angles
of Complex C
Bond
Distances (Å)
Ru1–Cl1
2.4309(16)
N2–C18
1.377(8)
Ru1–N1
2.106(4)
N3–C18
1.308(8)
Ru1–N3
2.074(5)
N4–C18
1.340(9)
N1–C16
1.406(8)
S1–C15
1.734(7)
N1–C17
1.316(7)
S1–C17
1.739(6)
N2–C17
1.356(8)
Bond
Angles (deg)
N1–Ru1–N3
82.83(19)
Cl1–Ru1–N1
86.13(16)
Cl1–Ru1–N3
88.50(17)
Crystal structure of the [Ru(p-cymene)(L)Cl]+ complex.
Spectral
Studies
Magnetic susceptibility
measurements show that the C–Ccomplexesare diamagnetic,
which is consistent with the +2 oxidation state of ruthenium (low-spin
d6S = 0) in them. In the 1H NMR spectra of the complexes, signals from both the coordinated p-cymene and guanidine-based ligand were expected, majority
of which could be identified. For example, all the signals for the p-cymene ligand could be distinctly observed in all three
complexes. Three signals from the alkyl groups are observed within
1.1–3.5 ppm and two signals from the aromatic fragment are
observed within 5.2–5.8 ppm. From the guanidine-based ligands,
the NH and NH2 signals appeared within 6.5–8.2 ppm,
while signals from the aromatic protons are observed around 7.5 ppm.
In complex C, a broad signal
is observed at 10.50 ppm, which is absent in the spectra of the other
two complexes, and hence, it is attributable to the benzimidazole-NH
in metal-bound L. 13C NMR spectra of the complexesare also found to be consistent with
their compositions. For the p-cymene ligand, three
signals from the alkyl carbonsare found below 40 ppm and four from
the aromaticcarbons appear within 70–90 ppm. For the guanidine-based
ligands, two signals are observed above 150 ppm and signals from the
phenyl ring are found within 110–130 ppm.The mass spectra
of complexesC–C, recorded in the positive ion mode,
provide proof of coordination of the guanidine-based ligands. Each
complex shows a peak at a m/z value
that corresponds to the [Ru(p-cymene)(L-H)]+ fragment, which is generated via loss of HCl from the cationic [Ru(p-cymene)(L)Cl]+ unit. Associated with loss of
the coordinated chloride ion, loss of proton is believed to take place
from the central −NH– moiety of the guanidine-based
ligand. This particular proton in such ligands is known to undergo
facile dissociation.[22] Elimination of one
equivalent HCl from compositionally similararene–rutheniumcomplexes is precedent.[22] A similar mass
spectral behavior of complexesC–C supports their similarcomposition and structure.Infrared spectra of complexesC–C, recorded in
the 450–4000 cm–1 region, exhibit several
bands. Upon comparison of the spectrum of each complex with that of
the starting [{Ru(p-cymene)Cl2}2]complex reveals the presence of several new bands (near 3400, 3180,
1680, 1615, 1256, and 752 cm–1) in the spectra of
the complexes, which are attributable to the coordinated guanidine-based
ligand. Among these bands, the two near 3400 and 3180 cm–1 are attributable to the −NH– and −NH2 fragments, respectively. The NMR and IR data are therefore in good
agreement with the composition of the complexes.The C–Ccomplexesare soluble in polar solvents,
such as water, methanol, ethanol, dimethylformamide, and dimethylsulfoxide,
producing yellow solutions. Electronic spectra of the complexes were
recorded in methanol solutions. Spectral data are presented in Table . Each complex shows
four absorptions spanning the visible and ultraviolet regions. To
have an insight into the nature of these absorptions, TDDFT calculations
have been performed on the C–Ccomplexes, using
the Gaussian 09 package,[11] and the results
are found to be similar for all the complexes. The DFT-optimized structures
of the complexesare shown in Figure S1 (Supporting Information) and some computed bond parameters are listed in
Table S1 (Supporting Information). The
computed bond parameters in the DFT-optimized structure of Care comparable with those found in its
crystal structure. The main calculated transitions for the C–Ccomplexes and compositions of the molecular orbitals associated
with the transitions are presented in Tables S2–S7 (Supporting Information), and contour plots of
selected molecular orbitals are shown in Figure S2 (Supporting Information). As the computed optical transitions
and compositions of the participating orbitals are similar for all
three complexes, the case of C is described here as representative. Plots of experimental and theoretical
spectra for Care deposited
in Figure S3 (Supporting Information).
The close match between each set of experimental and theoretical spectra
testifies validity of the optimized structures of the complexes, particularly
of complexesC and C, for which crystal structures remained
elusive. The lowest energy absorption at 451 nm is attributable primarily
to a HOMO – 1 → LUMO transition, with much less HOMO
– 3 → LUMO, HOMO – 2 → LUMO and HOMO →
LUMO character. Additionally, based on the nature of the participating
orbitals, the electronic excitation is best described as a MMCT transition
mixed with some MLCT, LMCT, and LLCT character. The next absorption
at 294 nm is mostly due to a HOMO – 1 → LUMO + 2 transition
and assignable primarily to a MLCT transition with much less LLCT
and ILCT character. The third absorption at 237 nm is largely due
to a HOMO – 3 → LUMO + 2 transition and has a dominant
MLCT character. The fourth absorption at 211 nm has a dominant HOMO
– 3 → LUMO + 4 character and is assignable to a MLCT
transition with some LLCT character.
Table 3
Electronic
Absorption and Emission
Spectral Data of the Complexes
absorption spectral dataa
emission spectral dataa
complex
λmax, nm (ε, M–1 cm–1)
λF, nm [ΦF × 10–3]c
life
time
(τ)
C1
451(400), 294 (5940), 237b, 211 (21,970)
518 [8.6]
τ1 = 0.26 ns, τ2 = 4.10 ns
C2
449 (530), 282 (9190), 239b, 209 (22,770)
558
[6.5]
τ1 = 1.75 ns
C3
429 (490),
289 (7780), 255b, 220 (28,780)
486 [17.3]
τ1 = 0.36 ns, τ2 = 5.02 ns
In methanol.
Shoulder.
Quantum yield
was calculated with
reference to [Ru(bpy)3]2+ (ΦF = 0.09).
In methanol.Shoulder.Quantum yield
was calculated with
reference to [Ru(bpy)3]2+ (ΦF = 0.09).Luminescence
property of the complexes was examined in methanol
solution. All three complexes were found to show prominent emission
when excited near 400 nm (Table ). It is interesting to note that the complexes absorb
and emit in the visible region, a property much sought after in an
antitumor agent for its easy identification in a biological matrix.
DNA Binding Studies
The interaction
of the calf thymus DNA with complex C, a representative of this family of complexes, was studied
in detail to assess its potential as an antitumor agent. Initially,
titration of the uncoordinated guanidine-based ligand L was carried out with the calf thymus DNA. Figure A shows a typical
plot of this titration at the ionic strength and pH mentioned. Figure B is a similar plot
for a titration performed with complex C. The plots show that saturation is achieved in the binding
of the compounds with DNA. Representative plots based on eq are shown in Figure S4 (Supporting Information), from which the apparent
binding constant (Kapp) was evaluated
(Table ).[13−18] Plots in Figure A,B were fitted by the non-linear square fit analysis that also helps
to evaluate Kapp. Different binding parameters
are shown in Table . The inset of Figure A,B provides nb, the number of nucleotides
bound to each compound (Table ).[13−18] It is worth mentioning that the value for nb obtained in the case of the complex binding to the calf thymus
DNA was approximately 1.5 times greater than that obtained when the
guanidine-based ligand L binds
to the same DNA, suggesting that the complex engages more nucleic
acid bases when it interacts with DNA, thus being able to bring about
more distortion in DNA, an outcome of enforced planarity of the guanidine
portion of the ligand following chelation to ruthenium (Table ).[15,16,18]
Figure 2
Fluorescence emission spectra of (A) 50 μM L and (B) 50 μM C in aqueous solution in the presence of
0.12
M NaCl and 30 mM Tris buffer (pH 7.4) in the absence (1) and presence
of different concentrations of calf thymus DNA; temperature = 300
K.
Table 4
Binding Constant Values Obtained for
the Interaction of Ligand L and
Complex C1 with the Calf Thymus DNA that was
Followed by Fluorescence Spectroscopy
Kapp
K* = Kapp × nb
compound
expt
from double-reciprocal plot (a)
from non-linear plot
(b)
from double-reciprocal plot with y-intercept = 1(c)
site size nb
from double-reciprocal plot (a)
from non-linear plot
(b)
from double-reciprocal plot with y-intercept = 1(c)
K* from Scatchard
nb from Scatchard as nb = (n–1)
L1
1
2.80 × 103
3.30 × 103
1.60 × 103
8
2.20 × 104
2.60 × 104
1.28 × 104
1.86 × 104
8
2
2.20 × 103
2.80 × 103
2.20 × 103
9
1.98 × 104
2.50 × 104
1.98 × 104
3.60 × 104
7
C1
1
0.94 × 104
0.80 × 104
0.70 × 104
13
1.20 × 104
1.04 × 104
0.91 × 104
1.53 × 104
15
2
1.60 × 103
1.04 × 103
0.78 × 103
14
2.20 × 104
1.40 × 104
1.09 × 104
2.15 × 104
14
Figure 3
Binding isotherms for (A) ligand L and (B) complex C binding
to the calf thymus DNA at pH ∼ 7.4 and an ionic strength of
0.12 M. Corresponding non-linear fits are shown for these titrations
that evaluate Kapp. Inset: plot of the
normalized increase in fluorescence as a function of the ratio of
the calf thymus DNA to (A) ligand L and (B) complex C. [L] = [C] = 50 μM, pH = 7.40, T = 300 K.
Fluorescence emission spectra of (A) 50 μM L and (B) 50 μM C in aqueous solution in the presence of
0.12
M NaCl and 30 mM Tris buffer (pH 7.4) in the absence (1) and presence
of different concentrations of calf thymus DNA; temperature = 300
K.Binding isotherms for (A) ligand L and (B) complex C binding
to the calf thymus DNA at pH ∼ 7.4 and an ionic strength of
0.12 M. Corresponding non-linearfitsare shown for these titrations
that evaluate Kapp. Inset: plot of the
normalized increase in fluorescence as a function of the ratio of
the calf thymus DNA to (A) ligand L and (B) complex C. [L] = [C] = 50 μM, pH = 7.40, T = 300 K.Utilizing Kapp and nb from Table and the relation Kapp × nb = K′, the overall
binding constants could be evaluated for the uncoordinated guanidine-based
ligand L and complex C.[13−18] Overall binding constants were also obtained from a modified version
of the original Scatchard equation (eq ),[19] and plots obeying this
equation are shown in Figure S5 (Supporting Information). The overall binding constant values from the Scatchard equation
were strikingly similar to those evaluated by multiplying Kapp with nb (eqs and 4).An interesting aspect regarding titration of the complex
with the
calf thymus DNA, followed by fluorescence spectroscopy, was that with
an increase in the concentration of DNA (Figure ), there was a gradual increase in fluorescence
similar to that observed for compounds including ethidium bromide
that are known to intercalate DNA.[23−26] Hence, a logical conclusion is
that the complex too is able to intercalate DNA, registering an increase
in fluorescence. Such an increase in fluorescence upon interactions
is an important attribute of the complex that may be utilized in a
number of biological experiments to realize possible interactions
of the compound with a biological target.
Figure 4
Gradual variation in
fluorescence observed for the compounds as
the calf thymus DNA was added during titration; (■) complex C, (•) ligand L. Ionic strength of medium = 0.12 M; pH
∼ 7.4; [L] = [C] = 50 μM; temperature = 300 K.
Gradual variation in
fluorescence observed for the compounds as
the calf thymus DNA was added during titration; (■) complex C, (•) ligand L. Ionic strength of medium = 0.12 M; pH
∼ 7.4; [L] = [C] = 50 μM; temperature = 300 K.
Molecular Docking with
DNA
To elucidate
the mode of interaction and binding affinity, docking studies were
performed on B-DNA (PDB ID: 1BNA) in the presence of all the three complexes. The results
show that the complexes interact with DNA quite similarly via the
electrostatic mode. This is illustrated in Figure for complex C and in Figures S6 and S7 (Supporting Information) for complexesC and C, respectively. In each case, the
guanidine-based ligand is observed to form H-bonds with oxygen atoms
of the phosphate backbone and also with the N3 atom of adenine of
a single DNA strand. Additional H-bonding is observed in complex C due to the presence of an NH fragment
in L that is absent in L or L. The coordination-induced planarity of the guanidine-based
ligands is found to favor strong H-bonding interactions, with better
match of the complexes inside DNA strands allowing partial intercalation.
Due to the combined effect of the van der Waals and H-bonding interactions,
the complexes fit comfortably into the minor groove of the targeted
DNA near the A–T rich regions.
Figure 5
(a) Complex C interacted
with the DNA strand and (b) core view of the interaction (ball and
stick model).
(a) Complex C interacted
with the DNA strand and (b) core view of the interaction (ball and
stick model).Docking of the individual guanidine-based
ligands with DNA has
also been looked into. From the docked structures (Figure S8; Supporting Information), it is observed that
the L ligand shows the highest
binding affinity to DNA, which is attributable to additional H-bonding
possible due to the presence of an NH fragment in L, instead of oxygen (in L) or sulfur (in L). The same trend is observed in the complexes, which is also manifested
in the biological studies. It is interesting to note that while all
the uncoordinated guanidine-based ligands preferred to approach the
G–C base pairs, upon binding to the metalcenter, the A–T
base pairs have become their preferred binding location. Planarity
of the guanidine-based ligands in the complexes and the presence of
the Ru-coordinated p-cymene probably have caused
this observed variation in their binding preference.
Cytotoxicity Studies
Cisplatin shows
a remarkable efficacy in treating prostate cancer and has been quite
successfully and extensively used in the last few decades.[27−29] However, as delineated in the Introduction, ruthenium-based molecular species, particularly the half-sandwich
ruthenium–arenecomplexes, are also attracting attention owing
to their demonstrated anticancer activities with minimal side effects.
Encouraged by the prominent DNA binding properties of our three complexes
(C, C, and C), we
also determined the potency of these three complexes and cisplatin
on the humanprostate cancercell line PC3 and the human benign prostate
tumor cell line BPH1. Similarly, we have evaluated the toxic effect
of these three complexes and cisplatin on the lung cancercell line
A549 and the normal lung fibroblast cell line WI-38. Cells were treated
with three complexes (C, C, and C) in the concentration range of 0–100 μM for 24
h, followed by MTT assay. The results are displayed in Figure . Complex C was found to be the most cytotoxic to PC3cells
(IC50 = 39.5 ± 1.57 μM) among the three complexes
(Figure A). Complex C was found to be non-toxic to the
human prostate benign tumorcell line BPH1 even after 24 h of treatment,
which suggests no side effects of it on non-carcinomacells in our
body. In this context, it is worth mentioning that cisplatin shows
comparable cytotoxicity toward both PC3 and BPH1 cell lines. In A549
and WI-38cell lines, C shows
moderate (IC50 values 69.4 ± 1.2 and 69.6 ± 3.45)
and almost comparable cytotoxicity like cisplatin (IC50 values 60.1 ± 2.43 and 66.5 ± 2.12). Guanidine-based ligands
(L, L, and L) have
no cytotoxicity toward any type of cell lines, which signifies the
effect of their coordination to ruthenium in antiproliferative activity. Figure B shows the % cell
survival comparison between BPH1 and PC3 with Ccomplex for 24 h. The IC50 doses for all
the complexes and cisplatinare summarized in Table . Interestingly, it was observed that the
complexescould not exert significant toxicity toward A549 and WI-38.
Among the different cells we have tested, only PC3 is PTEN-negative.
Thus, it is reasonable to state that cytotoxicity generated through
the compounds is presumably governed by nonfunctional PTEN. Similar
results are also found in the case of curcumin.[30] All the other biological studies were done on two cell
lines, PC3 and BPH1, taking C as a model complex.
Figure 6
(A) MTT assay on the PC3 cell line after 24 h of treatment
with
three separate complexes C–C. (B) MTT assay on PC3 and BPH1
cell lines after 24 h of treatment with the C complex. Data are presented as % survival relative
to the untreated control. They are the mean ± SD of three independent
experiments.
Table 5
IC50 Values
of Ligands
(L) and Complexes toward Different Cell Linesa
complex
PC3
BPH1
A549
WI-38
R1b
R2c
L1
121.0 ± 1.57
>500 ± 4.5
>200 ± 2.2
>500 ± 4.67
nd
nd
L2
168.8 ± 1.9
>500 ± 4.09
167.4 ± 1.77
>500 ± 5.03
nd
nd
L3
446.4 ± 1.05
>500 ± 2.32
>200 ± 1.86
>500 ± 4.41
nd
nd
C1
39.5 ± 1.57
263.0 ± 1.87
69.4 ± 1.2
69.6 ± 3.45
6.6
1.0
C2
267.3 ± 2.01
443.9 ± 1.04
168.6 ± 1.5
135.9 ± 3.21
1.6
0.8
C3
125.0 ± 1.43
175.2 ± 1.88
112.2 ± 1.9
135.7 ± 4.89
1.4
1.2
cisplatin
5.4 ± 1.93
8.0 ± 1.03
60.1 ± 2.43
66.5 ± 2.12
1.4
1.1
The drug treatment period was 24
h.
R1 = IC50 ratio
of BPH1 cells to PC3 cells.
R2 = IC50 ratio
of WI-38 cells to A549 cells.
(A) MTT assay on the PC3cell line after 24 h of treatment
with
three separate complexesC–C. (B) MTT assay on PC3 and BPH1
cell lines after 24 h of treatment with the Ccomplex. Data are presented as % survival relative
to the untreated control. They are the mean ± SD of three independent
experiments.The drug treatment period was 24
h.R1 = IC50 ratio
of BPH1 cells to PC3cells.R2 = IC50 ratio
of WI-38cells to A549cells.The cytotoxic effect of the Ccomplex is most likely linked to the DNA-damaging effects of the
compound, and hence, we performed comet assay, a very useful and sensitive
experiment for elucidating single- or double-strand DNA damage caused
by any exogenous or endogenous species.[31] A small amount of nuclides/cells is required to perform this experiment
and the tail length is considered to represent the level of DNA damage.[32] As shown in Figure A, the Ccomplex caused significantly (p < 0.05) more
DNA damage at a concentration C50 or >IC50 dose in the case of PC3. However, at the same concentrations, no
significant (p > 0.05) DNA damage was observed
when
BPH1 is used. Consistent with these data, the percentage of tail DNA
increased significantly (p < 0.05) for Ccomplex-treated PC3cells after 24 h of
treatment (Figure B).
Figure 7
(A) Representative images of the comet assay of PC3 and BPH1 cell
lines treated with the C complex
with respect to the untreated control. (B) Histogram shows % of comet
tail DNA for PC3 and BPH1 cells treated with the C complex for 24 h with respect to their untreated
control at two different doses (IC50 doses of PC3). Values are the mean ± SD of three independent
experiments. *(p < 0.005) and **(p < 0.005) denote the statistically significant difference compared
to the untreated control.
(A) Representative images of the comet assay of PC3 and BPH1 cell
lines treated with the Ccomplex
with respect to the untreated control. (B) Histogram shows % of comet
tail DNA for PC3 and BPH1 cells treated with the Ccomplex for 24 h with respect to their untreated
control at two different doses (C50 and >IC50 doses of PC3). Valuesare the mean ± SD of three independent
experiments. *(p < 0.005) and **(p < 0.005) denote the statistically significant difference compared
to the untreated control.
Nuclear morphology and the nature of cell death were studied by
DAPI staining. The fluorescence micrographs of DAPI-stained PC3 and
BPH1 cell linesare shown in Figure A, and the percentage of apoptoticcells is presented
graphically in Figure B. When PC3cells were treated with the Ccomplex at a concentration >IC50 dose
for
24 h, we found that 80–85% cells were apoptotic in nature,
whereas for BPH1, the amount of apoptoticcells reduced to only 30–35%.
Figure 8
(A) Fluorescence
micrographs of DAPI-stained PC3 and BPH1 cell
lines under 40× magnification. Both the cells are treated with
the C complex for 24 h at two
different doses (IC50 doses
of
PC3). The arrow represents the decondensed nucleus of the apoptotic
cells. (B) % of apoptotic cells as determined by DAPI staining followed
by fluorescence microscopic observations. Each value represents the
mean ± S.D. of three independent experiments. *(p < 0.05) and **(p < 0.005) denote the statistically
significant difference compared to the untreated control.
(A) Fluorescence
micrographs of DAPI-stained PC3 and BPH1 cell
lines under 40× magnification. Both the cells are treated with
the Ccomplex for 24 h at two
different doses (C50 and >IC50 doses
of
PC3). The arrow represents the decondensed nucleus of the apoptoticcells. (B) % of apoptoticcells as determined by DAPI staining followed
by fluorescence microscopic observations. Each value represents the
mean ± S.D. of three independent experiments. *(p < 0.05) and **(p < 0.005) denote the statistically
significant difference compared to the untreated control.
Ru(II)–arenecomplexesare well known to bring about
cell
damage via production of ROS within the cells.[33] It is interesting to note that among the other metals (such
as Pt, Pd, and Au) used as therapeutics, only Ru shows higher antitumor
activity mediated by an enhanced ROS production.[34] Apoptoticcell death and DNA damage are connected with
ROS production, and we also estimated ROS production induced by the Ccomplex in PC3 and BPH1, where
we have used a fluorescent dye, DCFDA (2′,7′-dichlorofluorescin
diacetate), for indicating oxidative stress and hydroxyl and peroxyl
radical generation.[35] The ROS generation
in PC3 and BPH1 cell lines after treating with the Ccomplex for 24 h is shown in Figure . It was observed that PC3cells exposed to the Ccomplex
produced a significantly high amount (p < 0.005)
of ROS, compared to BPH1 cells. The intracellular imaging behavior
of the Ccomplex was studied
in both PC3 and BPH1 cell lines using fluorescence microscopy, and
the results obtained are illustrated in Figure . After incubation with the Ccomplex, BPH1 cells display no intracellular
fluorescence. However, PC3 shows green fluorescence both in the cytoplasm
and nuclei, suggesting that the Ccomplex was distributed both in the cytosol and nucleus in the proliferating
cancercell line.
Figure 9
Intracellular ROS generation of PC3 and BPH1 cell lines
treated
with the C complex for 24 h
at two different doses (IC50 doses
of PC3). Data are presented as % increase in ROS relative to untreated
controls. Values are the mean ± SD of three independent experiments.
***(p < 0.0005) and ** (p <
005) denote the statistically significant difference compared to untreated
controls.
Figure 10
Cellular imaging of the C complex in both PC3 and BPH1 cell lines
under a fluorescence microscope.
DAPI and the C complex are visible
as blue and green fluorescence, respectively. BPH1 cells display no
intracellular fluorescence but PC3 shows green fluorescence both in
the cytoplasm and nuclei.
IntracellularROS generation of PC3 and BPH1 cell lines
treated
with the Ccomplex for 24 h
at two different doses (C50 and >IC50 doses
of PC3). Data are presented as % increase in ROS relative to untreated
controls. Valuesare the mean ± SD of three independent experiments.
***(p < 0.0005) and ** (p <
005) denote the statistically significant difference compared to untreated
controls.
Cellular imaging of the Ccomplex in both PC3 and BPH1 cell lines
under a fluorescence microscope.
DAPI and the Ccomplex are visible
as blue and green fluorescence, respectively. BPH1 cells display no
intracellular fluorescence but PC3 shows green fluorescence both in
the cytoplasm and nuclei.
Conclusions
The present study shows that
the guanidine-based ligands (L) undergo
facile reaction with [{Ru(p-cymene)Cl2}2] to furnish cationic half-sandwich complexes of type
[Ru(p-cymene)(L)Cl]+. This study also
reveals that the complexesare better DNA binders than the corresponding
uncoordinated guanidine-based ligands, and the observed enhancement
in DNA binding is attributable to the imposed planarity of the guanidine-based
ligand upon coordination to ruthenium that enabled it to serve as
a better intercalator. Cytotoxicity studies also show a similar trend,
the complexes being more cytotoxic than the uncoordinated guanidine-based
ligands, presumably because complex formation leads to an improvement
in cellular uptake that permits more molecules to enter cells, showing
greater cytotoxicity. The other important aspect is that compared
to cancercells, the complexes were found to be significantly less
toxic to normal cells, and this is most prominent in the Ccomplex. This is probably due to the increased
uptake of the complex molecules in cancercells than normal cells,
as the membrane transport system of cancercells is more active than
that of the normal cells or benign cells. Additionally, more uptake
of complex molecules generates more reactive oxygen species that lead
to more oxidative DNA damage as observed by the comet assay.[36] This study also demonstrates that inclusion
of the guanidine-based ligands in the half-sandwich ruthenium–arenecomplexes, particularly in the Ccomplex, has been useful for exhibition of remarkable antiproliferative
activity against cancercells with high selectivity and also for convenient
tracing of the complexes in cells due to their prominent emissive
nature. It is worth highlighting that such studies involving ligand
modification at a single point (NH vs O vs S) are rare in the literature.[37]
Authors: Fengrui Qu; Seungjo Park; Kristina Martinez; Jessica L Gray; Fathima Shazna Thowfeik; John A Lundeen; Ashley E Kuhn; David J Charboneau; Deidra L Gerlach; Molly M Lockart; James A Law; Katherine L Jernigan; Nicole Chambers; Matthias Zeller; Nicholas A Piro; W Scott Kassel; Russell H Schmehl; Jared J Paul; Edward J Merino; Yonghyun Kim; Elizabeth T Papish Journal: Inorg Chem Date: 2017-06-21 Impact factor: 5.165
Authors: Laura A Hager; Stephan Mokesch; Claudia Kieler; Silvia Alonso-de Castro; Dina Baier; Alexander Roller; Wolfgang Kandioller; Bernhard K Keppler; Walter Berger; Luca Salassa; Alessio Terenzi Journal: Dalton Trans Date: 2019-07-11 Impact factor: 4.390