2-Aminopyrimidinium Decavanadate: Experimental and Theoretical Characterization, Molecular Docking, and Potential Antineoplastic Activity

: The interest in decavanadate anions has increased in recent decades, since these clusters show interesting applications as varied as sensors, batteries, catalysts, or new drugs in medicine. Due to the capacity of the interaction of decavanadate with a variety of biological molecules because of its high negative charge and oxygen-rich surface, this cluster is being widely studied both in vitro and in vivo as a treatment for several global health problems such as diabetes mellitus, cancer, and Alzheimer’s disease. Here, we report a new decavanadate compound with organic molecules synthesized in an aqueous solution and structurally characterized by elemental analysis, infrared spectroscopy, thermogravimetric analysis, and single-crystal X-ray diffraction. The decavanadate anion was combined with 2-aminopyrimidine to form the compound [2-ampymH]6[V10O28] · 5H2O ( 1 ). In the crystal lattice, organic molecules are stacked by π – π interactions, with a centroid-to-centroid distance similar to that shown in DNA or RNA molecules. Furthermore, computational DFT calculations of Compound 1 corroborate the hydrogen bond interaction between pyrimidine molecules and decavanadate anions, as well as the π – π stacking interactions between the central pyrimidine molecules. Finally, docking studies with test RNA molecules indicate that they could serve as other potential targets for the anticancer activity of decavanadate anion. zones). It was observed that the areas with the highest charge density were located on the decavanadate anion, while the positively charged areas were located in the pyrimidine molecules. The non-covalent interactions between the O atoms of decavanadate anion and the H atoms of pyrimidine molecules were found in intermediate electron density zones, represented in yellow-green colors.


Introduction
Polyoxidometalates (POMs) are defined as clusters made from early transition-metals, typically d 0 species V(V), Nb(V), Ta(V), Mo(VI), and W(VI), bridged by oxide anions. This class of compounds is highly interesting in molecular structural variety, reactivity, and applications in analytical chemistry, catalysis, medicine, and materials research [1]. POMs have great potential in biological applications, since every aspect that involves the interaction of POM with biological target macromolecules could be modified to improve their beneficial effects on a biological system. Thus, interesting POMs with anticancer [2] and antibiotic activities [3], among others, have been obtained to date. [2-ampymH] 6

