Effect of sulfonamidoethylenediamine substituents in Ru arene anticancer catalysts on transfer hydrogenation of coenzyme NAD by formate†

A series of neutral pseudo-octahedral Ru sulfonamidoethylenediamine complexes [(η-p-cym)Ru(N,N’) Cl] where N,N’ is N-(2-(R,R-amino)ethyl)-4-toluenesulfonamide (TsEn(R,R)) R,R = Me,H (1); Me,Me (2); Et,H (3); benzyl,H (Bz, 4); 4-fluorobenzyl,H (4-F-Bz, 5) or naphthalen-2-ylmethyl,H (Naph, 6), were synthesised and characterised including the X-ray crystal structure of 3. These complexes catalyse the reduction of NAD regioselectively to 1,4-NADH by using formate as the hydride source. The catalytic efficiency depends markedly on the steric and electronic effects of the N-substitutent, with turnover frequencies (TOFs) increasing in the order: 1 < 2 < 3, 6 < 4, 5, achieving a TOF of 7.7 h for 4 with a 95% yield of 1,4-NADH. The reduction rate was highest between pH* (deuterated solvent) 6 and 7.5 and improved with an increase in formate concentration (TOF of 18.8 h, 140 mM formate). The calculations suggested initial substitution of an aqua ligand by formate, followed by hydride transfer to Ru and then to NAD, and indicated specific interactions between the aqua complex and both NAD and NADH, the former allowing a preorganisation involving interaction between the aqua ligand, formate anion and the pyridine ring of NAD. The complexes exhibited antiproliferative activity towards A2780 human ovarian cancer cells with IC50 values ranging from 1 to 31 μM, the most potent complex, [(η-p-cym)Ru(TsEn(Bz, H))Cl] (4, IC50 = 1.0 ± 0.1 μM), having a potency similar to the anticancer drug cisplatin. Co-administration with sodium formate (2 mM), increased the potency of all complexes towards A2780 cells by 20–36%, with the greatest effect seen for complex 6.


Introduction
Nicotinamide adenine dinucleotide (NAD + ) and its reduced form (NADH), as well as their phosphorylated derivatives, NADP + and NADPH, play a vital role in biological systems as redox coenzymes. 13][4] The reduction of pyridinium salts (e.g.NAD + ) to dihydropyridine com-pounds (e.g.NADH) is of critical importance for energy storage and release in cell metabolism. 3,4][11][12][13] Compared to reduction with H 2 (hydrogenation), transfer hydrogenation (TH) reactions have the advantage of being simpler, without the need for any high external pressure and use readily available, safer-to-handle, hydride sources. 14Also, TH reduction of NAD(P) + artificially has attracted wide interest as an in vitro mimic for enzymatic reactions performed under biologically relevant conditions. 15he pathways of hydride transfer between pyridinium salts and dihydropyridine compounds are also of interest.16,17 They proposed a catalytic cycle involving a ring-slippage η 4 -Cp* intermediate with Rh coordinated to the amide of the pyridine ring. 18Knör et al. reported a Rhcoordinated poly(arylene-ethynylene)-alt-poly(arylene-vinylene) polymer as photocatalyst for the reduction of NAD + ; involving a possible photoexcited polymer chain being quenched and transferring an electron to the Rh III active centre. 19More recently, Yoon et al. described a mechanism involving hydride transfer to Cp* and formation of the Rh I intermediate [(η 4 -Cp*-H)Rh((CH 2 OH) 2 -bipy)] + followed by hydride transfer from the endo orientation of the C-H bond to maintain the 1,4-regioselectivity of NADH. 20he half-sandwich ruthenium complex [(η 6 -p-cym)Ru (TsDPEN)Cl] (TsDPEN: N-((1S,2S)-2-amino-1,2-diphenylethyl)-4-methylbenzenesulfonamide) was first reported by Noyori and coworkers in 1995. 21,22Potent catalytic activity has been shown in asymmetric TH reduction of aromatic ketones.Most recently, the 16-electron Os analogues [(η 6 -arene)Os(TsDPEN)] of Noyori type complexes were reported to reduce pyruvate enantioselectively to (Dor L-) lactate via asymmetric transfer hydrogenation in human cancer cells. 23Nonetheless, the hydrophobic nature of the two phenyl groups on the ethylene backbone limits its application as a possible catalyst for TH reduction of NAD + under biologically relevant conditions.Complexes with chelating diamine ligands such as complex 7 in Fig. 1, display good aqueous solubility but poor catalytic activity in TH reduction of NAD + . 24However, p-cymene ( p-cym) complexes with functional sulfonyl substituents such as [(η 6 p-cym)Ru(TsEn)Cl] (e.g.complex 8 in Fig. 1), 25 exhibit good solubility in water and improved catalysis for NAD + reduction to NADH in aqueous media.Moreover, co-administration of [(η 6 -p-cym)Ru(TsEn)Cl] with low non-cytotoxic doses of sodium formate led to an enhancement of the antiproliferative activity against A2780 human ovarian cancer cells by up to 50×. 15,26ere we investigate the effect on catalytic reduction of NAD + using formate as a hydride source upon variation of substituents on the amino group of the N,N-chelating TsEn ligand in Ru II 1).In addition, the catalytic mechanism was investigated both experimentally and by density functional theory (DFT) calculations.We also explored the effect of co-administration of formate on the antiproliferative activity of these complexes against A2780 human ovarian cancer cells.
A2780 human ovarian carcinoma cells were obtained from the European Collection of Cell Cultures.The cell line was grown in Roswell Park Memorial Institute medium (RPMI-1640) supplemented with 10% of foetal calf serum, 1% v/v of 2 mM glutamine and 1% v/v penicillin/streptomycin (10 000 units).All cells were grown as adherent monolayers at 310 K in a 5% CO 2 -humidified atmosphere and passaged at ca. 70-80% confluency.

