Ternary Copper Complex of L-Glutamine and Phenanthroline as Counterions of Cyclo-Tetravanadate Anion: Experimental–Theoretical Characterization and Potential Antineoplastic Activity

: Over the last decade, therapeutic metallodrugs have become substantially effective in the treatment of cancer. Thus, developing new effective anticancer drugs is a significant research area against the continuing increase in cancers worldwide. In the search for heterobimetallic prodrugs containing V/Cu, a new cyclo-tetravanadate was synthesized and characterized by UV-visible and FTIR spectroscopies and single-crystal X-ray diffraction. L-Glutamine and 1,10-phenanthroline al-low the crystallization of [Cu(L-Gln)(phen)(H 2 O)] 4 [V 4 O 12 ]∙8(H 2 O) ( 1 ), in which the cyclo-tetra-vanadate acts as a free anion. Density functional theory (DFT) calculations were carried out to characterize the frontier molecular orbitals and molecular electrostatic potential. Global reactivity indexes were calculated and analyzed to give insight into the cyclo-tetravanadate anion and complex counterions interactions. Also, using Bader’s theory of atoms in molecules (AIM), non-covalent interactions were analyzed. Docking analysis with the Casiopeina-like complex resulting from the hydrolysis of compound 1 provided insights into these complex potential anticancer activities by interacting with DNA/tRNA via H-bonds and hydrophobic interactions. The release of both components could act together or separately, acting as prodrugs with potential dual antineoplastic activities.


