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

Activation and friction in enzymatic loop opening and closing dynamics

Protein loop dynamics have recently been recognized as central to enzymatic activity, specificity and stability. However, the factors controlling loop opening and closing kinetics have remained elusive. Here, we combine molecular dynamics simulations with string-method determination of complex reaction coordinates to elucidate the molecular mechanism and rate-limiting step for WPD-loop dynamics in the PTP1B enzyme. While protein conformational dynamics is often represented as diffusive motion hindered by solvent viscosity and internal friction, we demonstrate that loop opening and closing is activated. It is governed by torsional rearrangement around a single loop peptide group and by significant friction caused by backbone adjustments, which can dynamically trap the loop. Considering both torsional barrier and time-dependent friction, our calculated rate constants exhibit very good agreement with experimental measurements, reproducing the change in loop opening kinetics between proteins. Furthermore, we demonstrate the applicability of our results to other enzymatic loops, including the M20 DHFR loop, thereby offering prospects for loop engineering potentially leading to enhanced designs.

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.

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.

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.

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.

Protein Dynamics and Enzymatic Catalysis

This Perspective presents a review of our work and that of others in the highly controversial topic of the coupling of protein dynamics to reaction in enzymes. We have been involved in studying this topic for many years. Thus, this perspective will naturally present our own views, but it also is designed to present an overview of the variety of viewpoints of this topic, both experimental and theoretical. This is obviously a large and contentious topic.

Loop Dynamics and Conformational Flexibility in Dengue Serine Protease Activity: Noninvasive Perturbation by Solvent Exchange

Molecular mechanics play an important role in enzyme action and understanding the dynamics of loop motion is key for designing inhibitors of an enzyme, particularly targeting the allosteric sites. For the successful creation of new protease inhibitors targeting the dengue serine protease, our current investigation detailed the intricate structural dynamics of NS2B/NS3 dengue protease. This enzyme is one of the most essential enzymes in the life cycle of the dengue virus, which is responsible for the activation/processing of viral polyprotein, thus making it a potential target for drug discovery. We showed that the internal dynamics of two regions, fingers 1 and 2 (R24-G39 and L149-A164, respectively) adjacent to the active site triad of this protease, control the enzyme action. Each of these regions is composed of two antiparallel β-strands connected by β-turn/hairpin loops. The correlated bending and rocking motions in the two β-turns on either side of the active site were found to modulate the activity of the enzyme to a large extent. With increasing concentration of cosolvent dimethyl sulfoxide, correlated motions in the finger 2 region get diminished and bending of finger 1 increases, which are also reflected in the loss of enzyme activity. Decreasing temperature and mutations in neighboring nonsubstrate binding residues show similar effects on loop motion and enzyme kinetics. Therefore, in vitro noninvasive perturbation of these motions by the solvent exchange as well as cold stress in combination with in silico molecular dynamics simulations established the importance of the two β-turns in the functioning of dengue virus serotype 2 NS2B/NS3 serine protease.

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.

Increasing the O2 Resistance of the [FeFe]-Hydrogenase CbA5H through Enhanced Protein Flexibility

The high turnover rates of [FeFe]-hydrogenases under mild conditions and at low overpotentials provide a natural blueprint for the design of hydrogen catalysts. However, the unique active site (H-cluster) degrades upon contact with oxygen. The [FeFe]-hydrogenase fromClostridium beijerinckii (CbA5H) is characterized by the flexibility of its protein structure, which allows a conserved cysteine to coordinate to the active site under oxidative conditions. Thereby, intrinsic cofactor degradation induced by dioxygen is minimized. However, the protection from O₂ is only partial, and the activity of the enzyme decreases upon each exposure to O₂. By using site-directed mutagenesis in combination with electrochemistry, ATR-FTIR spectroscopy, and molecular dynamics simulations, we show that the kinetics of the conversion between the oxygen-protected inactive state (cysteine-bound) and the oxygen-sensitive active state can be accelerated by replacing a surface residue that is very distant from the active site. This sole exchange of methionine for a glutamate residue leads to an increased resistance of the hydrogenase to dioxygen. With our study, we aim to understand how local modifications of the protein structure can have a crucial impact on protein dynamics and how they can control the reactivity of inorganic active sites through outer sphere effects.

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.

