Dissection of MKK6 and p38 Signaling Using Light-Activated Protein Kinases.

Stress-activated signaling pathways orchestrate cellular behaviors and fates. Studying the precise role(s) of stress-activated protein kinases is challenging, because stress conditions induce adaptation and impose selection pressure. To meet this challenge, we have applied an optogenetic system with a single plasmid to express light-activated p38α or its upstream activator, MKK6, in conjunction with live-cell fluorescence microscopy. In starved cells, decaging of constitutively active p38α or MKK6 by brief exposure to UV light elicits rapid p38-mediated signaling, release of cytochrome c from mitochondria, and apoptosis with different kinetics. In parallel, light activation of p38α also suppresses autophagosome formation, similarly to stimulation with growth factors that activate PI3K/Akt/mTORC1 signaling. Active MKK6 negatively regulates serum-induced ERK activity, which is p38-independent as previously reported. Here, we reproduce that result with the one plasmid system and show that although decaging active p38α does not reduce basal ERK activity in our cells, it can block growth factor-stimulated ERK signaling in serum-starved cells. These results clarify the roles of MKK6 and p38α in dynamic signaling programs, which act in concert to actuate apoptotic death while suppressing cell survival mechanisms.

G-actin diffusion is insufficient to achieve F-actin assembly in fast-treadmilling protrusions.

To generate forces that drive migration of a eukaryotic cell, arrays of actin filaments (F-actin) are assembled at the cell's leading membrane edge. To maintain cell propulsion and respond to dynamic external cues, actin filaments must be disassembled to regenerate the actin monomers (G-actin), and transport of G-actin from sites of disassembly back to the leading edge completes the treadmilling cycle and limits the flux of F-actin assembly. Whether or not molecular diffusion is sufficient for G-actin transport has been a long-standing topic of debate, in part because the dynamic nature of cell motility and migration hinders the estimation of transport parameters. In this work, we applied an experimental system in which cells adopt an approximately constant and symmetrical shape; they cannot migrate but exhibit fast, steady treadmilling in the thin region protruding from the cell. Using fluorescence recovery after photobleaching, we quantified the relative concentrations and corresponding fluxes of F- and G-actin in this system. In conjunction with mathematical modeling, constrained by measured features of each region of interest, this approach revealed that diffusion alone cannot account for the transport of G-actin to the leading edge. Although G-actin diffusion and vectorial transport might vary with position in the protruding region, good agreement with the fluorescence recovery after photobleaching measurements was achieved by a model with constant G-actin diffusivity ∼2 µm2/s and anterograde G-actin velocity less than 1 µm/s.

A kinetic model of phospholipase C-γ1 linking structure-based insights to dynamics of enzyme autoinhibition and activation.

Phospholipase C-γ1 (PLC-γ1) is a receptor-proximal enzyme that promotes signal transduction through PKC in mammalian cells. Because of the complexity of PLC-γ1 regulation, a two-state (inactive/active) model does not account for the intricacy of activation and inactivation steps at the plasma membrane. Here, we introduce a structure-based kinetic model of PLC-γ1, considering interactions of its regulatory Src homology 2 (SH2) domains and perturbation of those dynamics upon phosphorylation of Tyr783, a hallmark of activation. For PLC-γ1 phosphorylation to dramatically enhance enzyme activation as observed, we found that high intramolecular affinity of the C-terminal SH2 (cSH2) domain-pTyr783 interaction is critical, but this affinity need not outcompete the autoinhibitory interaction of the cSH2 domain. Under conditions for which steady-state PLC-γ1 activity is sensitive to the rate of Tyr783 phosphorylation, maintenance of the active state is surprisingly insensitive to the phosphorylation rate, since pTyr783 is well protected by the cSH2 domain while the enzyme is active. In contrast, maintenance of enzyme activity is sensitive to the rate of PLC-γ1 membrane (re)binding. Accordingly, we found that hypothetical PLC-γ1 mutations that either weaken autoinhibition or strengthen membrane binding influence the activation kinetics differently, which could inform the characterization of oncogenic variants. Finally, we used this newly informed kinetic scheme to refine a spatial model of PLC/PKC polarization during chemotaxis. The refined model showed improved stability of the polarized pattern while corroborating previous qualitative predictions. As demonstrated here for PLC-γ1, this approach may be adapted to model the dynamics of other receptor- and membrane-proximal enzymes.

