Dissociative Transition State Research Articles

Density functional theory study of sulfur tolerance of CO adsorption and dissociation on Rh–Ni binary metals

The effect of Ni addition to improve the sulfur tolerance of a Rh catalyst for CO dissociation was studied using density functional theory (DFT) methods. Adsorption and dissociation were considered over the (1 1 1) surfaces of binary Rh 1Ni 2 and Rh 2Ni 1 metals with comparison to pure Rh and Ni surfaces. Sulfur adsorption on the Rh 1Ni 2(1 1 1) surface is 0.21 eV more endothermic than on the Rh(1 1 1) surface, suggesting that a Rh 1Ni 2 bimetallic catalyst has a higher sulfur tolerance than pure Rh catalysts due to a lower surface coverage of the sulfur poison. To compare catalytic activity in the presence of adsorbed sulfur, the CO dissociation rates over the binary and pure metals were calculated with 1/9 sulfur coverage. CO dissociation is fastest on the pure Rh surface under sulfur-free conditions, whereas among sulfur poisoned surfaces, the Rh 1Ni 2 surface shows the fastest CO dissociation rate. The CO dissociation barrier on Rh 1Ni 2 is destabilized less by a S coadsorbate than for the other metals. The addition of Ni atoms to a Rh catalyst improves the sulfur tolerance of the catalyst for CO dissociation by minimizing the repulsion between the adsorbed S atom and the CO dissociation transition state, as evidenced through a projected density of states analysis. The Rh 1Ni 2(2 2 1) stepped surface also shows a lower activation barrier and higher CO dissociation rate in the presence of sulfur than the Rh(2 2 1) stepped surface.

ChemInform Abstract: Synthesis and Hydrolysis Studies of Phosphonopyruvate.

Phosphonopyruvate (1) is prepared in 56% yield from triethyl phosphonopyruvate (3) and pKa values of 1.63, 2.40 and 7.41 determined. The stability of phosphonopyruvate is monitored at 75 °C in aqueous buffers over the pH range 0.6 to 8.3. The only products detected are pyruvate and inorganic phosphate. The pH–rate profile shows that the C–P bond is hydrolysed fastest in the monoanion and dianion of phosphonopyruvate, with rate constants of 1.94 × 10–4 s–1 for the monoanion and 1.05 × 10–4 s–1 for the dianion. For the reaction at pH 4.75, where the dianion predominates, ΔH‡ is 114.0 kJ mor–1 and ΔS‡ is 5.4 J mol–1 K–1, the deuterium isotope effect, kH/kD is 1.08 ± 0.04, and there is non-selective phosphorylation of methanol and H2O in an equimolar solution of these solvents. These results are consistent with a mechanism of hydrolysis of both the monoanion and dianion that involves a very largely dissociative transition state with monomeric metaphosphate character.

Phosphoenolpyruvate and Mg2+ Binding to Pyruvate Kinase Monitored by Infrared Spectroscopy

Structural changes in rabbit muscle pyruvate kinase (PK) induced by phosphoenolpyruvate (PEP) and Mg2+ binding were studied by attenuated total reflection Fourier transform infrared spectroscopy in combination with a dialysis accessory. The experiments indicated a largely preserved secondary structure upon PEP and Mg2+ binding but also revealed small backbone conformational changes of PK involving all types of secondary structure. To assess the effect of the protein environment on the bound PEP, we assigned and evaluated the infrared absorption bands of bound PEP. These were identified using 2,3-13C2-labeled PEP. We obtained the following assignments: 1589 cm−1 (antisymmetric carboxylate stretching vibration); 1415 cm−1 (symmetric carboxylate stretching vibration); 1214 cm−1 (C-O stretching vibration); 1124 and 1110 cm−1 (asymmetric PO32- stretching vibrations); and 967 cm−1 (symmetric PO32- stretching vibration). The corresponding band positions in solution are 1567, 1407, 1229, 1107, and 974 cm−1. The differences for bound and free PEP indicate specific interactions between ligand and protein. Quantification of the interactions with the phosphate group indicated that the enzyme environment has little influence on the P-O bond strengths, and that the bridging P-O bond, which is broken in the catalytic reaction, is weakened by <3%. Thus, there is only little distortion toward a dissociative transition state of the phosphate transfer reaction when PEP binds to PK. Therefore, our results are in line with an associative transition state. Carboxylate absorption bands indicated a maximal shortening of the length of the shorter C-O bond by 1.3 pm. PEP bound to PK in the presence of the monovalent ion Na+ exhibited the same band positions as in the presence of K+, indicating very similar interaction strengths between ligand and protein in both cases.