[V 10 O 28 ]·5H 2 O (1)
The decavanadate anion is a very useful cluster in coordination chemistry as a building block for various structures. Since they have a high negative charge and oxygen-rich surfaces, interactions between this anion and metals or organic ligands are easily formed by coordination or hydrogen bonds [39]. Based on this cluster, the crystal structure of a new compound with a cationic ligand was determined in the present study. Table 1 presents the corresponding crystal data, and Tables S1-S7 in the supplementary section contain additional crystallographic information. Compound 1 crystallizes in the P1 space group of the triclinic system and consists of a three-dimensional supramolecular structure, where units of [V 10 O 28 ] 6− interact with the cationic ligand 2-aminopyrimidinium, [2-ampymH] + and crystallization water molecules through hydrogen bonds. The asymmetric unit contains half of a decavanadate anion arranged in an inversion center, three independent molecules of [2-ampymH] + (named as A, B, and C) protonated in the N3 position, and two and a half crystallization water molecules (two of them with 100% occupancy and one with occupancy set at 50%). Therefore, the (6-) charge of the decavanadate cluster is stabilized by six [2-ampymH] + cations, which are arranged on both sides of the cluster interacting by hydrogen bonds. Thus, the formula established by monocrystal X-ray diffraction was [2-ampymH] 6 [V 10 O 28 ]·5H 2 O (Figure 1).
The protonation effect of the 2-ampym ligand mainly results in a variation in the internal angles N2-C1-N3 and, to a lesser extent, in the angles C1-N2-C2 and C1-N3-C4. Table 2 shows these angles compared with the values given by J. Scheinbeim et al. for a non-protonated 2-ampym molecule [40]. Protonation also affects the distances between the atoms within the cation. Specifically, there is a shortening of the N1-C1 bond, from 1.342 Å for the ligand to 1.317 Å on average for the cation. In addition, in the 2-ampym ligand, the N2-C2 and N3-C4 bonds are symmetrical, with a value of 1.331 Å, while, N3-C4 bond is greater for the protonated form than the N2-C2 one (1.351 Å versus 1.323 Å on average, respectively). The protonation effect of the 2-ampym ligand mainly results in a variation in the internal angles N2-C1-N3 and, to a lesser extent, in the angles C1-N2-C2 and C1-N3-C4. Table 2 shows these angles compared with the values given by J. Scheinbeim et al. for a non-protonated 2-ampym molecule [40]. Protonation also affects the distances between the atoms within the cation. Specifically, there is a shortening of the N1-C1 bond, from 1.342 Å for the ligand to 1.317 Å on average for the cation. In addition, in the 2-ampym ligand, the N2-C2 and N3-C4 bonds are symmetrical, with a value of 1.331 Å, while, N3-C4 bond is greater for the protonated form than the N2-C2 one (1.351 Å versus 1.323 Å on average, respectively). The decavanadate cluster comprises ten edge-sharing VO6 octahedra, containing ten vanadium atoms and twenty-eight oxygen atoms. In total, the cluster possesses sixty V-O bonds which can be differentiated by their bond distances. First, two oxygen atoms are located inside the polyanion (Oc) bonded to six vanadium atoms, each one with the largest bond distance in the range of 2.106-2.340 Å. Four other oxygen atoms are arranged on the surface (Ob1) and coordinate to three vanadium atoms each, with bond lengths of V-Ob1 in the range of 1.900-2.084 Å. Fourteen oxygen atoms at the corners (Ob2) coordinate to two vanadium atoms each, where V-Ob2 bond length distances are ranged from 1.684 to 2.069 Å. Lastly, eight terminal oxygen atoms (Ot) are coordinated to only one vanadium each,  The decavanadate cluster comprises ten edge-sharing VO 6 octahedra, containing ten vanadium atoms and twenty-eight oxygen atoms. In total, the cluster possesses sixty V-O bonds which can be differentiated by their bond distances. First, two oxygen atoms are located inside the polyanion (O c ) bonded to six vanadium atoms, each one with the largest bond distance in the range of 2.106-2.340 Å. Four other oxygen atoms are arranged on the surface (O b1 ) and coordinate to three vanadium atoms each, with bond lengths of V-O b1 in the range of 1.900-2.084 Å. Fourteen oxygen atoms at the corners (O b2 ) coordinate to two vanadium atoms each, where V-O b2 bond length distances are ranged from 1.684 to 2.069 Å. Lastly, eight terminal oxygen atoms (O t ) are coordinated to only one vanadium each, and the V = O t distance is the smallest one with values in the range of 1.595-1.611 Å. These bond distances agree with other [V 10 O 28 ] 6− anions reported in the literature [41][42][43][44]. Bond lengths are listed in Table S3.
The supramolecular structure is dominated by hydrogen bonds of type N-H···N, N-H···O, and O-H···O, in addition to π-stacking interactions between pyrimidine rings. On the one hand, as shown in Figure 2, the decavanadate anion interacts by hydrogen bonds of type N-H···O with six [2-ampymH] + cations: two hydrogen bonds are formed thanks to N1 and N3 atoms of molecules A and C, giving rise to a R 2 2 (8) ring motif, according to the nomenclature of Etter et al. [45], while molecule B only forms one hydrogen bond through the N3 atom. Furthermore, [2-ampymH]       On the other hand, due to the arrangement of [2-ampymH] + cations, π-π stacking interactions are formed between pyrimidine rings, by forming a chain of six rings that repeats infinitely in space ( Figure 3), with an average distance between centroids of 3.617 Å, which falls within the range of a π-π interaction. Some of the rings are stacked in an aligned way, while others form a π-π stacking shifted parallel about 22 • , also producing a shorter distance between centroids because π-σ attraction forces predominate in this arrangement [46]. A similar structure has been published by S. Sedghiniya et al. [6], in which decavanadate anion interacts by hydrogen bonds with cations of Adenine, one of the five nitrogenous bases that form DNA and RNA. Similarly, in this structure, the cations interact with each other by hydrogen bonds and π-π stacking interactions, with a centroid-to-centroid distance of 3.5 Å and a displacement angle of approximately 24 • . aligned way, while others form a π-π stacking shifted parallel about 22°, also producing a shorter distance between centroids because π-σ attraction forces predominate in this arrangement [46]. A similar structure has been published by S. Sedghiniya et al. [6], in which decavanadate anion interacts by hydrogen bonds with cations of Adenine, one of the five nitrogenous bases that form DNA and RNA. Similarly, in this structure, the cations interact with each other by hydrogen bonds and π-π stacking interactions, with a centroid-to-centroid distance of 3.5 Å and a displacement angle of approximately 24°. This arrangement is clearly reminiscent of DNA structure, a polymer formed by base pairs that interact with hydrogen bonds and are stacked by π-π stacking interactions with a distance between centroids of 3.4 Å. Similarly, in Compound 1, [2-ampymH] + molecules interact by hydrogen bonds and are stacked by π-π stacking interactions with a distance between centroids of approximately 3.6 Å ( Figure 4). This arrangement is clearly reminiscent of DNA structure, a polymer formed by base pairs that interact with hydrogen bonds and are stacked by π-π stacking interactions with a distance between centroids of 3.4 Å. Similarly, in Compound 1, [2-ampymH] + molecules interact by hydrogen bonds and are stacked by π-π stacking interactions with a distance between centroids of approximately 3.6 Å ( Figure 4).  Lastly, this structure contains water molecules of crystallization ( Figure 5). Two of the water molecules of the unit cell (O1W and O3W) are arranged between the decavanadate and the cations in general positions, while the third one (O2W) is located at Lastly, this structure contains water molecules of crystallization ( Figure 5). Two of the water molecules of the unit cell (O1W and O3W) are arranged between the decavanadate and the cations in general positions, while the third one (O2W) is located at the apex, making its contribution to the unit cell 50%. These water molecules form hydrogen bonds within the structure, helping to stabilize it. All distances of hydrogen bonds along with their angles are found in Table 3. Lastly, this structure contains water molecules of crystalliz the water molecules of the unit cell (O1W and O3W) are decavanadate and the cations in general positions, while the third the apex, making its contribution to the unit cell 50%. These wate gen bonds within the structure, helping to stabilize it. All dista along with their angles are found in Table 3.