Introduction
The synthesis of metallodrugs for cancer treatment has been rising due to its substantial success in the last decade [1,2]. Cancer is a multi-factorial, multibillion-dollar public health epidemic around the world. Considering that in 2020 there were approximately 18.1 million new cases and 9.6 million deaths worldwide, it is important to find safe and low-cost ways to fight it [3]. In the Region of the Americas, it is the second leading cause of death. If no steps are taken to prevent and control cancer, the number of persons diagnosed with cancer will rise by 55 percent by 2040, totaling roughly 6.23 million people. [4]. About 195,499 new cases and 90,222 deaths were estimated in Mexico in 2020, with a higher prevalence of breast cancer, followed by prostate cancer with 15.3% and 13.7%, respectively [3]. While several drugs containing metals based on gold, ruthenium, gallium, titanium, iron, and copper are in preclinical and clinical phase I and II trials [5,6], the most effective antitumor agents used in clinical practice remain cis-platinum and also second-and third-generation platinum coordination compounds (carboplatin, oxaliplatin, and picoplatin). Nevertheless, platinum-based products' clinic use carries significant side effects, including nephrotoxicity, neurotoxicity, ototoxicity, hepatotoxicity, and myelosuppression [7,8]. Endogenous metal-based antitumor drugs (Co, Cu, Zn, and Fe) are less toxic than platinum analogs and display promising pharmacological properties [9,10].
In the specific case of copper, it provides an ideal geometry to interact with DNA molecules, making them a viable alternative to platinum-based anticancer drugs, with the advantage that copper is better tolerated and can be more easily managed in vivo than other transition metals [11,12]. Complexes containing symmetric aromatic binders, such as 1,10-phenanthroline and 2,2'-bipyridine, and how they interact with DNA have received much attention [13][14][15][16]. In 1979, Sigman et al. discovered the first metallonuclease by demonstrating that [Cu(phen)2] 2+ complex inhibits DNA and RNA polymerase activity and can induce DNA strand scission [17]. Numerous studies since then have demonstrated that the geometry of the metal center, combined with planar bidentate ligands, contains an optimal structure for interacting with a wide variety of biological molecules and possesses antitumoral and antiviral properties [18][19][20][21][22]. The synthesis, design, and development of complexes as anticancer agents have been presented in several reviews over the last decade, indicating an exponential increase in the field [23][24][25][26][27][28][29][30][31][32][33][34][35]. Ng et al. (2014) showed that copper(II) complexes containing 1,10-phenanthroline and amino acids were significantly more antiproliferative and induced more apoptotic cell death in the cisplatinresistant MDA-MB-231 cells than in the non-cancerous MCF10A cells, and only induced cell cycle arrest in cancer cells. Additionally, these compounds were found to be effective against cervical (HeLa), ovarian (SKOV3), lung (A549, PC9), non-small cell lung (NPC) (Hone1, HK1, C666-1), breast (MCF7, T47D), lymphoma and leukemia (Nalmawa, HL60), and colorectal (SW480, SW48, HCT118) cancer cell lines with IC50 values (24 h) ranging from 1.7 to 19.0 μM. Importantly, these compounds were even more effective than cisplatin in some cases [36]. A series of ternary copper(II)-L-dipeptide-neocuproine complexes demonstrated cytotoxicity against cancer cells, including the triple-negative breast cancer cell line MDA-MB-231 [37]. The complex bis[(2-chloro)chloro(1,10-phenanthroline)copper(II)] exhibited a high degree of reactivity. Anticancer activity was observed against the B16, MDA-MB-32, A549, HT-29, and SF cell lines, with an average of IC50 of 0.726 g/mL (1.15 μM), compared to cisplatin's 4.88 g/mL (16.3 μM). Furthermore, it demonstrated a higher selectivity for cancer cells than human bone marrow stem cells, and in rats, it is less toxic than cisplatin [38]. Copper-based metallodrugs loaded in nanoparticles have also shown promising results. The survival rate of a murine B16 melanoma model, when given chitosan nanoparticles loaded with Cassiopeina, was higher than those that received the drug as such. This was achieved due to a longer residence time at the site of action because the release of the substance from the nanoparticles was achieved at the acidic pH of the tumor [39]. Reina [40] by substituting acetylacetonate for 2-aminomethyl benzimidazole (2AMB) as a bidentate chelate ligand. The change increases hydrophilicity and antiproliferative activity, resulting in two active complexes against the human tumoral cell lines, HeLa and MCF-7. A small group of copper complex exerts an effective inhibitory action on topoisomerases, which participate in the regulation of DNA topology. Copper complexes of topoisomerase inhibitors work via different molecular mechanisms that affect the cell cycle checkpoints and death effectors. Expanding this family of highly active anticancer drugs and their use with other emerging cancer therapies opens new avenues for treating cancers [28]. Complexes of Cu(II)-phenanthroline with L-asparagine and L-methionine were highly toxic to the prostate cancer cell lines (DU-145 and PC3), breast cancer cell lines (MDA-MB-231 and MCF7), and MV3 (melanoma) [41]. The mechanism is a redox reaction of Cu (II)/Cu (I) that catalyzes the formation of reactive oxygen species (ROS). The structure of this type of complex consists of a copper(II) core with a distorted square pyramidal geometry, which, in the presence of ascorbate, exhibits effective DNA cleavage activity at micromolar concentrations, producing hydroxyl radicals as the active species [42]. Recently, it has been demonstrated that sulfonamide copper(II) complexes with 1,10-phenanthroline and 2,2′-bipyridine possess potent nuclease activities with antiproliferative and anti-M. tuberculosis actions [43]. Casiopeinas ® have gained an important place in the efforts to generate antineoplastic drugs with copper, as expressed in the recent invention of a safe parental solution with Cassiopeia III-ia [44]. In the review of copper coordination compounds as biologically active agents," it is stated that the redox activity of copper ions, as well as the stability of copper coordination compounds in the bloodstream, and the extremely encouraging therapeutic outcomes both in vitro and in vivo, show the potential of copper coordination compounds, which is expected to gain widespread clinical use" [45].
Vanadium, on the other hand, has therapeutic properties that are well known. Its pharmacological properties could be used to treat diabetes, however, concern about toxicity is still an issue [46,47]. However, some studies have also shown that vanadium could be used to treat cancer since it causes cell apoptosis, displaying cytotoxic and antiproliferative effects [48]. In fact, several of the metabolic pathways are used by both diabetes mellitus and cancer [49][50][51][52]. Recent studies of pancreatic cancer, malignant melanoma, glioma multiforme, and an aggressive brain cancer have shown that vanadates and oligovanadates function as anticancer drugs [53][54][55]. Also, in cytotoxicity experiments with human osteosarcoma (MG-63) and human colorectal adenocarcinoma (HT-29) cell lines, the compound bis (4,7-dimethyl-1,10-phenanthroline) sulfatooxidovanadium (IV), known as Metvan, showed altered cell viability of both cancer cell lines in the low concentration range (0.25-5.0 μM) [56]. Another recent study showed the importance of the interaction of polyoxovanadates with the plasma membrane to subsequently hydrolyze and internalize the compound, thus allowing partial participation in the generation of ROS and electron transfer processes [57,58]. In addition, vanadium-based nanomaterials have been found to be helpful for theranostic applications. Vanadium nitride nanosheets are readily convertible into biocompatible materials that display near-infrared (NIR) absorption, producing ROS and effectively killing cancerous tissue [59]; vanadium oxide has an activity similar to peroxide, which would be helpful in tumor-specific chemodynamic therapy [60]. In addition, it is well known that oxidovanadium complexes with organic ligands have cytotoxic or differentiating properties against various cancer cell types [61][62][63][64][65][66]. However, their limited use in clinical trials stems primarily from concerns about the long-term toxicity of such complexes due to a lack of data. However, recent administration of [VO(HSHED)dtb] complex in mice did not show any signs of toxicity up to a dose of 300 mg/kg. The complex was less toxic than orthovanadate salt, consistent with the compound being partially intact during the administration [67]. The low toxicity is attributed to the redox properties obtained when combining the redox-active ligand 3,5-di(tert-butyl)catechol with the hydrolytic stability of the [VO(HSHED)dtb], which prevent the formation of the vanadate and catechol ligand. These results point out that the protection of even anionic vanadium species may help increase activity and reduce toxicity. Although anionic compounds of vanadium can survive the transit in the bloodstream and enter the cell by several mechanisms, it is possible to suggest that they generate mainly V1, V2, and V4 vanadates through speciation [63,64]. Thus, it is crucial to investigate the cyclo-tetravanadate anion since it will be the primary species at physiological pH [68]. Anionic compounds of vanadium(V) have shown interesting anticancer properties [69,70]. Recently, the compounds (NH4)4Li2V10O28•10H2O and Mg(H2O)6(C4N2H7)4V10O28•4H2O demonstrated dose-dependent antiproliferative activity on human cancer cells U87, IGR39, and MDA-MB-231 [71,72]. Again, nanotechnology provides new possibilities. Cationic liposomes have been used to encapsulate a hydrophobic oxidovanadium complex that has shown cytotoxic properties against neuroblastoma cell lines while also improving bioavailability [73]. V(IV)-curcumin-bipyridine (VCur) was encapsulated in vitro using magnetic cationic liposomes, which gives the compound excellent stability and solubility in physiological media [74]. Biocompatible vanadium nitride nanosheets have shown relevant near-infrared (NIR) absorption that can be used for imaging and produce reactive oxygen species upon NIR excitation, effectively killing cancer tissue. Recent studies of eleven vanadium species, compounds, and materials, in several human melanoma cells lines, have shown interesting antitumor capabilities to a variety of effects, including: (1) cell viability, (2) cell morphology changes and apoptosis, (3) cell-cycle arrest, (4) ROS production, (5) mitochondrial dysfunction, (6) protein expression, and (7) in vivo tumor regression and survival rates. Therefore, it was concluded that essential questions would be answered this decade to push forward the practice of metals compounds-particularly vanadium-to treat cancer [75].
As for copper, the new results and new strategies for delivery point to significant advances in this decade. It is essential to consider that compounds of vanadium and copper will act as prodrugs since they will undergo interactions with proteins such as human serum albumin and transferrin, and these will release their ligands, or in some robust cases, transport of the intact prodrug in the bloodstream could take place [76][77][78][79][80][81]. Speciation outside or inside the cells is most probable, but regardless of this, important anticancer actions are provided by polyoxidovanadates (POVs) and Casiopeinas ® , so we expect dual anticancer activities of the compounds reported here. It is now clear that both metals will provide new metallodrugs to help in the global problem of cancer treatment in the near future. Also, POVs have a diverse spectrum of structural, biological, and pharmacological properties that are intriguing and potentially useful. Recent research has indicated that these compounds may effectively treat solid tumors, DNA and RNA viruses, and drug-resistant bacteria. Changes in cell structure, interference with the ion transport system, inhibition of mRNA synthesis, and disruption of metabolic pathways and communication mechanisms are all known antibacterial modes of action for various POVs. Additionally, POVs have been found to reduce mitochondrial respiration in vitro. However, in the future, researchers need to explore vanadium-containing POV-based nanohybrids rather than pure POVs to reduce toxicity and boost efficacy [82].
As the first collection of V/Cu heterobimetallic compounds with the potential to be used as metallodrugs in cancer treatment, our group has identified several cyclo-tetravanadates [83][84][85]. We recently reported the synthesis and theoretical-experimental characterization of two compounds resulting from our quest to complete the family of copper cyclo-tetravanadate complexes. Interestingly, the compounds obtained were different, although the same synthetic technique was used; however, it is possible to consider them Casiopeina ® analogs [86]. Here, a new cyclo-tetravanadate is presented. This compound was characterized by visible and IR spectroscopies as well as single-crystal X-ray diffraction. In addition, the compound was studied using DFT computational methods. The frontier molecular orbitals and global reactivity indexes were analyzed to show interesting characteristics about the donor-acceptor interactions. These insights about the compound's reactivity were corroborated by analyzing the non-covalent interactions using the AIM approach. The Casiopeina-like complex in docking studies shows the potential to interact with DNA/RNA, thus providing a new compound with potential dual anticancer activity in various cell lines.

