The molecular details of WPD-loop movement differ in the protein-tyrosine phosphatases YopH and PTP1B

Native dynamics and allosteric responses in PTP1B probed by high-resolution HDX-MS

Protein tyrosine phosphatase 1B (PTP1B) is a validated therapeutic target for obesity, diabetes, and certain types of cancer. In particular, allosteric inhibitors hold potential for therapeutic use, but an incomplete understanding of conformational dynamics and allostery in this protein has hindered their development. Here, we interrogate solution dynamics and allosteric responses in PTP1B using high-resolution hydrogen-deuterium exchange mass spectrometry (HDX-MS), an emerging and powerful biophysical technique. Using HDX-MS, we obtain a detailed map of backbone amide exchange that serves as a proxy for the solution dynamics of apo PTP1B, revealing several flexible loops interspersed among more constrained and rigid regions within the protein structure, as well as local regions that exchange faster than expected from their secondary structure and solvent accessibility. We demonstrate that our HDX rate data obtained in solution adds value to estimates of conformational heterogeneity derived from a pseudo-ensemble constructed from ~200 crystal structures of PTP1B. Furthermore, we report HDX-MS maps for PTP1B with active-site versus allosteric small-molecule inhibitors. These maps suggest distinct and widespread effects on protein dynamics relative to the apo form, including changes in locations distal (>35 Å) from the respective ligand binding sites. These results illuminate that allosteric inhibitors of PTP1B can induce unexpected changes in dynamics that extend beyond the previously understood allosteric network. Together, our data suggest a model of BB3 allostery in PTP1B that combines conformational restriction of active-site residues with compensatory liberation of distal residues that aid in entropic balancing. Overall, our work showcases the potential of HDX-MS for elucidating aspects of protein conformational dynamics and allosteric effects of small-molecule ligands and highlights the potential of integrating HDX-MS alongside other complementary methods, such as room-temperature X-ray crystallography, NMR spectroscopy, and molecular dynamics simulations, to guide the development of new therapeutics.

Allostery in Protein Tyrosine Phosphatases is Enabled by Divergent Dynamics

Dynamics-driven allostery provides important insights into the working mechanics of proteins, especially enzymes. In this study, we employ this paradigm to answer a basic question: in enzyme superfamilies, where the catalytic mechanism, active sites, and protein fold are conserved, what accounts for the difference in the catalytic prowess of the individual members? We show that when subtle changes in sequence do not translate to changes in structure, they do translate to changes in dynamics. We use sequentially diverse PTP1B, TbPTP1, and YopH as representatives of the conserved protein tyrosine phosphatase (PTP) superfamily. Using amino acid network analysis of group behavior (community analysis) and influential node dominance on networks (eigenvector centrality), we explain the dynamic basis of the catalytic variations seen between the three proteins. Importantly, we explain how a dynamics-based blueprint makes PTP1B amenable to allosteric control and how the same is abstracted in TbPTP1 and YopH.

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.

Cancer and Autism: How PTEN Mutations Degrade Function at the Membrane and Isoform Expression in the Human Brain

Mutations causing loss of PTEN lipid phosphatase activity can promote cancer, benign tumors (PHTS), and neurodevelopmental disorders (NDDs). Exactly how they preferentially trigger distinct phenotypic outcomes has been puzzling. Here, we demonstrate that PTEN mutations differentially allosterically bias P loop dynamics and its connection to the catalytic site, affecting catalytic activity. NDD-related mutations are likely to sample conformations of the functional wild-type state, while sampled conformations for the strong, cancer-related driver mutation hotspots favor catalysis-primed conformations, suggesting that NDD mutations are likely to be weaker, and our large-scale simulations show why. Prenatal PTEN isoform expression data suggest exons 5 and 7, which harbor NDD mutations, as cancer-risk carriers. Since cancer requires more than a single mutation, our conformational and genomic analysis helps discover how same protein mutations can foster different clinical manifestations, articulates a role for co-occurring background latent driver mutations, and uncovers relationships of splicing isoform expression to life expectancy.

Native dynamics and allosteric responses in PTP1B probed by high-resolution HDX-MS