Infrared (IR) Spectroscopy
At the top of the 2-ampym ligand spectrum ( Figure 6) were the characteristic bands at 3331 and 3165 cm −1 that belong to the stretching vibration modes of the N-H and C-H bonds, respectively. In addition, the bending vibration modes appeared at 1645 and 1128 cm −1 , respectively. Lower wavenumbers show the characteristic bands of the aromatic pyrimidine ring with peaks at 1556, 1469, 1356, 794, and 555 cm −1 , belonging to the C-N, C=C, C-C, and CCC bonds, respectively [47][48][49]. As for the spectrum of Compound 1, it was possible to see both characteristic bands of the ligand 2-aminopyrimidine and the bands belonging to the different vibration modes of the V-O bond of decavanadate anion. Although the band belonging to the vibration of the O-H bond was not distinguished precisely, there was a wideband between 3400 and 2900 cm −1 , which could be assigned to water molecules in the structure. In addition, there was a slight displacement of the bands belonging to the V=O and N-H groups, so the interaction between both molecules was weak, mainly due to hydrogen bonds. Main IR bands are collected in Table 4.
belonging to the different vibration modes of the V-O bond of decavanadate anion. Although the band belonging to the vibration of the O-H bond was not distinguished precisely, there was a wideband between 3400 and 2900 cm −1 , which could be assigned to water molecules in the structure. In addition, there was a slight displacement of the bands belonging to the V=O and N-H groups, so the interaction between both molecules was weak, mainly due to hydrogen bonds. Main IR bands are collected in Table 4.

