What roles do the residue Asp229 and the coordination variation of calcium play of the reaction mechanism of the diisopropyl-fluorophosphatase? A DFT investigation

Eliminating Imaginary Vibrational Frequencies in Quantum-Chemical Cluster Models of Enzymatic Active Sites

In constructing finite models of enzyme active sites for quantum-chemical calculations, atoms at the periphery of the model must be constrained to prevent unphysical rearrangements during geometry relaxation. A simple fixed-atom or "coordinate-lock" approach is commonly employed but leads to undesirable artifacts in the form of small imaginary frequencies. These preclude evaluation of finite-temperature free-energy corrections, limiting thermochemical calculations to enthalpies only. Full-dimensional vibrational frequency calculations are possible by replacing the fixed-atom constraints with harmonic confining potentials. Here, we compare that approach to an alternative strategy in which fixed-atom contributions to the Hessian are simply omitted. While the latter strategy does eliminate imaginary frequencies, it tends to underestimate both the zero-point energy and the vibrational entropy while introducing artificial rigidity. Harmonic confining potentials eliminate imaginary frequencies and provide a flexible means to construct active-site models that can be used in unconstrained geometry relaxations, affording better convergence of reaction energies and barrier heights with respect to the model size, as compared to models with fixed-atom constraints.

The His23 and Lys79 pair determines the high catalytic efficiency of the inorganic pyrophosphatase of the haloacid dehalogenase superfamily

Haloacid dehalogenase (HAD) superfamily members are mainly phosphomonoesterases, while BT2127 from Bacteroides thetaiotaomicron of the HAD superfamily is identified as an inorganic pyrophosphatase. In this study, to explore the roles of the Lys79 and His23 pair in the hydrolysis reaction of inorganic pyrophosphate (PP_i) catalyzed by BT2127, a series of models were designed. Calculations were performed by using the density functional theory (DFT) method with the dispersion energy D3-B3LYP. The His23 and Lys79 pair plays a key role in the high catalytic efficiency of BT2127 with PPi. First, the His23 and Lys79 pair prompts Asp13 to easily provide a proton to the leaving group, which remarkably reduces the energy barrier of the phospho-transfer step; then, Lys79 provides a proton to the first leaving phosphate group via His23, produces a more electrically stabilized phosphate (H₃PO₄), makes this step exothermal, and further promotes the subsequent phospho-enzyme intermediate hydrolysis. The results suggest that the Lys79-His23 pair helps BT2127 reach high catalytic efficiency by strengthening the acid catalysis. Our study provides detailed chemical insights into the evolution of the inorganic pyrophosphatase function of BT2127 from the phosphomonoesterase of the HAD superfamily and the biomimetic enzyme design.

Probing the Suitability of Different Ca2+ Parameters for Long Simulations of Diisopropyl Fluorophosphatase

Organophosphate hydrolases are promising as potential biotherapeutic agents to treat poisoning with pesticides or nerve gases. However, these enzymes often need to be further engineered in order to become useful in practice. One example of such enhancement is the alteration of enantioselectivity of diisopropyl fluorophosphatase (DFPase). Molecular modeling techniques offer a unique opportunity to address this task rationally by providing a physical description of the substrate-binding process. However, DFPase is a metalloenzyme, and correct modeling of metal cations is a challenging task generally coming with a tradeoff between simulation speed and accuracy. Here, we probe several molecular mechanical parameter combinations for their ability to empower long simulations needed to achieve a quantitative description of substrate binding. We demonstrate that a combination of the Amber19sb force field with the recently developed 12-6 Ca2+ models allows us to both correctly model DFPase and obtain new insights into the DFP binding process.

Enzymes, Reacting with Organophosphorus Compounds as Detoxifiers: Diversity and Functions

Organophosphorus compounds (OPCs) are able to interact with various biological targets in living organisms, including enzymes. The binding of OPCs to enzymes does not always lead to negative consequences for the body itself, since there are a lot of natural biocatalysts that can catalyze the chemical transformations of the OPCs via hydrolysis or oxidation/reduction and thereby provide their detoxification. Some of these enzymes, their structural differences and identity, mechanisms, and specificity of catalytic action are discussed in this work, including results of computational modeling. Phylogenetic analysis of these diverse enzymes was specially realized for this review to emphasize a great area for future development(s) and applications.

Enzymatic Degradation of Organophosphorus Pesticides and Nerve Agents by EC: 3.1.8.2

The organophosphorus substances, including pesticides and nerve agents (NAs), represent highly toxic compounds. Standard decontamination procedures place a heavy burden on the environment. Given their continued utilization or existence, considerable efforts are being made to develop environmentally friendly methods of decontamination and medical countermeasures against their intoxication. Enzymes can offer both environmental and medical applications. One of the most promising enzymes cleaving organophosphorus compounds is the enzyme with enzyme commission number (EC): 3.1.8.2, called diisopropyl fluorophosphatase (DFPase) or organophosphorus acid anhydrolase from Loligo Vulgaris or Alteromonas sp. JD6.5, respectively. Structure, mechanisms of action and substrate profiles are described for both enzymes. Wild-type (WT) enzymes have a catalytic activity against organophosphorus compounds, including G-type nerve agents. Their stereochemical preference aims their activity towards less toxic enantiomers of the chiral phosphorus center found in most chemical warfare agents. Site-direct mutagenesis has systematically improved the active site of the enzyme. These efforts have resulted in the improvement of catalytic activity and have led to the identification of variants that are more effective at detoxifying both G-type and V-type nerve agents. Some of these variants have become part of commercially available decontamination mixtures.

