Sequence - dynamics - function relationships in protein tyrosine phosphatases.

Protein tyrosine phosphatases (PTPs) are crucial regulators of cellular signaling. Their activity is regulated by the motion of a conserved loop, the WPD-loop, from a catalytically inactive open to a catalytically active closed conformation. WPD-loop motion optimally positions a catalytically critical residue into the active site, and is directly linked to the turnover number of these enzymes. Crystal structures of chimeric PTPs constructed by grafting parts of the WPD-loop sequence of PTP1B onto the scaffold of YopH showed WPD-loops in a wide-open conformation never previously observed in either parent enzyme. This wide-open conformation has, however, been observed upon binding of small molecule inhibitors to other PTPs, suggesting the potential of targeting it for drug discovery efforts. Here, we have performed simulations of both enzymes and show that there are negligible energetic differences in the chemical step of catalysis, but significant differences in the dynamical properties of the WPD-loop. Detailed interaction network analysis provides insight into the molecular basis for this population shift to a wide-open conformation. Taken together, our study provides insight into the links between loop dynamics and chemistry in these YopH variants specifically, and how WPD-loop dynamic can be engineered through modification of the internal protein interaction network.

Insights into the importance of WPD-loop sequence for activity and structure in protein tyrosine phosphatases.

Protein tyrosine phosphatases (PTPs) possess a conserved mobile catalytic loop, the WPD-loop, which brings an aspartic acid into the active site where it acts as an acid/base catalyst. Prior experimental and computational studies, focused on the human enzyme PTP1B and the PTP from Yersinia pestis, YopH, suggested that loop conformational dynamics are important in regulating both catalysis and evolvability. We have generated a chimeric protein in which the WPD-loop of YopH is transposed into PTP1B, and eight chimeras that systematically restored the loop sequence back to native PTP1B. Of these, four chimeras were soluble and were subjected to detailed biochemical and structural characterization, and a computational analysis of their WPD-loop dynamics. The chimeras maintain backbone structural integrity, with somewhat slower rates than either wild-type parent, and show differences in the pH dependency of catalysis, and changes in the effect of Mg2+. The chimeric proteins' WPD-loops differ significantly in their relative stability and rigidity. The time required for interconversion, coupled with electrostatic effects revealed by simulations, likely accounts for the activity differences between chimeras, and relative to the native enzymes. Our results further the understanding of connections between enzyme activity and the dynamics of catalytically important groups, particularly the effects of non-catalytic residues on key conformational equilibria.

Competitive measurement of ß/α naphthyl phosphate catalytic efficiency by phosphatases utilizing quantitative NMR.

The two constitutional isomers of naphthyl phosphate have different steric properties, analogous to those of phosphotyrosine versus phosphoserine/threonine within a peptide or protein. The ratios of their respective rates of hydrolysis, assayed by measuring rates of inorganic phosphate release, have been used to probe the steric requirements around the active sites of many phosphatases in the literature. We report an NMR-based competitive method that is simpler to execute and has other advantages. It directly yields the ratio of catalytic efficiencies (V/K) of the two substrates, a more biologically relevant comparison than the ratio of initial rates (vo) or maximal rates (Vmax). The competitive method ensures that temperature, pH, enzyme and substrate concentrations, and the presence of any potential inhibitors are identical and will not skew the results. The method can be easily applied at any chosen temperature or pH, and to mutants, or under any other condition that might influence protein conformation and, thus, substrate specificity. It provides a facile screening method to select conditions for a detailed phosphopeptide screen to provide deeper insight into substrate preference.

Significant Loop Motions in the SsoPTP Protein Tyrosine Phosphatase Allow for Dual General Acid Functionality.