Semi-autonomous wound invasion via matrix-deposited, haptotactic cues.

Proper wound healing relies on invasion of fibroblasts via directed migration. While the related experimental and mathematical modeling literature has mainly focused on cell migration directed by soluble cues (chemotaxis), there is ample evidence that fibroblast migration is also directed by insoluble, matrix-bound cues (haptotaxis). Furthermore, numerous studies indicate that fibronectin (FN), a haptotactic ligand for fibroblasts, is present and dynamic in the provisional matrix throughout the proliferative phase of wound healing. In the present work, we show the plausibility of a hypothesis that fibroblasts themselves form and maintain haptotactic gradients in a semi-autonomous fashion. As a precursor to this, we examine the positive control scenario where FN is pre-deposited in the wound matrix, and fibroblasts maintain haptotaxis by removing FN at an appropriate rate. After developing conceptual and quantitative understanding of this scenario, we consider two cases in which fibroblasts activate the latent form of a matrix-loaded cytokine, TGFß, which upregulates the fibroblasts' own secretion of FN. In the first of these, the latent cytokine is pre-patterned and released by the fibroblasts. In the second, fibroblasts in the wound produce the latent TGFß, with the presence of the wound providing the only instruction. In all cases, wound invasion is more effective than a negative control model with haptotaxis disabled; however, there is a trade-off between the degree of fibroblast autonomy and the rate of invasion.

Modeling cell protrusion predicts how myosin II and actin turnover affect adhesion-based signaling.

Orchestration of cell migration is essential for development, tissue regeneration, and the immune response. This dynamic process integrates adhesion, signaling, and cytoskeletal subprocesses across spatial and temporal scales. In mesenchymal cells, adhesion complexes bound to extracellular matrix mediate both biochemical signal transduction and physical interaction with the F-actin cytoskeleton. Here, we present a mathematical model that offers insight into both aspects, considering spatiotemporal dynamics of nascent adhesions, active signaling molecules, mechanical clutching, actin treadmilling, and nonmuscle myosin II contractility. At the core of the model is a positive feedback loop, whereby adhesion-based signaling promotes generation of barbed ends at, and protrusion of, the cell's leading edge, which in turn promotes formation and stabilization of nascent adhesions. The model predicts a switch-like transition and optimality of membrane protrusion, determined by the balance of actin polymerization and retrograde flow, with respect to extracellular matrix density. The model, together with new experimental measurements, explains how protrusion can be modulated by mechanical effects (nonmuscle myosin II contractility and adhesive bond stiffness) and F-actin turnover.

Optical control of MAP kinase kinase 6 (MKK6) reveals that it has divergent roles in pro-apoptotic and anti-proliferative signaling.

The selective pressure imposed by extrinsic death signals and stressors adds to the challenge of isolating and interpreting the roles of proteins in stress-activated signaling networks. By expressing a kinase with activating mutations and a caged lysine blocking the active site, we can rapidly switch on catalytic activity with light and monitor the ensuing dynamics. Applying this approach to MAP kinase 6 (MKK6), which activates the p38 subfamily of MAPKs, we found that decaging active MKK6 in fibroblasts is sufficient to trigger apoptosis in a p38-dependent manner. Both in fibroblasts and in a murine melanoma cell line expressing mutant B-Raf, MKK6 activation rapidly and potently inhibited the pro-proliferative extracellular signal-regulated kinase (ERK) pathway; to our surprise, this negative cross-regulation was equally robust when all p38 isoforms were inhibited. These results position MKK6 as a new pleiotropic signal transducer that promotes both pro-apoptotic and anti-proliferative signaling, and they highlight the utility of caged, light-activated kinases for dissecting stress-activated signaling networks.

Mechanistic models of PLC/PKC signaling implicate phosphatidic acid as a key amplifier of chemotactic gradient sensing.