Methane dissociation and adsorption on Ni(111), Pt(111), Ni(100), Pt(100), and Pt(110)-(1×2): Energetic study

We use density functional theory to examine 24 transition states for methane dissociation on five different metal surfaces. In our calculations, the nonlocal exchange-correlation effects are treated within the generalized gradient approximation using the Perdew-Burke-Ernzerhof functional. In all cases, the minimum energy path for dissociation is over a top site. The barriers are large, 0.66-1.12 eV, and relatively insensitive to the rotational orientation of the (nonreacting) methyl group and the azimuthal orientation of the reactive C-H bond. There is a strong preference on the Pt surfaces for the methyl fragment to bond on the top site, while on the Ni surfaces there is a preference for the hollow or bridge sites. Thus, during the dissociation on Pt, only the low mass H atom needs to significantly move or tunnel, while on Ni, both the dissociating H and the methyl fragment move away from the top site. For all 24 configurations there is a strong force at the transition state to pucker the metal atom over which the reaction occurs. The resulting magnitude of the variation in the barrier height with the motion of this atom varies a bit from surface-to-surface, but is of the order of 1 eV/A. We derive a model for the effective reaction barrier height that includes the effects of lattice motion and substrate temperature and compare with recent experiments and other theoretical studies.

A Substructure Combination Strategy To Create Potent and Selective Transthyretin Kinetic Stabilizers That Prevent Amyloidogenesis and Cytotoxicity

Transthyretin aggregation-associated proteotoxicity appears to cause several human amyloid diseases. Rate-limiting tetramer dissociation and monomer misfolding of transthyretin (TTR) occur before its aggregation into cross-beta-sheet amyloid fibrils. Small molecule binding to and preferential stabilization of the tetrameric state of TTR over the dissociative transition state raises the kinetic barrier for dissociation, imposing kinetic stabilization on TTR and preventing aggregation. This is an effective strategy to halt neurodegeneration associated with polyneuropathy, according to recent placebo-controlled clinical trial results. In three recent papers, we systematically ranked possibilities for the three substructures composing a typical TTR kinetic stabilizer, using fibril inhibition potency and plasma TTR binding selectivity data. Herein, we have successfully employed a substructure combination strategy to use these data to develop potent and selective TTR kinetic stabilizers that rescue cells from the cytotoxic effects of TTR amyloidogenesis. Of the 92 stilbene and dihydrostilbene analogues synthesized, nearly all potently inhibit TTR fibril formation. Seventeen of these exhibit a binding stoichiometry of >1.5 of a maximum of 2 to plasma TTR, while displaying minimal binding to the thyroid hormone receptor (<20%). Six analogues were definitively categorized as kinetic stabilizers by evaluating dissociation time-courses. High-resolution TTR.(kinetic stabilizer)(2) crystal structures (1.31-1.70 A) confirmed the anticipated binding orientation of the 3,5-dibromo-4-hydroxyphenyl substructure and revealed a strong preference of the isosteric 3,5-dibromo-4-aminophenyl substructure to bind to the inner thyroxine binding pocket of TTR.

Transition state analogues in structures of ricin and saporin ribosome-inactivating proteins

Ricin A-chain (RTA) and saporin-L1 (SAP) catalyze adenosine depurination of 28S rRNA to inhibit protein synthesis and cause cell death. We present the crystal structures of RTA and SAP in complex with transition state analogue inhibitors. These tight-binding inhibitors mimic the sarcin-ricin recognition loop of 28S rRNA and the dissociative ribocation transition state established for RTA catalysis. RTA and SAP share unique purine-binding geometry with quadruple pi-stacking interactions between adjacent adenine and guanine bases and 2 conserved tyrosines. An arginine at one end of the pi-stack provides cationic polarization and enhanced leaving group ability to the susceptible adenine. Common features of these ribosome-inactivating proteins include adenine leaving group activation, a remarkable lack of ribocation stabilization, and conserved glutamates as general bases for activation of the H(2)O nucleophile. Catalytic forces originate primarily from leaving group activation evident in both RTA and SAP in complex with transition state analogues.

