- 1Center for Structure-based Drug Design and Development, Department of Pharmaceutical Sciences, Concordia University Wisconsin, Mequon, WI, United States
- 2Department of Chemistry and Biochemistry, New Mexico State University, MSC 3C, Las Cruces, NM, United States
- 3Department of Pediatrics, Division of Neonatology, Children’s Research Institute (CRI) Developmental Vascular Biology Program, Translational and Biomedical Research Center, Milwaukee, WI, United States
- 4Department of Chemistry, Colorado State University, Fort Collins, CO, United States
- 5Department of Pharmacology and Toxicology, Medical College of Wisconsin, Milwaukee, WI, United States
- 6Cell and Molecular Biology Program, Colorado State University, Fort Collins, CO, United States
Introduction: Dual specific phosphatases (DUSPs) are mitogen-activated protein kinase (MAPK) regulators, which also serve as drug targets for treating various vascular diseases. Previously, we have presented mechanistic characterizations of DUSP5 and its interaction with pERK, proposing a dual active site.
Methods: Herein, we characterize the interactions between the DUSP5 phosphatase domain and the pT-E-pY activation loop of ERK2, with specific active site assignments. We also report the full NMR chemical shift assignments of DUSP5 that now enable chemical shift perturbation and dynamics studies.
Results and Discussion: Both phosphates of the pT-E-pY tripeptide are dephosphorylated, based on 31P NMR; but, steady state kinetic studies of the tripeptide both as a substrate and as an inhibitor indicate a preference for binding and dephosphorylation of the phospho-tyrosine before the phospho-threonine. Catalytic efficiency (kcat/Km) is 3.7 M−1S−1 for T-E-pY vs 1.3 M−1S−1 for pT-E-Y, although the diphosphorylated peptide (pT-E-pY) is a better substrate than both, with kcat/Km = 18.2 M−1S−1 . Steady state inhibition studies with the pNPP substrate yields Kis values for the peptide inhibitors of: 15.82 mM (pT-E-Y), 4.932 mM (T-E-pY), 1.672 mM (pT-E-pY). Steady state inhibition studies with pNPP substrate and with vanadate or phosphate inhibitors indicated competitive inhibition with Kis values of 0.0006122 mM (sodium vanadate) and 17.32 mM (sodium phosphate), similar to other Protein Tyrosine Phosphatases with an active site cysteine nucleophile that go through a five-coordinate high energy transition state or intermediate. Molecular dynamics (MD) studies confirm preferential binding of the diphosphorylated peptide, but with preference for binding the pY over the pT reside in the catalytic site proximal to the Cys263 nucleophile. Based on MD, the monophosphorylated peptide binds tighter if phosphorylated on the Tyr vs the Thr. And, if the starting pose of the docked diphosphorylated peptide has pT in the catalytic site, it will adjust to have the pY in the catalytic site, suggesting a dynamic shifting of the peptide orientation. 2D 1H-15N HSQC chemical shift perturbation studies confirm that DUSP5 with tripeptide bound is in a dynamic state, with extensive exchange broadening observed—especially of catalytic site residues. The availability of NMR chemical shift assignments enables additional future studies of DUSP5 binding to the ERK2 diphosphorylated activation loop.
Summary: These studies indicate a preference for pY before pT binding, but with ability to bind and dephosphorylate both residues, and with a dynamic active site pocket that accommodates multiple tripeptide orientations.
1 Introduction
In eukaryotes, the prototypical mitogen-activated protein kinases (MAPK) have been widely studied and are comprised of a three-tier protein phosphorylation system (Cargnello and Roux, 2011; Raman et al., 2007). In this system, a series of proteins get phosphorylated in response to a ligand-receptor interaction and the signal leads into the nucleus to dictate gene expression. Extracellular regulated kinase (ERK), c-Jun N-terminal kinase (JNK) and p38 kinase are the terminal kinases in this prototypical pathway (Figure 1). As one would expect, exquisite control of the MAPK pathway is established by a second group of proteins called phosphatases (Caunt and Keyse, 2012). Dual specificity phosphatases (DUSPs) are a unique class of mitogen-activated kinase phosphatases (MKPs) that remove phosphates from both tyrosine and threonine residues of proteins (Caunt et al., 2008). The removal of phosphate groups renders the protein active or inactive based on the physiological status of the cell. Our laboratory studies the fifth member of DUSP family, DUSP5, which is responsible for dephosphorylating ERK in the nucleus (Kutty et al., 2017). ERK has two isoforms, ERK1 and ERK2, and they become fully active upon phosphorylation at two sites, namely, T202/185 & Y204/187 (Canagarajah, et al., 1997). Therefore, at any given time in a cell, there are four isoforms present, namely: nonphosphorylated ERK (np-ERK), monophosphorylated ERKs (pT-ERK or pY-ERK) or dual phosphorylated ERK (pTpY-ERK). ERKs are phosphorylated in the cytoplasm and the phosphorylated ERKs are transported into the nucleus where they exert their biological effect, such as gene expression. Maintaining pERK concentration in the nucleus is partly dictated by available pools of ERK isoforms for phosphorylation. Contribution to this ERK isoform pool can also arise from dephosphorylation of ERK. DUSP5 is responsible for dephosphorylation of both ERK1 and ERK2 in the nucleus (Kidger et al., 2017). DUSP5 removes phosphate groups from both T and Y residues, thereby inactivating ERK1/2. Therefore, the mechanism of dephosphorylation of one or both these ERK phosphates is important to understanding the cellular biology of ERK. Such a mechanistic understanding would also facilitate inhibitor development, with some inhibitors possibly having therapeutic applications.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1385560/fchbi-03-1385560-HTML/image_m/fchbi-03-1385560-g001.jpg)
Figure 1. Model of the DUSP5-ERK2 mechanism. (A, B) Schematic representations of pERK binding to the DUSP5 phosphatase domain (PD) and ERK binding domain (EBD), leading to the near attack conformation. The corresponding kinetic scheme (C) shows that kcat/Km reflects steps up to formation of the covalent adduct formed with the cysteine thiol nucleophile. Adapted from our prior studies (Talipov et al., 2016; Gupta et al., 2019). (D) Model illustrating the cellular-level impact of the DUSP5 - pERK interaction on biological function.
Developments in the field of genetics and epigenetics have enabled the first wave of genome-wide analyses on people diagnosed with specific disorders. These analyses have provided critical insights, suggesting the involvement of distinct genes and biological processes in the onset of mental, neurological and vascular disorders. Increasing evidence has identified associations in the DUSP family with substance use disorder, cancer and vascular disorders (Ueda et al., 2003; Kuntz-Melcavage et al., 2009; Pramanik et al., 2009; Seleman et al., 2014; An et al., 2021). Our initial interest in investigating DUSP5:ERK2 dephosphorylation mechanism stemmed largely from the identification of a somatic mutation in DUSP5 (S147P) in patients diagnosed with vascular anomalies (Pramanik et al., 2009). The S147P mutation in DUSP5 causes hypoactivity of DUSP5 towards ERK2 (Nayak et al., 2014). To determine the molecular mechanism of DUSP5-mediated dephosphorylation of ERK2, we constructed a model of DUSP5/pT185pY187-ERK2 complex and noted in the X-ray structure of the DUSP5 phosphatase domain two phosphate binding sites, one in the active site where dephosphorylation occurs and the catalytic Cys263 resides, and one adjacent to the active site (Talipov et al., 2016). The presence of two phosphate pockets distinguishes DUSPs from other phosphatases which typically contain only the active site pocket. DUSP5 is also unique from other DUSPs in the presence of a catalytic Cys nucleophile.
Of mechanistic relevance for DUSP5, Protein Tyrosine Phosphatase 1B (PTP1B) also has a second site adjacent to the active site, which was found to be important in studies of vanadium-PTP1B complexes. The H-bonding to the five-coordinate protein adduct between the Cys thiol and the vanadium is stabilized by the many H-bonds to the vanadium and ligands in the first coordination sphere, as well at several H-bonds which extend into the adjacent site, further stabilizing the protein-inhibitor complex. The nature of the interplay between two adjacent sites in DUSP5 and PTP1B (Huyer et al., 1997; Zhang. and Zhang, 1998; Heo et al., 2002; Li et al., 2008; Feng et al., 2021) may provide a model for future DUSP5 inhibitor design, which has not yet been explored with all the phosphatases containing Cys in the active site (Lawrence and Van Etten, 1981; Zhou et al., 1993; Huyer et al., 1997; Zhang et al., 1997; Pannifer et al., 1998; Zhang. and Zhang, 1998; Reiter et al., 2001; Deng et al., 2002; Heo et al., 2002). This mechanism, where a Cys thiol attacks the phosphorylated peptide forming a protein-phosphate covalent adduct, was proposed for DUSP5 in our prior studies (Figure 1).
Our hypothesis in this work is that the two phosphate binding sites will be occupied by the phosphate groups (pT185& pY187) of ERK2. Indeed, molecular modeling studies confirmed this and further suggested the entry of the T185 residue present on the activation loop of ERK2 into the secondary DUSP5 binding site, which disrupts a salt bridge formed across the primary active site, thus facilitating entry of Y187 into the active site for subsequent phosphate hydrolysis (Figure 1). It is known that ERK2 is fully active when it is doubly phosphorylated (Canagarajah et al., 1997) and sequential dephosphorylation occurs with pY dephosphorylated first followed by pT. Whether the dephosphorylation events mediated by DUSP5 also follow a similar preferential sequence of events was not known. The studies performed herein partly address the mechanism and preference for the dephosphorylation reaction of ERK pT or pY residues in the DUSP5 pocket. The DUSP5-catalyzed ERK2 dephosphorylation mechanism, based on our prior studies, is summarized in Figure 1A–C (Talipov et al., 2016; Gupta et al., 2019). The biological role of DUSP5 is illustrated in Figure 1D, where it is shown that there are 2 isoforms of ERK (1 and 2); and, upon phosphorylation at two different sites (tyrosine Y187 and threonine T185), the ERKs become fully active. ERK is phosphorylated in the cytoplasm, crosses the nuclear membrane and exerts its biological effect via gene expression. DUSP5 is responsible for dephosphorylating ERK in the nucleus, controlling p-ERK activity by removing the phosphates from one or both sites. Dephosphorylated ERK 1 and 2 are inactive kinases. Therefore, it is critical to understand the mechanisms of phosphorylation and dephosphorylation of ERKs. Moreover, clinical research from genetic and epigenetics studies have identified associations between DUSP family genes and different mental, neurological and vascular disorders. Specifically, DUSP5’s role in substance use disorder, cancer, autoimmune diseases and vascular development. The mechanistic studies presented herein therefore provide a foundation for a better understanding of these downstream effects of ERK phosphorylation state.
2 Materials and methods
2.1 Protein expression and purification
The Dusp5 phosphatase domain gene, both wild type (WT; DUSP5pdWT; used for kinetic studies) and the RDB version (DUSP5pdRDB; used for NMR studies), was synthesized by Blue Heron (Bothell, WA); and WT protein was expressed and purified as previously described (Nayak et al., 2014; Neumann et al., 2015). The DUSP5pdRDB construct used for the studies herein differs from the original wild type construct due to the addition of a thrombin cleavage site and the mutation of cysteine at position 263 to serine, to avoid catalytic turnover. This construct is identical to that which has been crystalized (2g6z), thereby permitting better comparison of NMR and x-ray data. The DUSP5pdRDB thrombin cleavage site also avoids additional N-terminal amino acids, which affected protein stability and quality of NMR spectra (Supplementary Figure S1).
Dual labeled (13C/15N) DUSP5pdRDB protein with an N-terminal hexa-histidine tag was expressed in BL21(DE3) cells grown in minimal media containing ammonium-15N chloride and 13C6 D-glucose, using a New Brunswick™ BioFlo®/Celligen® 115 Fermentor/Bioreactor (Eppendorf). The vessel media temperature, dissolved oxygen, agitation rate and pH levels were maintained at 25°C, 100% room air oxygen, 400 rpm and pH 7.2, respectively. Protein expression was induced with 0.1 mM IPTG once the media O.D. at 600 nm reached 1.5. The cells were incubated an additional 14 h (h) and harvested by centrifugation. The resulting cell pellets were stored at −80°C. Dual labeled DUSP5pdRDB was isolated using a multi-step purification. The harvested cell pellets were lysed using sonication and the supernatant was passed through a Ni-sepharose Fast Flow resin column (GE Healthcare Life Sciences) to isolate the His-tagged protein. The histidine tag was cleaved with thrombin (1 unit thrombin per 100 µg protein), and a second Ni-sepharose column was used to isolate the protein without the histidine tag. Thrombin removal was accomplished by passage of the His-cleaved protein through a benzamidine column. The resulting sample was dialyzed (20 mM KH2PO4(pH 5.0), 75 mM KCl and 2 mM DTT) and concentrated to 1 mM by ultrafiltration in preparation for NMR spectroscopy analysis. Histidine tag cleavage and protein purification was monitored by SDS PAGE, with stacking and resolving gel compositions of 5% and 10%, respectively. See Supplementary Figure S2. Unlabeled protein was expressed and purified in a like manner, but using LB media rather than minimal media. NMR studies were performed using protein with a serine in the active site, whereas the wild type (WT) protein for enzyme kinetic studies had the cysteine (Cys263).
2.2 Multinuclear NMR studies of tripeptide reaction with DUSPpdWT
2.2.1 Preparation of NMR samples and optimization of DUSP5RDB buffer conditions
Various NMR buffer conditions were explored, to determine conditions which gave maximum stability for NMR studies, and highest quality 2D HSQC spectra. The conditions examined were deuterated Tris buffer and potassium phosphate buffer, with pH values of 5.0 and 6.8. The procedure outlined below is for the 50 mM Tris, 100 mM KCl, and 2 mM DTT buffer samples at pH 5.0 and 6.8. This same procedure was followed for the investigation in buffer containing 20 mM potassium phosphate, 75 mM KCl, and 2 mM DTT.
An 15N labeled DUSP5pdRDB sample with a concentration of 500 µM was thawed on ice from −80°C storage. Samples without glycerol-d8 were prepared by transferring 300 μL of the 15N labeled protein and 30 μL D2O to a 1.5 mL Eppendorf tube. Glycerol-d8 containing samples were prepared in a similar fashion by transferring 300 μL of 15N labeled protein, 30 μL D2O, and 30 μL glycerol-d8 for a final concentration of 10% D2O and 10% glycerol-d8. The samples were spun down in a centrifuge at 14,000 rpm for 2 min to remove any particulate. The supernatant was transferred to a Shigemi NMR tube.
2.2.2 31P Multinuclear NMR studies of tripeptide reaction with DUSPpdWT
Hydrolysis of the pThr-Glu-pTyr peptide was monitored using 31P 1D NMR and a Varian 500 MHz spectrometer equipped with a broad band probe. 1 mM of peptide was incubated with 100 μM of DUSPpdWT protein for 3 min, and spectra were collected at 25oC. See Figure 2 and, for a time course, Supplementary Figure S3.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1385560/fchbi-03-1385560-HTML/image_m/fchbi-03-1385560-g002.jpg)
Figure 2. 31P NMR spectra of the DUSP5-catalyzed hydrolysis of the pT-E-pY tripeptide. (A) The DUSP5-catalyzed hydrolysis reaction. (B) 31P NMR spectrum of 1 mM pT-E-pY tripeptide before adding enzyme, showing signals for phospho-threonine (doublet) and phospho-tyrosine (singlet). (C) 31P NMR spectrum 3 min after addition of 100 μM DUSP5pdWT, showing 100% conversion to orthophosphate. Reactions with lower enzyme concentration showed no clear preference for hydrolysis of pThr or pTyr (Supplementary Figure S3).
2.2.3 2D 1H-15N HSQC spectra of DUSP5pdRDB protein labeled with 15N and 13C
2D 1H-15N HSQC spectra were collected using a 500 MHz Varian NMR spectrometer of a sample of expressed and purified DUSP5pdRDB protein labeled with 15N and 13C. The temperature for the experiments was adjusted to either 25°C, 30°C, or 37°C. After the sample was done running on the NMR spectrometer it was transferred to an oven set to either 25°C, 30°C, or 37°C. The samples were held at their respective temperatures for a total of 72 h 2D 1H-15N HSQC experiments were run on each sample after the completion of each 24 h period until 72 h. Final optimized conditions, based on quality of 2D 1H-15N HSQC spectra, were determined to be 20 mM potassium phosphate (pH 5), 75 mM KCl, 2 mM DTT with 10% glycerol, with NMR experiments run at 25oC. These are the conditions used for the 3D NMR experiments, and for chemical assignments, using the DUSP5pdRGB construct (Supplementary Figure S1).
2.2.4 2D and 3D NMR and chemical shift assignments
Using these optimized conditions, protein chemical shift assignments were made for 15N and 13C labeled and cleaved DUSP5pdRDB protein (summarized in Supplementary Table S1). The multidimensional NMR experiments were performed on a 600 MHz Bruker NMR spectrometer equipped with a cryoprobe, using instrumentation at Medical College of Wisconsin.
The sample was concentrated to 1 mM in a total volume of 0.3 mL. 30µL glycerol-d8 was added to the sample and it was flash frozen and stored in the −80°C freezer until ready to analyze with the NMR spectrometer. 30µL D2O and 0.02% sodium azide were added to the sample. The sample was transferred to a Shigemi NMR tube and a standard suite of 2D and 3D experiments were collected for structure determination: 2D 1H-15N HSQC, 3D HNCO, HN(CO)CA, HN(CO)CACB, CCONH, HN(CA)CO, HNCA, and HNCACB.
Spectral data were processed using NMRpipe, then experiments were viewed in Xeasy and were peak picked. Once all the peaks were picked in all of the experiments, the data were run through Garant, which was able to auto assign roughly 85% of the chemical shifts. Next, the Garant assignments were opened into the HSQC experiment in Xeasy. The assignments were checked and corrected if necessary. This procedure was repeated for the HNCO experiment. The correct peak list of the HNCO experiment were loaded into the HN(CO)CA. The corrected HN(CO)CA list was loaded into the HNCACB, and the HN(CO)CACB experiments. Lastly, the HN(CO)CACB peak list was loaded into the CCONH experiment. The completed peak list was once again opened into the HSQC file for a final check of chemical shifts. These chemical shifts were saved to the prot file, with final assignments provided in Supplementary Table S1. Assignments mapped onto 2D 1H-15N HSQC spectra are in Supplementary Figure S4, S5.
2.2.5 Secondary structure assignments
The chemical shifts that were assigned for the primary structure were re-formatted for the TALOS+ web based secondary structure assignment (National Institute of Diabetes and Digestive and Kidney Diseases, 2014). Once the file was reformatted, the data file was uploaded to the TALOS+ server (Shen et al., 2009). This produced the analysis files. These files were then read through the jRAMA+ web-based java viewer for TALOS+ torsion angle predictions. These predictions were compared to the crystal structure, with pdb code 2G6Z (See Figure 3).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1385560/fchbi-03-1385560-HTML/image_m/fchbi-03-1385560-g003.jpg)
Figure 3. Comparison of crystal structure with NMR-based predictions. (A) Crystal structure of the DUSP5pd-WT protein (pdb code 2G6Z). (B) TALOS+ secondary structure predictions superimposed on the crystal structure. Red is for helix, Blue is for beta sheet and Green is for loop. (C) TALOS+ secondary structure and order parameter (flexibility) predictions. Two highly flexible regions (A, B) are identified, based on low predicted order parameter (S2) values.
2.2.6 DUSP5-tripeptide (pT-E-pY) NMR binding studies
After dialysis, two NMR sample tubes were prepared by adding 275 µL of 564 µM DUSP5pdRDB protein, 30 µL glycerol-d8, 30 µL D2O, and varying concentrations of the Ac-pThr-Glu-pTyr-NH2 (pT-E-pY) peptide (N-terminus acetylated and C-terminus amide) derived from pERK. Peptides were synthesized (Biomatik LLC, Wilmington, Delaware) to greater than 98% purity as determined by HPLC. The concentration of the tripeptide added was 0.250 mM. The Eppendorf tubes were centrifuged at 14,000 rpm for 2 min. After centrifugation the supernatant of each Eppendorf was transferred to each of two Shigemi NMR tubes. 2D 1H-15N HSQC spectra were collected on a 500 MHz Varian NMR spectrometer, at 25°C. Spectra were processed then overlaid using vNMRDraw software (shown in Figure 4). Exchange-broadened active site residues (identified in Figure 4) as well as residues with low predicted order parameters (low S2; from Figure 3) are indicated in our previously reported model of the DUSP5-ERK2 complex (Figure 5).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1385560/fchbi-03-1385560-HTML/image_m/fchbi-03-1385560-g004.jpg)
Figure 4. 2D 1H-15N HSQC spectrum of DUSP5pdRDB ± tripeptide. 1H-15N HSQC spectrum for DUSP5pdWT protein alone (blue) or in the presence of the pT-E-pY tripeptide (red). Exchange broadening due to peptide binding occurred for many residues, including active site residues A265, S268 and I267. Full chemical shift assignments are provided in Supplementary Figure S4, S5, Supplementary Table S1).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1385560/fchbi-03-1385560-HTML/image_m/fchbi-03-1385560-g005.jpg)
Figure 5. Model of the DUSP5-ERK2 complex. Results of NMR chemical shift analysis (Figures 3, 4) identifying exchange broadened active site residues (A) and dynamic regions with low predicted S2 values (B) have been mapped onto our previously reported model of the DUSP5-ERK2 complex.
2.3 Steady state kinetic analysis of phosphate hydrolysis
2.3.1 Steady state kinetic analysis of phosphate hydrolysis using a continuous Biomol Green assay
Peptide stock solutions were kept in buffer. Steady state kinetics were performed as described for the peptide inhibition studies, but monitoring actual dephosphorylation of the three tripeptides by DUSP5pdWT: pT-E-Y, T-E-pY, pT-E-pY. 1 mM of each tripeptide was incubated with 0.18 mM DUSP5pdWT enzyme at 25oC, and at timed intervals samples were analyzed for inorganic phosphate using the Biomol Green reagent (at 620 nm) (Supplementary Figure S6). Lack of protein aggregation at these concentrations is evidenced by consistently narrow linewidths in 2D HSQC spectra at various enzyme concentrations (i.e., minimal T2 broadening in Figs. 4, S4 and S5) as well as our previous data on dynamic light scattering of DUSP5 (Gupta et al., 2019). Kinetic studies were also repeated with varied tripeptide substrate concentration, to generate Michaelis-Menten curves. Results of fitting to Equation (1) are shown in Table 1.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1385560/fchbi-03-1385560-HTML/image_m/fchbi-03-1385560-t001.jpg)
Table 1. DUSP5pdWT kinetic parameters (±SE) estimated from Michaelis-Menten fits (Eq. 1) of substrate velocity curves with the indicated tripeptides as substrates. Velocities were determined from the rates of inorganic phosphate generation, monitored with Biomol green reagent.
2.3.2 Steady state kinetic analysis of phosphate hydrolysis and inhibition by tripeptides
Recombinant DUSP5 PD (DUSP5pdWT) phosphatase activity was assayed as previously described, with minor modifications (Nayak et al., 2014; Neumann et al., 2015). All assays were performed in clear bottom 96-well microplates (Greiner 655001) with a total assay volume of 210 μL, using a SpectraMax M5 microplate reader (Molecular Devices). Assay buffer contained 100 mM Tris, 100 mM NaCl, 5 mM MgCl2·6H2O with 1 mM dithiothreitol (DTT) at pH 7.5. p-nitrophenol phosphate (pNPP, Sigma Aldrich) at 5 mM, a concentration near the established Km, was used for substrate. DUSP5 PD-mediated hydrolysis of pNPP to p-nitrophenolate, a chromogenic product which absorbs at 405 nm with a molar extinction coefficient of 18,000 M−1 cm−1, was monitored over time using the Spectramax M5 microplate reader. Wells that contained assay buffer with pNPP were used as blanks. Stock solutions of peptides were kept in buffer. Serial dilutions of each stock solution were prepared to generate a range of inhibitor concentrations for the steady state inhibition studies. Absorbance values were converted to p-nitrophenolate concentration using the product extinction coefficient and a micro well path length of 0.545 cm DUSP5 PD (DUSP5pdWT) enzyme progress curves over a range of substrate and peptide inhibitor concentrations were fitted globally using GraphPad Prism 6 to Equation (2) for steady-state competitive inhibition:
where
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1385560/fchbi-03-1385560-HTML/image_m/fchbi-03-1385560-t002.jpg)
Table 2. Kinetic parameters for inhibition of DUSP5pdWT hydrolysis of para-nitrophenol phosphate (pNPP) in the presence of pERK activation loop tripeptide mimetics as inhibitors, having different phosphorylation states. Kinetic parameters estimated from global non-linear regression competitive inhibition model fits (Michaelis-Menten) to Eq. 2.
2.3.3 Steady state kinetic analysis of phosphate hydrolysis and inhibition by vanadate and phosphate
All assays were carried out using pNPP as a substrate in the presence of DTT using a microplate reader and a total volume of 100 μL. Each assay contained an enzyme concentration of approximately 0.072 mg/mL (appropriately 3.13 μM DUSP5pdWT) which was left in a stock solution at about 1 h at 4°C before use. Each assay consisted of a reaction mixture in 100 mM tris-hydrochloride buffer containing 100 mM NaCl and 5.0 mM magnesium chloride (MgCl2) at pH 7.5. The concentration of pNPP was varied from 2 to 25 mM (2, 4, 8, 15 and 25 mM) which was determined after measuring a Km value of 8.02 mM using these conditions. The kinetic parameters were measured in the presence of 0.045 mM DTT in the assay solution and the assay was started with the addition of pNPP.
The inhibition by vanadate was carried out using inhibitor concentrations of 100, 200, 400, 600 and 800 nM under the assay conditions described above and in an assay solution containing 0.045 mM DTT. For phosphate inhibition studies, the phosphate inhibitor concentrations used were 0.5; 1; 2; 5; 10 mM. Both sets of inhibition studies are shown in Supplementary Figure S8, with both phosphate and vanadate behaving as competitive inhibitors. The kinetic runs were repeated three times and the data were initial processed using excel and the global kinetics analysis was performed using GraphPad Prism 10.1.2. The errors provided in the analysis were calculated by the program.
2.4 Docking and molecular dynamics
2.4.1 Docking and molecular dynamics
Three series of molecular dynamics simulations with varied structural features of the molecules were conducted: neutral cysteine (Cys263) in the active site of DUSP5 and phosphorylated residues of the tripeptides in dianionic PO32- form (Series 1), deprotonated Cys263 and phosphorylated residues of the tripeptides in dianionic PO32- form (Series 2), deprotonated Cys263 and phosphorylated residues of the tripeptides in monoanionic PO3H− form (Series 3). Results are summarized in Table 4 and Supplementary Figure S9, S10.
Calculations were done using methods described previously (Talipov et al., 2016; Gupta et al., 2019). Briefly, MD simulations of DUSP5pdWT were performed with the secondary binding site and active/catalytic binding sites loaded or unloaded with mono-phosphorylated ERK2 (i.e., pT-ERK2 or pY-ERK2) tripeptides described above (T-E-pY; pT-E-Y; also T-E-Y), while the remaining site was vacant and conformationally flexible in both cases. These simulations of the PD with the loaded secondary binding site or active site serve to probe the structural rearrangement of the active and secondary binding sites as compared with the original X-ray structure of PD domain. MD simulations were performed as before (Gupta et al., 2019) using the Amber MD software (Case et al., 2005) with the ff14SB force field with a 10.0 Å force cutoff and the Particle Mesh Ewald algorithm (Duan et al., 2003) to treat long-range electrostatic interactions. A triplicate 8 microsecond MD simulation (24 μs in total) was performed with each of the tripeptides: pT-E-pY, pT-E-Y, T-E-pY, and T-E-Y. Solvent was modeled explicitly with the TIP3P water model (Jorgensen et al., 1983). For each simulation, the initial structure (prepared as described below) was subjected to 500 steepest descent (SD) + 500 conjugate gradient (CG) minimization, in which protein structure was frozen, and subsequent 1000 SD + 1500 CG full minimization steps performed. The resulting structure was subjected to a warm-up NVT simulation with temperature being increased from 0 to 300K during 40 ps and subsequent 3-ns NPT equilibration. The sets of coordinates and velocities after each nanosecond of simulation were used as a starting point for three independent production run simulations. Production runs were performed using an NPT ensemble, with the time constant for heat bath coupling for the system equal to 10 ps and a Monte Carlo barostat. For each production run simulation, it was confirmed that temperature, mass density, and potential energy were nearly constant and did not show any noticeable drift.
A docking study with molecular dynamics, as described above, was performed also using the fully phosphorylated peptide (pT-E-pY), and a population analysis was used to assess which of the two phosphate binding pockets was preferentially occupied by pT or pY, and at what energy (Figure 6).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1385560/fchbi-03-1385560-HTML/image_m/fchbi-03-1385560-g006.jpg)
Figure 6. Molecular dynamics-docking population analysis. (A) Structural model of DUSP5 with the pT-E-pY (pThr-Glu-pTyr) peptide docked into the active site. (B) Relative population and total potential energy of binding modes with pY (blue) vs pT (gold) in the active site pocket, indicating that after MD, the preferred binding mode has pY in the active site pocket rather than pT. (C) A plot of distance from the active site serine to phosphate of pY vs phosphate of pT, indicating multiple low energy binding modes and preferential binding of pY to the active site pocket.
The starting structure for simulations was the 2.7 Å resolution X-ray structure (PDB code 2G6Z) (Jeong et al., 2006), with Ser263 typically replaced with Cys, which is the WT (catalytically active) form. Structures were prepared using UCSF Chimera (Pettersen et al., 2004). Statistical analysis of the conformational dynamics was performed using R software version 3.3.1. For the cluster analysis, a hierarchical cluster analysis technique was used for the combined set of MD trajectories.
3 Results
3.1 Protein expression and purification for NMR studies
The DUSP5pdRDB construct (Supplementary Figure S1) was created with a thrombin cleavage site and the active site Cysteine 263 mutated to Serine to avoid catalytic turnover and oxidation during NMR experiments. Purification with a Nickel NTA column yielded pure protein for both steady state kinetic (DUSP5pdWT; with Cysteine 263) and NMR (DUSP5pdRDB; with Serine 263) studies (Supplementary Figure S2). The 15N-13C double labeled protein was well-behaved for >72 h, enabling collection of the suite of 3D NMR experiments needed to determine full chemical shift assignments (HNCO, HN(CO)CA, HN(CO)CACB, CCONH, HN(CA)CO, HNCA, and HNCACB); and, it gave well-dispersed 1H-15N HSQC spectra, indicative of a well-folded and stable protein (Supplementary Figure S4, S5). Chemical shift assignments were obtained and are summarized in Supplementary Table S1.
3.1.1 31P NMR-based characterization of phosphate hydrolysis of the pT-E-pY tripeptide
DUSP5 is a dual specificity phosphatase that cleaves both phospho-threonine (pT) and phospho-tyrosine (pY) residues, and it does this via a two-pocket gated mechanism that we have previously described (Figure 1). The peptide from the ERK activation loop that is dephosphorylated by DUSP5 actually contains proximal phosphates on both of these residues, separated by a glutamate: pT-E-pY. Studies presented herein characterize the binding and reaction of this tripeptide, derived from ERK. This pT-E-pY tripeptide is a substrate for DUSPpdWT with both phosphates cleaved 100%, based on 31P NMR analysis of the reaction progress (Figure 2). This result demonstrates that both the phosphate on the pY and the pT residues can be hydrolyzed by DUSPpdWT.
Hydrolysis was also examined at 10-fold lower enzyme concentration (Supplementary Figure S3) with the objective of observing a preference for either the pT or the pY residue. At the first timepoint after 3 min incubation, about 50% conversion was observed. Thus, about 25% each of the pT and pY peptides remained even after prolonged incubation. While both phosphates (pY and pT) are hydrolyzed efficiently in 31P NMR studies at 10-fold higher enzyme concentration, the fact that no more hydrolysis is occurring at lower enzyme concentration suggests either the enzyme is decomposed, or the enzyme is inhibited. In addition, the question remains as to whether there is differential affinity or binding order for pT or pY in each of the two binding pockets of DUSP5 (Figure 1), and how the enzyme catalyzed reaction occurs under steady state conditions.
3.1.2 Protein NMR-based secondary structure prediction and pT-E-pY tripeptide binding study
Based on DUSP5 chemical shift values, secondary structure was predicted, with these predictions mapped onto the known DUSP5 x-ray crystal structure in Figure 3B. NMR data confirm that the solution structure is the same or at least consistent with the crystal structure, with perhaps more flexibility (loop-like structure) observed in some regions at the end of helices (Figure 3). Furthermore, the NMR data yielded a prediction of high flexibility (low S2 values) in regions 215-220 and 232-237, which we term “Flexible Regions A and B.” Extensive exchange broadening in 2D 1H-15N HSQC spectra (Figure 4) is consistent with protein conformational mobility in general. Mapping the Flexible Regions onto our modeled DUSP5-ERK complex reveals that these two loop regions define the binding interface with ERK (Figure 5), perhaps positioning the DUSP5 phosphatase domain relative to ERK so that the ERK activation loop (containing the pT-E-pY peptide) fits into the DUSP5 active site pocket for hydrolysis to occur.
Chemical shift perturbation studies, comparing 2D 1H-15N HSQC spectra before and after binding of the pT-E-pY tripeptide, indicate there is significant exchange broadening of cross peaks (Figure 4), including primary active site residues Ala265, Ile267 and Ser268. This suggests that the DUSP5-(pT-E-pY) complex is in a dynamic rather than rigid state, consistent with our previously published mechanism (Figure 1). These residues are shown on the DUSP5 phosphatase domain structure in Figure 6A and are close in space to the catalytic Cys163 residue in the primary active site pocket. Their flexibility on the NMR exchange timescale (typically msec) is consistent with motions and flexibility in the primary active site pocket that may be necessary to transfer pY or pT from binding in one pocket to the next.
3.1.3 Steady state kinetic analysis of phosphate hydrolysis of phosphorylated tripeptides
Steady state kinetic studies were performed with the tripeptides acting as substrates (vs inhibitors), with progress curves and Michaelis-Menten plots shown in Supplementary Figure S6. Fitted Michaelis-Menten kinetic profiles yielded a relative activity as a substrate, expressed as kcat/Km, of: pT-E-pY > T-E-pY > pT-E-Y (Table 1). This clear trend is only with the kinetic parameter kcat/Km, the “catalytic efficiency,” that reflects all steps from (and including) initial binding up to the first irreversible step (Figure 1C), which is most likely the (initial) phosphate transfer that creates the phosphor-enzyme adduct (as also reported by Brandão, et al. (2010), for PTP1B). What is not clear from these results is whether the pT-E-pY peptide binds sequentially in one site then the next, and—even though there is higher overall affinity for pY over pT—what the preference is in each of the two binding pockets. In other words, the binding mechanism needs further exploration, which we performed below using kinetic studies with pNPP as substrate and phosphorylated peptides and with molecular modeling.
3.1.4 Steady state kinetic analysis of pNPP hydrolysis inhibited by phosphorylated tripeptides
Steady state inhibition studies were performed with pNPP as substrate and the fully phosphorylated, nonphosphorylated, and both mono-phosphorylated T-E-Y tripeptides as inhibitors. The absorbances at 405 nm, measuring the dephosphorylation of pNPP, was analyzed with global fit analysis fitting to the equation for competitive inhibition (Supplementary Figure S7). The resulting kinetic parameters (Km, kcat, Ki and kcat/Km) summarized in Table 2 indicate an overall binding affinity preference and relative phospho-tripeptide inhibition potency of: pT-E-pY > T-E-pY > pT-E-Y > T-E-Y. The diphosphorylated peptide is a roughly 3x more potent binder than the pY monophosphoryled peptide, which is in turn 3x more potent than the pT monophosphorylated peptide. But the greatest affinity boost comes from the addition of the initial phosphate group, since the pT monophosphoryleted peptide binds 30x stronger than nonphosphorylated peptide.
These data demonstrate that the phosphorylated peptides are inhibitory. While both pNPP and the phosphorylated peptide inhibitor can hydrolyze during the course of the assay, the rate of hydrolysis of pNPP is much faster than of the phosphorylated peptide. This can be rationalized because the pNPP concentration range is 2–20 mM whereas the peptides are present at much less than 1 mM, and the amount of the latter hydrolyzed is not contributing significantly to the overall level and rate of phosphate hydrolysis. Hence incubation of phosphorylated peptides with DUSP5pdWT will inhibit this enzyme and the observed 31P NMR results with low enzyme concentration where the reaction stops at 50% conversion can be explained at least in part by incubation of the phosphorylated peptide. Furthermore, these data show that the pY binds tighter than pT, consistent with the possibility that the initial binding of the diphosphorylated peptide will be through binding of pY. This hypothesis is explored further in molecule docking and dynamics experiments described below.
3.1.5 Steady state kinetic analysis of pNPP hydrolysis inhibited by vanadate and phosphate
Since DUSP5 has an active site cysteine that is thought to act as a nucleophile (Figure 1), forming a Cys-phospho intermediate analogous to PTP1B, it might be expected DUSP5 is potently inhibited by vanadate, since Brandão, et al. (2009 and Brandao et al., 2010) had previously reported such a potent inhibitory adduct formed between vanadate and the cysteine thiol of PTP1B. Vanadate readily forms trigonal bipyramidal structures and are widely accepted as mimicking transition state or high energy intermediate structures for phosphate ester hydrolysis (Crans et al., 2014; Crans, 2015; McLauchlan et al., 2015). The structural advantage that vanadate has over other analogs including tungstate (WO42-), molybdate (MoO42-), magnesium fluoride (MgF3−), and nitrate (NO3−) is that it has the same overall charge as phosphate. Herein we explored the inhibition by vanadate and compared it with inhibition of phosphate, as our lab has done previously (Crans, 2015; McLauchlan et al., 2015). However, because of the redox properties of vanadate (Crans et al., 2004), it can be challenging to use in the presence of a reducing agent (DTT) that is needed to maintain the stability of protein phosphatase, with a readily oxidized cysteine in the active site (Crans et al., 2014; Crans, 2015; McLauchlan et al., 2015). We carried out preliminary studies with DUSP5 using varying concentrations of DTT. Tris forms a weak complex with vanadate, and therefore was compared to Hepes, which is a buffer with similar pKa value that interacts minimally with vanadate (Crans et al., 1989). Preliminary studies in both buffers showed that some variations in Km values for pNPP and Ki values for vanadate were observed. These preliminary studies allowed us to choose Tris buffer with minimum DTT concentration for studies with the DUSP5pdWT enzyme for comparison with the kinetic analysis carried out in this work.
Using the conditions determined by the pNPP Km of 8.02 mM that was measured and calculated using global kinetics analysis, with vanadate at concentrations at 100, 200, 400, 600 and 800 nM, a Ki value of 638 nM was obtained for vanadate (Supplementary Figure S8). Considering that the active site of DUSP5 contains a cysteine residue and it has a secondary binding site which does not, DUSP5 is distinct from than other cysteine containing phosphatases. However, as shown by our data, we find that the inhibition is competitive and the Ki value is very potent, submicromolar. A summary of Ki values surveyed from the literature is listed in Table 3, and demonstrates the variations in Ki values reported for vanadate, and for comparison phosphate, considering the differences in phosphatases with cysteine in the active site, substrate, and pH values in the assay (Lawrence and Van Etten, 1981; Zhou et al., 1993; Huyer et al., 1997; Zhang et al., 1997; Zhang. and Zhang, 1998; Deng et al., 2002; Heo et al., 2002; Lu et al., 2012). Except for one bovine liver enzyme (Lawrence and Van Etten, 1981) these phosphatases all have micromolar or submicromolar Ki values. The most well studied vanadate-phosphatase system with a cysteine residue in the active site is protein phosphatase 1B (PTP1B) (Brandão et al., 2009; Brandão, et al., 2010; Crans, 2015). There are multiple X-ray structures of the vanadate-PTP1B complex showing a range of different complexes. One complex in particular has part of the peptide tyrosine-phosphate substrate bound to the PTP1B. The substrate is a hexapeptide with a sequence DADEYL which is a fragment of the epidermal growth factor receptor, which, in phosphorylated form (DADEpYL) exhibits strong affinity for PTP1B (Pedersen et al., 2004). A previous crystallographic study shows a Michaelis complex between a catalytically inactive C215S mutant of PTP1B and the phosphorylated peptide (Jia et al., 1995) but the second transition state of this protein with vanadate (Pedersen et al., 2004; Brandão, et al., 2010) has lost the tyrosine peptide. The crystal structure of the C215S mutant of PTP1B and the phosphorylated peptide shows the peptide extending from the active site to the adjacent site which provides a number additional stabilizing H-bonds and as such a system to which there is more direct comparison to DUSP5 than the other Cys-containing phosphatases.
For comparison, the inhibitory effect of phosphate on DUSP5 was also determined (Supplementary Figure S8B). Phosphate also exhibits a competitive inhibition pattern with DUSP5 and yields a Ki of 17.3 mM, similar to those reported previously (Lawrence and Van Etten, 1981; Zhang et al., 1997; Zhang. and Zhang, 1998). The considerably lower inhibitory effects of phosphate (5 orders of magnitude higher Ki) compared to vanadate (Ki = 638 nM) is due to the fact that phosphate is not able to form the stable pentacoordinate complex in the trigonal bipyramidal structure that vanadate can adopt, resembling the transition state or high energy intermediate as reported for other protein phosphatases (McLauchlan et al., 2015). Our studies demonstrate that DUSP5 is subjected to the same effects during its catalytic cycle, and in this regard resemble PTP1B.
3.2 Molecular docking and dynamics of tripeptide binding
3.2.1 Molecular docking and dynamics of tripeptide binding
Binding of the pT-E-pY tripeptide, and its monophosphorylated analogs, was explored using docking and molecular dynamics to assess preference for the two DUSP5 phosphate binding pockets (catalytic and the secondary site; Figure 1) and to verify if the pY residue binds tighter than the pT residue. A 24-microsecond molecular dynamics (MD) simulation of the DUSP5 complex with the pT-E-pY tripeptide was performed, and which showed that (a) the complex is stable and maintains pY in the catalytic site and pT in the secondary site (Figure 6A), (b) there is preferential binding of the pY group of pT-E-pY before pT (Figure 6B), (c) distance between the active site residue (Cys263) is shorter to pY than to pT (Figure 6C), and (d) the pY and pT groups bind in the catalytic (active) and the secondary binding site, respectively, with binding of pT in the secondary pocket permitting binding of pY in the active site, consistent with our published mechanism (Figure 1; Talipov et al., 2016). The results from the accelerated sampling of free energy profile for the complex formation (replica exchange umbrella sampling technique) are also consistent with the catalytic efficiency measurements: pT-E-pY (∆G = −6.8 kcal/mol) > T-E-pY (−6.5 kcal/mol) > pT-E-Y (−6.3 kcal/mol) > T-E-Y (−5.9 kcal/mol).
Our calculations show that the active site (containing the catalytic residue Cys263) is preferred over the secondary active site for binding phosphorylated residues. Tripeptides pT-E-Y and T-E-pY bind to the active site (∼5 Å) by their single phosphate groups in more than 50% of frames in all series in the MD simulation. pY is more likely to be involved in this interaction than pT. Numbers of frames with corresponding interactions for T-E-pY and pT-E-Y are about 80–90% and 50–60% respectively. Among all the series of simulations, the strongest interactions are observed in Series 3 with deprotonated Cys263 and the phosphate groups in monoanionic PO3H− form; and, Series 2 with deprotonated Cys263 and phosphate groups in dianionic PO32- form showed poorest binding between modeled fragments, due to high negative charge of the ligands (Table 4; Supplementary Figure S9, S10).
In the case of a tripeptide with two phosphate groups (pT-E-pY), the main active site is occupied by the pY residue most of the time. The most illustrious results in this connection are from Series 3: the distribution between frames with Cys-pY and Cys-pT interactions are about 90% and 15% respectively (in 5% of cases both phosphorylated groups are located near Cys). According to the Boltzmann distribution, this ratio between conformers corresponds to 1.3 kcal/mol of difference in free energy. Calculations from Series 2 show that two dianions make the ligand too negative for interaction with Cys-. In this case the active site in most of the frames is empty, while the occupation of the secondary binding site is considerably higher than in other sets of calculations. In simulations involving neutral cysteine, pT did not show interactions with the main active site, while pY is bound proximal to Cys in 40% of frames.
Since the single phosphate groups of pT-E-Y and T-E-pY tripeptides are tending to occupy the main active site, the corresponding secondary binding site is mostly empty, with the average distance between a center of mass of arginines ARG213-214 forming the second active site and phosphorylated residues of the tripeptides being about 10 Å (Supplementary Figure S10). The other reason for these long distances is that in most of the frames these arginines are not close to each other and do not form a tight active site. Residues interacting with only one arginine are located some distance from the center of mass.
4 Conclusion
In summary, based on NMR binding studies and the observation of competitive inhibition, the ERK tripeptides–in all phosphorylation states—appear to bind preferentially in the expected primary active site pocket, comprised of residues that include Cys263 (nucleophilic catalyst), Ser268, A265 and I267 (Figure 5). But the observation of significant NMR exchange broadening (Figure 4) suggests the DUSP5 protein undergoes dynamic motions on an intermediate timescale (probably msec), likely due to local conformational changes indicative of a somewhat dynamic protein capable of induced fit binding. This mobility is especially pronounced in loop regions we have defined as Flexible Regions A and B (Figure 3C). Such flexibility may be needed to accommodate binding of the ERK activation loop tripeptide in one of two potential catalytically competent orientations, with pY or pT in the active site pocket (with pY being favored). Surface loops on the DUSP5 phosphatase domain that interface with ERK2 appear especially flexible, based on low predicted S2 values, consistent with a flexible ERK-DUSP5 binding interface, to go along with a malleable DUSP5 active site.
The inhibition studies with vanadate demonstrate that the transition states and catalytic cycle of DUSP5 is likely to be affected similarly to that of other Cys-containing phosphatases, such as PTP1B. We therefore explored the interaction of vanadate with DUSP5 and compared with phosphate. Although DUSP5 is unique in that it has two phosphate binding pockets, an active site and a secondary binding pocket, one for the phosphorylated tyrosine and one for the phosphorylated threonine. These sites likely flex upon peptide binding, with associated conformational changes. The crystal structure of the phosphotyrosine complex of PTP1B has provided a snapshot of its catalytic cycle; but, that phosphatase has only one active site. DUSP5’s substrate extends into an adjacent binding pocket. These interesting similarities with PTP1B will be explored further in future work.
Based on docking and molecular dynamics (MD) studies, the phosphorylated residues (pY; pT) of the ERK activation loop tripeptides can readily bind to the active site of DUSP5pdWT, with a preference for the pY residue binding in the active site pocket relative to the secondary binding site. Still, the MD simulations indicate that both binding orientations are possible, even if pY in the main active site pocket is preferred. This preference for pY is consistent with the higher affinity for the T-E-pY peptide relative to the pT-E-Y peptide in steady state kinetic inhibition studies; and, with the observation that the T-E-pY peptide is more readily hydrolyzed based on the steady state kcat/Km values. Thus, pY binding in the active site pocket, proximal to the catalytic Cys163 residue, is more favorable. But, pT can occupy this site as well and is also readily hydrolyzed, with only 3-fold lower kcat/Km relative to pY. Importantly, DUSP5 is behaving as other Cys-containing protein phosphatases, which are all potently inhibited by vanadate, due to vanadate’s ability to form a stable five-coordinate transition state protein complex.
Future studies will be directed towards further characterizing the structural and dynamics changes associated with DUSP5 binding to ERK and ERK-derived phosphorylated peptides, enabled by the presentation herein of full NMR chemical shift assignments. These studies will enable the design of inhibitors as mechanistic probes, and potential drug leads. Potential candidate inhibitors could include transition state analogs such as vanadate or derivatives obtained from information gained in studies using a vanadate-containing potent transition state analog.
The literature has shown the critical role that DUSP5 plays in disorders that include vascular disease, cancer, autoimmune disease and substance use disorder (Figure 1D). In particular, increases or decreases in DUSP5 mRNA expression were observed in different brain structures involved in reward circuits like the prefrontal cortex and nucleus accumbens, respectively. Moreover, studies suggest that DUSP5 has a role in drug exposure as well as synaptic plasticity during drug-seeking behavior (relapse) (Kuntz-Melcavage et al., 2009; Seleman et al., 2014; An et al., 2021). Future studies in our laboratories will be directed to understanding these biological processes underlying or mediated by DUSP5, and their biological relevance in the context of substance use disorder.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding authors.
Author contributions
AI: Investigation, Methodology, Writing–original draft, Writing–review and editing. NS: Investigation, Methodology, Writing–original draft, Writing–review and editing. RB: Formal Analysis, Investigation, Methodology, Writing–original draft, Writing–review and editing. RS: Investigation, Writing–original draft, Writing–review and editing. AG: Formal Analysis, Writing–original draft, Writing–review and editing. ED: Investigation, Methodology, Writing–original draft, Writing–review and editing. MT: Conceptualization, Formal Analysis, Investigation, Methodology, Writing–original draft, Writing–review and editing. CG-K: Conceptualization, Project administration, Writing–original draft, Writing–review and editing. DC: Conceptualization, Formal Analysis, Project administration, Writing–original draft, Writing–review and editing, Supervision. RR: Conceptualization, Formal Analysis, Writing–original draft, Writing–review and editing, Project administration. DS: Conceptualization, Formal Analysis, Funding acquisition, Investigation, Methodology, Supervision, Writing–original draft, Writing–review and editing.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This work was supported by grant R15GM118304 from the National Institute of General Medical Sciences to DSS and an Institutional Development Award (IDeA) from the National Institute of General Medical Sciences of the National Institutes of Health under grant number P20GM103451 to RS and MRT. RR and AG were partly supported by funds from R33HL154254 and programmatic funds provided by the Department of Pediatrics, Medical College of Wisconsin and Children’s Research Institute, Milwaukee WI.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The author(s) declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fchbi.2024.1385560/full#supplementary-material
References
An, N., Bassil, K., Al Jowf, G. I., Steinbusch, H. W. M., Rothermel, M., de Nijs, L., et al. (2021). Dual-specificity phosphatases in mental and neurological disorders. Prog. Neurobiol. 198, 101906. doi:10.1016/j.pneurobio.2020.101906
Brandão, T. A. S., Hengge, A. C., and Johnson, S. J. (2010). Insights into the reaction of protein-tyrosine phosphatase 1B. J. Biol. Chem. 285 (21), 15874–15883. doi:10.1074/jbc.M109.066951
Brandão, T. A. S., Robinson, H., Johnson, S. J., and Hengge, A. C. (2009). Impaired acid catalysis by mutation of a protein loop hinge residue in a YopH mutant revealed by crystal structures. J. Am. Chem. Soc. 131 (2), 778–786. doi:10.1021/ja807418b
Canagarajah, B. J., Khokhlatchev, A., Cobb, M. H., and Goldsmith, E. J. (1997). Activation mechanism of the MAP kinase ERK2 by dual phosphorylation. Cell. 90 (5), 859–869. doi:10.1016/s0092-8674(00)80351-7
Cargnello, M., and Roux, P. P. (2011). Activation and function of the MAPKs and their substrates, the MAPK-activated protein kinases. Microbiol. Mol. Biol. Rev. 75 (1), 50–83. doi:10.1128/mmbr.00031-10
Case, D. A., Cheatham, T. E., Darden, T., Gohlke, H., Luo, R., Merz, K. M., et al. (2005). The Amber biomolecular simulation programs. J. Comput. Chem. 26 (16), 1668–1688. doi:10.1002/jcc.20290
Caunt, C. J., Armstrong, S. C., Rivers, C. A., Norman, M. L., and McArdle, C. A. (2008). Spatiotemporal regulation of ERK2 by dual specificity phosphatases. J. Biol. Chem. 283 (39), 26612–26623. doi:10.1074/jbc.m801500200
Caunt, C. J., and Keyse, S. M. (2012). Dual-specificity MAP kinase phosphatases (MKPs). FEBS J. 280 (2), 489–504. doi:10.1111/j.1742-4658.2012.08716.x
Crans, D. C. (2015). Antidiabetic, chemical, and physical properties of organic vanadates as presumed transition-state inhibitors for phosphatases. J. Org. Chem. 80 (24), 11899–11915. doi:10.1021/acs.joc.5b02229
Crans, D. C., Bunch, R. L., and Theisen, L. A. (1989). Interaction of trace levels of vanadium(IV) and vanadium(V) in biological systems. J. Am. Chem. Soc. 111 (19), 7597–7607. doi:10.1021/ja00201a049
Crans, D. C., Smee, J. J., Gaidamauskas, E., and Yang, L. (2004). The chemistry and biochemistry of vanadium and the biological activities exerted by vanadium compounds. Chem. Rev. 104 (2), 849–902. doi:10.1021/cr020607t
Crans, D. C., Tarlton, M. L., and McLauchlan, C. C. (2014). Trigonal bipyramidal or square pyramidal coordination geometry? Investigating the most potent geometry for vanadium phosphatase inhibitors. Eur. J. Inorg. Chem., 2014(27), 4450–4468. doi:10.1002/ejic.201402306
Deng, H., Callender, R., Huang, Z., and Zhang, Z. Y. (2002). Is the PTPase−Vanadate complex a true transition state analogue? Biochemistry 41 (18), 5865–5872. doi:10.1021/bi016097z
Duan, Y., Wu, C., Chowdhury, S., Lee, M. C., Xiong, G., Zhang, W., et al. (2003). A point-charge force field for molecular mechanics simulations of proteins based on condensed-phase quantum mechanical calculations. J. Comput. Chem. 24 (16), 1999–2012. doi:10.1002/jcc.10349
Feng, B., Dong, Y., Shang, B., Zhang, B., Crans, D. C., and Yang, X. (2021). Convergent protein phosphatase inhibitor design for PTP1B and TCPTP: exchangeable vanadium coordination complexes on graphene quantum dots. Adv. Funct. Mater. 32 (5). doi:10.1002/adfm.202108645
Gupta, A., Brahmbhatt, J., Raulia, S., Bodnar, C., Bodnar, N., Bongard, R., et al. (2019). Role of conserved histidine and serine in the HCXXXXXRS motif of human dual-specificity phosphatase 5. J. Chem. Inf. Model. 59 (4), 1563–1574. doi:10.1021/acs.jcim.8b00919
Heo, Y. S., Ryu, J. M., Park, S. M., Park, J. H., Lee, H. C., Hwang, K. Y., et al. (2002). Structural basis for inhibition of protein tyrosine phosphatases by Keggin compounds phosphomolybdate and phosphotungstate. Exp. Mol. Med. 34 (3), 211–223. doi:10.1038/emm.2002.30
Huyer, G., Liu, S., Kelly, J., Moffat, J., Payette, P., Kennedy, B., et al. (1997). Mechanism of inhibition of protein-tyrosine phosphatases by vanadate and pervanadate. J. Biol. Chem. 272 (2), 843–851. doi:10.1074/jbc.272.2.843
Jeong, D. G., Cho, Y. H., Yoon, T.-S., Kim, J. H., Ryu, S. E., and Kim, S. J. (2006). Crystal structure of the catalytic domain of human DUSP5, a dual specificity MAP kinase protein phosphatase. Proteins 66, 253–258. doi:10.1002/prot.21224
Jia, Z., Barford, D., Flint, A., and Tonks, N. (1995). Structural basis for phosphotyrosine peptide recognition by protein tyrosine phosphatase 1B. Science 268 (5218), 1754–1758. doi:10.1126/science.7540771
Jorgensen, W. L., Chandrasekhar, J., Madura, J. D., Impey, R. W., and Klein, M. L. (1983). Comparison of simple potential functions for simulating liquid water. J. Chem. Phys. 79 (2), 926–935. doi:10.1063/1.445869
Kidger, A. M., Rushworth, L. K., Stellzig, J., Davidson, J., Bryant, C. J., Bayley, C., et al. (2017). Dual-specificity phosphatase 5 controls the localized inhibition, propagation, and transforming potential of ERK signaling. Proc. Natl. Acad. Sci. U. S. A. 114 (3), E317. doi:10.1073/pnas.1614684114
Kuntz-Melcavage, K. L., Brucklacher, R. M., Grigson, P. S., Freeman, W. M., and Vrana, K. E. (2009). Gene expression changes following extinction testing in a heroin behavioral incubation model. BMC Neurosci. 10, 95. doi:10.1186/1471-2202-10-95
Kutty, R. G., Talipov, M. R., Bongard, R. D., Jones, R. A., Sweeney, N. L., Sem, D. S., et al. (2017). Dual specificity phosphatase 5-substrate interaction: a mechanistic perspective. Compr. Physiol. 7 (4), 1449–1461. doi:10.1002/cphy.c170007
Lawrence, G. L., and Van Etten, R. L. (1981). The low-molecular-weight acid phosphatase from bovine liver: isolation, amino acid composition, and chemical modification studies. Archives Biochem. Biophysics 206 (1), 122–131. doi:10.1016/0003-9861(81)90073-4
Li, M., Ding, W., Baruah, B., Crans, D. C., and Wang, R. (2008). Inhibition of protein tyrosine phosphatase 1B and alkaline phosphatase by bis(maltolato)oxovanadium (IV). J. Inorg. Biochem. 102 (10), 1846–1853. doi:10.1016/j.jinorgbio.2008.06.007
Lu, L., Gao, X., Zhu, M., Wang, S., Wu, Q., Xing, S., et al. (2012). Exploration of biguanido–oxovanadium complexes as potent and selective inhibitors of protein tyrosine phosphatases. BioMetals 25 (3), 599–610. doi:10.1007/s10534-012-9548-4
McLauchlan, C. C., Peters, B. J., Willsky, G. R., and Crans, D. C. (2015). Vanadium–phosphatase complexes: phosphatase inhibitors favor the trigonal bipyramidal transition state geometries. Coord. Chem. Rev. 301-302, 163–199. doi:10.1016/j.ccr.2014.12.012
National Institute of Diabetes and Digestive and Kidney Diseases (2024). Ad bax group software - biophysical nuclear magnetic resonance spectroscopy section - NIDDK. Available at: https://spin.niddk.nih.gov/bax/software/TALOS/(Accessed November 13, 2017).
Nayak, J., Gastonguay, A. J., Talipov, M. R., Padmanabhan, V., Span, E. A., Kalous, K. S., et al. (2014). Protein expression, characterization and activity comparisons of wild type and mutant DUSP5 proteins. BMC Biochem. 15, 27. doi:10.1186/s12858-014-0027-0
Neumann, T. S., Span, E. A., Kalous, K. S., Bongard, R., Gastonguay, A., Lepley, M. A., et al. (2015). Identification of inhibitors that target dual-specificity phosphatase 5 provide new insights into the binding requirements for the two phosphate pockets. BMC Biochem. 16, 19. doi:10.1186/s12858-015-0048-3
Pannifer, A. D. B., Flint, A. J., Tonks, N. K., and Barford, D. (1998). Visualization of the cysteinyl-phosphate intermediate of a protein-tyrosine phosphatase by X-ray crystallography. J. Biol. Chem. 273 (17), 10454–10462. doi:10.1074/jbc.273.17.10454
Pedersen, A. K., Guo, X. L., Moller, K. B., Peters, G. H., Andersen, H. S., Kastrup, J. S., et al. (2004). Residue 182 influences the second step of protein-tyrosine phosphatase-mediated catalysis. Biochem. J. 378(2), 421–433. doi:10.1042/bj20030565
Pettersen, E. F., Goddard, T. D., Huang, C. C., Couch, G. S., Greenblatt, D. M., Meng, E. C., et al. (2004). UCSF Chimera--A visualization system for exploratory research and analysis. J. Comput. Chem. 25 (13), 1605–1612. doi:10.1002/jcc.20084
Pramanik, K., Chun, C. Z., Garnaas, M. K., Samant, G. V., Li, K., Horswill, M. A., et al. (2009). Dusp-5 and Snrk-1 coordinately function during vascular development and disease. Blood 113 (5), 1184–1191. doi:10.1182/blood-2008-06-162180
Raman, M., Chen, W., and Cobb, M. H. (2007). Differential regulation and properties of MAPKs. Oncogene 26 (22), 3100–3112. doi:10.1038/sj.onc.1210392
Reiter, N. J., White, D. J., and Rusnak, F. (2001). Inhibition of bacteriophage λ protein phosphatase by organic and oxoanion inhibitors. Biochemistry 41 (3), 1051–1059. doi:10.1021/bi011577b
Seleman, M., Chapy, H., Cisternino, S., Courtin, C., Smirnova, M., Schlatter, J., et al. (2014). Impact of P-glycoprotein at the blood-brain barrier on the uptake of heroin and its main metabolites: behavioral effects and consequences on the transcriptional responses and reinforcing properties. Psychopharmacology 231 (16), 3139–3149. doi:10.1007/s00213-014-3490-9
Shen, Y., Delaglio, F., Cornilescu, G., and Bax, A. (2009). TALOS+: a hybrid method for predicting protein backbone torsion angles from NMR chemical shifts. J. Biomol. NMR 44 (4), 213–223. doi:10.1007/s10858-009-9333-z
Talipov, M. R., Nayak, J., Lepley, M., Bongard, R. D., Sem, D. S., Ramani, R., et al. (2016). Critical role of the secondary binding pocket in modulating the enzymatic activity of DUSP5 toward phosphorylated ERKs. Biochemistry 55 (44), 6187–6195. doi:10.1021/acs.biochem.6b00498
Ueda, K., Arakawa, H., and Nakamura, Y. (2003). Dual-specificity phosphatase 5 (DUSP5) as a direct transcriptional target of tumor suppressor p53. Oncogene 22 (36), 5586–5591. doi:10.1038/sj.onc.1206845
Zhang, M., Zhou, M., Van, R. L., and Stauffacher, C. V. (1997). Crystal structure of bovine low molecular weight phosphotyrosyl phosphatase complexed with the transition state analog vanadate. Biochemistry 36 (1), 15–23. doi:10.1021/bi961804n
Zhang, Y. L., and Zhang, Z. Y. (1998). Low-affinity binding determined by titration calorimetry using a high-affinity coupling ligand: a thermodynamic study of ligand binding to protein tyrosine phosphatase 1B. Anal. Biochem. 261 (2), 139–148. doi:10.1006/abio.1998.2738
Zhou, M. M., Davis, J. P., and Van Etten, R. L. (1993). Identification and pKa determination of the histidine residues of human low-molecular-weight phosphotyrosyl protein phosphatases: a convenient approach using MLEV-17 spectral editing scheme. Biochemistry 32 (33), 8479–8486. doi:10.1021/bi00084a012
Keywords: DUSP5, pERK, phosphatase, regulation, enzyme kinetics, phosphorylated peptide substrates, dual active site, vanadate
Citation: Imhoff A, Sweeney NL, Bongard RD, Syrlybaeva R, Gupta A, Del Carpio E, Talipov MR, Garcia-Keller C, Crans DC, Ramchandran R and Sem DS (2024) Structural and kinetic characterization of DUSP5 with a Di-phosphorylated tripeptide substrate from the ERK activation loop. Front. Chem. Biol 3:1385560. doi: 10.3389/fchbi.2024.1385560
Received: 13 February 2024; Accepted: 15 July 2024;
Published: 05 August 2024.
Edited by:
Adam Offenbacher, East Carolina University, United StatesReviewed by:
Sourav Roy, St. Jude Children’s Research Hospital, United StatesMing Dong, University of North Carolina Wilmington, United States
Copyright © 2024 Imhoff, Sweeney, Bongard, Syrlybaeva, Gupta, Del Carpio, Talipov, Garcia-Keller, Crans, Ramchandran and Sem. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Daniel S. Sem, Daniel.Sem@cuw.edu; Ramani Ramchandran, rramchan@mcw.edu; Debbie C. Crans, debbie.crans@colostate.edu; Marat R. Talipov, talipovm@nmsu.edu; Costanza Garcia-Keller, cgkeller@mcw
†Present address: Raulia Syrlybaeva, Washington University School of Medicine, Division of Infectious Diseases, St. Louis, MO, United StatesEdgar Del Carpio, Unidad de Química Medicinal, Facultad de Farmacia, Universidad Central de Venezuela, CCS 1050, VenezuelaDaniel S. Sem, Walsh College, Troy, Michigan, United States