Thermal Study
The thermal decomposition behavior of [2-ampymH] 6 [V 10 O 28 ]·5H 2 O (1) in solid-state was studied by TG analysis in the range of 35-950 • C ( Figure S1). The thermogram of Compound 1 revealed a mass loss of around 5% in the range of 35-180 • C, which corresponded to five water molecules of crystallization (% mass, calc. (found) for 5 × H 2 O: 5.5% (5.4%)). Then, a continuous mass loss occurred in the range from 180 to 450 • C, which may correspond to the thermal degradation of six 2-aminopyridinium cations (% mass, calc. (found) for 6 × C 4 H 6 N 3 : 35.5% (36.3%)). The total mass loss of Compound 1 was 41.85%, which was in good agreement with the calculated value of 41%. No further losses in mass were observed above 460 • C, which can be attributed to the remaining [V 10 O 28 ] 6− fragment that may be thermally stable up to 600 • C, as previously reported [50,51]. Nevertheless, some authors point out that further thermal treatment of decavanadate anion could result in the formation of other vanadium oxides, such as V 2 O 5 or other mixed-valence oxides [52][53][54]. Figure 7 shows the molecular structure and isosurfaces of the Highest Occupied Molecular Orbital (HOMO) and Lowest Unoccupied Molecular Orbital (LUMO) molecular orbitals of Compound 1 plotted with an isovalue of 0.02 a.u. The frontier orbitals HOMO and LUMO can be related to donating or accepting electrons, respectively. The higher energy of the HOMO orbital indicates a more pronounced behavior as an electron donor, while the lower energy of the LUMO orbital indicates a higher electron affinity. The HOMO orbital, with energy of −7.4671 eV, is located on the decavanadate anion, mainly on the Obridge atoms, as shown in Figure 7b, while the LUMO orbital, with energy of −5.7381 eV, is located on the organic counterions of pyrimidine, mainly on the C and N atoms of the rings, as shown in Figure 7c.      The main non-covalent bonds between the decavanadate cluster a rings were characterized by topological electron density parameters, The main non-covalent bonds between the decavanadate cluster and pyrimidine rings were characterized by topological electron density parameters, ρ(r), the Laplacian of density, ∇ 2 ρ(r), the energy of interaction, E H···Y , and the interatomic distance, D int . The results are shown in Table 5. In addition, Figure 9 shows the molecular graphs of Compound 1. From the results, it can be seen that the positive values of ∇ 2 ρ(r) indicates that non-covalent interactions are hydrogen bonds. The interaction energy is calculated using the equation E H···Y = 1 2 |V(r)|. The maximum value was found for the interaction O10· · · H3CA with a value of 0.0788 a.u. and the highest interaction energy of 23.02 kcal mol −1 . In addition, the interaction O3· · · H3AA is strong with an energy of 11.26 kcal mol −1 . The interaction energies between pyrimidine molecules are in the range of 4.54-5.05 kcal mol −1 (N2A· · · H1BA, N2B· · · H1AA, and N2C· · · H1CA). Figure 7a,b shows the π-π stacking interactions between pyrimidine molecules. Many rings and cage critical points were also observed, which indicate the formation of a stable ring and cage-like structures that form between pyrimidine molecules and decavanadate, as well as between the central pyrimidine rings.

