ERK as a model for systems biology of enzyme kinetics in cells

Determining the ERK-regulated phosphoproteome driving KRAS-mutant cancer

To delineate the mechanisms by which the ERK1 and ERK2 mitogen-activated protein kinases support mutant KRAS-driven cancer growth, we determined the ERK-dependent phosphoproteome in KRAS-mutant pancreatic cancer. We determined that ERK1 and ERK2 share near-identical signaling and transforming outputs and that the KRAS-regulated phosphoproteome is driven nearly completely by ERK. We identified 4666 ERK-dependent phosphosites on 2123 proteins, of which 79 and 66%, respectively, were not previously associated with ERK, substantially expanding the depth and breadth of ERK-dependent phosphorylation events and revealing a considerably more complex function for ERK in cancer. We established that ERK controls a highly dynamic and complex phosphoproteome that converges on cyclin-dependent kinase regulation and RAS homolog guanosine triphosphatase function (RHO GTPase). Our findings establish the most comprehensive molecular portrait and mechanisms by which ERK drives KRAS-dependent pancreatic cancer growth.

An Affordable Topography-Based Protocol for Assigning a Residue's Character on a Hydropathy (PARCH) Scale

The hydropathy of proteins or quantitative assessment of protein-water interactions has been a topic of interest for decades. Most hydropathy scales use a residue-based or atom-based approach to assign fixed numerical values to the 20 amino acids and categorize them as hydrophilic, hydroneutral, or hydrophobic. These scales overlook the protein's nanoscale topography, such as bumps, crevices, cavities, clefts, pockets, and channels, in calculating the hydropathy of the residues. Some recent studies have included protein topography in determining hydrophobic patches on protein surfaces, but these methods do not provide a hydropathy scale. To overcome the limitations in the existing methods, we have developed a Protocol for Assigning a Residue's Character on the Hydropathy (PARCH) scale that adopts a holistic approach to assigning the hydropathy of a residue. The parch scale evaluates the collective response of the water molecules in the protein's first hydration shell to increasing temperatures. We performed the parch analysis of a set of well-studied proteins that include the following─enzymes, immune proteins, and integral membrane proteins, as well as fungal and virus capsid proteins. Since the parch scale evaluates every residue based on its location, a residue may have very different parch values inside a crevice versus a surface bump. Thus, a residue can have a range of parch values (or hydropathies) dictated by the local geometry. The parch scale calculations are computationally inexpensive and can compare hydropathies of different proteins. The parch analysis can affordably and reliably aid in designing nanostructured surfaces, identifying hydrophilic and hydrophobic patches, and drug discovery.

Phase plane dynamics of ERK phosphorylation

The extracellular signal-regulated kinase (ERK) controls multiple critical processes in the cell and is deregulated in human cancers, congenital abnormalities, immune diseases, and neurodevelopmental syndromes. Catalytic activity of ERK requires dual phosphorylation by an upstream kinase, in a mechanism that can be described by two sequential Michaelis-Menten steps. The estimation of individual reaction rate constants from kinetic data in the full mechanism has proved challenging. Here, we present an analytically tractable approach to parameter estimation that is based on the phase plane representation of ERK activation and yields two combinations of six reaction rate constants in the detailed mechanism. These combinations correspond to the ratio of the specificities of two consecutive phosphorylations and the probability that monophosphorylated substrate does not dissociate from the enzyme before the second phosphorylation. The presented approach offers a language for comparing the effects of mutations that disrupt ERK activation and function in vivo. As an illustration, we use phase plane representation to analyze dual phosphorylation under heterozygous conditions, when two enzyme variants compete for the same substrate.

Stochastic phenotypes in RAS-dependent developmental diseases

Germline mutations upregulating RAS signaling are associated with multiple developmental disorders. A hallmark of these conditions is that the same mutation may present vastly different phenotypes in different individuals, even in monozygotic twins. Here, we demonstrate how the origins of such largely unexplained phenotypic variations may be dissected using highly controlled studies in Drosophila that have been gene edited to carry activating variants of MEK, a core enzyme in the RAS pathway. This allowed us to measure the small but consistent increase in signaling output of such alleles in vivo. The fraction of mutation carriers reaching adulthood was strongly reduced, but most surviving animals had normal RAS-dependent structures. We rationalize these results using a stochastic signaling model and support it by quantifying cell fate specification errors in bilaterally symmetric larval trachea, a RAS-dependent structure that allows us to isolate the effects of mutations from potential contributions of genetic modifiers and environmental differences. We propose that the small increase in signaling output shifts the distribution of phenotypes into a regime, where stochastic variation causes defects in some individuals, but not in others. Our findings shed light on phenotypic heterogeneity of developmental disorders caused by deregulated RAS signaling and offer a framework for investigating causal effects of other pathogenic alleles and mild mutations in general.

hucMSC-Ex Alleviates IBD-Associated Intestinal Fibrosis by Inhibiting ERK Phosphorylation in Intestinal Fibroblasts

Background

Intestinal fibrosis, one of the complications of inflammatory bowel disease (IBD), is associated with fistula and intestinal stricture formation. There are currently no treatments for fibrosis. Mesenchymal stem cell-derived exosomes have been proven to exert inhibitory and reversal effects in IBD and other organ fibrosis. In this study, we explored the role of human umbilical cord mesenchymal stem cell-derived exosomes (hucMSC-Ex) in IBD-related fibrosis and its associated mechanism to provide new ideas for the prevention and treatment of IBD-related intestinal fibrosis.

