Atomic structure of the MAP kinase ERK2 at 2.3 A resolution

Apprehensions and emerging solutions in ML-based protein structure prediction

The three-dimensional structure of proteins determines their function in vital biological processes. Thus, when the structure is known, the molecular mechanism of protein function can be understood in more detail and obtained information utilized in biotechnological, diagnostics, and therapeutic applications. Over the past five years, machine learning (ML)-based modeling has pushed protein structure prediction to the next level with AlphaFold in the front line, predicting the structure for hundreds of millions of proteins. Further advances recently report promising ML-based approaches for solving remaining challenges by incorporating functionally important metals, co-factors, post-translational modifications, structural dynamics, and interdomain and multimer interactions in the structure prediction process.

Is allostery a fuzzy concept?

Allostery is an important property of biological macromolecules which regulates diverse biological functions such as catalysis, signal transduction, transport, and molecular recognition. However, the concept was expressed using two different definitions by J. Monod and, over time, more have been added by different authors, making it fuzzy. Here, we reviewed the different meanings of allostery in the current literature and found that it has been used to indicate that the function of a protein is regulated by heterotropic ligands, and/or that the binding of ligands and substrates presents homotropic positive or negative cooperativity, whatever the hypothesized or demonstrated reaction mechanism might be. Thus, proteins defined to be allosteric include not only those that obey the two-state concerted model, but also those that obey different reaction mechanisms such as ligand-induced fit, possibly coupled to sequential structure changes, and ligand-linked dissociation-association. Since each reaction mechanism requires its own mathematical description and is defined by it, there are many possible 'allosteries'. This lack of clarity is made even fuzzier by the fact that the reaction mechanism is often assigned imprecisely and/or implicitly in the absence of the necessary experimental evidence. In this review, we examine a list of proteins that have been defined to be allosteric and attempt to assign a reaction mechanism to as many as possible.

Conformation selection by ATP-competitive inhibitors and allosteric communication in ERK2

Activation of the extracellular signal-regulated kinase-2 (ERK2) by phosphorylation has been shown to involve changes in protein dynamics, as determined by hydrogen-deuterium exchange mass spectrometry (HDX-MS) and NMR relaxation dispersion measurements. These can be described by a global exchange between two conformational states of the active kinase, named 'L' and 'R,' where R is associated with a catalytically productive ATP-binding mode. An ATP-competitive ERK1/2 inhibitor, Vertex-11e, has properties of conformation selection for the R-state, revealing movements of the activation loop that are allosterically coupled to the kinase active site. However, the features of inhibitors important for R-state selection are unknown. Here, we survey a panel of ATP-competitive ERK inhibitors using HDX-MS and NMR and identify 14 new molecules with properties of R-state selection. They reveal effects propagated to distal regions in the P+1 and helix αF segments surrounding the activation loop, as well as helix αL16. Crystal structures of inhibitor complexes with ERK2 reveal systematic shifts in the Gly loop and helix αC, mediated by a Tyr-Tyr ring stacking interaction and the conserved Lys-Glu salt bridge. The findings suggest a model for the R-state involving small movements in the N-lobe that promote compactness within the kinase active site and alter mobility surrounding the activation loop. Such properties of conformation selection might be exploited to modulate the protein docking interface used by ERK substrates and effectors.

Synergistic anticancer effect by targeting CDK2 and EGFR-ERK signaling

The EGFR-RAS-ERK pathway is one of the most important signaling cascades in cell survival, growth, and proliferation. Aberrant activation of this pathway is a common mechanism in various cancers. Here, we report that CDK2 is a novel regulator of the ERK pathway via USP37 deubiquitinase (DUB). Mechanistically, CDK2 phosphorylates USP37, which is required for USP37 DUB activity. Further, USP37 deubiquitinates and stabilizes ERK1/2, thereby enhancing cancer cell proliferation. Thus, CDK2 is able to promote cell proliferation by activating USP37 and, in turn, stabilizing ERK1/2. Importantly, combined CDK1/2 and EGFR inhibitors have a synergetic anticancer effect through the downregulation of ERK1/2 stability and activity. Indeed, our patient-derived xenograft (PDX) results suggest that targeting both ERK1/2 stability and activity kills cancer cells more efficiently even at lower doses of these two inhibitors, which may reduce their associated side effects and indicate a potential new combination strategy for cancer therapy.

Activation Loop Plasticity and Active Site Coupling in the MAP Kinase, ERK2

Previous studies of the protein kinase, ERK2, using NMR and hydrogen-exchange measurements have shown changes in dynamics accompanying its activation by phosphorylation. However, knowledge about the conformational motions involved is incomplete. Here, we examined ERK2 using long conventional molecular dynamics (MD) simulations starting from crystal structures of phosphorylated (2P) and unphosphorylated (0P) forms. Individual trajectories were run for (5 to 25) μs, totaling 727 μs. The results show unexpected flexibility of the A-loop, with multiple long-lived (>5 μs) conformational states in both 2P- and 0P-ERK2. Differential contact network and principal component analyses reveal coupling between the A-loop fold and active site dynamics, with evidence for conformational selection in the kinase core of 2P-ERK2 but not 0P-ERK2. Simulations of 2P-ERK2 show A-loop states corresponding to restrained dynamics within the N-lobe, including regions around catalytic residues. One A-loop conformer forms lasting interactions with the L16 segment, leading to reduced RMSF and greater compaction in the active site. By contrast, simulations of 0P-ERK2 reveal excursions of A-loop residues away from the C-lobe, leading to greater active site mobility. Thus, the A-loop in ERK2 switches between distinct conformations that reflect coupling with the active site, possibly via the L16 segment. Crystal packing interactions suggest that lattice contacts with the A-loop may restrain its structural variation in X-ray structures of ERK2. The novel conformational states identified by MD expand our understanding of ERK2 regulation, by linking the activated state of the kinase to reduced dynamics and greater compaction surrounding the catalytic site.

Uncoupling p38α nuclear and cytoplasmic functions and identification of two p38α phosphorylation sites on β-catenin: implications for the Wnt signaling pathway in CRC models

BACKGROUND

Activation of the Wnt pathway has been linked to colorectal cancer (CRC). Previous reports suggest that Wnt3a can activate p38. Besides, p38α feeds into the canonical Wnt/β-catenin pathway by inhibiting GSK3β through phosphorylation. Recently, we identified p38α as a new druggable member of β-catenin chromatin-associated kinase complexes in CRC.

METHODS

The functional relationship between p38α and β-catenin was characterized in CRC cells, patient-derived CRC stem cells, patient-derived tumor intestinal organoids, and in vivo models (C57BL/6-APCMin/+ mice). The role of p38α in β-catenin transcriptional activity was assessed by pharmacological inhibition with ralimetinib.

RESULTS

We used the GSK3β inhibitor TWS-119, which promotes the activation of Wnt signaling, to uncouple p38α nuclear/cytoplasmatic functions in the Wnt pathway. Upon GSK3β inhibition, nuclear p38α phosphorylates β-catenin at residues S111 and T112, allowing its binding to promoter regions of Wnt target genes and the activation of a transcriptional program implicated in cancer progression. If p38α is pharmacologically inhibited in addition to GSK3β, β-catenin is prevented from promoting target gene transcription, which is expected to impair carcinogenesis.

CONCLUSIONS

p38α seems to play a dual role as a member of the β-catenin destruction complex and as a β-catenin chromatin-associated kinase in CRC. This finding may help elucidate mechanisms contributing to human colon tumor pathogenesis and devise new strategies for personalized CRC treatment.

Conformation Selection by ATP-competitive Inhibitors and Allosteric Communication in ERK2