Results
Compound 1 crystallizes in the triclinic space group P1, with only one asymmetric unit filling the unit cell. Figure 1 shows a representation of the molecule and Table 1 presents the corresponding crystal data. The compound contains four [Cu(L-Gln)(phen)(H2O)] + cations displaying a squared-pyramidal geometry (τ5 = 0.017), of which the two N donors of phenanthroline and the amino and carboxyl groups of glutamine are coordinated in basal positions. As shown in Figure 1, four water molecules occupy the apical positions of all the square pyramids. In this case, the glutamine ligands are L isomers. However, two hydrogen atoms in the alpha carbon point in the opposite direction of the coordinated water, so they are considered anti-conformers; the other two pointing in the same direction are considered syn-conformers, as seen in Figures 2 and 3. It can be observed that all copper coordination compounds interact non-covalently with the cyclotetravanadate structure through different hydrogen bonds. Figure 2 shows the hydrogen bonds surrounding the cyclo-tetravanadate anion. Eight H2O molecules stabilize the structure, five of them (O33 to O37) in a chain array connecting one complex molecule with the [V4O12] 4-anion. Tables S1-S5 in the supplementary section contain more details of structural data.  L-glutamine is present as two syn and two anti-conformers, data from [87].

Infrared Spectroscopy
Due to the relatively rigid structure of the polyoxovanadate ions, the IR spectra of the individual types of polyoxovanadates have characteristic bands in the region of asymmetric and symmetric VO stretching vibrations. As a result, IR spectroscopy can be used to characterize such compounds [88][89][90]. The spectrum of Compound 1 in the region 1000-500 cm -1 is presented in Figure 4. The complete spectrum is presented in Figure S1 in the Supplementary Materials. Bands resulting from the symmetric and asymmetric vibrations of terminal V(Ot)2 commonly occur between 910-1000 cm −1 for νas and 860-900 cm −1 for νs. On the other hand, bands below 850 cm −1 could be assigned to the asymmetric vibrations of the bridging V-Ob-V groups, and the small band at 530 cm −1 could be assigned to νs V-Ob-V [91][92][93][94][95][96][97]. The positions of absorption bands in the region of VO bond stretching vibrations (500-1000 cm −1 ) will make it possible to determine whether a compound's structure is made up of cyclic [V4O12] 4-as a free anion or is involved in several hydrogen bond interactions, as shown in Figure 4. Compound 1 shows that cyclo-tetravanadate ions are involved in different hydrogen bond interactions, causing the splitting of the typical bands. Several shoulders appear in the main bands, indicating different V(Ot)2 and V-Ob-V vibration strengths. These are seen in the inset, where the second derivative spectrum is depicted. The presence of several shoulders can be evidenced by second derivative spectroscopy [98][99][100]. The shoulders are now presented as peaks and show how different environments are experienced by each VO bond, consistent with the hydrogen bond interactions presented in the solid-state structures. As for the bands of phenanthroline, the only bands seen in this region are located at 727 cm −1 and 736 cm −1, respectively, and are assigned according to Schilt and Taylor [101] and Baskaran et al. [102], who observed two strong bands at 850 cm −1 and 725 cm −1 in the free ligand. The band at 850 cm −1 was assigned to the out-of-plane motion of the carbocyclic ring's hydrogen atoms and the 725 cm −1 to those of the two heterocyclic rings' hydrogen atoms. The band at 850 cm −1 is overlapped with the V-Ob-V antisymmetric bands. As for the amino acid moiety, the strong bands of the amino and carboxylic groups are outside this region [103].  Figure S2 shows the visible spectrum of Compound 1 after it had been dissolved in phosphate-buffered saline (PBS) solution at pH 7.4. Compound 1 showed an asymmetric absorption band at λmax value 627 nm, and a molar extinction coefficient ε = 45 L mol −1 cm −1 . Valora et al., has shown that several ternary complexes containing 1,10-phenanthroline and amioacidates have d-d bands with λmax value at 618 ± 6 nm and ε of 57 ± 5 due to the promotion of low energy d electrons to the dx 2 −y 2 hole, which is the value expected for a CuN3O in-plane chromophore [104,105]. Thus, although in the solid-state, the geometry of the compound was a square pyramid, it is likely that in aqueous solution, the compound exhibits a unique molecular geometry, which is a tetragonally elongated octahedron.