Theoretical Calculations
Lastly, the Hirshfeld surface of Compound 1, as shown in Figure 10a, was mapped with the normalized contact distance, d norm . It was observed that the red spots on the Hirshfeld surface were due to close intermolecular interactions of O· · · H, between oxygen atoms of [V 10 O 28 ] 6− ion and hydrogen atoms of 2-aminopyrimidine. In Figure 10b, in the fingerprint plot, d i indicates the distance from the surface to the nearest nucleus inside the surface, and d e is the distance from the surface to the nearest nucleus outside the surface. The interaction O· · · H had the most significant contribution to the overall surface with 96.5%, blue-green region of the surface. Other interactions with minor contributions were O· · · N (1.4%) and O· · · O (1.3%), which corresponded to the small grey zone of the whole surface.  Table 5. In addition, Figure 9 shows the molecular graphs of Compound 1.   surface, and de is the distance from the surface to the nearest nucleus outs The interaction O⋯H had the most significant contribution to the over 96.5%, blue-green region of the surface. Other interactions with minor con O⋯N (1.4%) and O⋯O (1.3%), which corresponded to the small grey zo surface.

Docking Analysis
Molecular docking analysis was performed to test DNA/RNA interactions with Compound 1 using DNA and several microRNAs (miRNA), such as lncRNA, miR-21, and let-7 miRNA (Figure 11). The docked binding energies and interactions with the miRNAs structures were between −8.91 and −8.39 Kcal/mol ( Table 6). The interactions involve different positions between decavanadate and miRNAs that comprise hydrogen bonds and hydrophobic interactions mainly with the nitrogen bases Guanine and Adenine followed by Cytosine and Uracil (Figure 11a-c).   The results obtained when docking was carried out considering DNA as a macromolecule showed binding energies of −10.79 kcal/mol and the interaction involved hydrogen bonds with Guanine and Adenine (Figure 11d).

