Mechanism of substrate dephosphorylation in lowMr protein tyrosine phosphatase

Evaluating transition state structures of vanadium-phosphatase protein complexes using shape analysis

Shape analysis of coordination complexes is well-suited to evaluate the subtle distortions in the trigonal bipyramidal (TBPY-5) geometry of vanadium coordinated in the active site of phosphatases and characterized by X-ray crystallography. Recent studies using the tau (τ) analysis support the assertion that vanadium is best described as a trigonal bipyramid, because this geometry is the ideal transition state geometry of the phosphate ester substrate hydrolysis (C.C. McLauchlan, B.J. Peters, G.R. Willsky, D.C. Crans, Coord. Chem. Rev. http://dx.doi.org/10.1016/j.ccr.2014.12.012 ; D.C. Crans, M.L. Tarlton, C.C. McLauchlan, Eur. J. Inorg. Chem. 2014, 4450-4468). Here we use continuous shape measures (CShM) analysis to investigate the structural space of the five-coordinate vanadium-phosphatase complexes associated with mechanistic transformations between the tetrahedral geometry and the five-coordinate high energy TBPY-5 geometry was discussed focusing on the protein tyrosine phosphatase 1B (PTP1B) enzyme. No evidence for square pyramidal geometries was observed in any vanadium-protein complexes. The shape analysis positioned the metal ion and the ligands in the active site reflecting the mechanism of the cleavage of the organic phosphate in a phosphatase. We identified the umbrella distortions to be directly on the reaction path between tetrahedral phosphate and the TBPY-5-types of high-energy species. The umbrella distortions of the trigonal bipyramid are therefore identified as being the most relevant types of transition state structures for the phosphoryl group transfer reactions for phosphatases and this may be related to the possibility that vanadium is an inhibitor for enzymes that support both exploded and five-coordinate transition states.

Why nature really chose phosphate

Phosphoryl transfer plays key roles in signaling, energy transduction, protein synthesis, and maintaining the integrity of the genetic material. On the surface, it would appear to be a simple nucleophile displacement reaction. However, this simplicity is deceptive, as, even in aqueous solution, the low-lying d-orbitals on the phosphorus atom allow for eight distinct mechanistic possibilities, before even introducing the complexities of the enzyme catalyzed reactions. To further complicate matters, while powerful, traditional experimental techniques such as the use of linear free-energy relationships (LFER) or measuring isotope effects cannot make unique distinctions between different potential mechanisms. A quarter of a century has passed since Westheimer wrote his seminal review, 'Why Nature Chose Phosphate' (Science 235 (1987), 1173), and a lot has changed in the field since then. The present review revisits this biologically crucial issue, exploring both relevant enzymatic systems as well as the corresponding chemistry in aqueous solution, and demonstrating that the only way key questions in this field are likely to be resolved is through careful theoretical studies (which of course should be able to reproduce all relevant experimental data). Finally, we demonstrate that the reason that nature really chose phosphate is due to interplay between two counteracting effects: on the one hand, phosphates are negatively charged and the resulting charge-charge repulsion with the attacking nucleophile contributes to the very high barrier for hydrolysis, making phosphate esters among the most inert compounds known. However, biology is not only about reducing the barrier to unfavorable chemical reactions. That is, the same charge-charge repulsion that makes phosphate ester hydrolysis so unfavorable also makes it possible to regulate, by exploiting the electrostatics. This means that phosphate ester hydrolysis can not only be turned on, but also be turned off, by fine tuning the electrostatic environment and the present review demonstrates numerous examples where this is the case. Without this capacity for regulation, it would be impossible to have for instance a signaling or metabolic cascade, where the action of each participant is determined by the fine-tuned activity of the previous piece in the production line. This makes phosphate esters the ideal compounds to facilitate life as we know it.

Mechanism of Cdc25B phosphatase with the small molecule substrate p-nitrophenyl phosphate from QM/MM-MFEP calculations

Cdc25B is a dual-specificity phosphatase that catalyzes the dephosphorylation of the Cdk2/CycA protein complex. This enzyme is an important regulator of the human cell cycle and has been identified as a potential anticancer target. In general, protein tyrosine phosphatases are thought to bind the dianionic form of the phosphate and employ general acid catalysis via the Asp residue in the highly conserved WPD-loop. However, the Cdc25 phosphatases form a special subfamily based on their distinct differences from other protein tyrosine phosphatases. Although Cdc25B contains the (H/V)CX(5)R catalytic motif present in all other protein tyrosine phosphatases, it lacks an analogous catalytic acid residue. No crystallographic data currently exist for the complex of Cdc25B with Cdk2/CycA, so in addition to its natural protein substrate, experimental and theoretical studies are often carried out with small molecule substrates. In an effort to gain understanding of the dephosphorylation mechanism of Cdc25B with a commonly used small molecule substrate, we have performed simulations of the rate-limiting step of the reaction catalyzed by Cdc25B with the substrate p-nitrophenyl phosphate using the recently developed QM/MM Minimum Free Energy Path method (Hu et al. J. Chem. Phys. 2008, 034105). We have simulated the first step of the reaction with both the monoanionic and the dianionic forms of the substrate, and our calculations favor a mechanism involving the monoanionic form. Thus, Cdc25 may employ a unique dephosphorylation mechanism among protein tyrosine phosphatases, at least in the case of the small molecule substrate p-nitrophenyl phosphate.

Theoretical study of specific hydrogen-bonding effects on the bridging P-OR bond strength of phosphate monoester dianions

It has been proposed that the driving force for the initial phosphoryl transfer step of protein tyrosine phosphatases (PTPases) could be activation of the substrate ROPO32- by means of an enforced hydrogen-bonding interaction between an aspartic general acid and the bridging oxygen atom O (Zhang et al. Biochemistry 1995, 34, 16088-16096). The potential catalytic effect of this type of interaction, with regard to P-OR bond cleavage, was investigated computationally through simple model systems in which an efficient intramolecular hydrogen bond can take place between a H-bond donor group and the bridging oxygen atom of the dianionic phosphate. The dielectric effect of the environment (epsilon = 1, 4, and 78) was also explored. The results indicate that this interaction causes significant lengthenings of the scissile P-OR bond in all media but with more extreme effects observed in the low dielectric fields epsilon = 1 and epsilon = 4. It is interesting that, in all cases examined, this interaction actually contributes to stabilize the reactant state while causing its P-OR bond to lengthen. Overall, our results support the idea that this specific hydrogen-bonding situation might well be used by PTPases as an important driving force for promoting phosphoryl transfer reactions through highly dissociative transition states.

The role of metal ions in phosphate ester hydrolysis

Many phosphatases make use of metal ions to aid catalysis of phosphate ester hydrolysis. Here, we investigate the impact of metal ions on the potential energy surface (PES), and hence the preferred reaction mechanism, for a simple model for hydrolysis of phosphate ester monoanions. We show that, while both associative (A(N) + D(N)) and dissociative (D(N) + A(N)) mechanisms are represented on the potential energy surfaces both in the presence and absence of metal ions, the D(N) + A(N) process is favoured when there are no metal ions present and the A(N) + D(N) process is favoured in the presence of two metal ions. A concerted (A(N)D(N)) process is also available in the presence of two metal ions, but proceeds via a high-energy transition state. In the presence of only a single metal ion the A(N)D(N) process is the most favoured, but still proceeds via a high-energy transition state. Thus, we conclude that metallo-enzyme phosphatases are likely to utilise an associative process, while those that function without metal ions may well follow a dissociative process.

Computer simulation of the chemical catalysis of DNA polymerases: discriminating between alternative nucleotide insertion mechanisms for T7 DNA polymerase

Understanding the chemical step in the catalytic reaction of DNA polymerases is essential for elucidating the molecular basis of the fidelity of DNA replication. The present work evaluates the free energy surface for the nucleotide transfer reaction of T7 polymerase by free energy perturbation/empirical valence bond (FEP/EVB) calculations. A key aspect of the enzyme simulation is a comparison of enzymatic free energy profiles with the corresponding reference reactions in water using the same computational methodology, thereby enabling a quantitative estimate for the free energy of the nucleotide insertion reaction. The reaction is driven by the FEP/EVB methodology between valence bond structures representing the reactant, pentacovalent intermediate, and the product states. This pathway corresponds to three microscopic chemical steps, deprotonation of the attacking group, a nucleophilic attack on the P(alpha) atom of the dNTP substrate, and departure of the leaving group. Three different mechanisms for the first microscopic step, the generation of the RO(-) nucleophile from the 3'-OH hydroxyl of the primer, are examined: (i) proton transfer to the bulk solvent, (ii) proton transfer to one of the ionic oxygens of the P(alpha) phosphate group, and (iii) proton transfer to the ionized Asp654 residue. The most favorable reaction mechanism in T7 pol is predicted to involve the proton transfer to Asp654. This finding sheds light on the long standing issue of the actual role of conserved aspartates. The structural preorganization that helps to catalyze the reaction is also considered and analyzed. The overall calculated mechanism consists of three subsequent steps with a similar activation free energy of about 12 kcal/mol. The similarity of the activation barriers of the three microscopic chemical steps indicates that the T7 polymerase may select against the incorrect dNTP substrate by raising any of these barriers. The relative height of these barriers comparing right and wrong dNTP substrates should therefore be a primary focus of future computational studies of the fidelity of DNA polymerases.

Density functional study of the mechanism of a tyrosine phosphatase: I. Intermediate formation

The first step in the catalytic mechanism of a protein tyrosine phosphatase, the transfer of a phosphate group from the phosphotyrosine substrate to a cysteine side chain of the protein to form a phosphoenzyme intermediate, has been studied by combining density functional calculations of an active-site cluster with continuum electrostatic descriptions of the solvent and the remainder of the protein. This approach provides the high level of quantum chemical methodology needed to adequately model phosphotransfer reactions with a reasonable description of the environment around the active site. Energy barriers and geometries along a reaction pathway are calculated. In the literature, mechanisms assuming both a monoanionic and a dianionic substrate have been proposed; this disagreement is addressed by performing calculations for both possibilities. For the dianionic substrate, a dissociative reaction pathway with early proton transfer to the leaving group and a 9 kcal/mol energy barrier is predicted (the experimental estimate is ca. 14 kcal/mol), while for the monoanionic substrate, an associative pathway with late proton transfer and a 22 kcal/mol energy barrier is predicted. These results, together with a review of experimental evidence, support the dianionic-substrate/dissociative-pathway alternative. The relationship between a dianionic or monoanionic substrate and a dissociative or associative pathway, respectively, can be understood in terms of classical organic chemical reaction pathways.

Simulating enzyme reactions: challenges and perspectives

Elucidating how enzymes enhance the rates of the reactions that they catalyze is a major goal of contemporary biochemistry, and it is an area in which computational and theoretical techniques can make a major contribution. This article outlines some of the processes that need to be investigated if enzyme catalysis is to be understood, reviews the current state-of-the-art in enzyme simulation work, and highlights challenges for the future.

Insights into the structure, solvation, and mechanism of ArsC arsenate reductase, a novel arsenic detoxification enzyme

BACKGROUND

In Escherichia coli bearing the plasmid R773, resistance to arsenite, arsenate, antimonite, and tellurite is conferred by the arsRDABC plasmid operon that codes for an ATP-dependent anion pump. The product of the arsC gene, arsenate reductase (ArsC), is required to efficiently catalyze the reduction of arsenate to arsenite prior to extrusion.

RESULTS

Here, we report the first X-ray crystal structures of ArsC at 1.65 A and of ArsC complexed with arsenate and arsenite at 1.26 A resolution. The overall fold is unique. The native structure shows sulfate and sulfite ions binding in the active site as analogs of arsenate and arsenite. The covalent adduct of arsenate with Cys-12 in the active site of ArsC, which was analyzed in a difference map, shows tetrahedral geometry with a sulfur-arsenic distance of 2.18 A. However, the corresponding adduct with arsenite binds as a hitherto unseen thiarsahydroxy adduct. Finally, the number of bound waters (385) in this highly ordered crystal structure approaches twice the number expected at this resolution for a structure of 138 ordered residues.

CONCLUSIONS

Structural information from the adduct of ArsC with its substrate (arsenate) and with its product (arsenite) together with functional information from mutational and biochemical studies on ArsC suggest a plausible mechanism for the reaction. The exceptionally well-defined water structure indicates that this crystal system has precise long-range order within the crystal and that the upper limit for the number of bound waters in crystal structures is underestimated by the structures in the Protein Data Bank.

The mechanism of the phosphoryl transfer catalyzed by Yersinia protein-tyrosine phosphatase: a computational and isotope effect study

In order to evaluate various mechanistic proposals that have been made regarding the mechanism of the first step of the reaction catalyzed by protein-tyrosine phosphatases, new experimental data have been obtained, and some existing data have been carefully reevaluated. New kinetic isotope effect data for the uncatalyzed hydrolysis of p-nitrophenyl phosphate allow a better evaluation of previously reported data from enzymatic reactions with this substrate. The interpretation, and misinterpretation, of pH rate studies is considered. The pathway of this reaction has been modeled computationally and is found to be generally consistent with experimental studies, except for the extent of proton transfer to the leaving group.

The catalytic mechanism of protein tyrosine phosphatases revisited

Experimental and theoretical studies of the catalytic mechanism in protein tyrosine phosphatases and dual specific phosphatases are reviewed. The structural properties of these enzymes contributing to the efficient rate enhancement of phosphate monoester hydrolysis have been established during the last decade. There are, however, uncertainties in the interpretation of available experimental data that make the commonly assumed reaction mechanism somewhat doubtful. Theoretical calculations as well as analysis of crystal structures point towards an alternative interpretation of the ionisation state in the reactive complex.

Prediction of a ligand-induced conformational change in the catalytic core of Cdc25A

The cell cycle control phosphatases Cdc25 are dual specificity phosphatases that dephosphorylate both phosphothreonine and phosphotyrosine residues on their substrate proteins. The determination of the apo-protein structure of Cdc25A revealed that this enzyme has a completely different fold compared to all other phosphatases crystallised to date. The conformation of the active site residues does not seem very suitable for catalysis in this unliganded structure. We have studied some structural features of the Cdc25A apo-structure and a modelled Cdc25A-ligand complex by molecular dynamics simulations. The simulations predict a conformational change in the peptide backbone of the complex, which is not observed in the apo-structure. This ligand-induced conformational change yields a structure that is similar to other protein tyrosine phosphatase-ligand complexes that have been crystallised. The change in conformation takes place in the position between a serine and a glutamic acid residue in the phosphate binding loop. We suggest that this type of conformational change is an important molecular switch in the catalytic process.