Hydrolysis of chemically distinct sites of human serum albumin by polyoxometalate: A hybrid QM/MM (ONIOM) study

Mechanism of Nucleic Acid Phosphodiester Bond Cleavage by Human Endonuclease V: MD and QM/MM Calculations Reveal a Versatile Metal Dependence

Human endonuclease V (EndoV) catalytically removes deaminated nucleobases by cleaving the phosphodiester bond as part of RNA metabolism. Despite being implicated in several diseases (cancers, cardiovascular diseases, and neurological disorders) and potentially being a useful tool in biotechnology, details of the human EndoV catalytic pathway remain unclear due to limited experimental information beyond a crystal structure of the apoenzyme and select mutational data. Since a mechanistic understanding is critical for further deciphering the central roles and expanding applications of human EndoV in medicine and biotechnology, molecular dynamics (MD) simulations and quantum mechanics/molecular mechanics (QM/MM) calculations were used to unveil the atomistic details of the catalytic pathway. Due to controversies surrounding the number of metals required for nuclease activity, enzyme-substrate models with different numbers of active site metals and various metal-substrate binding configurations were built based on structural data for other nucleases. Subsequent MD simulations revealed the structure and stability of the human EndoV-substrate complex for a range of active site metal binding architectures. Four unique pathways were then characterized using QM/MM that vary in metal number (one versus two) and modes of substrate coordination [direct versus indirect (water-mediated)], with several mechanisms being fully consistent with experimental structural, kinetic, and mutational data for related nucleases, including members of the EndoV family. Beyond uncovering key roles for several active site amino acids (D240 and K155), our calculations highlight that while one metal is essential for human EndoV activity, the enzyme can benefit from using two metals due to the presence of two suitable metal binding sites. By directly comparing one- versus two-metal-mediated P-O bond cleavage reactions within the confines of the same active site, our work brings a fresh perspective to the "number of metals" controversy.

Nanoproteases: Alternatives to Natural Protease for Biotechnological Applications

Some nanomaterials with intrinsic protease-like activity have the advantages of good stability, biosafety, low price, large-scale preparation and unique property of nanomaterials, which are promising alternatives for natural proteases in various applications. An especial term, "nanoprotease", has been coined to stress the intrinsic proteolytic property of these nanomaterials. As a new generation of artificial proteases, they have become a burgeoning field, attracting many researchers to design and synthesize high performance nanoproteases. In this review, we summarize recent progress on all types of nanoproteases with regard of their activity, mechanism and application and introduce a new and effective strategy for engineering high-performance nanoproteases. In addition, we discuss the challenges and opportunities of nanoprotease research in the future.

Implications of albumin in cell culture media on the biological action of vanadates(V)

In this article, the implications of binding competition of vanadates(V) with dodecyl sulfates for bovine serum albumin on cytotoxicity of vanadium(V) species against prostate cancer cells have been investigated. The pH- and SDS-dependent vanadate(V)-BSA interactions were observed. At pH 5, there is only one site capable of binding ten vanadates(V) ions (logK(ITC)1 = 4.96 ± 0.06; ΔH(ITC)1 = -1.04 ± 0.03 kcal mol-1), whereas at pH 7 two distinctive binding sites on protein were found, saturated with two and seven V(V) ions, respectively (logK(ITC)1 = 6.11 ± 0.06; ΔH(ITC)1 = 0.78 ± 0.12 kcal mol-1; logK(ITC)2 = 4.80 ± 0.02; ΔH(ITC)2 = - 4.95 ± 0.14 kcal mol-1). SDS influences the stoichiometry and the stability of the resulting V(V)-BSA complexes. Finally, the cytotoxicity of vanadates(V) against prostate cancer cells (PC3 line) was examined in the presence and absence of SDS in the culture medium. In the case of a 24-h incubation with 100 μM vanadate(V), a ca. 20 % reduction in viability of PC3 cells was observed in the presence of SDS. However, in other considered cases (various concentrations and time of incubation) SDS does not affect the dose-dependent action of vanadates(V) on the investigated prostate cancer cells.

Distinct chemical factors in hydrolytic reactions catalyzed by metalloenzymes and metal complexes

The selective hydrolysis of the extremely stable phosphoester, peptide and ester bonds of molecules by bio-inspired metal-based catalysts (metallohydrolases) is required in a wide range of biological, biotechnological and industrial applications. Despite the impressive advances made in the field, the ultimate goal of designing efficient enzyme mimics for these reactions is still elusive. Its realization will require a deeper understanding of the diverse chemical factors that influence the activities of both natural and synthetic catalysts. They include catalyst-substrate complexation, non-covalent interactions and the electronic nature of the metal ion, ligand environment and nucleophile. Based on our computational studies, their roles are discussed for several mono- and binuclear metallohydrolases and their synthetic analogues. Hydrolysis by natural metallohydrolases is found to be promoted by a ligand environment with low basicity, a metal bound water and a heterobinuclear metal center (in binuclear enzymes). Additionally, peptide and phosphoester hydrolysis is dominated by two competing effects, i.e. nucleophilicity and Lewis acid activation, respectively. In synthetic analogues, hydrolysis is facilitated by the inclusion of a second metal center, hydrophobic effects, a biological metal (Zn, Cu and Co) and a terminal hydroxyl nucleophile. Due to the absence of the protein environment, hydrolysis by these small molecules is exclusively influenced by nucleophile activation. The results gleaned from these studies will enhance the understanding of fundamental principles of multiple hydrolytic reactions. They will also advance the development of computational methods as a predictive tool to design more efficient catalysts for hydrolysis, Diels-Alder reaction, Michael addition, epoxide opening and aldol condensation.

Computational Modelling of the Interactions Between Polyoxometalates and Biological Systems

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.