Methods

We established a DSS-induced mouse IBD-related intestinal fibrosis model and observed the effect of hucMSC-Ex on the mouse model. We also used the TGF-induced human intestinal fibroblast CCD-18Co to observe the role of hucMSC-Ex in the proliferation, migration, and activation of intestinal fibroblasts. Having observed that the extracellular-signal-regulated kinase (ERK) pathway in intestinal fibrosis can be inhibited by hucMSC-Ex, we treated intestinal fibroblasts with an ERK inhibitor to emphasize the potential target of ERK phosphorylation in the treatment of IBD-associated intestinal fibrosis.

Results

In the animal model of IBD-related fibrosis, hucMSC-Ex alleviated inflammation-related fibrosis as evident in the thinning of the mice's intestinal wall and decreased expression of related molecules. Moreover, hucMSC-Ex inhibited TGF-β-induced proliferation, migration, and activation of human intestinal fibroblasts, and ERK phosphorylation played a key role in IBD-associated fibrosis. The inhibition of ERK decreased the expression of fibrosis-related indicators such as α-SMA, fibronectin, and collagen I.

Conclusion

hucMSC-Ex alleviates DSS-induced IBD-related intestinal fibrosis by inhibiting profibrotic molecules and intestinal fibroblast proliferation and migration by decreasing ERK phosphorylation.

ERK2 MAP kinase regulates SUFU binding by multisite phosphorylation of GLI1

Crosstalk between the Hedgehog and MAPK signaling pathways occurs in several types of cancer and contributes to clinical resistance to Hedgehog pathway inhibitors. Here we show that MAP kinase-mediated phosphorylation weakens the binding of the GLI1 transcription factor to its negative regulator SUFU. ERK2 phosphorylates GLI1 on three evolutionarily conserved target sites (S102, S116, and S130) located near the high-affinity binding site for SUFU; these phosphorylations cooperate to weaken the affinity of GLI1-SUFU binding by over 25-fold. Phosphorylation of any one, or even any two, of the three sites does not result in the level of SUFU release seen when all three sites are phosphorylated. Tumor-derived mutations in R100 and S105, residues bordering S102, also diminish SUFU binding, collectively defining a novel evolutionarily conserved SUFU affinity-modulating region. In cultured mammalian cells, GLI1 variants containing phosphomimetic substitutions of S102, S116, and S130 displayed an increased ability to drive transcription. We conclude that multisite phosphorylation of GLI1 by ERK2 or other MAP kinases weakens GLI1-SUFU binding, thereby facilitating GLI1 activation and contributing to both physiological and pathological crosstalk.

Association of two polyelectrolytes in salt solutions

The association of polyelectrolytes (PEs) in solution affects a wealth of structural and dynamic behaviors, and is also fundamentally important for an understanding of protein association and aggregation. Here, we theoretically study the association of two PE chains by addressing the stability and morphology of the non-spherical associates. Our theory predicts that an elongated pearl-necklace (PN) associate can be stable at high salt concentrations due to the screened electrostatic repulsion. This contradicts the implication of scaling theory. In addition, there is no one-to-one correspondence between the morphology of the associate and its constituting unimers, which is demonstrated by the existence of different association modes.

Identifying hydrophobic protein patches to inform protein interaction interfaces

Interactions between proteins lie at the heart of numerous biological processes and are essential for the proper functioning of the cell. Although the importance of hydrophobic residues in driving protein interactions is universally accepted, a characterization of protein hydrophobicity, which informs its interactions, has remained elusive. The challenge lies in capturing the collective response of the protein hydration waters to the nanoscale chemical and topographical protein patterns, which determine protein hydrophobicity. To address this challenge, here, we employ specialized molecular simulations wherein water molecules are systematically displaced from the protein hydration shell; by identifying protein regions that relinquish their waters more readily than others, we are then able to uncover the most hydrophobic protein patches. Surprisingly, such patches contain a large fraction of polar/charged atoms and have chemical compositions that are similar to the more hydrophilic protein patches. Importantly, we also find a striking correspondence between the most hydrophobic protein patches and regions that mediate protein interactions. Our work thus establishes a computational framework for characterizing the emergent hydrophobicity of amphiphilic solutes, such as proteins, which display nanoscale heterogeneity, and for uncovering their interaction interfaces.

Dynamics of ERK regulation in the processive limit

We consider a model of extracellular signal-regulated kinase regulation by dual-site phosphorylation and dephosphorylation, which exhibits bistability and oscillations, but loses these properties in the limit in which the mechanisms underlying phosphorylation and dephosphorylation become processive. Our results suggest that anywhere along the way to becoming processive, the model remains bistable and oscillatory. More precisely, in simplified versions of the model, precursors to bistability and oscillations (specifically, multistationarity and Hopf bifurcations, respectively) exist at all "processivity levels". Finally, we investigate whether bistability and oscillations can exist together.

Molecular mechanisms underlying cellular effects of human MEK1 mutations