Side chain specificity of ADP-ribosylation by a sirtuin.

Endogenous mono-ADP-ribosylation in eukaryotes is involved in regulating protein synthesis, signal transduction, cytoskeletal integrity, and cell proliferation, although few cellular ADP-ribosyltransferases have been identified. The sirtuins constitute a highly conserved family of protein deacetylases, and several family members have also been reported to perform protein ADP-ribosylation. We characterized the ADP-ribosylation reaction of the nuclear sirtuin homolog Trypanosoma brucei SIR2-related protein 1 (TbSIR2RP1) on both acetylated and unacetylated substrates. We demonstrated that an acetylated substrate is not required for ADP-ribosylation to occur, indicating that the reaction performed by TbSIR2RP1 is a genuine enzymatic reaction and not a side reaction of deacetylation. Biochemical and MS data showed that arginine is the major ADP-ribose acceptor for unacetylated substrates, whereas arginine does not appear to be the major ADP-ribose acceptor in reactions with acetylated histone H1.1. We performed combined ab initio quantum mechanical/molecular mechanical molecular dynamics simulations, which indicated that sirtuin ADP-ribosylation at arginine is energetically feasible, and involves a concerted mechanism with a highly dissociative transition state. In comparison with the corresponding nicotinamide cleavage in the deacetylation reaction, the simulations suggest that sirtuin ADP-ribosylation would be several orders slower but less sensitive to nicotinamide inhibition, which is consistent with experimental results. These results suggest that TbSIR2RP1 can perform ADP-ribosylation using two distinct mechanisms, depending on whether or not the substrate is acetylated.

Transition States

Transition States

A Computational Study of the Phosphorylation Mechanism of the Insulin Receptor Tyrosine Kinase

Although various groups have studied the phosphorylation mechanism of the insulin receptor tyrosine kinase (IRK), an unanimous picture has not yet emerged. In this work, we performed a computational study to gain further insights. We first built a structural model of the reactant complex with the guide of several crystal structures and previous computational studies of the cyclic AMP-dependent protein kinase. We then optimized the structure by performing geometry optimization using a quantum mechanical model containing nearly 300 atoms. A reaction path was then traced between the reactant and the product by using a multiple coordinate-driven method. The calculations mapped out a sequence of structural changes depicting the conversion of the reactant to the product. Analysis of the structural changes revealed the formation of a dissociative transition state and the involvement of a proton transfer from the hydroxyl group of the tyrosyl residue of the peptide substrate to a conserved aspartate in the active site of the enzyme. The proton transfer began well before the transition state was reached and finished only shortly before the product was completely formed. In addition, the formation of a hydrogen bonding network among Arg1136, Asp1132, the gamma-phosphate of ATP, and the tyrosine residue of the substrate appeared to hold the latter two in a near-attack position for reaction. The model estimated a reaction barrier of 14 kcal/mol, semiquantitatively in accord with experiment.

Binding Mechanism of an SH3 Domain Studied by NMR and ITC