Activation of the extracellular signal regulated kinase-2 (ERK2) by phosphorylation has been shown to involve changes in protein dynamics, as determined by hydrogen-deuterium exchange mass spectrometry (HDX-MS) and NMR relaxation dispersion measurements. These can be described by a global exchange between two conformational states of the active kinase, named "L" and "R", where R is associated with a catalytically productive ATP-binding mode. An ATP-competitive ERK1/2 inhibitor, Vertex-11e, has properties of conformation selection for the R-state, revealing movements of the activation loop that are allosterically coupled to the kinase active site. However, the features of inhibitors important for R-state selection are unknown. Here we survey a panel of ATP-competitive ERK inhibitors using HDX-MS and NMR and identify 14 new molecules with properties of R-state selection. They reveal effects propagated to distal regions in the P+1 and helix αF segments surrounding the activation loop, as well as helix αL16. Crystal structures of inhibitor complexes with ERK2 reveal systematic shifts in the Gly loop and helix αC, mediated by a Tyr-Tyr ring stacking interaction and the conserved Lys-Glu salt bridge. The findings suggest a model for the R-state involving small movements in the N-lobe that promote compactness within the kinase active site and alter mobility surrounding the activation loop. Such properties of conformation selection might be exploited to modulate the protein docking interface used by ERK substrates and effectors.

Signal transduction at GPCRs: Allosteric activation of the ERK MAPK by β-arrestin

β-arrestins are multivalent adaptor proteins that bind active phosphorylated G protein-coupled receptors (GPCRs) to inhibit G protein signaling, mediate receptor internalization, and initiate alternative signaling events. β-arrestins link agonist-stimulated GPCRs to downstream signaling partners, such as the c-Raf-MEK1-ERK1/2 cascade leading to ERK1/2 activation. β-arrestins have been thought to transduce signals solely via passive scaffolding by facilitating the assembly of multiprotein signaling complexes. Recently, however, β-arrestin 1 and 2 were shown to activate two downstream signaling effectors, c-Src and c-Raf, allosterically. Over the last two decades, ERK1/2 have been the most intensely studied signaling proteins scaffolded by β-arrestins. Here, we demonstrate that β-arrestins play an active role in allosterically modulating ERK kinase activity in vitro and within intact cells. Specifically, we show that β-arrestins and their GPCR-mediated active states allosterically enhance ERK2 autophosphorylation and phosphorylation of a downstream ERK2 substrate, and we elucidate the mechanism by which β-arrestins do so. Furthermore, we find that allosteric stimulation of dually phosphorylated ERK2 by active-state β-arrestin 2 is more robust than by active-state β-arrestin 1, highlighting differential capacities of β-arrestin isoforms to regulate effector signaling pathways downstream of GPCRs. In summary, our study provides strong evidence for a new paradigm in which β-arrestins function as active "catalytic" scaffolds to allosterically unlock the enzymatic activity of signaling components downstream of GPCR activation.

Navigating the ERK1/2 MAPK Cascade

The RAS-ERK pathway is a fundamental signaling cascade crucial for many biological processes including proliferation, cell cycle control, growth, and survival; common across all cell types. Notably, ERK1/2 are implicated in specific processes in a context-dependent manner as in stem cells and pancreatic β-cells. Alterations in the different components of this cascade result in dysregulation of the effector kinases ERK1/2 which communicate with hundreds of substrates. Aberrant activation of the pathway contributes to a range of disorders, including cancer. This review provides an overview of the structure, activation, regulation, and mutational frequency of the different tiers of the cascade; with a particular focus on ERK1/2. We highlight the importance of scaffold proteins that contribute to kinase localization and coordinate interaction dynamics of the kinases with substrates, activators, and inhibitors. Additionally, we explore innovative therapeutic approaches emphasizing promising avenues in this field.

Probing conformational dynamics to understand kinase inhibition

Why do some inhibitors select the on-state in ERK2, a kinase that is involved in many signaling pathways in cells, whereas others bind to more than one conformation?

Architecture of the MKK6-p38α complex defines the basis of MAPK specificity and activation

The mitogen-activated protein kinase (MAPK) p38α is a central component of signaling in inflammation and the immune response and is, therefore, an important drug target. Little is known about the molecular mechanism of its activation by double phosphorylation from MAPK kinases (MAP2Ks), because of the challenge of trapping a transient and dynamic heterokinase complex. We applied a multidisciplinary approach to generate a structural model of p38α in complex with its MAP2K, MKK6, and to understand the activation mechanism. Integrating cryo-electron microscopy with molecular dynamics simulations, hydrogen-deuterium exchange mass spectrometry, and experiments in cells, we demonstrate a dynamic, multistep phosphorylation mechanism, identify catalytically relevant interactions, and show that MAP2K-disordered amino termini determine pathway specificity. Our work captures a fundamental step of cell signaling: a kinase phosphorylating its downstream target kinase.

Mitogen-activated protein kinase (MAPK) cascades-A yeast perspective

Discovery of the class of protein kinase now dubbed a mitogen (or messenger)-activated protein kinase (MAPK) is an illustrative example of how disparate lines of investigation can converge and reveal an enzyme family universally conserved among eukaryotes, from single-celled microbes to humans. Moreover, elucidation of the circuitry controlling MAPK function defined a now overarching principle in enzyme regulation-the concept of an activation cascade mediated by sequential phosphorylation events. Particularly ground-breaking for this field of exploration were the contributions of genetic approaches conducted using several model organisms, but especially the budding yeast Saccharomyces cerevisiae. Notably, examination of how haploid yeast cells respond to their secreted peptide mating pheromones was crucial in pinpointing genes encoding MAPKs and their upstream activators. Fully contemporaneous biochemical analysis of the activities elicited upon stimulation of mammalian cells by insulin and other growth- and differentiation-inducing factors lead eventually to the demonstration that components homologous to those in yeast were involved. Continued studies of these pathways in yeast were integral to other foundational discoveries in MAPK signaling, including the roles of tethering, scaffolding and docking interactions.

Activation loop plasticity and active site coupling in the MAP kinase, ERK2

Changes in the dynamics of the protein kinase, ERK2, have been shown to accompany its activation by dual phosphorylation. However, our knowledge about the conformational changes represented by these motions is incomplete. Previous NMR relaxation dispersion studies showed that active, dual-phosphorylated ERK2 undergoes global exchange between at least two energetically similar conformations. These findings, combined with measurements by hydrogen exchange mass spectrometry (HX-MS), suggested that the global conformational exchange involves motions of the activation loop (A-loop) that are coupled to regions surrounding the kinase active site. In order to better understand the contribution of dynamics to the activation of ERK2, we applied long conventional molecular dynamics (MD) simulations starting from crystal structures of active, phosphorylated (2P), and inactive, unphosphorylated (0P) ERK2. Individual trajectories were run for (5 to 25) µ s and totaled 727 µ s. The results showed that the A-loop is unexpectedly flexible in both 2P- and 0P-ERK2, and able to adopt multiple long-lived (>5 µ s) conformational states. Simulations starting from the X-ray structure of 2P-ERK2 (2ERK) revealed A-loop states corresponding to restrained dynamics within the N-lobe, including regions surrounding catalytic residues. One A-loop conformer forms lasting interactions with the C-terminal L16 segment and shows reduced RMSF and greater compaction in the active site. By contrast, simulations starting from the most common X-ray conformation of 0P-ERK2 (5UMO) reveal frequent excursions of A-loop residues away from a C-lobe docking site pocket and towards a new state that shows greater dynamics in the N-lobe and disorganization around the active site. Thus, the A-loop in ERK2 appears to switch between distinct conformational states that reflect allosteric coupling with the active site, likely occurring via the L16 segment. Analyses of crystal packing interactions across many structural datasets suggest that the A-loop observed in X-ray structures of ERK2 may be driven by lattice contacts and less representative of the solution structure. The novel conformational states identified by MD expand our understanding of ERK2 regulation, by linking the activated state of the kinase to reduced dynamics and greater compaction surrounding the catalytic site.