Terminal regions of Drosophila embryos are patterned by signaling through ERK, which is genetically deregulated in multiple human diseases. Quantitative studies of terminal patterning have been recently used to investigate gain-of-function variants of human MEK1, encoding the MEK kinase that directly activates ERK by dual phosphorylation. Unexpectedly, several mutations reduced ERK activation by extracellular signals, possibly through a negative feedback triggered by signal-independent activity of the mutant variants. Here we present experimental evidence supporting this model. Using a MEK variant that combines a mutation within the negative regulatory region with alanine substitutions in the activation loop, we prove that pathogenic variants indeed acquire signal-independent kinase activity. We also demonstrate that signal-dependent activation of these variants is independent of kinase suppressor of Ras, a conserved adaptor that is indispensable for activation of normal MEK. Finally, we show that attenuation of ERK activation by extracellular signals stems from transcriptional induction of Mkp3, a dual specificity phosphatase that deactivates ERK by dephosphorylation. These findings in the Drosophila embryo highlight its power for investigating diverse effects of human disease mutations.

Inference of Multisite Phosphorylation Rate Constants and Their Modulation by Pathogenic Mutations

Multisite protein phosphorylation plays a critical role in cell regulation [1-3]. It is widely appreciated that the functional capabilities of multisite phosphorylation depend on the order and kinetics of phosphorylation steps, but kinetic aspects of multisite phosphorylation remain poorly understood [4-6]. Here, we focus on what appears to be the simplest scenario, when a protein is phosphorylated on only two sites in a strict, well-defined order. This scenario describes the activation of ERK, a highly conserved cell-signaling enzyme. We use Bayesian parameter inference in a structurally identifiable kinetic model to dissect dual phosphorylation of ERK by MEK, a kinase that is mutated in a large number of human diseases [7-12]. Our results reveal how enzyme processivity and efficiencies of individual phosphorylation steps are altered by pathogenic mutations. The presented approach, which connects specific mutations to kinetic parameters of multisite phosphorylation mechanisms, provides a systematic framework for closing the gap between studies with purified enzymes and their effects in the living organism.

The design and logic of terminal patterning in Drosophila

Terminal regions of the early Drosophila embryo are patterned by the highly conserved ERK cascade, giving rise to the nonsegmented terminal structures of the future larva. In less than an hour, this signaling event establishes several gene expression boundaries and sets in motion a sequence of elaborate morphogenetic events. Genetic studies of terminal patterning discovered signaling components and transcription factors that are involved in numerous developmental contexts and deregulated in human diseases. This review summarizes current understanding of signaling and morphogenesis during terminal patterning and discusses several open questions that can now be rigorously investigated using live imaging, omics, and optogenetic approaches. The anatomical simplicity of the terminal patterning system and its amenability to a broad range of increasingly sophisticated genetic perturbations will continue to make it a premier quantitative model for studying multiple aspects of tissue patterning by dynamically controlled cell signaling pathways.

Activation-induced substrate engagement in ERK signaling

The extracellular signal-regulated kinase (ERK) pathway is an essential component of developmental signaling in metazoans. Previous models of pathway activation suggested that dissociation of activated dually phosphorylated ERK (dpERK) from MAPK/ERK kinase (MEK), a kinase that phosphorylates ERK, and other cytoplasmic anchors, is sufficient for allowing ERK interactions with its substrates. Here, we provide evidence for an additional step controlling ERK’s access to substrates. Specifically, we demonstrate that interaction of ERK with its substrate Capicua (Cic) is controlled at the level of ERK phosphorylation, whereby Cic binds to dpERK much stronger than to unphosphorylated ERK, both in vitro and in vivo. Mathematical modeling suggests that the differential affinity of Cic for dpERK versus ERK is required for both down-regulation of Cic and stabilizing phosphorylated ERK. Preferential association of Cic with dpERK serves two functions: it prevents unproductive competition of Cic with unphosphorylated ERK and contributes to efficient signal propagation. We propose that high-affinity substrate binding increases the specificity and efficiency of signal transduction through the ERK pathway.

Oscillations and bistability in a model of ERK regulation

This work concerns the question of how two important dynamical properties, oscillations and bistability, emerge in an important biological signaling network. Specifically, we consider a model for dual-site phosphorylation and dephosphorylation of extracellular signal-regulated kinase (ERK). We prove that oscillations persist even as the model is greatly simplified (reactions are made irreversible and intermediates are removed). Bistability, however, is much less robust-this property is lost when intermediates are removed or even when all reactions are made irreversible. Moreover, bistability is characterized by the presence of two reversible, catalytic reactions: as other reactions are made irreversible, bistability persists as long as one or both of the specified reactions is preserved. Finally, we investigate the maximum number of steady states, aided by a network's "mixed volume" (a concept from convex geometry). Taken together, our results shed light on the question of how oscillations and bistability emerge from a limiting network of the ERK network-namely, the fully processive dual-site network-which is known to be globally stable and therefore lack both oscillations and bistability. Our proofs are enabled by a Hopf bifurcation criterion due to Yang, analyses of Newton polytopes arising from Hurwitz determinants, and recent characterizations of multistationarity for networks having a steady-state parametrization.

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.

Raman micro-spectroscopy monitors acquired resistance to targeted cancer therapy at the cellular level

Monitoring the drug efficacy or resistance in vitro is usually carried out by measuring the response of single few proteins. However, observation of single proteins instead of an integral cell response may lead to results that are not consistent with patient’s response to a drug. We present a Raman spectroscopic method that detects the integral cell response to drugs such as tyrosine kinase inhibitors (TKIs). Non-small cell lung cancer (NSCLC) patients with EGFR mutations develop acquired resistance to first (erlotinib)- and third (osimertinib)-generation TKIs. Large erlotinib-induced differences were detected by Raman micro-spectroscopy in NSCLC cells without T790M EGFR mutation but not in cells with this mutation. Additionally, Raman difference spectra detected the response of NSCLC cells with T790M EGFR mutation to second- (neratinib) and third-generation (osimertinib) TKIs, and the resistance of cells with T790M/C797S EGFR mutation to osimertinib. Thus, the in vitro Raman results indicated that NSCLC cells with T790M and T790M/C797S EGFR mutations are resistant to erlotinib- and osimertinib, respectively, consistent with the observed responses of patients. This study shows the potential of Raman micro-spectroscopy to monitor drug resistance and opens a new door to in vitro companion diagnostics for screening personalized therapies.