Instruments
NMR spectra were acquired on Bruker HD-300, HD-400, HD-500, and AV III 600 spectrometers. 1 H NMR chemical shifts were internally referenced to TMS via 1,4-dioxane in D 2 O (δ = 3.75 ppm) or residual protiated d 4 -MeOD (δ = 3.31 ppm), or CDCl 3 (δ = 7.26 ppm).1D spectra were recorded using standard pulse sequences.Typically, data were acquired with 16 transients into 32k data points over a spectral width of 14 ppm and, for the kinetic experiment, 32 transients into 32k data points over a spectral width of 30 ppm using a relaxation delay of 2 s.Elemental analysis were performed by Warwick Analytical using an Exeter Analytical elemental analyzer (CE440).
Positive ion electrospray mass spectra were obtained on an Agilent 6130B ion mass spectrometer.High resolution mass spectrometry data were obtained on a Bruker Maxis Plus Q-TOF instrument.
X-ray crystallographic diffraction data were collected on an Oxford Diffraction Gemini four-circle system with a Ruby CCD area detector.The structure was refined by full-matrix leastsquares against F 2 using SHELXL 9735 and solved by direct methods using SHELXS36 (TREF) with additional light atoms found by Fourier methods.The atoms from the sulfonamide nitrogen to the end of the chain (C10 C11 N12 C13) were modelled as disordered over two positions related by a small ruffle in the chain.The occupancy of the two positions was linked to a free variable which refined to 86 : 14.The minor component was refined isotropically.X-ray crystallographic data for complex 3 has been deposited in the Cambridge Crystallographic Data Center (CCDC) under the accession number CCDC 1571331.† ICP-OES analysis were carried out on a PerkinElmer Optima 5300 DV series Optical Emission Spectrophotometer.The water used for ICP-OES analysis was doubly deionized (DDW) using a Millipore Milli-Q water purification system and a USF Elga UHQ water deionizer.The ruthenium Specupure plasma standard (ruthenium chloride, 1004 ± 5 µg mL −1 in 10% v/v hydrochloric acid) was diluted with 3.6% v/v HNO 3 to freshly prepare calibrants at concentrations of 50-700 ppb.Calibration standards were adjusted to match the sample matrix by standard addition of sodium chloride (TraceSELECT®).Total dissolved solids did not exceed 0.2% w/v.Data were acquired and processed using WinLab32 V3.4.1 for Windows.
ICP-MS analysis were carried out on an Agilent Technologies 7500 series ICP-MS instrument.The water used for ICP-MS analysis was double-deionized (DDW) using a Millipore Milli-Q water purification system and a USF Elga UHQ water deionizer.The Ruthenium Specpure plasma standard (ruthenium chloride, 1004 ± 5 µg mL −1 in 10% v/v hydrochloric acid) was diluted with 3.6% v/v HNO 3 to prepare calibrants freshly at concentrations of 0.1-1000 ppb.The ICP-MS instrument was set to detect 101 Ru in no gas mode.Total dissolved solids did not exceed 0.1% w/v.An internal standard of 166 Er (50 ppb) was used.Data were acquired using ICP-MS-TOP and proceeded using Offline Data Analysis (ChemStation version B.03.05, Agilent Technologies, Inc.).pH values were measured using a Minilab IQ125 pH meter equipped with a ISFET silicon chip pH sensor and referenced in KCl gel.pH* values ( pH meter reading without correction for the effect of deuterium on the sensor) of NMR samples in D 2 O were measured at 310 K. Relative hydrophobicity measurements were performed utilising the Agilent 1200 HPLC system with a VWD and 50 µL loop.The column was an Agilent Zorbax 300SB C18, 150 × 4.6 mm with a 5 µm pore size.The mobile phase was H 2 O (50 mM NaCl)/H 2 O/CH 3 CN 1 : 1 (50 mM NaCl), with a flow of 1 mL min −1 .The detection wavelength was set at 254 nm with the reference wavelength at 360 nm.