Discussion
Every year, the global prevalence of cancer rises, as does the resistance to current chemotherapeutic agents like cisplatin; nowadays, one of the main goals of the pharmaceutical industry is the development of more effective drugs for the treatment of cancer [55,56].
This work presents a new structure based on the decavanadate cluster and the organic ligand, 2-aminopyrimidine. In general, pyrimidines are one of the most bioactive classes of compounds, with a wide range of biological activities, including in vitro antiviral, diuretic, antitumor, anti-HIV, and cardiovascular effects. Furthermore, pyrimidine moiety is present in the nucleobases that act as building blocks of nucleic acid, DNA or RNA, thymine, cytosine, and uracil, which might be one possible reason for the vast medicinal applications of pyrimidine-derived compounds [57]. Within the pyrimidines, the aminopyrimidine scaffold is present in the structure of a wide variety of natural and synthetic products, such as Thiamine (vitamin B1), Meridianins (a class of marine alkaloids), or Imatinib, a drug used against leukemia, which displays biological activities such as neuroprotection, antibiotic, antidiabetic, anti-Alzheimer, and anticancer. Among the large family of aminopyrimidines, the 2-aminopyrimidine isomer is the most studied, mainly because of its versatility as a starting material for synthesizing many other bioactive compounds [55].
Vanadium is an essential element, and it is important to keep the average vanadium concentration in the human body at around 0.3 µM by supplementation via food and drinking water. In the organism, vanadium is found primarily in the form of vanadate H 2 VO 4 − that, due to the structural and chemical similarity with phosphate, is most likely involved in the regulation of phosphate-dependent processes, such as metabolic pathways involving phosphatases and kinases, as well as phosphate metabolism in general. In addition, organic ligands could aid in the modulation of vanadium's bioavailability, transport, and targeting mechanism, so, nowadays, coordination compounds containing vanadium are gaining in popularity because of their potential in the treatment of diabetes and cancer, leishmaniasis, and HIV [58]. In light of this, among vanadate compounds, some researchers have pointed out decavanadates as alternative antitumor agents with promising findings in tumor growth inhibition. Although the anticancer activity of decavanadate is more recent and not yet fully understood, it is probably due to the inhibition of different enzymes such as alkaline phosphatases, ectonucleotidases, and P-type ATPases [56]. One of the first articles published about decavanadate compounds with potential antitumor activity was Na 4 Co(H 2 O) 6 V 10 O 28 ·18H 2 O reported by Zhai et al. [15]. This compound in vitro displays higher inhibitory activity to human liver cancer (SMMC-7721) and ovary cancer (SK-OV-3) cell lines than 5-fluorouracil, the antitumoral drug clinically used, while in vivo, it can decrease liver tumor mass in rats. Shortly after, in 2010, Li et al. synthesized two decavanadates compounds with organic ligands that can inhibit proliferation of human lung (A549) and murine leukemia (P388) tumor cell lines in vitro: (H 2 tmen) 3 V 10 O 28 ·6H 2 O and (H 2 en) 3 V 10 O 28 ·2H 2 O [27]. It is important to highlight that the first compound contains four methyl (-CH 3 ) substituents in the cation moiety. Their presence may enhance the lipophilic effect of this compound, which increases its penetration through the lipid bilayer of the cell membrane, thus showing higher inhibitory activity than the second compound. After that, a few more articles based on decavanadates compounds with potential antitumor activity have been published up to the present [29,53,[59][60][61].
Since hydrogen bonds naturally occur between nucleobases and are of great biological importance for DNA and RNA structures, we chose the 2-aminopyrimidine molecule to obtain a new decavanadate compound with potential anticancer activity. Similar to this organic molecule, some structures with cytosine and decavanadate have been published. In the first compound published by Bošnjaković-Pavlović et al. [36], cytosine forms dimeric cytosine-cytosinium cations that stabilize the charge of the decavanadate anion. A few years later, the second material of decavanadate with cytosine in the literature was published [37]. In this compound, all cytosine (C) molecules are protonated, and the supramolecular structure is dominated by C-C and C-Decavanadate hydrogen bonds, as well as π-stacking interaction among heterocyclic rings with a centroid-to-centroid distance similar to that found in DNA and RNA structures. The fact that in both structures, cytosines interact with one another by hydrogen bonds and π-π stacking similarly as in DNA and RNA polymers lead to two hypotheses: (i) polyoxidometalates could be employed as templates or catalysts for base-base linkages, and (ii) base-base pairing was crucial in the early stages of life [37]. Why nucleobases tend to form this arrangement in DNA and RNA has been studied by Francés-Monerris et al. [62]. The intra-and inter-strand interactions in Watson-Crick base pairs are of great importance for the thermal stability of the double-helix structure. In fact, in their experiments, upon UV light exposure at a B-DNA arrangement, C-C dimers twist towards a face-to-face arrangement to increase the π-stacking interaction and further promote the photostability of the genomic material. On the other hand, structures containing purine-based ligands show the same fact, as in the (NH 4 ) 2 (C 8 H 10 N 4 O 2 ) 4 [H 4 V 10 O 28 ]·2H 2 O compound, π-π stacking interactions with an interplanar distance of 3.38 Å exist in the purine rings of caffeine [63], and in [AdH] 6 [V 10 O 38 ]·4H 2 O, Adeninium cations form simultaneous hydrogen bond interactions in a ribbon-like geometry as well as π-π interaction between cationic ribbons [6].
All these facts support the structure of Compound 1 in this paper, since 2-aminopyrimidinium cations are disposed of in the form that they interact with each other by hydrogen bonds and π-π stacking interactions with a centroid-to centroid distance of around 3.6 Å.
Furthermore, to study the potential anticarcinogenic activity of Compound 1, docking tests were performed to analyze the interaction of the decavanadate anion with some types of microRNAs (miRNAs) and DNA molecules. MiRNAs are small, non-coding RNA molecules with a length of 19-25 nucleotides that regulate various target genes. MiRNAs play a role in the cell cycle, differentiation, proliferation, energy metabolism, and immunological response, among other biological activities, regulating around 30% of human genes, with half of these genes being tumor-related. Recent studies have found that miRNAs play an important role in cancer progression, including tumor growth, differentiation, adhesion, apoptosis, invasion, and metastasis [64,65]. Within the miRNA family, the upregulations of miR-21, as well as lncRNA (long non-coding RNA), are linked to several types of cancer such as malignant B-cell lymphoma or breast cancer, respectively [66,67]. On the other hand, let-7 miRNA (member of the family of let-7 RNA) is known as the keeper of differentiation and has also emerged as a promising therapeutic agent to treat cancer and immune responses [68,69]. Thus, it is interesting to carry out experiments to prove this hypothesis, since our preliminary analyses show great affinity of decavanadate for miRNA fragments.
In addition to the possible antitumoral activity that Compound 1 could exhibit, the interesting structure that it shows, along with all the above examples about structures with decavanadate clusters and organic ligands that arrange remembering the DNA or RNA structures, highlights the idea that these polyoxidovanadates could act as templates or catalysts for base-base pairing. Bernal, in 1949, was the first who proposed the important role of clay minerals in the origin of life [70]. The advantage of these clays could include ordered arrangement, substantial adsorption capacity, UV protection, ability to concentrate organic compounds, and potential to serve as polymerization templates so that clay minerals could have played a key role in chemical evolution and the origin of life [71]. Intercalation of decavanadate into laminar minerals has already been achieved [72], so it is possible to concentrate organic compounds with a high capacity for hydrogen bond formation and catalyze their polymerization. Therefore, it will be worthwhile to explore this idea.
Additionally, it is important to mention that recent work from Aureliano's group about melanoma anticancer activity of a variety of vanadium compounds [73] and potential anti-SARS-CoV-2 activity of vanadium compounds by Scior et al. [74], present a promising future for vanadium compounds as metallodrugs.