Protein tyrosine phosphatase 1B (PTP1B) is a validated therapeutic target for obesity, diabetes, and certain types of cancer. In particular, allosteric inhibitors hold potential for therapeutic use, but an incomplete understanding of conformational dynamics and allostery in this protein has hindered their development. Here, we interrogate solution dynamics and allosteric responses in PTP1B using high-resolution hydrogen-deuterium exchange mass spectrometry (HDX-MS), an emerging and powerful biophysical technique. Using HDX-MS, we obtain a detailed map of the solution dynamics of apo PTP1B, revealing several flexible loops interspersed among more constrained and rigid regions within the protein structure, as well as local regions that exchange faster than expected from their secondary structure and buriedness. We demonstrate that our HDX rate data obtained in solution adds value to predictions of dynamics derived from a pseudo-ensemble constructed from ~200 crystal structures of PTP1B. Furthermore, we report HDX-MS maps for PTP1B with active-site vs. allosteric small-molecule inhibitors. These maps reveal distinct, dramatic, and widespread effects on protein dynamics relative to the apo form, including changes to dynamics in locations distal (>35 Å) from the respective ligand binding sites. These results help shed light on the allosteric nature of PTP1B and the surprisingly far-reaching consequences of inhibitor binding in this important protein. Overall, our work showcases the potential of HDX-MS for elucidating protein conformational dynamics and allosteric effects of small-molecule ligands, and highlights the potential of integrating HDX-MS alongside other complementary methods to guide the development of new therapeutics.

Allostery in Protein Tyrosine Phosphatases is Enabled by Divergent Dynamics

Dynamics-driven allostery provides important insights into the working mechanics of proteins, especially enzymes. In this study we employ this paradigm to answer a basic question: in enzyme superfamilies where the catalytic mechanism, active sites and protein fold are conserved, what accounts for the difference in the catalytic prowess of the individual members? We show that when subtle changes in sequence do not translate to changes in structure, they do translate to changes in dynamics. We use sequentially diverse PTP1B, TbPTP1, and YopH as the representatives of the conserved Protein Tyrosine Phosphatase (PTP) superfamily. Using amino acid network analysis of group behavior (community analysis) and influential node dominance on networks (eigenvector centrality), we explain the dynamic basis of catalytic variations seen between the three proteins. Importantly, we explain how a dynamics-based blueprint makes PTP1B amenable to allosteric control and how the same is abstracted in TbPTP1 and YopH.

Q-RepEx: A Python pipeline to increase the sampling of empirical valence bond simulations

The exploration of chemical systems occurs on complex energy landscapes. Comprehensively sampling rugged energy landscapes with many local minima is a common problem for molecular dynamics simulations. These multiple local minima trap the dynamic system, preventing efficient sampling. This is a particular challenge for large biochemical systems with many degrees of freedom. Replica exchange molecular dynamics (REMD) is an approach that accelerates the exploration of the conformational space of a system, and thus can be used to enhance the sampling of complex biomolecular processes. In parallel, the empirical valence bond (EVB) approach is a powerful approach for modeling chemical reactivity in biomolecular systems. Here, we present an open-source Python-based tool that interfaces with the Q simulation package, and increases the sampling efficiency of the EVB free energy perturbation/umbrella sampling approach by means of REMD. This approach, Q-RepEx, both decreases the computational cost of the associated REMD-EVB simulations, and opens the door to more efficient studies of biochemical reactivity in systems with significant conformational fluctuations along the chemical reaction coordinate.

How PTEN mutations degrade function at the membrane and life expectancy of carriers of mutations in the human brain

PTEN dysfunction, caused by loss of lipid phosphatase activity or deletion, promotes pathologies, cancer, benign tumors, and neurodevelopmental disorders (NDDs). Despite efforts, exactly how the mutations trigger distinct phenotypic outcomes, cancer or NDD, has been puzzling. It has also been unclear how to distinguish between mutations harbored by isoforms, are they cancer or NDDs-related. Here we address both. We demonstrate that PTEN mutations differentially allosterically bias P-loop dynamics and its connection to the catalytic site, affecting catalytic activity. NDD-related mutations are likely to sample conformations present in the wild-type, while sampled conformations sheltering cancer-related hotspots favor catalysis-prone conformations, suggesting that NDD mutations are weaker. Analysis of isoform expression data indicates that if the transcript has NDD-related mutations, alone or in combination with cancer hotspots, there is high prenatal expression. If no mutations within the measured days, low expression levels. Cancer mutations promote stronger signaling and cell proliferation; NDDs' are weaker, influencing brain cell differentiation. Further, exon 5 is impacted by NDD or non-NDD mutations, while exon 7 is exclusively impacted by NDD mutations. Our comprehensive conformational and genomic analysis helps discover how same allele mutations can foster different clinical manifestations and uncovers correlations of splicing isoform expression to life expectancy.

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.