When More Is Less: Dual Phosphorylation Protects Signaling Off State against Overexpression

Kinases in signaling pathways are commonly activated by multisite phosphorylation. For example, the mitogen-activated protein kinase Erk is activated by its kinase Mek by two consecutive phosphorylations within its activation loop. In this article, we use kinetic models to study how the activation of Erk is coupled to its abundance. Intuitively, Erk activity should rise with increasing amounts of Erk protein. However, a mathematical model shows that the signaling off state is robust to increasing amounts of Erk, and Erk activity may even decline with increasing amounts of Erk. This counterintuitive, bell-shaped response of Erk activity to increasing amounts of Erk arises from the competition of the unmodified and single phosphorylated form of Erk for access to its kinase Mek. This shows that phosphorylation cycles can contain an intrinsic robustness mechanism that protects signaling from aberrant activation e.g., by gene expression noise or kinase overexpression after gene duplication events in diseases like cancer.

A quantitative model of developmental RTK signaling

Receptor tyrosine kinases (RTKs) control a wide range of developmental processes, from the first stages of embryogenesis to postnatal growth and neurocognitive development in the adult. A significant share of our knowledge about RTKs comes from genetic screens in model organisms, which provided numerous examples demonstrating how specific cell fates and morphologies are abolished when RTK activation is either abrogated or significantly reduced. Aberrant activation of such pathways has also been recognized in many forms of cancer. More recently, studies of human developmental syndromes established that excessive activation of RTKs and their downstream signaling effectors, most notably the Ras signaling pathway, can also lead to structural and functional defects. Given that both insufficient and excessive pathway activation can lead to abnormalities, mechanistic analysis of developmental RTK signaling must address quantitative questions about its regulation and function. Patterning events controlled by the RTK Torso in the early Drosophila embryo are well-suited for this purpose. This mini review summarizes current state of knowledge about Torso-dependent Ras activation and discusses its potential to serve as a quantitative model for studying the general principles of Ras signaling in development and disease.

Numerical algebraic geometry for model selection and its application to the life sciences

Researchers working with mathematical models are often confronted by the related problems of parameter estimation, model validation and model selection. These are all optimization problems, well known to be challenging due to nonlinearity, non-convexity and multiple local optima. Furthermore, the challenges are compounded when only partial data are available. Here, we consider polynomial models (e.g. mass-action chemical reaction networks at steady state) and describe a framework for their analysis based on optimization using numerical algebraic geometry. Specifically, we use probability-one polynomial homotopy continuation methods to compute all critical points of the objective function, then filter to recover the global optima. Our approach exploits the geometrical structures relating models and data, and we demonstrate its utility on examples from cell signalling, synthetic biology and epidemiology.

[The role of extracellular signal regulated kinase 1/2 in mediating osteodifferentiation of human periodontal ligament cells induced by cyclic stretch]

OBJECTIVE

This study aimed to investigate the mechanism of cyclic stretch that promotesthe osteogenic differentiation of human periodontal ligament cells (hPDLCs) through the mediation of extracellular-signal-regulated kinase 1/2 (ERK1/2).

METHODS

hPDLCs were isolated through the explant method and cultured in vitro. hPDLCs were mechanically stimulated by a multi-channel cell-stress-loading system for 1, 3, 6, 12, and 24 h. The magnitude of stretch was 10% deformation, and the frequency was 0.5 Hz. Nonloaded cells were used as control group. ERK1/2 activation was blocked by U0126, a specific ERK1/2 pathway inhibitor. Additionally, hPDLCs were transfected with adenoviral vector encoding dominant negative ERK1/2 (DN-ERK1/2) to continuouslyinhibit ERK1/2 activation. The mRNA and protein levels of target geneswere detected through real-time polymerase chain reaction and Western blot.

RESULTS

Cyclic stretching promoted the expression of ERK1/2, osteocalcin (OCN) mRNA, and bone sialoprotein (BSP) mRNA. The expression of runt-related transcription factor (Runx) 2 protein and mRNA also increased at 3 and 6 h of cyclic stretching. The inhibition of ERK1/2 by U0126 and DN-ERK1/2 suppressed the expressionof Runx2 mRNA, OCN mRNA, BSP mRNA, Runx 2 protein, and p-ERK1/2 protein relative to that in stretched cells without the ERK1/2 inhibitor.

CONCLUSIONS

ERK1/2 is a critical molecule in the mediation ofthe osteogenic differentiation of hPDLCs under mechanical stimulation. ERK1/2 activation induced the elevation of Runx2 protein levels, which may be involved in the stretch-induced expressions of OCN and BSP.

Two faces of competition: target-mediated reverse signalling in microRNA and mitogen-activated protein kinase regulatory networks

Biomolecular regulatory networks are organised around hubs, which can interact with a large number of targets. These targets compete with each other for access to their common hubs, but whether the effect of this competition would be significant in magnitude and in function is not clear. In this review, the authors discuss recent in vivo studies that analysed the system level retroactive effects induced by target competition in microRNA and mitogen-activated protein kinase regulatory networks. The results of these studies suggest that downstream targets can regulate the overall state of their upstream regulators, and thus cannot be ignored in analysing biomolecular networks.