Turnover frequency determination
UV-vis spectroscopy.In a typical experiment, 330 μL of each solution (84 µM complex in MeOH/H 2 O 1 : 9 v/v, 102 mM sodium formate and 510 µM NAD + in H 2 O) was added to a 1 mL cuvette, and the pH adjusted to 7.2, bringing the total volume to 1 mL (final concentrations: Ru complex 28 µM; NAD + 170 µM; NaHCO 2 34 mM; molar ratio 1 : 6 : 1200).UV spectra were recorded and the absorbance at 340 nm was monitored every 5 min until completion of the reaction.
NMR spectroscopy.Complexes were dissolved in d 4 -MeOD/ D 2 O (1 : 4, v/v) (1.4 mM) in a glass vial.Solutions of sodium formate (35 mM) and NAD + (5.6 mM) in D 2 O were also prepared and then incubated at 310 K, pH* 7.2 ± 0.1.An aliquot of 200 μL from each solution was added to a 5 mm NMR tube, giving a final volume of 0.64 mL (Ru complex 0.44 mM; NAD + 1.75 mM; NaHCO 2 10.94 mM; molar ratio 1 : 4 : 25).A 1 H NMR spectrum was recorded at 310 K every 162 s until the completion of the reaction.Further experiments under similar conditions using different concentrations of sodium formate (complex 4, NAD + and sodium formate in ratio of 1: 4: X, where X = 10, 25, 50, and 100 mol equiv.)and different concentrations of NAD + (complex 4, NAD + and sodium formate in ratio of 1 : Y : 25, where Y = 2, 4, 6 and 10) were also studied.Another series of experiments using different pH* values of the reaction solutions (5, 6, 7, 8 and 9) were also performed.Molar ratios of NAD + and NADH were determined by integrating 1 H NMR peaks corresponding to NAD + (9.33 ppm) and 1,4-NADH (6.96 ppm).The turnover number (TON) for the reaction was calculated as follows: where I n is the integral of the signal at n ppm and [NAD + ] is the concentration of NAD + at the start of the reaction.

Cell growth inhibition assays
The antiproliferative activity of complexes 1-6 was determined in A2780 ovarian cancer cells.Briefly, 96-well plates were used to seed 5000 cells per well.Cells were incubated in drug-free medium at 310 K for 48 h before addition of tested compounds ( prepared by serial dilution in culture medium containing 5% DMSO, typically 6 concentrations in the range: 0.01-100 μM).Exact Ru concentrations were determined by ICP-OES and the maximum concentration of DMSO to which cells were exposed never exceeded 0.5% v/v.A drug exposure period of 24 h was allowed.After this, supernatants were removed by suction and each well was washed with PBS.A further 72 h were allowed for the cells to recover in drug-free medium at 310 K.The sulforhodamine B (SRB) assay was used to determine cell viability. 30IC 50 values, as the concentration that causes 50% cell death, were determined as duplicates of triplicates in two independent sets of experiments and their standard deviation were calculated.Data were processed using Microsoft Excel and sigmoidal curves fitted using Origin 9.1.

Co-administration of Ru complexes with formate
Cell viability assays were carried out with complexes 1-6 with co-administration of sodium formate in A2780 ovarian cancer cells.These experiments were carried out as above (in vitro growth inhabitation assay) with the following modifications: a fixed equipotent concentration of each Ru complex equal to 1/3 × IC 50 in that cell line was used in coadministration with three different concentrations of sodium formate (0.5, 1.0 and 2.0 mM).Drug stock solutions (ca. 100 µM) were prepared and they were further diluted using media until working concentrations were achieved.Separately, a stock solution of sodium formate was prepared in saline.The complex and formate solutions were added to each well independently, but within 5 minutes of each other.All other experiment conditions were kept unchanged (drug exposure and cell recovery time, as well as, end point assay used).

Cellular accumulation
The accumulation studies for Ru complexes 1-6 were performed on A2780 ovarian cancer cells.1.5 × 10 6 cells were seeded on a six-well plate using 2 mL of cell culture medium.After 24 h of pre-incubation in drug-free medium at 310 K, cells were exposed to complexes at equipotent IC 50 concentrations for 24 h ( prepared by serial dilution of a ca. 100 μM stock solution, prepared using culture medium containing 5% DMSO.This solution was analysed by ICP-OES to determine Ru concentration before treatment of cells with Ru complex).After this time, drug solutions were removed by suction, cells were washed with PBS and then treated with trypsin-EDTA.A suspension of single cells was counted, and cell pellets were collected.Each pellet was digested overnight in freshly-distilled concentrated nitric acid (200 μL, 72% v/v) at 353 K; the resulting solutions were diluted with double-distilled water to a final concentration of 3.6% v/v HNO 3 , and the amount of Ru in A2780 ovarian cells was determined by ICP-MS.These experiments did not include any cell recovery time in drug-free media; they were carried out in triplicate, and the standard deviations were calculated.Data were processed using Microsoft Excel and reported as ng Ru × 10 6 cells.

ROS determination
Flow cytometry analysis of ROS/superoxide induction in A2780 cells caused by exposure to complexes 1 and 4 was carried out using the Total ROS/Superoxide detection kit (Enzo-Life Sciences) according to the manufacturer's instructions.Briefly, 1.0 × 10 6 A2780 cells per well were seeded in a six-well plate.Cells were preincubated in drug-free media at 310 K for 24 h in a 5% CO 2 humidified atmosphere, and then drugs were added to triplicates wells at IC 50 concentration.After 24 h of drug exposure, supernatants were removed by suction and cells were washed with PBS and harvested.Staining was achieved by resuspending the cell pellets in buffer containing the orange/ green fluorescent reagents.Cells were analysed in a Becton Dickinson FACScan flow cytometer using FL1 channel Ex/Em: 490/525 nm for the oxidative stress and FL2 channel Ex/Em: 550/620 nm for superoxide detection.Data were gated using positive-stained ( pyocyanin positive control), untreated-stained and untreated-unstained control samples, acquired as instrumental triplicates, and were processed using FlowJo V10 for Windows.At all times, samples were kept under dark conditions to avoid light-induced ROS production.