Drug Discovery of Plausible Lead Natural Compounds That Target the Insulin Signaling Pathway: Bioinformatics Approaches

The growing smooth talk in the field of natural compounds is due to the ancient and current interest in herbal medicine and their potentially positive effects on health. Dozens of antidiabetic natural compounds were reported and tested in vivo, in silico, and in vitro. The role of these natural compounds, their actions on the insulin signaling pathway, and the stimulation of the glucose transporter-4 (GLUT4) insulin-responsive translocation to the plasma membrane (PM) are all crucial in the treatment of diabetes and insulin resistance. In this review, we collected and summarized a group of available in vivo and in vitro studies which targeted isolated phytochemicals with possible antidiabetic activity. Moreover, the in silico docking of natural compounds with some of the insulin signaling cascade key proteins is also summarized based on the current literature. In this review, hundreds of recent studies on pure natural compounds that alleviate type II diabetes mellitus (type II DM) were revised. We focused on natural compounds that could potentially regulate blood glucose and stimulate GLUT4 translocation through the phosphoinositide 3-kinase (PI3K)/protein kinase B (Akt) pathway. On attempt to point out potential new natural antidiabetic compounds, this review also focuses on natural ingredients that were shown to interact with proteins in the insulin signaling pathway in silico, regardless of their in vitro/in vivo antidiabetic activity. We invite interested researchers to test these compounds as potential novel type II DM drugs and explore their therapeutic mechanisms.

Crystal structure of the Cys-NO modified YopH tyrosine phosphatase

Protein tyrosine phosphatases (PTPs) are key virulence factors in pathogenic bacteria, consequently, they have become important targets for new approaches against these pathogens, especially in the fight against antibiotic resistance. Among these targets of interest YopH (Yersinia outer protein H) from virulent species of Yersinia is an example. PTPs can be reversibly inhibited by nitric oxide (NO) since the oxidative modification of cysteine residues may influence the protein structure and catalytic activity. We therefore investigated the effects of NO on the structure and enzymatic activity of Yersinia enterocolitica YopH in vitro. Through phosphatase activity assays, we observe that in the presence of NO YopH activity was inhibited by 50%, and that this oxidative modification is partially reversible in the presence of DTT. Furthermore, YopH S-nitrosylation was clearly confirmed by a biotin switch assay, high resolution mass spectrometry (MS) and X-ray crystallography approaches. The crystal structure confirmed the S-nitrosylation of the catalytic cysteine residue, Cys403, while the MS data provide evidence that Cys221 and Cys234 might also be modified by NO. Interestingly, circular dichroism spectroscopy shows that the S-nitrosylation affects secondary structure of wild type YopH, though to a lesser extent on the catalytic cysteine to serine YopH mutant. The data obtained demonstrate that S-nitrosylation inhibits the catalytic activity of YopH, with effects beyond the catalytic cysteine. These findings are helpful for designing effective YopH inhibitors and potential therapeutic strategies to fight this pathogen or others that use similar mechanisms to interfere in the signal transduction pathways of their hosts.

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.

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.

Are vanadium complexes druggable against the main protease Mpro of SARS-CoV-2? - A computational approach