Conformational dynamics are important factors in the function of enzymes, including protein tyrosine phosphatases (PTPs). Crystal structures of PTPs first revealed the motion of a protein loop bearing a conserved catalytic aspartic acid, and subsequent nuclear magnetic resonance and computational analyses have shown the presence of motions, involved in catalysis and allostery, within and beyond the active site. The tyrosine phosphatase from the thermophilic and acidophilic Sulfolobus solfataricus (SsoPTP) displays motions of its acid loop together with dynamics of its phosphoryl-binding P-loop and the Q-loop, the first instance of such motions in a PTP. All three loops share the same exchange rate, implying their motions are coupled. Further evidence of conformational flexibility comes from mutagenesis, kinetics, and isotope effect data showing that E40 can function as an alternate general acid to protonate the leaving group when the conserved acid, D69, is mutated to asparagine. SsoPTP is not the first PTP to exhibit an alternate general acid (after VHZ and TkPTP), but E40 does not correspond to the sequence or structural location of the alternate general acids in those precedents. A high-resolution X-ray structure with the transition state analogue vanadate clarifies the role of the active site arginine R102, which varied in structures of substrates bound to a catalytically inactive mutant. The coordinated motions of all three functional loops in SsoPTP, together with the function of an alternate general acid, suggest that catalytically competent conformations are present in solution that have not yet been observed in crystal structures.

Plasmodium falciparum finds the winning combination.

After 3 years of laboratory drug pressure in the presence of a picomolar inhibitor, the parasite Plasmodium falciparum developed a combination strategy of gene amplification and mutation to regain viability. The mutation observed led to a dysfunctional enzyme, but new research reveals the clever mechanism behind its success. Not that we needed a reminder of nature's creativity in the time of a pandemic.

Single Residue on the WPD-Loop Affects the pH Dependency of Catalysis in Protein Tyrosine Phosphatases.

Catalysis by protein tyrosine phosphatases (PTPs) relies on the motion of a flexible protein loop (the WPD-loop) that carries a residue acting as a general acid/base catalyst during the PTP-catalyzed reaction. The orthogonal substitutions of a noncatalytic residue in the WPD-loops of YopH and PTP1B result in shifted pH-rate profiles from an altered kinetic pK a of the nucleophilic cysteine. Compared to wild type, the G352T YopH variant has a broadened pH-rate profile, similar activity at optimal pH, but significantly higher activity at low pH. Changes in the corresponding PTP1B T177G variant are more modest and in the opposite direction, with a narrowed pH profile and less activity in the most acidic range. Crystal structures of the variants show no structural perturbations but suggest an increased preference for the WPD-loop-closed conformation. Computational analysis confirms a shift in loop conformational equilibrium in favor of the closed conformation, arising from a combination of increased stability of the closed state and destabilization of the loop-open state. Simulations identify the origins of this population shift, revealing differences in the flexibility of the WPD-loop and neighboring regions. Our results demonstrate that changes to the pH dependency of catalysis by PTPs can result from small changes in amino acid composition in their WPD-loops affecting only loop dynamics and conformational equilibrium. The perturbation of kinetic pK a values of catalytic residues by nonchemical processes affords a means for nature to alter an enzyme's pH dependency by a less disruptive path than altering electrostatic networks around catalytic residues themselves.

Loop Dynamics and Enzyme Catalysis in Protein Tyrosine Phosphatases.

Protein tyrosine phosphatases (PTPs) play an important role in cellular signaling and have been implicated in human cancers, diabetes, and obesity. Despite shared catalytic mechanisms and transition states for the chemical steps of catalysis, catalytic rates within the PTP family vary over several orders of magnitude. These rate differences have been implied to arise from differing conformational dynamics of the closure of a protein loop, the WPD-loop, which carries a catalytically critical residue. The present work reports computational studies of the human protein tyrosine phosphatase 1B (PTP1B) and YopH from Yersinia pestis, for which NMR has demonstrated a link between their respective rates of WPD-loop motion and catalysis rates, which differ by an order of magnitude. We have performed detailed structural analysis, both conventional and enhanced sampling simulations of their loop dynamics, as well as empirical valence bond simulations of the chemical step of catalysis. These analyses revealed the key residues and structural features responsible for these differences, as well as the residues and pathways that facilitate allosteric communication in these enzymes. Curiously, our wild-type YopH simulations also identify a catalytically incompetent hyper-open conformation of its WPD-loop, sampled as a rare event, previously only experimentally observed in YopH-based chimeras. The effect of differences within the WPD-loop and its neighboring loops on the modulation of loop dynamics, as revealed in this work, may provide a facile means for the family of PTP enzymes to respond to environmental changes and regulate their catalytic activities.