17β-Estradiol Concentration and Direct β2-Adrenoceptor Inhibition Determine Estrogen-Mediated Reversal of Adrenergic Immunosuppression

Background:

Sympathetic innervation of lymphoid organs, and the presence of 17β-estradiol (estrogen or E2) and adrenergic receptors (ARs) on lymphocytes, suggests that sympathetic stimulation and hormonal activation may influence immune functions.

Purpose:

Modeling and simulating these pathways may help to understand the dynamics of neuroendocrine-immune modulation at the cellular and molecular levels.

Methods:

Dose- and receptor-dependent effects of E2 and AR subtype-specific agonists were established in vitro on lymphocytes from young male Sprague-Dawley rats and were modeled in silico using the MATLAB Simbiology toolbox. Kinetic principles were assigned to define receptor–ligand dynamics, and concentration/time plots were obtained using Ode15s solvers at different time intervals for key regulatory molecules. Comparisons were drawn between in silico and in vitro data for validating the constructed model with sensitivity analysis of key regulatory molecules to assess their individual impacts on the dynamics of the system. Finally, docking studies were conducted with key ligands E2 and norepinephrine (NE) to understand the mechanistic principles underlying their interactions.

Results:

Adrenergic activation triggered proapoptotic signals, while E2 enhanced survival signals, showing opposing effects as observed in vitro. Treatment of lymphocytes with E2 shows a 10-fold increase in survival signals in a dose-dependent manner. Cyclic adenosine monophosphate (cAMP) activation is crucial for the activation of survival signals through extracellular signal-regulated kinase (p-ERK) and cAMP responsive element binding (p-CREB) protein. Docking studies showed the direct inhibition of ERK by NE and β2-AR by E2 explaining how estrogen signaling overrides NE-mediated immunosuppression in vitro.

Conclusion:

The cross-talk between E2 and adrenergic signaling pathways determines lymphocyte functions in a receptor subtype and coactivation-dependent manner in health and disease.

A non-catalytic herpesviral protein reconfigures ERK-RSK signaling by targeting kinase docking systems in the host

The Kaposi's sarcoma associated herpesvirus protein ORF45 binds the extracellular signal-regulated kinase (ERK) and the p90 Ribosomal S6 kinase (RSK). ORF45 was shown to be a kinase activator in cells but a kinase inhibitor in vitro, and its effects on the ERK-RSK complex are unknown. Here, we demonstrate that ORF45 binds ERK and RSK using optimized linear binding motifs. The crystal structure of the ORF45-ERK2 complex shows how kinase docking motifs recognize the activated form of ERK. The crystal structure of the ORF45-RSK2 complex reveals an AGC kinase docking system, for which we provide evidence that it is functional in the host. We find that ORF45 manipulates ERK-RSK signaling by favoring the formation of a complex, in which activated kinases are better protected from phosphatases and docking motif-independent RSK substrate phosphorylation is selectively up-regulated. As such, our data suggest that ORF45 interferes with the natural design of kinase docking systems in the host.

The two non-visual arrestins engage ERK2 differently

Arrestin binding to active phosphorylated G protein-coupled receptors terminates G protein coupling and initiates another wave of signaling. Among the effectors that bind directly to receptor-associated arrestins are extracellular signal-regulated kinases 1/2 (ERK1/2), which promote cellular proliferation and survival. Arrestins may also engage ERK1/2 in isolation in a pre- or post-signaling complex that is likely in equilibrium with the full signal initiation complex. Molecular details of these binary complexes remain unknown. Here, we investigate the molecular mechanisms whereby arrestin-2 and arrestin-3 (a.k.a. β-arrestin1 and β-arrestin2, respectively) engage ERK1/2 in pairwise interactions. We find that purified arrestin-3 binds ERK2 more avidly than arrestin-2. A combination of biophysical techniques and peptide array analysis demonstrates that the molecular basis in this difference of binding strength is that the two non-visual arrestins bind ERK2 via different parts of the molecule. We propose a structural model of the ERK2-arrestin-3 complex in solution using size-exclusion chromatography coupled to small angle X-ray scattering (SEC-SAXS). This binary complex exhibits conformational heterogeneity. We speculate that this drives the equilibrium either toward the full signaling complex with receptor-bound arrestin at the membrane or toward full dissociation in the cytoplasm. As ERK1/2 regulates cell migration, proliferation, and survival, understanding complexes that relate to its activation could be exploited to control cell fate.

Kinobead Profiling Reveals Reprogramming of BCR Signaling in Response to Therapy within Primary CLL Cells

PURPOSE

B-cell receptor (BCR) signaling is critical for the pathogenesis of chronic lymphocytic leukemia (CLL), promoting both malignant cell survival and disease progression. Although vital, understanding of the wider signaling network associated with malignant BCR stimulation is poor. This is relevant with respect to potential changes in response to therapy, particularly involving kinase inhibitors. In the current study, we describe a novel high-resolution approach to investigate BCR signaling in primary CLL cells and track the influence of therapy on signaling response.

EXPERIMENTAL DESIGN

A kinobead/mass spectrometry-based protocol was used to study BCR signaling in primary CLL cells. Longitudinal analysis of samples donated by clinical trial patients was used to investigate the impact of chemoimmunotherapy and ibrutinib on signaling following surface IgM engagement. Complementary Nanostring and immunoblotting analysis was used to verify our findings.

RESULTS

Our protocol isolated a unique, patient-specific signature of over 30 kinases from BCR-stimulated CLL cells. This signature was associated with 13 distinct Kyoto Encyclopedia of Genes and Genomes pathways and showed significant change in cells from treatment-naïve patients compared with those from patients who had previously undergone therapy. This change was validated by longitudinal analysis of clinical trials samples where BCR-induced kinome responses in CLL cells altered between baseline and disease progression in patients failing chemoimmunotherapy and between baseline and treatment in patients taking ibrutinib.

CONCLUSIONS

These data comprise the first comprehensive proteomic investigation of the BCR signaling response within CLL cells and reveal unique evidence that these cells undergo adaptive reprogramming of this signaling in response to therapy.

c-Jun N-terminal kinase (JNK) in Procambarus clarkii: Molecular characterization and involvement in oxidative stress-induced apoptosis during molting cycle

The present study was conducted to characterize the full-length cDNA of c-Jun N-terminal kinase (JNK) in Procambarus clarkii (Pcjnk) and evaluate its potential function under different molt cycle. The full-length cDNA of Pcjnk covered 2937 bp with an open reading frame of 1320 bp, encoding 439 amino acids. A typical conserved TPY motif (118Thr-Pro-120Tyr) was found in Pcjnk. Quantitative real-time PCR (qRT-PCR) analysis revealed a constitutive expression of Pcjnk in the tested tissue, with the highest expression occurring in the hepatopancreas. Additionally, the present study initially revealed that relative mRNA expression of Pcjnk and apoptosis level were significantly higher in the premolt stage (D1/D2 and D3/D4 stage) as compared to other molt stages. In contrast to the levels of superoxide dismutase (SOD) and malondialdehyde (MDA), catalase (CAT) and glutathione peroxidase (GPX) level decreased significantly from the intermolt stage (C stage) to the premolt stage (D1/D2 and D3/D4 stage), then increased from the premolt stage to the postmolt stage (A and B stage). The results obtained in the present study indicated that molt could cause apoptosis induced by oxidative stress through the activation of JNK in Procambarus clarkii.

Predicting aptamer sequences that interact with target proteins using an aptamer-protein interaction classifier and a Monte Carlo tree search approach