In silico techniques helped explore the binding capacities of the SARS-CoV-2 main protease (M^pro) for a series of metalloorganic compounds. Along with small size vanadium complexes a vanadium-containing derivative of the peptide-like inhibitor N3 (N-[(5-methylisoxazol-3-yl)carbonyl]alanyl-l-valyl-N1-((1R,2Z)-4-(benzyloxy)-4-oxo-1-{[(3R)-2-oxopyrrolidin-3-yl] methyl }but-2-enyl)-l-leucinamide) was designed from the crystal structure with PDB entry code 6LU7. On theoretical grounds our consensus docking studies evaluated the binding affinities at the hitherto known binding site of Chymotrypsin-like protease (3CLpro) of SARS-CoV-2 for existing and designed vanadium complexes. This main virus protease (M^pro) has a Cys-His dyad at the catalytic site that is characteristic of metal-dependent or metal-inhibited hydrolases. M^pro was compared to the human protein-tyrosine phosphatase 1B (hPTP1B) with a comparable catalytic dyad. HPTP1B is a key regulator at an early stage in the signalling cascade of the insulin hormone for glucose uptake into cells. The vanadium-ligand binding site of hPTP1B is located in a larger groove on the surface of M^pro. Vanadium constitutes a well-known phosphate analogue. Hence, its study offers possibilities to design promising vanadium-containing binders to SARS-CoV-2. Given the favourable physicochemical properties of vanadium nuclei, such organic vanadium complexes could become drugs not only for pharmacotherapy but also diagnostic tools for early infection detection in patients. This work presents the in silico design of a potential lead vanadium compound. It was tested along with 20 other vanadium-containing complexes from the literature in a virtual screening test by docking to inhibit M^pro of SARS-CoV-2.

The mechanism of full activation of tumor suppressor PTEN at the phosphoinositide-enriched membrane

Tumor suppressor PTEN, the second most highly mutated protein in cancer, dephosphorylates signaling lipid PIP₃ produced by PI3Ks. Excess PIP₃ promotes cell proliferation. The mechanism at the membrane of this pivotal phosphatase is unknown hindering drug discovery. Exploiting explicit solvent simulations, we tracked full-length PTEN trafficking from the cytosol to the membrane. We observed its interaction with membranes composed of zwitterionic phosphatidylcholine, anionic phosphatidylserine, and phosphoinositides, including signaling lipids PIP₂ and PIP₃. We tracked its moving away from the zwitterionic and getting absorbed onto anionic membrane that harbors PIP₃. We followed it localizing on microdomains enriched in signaling lipids, as PI3K does, and observed PIP₃ allosterically unfolding the N-terminal PIP₂ binding domain, positioning it favorably for the polybasic motif interaction with PIP₂. Finally, we determined PTEN catalytic action at the membrane, all in line with experimental observations, deciphering the mechanisms of how PTEN anchors to the membrane and restrains cancer.

•

PTEN localizes on membrane microdomains enriched in phosphoinositides, as PI3K does

•

Full PTEN activation requires both signaling lipids, PIP₂ and PIP₃

•

Strong salt bridge interactions sustain stable PTEN membrane localization

•

Substrate-induced P loop conformational change implicates PTEN catalytic activity

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.

Ligand-Binding-Site Refinement to Generate Reliable Holo Protein Structure Conformations from Apo Structures

The first important step in a structure-based virtual screening is the judicious selection of a receptor protein. In cases where the holo protein receptor structure is unavailable, significant reduction in virtual screening performance has been reported. In this work, we present a robust method to generate reliable holo protein structure conformations from apo structures using molecular dynamics (MD) simulation with restraints derived from holo structure binding-site templates. We perform benchmark tests on two different datasets: 40 structures from a directory of useful decoy-enhanced (DUD-E) and 84 structures from the Gunasekaran dataset. Our results show successful refinement of apo binding-site structures toward holo conformations in 82% of the test cases. In addition, virtual screening performance of 40 DUD-E structures is significantly improved using our MD-refined structures as receptors with an average enrichment factor (EF), an EF_1% value of 6.2 compared to apo structures with 3.5. Docking of native ligands to the refined structures shows an average ligand root mean square deviation (RMSD) of 1.97 Å (DUD-E dataset and Gunasekaran dataset) relative to ligands in the holo crystal structures, which is comparable to the self-docking (i.e., docking of the native ligand back to its crystal structure receptor) average, 1.34 Å (DUD-E dataset) and 1.36 Å (Gunasekaran dataset). On the other hand, docking to the apo structures yields an average ligand RMSD of 3.65 Å (DUD-E) and 2.90 Å (Gunasekaran). These results indicate that our method is robust and can be useful to improve virtual screening performance of apo structures.