DFT calculations
The DFT calculations of electronic energy levels of the catalytic cycle were based on the crystal structure of complex 3.The method of the calculation was functional CAM-B3LYP 31 with basis set CEP-31G, [32][33][34] using Gaussian 16 software. 35Ultrafine grid of integration was used in each case.The starting geometry was taken from X-ray data for 3, with an appropriate change of substituents for other systems.All given energy values are the result of the full geometry optimisation with subsequent frequency calculations.Optimisations were performed with modelling of water as solvent, within the continuous polarisation model with integral equation formalism variant (IEFPCM keyword of Gaussian).Grimme empirical corrections for dispersion were applied (keyword GD3).The optimisations were performed using the solvent-accessible surface option and the final energy was calculated with using solvent-excluding surface options (keywords surface = sas and ses, respectively).NAD + was modelled with an effective charge of −1, with two deprotonated phosphate groups; the same protonation of phosphate was used for NADH, giving an effective charge of −2.

DNA binding
The reactions of complex 4 (ca. 2 mM) with nucleobases (9-EtG and 5′-AMP) were studied typically by addition of an aqueous solution of nucleobase (3 mM, 1.5 mol equiv.) in 10% of d 4 -MeOD and 90% of D 2 O, pH* 7.2 ± 0.1, and monitored by 1 H NMR at 310 K. Solutions of double-helical calf thymus DNA (ct-DNA) at a concentration of 32 µg mL −1 were incubated with complex 4 at r i value of 0.1 in 10 mM NaClO 4 at 310 K (r i is defined as the molar ratio of free ruthenium complex to nucleotide phosphates at the onset of incubation with DNA).The concentration of ruthenium associated with DNA in these samples was determined by flameless atomic absorption spectrometry (FAAS).The concentrations of DNA were determined by absorption spectrophotometry.Plasmid DNA pBR322 (28 µg mL −1 ) and complex 4 in various molar ratios (r i = 0.05-1) were incubated in 0.01 M NaClO 4 at 310 K for 24 h in the dark.
Then the samples were directly mixed with the loading buffer and loaded onto a 1% agarose gel running at 298 K in the dark with Tris-acetate-EDTA (TAE) buffer and the voltage set at 25 V.No separation step was included before loading the samples into the gel to allow detection of potential noncovalent binding (if any).The gels were then stained with EtBr, followed by photography with a transilluminator.

Synthesis and characterisation
Ru II complexes 1-6 were synthesised using a similar procedure to that reported for related complexes (Scheme 1). 25 Typically, triethylamine (4 mol equiv.)and ligands (ca. 2 mol equiv.)were added to a solution of [(η 6 -p-cym)RuCl 2 ] 2 in degassed isopropanol, and the reaction was stirred under a N 2 atmosphere at 365 K for 12 h.All synthesised complexes were characterised by 1 H and 13 C NMR spectroscopy, mass spectrometry (ESI-MS) and elemental analysis (CHN).A crystal of complex 3 suitable for X-ray analysis was obtained by diffusion of diethyl ether into a solution of the complex in methanol at ambient temperature.Selected bond lengths and angles for complex 3 are listed in Table 2. Crystallographic data are presented in Table S1, † and the structure of complex 3 is shown in Fig. 2. Complex 3 adopts a pseudo-octahedral geometry with the η 6bonded aromatic ring occupying 3 coordination sites.The chelating ligand is deprotonated and bonded as a monoanionic bidentate ligand.The CH 2 CH 2 N-Et atoms from N,N′ chelated ligand (C10 C11 N12 C13) were modelled as disordered over two positions whose occupancy refined to 86 : 14.Compared to reported ruthenium ethylenediamine complexes (either neutral or +1 charge), 25,36,37 the Ru-N − bond length (N9, 2.126( 9)) is within the expected range of 2.11-2.14Å, 37 but the Ru-N12 length (2.1702(11) Å) is longer than the neutral analogue [(η 6 -biph)Ru(TsEn)Cl] (2.122(3) Å), 25 suggesting that the presence of N-ethyl substituent causes a slight weakening of this Ru-N bond.The remaining bond length and angles show no significant difference.

Hydrolysis and pK *
a determination The hydrolysis of complex 4 was studied by dissolving the Ru II complex in d 4 -MeOD/D 2 O (1.4 mM, 1 : 9 (v/v)).The 1 H NMR spectrum remained unchanged after 24 h and the hydrolysis was assumed to be rapid since the peaks could be assigned to the aqua Ru II species (4a) by comparison to those from the aqua species generated in a reaction with silver nitrate in D 2 O (1 mol equiv.).The pK * a (pK a value determined in deuterated solvent) of complex 4a was determined by a pH* (meter reading) titration ranging from 2 to 12 by addition of NaOD or DNO 3 solutions as appropriate.Changes in the chemical shift of a tosyl 1 H NMR resonance were followed and the data were fitted to the Henderson-Hasselbalch equation, giving a pK * a value of 9.73 ± 0.06 (Fig. 3).

