Biospeciation of Potential Vanadium Drugs of Acetylacetonate in the Presence of Proteins
Giuseppe Sciortino1,2 
Valeria Ugone1 
Daniele Sanna3 
Giuseppe Lubinu1 
Simone Ruggiu1 
Jean-Didier Maréchal2 
Eugenio Garribba1*- 1Dipartimento di Chimica e Farmacia, Università di Sassari, Sassari, Italy
- 2Departament de Química, Universitat Autònoma de Barcelona, Cerdanyola del Vallés, Barcelona, Spain
- 3Istituto di Chimica Biomolecolare, Consiglio Nazionale delle Ricerche, Sassari, Italy
Among vanadium compounds with potential medicinal applications, [VIVO(acac)2] is one of the most promising for its antidiabetic and anticancer activity. In the organism, however, interconversion of the oxidation state to +III and +V and binding to proteins are possible. In this report, the transformation of VIII(acac)3, VIVO(acac)2, and VVO2(acac) after the interaction with two model proteins, lysozyme (Lyz) and ubiquitin (Ub), was studied with ESI-MS (ElectroSpray Ionization-Mass Spectroscopy), EPR (Electron Paramagnetic Resonance), and computational (docking) techniques. It was shown that, in the metal concentration range close to that found in the organism (15–250 μM), VIII(acac)3 is oxidized to VIVO(acac)+ and VIVO(acac)2, which—in their turn—interact with proteins to give n[VIVO(acac)]–Protein and n[VIVO(acac)2]–Protein adducts. Similarly, the complex in the +IV oxidation state, VIVO(acac)2, dissociates to the mono-chelated species VIVO(acac)+ which binds to Lyz and Ub. Finally, VVO2(acac) undergoes complete dissociation to give the 'bare' VVO ion that forms adducts n[VVO2]–Protein with n = 1–3. Docking calculations allowed the prediction of the residues involved in the metal binding. The results suggest that only the VIVO complex of acetylacetonate survives in the presence of proteins and that its adducts could be the species responsible of the observed pharmacological activity, suggesting that in these systems VIVO2+ ion should be used in the design of potential vanadium drugs. If VIII or VVO2 potential active complexes had to be designed, the features of the organic ligand must be adequately modulated to obtain species with high redox and thermodynamic stability to prevent oxidation and dissociation.
Introduction
Among the metal-based drugs, vanadium (V) compounds have drawn a growing interest for their pharmacological properties shown both in vivo and in vitro studies. In particular, several V complexes have been proposed as antidiabetic and anticancer potential drugs (Sakurai et al., 2008; Costa Pessoa et al., 2015a; Kioseoglou et al., 2015; Rehder, 2016) and some of them reached the clinical trials, for example bis(ethylmaltolato)oxidovanadium(IV) or BEOV (Thompson and Orvig, 2006; Thompson et al., 2009). Over the last years, many complexes were studied with vanadium in one of the three oxidation states found in biological systems, +III, +IV, and +V. However, the mechanism of action and active oxidation state in the organism remain unclear and biospeciation has not been completely ascertained. This lack of knowledge is one of the reasons why V compounds are not considered for clinical tests by pharmaceutical companies and limits the possibility of a rational design of new potential drugs with more activity and lower toxicity than those proposed up to now (Scior et al., 2016).
Among these compounds, bis(acetylacetonato)-oxidovanadium(IV) or VIVO(acac)2 is one of the most studied (Makinen and Brady, 2002; Crans et al., 2018). It favors the decrease of blood glucose concentration and is more effective than the inorganic salt VIVOSO4 in the suppression of hepatic glycolysis in streptozotocin-induced diabetic rats, and this effect cannot be related only to an improved intestinal absorption (Reul et al., 1999; Amin et al., 2000). VIVO(acac)2 enhances in cells the IR (insulin receptor) kinase activity more than other VIVO chelate species (Makinen and Salehitazangi, 2014), even though it has not been proven if its action is direct through the stimulation of the IR kinase activity or indirect through the activation of other tyrosine kinases or inhibition of tyrosine phosphatases, for example PTP-1B (Ou et al., 2005; Makinen et al., 2007). VIVO(acac)2 exhibits also important anti-proliferative effects: it blocks the G1/S phase of the cell cycle progression via an activated ERK (extracellular signal-regulated kinase) signal in human hepatoma HepG2 cells (Fu et al., 2008), while in human pancreatic cancer cell line AsPC-1 it induces cell cycle arrest at the G2/M phase and causes an increase of ROS (reactive oxygen species) concentration (Wu et al., 2016). These results indicate that V-based complexes may be potentially applied as antidiabetic and anticancer agents to treat patients suffering from both the diseases (Liu et al., 2013).
Species related to VIVO(acac)2 in the organism are VIII(acac)3 and VVO2(acac). In fact, on one hand several experimental evidences showed that VIV complexes can be rapidly oxidized to VV in the cell culture medium and, on the other, that +IV state can converge to +III in the gastrointestinal medium, in the blood and reducing cellular environment; for these reasons, more than one species could contribute to the global biological action. The biotransformation of VIVO(acac)2 is obviously related to the pharmacological efficacy, since it determines the gastrointestinal absorption, the transport in blood and—overall—the availability of the active species. The interaction with proteins plays a key role in the biospeciation of a potential metal-based drug in the organism due to their high concentration and/or high affinity toward the metals (Costa Pessoa et al., 2015b).
In this paper, the transformation of VIII(acac)3, VIVO(acac)2, and VVO2(acac) in the presence of two model proteins, lysozyme (Lyz) and ubiquitin (Ub), was studied by ESI-MS (ElectroSpray Ionization-Mass Spectrometry), EPR (Electron Paramagnetic Resonance), and computational (docking and QM) techniques. Both ESI-MS and EPR techniques present, for what concerns this work, pros and cons. The main advantage of ESI-MS is that the metal concentrations used for experiments are in the range 1–100 μM, i.e., that found for V at physiological conditions (Jakusch and Kiss, 2017); the limitation is that it does not give any information on the structure of the revealed vanadium adducts and that the intensity of the peaks cannot be related to the amount of the species in solution. With regard to EPR spectroscopy, it allows only the recognition of the type of amino acid donors involved in the metal coordination (His-N, Asp/Glu-COO, or carbonyl-CO, etc.), with the limitation that the specific site and its three-dimensional structure cannot be identified (Chasteen, 1981; Smith et al., 2002).
The results of this study on the transformation of VIVO(acac)2 and—in general—of bis-chelated complexes with stoichiometry VIVOL2, where L is a bidentate monoanionic organic ligand, could be very useful to ascertain the active species in the organism and be used as a guide for a rational design of new potential V drugs.
Materials and Methods
Chemicals
Water was deionized through the purification system Millipore MilliQ Academic or purchased from Sigma-Aldrich (LC-MS grade). VIII(acac)3 and VIVO(acac)2, NH4VVO3, acetylacetone (Hacac), and 4-(2-hydroxyethyl)piperazine-1-ethanesulfonic acid (HEPES) were Sigma-Aldrich products of the highest grade available and used as received.
Chicken egg white lysozyme (Lyz; molecular mass of 14.3 kDa, code 62970) and ubiquitin from bovine erythrocytes (Ub; molecular mass of 8.6 kDa, code U6253) were purchased from Sigma-Aldrich.
ESI-MS Measurements
The solutions for ESI-MS measurements were prepared as follows. Briefly, the solid complexes VIII(acac)3 and VIVO(acac)2 or NH4VVO3 plus Hacac in molar ratio 1/2 were dissolved in LC-MS grade water to obtain a V concentration of 1.0–2.0 mM. Subsequently, all the solutions were diluted in LC-MS grade water with an aliquot of a stock protein solution (500 μM) to reach (V complex)/Protein ratios of 3/1, 5/1, and 10/1 and a protein concentration of 5 μM. In all the solutions pH was in the range 6–7. Immediately after the preparation of the solutions, ESI-MS spectra were acquired.
Mass spectra in the positive-ion or negative-ion mode were recorded on a Q Exactive™ Plus Hybrid Quadrupole-Orbitrap™ (Thermo Fisher Scientific) mass spectrometer. The solutions were infused at a flow rate of 5.00 μL/min into the ESI chamber. The spectra were obtained in the m/z range 300–4,500 at a resolution of 140,000 and accumulated for at least 5 min to increase the signal-to-noise ratio. The instrumental setting for the measurement of the spectra in positive-ion mode were: spray voltage 2,300 V, capillary temperature 250°C, sheath gas 10 (arbitrary units), auxiliary gas 3 (arbitrary units), sweep gas 0 (arbitrary units), probe heater temperature 50°C. The setting used for negative-ion mode spectra were: spray voltage −1,900 V, capillary temperature 250°C, sheath gas 20 (arbitrary units), auxiliary gas 5 (arbitrary units), sweep gas 0 (arbitrary units), probe heater temperature 14°C. ESI-MS spectra were analyzed by using Thermo Xcalibur 3.0.63 software (Thermo Fisher Scientific) and the average deconvoluted monoisotopic masses were obtained through the Xtract tool integrated in the software.
EPR Measurements
EPR spectra were acquired by dissolving in ultra-pure water VIII(acac)3 to obtain a VIII concentration of 1 mM. To the solutions, HEPES buffer of 0.1 M concentration was added. The value of pH was brought in the pH range 5–7 and a protein (Lyz or Ub) was added to 1 mL of the solution with the V species, to reach a concentration of 0.5 mM and a molar ratio (V complex)/Protein of 2/1. Subsequently, an aliquot of the solution with V complex was added to bring the ratio (V complex)/Protein to 3/1.
EPR spectra were recorded at 120 K with an X-band Bruker EMX spectrometer equipped with a HP 53150A microwave frequency counter using a microwave frequency of 9.40–9.41 GHz, microwave power of 20 mW, time constant of 81.92 ms, modulation frequency of 100 kHz, modulation amplitude of 0.4 mT, and resolution of 4,096 points. The spectra were immediately recorded after the preparation of the samples. To increase the signal to noise ratio, signal averaging was employed (Sanna et al., 2009).
DFT and Docking Calculations
All the DFT calculations were performed with Gaussian 09 (revision D.01) (Frisch et al., 2010). The geometries and harmonic frequencies of VIVO and VVO2 complexes were computed at the level of theory B3P86/6-311++g(d,p) with the SMD model (Marenich et al., 2009) for water.
Docking calculations were carried out with GOLD 5.8 software (Jones et al., 1997) on the X-ray structure available in the Protein Data Bank (PDB) of lysozyme [PDB code: 2LYZ (Diamond, 1974)] and ubiquitin [3H1U (Qureshi et al., 2009)]. The PDB structure was cleaned removing all the small molecules and crystallographic waters, and hydrogen atoms were added with the UCSF Chimera (Pettersen et al., 2004).
The docking simulations were carried out constructing in the region of interest an evaluation sphere of 12 Å. The coordination docking calculations (with an active binding V–donor) were performed on the DFT optimized moieties VIVO(acac)+ and cis-VVO, activating one, two or three equatorial coordination positions, respectively, with dummy hydrogen atoms (Sciortino et al., 2017, 2018a,b, 2019a). All the possible moieties VVO2(H2O)n+ with n = 0–2 were taken into account during dockings and the proposed solutions in the text are those with higher scoring. For coordination dockings the regions with the potential coordinating residues were examined with the GOLD rotamers libraries (Lovell et al., 2000). For the non-coordination dockings (inert binding V···donor) with VIVO(acac)2, the whole rigid protein was considered. The solutions were analyzed by means of GaudiView (Rodríguez-Guerra Pedregal et al., 2017).
To identify potential binding sites, the protein space was first probed with multisite.py for zones where at least two potential coordinating residues (Asp, Glu, Asn, Gln, and His) featured an α- and a β-carbon within a range of distances from each queried grid point of 2.6–6.4 and 3.4–7.2 Å, respectively (Sciortino et al., 2019b). Both the protonation states at δ and ε nitrogens of His imidazole ring were taken into account during the simulations.
Genetic algorithm (GA) parameters have been set to 50 GA runs and a minimum of 100,000 operations. The other parameters of GA were set to default.
The scoring (Fitness of GoldScore) was calculated with the modified version of GoldScore scoring function, which was validated in previously published papers (Sciortino et al., 2017, 2018a, 2019a). The best solutions (binding poses) were evaluated through the mean (Fmean) and the highest value (Fmax) of the Fitness associated with each pose, the population of the cluster containing the best pose and the position in the Fitness ranking of the computed cluster.
Results and Discussion
Interaction of VIII(acac)3 and VIVO(acac)2 With Proteins
The species distribution diagram of the VIII/acac system with molar ratio 1/3 and V concentration of 150 μM, that does not account for the oxidation process to VIV, can be built using the data of the thermodynamic stability constants reported in the literature (Brito et al., 2009). The diagram is shown in Figure S1 of the Supplementary Material and suggests that the main species in the pH range 6.0–7.0 is VIII(acac)3. ESI-MS spectrum was recorded on VIII(acac)3 at pH 6.2 in a degassed aqueous solution to evaluate its stability in the absence of a protein (Figure S2). The results indicate that the oxidation to VIVO is only partial: even though the intensity of the MS signal cannot be directly related to the concentration in solution, the peaks of the adducts [VIII(acac)3+H+] and [VIII(acac)3+Na+] (m/z = 349.09 and 371.07, respectively) are much more intense than those assigned to [VIVO(acac)2+H+] and [VIVO(acac)2+Na+] (m/z = 266.04 and 288.02). The simulation of the two peaks of VIII(acac)3 is shown in Figures S3, S4 and confirms the attribution.
When the spectra were measured on the system VIII(acac)3/Lyz the transformation to VIVO is complete and no peaks assignable to the active or inert binding of VIII(acac)3 to the protein could be revealed. With a molar ratio of 3/1 and a protein concentration of 5 μM, the peaks of the free protein at 14304 Da plus the signals of the adducts n[VIVO(acac)]–Lyz with n = 1–4 (with mass of 14469, 14635, 14800 and 14965 Da, respectively) were observed (Figure S5). For these species the equatorial interaction with one or two protein residues is possible. Similar results were obtained at (V complex)/Protein ratio of 5. When higher protein concentration was used (50 μM), the spectra show a series of peaks belonging to n[VIVO(acac)]–Lyz, with n = 1–4, while the signals at 14570 and 14836 Da can be attributed to [VIVO(acac)2]–Lyz and 2[VIVO(acac)2]–Lyz (Figure 1A). Moreover, the adduct [VIVO(acac)+VIVO(acac)2]–Lyz at 14735 Da is detected with the contemporaneous coordination of mono- and bis-chelated moieties (Figure 1A). Therefore, the increase of the V concentration favors the interaction of the bis-chelated complex VIVO(acac)2 with lysozyme, even though the binding of the fragment VO(acac)+ appears to be favored compared to VIVO(acac)2. Concerning the interaction between [VO(acac)2] with lysozyme, no equatorial sites are accessible and so only two possibilities exist: an inert binding with the protein surface or a weak axial active binding through an amino acid residue.
Figure 1. Deconvoluted ESI-MS spectrum recorded on the systems: (A) VIII(acac)3 and lysozyme, pH 6.2, and (B) VIII(acac)3 and ubiquitin, pH 6.3. The molar ratio VIII/Protein was 3/1 and Protein concentration 50 μM. Under these experimental conditions, VIII species fully oxidizes to VIVO.
In the analysis of the results on the system VIVO(acac)2/Lyz, it must be considered that in the examined pH range VIVO(acac)2 is the major species in aqueous solution and that it coexists with 1:1 complex VIVO(acac)+. The diagram with the distribution of the species is represented in Figure 2. The ESI-MS spectra are in agreement with these data and show similar results with formation of the adducts n[VIVO(acac)]–Lyz and n[VIVO(acac)2]–Lyz (Figure S6).
Figure 2. Distribution diagram of the species as a function of pH in a solution containing VIVO2+ ion and acac (VIVO/acac = 1/2, VIVO = 1.50 ×10−4 M). The stability constants of the VIVO–acac complexes were taken from Crans et al. (2001) and of the hydrolytic species from Buglyo et al. (2000).
EPR spectroscopy confirms the oxidation of VIII and the resonances due to the VIVO ion (3d1 electronic configuration) are revealed (Figure 3A). Beside the absorption due to VIVO(acac)2 not covalently bound to lysozyme, the signals were attributed to [VIVO(acac)]–Lyz adducts with the equatorial binding of two protein residues (gz ~ 1.943 and Az ~ 171–172 × 10−4 cm−1). This is illustrated in Figures S7A,B.
Figure 3. EPR spectra recorded on frozen solutions (120 K) of the systems: (A) VIII(acac)3 and lysozyme and (B) VIII(acac)3 and ubiquitin. VIII/Protein molar ratio was 3/1 and V concentration 1 mM. pH was 5.70 and 5.15, respectively. The region at high field was amplified 10 times. At these conditions, VIII species fully oxidizes to VIVO.
The deconvoluted ESI-MS spectrum of the system VIII(acac)3/Ub 3/1 using an Ub concentration of 5 or 50 μM (Figure 1B) showed the presence of [VIVO(acac)]–Ub, 2[VIVO(acac)]–Ub and 3[VIVO(acac)]–Ub. This means that three moieties VIVO(acac)+ bind to ubiquitin, which replaces two weak water ligands in the two adjacent equatorial positions of VIVO(acac)(H2O). When the interaction between VIVO(acac)2 and ubiquitin is considered, the results are similar to those revealed with VIII and recently published, that indicate the formation of the adducts n[VIVO(acac)]–Ub, with n = 1–3 (Ugone et al., 2019). For the system containing Ub as well, EPR measurements indicated the complete transformation of VIII(acac)3 to VIVO(acac)+ moiety, which in solution interacts with ubiquitin (Figure 3B).
Overall, the behavior of the VIII and VIVO systems with the two model proteins is very similar. The results indicate that, in the metal concentration range 15–250 μM, i.e., that found in the organism (Jakusch and Kiss, 2017), VIII(acac)3 is quantitatively oxidized to VIVO(acac)+ and VIVO(acac)2, which—in their turn—interact with proteins. Compared to the binary systems without proteins, the oxidation is favored probably due to the stabilization of the VIVO adducts upon the binding with the amino acid residues. So, the results obtained with VIII can be interpreted considering the species distribution diagram shown in Figure 2 for VIVO–acac system, postulating that the oxidation process in the presence of the proteins is quantitative. The reactions which account for the formation of the adducts [VIVO(acac)+]–Protein and [VIVO(acac)2]–Protein are as follows:
The number of VIVO(acac)+ and VIVO(acac)2 moieties bound depends on the vanadium concentration and type of protein; in particular, the binding of VIVO(acac)2 is favored with increasing V concentration and going from Ub to Lyz.
Docking calculations were performed to confirm ESI-MS and EPR findings and get information on the residues involved in the VIV binding. For lysozyme, the results indicate that four-five sites are available for VIVO(acac)+ interaction, in agreement with ESI-MS measurements (Table 1 and Figure 4). In the region 44–52 one or two VIVO(acac)+ can bind to protein with the simultaneous coordination of two residues. The affinity order is as follows: (Asn46, Asp48) > (Asn46, Asp52) ~ (Asn44, Asp52), the Fmean values being in the range 35.9–40.8 with a population between 41/50 and 50/50. In all the three sites, named A, the donor set is (NCO, COO−) and this accounts for the Az value measured in the EPR spectra (Az ~ 171–172 × 10−4 cm−1, see above). It must be observed that the contemporaneous coordination to both (Asn46, Asp52) and (Asn44, Asp52) donor set is not compatible and the binding to the first site excludes the second one. The second site (site B) is located in the N-terminal region, where His15 plus a water ligand are coordinated to the VIVO2+ ion, and a hydrogen bond between V=O group and Thr89 stabilizes the binding. Interestingly, the His15 equatorial binding is not expected for bis-chelated cis moieties such as cis-VIVO(picolinato)2 or cis-VIVO(maltolato)2 due to the steric hindrance of the second ligand with an (equatorial-axial) arrangement (Sciortino et al., 2017). The third and fourth sites are weaker than these ones and are based on the coordination of (Asp18, Asn19) with Fmean = 30.9 (site C) and (Asp119, Gln121) with Fmean = 30.7 (site D). The four identified sites are represented in Figure 4. For the interaction of VIVO(acac)2, only inert binding must be considered: the best site is based on a hydrogen bond stabilization between the VIV=O and NH group of Trp63 residue and a weak axial interaction of NCO of Ala107 with vanadium. The two solutions are displayed in Figures 5A,B, showing that with increasing the strength of the axial interaction with Ala107 decreases the strength of the contact of V=O with Trp63 and viceversa; in Figure 5A V=O···HN–Trp63 is 2.029 Å and V···OCN–Ala107 is 3.370 Å, while in Figure 5B V=O···HN–Trp63 is 1.644 Å and V···OCN–Ala107 is 3.875 Å.
Table 1. Docking solutions for the active binding of the VIVO(acac)+ moiety with Lyz and Ub.
Figure 4. Most stable structures obtained by docking methods for the active binding of VIVO(acac)+ to Lyz: (A) site A with (Asn44, Asp52) and (Asn46, Asp52); (B) site B with His15; (C) site C with (Asp18, Asn19); and (D) site D with (Asp119, Gln121). The hydrogen bonds are indicated with the full blue lines.
Figure 5. Most stable site obtained by docking methods for the inert binding of VIVO(acac)2 to Lyz: (A) V=O···HN–Trp63 is 2.029 Å and V···OCN–Ala107 3.370 Å and (B) V=O···HN–Trp63 is 1.644 Å and V···OCN–Ala107 3.875 Å. The weak bonds with Trp63 and Ala107 are indicated with the full blue lines.
When the interaction of VIVO species and ubiquitin is examined, the data indicated that two main sites exist. The first is based on the coordination of (Glu 16, Glu18; site 1) or (Glu18, Asp21; site 1') (Ugone et al., 2019); the two sites are mutually exclusive and have Fmean in the range 39.3–44.3 with population between 14 and 20%. The second site is due to the equatorial binding of the couple (Glu24, Asp52; site 2) or (Glu51, Asp52; site 2') (Ugone et al., 2019); in this case too, these sites are not independent with Fmean = 34.9–37.2 and population = 12–88%. The two sites are shown in Figure 6 and the values of Fmax, Fmean and population in Table 1.
Figure 6. Most stable structures obtained by docking methods for the active binding of VIVO(acac)+ to Ub: (A) site 1 with (Glu16, Glu18) and (B) site 2 with (Glu24, Asp52). The hydrogen bonds are indicated with the full blue lines.
Interaction of VVO2(acac) With Proteins
Positive-ion mode and negative-ion mode ESI-MS spectra were recorded on the system H2VVO/Hacac 1/2 with V concentration of 150 μM at pH 7.0 (Figure S8). In the literature, the thermodynamic stability constants on the system VV/acac are lacking. The data for maltol, a (O,O) ligand with similar basicity [pKa is 8.44 for maltol (Buglyo et al., 2002) and 8.76 for acetylacetone (Crans et al., 2001)], indicate that at these experimental conditions, the percentage of V in solution as VVO2(maltolato) in the pH range 6–7 is <30% (Elvingson et al., 1996); the amount of VVO2(maltolato) increases significantly with V concentration and overcomes 70% when it is ten times larger, i.e., 1.5 mM. This demonstrates that the hydrolytic processes cannot be neglected with decreasing V concentration. A similar behavior is expected for VV-acac system: specifically, the complex VVO2(acac) is formed only at high V concentration (around some mM), but it is not stable when the metal concentration is in the order of μM.
ESI-MS results indicated the presence of [Hacac+H+] and [Hacac+Na+] in the positive-ion mode (Figure S8A), and of [H2VVO] in the negative-ion mode (Figure S8B), derived from the hydrolysis of VVO2(acac). In the deconvoluted spectrum recorded in the system containing VVO2(acac) and lysozyme (Figure 7A), the signal corresponding to the mass of Lyz at 14304 Da and the series of signals of the adducts between the protein and one, two or three VVO ions are revealed. In all these adducts the acetylacetonato ligand is not coordinated to vanadium(V), confirming that the hydrolysis causes the removal of acac− from the first coordination sphere of the metal; interestingly, the free VVO ion interacts with lysozyme which behaves as a polydentate ligand. Moreover, the adducts with one or two water ligands were revealed (peaks indicated by the asterisks in Figure 7).
Figure 7. Deconvoluted ESI-MS spectrum recorded on the systems containing: (A) VVO2(acac) and lysozyme, pH 6.4, and (B) VVO2(acac) and ubiquitin, pH 6.4. The molar ratio VIII/Protein was 3/1 and Protein concentration 50 μM. With the asterisks the adducts with the binding of one and two water molecules are indicated.
The behavior of the system with ubiquitin is similar and the adducts n[VVO2]–Ub are detected with n = 1–3 (Figure 7B). This means that three sites of Ub are available for the metal binding.
The docking results confirmed the suggestions of ESI-MS data and indicate that three sites are available both for Lyz and Ub. One, two or three donors are coordinated to V, which shows a slightly distorted trigonal bipyramidal environment (Figure 8). With lysozyme the binding occurs at site A (Asn46, Asp52, Asn59), B (His15, Asp87), C (Asp18, Asn 19) with different affinity and Fmax of 48.5, 41.0, and 22.4 (Table 2). With ubiquitin the coordination to sites 1 (Glu16, Glu18) and 2′ (Glu24, Asp52) plus to His68 with Fmax in the range 22.3–45.2 is predicted (Table 2). The monodentate coordination of Ub at His68 is stabilized by a V=O···H3N–Lys6 hydrogen bond and reminds that of His12 in acid phosphatase (Lindqvist et al., 1994).
Figure 8. Most stable structures obtained by docking methods for the active binding of VVO to Lyz (A–C) and Ub (D–F). The hydrogen bonds are highlighted with the full blue lines.
Table 2. Docking solutions for the active binding of the VVO ion with Lyz and Ub.
Conclusions
Even though over the recent years many studies in bioinorganic and medicinal inorganic chemistry were dedicated to the development of vanadium-based potential drugs, the active oxidation state in the organism and mechanism of action are still elusive and their biospeciation has not completely explained. This lack of knowledge is one of the reasons why the number of V complexes under clinical tests is significantly lower than those of other metal ions. For example, while the oxidation state of Pt remains stable under physiological conditions (+II), for V three oxidation states are possible (+III, +IV, and +V) with an interconversion among them in the body. Therefore, in the design of vanadium drugs, such biotransformation and active state in the organism must be taken into account.
In this study, the V complexes with acetylacetonate, which displayed very promising antidiabetic and anticancer activity, were examined. It may be considered just as an example for the behavior of other V compounds. It was shown that, in the metal concentration range close to that found in the organism (15–250 μM), the species with +III state, VIII(acac)3, is completely oxidized to VIVO(acac)+ and VIVO(acac)2 and this process is favored—with respect to the binary system without protein—by the binding of one or more amino acid side-chains. VIVO(acac)+ and VIVO(acac)2, in their turn, interact with proteins: up to four VIVO(acac)+ can bind to proteins to yield n[VIVO(acac)]–Ub/Lyz adducts, the binding of VIVO(acac)+ being favored with decreasing the vanadium concentration. At high concentration, also n[VIVO(acac)2]–Protein adduct may be revealed with an axial active or a surface inert interaction with the protein. The results with the complex in the +IV oxidation state are the same, with VIVO(acac)2 undergoing dissociation to the mono-chelated species VIVO(acac)+ moiety which binds to the proteins. Finally, VVO2(acac) undergoes complete dissociation to give the “bare” VVO ion. Docking calculations allowed us to predict the residues involved in the metal binding and three-dimensional structure of the formed adducts. Three important sites were predicted both for lysozyme and ubiquitin, A (with Asn44, Asn46, Asp48, Asp52, and Asn59), B (His15), and C (Asp18, Asn 19) for Lyz, and 1/1' (Glu16, Glu18, Asp21), 2/2' (Glu24, Glu51, Asp52), and His68 for Ub.
The results suggest that, at the studied conditions, only VIVO species of acetylacetonate survive in the presence of proteins and these could be the species responsible of the observed pharmacological activity. Starting from VIII or VV may be useless, since VIII gives quantitatively VIVO, while VV undergoes hydrolysis to the inorganic ions and the active species should be the same as when inorganic vanadium(V) was administered. This finding allow us to suggest that in this system and in those containing a bidentate organic ligand, VIVO ion should be used in the design of potential vanadium drugs, and that, if VIII or VVO2 potential active complexes had to be designed, species with high redox and thermodynamic stability should be synthesized to prevent oxidation and dissociation processes. To reach this aim, it is necessary to modulate adequately the features of the organic ligand and ligands with weak or intermediate strength such as acetylacetonate or maltolate should not be associated to VIII and VV but only with VIV.
Data Availability Statement
All relevant data is contained within the article. All datasets analyzed in this study are included in the article and the Supplementary Material.
Author Contributions
EG conceived this study and designed the experiments and calculations. GS, SR, and J-DM performed the computational calculations. DS and GL performed the EPR measurements and VU the ESI-MS experiments. All authors contributed to the manuscript revision, read, and approved the submitted version.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
EG, GS, VU, DS, GL, and SR thank Regione Autonoma della Sardegna (grant RASSR79857) and Università di Sassari (fondo di Ateneo per la ricerca 2019) for financial support. GS and J-DM thank also Spanish MINECO (grant CTQ2017-87889-P) and Generalitat de Catalunya (2017SGR1323).
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fchem.2020.00345/full#supplementary-material
References
Amin, S. S., Cryer, K., Zhang, B., Dutta, S. K., Eaton, S. S., Anderson, O. P., et al. (2000). Chemistry and insulin-mimetic properties of Bis(acetylacetonate)oxovanadium(IV) and derivatives. Inorg. Chem. 39, 406–416. doi: 10.1021/ic9905897
Brito, F., Araujo, M. L., Martínez, J. D., Hernández, Y., Moh, A., and Lubes, V. (2009). Speciation of the vanadium(III)–acetylacetone system in 3.0 M KCl ionic medium at 25°C. J. Coord. Chem. 62, 52–62. doi: 10.1080/00958970802474763
Buglyo, P., Kiss, E., Fabian, I., Kiss, T., Sanna, D., Garribba, E., et al. (2000). Speciation and NMR relaxation studies of VO(IV) complexes with several O-donor containing ligands: oxalate, malonate, maltolate and kojate. Inorg. Chim. Acta 306, 174–183. doi: 10.1016/S0020-1693(00)00168-7
Buglyo, P., Kiss, T., Kiss, E., Sanna, D., Garribba, E., and Micera, G. (2002). Interaction between the low molecular mass components of blood serum and the VO(IV)-DHP system (DHP=1,2-dimethyl-3-hydroxy-4(1H)-pyridinone). J. Chem. Soc. Dalton Trans. 11, 2275–2282. doi: 10.1039/b200688j
Chasteen, D. N. (1981). “Vanadyl(IV) EPR spin probe. Inorganic and biochemical aspects,” in Biological Magnetic Resonance, eds L. J. J. Berliner and J. Reuben (New York, NY: Plenum Press), 53–119. doi: 10.1007/978-1-4613-3201-5_2
Costa Pessoa, J., Etcheverry, S., and Gambino, D. (2015a). Vanadium compounds in medicine. Coord. Chem. Rev. 30, 24–48. doi: 10.1016/j.ccr.2014.12.002
Costa Pessoa, J., Garribba, E., Santos, M. F. A., and Santos-Silva, T. (2015b). Vanadium and proteins: uptake, transport, structure, activity and function. Coord. Chem. Rev. 30, 49–86. doi: 10.1016/j.ccr.2015.03.016
Crans, D. C., Khan, A. R., Mahroof-Tahir, M., Mondal, S., Miller, S. M., la Cour, A., et al. (2001). Bis(acetylamido)oxovanadium(IV) complexes: solid state and solution studies. J. Chem. Soc. Dalton Trans. 3337–3345. doi: 10.1039/b101718g
Crans, D. C., Yang, L., Haase, A., and Yang, X. (2018). “Health benefits of vanadium and its potential as an anticancer agent,” in Metallo-Drugs: Development and Action of Anticancer Agents, eds A. Sigel, H. Sigel, E. Freisinger, and R. K. O. Sigel (Berlin: De Gruyter GmbH), 251–280. doi: 10.1515/9783110470734-009
Diamond, R. (1974). Real-space refinement of the structure of hen egg-white lysozyme. J. Mol. Biol. 82, 371–391. doi: 10.1016/0022-2836(74)90598-1
Elvingson, K., González Baró, A., and Pettersson, L. (1996). Speciation in vanadium bioinorganic systems. 2. An NMR, ESR, and potentiometric study of the aqueous H+-vanadate–maltol system. Inorg. Chem. 35, 3388–3393. doi: 10.1021/ic951195s
Frisch, M. J., Trucks, G. W., Schlegel, H. B., Scuseria, G. E., Robb, M. A., Cheeseman, J. R., et al. (2010). Gaussian 09, Revision D.01. Wallingford, CT: Gaussian, Inc.
Fu, Y., Wang, Q., Yang, X.-G., Yang, X.-D., and Wang, K. (2008). Vanadyl bisacetylacetonate induced G1/S cell cycle arrest via high-intensity ERK phosphorylation in HepG2 cells. J. Biol. Inorg. Chem. 13, 1001–1009. doi: 10.1007/s00775-008-0387-2
Jakusch, T., and Kiss, T. (2017). In vitro study of the antidiabetic behavior of vanadium compounds. Coord. Chem. Rev. 351, 118–126. doi: 10.1016/j.ccr.2017.04.007
Jones, G., Willett, P., Glen, R. C., Leach, A. R., and Taylor, R. (1997). Development and validation of a genetic algorithm for flexible docking1. J. Mol. Biol. 267, 727–748. doi: 10.1006/jmbi.1996.0897
Kioseoglou, E., Petanidis, S., Gabriel, C., and Salifoglou, A. (2015). The chemistry and biology of vanadium compounds in cancer therapeutics. Coord. Chem. Rev. 87–105. doi: 10.1016/j.ccr.2015.03.010
Lindqvist, Y., Schneider, G., and Vihko, P. (1994). Crystal structures of rat acid phosphatase complexed with the transition-state analogs vanadate and molybdate. Eur. J. Biochem. 221, 139–142. doi: 10.1111/j.1432-1033.1994.tb18722.x
Liu, J.-C., Yu, Y., Wang, G., Wang, K., and Yang, X.-G. (2013). Bis(acetylacetonato)-oxovanadium(IV), bis(maltolato)-oxovanadium(IV) and sodium metavanadate induce antilipolytic effects by regulating hormone-sensitive lipase and perilipin via activation of Akt. Metallomics 5, 813–820. doi: 10.1039/c3mt00001j
Lovell, S. C., Word, J. M., Richardson, J. S., and Richardson, D. C. (2000). The penultimate rotamer library. Proteins 40, 389–408. doi: 10.1002/1097-0134(20000815)40:3<389::AID-PROT50>3.0.CO;2-2
Makinen, M. W., and Brady, M. J. (2002). Structural origins of the insulin-mimetic activity of Bis(acetylacetonato)oxovanadium(IV). J. Biol. Chem. 277, 12215–12220. doi: 10.1074/jbc.M110798200
Makinen, M. W., Rivera, S. E., Zhou, K. I., and Brady, M. J. (2007). “Enhancement of insulin action by bis(Acetylacetonato)oxovanadium(IV) occurs through uncompetitive inhibition of protein Tyrosine phosphatase-1B,” in Vanadium: The Versatile Metal, eds K. Kustin, D. C. Crans and J. Costa Pessoa (American Chemical Society, 82–92.
Makinen, M. W., and Salehitazangi, M. (2014). The structural basis of action of vanadyl (VO2+) chelates in cells. Coord. Chem. Rev. 279, 1–22. doi: 10.1016/j.ccr.2014.07.003
Marenich, A. V., Cramer, C. J., and Truhlar, D. G. (2009). Universal solvation model based on solute electron density and on a continuum model of the solvent defined by the bulk dielectric constant and atomic surface tensions. J. Phys. Chem. B 113, 6378–6396. doi: 10.1021/jp810292n
Ou, H., Yan, L., Mustafi, D., Makinen, M. W., and Brady, M. J. (2005). The vanadyl (VO2+) chelate bis(acetylacetonato)oxovanadium(IV) potentiates tyrosine phosphorylation of the insulin receptor. J. Biol. Inorg. Chem. 10, 874–886. doi: 10.1007/s00775-005-0037-x
Pettersen, E. F., Goddard, T. D., Huang, C. C., Couch, G. S., Greenblatt, D. M., Meng, E. C., et al. (2004). UCSF Chimera-A visualization system for exploratory research and analysis. J. Comput. Chem. 25, 1605–1612. doi: 10.1002/jcc.20084
Qureshi, I. A., Ferron, F., Seh, C. C., Cheung, P., and Lescar, J. (2009). Crystallographic structure of ubiquitin in complex with cadmium ions. BMC Res. Notes 2:251. doi: 10.1186/1756-0500-2-251
Rehder, D. (2016). Perspectives for vanadium in health issues. Future Med. Chem. 8, 325–338. doi: 10.4155/fmc.15.187
Reul, B. A., Amin, S. S., Buchet, J.-P., Ongemba, L. N., Crans, D. C., and Brichard, S. M. (1999). Effects of vanadium complexes with organic ligands on glucose metabolism: a comparison study in diabetic rats. Br. J. Pharmacol. 126, 467–477. doi: 10.1038/sj.bjp.0702311
Rodríguez-Guerra Pedregal, J., Sciortino, G., Guasp, J., Municoy, M., and Maréchal, J.-D. (2017). GaudiMM: a modular multi-objective platform for molecular modeling. J. Comput. Chem. 38, 2118–2126. doi: 10.1002/jcc.24847
Sakurai, H., Yoshikawa, Y., and Yasui, H. (2008). Current state for the development of metallopharmaceutics and anti-diabetic metal complexes. Chem. Soc. Rev. 37, 2383–2392. doi: 10.1039/b710347f
Sanna, D., Garribba, E., and Micera, G. (2009). Interaction of VO2+ ion with human serum transferrin and albumin. J. Inorg. Biochem. 103, 648–655. doi: 10.1016/j.jinorgbio.2009.01.002
Scior, T., Guevara-Garcia, J. A., Do, Q.-T., Bernard, P., and Laufer, S. (2016). Why antidiabetic vanadium complexes are not in the pipeline of “big pharma” drug research? A critical review. Curr. Med. Chem. 23, 2874–2891. doi: 10.2174/0929867323666160321121138
Sciortino, G., Garribba, E., and Maréchal, J.-D. (2019a). Validation and applications of protein–ligand docking approaches improved for metalloligands with multiple vacant sites. Inorg. Chem. 58, 294–306. doi: 10.1021/acs.inorgchem.8b02374
Sciortino, G., Garribba, E., Rodríguez-Guerra Pedregal, J., and Marechal, J.-D. (2019b). Simple coordination geometry descriptors allow to accurately predict metal-binding sites in proteins. ACS Omega 4, 3726–3731. doi: 10.1021/acsomega.8b03457
Sciortino, G., Rodríguez-Guerra Pedregal, J., Lledós, A., Garribba, E., and Maréchal, J.-D. (2018a). Prediction of the interaction of metallic moieties with proteins: an update for protein-ligand docking techniques. J. Comput. Chem. 39, 42–51. doi: 10.1002/jcc.25080
Sciortino, G., Sanna, D., Ugone, V., Lledós, A., Maréchal, J.-D., and Garribba, E. (2018b). Decoding surface interaction of VIVO metallodrug candidates with lysozyme. Inorg. Chem. 57, 4456–4469. doi: 10.1021/acs.inorgchem.8b00134
Sciortino, G., Sanna, D., Ugone, V., Micera, G., Lledós, A., Maréchal, J.-D., et al. (2017). Elucidation of binding site and chiral specificity of oxidovanadium drugs with lysozyme through theoretical calculations. Inorg. Chem. 56, 12938–12951. doi: 10.1021/acs.inorgchem.7b01732
Smith, T. S. II., LoBrutto, R., and Pecoraro, V. L. (2002). Paramagnetic spectroscopy of vanadyl complexes and its applications to biological systems. Coord. Chem. Rev. 228, 1–18. doi: 10.1016/S0010-8545(01)00437-4
Thompson, K. H., Lichter, J., LeBel, C., Scaife, M. C., McNeill, J. H., and Orvig, C. (2009). Vanadium treatment of type 2 diabetes: a view to the future. J. Inorg. Biochem. 103, 554–558. doi: 10.1016/j.jinorgbio.2008.12.003
Thompson, K. H., and Orvig, C. (2006). Vanadium in diabetes: 100 years from Phase 0 to Phase I. J. Inorg. Biochem. 100, 1925–1935. doi: 10.1016/j.jinorgbio.2006.08.016
Ugone, V., Sanna, D., Sciortino, G., Marechal, J. D., and Garribba, E. (2019). Interaction of vanadium(IV) species with Ubiquitin: a combined instrumental and computational approach. Inorg. Chem. 58, 8064–8078. doi: 10.1021/acs.inorgchem.9b00807
Wu, J.-X., Hong, Y.-H., and Yang, X.-G. (2016). Bis(acetylacetonato)-oxidovanadium(IV) and sodium metavanadate inhibit cell proliferation via ROS-induced sustained MAPK/ERK activation but with elevated AKT activity in human pancreatic cancer AsPC-1 cells. J. Biol. Inorg. Chem. 21, 919–929. doi: 10.1007/s00775-016-1389-0
Keywords: metal drugs, vanadium, anticancer action, antidiabetic action, proteins, transport in the organism, drug design
Citation: Sciortino G, Ugone V, Sanna D, Lubinu G, Ruggiu S, Maréchal J-D and Garribba E (2020) Biospeciation of Potential Vanadium Drugs of Acetylacetonate in the Presence of Proteins. Front. Chem. 8:345. doi: 10.3389/fchem.2020.00345
Received: 11 February 2020; Accepted: 02 April 2020;
Published: 07 May 2020.
Edited by:
Arthur David Tinoco, University of Puerto Rico, Puerto RicoReviewed by:
Sergio Armando Loza Rosas, University of Boyaca, ColombiaRitika Gautam, Indian Institute of Technology Kanpur, India
Copyright © 2020 Sciortino, Ugone, Sanna, Lubinu, Ruggiu, Maréchal and Garribba. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Eugenio Garribba, garribba@uniss.it