Loop Dynamics and Enzyme Catalysis in Protein Tyrosine Phosphatases

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.

Rational Design of Loop Dynamics for a Barrel-Shaped Enzyme by Introducing Disulfide Bonds

Loop dynamics redesign is an important strategy to manipulate protein function. Cellobiose 2-epimerase (CE) and other members of its superfamily are widely used for diverse industrial applications. The structural feature of the loops connecting barrel helices contributes greatly to the differences in their functional characteristics. Inspired by the in-silico mutation with molecular dynamics (MD) simulation analysis, we propose a strategy for identifying disulfide bond mutation candidates based on the prediction of protein flexibility and residue-residue interaction. The most beneficial mutant with the newly introduced disulfide bond would simultaneously improve both its thermostability and its reaction propensity to the targeting isomerization product. The ratio of the isomerization/epimerization catalytic rate was improved from 4:103 to 9:22. MD simulation and binding free energy calculations were applied to provide insights into molecular recognition upon mutations. The comparative analysis of enzyme/substrate binding modes indicates that the altered catalytic reaction pathway is due to less efficient binding of the native product. The key residue responsible for the observed phenotype was identified by energy decomposition and was further confirmed by the mutation experiment. The rational design of the key loop region might be a promising strategy to alter the catalytic behavior of all (α/α)6-barrel-like proteins.

A novel gain-of-function phosphorylation site modulates PTPN22 inhibition of TCR signaling

Protein tyrosine phosphatase nonreceptor type 22 (PTPN22) is encoded by a major autoimmunity gene and is a known inhibitor of T cell receptor (TCR) signaling and drug target for cancer immunotherapy. However, little is known about PTPN22 posttranslational regulation. Here, we characterize a phosphorylation site at Ser325 situated C terminal to the catalytic domain of PTPN22 and its roles in altering protein function. In human T cells, Ser325 is phosphorylated by glycogen synthase kinase-3 (GSK3) following TCR stimulation, which promotes its TCR-inhibitory activity. Signaling through the major TCR-dependent pathway under PTPN22 control was enhanced by CRISPR/Cas9-mediated suppression of Ser325 phosphorylation and inhibited by mimicking it via glutamic acid substitution. Global phospho-mass spectrometry showed Ser325 phosphorylation state alters downstream transcriptional activity through enrichment of Swi3p, Rsc8p, and Moira domain binding proteins, and next-generation sequencing revealed it differentially regulates the expression of chemokines and T cell activation pathways. Moreover, in vitro kinetic data suggest the modulation of activity depends on a cellular context. Finally, we begin to address the structural and mechanistic basis for the influence of Ser325 phosphorylation on the protein's properties by deuterium exchange mass spectrometry and NMR spectroscopy. In conclusion, this study explores the function of a novel phosphorylation site of PTPN22 that is involved in complex regulation of TCR signaling and provides details that might inform the future development of allosteric modulators of PTPN22.

Protein tyrosine phosphatase inhibitors: a patent review and update (2012-2023)

INTRODUCTION

Protein tyrosine phosphatases (PTPs), essential and evolutionarily highly conserved enzymes, govern cellular functions by modulating tyrosine phosphorylation, a pivotal post-translational modification for signal transduction. The recent strides in phosphatase drug discovery, leading to the identification of selective modulators for enzymes, restoring interest in the therapeutic targeting of protein phosphatases.

AREAS COVERED

The compilation of patents up to the year 2023 focuses on the efficacy of various classes of Tyrosine phosphatases and their inhibitors, detailing their chemical structure and biochemical characteristics. These findings have broad implications, as they can be applied to treating diverse conditions like cancer, diabetes, autoimmune disorders, and neurological diseases. The search for scientific articles and patent literature was conducted using well known different platforms to gather information up to 2023.

EXPERT OPINION

The latest improvements in protein tyrosine phosphatase (PTP) research include the discovery of new inhibitors targeting specific PTP enzymes, with a focus on developing allosteric site covalent inhibitors for enhanced efficacy and specificity. These advancements have not only opened up new possibilities for therapeutic interventions in various disease conditions but also hold the potential for innovative treatments. PTPs offer promising avenues for drug discovery efforts and innovative treatments across a spectrum of health conditions.

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.