Complex Loop Dynamics Underpin Activity, Specificity, and Evolvability in the (βα)8 Barrel Enzymes of Histidine and Tryptophan Biosynthesis

Enzymes are conformationally dynamic, and their dynamical properties play an important role in regulating their specificity and evolvability. In this context, substantial attention has been paid to the role of ligand-gated conformational changes in enzyme catalysis; however, such studies have focused on tremendously proficient enzymes such as triosephosphate isomerase and orotidine 5'-monophosphate decarboxylase, where the rapid (μs timescale) motion of a single loop dominates the transition between catalytically inactive and active conformations. In contrast, the (βα)₈-barrels of tryptophan and histidine biosynthesis, such as the specialist isomerase enzymes HisA and TrpF, and the bifunctional isomerase PriA, are decorated by multiple long loops that undergo conformational transitions on the ms (or slower) timescale. Studying the interdependent motions of multiple slow loops, and their role in catalysis, poses a significant computational challenge. This work combines conventional and enhanced molecular dynamics simulations with empirical valence bond simulations to provide rich details of the conformational behavior of the catalytic loops in HisA, PriA, and TrpF, and the role of their plasticity in facilitating bifunctionality in PriA and evolved HisA variants. In addition, we demonstrate that, similar to other enzymes activated by ligand-gated conformational changes, loops 3 and 4 of HisA and PriA act as gripper loops, facilitating the isomerization of the large bulky substrate ProFAR, albeit now on much slower timescales. This hints at convergent evolution on these different (βα)₈-barrel scaffolds. Finally, our work reemphasizes the potential of engineering loop dynamics as a tool to artificially manipulate the catalytic repertoire of TIM-barrel proteins.

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.

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.

Harnessing Conformational Plasticity to Generate Designer Enzymes

Recent years have witnessed an explosion of interest in understanding the role of conformational dynamics both in the evolution of new enzymatic activities from existing enzymes and in facilitating the emergence of enzymatic activity de novo on scaffolds that were previously non-catalytic. There are also an increasing number of examples in the literature of targeted engineering of conformational dynamics being successfully used to alter enzyme selectivity and activity. Despite the obvious importance of conformational dynamics to both enzyme function and evolvability, many (although not all) computational design approaches still focus either on pure sequence-based approaches or on using structures with limited flexibility to guide the design. However, there exist a wide variety of computational approaches that can be (re)purposed to introduce conformational dynamics as a key consideration in the design process. Coupled with laboratory evolution and more conventional existing sequence- and structure-based approaches, these techniques provide powerful tools for greatly expanding the protein engineering toolkit. This Perspective provides an overview of evolutionary studies that have dissected the role of conformational dynamics in facilitating the emergence of novel enzymes, as well as advances in computational approaches that allow one to target conformational dynamics as part of enzyme design. Harnessing conformational dynamics in engineering studies is a powerful paradigm with which to engineer the next generation of designer biocatalysts.

Hybrid Thermophilic/Mesophilic Enzymes Reveal a Role for Conformational Disorder in Regulation of Bacterial Enzyme I

Conformational disorder is emerging as an important feature of biopolymers, regulating a vast array of cellular functions, including signaling, phase separation, and enzyme catalysis. Here we combine NMR, crystallography, computer simulations, protein engineering, and functional assays to investigate the role played by conformational heterogeneity in determining the activity of the C-terminal domain of bacterial Enzyme I (EIC). In particular, we design chimeric proteins by hybridizing EIC from thermophilic and mesophilic organisms, and we characterize the resulting constructs for structure, dynamics, and biological function. We show that EIC exists as a mixture of active and inactive conformations and that functional regulation is achieved by tuning the thermodynamic balance between active and inactive states. Interestingly, we also present a hybrid thermophilic/mesophilic enzyme that is thermostable and more active than the wild-type thermophilic enzyme, suggesting that hybridizing thermophilic and mesophilic proteins is a valid strategy to engineer thermostable enzymes with significant low-temperature activity.