Next Article in Journal
Assistance of a Person with Muscular Weakness Using a Joint-Torque-Assisting Exoskeletal Robot
Next Article in Special Issue
The Cytotoxic Activity of Diiron Bis-Cyclopentadienyl Complexes with Bridging C3-Ligands
Previous Article in Journal
Hierarchical Task Assignment and Path Finding with Limited Communication for Robot Swarms
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Ruthenium–Thymine Acetate Binding Modes: Experimental and Theoretical Studies

1
Dipartimento di Chimica Industriale “Toso Montanari”, Università di Bologna, Viale Risorgimento 4, 40136 Bologna, Italy
2
Dipartimento di Chimica “Giacomo Ciamician”, Università di Bologna, Via Selmi 2, 40126 Bologna, Italy
*
Author to whom correspondence should be addressed.
Appl. Sci. 2021, 11(7), 3113; https://doi.org/10.3390/app11073113
Submission received: 8 February 2021 / Revised: 26 March 2021 / Accepted: 27 March 2021 / Published: 31 March 2021
(This article belongs to the Special Issue Design of New Organometallic Bioactive Molecules)

Abstract

:
Ruthenium complexes have proved to exhibit antineoplastic activity, related to the interaction of the metal ion with DNA. In this context, synthetic and theoretical studies on ruthenium binding modes of thymine acetate (THAc) have been focused to shed light on the structure-activity relationship. This report deals with the reaction between dihydride ruthenium mer-[Ru(H)2(CO)(PPh3)3], 1 and the thymine acetic acid (THAcOH) selected as model for nucleobase derivatives. The reaction in refluxing toluene between 1 and THAcOH excess, by H2 release affords the double coordinating species κ1-(O)THAc-, κ2-(O,O)THAc-[Ru(CO)(PPh3)2], 2. The X-ray crystal structure confirms a simultaneous monohapto, dihapto- THAc coordination in a reciprocal facial disposition. Stepwise additions of THAcOH allowed to intercept the monohapto mer-κ1(O)THAc-Ru(CO)H(PPh3)3] 3 and dihapto trans(P,P)2(O,O)THAc-[Ru(CO)H(PPh3)2] 4 species. Nuclear magnetic resonance (NMR) studies, associated with DFT (Density Function Theory)-calculations energies and analogous reactions with acetic acid, supported the proposed reaction path. As evidenced by the crystal supramolecular hydrogen-binding packing and 1H NMR spectra, metal coordination seems to play a pivotal role in stabilizing the minor [(N=C(OH)] lactim tautomers, which may promote mismatching to DNA nucleobase pairs as a clue for its anticancer activity.

Graphical Abstract

1. Introduction

In the last decade, a plethora of novel developments on ruthenium complexes exhibited a remarkable efficiency as potential metal-drugs in anticancer treatments together with a consistent reduction of the related toxic issues exhibited by the former Pt-agents. The ruthenium complexes possessing anticancer activity exhibit a plethora of different structures, mainly belonging to coordination chemistry with labile chloride groups [Ru(II)(biphenyl)Cl(en)]+, where en = 1,2-ethylenediamine [1] or the efficient solvato-complexes [Ru(III)Cl4(DMSO)(imid)][imidH] [imid=imidazole] [2]; at the same time, organometallic piano stool compounds bearing the ubiquitous p-cymene [Ru(p-cymene)Cl(H2O)(PTA)]+, arene or cyclopentadienyl groups and a variety of chelate donor ligands as bipyridine moieties as shown by [Ru(η5-C5H5)(PPh3)bipy)][CF3SO3][bipy=bypyridyl][3]. Therefore, it is of fundamental importance to disclose more evidence on the structure–effect relationship. In this respect, great efforts have been made in ruthenium novel drugs [1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16] exploration, focused on the interaction modes towards DNA nucleic bases [17,18,19,20,21,22,23]. A useful approach to shed some light on the mechanism related to biological action effectiveness [24] is to study the coordination modes of a metal-anchored thymine acetic acid (THAcOH), proposed here as ligand to design potential bioactive ruthenium complexes. Due to chelate properties, the side arm carboxylate group may induce further stability, analogously to transferrin behavior in the inter-strand crosslinks. The reaction course will be explored by spectroscopic and structural investigations, DFT-calculated energy profiles and by comparing analogous reactivity with acetic acid. The intermolecular H-binding network in THAc-Ru systems, which is remarkably evident in the X-ray structure, prompt us to investigate the occurrence in solution. Evaluation of the 1H nuclear magnetic resonance (NMR) intermolecular H-bonding [25] signals has been suggested by the help of DFT-calculated energies, through simulation of the interactions of two THAc- anions, selected as model. By dealing with the mutagenic properties of metallodrugs in anticancer activity, Lippert’s pioneering work [21,22,23,26,27,28,29,30] suggested that proton relocation, obtained by acid treatment, deprotonation, or H-bond interactions is responsible for the stabilization of rare tautomeric forms. The proton relocation can be a crucial key in mismatching DNA nucleobase pairs, by impeding cell replication.

2. Materials and Methods

2.1. General

All reactions were routinely carried out under argon atmosphere, using standard Schlenk techniques. Solvents were distilled immediately before use under nitrogen from appropriate drying agents. Chromatography separations were carried out on columns of dried celite. Glassware was oven-dried before use. Infrared spectra were recorded at 298 K on a PerkinElmer Spectrum 2000 FT-IR (Fourier transform infrared) spectrophotometer (Waltham, MA, USA), electrospray ionization mass spectrometry (ESI MS) spectra were recorded on a Waters Micromass ZQ 4000 (Milford, MA, USA), with samples dissolved in CH3OH. All deuterated solvents were degassed before use. All NMR measurements were performed on a Mercury Plus 400 instrument (Oxford Instruments, Abingdon-on-Thames, UK). The chemical shifts for 1H and 13C were referenced to internal TMS. 1H, 13C correlation measured using gs-HSQC (gradient selected Heteronuclear Single Quantum Coherence) and gs-HMBC (gradient selected Heteronuclear Multiple Bond Correlation) experiments. All NMR spectra were recorded at 298 K. All the reagents were commercial products (Aldrich) of the highest purity available and used as received. [RuCl3∙xH2O] was purchased from Strem and used as received. Compound [Ru(H)2(CO)(PPh3)3] (1) was prepared by published methods [31].

2.2. Experimental Procedure

2.2.1. Synthesis of Complex. κ1O)-THA κ2(O,O)-THAc-[Ru(CO)(PPh3)2], 2