Druggable targets of protein tyrosine phosphatase Family, viz. PTP1B, SHP2, Cdc25, and LMW-PTP: Current scenario on medicinal Attributes, and SAR insights

Protein tyrosine phosphatases (PTPs) are the class of dephosphorylation enzymes that catalyze the removal of phosphate groups from tyrosine residues on proteins responsible for various cellular processes. Any disbalance in signal pathways mediated by PTPs leads to various disease conditions like diabetes, obesity, cancers, and autoimmune disorders. Amongst the PTP superfamily, PTP1B, SHP2, Cdc25, and LMW-PTP have been prioritized as druggable targets for developing medicinal agents. PTP1B is an intracellular PTP enzyme that downregulates insulin and leptin signaling pathways and is involved in insulin resistance and glucose homeostasis. SHP2 is involved in the RAS-MAPK pathway and T cell immunity. Cdk-cyclin complex activation occurs by Cdc25-PTPs involved in cell cycle regulation. LMW-PTPs are involved in PDGF/PDGFR, Eph/ephrin, and insulin signaling pathways, resulting in certain diseases like diabetes mellitus, obesity, and cancer. The signaling cascades of PTP1B, SHP2, Cdc25, and LMW-PTPs have been described to rationalize their medicinal importance in the pathophysiology of diabetes, obesity, and cancer. Their binding sites have been explored to overcome the hurdles in discovering target selective molecules with optimum potency. Recent developments in the synthetic molecules bearing heterocyclic moieties against these targets have been explored to gain insight into structural features. The elaborated SAR investigation revealed the effect of substituents on the potency and target selectivity, which can be implicated in the further discovery of newer medicinal agents targeting the druggable members of the PTP superfamily.

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.

Comprehensive assessment of protein loop modeling programs on large-scale datasets: prediction accuracy and efficiency

Protein loops play a critical role in the dynamics of proteins and are essential for numerous biological functions, and various computational approaches to loop modeling have been proposed over the past decades. However, a comprehensive understanding of the strengths and weaknesses of each method is lacking. In this work, we constructed two high-quality datasets (i.e. the General dataset and the CASP dataset) and systematically evaluated the accuracy and efficiency of 13 commonly used loop modeling approaches from the perspective of loop lengths, protein classes and residue types. The results indicate that the knowledge-based method FREAD generally outperforms the other tested programs in most cases, but encountered challenges when predicting loops longer than 15 and 30 residues on the CASP and General datasets, respectively. The ab initio method Rosetta NGK demonstrated exceptional modeling accuracy for short loops with four to eight residues and achieved the highest success rate on the CASP dataset. The well-known AlphaFold2 and RoseTTAFold require more resources for better performance, but they exhibit promise for predicting loops longer than 16 and 30 residues in the CASP and General datasets. These observations can provide valuable insights for selecting suitable methods for specific loop modeling tasks and contribute to future advancements in the field.

Loop Evolutionary Patterns Shape Catalytic Efficiency of TRI101/201 for Trichothecenes: Insights into Protein-Substrate Interactions

Trichothecenes are highly toxic mycotoxins produced by Fusarium fungi, while TRI101/201 family enzymes play a crucial role in detoxification through acetylation. Studies on the substrate specificity and catalytic kinetics of TRI101/201 have revealed distinct kinetic characteristics, with significant differences observed in catalytic efficiency toward deoxynivalenol, while the catalytic efficiency for T-2 toxin remains relatively consistent. In this study, we used structural bioinformatics analysis and a molecular dynamics simulation workflow to investigate the mechanism underlying the differential catalytic activity of TRI101/201. The findings revealed that the binding stability between trichothecenes and TRI101/201 hinges primarily on a hydrophobic cage structure within the binding site. An intrinsic disordered loop, termed loop cover, defined the evolutionary patterns of the TRI101/201 protein family that are categorized into four subfamilies (V1/V2/V3/M). Furthermore, the unique loop displayed different conformations among these subfamilies' structures, which served to disrupt (V1/V2/V3) or reinforce (M) the hydrophobic cages. The disrupted cages enhanced the water exposure of the hydrophilic moieties of substrates like deoxynivalenol and thereby hindered their binding to the catalytic sites of V-type enzymes. In contrast, this water exposure does not affect substrates like T-2 toxin, which have more hydrophobic substituents, resulting in a comparable catalytic efficiency of both V- and M-type enzymes. Overall, our studies provide theoretical support for understanding the catalytic mechanism of TRI101/201, which shows how an intrinsic disordered loop could impact the protein-ligand binding and suggests a direction for rational protein design in the future.

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.