Computational Calculations
The molecular structure, molecular electrostatic potential (MEP) map, and isosurfaces of the highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) are shown in Figure 5. Global reactivity descriptors, such as chemical potential (μ), electronegativity (χ), hardness (η), softness (s), and electrophilicity index (ω), were evaluated based on a conceptual DFT approach using the following equations: μ = (EHOMO+ELUMO) ⁄ 2; χ = −(EHOMO+ELUMO) ⁄ 2; η = ELUMO−EHOMO; s = 1 ⁄ η, ω = μ 2 ⁄ 2η, where EHOMO and ELUMO are the frontier molecular orbitals energies. Table 2 shows the results. The electron charge transfer process can be analyzed in terms of chemical potential (μ) (or electronegativity (χ = −μ)), hardness (η) (or softness (s = η −1 )), and electrophilicity (ω). The results show that the chemical potential value, μ, of Compound 1, indicates the capacity for exchanging electrons. The chemical hardness value, η, indicates the resistance for changing its electron distribution. On the other hand, the behavior of the electrophile or nucleophile can be related to a higher or lower ω index. In conjunction with the MEP analysis, global reactivity indexes give insights regarding the electrophilic or nucleophilic power of the compounds' sites, which is important for analyzing the interaction. Topological parameters were used to characterize intramolecular interactions, such as electron density ( ), Laplacian ∇ 2 ( ), Lagrangian kinetic energy G, potential energy density V, Hamiltonian kinetic energy H, interaction energy EH…Y, and interatomic distance . Results are shown in Table S6 | ( )|. The molecular graph for Compound 1 is shown in Figure 6. Cyan dots represent the bond critical points (BCPs), yellow dots represent the ring critical points (RCPs), and orange dots represent the cage critical points (CCPs). In Figure 6, the main H-bond interactions between O atoms of the cyclo-tetravanadate anion with water molecules and Hatoms of glutamine are shown. Table S6