Why and how to investigate the role of protein phosphorylation in ZIP and ZnT zinc transporter activity and regulation

Zinc is required for the regulation of proliferation, metabolism, and cell signaling. It is an intracellular second messenger, and the cellular level of ionic, mobile zinc is strictly controlled by zinc transporters. In mammals, zinc homeostasis is primarily regulated by ZIP and ZnT zinc transporters. The importance of these transporters is underscored by the list of diseases resulting from changes in transporter expression and activity. However, despite numerous structural studies of the transporters revealing both zinc binding sites and motifs important for transporter function, the exact molecular mechanisms regulating ZIP and ZnT activities are still not clear. For example, protein phosphorylation was found to regulate ZIP7 activity resulting in the release of Zn2+ from intracellular stores leading to phosphorylation of tyrosine kinases and activation of signaling pathways. In addition, sequence analyses predict all 24 human zinc transporters to be phosphorylated suggesting that protein phosphorylation is important for regulation of transporter function. This review describes how zinc transporters are implicated in a number of important human diseases. It summarizes the current knowledge regarding ZIP and ZnT transporter structures and points to how protein phosphorylation seems to be important for the regulation of zinc transporter activity. The review addresses the need to investigate the role of protein phosphorylation in zinc transporter function and regulation, and argues for a pressing need to introduce quantitative phosphoproteomics to specifically target zinc transporters and proteins involved in zinc signaling. Finally, different quantitative phosphoproteomic strategies are suggested.

Allosteric Impact of the Variable Insert Loop in Vaccinia H1-Related (VHR) Phosphatase

Dynamics and conformational motions are important to the activity of enzymes, including protein tyrosine phosphatases. These motions often extend to regions outside the active site, called allosteric regions. In the tyrosine phosphatase Vaccinia H1-related (VHR) enzyme, we demonstrate the importance of the allosteric interaction between the variable insert region and the active-site loops in VHR. These studies include solution nuclear magnetic resonance, computation, steady-state, and rapid kinetic measurements. Overall, the data indicate concerted millisecond motions exist between the variable insert and the catalytic acid loop in wild-type (WT) VHR. The 150 ns computation studies show a flexible acid loop in WT VHR that opens during the simulation from its initial closed structure. Mutation of the variable insert residue, asparagine 74, to alanine results in a rigidification of the acid loop as observed by molecular dynamics simulations and a disruption of crucial active-site hydrogen bonds. Moreover, enzyme kinetic analysis shows a weakening of substrate affinity in the N74A mutant and a >2-fold decrease in substrate cleavage and hydrolysis rates. These data show that despite being nearly 20 Å from the active site, the variable insert region is linked to the acid loop by coupled millisecond motions, and that disruption of the communication between the variable insert and active site alters the normal catalytic function of VHR and perturbs the active-site environment.

Uncovering the Molecular Interactions in the Catalytic Loop That Modulate the Conformational Dynamics in Protein Tyrosine Phosphatase 1B

Active-site loops are integral to the function of numerous enzymes. They enable substrate and product binding and release, sequester reaction intermediates, and recruit catalytic groups. Here, we examine the catalytic loop in the enzyme protein tyrosine phosphatase 1B (PTP1B). PTP1B has a mobile so-called WPD loop (named for its three N-terminal residues) that initiates the dephosphorylation of phosphor-tyrosine substrates upon loop closure. We have combined X-ray crystallography, solution NMR, and pre-steady-state kinetics experiments on wild-type and five WPD loop mutants to identify the relationships between the loop structure, dynamics, and function. The motions of the WPD loop are modulated by the formation of weak molecular interactions, where perturbations of these interactions modulate the conformational equilibrium landscape. The point mutants in the WPD loop alter the loop equilibrium position from a predominantly open state (P185A) to 50:50 (F182A), 35:65 (P188A), and predominantly closed states (T177A and P188A). Surprisingly, there is no correlation between the observed catalytic rates in the loop mutants and changes to the WPD loop equilibrium position. Rather, we observe a strong correlation between the rate of dephosphorylation of the phosphocysteine enzyme intermediate and uniform millisecond motions, not only within the loop but also in the adjacent α-helical domain of PTP1B. Thus, the control of loop motion and thereby catalytic activity is dispersed and resides within not only the loop sequence but also the surrounding protein architecture. This has broad implications for the general mechanistic understanding of enzyme reactions and the role that flexible loops play in the catalytic cycle.

