Xin
Zhang‡
a,
Fortuna
Ponte‡
b,
Elisa
Borfecchia
c,
Andrea
Martini
c,
Carlos
Sanchez-Cano§
*a,
Emilia
Sicilia
*b and
Peter J.
Sadler
*a
aDepartment of Chemistry, University of Warwick, Coventry, CV4 7AL, UK. E-mail: p.j.sadler@warwick.ac.uk
bDipartimento di Chimica e Tecnologie Chimiche, Università della Calabria, Ponte P. Bucci, Cubo 14cI-87030, Arcavacata di Rende, Italy. E-mail: emilia.sicilia@unical.it
cDepartment of Chemistry, NIS Center and INSTM Reference Center, University of Turin, via Giuria 7, 10125, Turin, Italy
First published on 19th November 2019
In contrast to the clinical drug cisplatin, the anticancer complex [Os(η6-p-cymene)(4-(2-pyridylazo)-N,N-dimethylaniline)I]+ [1-I] is inert towards hydrolysis and targets cancer cell metabolism rather than DNA. A combination of DFT calculations and X-ray absorption spectroscopy (XAS) suggests that hydrolytic activation of 1-I involves catalytic attack by the intracellular tripeptide glutathione (GSH) on the azo bond of the chelating ligand in the complex.
Prodrugs hypothesised to be activated by intracellular reduction include PtIV complexes which are reduced to active PtII species.3 A number of RuIII anticancer complexes are also thought to be prodrugs activated in a similar way,4 but recent experiments have not detected reduction of RuIII compounds in tumours.5 However, activation of 1-I and its analogues does not appear to involve changes in the oxidation state of OsII in cells. On the contrary, initial experiments suggest that hydrolysis of the Os–I bond is promoted by the tripeptide glutathione (GSH). This generates the more reactive hydroxido adduct 1-OH, which can react with other molecules in cells, such as H2O2 (generating hydroxyl radicals).1 The aim of the present work was to establish both experimental and theoretical (density functional theory, DFT) evidence for this proposed highly unusual mechanism of activation of a metallo-prodrug.
DFT computational analysis can provide insights into reaction mechanisms for metal coordination compounds.6 Hence we probed the most favourable mechanism, involving deprotonated glutathione (GS−), for the hydrolysis of 1-I using quantum-mechanical DFT calculations. First, the direct replacement of the iodido ligand by water (pathway a), a process observed under physiological conditions for many OsII, RuII, or IrIII half-sandwich complexes [M(arene)(L)X]+,7 was explored. Our calculations suggested an overall energy barrier of 25.4 kcal mol−1 and an endergonic (by 9.8 kcal mol−1) exchange by second-order nucleophilic substitution (SN2) (pathway a; Fig. S1a, ESI†). An alternative, direct replacement of the iodide ligand by GS−, to form an Os-SG adducts, followed by replacement of the thiolate ligand by water was also studied (pathway b). Such a reaction has been reported for other complexes, e.g. [Ru(arene)(en)Cl]+ and its analogues.8 In this case, the energy barrier to generate the stable 1-SG complex was 23.2 kcal mol−1, and the substitution product lies 18.9 kcal mol−1 below the entrance channel. However, hydrolysis of 1-SG to form the experimentally observed 1-OH required surmounting a very high energy barrier of 30.3 kcal mol−1 (pathway b, Fig. S1b, ESI†). Finally, we explored ligand-based attack on the complex by GSH resulting in hydrolysis of the Os–I bond through redox mediation involving the NN double bond (pathway c).
Some organic azo compounds are known to catalyse oxidation of GSH in presence of O2,9 reactions also possible for azo-containing ligands in RuII or IrIII complexes.10 The reported mechanism can involve addition of GSH to the diazene double bond and subsequent reduction to an hydrazo group.9 However, the free azo ligand in 1-I does not react with GSH.11
We probed azo bond attack by GSH using DFT, but the most favourable energy barrier for the only viable reaction mechanism was high (40.2 kcal mol−1 Fig. S2, ESI†). Nevertheless, the attack by GS− on the azo N atom bound to the pyridine ring, preluding hydrolysis of 1-I, yielded the lowest barrier (19.0 kcal mol−1). This led to intermediate 1-GS/NN (pathway c, Fig. 1), that lies 11.8 kcal mol−1 above the reactant reference energy, and accumulated electronic density near Os, causing the weakening of the Os–I bond. This facilitates displacement of the I− ligand by water (13.9 kcal mol−1), which occurred simultaneously with detachment of the GS− and regeneration of the azo bond. Further deprotonation of the coordinated water (requiring 5.5 kcal mol−1) yielded 1-OH, which was stabilised by 7.3 kcal mol−1 (with respect to entrance channel). Computational analysis of the activation of 2-I, an analogue with an electron donating ethoxy group on the pyridine ring gave similar results, but with lower energy barriers (Fig. S3, ESI†). This paralleled experimental data showing faster hydrolysis of 2-I compared to 1-I.1
Fig. 1 DFT free energy profile describing the most favourable pathway for the reaction of 1-I with GS− (pathway c). Calculated energies are in kcal mol−1 and relative to reference energy of separated reactants. Alternative higher energy pathways a and b are in Fig. S1 (ESI†). |
Interestingly, our calculations proposed thermodynamically favoured hydrolysis pathways in the presence of glutathione, and kinetic control favouring the azo-based mechanism. This was confirmed by following the hydrolysis of chlorido [Os(η6-p-cymene)(4-(2-pyridylazo)-N,N-dimethylaniline)Cl]+ [1-Cl] and imino [Os(η6-p-cymene)(4-(2-iminopyridine)-N,N-dimethyl-aniline)I]+ [1-impy-I] analogues in presence or absence of GSH using HPLC (the solvent gradient used is on page S11 of the ESI†). As for 1-I, 1-Cl is stable in aqueous solution, but treatment with GSH promotes faster hydrolysis than 1-I (Fig. S18, ESI†). However, both 1-Cl and 1-I hydrolysed in the presence of GSH following the same mechanism, generating first 1-OH (Fig. S19, S20 and Tables S3, S4, ESI†), which can be attacked by GSH or Cl− to form 1-Cl or 1-SG. 1-SG can be further oxidised to the sulfenate complex 1-SOG (Fig. S21 and S22, ESI†). Notably, similar pathways were observed when different thiols induce hydrolysis of 1-I (Fig. S23, S24 and Table S5, ESI†). In contrast, the iminopyridine complex 1-impy-I readily hydrolysed in the absence of GSH, generating both hydroxido 1-impy-OH and chlorido 1-impy-Cl species. Incubation with GSH did not increase the extent of hydrolysis, but generated the thiolato adduct 1-impy-SG (Fig. S18 and S25, ESI†), even when 1-impy-I was in large excess. This confirmed that the azo bond in the azopyridine complexes is the likely target site for GS− (as it is the only difference between the chemical structure of 1-I and 1-impy-I). Hydrolysis of 1-I was ca. 3× faster at pH 10 compared to pH 7, attributable to the higher proportion of GS− present (pKa(SH): 9.42),12 whereas hydrolysis was inhibited at pH 2.4 (Fig. S26 and Table S6, ESI†). No 1-OH was observed in the absence of thiol at basic pH (Fig. S27, ESI†).
The proposed attack of GSH on the –NN– bond in 1-I during the hydrolysis was further probed using X-ray absorption spectroscopy (XAS), to provide structural information.13 We acquired XAS spectra at the Os L3-edge for solid reference compounds 1-I, 1-Cl, 1-OH and 1-SG (for synthesis and characterisation see ESI†), and an aqueous solution (pH 7 buffer) of 1-I on beamline BM23 at ESRF.14 Discrimination between the different references could not rely on analysis of the XANES region (Fig. 2a and Fig. S28, ESI†).15 However, there were distinguishing differences in the EXAFS oscillations for 1-I, when compared with 1-Cl, 1-OH and 1-SG especially from k = 10 Å−1 upward (Fig. 2b–d and Fig. S29, S30, ESI†). In contrast, the spectra of 1-Cl and 1-SG were similar (due to similar size and atomic number of Cl and S). XAS spectra from the reference compounds were also used to validate structural models generated by our DFT calculations (Fig. S31, ESI†). Based on these data, we simulated an EXAFS signature corresponding to the DFT model of the intermediate 1-GS/NN (Fig. S31, see ESI† for details). This allowed us to probe the presence of 1-GS/NN in reactions between 1-I and GSH followed in situ by EXAFS.
The kinetics for the hydrolysis of 1-I in presence of GSH was further monitored in situ by acquiring consecutive EXAFS spectra (Fig. 3). Data analysis used the structural models obtained from EXAFS spectra of the experimental standards and DFT calculations for the 1-GS/NN intermediate. Initially, an in situ XAS experiment was carried at pH 7, at an Os:GSH mol ratio of 1:2 (Fig. 3a). Under these conditions (30 mM 1-I, 60 mM GSH; pH 7) Os–I hydrolysis was much faster than in HPLC experiments, and the reaction was complete after ca. 40 min, with only minor further spectral modifications up to 220 min. Linear combination fits (LCF) indicated the dominant presence of 1-SG and 1-Cl after 80 min, accounting for >80% of the Os species in the reaction mixture, together with contributions of <10% from 1-I and 1-OH (Table S9 and Fig. S33, ESI†). LCF provided estimates of the combined % of chlorido and thiolato adducts (Table S9, ESI†) since 1-Cl and 1-SG could not be readily discriminated by EXAFS analysis. In contrast, when the same experiment (30 mM 1-I, 60 mM GSH) was repeated at pH 3, EXAFS spectra were almost identical to 1-I standards, and remained unchanged over the duration of study (Fig. 3b). This suggested that it is possible to capture the dynamics of some steps in these reactions from XAS measurements by careful tuning of the reaction conditions.
The reaction of 1-I (30 mM) with GSH was further studied at pH = 5, at an Os:GSH ratio of 2:1, in an attempt to moderate the reaction rate (as the [GS−] at pH 5 is two orders of magnitude lower than at pH 7, according to the Henderson–Hasselbalch equation;16Fig. 3c and Fig. S34, ESI†). Remarkably, EXAFS spectra in the 40–220 min reaction time range, showed an oscillatory behaviour. Periodic and opposite changes in the intensities of the peaks at ca. 2.5 Å (Os–I single scattering [SS] contribution from 1-I), and the shoulder at ca. 2.0 Å (Os–S SS contribution from 1-SG) in the phase-uncorrected spectra were observed (Fig. 3c). Yet, certain time points showed dampened contributions from Os–I SS without an increase in the characteristic Os–S SS (e.g. after 200 min; red curve in Fig. 3c), suggesting the presence within the reaction of additional unknown Os species. Furthermore, the cyclic behaviour might arise from alternating periods dominated by cleavage or formation of Os–I bonds. To check this, Os-speciation as a function of the reaction time was evaluated by applying EXAFS LCF analysis in k-space (Fig. 3d). Based on our observations, three different components were used to generate the LCF model: the reference compounds 1-I and 1-SG and the simulated EXAFS signature for the 1-GS/NN intermediate predicted by theory. The model adopted reproduced the experimental data well (Fig. S35, ESI†). Inclusion in the LCF model of the 1-GS/NN component resulted in a significant improvement of the fitting quality for the time points with high concentrations of such species (coloured arrows in Fig. 3d) with respect to a two-component model of only 1-I and 1-SG (Fig. S36, ESI†). Crucially, the maxima in the relative abundance of 1-GS/NN (16–24% total Os) systematically preluded increases in the concentration of 1-SG. These data supported the presence of the 1-GS/NN species during the reaction, in line with its intermediate character defined by DFT. LCF results highlighted the oscillating character of the reaction. 1-I was cyclically converted into 1-SG (with a period of 60–80 min), which then was transformed again to 1-I. This may be related to fluctuations in the concentrations of free I− ions and free GSH as a result of the hydrolysis of Os–I bonds, and the formation of new Os–SG species. Indeed, our chromatographic experiments have shown previously that the presence of high concentrations of Cl− ions in the medium can alter the amount of 1-Cl detected during the hydrolysis of 1-I in the presence of GSH. Moreover, a more complex four-component LCF model (with 1-I, 1-SG and 1-OH and simulated 1-GS/NN references) suggested negligible presence of 1-OH in the reaction mixture (Fig. S37, ESI†). This might indicate the presence of high concentrations of free I− ions and free GSH. By calculating the concentrations of free I− ions and free GSH (from amounts of 1-I and 1-SG species indicated by LCF), we observed an alternation between long periods of low I−/GSH ratios and short periods of high I−/GSH ratios (every 60–80 min). Remarkably, the highest concentrations of free I− ions were always just before periods dominated by the formation of Os–I bonds (Fig. S38, ESI†), providing a possible explanation for the regeneration of 1-I. The presence of fluctuations in GSH, together with complete hydrolysis of Os–I bonds (31.2 mM in total, calculated from amounts of 1-I and 1-SG indicated by LCF) by 15 mM of GSH suggested a catalytic reaction, with the thiol being regenerated. This was confirmed by the almost total conversion of 1-Cl to 1-OH after 12 h incubation with sub-stoichiometric quantities of GSH (0.5, 0.25 or 0.125 mol equiv.; Fig. S39, ESI†). Furthermore, significant amounts of 1-OH (525 μM) were generated with little consumption of GSH when 1-I (2 mM) reacted with 0.1 mol equiv. of the thiol (200 μM) for 2 h (Fig. 4) under N2 (protecting GSH from oxidation). 1H NMR spectra provided a turnover frequency (TOF) number of 1.97 h−1 for the reaction under such conditions.
Fig. 4 Time-dependent concentration changes for GSH, 1-I and 1-OH in the reaction of 1-I (2 mM) with GSH (200 μM), 75 mM phosphate buffer (pH 7.4), 298 K, for 2 h based on integration of 1H NMR peak areas (Table S10 and Fig. S40, ESI†). |
In summary, DFT calculations suggested that substitution of bound iodide by water in the inert organo-osmium anticancer complex 1-I is dependent on high energy kinetic barriers, but revealed a favourable reaction pathway involving attack of GS− on the –NN– azo bond in the chelating ligand, giving a 1-GS/NN intermediate (Fig. S31, ESI†) containing an N–SG bond.
Experimental chromatographic analysis of the GSH-induced hydrolysis then confirmed the importance of the azo ligand in the mechanism of the hydrolysis, and demonstrated enhanced rates at high pH and inhibition at low pH (importance of GS− in the attack). Time resolved XAS data provided evidence for the proposed 1-GS/NN intermediate during the hydrolysis reaction. Our experiments also revealed the possible catalytic nature of the reaction (TOF ca. 2 h−1). Overall, these investigations reveal a novel activation pathway for a metallodrug in which an inert metal–ligand bond is activated towards hydrolysis, and subsequent interactions with biomolecules, by the attack of the abundant intracellular thiol glutathione on a chelating ligand. Such findings have implications for the design of new generations of anticancer metallodrugs.
XAS spectra were acquired in the frame of the ESRF proposal CH-5212. We thank M. Monte Caballero and E. M. Bolitho for the support during experiments at BM23. We thank CRUK/EPSRC (Grant No. C53561/A19933), EPSRC (Grant No. EP/F034210/1) and Wellcome Trust (Grant No. 107691/Z/15/Z) for financial support.
Footnotes |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/c9cc06725f |
‡ These authors contributed equally to the manuscript. |
§ Current address: CIC biomaGUNE, Parque Científico y Tecnológico de Gipuzkoa, Paseo Miramón 182, 20014 San Sebastian, Spain. E-mail: csanchez@cicbiomagune.es |
This journal is © The Royal Society of Chemistry 2019 |