Oligonucleotide-based aptamers, which have a three-dimensional structure with a single-stranded fragment, feature various characteristics with respect to size, toxicity, and permeability. Accordingly, aptamers are advantageous in terms of diagnosis and treatment and are materials that can be produced through relatively simple experiments. Systematic evolution of ligands by exponential enrichment (SELEX) is one of the most widely used experimental methods for generating aptamers; however, it is highly expensive and time-consuming. To reduce the related costs, recent studies have used in silico approaches, such as aptamer-protein interaction (API) classifiers that use sequence patterns to determine the binding affinity between RNA aptamers and proteins. Some of these methods generate candidate RNA aptamer sequences that bind to a target protein, but they are limited to producing candidates of a specific size. In this study, we present a machine learning approach for selecting candidate sequences of various sizes that have a high binding affinity for a specific sequence of a target protein. We applied the Monte Carlo tree search (MCTS) algorithm for generating the candidate sequences using a score function based on an API classifier. The tree structure that we designed with MCTS enables nucleotide sequence sampling, and the obtained sequences are potential aptamer candidates. We performed a quality assessment using the scores of docking simulations. Our validation datasets revealed that our model showed similar or better docking scores in ZDOCK docking simulations than the known aptamers. We expect that our method, which is size-independent and easy to use, can provide insights into searching for an appropriate aptamer sequence for a target protein during the simulation step of SELEX.

Calculating Transfer Entropy from Variance-Covariance Matrices Provides Insight into Allosteric Communication in ERK2

We develop an approach by which reliable estimates of the transfer entropy can be obtained from the variance-covariance matrix of atomic fluctuations, which converges quickly and retains sensitivity to the full chemical profile of the biomolecular system. We validate our method on ERK2, a well-studied kinase involved in the MAPK signaling cascade for which considerable computational, experimental, and mutation data are available. We present the results of transfer entropy analysis on data obtained from molecular dynamics simulations of wild-type active and inactive ERK2, along with mutants Q103A, I84A, L73P, and G83A. We show that our method is systematically consistent within the context of other approaches for calculating transfer entropy, and we provide a method for interpreting networks of interconnected residues in the protein from a perspective of allosteric coupling. We introduce new insights about possible allosteric activity of the extreme N-terminal region of the kinase, and we describe evidence that suggests that activation may occur by different paths or routes in different mutants. Our results highlight systematic advantages and disadvantages of each method for calculating transfer entropy and show the important role of transfer entropy analysis for understanding allosteric behavior in biomolecular systems.

6-(methylsulfinyl)hexyl isothiocyanate (6-MITC) from Wasabia japonica alleviates inflammatory bowel disease (IBD) by potential inhibition of glycogen synthase kinase 3 beta (GSK-3β)

Inflammatory bowel disease (IBD) describes a set of disorders involving alterations to gastrointestinal physiology and mucosal immunity. Unravelling its complex pathophysiology is important since many IBD patients are refractory to or suffer adverse side effects from current treatments. Isothiocyanates (ITCs), such as 6-(methylsulfinyl)hexyl ITC (6-MITC) in Wasabia japonica, have potential anti-inflammatory activity. We aimed to elucidate the pathways through which 6-MITC alleviates inflammation by examining its role in the nuclear factor-kappa B (NF-κB) pathway through inhibition of glycogen synthase kinase 3 beta (GSK-3β) using a chemically induced murine model of IBD, cell-based and in silico techniques. The effects of 6-MITC and two NF-κB inhibitors, sulfasalazine (SS), pyrrolidine dithiolcarbamate (PDTC) were investigated on a dextran sulfate sodium (DSS)-induced murine mouse model of acute and chronic colitis using macroscopic measurements and pro-inflammatory markers. The effect of 6-MITC on NF-κB induction was assessed using a murine macrophage cell line. Complexes of GSK-3β-6-MITC and GSK-3β-ATP were generated in silico to elucidate the mechanism of 6-MITC's direct inhibition of GSK-3β. Changes in pro-inflammatory markers, inducible nitric oxide synthase (iNOS) (increased) and interleukin-6 (IL-6) (decreased) demonstrated that iNOS regulation occurred at the translational level. Intraperitoneal (ip) injection of 6-MITC to the colitis-induced mice ameliorated weight loss whereas oral administration had negligible effect. Fecal blood and colon weight/length ratio parameters improved on treatment with 6-MITC and the other NF-κB inhibitors. Levels of NF-κB decreased upon addition of 6-MITC in vitro while structural studies showed 6-MITC acts competitively to inhibit GSK-3β at the ATP binding site. In this study we demonstrated that 6-MITC inhibits NF-κB signaling via GSK-3β inhibition ameliorating fecal blood, colonic alterations and DSS-induced weight loss indirectly indicating reduced intestinal stress. Taken together these results suggest a role for 6-MITC in the treatment of IBD acting to alleviate inflammation through the GSK-3β/NF-κB pathway. Furthermore, the GSK-3β-6-MITC model can be utilized as a basis for development of novel therapeutics targeting GSK-3β for use in other disorders including cancer.

Crystal Structure and Inhibitor Identifications Reveal Targeting Opportunity for the Atypical MAPK Kinase ERK3

Extracellular signal-regulated kinase 3 (ERK3), known also as mitogen-activated protein kinase 6 (MAPK6), is an atypical member of MAPK kinase family, which has been poorly studied. Little is known regarding its function in biological processes, yet this atypical kinase has been suggested to play important roles in the migration and invasiveness of certain cancers. The lack of tools, such as a selective inhibitor, hampers the study of ERK3 biology. Here, we report the crystal structure of the kinase domain of this atypical MAPK kinase, providing molecular insights into its distinct ATP binding pocket compared to the classical MAPK ERK2, explaining differences in their inhibitor binding properties. Medium-scale small molecule screening identified a number of inhibitors, several of which unexpectedly exhibited remarkably high inhibitory potencies. The crystal structure of CLK1 in complex with CAF052, one of the most potent inhibitors identified for ERK3, revealed typical type-I binding mode of the inhibitor, which by structural comparison could likely be maintained in ERK3. Together with the presented structural insights, these diverse chemical scaffolds displaying both reversible and irreversible modes of action, will serve as a starting point for the development of selective inhibitors for ERK3, which will be beneficial for elucidating the important functions of this understudied kinase.

Enhanced MAPK1 Function Causes a Neurodevelopmental Disorder within the RASopathy Clinical Spectrum

Signal transduction through the RAF-MEK-ERK pathway, the first described mitogen-associated protein kinase (MAPK) cascade, mediates multiple cellular processes and participates in early and late developmental programs. Aberrant signaling through this cascade contributes to oncogenesis and underlies the RASopathies, a family of cancer-prone disorders. Here, we report that de novo missense variants in MAPK1, encoding the mitogen-activated protein kinase 1 (i.e., extracellular signal-regulated protein kinase 2, ERK2), cause a neurodevelopmental disease within the RASopathy phenotypic spectrum, reminiscent of Noonan syndrome in some subjects. Pathogenic variants promote increased phosphorylation of the kinase, which enhances translocation to the nucleus and boosts MAPK signaling in vitro and in vivo. Two variant classes are identified, one of which directly disrupts binding to MKP3, a dual-specificity protein phosphatase negatively regulating ERK function. Importantly, signal dysregulation driven by pathogenic MAPK1 variants is stimulus reliant and retains dependence on MEK activity. Our data support a model in which the identified pathogenic variants operate with counteracting effects on MAPK1 function by differentially impacting the ability of the kinase to interact with regulators and substrates, which likely explains the minor role of these variants as driver events contributing to oncogenesis. After nearly 20 years from the discovery of the first gene implicated in Noonan syndrome, PTPN11, the last tier of the MAPK cascade joins the group of genes mutated in RASopathies.

Growth Inhibitory Signaling of the Raf/MEK/ERK Pathway

