Theoretical study of the phosphotriesterase reaction mechanism

Catalytic mechanisms for phosphotriesterases

Phosphotriesters are one class of highly toxic synthetic compounds known as organophosphates. Wide spread usage of organophosphates as insecticides as well as nerve agents has lead to numerous efforts to identify enzymes capable of detoxifying them. A wide array of enzymes has been found to have phosphotriesterase activity including phosphotriesterase (PTE), methyl parathion hydrolase (MPH), organophosphorus acid anhydrolase (OPAA), diisopropylfluorophosphatase (DFP), and paraoxonase 1 (PON1). These enzymes differ widely in protein sequence and three-dimensional structure, as well as in catalytic mechanism, but they also share several common features. All of the enzymes identified as phosphotriesterases are metal-dependent hydrolases that contain a hydrophobic active site with three discrete binding pockets to accommodate the substrate ester groups. Activation of the substrate phosphorus center is achieved by a direct interaction between the phosphoryl oxygen and a divalent metal in the active site. The mechanistic details of the hydrolytic reaction differ among the various enzymes with both direct attack of a hydroxide as well as covalent catalysis being found. This article is part of a Special Issue entitled: Chemistry and mechanism of phosphatases, diesterases and triesterases.

Repair of DNA Dewar photoproduct to (6-4) photoproduct in (6-4) photolyase

Dewar photoproduct (Dewar PP) is the valence isomer of (6-4) photoproduct ((6-4)PP) in photodamaged DNA. Compared to the extensive studied CPD photoproducts, the underlying repair mechanisms for the (6-4)PP, and especially for the Dewar PP, are not well-established to date. In this paper, the repair mechanism of DNA Dewar photoproduct T(dew)C in (6-4) photolyase was elucidated using hybrid density functional theory. Our results showed that, during the repair process, the T(dew)C has to isomerize to T(6-4)C photolesion first via direct C6'-N3' bond cleavage facilitated by electron injection. This isomerization mechanism is energetically much more efficient than other possible rearrangement pathways. The calculations provide a theoretical interpretation to recent experimental observations.

The organophosphate-degrading enzyme from Agrobacterium radiobacter displays mechanistic flexibility for catalysis

The OP (organophosphate)-degrading enzyme from Agrobacterium radiobacter (OpdA) is a binuclear metallohydrolase able to degrade highly toxic OP pesticides and nerve agents into less or non-toxic compounds. In the present study, the effect of metal ion substitutions and site-directed mutations on the catalytic properties of OpdA are investigated. The study shows the importance of both the metal ion composition and a hydrogen-bond network that connects the metal ion centre with the substrate-binding pocket using residues Arg254 and Tyr257 in the mechanism and substrate specificity of this enzyme. For the Co(II) derivative of OpdA two protonation equilibria (pKa1 ~5; pKa2 ~10) have been identified as relevant for catalysis, and a terminal hydroxide acts as the likely hydrolysis-initiating nucleophile. In contrast, the Zn(II) and Cd(II) derivatives only have one relevant protonation equilibrium (pKa ~4-5), and the μOH is the proposed nucleophile. The observed mechanistic flexibility may reconcile contrasting reaction models that have been published previously and may be beneficial for the rapid adaptation of OP-degrading enzymes to changing environmental pressures.

Protonation of a hydroxide anion bridging two divalent magnesium cations in water probed by first-principles metadynamics simulation

The protonation of a hydroxide anion (OH(-)) located between two magnesium cations (Mg(2+)) in aqueous solution has been investigated by first-principles metadynamics simulation. We observe that the complex Mg(2+)-OH(-)-Mg(2+) is stabilized by the coparticipation of the hydroxide anion to the first hydration shells of both the Mg(2+) cations. Contrary to the cases of OH(-) in pure water, the transfer of protons in the presence of the divalent metal ions turns out to be a slow chemical event. This can be ascribed to the decreased proton affinity of the bridging OH(-). Metadynamics simulation, used to overcome the difficulty of the long time scale required by the protonation of the bridging OH(-), has shown that the system possesses a great stability on the reactant state, characterized by a bioctahedral (6,6) solvation structure around the two Mg(2+) cations. The exploration of the free energy landscape shows that this stable bioctahedral configuration converts into a lower coordinated (5,6) structure, leading to a proton transfer from a water molecule belonging to the first solvation shell of the Mg(2+) ion having the lower coordination to the bridging OH(-); the free energy barrier for the protonation reaction is 11 kcal/mol, meaning that the bridging hydroxide is a weak base. During the proton transfer, the bridging OH(-) reverts to an H(2)O molecule, and this breaks the electrostatic coupling of the two Mg(2+) ions, which depart independently with their own hydration shells, one of which is entirely formed by water molecules. The second one carries the newly created OH(-). Our results show that the flexibility in the metal coordination plays a crucial role in both the protonation process of the bridging OH(-) and the separation of the metal cations, providing useful insight into the nature of proton transfer in binuclear divalent metal ions, with several biological implications, such as, for instance, transesterification of catalytic RNA.