To a 100 mg of 1 were added upon stirring in refluxing toluene two equivalents of THAcH (40 mg, 0.218 mmol) for 240 min. After 15 min evolution of a gas was observed until the Ru-CO absorption at 1940 cm−1 disappeared. After drying under vacuum, the light brown solid was washed with Et2O (3 × 10 mL) to remove the released PPh3. The crystallization occurred in CDCl3 gave white crystals of 2 (90%, 97 mg). The same reaction was carried out by adding an excess of THAcH or two equivalents consecutively, until the starting material infrared (IR) absorption bands disappeared.
[RuC53H44O9P2; 988.15 g mol−1]; IR (KBr, cm−1) trans-2 ν: Ru-CO 1969 vs; κ2THAc, νasym C(O)O 1653s, νsym C(O)O 1630s; κ2THAc ν C(O)O 1521 vs. κ1THAc ν C(O)O 1361 vs, (Δν = 280 cm−1); 1H-NMR (400, CDCl3): δ 1.74 (s, 6H, Me κ1THAc + Me κ2THAc), 1.82 (aq, ∑2JHP = 13 Hz CHαHβ, κ1THAc); 4.29 (aq, 5JHP = 89 Hz, 2JHH = 17 Hz, CHαHβ, κ2THAc [C(O)Oκ1THAc]; 5.85 (s, 1H,CH-methyne); 6.80–7.70 (m, 30H, 2 PPh3); 7.90 (s, br,1H, NH), 8.10 (s, br, 1H, NH). 13C{1H}(100 MHz, CDCl3) δ 205.6 (t, CO 2JCP = 15Hz,); 177.4 (COO THAc); 164.6 150.5 (CO, THAc); 141.7 (Cq THAc); 126–138 (PPh3); 109.7 (CH); 49.5 (CH2); 14.7 (Me) cis-2 ν: Ru-CO 1969 vs; κ2THAc, νasym C(O)O 1658s, νsym C(O)O 1635s; κ2THAc ν C(O)O 1536 vs (Δν = 111 cm−1) κ1THAc ν C(O)O 1358m, (Δν = 289 cm−1); 1H-NMR (400, CDCl3): δ 1.05 (s, 6H, Me κ1THAc + Me κ2THAc), 1.43 (s, 2H, CH2N); 3.27 (aq, 5JHP = 828 Hz, 2JHH = 16 Hz, CHαHβ, κ2THAc [C(O)Oκ1THAc]; 4.43 (s, 1H, CH-methyne); 6.80–7.70 (m, 30H, 2 PPh3); 7.67 (m, 2H, 2NH). The κ1- and κ2- coordination are undistinguishable because of the mutual scrambling at room temperature that affect either 1H or 13C NMR spectra. 31P{1H}-NMR (161.9 MHz, CDCl3): cis 57.0 (dd, 2JPa-Pb = 15 Hz, 2P cis). 177,168 (COTH)TRANS 170, 173 (COTH)CIS 163 (COO)TRANS, 162 (COO)CIS, 150, 141, 125, 112 (J2CP = 4 Hz) PPh3 109.5 (CH) TRANS 109 (CH) CIS 110 (CMe)TRANS 49 (broad NCH2)TRANS 52 (broad NCH2)CIS 39 (Me)TRANS 32 (Me)CIS.

2.2.2. Synthesis of Complexes mer–κ1(O)THAc-[RuH(CO)(PPh3)3] (3a, 3b) + trans(P,P)2(O,O)THAc -[RuH(CO)(PPh3)2], 4

To a 100 mg of 1 was added upon stirring in refluxing toluene an equimolar amount of THAcH (20 mg, 0.109 mmol). After 10 min evolution of a gas was observed until the Ru-CO absorption at 1940 cm−1 disappeared. After 30 min the solid was dried under vacuum, the light brown solid was washed with Etp (3 × 10 mL) to remove the released PPh3 and chromatographed on a celite column (7 × 1.5 cm). By elution with neat Et2O a first pale green fraction was collected giving the monohapto acetate derivatives 3 (72 mg, 60%), whereas the elution with CH2Cl2 gave a purple fraction prevalently constituted by the dihapto species 4 (27 mg, 30%).
By changing the reaction times to 60 min the major fraction mainly afforded the purple fraction constituted by the chelate derivative 4 (78 mg, 85 %); merκ1(O)THAc-[Ru (CO)H(PPh3)3] (3a,3b), green powdered solid; IR (KBr, cm−1) ν: Ru-CO 1918 vs; νas C(O)O κ1THAc 1683, νs C(O)Oκ1THAc 1654s (Δν = 175 cm−1), 1480 m (P-C) 1434 vs. (P-C); 1361m C(O)O; 1H-NMR (400.0 MHz, CDCl3: δ 6.8–7.7 (m, 45H, 3 PPh3); 5.7 (s br, 1H, CH THAc), 5.4 (s br, 1H, CH THAc); 3.9 (s br, 2H, CH2); 3.1 (s br, 2H, CH2); 1.6 (s, 3H, Me); −6.1 (dt, JP-H = 112 Hz, JP-H = 24Hz, 1H); −7.2 (dt, JP-H = 112Hz, JP-H = 24 Hz, 1H); 13C{1H}-NMR (100.6 MHz, CDCl3: δ 205.6 (s br, CO); 203.9 (s br, CO); 177.3 (COO); 170.9 (COO); 167.9 (CO(1) THAc); 164.7 (CO(1) THAc); 151.5 (CO(2) THAc); 150.5 (CO(2) THAc); 142.8 (Cq THAc); 141.6 (Cq THAc); 135.3 (t, JCP = 5 Hz, PPh3); 134.5 (t, JPC = 22 Hz, Cq PPh3); 130.7 (s, CHpara PPh3); 128.0 (t, JPC = 4 Hz, CH PPh3); 108.9 (CH THAc); 108.5 (CH THAc); 50.5 (CH2); 49.5 (CH2); 12.8 (Me); 12.7 (Me); 13C{1H}-NMR (100.6 MHz, CDCl3): δ 31P{1H}-NMR (161.9 MHz, CDCl3: δ 41.2 (s br, 1P); 36.5 (s br, 2P).
trans(P,P)2(O,O)-[Ru(CO)H(PPh3)2THAc], 4 waxy reddish solid; IR (KBr, cm−1) ν 1963 vs. (RuCO); 1683 vs, b; 1443, 1434 (P-C); 1H-NMR (400,0MHz,CDCl3: δ 6.8–7.7 (m, 30H, 2 PPh3); 5.4 (s br, 1H, CH THAc); 3.9 (s br, 2H, CH2); 1.6 (s, 3H, Me); −17.0, (t, 2JHP = 20.0 Hz);13C{1H}-NMR (100.6 MHz: δ Ru-CO not observed, 179.8 (s, COO κ2-THAc); 163.9 (s, CO(1) κ2-THAc); 149.8 (s, CO(2) κ2-THAc); 141.4 (s, Cq κ2-THAc); 124.0–136.0 (PPh3) 112.0 (s, CH κ2-THAc); 64.4 (s, CH2 κ2-THAc); 13.4 (s, Me κ2-THAc); 31P{1H}-NMR (161.9 MHz, CDCl3: δ 38.7 (d, ∑2JPH = 14 Hz, 2P); Ms-ESI (MeOH, m/z): 1043 {20%, M=[Ru(CO)(THAc)2(PPh3)2]++Na} 859 {25%, [Ru(CO)(THAc)(PPh3)2]++Na} 837 {100%, [Ru(CO)(THAc)(PPh3)2]++H}; 653 {15%, [Ru(CO)(PPH3)2]++H}.

2.3. Reactions of 1 with Acetic Acid

For optimizing the sake of product selectivity, we run the reactions by refluxing in different solvents and by using different quantities of reactant (Table 1).
A 200 mg of 1 has added upon stirring in refluxing solvents with variable quantities of acetic acid (25 μL ~2 eq). After 5 min, evolution of a gas was observed until the Ru-CO absorption at 1940 cm−1 disappeared. After the reaction times displayed in Table 1 the solid was dried under vacuum and the light brown solid was washed with Etp (2 × 10 mL) to remove the released PPh3 and filtered on a celite column (8 × 1.5 cm). By elution with neat Et2O two fractions were collected giving a mixture of acetate adducts with prevalence of the desired product. Table 1 summarizes the reaction conditions of different reactions run by seeking the suitable conditions to reach the maximum conversion and selectivity.
mer-κ1(O)Ac-[Ru(CO)H(PPh3)3] 5 pale green oily solid; IR (KBr, cm−1) ν: 1923 vs. (CO); 1H-NMR (400.0 MHz, CDCl3): δ 6.80–780 (m, 45H; PPh3); 1.28 (s, br 3H, CH3COO); −7.13 (dt, 2JPPcis = 24 Hz, 2JPPtrans = 105 Hz, Ru-H); 13C{1H}-NMR (100,6MHz): δ Ru-CO 204.3 (t, 2JPC = 13.9 Hz, COtrans), 203.0 (dd, ∑S2JPaC = 15.2 Hz, COcis), 189.3 (COO), 24.6 (Me); 31P{1H}-NMR (161.9 MHz, CDCl3): δ 36.3.
trans(P,P)-κ2(O,O)Ac-[Ru(CO)H(PPh3)2] 6, deep purple waxy solid; IR (KBr, cm−1) ν: 1929 vs. (CO); 1H-NMR (400.0 MHz, CDCl3): δ 6.80–780 (m, 30H; PPh3); 0.62 (s, br 3H, CH3COO);-16.50 ppm (2JHP = 20 Hz); 13C{1H}-NMR (100,6MHz): δ Ru-CO, 186.1 (COO), 135.5–126.9 (PPh3); 30.3 (Me).31P{1H}-NMR (161.9 MHz, CDCl3): δ 38.4.
trans(P,P)1(O)Ac,-κ2(O,O)Ac-[Ru(CO)(PPh3)2] 7a, pink solid; IR (KBr, cm−1) ν: 1940 vs. (CO), 1625 COO, 1465 (κ2) 1434sy (κ1) (COO); 1H-NMR (400.0 MHz, CDCl3): δ 7.80–6.80 (m, 30H; PPh3); 0.69 (s, br 6H, CH3COO); 13C{1H}-NMR (100,6 MHz): δ 206.0 (t, 2JPC = 14 Hz, COtrans), 186.0 (COO), 133.5 (t, JPC=8Hz, Cq PPh3); 133.3 (t, JPC = 5 Hz, CH PPh3); 129.5 (s, CHpara PPh3); 127.2 (t, 2JPC = 5 Hz, CH PPh3); 21.9 (Me). 31P{1H}-NMR (161,9MHz, CDCl3): δ 39.1.
cis(P,P)1(O)Ac,-κ2(O,O)Ac-[Ru(CO)(PPh3)2] 7b, violet oily solid; IR (KBr, cm−1) ν: 1943 vs. (CO), 1515, 1434 (COO); 1482, 1466 (P-C) 1H-NMR (400.0 MHz, CDCl3): δ 7.84–6.95 (m, 30H; PPh3); 0.52 (s, br 6H, CH3COO); 13C{1H}-NMR (100.6 MHz): 204.7 (dd, 2JPaC = 19 Hz, 2JPbC = 16 Hz COcis), 189.30, (COO), 189.27 (COO), 133.5 (t, JPC = 8 Hz, Cq PPh3); 133.3 (t, JPC = 5 Hz, CH PPh3); 129.5 (s, CHpara PPh3); 127.2 (t, 2JPC = 5 Hz, CH PPh3); 22.6 (Me), 22.3 (Me). 31P{1H}-NMR (161.9 MHz, CDCl3): δ 35.5.

2.4. X-ray Crystallography

The X-ray intensity data were measured on a Bruker Apex II CCD diffractometer. Cell dimensions, and the orientation matrix was initially determined from a least-squares refinement on reflections measured in three sets of 20 exposures, collected in three different ω regions, and eventually refined against all data. A full sphere of reciprocal space was scanned by 0.3° ω steps. The software SMART was used for collecting frames of data, indexing reflections and determination of lattice parameters. The collected frames were then processed for integration by the SAINT program, and an empirical absorption correction was applied using SADABS. The structures were solved by direct methods (SIR 97) and subsequent Fourier syntheses and refined by full-matrix least-squares on F2 (SHELXTL) [32,33,34,35] using anisotropic thermal parameters for all non-hydrogen atoms. The structure contained significant void space (527 Å3) and residual electron density that could not be meaningfully modelled; hence the SQUEEZE routine of PLATON was employed.
Cambridge Crystallographic Data Centre (CCDC)-1,868,289 contains the supplementary crystallographic data for this paper. These data can be obtained free of charge at www.ccdc.cam.ac.uk/conts/retrieving.html, accessed on 17 September 2018 (or from the Cambridge Crystallographic Data Centre, 12 Union Road, Cambridge CB2 1EZ, UK; Fax: +44-1223-336-033; e-mail: [email protected]).

2.5. Computational Details

DFT calculations were performed using the Molpro2010 software [36]. Due to the size of the involved species, geometry optimizations and free energy calculations (at 298 K) were undertaken with the small LANL2DZ [37] basis and the B3LYP [38] functional. Additional single point energy calculations with the larger def2-TZP [39] basis and the same B3LYP functional were performed at the previously optimized geometries.
The final energy of each structure, used to evaluate the relative free energies of the various products and intermediates, was built by summing the difference between the LANL2DZ electronic and free energies to the def2-TZP electronic energy. All calculations were performed in vacuo, i.e., commonly without considering the effect of the dielectric constant of the solvent unless otherwise stated.

3. Results

3.1. Reactions of [Ru(CO)H2(PPh3)3] 1 with THAcOH

The Ru-hydride mer-[Ru(CO)H2(PPh3)3] [40], 1 reacted in refluxing toluene within 240 min with thymine-acetic acid (THAcH) excess, by molecular hydrogen release until no 1H NMR signals for residual Ru-hydride in the crude solution were observed. Attribution to a structure, showing IR frequency at 1970 cm−1 for Ru-carbonyl function and simultaneous mono- and di-hapto-acetate-Ruthenium coordination like κ1(O)-THAc, κ2(O,O)-THAc[Ru(CO)(PPh3)2], 2 has been suggested on the basis of the IR frequency pattern. Namely, antisymmetric, and symmetric stretching vibrations of the COO- group of acetate at 1653 and 1363 (Δν = 290 cm−1) indicate κ1-THAc coordination, whereas narrower frequency differences (between ν = 1630 and 1532, Δν = 98 cm−1) are typical for κ2-THAc chelate mode. Species 2 in CDCl3 showing a single sets of 1H NMR resonances has been assigned to the less sterically hindered transP,P-configuration on the basis of stability considerations reported for similar complexes [41,42]. The 1H NMR signals at δ 1.74, and 5.85 ppm, are respectively attributed to methyl and methyne thymine substituents, suggesting facile fluxional exchange at room temperature between κ1- and κ2-THAc rings, which may plausibly occur through a penta-coordinated transient species, in which the THAc ligands coordinatively exchange by a merry-go-round (THAc) mutual interconversion. Then, the merTHAc- stereogeometry designation is independent from their hapticity mode. The 1H NMR signal at δ 2.35 is associated to the κ2-NCH2 group in transP-P-merTHAc-form, whereas a multiplet centered at 4.29 ppm showing a typical roof effect (doublet of doublet CHαHβ, 2JHH = 17; 2JHH = 90 Hz), is attributed to the NCHαHβ methylene moiety of the κ1-carboxy moiety. DFT molecular modelling calculations reveal that the Ru-OC(O) carboxy rotation is hampered by the axial phosphine steric crowding. By contrast, the absence of detectable JH-H coupling in the case of the dihapto-coordination may be attributed to the methylene group of chelate carboxy function, lying in planes that are mutually perpendicular (Karplus effect) [43]. This reciprocal disposition was evidenced both by the DFT transP-P- calculated isomer (94° in CDCl3 solution) and by the X-ray crystal structure of the cis isomer (87° A) (Figure 1). DFT calculations predicted also a low energy interconversion process from transP,P-merTHAc to cisP,P-facTHAc species (ΔE = 3.7 KJ mol−1 in vacuo), which decreased to 0.8 kJ mol−1, if the calculations take into account of the solvent effect (CDCl3). We supposed that Ru-O cleavage of the chelate-THAc could promote PPh3 migration from axial to equatorial position in a pentacoordinated intermediate species, as observed in the case of the precursor 1 solicited by light [44,45]. Accordingly, complete interconversion from transP,P-merTHAc to cisP,P-facTHAc isomer occurs at r.t. in CDCl3 within 3 h. Two 31P NMR doublets centered at δ 46.5, δ 43.2; (JPP = 25 Hz) are representative for cis-fac 2 form, while a broad doublet detected at 35.9 (2JP,P = 49 Hz) is the only observed residual signal due to the trans-mer form. The rapid mutual room temperature scrambling of the κ1- and κ2- coordinated THAc fragments made undistinguishable both the Me and the CH methyne signals in 1H or 13C NMR spectra. On the contrary, the THAc-methylene moiety and the phosphine 31P NMR signals resulted both in being non-equivalent, since they were affected by punctual different molecular environments. No-planar disposition of the thymine rings altered the global symmetry, therefore was responsible for doublets related to 1H NMR mutual coupling even in trans-configuration. As expected, the IR frequency pattern of the cis form was totally analogous to the trans-, but slightly lower-shifted. Due to the scarce solubility in less polar solvent, the 13C NMR spectrum, which was recorded within 24 h for resolution requirements, exhibited a broad triplet signal at δ 205.6 for the Ru-CO of the former transP,P-merTHAc isomer next to a multiplet at δ 203.9 attributable to the growing cisP,P-facTHAc species (Scheme 1). The 13C signals were attributed in a spectrum of trans:cis = 2:1 mixture, where the growing signals, related to the cis form, were adjacent to the trans one. Both the isomers were unstable in solution, and exhibited complete decomposition within 36 h at room temperature in polar coordinating solvents as DMSO (dimethylsulphoxide).
The cisP,P-facTHAc configuration in CDCl3 is calculated to be almost isoenergetic and stabilized by a multi-intermolecular H-binding network, which is clear in the X-ray crystal packing. In the crystal environment, the κ1-THAc ligand is able to intermolecularly couple through Watson–Crick N-H...O interactions with the chelate κ2 (O,O)-THAc counterpart of the adjacent molecule, forming a supramolecular metallacycle, reinforced by antiparallel thymine ring π–π stacking together with a phosphine phenyl C-H...O interaction [12]. In the crystal packing (Figure S2), symmetry related molecules of 2 are connected via two pairs of N-H...O hydrogen bonds between κ1(O)THAc -and κ2 (O,O)THAc ligands of neighboring molecules thus forming dimers centered about inversion centres. Further stabilization of the dimer is achieved by π–π stacking of the thymine rings (centroid-centroid distance 3.65 Å) (Figure 2). In this way the thymine-1-acetic ligand preserves the N-H–O interactions observed in the crystal packing of thymine-1-acetic acid. A similar self-assembly showing intermolecular N-H–O hydrogen bonds and the π–π stacking of the thymine residue has been previously observed, for example, in the paddle-wheel dicopper complex [Cu2(O(O)CCH2-T)4(DMSO)2]. The supramolecular architecture is completed by an intricate intermolecular C-H–O network involving the hydrogen atoms of two phenyl rings belonging to two PPh3 groups of two different molecules and the uncoordinated O atom (O6) of the monodentate THAc ligand. An analogous self-assembly interaction was observed in the THAcOH crystal structure [46,47,48,49]. Therefore, the cis-geometry of 2 could also be ascribed to H-bonding interactions between the THAc-rings and the low solubility exhibited by the cis-form attributed to the inter-molecular H-bonding network in the solid packing as inferred by the small energy gap (3.7 kJ mol−1) between the trans and cis configurations, calculated by DFT. It is worth noting that the Ru-CO function lies on the same side (syn) of the C(O)O carboxy moiety, in accord with the DFT calculations, which designate this stereogeometry as more stable in the observed cis-fac and in the unstable trans-mer form (See Scheme 4 for calculations). The NMR spectra relative to trans-2 and cis-2 isomers both show the presence of many solvents, which have been used during the synthesis (toluene), in chromatographic processes (CH2Cl2 and Et2O), or in multi-solvent stratification procedures for crystallization (CH2Cl2, Et2O and Etp). Fluxional behaviour, that plausibly promotes interconversion from trans to cis forms, implies generation of pentacoordinated intermediates, which in turn are stabilized by coordinating solvents, as diethylic or methyl-cyclopentyl ether. As shown by 1H NMR spectra, species 2 strongly attract solvent molecules up to the second coordination sphere, exhibiting a solvatochromic effect in solution, and solvent inclusion in the crystals, selected for X-ray diffraction studies.

3.2. Studies on the Reaction Path: Isolation of the Intermediates mer–κ1(O)-THAc [Ru(CO)H(PPh3)3] 3a, 3b and trans(P,P)-[κ2(O,O)-THAc-[Ru(CO)H(PPh3)2] 4

With the purpose to investigate the reaction course of 2, a stepwise addition of two subsequent equivalents of THAcH intercepted two intermediates, which were spectroscopically characterized and analyzed by DFT calculations. The reaction (1:1 molar ratio) were stopped after 30 min. and the dried residue was eluted by Et2O on a celite pad giving two distinct fractions. The nature of the more abundant pale green band (60%) has been assigned to a monohapto-coordinated mer–κ1(O)THAc-[Ru(CO)H(PPh3)3] 3a, 3b (Scheme 1) on the basis of a single Ru-CO IR absorption at ν = 1918 cm−1 and the characteristic acetate band pattern. The 31P{1H} NMR broad signals at δ 41.2 and 36.5 in reciprocal 2:1 ratio is attributed to the phosphine ligands by the hindered rotation about the CO-Ru rotation. The broad growing doublet centered at δ 38.7, with residual Ru-H coupling (∑JPH = 15 Hz), supports the formation of the transP,P- κ2-(O,O)-THAc, 4 as the more stable species (Scheme 1). The 1H NMR Ru-H triplets centered at −6.1 and −7.2 ppm (Figure 3), (2JH-Ptrans = 111 Hz and 2JH-Pcis = 24 Hz) indicate that hydride substitution mainly occurred trans to carbonyl ligand, which is responsible for stabilizing the donor THAc ligand. The second eluted reddish fraction is the result of acetate chelation by concomitant release of PPh3. THAc-chelation has confirmed by a strong Ru-carbonyl IR band at ν 1960 cm−1 and by a narrower separation of dihapto-acetate IR bands at 1682 s, 1656 s cm−1 (Δν = 26 cm−1). The highfield-shifted 1H NMR hydride triplet at –17.0 ppm (t, 2JHP = 20.0 Hz) (Figure 3), due to the notable electron-donor dihapto ligand, is compatible with transP,P-phosphine location, analogously to the κ2-O(O)-[Ru(CO)H(COH)(PCy3)2] complex (−18.5 ppm) [50].
Although the first reaction step, leading to κ1 intermediates, was not energetically favored, the energies required were still accessible at the experimental conditions. The distinct cis conformers 3a, 3b exhibited comparable energy, but were separated by a quite high barrier, as predicted by DFT calculations (Scheme 2). A third isomer, having κ1-(O) trans-located with respect the PPh3, was predicted at too high energy to be observed. The further evolution by PPh3 release, forming κ2-intermediates 4, exhibited energies lower than reagents only in the case of trans-PPh3, whereas cis-isomers always fall at higher energies. The last steps, leading to κ1-, κ2- species both in cis and trans- forms, are almost isoenergetic, as previously pointed out and fully supported by the experimental findings.
With the purpose to further investigate the reactivity path of 2, a reaction involving subsequent addition of THAcOH was conducted under the same conditions. The addition of the first equivalent of THAc promptly formed the pale green monohapto-3 complexes, which in turn converted the purple chelate 4. The transformation occurred spontaneously, even in the solid state. By contrast, addition of a further equivalent of THAcH to the isolated pale green monohapto-3 species, along the expected growing triplet at −17.0 of 4, showed the formation of two 1H NMR triplets in the interval δ −10.8–11.3 (2JHH = 11.7), likely attributable to transient bis-monohapto species. The related structure was suggested by DFT calculations as κ1(O, κ2(O)-[TH(CH2)C(O)O···H···OC(O)(CH2)TH][RuH(CO)(PPh3)2] [TH=thymine ring], in which the carboxy-carboxylic -(O)C(O)···HO(O)C- moieties of the two κ1-functionalized arms resulted intramolecularly stabilized by a pentametallacycle, showing very similar H-bonds (1.22 vs. 1.18 A). On the energetic scale, the structure appeared almost isoenergetic (−28.1 kJmol−1) with the trans form, used as reference, (Figure S18a). By further addition of THAc, the monohapto-, dihapto- species 2 was exclusively formed. The reaction requires one more hour at room temperature by using low-polar, non-coordinating solvents such as CHCl3 to avoid further evolution to new pentacoordinate species, formed by facile releasing of phosphine ligand, as observed in DMSO.

3.3. DFT Theoretical Calculations

The DFT calculations, [37] performed on both the monohapto-3 (blue line) and dihapto-4 (red line) rotamers, indicate a high torsional barrier (50 kJ mol−1) along the C-O bond for κ1(O)-3, forming distinct downward or upward THAc- conformations (Scheme 2). In the case of the favorite dihapto κ2-(O,O)-4 species a much lower energy barrier (10 kJ mol−1) was evaluated for the torsion along (carboxy)C-C(methylene) bond, confirming the equatorial THAc location as preferred for the more stable isomers, whatever coordinative mode would be adopted. Species κ1(O)-mer-3 was unstable and promptly converted to the entropically-driven chelate κ2(O,O)-trans-4. Although not totally hindered cisP,P- configurations of 3 and 4 were less favored because of steric contacts. Mixtures with different ratio of 3 and 4 intermediates were observed depending on the reaction time and the solvent nature. However, evident broadening of the 1H NMR triplet at −17.0 ppm, observed by keeping the solution of (3 + 4) in chlorinated solvent, may indicate a rearrangement by reposition of the phosphine ligands. This process seemed to be favored by the reduction of steric congestion then promoting complete isomerization. However, DFT calculations in the vacuo (Scheme 3) indicated the cisP,P-(−20.7 KJmol−1) and transP,P structures (−24.4 KJmol−1) to be almost isoenergetic, so that rearrangement to cisP,P may result prevalently by inter-molecular H-interactions as shown by the strongly stabilized Watson–Crick H-binding contacts in the X-ray crystal packing.

3.4. Reactions of 1 with AcOH and DFT Calculations

With the purpose to validate the proposed mechanism (Scheme 2), analogous reactions between 1 and acetic acid (AcOH) has been exploited, by changing reactant ratio, solvent polarity, and reaction time. The 1H-NMR spectrum of the dried solution displays two set of multiplets in 1:1.3 ratio, centered at −7.13 ppm (dt, 2JPPcis = 24 Hz, 2JPPtrans = 105 Hz) and a triplet at −16.50 ppm (2JHP = 20 Hz), which are attributed to mer1(O)Ac-[RuH(CO)(PPh3)3], 5 and trans(P,P)-κ2(O,O)Ac-[RuH(CO)(PPh3)2], 6 respectively (Scheme 4). The selectivity results strongly dependent on the temperature and solvent polarity. The DFT-calculated energy scale (Scheme 5) displays which are the favorite species. By adding two equivalent of acetic acid, the double coordinated trans(P,P)-κ1(O)Ac, κ2(O,O)Ac-[Ru(CO(PPh3)2], 7 appears as the main product. The presence of 13C NMR Ru-carbonyl triplet at 205.3 ppm, the 31P{1H} singlet at δ 39.1 and the resonance at δ 0.64 for both the 1H NMR methyl groups, strongly support acetate interconversion at room temperature, in accord with the analogous fluxional behavior reported for [κ1-(O)Ac-, κ2-(O,O)Ac-[Ru(CO)(PPh3)2] [51] by using a different reaction procedure.
In accord with what has been described for THAcH, subsequent additions of acetic acid, prior to giving birth to the signals relative to complex 2, unexpectedly lead to a 1H NMR signal at δ 14.0 ppm, likely due to a novel transient penta-coordinated species (Figure S18b). The structure of the transient intermediate is assigned to a bis-monohapto species of the type κ1(O), κ2(O)-[MeC(O)O···H···OC(O)Me][RuH(CO)(PPh3)2], showing the functionalized opened side arms interconnected by a H-bond, as suggested by DFT calculations (Scheme 6).
The following scheme, which demonstrates the intercepted intermediates in the sequenced reactions with AcOH, is supposed to simulate the reaction course with the THAcOH.

3.5. Intra or Inter-Molecular Bonding Network?

The observation of the unusually crowded 1H NMR spectral interval in the 8–12 ppm region suggested the need to check if the supramolecular H-binding network, exhibited by the X-ray packing of 2 occurred in solution also. The idea is to attribute the signals to the corresponding structures by the help of DFT free energies, to evaluate the related stabilization scale. The DFT energies (Scheme 7) for the various THAcH keto-enol monomers indicated that the -[N(H)-C(O)]- lactam-form (A, also referred as keto2 hereafter), was more stable than relative -[N=C(OH)]- iminol-lactim tautomers (C + D, referred to as enol4 and enol2, respectively) [26,27,28,29,30]. This order does not change for THAcK salts.
The energy of (A) shows a remarkable decrease DG = −16.2 kJ/mol for B) in the case of intra-molecular H-bond interaction between carboxy proton and the CO(2).
To investigate supra-molecular H-binding networks, a model by using κ2(O,O)-THAc-K+ to simulate the Ru-coordination sphere was adopted, since the full Ru-coordinated molecules were too large to be handled by computations.
In spite of the bonding nature, we proposed substituting κ2-(O,O)THAc-[Ru(CO)(PPh3)2] with the κ2(O,O)-THAcK+ thymine -acetate ligand, with the purpose of achieving feasible DFT calculations. In Figure 4 we compare the relative space filling models of the THAc-Ru fragment belonging to species 4 with THAcK+. We chose to freeze the rotamers into two distinct perpendicular and parallel limit conformations (Figure 2). The ability of the thymine ligands to interact each other is related to the external disposition of the THAc rings with respect to the crowded volume spanned by the Ru-coordination sphere. Although the steric encumbrance of THAcK compared to κ2(O)THAc-[RuH(CO)(PPh3)2 ] were very different, the molecular modelling evidenced that the thymine rings were external enough to reciprocally interact (Figure 4).
DFT calculations suggested a lot of dimerization possibilities, which are sketched in Scheme 8. All but the dimeric structures involving two enol forms (enol4 and enol2), which are found at energies larger than ~ 60 kJmol−1 were not considered for the assignment in 1H-NMR spectrum.
We selected the space filling figure of trans- κ2(O,O)-[RuH(CO)(PPh3)2], 4 as model of the (3+4) mixture solution which exhibits various 1H NMR signals in the appropriate range for H-bonding interactions. We analyzed the 1H NMR spectra (Figure 5), which contains different reaction mixtures after approximately after 45 min showing a (κ1- + κ2-) ratio rather equivalent (κ1-3: κ2-4 = 1.3:1). The weak sharp 1H NMR signal at δ 11.1 was attributed to the inter-molecular anti- enol-carboxy dimer (entry m, −20.4 kJmol−1, Scheme 8). Similarly, the three small broad signals in the interval of δ 9.2–8.8 ppm were assigned to the hydroxy function of the enol-lactim -N=C2,4(OH)- species, mutually H-bound to the keto- tautomers (entries cf). The 1H NMR signal at 8.3 ppm, which appears significantly more intense (5:1 with respect the C-H of thymine ring), likely belonged to the [NH-(O)C] Watson–Crick H-bonds in Scheme 6. We assumed the four overlapped N-H-N signals (7.72–7.93 ppm) to be a measure of the concomitant enol- species H-bound to the more stable keto-forms. The signals attributed to hydroxy functions were indicative for the lactim tautomers, including those involved into intramolecular H-bonding structures (entries F, G Scheme 7).
In the case of the species κ1(O)THAc-, κ2(O,O)THAc-[Ru(CO)(PPh3)2], 2 no evidence of tautomerization has been observed in the 1H NMR spectra (Figures S5a,b and S7a,b), since the NH-(O)C signals at δ 7.8 and 7.9 ppm indicate exclusively the formation of amino-carbonyl intermolecular bonds. To carry out this investigation we selected CDCl3, which concomitantly showed low polarity (χ = 4.8) and scarce coordination ability [52]. Considering the low solubility of the observed mixture and the faint response of the X···H–O interactions, by comparison with the intensities of the other signal, the 1H NMR spectra were recorded from freshly prepared samples (within 15 min) to exclude facile isotopic exchange processes with deuterated solvents. We are convinced that Ru-carboxy coordination played a pivotal role to induce proton relocation, boosting the tautomerization from -[N(H)C(O)]-amido to –[N=C(OH)]-iminol form, which finally resulted in being energetically stabilized by the inter-molecular H-binding network. The prevalence keto distribution (3.2:1) is correlated to inter-connections shown by the dimeric species, proposed as case-study model, and reported in Scheme 7. The NMR signal attribution to the X-H-O (X=N, O) intermolecular interactions was based on the electron withdrawing ability, the shielding effects of O–H–O (c–g, m and intramolecularly G) < N–H–O (h–i, n) < N–H–N (concomitant c–f) chemical shifts scale, the DFT-calculated energy scale and correlated to the resonance intensities. The carboxy-enol signal interaction (δ or 11.2) and the minor hydroxy signals around δ 8, which all belong to the minor enol-forms, were compared to the intense resonance at δ 8.3 which represented the prevalent keto isomers. The presence of iminol forms was also responsible for the concomitant growing of the NH-N signals (δ 7.8–8.0). As suggested by DFT-space filling evaluations, the H-bond interactions have merely been considered in the case of thymine rings anti-disposed to minimize the steric interference with the crowded trans-phosphines, therefore allowing the remarkable H-bond network. (Figure 5). A variety of 1H NMR spectra of intermediate mixtures, because they were analyzed at different reaction times, and although composed by different proportion of chelate-4 and monohapto-3 structures, showed a similar low-shifted signal pattern (two distinct examples are reported in Figure S11b,c). Considering the uncertainty of 10–15% due to integral determination and the different proportion of κ12 coordination of the non-conjugated metal fragments, the analyzed spectra exhibited an enol/keto ratio of 0.31 (Figure 5), compared to the enol/keto ratio as 0.13, reported in the literature [53].
The assignments were accomplished by the help of DFT-calculated energies for THAc anionic dimers selected as models, supporting the belief that the intermolecular interactions, promoted by H-binding network, played a key role in the solution also. The larger line widths, which were observed for the intermolecular O···H···O bonds (violet entries c–f) may have been due to the higher degree of freedom shown by OH···O(C) functions, which presented two O-lone pair able to generate in turn more than one H-bonded conformation. On the contrary, NH···O bond type (green entries n–l and n) or the enol-carboxy C(O)O···HO interactions at δ 11. 1 (m), which can occur intra- (G entry in Scheme 7) or intra-molecularly (entries h–l and n (Figure 5), resulted in being sharper since they are affected by a stronger H-bond strings network as in the case of N–H–N interactions (blue entries c–f).

4. Conclusions

The report deals with the reaction of 1 with THAcH excess, forming double coordinated mononuclear κ1(O)THAc-, κ2(O,O)THAc-[Ru(CO)(PPh3)2] species 2, which indicates the versatility of the coordination modes exhibited by thymine acetate. However, the X-ray characterized cisP,P-facTHAc isomeric form of 2 did not correspond to the transP,P-merTHAc structure observed in solution. A single 1H NMR signal for the THAc methyl moieties indicated a low-energy monohapto–dihapto interconversion at room temperature. DFT energetic studies confirmed the stabilization in solution of the cis configuration, governed by the reduced steric hindrance and by maximizing either intra- and inter-molecular H-bonding or π–interactions. Metal hydride upfield-shifted 1H NMR signals indicated the nature of the isolated intermediates 3 and 4, the rotamer energies of which were studied also by DFT-calculations. To elucidate the reaction course, similar reactions were run between 1 and acetic acid, confirming the proposed mechanism of cis-trans interconversion by the help of DFT calculations, despite the remarkable reduced steric requirements and the limited H-binding features. Upon ruthenium complexation, due to the further stabilization imparted by X–H–O (X=O, N) supramolecular H-binding network, the rare -(O)C-N=C(OH)- THAc-lactim tautomer results increased remarkably (0.31) compared to the predominant (O)C-NH-C(O) lactam form.
Lippert’s pioneering work, by dealing with the mutagenic properties of metallodrugs in anticancer activity, suggested proton relocation was responsible for the stabilization of rare tautomeric forms, being promoted by direct (N,O) metal coordination of the thyminate nucleobase. Herein, we further support the evidence that even non-conjugated thymine derivatives, once coordinated to ruthenium, may notably influence the stability of the minor iminol (0.31 ratio) compared to the tautomeric thymine equilibrium (pka = 10), where the formation of iminol species is promoted by acidic treatment. The occurrence of the latter tautomer, further stabilized by the preponderant H-binding network, is recognized to be active in DNA mismatching, limiting the replication of cancer cells [54,55,56,57].

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/app11073113/s1, Table S1: Crystal data and structure refinement for compound cis-2, Table S2: Most relevant hydrogen bonds for cis-2 [Å and °]; Figure S1: Molecular structure of cis-2 with the atom labelling, Figure S2: Molecular structure of cis-2 with π-π stacking and Watson-Creek intermolecular interaction, Figure S3: View down. Fragment of the crystal packing of complex 2 illustrating the intermolecular H bonding pattern. For the sake of clarity only the H atoms engaged in H bonds are shown. Ball and stick representation is used for the dimer arising by the strong N-H···O hydrogen bonding and for the atoms connected to it. H bonds are shown with blue dashed lines. The a axis of the crystal packing of cis-2. Black dotted lines indicate the intermolecular N-H···O hydrogen bonds, Figure S4: Fourier transform infrared (FT-IR) spectra of trans-κ1(O)THAc-, κ2(O,O)THAc-[Ru(CO)(PPh3)2] 2, Figure S5a: 1H-NMR (400 MHz) of cis-κ1(O)THAc-, κ2(O,O)THAc-[Ru(CO)(PPh3)2] 2 in CDCl3, Figure S5b: 1H-NMR (400 MHz) of cis-κ1(O)THAc-, κ2(O,O)THAc-[Ru(CO)(PPh3)2] 2 full spectrum in CDCl3, Figure S6: ESI-MS of cis-κ1(O)THAc-, κ2(O,O)THAc-[Ru(CO)(PPh3)2] 2, Figure S7a: 1H-NMR (400 MHz) of cis-κ1(O)THAc-, κ2(O,O)THAc-[Ru(CO)(PPh3)2] 2 in CDCl3, Figure S7b: 1H-NMR (400 MHz) of cis-κ1(O)THAc-, κ2(O,O)THAc-[Ru(CO)(PPh3)2] 2 aromatic region in CDCl3, Figure S8: FT-IR (KBr) κ1(O)THAc-[RuH(CO)(PPh3)3] 3a, b, Figure S9: 31P{1H} NMR (161 MHz) of 2 in CDCl3, Figure S10: FT-IR (KBr) of κ2(O,O)THAc-[RuH(CO)(PPh3)2] 4, Figure S11a: 31P{1H} NMR (161,9 MHz) κ1(O)-[RuH(CO)(THAc)(PPh3)3 3a,b + κ2(O,O)-[RuH(CO)(THAc)(PPh3)2] 4 in CDCl3, Figure S11b: 1H-NMR spectrum of 3+4 (ratio 3:2) showing iminol/keto ratio = 0.31 in CDCl3, Figure S11c: 1H-NMR spectrum of 3+4 (ratio 3:7) in CDCl3 showing iminol/keto ratio = 0.25, Figure S11d: Complexes formed as impurities during the preparative reduction of starting material 1, Figure S12: Electrospray ionization mass spectrometry (ESI MS) (MeOH) of κ2(O,O)THAc-[Ru(CO)H(PPh3)2] 4, Figure S13: 1H-NMR hydride region spectrum of 1+ κ1(O)Ac-[RuH(CO)(PPh3)2] 5+ κ2(O,O)Ac-[RuH(CO)(PPh3)2] 6 in CDCl3, Figure S14a: FT-IR of compound 5, Figure S14b: 31P{1H}-NMR spectrum of κ2(O,O)Ac-[RuH(CO)(PPh3)2] 6 species, Figure S15: 13C{1H}-NMR spectrum of 6 species, Figure S16: FT-IR (KBr) κ1(O)Ac-, κ2(O,O)Ac- [Ru(CO)(PPh3)2], 7, Figure S17: Carbonyl region in 13C{1H}-NMR spectrum of (trans+cis)(P,P)- κ1(O)Ac,- κ2(O,O)Ac-[Ru(CO)(PPh3)2] 7 in CDCl3, Figure S18a: portion of 1HNMR spectrum in hydride region of intermediate κ1(O)-, κ2(O)-[THCH2C(O)O···H···OC(O) CH2TH][RuH(CO)(PPh3)2][TH=thymine ring] A, Figure S18b: portion of 1HNMR spectrum in hydride region of intermediate κ1(O), κ2(O)-[MeC(O)O···H···OC(O)Me][RuH(CO)(PPh3)2] B, Figure S18c: 1H-NMR portion of the spectrum of A in CDCl3 displaying the CH2 moiety of the Et2O solvent, remarkably shifted upon Ru-coordination, Figure S18d: 1H-NMR portion of the spectrum of A in CDCl3 displaying the CH3 moiety of the Et2O solvent, shifted upon Ru-coordination

Author Contributions

Conceptualization, S.B.; methodology, S.C. and G.M.; software, R.T. and M.M.; formal analysis, S.C., M.M., C.B. and R.T.; investigation, S.B., S.C. and G.M.; resources, S.B.; data curation, S.C. and C.B.; writing—original draft preparation, S.B.; funding acquisition, S.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding. We thank the University of Bologna.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Acknowledgments

We thank the University of Bologna. Pietro Paolo Cristallini, Fabio Battaglia and Federica Zizzi are also gratefully thanked for the relevant and accurate experimental work carried out during their graduation internship.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bergamo, A.; Masi, A.; Peacock, A.F.A.; Habtemariam, A.; Sadler, P.J.; Sava, G. In vivo tumour and metastasis reduction and in vitro effects on invasion assays of the ruthenium RM175 and osmium AFAP51 organometallics in the mammary cancer model. J. Inorg. Biochem. 2010, 104, 79–86. [Google Scholar] [CrossRef]
  2. Alessio, E. Thirty Years of the Drug Candidate NAMI-A and the Myths in the Field of Ruthenium Anticancer Compounds: A Personal Perspective. Eur. J. Inorg. Chem. 2017, 12, 1549–1560. [Google Scholar] [CrossRef]
  3. Morais, T.S.; Santos, F.; Côrte-Real, L.; Marques, F.; Robalo, M.P.; Amorim Madeira, P.J.; Garcia, M.H. Biological activity and cellular uptake of [Ru(η5-C5H5)(PPh3)(Me2bpy)][CF3SO3] complex. J. Inorg. Biochem. 2013, 122, 8–17. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Allardyce, C.S.; Dyson, P.J. Ruthenium in Medicine: Current Clinical Uses and Future Prospects. Platin. Metals Rev. 2001, 45, 62–69. [Google Scholar]
  5. Gasser, G.; Ott, I.; Metzler-Nolte, N. Organometallic anticancer compounds. J. Med. Chem. 2011, 54, 3–25. [Google Scholar] [CrossRef] [PubMed]
  6. Frezza, M.; Hindo, S.; Chen, D.; Davenport, A.; Schmitt, S.; Tomco, D.; Ping Dou, Q. Novel metals and metal complexes as platforms for cancer therapy. Curr. Pharm. Des. 2010, 16, 1813–1825. [Google Scholar] [CrossRef] [Green Version]
  7. de Silva, D.O.; Pinto, D. Anticancer Natural Coumarins as Lead Compounds for the Discovery of New Drugs. Anticancer Agents Med. Chem. 2010, 10, 312–323. [Google Scholar]
  8. Bergamo, A.; Gaiddon, C.; Schellens, J.H.; Beijnen, J.H.; Sava, G. Approaching tumour therapy beyond platinum drugs: Status of the art and perspectives of ruthenium drug candidates. J. Inorg. Biochem. 2012, 106, 90–99. [Google Scholar] [CrossRef] [PubMed]
  9. Messori, L.; Merlino, A. Ruthenium metalation of proteins: The X-ray structure of the complexformed between NAMI-A and hen egg white lysozyme. Dalton Trans. 2014, 43, 6128–6131. [Google Scholar] [CrossRef]
  10. Motswainyana, W.M.; Ajibade, P.A. Anticancer Activities of Mononuclear Ruthenium(II) Coordination Complexes. Adv. Chem. 2015, 2015, 1–21. [Google Scholar] [CrossRef] [Green Version]
  11. Dragutan, I.; Dragutan, V.; Demonceau, A. The Expanding Chemistry of the Ruthenium Complexes. Molecules 2015, 20, 17244–17274. [Google Scholar] [CrossRef] [Green Version]
  12. Lin, Y.; Huang, Y.; Zheng, W.; Wang, F.; Habtemariam, A.; Luo, Q.; Li, X.; Wu, K.; Sadler, P.J.; Xiong, S. Organometallic ruthenium anticancer complexes inhibit human glutathione-S-transferase π. J. Inorg. Biochem. 2013, 128, 77–84. [Google Scholar] [CrossRef]
  13. Abid, M.; Shamsi, F.; Azam, A. Editorial of Special Issue Ruthenium Complex: The Expanding Chemistry of the Ruthenium Complexes. Mini-Rev. Med. Chem. 2016, 16, 772–786. [Google Scholar] [CrossRef] [PubMed]
  14. Kostova, I. Ruthenium complexes as anticancer agents. Curr. Med. Chem. 2006, 13, 1085–1107. [Google Scholar] [CrossRef] [PubMed]
  15. Sadler, P.J.; Guo, Z. Metal complexes in medicine: Design and mechanism of action. Pure Appl. Chem. 1998, 70, 863–871. [Google Scholar] [CrossRef] [Green Version]
  16. Ang, W.H.; Dyson, P.J. Classical and Non-Classical Ruthenium-Based Anticancer Drugs: Towards Targeted Chemotherapy. Eur. J. Inorg. Chem. 2006, 2006, 4003–4018. [Google Scholar] [CrossRef]
  17. Clarke, M.J. Ruthenium metallopharmaceuticals. Coord. Chem. Rev. 2003, 236, 209–233. [Google Scholar] [CrossRef]
  18. Goncharov, N.V. On the enzymatic activity of albumin. Russ. J. Bioorg. Chem. 2015, 41, 113–124. [Google Scholar] [CrossRef] [PubMed]
  19. Ruiz, J.; Villa, M.D.; Rodriguez, V.; Cutillas, N.; Vicente, C.; Lopez, G.; Bautista, D. A novel Metal-binding mode of thymine nucleobases: N(3) and O(4) chelation. Inorg. Chem. 2007, 46, 548–549. [Google Scholar] [CrossRef]
  20. Sletten, E. Metal binding to nucleic acid—A journey from the beginning. Inorg. Chim. Acta 2016, 452, 273–278. [Google Scholar] [CrossRef]
  21. Lippert, B. Multiplicity of metal ion binding patterns to nucleobases. Coord. Chem. Rev. 2000, 200–202, 487–516. [Google Scholar] [CrossRef]
  22. Pereira Lima, A.; Castro Pereira, F.; Pinheiro Almeida, M.A.; Santos Mello, F.M.; Carvalho Pires, W.; Monteiro Pinto, T.; Delella, F.K.; Felisbino, S.L.; Moreno, V.; Azevedo Batista, A.; et al. Novel piplartine-containing ruthenium complexes: Synthesis, cell growth inhibition, apoptosis induction and ROS production on HCT116 cells. PLoS ONE 2014, 9, e105865. [Google Scholar]
  23. Lippert, B.; Sanz Miguel, P.J. The Renaissance of Metal−Pyrimidine Nucleobase Coordination Chemistry. Acc. Chem. Res. 2016, 49, 1537–1545. [Google Scholar] [CrossRef] [PubMed]
  24. Nowak-Sliwinska, P.; van Beijnum, J.R.; Casini, A.; Nazarov, A.A.; Wagnieres, G.; van den Bergh, H.; Dyson, P.J.; Griffioen, A.W. Organometallic ruthenium(II) arene compounds with antiangiogenic activity. J. Med. Chem. 2011, 54, 3895–3902. [Google Scholar] [CrossRef] [PubMed]
  25. Steiner, T. The hydrogen bond in the solid state. Angew. Chem. Int. Ed. 2002, 41, 48–76. [Google Scholar] [CrossRef]
  26. Lippert, B. Rare iminol tautomer of 1-methylthymine through metal coordination at N(3). Inorg. Chim. Acta 1981, 55, 5–14. [Google Scholar] [CrossRef]
  27. Schollhorn, H.; Thewalt, U.; Lippert, B. Metal-Stabilized Rare Tautomers of Nucleobases. 2.′2-Oxo-4-hydroxo Form of Uracil: Crystal Structures and Solution Behavior of Two Platinum (II) Complexes Containing Imino1 Tautomers of l-Methyluracile. J. Am. Chem. Soc. 1989, 111, 7213–7221. [Google Scholar] [CrossRef]
  28. Lippert, B.; Neugebauer, D. Simultaneous Binding of Two Different Transition Metals to the DNA Model Base I-Methylthymine: The X-ray Structure of Bis[bis(p-1-methylthyminato-N3,04)-cis-diammineplatinum (II)] Silver Nitrate Pentahydrate. Inorg. Chim. Acta 1980, 46, 171–179. [Google Scholar] [CrossRef]
  29. Lippert, B. Effects of metal-ion binding on nucleobase pairing: Stabilization, prevention and mismatch formation. J. Chem. Soc. Dalton Trans. 1997, 21, 3971–3976. [Google Scholar] [CrossRef]
  30. Bruning, W.; Ascaso, I.; Freisinger, E.; Sabat, M.; Lippert, B. Metal-stabilized rare tautomers of nucleobases. 8. Promotion of rare cytosine tautomer upon complex formation with (dien)M2_ (M_/Pt, Pd). Inorg. Chim. Acta 2002, 339, 400–410. [Google Scholar] [CrossRef]
  31. Samouei, H.; Miloserdov, F.M.; Escudero-Adan, E.C.; Grushin, V.V. Solid-state structure and solution reactivity of [(Ph3P)4Ru(H)2] and related Ru(II) complexes used in catalysis: A reinvestigation. Oragnometallics 2014, 33, 7279–7283. [Google Scholar] [CrossRef]
  32. SMART & SAINT Software Reference Manuals, Version 5.051 (Windows NT Version); Bruker Analytical X-ray Instruments Inc.: Madison, WI, USA, 1998.
  33. Sheldrick, G.M. SADABS Program for Empirical Absorption Correction; University of Göttingen: Göttingen, Germany, 1996. [Google Scholar]
  34. Altomare, A.; Cascarano, G.; Giacovazzo, C.; Guagliardi, A.; Moliterni, A.G.G.; Burla, M.C.; Polidori, G.; Camalli, M.; Siliqi, D. SIRWARE: A universal data exchange format for crystallography. Acta Crystallogr. Sect. A 1996, 52, C79. [Google Scholar]
  35. Sheldrick, G.M. SHELXTL Plus (Windows NT Version) Structure Determination Package, Version 5.1; Bruker Analytical X-ray Instruments Inc.: Madison, WI, USA, 1998. [Google Scholar]
  36. Werner, H.-J.; Knowles, P.J.; Knizia, G.; Manby, F.R.; Schütz, M. Molpro: A general-purpose quantum chemistry program package. WIREs Comput. Mol. Sci. 2012, 2, 242–253. [Google Scholar] [CrossRef]
  37. Hay, P.J.; Wadt, W.R. Ab initio effective core potentials for molecular calculations. Potentials for the transition-metal atoms Sc to Hg. J. Chem. Phys. 1985, 82, 270–283. [Google Scholar] [CrossRef]
  38. Becke, A.D. Density–functional thermochemistry III: The role of exact exchange. J. Chem. Phys. 1993, 98, 5648–5652. [Google Scholar] [CrossRef] [Green Version]
  39. Weigend, F.; Ahlrichs, R. Balanced basis sets of split valence, triple zeta valence and quadruple zeta valence quality for H to Rn: Design and assessment of accuracy. Phys. Chem. Chem. Phys. 2005, 7, 3297–3305. [Google Scholar] [CrossRef] [PubMed]
  40. Samouei, H.; Grushin, V.V. New, Highly Efficient, Simple, Safe, and Scalable Synthesis of [1-(Carboxy-methyl)-thymine; (Ph3P)3Ru(CO)(H)2]. Organometallics 2013, 32, 4440–4443. [Google Scholar] [CrossRef]
  41. Lynam, J.M.; Welby, C.E.; Whitwood, A.C. Exploitation of a Chemically Non-innocent Acetate Ligand in the Synthesis and Reactivity of Ruthenium Vinylidene Complexes. Organometallics 2009, 28, 1320–1328. [Google Scholar] [CrossRef]
  42. Welby, C.E.; Eschemann, T.O.; Unsworth, C.A.; Smith, E.J.; Thatcher, R.J.; Whitwood, A.C.; Lynam, J.M. Ruthenium Acetate Complexes as Versatile Probes of Metal-Ligand Interactions: Insight into the Ligand Effects of Vinylidene, Carbene, Carbonyl, Nitrosyl and Isocyanide. Eur. J. Inorg. Chem. 2012, 2012, 1493–1506. [Google Scholar] [CrossRef]
  43. Akitt, J.W.; Mann, B.E. NMR and Chemistry. An Introduction to Modern NMR Spectroscopy, 4th ed.; Stamley Thornes: Cheltenham, UK, 2000. [Google Scholar]
  44. Colombo, M.; George, M.W.; Moore, J.N.; Pattison, D.I.; Perutz, R.N.; Virrels, I.G.; Ye, T.-Q. Ultrafast reductive elimination of Hydrogen from a metal carbonyl dihydride complex; A study by time-resolved and visible spectroscopy. J. Chem. Soc. Dalton Trans. 1997, 17, 2857–2859. [Google Scholar] [CrossRef]
  45. Procacci, B.; Duckett, S.B.; George, M.W.; Hanson-Heine, M.W.D.; Horwath, R.; Perutz, R.N.; Sun, X.-Z.; Vuong, K.Q.; Welch, J.A. Competing Pathways in the Photochemistry of Ru(H)2(CO)(PPh3)3. Organometallics 2018, 37, 865–868. [Google Scholar] [CrossRef] [Green Version]
  46. Miao-Chang, L.; Feng, W.; Ding, J.C. 1-(Carboxy-methyl)-hymine. Acta Crystallogr. 2004, 60, 1611. [Google Scholar]
  47. Hassanein, K.; Zamora, F.; Castillo, O.; Amo-Ochoa, P. Supramolecular interactions in Cobalt (II)–nucleobases complexes: A methyl matter. Inorg. Chim. Acta 2016, 452, 251–257. [Google Scholar] [CrossRef]
  48. Gawinecki, R.; Kuczek, A.; Kolehmainen, E.; Osmialowski, B.; Krygowski, T.M.; Kauppinen, R. Influence of bond fixation in Benzo-annulated N-salicyldeneanilines and their ortho-C(=O)X derivatives (X=CH3, NH2, OCH3) on tautomeric equilibria in solution. J. Org. Chem. 2007, 72, 4498–5607. [Google Scholar] [CrossRef] [PubMed]
  49. Hassanein, K.; Castillo, O.; Gómez-García, C.J.; Zamora, F.; Amo-Ochoa, P. Asymmetric and Symmetric Dicopper(II) Paddle-Wheel Units with Modified Nucleobases. Cryst. Growth Des. 2015, 15, 5485–5494. [Google Scholar] [CrossRef]
  50. Vendier, L.; Sabo-Etienne, S. Borane–Mediated Carbon Dioxide Reduction at Ruthenium: Formation of C1 and C2 Compounds. Angew. Chem. Int. Ed. 2012, 51, 1671–1674. [Google Scholar]
  51. Robinson, S.D.; Uttley, M.F. Complexes of the platinum metals. Part II. Carboxylato (triphenylphosphine) derivatives of ruthenium, osmium, rhodium, and iridium. Dalton Trans. 1973, 18, 1912–1920. [Google Scholar] [CrossRef]
  52. Emsley, J.; Freeman, N.J. β-diketone interaction. Part 5: Solvent effect on the keto-enol equilibrium. J. Mol. Struct. 1987, 161, 193–204. [Google Scholar] [CrossRef]
  53. Kawaguchi, S. Variety in the coordination modes of β-dicarbonyl compounds in metal complexes. Coord. Chem. Rev. 1986, 70, 51–84. [Google Scholar] [CrossRef]
  54. Correa, V.R.S.; Freire, M.; Barbosa, D.P.; Bezerra, L.; Bomfim, D.R.; Moreira, M.B.; Soares, J.A.; Ellena, J.; Batista, A.A. Ru(II)-thyminate complexes: New metallodrug candidates against tumor cells. New J. Chem. 2018, 42, 6794–6802. [Google Scholar] [CrossRef]
  55. Correa, R.S.; Bomfim, L.M.; Oliveira, K.M.; Moreirab, D.R.M.; Soares, M.B.P.; Ellena, J.; Bezerra, D.P.; Batista, A.A. Ru(II) complexes containing uracil nucleobase analogs with cytotoxicity against tumor cells. J. Inorg. Biochem. 2019, 198, 110751. [Google Scholar] [CrossRef] [PubMed]
  56. de Souza Oliveira, M.; Dantas de Santana, A.A.; Correa, R.S.; Botelho, M.; Soares, P.; Batista, A.A.; Bezerra, D.P. Ru(II)-Thymine Complex Causes Cell Growth Inhibition and Induction of Caspase-Mediated Apoptosis in Human Promyelocytic Leukemia HL-60 Cells. Int. J. Mol. Sci. 2018, 19, 1609–1622. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Peng, C.S.; Tokmakoff, A. Identification of Lactam-Lactim Tautomers of Aromatic Heterocycles in Aqueous Solution Using 2D IR Spectroscopy. J. Phys. Chem. Lett. 2013, 3, 3302–3306. [Google Scholar]
Scheme 1. The fluxional interconversion between trans-, mer-2 isomer and the related cis,fac-form envisages the formation of a transient penta-coordinated species, bearing double κ1-THAc coordination mode which affords the complete transformation to the cis-2 isomer.
Scheme 1. The fluxional interconversion between trans-, mer-2 isomer and the related cis,fac-form envisages the formation of a transient penta-coordinated species, bearing double κ1-THAc coordination mode which affords the complete transformation to the cis-2 isomer.
Applsci 11 03113 sch001
Figure 1. X-ray structure of κ1(O)-THAc, κ2(O,O)-THAc [Ru(CO)(PPh3)2] cis-2. Hydrogen atoms have been omitted for clarity.
Figure 1. X-ray structure of κ1(O)-THAc, κ2(O,O)-THAc [Ru(CO)(PPh3)2] cis-2. Hydrogen atoms have been omitted for clarity.
Applsci 11 03113 g001
Figure 2. Distinct H-bond interactions have been evidenced by different color notifications. Applsci 11 03113 i001 Expected intermolecular Watson-Crick H-interactions Å NH-OC between two adjacent THAc rings. Applsci 11 03113 i002 OC(O)-H (C5H5-P) between the κ1-THAc C(O)O and the proton belonging to a phenyl of PPh3 of the adjacent molecule O6-H45 (2.51) O-H31 (2.44) N2H-O7 (2.12) N4H-O4 (2.05). Applsci 11 03113 i003 π–π stacking between two aromatic rings (3.6 Å) of THAc ring belonging to opposite molecules.
Figure 2. Distinct H-bond interactions have been evidenced by different color notifications. Applsci 11 03113 i001 Expected intermolecular Watson-Crick H-interactions Å NH-OC between two adjacent THAc rings. Applsci 11 03113 i002 OC(O)-H (C5H5-P) between the κ1-THAc C(O)O and the proton belonging to a phenyl of PPh3 of the adjacent molecule O6-H45 (2.51) O-H31 (2.44) N2H-O7 (2.12) N4H-O4 (2.05). Applsci 11 03113 i003 π–π stacking between two aromatic rings (3.6 Å) of THAc ring belonging to opposite molecules.
Applsci 11 03113 g002
Figure 3. 1H NMR (nuclear magnetic resonance) spectrum of the two 3a, 3bmer- κ1(O)THAc-[Ru(CO)H (PPh3)3] rotamers evidences two doublets of triplets (dt, −6.12 ppm; −7.27ppm, in red) and a triplet for the chelate trans(P,P)2(O,O)THAc-[Ru(CO)H(PPh3)2] 4 (t, −17.01 ppm in blue).
Figure 3. 1H NMR (nuclear magnetic resonance) spectrum of the two 3a, 3bmer- κ1(O)THAc-[Ru(CO)H (PPh3)3] rotamers evidences two doublets of triplets (dt, −6.12 ppm; −7.27ppm, in red) and a triplet for the chelate trans(P,P)2(O,O)THAc-[Ru(CO)H(PPh3)2] 4 (t, −17.01 ppm in blue).
Applsci 11 03113 g003
Scheme 2. Reaction key steps sequence to transform 1 into 2. Monohapto intermediates 3a, 3b have generated by restricted C-O rotation of the THAc fragment about the phosphine ligands.
Scheme 2. Reaction key steps sequence to transform 1 into 2. Monohapto intermediates 3a, 3b have generated by restricted C-O rotation of the THAc fragment about the phosphine ligands.
Applsci 11 03113 sch002
Scheme 3. Calculations on torsional barriers. κ1-THAc rotation in 3 takes place along the C-O bond of the acetate group, whereas κ2-THAc rotation in 4 is along C-C bond of the THAc-side arm. Only half of the torsion has displayed for clarity.
Scheme 3. Calculations on torsional barriers. κ1-THAc rotation in 3 takes place along the C-O bond of the acetate group, whereas κ2-THAc rotation in 4 is along C-C bond of the THAc-side arm. Only half of the torsion has displayed for clarity.
Applsci 11 03113 sch003
Scheme 4. Calculated DFT free energies (relative to the reaction of 1 with one equivalent of THAcOH (blue dash line) for the intermediates shown in Scheme 2 with other transient species intercepted spectroscopically only. It is noteworthy that the trans conformations are the more stable in all the cases, but the related cis forms are notably stabilized by intermolecular H-bonding interactions.
Scheme 4. Calculated DFT free energies (relative to the reaction of 1 with one equivalent of THAcOH (blue dash line) for the intermediates shown in Scheme 2 with other transient species intercepted spectroscopically only. It is noteworthy that the trans conformations are the more stable in all the cases, but the related cis forms are notably stabilized by intermolecular H-bonding interactions.
Applsci 11 03113 sch004
Scheme 5. Subsequent addition reactions of acetic acid (AcOH) to 1. The reactions have been performed by adopting different conditions of temperature, time duration, reactant ratio and solvent polarity. All the optimized tests have been reported in the experimental section.
Scheme 5. Subsequent addition reactions of acetic acid (AcOH) to 1. The reactions have been performed by adopting different conditions of temperature, time duration, reactant ratio and solvent polarity. All the optimized tests have been reported in the experimental section.
Applsci 11 03113 sch005
Scheme 6. Calculated DFT free energies (relative to 1 + 2 AcH, blue dashed line) for the reaction intermediates shown above and H-bonded postulated isomers of similar energy.
Scheme 6. Calculated DFT free energies (relative to 1 + 2 AcH, blue dashed line) for the reaction intermediates shown above and H-bonded postulated isomers of similar energy.
Applsci 11 03113 sch006
Scheme 7. Calculated free energies (kJmol−1) for the equilibrium between THAcH (AD) and anionic THAcK (EG) amido and imino-enol tautomers.
Scheme 7. Calculated free energies (kJmol−1) for the equilibrium between THAcH (AD) and anionic THAcK (EG) amido and imino-enol tautomers.
Applsci 11 03113 sch007
Figure 4. Steric comparison of the hindrance between THAcK and κ2-(O,O)[RuH(CO)(THAc)(PPh3)2].
Figure 4. Steric comparison of the hindrance between THAcK and κ2-(O,O)[RuH(CO)(THAc)(PPh3)2].
Applsci 11 03113 g004
Scheme 8. Relative free energies (kJmol−1) of THAc-K intermolecularly bound dimeric forms. The relative energies have been evaluated has half of the difference between the energy of the dimer and twice the energy of keto monomer (Structure A, Scheme 7). The O···H-N interactions are displayed in green, N-H···N in blue and O-H···O in purple and considered to have similar energy to generate coincident H-bonds in the 1H NMR spectrum. The designation of carboxy-anti has been referred to the reciprocal disposition of the C(O)O functions in the dimer. The same color notification is reported in the 1H NMR spectrum (12–7 ppm interval) (Figure 5).
Scheme 8. Relative free energies (kJmol−1) of THAc-K intermolecularly bound dimeric forms. The relative energies have been evaluated has half of the difference between the energy of the dimer and twice the energy of keto monomer (Structure A, Scheme 7). The O···H-N interactions are displayed in green, N-H···N in blue and O-H···O in purple and considered to have similar energy to generate coincident H-bonds in the 1H NMR spectrum. The designation of carboxy-anti has been referred to the reciprocal disposition of the C(O)O functions in the dimer. The same color notification is reported in the 1H NMR spectrum (12–7 ppm interval) (Figure 5).
Applsci 11 03113 sch008
Figure 5. 1H spectra relative to (3 + 4) mixture showing the intermediate trapping during the reaction course. The integrals evaluation for the tautomeric ratio results to be 1:3.2 with the prevalence of the stable amido conformation. The attributions of the H-bond functions of the type X–H–O (X=N, O) have been assigned based on the shielding effect (δ) O–H–O (11) < N–H–O (10-8) < N–H–N (8), related to the calculated DFT energy scale of the dimers sketched in Scheme 8.
Figure 5. 1H spectra relative to (3 + 4) mixture showing the intermediate trapping during the reaction course. The integrals evaluation for the tautomeric ratio results to be 1:3.2 with the prevalence of the stable amido conformation. The attributions of the H-bond functions of the type X–H–O (X=N, O) have been assigned based on the shielding effect (δ) O–H–O (11) < N–H–O (10-8) < N–H–N (8), related to the calculated DFT energy scale of the dimers sketched in Scheme 8.
Applsci 11 03113 g005
Table 1. Selected conditions for the reaction of 1 with AcOH.
Table 1. Selected conditions for the reaction of 1 with AcOH.
Solvent, rfxRu/AcH RatioTime (min)SpeciesYield (%)
CHCl31:2240765
CDCl3 r.t.1:114 h (r.t);5 + 650 (6) + 38 (5)
CDCl31:230750 (7)
CH2Cl21:240570
toluene1:12405 + 670 (5) + 30 (6)
CPME1:240770 (7)
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Bordoni, S.; Cerini, S.; Tarroni, R.; Monari, M.; Micheletti, G.; Boga, C. Ruthenium–Thymine Acetate Binding Modes: Experimental and Theoretical Studies. Appl. Sci. 2021, 11, 3113. https://doi.org/10.3390/app11073113

AMA Style

Bordoni S, Cerini S, Tarroni R, Monari M, Micheletti G, Boga C. Ruthenium–Thymine Acetate Binding Modes: Experimental and Theoretical Studies. Applied Sciences. 2021; 11(7):3113. https://doi.org/10.3390/app11073113

Chicago/Turabian Style

Bordoni, Silvia, Stefano Cerini, Riccardo Tarroni, Magda Monari, Gabriele Micheletti, and Carla Boga. 2021. "Ruthenium–Thymine Acetate Binding Modes: Experimental and Theoretical Studies" Applied Sciences 11, no. 7: 3113. https://doi.org/10.3390/app11073113

APA Style

Bordoni, S., Cerini, S., Tarroni, R., Monari, M., Micheletti, G., & Boga, C. (2021). Ruthenium–Thymine Acetate Binding Modes: Experimental and Theoretical Studies. Applied Sciences, 11(7), 3113. https://doi.org/10.3390/app11073113

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop