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

Revealing the principles of inter- and intra-domain regulation in a signaling enzyme via scanning mutagenesis

Multi-domain enzymes can be regulated by both inter-domain interactions and structural features intrinsic to the catalytic domain. The tyrosine phosphatase SHP2 is a quintessential example of a multi-domain protein that is regulated by inter-domain interactions. This enzyme has a protein tyrosine phosphatase (PTP) domain and two phosphotyrosine-recognition domains (N-SH2 and C-SH2) that regulate phosphatase activity through autoinhibitory interactions. SHP2 is canonically activated by phosphoprotein binding to the SH2 domains, which causes large inter-domain rearrangements, but autoinhibition can also be disrupted by disease-associated mutations. Many details of the SHP2 activation mechanism are still unclear, the physiologically-relevant active conformations remain elusive, and hundreds of human variants of SHP2 have not been functionally characterized. Here, we perform deep mutational scanning on both full-length SHP2 and its isolated PTP domain to examine mutational effects on inter-domain regulation and catalytic activity. Our experiments provide a comprehensive map of SHP2 mutational sensitivity, both in the presence and absence of inter-domain regulation. Coupled with molecular dynamics simulations, our investigation reveals novel structural features that govern the stability of the autoinhibited and active states of SHP2. Our analysis also identifies key residues beyond the SHP2 active site that control PTP domain dynamics and intrinsic catalytic activity. This work expands our understanding of SHP2 regulation and provides new insights into SHP2 pathogenicity.

Motif-VI loop acts as a nucleotide valve in the West Nile Virus NS3 Helicase

The Orthoflavivirus NS3 helicase (NS3h) is crucial in virus replication, representing a potential drug target for pathogenesis. NS3h utilizes nucleotide triphosphate (ATP) for hydrolysis energy to translocate on single-stranded nucleic acids, which is an important step in the unwinding of double-stranded nucleic acids. Intermediate states along the ATP hydrolysis cycle and conformational changes between these states, represent important yet difficult-to-identify targets for potential inhibitors. Extensive molecular dynamics simulations of West Nile virus NS3h+ssRNA in the apo, ATP, ADP+Pi and ADP bound states were used to model the conformational ensembles along this cycle. Energetic and structural clustering analyses depict a clear trend of differential enthalpic affinity of NS3h with ADP, demonstrating a probable mechanism of hydrolysis turnover regulated by the motif-VI loop (MVIL). Based on these results, MVIL mutants (D471L, D471N and D471E) were found to have a substantial reduction in ATPase activity and RNA replication compared to the wild-type. Simulations of the mutants in the apo state indicate a shift in MVIL populations favoring either a closed or open 'valve' conformation, affecting ATP entry or stabilization, respectively. Combining our molecular modeling with experimental evidence highlights a conformation-dependent role for MVIL as a 'valve' for the ATP-pocket, presenting a promising target for antiviral development.

NMR and Single-Molecule FRET Insights into Fast Protein Motions and Their Relation to Function

Proteins often undergo large-scale conformational transitions, in which secondary and tertiary structure elements (loops, helices, and domains) change their structures or their positions with respect to each other. Simple considerations suggest that such dynamics should be relatively fast, but the functional cycles of many proteins are often relatively slow. Sophisticated experimental methods are starting to tackle this dichotomy and shed light on the contribution of large-scale conformational dynamics to protein function. In this review, we focus on the contribution of single-molecule Förster resonance energy transfer and nuclear magnetic resonance (NMR) spectroscopies to the study of conformational dynamics. We briefly describe the state of the art in each of these techniques and then point out their similarities and differences, as well as the relative strengths and weaknesses of each. Several case studies, in which the connection between fast conformational dynamics and slower function has been demonstrated, are then introduced and discussed. These examples include both enzymes and large protein machines, some of which have been studied by both NMR and fluorescence spectroscopies.

Effects of Xylanase A double mutation on substrate specificity and structural dynamics

While protein activity is traditionally studied with a major focus on the active site, the activity of enzymes has been hypothesized to be linked to the flexibility of adjacent regions, warranting more exploration into how the dynamics in these regions affects catalytic turnover. One such enzyme is Xylanase A (XylA), which cleaves hemicellulose xylan polymers by hydrolysis at internal β-1,4-xylosidic linkages. It contains a "thumb" region whose flexibility has been suggested to affect the activity. The double mutation D11F/R122D was previously found to affect activity and potentially bias the thumb region to a more open conformation. We find that the D11F/R122D double mutation shows substrate-dependent effects, increasing activity on the non-native substrate ONPX2 but decreasing activity on its native xylan substrate. To characterize how the double mutant causes these kinetics changes, nuclear magnetic resonance (NMR) and molecular dynamics (MD) simulations were used to probe structural and flexibility changes. NMR chemical shift perturbations revealed structural changes in the double mutant relative to the wild-type, specifically in the thumb and fingers regions. Increased slow-timescale dynamics in the fingers region was observed as intermediate-exchange line broadening. Lipari-Szabo order parameters show negligible changes in flexibility in the thumb region in the presence of the double mutation. To help understand if there is increased energetic accessibility to the open state upon mutation, alchemical free energy simulations were employed that indicated thumb opening is more favorable in the double mutant. These studies aid in further characterizing how flexibility in adjacent regions affects the function of XylA.

Protein Tyrosine Phosphatase regulation by Reactive Oxygen Species

Protein Tyrosine Phosphatases (PTPs) help to maintain the balance of protein phosphorylation signals that drive cell division, proliferation, and differentiation. These enzymes are also well-suited to redox-dependent signaling and oxidative stress response due to their cysteine-based catalytic mechanism, which requires a deprotonated thiol group at the active site. This review focuses on PTP structural characteristics, active site chemical properties, and vulnerability to change by reactive oxygen species (ROS). PTPs can be oxidized and inactivated by H2O2 through three non-exclusive mechanisms. These pathways are dependent on the coordinated actions of other H2O2-sensitive proteins, such as peroxidases like Peroxiredoxins (Prx) and Thioredoxins (Trx). PTPs undergo reversible oxidation by converting their active site cysteine from thiol to sulfenic acid. This sulfenic acid can then react with adjacent cysteines to form disulfide bonds or with nearby amides to form sulfenyl-amide linkages. Further oxidation of the sulfenic acid form to the sulfonic or sulfinic acid forms causes irreversible deactivation. Understanding the structural changes involved in both reversible and irreversible PTP oxidation can help with their chemical manipulation for therapeutic intervention. Nonetheless, more information remains unidentified than is presently known about the precise dynamics of proteins participating in oxidation events, as well as the specific oxidation states that can be targeted for PTPs. This review summarizes current information on PTP-specific oxidation patterns and explains how ROS-mediated signal transmission interacts with phosphorylation-based signaling machinery controlled by growth factor receptors and PTPs.

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.

Multiscale In Silico Study of the Mechanism of Activation of the RtcB Ligase by the PTP1B Phosphatase

Inositol-requiring enzyme 1 (IRE1) is a transmembrane sensor that is part of a trio of sensors responsible for controlling the unfolded protein response within the endoplasmic reticulum (ER). Upon the accumulation of unfolded or misfolded proteins in the ER, IRE1 becomes activated and initiates the cleavage of a 26-nucleotide intron from human X-box-containing protein 1 (XBP1). The cleavage is mediated by the RtcB ligase enzyme, which splices together two exons, resulting in the formation of the spliced isoform XBP1s. The XBP1s isoform activates the transcription of genes involved in ER-associated degradation to maintain cellular homeostasis. The catalytic activity of RtcB is controlled by the phosphorylation and dephosphorylation of three tyrosine residues (Y306, Y316, and Y475), which are regulated by the ABL1 tyrosine kinase and PTP1B phosphatase, respectively. This study focuses on investigating the mechanism by which the PTP1B phosphatase activates the RtcB ligase using a range of advanced in silico methods. Protein-protein docking identified key interacting residues between RtcB and PTP1B. Notably, the phosphorylated Tyr306 formed hydrogen bonds and salt bridge interactions with the "gatekeeper" residues Arg47 and Lys120 of the inactive PTP1B. Classical molecular dynamics simulation emphasized the crucial role of Asp181 in the activation of PTP1B, driving the conformational change from an open to a closed state of the WPD-loop. Furthermore, QM/MM-MD simulations provided insights into the free energy landscape of the dephosphorylation reaction mechanism of RtcB, which is mediated by the PTP1B phosphatase.

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.

Motif-VI Loop Acts as a Nucleotide Valve in the West Nile Virus NS3 Helicase

The flavivirus NS3 helicase (NS3h), a highly conserved protein, plays a pivotal role in virus replication and thus represents a potential drug target for flavivirus pathogenesis. NS3h utilizes nucleotide triphosphate, such as ATP, for hydrolysis energy (ATPase) to translocate on single-stranded nucleic acids, which is an important step in the unwinding of double-stranded nucleic acids. The intermediate states along the ATP binding and hydrolysis cycle, as well as the conformational changes between these states, represent important yet difficult-to-identify targets for potential inhibitors. We use extensive molecular dynamics simulations of apo, ATP, ADP+Pi, and ADP bound to WNV NS3h+ssRNA to model the conformational ensembles along this cycle. Energetic and structural clustering analyses on these trajectories depict a clear trend of differential enthalpic affinity of NS3h with ADP, demonstrating a probable mechanism of hydrolysis turnover regulated by the motif-VI loop (MVIL). These findings were experimentally corroborated using viral replicons encoding three mutations at the D471 position. Replication assays using these mutants demonstrated a substantial reduction in viral replication compared to the wild-type. Molecular simulations of the D471 mutants in the apo state indicate a shift in MVIL populations favoring either a closed or open 'valve' conformation, affecting ATP entry or stabilization, respectively. Combining our molecular modeling with experimental evidence highlights a conformation-dependent role for MVIL as a 'valve' for the ATP-pocket, presenting a promising target for antiviral development.

Enhanced Thermostability of an l-Rhamnose Isomerase for d-Allose Synthesis by Computation-Based Rational Redesign of Flexible Regions

d-Allose is a low-calorie rare sugar with great application potential in the food and pharmaceutical industries. The production of d-allose has been accomplished using l-rhamnose isomerase (L-RI), but concomitantly increasing the enzyme's stability and activity remains challenging. Here, we rationally engineered an L-RI from Clostridium stercorarium to enhance its stability by comprehensive computation-aided redesign of its flexible regions, which were successively identified using molecular dynamics simulations. The resulting combinatorial mutant M2-4 exhibited a 5.7-fold increased half-life at 75 °C while also exhibiting improved catalytic efficiency. Especially, by combining structure modeling and multiple sequence alignment, we identified an α0 region that was universal in the L-RI family and likely acted as a "helix-breaker". Truncating this region is crucial for improving the thermostability of related enzymes. Our work provides a significantly stable biocatalyst with potential for the industrial production of d-allose.

Conformational Switch of the 250s Loop Enables the Efficient Transglycosylation in GH Family 77

Amylomaltases (AMs) play important roles in glycogen and maltose metabolism. However, the molecular mechanism is elusive. Here, we investigated the conformational dynamics of the 250s loop and catalytic mechanism of Thermus aquaticus TaAM using path-metadynamics and QM/MM MD simulations. The results demonstrate that the transition of the 250s loop from an open to closed conformation promotes polysaccharide sliding, leading to the ideal positioning of the acid/base. Furthermore, the conformational dynamics can also modulate the selectivity of hydrolysis and transglycosylation. The closed conformation of the 250s loop enables the tight packing of the active site for transglycosylation, reducing the energy penalty and efficiently preventing the penetration of water into the active site. Conversely, the partially closed conformation for hydrolysis results in a loosely packed active site, destabilizing the transition state. These computational findings guide mutation experiments and enable the identification of mutants with an improved disproportionation/hydrolysis ratio. The present mechanism is in line with experimental data, highlighting the critical role of conformational dynamics in regulating the catalytic reactivity of GHs.

Interaction of Vanadium Complexes with Proteins: Revisiting the Reported Structures in the Protein Data Bank (PDB) since 2015

The structural determination and characterization of molecules, namely proteins and enzymes, is crucial to gaining a better understanding of their role in different chemical and biological processes. The continuous technical developments in the experimental and computational resources of X-ray diffraction (XRD) and, more recently, cryogenic Electron Microscopy (cryo-EM) led to an enormous growth in the number of structures deposited in the Protein Data Bank (PDB). Bioinorganic chemistry arose as a relevant discipline in biology and therapeutics, with a massive number of studies reporting the effects of metal complexes on biological systems, with vanadium complexes being one of the relevant systems addressed. In this review, we focus on the interactions of vanadium compounds (VCs) with proteins. Several types of binding are established between VCs and proteins/enzymes. Considering that the V-species that bind may differ from those initially added, the mentioned structural techniques are pivotal to clarifying the nature and variety of interactions of VCs with proteins and to proposing the mechanisms involved either in enzymatic inhibition or catalysis. As such, we provide an account of the available structural information of VCs bound to proteins obtained by both XRD and/or cryo-EM, mainly exploring the more recent structures, particularly those containing organic-based vanadium complexes.

Identifying structural and dynamic changes during the Biliverdin Reductase B catalytic cycle

Biliverdin Reductase B (BLVRB) is an NADPH-dependent reductase that catalyzes the reduction of multiple substrates and is therefore considered a critical cellular redox regulator. In this study, we sought to address whether both structural and dynamics changes occur between different intermediates of the catalytic cycle and whether these were relegated to just the active site or the entirety of the enzyme. Through X-ray crystallography, we determined the apo BLVRB structure for the first time, revealing subtle global changes compared to the holo structure and identifying the loss of a critical hydrogen bond that "clamps" the R78-loop over the coenzyme. Amide and Cα chemical shift perturbations were used to identify environmental and secondary structural changes between intermediates, with more distant global changes observed upon coenzyme binding compared to substrate interactions. NMR relaxation rate measurements provided insights into the dynamic behavior of BLVRB during the catalytic cycle. Specifically, the inherently dynamic R78-loop that becomes ordered upon coenzyme binding persists through the catalytic cycle while similar regions experience dynamic exchange. However, the dynamic exchange processes were found to differ through the catalytic cycle with several groups of residues exhibiting similar dynamic responses. Finally, both local and distal structural and dynamic changes occur within BLVRB that are dependent solely on the oxidative state of the coenzyme. Thus, through a comprehensive analysis here, this study revealed structural and dynamic alterations in BLVRB through its catalytic cycle that are not simply relegated to the active site, but instead, are allosterically coupled throughout the enzyme.

Analysis of neutral mutational drift in an allosteric enzyme

Neutral mutational drift is an important source of biological diversity that remains underexploited in fundamental studies of protein biophysics. This study uses a synthetic transcriptional circuit to study neutral drift in protein tyrosine phosphatase 1B (PTP1B), a mammalian signaling enzyme for which conformational changes are rate limiting. Kinetic assays of purified mutants indicate that catalytic activity, rather than thermodynamic stability, guides enrichment under neutral drift, where neutral or mildly activating mutations can mitigate the effects of deleterious ones. In general, mutants show a moderate activity-stability tradeoff, an indication that minor improvements in the activity of PTP1B do not require concomitant losses in its stability. Multiplexed sequencing of large mutant pools suggests that substitutions at allosterically influential sites are purged under biological selection, which enriches for mutations located outside of the active site. Findings indicate that the positional dependence of neutral mutations within drifting populations can reveal the presence of allosteric networks and illustrate an approach for using synthetic transcriptional systems to explore these mutations in regulatory enzymes.

Loop dynamics and the evolution of enzyme activity

In the early 2000s, Tawfik presented his 'New View' on enzyme evolution, highlighting the role of conformational plasticity in expanding the functional diversity of limited repertoires of sequences. This view is gaining increasing traction with increasing evidence of the importance of conformational dynamics in both natural and laboratory evolution of enzymes. The past years have seen several elegant examples of harnessing conformational (particularly loop) dynamics to successfully manipulate protein function. This Review revisits flexible loops as critical participants in regulating enzyme activity. We showcase several systems of particular interest: triosephosphate isomerase barrel proteins, protein tyrosine phosphatases and β-lactamases, while briefly discussing other systems in which loop dynamics are important for selectivity and turnover. We then discuss the implications for engineering, presenting examples of successful loop manipulation in either improving catalytic efficiency, or changing selectivity completely. Overall, it is becoming clearer that mimicking nature by manipulating the conformational dynamics of key protein loops is a powerful method of tailoring enzyme activity, without needing to target active-site residues.

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.

A Conserved Local Structural Motif Controls the Kinetics of PTP1B Catalysis

Protein tyrosine phosphatase 1B (PTP1B) is a negative regulator of the insulin and leptin signaling pathways, making it a highly attractive target for the treatment of type II diabetes. For PTP1B to perform its enzymatic function, a loop referred to as the "WPD loop" must transition between open (catalytically incompetent) and closed (catalytically competent) conformations, which have both been resolved by X-ray crystallography. Although prior studies have established this transition as the rate-limiting step for catalysis, the transition mechanism for PTP1B and other PTPs has been unclear. Here we present an atomically detailed model of WPD loop transitions in PTP1B based on unbiased, long-timescale molecular dynamics simulations and weighted ensemble simulations. We found that a specific WPD loop region─the PDFG motif─acted as the key conformational switch, with structural changes to the motif being necessary and sufficient for transitions between long-lived open and closed states of the loop. Simulations starting from the closed state repeatedly visited open states of the loop that quickly closed again unless the infrequent conformational switching of the motif stabilized the open state. The functional importance of the PDFG motif is supported by the fact that it is well conserved across PTPs. Bioinformatic analysis shows that the PDFG motif is also conserved, and adopts two distinct conformations, in deiminases, and the related DFG motif is known to function as a conformational switch in many kinases, suggesting that PDFG-like motifs may control transitions between structurally distinct, long-lived conformational states in multiple protein families.

Deconstructing allostery by computational assessment of the binding determinants of allosteric PTP1B modulators

Fragment-based drug discovery is an established methodology for finding hit molecules that can be elaborated into lead compounds. However it is currently challenging to predict whether fragment hits that do not bind to an orthosteric site could be elaborated into allosteric modulators, as in these cases binding does not necessarily translate into a functional effect. We propose a workflow using Markov State Models (MSMs) with steered molecular dynamics (sMD) to assess the allosteric potential of known binders. sMD simulations are employed to sample protein conformational space inaccessible to routine equilibrium MD timescales. Protein conformations sampled by sMD provide starting points for seeded MD simulations, which are combined into MSMs. The methodology is demonstrated on a dataset of protein tyrosine phosphatase 1B ligands. Experimentally confirmed allosteric inhibitors are correctly classified as inhibitors, whereas the deconstructed analogues show reduced inhibitory activity. Analysis of the MSMs provide insights into preferred protein-ligand arrangements that correlate with functional outcomes. The present methodology may find applications for progressing fragments towards lead molecules in FBDD campaigns.

Mapping the Chemical Space of Active-Site Targeted Covalent Ligands for Protein Tyrosine Phosphatases

Protein tyrosine phosphatases (PTPs) are an important class of enzymes that modulate essential cellular processes through protein dephosphorylation and are dysregulated in various disease states. There is demand for new compounds that target the active sites of these enzymes, for use as chemical tools to dissect their biological roles or as leads for the development of new therapeutics. In this study, we explore an array of electrophiles and fragment scaffolds to investigate the required chemical parameters for covalent inhibition of tyrosine phosphatases. Our analysis juxtaposes the intrinsic electrophilicity of these compounds with their potency against several classical PTPs, revealing chemotypes that inhibit tyrosine phosphatases while minimizing excessive, potentially non-specific reactivity. We also assess sequence divergence at key residues in PTPs to explain their differential susceptibility to covalent inhibition. We anticipate that our study will inspire new strategies to develop covalent probes and inhibitors for tyrosine phosphatases.

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.

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 Mg²⁺. 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.

Allosteric Inhibition of PTP1B by a Nonpolar Terpenoid

Protein tyrosine phosphatases (PTPs) are promising drug targets for treating a wide range of diseases such as diabetes, cancer, and neurological disorders, but their conserved active sites have complicated the design of selective therapeutics. This study examines the allosteric inhibition of PTP1B by amorphadiene (AD), a terpenoid hydrocarbon that is an unusually selective inhibitor. Molecular dynamics (MD) simulations carried out in this study suggest that AD can stably sample multiple neighboring sites on the allosterically influential C-terminus of the catalytic domain. Binding to these sites requires a disordered α7 helix, which stabilizes the PTP1B-AD complex and may contribute to the selectivity of AD for PTP1B over TCPTP. Intriguingly, the binding mode of AD differs from that of the most well-studied allosteric inhibitor of PTP1B. Indeed, biophysical measurements and MD simulations indicate that the two molecules can bind simultaneously. Upon binding, both inhibitors destabilize the α7 helix by disrupting interactions at the α3-α7 interface and prevent the formation of hydrogen bonds that facilitate closure of the catalytically essential WPD loop. These findings indicate that AD is a promising scaffold for building allosteric inhibitors of PTP1B and illustrate, more broadly, how unfunctionalized terpenoids can engage in specific interactions with protein surfaces.

Conserved conformational dynamics determine enzyme activity

Homologous enzymes often exhibit different catalytic rates despite a fully conserved active site. The canonical view is that an enzyme sequence defines its structure and function and, more recently, that intrinsic protein dynamics at different time scales enable and/or promote catalytic activity. Here, we show that, using the protein tyrosine phosphatase PTP1B, residues surrounding the PTP1B active site promote dynamically coordinated chemistry necessary for PTP1B function. However, residues distant to the active site also undergo distinct intermediate time scale dynamics and these dynamics are correlated with its catalytic activity and thus allow for different catalytic rates in this enzyme family. We identify these previously undetected motions using coevolutionary coupling analysis and nuclear magnetic resonance spectroscopy. Our findings strongly indicate that conserved dynamics drives the enzymatic activity of the PTP family. Characterization of these conserved dynamics allows for the identification of novel regulatory elements (therapeutic binding pockets) that can be leveraged for the control of enzymes.

Biochemical, Enzymatic, and Computational Characterization of Recurrent Somatic Mutations of the Human Protein Tyrosine Phosphatase PTP1B in Primary Mediastinal B Cell Lymphoma

Human protein tyrosine phosphatase 1B (PTP1B) is a ubiquitous non-receptor tyrosine phosphatase that serves as a major negative regulator of tyrosine phosphorylation cascades of metabolic and oncogenic importance such as the insulin, epidermal growth factor receptor (EGFR), and JAK/STAT pathways. Increasing evidence point to a key role of PTP1B-dependent signaling in cancer. Interestingly, genetic defects in PTP1B have been found in different human malignancies. Notably, recurrent somatic mutations and splice variants of PTP1B were identified in human B cell and Hodgkin lymphomas. In this work, we analyzed the molecular and functional levels of three PTP1B mutations identified in primary mediastinal B cell lymphoma (PMBCL) patients and located in the WPD-loop (V184D), P-loop (R221G), and Q-loop (G259V). Using biochemical, enzymatic, and molecular dynamics approaches, we show that these mutations lead to PTP1B mutants with extremely low intrinsic tyrosine phosphatase activity that display alterations in overall protein stability and in the flexibility of the active site loops of the enzyme. This is in agreement with the key role of the active site loop regions, which are preorganized to interact with the substrate and to enable catalysis. Our study provides molecular and enzymatic evidence for the loss of protein tyrosine phosphatase activity of PTP1B active-site loop mutants identified in human lymphoma.