Setting sail: Maneuvering SHP2 activity and its effects in cancer

Since the discovery of tyrosine phosphorylation being a critical modulator of cancer signaling, proteins regulating phosphotyrosine levels in cells have fast become targets of therapeutic intervention. The nonreceptor protein tyrosine phosphatase (PTP) coded by the PTPN11 gene "SHP2" integrates phosphotyrosine signaling from growth factor receptors into the RAS/RAF/ERK pathway and is centrally positioned in processes regulating cell development and oncogenic transformation. Dysregulation of SHP2 expression or activity is linked to tumorigenesis and developmental defects. Even as a compelling anti-cancer target, SHP2 was considered "undruggable" for a long time owing to its conserved catalytic PTP domain that evaded drug development. Recently, SHP2 has risen from the "undruggable curse" with the discovery of small molecules that manipulate its intrinsic allostery for effective inhibition. SHP2's unique domain arrangement and conformation(s) allow for a truly novel paradigm of inhibitor development relying on skillful targeting of noncatalytic sites on proteins. In this review we summarize the biological functions, signaling properties, structural attributes, allostery and inhibitors of SHP2.

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.

Allosteric regulation of β-reaction stage I in tryptophan synthase upon the α-ligand binding

Tryptophan synthase (TRPS) is a bifunctional enzyme consisting of α- and β-subunits that catalyzes the last two steps of L-tryptophan (L-Trp) biosynthesis. The first stage of the reaction at the β-subunit is called β-reaction stage I, which converts the β-ligand from an internal aldimine [E(Ain)] to an α-aminoacrylate [E(A-A)] intermediate. The activity is known to increase 3-10-fold upon the binding of 3-indole-D-glycerol-3'-phosphate (IGP) at the α-subunit. The effect of α-ligand binding on β-reaction stage I at the distal β-active site is not well understood despite the abundant structural information available for TRPS. Here, we investigate the β-reaction stage I by carrying out minimum-energy pathway searches based on a hybrid quantum mechanics/molecular mechanics (QM/MM) model. The free-energy differences along the pathway are also examined using QM/MM umbrella sampling simulations with QM calculations at the B3LYP-D3/aug-cc-pVDZ level of theory. Our simulations suggest that the sidechain orientation of βD305 near the β-ligand likely plays an essential role in the allosteric regulation: a hydrogen bond is formed between βD305 and the β-ligand in the absence of the α-ligand, prohibiting a smooth rotation of the hydroxyl group in the quinonoid intermediate, whereas the dihedral angle rotates smoothly after the hydrogen bond is switched from βD305-β-ligand to βD305-βR141. This switch could occur upon the IGP-binding at the α-subunit, as evidenced by the existing TRPS crystal structures.

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.

Substantial Improvement of an Epimerase for the Synthesis of D-Allulose by Biosensor-Based High-Throughput Microdroplet Screening

Biosynthesis of D-allulose has been achieved using ketose 3-epimerases (KEases), but its application is limited by poor catalytic performance. In this study, we redesigned a genetically encoded biosensor based on a D-allulose-responsive transcriptional regulator for real-time monitoring of D-allulose. An ultrahigh-throughput droplet-based microfluidic screening platform was further constructed by coupling with this D-allulose-detecting biosensor for the directed evolution of the KEases. Structural analysis of Sinorhizobium fredii D-allulose 3-epimerase (SfDAE) revealed that a highly flexible helix/loop region exposes or occludes the catalytic center as an essential lid conformation regulating substrate recognition. We reprogrammed SfDAE using structure-guided rational design and directed evolution, in which a mutant M3-2 was identified with 17-fold enhanced catalytic efficiency. Our research offers a paradigm for the design and optimization of a biosensor-based microdroplet screening platform.