Vanadium Compounds as Enzyme Inhibitors with a Focus on Anticancer Effects

Vanadium salts and coordination compounds have desirable cellular anticancer effects, and although they have been investigated in detail as a potential treatment for diabetes, less attention has been given to the anticancer effects. The inhibition of some signal transduction enzymes is known, and studies of the metabolism and activation pathways both in vitro and in vivo are important for future investigations and development of vanadium's role as a new potential drug. In addition, a new approach has demonstrated that the enhancement of oncolytic viruses using vanadium salts and coordination complexes for immunotherapy is very promising. Some differences exist between this approach and current antidiabetic and anticancer studies because vanadium(iv) complexes have been found to be most potent in the latter approach, but the few compounds investigated with oncolytic viruses show that vanadium(v) systems are more effective. We conclude that recent studies demonstrate effects on signal transduction enzymes and anticancer pathways, thus suggesting potential applications of vanadium as anticancer agents in the future both as standalone treatments as well as combination therapies.

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.

Coordination environment changes of the vanadium in vanadium-dependent haloperoxidase enzymes

Vanadium-dependent haloperoxidases are a class of enzymes that catalyze oxidation reactions with halides to form halogenated organic products and water. These enzymes include chloroperoxidase and bromoperoxidase, which have very different protein sequences and sizes, but regardless the coordination environment of the active sites is surprisingly constant. In this manuscript, the comparison of the coordination chemistry of V-containing-haloperoxidases of the trigonal bipyramidal geometry was done by data mining. The catalytic cycle imposes changes in the coordination geometry of the vanadium to accommodate the peroxidovanadium(V) intermediate in an environment we describe as a distorted square pyramidal geometry. During the catalytic cycle, this intermediate converts to a trigonal bipyramidal intermediate before losing the halogen and forming a tetrahedral vanadium-protein intermediate. Importantly, the catalysis is facilitated by a proton-relay system supplied by the second sphere coordination environment and the changes in the coordination environment of the vanadium(V) making this process unique among protein catalyzed processes. The analysis of the coordination chemistry shows that the active site is very tightly regulated with only minor changes in the coordination geometry. The coordination geometry in the protein structures deviates from that found for both small molecules crystalized in the absence of protein and the reported functional small molecule model compounds. At this time there are no examples reported of a structurally similar small molecule with the geometry observed for the peroxidovanadium(V) in the active site of the vanadium-containing haloperoxidases.

Structural analysis of protein tyrosine phosphatase 1B reveals potentially druggable allosteric binding sites

Catalytic proteins such as human protein tyrosine phosphatase 1B (PTP1B), with conserved and highly polar active sites, warrant the discovery of druggable nonactive sites, such as allosteric sites, and potentially, therapeutic small molecules that can bind to these sites. Catalyzing the dephosphorylation of numerous substrates, PTP1B is physiologically important in intracellular signal transduction pathways in diverse cell types and tissues. Aberrant PTP1B is associated with obesity, diabetes, cancers, and neurodegenerative disorders. Utilizing clustering methods (based on root mean square deviation, principal component analysis, nonnegative matrix factorization, and independent component analysis), we have examined multiple PTP1B structures. Using the resulting representative structures in different conformational states, we determined consensus clustroids and used them to identify both known and novel binding sites, some of which are potentially allosteric. We report several lead compounds that could potentially bind to the novel PTP1B binding sites and can be further optimized. Considering the possibility for drug repurposing, we discovered homologous binding sites in other proteins, with ligands that could potentially bind to the novel PTP1B binding sites.

Conformational Rigidity and Protein Dynamics at Distinct Timescales Regulate PTP1B Activity and Allostery