Molecular Docking
It was hypothesized that in an aqueous solution, because of reduction, re-oxidation, and hydrolysis, the vanadium and copper compounds would be released and act separately; thus, the copper Casiopeina like complex was studied as the corresponding aquo complex. In Table 3, the docked binding energies of the interaction with DNA/RNA test molecules, corresponding to the top molecular pose (lowest energy) for complex ions already reported, are presented for comparison.  [83,86,[106][107][108].
In Table 4, the docked binding energies and the interaction with DNA/RNA corresponding to the top molecular pose (lowest energy) for complex ions [Cu(L-Gln)(phen)(H2O)] + are shown. As was expected from previous studies, all docked structures occupy similar positions in the minor groove of the 1BNA/DNA and intercalate in the 151D/DNA fragment structures (Figures 7 and 8). The compounds interact well with the DNA test molecules, yet considering the binding energies among these complexes, the best affinity energy is found in the complex ion anti-[Cu(L-Gln)(phen)(H2O)] + . Also, it can be implied that ligand phenanthroline has an essential role in the interaction; this result agrees with our previous results [83,84,86]. Regarding the interaction with an RNA test molecule (Figure 9), the hydrogen bonds between L-syn and tRNA involve C24, C40, and PSU39 (pseudouridine). Salt bridge interactions involve A38 and G26 and the van der Waals interactions include A29, G42, and U41. On the other hand, L-anti interactions consist of hydrogen bonds with C27 and G26 and van der Waals interactions with C25, C28, A29, U41, G30, PSU39, A38, and 5MC40 (5-methyl-cytosine).   However, the combination of interactions in this study favors the syn-complex. The interactions involve several H-bonds, hydrophobic interactions, salt bridges, attractive charges, and π-anions, as shown in Figure 10.

Discussion
In our search to provide vanadium complexes that mimic the interactions of decavanadates and cyclo-tetravanadates with side-chain groups of proteins with positive or hydrogen-bond donor amino acids, the series of amino acids studied has been expanded to L-glutamine. A new heterobimetallic cyclo-tetravanadate was synthesized and characterized. L-glutamine and phenanthroline mixed copper complex was used to provide potential hydrogen bonds to anchor the complexes of planar di-imines, which are well known to interact with DNA and RNA test molecules [109][110][111]. Compound 1 crystallizes in the non-centrosymmetric space group P1. Interestingly L-glutamine is present in two conformations, where two of the copper complexes show a syn-conformation, while the other two show an anti-conformation.
The existence of [V4O12] 4-was confirmed by FTIR spectroscopy analysis, where the characteristics bands are located between 910-1000 cm −1 for asymmetric and symmetric vibrations of terminal V(Ot)2. The asymmetric vibrations of the bridge V-Ob-V groups could be allocated to bands below 850 cm −1 . It is feasible to establish that cyclo-tetravanadate ions are involved in several hydrogen-bond interactions, causing the usual bands to divide; as a result, there are numerous shoulders, each indicating different V (Ot)2 and V-Ob-V vibration strengths. Although second derivative spectra are rare, these are very informative since shoulders end up becoming peaks. On the other hand, the visible spectrum shows an asymmetric absorption band at λmax value 627 nm because there is the promotion of d-d electrons from (dxz, dyz), dxy, and dz 2 to the orbital dx 2 -y 2 . Therefore, although in solid-state, the geometry corresponds to a square-based pyramid, it may change in aqueous solution to elongated tetragonal octahedron.
Using DFT calculations, the structure was minimized from the crystallographic data, allowing the frontier molecular orbitals and the map of electrostatic potential to be calculated. Global reactivity descriptors, such as chemical potential (μ), electronegativity (χ), hardness (η), softness (s), and electrophilicity index (ω), were evaluated based on conceptual DFT approach. The presence of hydrogen bonds dominates all the interactions. Thus an AIM study was conducted to characterize the compound.
It is important to mention that there have only been a few reports on cyclo-tetravanadates: two of them were conducted with organic molecules (C6H9N2)4(V4O12)•4H2O and [(CH3)3-CNH3]4(V4O12) [112,113], vanadium [114,115], iron [116], cobalt [117], nickel [118], zinc [119,120], and Zr [121]. In the particular case of copper [96,122,123], four complexes have been reported; in all those cases, copper ions were in a trigonal bipyramidal geometry. There have also been also four cyclo-tetravanadates containing three metals [124]. , where the cyclo-tetravanadate acts as a bridge and the copper atoms are in a square pyramidal geometry [83][84][85]. In this paper, we present [Cu(L-Gln)(phen)(H2O)]4(V4O12), the first anionic cyclo-tetravanadate of this type in which the anion is not coordinated to the contraionic copper complexes. Concerning the planarity of the cyclo-tetravanadates, there are many possibilities, ranging from completely flat (D4h symmetry) to Cs symmetry; the compound reported here has Cs symmetry. Although other cyclo-tetravanadates, such as [Fe(bipy)3]2[V4O12]•10H2O, Zn(2,2′-bipy)3] 2V4O12•8H2O, show uncoordinated (V4O12) anions, we expected a similar structure as our compound with glycinate. However, the introduction of the chiral carbon of the glutamine aminoacidate requires the use of DL-glutamine to surround the cyclo-tetravanadate. Progress has been made in that direction [125]. As pointed out before, the release of the components upon hydrolysis generates active species of vanadium and copper [83][84][85]. The cationic copper complexes generated have a strikingly similar molecular structure to Casiopeinas ® (CAS), a class of copper-based drugs developed by Ruiz-Azuara and colleagues [126][127][128]. CAS is the first family of antineoplastic copper compounds that have reached the clinical phase in the world [129]. has shown that Cas II-gly binds to low-molecular-mass bioligands found in blood serum (citric, L-lactic acid, and L-histidine) and cytosol (reduced glutathione (GSH), reduced nicotinamide adenine dinucleotide (NADH), adenosine triphosphate (ATP), and L-ascorbic acid). The bioligands replaced glycinato in mixed-species Cu II -Me2phen-bL/cL. The creation of these adducts may aid copper transport to the target organs [132]. Our computational analyses show that several compounds have better binding energies with nucleic acid test molecules than those already in phase 1 clinical trials. It can be observed that the ligand phenanthroline has an essential role in these interactions since it is the moiety that intercalates well into the DNA [133]. On the other hand, when the docking study was carried out considering 6TNA, the binding energies obtained were below -9.3 kcal/mol. The combination of interactions in this study favors the syn-complex. The interactions involved several H-bond, hydrophobic interactions, salt bridges, attractive charges, and π-anions. Therefore interactions with tRNA near the anticodon arm are needed to study experimentally in great detail.