In response to extracellular stimuli, the Raf/MEK/extracellular signal-regulated kinase (ERK) pathway regulates diverse cellular processes. While mainly known as a mitogenic signaling pathway, the Raf/MEK/ERK pathway can mediate not only cell proliferation and survival but also cell cycle arrest and death in different cell types. Growing evidence suggests that the cell fate toward these paradoxical physiological outputs may be determined not only at downstream effector levels but also at the pathway level, which involves the magnitude of pathway activity, spatial-temporal regulation, and non-canonical functions of the molecular switches in this pathway. This review discusses recent updates on the molecular mechanisms underlying the pathway-mediated growth inhibitory signaling, with a major focus on the regulation mediated at the pathway level.

Activation Loop Dynamics Are Coupled to Core Motions in Extracellular Signal-Regulated Kinase-2

The activation loop segment in protein kinases is a common site for regulatory phosphorylation. In extracellular signal-regulated kinase 2 (ERK2), dual phosphorylation and conformational rearrangement of the activation loop accompany enzyme activation. X-ray structures show the active conformation to be stabilized by multiple ion pair interactions between phosphorylated threonine and tyrosine residues in the loop and six arginine residues in the kinase core. Despite the extensive salt bridge network, nuclear magnetic resonance Carr-Purcell-Meiboom-Gill relaxation dispersion experiments show that the phosphorylated activation loop is conformationally mobile on a microsecond to millisecond time scale. The dynamics of the loop match those of previously reported global exchange within the kinase core region and surrounding the catalytic site that have been found to facilitate productive nucleotide binding. Mutations in the core region that alter these global motions also alter the dynamics of the activation loop. Conversely, mutations in the activation loop perturb the global exchange within the kinase core. Together, these findings provide evidence for coupling between motions in the activation loop and those surrounding the catalytic site in the active state of the kinase. Thus, the activation loop segment in dual-phosphorylated ERK2 is not held statically in the active X-ray conformation but instead undergoes exchange between conformers separated by a small energetic barrier, serving as part of a dynamic allosteric network controlling nucleotide binding and catalytic function.

ERK signalling: a master regulator of cell behaviour, life and fate

The proteins extracellular signal-regulated kinase 1 (ERK1) and ERK2 are the downstream components of a phosphorelay pathway that conveys growth and mitogenic signals largely channelled by the small RAS GTPases. By phosphorylating widely diverse substrates, ERK proteins govern a variety of evolutionarily conserved cellular processes in metazoans, the dysregulation of which contributes to the cause of distinct human diseases. The mechanisms underlying the regulation of ERK1 and ERK2, their mode of action and their impact on the development and homeostasis of various organisms have been the focus of much attention for nearly three decades. In this Review, we discuss the current understanding of this important class of kinases. We begin with a brief overview of the structure, regulation, substrate recognition and subcellular localization of ERK1 and ERK2. We then systematically discuss how ERK signalling regulates six fundamental cellular processes in response to extracellular cues. These processes are cell proliferation, cell survival, cell growth, cell metabolism, cell migration and cell differentiation.

Immunostimulatory efficacy and protective potential of putative TgERK7 protein in mice experimentally infected by Toxoplasma gondii

The extracellular signal-regulated kinases (ERKs) serve as important determinants of cellular signal transduction pathways, and hence may play important roles during infections. Previous work suggested that putative ERK7 of Toxoplasma gondii is required for efficient intracellular replication of the parasite. However, the antigenic and immunostimulatory properties of TgERK7 protein remain unknown. The objective of this study was to produce a recombinant TgERK7 protein in vitro and to evaluate its effect on the induction of humoral and T cell-mediated immune responses against T. gondii infection in BALB/c mice. Immunization using TgERK7 mixed with Freund's adjuvants significantly increased the ratio of CD3e⁺CD4⁺ T/CD3e⁺CD8a⁺ T lymphocytes in spleen and elevated serum cytokines (IFN-γ, IL-2, IL-4, IL-10, IL-12p70, IL-23, MCP-1, and TNF-α) in immunized mice compared to control mice. On the contrary, immunization did not induce high levels of serum IgG antibodies. Five predicted peptides of TgERK7 were synthesized and conjugated with KLH and used to analyze the antibody specificity in the sera of immunized mice. We detected a progressive increase in the antibody level only against TgERK7 peptide A (DEVDKHVLRKYD). Antibody raised against this peptide significantly decreased intracellular proliferation of T. gondii in vitro, suggesting that peptide A can potentially induce a protective antibody response. We also showed that immunization improved the survival rate of mice challenged with a virulent strain and significantly reduced the parasite cyst burden within the brains of chronically infected mice. Our data show that TgERK7-based immunization induced TgERK7 peptide A-specific immune responses that can impart protective immunity against T. gondii infection. The therapeutic potential of targeting ERK7 signaling pathway for future toxoplasmosis treatment is warranted.

The bacterial metalloprotease NleD selectively cleaves mitogen-activated protein kinases that have high flexibility in their activation loop

Microbial pathogens often target the host mitogen-activated protein kinase (MAPK) network to suppress host immune responses. We previously identified a bacterial type III secretion system effector, termed NleD, a metalloprotease that inactivates MAPKs by specifically cleaving their activation loop. Here, we show that NleDs form a growing family of virulence factors harbored by human and plant pathogens as well as insect symbionts. These NleDs disable specifically Jun N-terminal kinases (JNKs) and p38s that are required for host immune response, whereas extracellular signal-regulated kinase (ERK), which is essential for host cell viability, remains intact. We investigated the mechanism that makes ERK resistant to NleD cleavage. Biochemical and structural analyses revealed that NleD exclusively targets activation loops with high conformational flexibility. Accordingly, NleD cleaved the flexible loops of JNK and p38 but not the rigid loop of ERK. Our findings elucidate a compelling mechanism of native substrate proteolysis that is promoted by entropy-driven specificity. We propose that such entropy-based selectivity is a general attribute of proteolytic enzymes.

Spotlight on the Ballet of Proteins: The Structural Dynamic Properties of Proteins Illuminated by Solution NMR

Solution NMR spectroscopy is a unique and powerful technique that has the ability to directly connect the structural dynamics of proteins in physiological conditions to their activity and function. Here, we summarize recent studies in which solution NMR contributed to the discovery of relationships between key dynamic properties of proteins and functional mechanisms in important biological systems. The capacity of NMR to quantify the dynamics of proteins over a range of time scales and to detect lowly populated protein conformations plays a critical role in its power to unveil functional protein dynamics. This analysis of dynamics is not only important for the understanding of biological function, but also in the design of specific ligands for pharmacologically important proteins. Thus, the dynamic view of structure provided by NMR is of importance in both basic and applied biology.

Structural Consequences of Multisite Phosphorylation in the BAK1 Kinase Domain

Multisite phosphorylation is an important mechanism of post-translational control of protein kinases. The effects of combinations of possible phosphorylation states on protein kinase activity are difficult to study experimentally because of challenges in isolating a particular phosphorylation state; surprising little effort on this topic has been expended in computational studies. To understand the effects of multisite phosphorylation on the plant protein kinase brassinosteroid insensitive 1-associated kinase 1 (BAK1) conformational ensemble, we performed Gaussian accelerated molecular dynamics simulations on eight BAK1 mod-forms involving phosphorylation of the four activation-loop threonine residues and binding of ATP-Mg²⁺. We find that unphosphorylated BAK1 transitions into an inactive conformation with a "cracked" activation loop and with the αC helix swung away from the active site. T450 phosphorylation can prevent the activation loop from cracking and keep the αC helix in an active-like conformation, whereas phosphorylation of T455 only slightly stabilizes the activation loop. There is a general trend of reduced flexibility in interlobe motion with increased phosphorylation. Interestingly, the αC helix is destabilized when the activation loop is fully phosphorylated but is again stabilized with ATP-Mg²⁺ bound. Our results provide insight into the mechanism of phosphorylation-controlled BAK1 activation while at the same time represent the first, to our knowledge, comprehensive, comparative study of the effects of combinatorial phosphorylation states on protein kinase conformational dynamics.