Modeling the Alkaline Hydrolysis of Diaryl Sulfate Diesters: A Mechanistic Study.

Phosphate and sulfate esters have important roles in regulating cellular processes. However, while there has been substantial experimental and computational investigation of the mechanisms and the transition states involved in phosphate ester hydrolysis, there is far less work on sulfate ester hydrolysis. Here, we report a detailed computational study of the alkaline hydrolysis of diaryl sulfate diesters, using different DFT functionals as well as mixed implicit/explicit solvation with varying numbers of explicit water molecules. We consider the impact of the computational model on computed linear free-energy relationships (LFER) and the nature of the transition states (TS) involved. We obtain good qualitative agreement with experimental LFER data when using a pure implicit solvent model and excellent agreement with experimental kinetic isotope effects for all models used. Our calculations suggest that sulfate diester hydrolysis proceeds through loose transition states, with minimal bond formation to the nucleophile and bond cleavage to the leaving group already initiated. Comparison to prior work indicates that these TS are similar in nature to those for the alkaline hydrolysis of neutral arylsulfonate monoesters or charged phosphate diesters and fluorophosphates. Obtaining more detailed insights into the transition states involved assists in understanding the selectivity of enzymes that hydrolyze these reactions.

Transition-State Interactions in a Promiscuous Enzyme: Sulfate and Phosphate Monoester Hydrolysis by Pseudomonas aeruginosa Arylsulfatase.

Pseudomonas aeruginosa arylsulfatase (PAS) hydrolyzes sulfate and, promiscuously, phosphate monoesters. Enzyme-catalyzed sulfate transfer is crucial to a wide variety of biological processes, but detailed studies of the mechanistic contributions to its catalysis are lacking. We present linear free energy relationships (LFERs) and kinetic isotope effects (KIEs) of PAS and analyses of active site mutants that suggest a key role for leaving group (LG) stabilization. In LFERs PASWT has a much less negative Brønsted coefficient (ßleaving groupobs-Enz = -0.33) than the uncatalyzed reaction (ßleaving groupobs = -1.81). This situation is diminished when cationic active site groups are exchanged for alanine. The considerable degree of bond breaking during the transition state (TS) is evidenced by an 18Obridge KIE of 1.0088. LFER and KIE data for several active site mutants point to leaving group stabilization by active site K375, in cooperation with H211. 15N KIEs and the increased sensitivity to leaving group ability of the sulfatase activity in neat D2O (Δßleaving groupH-D = +0.06) suggest that the mechanism for S-Obridge bond fission shifts, with decreasing leaving group ability, from charge compensation via Lewis acid interactions toward direct proton donation. 18Ononbridge KIEs indicate that the TS for PAS-catalyzed sulfate monoester hydrolysis has a significantly more associative character compared to the uncatalyzed reaction, while PAS-catalyzed phosphate monoester hydrolysis does not show this shift. This difference in enzyme-catalyzed TSs appears to be the major factor favoring specificity toward sulfate over phosphate esters by this promiscuous hydrolase, since other features are either too similar (uncatalyzed TS) or inherently favor phosphate (charge).

Transition State Analysis of the Reaction Catalyzed by the Phosphotriesterase from Sphingobium sp. TCM1.

Organophosphorus flame retardants are stable toxic compounds used in nearly all durable plastic products and are considered major emerging pollutants. The phosphotriesterase from Sphingobium sp. TCM1 ( Sb-PTE) is one of the few enzymes known to be able to hydrolyze organophosphorus flame retardants such as triphenyl phosphate and tris(2-chloroethyl) phosphate. The effectiveness of Sb-PTE for the hydrolysis of these organophosphates appears to arise from its ability to hydrolyze unactivated alkyl and phenolic esters from the central phosphorus core. How Sb-PTE is able to catalyze the hydrolysis of the unactivated substituents is not known. To interrogate the catalytic hydrolysis mechanism of Sb-PTE, the pH dependence of the reaction and the effects of changing the solvent viscosity were determined. These experiments were complemented by measurement of the primary and secondary 18-oxygen isotope effects on substrate hydrolysis and a determination of the effects of changing the p Ka of the leaving group on the magnitude of the rate constants for hydrolysis. Collectively, the results indicated that a single group must be ionized for nucleophilic attack and that a separate general acid is not involved in protonation of the leaving group. The Brønsted analysis and the heavy atom kinetic isotope effects are consistent with an early associative transition state with subsequent proton transfers not being rate limiting. A novel binding mode of the substrate to the binuclear metal center and a catalytic mechanism are proposed to explain the unusual ability of Sb-PTE to hydrolyze unactivated esters from a wide range of organophosphate substrates.