Chemotaxis of fibroblasts and other mesenchymal cells is critical for embryonic development and wound healing. Fibroblast chemotaxis directed by a gradient of platelet-derived growth factor (PDGF) requires signaling through the phospholipase C (PLC)/protein kinase C (PKC) pathway. Diacylglycerol (DAG), the lipid product of PLC that activates conventional PKCs, is focally enriched at the up-gradient leading edge of fibroblasts responding to a shallow gradient of PDGF, signifying polarization. To explain the underlying mechanisms, we formulated reaction-diffusion models including as many as three putative feedback loops based on known biochemistry. These include the previously analyzed mechanism of substrate-buffering by myristoylated alanine-rich C kinase substrate (MARCKS) and two newly considered feedback loops involving the lipid, phosphatidic acid (PA). DAG kinases and phospholipase D, the enzymes that produce PA, are identified as key regulators in the models. Paradoxically, increasing DAG kinase activity can enhance the robustness of DAG/active PKC polarization with respect to chemoattractant concentration while decreasing their whole-cell levels. Finally, in simulations of wound invasion, efficient collective migration is achieved with thresholds for chemotaxis matching those of polarization in the reaction-diffusion models. This multi-scale modeling framework offers testable predictions to guide further study of signal transduction and cell behavior that affect mesenchymal chemotaxis.

Kinetic Modeling and Analysis of the Akt/Mechanistic Target of Rapamycin Complex 1 (mTORC1) Signaling Axis Reveals Cooperative, Feedforward Regulation.

Mechanistic target of rapamycin complex 1 (mTORC1) controls biosynthesis and has been implicated in uncontrolled cell growth in cancer. Although many details of mTORC1 regulation are well understood, a systems-level, predictive framework synthesizing those details is currently lacking. We constructed various mathematical models of mTORC1 activation mediated by Akt and aligned the model outputs to kinetic data acquired for growth factor-stimulated cells. A model based on a putative feedforward loop orchestrated by Akt consistently predicted how the pathway was altered by depletion of key regulatory proteins. Analysis of the successful model also elucidates two dynamical motifs: neutralization of a negative regulator, which characterizes how Akt indirectly activates mTORC1, and seesaw enzyme regulation, which describes how activated and inhibited states of mTORC1 are controlled in concert to produce a nonlinear, ultrasensitive response. Such insights lend quantitative understanding of signaling networks and their precise manipulation in various contexts.

Lamellipodia are crucial for haptotactic sensing and response.

Haptotaxis is the process by which cells respond to gradients of substrate-bound cues, such as extracellular matrix proteins (ECM); however, the cellular mechanism of this response remains poorly understood and has mainly been studied by comparing cell behavior on uniform ECMs with different concentrations of components. To study haptotaxis in response to gradients, we utilized microfluidic chambers to generate gradients of the ECM protein fibronectin, and imaged the cell migration response. Lamellipodia are fan-shaped protrusions that are common in migrating cells. Here, we define a new function for lamellipodia and the cellular mechanism required for haptotaxis - differential actin and lamellipodial protrusion dynamics lead to biased cell migration. Modest differences in lamellipodial dynamics occurring over time periods of seconds to minutes are summed over hours to produce differential whole cell movement towards higher concentrations of fibronectin. We identify a specific subset of lamellipodia regulators as being crucial for haptotaxis. Numerous studies have linked components of this pathway to cancer metastasis and, consistent with this, we find that expression of the oncogenic Rac1 P29S mutation abrogates haptotaxis. Finally, we show that haptotaxis also operates through this pathway in 3D environments.

A Reaction-Diffusion Model Explains Amplification of the PLC/PKC Pathway in Fibroblast Chemotaxis.