A Dynamic Switch in Inactive p38γ Leads to an Excited State on the Pathway to an Active Kinase

The inactive state of mitogen-activated protein kinases (MAPKs) adopts an open conformation while the active state exists in a compact form stabilized by phosphorylation. In the active state, eukaryotic kinases undergo breathing motions related to substrate binding and product release that have not previously been detected in the inactive state. However, docking interactions of partner proteins with inactive MAPK kinases exhibit allostery in binding of activating kinases. Interactions at a site distant from the activation loop are coupled to the configuration of the activation loop, suggesting that the inactive state may also undergo concerted dynamics. X-ray crystallographic studies of nonphosphorylated, inactive p38γ reveal differences in domain orientations and active site structure in the two molecules in the asymmetric unit. One molecule resembles an inactive kinase with an open active site. The second molecule has a rotation of the N-lobe that leads to partial compaction of the active site, resulting in a conformation that is intermediate between the inactive open state and the fully closed state of the activated kinase. Although the compact state of apo p38γ displays several of the features of the activated enzyme, it remains catalytically inert. In solution, the kinase fluctuates on a millisecond time scale between the open ground state and a weakly populated excited state that is similar in structure to the compact state observed in the crystal. The nuclear magnetic resonance and crystal structure data imply that interconversion between the open and compact states involves a molecular switch associated with the DFG loop.

Kinase Activation by Small Conformational Changes

Protein kinases (PKs) are allosteric enzymes that play an essential role in signal transduction by regulating a variety of key cellular processes. Most PKs suffer conformational rearrangements upon phosphorylation that strongly enhance the catalytic activity. Generally, it involves the movement of the phosphorylated loop toward the active site and the rotation of the whole C-terminal lobe. However, not all kinases undergo such a large configurational change: The MAPK extracellular signal-regulated protein kinases ERK1 and ERK2 achieve a 50 000 fold increase in kinase activity with only a small motion of the C-terminal region. In the present work, we used a combination of molecular simulation tools to characterize the conformational landscape of ERK2 in the active (phosphorylated) and inactive (unphosphorylated) states in solution in agreement with NMR experiments. We show that the chemical reaction barrier is strongly dependent on ATP conformation and that the "active" low-barrier configuration is subtly regulated by phosphorylation, which stabilizes a key salt bridge between the conserved Lys52 and Glu69 belonging to helix-C and promotes binding of a second Mg ion. Our study highlights that the on-off switch embedded in the kinase fold can be regulated by small, medium, and large conformational changes.

Ancestral reconstruction reveals mechanisms of ERK regulatory evolution

Protein kinases are crucial to coordinate cellular decisions and therefore their activities are strictly regulated. Previously we used ancestral reconstruction to determine how CMGC group kinase specificity evolved (Howard et al., 2014). In the present study, we reconstructed ancestral kinases to study the evolution of regulation, from the inferred ancestor of CDKs and MAPKs, to modern ERKs. Kinases switched from high to low autophosphorylation activity at the transition to the inferred ancestor of ERKs 1 and 2. Two synergistic amino acid changes were sufficient to induce this change: shortening of the β3-αC loop and mutation of the gatekeeper residue. Restoring these two mutations to their inferred ancestral state led to a loss of dependence of modern ERKs 1 and 2 on the upstream activating kinase MEK in human cells. Our results shed light on the evolutionary mechanisms that led to the tight regulation of a kinase that is central in development and disease.

Functional divergence caused by mutations in an energetic hotspot in ERK2

The most frequent extracellular signal-regulated kinase 2 (ERK2) mutation occurring in cancers is E322K (E-K). ERK2 E-K reverses a buried charge in the ERK2 common docking (CD) site, a region that binds activators, inhibitors, and substrates. Little is known about the cellular consequences associated with this mutation, other than apparent increases in tumor resistance to pathway inhibitors. ERK2 E-K, like the mutation of the preceding aspartate (ERK2 D321N [D-N]) known as the sevenmaker mutation, causes increased activity in cells and evades inactivation by dual-specificity phosphatases. As opposed to findings in cancer cells, in developmental assays in Drosophila, only ERK2 D-N displays a significant gain of function, revealing mutation-specific phenotypes. The crystal structure of ERK2 D-N is indistinguishable from that of wild-type protein, yet this mutant displays increased thermal stability. In contrast, the crystal structure of ERK2 E-K reveals profound structural changes, including disorder in the CD site and exposure of the activation loop phosphorylation sites, which likely account for the decreased thermal stability of the protein. These contiguous mutations in the CD site of ERK2 are both required for docking interactions but lead to unpredictably different functional outcomes. Our results suggest that the CD site is in an energetically strained configuration, and this helps drive conformational changes at distal sites on ERK2 during docking interactions.

Activation loop dynamics are controlled by conformation-selective inhibitors of ERK2

Conformational selection by small molecules expands inhibitory possibilities for protein kinases. Nuclear magnetic resonance (NMR) measurements of the mitogen-activated protein (MAP) kinase ERK2 have shown that activation by dual phosphorylation induces global motions involving exchange between two states, L and R. We show that ERK inhibitors Vertex-11e and SCH772984 exploit the small energetic difference between L and R to shift the equilibrium in opposing directions. An X-ray structure of active 2P-ERK2 complexed with AMP-PNP reveals a shift in the Gly-rich loop along with domain closure to position the nucleotide in a more catalytically productive conformation relative to inactive 0P-ERK2:ATP. X-ray structures of 2P-ERK2 complexed with Vertex-11e or GDC-0994 recapitulate this closure, which is blocked in a complex with a SCH772984 analog. Thus, the L→R shift in 2P-ERK2 is associated with movements needed to form a competent active site. Solution measurements by hydrogen-exchange mass spectrometry (HX-MS) reveal distinct binding interactions for Vertex-11e, GDC-0994, and AMP-PNP with active vs. inactive ERK2, where the extent of HX protection correlates with R state formation. Furthermore, Vertex-11e and SCH772984 show opposite effects on HX near the activation loop. Consequently, these inhibitors differentially affect MAP kinase phosphatase activity toward 2P-ERK2. We conclude that global motions in ERK2 reflect conformational changes at the active site that promote productive nucleotide binding and couple with changes at the activation loop to allow control of dephosphorylation by conformationally selective inhibitors.

A Novel Class of Common Docking Domain Inhibitors That Prevent ERK2 Activation and Substrate Phosphorylation

Extracellular signal-regulated kinases (ERK1/2) are mitogen-activated protein kinases (MAPKs) that play a pro-tumorigenic role in numerous cancers. ERK1/2 possess two protein-docking sites that are distinct from the active site: the D-recruitment site (DRS) and the F-recruitment site. These docking sites facilitate substrate recognition, intracellular localization, signaling specificity, and protein complex assembly. Targeting these sites on ERK in a therapeutic context may overcome many problems associated with traditional ATP-competitive inhibitors. Here, we identified a new class of inhibitors that target the ERK DRS by screening a synthetic combinatorial library of more than 30 million compounds. The screen detects the competitive displacement of a fluorescent peptide from the DRS of ERK2. The top molecular scaffold from the screen was optimized for structure-activity relationship by positional scanning of different functional groups. This resulted in 10 compounds with similar binding affinities and a shared core structure consisting of a tertiary amine hub with three functionalized cyclic guanidino branches. Compound 2507-1 inhibited ERK2 from phosphorylating a DRS-targeting substrate and prevented the phosphorylation of ERK2 by a constitutively active MEK1 (MAPK/ERK kinase 1) mutant. Interaction between an analogue, 2507-8, and the ERK2 DRS was confirmed by nuclear magnetic resonance and X-ray crystallography. 2507-8 forms critical interactions at the common docking domain residue Asp319 via an arginine-like moiety that is shared by all 10 hits, suggesting a common binding mode. The structural and biochemical insights reported here provide the basis for developing new ERK inhibitors that are not ATP-competitive but instead function by disrupting critical protein-protein interactions.