Complexes between Src-homology 3 domains and proline-rich target peptides can have lifetimes on the order of milliseconds, making them too short-lived for kinetic characterization by conventional methods. Nuclear magnetic resonance (NMR) dynamics experiments are ideally suited to study such rapid binding equilibria, and additionally provide information on partly bound intermediate states. We used NMR together with isothermal titration calorimetry (ITC) to characterize the interaction of the SH3 domain from the Fyn tyrosine kinase with a 12-residue peptide at temperatures between 10 and 50 degrees C. NMR data at all temperatures are consistent with an effectively two-state binding reaction, such that any intermediates are either very weakly populated or exchange extremely rapidly with the free or bound forms. Dissociation rate constants, determined by CPMG and ZZ-exchange NMR experiments, range from k(off)(10 degrees C) = 4.5 s(-1) to k(off)(50 degrees C) = 331 s(-1). ITC data at all temperatures follow a simple two-state interaction model. Binding is favored enthalpically, with a dissociation enthalpy, DeltaH(D)(30 degrees C) = 15.4 kcal mol(-1), and disfavored entropically, with a dissociation entropy, DeltaS(D)(30 degrees C) = 20.0 cal mol(-1) K(-1). The free protein and peptide have significantly higher heat capacity than the bound complex, DeltaC(p) = 352 cal mol(-1) K(-1), which is consistent with the largely hydrophobic character of the binding interface. An Eyring plot of k(off) values gives an activation enthalpy of dissociation, DeltaH(D)(double dagger)(30 degrees C) = 19.3 kcal mol(-1) and exhibits slight curvature consistent with the ITC-derived value of DeltaC(p). The curvature suggests that nonpolar residues of the hydrophobic interface are solvated in the transition state for dissociation. Association rate constants were calculated as k(on) = k(off)/K(D), and range from k(on)(10 degrees C) = 1.03 x 10(8) M(-1) s(-1) to k(on)(50 degrees C) = 2.0 x 10(8) M(-1) s(-1), with an apparent activation enthalpy, DeltaH(A)(double dagger) = 3.4 kcal mol(-1). Both the magnitudes and temperature dependence of k(on) values are consistent with a diffusion-limited association mechanism. The combination of NMR and ITC data sheds light on how the Fyn tyrosine kinase is activated by binding to proline-rich targets, and represents a powerful approach for characterizing transient protein/ligand interactions.

Molecular and Atomic Hydrogen Interactions with Au−Ir Near-Surface Alloys

Using a combination of density functional theory and experimental electrochemical methods, we have investigated Au−Ir near-surface alloys and their ability to take up hydrogen (H). Despite the relative instability of H in the subsurface of pure Au(111) and Ir(111), H in the subsurface of a near-surface alloy made of a monolayer of Au over Ir(111) (Au*/Ir) is stabilized compared to surface H. While dissociative adsorption of H2 is activated on this alloy surface, the presence of subsurface H stabilizes the transition state for further H2 dissociation. This is explained by the upshift in the d-band center of the surface Au atoms induced by the presence of subsurface H in Au*/Ir. Submonolayers of Au on Ir(111) stabilize H at the Au−Ir interface while allowing for nonactivated H2 dissociation at the exposed Ir atoms. Synthesis of submonolayer alloys of Au on Ir(111) shows that Au does not block hydrogen adsorption on the Ir surface, despite the inability of pure Au to take up hydrogen under these conditions. The possible existence of subsurface hydrogen in these alloys may explain this unexpected behavior.

Hydrolysis of 8-Quinolyl Phosphate Monoester: Kinetic and Theoretical Studies of the Effect of Lanthanide Ions

8-Quinolyl phosphate (8QP) in the presence of the trivalent lanthanide ions (Ln = La, Sm, Eu, Tb, and Er) forms a [Ln x 8QP]+ complex where the lanthanide ion catalyzes hydrolysis of 8QP. In reactions with Tb3+ or Er3+, there is evidence of limited intervention by a second lanthanide ion. Rate constants are increased by more than 10(7)-fold, and kinetic data and B3LYP/ECP calculations indicate that the effects are largely driven by leaving group and metaphosphate ion stabilization. The lanthanides favor a single-step D(N)A(N) mechanism with a dissociative transition state, with limited nucleophilic assistance, consistent with the low hydroxide ion dependence and the small kinetic effect of Ln3+ radii.

Mechanism of ATP turnover inhibition in the EJC