Catalytic mechanism for the conversion of salicylate into catechol by the flavin-dependent monooxygenase salicylate hydroxylase.

Salicylate hydroxylase (NahG) is a flavin-dependent monooxygenase that catalyzes the decarboxylative hydroxylation of salicylate into catechol in the naphthalene degradation pathway in Pseudomonas putida G7. We explored the mechanism of action of this enzyme in detail using a combination of structural and biophysical methods. NahG shares many structural and mechanistic features with other versatile flavin-dependent monooxygenases, with potential biocatalytic applications. The crystal structure at 2.0 Šresolution for the apo form of NahG adds a new snapshot preceding the FAD binding in flavin-dependent monooxygenases. The kcat/Km for the salicylate reaction catalyzed by the holo form is >105 M-1 s-1 at pH 8.5 and 25 °C. Hammett plots for Km and kcat using substituted salicylates indicate change in rate-limiting step. Electron-donating groups favor the hydroxylation of salicylate by a peroxyflavin to yield a Wheland-like intermediate, whereas the decarboxylation of this intermediate is faster for electron-withdrawing groups. The mechanism is supported by structural data and kinetic studies at different pHs. The salicylate carboxyl group lies near a hydrophobic region that aids decarboxylation. A conserved histidine residue is proposed to assist the reaction by general base/general acid catalysis.

A YopH PTP1B Chimera Shows the Importance of the WPD-Loop Sequence to the Activity, Structure, and Dynamics of Protein Tyrosine Phosphatases.

To study factors that affect WPD-loop motion in protein tyrosine phosphatases (PTPs), a chimera of PTP1B and YopH was created by transposing the WPD loop from PTP1B to YopH. Several subsequent mutations proved to be necessary to obtain a soluble, active enzyme. That chimera, termed chimera 3, retains productive WPD-loop motions and general acid catalysis with a pH dependency similar to that of the native enzymes. Kinetic isotope effects show the mechanism and transition state for phosphoryl transfer are unaltered. Catalysis of the chimera is slower than that of either of its parent enzymes, although its rate is comparable to those of most native PTPs. X-ray crystallography and nuclear magnetic resonance were used to probe the structure and dynamics of chimera 3. The chimera's structure was found to sample an unproductive hyper-open conformation of its WPD loop, a geometry that has not been observed in either of the parents or in other native PTPs. The reduced catalytic rate is attributed to the protein's sampling of this conformation in solution, reducing the fraction in the catalytically productive loop-closed conformation.

Transition States and Control of Substrate Preference in the Promiscuous Phosphatase PP1.

Catalytically promiscuous enzymes are an attractive frontier for biochemistry, because enzyme promiscuities not only plausibly explain enzyme evolution through the mechanism of gene duplication but also could provide an efficient route to changing the catalytic function of proteins by mimicking this evolutionary process. PP1γ is an effectively promiscuous phosphatase for the hydrolysis of both monoanionic and dianionic phosphate ester-based substrates. In addition to its native phosphate monoester substrate, PP1γ catalyzes the hydrolysis of aryl methylphosphonates, fluorophosphate esters, phosphorothioate esters, and phosphodiesters, with second-order rate accelerations that fall within the narrow range of 1011-1013. In contrast to the different transition states in the uncatalyzed hydrolysis reactions of these substrates, PP1γ catalyzes their hydrolysis through similar transition states. PP1γ does not catalyze the hydrolysis of a sulfate ester, which is unexpected. The PP1γ active site is tolerant of variations in the geometry of bound ligands, which permit the effective catalysis even of substrates whose steric requirements may result in perturbations to the positioning of the transferring group, both in the initial enzyme-substrate complex and in the transition state. The conservative mutation of arginine 221 to lysine results in a mutant that is a more effective catalyst toward monoanionic substrates. The surprising conversion of substrate preference lends support to the notion that mutations following gene duplication can result in an altered enzyme with different catalytic capabilities and preferences and may provide a pathway for the evolution of new enzymes.