Kinetics of transfer hydrogenation reactions
The ratio of coenzyme NAD + /NADH greatly influences the intracellular potential and can drive many reactions in vivo. 38he reduction of coenzyme nicotinamide adenine dinucleotide (NAD + ) to NADH was investigated in an aqueous medium using complexes 1-6 as catalysts and sodium formate as the hydride source.Initially, the TH reactions were studied by UVvisible spectroscopy under conditions of pH 7.2 ± 0.1, 310 K and MeOH/H 2 O (1 : 9, v/v, Table 3); in all the cases, an increase in intensity of the band at 340 nm was observed, which is assignable to formation of NADH (Fig. S1, ESI †).The kinetics of conversion were also monitored by 1 H NMR at 310 K and pH* 7.2 ± 0.1.The reactions were performed in a mixed solvent d 4 -MeOD/D 2 O (1 : 4, v/v), due to the poor aqueous solubility of complexes 5 and 6, although the presence of methanol in such reactions is known to enhance the reaction rate. 25heme 1 Synthetic routes for diamine ligands and Ru II complexes 1-6.In general, the introduction of substituents on the terminal nitrogen improved the catalytic activity.The bulkier the substituents on the terminal nitrogen, the higher the TH reaction rate becomes.The turnover frequency reaches a maximum (ca.7.54 h −1 ) when the substituent on the terminal N is benzyl (complex 4), making it as efficient as the Rh III complex [(η 5 -Cp*)Rh(bipy)Cl]PF 6 . 16Interestingly, the TOF decreases when the substituent is para-fluoro-benzyl (complex 5) or naphthalene (complex 6), probably, because these ligands hamper the approach of NAD + to the Ru centre.Compared to the en complex with unsubstituted nitrogens [(η 6 -biph)Ru(en)Cl]PF 6 , the turnover frequency of complex 4 is 41× higher, 24 and 2.7× higher compared to [(η 6 -p-cym)Ru(TsEn)Cl]. 25he NH proton of the chelated diamine ligand appears to be essential for the TH reduction of ketones to alcohols; 39 normally, Ru II catalysts for TH of ketones form 16-e intermediates. 40,41It has been reported that a Ru II complex with two N-alkyl groups (R,R)-[(η 6 -benzene)Ru(TsDPEN-Me 2 )Cl] exhibited poor catalytic reactivity in TH reaction of ketones. 40owever, complex 2 [(η 6 -p-cym)Ru(TsEn(Me,Me))Cl] exhibited good catalytic activity towards the TH reduction of NAD + to NADH (TOF = 4.1 h −1 , Table 3), despite not having an NH proton, which suggests, as expected, that an N-H is not essential in the transfer reduction of NAD + to NADH.
The dependence of the rate of catalysis on pH was determined.Six pH* values ranging from 5 to 9 were studied for complex 4 at a mol ratio complex 4 : NAD + : formate of 1 : 4 : 25, in the same mixed solvent at 310 K (Fig. S2, ESI †).The TOF was relatively insensitive to pH* over the range pH* 6-8 (ca.7.5 h −1 ), but decreased slightly at lower and higher pH* (5.6 h −1 at pH* 5, 6.6 h −1 at pH* 9).
The dependence of turnover frequency on the concentrations of sodium formate and NAD + was also investigated for complex 4 in d 4 -MeOD/D 2 O (1 : 4) at 310 K.The mol ratio of complex 4 : NAD + : formate was 1 : 4 : X, where X = 5, 10, 25, 50 and 100 (Fig. S3, ESI †).The TOF increased steadily from 2.2 h −1 to 18.8 h −1 as the concentration of formate was increased from 7 mM to 140 mM.Next the dependence of TOF on the NAD + concentration was studied for mol ratio complex 4 : NAD + : formate = 1 : Y : 25, where Y = 2, 6 and 10.The TOF was found to be independent of NAD + concentration (7.7 ± 0.5 h −1 ).