Materials and Methods
All starting reagents were purchased from Sigma-Aldrich (Merck, Naucalpan de Juarez, Mexico) and were used as received, without additional purification. [2-ampymH] 6

Synthesis of
A total of 0.4 mmol (0.047 gr) of ammonium metavanadate (NH 4 VO 3 ), in situ precursor of decavanadate anion, was dissolved in 10 mL of distilled water at 80 • C and stirred. After cooling to room temperature, the pH was lowered with a dissolution of concentrated HCl (37%) to reach pH = 3. In a separate glass vessel, 0.2 mmol (0.019 gr) of 2-aminopyrimidine ligand was dissolved in 4 mL of distilled water. This solution was gradually added to the orange solution of decavanadate, resulting in a pH close to 4.3. The reaction mixture was then allowed to evaporate slowly at room temperature. After a week, orange crystals of Compound 1 for single-crystal X-ray diffraction were obtained, filtered, and dried in air. Anal. Calcd. for C 12

Physicochemical Characterization
Elemental analyses (C, H, and N) were carried out on a THERMO SCIENTIFIC analyzer model Flash 2000 (Thermo Fisher Scientific S.p.A., Milan, Italy), while thermal analyses were examined with a Shimadzu TGA-50H thermogravimetric analyzer (Shimadzu Corporation, Kyoto, Japan), with a heating rate of 10 • C/min under air atmosphere in the range of 35-950 • C, both of them at the "Centro de Instrumentación Científica" (University of Granada, Granada, Spain). The IR spectra of powdered samples were recorded in the 4000-400 cm −1 region with a BRUKER TENSOR 27 FT-IR, using the OPUS program as a data collector (Bruker Corporation, Billerica, MA, USA).

Single-Crystal X-ray Diffraction
X-ray quality orange crystals of Compound 1 were obtained. The crystals were prepared under inert conditions immersed in perfluoropolyether as protecting oil for manipulation. Then, a suitable crystal was mounted on a MiTeGen Micromount TM (MiTeGen, Ithaca, NY, USA) for data collection at 300(2) K. For Compound 1, diffraction intensities were collected with a Bruker D8 Venture diffractometer (Bruker Corporation, Billerica, MA, USA), using a photon detector equipped with graphite monochromated MoKα radiation (λ = 0.71073 Å). Data reduction was performed with the APEX3 software [75] and corrected for absorption using SADABS [76]. The structure was solved by direct methods using the SHELXS-2013 program [77] and refined by full-matrix least-squares on F 2 with SHELXL-2018 [78]. The positional and anisotropic atomic displacement parameters were refined for all non-hydrogen atoms. Hydrogen atoms were located in difference Fourier maps and included as fixed contributions riding their parent atoms with isotropic thermal factors chosen as 1.2 times or 1.5 times those of their carrier atoms. One water molecule (O2W) is disordered into two dispositions, so it was refined with occupancy of 50% and anisotropic displacement parameters. OLEX2 software [79] was used as a graphical interface. The Addsym routine implemented in the program PLATON [80]