Modern non-polarizable force fields diverge in modeling the enzyme-substrate complex of a canonical serine protease

Classical molecular dynamics simulation is a powerful and established method of modern computational chemistry. Being able to obtain accurate information on molecular behavior is crucial to get valuable insights into structure-function relationships that translate into fundamental findings and practical applications. Active sites of enzymes are known to be particularly intricate, therefore, simpler non-polarizable force fields may provide an inaccurate description. In this work, we addressed this hypothesis in a case of a canonical serine triad protease trypsin in its complex with a substrate-mimicking inhibitor. We tested six modern and popular force fields to find that significantly diverging results may be obtained. Amber FB-15 and OPLS-AA/M turned out to model the active site incorrectly. Amber ff19sb and ff15ipq demonstrated mixed performance. The best performing force fields were CHARMM36m and Amber ff99sb-ildn, therefore, they are recommended for use with this and related systems. We speculate that a similar lack of cross-force field convergence may be characteristic of other enzymatic systems. Therefore, we advocate for careful consideration of different force fields in any study within the field of computational enzymology.

Cooperation of structural motifs controls drug selectivity in cyclin-dependent kinases: an advanced theoretical analysis

Understanding drug selectivity mechanism is a long-standing issue for helping design drugs with high specificity. Designing drugs targeting cyclin-dependent kinases (CDKs) with high selectivity is challenging because of their highly conserved binding pockets. To reveal the underlying general selectivity mechanism, we carried out comprehensive analyses from both the thermodynamics and kinetics points of view on a representative CDK12 inhibitor. To fully capture the binding features of the drug-target recognition process, we proposed to use kinetic residue energy analysis (KREA) in conjunction with the community network analysis (CNA) to reveal the underlying cooperation effect between individual residues/protein motifs to the binding/dissociating process of the ligand. The general mechanism of drug selectivity in CDKs can be summarized as that the difference of structural cooperation between the ligand and the protein motifs leads to the difference of the energetic contribution of the key residues to the ligand. The proposed mechanisms may be prevalent in drug selectivity issues, and the insights may help design new strategies to overcome/attenuate the drug selectivity associated problems.

Recent applications of computational methods to allosteric drug discovery

Interest in exploiting allosteric sites for the development of new therapeutics has grown considerably over the last two decades. The chief driving force behind the interest in allostery for drug discovery stems from the fact that in comparison to orthosteric sites, allosteric sites are less conserved across a protein family, thereby offering greater opportunity for selectivity and ultimately tolerability. While there is significant overlap between structure-based drug design for orthosteric and allosteric sites, allosteric sites offer additional challenges mostly involving the need to better understand protein flexibility and its relationship to protein function. Here we examine the extent to which structure-based drug design is impacting allosteric drug design by highlighting several targets across a variety of target classes.

Ligands-induced open-close conformational change during DapE catalysis: Insights from molecular dynamics simulations

The microbial enzyme DapE plays a critical role in the lysine biosynthetic pathway and is considered as a potentially safe antibiotic target. In this study, atomistic simulations are employed to identify the modes of essential dynamics that define the conformational response of the enzyme to ligand binding and unbinding. The binding modes and the binding affinities of the products to the DapE enzyme are estimated from the MM-PBSA method, and the residues contributing to the ligand binding are identified. Various structural analyses and the principal component analysis of the molecular dynamics trajectories reveal that the removal of products from the active site causes a significant change in the overall enzyme structure. Both Cartesian and dihedral principal component analyses are used to characterize the structural changes in terms of domain unfolding and domain twisting motions. In the most dominant mode, that is, the domain unfolding motion, the two catalytic domains move away from the two dimerization domains of the dimeric enzyme, representing a closed-to-open conformational change. The conformational changes are initiated by the coordinated movement of three loops (Asp75-Pro82, Gly240-Asn244, and Thr347-Glu353) that trigger a domain-level movement. From multiple short trajectories, the time constant associated with the domain opening motion is estimated as 43.6 ns. Physiologically, this close-to-open conformational change is essential for the regeneration of the initial state of the enzyme for the subsequent cycle of catalytic action and provides the apo enzyme enough flexibility for efficient substrate binding.

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.

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.