Parallel imaging of Drosophila embryos for quantitative analysis of genetic perturbations of the Ras pathway

The Ras pathway patterns the poles of the Drosophila embryo by downregulating the levels and activity of a DNA-binding transcriptional repressor Capicua (Cic). We demonstrate that the spatiotemporal pattern of Cic during this signaling event can be harnessed for functional studies of mutations in the Ras pathway in human diseases. Our approach relies on a new microfluidic device that enables parallel imaging of Cic dynamics in dozens of live embryos. We found that although the pattern of Cic in early embryos is complex, it can be accurately approximated by a product of one spatial profile and one time-dependent amplitude. Analysis of these functions of space and time alone reveals the differential effects of mutations within the Ras pathway. Given the highly conserved nature of Ras-dependent control of Cic, our approach provides new opportunities for functional analysis of multiple sequence variants from developmental abnormalities and cancers.

Mitogen-activated protein kinases participate in determination of apical-basal symmetry in Pisum sativum

Mitogen-activated protein kinases (MAPKs) are implicated in various processes in plants. Apart from response to biotic and abiotic stresses they are involved in regulation of embryo development. Although MAPKs were found to be indispensable during embryo development it has never been established at which stages of embryogenesis and in which regions of a plant embryo activated MAPKs can be observed. Here, we show that apical and basal regions display activation of the same types of MAPKs and the only difference concerns the level of their phosphorylation and cellular localization. Dually-phosphorylated MAPKs were found in nuclei of the apical region of an embryo both at the early and late cotyledonary stage while no immunofluorescence signals were detected in nuclei of the basal region. However, in this case phosphorylated MAPKs were immunodetected in cytoplasm in the apical domain of cortex cells, indicating their role in auxin transport from the basal to apical region of an embryo. Furthermore, obtained data indicate that nuclear localization of activated MAPKs may result from epigenetic modifications and polar auxin transport. The presented data and previous studies lead to the conclusion that activated MAPKs and their cellular localization define apical and basal regions during formation of an apical-basal axis.

In vivo severity ranking of Ras pathway mutations associated with developmental disorders

Germ-line mutations in components of the Ras/MAPK pathway result in developmental disorders called RASopathies, affecting about 1/1,000 human births. Rapid advances in genome sequencing make it possible to identify multiple disease-related mutations, but there is currently no systematic framework for translating this information into patient-specific predictions of disease progression. As a first step toward addressing this issue, we developed a quantitative, inexpensive, and rapid framework that relies on the early zebrafish embryo to assess mutational effects on a common scale. Using this assay, we assessed 16 mutations reported in MEK1, a MAPK kinase, and provide a robust ranking of these mutations. We find that mutations found in cancer are more severe than those found in both RASopathies and cancer, which, in turn, are generally more severe than those found only in RASopathies. Moreover, this rank is conserved in other zebrafish embryonic assays and Drosophila-specific embryonic and adult assays, suggesting that our ranking reflects the intrinsic property of the mutant molecule. Furthermore, this rank is predictive of the drug dose needed to correct the defects. This assay can be readily used to test the strengths of existing and newly found mutations in MEK1 and other pathway components, providing the first step in the development of rational guidelines for patient-specific diagnostics and treatment of RASopathies.

Single-Step Affinity Purification of ERK Signaling Complexes Using the Streptavidin-Binding Peptide (SBP) Tag

Elucidation of biological functions of signaling proteins is facilitated by studying their protein-protein interaction networks. Affinity purification combined with mass spectrometry (AP-MS) has become a favorite method to study protein complexes. Here we describe a procedure for single-step purification of ERK (Rolled) and associated proteins from Drosophila cultured cells. The use of the streptavidin-binding peptide (SBP) tag allows for a highly efficient isolation of native ERK signaling complexes, which are suitable for subsequent analysis by mass spectrometry. Our analysis of the ERK interactome has identified both known and novel signaling components. This method can be easily adapted for SBP-based purification of protein complexes in any expression system.

Co-culture Activation of MAP Kinase in Drosophila S2 Cells

Intercellular communication often involves phosphorylation of signal transduction proteins, including mitogen-activated protein kinases (MAPKs). Immunological detection of phosphorylated MAPK can be used to monitor signaling in vivo, identify novel pathway components, and assess ligand activity. In this chapter, I describe a cell co-culture method to assess activity of cell-bound extracellular ligands that result in phosphorylation of the ERK (extracellular signal-regulated kinase) MAPK in Drosophila. This protocol may be adaptable to other pathways and/or model systems.

Reconstructing ERK Signaling in the Drosophila Embryo from Fixed Images

The early Drosophila embryo provides unique opportunities for quantitative studies of ERK signaling. This system is characterized by simple anatomy, the ease of obtaining large numbers of staged embryos, and the availability of powerful tools for genetic manipulation of the ERK pathway. Here, we describe how these experimental advantages can be combined with recently developed microfluidic devices for high throughput imaging of ERK activation dynamics. We focus on the stage during the third hour of development, when ERK activation is essential for patterning of the future nerve cord. Our approach starts with an ensemble of fixed embryos stained with an antibody that recognizes the active, dually phosphorylated form of ERK. Each embryo in this ensemble provides a snapshot of the spatial and temporal pattern of ERK activation during development. We then quantitatively estimate the ages of fixed embryos using a model that links their morphology and developmental time. This model is learned based on live imaging of cellularization and gastrulation, two highly stereotyped morphogenetic processes at this stage of embryogenesis. Applying this approach, we can characterize ERK signaling at high spatial and temporal resolution. Our methodology can be readily extended to studies of ERK regulation and function in multiple mutant backgrounds, providing a versatile assay for quantitative studies of developmental ERK signaling.