Structural and Kinetic Properties of the Aldehyde Dehydrogenase NahF, a Broad Substrate Specificity Enzyme for Aldehyde Oxidation.

The salicylaldehyde dehydrogenase (NahF) catalyzes the oxidation of salicylaldehyde to salicylate using NAD(+) as a cofactor, the last reaction of the upper degradation pathway of naphthalene in Pseudomonas putida G7. The naphthalene is an abundant and toxic compound in oil and has been used as a model for bioremediation studies. The steady-state kinetic parameters for oxidation of aliphatic or aromatic aldehydes catalyzed by 6xHis-NahF are presented. The 6xHis-NahF catalyzes the oxidation of aromatic aldehydes with large kcat/Km values close to 10(6) M(-1) s(-1). The active site of NahF is highly hydrophobic, and the enzyme shows higher specificity for less polar substrates than for polar substrates, e.g., acetaldehyde. The enzyme shows α/ß folding with three well-defined domains: the oligomerization domain, which is responsible for the interlacement between the two monomers; the Rossmann-like fold domain, essential for nucleotide binding; and the catalytic domain. A salicylaldehyde molecule was observed in a deep pocket in the crystal structure of NahF where the catalytic C284 and E250 are present. Moreover, the residues G150, R157, W96, F99, F274, F279, and Y446 were thought to be important for catalysis and specificity for aromatic aldehydes. Understanding the molecular features responsible for NahF activity allows for comparisons with other aldehyde dehydrogenases and, together with structural information, provides the information needed for future mutational studies aimed to enhance its stability and specificity and further its use in biotechnological processes.

Mechanistic Aspects of Phosphate Diester Cleavage Assisted by Imidazole. A Template Reaction for Obtaining Aryl Phosphoimidazoles.

