Challenges and advances in the computational modeling of biological phosphate hydrolysis

Hydrolytic activity of new bioinspired MnIIIMnII and FeIIIMnII complexes as mimetics of PAPs: Biological and environmental interest

Coordination compounds that mimic Purple Acid Phosphatases (PAPs) have drawn attention in the bioinorganic field due to their capacity to cleave phosphodiester bonds. However, their catalytic activity upon phosphate triesters is still unexplored. Thus, we report the synthesis and characterization of two binuclear complexes, [Mn^IIMn^III(L₁)(OAc)₂]BF₄ (1) and [Mn^IIFe^III(L₁)(OAc)₂]BF₄ (2) (H₂L₁ = 2-[N,N-bis-(2- pyridilmethyl)aminomethyl]-4-methyl-6-[N-(2-hydroxy-3-formyl-5-methylbenzyl)-N-(2-pyridylmethyl)aminomethyl]phenol), their hydrolytic activity and antioxidant potential. The complexes were fully characterized, including the X-Ray diffraction (XRD) of 1. Density functional theory (DFT) calculations were performed to better understand their electronic and structural properties and phosphate conjugates. The catalytic activity was analyzed for two model substrates, a diester (BDNPP) and a triester phosphate (DEDNPP). The results suggest enhancement of the hydrolysis reaction by 170 to 1500 times, depending on the substrate and complex. It was possible to accompany the catalytic reaction of DEDNPP hydrolysis by phosphorus nuclear magnetic resonance (³¹P NMR), showing that both 1 and 2 are efficient catalysts. Moreover, we also addressed that 1 and 2 present a relevant antioxidant potential, protecting the yeast Saccharomyces cerevisiae, used as eukaryotic model of study, against the exposure of cells to acute oxidative stress.

Water-Regulated Mechanisms for Degradation of Pesticides Paraoxon and Parathion by Phosphotriesterase: Insight from QM/MM and MD Simulations

The enzymatic degradation of pesticides paraoxon (PON) and parathion (PIN) by phosphotriesterase (PTE) has been investigated by QM/MM calculations and MD simulations. In the PTE-PON complex, Zn α and Zn β in the active site are five- and six-coordinated, respectively, while both zinc ions are six coordinated with the Zn α -bound water molecule (WT1) for the PTE-PIN system. The hydrolytic reactions for PON and PIN are respectively driven by the nucleophilic attack of the bridging-OH - and the Zn α -bound water molecule on the phosphorus center of substrate, and the two-step hydrolytic process is predicted to be the rate-limiting step with the energy spans of 13.8 and 14.4 kcal/mol for PON and PIN, respectively. The computational studies reveal that the presence of the Zn α -bound water molecule depends on the structural feature of substrates characterized by P=O and P=S, which determines the hydrolytic mechanism and efficiency for the degradation of organophosphorus pesticides by PTE.

Origin of iodine preferential attack at sulfur in phosphorothioate and subsequent P-O or P-S bond dissociation

Iodine-induced cleavage at phosphorothioate DNA (PT-DNA) is characterized by extremely high sensitivity (∼1 phosphorothioate link per 106 nucleotides), which has been used for detecting and sequencing PT-DNA in bacteria. Despite its foreseeable potential for wide applications, the cleavage mechanism at the PT-modified site has not been well established, and it remains unknown as to whether or not cleavage of the bridging P-O occurs at every PT-modified site. In this work, we conducted accurate ωB97X-D calculations and high-performance liquid chromatography-mass spectrometry to investigate the process of PT-DNA cleavage at the atomic and molecular levels. We have found that iodine chemoselectively binds to the sulfur atom of the phosphorothioate link via a strong halogen-chalcogen interaction (a type of halogen bond, with binding affinity as high as 14.9 kcal/mol) and thus triggers P-O bond cleavage via phosphotriester-like hydrolysis. Additionally, aside from cleavage of the bridging P-O bond, the downstream hydrolyses lead to unwanted P-S/P-O conversions and a loss of the phosphorothioate handle. The mechanism we outline helps to explain specific selectivity at the PT-modified site but also predicts the dynamic stoichiometry of P-S and P-O bond breaking. For instance, Tris is involved in the cascade derivation of S-iodo-phosphorothioate to S-amino-phosphorothioate, suppressing the S-iodo-phosphorothioate hydrolysis to a phosphate diester. However, hydrolysis of one-third of the Tris-O-grafting phosphotriester results in unwanted P-S/P-O conversions. Our study suggests that bacterial DNA phosphorothioation may more frequently occur than previous bioinformatic estimations have predicted from iodine-induced deep sequencing data.

Degradation of pesticides diazinon and diazoxon by phosphotriesterase: insight into divergent mechanisms from QM/MM and MD simulations

