- 1Centro de Química del Instituto de Ciencias, Benemérita Universidad Autónoma de Puebla, Mexico
- 2Facultad de Ciencias Químicas, Benemérita Universidad Autónoma de Puebla, Mexico
- 3Instituto de Física “Luis Rivera Terrazas”, Benemérita Universidad Autónoma de Puebla, Mexico
- 4Departamento de Ciencias Químico-Biológicas, Universidad de las Américas Puebla, Mexico
Cytosine, a DNA and RNA building-block, and Metformin, the most widely prescribed drug for the treatment of Type 2 Diabetes mellitus were made to react separately with ammonium or sodium metavanadates in acidic aqueous solutions to obtain two polyoxovanadate salts with a 6:1 ratio of cation-anion. Thus, compounds [HCyt]6[V10O28]·4H2O, 1 and [HMetf]6[V10O28]·6H2O, 2 (where HCyt = Cytosinium cation, [C4H6N3O]+ and HMetf = Metforminium cation, [C4H12N5]+) were obtained and characterized by elemental analysis, single crystal X-ray diffraction, vibrational spectroscopy (IR and Raman), solution 51V-NMR, thermogravimetric analysis (TGA-DTGA), as well as, theoretical methods. Both compounds crystallized in P1¯ space group with Z' = 1/2, where the anionic charge of the centrosymmetric ion [V10O28]6− is balanced by six Cytosinium and six Metforminium counterions, respectively. Compound 1 is stabilized by π-π stacking interactions coming from the aromatic rings of HCyt cations, as denoted by close contacts of 3.63 Å. On the other hand, guanidinium moieties from the non-planar HMetf in Compound 2 interact with decavanadate μ2-O atoms via N−H···O hydrogen bonds. The vibrational spectroscopic data of both IR and Raman spectra show that the dominant bands in the 1000-450 cm−1 range are due to the symmetric and asymmetric ν(V−O) vibrational modes. In solution, 51V-NMR experiments of both compounds show that polyoxovanadate species are progressively transformed into the monomeric, dimeric and tetrameric oxovanadates. The thermal stability behavior suggests a similar molecular mechanism regarding the loss of water molecules and the decomposition of the organic counterions. Yet, no changes were observed in the TGA range of 540–580°C due to the stability of the [V10O28]6− fragment. Dispersion-corrected density functional theory (DFT-D) calculations were carried out to model the compounds in aqueous phase using a polarized continuum model calculation. Optimized structures were obtained and the main non-covalent interactions were characterized. Biological activities of these compounds are also under investigation. The combination of two therapeutic agents opens up a window toward the generation of potential metalopharmaceuticals with new and exciting pharmacological properties.
Introduction
Vanadium bioinorganic chemistry is becoming a very important topic in the last decades. The information about vanadium in biological systems is increasing in giant steps, so as a guide for our research we have considered many highlights in the history of the field (See Scheme 1). In this introduction, relevant information will be presented in a limited way.
Scheme 1. A brief timeline of the vanadium bioinorganic chemistry, highlighting the vanadium role as an antidiabetic, and anticancer agent (Cantley et al., 1977; Thompson et al., 1984; Vilter, 1984; Orvig et al., 1995; Elberg et al., 1997; Crans, 2000; Narla et al., 2001; Rehder, 2008; Aureliano and Crans, 2009; Steens et al., 2010; Aureliano, 2017; Bijelic et al., 2018).
Vanadium(IV) and (V) based-compounds have been tested for their potential biomedical use in the treatment of Diabetes mellitus, cancer, bacterial diseases and viral infections (Crans, 2000; Rehder, 2008, 2012, 2015a,b). Its therapeutic effects have been attributed mainly to the fact that vanadate has structural and electronic similarities to phosphate. In this sense, vanadium species like vanadate, (H2VO4)−, can adopt a stable trigonal-bipyramidal geometry similar to that of phosphate transition state in phosphate-metabolizing enzymes, and inhibit its biological activity (Costa-Pessoa et al., 2015; Crans, 2015; Dorsey et al., 2018). Since most of the enzymes inhibited by vanadate participate in key intracellular signaling processes, vanadium has been considered a transition element with relevant medicinal applications (Rehder, 2015a,b; Del Carpio et al., 2018; Selman et al., 2018).
Unlike phosphate, vanadate(V) solutions can form, under specific conditions, polyoxovanadate species such as the decavanadate anion, [V10O28]6− (Hayashi, 2011). The structure of this oxo-cluster is stable at acidic pH range and structurally contains ten vanadium atoms assembled into a compact structure with unit cell dimensions of 8.3 Å x 7.7 Å x 5.4 Å, where the V5+ metal ions occupy the octahedral interstices in the ten [VO6] units (Swallow and Barnes, 1964; Aureliano, 2011).
With respect to the biological activity of decameric species, it has been shown that decavanadate induce several changes in biological processes through the interaction with biological systems, such as myosin, actin, and ion pumps, which are major proteins implicated in muscle contraction and its regulation (Aureliano, 2011, 2016, 2017; Winkler et al., 2017). Furthermore, it has been reported that decavanadate has an insulin-enhancing activity in streptozotocin-induced diabetic rats (Pereira et al., 2009) and in murine diabetic models induced by hypercaloric diets or alloxan (Treviño et al., 2015, 2016, 2018). Recently it has been suggested that it's in vitro anticancer activity may be related to the inhibitory effects on specific enzymes involved in tumor proliferation (Aureliano, 2017; Bijelic et al., 2018). Also, the structure determination of the TRMP4 calcium channel has revealed two binding sites for the decavanadate anions, Thus, open up new ways for decavanadate biological actions (Winkler et al., 2017).
Because of its high anionic charge, decavanadate is mainly stabilized by counterions through electrostatic interactions and hydrogen bonds. In previous works, decavanadate has been combined with a variety of cations, including protons, alkali-metal ions, ammonium, alkylammonium, phosphonium, organic and organometallic cations (Chinea et al., 2000; Lee and Joo, 2003; Correia et al., 2004; Kojima et al., 2011; Bartošová et al., 2012; Kioseoglou et al., 2013; Mal et al., 2013; Sánchez-Lombardo et al., 2014; Crans, 2015; Sánchez-Lara et al., 2015, 2016b). In this work, we combined decavanadate with two biologically relevant molecules Cytosine and Metformin (Scheme 2) with the aim to synthesize compounds with potential pharmacological activity.
About the counterions used in this work, Cytosine is one of the naturally occurring nitrogenous bases found in DNA and RNA. Some pyrimidine derivatives exhibit anticancer activity, and some compounds have been used for the treatment of fungal infections (Vermes et al., 2000; Parker, 2009). On the other hand, N,N-dimethylbiguanide or Metformin, as it is known worldwide, is the favorite prescription drug to control the effects of type 2 D. mellitus (UK Prospective Diabetes Study (UKPDS) Group., 1998) and is known to affect the cellular housekeeping of copper (Repišcák et al., 2014). Although Metformin is a relatively simple molecule, it possess a versatile behavior as indicated by its capacity of forming salts and coordination compounds in which it participates as a dicationic, monocationic, neutral, and anionic or dianionic species (Zhu et al., 2002).
Besides its role as an insulin sensitizer, new evidence suggests that administration of Metformin may reduce the risk of several types of cancer, including breast, pancreatic and colon cancer (Evans et al., 2005; Libby et al., 2009; Kasznicki et al., 2014). Although the exact molecular mechanism of action for involving Metformin remains unclear, it is known that one of its primary targets is complex I of the mitochondrial electron transport system affecting ATP generation (El-Mir et al., 2000).
Following on from this, the primary focus of this paper is to present two polyoxovanadates-based compounds with potential pharmacological activity. Special emphasis was placed in the compounds' physical and chemical characterization and strengthened with the theoretical analysis. We thus think that the structural study here presented to justify an in-depth biological study to understand, for example, the outcome on the bioactivity of some decavanadate compounds with specific counterions (Yraola et al., 2007; Zorzano et al., 2009; Treviño et al., 2015, 2016, 2018).
Experimental Section
All manipulations were carried out at room temperature and with no special solvent and reagent purification. Infrared spectra were obtained in KBr pellets in the range from 400-4000 cm−1 by using an IR Digilab, Mod. Scmitar FT-IR spectrophotometer. Raman spectra were obtained at room temperature in backscattering configuration using the 633 nm line of a He-Ne laser as an excitation source by using a LabRAM HR-Olympus Micro Raman system. The 1H, 13C were obtained at 500 MHz and 125 MHz respectively, while the 51V-NMR spectra were recorded at 131.5 MHz with a Bruker AVANCE III 500 MHz spectrometer using deuterated water (D2O). Chemical shifts were referenced to VOCl3 as an external standard in the case of 51V-NMR. The elemental analysis was performed by Elemental Analyser CHNS/O Thermo Scientific Flash 2000. Thermogravimetric (TGA) and differential thermal (DTA) analysis were carried out under continuous N2 flow on an STA 2500 Regulus differential scanning calorimeter (Netzsch Instruments, Burlington, MA). Samples were weighed (~10 mg) and heated from 30 to 700 °C, at a heating rate of 25 °C/min in alumina pans. An empty pan was used as a reference. Single-crystal X-ray data were recorded with an Agilent Gemini A diffractometer, software SHELX-2014/7 (Sheldrick, 2015). Selected crystal data and details of the structure determination of the compounds are summarized in Table 1. CCDC numbers 1850706 (Compound 1) and 1850707 (Compound 2) contain the supplementary crystallographic data for this paper. These data can be obtained for free at http://www.ccdc.cam.ac.uk/conts/ retrieving.html (or from the CCDC, 12 Union Road, Cambridge CB2 1EZ, UK; Fax: +44-1223-336-033; e-mail address: deposit@ccdc.cam.ac.uk). The supramolecular networks were studied by using Mercury CSD (release 3.1.2) (Macrae et al., 2008), which, together with OLEX-2 (Dolomanov et al., 2009) were used to produce crystallographic artwork.
Table 1. Single crystal data and structure refinement details for [HCyt]6[V10O28]·4H2O, 1, and [HMetf]6[V10O28]·6H2O, 2.
Computational Methods
The structural and electronic properties of compounds [HCyt]6[V10O28]·4H2O, 1, and [HMet]6[V10O28]·6H2O, 2, were computed from theoretical calculations based on the density functional theory (DFT) (Hohenberg and Kohn, 1964). Geometry optimization of the asymmetric unit of Compound 1 and 2 was obtained using the pure functional B97-D3 (Grimme et al., 2011). B97-D3 is a Grimme's functional that includes the Becke and Johnson dispersion corrections (BJ-damping), recommended for non-covalent interactions (Grimme et al., 2011). The split-valence 6-31G(d) basis set (Rassolov et al., 1998) used in the C, H, O, and N atoms includes a single set of Gaussian polarization functions. A LanL2MB basis set (Hay and Wadt, 1985) and an effective core potential (ECP) which replaces the effects of the inner core electrons with a pseudopotential specific for transitions metal atoms, were used for the V atom. Aqueous solvent effects were calculated with the polarized continuum model using the conductor-like polarizable continuum model (CPCM) (Barone et al., 1998). Calculations were performed with the Gaussian16 program (Frisch et al., 2016) and visualization of the results was carried out with the GaussianView 6.0.16 program (Dennington et al., 2016).
Synthesis of [HCyt]6[V10O28]·4H2O (1)
First, 0.5 g of NH4VO3 was dissolved in 30 mL of distilled water and heated up to dissolution. Then, three drops of concentrated hydrochloric acid (37 %) were added at room temperature to allow the formation of decavanadate anions. After obtaining an orange solution, 0.111 g (1 mmol) of cytosine previously dissolved in 20 mL of distilled water was added dropwise under stirring, obtaining a crystalline orange precipitate, which was filtered from the mother solution (pH = 5.8). Recrystallization from 20 mL of hot water gave a low yield of block-shaped orange crystals, which have low solubility in water. 51V NMR (D2O): −515.84, −500.31, −423.55. IR (KBr, cm−1): 3383, 3321 [vas(N–H)], 3170 [vs(N–H)], 2984 [v(C–H)], 1718 [v(C = O)], 1686 [v(C = C)], 1656 [v(C = N)], 1545 [δ(N–H)], 987, 952 [v(V = O)], 820, 737 [vas(V–O–V)], 603, 571[vas(V–O–V)]. Anal. Calc. for C24H44N18O38V10 (MW = 1702.17 g/mol) C, 16.93%; H, 2.60%; N, 14.81%. Found: C, 17.23%; H, 2.76%; N, 14.94%.
Synthesis of [HMetf]6[V10O28]·6H2O (2)
Two 850-mg-Metformin hydrochloride tablets (Alpharma laboratories) were each crushed using a mortar and pestle and added to a solution containing 0.5 g of NaVO3, previously dissolved in 30 mL of distilled water. Three drops of concentrated hydrochloric acid (37%) were added by stirring and then filtered. Needle-shaped orange crystals of Compound 2, were isolated from their mother solution (pH = 6.3) and filtered after 3 days. It is important to note that the pH value was a key factor to obtain this compound since the metformin molecule can act as monocatioinic or dicationic species at different acidic pH values (Chatkon et al., 2014; Sánchez-Lombardo et al., 2014). The slightly acidic medium of the mother solution maintains the monoprotonation of the metformin molecule and eventually the formation of Compound 2. Second, the monoprotonated metformin is a good hydrogen-bonding donor leading to better solvation in the aqueous medium, therefore, allowing a good solubility of this compound. Yield (based on vanadium): 0.260 g, 35%. 1H NMR (D2O): δ = 3.027 (s, 6H, CH3), 13C (D2O): δ = 37.7 (CH3), 158.7 (C–NH2), 160.2 (C = ). 51V (D2O): −514.43, −500.02, −422.31. IR (KBr, cm−1): 3502 [v(O–H)], 3396–3321 [vas(N–H)], 3170 [vs(N–H)], 1624, 1569 [v(C = N)], 1481, 1419 [δ(CH3)], 1053 [δ(N–H)], 945 [vs(V = O)], 830, 736 [vas(V–O–V)], 603, 571[vs(V–O–V)]. Anal. Calc. for C24H84N30O34V10 (MW = 1846.61 g/mol) C, 15.61%; H, 4.60%; N, 22.75%. Found: C, 16%; H, 4.5%; N, 22.74%.
Results and Discussion
Synthetic Strategy
Compound [HCyt]6[V10O28]·4H2O, 1, was isolated from an acidic mother solution of NH4VO3 with a pH of around 6, which allowed the protonation of the cytosine molecules and the formation of decavanadate anions. The first attempt to obtain a compound based on [V10O28]6− units and cytosine was initially explored by our group using sodium metavanadate salt NaVO3 as a precursor of polyoxovanadate ions. However, as a result, we obtained a compound previously published, where the sodium ions are present in the crystal structure forming metal-complexes with the cytosine molecules and the [V10O28]6− ions (Bošnjakovic-Pavlovic et al., 2009). With the purpose of replacing the Na+ ions to get Compound 1, the starting material NaVO3 was replaced by NH4VO3. In this way, the cytosinium cations were incorporated in the resulting compound, without the presence of ions, and the water molecules were occluded during the crystallization process as lattice solvent.
On the other hand, crystals of [HMetf]6[V10O28]·6H2O, 2, were isolated in good yield (35%) from an aqueous solution of sodium metavanadate NaVO3 at pH 6.3. Again, the pH value and the vanadium source were the key factors to obtain this Compound. The synthetic strategy was, first, take into account that Metformin molecule is a moderately strong base with two pKa values of 2.8 and 11.5 (Kathuria et al., 2018) and therefore can act as monocationic or dicationic species at different acidic pH values. Second, we needed to avoid the formation of metforminium(2+) decavanadates salts, previously reported by us (Sánchez-Lombardo et al., 2014). From this, we decided to use a slightly acidic medium in order to maintain the monoprotonation of the Metformin molecule present in the Metformin hydrochloride tablets, and eventually to get the Compound 2. Is important to mention herein that the solubility of this compound is very high in contrast with the Compound 1, may be due to protonated form of Metformin (HMetf) is a good hydrogen-bonding donor, leading to better solvation in aqueous medium allows, therefore, a good solubility of this compound. This point is particularly relevant to carry out biological testing.
Crystal Structure Description
Single-crystal X-ray structural analysis shows that Compound 1 crystallizes in the centrosymmetric P space group (Table 1) with the asymmetric unit, containing one-half of decavanadate anion placed in an inversion center, three independent HCyt cations in the amino-oxo tautomeric form and two lattice water molecules, each located in general positions are depicted in Figure 1. The structure of the [V10O28]6− cluster is well known and consists of an arrangement of a 10-edge-shared [VO6] octahedron close to a local D2h symmetry, though the crystallographic symmetry is Ci.
Figure 1. ORTEP-like view of the asymmetric unit of Compound 1, with displacement ellipsoids for non-H atoms drawn at the 50% probability level.
In agreement with the Mogul geometry check (Macrae et al., 2008), the V–O bond lengths and the V–O–V angles are found in normal ranges; comparable with other structures containing this oxo-cluster and showed no unusual geometrical parameters. The backbone of the [V10O28]6− anion is a very rigid entity that does not present significant structural changes caused by non-covalent interactions. The (3-)anion charge in the asymmetric unit is therefore stabilized by three HCyt cations, protonated at N2, N12, and N22 sites. All the H-atoms of HCyt were found from the difference-Fourier map and refined isotropically. The residual electronic density is actually low for room temperature data of 0.80 e-/Å3. From this, Compound 1 can be formulated as [HCyt]6[V10O28]·4H2O.
Ionization of the cytosine molecule in an acid medium leads to the introduction of a third donor group. This protonation causes an increase in the internal C–N–C bond angles of the pyrimidine rings in the three independent organic molecules. For example, the C1–N2–C2 angle of 124.81(3)° is larger in cation N1 than in the corresponding value found in the 118.90°-non-protonated cytosine molecule (Lee and Wang, 2010). Also, the short bond lengths between C2–N3 and N1–C4 of 1.302(5) and 1.358(5) Å, respectively, show a resonance effect of π-electron density on the entire molecule. These values are very similar in the three independent cations.
According to the Cambridge Structural Database (CSD, version 5.39, last updated November 2017; Groom et al., 2016), there are no direct structure precedents of Compound 1 in the crystallographic literature, but some complexes containing HCyt or neutral cytosine interacting with polyoxometalate ions have been structurally characterized by X-ray diffraction. One of these structures is an inorganic hybrid polymer built from P2Mo5O23 clusters and Cu(II)–cytosine subunits (Weng et al., 2002). There are two structures based on hetero-polyoxometalates, [SeMo5O21]4− and [(HAsO4)Mo6O19]6− (Nagazi and Haddad, 2014; Ayed et al., 2015). The last and possibly the closest structure to Compound 1 is a compound co-crystallized with [V10O28]6− and Na+ ions (Bošnjakovic-Pavlovic et al., 2009). However, a duplex structure containing a cytosine/cytosinium dimer is present in this case. In Compound 1, dimer formation is disrupted despite a cooperative triple hydrogen bonding of the duplex structure, possibly owing to the low pKa value of the hydrochloric acid used to adjust the pH that produces the protonation of all the cytosine molecules, avoiding the formation of the hemicytosinium pair (Perumalla et al., 2013).
The supramolecular structure of Compound 1 is dominated by classic N−H···O, O–H···O, C–H···O and π-stacking interaction among heterocyclic rings. In this regard, the independent cations N1 and N11 of HCyt interact with [V10O28]6− through N−H···O hydrogen bonds, which involve protonated pyrimidine rings as donor groups and two μ2-O sites of the anion as acceptors, with a N1–H1···O9 distance of 2.77(4) Å and a N11–H11···O12 distance of 2.73(4) Å, as well as N–H···O angles of 169 and 172°, respectively (Table 2). These cations generate π-π contacts characterized by a centroid-to-centroid distance of 3.92(3) Å. The remaining N2 and N3 atoms from cation N1 comprise planar rings represented by (Etter et al., 1990) with μ2-O decavanadate atoms.
On the other hand, N11- and N21-based HCyt cations interact with each other via one hydrogen bond between the acceptor O carboxylic atom and the exocyclic amino group, with distance of 2.81(5) Å for O21C···H131–N13 and 2.79(5) Å for O11B···H231–N23, forming a C(12) (Etter et al., 1990) infinite chain extended in the [101] direction, which is stacked through the inversion center coming from the space-group symmetry. In this regard, the π-π parallel-displaced stacking interaction between N11 and N21 cations in the layer is characterized by a centroid-to-centroid distance of 3.63(3) Å and by a displacement angle of approximately 20°. These chains are hydrogen-bonded with decavanadate oxygen atoms forming and ring motifs (Figure 2). It is important to mention that distances of 3.6–3.8 Å between centroids of the π-π stacking interactions are also found in DNA and RNA structures. The relevance of this phenomena is unclear at the moment. However, it can be hypothezised that polyoxometalates could be used as templates for base-base linkages and base-base pairing in the early development of life (Rehder, 2011; Cafferty et al., 2016).
Figure 2. Partial crystal structure of Compound 1, showing the supramolecular network, which is based on the hydrogen bonds (dashed lines) formed between HCyt and decavanadates: Orange H-bonds are used to form the C(12) chain that is connected to [V10O28]6− anions by purple H-bonds in the crystal structure. The inset displays the π-π contacts between two non-symmetry related HCyt.
The molecular formula of Compound 2 has been determined by single-crystal X-ray diffraction as [HMetf]6[V10O28]·6H2O (Table 1). This compound also crystallizes in the triclinic P space group with half of the decavanadate anion (Z' = ), three independent monoprotonated HMetf cations and three lattice water molecules in the asymmetric unit (Figure 3). Therefore, six HMetf neutralize the [V10O28]6− anion charge in the unit cell.
Figure 3. ORTEP-like view of the asymmetric unit of Compound 2, with displacement ellipsoids for non-H atoms drawn at the 50% probability level. Only the major component of disorder for one HMetf is shown.
In the refinement process, one of the HMetf is split into two sites with refined occupancies of 0.616(5) in part A and 0.382(5) in part B. It is not surprising to detect a positional disorder for this cation because this limited disorder should be a stabilizing factor for the whole crystal structure. On the other hand, O52, O53, O54, and O55 from water molecules were refined more effectively with an occupancy factor of 1/2. The largest peak in the final difference electron density synthesis was 1.01 e-/Å3 and the largest hole was 0.52 e-/Å3.
The non-planar HMetf in Compound 2 shows the geometrical features of Metformin hydrochloride (Form B) (Childs et al., 2004), where the orientation of the –C(NH2) and –N(CH3)2 groups are found on the opposite side due to a decrease in van der Waals repulsion. The torsion angles for cations N10, N20 and N30A are 158.3(6)°, 137.2(4)°, and 144.7(13)°, respectively. This type of conformation is observed in several Metforminium(1+) based-salts containing chloride, nitrate, acetate, salicylate, squarate, and carbonate (Zhu et al., 2003; Childs et al., 2004; Olar et al., 2010; Pérez-Fernández et al., 2013; Serb et al., 2014; Dong et al., 2015).
There are some tautomeric forms reported for HMetf in the solid state, all of which are fully delocalized (Scheme 3) (Pérez-Fernández et al., 2013; Kathuria et al., 2018). Cation N10 present in Compound 2 corresponds to tautomer (c), while N30A and N20 cations correspond to tautomer (d). Although the central N atoms can act as acceptor groups, only N23 accepts one hydrogen bond from O55, with a distance O55–H55A···N23 of 2.84(11) Å (Table 3).
By analyzing the bond lengths in HMetf, we observed that the N–(CH3)2 methyl group bond lengths correspond to single bonds with a distance ranging from 1.44(7) to 1.47(9) Å. For the C–N bonds, the bond distances range from 1.31(7) to 1.35(7) Å and the central C–N–C angle is reduced to nearly 120°, all of which suggests a π-electron density delocalization across the biguanide group.
Regarding structures crystallized with decavanadate and the anti-diabetic drug Metformin, our research group has previously obtained three decavanadate Metforminium salts at low pH values, all of them characterized by X-ray diffraction (Sánchez-Lombardo et al., 2014). However, in these cases, the Metforminium ions are acting as a dicationic species. Debbie Crans group, on the other hand, has described a water-insoluble double salt containing decavanadate and a mixture of Metforminium(1+) and a protonated guanylurea counterion (Chatkon et al., 2014). The Metformin molecule in the reported compound presents equal proton distributions and similar structural features to those of HMetf found in Compound 2.
Further describing the structure of Compound 2, all guanidinium moieties of HMetf cations are involved in hydrogen bonds with decavanadate anion, through N–H···O with and ring motifs with distances ranging from 2.89(5) to 3.25(7) Å and characterized by angles D–H···A > 140° (Figure 4). These guanidinium-decavanadate interactions are also favored by Coulombic forces.
Figure 4. Part of the crystal structure of Compound 2, showing the interaction between HMetf and [V10O28]6−. The main hydrogen bonds between the cations and anion that form the (8) and (6) rings are shown as purple dashed lines.
A complete supramolecular analysis of this compound is difficult to attain because one HMetf is split into two fractional parts and the disorder in O51, O52, O53, and O54 affect the water H-bond positions. However, considering only the N10 and N20 cations, we observe the formation of a four-level ring involving two non-symmetry related HMetf, that interact with terminal and bridging μ2- and μ3-O from [V10O28]6−. This tetramer is formed by N11-H11A···O12 3.10(6) Å, N15-H15B···O13 2.93(8) Å, N20-H20B ⋯ O3 2.93(6) Å, and N21-H21A ⋯ O11 2.95(5) Å hydrogen bonds. Also, a centrosymmetric ring is present in this compound. In this case, two inversion-related HMetf cations interact with terminal and bridging μ3-O atoms from decavanadate as acceptor groups. This third-level motif with pattern R(> k < l > k < l) shows a chair-like conformation. Both supramolecular patterns are represented in Figure 5. It is important to notice the similarity between the binding sites of TRPM4 calcium channel which are rich in arginine and lysine residues and Compound 2. In this sense, the interaction of the biguanide Metformin, resembles the interactions of the guanidinium ions in the TRMP4 decavanadate binding sites. (see Figure 4) PDB: 5WP6.
Figure 5. Hydrogen-bond patterns in Compound 2, formed through hydrogen bonds (purple dashed lines) between two HMetf cations and two symmetry-related [V10O28]6−moities.
Solid-State Vibrational Spectroscopy
FTIR and Raman spectra of Compounds 1-2 are shown in Figure 6 and can be compared with those of neutral cytosine and Metformin hydrochloride, respectively. Furthermore, the Raman spectra of Compounds 1-2 are compared with that of the ammonium decavanadate salt, [NH4]6[V10O28]·6H2O, for a better assignment of the Raman bands.
Figure 6. Vibrational spectroscopy of Compounds 1-2 recorded in solid-state. (A) and (B) are the FT-IR and Raman spectra for Cytosine and Compound 1, while (C) and (D) are the FT-IR and Raman spectra for Metformin hydrochloride and Compound 2, respectively.
The middle IR (Figure 6A) of Compound 1 shows a set of peaks in the high-frequency region at 3486-2684 cm−1. νas(N–H) appears at 3383 and 3321 cm−1, while νs(N–H) appears at 3135 cm−1. v(Csp2−−H) is observed in the range of 2984–2694 cm−1 (Rozenberg et al., 2004). Strong IR vibrations seen at 1,718, 1,686, and 1,656 cm−1 have been assigned to v(C = O), v(C = C), and v(C = N), respectively. It is known that the position of these bands is shifted upon protonation of the cytosine ring by shortening bond lengths (Rasheed and Ahmad, 2010; Sridhar et al., 2012).
The bands observed at 987 and 952 cm−1 in the IR spectra and the corresponding Raman peaks (Figure 6B) at 992 and 958 cm−1 are assigned to symmetric stretching modes of the terminal V = O bonds of the decavanadate anion. However, the (V10O28)6− vibrations are mixed with the corresponding vibrations of the HCyt molecules (see Table 4).
The bands at the 833–737 cm−1 region in the infrared spectrum are attributed to the asymmetric stretching modes of (V–O–V) units, overlapped with the cytosine ring “breathing” vibration observed at 792 cm−1 in the neutral cytosine molecule. In Raman spectra, the asymmetric stretching mode of (V–O–V) is found at 835 cm−1 as a medium-intensity peak, and the intense band at 793 cm−1 is due to the ring-breathing mode of the HCyt (Mathlouthi et al., 1986; Frost et al., 2005; Frost and Palmer, 2011; Pavliuk et al., 2014; Madzharova et al., 2016; Sánchez-Lara et al., 2016b).
Two IR intense bands at 572 and 513 cm−1, and the corresponding to Raman band at 595 cm−1 are observed. These frequencies were assigned to νs(V–O–V). The Raman bands at 423–320 cm−1 may be attributed to the δ(V–O–V). Bands below 240 cm−1 are attributed to V−O bonds and lattice vibrations (Frost et al., 2005; Frost and Palmer, 2011).
For Compound 2, the medium-intensity bands in the high-frequency region in the IR spectra (Figure 6C) at 3,395–3,289 cm−1 are attributed to the vas(N–H) bond, while the band corresponding to vs(N–H) is observed in the region of 3,200–3,170 cm−1. On the other hand, the very weak bands at 2,972 and 2,936 cm−1 correspond to the symmetric and asymmetric stretching vibrations of the methyl groups of HMetf, respectively (Gunasekaran et al., 2006). No significant difference between the IR spectra of Compound (2) and Metformin hydrochloride were observed possibly because the protonation state remained unchanged.
For Compound 2, the strong intensity bands in the IR spectrum at 1,624 and 1,569 cm−1 could originate from v(C = N). The medium-intensity bands in the IR at 1,481 and 1,419 cm−1 could arise from asymmetric deformations of CH3 groups (Gunasekaran et al., 2006; Sheela et al., 2010; Ghasemi et al., 2018).
The strong IR vibration at 950 cm−1 corresponds to v(V = O). The corresponding vibrations appear in the Raman spectrum (Figure 6D) at 977 and 952 cm−1. In the same Raman spectrum, an intense peak is observed at 938 cm−1, which can be assigned to δ(N–H) with a contribution of vas(V–O–V) (Frost et al., 2005; Gunasekaran et al., 2006). The IR absorption at 830 and 736 cm−1 corresponds to the vas(V–O–V) bond. The Raman line at 833 cm−1 originated from the same vibration is very weak. Returning to the IR absorption spectrum, medium-intensity bands are observed at 603 and 571 cm−1, which correspond to vs(V–O–V). This band is observed also at 602 cm−1 in the Raman spectrum (Frost et al., 2005; Frost and Palmer, 2011; Pavliuk et al., 2014; Sánchez-Lara et al., 2016b). Vibrations below 450 cm−1 are quite similar to those found for Compound 1 crystallized with cytosine.
Solution NMR-Spectroscopy
It has been shown that vanadate(V) species in solution can occur simultaneously in equilibrium with a different state of protonation and, in some cases, with different conformations (Aureliano and Ohlin, 2014). These characteristics of the aqueous vanadium(V) chemistry are highly dependent on pH, vanadium concentration and ionic strength (Dorsey et al., 2018). Decavanadate structure in solution, for example, could be decomposed into three species: monomeric [H2VO4]−, dimeric [H2V2O7]2− and tetrameric [V4O12]4− vanadates. This speciation has been studied mainly by 51V-NMR spectroscopy and these findings have significant consequences for toxicology activities and pharmacological applications of decavanadate based-compounds (Soares et al., 2007; Aureliano et al., 2016). Due to this last, the 51V-NMR spectroscopy for both compounds was carried out. Aqueous solutions of Compound 1 at pH close to 7 display the classical 51V resonance signals ascribed to the three different vanadium atoms of the decavanadate structure, present at: V10A = −515.88 ppm, V10B = −500.31 ppm, and V10C = −423.55 ppm (Figure 7) (Rehder et al., 2007; Rehder, 2015a,b). However, it is known that decavanadate slowly decomposes when the pH > 6.5, where it is thermodynamically unstable and it is transformed into the labile monomers and cyclic vanadates (Soares et al., 2007; Aureliano and Crans, 2009). Thus, the sharp signal at δ = −560.04 ppm can be assigned to the biologically active species [H2VO4]− (V1), while the signals at δ = −572.52 and −577.06 ppm can be attributed to the divanadate [H2V2O7]2− (V2) and cyclic tetramer [V4O12]4− (V4), respectively (Rehder et al., 2007).
Figure 7. 51V-NMR spectra in D2O of 1-2. The signals correspond to sites VA (low field), VB and VC (high field), and the oligomeric vanadium species V1, V2, and V4. The inset displays the three different vanadium sites found in decavanadate.
For Compound 2, the shifts for [V10O28]6− structure are very similar to those for Compound 1, which are observed at: V10A = −514.43 ppm, V10B = −500, V10C = −422.31 ppm. We also observe the presence of other vanadate species as a consequence of decavanadate hydrolysis, with signals at −560.39, −572.19, and −576.80 ppm, which correspond to V1, V2, and V4, respectively (see Figure 7).
On the other hand, the 13C-NMR spectrum (Figure S1) shows one signal at 37.41 ppm corresponding to the equivalent methyl groups of the Metformin molecules, and the two signals observed at 158.43 and 160.12 ppm corresponding to the tertiary guanidyl carbons C–NH2 and C = , respectively (Chatkon et al., 2013; Ibrahim et al., 2015). The 1H-NMR in water-D2 shows only one resonance peak at 2.7 ppm due to two equivalent methyl groups (data not shown). The N–H protons were exchanged in D2O by eliminating the signals at 7.20 and 7.69 ppm, previously observed in DMSO-d6 (Gadape and Parikh, 2011).
DSC-TGA Analysis
The thermal behavior of [HCyt]6[V10O28]·4H2O (1), and [HMetf]6[V10O28]·6H2O (2), was assessed by TGA/DTA analysis (Figures S2, S3). For Compound 1, the thermogram shows a mass loss of around 3-4% corresponding to lattice water molecules in the range of 30–150°C, which fits with four water molecules [% mass, calc. (found): 4.23% (4%)]. An endothermic process occurs around 225°C which may be attributed to the melting of cytosine, followed by a continuous mass loss in the range from 220 to 580°C which may correspond to the thermal degradation of the six cytosinium ions [% mass, calc. (found) for 6 × C4H6N3O: 39.5% (40%)]. For Compound 2, the DTA thermogram shows three well defined endothermic peaks in a temperature range of 37–70°C, all related to the loss of lattice water molecules. This water loss is in agreement with the corresponding mass change observed in the TGA curve suggesting the loss of six uncoordinated water molecules [% mass, calc (found): 5.85% (6%)]. Completion of this stage leads to the anhydrous phase, which shows a significant range of thermal stability extending to ca. 200°C. This thermal behavior has been observed for other polyoxometalates ions with organic moieties (Iyer et al., 2014; Sánchez-Lombardo et al., 2014; Martin-Caballero et al., 2016; Dissem et al., 2018). Considering that the melting point of Metformin hydrochloride is 224°C (Benmessaoud et al., 2016), we can assign the endothermic peaks observed in the DTA thermogram in a temperature range of 190–210°C, to the melting of Metforminium cations (HMetf) present in Compound 2. Above this temperature, an exothermic mass loss takes place that corresponds to the combustion of the metforminium cations per decavanadate unit [%mass, calc. (found) for 6 × C4H12N5: 42.3% (40%)]. In both cases, above 540–580°C no further mass losses are observed, which can be attributed to the remaining inorganic fragment [V10O28]6− which is thermally stable event up to 600°C as observed in previously reported decavanadates compounds (Omri et al., 2015; Ortaboy et al., 2018). It may be possible that further thermal treatment of the decavanadate ion may yield other vanadium oxides, possibly V2O5 or other mixed-valence oxides (Riou et al., 1998). However, a more in-depth analysis such as thermodiffractometry is necessary to confirm the nature of the resulting metal oxides, which may probably consist of vanadium oxide nanoparticles produced by thermal decomposition of the polyoxovanadate systems (Martin-Caballero et al., 2016).
Theoretical Results
Optimized structures of the asymmetric unit in aqueous solvation phase of Compounds 1 and 2 are shown in Figures 8, 9. It can be seen that the theoretical structures are similar to those mapped from crystallographic data. Selected parameters for the compounds are shown in Tables S1, S2. Selected hydrogen bonds are also reported in Table S3. Atom labels are the same as those used in Figures 1, 3.
Figure 8. Solution-phase of the optimized structure of the asymmetric unit of Compound 1, using the B97-D3 functional with the 6-31G(d) basis set for the C, H, O, N, atoms and the LanL2MB basis set for the V atoms.
Figure 9. Solution-phase optimized structure of the asymmetric unit of Compound 2 using the B97-D3 functional with the 6-31G(d) basis set for the C, H, O, N, atoms and the LanL2MB basis set for V atoms.
Table S1 shows the structural parameters for Compound 1 obtained by the DFT-D methods in solution phase from X-ray crystallographic data used as an initial structure for calculation. In the aqueous phase, some variations were observed in the dihedral angle values calculated in the decavanadate anion [V10O28]6− due to the solvent effect. The group of cytosinium cations ([C4H6N3O]+) interacting with the decavanadate anion [V10O28]6− replicates the behavior in the crystallographic unit cell. Water molecules in Compound 1 form important hydrogen bond interactions.
For Compound 2, the full optimization of the asymmetric unit showed that the intramolecular parameters for the metforminium cations ([C4H12N5]+) had significant structural displacements relative to the position of the decavanadate. The decavanadate ion optimized parameters for Compound 2 were obtained with good precision using the crystallographic data as a reference. Moreover, the optimized parameters for the anion structure were similar in both compounds. Contrary to what is observed in 1, where one of the water molecules acts as a water bridge (μ2-O) between the cation and the decavanadate, the optimized water molecules in Compound 2 act as a cluster (O51–O55) without any apparent cation-anion interaction (Figure 9).
Including the solvent effect in both compounds was a decisive step for geometry optimization, as well as taking into consideration the non-covalent interactions N–H···O and O–H···O through the inclusion of the dispersion forces by the DFT-D methods.
The main hydrogen bonds were also theoretically characterized. A previous study about the electronic properties of a cytosine-decavanadate system was carried out by using Atoms in Molecules Topological Analysis (Bošnjakovic-Pavlovic et al., 2009). In this work, the density of D–H···A hydrogen bonds were analyzed, suggesting the N–H···O as strong hydrogen bonds and C–H···O as weak hydrogen bonds in the cytosine-decavanadate system containing sodium cations. In the present work, the non-covalent interactions N–H···O were also found using the DFT-D methodology and including a solvent effect by means of the CPCM approach. Values similar to crystallographic data, used as a reference, were obtained for interactions N–H···O for Compounds 1 and 2 as reported in Table 2. In addition, important interactions were calculated between water molecules in both compounds (Table S3).
Figures 10, 11 show the distribution for the frontier molecular orbitals [the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO)], and the molecular electrostatic potential (MEP) for Compounds 1 and 2. When the decavanadate anion [V10O28]6− is built together with cytosinium or metforminium cations and water molecules, significant symmetric non-covalent interactions are observed in the resulting compounds. For compound 1, the HOMO was mainly delocalized on all oxygen atoms of [V10O28]6− except O1, as shown in Figure 8. Furthermore, the nitrogen atoms and the partially delocalized C = C double bonds of Cytosinium cations provide a significant contribution to HOMO. On the other hand, LUMO was mainly delocalized on the nitrogen and carbon atoms of cations [C4H6N3O]+, (see Figures 10A,B). For compound 2 the HOMO distribution is mainly delocalized on the oxygen atoms of [V10O28]6−, with contributions from the O1, O3 and O5 atoms (and their symmetric equivalents in decavanadate anion), see Figure 9. Additionally, a significant contribution to HOMO is provided from the atoms forming C–N bonds in Metforminium cations. LUMO is mainly delocalized on nitrogen atoms of cations [C4H12N5]+ (see Figures 11A,B). The total electron density was mapped with the electrostatic potential surface (isovalue = 0.004) for both compounds, as shown in Figures 10C, 11C. The qualitative color code indicates red regions with a negative charge, while blue regions indicate positive charge. Yellow and green regions correspond to intermediate values toward negative or positive charges, respectively. The negative charge is gathered mainly on the anion [V10O28]6−, while counterions [C4H6N3O]+, and [C4H12N5]+ have a positive charge. Non-covalent interactions N–H···O and O–H···O are located in intermediate regions of electron density. The distribution of charges in both compounds clearly indicate the relevance of Coulombic interactions.
Figure 10. Isosurfaces of frontier molecular orbitals (A) HOMO and (B) LUMO, and (C) MEP of Compound 1 at B97-D3 functional with the 6-31G(d) basis set for C, H, O, N, atoms and LanL2MB basis set for V atoms in solution phase.
Figure 11. Isosurfaces of frontier molecular orbitals (A) HOMO and (B) LUMO, and (C) MEP of Compound 2 at B97-D3 functional with the 6-31G(d) basis set for C, H, O, N, atoms and LanL2MB basis set for V atoms in solution phase.
Conclusions
This paper focused on the structural description and characterization of two new decavanadate salts, obtained at acidic pH under mild conditions by a combination of ammonium and sodium metavanadates with pharmacologically active counterions from Cytosine and Metformin. In summary, the preparation and complete characterization by elemental analysis, XRD single crystal analysis, spectroscopy and thermal analysis of these new compounds was reported. Also, theoretical studies were carried out to analyze the molecular structure and characterize the intermolecular interactions of these hybrid compounds.
About the crystal structures, we observed that molecular packing in these compounds is predominantly driven by strong hydrogen bonds, as well as Coulombic interactions between cations and decavanadate cluster, building interesting supramolecular networks that have been described herein. The solid-state spectroscopy analyses (FT-IR and Raman) confirmed the presence of [V10O28]6− and the characteristic vibrations of the counterions. The terminal V = O bonds occur in the 900–1,000 cm−1 region; bridging V–O–V bonds vibrate in the 500–700 cm−1 region as symmetric and asymmetric stretching modes, respectively, while the V–O bending modes occur in the 300–400 cm−1 region. The 51V-NMR experiments showed that the equilibrium between decavanadate and vanadate species is progressively displaced toward the formation of V1, V2, and V4 species.
On the other hand, theoretical characterization of optimized structures, frontier molecular orbitals, and the electrostatic potential distribution were obtained. The main non-covalent interactions N–H···O and O–H···O were characterized using DFT-D methods, including dispersion forces that are essential for the treatment of this kind of complexes. The inclusion of solvent effect showed interesting results for correlating the electronic and electrostatic properties with the optimized structures, regardless of those obtained for the crystalline environment.
Currently, further calculations of the theoretical spectroscopic characterization of the IR, Raman, and 51V-NMR spectra, as well as a topological analysis of the complexes, are being carried out. Although the scope of this work is far from a biological analysis, we consider that the reported structures can provide interesting insights into their potential pharmacological effects, due to biological activities of the counterions used. First, Compound 2 combined with Metformin shows higher water solubility, becoming promising candidates for carrying out biological studies. Second, it has been reported that decavanadate has anticancer and antidiabetic properties and Metformin, widely recognized as one of the safest and most effective therapeutic compounds for the treatment of Type 2 D. mellitus is currently used in the treatment of different types of cancer. Also, deoxycytidine, cytosine arabinoside, and gemcitabine.6 which are derivatives of Cytosine are widely used for the treatment of various malignancies. Benchmark drugs are cytarabine for acute myeloid leukemia and gemcitabine for pancreatic and lung cancer (Muggia et al., 2012).
Biological studies are currently under investigation in our laboratory. However, preliminary results in a mouse colon adenocarcinoma model have shown similar IC50 value as CISPLATIN for Cytosinium decavanadate (Sánchez-Lara et al., 2016a). As for the Metforminium decavanadate 3:1, promising results have been obtained in type 1 and 2 murine models of D. mellitus. (Treviño et al., 2015, 2016, 2018). The combination of two therapeutic agents opens up a window toward the generation of potential metalopharmaceuticals with new and exciting pharmacological properties. The results of those investigations will be reported elsewhere in due course.
Author Contributions
All authors listed have made a substantial, direct and intellectual contribution to the work, and approved it for publication.
Conflict of Interest Statement
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
ES-L thanks CONACYT (Mexico) for the Ph.D. Scholarship (Support No. 293256). We thank Prof. S. Bernès for the structures solution and refinement, as well as for his comments to the final version of the draft. MEC and FJM-B would like to thank the Laboratorio Nacional de Supercómputo del Sureste de México (LNS-BUAP) of the CONACYT network of national laboratories, for the computer resources and support provided. FT-Raman measurements were carried out by Dra. L. Serrano, Laboratorio Central del Instituto de Física Luis Rivera Terrazas (IFUAP), Benemérita Universidad Autónoma de Puebla. TGA-DTA measurements were carried out by Rubi E. Gonzalez-Campos, UDLAP. The authors wish to thank Vicerrectoria de Investigación y Estudios de Posgrado (VIEP-BUAP) for the financial support obtained through projects 100517029-VIEP2018 and 100256733-VIEP2018.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fchem.2018.00402/full#supplementary-material
References
Aureliano, M. (2011). Recent perspectives into the biochemistry of decavanadate. World J. Biol. Chem. 2:215. doi: 10.4331/wjbc.v2.i10.215
Aureliano, M. (2016). Decavanadate toxicology and Pharmacological Activities: V10 or V1, Both or None? Oxid. Med. Cell Longev. 2016:6103457. doi: 10.1155/2016/6103457
Aureliano, M. (2017). The role of decavanadate in anti-tumor activity. Glob. J. Cancer Ther. 3, 12–14. doi: 10.17352/gjct.000015
Aureliano, M., and Crans, D. C. (2009). Decavanadate (V10O286-) and oxovanadates: oxometalates with many biological activities. J. Inorg. Biochem. 103, 536–546. doi: 10.1016/j.jinorgbio.2008.11.010
Aureliano, M., and Ohlin, C. A. (2014). Decavanadate in vitro and in vivo effects: facts and opinions. J. Inorg. Biochem. 137, 123–130. doi: 10.1016/j.jinorgbio.2014.05.002
Aureliano, M., Ohlin, C. A., Meira, M. O., Marques, M. P. M., Casey, W. H., and De Carvalho, L. A. B. (2016). Characterization of decavanadate and decaniobate solutions by Raman spectroscopy. Dalton Trans. 45, 7391–7399. doi: 10.1039/C5DT04176G
Ayed, M., Ayed, B., and Haddad, A. (2015). Synthesis and structural characterization of new inorganic–organic hybrid: arsenomolybdate compound with cytosinium cations. Bull. Mater. Sci. 38, 13–21. doi: 10.1007/s12034-014-0805-8
Barone, V., Cossi, M., and Tomasi, J. (1998). Geometry optimization of molecular structures in solution by the polarizable continuum model. J. Comput. Chem. 19, 404–417. doi: 10.1002/(SICI)1096-987X(199803)19:4 < 404::AID-JCC3>3.0.CO;2-W
Bartošová, L., Padělková, Z., Rakovský, E., and Schwendt, P. (2012). Synthesis and crystal structure of two copper (II) complexes with coordinated decavanadate ion. Polyhedron 31, 565–569. doi: 10.1016/j.poly.2011.10.042
Benmessaoud, I., Koutchoukali, O., Bouhelassa, M., Nouar, A., and Veesler, S. (2016). Solvent screening and crystal habit of metformin hydrochloride. J. Cryst. Growth. 451, 42–51. doi: 10.1016/j.jcrysgro.2016.07.001
Bijelic, A., Aureliano, M., and Rompel, A. (2018). Polyoxometalates as potential next-generation metallodrugs in the combat against cancer. Angewante Chem. Int Ed. doi: 10.1002/anie.201803868. [Epub ahead of print].
Bošnjakovic-Pavlovic, N., Spasojevic-de Biré, A., Tomaz, I., Bouhmaida, N., Avecilla, F., Mioc, U. B., et al. (2009). Electronic properties of a cytosine decavanadate: toward a better understanding of chemical and biological properties of decavanadates. Inorg. Chem. 48, 9742–9753. doi: 10.1021/ic9008575
Cafferty, B. J., Fialho, D. M., Khanam, J., Krishnamurthy, R., and Hud, N. V. (2016). Spontaneous formation and base pairing of plausible prebiotic nucleotides in water. Nat. Commun. 7:11328. doi: 10.1038/ncomms11328
Cantley, L. C. Jr., Josephson, L., Warner, R., Yanagisawa, M., Lechene, C., and Guidotti, G. (1977). Vanadate is a potent (Na,K)-ATPase inhibitor found in ATP derived from muscle. J Biol Chem. 252, 7421–7423.
Chatkon, A., Barres, A., Samart, N., Boyle, S. E., Haller, K. J., and Crans, D. C. (2014). Guanylurea metformium double salt of decavanadate,(HGU+)4(HMet+)2(V10O286–)·2H2O. Inorg. Chim. Acta 420, 85–91. doi: 10.1016/j.ica.2013.12.031
Chatkon, A., Chatterjee, P. B., Sedgwick, M. A., Haller, K. J., and Crans, D. C. (2013). Counterion affects interaction with interfaces: the antidiabetic drugs metformin and decavanadate. Eur. J. Inorg. Chem. 2013, 1859–1868. doi: 10.1002/ejic.201201345
Childs, S. L., Chyall, L. J., Dunlap, J. T., Coates, D. A., Stahly, B. C., and Stahly, G. P. (2004). A metastable polymorph of metformin hydrochloride: isolation and characterization using capillary crystallization and thermal microscopy techniques. Cryst. Growth Des. 4, 441–449. doi: 10.1021/cg034243p
Chinea, E., Dakternieks, D., Duthie, A., Ghilardi, C. A., Gili, P., Mederos, A., et al. (2000). Synthesis and characterization of (NH4)4[H4V10O28][NTPH2]2·4H2O (NTPH3 = N(CH2CH2COOH)3). Inorg. Chim. Acta 298, 172–177. doi: 10.1016/S0020-1693(99)00440-5
Correia, I., Avecilla, F., Marcão, S., and Costa-Pessoa, J. (2004). Structural studies of decavanadate compounds with organic molecules and inorganic ions in their crystal packing. Inorg. Chim. Acta 357, 4476–4487. doi: 10.1016/j.ica.2004.06.055
Costa-Pessoa, J., Etcheverry, S., and Gambino, D. (2015). adium compounds in medicine. Coord. Chem. Rev. 301, 24–48. doi: 10.1016/j.ccr.2014.12.002
Crans, D. C. (2000). Chemistry and insulin-like properties of vanadium(IV) and vanadium(V) compounds. J. Inorg. Biochem. 80, 123–131. doi: 10.1016/S0162-0134(00)00048-9
Crans, D. C. (2015). Antidiabetic, chemical, and physical properties of organic vanadates as presumed transition-state inhibitors for phosphatases. J. Org. Chem. 80, 11899–11915. doi: 10.1021/acs.joc.5b02229
Del Carpio, E., Hernández, L., Ciangherotti, C., Villalobos Coa, V., Jiménez, L., Lubes, V., et al. (2018). Vanadium: history, chemistry, interactions with α-amino acids and potential therapeutic applications. Coord. Chem. Rev. 372, 117–140. doi: 10.1016/j.ccr.2018.06.002
Dennington, R. D. II., Keith, T. A., and Millam, J. M. (2016). Gauss View, Version 6.0.16. Shawnee Mission: Semichem Inc.
Dissem, N., Artetxe, B., San Felices, L., Lezama, L., Haddad, A., and Guitierrez-Zorrilla, J. M. (2018). A robust framework based on polymeric octamolybdate anions and copper (ii) complexes of tetradentate N-donor ligands. Crystals 8:20. doi: 10.3390/cryst8010020
Dolomanov, O. V., Bourhis, L. J., Gildea, R. J., Howard, J. A., and Puschmann, H. (2009). OLEX2: a complete structure solution, refinement, and analysis program. J. Appl. Crys. 42, 339–341. doi: 10.1107/S0021889808042726
Dong, J., Liu, B., and Yang, B. (2015). Crystal structure of bis {2-[amino (iminiumyl) methyl]-1, 1-dimethylguanidine} carbonate methanol disolvate. Acta Crys. E 71, o747–o748. doi: 10.1107/S2056989015016771
Dorsey, B. M., McLauchlan, C. C., and Jones, M. A. (2018). Evidence that speciation of oxovanadium complexes does not solely account for inhibition of leishmania acid phosphatases. Front. Chemi. 6:109. doi: 10.3389/fchem.2018.00109
Elberg, G., He, Z., Li, J., Sekar, N., and Shechter, Y. (1997). Vanadate activates membranous nonreceptor protein tyrosine kinase in rat adipocytes. Diabetes 46:16841690.
El-Mir, M. Y., Nogueira, V., Fontaine, E., Avéret, N., Rigoulet, M., and Leverve, X. (2000). Dimethylbiguanide inhibits cell respiration via an indirect effect targeted on the respiratory chain complex I. J. Biol. Chem. 275, 223–228. doi: 10.1074/jbc.275.1.223
Etter, M. C., MacDonald, J. C., and Bernstein, J. (1990). Graph-set analysis of hydrogen-bond patterns in organic crystals. Acta Crys. B 46, 256–262. doi: 10.1107/S0108768189012929
Evans, J. M., Donnelly, L. A., Emslie-Smith, A. M., Alessi, D. R., and Morris, A. D. (2005). Metformin and reduced risk of cancer in diabetic patients. Brit. Med. J. 330, 1304–1305. doi: 10.1136/bmj.38415.708634.F7
Frisch, M. J., Trucks, G. W., Schlegel, H. B., Scuseria, G. E., Robb, M. A., Cheeseman, J. R., et al. (2016). Gaussian 16, Revision B.01. Wallingford, CT: Gaussian, Inc.
Frost, R. L., Erickson, K. L., Weier, M. L., and Carmody, O. (2005). Raman and infrared spectroscopy of selected vanadates. Spectrochim. Acta A 61, 829–834. doi: 10.1016/j.saa.2004.06.006
Frost, R. L., and Palmer, S. J. (2011). Raman spectroscopic study of pascoite Ca3V10O28·17H2O. Spectrochim. Acta A 78, 248–252. doi: 10.1016/j.saa.2010.10.002
Gadape, H. H., and Parikh, K. S. (2011). Quantitative determination and validation of metformin hydrochloride in pharmaceutical using quantitative nuclear magnetic resonance spectroscopy. J. Chem. 8, 767–781. doi: 10.1155/2011/461898
Ghasemi, F., Rezvani, A. R., Ghasemi, K., and Graiff, C. (2018). Glycine and metformin as new counter ions for mono and dinuclear vanadium (V)-dipicolinic acid complexes based on the insulin-enhancing anions: synthesis, spectroscopic characterization, and crystal structure. J. Mol. Struct. 1154, 319–326. doi: 10.1016/j.molstruc.2017.10.055
Grimme, S., Ehrlich, S., and Goerigk, L. (2011). Effect of the damping function in dispersion corrected density functional theory. J. Comput. Chem. 32, 1456–1465. doi: 10.1002/jcc.21759
Groom, C. R., Bruno, I. J., Lightfoot, M. P., and Ward, S. C. (2016). The Cambridge structural database. Acta Cryst. B 72, 171–179. doi: 10.1107/S2052520616003954
Gunasekaran, S., Natarajan, R. K., Renganayaki, V., and Natarajan, S. (2006). Vibrational Spectra and Thermodynamic Analysis of Metformin. Available online at: http://nopr.niscair.res.in/handle/123456789/8327
Hay, P. J., and Wadt, W. R. (1985). Ab initio effective core potentials for molecular calculations. Potentials for the transition metal atoms Sc to Hg. J. Chem. Phys. 82, 270–283. doi: 10.1063/1.448799
Hayashi, Y. (2011). Hetero and lacunary polyoxovanadate chemistry: synthesis, reactivity, and structural aspects. Coord. Chem. Rev. 255, 2270–2280. doi: 10.1016/j.ccr.2011.02.013
Hohenberg, P., and Kohn, W. (1964). Inhomogeneous electron gas. Phys. Rev. 136:B864. doi: 10.1103/PhysRev.136.B864
Ibrahim, M. M., Mohamed, M. A., Mersal, G. A., and Al-Juaid, S. (2015). Insulin-like action of novel metformin-containing vanadate as a new antidiabetic drug: synthesis, characterization, and crystal structure of [Metformin-H]2[V2O6]·H2O. J. Mol. Struct. 1098, 92–100. doi: 10.1016/j.molstruc.2015.05.050
Iyer, A. K., Roy, S., Haridasan, R., Sarkar, S., and Peter, S. C. (2014). Ligand-mediated valence fultuation of copper in new hybrid materials constructed from decavanadate and a Cu(1,10-phenanthroline) complex. Dalton Trans. 43, 2153–2160. doi: 10.1039/C3DT52484A
Kasznicki, J., Sliwinska, A., and Drzewoski, J. (2014). Metformin in cancer prevention and therapy. Ann. Transl. Med. 2:57. doi: 10.3978/j.issn.2305-5839.2014.06.01
Kathuria, D., Bankar, A. A., and Bharatam, P. V. (2018). “What's in a structure?” The story of biguanides. J. Mol. Struct. 1152, 61–78. doi: 10.1016/j.molstruc.2017.08.100
Kioseoglou, E., Gabriel, C., Petanidis, S., Psycharis, V., Raptopoulou, C. P., Terzis, A., et al. (2013). Binary decavanadate-betaine composite materials of potential anticarcinogenic activity. Z. Anorg. Allg. Chem. 639, 1407–1416. doi: 10.1002/zaac.201300144
Kojima, T., Antonio, M. R., and Ozeki, T. (2011). Solvent-driven association and dissociation of the hydrogen-bonded protonated decavanadates. J. Am. Chem. Soc. 133, 7248–7251. doi: 10.1021/ja110781f
Lee, T., and Wang, P. Y. (2010). Screening, manufacturing, photoluminescence, and molecular recognition of co-crystals: cytosine with dicarboxylic acids. Cryst. Growth Des. 10, 1419–1434. doi: 10.1021/cg901412d
Lee, U. K., and Joo, H. C. (2003). Potassium–sodium double salt of decavanadate, K4Na2 [V10O28]·10H2O. Acta Crys. E 59, i122–i124. doi: 10.1107/S1600536803016453
Libby, G., Donnelly, L. A., Donnan, P. T., Alessi, D. R., Morris, A. D., and Evans, J. M. (2009). New users of metformin are at low risk of incident cancer: a cohort study among people with type 2 diabetes. Diabetes care 32, 1620–1625. doi: 10.2337/dc08-2175
Macrae, C. F., Bruno, I. J., Chisholm, J. A., Edgington, P. R., McCabe, P., Pidcock, E., et al. (2008). Mercury CSD 2.0 - new features for the visualization and investigation of crystal structures. J. Appl. Cryst. 41, 466–470. doi: 10.1107/S0021889807067908
Madzharova, F., Heiner, Z., Gühlke, M., and Kneipp, J. (2016). Surface-enhanced hyper-Raman spectra of adenine, guanine, cytosine, thymine, and uracil. J. Phys. Chem. 120, 15415–15423. doi: 10.1021/acs.jpcc.6b02753
Mal, S. S., Tröppner, O., Ivanović-Burmazović, I., and Burger, P. (2013). Tetraalkylphosphonium decavanadates: synthesis, structures, and solution properties. Eur. J. Inorg. Chem. 2013, 1960–1967. doi: 10.1002/ejic.201201447
Martín-Caballero, J., San José Wéry, A., Reinoso, S., Artetxe, B., San Felices, L., Bakkali, E. I., et al. (2016). A robust open framework formed by decavanadate clusters and copper(ii) complexes of macrocyclic polyamines: permanent microporosity and catalytic oxidation of cycloalkanes. Inorg. Chem. 55, 4970–4979. doi: 10.1021/acs.inorgchem.6b00505
Mathlouthi, M., Seuvre, A. M., and Koenig, J. L. (1986). FT-IR and laser-Raman spectra of cytosine and cytidine. Carbohydr. Res. 146, 1–13. doi: 10.1016/0008-6215(86)85019-4
Muggia, F., Diaz, I., and Peters, G. J. (2012). Nucleoside and nucleobase analogs in cancer treatment: not only sapacitabine but also gemcitabine. Expert Opin. Investig. Drugs. 21, 403–4088. doi: 10.1517/13543784.2012.666236
Nagazi, I., and Haddad, A. (2014). Synthesis, crystal structure and physicochemical proprieties of a cytosine selenomolybdate based on strandberg-type [Se2Mo5O21]4− polyanion. J. Cluster Sci. 25, 627–638. doi: 10.1007/s10876-013-0658-0
Narla, R. K., Chen, C. L., Dong, Y., and Uckun, F. M. (2001). In vivo antitumor activity of bis(4,7-dimethyl-1,10-phenanthroline) sulfatooxovanadium(IV) (METVAN [VO(SO4)(Me2-Phen)2]). Clin Cancer Res. 7, 2124–2133.
Olar, R., Badea, M., Marinescu, D., Chifiriuc, C. M., Bleotu, C., Grecu, M. N., et al. (2010). Prospects for new antimicrobials based on N, N-dimethylbiguanide complexes as effective agents on both planktonic and adhered to microbial strains. Eur. J. Med. Chem. 45, 2868–2875. doi: 10.1016/j.ejmech.2010.03.009
Omri, I., Mhiri, T., and Graia, M. (2015). Novel decavanadate cluster complexes (HImz)12(V10O28)·2.3H2O: synthesis, characterization, crystal structure, optical and thermal properties. J. Mol. Struct. 1098, 324–331. doi: 10.1016/j.molstruc.2015.06.011
Ortaboy, S., Acar, E. T., and Atun, G. (2018). The removal of radioactive strontium ions from aqueous solutions by isotopic exchange using a strontium decavanadate and corresponding mixed oxides. Chem. Eng. J. 344, 194–205. doi: 10.1016/j.cej.2018.03.069
Orvig, C., Caravan, P., Gelmini, L., Glover, N., Herring, F. G., Li, H., et al. (1995). Reaction chemistry of BMOV, bis (maltolato) oxovanadium (IV), a potent insulin mimetic agent. J. Am. Chem. Soc. 117, 12759–12770 doi: 10.1021/ja00156a013
Parker, W. B. (2009). Enzymology of purine and pyrimidine antimetabolites used in the treatment of cancer. Chem. Rev. 109, 2880–2893. doi: 10.1021/cr900028p
Pavliuk, M. V., Makhankova, V. G., Khavryuchenko, O. V., Kokozay, V. N., Omelchenko, I. V., Shishkin, O. V., et al. (2014). Decavanadates decorated with [Cu(en)2]2+: convenient synthetic route, crystal structures and analysis of vibrational spectra. Polyhedron 81, 597–606. doi: 10.1016/j.poly.2014.06.044
Pereira, M. J., Carvalho, E., Eriksson, J. W., Crans, D. C., and Aureliano, M. (2009). Effects of decavanadate and insulin enhancing vanadium compounds on glucose uptake in isolated rat adipocytes. J. Inorg. Biochem. 103, 1687–1692. doi: 10.1016/j.jinorgbio.2009.09.015
Pérez-Fernández, R., Fresno, N., Goya, P., Elguero, J., Menéndez-Taboada, L., García-Granda, S., et al. (2013). Structure and thermodynamical properties of metformin salicylate. Cryst. Growth Des. 13, 1780–1785. doi: 10.1021/cg400150f
Perumalla, S. R., Pedireddi, V. R., and Sun, C. C. (2013). Protonation of cytosine: cytosinium vs. hemicytosinium duplexes. Cryst. Growth Des. 13, 429–432. doi: 10.1021/cg3014915
Rasheed, T., and Ahmad, S. (2010). Approximate solution of the mode–mode coupling integral: application to cytosine and its deuterated derivative. Spectrochim. Acta A 77, 446–456. doi: 10.1016/j.saa.2010.06.015
Rassolov, V. A., Pople, J. A., Ratner, M. A., and Windus, T. L. (1998). 6-31G* basis set for atoms K through Zn. J. Chem. Phys. 109, 1223–1229. doi: 10.1063/1.476673
Rehder, D. (2011). A possible role for extraterrestrial vanadium in the encounter of life. Coord. Chem. Rev. 255, 2227–2231. doi: 10.1016/j.ccr.2011.04.015
Rehder, D. (2012). The potentiality of vanadium in medicinal applications. Future Med. Chem. 4, 1823–1837. doi: 10.4155/fmc.12.103
Rehder, D. (2015a). The (biological) speciation of vanadate (V) as revealed by 51V NMR: a tribute on Lage Pettersson and his work. J. Inorg. Biochem. 147, 25–31. doi: 10.1016/j.jinorgbio.2014.12.014
Rehder, D. (2015b). The role of Vanadium in biology. Metallomics 7, 730–742. doi: 10.1039/C4MT00304G
Rehder, D., Polenova, T., and Bühl, M. (2007). Vanadium-51 NMR. Annu. Rep. NMR Spectrosc. 62, 49–114. doi: 10.1016/S0066-4103(07)62002-X
Repišcák, P., Erhardt, S., Rena, G., and Paterson, M. J. (2014). Biomolecular mode of action of metformin in relation to its copper binding properties. Biochemistry 53, 787–795. doi: 10.1021/bi401444n
Riou, D., Roubeau, O., and Férey, G. (1998). Evidence for the solid state structural transformation of the network-type decavanadate (NC7H14)4[H2V10O28] into a lamellar topology (NC7H14)[V4O10]. Z. Anorg. Allg. Chem. 624, 1021–1025. doi: 10.1002/(SICI)1521-3749(199806)624:63.0.CO;2-K
Rozenberg, M., Shoham, G., Reva, I., and Fausto, R. (2004). Low-temperature FTIR spectroscopy and hydrogen bonding in cytosine polycrystals. Spectrochim. Acta A 60, 463–470. doi: 10.1016/S1386-1425(03)00251-8
Sánchez-Lara, E., Pérez-Benítez, A., Treviño, S., Bernès, S., Delgado, R., Sánchez-Mora, E., et al. (2016a). “Metforminium decavanadates as a potential metallodrugs in the treatment of cancer and diabetes,” in V Latinamerican Metting on Biological Inorganic Chemistry. Santiago de Querétaro.
Sánchez-Lara, E., Pérez-Benítez, A., Treviño, S., Mendoza, A., Meléndez, F. J., Sánchez-Mora, E., et al. (2016b). Synthesis and 3D network architecture of 1-and 16-hydrated salts of 4-dimethylaminopyridinium decavanadate,(DMAPH)6[V10O28]· nH2O. Crystals 6:65. doi: 10.3390/cryst6060065
Sánchez-Lara, E., Sanchez-Lombardo, I., Pérez-Benitez, A., Mendoza, A., Flores-Alamo, M., and Gonzalez-Vergara, E. (2015). A new dicationic ring [(Water)6–(Ammonium)2] acts as a building block for a supramolecular 3D assembly of decavanadate clusters and 4-(N,Ndimethylamino) pyridinium Ions. J. Cluster Sci. 26, 901–912. doi: 10.1007/s10876-014-0779-0
Sánchez-Lombardo, I., Sánchez-Lara, E., Pérez-Benítez, A., Mendoza, Á., Bernès, S., and González-Vergara, E. (2014). Synthesis of metforminium (2+) decavanadates–crystal structures and solid-state characterization. Eur. J. Inorg. Chem. 2014, 4581–4588. doi: 10.1002/ejic.201402277
Selman, M., Rousso, C., Bergeron, A., Son, H. H., Krishnan, R., El-Sayes, N. A., et al. (2018). Multi-modal potentiation of oncolytic virotherapy by vanadium compounds. Mol. Ther. 26, 55–69. doi: 10.1016/j.ymthe.2017.10.014
Serb, M. D., Kalf, I., and Englert, U. (2014). Biguanide and squaric acid as pH-dependent building blocks in crystal engineering. Cryst. Eng. Comm. 16, 10631–10639. doi: 10.1039/C4CE01643B
Sheela, N. R., Muthu, S., and Krishnan, S. S. (2010). FTIR, FT Raman and UV-visible spectroscopic analysis on metformin hydrochloride. Asian J. Chem. 22, 5049–5056.
Sheldrick, G. M. (2015). Crystal structure refinement with SHELXL. Acta Cryst. C 71, 3–8. doi: 10.1107/S2053229614024218
Soares, S. S., Martins, H., Duarte, R. O., Moura, J. J., Coucelo, J., Gutiérrez-Merino, C., et al. (2007). Vanadium distribution, lipid peroxidation and oxidative stress markers upon decavanadate in vivo admistration. J. Inorg. Biochem. 101 80–88. doi: 10.1016/j.jinorgbio.2005.02.023
Sridhar, B., Nanubolu, J. B., and Ravikumar, K. (2012). Competition between the two-point and three-point synthon in cytosine–carboxylic acid complexes. Cryst. Eng. Comm. 14, 7065–7074. doi: 10.1039/c2ce26076j
Steens, N., Ramadan, A. M., Absillis, G., and Parac-Vogt, T. N. (2010). Hydrolytic cleavage of DNA-model substrates promoted by polyoxovanadates. Dalton Trans. 14, 585–592. doi: 10.1039/b913471a
Swallow, A. G., and Barnes, W. H. (1964). The structure of the decavanadate Ion. J. Am. Chem. Soc. 86, 4209–4210. doi: 10.1021/ja01073a070
Thompson, H. J., Chasteen, D. N., and Neeker, L. (1984). Dietary vanadyl (IV) sulfate inhibits chemically-induced mammary carcinogenesis.Carcinogenesis 5,849–851. doi: 10.1093/carcin/5.6.849
Treviño, S., Díaz, A., Sánchez-Lara, E., Sarmiento-Ortega, V. E., Flores-Hernandez, J. A., Brambila, E., et al. (2018). Pharmacological and toxicological threshold of bisammonium tetrakis 4-(N,N-Dimethylamino) pyridinium decavanadate in a rat model of metabolic syndrome and insulin resistance. Bioinorganic Chem. Appl. 2018:2151079. doi: 10.1155/2018/2151079
Treviño, S., Velázquez-Vázquez, D., Sánchez-Lara, E., Diaz-Fonseca, A., Flores-Hernandez, J. Á., Pérez-Benítez, A., et al. (2016). Metforminium decavanadate as a potential metallopharmaceutical drug for the treatment of Diabetes mellitus. Oxid. Med. Cell. Longev. 2016:6058705. doi: 10.1155/2016/6058705
Treviño, S., Sánchez-Lara, E., Sarmiento-Ortega, V. E., Sánchez Lombardo, I., Flores-Hernández, J. A., Pérez-Benítez, A., et al. (2015). Hypoglycemic, lipid-lowering and metabolic regulation activities of metforminium decavanadate (H2Metf)3 [V10O28]∙8H2O using hypercaloric-induced carbohydrate and lipid deregulation in wistar rats as biological model. J. Inorg. Biochem. 147, 85–92. doi: 10.1016/j.jinorgbio.2015.04.002
UK Prospective Diabetes Study (UKPDS) Group. (1998). Effect of intensive blood-glucose control with metformin on complications in overweight patients with type 2 diabetes (UKPDS 34). Lancet 352, 854–865.
Vermes, A., Guchelaar, H. J., and Dankert, J. (2000). Flucytosine: a review of its pharmacology, clinical indications, pharmacokinetics, toxicity and drug interactions. J. Antimicrob. Chemother. 46, 171–179. doi: 10.1093/jac/46.2.171
Vilter, H. (1984). Peroxidase from Phaeophyceae A Vanadium(V)-dependent peroxidase from Ascophyllum nodosum. Phytochemistry 23, 1387–1390. doi: 10.1016/S0031-9422(00)80471-9
Weng, J., Hong, M., Liang, Y., Shi, Q., and Cao, R. (2002). A nucleobase–inorganic hybrid polymer consisting of copper bis (phosphopentamolybdate) and cytosine. J. Chem. Soc. Dalton Trans. 3, 289–290. doi: 10.1039/b111457n
Winkler, P. A., Huang, Y., Sun, W., Du, J., and Lü, W. (2017). Electron cryo microscopy structure of a human TRPM4 channel. Nature 552, 200–204. doi: 10.1038/nature24674
Yraola, F., García-Vicente, S., Marti, L., Albericio, F., Zorzano, A., and Royo, M. (2007). Understanding the mechanism of action of the novel SSAO substrate (C7NH10) 6(V10O28)· 2H2O, a prodrug of peroxovanadate insulin mimetics. Chem. Biol. Drug.Des. 69, 423–428. doi: 10.1111/j.1747-0285.2007.00516.x
Zhu, M., Lu, L., and Yang, P. (2003). N, N-dimethylbiguanidium nitrate. Acta Cryst E 59, o586–o588. doi: 10.1107/S1600536803006639
Zhu, M., Lu, L., Yang, P., and Jin, X. (2002). Bis (1, 1-dimethylbiguanido) nickel (II). Acta Cryst. E 58, m272–m274. doi: 10.1107/S1600536802007869
Keywords: decavanadate, metformin, cytosine, X-ray crystal structure, vibrational spectroscopy, 51V-NMR, theoretical studies, polyoxovanadates
Citation: Sánchez-Lara E, Treviño S, Sánchez-Gaytán BL, Sánchez-Mora E, Castro ME, Meléndez-Bustamante FJ, Méndez-Rojas MA and González-Vergara E (2018) Decavanadate Salts of Cytosine and Metformin: A Combined Experimental-Theoretical Study of Potential Metallodrugs Against Diabetes and Cancer. Front. Chem. 6:402. doi: 10.3389/fchem.2018.00402
Received: 29 June 2018; Accepted: 20 August 2018;
Published: 02 October 2018.
Edited by:
Debbie C. Crans, Colorado State University, United StatesReviewed by:
Manuel Aureliano, University of Algarve, PortugalSamar Kumar Das, University of Hyderabad, India
Copyright © 2018 Sánchez-Lara, Treviño, Sánchez-Gaytán, Sánchez-Mora, Castro, Meléndez-Bustamante, Méndez-Rojas and González-Vergara. 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: Enrique González-Vergara, ZW5yaXF1ZS5nb256YWxlekBjb3JyZW8uYnVhcC5teA==