Targeting ERK beyond the boundaries of the kinase active site in melanoma

Extracellular signal-regulated kinase 1/2 (ERK1/2) constitute a point of convergence for complex signaling events that regulate essential cellular processes, including proliferation and survival. As such, dysregulation of the ERK signaling pathway is prevalent in many cancers. In the case of BRAF-V600E mutant melanoma, ERK inhibition has emerged as a viable clinical approach to abrogate signaling through the ERK pathway, even in cases where MEK and Raf inhibitor treatments fail to induce tumor regression due to resistance mechanisms. Several ERK inhibitors that target the active site of ERK have reached clinical trials, however, many critical ERK interactions occur at other potentially druggable sites on the protein. Here we discuss the role of ERK signaling in cell fate, in driving melanoma, and in resistance mechanisms to current BRAF-V600E melanoma treatments. We explore targeting ERK via a distinct site of protein-protein interaction, known as the D-recruitment site (DRS), as an alternative or supplementary mode of ERK pathway inhibition in BRAF-V600E melanoma. Targeting the DRS with inhibitors in melanoma has the potential to not only disrupt the catalytic apparatus of ERK but also its noncatalytic functions, which have significant impacts on spatiotemporal signaling dynamics and cell fate.

Kinetic network model to explain gain-of-function mutations in ERK2 enzyme

ERK2 is a kinase protein that belongs to a Ras/Raf/MEK/ERK signaling pathway, which is activated in response to a range of extracellular signals. Malfunctioning of this cascade leads to a variety of serious diseases, including cancers. This is often caused by mutations in proteins belonging to the cascade, frequently leading to abnormally high activity of the cascade even in the absence of an external signal. One such "gain-of-function" mutation in the ERK2 protein, called a "sevenmaker" mutation (D319N), was discovered in 1994 in Drosophila. The mutation leads to disruption of interactions of other proteins with the D-site of ERK2 and results, contrary to expectations, in an increase of its activity in vivo. However, no molecular mechanism to explain this effect has been presented so far. The difficulty is that this mutation should equally negatively affect interactions of ERK2 with all substrates, activators, and deactivators. In this paper, we present a semiquantitative kinetic network model that gives a possible explanation of the increased activity of mutant ERK2 species. A simplified biochemical network for ERK2, viewed as a system of coupled Michaelis-Menten processes, is presented. Its dynamic properties are calculated explicitly using the method of first-passage processes. The effect of mutation is associated with changes in the strength of interaction energy between the enzyme and the substrates. It is found that the dependence of kinetic properties of the protein on the interaction energy is nonmonotonic, suggesting that some mutations might lead to more efficient catalytic properties, despite weakening intermolecular interactions. Our theoretical predictions agree with experimental observations for the sevenmaker mutation in ERK2. It is also argued that the effect of mutations might depend on the concentrations of substrates.

Disordered Protein Kinase Regions in Regulation of Kinase Domain Cores

Since publication of the crystal structure of protein kinase (PK)A three decades ago, a structural portrait of the conserved kinase core has been drawn. The next challenge is to elucidate structures of full-length kinases and to address the intrinsically disordered regions (IDRs) that typically flank the core as well as the small linear motifs (SLiMs) that are embedded within the IDRs. It is increasingly apparent that unstructured regions integrate the kinase catalytic chassis into multienzyme-based regulatory networks. The extracellular signal-regulated kinase-ribosomal S6 PK-phosphoinositide-dependent kinase (ERK-RSK-PDK) complex is an excellent example to demonstrate how IDRs and SLiMs govern communication between four different kinase catalytic cores to mediate activation and how in molecular terms these promote the formation of kinase heterodimers in a context dependent fashion.

Targeting ERK1/2 protein-serine/threonine kinases in human cancers

ERK1 and ERK2 are key protein kinases that contribute to the Ras-Raf-MEK-ERK MAP kinase signalling module. This pathway participates in the control of numerous processes including apoptosis, cell proliferation, the immune response, nervous system function, and RNA synthesis and processing. MEK1/2 activate human ERK1/2 by first catalyzing the phosphorylation of Y204/187 and then T202/185, both residues of which occur within the activation segment. The phosphorylation of both residues is required for enzyme activation. The only Raf substrates are MEK1/2 and the only MEK1/2 substrates are ERK1/2. In contrast, ERK1/2 catalyze the phosphorylation of many cytoplasmic and nuclear substrates including transcription factors and regulatory molecules. The linear MAP kinase pathway branches extensively at the ERK1/2 node. ERK1/2 are proline-directed kinases that preferentially catalyze the phosphorylation of substrates containing a PxS/TP sequence. The dephosphorylation and inactivation of ERK1/2 is catalyzed by dual specificity phosphatases, protein-tyrosine specific phosphatases, and protein-serine/threonine phosphatases. The combined functions of kinases and phosphatases make the overall process reversible. To provide an idea of the complexities involved in these reactions, somatic cell cycle progression involves the strict timing of more than 32,000 phosphorylation and dephosphorylation events as determined by mass spectrometry. The MAP kinase cascade is perhaps the most important oncogenic driver of human cancers and the blockade of this signalling module by targeted inhibitors is an important anti-tumor strategy. Although numerous cancers are driven by MAP kinase pathway activation, thus far the only orally effective approved drugs that target this signaling module are used for the treatment of BRAF-mutant melanomas. The best treatments include the combination of B-Raf and MEK inhibitors (dabrafenib and trametinib, encorafenib and binimetinib, vemurafenib and cobimetanib). However, resistance to these antagonists occurs within one year and additional treatment options are necessary. Owing to the large variety of malignancies that are driven by dysregulation of the MAP kinase pathway, additional tumor types should be amenable to MAP kinase pathway inhibitor therapy. In addition to new B-Raf and MEK inhibitors, the addition of ERK inhibitors should prove helpful. Ulixertinib, MK-8353, and GDC-0994 are orally effective, potent, and specific inhibitors of ERK1/2 that are in early clinical trials for the treatment of various advanced/metastatic solid tumors. These agents are effective against cell lines that are resistant to B-Raf and MEK1/2 inhibitor therapy. Although MK-8353 does not directly inhibit MEK1/2, it decreases the phosphorylation of ERK1/2 as well as the phosphorylation of RSK, an ERK1/2 substrate. The decrease in RSK phosphorylation appears to be a result of ERK inhibition and the decrease in ERK1/2 phosphorylation is related to the inability of MEK to catalyze the phosphorylation of the ERK-MK-8353 complex; these decreases characterize the ERK dual mechanism inhibition paradigm. Additional work will be required to determine whether ERK inhibitors will be successful in the clinic and are able to forestall the development of drug resistance of the MAP kinase pathway.

Arrestin-3 scaffolding of the JNK3 cascade suggests a mechanism for signal amplification

Scaffold proteins tether and orient components of a signaling cascade to facilitate signaling. Although much is known about how scaffolds colocalize signaling proteins, it is unclear whether scaffolds promote signal amplification. Here, we used arrestin-3, a scaffold of the ASK1-MKK4/7-JNK3 cascade, as a model to understand signal amplification by a scaffold protein. We found that arrestin-3 exhibited >15-fold higher affinity for inactive JNK3 than for active JNK3, and this change involved a shift in the binding site following JNK3 activation. We used systems biochemistry modeling and Bayesian inference to evaluate how the activation of upstream kinases contributed to JNK3 phosphorylation. Our combined experimental and computational approach suggested that the catalytic phosphorylation rate of JNK3 at Thr-221 by MKK7 is two orders of magnitude faster than the corresponding phosphorylation of Tyr-223 by MKK4 with or without arrestin-3. Finally, we showed that the release of activated JNK3 was critical for signal amplification. Collectively, our data suggest a "conveyor belt" mechanism for signal amplification by scaffold proteins. This mechanism informs on a long-standing mystery for how few upstream kinase molecules activate numerous downstream kinases to amplify signaling.