General Considerations
The synthesis of Compound 1 was conducted using ammonium metavanadate (NH4VO3), L-glutamine, and 1,10-phenanthroline hydrochloride (phen) purchased from Merck Mexico. KOH was purchased from Fermont, and CuCl2•2H2O was purchased from Química Dinámica S. A. de C. V, Nuevo Leon, Mexico. All manipulations were conducted without specific purification of solvents and reagents. The complexes' electronic spectra were determined by UV-Vis spectroscopy with a xenon lamp Varian Cary 50 spectrophotometer(Varian, Palo Alto, CA, USA) and using a quartz cuvette of 1 cm path length. The infrared spectrum was obtained from 500 to 4000 cm −1 in KBr pellets using an IR Digilab, Mod. Scimitar Spectrophotometer with FTIR (Digilab, Hopkinton, MA, USA).

Crystallization and Synthesis
Compound 1 was prepared by adding 1.00 × 10 −3 mol (0.198 g) of 1,10-phenanthroline hydrochloride to 30 mL distilled water with stirring and moderate heat; once dissolved, 1.00 × 10 −3 mol (0.170 g) of CuCl2•2H2O was added, and the clear solution was adjusted to pH 9.5 by slow addition of (10 percent) KOH solution. Subsequently, when a clear blue solution was achieved, 2.00 × 10 −3 mol of L-glutamine were added while stirring. Then, NH4VO3, 1.00 × 10 −3 mol (0.116 g) in 15 mL H2O, was added dropwise to the solution. The final pH at room temperature was 9. The solution was filtered and left at room temperature. Blue prismatic crystals were recovered by filtration after six days.