Theoretical Investigation of the Mechanisms of ERK2 Enzymatic Catalysis

ERK2 are protein kinases that during the enzymatic catalysis, in contrast to traditional enzymes, utilize additional interactions with substrates outside of the active sites. It is widely believed that these docking interactions outside of the enzymatic pockets enhance the specificity of these proteins. However, the molecular mechanisms of how the docking interactions affect the catalysis remain not well understood. Here, we develop a simple theoretical approach to analyze the enzymatic catalysis in ERK2 proteins. Our method is based on first-passage process analysis, and it provides explicit expressions for all dynamic properties of the system. It is found that there are specific binding energies for substrates in docking and catalytic domains that lead to maximal enzymatic reaction rates. Thus, we propose that the role of the docking interactions is not only to increase the enzymatic specificity but also to optimize the dynamics of the catalytic process. Our theoretical results are utilized to describe experimental observations on ERK2 enzymatic activities.

A Transport Model for Estimating the Time Course of ERK Activation in the C. elegans Germline

The Caenorhabditis elegans germline is a well-studied model system for investigating the control of cell fate by signaling pathways. Cell signals at the distal tip of the germline promote cell proliferation; just before the loop, signals couple cell maturation to organism-level nutrient status; at the proximal end of the germline, signals coordinate oocyte maturation and fertilization in the presence of sperm. The latter two events require dual phosphorylation and activation of ERK, the effector molecule of the Ras/MAPK cascade. In C. elegans, ERK is known as MPK-1. At this point, none of today's methods for real-time monitoring of dually phosphorylated MPK-1 are working in the germline. Consequently, quantitative understanding of the MPK-1-dependent processes during germline development is limited. Here, we make a step toward advancing this understanding using a model-based framework that reconstructs the time course of MPK-1 activation from a snapshot of a fixed germline. Our approach builds on a number of recent studies for estimating temporal dynamics from fixed organisms, but takes advantage of the anatomy of the germline to simplify the analysis. Our model predicts that the MPK-1 signal turns on ∼30 h into germ cell progression and peaks ∼7 h later.

Sulforaphane-cysteine suppresses invasion via downregulation of galectin-1 in human prostate cancer DU145 and PC3 cells

Our previous study showed that sulforaphane (SFN) inhibits invasion in human prostate cancer DU145 cells; however, the underlying mechanisms were not profoundly investigated. In the present study, we found that sulforaphane-cysteine (SFN-Cys), as a metabolite of SFN, inhibits invasion and possesses a novel mechanism in prostate cancer DU145 and PC3 cells. The scratch and Transwell assays showed that SFN-Cys (15 µM) inhibited both migration and invasion, with cell morphological changes, such as cell shrinkage and pseudopodia shortening. The cell proliferation (MTS) assay indicated that cell viability was markedly suppressed with increasing concentrations of SFN‑Cys. Furthermore, the Transwell assay showed that inhibition of SFN‑Cys‑triggered invasion was tightly linked to the sustained extracellular signal-regulated kinase 1/2 (ERK1/2) phosphorylation. Western blot analysis revealed that SFN-Cys downregulated galectin-1 protein, an invasion‑related protein, and that the galectin‑1 reduction could be blocked by ERK1/2 inhibitor PD98059 (25 µM). Moreover, immunofluorescence staining showed that the expression level of galectin-1 protein was significantly reduced in the cells treated with SFN‑Cys. Hence, SFN‑Cys‑inhibited invasion resulted from the sustained ERK1/2 phosphorylation and ERK1/2‑triggered galectin-1 downregulation, suggesting that galectin-1 is a new SFN-Cys target inhibiting invasion apart from ERK1/2, in the treatment of prostate cancer.

Long-term dynamics of multisite phosphorylation

A systematic framework for exploring the long-term dynamics of a reaction network is applied to a minimal model of ERK regulation that distinguishes both monophosphorylated forms and allows for nonzero enzyme processivity. Bistability and oscillations can be observed at high levels of processivity.

Multisite phosphorylation cycles are ubiquitous in cell regulation systems and are studied at multiple levels of complexity, from molecules to organisms, with the ultimate goal of establishing predictive understanding of the effects of genetic and pharmacological perturbations of protein phosphorylation in vivo. Achieving this goal is essentially impossible without mathematical models, which provide a systematic framework for exploring dynamic interactions of multiple network components. Most of the models studied to date do not discriminate between the distinct partially phosphorylated forms and focus on two limiting reaction regimes, distributive and processive, which differ in the number of enzyme–substrate binding events needed for complete phosphorylation or dephosphorylation. Here we use a minimal model of extracellular signal-related kinase regulation to explore the dynamics of a reaction network that includes all essential phosphorylation forms and arbitrary levels of reaction processivity. In addition to bistability, which has been studied extensively in distributive mechanisms, this network can generate periodic oscillations. Both bistability and oscillations can be realized at high levels of reaction processivity. Our work provides a general framework for systematic analysis of dynamics in multisite phosphorylation systems.

RASopathies: unraveling mechanisms with animal models