Can Second Coordination Sphere and Long-Range Interactions Modulate Hydrogen Atom Transfer in a Non-Heme Fe(II)-Dependent Histone Demethylase?

Fe(II)-dependent oxygenases employ hydrogen atom transfer (HAT) to produce a myriad of products. Understanding how such enzymes use dynamic processes beyond the immediate vicinity of the active site to control the selectivity and efficiency of HAT is important for metalloenzyme engineering; however, obtaining such knowledge by experiments is challenging. This study develops a computational framework for identifying second coordination sphere (SCS) and especially long-range (LR) residues relevant for catalysis through dynamic cross-correlation analysis (DCCA) using the human histone demethylase PHF8 (KDM7B) as a model oxygenase. Furthermore, the study explores the mechanistic pathways of influence of the SCS and LR residues on the HAT reaction. To demonstrate the plausibility of the approach, we investigated the effect of a PHF8 F279S clinical mutation associated with X-linked intellectual disability, which has been experimentally shown to ablate PHF8-catalyzed demethylation. In agreement, the molecular dynamics (MD) and quantum mechanics/molecular mechanics (QM/MM) studies showed a change in the H3_1-14K9me2 substrate orientation and an increased HAT barrier. We systematically analyzed the pathways by which the identified SCS and LR residues may influence HAT by exploring changes in H3K9me2 substrate orientation, interdomain correlated motions, HAT transition state stabilization, reaction energetics, electron transfer mechanism, and alterations in the intrinsic electric field of PHF8. Importantly, SCS and LR variations decrease key motions of α9-α12 of the JmjC domain toward the Fe(IV)-center that are associated with tighter binding of the H3_1-14K9me2 substrate. SCS and LR residues alter the intrinsic electric field of the enzyme along the reaction coordinate and change the individual energetic contributions of residues toward TS stabilization. The overall results suggest that DCCA can indeed identify non-active-site residues relevant for catalysis. The substitutions of such dynamically correlated residues might be used as a tool to tune HAT in non-heme Fe(II)- and 2OG-dependent enzymes.

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.

Targeting protein tyrosine phosphatases for the development of antivirulence agents: Yersinia spp. and Mycobacterium tuberculosis as prototypes

Protein phosphorylation mediated by protein kinases and phosphatases has a central regulatory function in many cellular processes in eukaryotes and prokaryotes. As a result, several diseases caused by imbalance in phosphorylation levels are known, especially due to protein tyrosine phosphatases (PTPs) activity, an important family of signaling enzymes. Furthermore, over the last decades several studies have shown the main role of PTPs in pathogenic bacteria: they are associated with growth, cell division, cell wall biosynthesis, biofilm formation, metabolic processes, as well as virulence factor. In this way, PTPs have ascended as targets for antibacterial drug design, particularly in view of the antibiotic resistance in pathogenic bacteria, which demands novel therapeutics strategies. Targeting secreted PTPs is an antivirulence strategy to combat the emergence of antimicrobial resistance (AMR). This review focuses on the recent advances in understanding the role of PTPs and the approaches to target them, with an emphasis in Yersinia spp. and Mycobacterium tuberculosis pathogenesis.

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.

A structural exposé of noncanonical molecular reactivity within the protein tyrosine phosphatase WPD loop