Single Crystal X-ray Diffraction
Single crystal of [Cu(L-Gln)(phen)(H2O)]4[V4O12]•(H2O)8 (Compound 1), suitable for X-ray diffraction, was selected and mounted on an Xcalibur, Atlas, Gemini diffractometer. The crystal was kept at 293.15 K during data collection. Using Olex2 [134], the structure was solved with the SHELXD [135] structure solution program using dual space and refined with the SHELXL [136] refinement package using least squares minimization. Compound 1 crystallizes in the triclinic space group P1, with only one asymmetric unit filling the unit cell. Figure 1 shows the complete molecule crystal data. Selected crystal data and details are shown in Table 1 and Tables S1-S5

Computational Methods
Molecular structure, frontier molecular orbitals (FMO), and molecular electrostatic potential (MEP) of Compound 1 was carried out using the hybrid functional B3LYP [137], including the empirical dispersion method of Grimme with Becke-Johnson damping, GD3BJ [138]. The double zeta split-valence Def2SVP basis set [139] for C, H, O, and N atoms and LanL2DZ basis set [140] for Cu and V atoms were used. In addition, LanL2DZ was used as an effective core potential (ECP) for Cu and V atoms [141]. The effect of the implicit solvation was included with the polarizable continuum model (PCM), using the integral equation formalism variant (IEFPCM) with water as the solvent [142]. Calculations were performed with the Gaussian16 and GaussView programs [143,144]. Global reactivity descriptors were calculated from the FMO energies based on the conceptual DFT approach [145]. Atoms in molecules analysis (AIM) using AIMAll software (Version 14.11.23, TK Gristmill Software, Cambridge, MA, USA) was performed for analyzing the non-covalent interactions [146].