Differences in ligand-induced protein dynamics extracted from an unsupervised deep learning approach correlate with protein-ligand binding affinities

Prediction of protein–ligand binding affinity is a major goal in drug discovery. Generally, free energy gap is calculated between two states (e.g., ligand binding and unbinding). The energy gap implicitly includes the effects of changes in protein dynamics induced by ligand binding. However, the relationship between protein dynamics and binding affinity remains unclear. Here, we propose a method that represents ligand-binding-induced protein behavioral change with a simple feature that can be used to predict protein–ligand affinity. From unbiased molecular simulation data, an unsupervised deep learning method measures the differences in protein dynamics at a ligand-binding site depending on the bound ligands. A dimension reduction method extracts a dynamic feature that strongly correlates to the binding affinities. Moreover, the residues that play important roles in protein–ligand interactions are specified based on their contribution to the differences. These results indicate the potential for binding dynamics-based drug discovery.

Differences in ligand-induced protein dynamics extracted as a single feature from a deep learning-based analysis of MD simulations correlate with ligand binding affinity.

Observation of Arginine Side-Chain Motions Coupled to the Global Conformational Exchange Process in Deubiquitinase A

Coupled motions have been demonstrated to be functionally important in a number of enzymes. Noncovalent side-chain interactions play essential roles in coordinating the motions across different structural elements in a protein. However, most of the dynamic studies of proteins are focused on backbone amides or methyl groups in the side chains and little is known about the polar and charged side chains. We have previously characterized the conformational dynamics of deubiquitinase A (DUBA), an isopeptidase, on the microsecond-to-millisecond (μs-ms) time scales with the amide ¹H Carr-Purcell-Meiboom-Gill (CPMG) experiment. We detected a global conformational exchange process on a time scale of approximately 200 μs, which involves most of the structural elements in DUBA, including the active site and the substrate binding interface. Here, we extend our previous study on backbone amides to the arginine side-chain N^ε-H^ε groups using a modified ¹H CPMG pulse sequence that can efficiently detect both backbone amide and arginine side-chain N^ε-H^ε signals in a single experiment. We found that the side chains of three arginines display motions on the same time scale as the backbone amides. Mutations of two of the three arginines to alanines result in a decrease in enzyme activity. One of these two arginines is located in a loop involved in substrate binding. This loop is not visible in the backbone amide-detected experiments due to excess line broadening induced by motions on the μs-ms time scales. These results clearly demonstrate that the motions of some arginine side chains are coupled to the global conformational exchange process and provide an additional probe for motions in a functionally important loop that did not yield visible backbone amide signals, suggesting the value of side-chain experiments on DUBA. The modified ¹H CPMG pulse sequence allows the simultaneous characterization of backbone and arginine side-chain dynamics without any increase in data acquisition time and can be applied to the dynamic studies of any protein that displays measurable amide ¹H relaxation dispersion.

Direct Intracellular Delivery of Benzene Diazonium Ions As Observed by Increased Tyrosine Phosphorylation

A challenge within the field of bioconjugation is developing probes to uncover novel information on proteins and other biomolecules. Intracellular delivery of these probes offers the promise of giving relevant context to this information, and these probes can serve as hypothesis-generating tools within complex systems. Leveraging the utility of triazabutadiene chemistry, herein, we discuss the development of a probe that undergoes reduction-mediated deprotection to rapidly deliver a benzene diazonium ion (BDI) into cells. The intracellular BDI resulted in an increase in global tyrosine phosphorylation levels. Seeing phosphatase dysregulation as a potential source of this increase, a tyrosine phosphatase (PTP1B) was tested and shown to be both inhibited and covalently modified by the BDI. In addition to the expected azobenzene formation at tyrosine side chains, key reactive histidine residues were also modified.

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.

A Computer-Driven Scaffold-Hopping Approach Generating New PTP1B Inhibitors from the Pyrrolo[1,2-a]quinoxaline Core