Peptide Hydrolysis by Metal (Oxa)cyclen Complexes: Revisiting the Mechanism and Assessing Ligand Effects

The mechanism responsible for peptide bond hydrolysis by Co(III) and Cu(II) complexes with (oxa)cyclen ligands has been revisited by means of computational tools. We propose that the mechanism starts by substrate coordination and an outer-sphere attack on the amide C atom of a solvent water molecule assisted by the metal hydroxo moiety as a general base, which occurs through six-membered ring transition states. This new mechanism represents a more likely scenario than the previously proposed mechanisms that involved an inner-sphere nucleophilic attack through more strained four-membered rings transition states. The corresponding computed overall free-energy barrier of 25.2 kcal mol^-1 for hydrolysis of the peptide bond in Phe-Ala by a cobalt(III) oxacyclen catalyst (1) is consistent with the experimental values obtained from rate constants. Also, we assessed the influence of the nature of the ligand throughout a systematic replacement of N by O atoms in the (oxa)cyclen ligand. Increasing the number of coordinating O atoms accelerates the reaction by increasing the Lewis acidity of the metal ion. On the other hand, the higher reactivity observed for the copper(II) oxacyclen catalyst with respect to the analogous Co(III) complex can be attributed to the larger Brönsted basicity of the copper(II) hydroxo ligand. Ultimately, the detailed understanding of the ligand and metal nature effects allowed us to identify the double role of the metal hydroxo complexes as Lewis acids and Brönsted bases and to rationalize the observed reactivity trends.

Mechanism of the highly effective peptide bond hydrolysis by MOF-808 catalyst under biologically relevant conditions

Efficient and selective hydrolysis of inert peptide bonds is of paramount importance. MOF-808, a metal-organic framework based on Zr₆ nodes, can hydrolyze peptide bonds efficiently under biologically relevant conditions. However, the details of the catalyst structure and of the underlying catalytic reaction mechanism are challenging to establish. By means of DFT calculations we first investigate the speciation of the Zr₆ nodes and identify the nature of ligands that bind to the Zr₆O₈H_4-x core in aqueous conditions. The core is predicted to strongly prefer a Zr₆O₈H₄ protonation state and to be predominantly decorated by bridging formate ligands, giving Zr₆(μ₃-O)₄(μ₃-OH)₄(BTC)₂(HCOO)₆ and Zr₆(μ₃-O)₄(μ₃-OH)₄(BTC)₂(HCOO)₅(OH)(H₂O) as the most favorable structures at physiological pH. The GlyGly peptide can bind MOF in several different ways, with the preferred structure involving coordination through the terminal carboxylate analogously to the binding mode of formate ligand. The pre-reactive binding mode in which the amide carbonyl oxygen coordinates the metal core lies 7 kcal higher in free energy. The preferred reaction pathway is predicted to have two close-lying transition states, either of which could be the rate-determining step: nucleophilic attack on the amide carbon atom and C-N bond breaking, with calculated relative free energies of 31 and 32 kcal mol^-1, respectively. Replacement of formate by water and hydroxide at the Zr₆ node is predicted to be possible, but does not appear to play a role in the hydrolysis mechanism.

Structure-Activity Relationships for the Affinity of Chaotropic Polyoxometalate Anions towards Proteins

The influence of the composition of chaotropic polyoxometalate (POM) anions on their affinity to biological systems was studied by means of atomistic molecular dynamics (MD) simulations. The variations in the affinity to hen egg-white lysozyme (HEWL) were analyzed along two series of POMs whereby the charge or the size and shape of the metal cluster are modified systematically. Our simulations revealed a quadratic relationship between the charge of the POM and its affinity to HEWL as a consequence of the parabolic growth of POM⋅⋅⋅water interaction with the charge. As the charge increases, POMs become less chaotropic (more kosmotropic) increasing the number and the strength of POM-water hydrogen bonds and structuring the solvation shell around the POM. This atomistic description explains the proportionally larger desolvation energies and less protein affinity for highly charged POMs, and consequently, the preference for moderate charge densities (q/M=0.33). Also, our simulations suggest that POM⋅⋅⋅protein interactions are size-specific. The cationic pockets of HEWL protein show a preference for Keggin-like structures, which display the optimal dimensions (≈1 nm). Finally, we developed a quantitative multidimensional model for protein affinity with predictive ability (r² =0.97; q² =0.88) 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; V_S ).

Effects of the Metal Ion on the Mechanism of Phosphodiester Hydrolysis Catalyzed by Metal-Cyclen Complexes

In this study, mechanisms of phosphodiester hydrolysis catalyzed by six di- and tetravalent metal-cyclen (M-C) complexes (Zn-C, Cu-C, Co-C, Ce-C, Zr-C and Ti-C) have been investigated using DFT calculations. The activities of these complexes were studied using three distinct mechanisms: (1) direct attack ( DA ), (2) catalyst-assisted ( CA ), and (3) water-assisted ( WA ). All divalent metal complexes (Zn-C, Cu-C and Co-C) coordinated to the BNPP substrate in a monodentate fashion and activated its scissile phosphoester bond. However, all tetravalent metal complexes (Ce-C, Zr-C, and Ti-C) interacted with BNPP in a bidentate manner and strengthened this bond. The DA mechanism was energetically the most feasible for all divalent M-C complexes, while the WA mechanism was favored by the tetravalent complexes, except Ce-C. The divalent complexes were found to be more reactive than their tetravalent counterparts. Zn-C catalyzed the hydrolysis with the lowest barrier among all M-C complexes, while Ti-C was the most reactive tetravalent complex. The activities of Ce-C and Zr-C, except Ti-C, were improved with an increase in the coordination number of the metal ion. The structural and mechanistic information provided in this study will be very helpful in the development of more efficient metal complexes for this critical reaction.