Abstract
Polyoxometalates (POMs) structures have raised considerable interest for the last years in their application to biological processes and medicine. Within this area, our mini-review shows that computational modelling is an emerging tool, which can play an important role in understanding the interaction of POMs with biological systems and the mechanisms responsible of their activity, otherwise difficult to achieve experimentally. During recent years, computational studies have mainly focused on the analysis of POM binding to proteins and other systems such as lipid bilayers and nucleic acids, and on the characterization of reaction mechanisms of POMs acting as artificial metalloproteases and phosphoesterases. From early docking studies locating binding sites, molecular dynamics (MD) simulations have allowed to characterize the nature of POM···protein interactions, and to evaluate the effect of the charge, size, and shape of the POM on protein affinity, including also, the atomistic description of chaotropic character of POM anions. Although these studies rely on the interaction with proteins and nucleic acid models, the results could be extrapolated to other biomolecules such as carbohydrates, triglycerides, steroids, terpenes, etc. Combining MD simulations with quantum mechanics/molecular mechanics (QM/MM) methods and DFT calculations on cluster models, computational studies are starting to shed light on the factors governing the activity and selectivity for the hydrolysis of peptide and phosphoester bonds catalysed by POMs.
Introduction
Polyoxometalates (POMs) are a unique class of well−defined polynuclear metal oxide clusters that are usually built up from early transition metal ions such as W, Mo or V in their highest oxidation state with an overwhelming diversity in size, composition and structures. They have been widely applied in catalysis due to their tunable Brønsted acidity combined with redox properties, as well as, their ability to accommodate other transition metals (TMs) in their structure. (Wang and Yang, 2015). In addition, polyoxometalates are water-soluble, large anions that give rise to unusual solution behavior that have found applications in biotechnology and supramolecular chemistry. (; ; ). Recently, this behavior has been attributed to the (super)chaotropic character of the POMs anions, whose low charge density make them to be weakly hydrated, and consequently, exhibit propensity to assemble with organic moieties and biomolecules (see below for a deeper description). (). Among the applications in biotechnology and medicine, POMs have shown in vitro and in vivo antiviral, antibacterial or antitumor properties; utility in protein crystallography; or activity as artificial metalloproteases and phosphoesterases. (), (; ; ; ) (; ; ; ).
The bioactivity of POMs depends largely on their ability to establish specific interactions with biomolecules, and therefore, the precise understanding of these interactions is crucial for further developments. Moreover, to exploit the use of metal-substituted POMs as a novel class of artificial enzymes, it is necessary to progress in characterizing and rationalizing their mechanism of action. Previous reviews in computational POM chemistry have mainly collected studies on structure, electronic properties, spectroscopy, and reactivity based on Density Functional Theory methods (DFT) in conjuction with continuous solvent models. (; ). More recently, the development of tailor-made, classical potentials for POMs has allowed to perform dynamic simulations, such as Molecular Dynamics (MD), on complex (bio)molecular systems including dynamic properties and explicit solvent effects. Here, we focus on how recent computational studies have contributed to increase understanding of the physicochemical foundations underlying the basic principles of the observed biological activity of POMs.
Interaction Between Polyoxometalates and Biomolecules
The interaction of POMs with proteins have attracted attention of many scientists because it plays a crucial role in mechanistic pathways governing the biological activity exhibited by POMs. (). Although covalent interaction between biomolecules and POMs have been reported, (), most of these interactions are of non-bonding nature. Experimentally, there are limitations to identify the binding modes and characterise the nature of these interaction, as well as to assess the competition for binding sites between POMs and other species in the media. Computationally, early docking studies explored the binding locations showing that POMs interact mainly at positively charged patches of the protein, where cationic- and polar-type amino acids predominate. (; ; ; ; ; ; ; ). Owing to the intrinsic limitations of docking methods, atomistic molecular dynamics (MD) simulations have been more recently performed to reveal the driving forces that are responsible for the specific interactions (Figure 1A)[ (; ; ; ; )] Pioneer MD simulations analysed the interaction between model protein hen egg-white lysoszyme (HEWL) and three different POMs, the Ce-substituted Keggin-type anion [PW11O39Ce(OH2)4]3− the corresponding 1:2 dimer [Ce(PW11O39)2]10− and the Zr-substituted Lindqvist-type anion [W5O18Zr(OH2) (OH)]3−, which differ in the overall charge, the size, the shape and the type of substituted metal. (). All POM structures interacted preferentially with positively charged Lys- and Arg-rich patches on the protein surface. The nature of these interactions comprises mainly electrostatic attraction, hydrogen bonding and water-mediated interactions, not only with positively-charged amino acids such as lysine and arginine, but also with uncharged polar amino acids such as tyrosine, serine and asparagine (Figure 1B). The basic oxygen atoms of POM framework interact with the side chains of the amino acids and, in lesser extent, with the N-H amide group of the main protein chain. Moreover, depending on the size and shape of the POM, several amino acids can interact simultaneously with the oxide framework, favoring the formation of POM-protein complexes. These results were subsequently backed by X-ray structural characterization of non-covalent complexes between HEWL and several transition metal-substituted tungstates. This also suggests that the interaction is largely independent on the nature of substituted metal within the same polyoxometalate structure type because the interaction occurs through the POM oxide framework. (; Vandebroek et al., 2019). Interactions of the same nature were computationally characterised for the binding of Zr-substituted Keggin tungstates to human serum albumin (HSA) protein. ().
FIGURE 1
The attempts to set structure-activity relationships for the affinity of POMs towards proteins, and then relate them to their biological activity, have some limitations from the experimental point of view. For example, they include the lack of stability under experimental conditions, or the limited number of structures available making them to differ in more than one feature at a time. On the other hand, computational modelling allows performing systematic variations on single parameters of well-defined POM structures. Thus, a fundamental study based on MD simulations and descriptor-based modeling has been recently reported, in which the charge or the size and shape of a series of tungstate POMs were modified systematically and the affinity towards HEWL protein evaluated. (). Using two molecular descriptors that account for the charge density (charge per metal atom ratio; q/M) and the size and shape (shape-weighted volume; VS), it was possible to build quantitative multidimensional regression models for protein affinity with predictive ability. Interestingly, the model reveals the non-linear relationship between protein affinity and both the charge density and the size of the POM, as a result of a delicate balance between POM···protein and POM···solvent interactions. Atomistic simulations revealed the variation of hydrogen bonding patterns for POM···water interactions depending on POM charge density. Moderately charged POMs anions have chaotropic character (water-structure-breaking), while highly charged POMs anions have kosmotropic properties (water-structure-forming) which results in proportionally larger desolvation energies, and consequently, less affinity towards proteins (Figure 1C). In fact, Nau and co-workers have defined POMs as superchaotropic anions that extent beyond the classical Hofmeister scale; that is, low charge density anions that are weakly hydrated and consequently have propensity to assemble with organic moieties and biomolecules. (). These simulations also showed that POM···protein interactions are size-specific, and that Keggin-type anions have the optimal size and shape to fit the cationic sites of HEWL. ().
More recently, computational studies based on MD simulations have also analysed the interaction of polyoxoniobates and polyoxovanadates with different proteins. (), () The nature of the characterised interactions are similar to those of tungstates and transition metal-substituted tungstates with the POM···protein binding dominated by electrostatic and hydrogen bonding forces. Experimentally, the analysis of X-ray polyoxovanadate-protein structures show that the binding sites include a variety of positively charged amino acids such as Arg, His and Lys. (). Computationally, a comparison of decavanadate anion (V10O286−) with the isostructural and equally charged decaniobate (Nb10O286−) was carried out in their interplay with the globular actin protein. (). Interestingly, they prefer different binding sites of the protein (the catalytic nucleotide site α for V10O286− and the β site for Nb10O286−), both inducing conformational arrangements in the protein, suggesting that biological activity could be synergistic. MD calculations confirmed that polyoxoniobates [Nb10O28]6− and [TiNb9O28]7− interact with positively charged sites of the surface of native S100A9 protein, which is a pro-inflammatory and amyloidogenic protein involved in neurodegenerative diseases, interacting simultaneously with several amino acids at a highly dynamic part of the protein. (). Experimentally, these two POMs acted as inhibitors of S100A9 amyloid assembly. (). Additionally, the same authors assessed how the ionic strength of the media influences the complex formation. Increasing the NaCl salt concentration from 20 to 150 mM in the simulated system reduced the formation of POM-protein complexes.
The interaction of POMs with protein surfaces can induce changes on its structure. For example the dimeric 1:2 POM K16 [Hf(α2-P2W17O61)2] co-crystalise with HEWL in its monomeric form, which was never observed in water solution. (Vandebroek et al., 2018). This indicates that the dimeric POM dissociates upon binding because it moves from the highly polar bulk water to the protein surface which has a lower polarity (lower dielectric constant). DFT calculations with continuum solvent model evaluated the free energy cost of dissociation process at different dielectric constants that is unfavorable at the dielectric constant of bulk water (εr = 80) and it becomes favorable at lower dielectric constants (εr = 20–50). (Vandebroek et al., 2018). Thus, theoretical calculations support the protein-assisted dissociation of group IV transition metal-substituted dimeric structures. It is to note that computational methods have studied other factors influencing the dissociation of these dimers such as the pH or concentration (; ). There have been also several molecular modeling studies on POM-oligopeptide hybrids, in which the polypeptide is covalently linked to the polyoxometalate. (Vilona et al. 2018; ). Besides their own interest as potential drugs or as building blocks for self-assembled materials, these compounds can serve as more tractable structures, allowing to combine classical MD with quantum mechanics DFT calculations, in order to understand the specific interactions of peptides with the metal-oxide surface of the POM. For example, in tin-substituted Dawson polyoxotungstates, the polyglycine side chains folds towards the metal-oxide surface forming zipper hydrogen bond networks. (Vilona et al., 2018). Interestingly, the intramolecular hydrogen bonds are formed preferentially with the terminal oxygens (W=O), even that they are less basic than the W-O-W bridging μ2 oxygens. Even classical molecular dynamics simulations present limitations since they cannot handle very big and complex systems such as lipid bilayers, and therefore this type of systems have been scarcely studied. Nevertheless, by sacrificing molecular details, coarse-grained MD simulations analysed the embedding of a giant POM (Mo132 type Keplerate capsule) in lipid bilayer membranes. (). The simulated system remained stable, and water was observed to flow into and out of the capsule as well as Na+ cations, suggesting that Mo132 can form a functional synthetic ion channel.
Polyoxometalates as Artificial Enzymes
The ability of POMs to form specific interactions with biomolecules combined with their capacity to catalyse reactions have prompted their use as artificial metalloenzymes. Computationally, two main processes have been analysed: 1) the peptide bond hydrolysis in di-, oligopeptides and proteins with potential applicability to proteomics, and 2) the phosphoester bond hydrolysis.
Mechanism of Peptide Bond Hydrolysis
The group of Parac−Vogt tested a series of TM-substituted POMs as catalysts in the hydrolysis of peptide bonds in small dipeptides, and among them, only ZrIV−, CeIV− and HfIV−substituted POMs showed hydrolytic activity. (). Since then, these type of POMs have been successfully applied to the hydrolysis of dipeptides and oligopeptides, as well as, to the selective hydrolysis of a wide range of proteins (for a recent review see ref ()). Interestingly, different patterns of hydrolysed sites were observed when moving from one protein to another. Computational studies have focus on the characterization of the reaction mechanism and the rationalization of the observed selectivity. (; ; ; ). In all cases, the reaction initiates by coordinating the hydrolytically active metal ion to the amide oxygen atom of the peptide bond. This coordination polarizes carbonyl group, owing the Lewis acid nature of these metal ions and renders the carbon atom more susceptible to suffer a nucleophilic attack. Then, the mechanism can be classified depending on which oxygen fragment is responsible of the nucleophilic attack (Figure 2A): 1) a hydroxo ligand of the TM via inner-sphere attack, 2) an external water molecule via outer-sphere attack assisted by the TM-OH moiety, or 3) by the carboxylate group of a neighbour amino acid (mechanism iii), and 4) a carboxylate group nearby via direct attack. Finally, the protonation of the amide nitrogen induces the C-N bond cleavage.
FIGURE 2
Pierloot and coworkers characterised the mechanism for the hydrolysis of several dipeptides catalysed by a tetrazirconium−containing sandwich POM by means of DFT calculations. (). The authors found that the most likely mechanism involves the nucleophilic attack of a solvent water molecule to the carbon of the peptide bond assisted by the carboxylate group of the C−terminus end acting as a Brönsted base (mechanism type iii). Note that when the dipeptide is embedded in the protein main chain, the carboxyl terminus group is not necessarily available to assist the outer-sphere mechanism. Prabhakar and coworkers analysed the hydrolysis mechanism for different peptide bonds in human serum albumin (HSA) protein catalysed by [Zr(OH) (PW11O39)]4− anion using static quantum mechanics/molecular mechanics (QM/MM) calculations in the full molecular system. (). The authors found that the inner-sphere nucleophilic attack of the Zr-OH ligand (mechanism i) is somewhat preferred with respect to the outer-sphere attack of an external water molecule assisted by the Zr-OH moiety (mechanism ii) for the Cys392-Glu393 bond. Later, the observed hydrolysis at the Asn44-Arg45 site in the HEWL protein by [W5O18Zr(H2O) (OH)]3− Lindqvist anion was analysed in detail by means of a combination of DFT calculations on cluster models obtained from molecular dynamics (MD) simulations, QM/MM calculations, and metadynamics simulations. (). Calculations on solvated full-protein model suggests that both mechanism i and ii are viable under experimental conditions, since the large radius of the Zr atom allows the strained four-membered ring transition state of mechanism i to be close in energy to the six-membered ring transition state of mechanism ii. Constrained MD simulations on the coordinated POM-Zr···HEWL complex at Asn44-Arg45 site can discard the carboxylate-assisted mechanisms iii and iv because the closest COO− moieties are far away from the hydrolysed peptide bond. Conversely, for hemoglobin (Hb), calculations support mechanism iii at peptide bonds adjacent to an aspartate residue (Asp—X). (). Interestingly, in protein hydrolysis, the C-N scission step governs the overall catalytic activity, (; ; ), while in dipeptides, the initial nucleophilic attack was found to be the rate−determining step. (). Probably, this is a consequence of the constrains imposed by the main chain of the protein to the peptide bond.
Early attempts to rationalise the observed selectivity pointed out that the cleavage bonds were located at the vicinity of positively charged patches on the protein surface, and that co-crystals confirmed that POM structures bind preferentially to these sites. (Vandebroek et al., 2018). Computational studies have helped to understand the correlation between the POM···protein interactions and the selectivity exhibited by Zr-substituted POMs, as well as, to identify other important factors. (; ). For Hb protein, the key factor to explain the selectivity towards Asp-X sites is the mechanism, in which direct intramolecular attack of the Asp COO− group on the amide carbon atom is optimal for Asp-X as compared to related X—Asp, X—Glu and Glu—X bonds. (). To explain that only the Asp-X bonds located at random coil fragments of the protein are hydrolysed, the authors suggested that less rigid parts of the protein backbone should be more amenable to electrophilic activation by the Zr ion. For HEWL protein, constrained MD simulations on POM coordination to reactive and non-reactive sites indicate that the origin of the selectivity can be ascribed to an enzyme-like recognition process. (). At the reactive site, the positively charged and polar amino acids traps the POM more strongly reducing the energy expense for reorientation and coordination, shifting down the energy profile, and accelerating the hydrolysis rate (Figure 2B). For HSA protein, the comparison of chemically equivalent sites (Arg-Leu) shows strong non-bonding POM···protein interactions in all of them, while the secondary structure of the protein chain favours the hydrolysis at less rigid, coil regions, where lower overall free-energy barriers were calculated using DFT calculations in cluster models ().
Mechanism of Phosphoester Bond Hydrolysis
The phosphoester bond hydrolysis catalysed by molybdate anions as artificial phosphoesterases has been also studied computationally. (; ; ). In combination with electrospray ionization mass spectroscopy (ESI-MS) experiments, the obtained results provided a different picture of the mechanism compared to previous experimental studies based on NMR, RAMAN and UV-Vis spectroscopy, obtained for [Mo7O24]6− species. (). It was found that a Mo-oxo binuclear species, generated in situ, is the promoter of the important catalytic effect when using the pNPP model molecule as substrate (see Figure 2C). Another interesting finding was that starting either from any simple, mononuclear Mo oxide or from polyoxmetalate [Mo7O24]6− anion, the system can converge to the same binuclear species formed in situ [Mo2O8H4]0, which promotes reduction of the energy barrier for the phosphoester bond hydrolysis. This barrier is presumably enhanced by the Coulombic repulsion between the negative phosphate environment of the substrate and pair of electrons or the eventual negative charge of the nucleophile, a water molecule or OH−. In binuclear [Mo2O8H4]0 species the charge is null what would facilitate the approximation of the catalyst to the negative charged phosphate fragment, reducing the energy expense to reach the transition state. This catalytic process shows a low energy barrier as compared to any of the proposed classical mechanisms, substate-assisted and solvent-assisted, for the non-catalysed mechanism. (; ). At this point, it must be mentioned that some hydrogen bond of the bridge OH of the [Mo2O8H4]0 species with the pNPP substrate seems to be involved in the stabilization of the transition state.
The mechanism characterised for the binuclear [Mo2O8H4]0 species generated in situ from the [Mo7O24]6− was also investigated for the [W7O24]6− species acting as catalyst. For the tungstate, it was observed that not only the activation barriers for the hydrolytic process are higher than those for the Mo counterparts but also the obtained products are less stable from a thermodynamic point of view. This behaviour for the hydrolysis of the phosphoester bond in the presence of the W-oxo species comes in hand with the previous experimental results () in which the authors did not observe any catalytic activity in the presence of the [W7O24]6− analogue.
Conclusion
This comprehensive review shows that since the early docking studies locating the binding sites of POMs on the surface of proteins, there has been a considerable progress in the computational analysis of the interactions between POMs and biological systems. Incorporating a range of computational tools such as MD simulations, QM/MM and QM/MD methods including metadynamics simulations, or DFT calculations on cluster models, researchers have made possible to provide atomistic description of the binding of POMs to biomolecules, and mechanistic insight into the hydrolysis of peptide and phosphoester bonds with POMs acting as metalloenzymes. Thus, the nature of non-bonding POM···protein interactions has been characterised, showing that the protein affinity depends on charge, size, shape of the POM, as a result of a delicate balance with POM···solvent interactions. Simulations have identified the preferred binding sites for several proteins, reveling in some cases that the specific interaction can be a function of POM composition. By sacrificing molecular details, coarse-grained MD simulations were able to analyse the interaction of POMs with more complex biological systems such as lipid bilayers.
Understanding the factors which govern the activity and selectivity of processes with POMs acting as metalloenzymes is more challenging. However, plausible reaction mechanisms have been proposed for the hydrolysis of peptide and phosphoester bonds catalysed by Zr-substituted POMs and molybdates, respectively. Moreover, in the former case, three factors influencing the selectivity were reported: the specific nature of the dipeptide bond, the secondary structure of peptide chain, and the strong electrostatic-type POM···protein interactions. Although the number of examples is still limited, we are confident that in the coming years the computational studies on the biological activity of POMs will grow significantly becoming an important subarea of computational bioinorganic chemistry. We expect that it will grow the interest for exploring the interaction with biomolecules not only of classical early-transition metal based POMs (Mo, W, or V) but also of the so-called noble POMs including less toxic metals such as Au or Pt. Computational studies will expand to analyse the interplay with other biomolecules such as carbohydrates, steroids, triglicerides, etc. Another important topic for biomedicine is the selective affinity of POMs towards one or another biomolecule that would have consequences on their use as eventual drugs targeting specific molecules.
Statements
Author contributions
All authors listed have made a substantial, direct, and intellectual contribution to the work and approved it for publication.
Funding
Grant PGC 2018-100780-B-I00 funded by MCIN/AEI/ 10.13039/501100011033 and by “ERDF A way of making Europe.” We also thank the Generalitat de Catalunya (No. 2017SGR629), and the Fundação para a Ciência e a Tecnologia (Nos. PTDC/QUI-QFI/29236/2017, UIDB/04046/2020, and UIDP/04046/2020) for generous support. AG is thankful to ARAID—Fundación Agencia Aragonesa para la Investigación y el Desarrollo for current funding in the frame of ARAID researcher.
Acknowledgments
AG is very grateful to Frederico F. Martins, Jose Lanuza, Nuno A. B. Bandeira, Ángel Sánchez-González, Haralampos N. Miras and Prof. Xabier Lopez for fruitful discussions on the topic and the efforts working together within the project.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1
AbsillisG.Parac-VogtT. N. (2012). Peptide Bond Hydrolysis Catalyzed by the Wells-Dawson Zr(α2-P2w17o61)2 Polyoxometalate. Inorg. Chem.51, 9902–9910. 10.1021/ic301364n
2
ArefianM.MirzaeiM.Eshtiagh-HosseiniH.FronteraA. (2017). A Survey of the Different Roles of Polyoxometalates in Their Interaction with Amino Acids, Peptides and Proteins. Dalton Trans.46, 6812–6829. 10.1039/c7dt00894e
3
AssafK. I.NauW. M. (2018). The Chaotropic Effect as an Assembly Motif in Chemistry. Angew. Chem. Int. Ed.57, 13968–13981. 10.1002/anie.201804597
4
AurelianoM.GumerovaN. I.SciortinoG.GarribbaE.McLauchlanC. C.RompelA.et al (2022). Polyoxidovanadates' Interactions with Proteins: An Overview. Coord. Chem. Rev.454, 214344. 10.1016/j.ccr.2021.214344
5
AurelianoM.GumerovaN. I.SciortinoG.GarribbaE.RompelA.CransD. C. (2021). Polyoxovanadates with Emerging Biomedical Activities. Coord. Chem. Rev.447, 214143. 10.1016/j.ccr.2021.214143
6
BijelicA.AurelianoM.RompelA. (2019). Polyoxometalates as Potential Next‐Generation Metallodrugs in the Combat against Cancer. Angew. Chem. Int. Ed.58, 2980–2999. 10.1002/anie.201803868
7
BijelicA.AurelianoM.RompelA. (2018). The Antibacterial Activity of Polyoxometalates: Structures, Antibiotic Effects and Future Perspectives. Chem. Commun.54, 1153–1169. 10.1039/c7cc07549a
8
CarrR.WeinstockI. A.SivaprasadaraoA.MüllerA.AksimentievA. (2008). Synthetic Ion Channels via Self-Assembly: A Route for Embedding Porous Polyoxometalate Nanocapsules in Lipid Bilayer Membranes. Nano Lett.8, 3916–3921. 10.1021/nl802366k
9
CartuyvelsE. (2008). Chemical Studies of Anti-tumor Active Polyoxomolybdate Complexes (Leuven: KU). PhD Thesis.
10
ChaudharyH.IashchishynI. A.RomanovaN. V.RambaranM. A.MusteikyteG.SmirnovasV.et al (2021). Polyoxometalates as Effective Nano-Inhibitors of Amyloid Aggregation of Pro-inflammatory S100A9 Protein Involved in Neurodegenerative Diseases. ACS Appl. Mater. Inter.13, 26721–26734. 10.1021/acsami.1c04163
11
DuarteF.ÅqvistJ.WilliamsN. H.KamerlinS. C. L. (2015). Resolving Apparent Conflicts between Theoretical and Experimental Models of Phosphate Monoester Hydrolysis. J. Am. Chem. Soc.137, 1081–1093. 10.1021/ja5082712
12
DuarteF.BarrozoA.ÅqvistJ.WilliamsN. H.KamerlinS. C. L. (2016). The Competing Mechanisms of Phosphate Monoester Dianion Hydrolysis. J. Am. Chem. Soc.138, 10664–10673. 10.1021/jacs.6b06277
13
HuD.ShaoC.GuanW.SuZ.SunJ. (2007). Studies on the Interactions of Ti-Containing Polyoxometalates (POMs) with SARS-CoV 3CLpro by Molecular Modeling. J. Inorg. Biochem.101, 89–94. 10.1016/j.jinorgbio.2006.08.013
14
Jayasinghe-ArachchigeV. M.HuQ.SharmaG.PaulT. J.LundbergM.QuiñoneroD.et al (2019). Hydrolysis of Chemically Distinct Sites of Human Serum Albumin by Polyoxometalate: A Hybrid QM/MM (ONIOM) Study. J. Comput. Chem.40, 51–61. 10.1002/jcc.25528
15
Jiménez-LozanoP.CarbóJ. J.ChaumontA.PobletJ. M.Rodríguez-ForteaA.WipffG. (2014). Nature of Zr-Monosubstituted Monomeric and Dimeric Polyoxometalates in Water Solution at Different pH Conditions: Static Density Functional Theory Calculations and Dynamic Simulations. Inorg. Chem.53, 778–786. 10.1021/ic401999r
16
Jiménez-LozanoP.Solé-DauraA.WipffG.PobletJ. M.ChaumontA.CarbóJ. J. (2017). Assembly Mechanism of Zr-Containing and Other TM-Containing Polyoxometalates. Inorg. Chem.56, 4148–4156. 10.1021/acs.inorgchem.7b00096
17
JuddD. A.NettlesJ. H.NevinsN.SnyderJ. P.LiottaD. C.TangJ.et al (2001). Polyoxometalate HIV-1 Protease Inhibitors. A New Mode of Protease Inhibition. J. Am. Chem. Soc.123, 886–897. 10.1021/ja001809e
18
LanuzaJ.Sánchez−GonzálezÁ.BandeiraN. A. G.LopezX.GilA. (2021). Mechanistic Insights into Promoted Hydrolysis of Phosphoester Bonds by MoO2Cl2(DMF)2. Inorg. Chem.60, 11177–11191. 10.1021/acs.inorgchem.1c01088
19
LokerenL. V.CartuyvelsE.AbsillisG.WillemR.Parac-VogtT. N. (2008). Phosphoesterase Activity of Polyoxomolybdates: Diffusion Ordered NMR Spectroscopy as a Tool for Obtaining Insights into the Reactivity of Polyoxometalate Clusters. Chem. Commun.2008, 2774–2776. 10.1039/b802671h
20
LópezX.CarbóJ. J.BoC.PobletJ. M. (2012). Structure, Properties and Reactivity of Polyoxometalates: a Theoretical Perspective. Chem. Soc. Rev.41, 7537–7571. 10.1039/C2CS35168D
21
LópezX.MiróP.Rodríguez-ForteaA.BoC.PobletJ. M. (2011). Current Trends in the Computational Modelling of Polyoxometalates. Theor. Chem. Acc.128, 393–404. 10.1007/s00214-010-0820-9
22
LyH. G. T.MihaylovT.AbsillisG.PierlootK.Parac-VogtT. N. (2015). Reactivity of Dimeric Tetrazirconium(IV) Wells-Dawson Polyoxometalate toward Dipeptide Hydrolysis Studied by a Combined Experimental and Density Functional Theory Approach. Inorg. Chem.54, 11477–11492. 10.1021/acs.inorgchem.5b02122
23
LyH. G. T.MihaylovT. T.ProostP.PierlootK.HarveyJ. N.Parac‐VogtT. N. (2019). Chemical Mimics of Aspartate‐Directed Proteases: Predictive and Strictly Specific Hydrolysis of a Globular Protein at Asp−X Sequence Promoted by Polyoxometalate Complexes Rationalized by a Combined Experimental and Theoretical Approach. Chem. Eur. J.25, 14370–14381. 10.1002/chem.201902675
24
MartinsF. F.Sánchez‐GonzálezÁ.LanuzaJ.MirasH. N.LopezX.BandeiraN. A.et al (2021). Probing the Catalytically Active Species in POM‐Catalysed DNA‐Model Hydrolysis. Chem. Eur. J.27, 8977–8984. 10.1002/chem.202004989
25
MolitorC.BijelicA.RompelA. (2016). In Situ formation of the First Proteinogenically Functionalized [TeW6O24O2(Glu)]7− Structure Reveals Unprecedented Chemical and Geometrical Features of the Anderson-type Cluster. Chem. Commun.52, 12286–12289. 10.1039/c6cc07004c
26
NarasimhanK.PillayS.Bin AhmadN. R.BikadiZ.HazaiE.YanL.et al (2011). Identification of a Polyoxometalate Inhibitor of the DNA Binding Activity of Sox2. ACS Chem. Biol.6, 573–581. 10.1021/cb100432x
27
NikoloudakisE.KarikisK.LauransM.KokotidouC.Solé-DauraA.CarbóJ. J.et al (2018). Self-assembly Study of Nanometric Spheres from Polyoxometalate-Phenylalanine Hybrids, an Experimental and Theoretical Approach. Dalton Trans.47, 6304–6313. 10.1039/c8dt00380g
28
PaulT. J.Parac-VogtT. N.QuiñoneroD.PrabhakarR. (2018). Investigating Polyoxometalate-Protein Interactions at Chemically Distinct Binding Sites. J. Phys. Chem. B122, 7219–7232. 10.1021/acs.jpcb.8b02931
29
PessoaJ. C.SantosM. F. A.CorreiaI.SannaD.SciortinoG.GarribbaE. (2021). Binding of Vanadium Ions and Complexes to Proteins and Enzymes in Aqueous Solution. Coord. Chem. Rev.449, 214192. 10.1016/j.ccr.2021.214192
30
PezzaR. J.VillarrealM. A.MontichG. G.ArgañaC. E. (2002). Vanadate Inhibits the ATPase Activity and DNA Binding Capability of Bacterial MutS. A Structural Model for the Vanadate-MutS Interaction at the Walker A Motif. Nucleic Acids Res.30, 4700–4708. 10.1093/nar/gkf606
31
PrudentR.MoucadelV.LaudetB.BaretteC.LafanechèreL.HasenknopfB.et al (2008). Identification of Polyoxometalates as Nanomolar Noncompetitive Inhibitors of Protein Kinase CK2. Chem. Biol.15, 683–692. 10.1016/j.chembiol.2008.05.018
32
PrudentR.SautelC. F.CochetC. (2010). Structure-based Discovery of Small Molecules Targeting Different Surfaces of Protein-Kinase CK2. Biochim. Biophys. Acta (Bba) - Proteins Proteomics1804, 493–498. 10.1016/j.bbapap.2009.09.003
33
Sánchez-GonzálezA.BandeiraN. A. G.Ortiz de LuzuriagaI.MartinsF. F.ElleuchiS.JarrayaK.et al (2021). New Insights on the Interaction of Phenanthroline Based Ligands and Metal Complexes and Polyoxometalates with Duplex DNA and G-Quadruplexes. Molecules26, 4737. 10.3390/molecules26164737
34
SapA.De ZitterE.Van MeerveltL.Parac-VogtT. N. (2015). Structural Characterization of the Complex between Hen Egg-White Lysozyme and ZrIV-Substituted Keggin Polyoxometalate as Artificial Protease. Chem. Eur. J.21, 11692–11695. 10.1002/chem.201501998
35
SarafianosS. G.KortzU.PopeM. T.ModakM. J. (1996). Mechanism of Polyoxometalate-Mediated Inactivation of DNA Polymerases: an Analysis with HIV-1 Reverse Transcriptase Indicates Specificity for the DNA-Binding Cleft. Biochem. J.319, 619–626. 10.1042/bj3190619
36
SciortinoG.AurelianoM.GarribbaE. (2021). Rationalizing the Decavanadate(V) and Oxidovanadium(IV) Binding to G-Actin and the Competition with Decaniobate(V) and ATP. Inorg. Chem.60, 334–344. 10.1021/acs.inorgchem.0c02971
37
Solé-DauraA.GoovaertsV.StroobantsK.AbsillisG.Jiménez-LozanoP.PobletJ. M.et al (2016). Probing Polyoxometalate–Protein Interactions Using Molecular Dynamics Simulations. Chem. Eur. J.22, 15280–15289. 10.1002/chem.201602263
38
Solé-DauraA.PobletJ. M.CarbóJ. J. (2020a). Structure–Activity Relationships for the Affinity of Chaotropic Polyoxometalate Anions towards Proteins. Chem. Eur. J.26, 5799–5809. 10.1002/chem.201905533
39
Solé-DauraA.Rodríguez-ForteaA.PobletJ. M.RobinsonD.HirstJ. D.CarbóJ. J. (2020b). Origin of Selectivity in Protein Hydrolysis by Zr(IV)-Containing Metal Oxides as Artificial Proteases. ACS Catal.10, 13455–13467. 10.1021/acscatal.0c02848
40
TiagoT.MartelP.Gutiérrez-MerinoC.AurelianoM. (2007). Binding Modes of Decavanadate to Myosin and Inhibition of the Actomyosin ATPase Activity. Biochim. Biophys. Acta (Bba) - Proteins Proteomics1774, 474–480. 10.1016/j.bbapap.2007.02.004
41
van RompuyL. S.Parac-VogtT. N. (2019). Interactions between Polyoxometalates and Biological Systems: from Drug Design to Artificial Enzymes. Curr. Opin. Biotechnol.58, 92–99. 10.1016/j.copbio.2018.11.013
42
VandebroekL.De ZitterE.LyH. G. T.ConićD.MihaylovT.SapA.et al (2018). Protein‐Assisted Formation and Stabilization of Catalytically Active Polyoxometalate Species. Chem. A Eur. J24, 10099–10108. 10.1002/chem.201802052
43
VandebroekL.MampaeyY.AntonyukS.Van MeerveltL.Parac‐VogtT. N. (2019). Noncovalent Complexes Formed between Metal‐Substituted Polyoxometalates and Hen Egg White Lysozyme. Eur. J. Inorg. Chem.2019, 506–511. 10.1002/ejic.201801113
44
VilonaD.LachkarD.DumontE.LelliM.LacôteE. (2018). Elucidation of the Conformation of Polyglycine Organo-Polyoxotungstates: Evidence for Zipper Folding. Chem. Eur. J.23, 13323–13327. 10.1002/chem.201703509
45
WangS.-S.YangG.-Y. (2015). Recent Advances in Polyoxometalate-Catalyzed Reactions. Chem. Rev.115, 4893–4962. 10.1021/cr500390v
Summary
Keywords
polyoxometalate, computational chemistry, protein, artificial enzymes, molecular dynamics, DFT, peptide hydrolysis, phosphoester hydrolysis
Citation
Gil A and Carbó JJ (2022) Computational Modelling of the Interactions Between Polyoxometalates and Biological Systems. Front. Chem. 10:876630. doi: 10.3389/fchem.2022.876630
Received
15 February 2022
Accepted
21 March 2022
Published
14 April 2022
Volume
10 - 2022
Edited by
Manuel Aureliano, University of Algarve, Portugal
Reviewed by
Eugenio Garribba, University of Sassari, Italy
Nadiia I. Gumerova, University of Vienna, Austria
Updates
Copyright
© 2022 Gil and Carbó.
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: Adrià Gil, adria.gil.mestres@csic.es; Jorge J. Carbó, j.carbo@urv.cat
† These authors have contributed equally to this work
This article was submitted to Inorganic Chemistry, a section of the journal Frontiers in Chemistry
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.