During the proliferative phase of cutaneous wound healing, dermal fibroblasts are recruited into the clotted wound by a concentration gradient of platelet-derived growth factor (PDGF), together with other spatial cues. Despite the importance of this chemotactic process, the mechanisms controlling the directed migration of slow-moving mesenchymal cells such as fibroblasts are not well understood. Here, we develop and analyze a reaction-diffusion model of phospholipase C/protein kinase C (PKC) signaling, which was recently identified as a requisite PDGF-gradient-sensing pathway, with the goal of identifying mechanisms that can amplify its sensitivity in the shallow external gradients typical of chemotaxis experiments. We show that phosphorylation of myristoylated alanine-rich C kinase substrate by membrane-localized PKC constitutes a positive feedback that is sufficient for local pathway amplification. The release of phosphorylated myristoylated alanine-rich C kinase substrate and its subsequent diffusion and dephosphorylation in the cytosol also serves to suppress the pathway in down-gradient regions of the cell. By itself, this mechanism only weakly amplifies signaling in a shallow PDGF gradient, but it synergizes with other feedback mechanisms to enhance amplification. This model offers a framework for a mechanistic understanding of phospholipase C/PKC signaling in chemotactic gradient sensing and can guide the design of experiments to assess the roles of putative feedback loops.

Are Filopodia Privileged Signaling Structures in Migrating Cells?

Filopodia are thin, fingerlike structures that contain bundled actin filaments and project from the cell periphery. These structures are dogmatically endowed with the ability to sense cues in the microenvironment, implying that filopodia foster local signal transduction, yet their small diameter hampers the imaging of dynamic processes therein. To overcome this challenge, we analyzed total internal reflection fluorescence images of migrating fibroblasts coexpressing either a plasma membrane marker or tagged AktPH domain, a translocation biosensor for signaling through the phosphoinositide 3-kinase pathway, along with a cytosolic volume marker. We devised a scheme to estimate the radii of filopodia using either the membrane marker or volume marker data, and we used that information to account for geometry effects in the biosensor data. With conservative estimates of relative target molecule abundance, it is revealed that filopodia typically harbor higher densities of 3' phosphoinositides than adjacent regions at the cell periphery. In this context at least, the analysis supports the filopodial signaling hypothesis.

Data-driven modeling reconciles kinetics of ERK phosphorylation, localization, and activity states.

The extracellular signal-regulated kinase (ERK) signaling pathway controls cell proliferation and differentiation in metazoans. Two hallmarks of its dynamics are adaptation of ERK phosphorylation, which has been linked to negative feedback, and nucleocytoplasmic shuttling, which allows active ERK to phosphorylate protein substrates in the nucleus and cytosol. To integrate these complex features, we acquired quantitative biochemical and live-cell microscopy data to reconcile phosphorylation, localization, and activity states of ERK. While maximal growth factor stimulation elicits transient ERK phosphorylation and nuclear translocation responses, ERK activities available to phosphorylate substrates in the cytosol and nuclei show relatively little or no adaptation. Free ERK activity in the nucleus temporally lags the peak in nuclear translocation, indicating a slow process. Additional experiments, guided by kinetic modeling, show that this process is consistent with ERK's modification of and release from nuclear substrate anchors. Thus, adaptation of whole-cell ERK phosphorylation is a by-product of transient protection from phosphatases. Consistent with this interpretation, predictions concerning the dose-dependence of the pathway response and its interruption by inhibition of MEK were experimentally confirmed.

Deactivation of a negative regulator: a distinct signal transduction mechanism, pronounced in Akt signaling.

Kinase cascades, in which enzymes are sequentially activated by phosphorylation, are quintessential signaling pathways. Signal transduction is not always achieved by direct activation, however. Often, kinases activate pathways by deactivation of a negative regulator; this indirect mechanism, pervasive in Akt signaling, has yet to be systematically explored. Here, we show that the indirect mechanism has properties that are distinct from direct activation. With comparable parameters, the indirect mechanism yields a broader range of sensitivity to the input, beyond saturation of regulator phosphorylation, and kinetics that become progressively slower, not faster, with increasing input strength. These properties can be integrated in network motifs to produce desired responses, as in the case of feedforward loops.

Quantifying collective motion patterns in mesenchymal cell populations using topological data analysis and agent-based modeling.

Fibroblasts in a confluent monolayer are known to adopt elongated morphologies in which cells are oriented parallel to their neighbors. We collected and analyzed new microscopy movies to show that confluent fibroblasts are motile and that neighboring cells often move in anti-parallel directions in a collective motion phenomenon we refer to as "fluidization" of the cell population. We used machine learning to perform cell tracking for each movie and then leveraged topological data analysis (TDA) to show that time-varying point-clouds generated by the tracks contain significant topological information content that is driven by fluidization, i.e., the anti-parallel movement of individual neighboring cells and neighboring groups of cells over long distances. We then utilized the TDA summaries extracted from each movie to perform Bayesian parameter estimation for the D'Orsgona model, an agent-based model (ABM) known to produce a wide array of different patterns, including patterns that are qualitatively similar to fluidization. Although the D'Orsgona ABM is a phenomenological model that only describes inter-cellular attraction and repulsion, the estimated region of D'Orsogna model parameter space was consistent across all movies, suggesting that a specific level of inter-cellular repulsion force at close range may be a mechanism that helps drive fluidization patterns in confluent mesenchymal cell populations.

Data-driven modelling of receptor tyrosine kinase signalling networks quantifies receptor-specific potencies of PI3K- and Ras-dependent ERK activation.

Signal transduction networks in mammalian cells, comprising a limited set of interacting biochemical pathways, are accessed by various growth factor and cytokine receptors to elicit distinct cell responses. This raises the question as to how specificity of the stimulus-response relationship is encoded at the molecular level. It has been proposed that specificity arises not only from the activation of unique signalling pathways, but also from quantitative differences in the activation and regulation of shared receptor-proximal signalling proteins. To address such hypotheses, data sets with greater precision and coverage of experimental conditions will need to be acquired, and rigorous frameworks that codify and parameterize the inherently non-linear relationships among signalling activities will need to be developed. In the present study we apply a systematic approach combining quantitative measurements and mathematical modelling to compare the signalling networks accessed by FGF (fibroblast growth factor) and PDGF (platelet-derived growth factor) receptors in mouse fibroblasts, in which the ERK (extracellular-signal-regulated kinase) cascade is activated by Ras- and PI3K (phosphoinositide 3-kinase)-dependent pathways. We show that, whereas the FGF stimulation of PI3K signalling is relatively weak, this deficiency is compensated for by a more potent Ras-dependent activation of ERK. Thus, as the modelling would predict, the ERK pathway is activated to a greater extent in cells co-stimulated with FGF and PDGF, relative to the saturated levels achieved with either ligand alone. It is envisaged that similar approaches will prove valuable in the elucidation of quantitative differences among other closely related receptor signalling networks.

On the inference of ERK signaling dynamics from protein biosensor measurements.

The extracellular signal-regulated kinase (ERK) signaling pathway plays prominent roles in cell growth, proliferation, and differentiation. ERK signaling is dynamic, involving phosphorylation/dephosphorylation, nucleocytoplasmic shuttling, and interactions with scores of protein substrates in the cytosol and in the nucleus. Live-cell fluorescence microscopy using genetically encoded ERK biosensors offers the potential to infer those dynamics in individual cells. In this study, we have monitored ERK signaling using four commonly used translocation- and Förster resonance energy transfer-based biosensors in a common cell stimulation context. Consistent with previous reports, we found that each biosensor responds with unique kinetics; it is clear that there is not a single dynamic signature characterizing the complexity of ERK phosphorylation, translocation, and kinase activity. In particular, the widely adopted ERK Kinase Translocation Reporter (ERKKTR) gives a readout that reflects ERK activity in both compartments. Mathematical modeling offers an interpretation of the measured ERKKTR kinetics, in relation to cytosolic and nuclear ERK activity, and suggests that biosensor-specific dynamics substantially influence the measured output.

Isoform-specific optical activation of kinase function reveals p38-ERK signaling crosstalk.

Evolution has diversified the mammalian proteome by the generation of protein isoforms that originate from identical genes, e.g., through alternative gene splicing or post-translational modifications, or very similar genes found in gene families. Protein isoforms can have either overlapping or unique functions and traditional chemical, biochemical, and genetic techniques are often limited in their ability to differentiate between isoforms due to their high similarity. This is particularly true in the context of highly dynamic cell signaling cascades, which often require acute spatiotemporal perturbation to assess mechanistic details. To that end, we describe a method for the selective perturbation of the individual protein isoforms of the mitogen-activated protein kinase (MAPK) p38. The genetic installation of a photocaging group at a conserved active site lysine enables the precise light-controlled initiation of kinase signaling, followed by investigation of downstream events. Through optical control, we have identified a novel point of crosstalk between two major signaling cascades: the p38/MAPK pathway and the extracellular signal-regulated kinase (ERK)/MAPK pathway. Specifically, using the photoactivated p38 isoforms, we have found the p38γ and p38δ variants to be positive regulators of the ERK signaling cascade, while confirming the p38α and p38ß variants as negative regulators.

Live-cell fluorescence microscopy with molecular biosensors: what are we really measuring?

Engineered protein biosensors, such as those based on Förster resonance energy transfer, membrane translocation, or solvatochromic shift, are being used in combination with live-cell fluorescence microscopy to reveal kinetics and spatial localization of intracellular processes as they occur. Progress in the application of this approach has been steady, yet its general suitability for quantitative measurements remains unclear. To address the pitfalls of the biosensor approach in quantitative terms, simple reaction-diffusion models were analyzed. The analysis shows that although high-affinity molecular recognition allows robust detection of the fluorescence readout, either of two detrimental effects is fostered. Binding of an intramolecular biosensor or of a relatively abundant intermolecular biosensor introduces observer effects in which the dynamics of the target molecule under study are significantly perturbed, whereas binding of a sparingly expressed intermolecular biosensor is subject to a saturation effect, where the pool of unbound biosensor is significantly depleted. The analysis explores how these effects are manifest in the kinetics and spatial gradients of the biosensor-target complex. A sobering insight emerges: the observer or saturation effect is always significant; the question is whether or not it can be tolerated or accounted for. The challenge in managing the adverse effects is that specification of the biosensor-target affinity to within a certain order of magnitude is required.

Allosteric modulation of Ras-GTP is linked to signal transduction through RAF kinase.

Ras is a key signal transduction protein in the cell. Mutants of Gly(12) and Gln(61) impair GTPase activity and are found prominently in cancers. In wild type Ras-GTP, an allosteric switch promotes disorder to order transition in switch II, placing Gln(61) in the active site. We show that the "on" and "off" conformations of the allosteric switch can also be attained in RasG12V and RasQ61L. Although both mutants have similarly impaired active sites in the on state, RasQ61L stabilizes an anti-catalytic conformation of switch II in the off state of the allosteric switch when bound to Raf. This translates into more potent activation of the MAPK pathway involving Ras, Raf kinase, MEK, and ERK (Ras/Raf/MEK/ERK) in cells transfected with RasQ61L relative to RasG12V. This differential is not observed in the Raf-independent pathway involving Ras, phosphoinositide 3-kinase (PI3K), and Akt (Ras/PI3K/Akt). Using a combination of structural analysis, hydrolysis rates, and experiments in NIH-3T3 cells, we link the allosteric switch to the control of signaling in the Ras/Raf/MEK/ERK pathway, supporting a GTPase-activating protein-independent model for duration of the Ras-Raf complex.

In chemotaxing fibroblasts, both high-fidelity and weakly biased cell movements track the localization of PI3K signaling.

Cell movement biased by a chemical gradient, or chemotaxis, coordinates the recruitment of cells and collective migration of cell populations. During wound healing, chemotaxis of fibroblasts is stimulated by platelet-derived growth factor (PDGF) and certain other chemoattractants. Whereas the immediate PDGF gradient sensing response has been characterized previously at the level of phosphoinositide 3-kinase (PI3K) signaling, the sensitivity of the response at the level of cell migration bias has not yet been studied quantitatively. In this work, we used live-cell total internal reflection fluorescence microscopy to monitor PI3K signaling dynamics and cell movements for extended periods. We show that persistent and properly aligned (i.e., high-fidelity) fibroblast migration does indeed correlate with polarized PI3K signaling; accordingly, this behavior is seen only under conditions of high gradient steepness (>10% across a typical cell length of 50 µm) and a certain range of PDGF concentrations. Under suboptimal conditions, cells execute a random or biased random walk, but nonetheless move in a predictable fashion according to the changing pattern of PI3K signaling. Inhibition of PI3K during chemotaxis is accompanied by loss of both cell-substratum contact and morphological polarity, but after a recovery period, PI3K-inhibited fibroblasts often regain the ability to orient toward the PDGF gradient.