Phosphoimidazole-containing compounds are versatile players in biological and chemical processes. We explore catalytic and mechanistic criteria for the efficient formation of cyclic aryl phosphoimidazoles in aqueous solution, viewed as a template reaction for the in situ synthesis of related compounds. To provide a detailed analysis for this reaction a series of o-(2'-imidazolyl)naphthyl (4-nitrophenyl) phosphate isomers were examined to provide a basis for analysis of both mechanism and the influence of structural factors affecting the nucleophilic attack of the imidazolyl group on the phosphorus center of the substrate. Formation of the cyclic aryl phosphoimidazoles was probed by NMR and ESI-MS techniques. Kinetic experiments show that cyclization is faster under alkaline conditions, with an effective molarity up to 2900 M for the imidazolyl group, ruling out competition from external nucleophiles. Heavy atom isotope effect and computational studies show that the reaction occurs through a SN2(P)-type mechanism involving a pentacoordinated phosphorus TS, with apical positions occupied by the incoming imidazolyl nucleophile and the p-nitrophenolate leaving group. The P-O bond to the leaving group is about 50-60% broken in the transition state.

Conservative tryptophan mutants of the protein tyrosine phosphatase YopH exhibit impaired WPD-loop function and crystallize with divanadate esters in their active sites.

Catalysis in protein tyrosine phosphatases (PTPs) involves movement of a protein loop called the WPD loop that brings a conserved aspartic acid into the active site to function as a general acid. Mutation of the tryptophan in the WPD loop of the PTP YopH to any other residue with a planar, aromatic side chain (phenylalanine, tyrosine, or histidine) disables general acid catalysis. Crystal structures reveal these conservative mutations leave this critical loop in a catalytically unproductive, quasi-open position. Although the loop positions in crystal structures are similar for all three conservative mutants, the reasons inhibiting normal loop closure differ for each mutant. In the W354F and W354Y mutants, steric clashes result from six-membered rings occupying the position of the five-membered ring of the native indole side chain. The histidine mutant dysfunction results from new hydrogen bonds stabilizing the unproductive position. The results demonstrate how even modest modifications can disrupt catalytically important protein dynamics. Crystallization of all the catalytically compromised mutants in the presence of vanadate gave rise to vanadate dimers at the active site. In W354Y and W354H, a divanadate ester with glycerol is observed. Such species have precedence in solution and are known from the small molecule crystal database. Such species have not been observed in the active site of a phosphatase, as a functional phosphatase would rapidly catalyze their decomposition. The compromised functionality of the mutants allows the trapping of species that undoubtedly form in solution and are capable of binding at the active sites of PTPs, and, presumably, other phosphatases. In addition to monomeric vanadate, such higher-order vanadium-based molecules are likely involved in the interaction of vanadate with PTPs in solution.

Kinetic isotope effects in the characterization of catalysis by protein tyrosine phosphatases.

Although thermodynamically favorable, the uncatalyzed hydrolysis of phosphate monoesters is extraordinarily slow, making phosphatases among the most catalytically efficient enzymes known. Protein-tyrosine phosphatases (PTPs) are ubiquitous in biology, and kinetic isotope effects were one of the key mechanistic tools used to discern molecular details of their catalytic mechanism and the transition state for phosphoryl transfer. Later, the unique level of detail KIEs provided led to deeper questions about the potential role of protein motions in PTP catalysis. The recent discovery that such motions are responsible for different catalytic rates between PTPs arose from questions originating from KIE data showing that the transition states and chemical mechanisms are identical, combined with structural data demonstrating superimposable active sites. KIEs also reveal perturbations to the transition state as mutations are made to residues directly involved in chemistry, and to residues that affect protein motions essential for catalysis. This article is part of a Special Issue entitled: Enzyme Transition States from Theory and Experiment.

New functional aspects of the atypical protein tyrosine phosphatase VHZ.

LDP3 (VHZ) is the smallest classical protein tyrosine phosphatase (PTP) known to date and was originally misclassified as an atypical dual-specificity phosphatase. Kinetic isotope effects with steady-state and pre-steady-state kinetics of VHZ and mutants with p-nitrophenol phosphate have revealed several unusual properties. VHZ is significantly more active than previously reported but remains one of the least active PTPs. Highly unusual for a PTP, VHZ possesses two acidic residues (E134 and D65) in the active site. D65 occupies the position corresponding to the typical general acid in the PTP family. However, VHZ primarily utilizes E134 as the general acid, with D65 taking over this role when E134 is mutated. This unusual behavior is facilitated by two coexisting, but unequally populated, substrate binding modes. Unlike most classical PTPs, VHZ exhibits phosphotransferase activity. Despite the presence of the Q-loop that normally prevents alcoholysis of the phosphoenzyme intermediate in other classical PTPs, VHZ readily phosphorylates ethylene glycol. Although mutations of Q-loop residues affect this phosphotransferase activity, mutations on the IPD loop that contains the general acid exert more control over this process. A single P68V substitution on this loop completely abolishes phosphotransferase activity. The ability of native VHZ to catalyze transphosphorylation may lead to an imbalance of intracellular phosphorylation, which could explain the correlation of its overexpression with several types of cancer.

Conformational motions regulate phosphoryl transfer in related protein tyrosine phosphatases.

Many studies have implicated a role for conformational motions during the catalytic cycle, acting to optimize the binding pocket or facilitate product release, but a more intimate role in the chemical reaction has not been described. We address this by monitoring active-site loop motion in two protein tyrosine phosphatases (PTPs) using nuclear magnetic resonance spectroscopy. The PTPs, YopH and PTP1B, have very different catalytic rates; however, we find in both that the active-site loop closes to its catalytically competent position at rates that mirror the phosphotyrosine cleavage kinetics. This loop contains the catalytic acid, suggesting that loop closure occurs concomitantly with the protonation of the leaving group tyrosine and explains the different kinetics of two otherwise chemically and mechanistically indistinguishable enzymes.