Improving a natural enzyme activity through incorporation of unnatural amino acids

The bacterial phosphotriesterases catalyze hydrolysis of the pesticide paraoxon with very fast turnover rates and are thought to be near to their evolutionary limit for this activity. To test whether the naturally evolved turnover rate could be improved through the incorporation of unnatural amino acids and to probe the role of peripheral active site residues in nonchemical steps of the catalytic cycle (substrate binding and product release), we replaced the naturally occurring tyrosine amino acid at position 309 with unnatural L-(7-hydroxycoumarin-4-yl)ethylglycine (Hco) and L-(7-methylcoumarin-4-yl)ethylglycine amino acids, as well as leucine, phenylalanine, and tryptophan. Kinetic analysis suggests that the 7-hydroxyl group of Hco, particularly in its deprotonated state, contributes to an increase in the rate-limiting product release step of substrate turnover as a result of its electrostatic repulsion of the negatively charged 4-nitrophenolate product of paraoxon hydrolysis. The 8-11-fold improvement of this already highly efficient catalyst through a single rationally designed mutation using an unnatural amino acid stands in contrast to the difficulty in improving this native activity through screening hundreds of thousands of mutants with natural amino acids. These results demonstrate that designer amino acids provide easy access to new and valuable sequence and functional space for the engineering and evolution of existing enzyme functions.

The mechanism of cyclization in chromophore maturation of green fluorescent protein: a theoretical study

An intriguing aspect of the green fluorescent protein (GFP) is the autocatalytic post-translational modification that results in the formation of its chromophore. Numerous experimental and theoretical studies indicate that cyclization is the first and the most important step in the maturation process. In this work, two proposed mechanisms for the cyclization were investigated by using the hybrid density functional theory method B3LYP. Cluster models corresponding to the two mechanisms proposed by Wachter et al. [J. Biol. Chem. 2005, 280, 26248-26255] are constructed on the basis of the X-ray crystal structure (PDB entry 2AWJ) and corresponding reaction path potential energy profiles for the two cyclization mechanisms are presented. Our results suggest that the backbone condensation initiated by deprotonation of the Gly67 amide nitrogen is easier than deprotonation of the Tyr66 alpha-carbon. Moreover, Arg96 fulfills the role of stabilizing the enolate moiety, and Glu222 plays the role of a general base. The formation of the cyclized product is found to be 16.0 and 18.6 kcal/mol endothermic with respect to the two models, which is in agreement with experimental observation.

Phosphate mono- and diesterase activities of the trinuclear zinc enzyme nuclease P1--insights from quantum chemical calculations

Nuclease P1 is a trinuclear zinc enzyme that catalyzes the hydrolysis of single-stranded DNA and RNA. Density functional calculations are used to elucidate the reaction mechanism of this enzyme with a model of the active site designed on the basis of the X-ray crystal structure. 2-Tetrahydrofuranyl phosphate and methyl 2-tetrahydrofuranyl phosphate substrates are used to explore the phosphomonoesterase and phosphodiesterase activities of this enzyme, respectively. The calculations reveal that for both activities, a bridging hydroxide performs an in-line attack on the phosphorus center, resulting in inversion of the configuration. Simultaneously, the P-O bond is cleaved, and Zn2 stabilizes the negative charge of the leaving alkoxide anion and assists its departure. All three zinc ions, together with Arg48, provide electrostatic stabilization to the penta-coordinated transition state, thereby lowering the reaction barrier.

Quantum chemical modeling of enzymatic reactions: the case of histone lysine methyltransferase