RASopathies are developmental disorders caused by germline mutations in the Ras-MAPK pathway, and are characterized by a broad spectrum of functional and morphological abnormalities. The high incidence of these disorders (∼1/1000 births) motivates the development of systematic approaches for their efficient diagnosis and potential treatment. Recent advances in genome sequencing have greatly facilitated the genotyping and discovery of mutations in affected individuals, but establishing the causal relationships between molecules and disease phenotypes is non-trivial and presents both technical and conceptual challenges. Here, we discuss how these challenges could be addressed using genetically modified model organisms that have been instrumental in delineating the Ras-MAPK pathway and its roles during development. Focusing on studies in mice, zebrafish and Drosophila, we provide an up-to-date review of animal models of RASopathies at the molecular and functional level. We also discuss how increasingly sophisticated techniques of genetic engineering can be used to rigorously connect changes in specific components of the Ras-MAPK pathway with observed functional and morphological phenotypes. Establishing these connections is essential for advancing our understanding of RASopathies and for devising rational strategies for their management and treatment.

Summary: Developmental disorders caused by germline mutations in the Ras-MAPK pathway are called RASopathies. Studies with animal models, including mice, zebrafish and Drosophila, continue to enhance our understanding of these diseases.

Kinetic Adaptations of Myosins for Their Diverse Cellular Functions

Members of the myosin superfamily are involved in all aspects of eukaryotic life. Their function ranges from the transport of organelles and cargos to the generation of membrane tension, and the contraction of muscle. The diversity of physiological functions is remarkable, given that all enzymatically active myosins follow a conserved mechanoenzymatic cycle in which the hydrolysis of ATP to ADP and inorganic phosphate is coupled to either actin-based transport or tethering of actin to defined cellular compartments. Kinetic capacities and limitations of a myosin are determined by the extent to which actin can accelerate the hydrolysis of ATP and the release of the hydrolysis products and are indispensably linked to its physiological tasks. This review focuses on kinetic competencies that - together with structural adaptations - result in myosins with unique mechanoenzymatic properties targeted to their diverse cellular functions.

Impact of ERK activation on fly survival and Wolbachia-mediated protection during virus infection

Elevated levels of reactive oxygen species (ROS) provide protection against virus-induced mortality in Drosophila. In addition to contributing to oxidative stress, ROS are known to activate a number of signalling pathways including the extracellular signal-regulated kinases (ERK) signalling cascade. It was recently shown that ERK signalling is important for resistance against viral replication and invasion in cultured Drosophila cells and the gut epithelium of adult flies. Here, using a Drosophila loss-of-function ERK (rolled) mutant we demonstrated that ERK is important for fly survival during virus infection. ERK mutant flies subjected to Drosophila C virus (DCV) oral and systemic infection were more susceptible to virus-induced mortality as compared with wild-type flies. We have demonstrated experimentally that ERK activation is important for fly survival during oral and systemic virus infection. Given that elevated ROS correlates with Wolbachia-mediated antiviral protection, we also investigated the involvement of ERK in antiviral protection in flies infected by Wolbachia. The results indicate that ERK activation is increased in the presence of Wolbachia but this does not appear to influence Wolbachia-mediated antiviral protection, at least during systemic infection.

Two hydrophobic residues can determine the specificity of mitogen-activated protein kinase docking interactions

MAPKs bind to many of their upstream regulators and downstream substrates via a short docking motif (the D-site) on their binding partner. MAPKs that are in different families (e.g. ERK, JNK, and p38) can bind selectively to D-sites in their authentic substrates and regulators while discriminating against D-sites in other pathways. Here we demonstrate that the short hydrophobic region at the distal end of the D-site plays a critical role in determining the high selectivity of JNK MAPKs for docking sites in their cognate MAPK kinases. Changing just 1 or 2 key hydrophobic residues in this submotif is sufficient to turn a weak JNK-binding D-site into a strong one, or vice versa. These specificity-determining differences are also found in the D-sites of the ETS family transcription factors Elk-1 and Net. Moreover, swapping two hydrophobic residues between these D-sites switches the relative efficiency of Elk-1 and Net as substrates for ERK versus JNK, as predicted. These results provide new insights into docking specificity and suggest that this specificity can evolve rapidly by changes to just 1 or 2 amino acids.

Selected Reaction Monitoring Mass Spectrometry for Absolute Protein Quantification

Absolute quantification of target proteins within complex biological samples is critical to a wide range of research and clinical applications. This protocol provides step-by-step instructions for the development and application of quantitative assays using selected reaction monitoring (SRM) mass spectrometry (MS). First, likely quantotypic target peptides are identified based on numerous criteria. This includes identifying proteotypic peptides, avoiding sites of posttranslational modification, and analyzing the uniqueness of the target peptide to the target protein. Next, crude external peptide standards are synthesized and used to develop SRM assays, and the resulting assays are used to perform qualitative analyses of the biological samples. Finally, purified, quantified, heavy isotope labeled internal peptide standards are prepared and used to perform isotope dilution series SRM assays. Analysis of all of the resulting MS data is presented. This protocol was used to accurately assay the absolute abundance of proteins of the chemotaxis signaling pathway within RAW 264.7 cells (a mouse monocyte/macrophage cell line). The quantification of Gi2 (a heterotrimeric G-protein α-subunit) is described in detail.

MAP kinase modules: the excursion model and the steps that count