Muscarinic receptor antagonists activate ERK-CREB signaling to augment neurite outgrowth of adult sensory neurons

A major cellular effector activated by G protein coupled receptors is extracellular signal-regulated kinase (ERK). The ERK signaling cascade regulates a variety of cellular processes including growth and proliferation. Both G protein and β-arrestin-mediated signaling lead to ERK activation by phosphorylation through different kinases. Recently, we have shown muscarinic acetylcholine type 1 receptor (M₁R) antagonists, muscarinic toxin 7 (MT7) and pirenzepine, elevated neurite outgrowth and protected from small and large fiber neuropathy in adult sensory neurons in various animal models. Thus, we tested the novel hypothesis that muscarinic antagonists could drive neurite outgrowth through altered M₁R-ERK signaling. We have used two dimensional isoelectric focusing/SDS-PAGE combined with analysis using multiple phospho-epitope specific antibodies to study ERK1/2 phosphorylation and activation of its downstream nuclear effector cyclic response element binding protein (CREB). Activated CREB is known to exhibit neuroprotective and growth promoting effects. One hour of treatment with MT7 and pirenzepine activated ERK through M₁R and induced a significant increase in levels of pCREB(S133) in cultured sensory neurons. Further, pharmacological blockade or siRNA based knockdown of ERK abolished the MT7 and pirenzepine mediated neuritogenic effect. In addition, we have shown drug-induced alterations of charged protein fractions that may possess additional post-translationally modified forms of ERK and CREB. For the first time we show that long-term treatment, e.g. 1 h, with muscarinic antagonists selective or specific for M₁R can activate a biased β-arrestin dependent ERK-CREB signal cascade. Our study gives novel insight into muscarinic antagonist-mediated modulation of M₁R-ERK-CREB signaling which could be exploited for therapy in neuropathic diseases.

H2O2 oxidation of cysteine residues in c-Jun N-terminal kinase 2 (JNK2) contributes to redox regulation in human articular chondrocytes

Reactive oxygen species (ROS), in particular H₂O₂, regulate intracellular signaling through reversible oxidation of reactive protein thiols present in a number of kinases and phosphatases. H₂O₂ has been shown to regulate mitogen-activated protein kinase (MAPK) signaling depending on the cellular context. We report here that in human articular chondrocytes, the MAPK family member c-Jun N-terminal kinase 2 (JNK2) is activated by fibronectin fragments and low physiological levels of H₂O₂ and inhibited by oxidation due to elevated levels of H₂O₂ The kinase activity of affinity-purified, phosphorylated JNK2 from cultured chondrocytes was reversibly inhibited by 5-20 μm H₂O₂ Using dimedone-based chemical probes that react specifically with sulfenylated cysteines (RSOH), we identified Cys-222 in JNK2, a residue not conserved in JNK1 or JNK3, as a redox-reactive site. MS analysis of human recombinant JNK2 also detected further oxidation at Cys-222 and other cysteines to sulfinic (RSO₂H) or sulfonic (RSO₃H) acid. H₂O₂ treatment of JNK2 resulted in detectable levels of peptides containing intramolecular disulfides between Cys-222 and either Cys-213 or Cys-177, without evidence of dimer formation. Substitution of Cys-222 to alanine rendered JNK2 insensitive to H₂O₂ inhibition, unlike C177A and C213A variants. Two other JNK2 variants, C116A and C163A, were also resistant to oxidative inhibition. Cumulatively, these findings indicate differential regulation of JNK2 signaling dependent on H₂O₂ levels and point to key cysteine residues regulating JNK2 activity. As levels of intracellular H₂O₂ rise, a switch occurs from activation to inhibition of JNK2 activity, linking JNK2 regulation to the redox status of the cell.

Endogenous, regulatory cysteine sulfenylation of ERK kinases in response to proliferative signals

ERK-dependent signaling is key to many pathways through which extracellular signals are transduced into cell-fate decisions. One conundrum is the way in which disparate signals induce specific responses through a common, ERK-dependent kinase cascade. While studies have revealed intricate ways of controlling ERK signaling through spatiotemporal localization and phosphorylation dynamics, additional modes of ERK regulation undoubtedly remain to be discovered. We hypothesized that fine-tuning of ERK signaling could occur by cysteine oxidation. We report that ERK is actively and directly oxidized by signal-generated H₂O₂ during proliferative signaling, and that ERK oxidation occurs downstream of a variety of receptor classes tested in four cell lines. Furthermore, within the tested cell lines and proliferative signals, we observed that both activation loop-phosphorylated and non-phosphorylated ERK undergo sulfenylation in cells and that dynamics of ERK sulfenylation is dependent on the cell growth conditions prior to stimulation. We also tested the effect of endogenous ERK oxidation on kinase activity and report that phosphotransfer reactions are reversibly inhibited by oxidation by as much as 80-90%, underscoring the importance of considering this additional modification when assessing ERK activation in response to extracellular signals.

Asp263 missense variants perturb the active site of human phosphoglucomutase 1

The enzyme phosphoglucomutase 1 (PGM1) plays a central role in glucose homeostasis. Clinical studies have identified mutations in human PGM1 as the cause of PGM1 deficiency, an inherited metabolic disease. One residue, Asp263, has two known variants associated with disease: D263G and D263Y. Biochemical studies have shown that these mutants are soluble and well folded, but have significant catalytic impairment. To better understand this catalytic defect, we determined crystal structures of these two missense variants, both of which reveal a similar and indirect structural change due to the loss of a conserved salt bridge between Asp263 and Arg293. The arginine reorients into the active site, making interactions with residues responsible for substrate binding. Biochemical studies also show that the catalytic phosphoserine of the missense variants is more stable to hydrolysis relative to wild-type enzyme. The structural perturbation resulting from mutation of this single amino acid reveals the molecular mechanism underlying PGM1 deficiency in these missense variants.

DATABASE

Structural data are available in the PDB under the accession numbers 5JN5 and 5TR2.

Determining time-weighted average concentrations of nitrate and ammonium in freshwaters using DGT with ion exchange membrane-based binding layers

Commercially-available AMI-7001 anion exchange and CMI-7000 cation exchange membranes were utilised as binding layers for DGT measurements of NO₃-N and NH₄-N in freshwaters. These ion exchange membranes are easier to prepare and handle than DGT binding layers consisting of hydrogels cast with ion exchange resins. The membranes showed good uptake and elution efficiencies for both NO₃-N and NH₄-N. The membrane-based DGTs are suitable for pH 3.5-8.5 and ionic strength ranges (0.0001-0.014 and 0.0003-0.012 mol L^-1 as NaCl for the AMI-7001 and CMI-7000 membrane, respectively) typical of most natural freshwaters. The binding membranes had high intrinsic binding capacities for NO₃-N and NH₄-N of 911 ± 88 μg and 3512 ± 51 μg, respectively. Interferences from the major competing ions for membrane-based DGTs are similar to DGTs employing resin-based binding layers but with slightly different selectivity. This different selectivity means that the two DGT types can be used in different types of freshwaters. The laboratory and field experiments demonstrated that AMI-DGT and CMI-DGT can be an alternative to A520E-DGT and PrCH-DGT for measuring NO₃-N and NH₄-N, respectively, as (i) membrane-based DGT have a consistent composition, (ii) avoid the use of toxic chemicals, (iii) provided highly representative results (C_DGT : C_SOLN between 0.81 and 1.3), and (iv) agreed with resin-based DGTs to within 85-120%.