Theoretical Methodology
The structural and electronic structure of Compound 1 were calculated using the density functional theory, DFT [81]. The hybrid functional B3LYP [82] was used, using the Def2SVP basis set [83] for the atoms of C, H, O, and N; and LANL2DZ [84] for vanadium atom with an effective core potential (ECP). The frontier molecular orbitals (highest occupied molecular orbital, HOMO, and lowest unoccupied molecular orbital, LUMO) and the molecular electrostatic potential (MEP) map were analyzed. The calculations were performed with the Gaussian program16 [85], and the visualization of the results was performed with the Gaussian View 6.0.16 program [86]. Additionally, the main non-covalent interactions in Compound 1 were characterized using the atoms in molecules (AIM) approach with AIMAll software [87] and Hirshfeld surface analysis with CrystalExplorer 17.5 software [88].

Docking Analysis
Molecular docking analysis was performed with the semi-flexible methodology, where the RNAs fragments were considered as a rigid entity, while flexibility was allowed for the decavanadate. The preparation of the macromolecule and the ligand was performed through the Autodock Tools 1.5.6 software [89], which includes polar hydrogens and empirical particles of atomic charges (Gasteiger-Marsili method). Different grid box sizes were used for each RNA molecule that encloses the entire fragment: for 6PK9, 60, 106, and 60 Å were used; 70, 126, and 70 Å for 2JXV fragment; and 94, 76, and 66 Å were used for 2MNC. In addition to those fragments, one DNA structure was considered in the docking study: 1BNA, with sizes of 70, 70, and 120 Å. The grid spacing for all the docking calculations was set to the default 0.375 Å value, using the Lamarckian genetic algorithm (LGA) searching methods. The parameters for the vanadium atom were the sum of VDW radii of two similar atoms (3.14 Å), plus the VDW well depth (0.016 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 indicate the type of H-bonding atom and indexes for the generation of the autogrid map (0, -1, -1, 1, respectively).

Conclusions
A new decavanadate compound containing 2-aminopyrimidine was synthesized and characterized experimentally through elemental analysis, infrared spectroscopy, thermogravimetric analyses, and single-crystal X-ray diffraction. Using theoretical studies based on DFT calculations and AIM, the non-covalent interactions present in the compound were also studied. Frontier molecular orbital, molecular electrostatic potential, and non-covalent interactions were completely characterized. The hydrogen bond patterns of this molecule, although similar to the already reported Adenine-decavanadate compound, has the interesting feature of resembling the ladder-like structure of DNA and RNA molecules, even with similar distances from the centroid-to-centroid and π-π interactions, it also occurs in triads, which could remind codon an anticodon interactions. A docking study was set to explore more of the resemblance with DNA and RNA to see the possibility of decavanadate interacting with RNA molecules. Thus, a set of test miRNA was considered for docking. The latest development in understanding the roles of non-coding RNAs and microRNAs is worthwhile to try some fitting. Although the negative charge of these molecules could prevent interactions with the anionic decavanadate, the formation of several hydrogen bonds with the nucleobases appears to counteract the charge problem, and relatively good interaction energies were observed. At the moment, this only suggests that experiments should be designed to test whether these interactions are worthwhile to pursue. The catalytic properties to build blocks with the correct orientation and similar types of interactions could also be explored since it could shed some light on the nature of the first replicator. On the other hand, since cancer cells exhibit increased levels of mRNA translation to meet tumor growth requirements [90], it will be interesting to explore the interactions of polyoxidovanadates with mRNA. Therefore, new and fascinating decavanadate chemistry is around the corner.