Protein function originates from a cooperation of structural rigidity, dynamics at different timescales, and allostery. However, how these three pillars of protein function are integrated is still only poorly understood. Here we show how these pillars are connected in Protein Tyrosine Phosphatase 1B (PTP1B), a drug target for diabetes and cancer that catalyzes the dephosphorylation of numerous substrates in essential signaling pathways. By combining new experimental and computational data on WT-PTP1B and ≥10 PTP1B variants in multiple states, we discovered a fundamental and evolutionarily conserved CH/π switch that is critical for positioning the catalytically important WPD loop. Furthermore, our data show that PTP1B uses conformational and dynamic allostery to regulate its activity. This shows that both conformational rigidity and dynamics are essential for controlling protein activity. This connection between rigidity and dynamics at different timescales is likely a hallmark of all enzyme function.

Leveraging Reciprocity to Identify and Characterize Unknown Allosteric Sites in Protein Tyrosine Phosphatases

Drug-like molecules targeting allosteric sites in proteins are of great therapeutic interest; however, identification of potential sites is not trivial. A straightforward approach to identify hidden allosteric sites is demonstrated in protein tyrosine phosphatases (PTP) by creation of single alanine mutations in the catalytic acid loop of PTP1B and VHR. This approach relies on the reciprocal interactions between an allosteric site and its coupled orthosteric site. The resulting NMR chemical shift perturbations (CSPs) of each mutant reveal clusters of distal residues affected by acid loop mutation. In PTP1B and VHR, two new allosteric clusters were identified in each enzyme. Mutations in these allosteric clusters altered phosphatase activity with changes in k_cat/K_M ranging from 30% to nearly 100-fold. This work outlines a simple method for identification of new allosteric sites in PTP, and given the basis of this method in thermodynamics, it is expected to be generally useful in other systems.

Recent advances in the development of protein tyrosine phosphatase 1B inhibitors for Type 2 diabetes

Diabetes mellitus is the most serious and prevalent metabolic disorders worldwide, complications of which can decrease significantly the quality of life and contribute to premature death. Resistance to insulin is a predominant pathophysiological factor of Type 2 diabetes (T2D). Protein tyrosine phosphatase 1B (PTP1B) is an important negative factor of insulin signal and a potent therapeutic target in T2D patients. This review highlights recent advances (2012-2015) in research related to the role of PTP1B in signal transduction processes implicated in pathophysiology of T2D, and novel PTP1B inhibitors with an emphasis on their chemical structures and modes of action.

A computational study on role of 6-(hydroxymethyl)-3-[3,4,5-trihydroxy-6-[(3,4,5-trihydroxyoxan-2-yl)oxymethyl]oxan-2-yl]oxyoxane-2,4,5-triol in the regulation of blood glucose level

6-(hydroxymethyl)-3-[3,4,5-trihydroxy-6-[(3,4,5-trihydroxyoxan-2-yl)oxymethyl]oxan-2-yl]oxyoxane-2,4,5-triol (SID 242078875) was isolated from the fruits of Syzygium densiflorum Wall. ex Wight & Arn (Myrtaceae), which has been traditionally used in the treatment of diabetes by the tribes of The Nilgiris, Tamil Nadu, India. In this study, reverse pharmacophore mapping approach and text-based database search identified the dipeptidyl peptidase-IV, protein-tyrosine phosphatase 1B, phosphoenolpyruvate carboxykinase, glycogen synthase kinase-3β and glucokinase as potential targets of SID 242078875 in diabetes management. Further, molecular docking was performed to predict the binding pose of SID 242078875 in the active site region of the target protein. In addition, dynamic behaviour and stability of protein-ligand complexes were observed for a period of 50 ns through molecular dynamics simulation.

Antidiabetic, Chemical, and Physical Properties of Organic Vanadates as Presumed Transition-State Inhibitors for Phosphatases