Molecular Docking
Molecular docking analysis was performed with the semi-flexible method. The DNA fragments used were Dickerson-Drew dodecamer (DDD) with the sequence d(CGCGAATTCGCG) (PDB ID: 1BNA [147]) and DNA fragment with the sequence d(CGATCG) (PDB ID:151D [148]. The RNA docking was carried out using the yeast tRNA (PDB: 6TNA) [149] and considered a rigid entity. Complete flexibility was allowed for the coordination compounds [150], and it was performed using the Autodock Tools 1.5.6 software (accessed on 1 September 2021) [151], which includes the addition of polar hydrogens and empirical particles of atomic charges (Gasteiger-Marsili method). A grid box that encloses the entire DNA fragment was used with sizes 70, 70, and 120 Å for the 1BNA DNA fragment and 60, 60, and 76 Å for the 151D DNA fragment. For tRNA, blind docking was carried out using three different boxes that enclosed the entire tRNA with the sizes 96, 80, 108 Å, followed by a re-docking of the compounds in their docked pose with potential and preferred minimized energy using a box centered on the compounds with the sizes 40, 40, and 40 Å. The grid spacing for all the docking calculations was set to the 0.375 Å default value, using the Lamarckian genetic algorithm (LGA) searching methods. The parameters for the copper atom were the sum of VDW radii of two similar atoms (3.50 Å), plus the VDW well depth (0.005 kcal/mol), plus the atomic solvation volume (12.0 Å3), plus the atomic solvation parameter (−0.00110). The H-bond radius of the heteroatom in contact with hydrogen (0.0 Å), the well depth of the H-bond (0.0 kcal/mol), and different integers showed the type of H-bonding atom and indexes for the generation of the auto grid map (0, −1, −1, 1, respectively). The corresponding figures were prepared using Biovia Discovery Studio [152].

Conclusions
A new V/Cu heterobimetallic compound was synthesized and characterized by visible and FTIR spectroscopies. Compound 1 crystallizes in the triclinic P1 space group with the cyclo-tetravanadate anion acting as a free ion, only bound to the copper moieties through hydrogen bonds and water molecules. The cyclo-tetravanadate anion provides a hydrogen-bonded bridge between two conformers syn and two anti of the [Cu(L-Gln)(phen)(H2O)] cationic complexes. The cyclo-tetravanadate is not coordinated to the copper, and thus, it can be easily separated in an aqueous solution to act independently. The AIM analysis indicated that the hydrogen-bond interactions govern the structural arrangement in the compound. Additionally, the Cu-O interactions (between Cu with O of apical water or with O of cyclo-tetravanadate) have high interaction energy values, indicating that they are responsible for the supramolecular arrangement. Also, the hydrogen bonds surrounding the cyclo-tetravanadate could serve as models for the interaction of side-chain hydrogen bond donors in hydrophilic pockets in proteins.
The copper moieties are reminiscent of Casiopeinas IIa y IIIia, which are in phase I clinical trials in Mexico. Docking studies indicated that both compounds can interact with DNA and RNA test molecules, acting as groove binders and intercalating agents. It is important to point out that a combination of electrostatic and hydrophobic interactions, ππ, and hydrogen bonds, are responsible for binding DNA/RNA molecules. The energies of interaction are comparable with other well-known and proven agents with anticancer activities. The binding energies of the RNA interactions are more favorable than other compounds already reported and point to RNA potential metallodrugs, which has become relevant as an easy-to-reach target since other important RNA molecules have been discovered to be involved in specific cancers [153]. The recently explored direct injection into the tumors or the encapsulation in liposomes or nisomes opens an excellent opportunity to improve the potential antitumor activities of the compound here presented. The recent reports in cancer therapy and the new strategies for delivering vanadium and copper point to significant advances in the near future.

Supplementary Materials:
The following are available online at www.mdpi.com/article/10.3390/met11101541/s1. Figure S1: FTIR spectrum of Compound 1 in the 500-3500 cm -1 region. Figure S2: The visible spectrum of Compound 1. Crystallographic data are presented in Table S1-S5. Topological parameters (a.u.), interaction energies EH…Y (kcal mol −1 ), and interatomic distances Dint (Å) are presented in Table S6. The check CIF of Compound 1 is also included in the Supplementary Materials.