Calcium Valence-to-Core X-ray Emission Spectroscopy: A Sensitive Probe of Oxo Protonation in Structural Models of the Oxygen-Evolving Complex

Calcium is an abundant, nontoxic metal that finds many roles in synthetic and biological systems including the oxygen-evolving complex (OEC) of photosystem II. Characterization methods for calcium centers, however, are underdeveloped compared to those available for transition metals. Valence-to-core X-ray emission spectroscopy (VtC XES) selectively probes the electronic structure of an element's chemical environment, providing insight that complements the geometric information available from other techniques. Here, the utility of calcium VtC XES is established using an in-house dispersive spectrometer in combination with density functional theory. Spectral trends are rationalized within a molecular orbital framework, and Kβ2,5 transitions, derived from molecular orbitals with primarily ligand p character, are found to be a promising probe of the calcium coordination environment. In particular, it is shown that calcium VtC XES is sensitive to the electronic structure changes that accompany oxo protonation in Mn3CaO4-based molecular mimics of the OEC. Through correlation to calculations, the potential of calcium VtC XES to address unresolved questions regarding the mechanism of biological water oxidation is highlighted.

Theoretical Studies on Catalysis Mechanisms of Serum Paraoxonase 1 and Phosphotriesterase Diisopropyl Fluorophosphatase Suggest the Alteration of Substrate Preference from Paraoxonase to DFP

The calcium-dependent β-propeller proteins mammalian serum paraoxonase 1 (PON1) and phosphotriesterase diisopropyl fluorophosphatase (DFPase) catalyze the hydrolysis of organophosphorus compounds and enhance hydrolysis of various nerve agents. In the present work, the phosphotriesterase activity development between PON1 and DFPase was investigated by using the hybrid density functional theory method B3LYP. Based on the active-site difference between PON1 and DFPase, both the wild type and the mutant (a water molecule replacing Asn270 in PON1) models were designed. The results indicated that the substitution of a water molecule for Asn270 in PON1 had little effect on the enzyme activity in kinetics, while being more efficient in thermodynamics, which is essential for DFP hydrolysis. Structure comparisons of evolutionarily related enzymes show that the mutation of Asn270 leads to the catalytic Ca²⁺ ion indirectly connecting the buried structural Ca²⁺ ion via hydrogen bonds in DFPase. It can reduce the plasticity of enzymatic structure, and possibly change the substrate preference from paraoxon to DFP, which implies an evolutionary transition from mono- to dinuclear catalytic centers. Our studies shed light on the investigation of enzyme catalysis mechanism from an evolutionary perspective.

Theoretical Studies Applied to the Evaluation of the DFPase Bioremediation Potential against Chemical Warfare Agents Intoxication

Organophosphorus compounds (OP) are part of a group of compounds that may be hazardous to health. They are called neurotoxic agents because of their action on the nervous system, inhibiting the acetylcholinesterase (AChE) enzyme and resulting in a cholinergic crisis. Their high toxicity and rapid action lead to irreversible damage to the nervous system, drawing attention to developing new treatment methods. The diisopropyl fluorophosphatase (DFPase) enzyme has been considered as a potent biocatalyst for the hydrolysis of toxic OP and has potential for bioremediation of this kind of intoxication. In order to investigate the degradation process of the nerve agents Tabun, Cyclosarin and Soman through the wild-type DFPase, and taking into account their stereochemistry, theoretical studies were carried out. The intermolecular interaction energy and other parameters obtained from the molecular docking calculations were used to construct a data matrix, which were posteriorly treated by statistical analyzes of chemometrics, using the PCA (Principal Components Analysis) multivariate analysis. The analyzed parameters seem to be quite important for the reaction mechanisms simulation (QM/MM). Our findings showed that the wild-type DFPase enzyme is stereoselective in hydrolysis, showing promising results for the catalytic degradation of the neurotoxic agents under study, with the degradation mechanism performed through two proposed pathways.

Similar Active Sites and Mechanisms Do Not Lead to Cross-Promiscuity in Organophosphate Hydrolysis: Implications for Biotherapeutic Engineering

Organophosphate hydrolases are proficient catalysts of the breakdown of neurotoxic organophosphates and have great potential as both biotherapeutics for treating acute organophosphate toxicity and as bioremediation agents. However, proficient organophosphatases such as serum paraoxonase 1 (PON1) and the organophosphate-hydrolyzing lactonase SsoPox are unable to hydrolyze bulkyorganophosphates with challenging leaving groups such as diisopropyl fluorophosphate (DFP) or venomous agent X, creating a major challenge for enzyme design. Curiously, despite their mutually exclusive substrate specificities, PON1 and diisopropyl fluorophosphatase (DFPase) have essentially identical active sites and tertiary structures. In the present work, we use empirical valence bond simulations to probe the catalytic mechanism of DFPase as well as temperature, pH, and mutational effects, demonstrating that DFPase and PON1 also likely utilize identical catalytic mechanisms to hydrolyze their respective substrates. However, detailed examination of both static structures and dynamical simulations demonstrates subtle but significant differences in the electrostatic properties and solvent penetration of the two active sites and, most critically, the role of residues that make no direct contact with either substrate in acting as "specificity switches" between the two enzymes. Specifically, we demonstrate that key residues that are structurally and functionally critical for the paraoxonase activity of PON1 prevent it from being able to hydrolyze DFP with its fluoride leaving group. These insights expand our understanding of the drivers of the evolution of divergent substrate specificity in enzymes with identical active sites and guide the future design of organophosphate hydrolases that hydrolyze compounds with challenging leaving groups.