Structural snapshots of protein/ligand complexes are a prerequisite for gaining atomic level insight into enzymatic reaction mechanisms. An important group of enzymes has been deprived of this analytical privilege: members of the protein tyrosine phosphatase (PTP) superfamily with catalytic WPD-loops lacking the indispensable general-acid/base within a tryptophan-proline-aspartate/glutamate context. Here, we provide the ligand/enzyme crystal complexes for one such PTP outlier: Arabidopsis thaliana Plant and Fungi Atypical Dual Specificity Phosphatase 1 (AtPFA-DSP1), herein unveiled as a regioselective and efficient phosphatase towards inositol pyrophosphate (PP-InsP) signaling molecules. Although the WPD loop is missing its canonical tripeptide motif, this structural element contributes to catalysis by assisting PP-InsP delivery into the catalytic pocket, for a choreographed exchange with phosphate reaction product. Subsequently, an intramolecular proton donation by PP-InsP substrate is posited to substitute functionally for the absent aspartate/glutamate general-acid. Overall, we expand mechanistic insight into adaptability of the conserved PTP structural elements.

LoopGrafter: a web tool for transplanting dynamical loops for protein engineering

The transplantation of loops between structurally related proteins is a compelling method to improve the activity, specificity and stability of enzymes. However, despite the interest of loop regions in protein engineering, the available methods of loop-based rational protein design are scarce. One particular difficulty related to loop engineering is the unique dynamism that enables them to exert allosteric control over the catalytic function of enzymes. Thus, when engaging in a transplantation effort, such dynamics in the context of protein structure need consideration. A second practical challenge is identifying successful excision points for the transplantation or grafting. Here, we present LoopGrafter (https://loschmidt.chemi.muni.cz/loopgrafter/), a web server that specifically guides in the loop grafting process between structurally related proteins. The server provides a step-by-step interactive procedure in which the user can successively identify loops in the two input proteins, calculate their geometries, assess their similarities and dynamics, and select a number of loops to be transplanted. All possible different chimeric proteins derived from any existing recombination point are calculated, and 3D models for each of them are constructed and energetically evaluated. The obtained results can be interactively visualized in a user-friendly graphical interface and downloaded for detailed structural analyses.

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.

Design and engineering of allosteric communications in proteins

Allostery in proteins plays an important role in regulating protein activities and influencing many biological processes such as gene expression, enzyme catalysis, and cell signaling. The process of allostery takes place when a signal detected at a site on a protein is transmitted via a mechanical pathway to a functional site and, thus, influences its activity. The pathway of allosteric communication consists of amino acids that form a network with covalent and non-covalent bonds. By mutating residues in this allosteric network, protein engineers have successfully established novel allosteric pathways to achieve desired properties in the target protein. In this review, we highlight the most recent and state-of-the-art techniques for allosteric communication engineering. We also discuss the challenges that need to be overcome and future directions for engineering protein allostery.

A distal regulatory strategy of enzymes: from local to global conformational dynamics

Modulating the distribution of various states in protein ensembles through distal sites may be promising in the evolution of enzymes in desired directions. However, the prediction of distal mutation hotspots that stabilize the favoured states from a computational perspective remains challenging. Here, we presented a strategy based on molecular dynamics (MD) and Markov state models (MSM) to predict distal mutation sites. Extensive MD combined with MSM was applied to determine the principally distributed metastable states interconverting at a slow timescale. Then, molecular docking was used to classify these states into active states and inactive ones. Distal mutation hotspots were targeted based on comparing the conformational features between active and inactive states, where mutations destabilize the inactive states and show little influence on the active state. The proposed strategy was used to explore the highly dynamic MHETase, which shows a potential application in the biodegradation of poly(ethylene terephthalate) (PET). Seven principally populated interrelated metastable states were identified, and the atomistic picture of their conformational changes was unveiled. Several residues at distal positions were found to adopt more H-bond occupancies in inactive states than active states, making them potential mutation hotspots for stabilizing the favoured conformations. In addition, the detailed mechanism revealed the significance of calcium ions at a distance from the catalytic centre in reshaping the free energy landscape. This study deepens the understanding of the conformational dynamics of α/β hydrolases containing a lid domain and advances the study of enzymatic plastic degradation.

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.

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.