Enzymatic hydrolysis by phosphotriesterase (PTE) is one of the most effective ways of degrading organophosphorus pesticides, but the catalytic efficiency depends on the structural features of substrates. Here the enzymatic degradation of diazinon (DIN) and diazoxon (DON), characterized by PS and PO, respectively, have been investigated by QM/MM calculations and MM MD simulations. Our calculations demonstrate that the hydrolysis of DON (with PO) is inevitably initiated by the nucleophilic attack of the bridging-OH^- on the phosphorus center, while for DIN (with PS), we proposed a new degradation mechanism, initiated by the nucleophilic attack of the Zn_α-bound water molecule, for its low-energy pathway. For both DIN and DON, the hydrolytic reaction is predicted to be the rate-limiting step, with energy barriers of 18.5 and 17.7 kcal mol^-1, respectively. The transportation of substrates to the active site, the release of the leaving group and the degraded product are generally verified to be favorable by MD simulations via umbrella sampling, both thermodynamically and dynamically. The side-chain residues Phe132, Leu271 and Tyr309 play the gate-switching role to manipulate substrate delivery and product release. In comparison with the DON-enzyme system, the degraded product of DIN is more easily released from the active site. These new findings will contribute to the comprehensive understanding of the enzymatic degradation of toxic organophosphorus compounds by PTE.

New Insights on the Interaction of Phenanthroline Based Ligands and Metal Complexes and Polyoxometalates with Duplex DNA and G-Quadruplexes

This work provides new insights from our team regarding advances in targeting canonical and non-canonical nucleic acid structures. This modality of medical treatment is used as a form of molecular medicine specifically against the growth of cancer cells. Nevertheless, because of increasing concerns about bacterial antibiotic resistance, this medical strategy is also being explored in this field. Up to three strategies for the use of DNA as target have been studied in our research lines during the last few years: (1) the intercalation of phenanthroline derivatives with duplex DNA; (2) the interaction of metal complexes containing phenanthroline with G-quadruplexes; and (3) the activity of Mo polyoxometalates and other Mo-oxo species as artificial phosphoesterases to catalyze the hydrolysis of phosphoester bonds in DNA. We demonstrate some promising computational results concerning the favorable interaction of these small molecules with DNA that could correspond to cytotoxic effects against tumoral cells and microorganisms. Therefore, our results open the door for the pharmaceutical and medical applications of the compounds we propose.

Mechanistic Insights into Promoted Hydrolysis of Phosphoester Bonds by MoO2Cl2(DMF)2

A phosphoester bond is a crucial structural block in biological systems, whose occurrence is regulated by phosphatases. Molybdenum compounds have been reported to be active in phosphate ester hydrolysis of model phosphates. Specifically, MoO₂Cl₂(DMF)₂ is active in the hydrolysis of para-nitrophenyl phosphate (pNPP), leading to heteropolyoxometalate structures. We use density functional theory (DFT) to clarify the mechanism by which these species promote the hydrolysis of the phosphoester bond. The present calculations give insight into several key aspects of this reaction: (i) the speciation of this complex prior to interaction with the phosphate (DMF release, Mo-Cl hydrolysis, and pH influence on the speciation), (ii) the competition between phosphate addition and the molybdate nucleation process, (iii) and the mechanisms by which some plausible active species promote this hydrolysis in different conditions. We described thoroughly two different pathways depending on the nucleation possibilities of the molybdenum complex: one mononuclear mechanism, which is preferred in conditions in which very low complex concentrations are used, and another dinuclear mechanism, which is preferred at higher concentrations.

Mechanisms of a Cyclobutane-Fused Lactone Hydrolysis in Alkaline and Acidic Conditions

Searching for functional polyesters with stability and degradability is important due to their potential applications in biomedical supplies, biomass fuel, and environmental protection. Recently, a cyclobutane-fused lactone (CBL) polymer was experimentally found to have superior stability and controllable degradability through hydrolysis reactions after activation by mechanical force. In order to provide a theoretical basis for developing new functional degradable polyesters, in this work, we performed a detailed quantum chemical study of the alkaline and acidic hydrolysis of CBL using dispersion-corrected density functional theory (DFT-D3) and mixed implicit/explicit solvent models. Various possible hydrolysis mechanisms were found: B_AC2 and B_AL2 in the alkaline condition and A_AC2, A_AL2, and A_AL1 in the acidic condition. Our calculations indicated that CBL favors the B_AC2 and A_AC2 mechanisms in alkaline and acidic conditions, respectively. In addition, we found that incorporating explicit water solvent molecules is highly necessary because of their strong hydrogen-bonding with reactant/intermediate/product molecules.

The evolution of nucleosidic analogues: self-assembly of prodrugs into nanoparticles for cancer drug delivery

Nucleoside and nucleotide analogs are essential tools in our limited arsenal in the fight against cancer. However, these structures face severe drawbacks such as rapid plasma degradation or hydrophilicity, limiting their clinical application. Here, different aspects of nucleoside and nucleotide analogs have been exposed, while providing their shortcomings. Aiming to improve their fate in the body and combating their drawbacks, two different approaches have been discussed, the prodrug and nanocarrier technologies. Finally, a novel approach called "PUFAylation" based on both the prodrug and nanocarrier technologies has been introduced, promising to be the supreme method to create a novel nucleoside or nucleotide analog based formulation, with enhanced efficacy and highly reduced toxicity.

Examining the Mechanism of Phosphite Dehydrogenase with Quantum Mechanical/Molecular Mechanical Free Energy Simulations