Quantum chemical cluster models of enzyme active sites are today an important and powerful tool in the study of various aspects of enzymatic reactivity. This methodology has been applied to a wide spectrum of reactions and many important mechanistic problems have been solved. Herein, we report a systematic study of the reaction mechanism of the histone lysine methyltransferase (HKMT) SET7/9 enzyme, which catalyzes the methylation of the N-terminal histone tail of the chromatin structure. In this study, HKMT SET7/9 serves as a representative case to examine the modeling approach for the important class of methyl transfer enzymes. Active site models of different sizes are used to evaluate the methodology. In particular, the dependence of the calculated energies on the model size, the influence of the dielectric medium, and the particular choice of the dielectric constant are discussed. In addition, we examine the validity of some technical aspects, such as geometry optimization in solvent or with a large basis set, and the use of different density functional methods.

Reaction mechanism of cis-3-chloroacrylic acid dehalogenase: a theoretical study

The reaction mechanism of cis-3-chloroacrylic acid dehalogenase (cis-CaaD) is studied using the B3LYP density functional theory method. This enzyme catalyzes the hydrolytic dehalogenation of cis-3-chloroacrylic acid to yield malonate semialdehyde and HCl. The uncatalyzed reaction is first considered, and excellent agreement is found between the calculated barrier and the measured rate constant. The enzymatic reaction is then studied with an active site model consisting of 159 atoms. The results suggest an alternative mechanism for cis-CaaD catalysis and different roles for some active site residues in this mechanism.

Molecular dynamics simulations of the detoxification of paraoxon catalyzed by phosphotriesterase

Combined QM(PM3)/MM molecular dynamics simulations together with QM(DFT)/MM optimizations for key configurations have been performed to elucidate the enzymatic catalysis mechanism on the detoxification of paraoxon by phosphotriesterase (PTE). In the simulations, the PM3 parameters for the phosphorous atom were reoptimized. The equilibrated configuration of the enzyme/substrate complex showed that paraoxon can strongly bind to the more solvent-exposed metal ion Zn(beta), but the free energy profile along the binding path demonstrated that the binding is thermodynamically unfavorable. This explains why the crystal structures of PTE with substrate analogues often exhibit long distances between the phosphoral oxygen and Zn(beta). The subsequent SN2 reaction plays the key role in the whole process, but controversies exist over the identity of the nucleophilic species, which could be either a hydroxide ion terminally coordinated to Zn(alpha) or the micro-hydroxo bridge between the alpha- and beta-metals. Our simulations supported the latter and showed that the rate-limiting step is the distortion of the bound paraoxon to approach the bridging hydroxide. After this preparation step, the bridging hydroxide ion attacks the phosphorous center and replaces the diethyl phosphate with a low barrier. Thus, a plausible way to engineer PTE with enhanced catalytic activity is to stabilize the deformed paraoxon. Conformational analyses indicate that Trp131 is the closest residue to the phosphoryl oxygen, and mutations to Arg or Gln or even Lys, which can shorten the hydrogen bond distance with the phosphoryl oxygen, could potentially lead to a mutant with enhanced activity for the detoxification of organophosphates.

DFT study on the mechanism of Escherichia coli inorganic pyrophosphatase

Escherichia coli inorganic pyrophosphatase (E-PPase) is a tetranuclear divalent metal dependent enzyme that catalyzes the reversible interconversion of pyrophosphate (PPi) and orthophosphate (Pi), with Mg(2+) conferring the highest activity. In the present work, the reaction mechanism of E-PPase is investigated using the hybrid density functional theory (DFT) method B3LYP with a large model of the active site. Our calculated results shed further light on the detailed reaction mechanism. In particular, the important residue Asp67, either protonated or unprotonated, was taken into account in the present calculations. Our calculations indicated that a protonated Asp67 is crucial for the reverse reaction to take place; however, it is lost sight of in the forward reaction. The bridging hydroxide is shown to be capable of performing nucleophilic in-line attack on the substrate from its bridging position in the presence of four Mg(2+) ions. During the catalysis, the roles of the four magnesium ions are suggested to provide a necessary conformation of the active site, facilitate the nucleophile formation and substrate orientation, and stabilize the trigonal bipyramid transition state, thereby lowering the barrier for the nucleophilic attack.

Structure of diethyl phosphate bound to the binuclear metal center of phosphotriesterase