Antiproliferative activity and anticancer activity with formate
Ruthenium complexes have shown promise for their activity against various types of cancer cells. 42The antiproliferative activity of complexes 1-6 towards A2780 human ovarian cancer cells was determined and compared with the clinically approved drug cisplatin, Fig. 4. The IC 50 values (50% inhibition of cell growth) range from 1 to 6.5 μM for complexes containing aromatic R substituents (4-6), whereas those containing aliphatic R substituents were less potent with IC 50 values of 12-31 µM.The complex [(η 6 -p-cym)Ru(TsEn(Bz,H))Cl] (4) (IC 50 , 1.0 μM) has a potency similar to cisplatin in this cell line (CDDP, 1.20 ± 0.02 µM).It is apparent that the presence of aromatic substituents on the chelated ligands of complexes 4-6 give rise to more potent cytotoxicity than aliphatic substituents in complexes 1-3, most probably due to their higher lipophilicity.
Combination treatment with formate can greatly increase the antiproliferative activity of Ru II arene sulfonyl diamine complexes, which offers a potential new strategy for cancer  Naph,H 6.1 ± 0.9 4.9 ± 0.5 a By UV-vis spectroscopy.b By NMR spectroscopy.
treatment. 15In this work, the antiproliferative activity of Ru II complexes in A2780 human ovarian cancer cells in the presence of sodium formate was studied (Fig. 5).Firstly, the cytotoxicity of sodium formate alone towards A2780 human ovarian cancer cells was investigated.No significant toxicity was found up to formate concentrations of 2 mM which is in agreement with the previous report. 15Then, A2780 human ovarian cancer cells were coincubated with equipotent concentrations of complexes 1-6 (1/3 × IC 50 ) and three different concentrations of sodium formate (0.5, 1 and 2 mM) in order to observe the formate-concentration dependence of the cell viability.The antiproliferative activity of complexes 1-6 increased significantly upon coincubation with 2 mM formate.The formate-induced decrease in viability of A2780 cells ranged from 20% to 36% in the presence of complexes 1-6.
Interestingly for complex 6, a 28% decrease in cell viability was observed with only 0.5 mM formate present (Fig. 5, for percentage of viability decrease see Table S3, ESI †).The largest decrease of cell survival was 31% for complex 6 in the presence of 2 mM sodium formate, followed by 29% and 32% for the other two complexes with aromatic substituents, complexes 4 and 5, respectively.Complexes 1-3 with aliphatic functional groups showed an increase in potency of 18%, 21% and 22%, respectively.

Cell accumulation and relative hydrophobicity
Hydrophobicity and cellular accumulation are often important factors that play key roles in the potency of organometallic and other anticancer drugs. 43The cellular accumulation, as an equilibrium between uptake and efflux, of ruthenium in A2780 human ovarian cancer cells after exposure to complexes 1-6 at their IC 50 equipotent concentrations was determined by inductively coupled plasma mass spectrometry (ICP-MS) and is shown in Fig. 6.Complex 4 gave the lowest cellular accumulation (0.52 ± 0.08 ng of Ru per 10 6 cells), while complex 6 with moderate anticancer activity, exhibited the highest extent of cell accumulation with 4.5 ± 0.2 ng of Ru per 10 6 cells at IC 50 concentration, 8.6× higher than complex 4. Complexes 1-3 and 5, gave rise to similar cell uptake 2.4 ± 0.3 ng, 1.2 ± 0.2 ng, 3.0 ± 0.2 ng and 1.3 ± 0.2 ng per 10 6 cells, respectively, following the order: 4 < 2, 5 < 1 < 3 < 6.
The relative hydrophobicity of complexes 1-6 was determined by RP-HPLC.The more hydrophobic complexes have longer retention times on a reverse-phase C 18 column. 44To ensure solubility of the Ru II complexes in water, methanol was used as co-solvent (MeOH/H 2 O, 1 : 9 v/v) together with NaCl (50 mM) to suppress hydrolysis of the complexes.The HPLC solvents were also prepared with 50 mM NaCl (measurements see Fig. S5, ESI †).The resulting retention times are shown in Table 4, and follow the order: 1, 2, 3 < 4, 5 < 6. Complex 3 shows the shortest retention time (least hydrophobic) of   14.0 min, while complex 6 shows the longest retention time (most hydrophobic), 20.9 min.
It is evident from Table 4 that the Ru II complexes with aromatic substituents (complexes 4-6) exhibit higher hydrophobicity than complexes with aliphatic substituents (complexes 1-3).The most hydrophobic complex (6) shows the highest cell accumulation.Nonetheless, there is no linear correlation between the hydrophobicity of complexes 1-6 and their cellular accumulation.This has been observed before. 45In these cases, the chemistry and the mechanism of action of each particular complex has a higher impact on the compound's anticancer activity than cellular accumulation per se.However, complex 4 has the lowest extent of cell uptake, but the most potent antiproliferative activity, suggesting that it is the chemical properties of the intracellular drug that are more important for activity than the total amount of Ru entering the cell.In general, a high hydrophobicity could facilitate interaction between the organometallic complex and cell membranes, and also correlate with the potency of the complex, but that is not always the case. 43,45

ROS induction
Reactive oxygen species (ROS) are metabolic byproducts of aerobic respiration and are responsible for maintaining redox homeostasis in cells. 46ROS also play a significant role in the mechanism of action of anticancer agents. 47,48Some organometallic complexes, e.g.Ir and Os, [49][50][51] can generate high ROS levels or bursts of superoxide in cancer cells to induce cell apoptosis, 49 but by comparison, other complexes are known to induce cell death by reductive stress. 15The levels of reactive oxygen species (ROS) were determined in A2780 human ovarian cancer cells for complexes 1 and 4 at IC 50 concentrations by flow cytometry fluorescence analysis (Fig. 7).This included the monitoring of H 2 O 2 , peroxy and hydroxyl radicals using a green probe, and superoxide levels using the orange channel.Induction of total ROS and superoxide were determined in A2780 cells after 24 h exposure to complexes 1 and 4 when compared to the negative untreated control.The populations of cells that show high fluorescence in both FL-1 and FL-2 channels (both high total ROS and high superoxide generation) for complexes 1 and 4 are 16.5 ± 1.0% and 31.3 ± 0.3%, respectively, which indicates a higher induction of superoxide by complex 4. Remarkably, the total increase of the population in the high FL-1 green channel shows that the levels of total ROS are induced in the majority, if not in all, of the cell population.These ROS may play a major role in killing the cancer cells (Table S4, ESI †).