The exon junction complex (EJC) is deposited onto spliced mRNAs and is involved in many aspects of mRNA function. We have recently reconstituted and solved the crystal structure of the EJC core made of MAGOH, Y14, the most conserved portion of MLN51, and the DEAD-box ATPase eIF4AIII bound to RNA in the presence of an ATP analog. The heterodimer MAGOH/Y14 inhibits ATP turnover by eIF4AIII, thereby trapping the EJC core onto RNA, but the exact mechanism behind this remains unclear. Here, we present the crystal structure of the EJC core bound to ADP-AIF(3), the first structure of a DEAD-box helicase in the transition-mimicking state during ATP hydrolysis. It reveals a dissociative transition state geometry and suggests that the locking of the EJC onto the RNA by MAGOH/Y14 is not caused by preventing ATP hydrolysis. We further show that ATP can be hydrolyzed inside the EJC, demonstrating that MAGOH/Y14 acts by locking the conformation of the EJC, so that the release of inorganic phosphate, ADP, and RNA is prevented. Unifying features of ATP hydrolysis are revealed by comparison of our structure with the EJC-ADPNP structure and other helicases. The reconstitution of a transition state mimicking complex is not limited to the EJC and eIF4AIII as we were also able to reconstitute the complex Dbp5-RNA-ADP-AlF(3), suggesting that the use of ADP-AlF(3) may be a valuable tool for examining DEAD-box ATPases in general.

Highly dissociative and concerted mechanism for the nicotinamide cleavage reaction in Sir2Tm enzyme suggested by ab initio QM/MM molecular dynamics simulations.

Sir2 enzymes catalyze the NAD+-dependent protein deacetylation and play critical roles in epigenetics, cell death, and lifespan regulation. In spite of a current flurry of experimental studies, the catalytic mechanism for this unique and important class of enzymes remains elusive. Employing on-the-fly Born-Oppenheimer molecular dynamics simulations with the B3LYP/6-31G(d) QM/MM potential and the umbrella sampling method, we have characterized the initial step of the Sir2Tm-catalyzed reaction, which is also the most controversial portion of its mechanism. Our results indicate that the nicotinamide cleavage reaction employs a highly dissociative and concerted displacement mechanism: the cleavage of the glycosidic bond is facilitated by the nucleophilic participation of the acetyl-lysine, and the dissociative transition state has a significant oxocarbenium ion character. During this step of the reaction, the Sir2Tm enzyme strongly stabilizes the covalent O-alkylamidate intermediate whereas its effect on the transition state is quite minimal. In addition, functional roles of key residues and motifs have been elucidated. This work further demonstrates the feasibility and applicability of the state-of-the-art ab initio QM/MM molecular dynamics approach in simulating enzyme reactions.

Theoretical study of adsorption and dissociation of NH3 on the Ir{110}(1×2) surface

The adsorption and dissociation of NH3 on Ir{110}(1×2) have been investigated using the density-functional calculations at a coverage of 0.25 ML. The adsorption sites, energy, and geometries were obtained for NH3, NH2, and H adsorptions on the surface. The transition state for NH3 dissociation on Ir{110}(1×2) was also identified. It was found that NH3 is adsorbed preferentially at the ridge atop site, while NH2 and H are adsorbed at the ridge bridge site. The activation barrier of NH3 dissociation is 78.4 kJ/mol, which is very close to the NH3 adsorption energy of 90.0 kJ/mol. This indicates that the desorption and dissociation of NH3 on Ir{110}(1×2) are very competitive, which is consistent with the recent experimental results.

Associative Versus Dissociative Mechanisms of Phosphate Monoester Hydrolysis: On the Interpretation of Activation Entropies

Phosphate monoester and anhydride hydrolysis is ubiquitous in biology, being involved in, amongst other things, signal transduction, energy production, and the regulation of protein function. Therefore, this reaction has understandably been the focus of intensive research. Nevertheless, the precise mechanism by which phosphate monoester hydrolysis proceeds remains controversial. Traditionally, it has been assumed and frequently implied that a near-zero activation entropy is indicative of a dissociative pathway. Herein, we examine free-energy surfaces for the hydrolysis of the methyl phosphate dianion and the methyl pyrophosphate trianion in aqueous solution. In both cases, the reaction can proceed through either compact or expansive concerted (A(N)D(N)) transition states, with fairly similar barriers. We have evaluated the activation entropies for each transition state and demonstrate that both associative and dissociative transition states have near-zero entropies of activation that are in good agreement with experimental values. Therefore, we believe that the activation entropy alone is not a useful diagnostic tool, as it depends not only on bond orders at the transition state, but also on other issues that include (but are not limited to) steric factors determining the configurational volumes available to reactants during the reaction, solvation and desolvation effects that may be associated with charge redistribution upon approaching the transition state and entropy changes associated with intramolecular degrees of freedom as the transition state is approached.