Protein tyrosine phosphatase 1B (PTP1B) is a very promising target for the treatment of metabolic disorders such as type II diabetes mellitus. Although it was validated as a promising target for this disease more than 30 years ago, as yet there is no drug in advanced clinical trials, and its biochemical mechanism and functions are still being studied. In the present study, based on our experience generating PTP1B inhibitors, we have developed and implemented a scaffold-hopping approach to vary the pyrrole ring of the pyrrolo[1,2-a]quinoxaline core, supported by extensive computational techniques aimed to explain the molecular interaction with PTP1B. Using a combination of docking, molecular dynamics and end-point free-energy calculations, we have rationally designed a hypothesis for new PTP1B inhibitors, supporting their recognition mechanism at a molecular level. After the design phase, we were able to easily synthesize proposed candidates and their evaluation against PTP1B was found to be in good concordance with our predictions. Moreover, the best candidates exhibited glucose uptake increments in cellulo model, thus confirming their utility for PTP1B inhibition and validating this approach for inhibitors design and molecules thus obtained.

Accurate Structure Prediction for Protein Loops Based on Molecular Dynamics Simulations with RSFF2C

Protein loops, connecting the α-helices and β-strands, are involved in many important biological processes. However, due to their conformational flexibility, it is still challenging to accurately determine three-dimensional (3D) structures of long loops experimentally and computationally. Herein, we present a systematic study of the protein loop structure prediction via a total of ∼850 μs molecular dynamics (MD) simulations. For a set of 15 long (10-16 residues) and solvent-exposed loops, we first evaluated the performance of four state-of-the-art loop modeling algorithms, DaReUS-Loop, Sphinx, Rosetta-NGK, and MODELLER, on each loop, and none of them could accurately predict the structures for most loops. Then, temperature replica exchange molecular dynamics (REMD) simulations were conducted with three recent force fields, RSFF2C with TIP3P water model, CHARMM36m with CHARMM-modified TIP3P, and AMBER ff19SB with OPC. We found that our recently developed residue-specific force field RSFF2C performed the best and successfully predicted 12 out of 15 loops with a root-mean-square deviation (RMSD) < 1.5 Å. As an alternative with lower computational cost, normal MD simulations at high temperatures (380, 500, and 620 K) were investigated. Temperature-dependent performance was observed for each force field, and, for RSFF2C+TIP3P, we found that three independent 100-ns MD simulations at 500 K gave comparable results with REMD simulations. These results suggest that MD simulations, especially with enhanced sampling techniques such as replica exchange, with the RSFF2C force field could be useful for accurate loop structure prediction.

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.

Optogenetic Analysis of Allosteric Control in Protein Tyrosine Phosphatases

Allosteric regulation enables dynamic adjustments to protein function that permit tight control over cellular biochemistry. Discrepancies in the allosteric systems of related proteins can thus reveal important differences in their susceptibilities to influential stimuli (e.g., allosteric ligands, mutations, or post-translational modifications). This study uses an optogenetic actuator as a tool to compare the allosteric systems of two structurally related regulatory proteins: protein tyrosine phosphatase 1B (PTP1B) and T-cell protein tyrosine phosphatase (TCPTP). It begins with an interesting observation: The fusion of a protein light switch to the allosterically influential α7 helix of PTP1B permits optical modulation of its catalytic activity, but a similar fusion to TCPTP does not. A subsequent analysis of different PTP chimeras shows that replacing regions of TCPTP with homologous regions from PTP1B can enhance photocontrol; as TCPTP becomes more "PTP1B-like", its photosensitivity increases. Interestingly, the structural changes required for photocontrol also enhance the sensitivity of TCPTP to other allosteric inputs, notably, an allosteric inhibitor and a newly reported activating mutation. Our findings indicate that the allosteric functionality of the α7 helix of PTP1B is not conserved across the PTP family and highlight residues necessary to transfer this functionality to other PTPs. More broadly, our results suggest that simple gene fusion events can strengthen allosteric communication within individual protein domains and describe an intriguing application for optogenetic actuators as structural probes-a sort of physically disruptive "ratchet"-for studying protein allostery.

Hydrogen-Deuterium Exchange within Adenosine Deaminase, a TIM Barrel Hydrolase, Identifies Networks for Thermal Activation of Catalysis