DNA related binding for complex 4
The interaction of complex 4 with DNA nucleobase models: 9-ethylguanine (9-EtG) and adenosine 5′-monophosphate (5′-AMP) was studied by 1 H NMR spectroscopy (Fig. S6, ESI †).The reactions were performed by adding nucleobase solution (3 mM in D 2 O) to Ru II complex solution (2 mM in 10% d 4 -MeOD/90% D 2 O) at 310 K, to give a final 1.5 : 1 mol ratio.The formation of adduct 4-9-EtG was confirmed by following the new set of peaks, and up to 90% yield of adduct was obtained when 1.5 mol equiv.nucleobase solution was added.However, no adduct was found when 1.5 mol equiv. of 5′-AMP was added to complex 4, even after 24 h incubation at 310 K. Reactions of double-helical calf thymus DNA (ct-DNA, 32 µg mL −1 ) and plasmid DNA pBR322 (28 µg mL −1 ) with complex 4 in various molar ratios (r i = 0.05-1, r i = the molar ratio of free Ru complex to nucleotide phosphates at the onset of incubation with DNA) were studied.Very low amounts of ruthenium (5-7% of initial Ru) were found in the samples of DNA treated with complex 4 for 24 h.No significant changes in the mobilities of supercoiled (sc) or open circular (oc) form of plasmid DNA were observed even when incubated with high concentration of complex 4 (r i = 1, Fig. S7, ESI †).DNA is thought to be a cellular target for the en complex 7 (Fig. 1). 52,53However, for the substituted-en complex studied here, no obvious unwinding of DNA was observed after coincubation of ct-DNA with complex 4, suggesting that binding is weak, nor changes in the ratio of sc and oc forms of plasmid DNA, suggesting that complex 4 does not cleave DNA.
For state 5, with ring-slipped coordinated η 2 -p-cymene, the introduction of water into coordination sphere was necessary, while highly distorted complexes without coordinated water were found ca. 100 kJ mol −1 higher in energy.It is notable that the Ru atoms of all complexes in state 5 are coordinated to the amide oxygen atom of NADH, while only weakly bound to the (hydridic) CH 2 of NADH, giving a Ru-H distance of 3.11-3.12Å for R 1 ,R 2 = Me,H (1); Et,H (3); Naph,H (6) and 3.06-3.07Å for R 1 ,R 2 = Bz,H (4) and 4-F-Bz,H (5).For R 1 ,R 2 = Me,Me (2) the calculations revealed a true bonding of the (hydridic) CH 2 , with a Ru-H distance of 1.99 Å.The results obtained are Fig. 8 (Top) Reduction cycle for conversion of NAD + to 1,4-NADH via transfer hydrogenation with formate as the hydride donor.(Bottom) DFT energy profile for the formation of Ru formate species, Ru hydride complex and hydride transfer from ruthenium; brown line, complex 1; red line, complex 2; blue line, complex 3; green line, complex 4; purple line, complex 5; black line, complex 6.Sets of calculated structures of states 1-7 are supplied in the ESI and illustrated graphically in Fig. S8 † for complex 2. To calculate the energy of states 1 and 7, the energies of the states represented in pdb files 1 and 7 were added to the energies calculated for NAD + and NADH.
shown in Fig. 8. Four general conclusions can be drawn from these data: (a) there is a strong interaction between the [(η 6 -p-cym)Ru(TsEn(R 1 ,R 2 ))(OH 2 )] + cation and NAD + and NADH molecules, leading to a stabilisation of the cationic form by 60-70 kJ mol −1 for NAD + (41 kJ mol −1 for R 1 ,R 2 = 4-F-Bz,H ( 5)) and 130-150 kJ mol −1 for NADH; (b) depending on the N-substituent, the species of the lowest energy is either [(η 6 -pcym)Ru(TsEn(R 1 ,R 2 ))(HCOO − )]•NAD + (R 1 ,R 2 = Me,H (1); Bz,H (4); 4-F-Bz,H (5) and Naph,H (6)) or [(η 6 -p-cym)Ru(TsEn(R 1 ,R 2 )) (OH 2 + (R 1 ,R 2 = Et,H (3) and Me,Me (2)), the difference between them being only 6-12 kJ mol −1 ; (c) the effective NADH-hydride coordination for bulky R 1 ,R 2 = Me,Me (2) lowers the energy, relative to the state of lowest energy, of the species with coordinating NADH by 30-40 kJ mol −1 , compared to other complexes; (d) the formation of the state with the Ru-H hydride bond, including the twist of formate and the elimination of carbon dioxide, corresponds to the highest energy step.These four factors seem all to influence the turnover.The energy barriers and optimized structures for the seven states of complex 2 in the cycle with NAD + are listed in Table S5 and illustrated graphically in Fig. S8.† The structure files for the remaining complexes are supplied as ESI.†

Conclusions
In this work, we have synthesised a series of new Ru II complexes of the type [(η 6 -p-cym)Ru(N,N′)Cl] where N,N′ are monosulfonamide chelating ligands derived from tosylethylenediamine, with either alkyl (Me,H (1); Me,Me (2); Et,H (3)) or aryl (Bz,H (4); 4-F-Bz,H (5); Naph,H (6)) substituents on the terminal N.These substituents have a significant effect on the rate of transfer hydrogenation of coenzyme NAD + with formate as hydride donor as determined by NMR and UV-vis spectroscopy.In general, the bulkier aromatic substituents gave rise to faster hydrogenation rates (Table 3).DFT calculations provided insight into the mechanism of hydride transfer form formate to NAD + involving initial coordination of formate followed by transfer of hydride to ruthenium and then to NAD + with release of CO 2 .The calculations suggested a preorganization of the initial aqua complex, formate and NAD + involving T-shaped adenosine NH-tosyl stacking, H-bonding of the NH of the chelated ligand and phosphate O of NAD + , H-bonding between formate and water, and between formate and the pyridine ring of NAD + .They also indicated strong interactions with NADH involving T-shaped adenosine NH-tosyl stacking, as well as H-bonds to phosphate and (hydridic) CH 2 -tosylate O (Fig. S8, ESI †).
To investigate the possibility of achieving transfer hydrogenation mediated by formate in cells, we investigated the effect of formate on the antiproliferative activity of these complexes towards human ovarian cancer cells.In each case a dose-dependent increase in potency of the complexes (20-36%) was observed with increasing formate concentration over a range of non-toxic formate concentrations (0-2 mM).The complexes with aromatic substituents were the most potent, the benzyl complex 4 being as potent as the anticancer drug cisplatin (Fig. 6).In general, the most hydrophobic complexes were found to be the most biologically active.However, the activity does not correlate closely with total cell accumulation of Ru or with hydrophobicity (Table 3).Although DNA can be a target for related arene Ru II diamine complexes, it does not appear to be a target for these sulfonyl-en Ru II catalysts since we observe very weak binding to both calf thymus and plasmid DNA (Fig. S7, ESI †).
We showed that complexes 1 and 4 can generate high levels of ROS in A2780 human ovarian cancer cells, especially 4, the most potent complex.This is consistent with interference in cellular redox pathways and possible attack on NAD + when sodium formate is present.The enhancement of anticancer activity by low non-toxic dose of formate might be clinically useful since it introduces a new mechanism of activity which does not involve DNA attack, unlike the clinical drug cisplatin.Such a regime might therefore avoid some unwanted sideeffects.Formate itself is a natural biochemical molecule enriched in some cancer cells. 54However, more work remains to be done to investigate possible intracellular catalysis, especially since a range of metabolites might readily poison these catalysts in cells.
Collaborative Postgraduate Research Scholarship (WCPRS, funding for J. P. C. C.); J. A. W. and V. S. acknowledge support of the research initiative NANOKAT and the German Federal Ministry of Education and Research (BMBF under 05K14UK1) and are grateful to the Allianz für Hochleistungsrechnen Rheinland-Pfalz (AHRP) for providing CPU-time within the project TUK-SPINPLUSVIB.
We also thank Dr Ivan Prokes, Dr Lijiang Song, and Mr Philip Aston (University of Warwick) for their excellent assistance with the NMR and MS measurements.

Fig. 2
Fig. 2 ORTEP diagrams for complex 3. Ellipsoids are shown at the 50% probability level.All hydrogen atoms have been omitted for clarity.

Fig. 3
Fig. 3 Dependence of the 1 H NMR chemical shift of a tosyl proton (red) on pH* of aqua complex 4a.The red curve is the best fit to the Henderson-Hasselbalch equation corresponding to a pK * a of 9.73 ± 0.06.

Fig. 6
Fig. 6 IC 50 values (µM) for complexes 1-6 against A2780 human ovarian cancer cells (orange bars) and cellular accumulation of Ru in A2780 cancer cells at equipotent IC 50 concentrations in the absence of sodium formate (in green).

Fig. 7
Fig. 7 ROS in A2780 cells induced by complexes 1 and 4, FL1 channel detects total oxidative stress, and FL2 channel detects superoxide production.(A) Induction of ROS by complexes 1 and 4. (B) Four different populations induced by complexes 1 and 4 at equipotent IC 50 concentrations.p-Values were calculated after a t-test against the negative control data, **p < 0.01, ***p < 0.001.

Table 1
Ru II complexes studied in this work

Table 3
Turnover frequencies for transfer hydrogenation reactions using Ru complexes 1-6 as catalysts

Table 4
Retention times (t R ) of Ru II complexes 1-6 by RP-HPLC and cellular accumulation (at equipotent of IC 50 concentrations) in A2780 cells