Studies of antidiabetic vanadium compounds, specifically the organic vanadate esters, are reviewed with regard to their chemistry and biological properties. The compounds are described from the perspective of how the fundamental chemistry and properties of organic vanadate esters impact their effects as inhibitors for phosphatases based on the structural information obtained from vanadium-phosphatase complexes. Vanadium compounds have been reported to have antidiabetic properties for more than a century. The structures and properties of organic vanadate complexes are reviewed, and the potency of such vanadium coordination complexes as antidiabetic agents is described. Because such compounds form spontaneously in aqueous environments, the reactions with most components in any assay or cellular environment has potential to be important and should be considered. Generally, the active form of vanadium remains elusive, although studies have been reported of a number of promising vanadium compounds. The description of the antidiabetic properties of vanadium compounds is described here in the context of recent characterization of vanadate-phosphatase protein structures by data mining. Organic vanadate ester compounds are generally four coordinate or five coordinate with the former being substrate analogues and the latter being transition-state analogue inhibitors. These studies demonstrated a framework for characterization of five-coordinate trigonal bipyramidal vanadium inhibitors by comparison with the reported vanadium-protein phosphatase complexes. The binding of the vanadium to the phosphatases is either as a five-coordinate exploded transition-state analogue or as a high energy intermediate, respectively. Even if potency as an inhibitor requires trigonal bipyramidal geometry of the vanadium when bound to the protein, such geometry can be achieved upon binding from compounds with other geometries. Desirable properties of ligands are identified and analyzed. Ligand interactions, as reported in one peptidic substrate, are favorable so that complementarity between phosphatase and coordinating ligand to the vanadium can be established resulting in a dramatic enhancement of the inhibitory potency. These considerations point to a frameshift in ligand design for vanadium complexes as phosphatase inhibitors and are consistent with other small molecule having much lower affinities. Combined, these studies do suggest that if effective delivery of potentially active antidiabetic compound such a the organic vanadate peptidic substrate was possible the toxicity problems currently reported for the salts and some of the complexes may be alleviated and dramatic enhancement of antidiabetic vanadium compounds may result.

The intrinsically disordered tails of PTEN and PTEN-L have distinct roles in regulating substrate specificity and membrane activity

Phosphatase and tensin homologue deleted on chromosome 10 (PTEN) is a lipid and protein phosphatase, and both activities are necessary for its role as a tumour suppressor. PTEN activity is controlled by phosphorylation of its intrinsically disordered C-terminal tail. A recently discovered variant of PTEN, PTEN-long (PTEN-L), has a 173-residue N-terminal extension that causes PTEN-L to exhibit unique behaviour, such as movement from one cell to another. Using hydrogen/deuterium exchange mass spectrometry (HDX-MS) and biophysical assays, we show that both the N-terminal extension of PTEN-L and C-terminal tail of PTEN affect the phosphatase activity using unique mechanisms. Phosphorylation of six residues in the C-terminal tail of PTEN results in auto-inhibitory interactions with the phosphatase and C2 domains, effectively blocking both the active site and the membrane-binding interface of PTEN. Partially dephosphorylating PTEN on pThr(366)/pSer(370) results in sufficient exposure of the active site to allow a selective activation for soluble substrates. Using HDX-MS, we identified a membrane-binding element in the N-terminal extension of PTEN-L, termed the membrane-binding helix (MBH). The MBH radically alters the membrane binding mechanism of PTEN-L compared with PTEN, switching PTEN-L to a 'scooting' mode of catalysis from the 'hopping' mode that is characteristic of PTEN.

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.

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.

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.

New aspects of the phosphatase VHZ revealed by a high-resolution structure with vanadate and substrate screening

The recently discovered 150-residue human VHZ (VH1-related protein, Z member) is one of the smallest protein tyrosine phosphatases (PTPs) known and contains only the minimal structural elements common to all PTPs. We report a substrate screening analysis and a crystal structure of the VHZ complex with vanadate at 1.1 Å resolution, with a detailed structural comparison with other members of the protein tyrosine phosphatase family, including classical tyrosine-specific protein tyrosine phosphatases (PTPs) and dual-specificity phosphatases (DSPs). A screen with 360 phosphorylated peptides shows VHZ efficiently catalyzes the hydrolysis of phosphotyrosine (pY)-containing peptides but exhibits no activity toward phosphoserine (pS) or phosphothreonine (pT) peptides. The new structure reveals a deep and narrow active site more typical of the classical tyrosine-specific PTPs. Despite the high degrees of structural and sequence similarity between VHZ and classical PTPs, its general acid IPD-loop is most likely conformationally rigid, in contrast to the flexible WPD counterpart of classical PTPs. VHZ also lacks substrate recognition domains and other domains typically found on classical PTPs. It is therefore proposed that VHZ is more properly classified as an atypical PTP rather than an atypical DSP, as has been suggested.