The projected decline of available phosphorus necessitates alternative methods to derive usable phosphate for fertilizer and other applications. Phosphite dehydrogenase oxidizes phosphite to phosphate with the cofactor NAD⁺ serving as the hydride acceptor. In addition to producing phosphate, this enzyme plays an important role in NADH cofactor regeneration processes. Mixed quantum mechanical/molecular mechanical free energy simulations were performed to elucidate the mechanism of this enzyme and to identify the protonation states of the substrate and product. Specifically, the finite temperature string method with umbrella sampling was used to generate the free energy surfaces and determine the minimum free energy paths for six different initial conditions that varied in the protonation state of the substrate and the position of the nucleophilic water molecule. In contrast to previous studies, the mechanism predicted by all six independent strings is a concerted but asynchronous dissociative mechanism in which hydride transfer from the phosphite substrate to NAD⁺ occurs prior to attack by the nucleophilic water molecule. His292 is identified as the most likely general base that deprotonates the attacking water molecule. However, Arg237 could also serve as this base if it were deprotonated and His292 were protonated prior to the main chemical transformation, although this scenario is less probable. The simulations indicate that the phosphite substrate is monoanionic in its active form and that the most likely product is dihydrogen phosphate. These mechanistic insights may be helpful for designing mutant enzymes or artificial constructs that convert phosphite to phosphate and NAD⁺ to NADH more effectively.

GTP Hydrolysis Without an Active Site Base: A Unifying Mechanism for Ras and Related GTPases

GTP hydrolysis is a biologically crucial reaction, being involved in regulating almost all cellular processes. As a result, the enzymes that catalyze this reaction are among the most important drug targets. Despite their vital importance and decades of substantial research effort, the fundamental mechanism of enzyme-catalyzed GTP hydrolysis by GTPases remains highly controversial. Specifically, how do these regulatory proteins hydrolyze GTP without an obvious general base in the active site to activate the water molecule for nucleophilic attack? To answer this question, we perform empirical valence bond simulations of GTPase-catalyzed GTP hydrolysis, comparing solvent- and substrate-assisted pathways in three distinct GTPases, Ras, Rab, and the G_αi subunit of a heterotrimeric G-protein, both in the presence and in the absence of the corresponding GTPase activating proteins. Our results demonstrate that a general base is not needed in the active site, as the preferred mechanism for GTP hydrolysis is a conserved solvent-assisted pathway. This pathway involves the rate-limiting nucleophilic attack of a water molecule, leading to a short-lived intermediate that tautomerizes to form H₂PO₄^- and GDP as the final products. Our fundamental biochemical insight into the enzymatic regulation of GTP hydrolysis not only resolves a decades-old mechanistic controversy but also has high relevance for drug discovery efforts. That is, revisiting the role of oncogenic mutants with respect to our mechanistic findings would pave the way for a new starting point to discover drugs for (so far) "undruggable" GTPases like Ras.

Assessing the Performance of the Nonbonded Mg2+ Models in a Two-Metal-Dependent Ribonuclease

Magnesium ions (Mg²⁺), abundant in living cells, are essential for biomolecular structure, dynamics, and function. The biological importance of Mg²⁺ has motivated continuous development and improvement of various Mg²⁺ models for molecular dynamics (MD) simulations during the last decades. There are four types of nonbonded Mg²⁺ models: the point charge models based on a 12-6 or 12-6-4 type Lennard-Jones (LJ) potential, and the multisite models based on a 12-6 or 12-6-4 LJ potential. Here, we systematically assessed the performance of these four types of nonbonded Mg²⁺ models (21 models in total) in terms of maintaining a challenging intermediate state configuration captured in the structure of a prototypical two-metal-ion RNase H complex with an RNA/DNA hybrid. Our data demonstrate that the 12-6-4 multisite models, which account for charge-induced dipole interactions, perform the best in reproducing all the unique coordination modes in this intermediate state and maintaining the correct carboxylate denticity. Our benchmark work provides a useful guideline for MD simulations and structural refinement of Mg²⁺-containing biomolecular systems.

Conformational dynamics and enzyme evolution

Enzymes are dynamic entities, and their dynamic properties are clearly linked to their biological function. It follows that dynamics ought to play an essential role in enzyme evolution. Indeed, a link between conformational diversity and the emergence of new enzyme functionalities has been recognized for many years. However, it is only recently that state-of-the-art computational and experimental approaches are revealing the crucial molecular details of this link. Specifically, evolutionary trajectories leading to functional optimization for a given host environment or to the emergence of a new function typically involve enriching catalytically competent conformations and/or the freezing out of non-competent conformations of an enzyme. In some cases, these evolutionary changes are achieved through distant mutations that shift the protein ensemble towards productive conformations. Multifunctional intermediates in evolutionary trajectories are probably multi-conformational, i.e. able to switch between different overall conformations, each competent for a given function. Conformational diversity can assist the emergence of a completely new active site through a single mutation by facilitating transition-state binding. We propose that this mechanism may have played a role in the emergence of enzymes at the primordial, progenote stage, where it was plausibly promoted by high environmental temperatures and the possibility of additional phenotypic mutations.