MAP kinase modules propagate diverse extracellular signals to downstream effectors. The two dual phosphorylation reactions catalyzed by the modules are thought to control the switch behavior of the pathway. Here we review recent approaches to understand these pathways through signal-to-response studies in cells and in vitro. These data are reconciled with physical models as well as predictions made on mathematical and theoretical grounds. Biochemical analysis has shown recently that the dual phosphorylation reactions catalyzed by MAP kinase modules are sequential at both levels of the cascade. The observed order of phosphorylation events suggests an excursion from the Ser/Thr kinase activity of the MAP3K into Tyr kinase activity of the central dual specificity MAP2K. How the order of events might be encoded in the structures and interactions is discussed. The ordered mechanism confirms predictions that reactions should be sequential to generate the steep signal-to-response curves and delayed responses observed in cells.

Mapping the binding interface of ERK and transcriptional repressor Capicua using photocrosslinking

Extracellular signal-regulated kinase (ERK) coordinates cellular responses to a range of stimuli by phosphorylating its numerous substrates. One of these substrates, Capicua (Cic), is a transcriptional repressor that was first identified in Drosophila and has been implicated in a number of human diseases. Here we use a chemical biology approach to map the binding interface of ERK and Cic. The noncanonical amino acid p-azidophenylalanine (AzF) was introduced into the ERK-binding region of Drosophila Cic, and photocrosslinking and tandem mass spectrometry were used to pinpoint its binding site on ERK. We also identified the ERK-binding region of human Cic and showed that it binds to the same site on ERK despite lacking conservation with the Drosophila Cic binding region. Finally, we mapped the amino acids involved in human Cic binding to ERK using AzF-labeled ERK. These results reveal the molecular details of the ERK-Cic interaction and demonstrate that the photocrosslinking approach is complementary to existing methods for mapping kinase-substrate binding interfaces.

Label-Free Raman Spectroscopic Imaging Monitors the Integral Physiologically Relevant Drug Responses in Cancer Cells

Predictions about the cellular efficacy of drugs tested in vitro are usually based on the measured responses of a few proteins or signal transduction pathways. However, cellular proteins are highly coupled in networks, and observations of single proteins may not adequately reflect the in vivo cellular response to drugs. This might explain some large discrepancies between in vitro drug studies and drug responses observed in patients. We present a novel in vitro marker-free approach that enables detection of cellular responses to a drug. We use Raman spectral imaging to measure the effect of the epidermal growth factor receptor (EGFR) inhibitor panitumumab on cell lines expressing wild-type Kirsten-Ras (K-Ras) and oncogenic K-Ras mutations. Oncogenic K-Ras mutation blocks the response to anti-EGFR therapy in patients, but this effect is not readily observed in vitro. The Raman studies detect large panitumumab-induced differences in vitro in cells harboring wild-type K-Ras as seen in A in red but not in cells with K-Ras mutations as seen in B; these studies reflect the observed patient outcomes. However, the effect is not observed when extracellular-signal-regulated kinase phosphorylation is monitored. The Raman spectra show for cells with wild-type K-Ras alterations based on the responses to panitumumab. The subcellular component with the largest spectral response to panitumumab was lipid droplets, but this effect was not observed when cells harbored K-Ras mutations. This study develops a noninvasive, label-free, in vitro vibrational spectroscopic test to determine the integral physiologically relevant drug response in cell lines. This approach opens a new field of patient-centered drug testing that could deliver superior patient therapies.

Dynamics of Inductive ERK Signaling in the Drosophila Embryo

Transient activation of the highly conserved extracellular-signal-regulated kinase (ERK) establishes precise patterns of cell fates in developing tissues. Quantitative parameters of these transients are essentially unknown, but a growing number of studies suggest that changes in these parameters can lead to a broad spectrum of developmental abnormalities. We provide a detailed quantitative picture of an ERK-dependent inductive signaling event in the early Drosophila embryo, an experimental system that offers unique opportunities for high-throughput studies of developmental signaling. Our analysis reveals a spatiotemporal pulse of ERK activation that is consistent with a model in which transient production of a short-ranged ligand feeds into a simple signal interpretation system. The pulse of ERK signaling acts as a switch in controlling the expression of the ERK target gene. The quantitative approach that led to this model, based on the integration of data from fixed embryos and live imaging, can be extended to other developmental systems patterned by transient inductive signals.

Paradoxical results in perturbation-based signaling network reconstruction

Mathematical models are extensively employed to understand physicochemical processes in biological systems. In the absence of detailed mechanistic knowledge, models are often based on network inference methods, which in turn rely upon perturbations to nodes by biochemical means. We have discovered a potential pitfall of the approach underpinning such methods when applied to signaling networks. We first show experimentally, and then explain mathematically, how even in the simplest signaling systems, perturbation methods may lead to paradoxical conclusions: for any given pair of two components X and Y, and depending upon the specific intervention on Y, either an activation or a repression of X could be inferred. This effect is of a different nature from incomplete network identification due to underdetermined data and is a phenomenon intrinsic to perturbations. Our experiments are performed in an in vitro minimal system, thus isolating the effect and showing that it cannot be explained by feedbacks due to unknown intermediates. Moreover, our in vitro system utilizes proteins from a pathway in mammalian (and other eukaryotic) cells that play a central role in proliferation, gene expression, differentiation, mitosis, cell survival, and apoptosis. This pathway is the perturbation target of contemporary therapies for various types of cancers. The results presented here show that the simplistic view of intracellular signaling networks being made up of activation and repression links is seriously misleading, and call for a fundamental rethinking of signaling network analysis and inference methods.