Femtosecond real-time probing of transition state dynamics in a surface photoreaction: Methyl desorption from CH3I on MgO(100)

A novel experimental approach to the investigation of surface adsorbate reaction dynamics is presented. The direct time-resolved monitoring of the surface reaction transition state and product formation dynamics were accomplished via pump-probe mass spectrometry. As an example, methyl iodide molecules adsorbed at submonolayer coverage on an ultrathin magnesia film on Mo(100) were photoexcited to the A-band by ultrafast laser pulse irradiation. Employing time-delayed multiphoton ionization the dynamics of the dissociative methyl iodide transition state and of the emerging methyl photoproduct could be detected with femtosecond resolution. The reaction times deduced from the temporal evolution of the methyl ion mass signal indicate a strong interaction of the methyl fragment with the substrate surface prior to desorption.

Calculating Geochemical Reaction Pathways - Exploration of the Inner-Sphere Water Exchange Mechanism in Al(H2O)63+(aq) + nH2O with ab Initio Calculations and Molecular Dynamics

We have simulated exchange of inner-sphere and bulk water molecules for different sizes of Al3+(aq) clusters, Al(H2O)63+ + nH2O for n = 0, 1, 6, or 12, with ab initio and molecular dynamics simulations, in order to understand how robust the ab initio method is for identifying hydrolytic reaction pathways of particular importance to geochemistry. In contrast to many interfacial reactions, this particular elementary reaction is particularly simple and well-constrained by experiment. Nevertheless, we find that a rich array of parallel reaction pathways depend sensitively on the details of the solvation sphere and structure and that larger clusters are not necessarily better. Inner-sphere water exchange in Al3+(aq) may occur through two Langford-Gray dissociative pathways, one in which the incoming and outgoing waters are cis, the other in which they are trans to one another. A large majority of exchanges in the molecular dynamics simulations occurred via the trans mechanism, in contrast to the predictions of the ab initio method. In Al(H2O)63+ + H2O, the cis mechanism has a transition state of 84.3 kJ/mol, which is in good agreement with previous experimental and ab initio results, while the trans mechanism has only a saddle point with two negative frequencies, not a transition state, at 89.7 kJ/mol. In addition to the exchange mechanisms, dissociation pathways could be identified that were considerably lower in energy than experiment and varied considerably between 60 and 100 kJ/mol, depending on the particular geometry and cluster size, with no clear relation between the two. Ab initio calculations using large clusters with full second coordination spheres (n = 12) were unable to find dissociation or exchange transition states because the network of hydrogen bonds in the second coordination sphere was too rigid to accommodate the outgoing inner-sphere water. Our results indicate that caution should surround ab initio simulation of complicated dynamic processes such as hydrolysis, ion exchange, and interfacial reactions that involve several steps. Dynamic methods of simulation need to accompany static methods such as ab initio calculation, and it is best to consider simulated pathways as hypotheses to be tested experimentally rather than definitive properties of the reaction.

Ab Initio QM/MM Studies of the Phosphoryl Transfer Reaction Catalyzed by PEP Mutase Suggest a Dissociative Metaphosphate Transition State

The interconversion between phosphoenolpyruvate (PEP) and phosphonopyruvate (P-pyr) catalyzed by PEP mutase is investigated using an ab initio QM/MM method with the QM region treated at the B3LYP/6-31G* level of theory. Two-dimensional minimum energy path calculations were carried out for both the wild-type enzyme and the N122A mutant. The calculations suggest a dissociative transition state featuring metaphosphate and Mg(2+)-coordinating pyruvate enolate, stabilized by an extensive hydrogen bond network involving Asn122, Ser123, Arg159, His190, Ser46, and Leu48. It is also found that a substantial conformational change in the pyruvyl group is required for the interconversion.

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.