Proteins are intrinsically flexible macromolecules that undergo internal motions with time scales spanning femtoseconds to milliseconds. These fluctuations are implicated in the optimization of reaction barriers for enzyme catalyzed reactions. Time, temperature, and mutation dependent hydrogen-deuterium exchange coupled to mass spectrometry (HDX-MS) has been previously employed to identify spatially resolved, catalysis-linked dynamical regions of enzymes. We now extend this technique to pursue the correlation of protein flexibility and chemical reactivity within the diverse and widespread TIM barrel proteins, targeting murine adenosine deaminase (mADA) that catalyzes the irreversible deamination of adenosine to inosine and ammonia. Following a structure-function analysis of rate and activation energy for a series of mutations at a second sphere phenylalanine positioned in proximity to the bound substrate, the catalytically impaired Phe61Ala with an elevated activation energy (Ea = 7.5 kcal/mol) and the wild type (WT) mADA (Ea = 5.0 kcal/mol) were selected for HDX-MS experiments. The rate constants and activation energies of HDX for peptide segments are quantified and used to assess mutation-dependent changes in local and distal motions. Analyses reveal that approximately 50% of the protein sequence of Phe61Ala displays significant changes in the temperature dependence of HDX behaviors, with the dominant change being an increase in protein flexibility. Utilizing Phe61Ile, which displays the same activation energy for k_cat as WT, as a control, we were able to further refine the HDX analysis, highlighting the regions of mADA that are altered in a functionally relevant manner. A map is constructed that illustrates the regions of protein that are proposed to be essential for the thermal optimization of active site configurations that dominate reaction barrier crossings in the native enzyme.

Cooperative dynamics across distinct structural elements regulate PTP1B activity

Protein-tyrosine phosphatase 1B (PTP1B) is the canonical enzyme for investigating how distinct structural elements influence enzyme catalytic activity. Although it is recognized that dynamics are essential for PTP1B function, the data collected thus far have not resolved whether distinct elements are dynamically coordinated or, alternatively, whether they fulfill their respective functions independently. To answer this question, we performed a comprehensive ¹³C-methyl relaxation study of Ile, Leu, and Val (ILV) residues of PTP1B, which, because of its substantially increased sensitivity, provides a comprehensive understanding of the influence of protein motions on different time scales for enzyme function. We discovered that PTP1B exhibits dynamics at three distinct time scales. First, it undergoes a distinctive slow motion that allows for the dynamic binding and release of its two most N-terminal helices from the catalytic core. Second, we showed that PTP1B ¹³C-methyl group side chain fast time-scale dynamics and ¹⁵N backbone fast time-scale dynamics are fully consistent, demonstrating that fast fluctuations are essential for the allosteric control of PTP1B activity. Third, and most importantly, using ¹³C ILV constant-time Carr-Purcell-Meiboom-Gill relaxation measurements experiments, we demonstrated that all four catalytically important loops-the WPD, Q, E, and substrate-binding loops-work in dynamic unity throughout the catalytic cycle of PTP1B. Thus, these data show that PTP1B activity is not controlled by a single functional element, but instead all key elements are dynamically coordinated. Together, these data provide the first fully comprehensive picture on how the validated drug target PTP1B functions.

Computational insight into the selective allosteric inhibition for PTP1B versus TCPTP: a molecular modelling study

All over the world, diabetes mellitus type 2 has spread as a problematic pandemic. Despite currently available treatments, approved drugs still show undesirable side effects and loss of efficacy or target symptoms instead of causes. Protein tyrosine phosphatase 1B (PTP1B), since its discovery, has emerged as a very promising target against this disease. Although the information regarding the enzyme is immense, little is known about the selectivity between this enzyme and its closest homologue, lymphocyte T tyrosine phosphatase (TCPTP), which is responsible for complicated side effects. In this study, on the basis of different computational approaches, we are able to highlight the importance of a phenylalanine residue located in PTP1B, but not in TCPTP, as a crucial hotspot that causes selectivity and stability for the whole ligand bound system. These results not only allow to explain the selectivity determinants of PTP1B but also provide a useful guide for the design of new allosteric inhibitors. Communicated by Ramaswamy H. Sarma.

Design and synthesis of tricyclic terpenoid derivatives as novel PTP1B inhibitors with improved pharmacological property and in vivo antihyperglycaemic efficacy

Overexpression of protein tyrosine phosphatase 1B (PTP1B) induces insulin resistance in various basic and clinical research. In our previous work, a synthetic oleanolic acid (OA) derivative C10a with PTP1B inhibitory activity has been reported. However, C10a has some pharmacological defects and cytotoxicity. Herein, a structure-based drug design approach was used based on the structure of C10a to elaborate the smaller tricyclic core. A series of tricyclic derivatives were synthesised and the compounds 15, 28 and 34 exhibited the most PTP1B enzymatic inhibitory potency. In the insulin-resistant human hepatoma HepG2 cells, compound 25 with the moderate PTP1B inhibition and preferable pharmaceutical properties can significantly increase insulin-stimulated glucose uptake and showed the insulin resistance ameliorating effect. Moreover, 25 showed the improved in vivo antihyperglycaemic potential in the nicotinamide–streptozotocin-induced T2D. Our study demonstrated that these tricyclic derivatives with improved molecular architectures and antihyperglycaemic activity could be developed in the treatment of T2D.