The bacterial phosphotriesterase (PTE) from Pseudomonas diminuta catalyzes the hydrolysis of organophosphate esters at rates close to the diffusion limit. X-ray diffraction studies have shown that a binuclear metal center is positioned in the active site of PTE and that this complex is responsible for the activation of the nucleophilic water from solvent. In this paper, the three-dimensional structure of PTE was determined in the presence of the hydrolysis product, diethyl phosphate (DEP), and a product analogue, cacodylate. In the structure of the PTE-diethyl phosphate complex, the DEP product is found symmetrically bridging the two divalent cations. The DEP displaces the hydroxide from solvent that normally bridges the two divalent cations in structures determined in the presence or absence of substrate analogues. One of the phosphoryl oxygen atoms in the PTE-DEP complex is 2.0 A from the alpha-metal ion, while the other oxygen is 2.2 A from the beta-metal ion. The two metal ions are separated by a distance of 4.0 A. A similar structure is observed in the presence of cacodylate. Analogous complexes have previously been observed for the product complexes of isoaspartyl dipeptidase, d-aminoacylase, and dihydroorotase from the amidohydrolase superfamily of enzymes. The experimentally determined structure of the PTE-diethyl phosphate product complex is inconsistent with a recent proposal based upon quantum mechanical/molecular mechanical simulations which postulated the formation of an asymmetrical product complex bound exclusively to the beta-metal ion with a metal-metal separation of 5.3 A. This structure is also inconsistent with a chemical mechanism for substrate hydrolysis that utilizes the bridging hydroxide as a base to abstract a proton from a water molecule loosely associated with the alpha-metal ion. Density functional theory (DFT) calculations support a reaction mechanism that utilizes the bridging hydroxide as the direct nucleophile in the hydrolysis of organophosphate esters by PTE.

Theoretical investigation of the reaction mechanism of the dinuclear zinc enzyme dihydroorotase

The reaction mechanism of the dinuclear zinc enzyme dihydroorotase was investigated by using hybrid density functional theory. This enzyme catalyzes the reversible interconversion of dihydroorotate and carbamoyl aspartate. Two reaction mechanisms in which the important active site residue Asp250 was either protonated or unprotonated were considered. The calculations establish that Asp250 must be unprotonated for the reaction to take place. The bridging hydroxide is shown to be capable of performing nucleophilic attack on the substrate from its bridging position and the role of Zn(beta) is argued to be the stabilization of the tetrahedral intermediate and the transition state leading to it, thereby lowering the barrier for the nucleophilic attack. It is furthermore concluded that the rate-limiting step is the protonation of the amide nitrogen by Asp250 coupled with C-N bond cleavage, which is consistent with previous experimental findings from isotope labeling studies.

The reaction mechanism of paraoxon hydrolysis by phosphotriesterase from combined QM/MM simulations

Molecular dynamics simulations employing combined quantum mechanical and molecular mechanical (QM/MM) potentials have been carried out to investigate the reaction mechanism of the hydrolysis of paraoxon by phosphotriesterase (PTE). We used a dual-level QM/MM approach that synthesizes accurate results from high-level electronic structure calculations with computational efficiency of semiempirical QM/MM potentials for free energy simulations. In particular, the intrinsic (gas-phase) energies of the active site in the QM region are determined by using density functional theory (B3LYP) and second-order Møller-Plesset perturbation theory (MP2) and the molecular dynamics free energy simulations are performed by using the mixed AM1:CHARMM potential. The simulation results suggest a revised mechanism for the phosphotriester hydrolysis mechanism by PTE. The reaction free energy profile is mirrored by structural motions of the binuclear metal center in the active site. The two zinc ions occupy a compact conformation with an average zinc-zinc distance of 3.5 +/- 0.1 A in the Michaelis complex, whereas it is elongated to 5.3 +/- 0.3 A at the transition state and product state. The substrate is loosely bound to the more exposed zinc ion (Znbeta2+) at an average distance of 3.8 A +/- 0.3 A. The P=O bond of the substrate paraoxon is activated by adopting a tight coordination to the Znbeta2+, releasing the coordinate to the bridging hydroxide ion and increasing its nucleophilicity. It was also found that a water molecule enters into the binding pocket of the loosely bound binuclear center, originally occupied by the nucleophilic hydroxide ion. We suggest that the proton of this water molecule is taken up by His254 at low pH or released to the solvent at high pH, resulting in a hydroxide ion that pulls the Znbeta2+ ion closer to form the compact configuration and restores the resting state of the enzyme.