1
|
Belliveau NM, Footer MJ, Akdoǧan E, van Loon AP, Collins SR, Theriot JA. Whole-genome screens reveal regulators of differentiation state and context-dependent migration in human neutrophils. Nat Commun 2023; 14:5770. [PMID: 37723145 PMCID: PMC10507112 DOI: 10.1038/s41467-023-41452-x] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/15/2023] [Accepted: 08/31/2023] [Indexed: 09/20/2023] Open
Abstract
Neutrophils are the most abundant leukocyte in humans and provide a critical early line of defense as part of our innate immune system. We perform a comprehensive, genome-wide assessment of the molecular factors critical to proliferation, differentiation, and cell migration in a neutrophil-like cell line. Through the development of multiple migration screen strategies, we specifically probe directed (chemotaxis), undirected (chemokinesis), and 3D amoeboid cell migration in these fast-moving cells. We identify a role for mTORC1 signaling in cell differentiation, which influences neutrophil abundance, survival, and migratory behavior. Across our individual migration screens, we identify genes involved in adhesion-dependent and adhesion-independent cell migration, protein trafficking, and regulation of the actomyosin cytoskeleton. This genome-wide screening strategy, therefore, provides an invaluable approach to the study of neutrophils and provides a resource that will inform future studies of cell migration in these and other rapidly migrating cells.
Collapse
Affiliation(s)
- Nathan M Belliveau
- Department of Biology and Howard Hughes Medical Institute, University of Washington, Seattle, WA, 98195, USA
| | - Matthew J Footer
- Department of Biology and Howard Hughes Medical Institute, University of Washington, Seattle, WA, 98195, USA
| | - Emel Akdoǧan
- Department of Microbiology and Molecular Genetics, University of California, Davis, Davis, CA, 95616, USA
| | - Aaron P van Loon
- Department of Biology and Howard Hughes Medical Institute, University of Washington, Seattle, WA, 98195, USA
| | - Sean R Collins
- Department of Microbiology and Molecular Genetics, University of California, Davis, Davis, CA, 95616, USA
| | - Julie A Theriot
- Department of Biology and Howard Hughes Medical Institute, University of Washington, Seattle, WA, 98195, USA.
| |
Collapse
|
2
|
Baldwin SA, Haugh JM. Semi-autonomous wound invasion via matrix-deposited, haptotactic cues. J Theor Biol 2023; 568:111506. [PMID: 37094713 PMCID: PMC10393182 DOI: 10.1016/j.jtbi.2023.111506] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/20/2022] [Revised: 04/18/2023] [Accepted: 04/19/2023] [Indexed: 04/26/2023]
Abstract
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.
Collapse
Affiliation(s)
- Scott A Baldwin
- Department of Chemical and Biomolecular Engineering, North Carolina State University, Campus Box 7905, Raleigh, NC 27695, USA
| | - Jason M Haugh
- Department of Chemical and Biomolecular Engineering, North Carolina State University, Campus Box 7905, Raleigh, NC 27695, USA.
| |
Collapse
|
3
|
Subic T, Sbalzarini IF. A Gaussian jump process formulation of the reaction–diffusion master equation enables faster exact stochastic simulations. J Chem Phys 2022; 157:194110. [DOI: 10.1063/5.0123073] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
We propose a Gaussian jump process model on a regular Cartesian lattice for the diffusion part of the Reaction–Diffusion Master Equation (RDME). We derive the resulting Gaussian RDME (GRDME) formulation from analogy with a kernel-based discretization scheme for continuous diffusion processes and quantify the limits of its validity relative to the classic RDME. We then present an exact stochastic simulation algorithm for the GRDME, showing that the accuracies of GRDME and RDME are comparable, but exact simulations of the GRDME require only a fraction of the computational cost of exact RDME simulations. We analyze the origin of this speedup and its scaling with problem dimension. The benchmarks suggest that the GRDME is a particularly beneficial model for diffusion-dominated systems in three dimensional spaces, often occurring in systems biology and cell biology.
Collapse
Affiliation(s)
- Tina Subic
- Technische Universität Dresden, Faculty of Computer Science, Dresden, Germany
- Max Planck Institute of Molecular Cell Biology and Genetics, Dresden, Germany
- Center for Systems Biology Dresden, Dresden, Germany
| | - Ivo F. Sbalzarini
- Technische Universität Dresden, Faculty of Computer Science, Dresden, Germany
- Max Planck Institute of Molecular Cell Biology and Genetics, Dresden, Germany
- Center for Systems Biology Dresden, Dresden, Germany
- Cluster of Excellence Physics of Life, TU Dresden, Dresden, Germany
| |
Collapse
|
4
|
Zhou L, Feng S, Li L, Lü S, Zhang Y, Long M. Two Complementary Signaling Pathways Depict Eukaryotic Chemotaxis: A Mechanochemical Coupling Model. Front Cell Dev Biol 2021; 9:786254. [PMID: 34869388 PMCID: PMC8635958 DOI: 10.3389/fcell.2021.786254] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/30/2021] [Accepted: 10/21/2021] [Indexed: 01/16/2023] Open
Abstract
Many eukaryotic cells, including neutrophils and Dictyostelium cells, are able to undergo correlated random migration in the absence of directional cues while reacting to shallow gradients of chemoattractants with exquisite precision. Although progress has been made with regard to molecular identities, it remains elusive how molecular mechanics are integrated with cell mechanics to initiate and manipulate cell motility. Here, we propose a two dimensional (2D) cell migration model wherein a multilayered dynamic seesaw mechanism is accompanied by a mechanical strain-based inhibition mechanism. In biology, these two mechanisms can be mapped onto the biochemical feedback between phosphoinositides (PIs) and Rho GTPase and the mechanical interplay between filamin A (FLNa) and FilGAP. Cell migration and the accompanying morphological changes are demonstrated in numerical simulations using a particle-spring model, and the diffusion in the cell membrane are simulations using a one dimensional (1D) finite differences method (FDM). The fine balance established between endogenous signaling and a mechanically governed inactivation scheme ensures the endogenous cycle of self-organizing pseudopods, accounting for the correlated random migration. Furthermore, this model cell manifests directional and adaptable responses to shallow graded signaling, depending on the overwhelming effect of the graded stimuli guidance on strain-based inhibition. Finally, the model cell becomes trapped within an obstacle-ridden spatial region, manifesting a shuttle run for local explorations and can chemotactically “escape”, illustrating again the balance required in the complementary signaling pathways.
Collapse
Affiliation(s)
- Lüwen Zhou
- Smart Materials and Advanced Structure Laboratory, School of Mechanical Engineering and Mechanics, Ningbo University, Ningbo Zhejiang, China.,Key Laboratory of Microgravity (National Microgravity Laboratory), and Beijing Key Laboratory of Engineered Construction and Mechanobiology, Center for Biomechanics and Bioengineering, Institute of Mechanics, Chinese Academy of Sciences, Beijing, China
| | - Shiliang Feng
- Smart Materials and Advanced Structure Laboratory, School of Mechanical Engineering and Mechanics, Ningbo University, Ningbo Zhejiang, China.,Key Laboratory of Microgravity (National Microgravity Laboratory), and Beijing Key Laboratory of Engineered Construction and Mechanobiology, Center for Biomechanics and Bioengineering, Institute of Mechanics, Chinese Academy of Sciences, Beijing, China
| | - Long Li
- State Key Laboratory of Nonlinear Mechanics (LNM) and Beijing Key Laboratory of Engineered Construction and Mechanobiology, Institute of Mechanics, Chinese Academy of Sciences, Beijing, China
| | - Shouqin Lü
- Key Laboratory of Microgravity (National Microgravity Laboratory), and Beijing Key Laboratory of Engineered Construction and Mechanobiology, Center for Biomechanics and Bioengineering, Institute of Mechanics, Chinese Academy of Sciences, Beijing, China.,School of Engineering Science, University of Chinese Academy of Sciences, Beijing, China
| | - Yan Zhang
- Key Laboratory of Microgravity (National Microgravity Laboratory), and Beijing Key Laboratory of Engineered Construction and Mechanobiology, Center for Biomechanics and Bioengineering, Institute of Mechanics, Chinese Academy of Sciences, Beijing, China.,School of Engineering Science, University of Chinese Academy of Sciences, Beijing, China
| | - Mian Long
- Key Laboratory of Microgravity (National Microgravity Laboratory), and Beijing Key Laboratory of Engineered Construction and Mechanobiology, Center for Biomechanics and Bioengineering, Institute of Mechanics, Chinese Academy of Sciences, Beijing, China.,School of Engineering Science, University of Chinese Academy of Sciences, Beijing, China
| |
Collapse
|
5
|
Three-dimensional stochastic simulation of chemoattractant-mediated excitability in cells. PLoS Comput Biol 2021; 17:e1008803. [PMID: 34260581 PMCID: PMC8330952 DOI: 10.1371/journal.pcbi.1008803] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/29/2021] [Revised: 08/03/2021] [Accepted: 06/08/2021] [Indexed: 01/21/2023] Open
Abstract
During the last decade, a consensus has emerged that the stochastic triggering of an excitable system drives pseudopod formation and subsequent migration of amoeboid cells. The presence of chemoattractant stimuli alters the threshold for triggering this activity and can bias the direction of migration. Though noise plays an important role in these behaviors, mathematical models have typically ignored its origin and merely introduced it as an external signal into a series of reaction-diffusion equations. Here we consider a more realistic description based on a reaction-diffusion master equation formalism to implement these networks. In this scheme, noise arises naturally from a stochastic description of the various reaction and diffusion terms. Working on a three-dimensional geometry in which separate compartments are divided into a tetrahedral mesh, we implement a modular description of the system, consisting of G-protein coupled receptor signaling (GPCR), a local excitation-global inhibition mechanism (LEGI), and signal transduction excitable network (STEN). Our models implement detailed biochemical descriptions whenever this information is available, such as in the GPCR and G-protein interactions. In contrast, where the biochemical entities are less certain, such as the LEGI mechanism, we consider various possible schemes and highlight the differences between them. Our simulations show that even when the LEGI mechanism displays perfect adaptation in terms of the mean level of proteins, the variance shows a dose-dependence. This differs between the various models considered, suggesting a possible means for determining experimentally among the various potential networks. Overall, our simulations recreate temporal and spatial patterns observed experimentally in both wild-type and perturbed cells, providing further evidence for the excitable system paradigm. Moreover, because of the overall importance and ubiquity of the modules we consider, including GPCR signaling and adaptation, our results will be of interest beyond the field of directed migration. Though the term noise usually carries negative connotations, it can also contribute positively to the characteristic dynamics of a system. In biological systems, where noise arises from the stochastic interactions between molecules, its study is usually confined to genetic regulatory systems in which copy numbers are small and fluctuations large. However, noise can have important roles when the number of signaling molecules is large. The extension of pseudopods and the subsequent motion of amoeboid cells arises from the noise-induced trigger of an excitable system. Chemoattractant signals bias this triggering thereby directing cell motion. To date, this paradigm has not been tested by mathematical models that account accurately for the noise that arises in the corresponding reactions. In this study, we employ a reaction-diffusion master equation approach to investigate the effects of noise. Using a modular approach and a three-dimensional cell model with specific subdomains attributed to the cell membrane and cortex, we explore the spatiotemporal dynamics of the system. Our simulations recreate many experimentally-observed cell behaviors thereby supporting the biased-excitable network hypothesis.
Collapse
|
6
|
Savage NS. Describing the movement of molecules in reduced-dimension models. Commun Biol 2021; 4:689. [PMID: 34099856 PMCID: PMC8184792 DOI: 10.1038/s42003-021-02200-3] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/13/2021] [Accepted: 05/10/2021] [Indexed: 12/05/2022] Open
Abstract
When addressing spatial biological questions using mathematical models, symmetries within the system are often exploited to simplify the problem by reducing its physical dimension. In a reduced-dimension model molecular movement is restricted to the reduced dimension, changing the nature of molecular movement. This change in molecular movement can lead to quantitatively and even qualitatively different results in the full and reduced systems. Within this manuscript we discuss the condition under which restricted molecular movement in reduced-dimension models accurately approximates molecular movement in the full system. For those systems which do not satisfy the condition, we present a general method for approximating unrestricted molecular movement in reduced-dimension models. We will derive a mathematically robust, finite difference method for solving the 2D diffusion equation within a 1D reduced-dimension model. The methods described here can be used to improve the accuracy of many reduced-dimension models while retaining benefits of system simplification.
Collapse
|
7
|
Floris E, Piras A, Dall’Asta L, Gamba A, Hirsch E, Campa CC. Physics of compartmentalization: How phase separation and signaling shape membrane and organelle identity. Comput Struct Biotechnol J 2021; 19:3225-3233. [PMID: 34141141 PMCID: PMC8190439 DOI: 10.1016/j.csbj.2021.05.029] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/31/2021] [Accepted: 05/15/2021] [Indexed: 11/29/2022] Open
Abstract
Compartmentalization of cellular functions is at the core of the physiology of eukaryotic cells. Recent evidences indicate that a universal organizing process - phase separation - supports the partitioning of biomolecules in distinct phases from a single homogeneous mixture, a landmark event in both the biogenesis and the maintenance of membrane and non-membrane-bound organelles. In the cell, 'passive' (non energy-consuming) mechanisms are flanked by 'active' mechanisms of separation into phases of distinct density and stoichiometry, that allow for increased partitioning flexibility and programmability. A convergence of physical and biological approaches is leading to new insights into the inner functioning of this driver of intracellular order, holding promises for future advances in both biological research and biotechnological applications.
Collapse
Affiliation(s)
- Elisa Floris
- Institute of Condensed Matter Physics and Complex Systems, Department of Applied Science and Technology, Politecnico di Torino, Corso Duca degli Abruzzi 24, 10129 Torino, Italy
| | - Andrea Piras
- Italian Institute for Genomic Medicine (IIGM), c/o IRCCS Candiolo, Str.Prov.le 142, km 3.95, Candiolo (TO) 10060, Italy
- Candiolo Cancer Institute, FPO - IRCCS, Str.Prov.le 142, km 3.95, Candiolo (TO) 10060, Italy
| | - Luca Dall’Asta
- Institute of Condensed Matter Physics and Complex Systems, Department of Applied Science and Technology, Politecnico di Torino, Corso Duca degli Abruzzi 24, 10129 Torino, Italy
- Collegio Carlo Alberto, Piazza Arbarello 8, 10122 Torino, Italy
| | - Andrea Gamba
- Institute of Condensed Matter Physics and Complex Systems, Department of Applied Science and Technology, Politecnico di Torino, Corso Duca degli Abruzzi 24, 10129 Torino, Italy
- Italian Institute for Genomic Medicine (IIGM), c/o IRCCS Candiolo, Str.Prov.le 142, km 3.95, Candiolo (TO) 10060, Italy
- Candiolo Cancer Institute, FPO - IRCCS, Str.Prov.le 142, km 3.95, Candiolo (TO) 10060, Italy
- Istituto Nazionale di Fisica Nucleare (INFN), sezione di Torino, Via Giuria 1, 10125 Torino, Italy
| | - Emilio Hirsch
- Department of Molecular Biotechnology and Health Sciences, Molecular Biotechnology Center, University of Torino, Via Nizza 52, 10126 Torino, Italy
| | - Carlo C. Campa
- Italian Institute for Genomic Medicine (IIGM), c/o IRCCS Candiolo, Str.Prov.le 142, km 3.95, Candiolo (TO) 10060, Italy
- Candiolo Cancer Institute, FPO - IRCCS, Str.Prov.le 142, km 3.95, Candiolo (TO) 10060, Italy
| |
Collapse
|
8
|
Dirks C, Striewski P, Wirth B, Aalto A, Olguin-Olguin A. A mathematical model for bleb regulation in zebrafish primordial germ cells. MATHEMATICAL MEDICINE AND BIOLOGY-A JOURNAL OF THE IMA 2021; 38:218-254. [PMID: 33601409 DOI: 10.1093/imammb/dqab002] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/19/2019] [Revised: 01/12/2021] [Accepted: 01/24/2021] [Indexed: 01/02/2023]
Abstract
Blebs are cell protrusions generated by local membrane-cortex detachments followed by expansion of the plasma membrane. Blebs are formed by some migrating cells, e.g. primordial germ cells of the zebrafish. While blebs occur randomly at each part of the membrane in unpolarized cells, a polarization process guarantees the occurrence of blebs at a preferential site and thereby facilitates migration toward a specified direction. Little is known about the factors involved in the controlled and directed bleb generation, yet recent studies revealed the influence of an intracellular flow and the stabilizing role of the membrane-cortex linker molecule Ezrin. Based on this information, we develop and analyse a coupled bulk-surface model describing a potential cellular mechanism by which a bleb could be induced at a controlled site. The model rests upon intracellular Darcy flow and a diffusion-advection-reaction system, describing the temporal evolution from a homogeneous to a strongly anisotropic Ezrin distribution. We prove the well-posedness of the mathematical model and show that simulations qualitatively correspond to experimental observations, suggesting that indeed the interaction of an intracellular flow with membrane proteins can be the cause of the Ezrin redistribution accompanying bleb formation.
Collapse
Affiliation(s)
- Carolin Dirks
- WWU Münster FB 10 Mathematik und Informatik, Institute for Analysis and Numerics, 48149 Münster, Germany
| | - Paul Striewski
- WWU Münster FB 10 Mathematik und Informatik, Institute for Analysis and Numerics, 48149 Münster, Germany
| | - Benedikt Wirth
- WWU Münster FB 10 Mathematik und Informatik, Institute for Analysis and Numerics, 48149 Münster, Germany
| | - Anne Aalto
- WWU Münster FB 13 Biologie, Institute of Cell Biology, Center for Molecular Biology of Inflammation, 48149 Münster, Germany
| | - Adan Olguin-Olguin
- WWU Münster FB 13 Biologie, Institute of Cell Biology, Center for Molecular Biology of Inflammation, 48149 Münster, Germany
| |
Collapse
|
9
|
Burbano-L. DA, Porfiri M. Modeling multi-sensory feedback control of zebrafish in a flow. PLoS Comput Biol 2021; 17:e1008644. [PMID: 33481795 PMCID: PMC7857640 DOI: 10.1371/journal.pcbi.1008644] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/17/2020] [Revised: 02/03/2021] [Accepted: 12/18/2020] [Indexed: 11/18/2022] Open
Abstract
Understanding how animals navigate complex environments is a fundamental challenge in biology and a source of inspiration for the design of autonomous systems in engineering. Animal orientation and navigation is a complex process that integrates multiple senses, whose function and contribution are yet to be fully clarified. Here, we propose a data-driven mathematical model of adult zebrafish engaging in counter-flow swimming, an innate behavior known as rheotaxis. Zebrafish locomotion in a two-dimensional fluid flow is described within the finite-dipole model, which consists of a pair of vortices separated by a constant distance. The strength of these vortices is adjusted in real time by the fish to afford orientation and navigation control, in response to of the multi-sensory input from vision, lateral line, and touch. Model parameters for the resulting stochastic differential equations are calibrated through a series of experiments, in which zebrafish swam in a water channel under different illumination conditions. The accuracy of the model is validated through the study of a series of measures of rheotactic behavior, contrasting results of real and in-silico experiments. Our results point at a critical role of hydromechanical feedback during rheotaxis, in the form of a gradient-following strategy.
Collapse
Affiliation(s)
- Daniel A. Burbano-L.
- Department of Mechanical and Aerospace Engineering, Tandon School of Engineering, New York University, New York City, New York, USA
| | - Maurizio Porfiri
- Department of Mechanical and Aerospace Engineering, Tandon School of Engineering, New York University, New York City, New York, USA
- Department of Biomedical Engineering, Tandon School of Engineering, New York University, New York City, New York, USA
- Center for Urban Sciences and Progress, Tandon School of Engineering, New York University, New York City, New York, USA
- * E-mail:
| |
Collapse
|
10
|
Johnson ME, Chen A, Faeder JR, Henning P, Moraru II, Meier-Schellersheim M, Murphy RF, Prüstel T, Theriot JA, Uhrmacher AM. Quantifying the roles of space and stochasticity in computer simulations for cell biology and cellular biochemistry. Mol Biol Cell 2021; 32:186-210. [PMID: 33237849 PMCID: PMC8120688 DOI: 10.1091/mbc.e20-08-0530] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/12/2020] [Revised: 10/13/2020] [Accepted: 11/17/2020] [Indexed: 12/29/2022] Open
Abstract
Most of the fascinating phenomena studied in cell biology emerge from interactions among highly organized multimolecular structures embedded into complex and frequently dynamic cellular morphologies. For the exploration of such systems, computer simulation has proved to be an invaluable tool, and many researchers in this field have developed sophisticated computational models for application to specific cell biological questions. However, it is often difficult to reconcile conflicting computational results that use different approaches to describe the same phenomenon. To address this issue systematically, we have defined a series of computational test cases ranging from very simple to moderately complex, varying key features of dimensionality, reaction type, reaction speed, crowding, and cell size. We then quantified how explicit spatial and/or stochastic implementations alter outcomes, even when all methods use the same reaction network, rates, and concentrations. For simple cases, we generally find minor differences in solutions of the same problem. However, we observe increasing discordance as the effects of localization, dimensionality reduction, and irreversible enzymatic reactions are combined. We discuss the strengths and limitations of commonly used computational approaches for exploring cell biological questions and provide a framework for decision making by researchers developing new models. As computational power and speed continue to increase at a remarkable rate, the dream of a fully comprehensive computational model of a living cell may be drawing closer to reality, but our analysis demonstrates that it will be crucial to evaluate the accuracy of such models critically and systematically.
Collapse
Affiliation(s)
- M. E. Johnson
- Thomas C. Jenkins Department of Biophysics, Johns Hopkins University, Baltimore, MD, 21218
| | - A. Chen
- Thomas C. Jenkins Department of Biophysics, Johns Hopkins University, Baltimore, MD, 21218
| | - J. R. Faeder
- Department of Computational and Systems Biology, University of Pittsburgh School of Medicine, Pittsburgh, PA, 15260
| | - P. Henning
- Institute for Visual and Analytic Computing, University of Rostock, 18055 Rostock, Germany
| | - I. I. Moraru
- Department of Cell Biology, Center for Cell Analysis and Modeling, University of Connecticut Health Center, Farmington, CT 06030
| | - M. Meier-Schellersheim
- Laboratory of Immune System Biology, National Institute of Allergy and Infectious Diseases, National Institutes of Health, Bethesda, MD 20892
| | - R. F. Murphy
- Computational Biology Department, Department of Biological Sciences, Department of Biomedical Engineering, Machine Learning Department, Carnegie Mellon University, Pittsburgh, PA 15289
| | - T. Prüstel
- Laboratory of Immune System Biology, National Institute of Allergy and Infectious Diseases, National Institutes of Health, Bethesda, MD 20892
| | - J. A. Theriot
- Department of Biology and Howard Hughes Medical Institute, University of Washington, Seattle, WA 98195
| | - A. M. Uhrmacher
- Institute for Visual and Analytic Computing, University of Rostock, 18055 Rostock, Germany
| |
Collapse
|
11
|
Kopfer KH, Jäger W, Matthäus F. A mechanochemical model for rho GTPase mediated cell polarization. J Theor Biol 2020; 504:110386. [DOI: 10.1016/j.jtbi.2020.110386] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/16/2020] [Revised: 06/15/2020] [Accepted: 06/17/2020] [Indexed: 01/13/2023]
|
12
|
Copos C, Mogilner A. A hybrid stochastic-deterministic mechanochemical model of cell polarization. Mol Biol Cell 2020; 31:1637-1649. [PMID: 32459563 PMCID: PMC7521800 DOI: 10.1091/mbc.e19-09-0549] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/27/2019] [Revised: 05/11/2020] [Accepted: 05/20/2020] [Indexed: 12/15/2022] Open
Abstract
Polarization is a crucial component in cell differentiation, development, and motility, but its details are not yet well understood. At the onset of cell locomotion, cells break symmetry to form well-defined cell fronts and rears. This polarity establishment varies across cell types: in Dictyostelium discoideum cells, it is mediated by biochemical signaling pathways and can function in the absence of a cytoskeleton, while in keratocytes, it is tightly connected to cytoskeletal dynamics and mechanics. Theoretical models that have been developed to understand the onset of polarization have explored either signaling or mechanical pathways, yet few have explored mechanochemical mechanisms. However, many motile cells rely on both signaling modules and actin cytoskeleton to break symmetry and achieve a stable polarized state. We propose a general mechanochemical polarization model based on coupling between a stochastic model for the segregation of signaling molecules and a simplified mechanical model for actin cytoskeleton network competition. We find that local linear coupling between minimally nonlinear signaling and cytoskeletal systems, separately not supporting stable polarization, yields a robustly polarized cell state. The model captures the essence of spontaneous polarization of neutrophils, which has been proposed to emerge due to the competition between frontness and backness pathways.
Collapse
Affiliation(s)
- Calina Copos
- Courant Institute, New York University, New York, NY 10012
| | - Alex Mogilner
- Courant Institute, New York University, New York, NY 10012
- Department of Biology, New York University, New York, NY 10012
| |
Collapse
|
13
|
Feng S, Zhou L, Zhang Y, Lü S, Long M. Mechanochemical modeling of neutrophil migration based on four signaling layers, integrin dynamics, and substrate stiffness. Biomech Model Mechanobiol 2018; 17:1611-1630. [PMID: 29968162 DOI: 10.1007/s10237-018-1047-2] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/21/2018] [Accepted: 06/24/2018] [Indexed: 01/09/2023]
Abstract
Directional neutrophil migration during human immune responses is a highly coordinated process regulated by both biochemical and biomechanical environments. In this paper, we developed an integrative mathematical model of neutrophil migration using a lattice Boltzmann-particle method built in-house to solve the moving boundary problem with spatiotemporal regulation of biochemical components. The mechanical features of the cell cortex are modeled by a series of spring-connected nodes representing discrete cell-substrate adhesive sites. The intracellular signaling cascades responsible for cytoskeletal remodeling [e.g., small GTPases, phosphoinositide-3-kinase (PI3K), and phosphatase and tensin homolog] are built based on our previous four-layered signaling model centered on the bidirectional molecular transport mechanism and implemented as reaction-diffusion equations. Focal adhesion dynamics are determined by force-dependent integrin-ligand binding kinetics and integrin recycling and are thus integrated with cell motion. Using numerical simulations, the model reproduces the major features of cell migration in response to uniform and gradient biochemical stimuli based on the quantitative spatiotemporal regulation of signaling molecules, which agree with experimental observations. The existence of multiple types of integrins with different binding kinetics could act as an adaptation mechanism for substrate stiffness. Moreover, cells can perform reversal, U-turn, or lock-on behaviors depending on the steepness of the reversal biochemical signals received. Finally, this model is also applied to predict the responses of mutants in which PTEN is overexpressed or disrupted.
Collapse
Affiliation(s)
- Shiliang Feng
- Center for Biomechanics and Bioengineering, Key Laboratory of Microgravity (National Microgravity Laboratory), and Beijing Key Laboratory of Engineered Construction and Mechanobiology, Institute of Mechanics, Chinese Academy of Sciences, Beijing, China
- School of Engineering Science, University of Chinese Academy of Sciences, Beijing, China
| | - Lüwen Zhou
- Center for Biomechanics and Bioengineering, Key Laboratory of Microgravity (National Microgravity Laboratory), and Beijing Key Laboratory of Engineered Construction and Mechanobiology, Institute of Mechanics, Chinese Academy of Sciences, Beijing, China
- School of Engineering Science, University of Chinese Academy of Sciences, Beijing, China
| | - Yan Zhang
- Center for Biomechanics and Bioengineering, Key Laboratory of Microgravity (National Microgravity Laboratory), and Beijing Key Laboratory of Engineered Construction and Mechanobiology, Institute of Mechanics, Chinese Academy of Sciences, Beijing, China
- School of Engineering Science, University of Chinese Academy of Sciences, Beijing, China
| | - Shouqin Lü
- Center for Biomechanics and Bioengineering, Key Laboratory of Microgravity (National Microgravity Laboratory), and Beijing Key Laboratory of Engineered Construction and Mechanobiology, Institute of Mechanics, Chinese Academy of Sciences, Beijing, China
- School of Engineering Science, University of Chinese Academy of Sciences, Beijing, China
| | - Mian Long
- Center for Biomechanics and Bioengineering, Key Laboratory of Microgravity (National Microgravity Laboratory), and Beijing Key Laboratory of Engineered Construction and Mechanobiology, Institute of Mechanics, Chinese Academy of Sciences, Beijing, China.
- School of Engineering Science, University of Chinese Academy of Sciences, Beijing, China.
| |
Collapse
|
14
|
Feng SL, Zhou LW, Lü SQ, Zhang Y. Dynamic seesaw model for rapid signaling responses in eukaryotic chemotaxis. Phys Biol 2018; 15:056004. [PMID: 29757152 DOI: 10.1088/1478-3975/aac45b] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023]
Abstract
Directed movement of eukaryotic cells toward spatiotemporally varied chemotactic stimuli enables rapid intracellular signaling responses. While macroscopic cellular manifestation is shaped by balancing external stimuli strength with finite internal delays, the organizing principles of the underlying molecular mechanisms remain to be clarified. Here, we developed a novel modeling framework based on a simple seesaw mechanism to elucidate how cells repeatedly reverse polarity. As a key feature of the modeling, the bottom module of bidirectional molecular transport is successively controlled by three upstream modules of signal reception, initial signal processing, and Rho GTPase regulation. Our simulations indicated that an isotropic cell is polarized in response to a graded input signal. By applying a reversal gradient to a chemoattractant signal, lamellipod-specific molecules (i.e. PIP3 and PI3K) disappear, first from the cell front, and then they redistribute at the opposite side, whereas functional molecules at the rear of the cell (i.e. PIP2 and PTEN) act oppositely. In particular, the model cell exhibits a seesaw-like spatiotemporal pattern for the establishment of front and rear and interconversion, consistent with those related experimental observations. Increasing the switching frequency of the chemotactic gradient causes the cell to stay in a trapped state, further supporting the proposed dynamics of eukaryotic chemotaxis with the underlying cytoskeletal remodeling.
Collapse
Affiliation(s)
- Shi Liang Feng
- Institute of mechanical engineering and mechanics, Ningbo University, Ningbo 315211, People's Republic of China. Center of Biomechanics and Bioengineering and Institute of Mechanics, Chinese Academy of Sciences, Beijing 100190, People's Republic of China
| | | | | | | |
Collapse
|
15
|
Sun M, Spill F, Zaman MH. A Computational Model of YAP/TAZ Mechanosensing. Biophys J 2017; 110:2540-2550. [PMID: 27276271 DOI: 10.1016/j.bpj.2016.04.040] [Citation(s) in RCA: 43] [Impact Index Per Article: 6.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/23/2015] [Revised: 03/24/2016] [Accepted: 04/12/2016] [Indexed: 12/15/2022] Open
Abstract
In cell proliferation, stem cell differentiation, chemoresistance, and tissue organization, the ubiquitous role of YAP/TAZ continues to impact our fundamental understanding in numerous physiological and disease systems. YAP/TAZ is an important signaling nexus integrating diverse mechanical and biochemical signals, such as ECM stiffness, adhesion ligand density, or cell-cell contacts, and thus strongly influences cell fate. Recent studies show that YAP/TAZ mechanical sensing is dependent on RhoA-regulated stress fibers. However, current understanding of YAP/TAZ remains limited due to the unknown interaction between the canonical Hippo pathway and cell tension. Furthermore, the multiscale relationship connecting adhesion signaling to YAP/TAZ activity through cytoskeleton dynamics remains poorly understood. To identify the roles of key signaling molecules in mechanical signal sensing and transduction, we present a, to our knowledge, novel computational model of the YAP/TAZ signaling pathway. This model converts extracellular-matrix mechanical properties to biochemical signals via adhesion, and integrates intracellular signaling cascades associated with cytoskeleton dynamics. We perform perturbations of molecular levels and sensitivity analyses to predict how various signaling molecules affect YAP/TAZ activity. Adhesion molecules, such as FAK, are predicted to rescue YAP/TAZ activity in soft environments via the RhoA pathway. We also found that changes of molecule concentrations result in different patterns of YAP/TAZ stiffness response. We also investigate the sensitivity of YAP/TAZ activity to ECM stiffness, and compare with that of SRF/MAL, which is another important regulator of differentiation. In addition, the model shows that the unresolved synergistic effect of YAP/TAZ activity between the mechanosensing and the Hippo pathways can be explained by the interaction of LIM-kinase and LATS. Overall, our model provides a, to our knowledge, novel platform for studying YAP/TAZ activity in the context of integrating different signaling pathways. This platform can be used to gain, to our knowledge, new fundamental insights into roles of key molecular and mechanical regulators on development, tissue engineering, or tumor progression.
Collapse
Affiliation(s)
- Meng Sun
- Department of Biomedical Engineering, Boston University, Boston, Massachusetts
| | - Fabian Spill
- Department of Biomedical Engineering, Boston University, Boston, Massachusetts; Department of Mechanical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts.
| | - Muhammad H Zaman
- Department of Biomedical Engineering, Boston University, Boston, Massachusetts; Howard Hughes Medical Institute, Boston University, Boston, Massachusetts.
| |
Collapse
|
16
|
Abstract
Cell polarization is a key step in the migration, development, and organization of eukaryotic cells, both at the single cell and multicellular level. Research on the mechanisms that give rise to polarization of a given cell, and organization of polarity within a tissue has led to new understanding across cellular and developmental biology. In this review, we describe some of the history of theoretical and experimental aspects of the field, as well as some interesting questions and challenges for the future.
Collapse
Affiliation(s)
- Wouter-Jan Rappel
- Department of Physics, University of California, San Diego, La Jolla, USA
| | | |
Collapse
|
17
|
Mathematical Models for Immunology: Current State of the Art and Future Research Directions. Bull Math Biol 2016; 78:2091-2134. [PMID: 27714570 PMCID: PMC5069344 DOI: 10.1007/s11538-016-0214-9] [Citation(s) in RCA: 79] [Impact Index Per Article: 9.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/03/2016] [Accepted: 09/26/2016] [Indexed: 01/01/2023]
Abstract
The advances in genetics and biochemistry that have taken place over the last 10 years led to significant advances in experimental and clinical immunology. In turn, this has led to the development of new mathematical models to investigate qualitatively and quantitatively various open questions in immunology. In this study we present a review of some research areas in mathematical immunology that evolved over the last 10 years. To this end, we take a step-by-step approach in discussing a range of models derived to study the dynamics of both the innate and immune responses at the molecular, cellular and tissue scales. To emphasise the use of mathematics in modelling in this area, we also review some of the mathematical tools used to investigate these models. Finally, we discuss some future trends in both experimental immunology and mathematical immunology for the upcoming years.
Collapse
|
18
|
Kim MO, Ryu JM, Suh HN, Park SH, Oh YM, Lee SH, Han HJ. cAMP Promotes Cell Migration Through Cell Junctional Complex Dynamics and Actin Cytoskeleton Remodeling: Implications in Skin Wound Healing. Stem Cells Dev 2015; 24:2513-24. [PMID: 26192163 DOI: 10.1089/scd.2015.0130] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023] Open
Abstract
Stem cells have attracted great interest for their therapeutic capacity in tissue regeneration. Cyclic adenosine 3',5'-monophosphate (cAMP), existing in high concentration at wound sites, mediated various signaling pathways such as cytoskeleton dynamics, cell adhesion, and cell migration in stem cells, which suggest the critical roles of cAMP in the wound healing process through functional regulation of stem cells. However, the mechanisms behind the effect of cAMP on mouse embryonic stem cell (mESC) motility and its roles on skin wound healing remain to be fully elucidated. In the present study, 8-Bromo cAMP-treated mESCs showed significant wound closure and improved neovascularization. Moreover, 8-Bromo cAMP stimulated mESC migration into the wound bed. 8-Bromo cAMP also increased ESC motility in in vitro migration assay. 8-Bromo cAMP induced myosin light chain phosphorylation through Rac1 and Cdc42 signaling, which were involved in 8-Bromo cAMP-induced decrease in expression of junction proteins (connexin 43, E-cadherin, and occludin) at the plasma membrane. Subsequently, 8-Bromo cAMP induced the disruption of cell junctions (including gap junctions, adherens junctions, and tight junctions), which reduced the function of the gap junctions and cell adhesion. In addition, 8-Bromo cAMP-induced Rac1 and Cdc42 activation increased Arp3, TOCA, PAK, and N-WASP expression, but decreased cofilin phosphorylation level, which elicited actin cytoskeleton remodeling. In contrast to the control, 8-Bromo cAMP evoked a substantial migration of cells into the denuded area, which was blocked by the small interfering RNAs of the signaling pathway-related molecules or by inhibitors. In conclusion, cAMP enhanced the migration of mESCs through effective coordination of junctional disruption and actin cytoskeleton remodeling, which increased the wound healing capacity of ESCs.
Collapse
Affiliation(s)
- Mi Ok Kim
- 1 Department of Veterinary Physiology, College of Veterinary Medicine and Research Institute for Veterinary Science, Seoul National University , Seoul, Republic of Korea.,2 BK21 PLUS Creative Veterinary Research Center, Seoul National University , Seoul, Republic of Korea
| | - Jung Min Ryu
- 1 Department of Veterinary Physiology, College of Veterinary Medicine and Research Institute for Veterinary Science, Seoul National University , Seoul, Republic of Korea.,2 BK21 PLUS Creative Veterinary Research Center, Seoul National University , Seoul, Republic of Korea
| | - Han Na Suh
- 1 Department of Veterinary Physiology, College of Veterinary Medicine and Research Institute for Veterinary Science, Seoul National University , Seoul, Republic of Korea.,2 BK21 PLUS Creative Veterinary Research Center, Seoul National University , Seoul, Republic of Korea
| | - Soo Hyun Park
- 3 College of Veterinary Medicine, Chonnam National University , Gwangju, Republic of Korea
| | - Yeon-Mok Oh
- 4 Department of Pulmonary and Critical Care Medicine, and Clinical Research Center for Chronic Obstructive Airway Diseases, Asan Medical Center, University of Ulsan College of Medicine , Seoul, Republic of Korea
| | - Sang Hun Lee
- 5 Medical Science Research Institute, Soonchunhyang University Seoul Hospital , Seoul, Republic of Korea
| | - Ho Jae Han
- 1 Department of Veterinary Physiology, College of Veterinary Medicine and Research Institute for Veterinary Science, Seoul National University , Seoul, Republic of Korea.,2 BK21 PLUS Creative Veterinary Research Center, Seoul National University , Seoul, Republic of Korea
| |
Collapse
|
19
|
On the Interaction between a Nanoparticulate System and the Human Body in Body Area Nanonetworks. MICROMACHINES 2015. [DOI: 10.3390/mi6091213] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
|
20
|
Zamparo M, Chianale F, Tebaldi C, Cosentino-Lagomarsino M, Nicodemi M, Gamba A. Dynamic membrane patterning, signal localization and polarity in living cells. SOFT MATTER 2015; 11:838-849. [PMID: 25563791 DOI: 10.1039/c4sm02157f] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/04/2023]
Abstract
We review the molecular and physical aspects of the dynamic localization of signaling molecules on the plasma membranes of living cells. At the nanoscale, clusters of receptors and signaling proteins play an essential role in the processing of extracellular signals. At the microscale, "soft" and highly dynamic signaling domains control the interaction of individual cells with their environment. At the multicellular scale, individual polarity patterns control the forces that shape multicellular aggregates and tissues.
Collapse
Affiliation(s)
- M Zamparo
- Human Genetics Foundation - Torino, Italy.
| | | | | | | | | | | |
Collapse
|
21
|
Liu X, Yang T, Suzuki K, Tsukita S, Ishii M, Zhou S, Wang G, Cao L, Qian F, Taylor S, Oh MJ, Levitan I, Ye RD, Carnegie GK, Zhao Y, Malik AB, Xu J. Moesin and myosin phosphatase confine neutrophil orientation in a chemotactic gradient. ACTA ACUST UNITED AC 2015; 212:267-80. [PMID: 25601651 PMCID: PMC4322047 DOI: 10.1084/jem.20140508] [Citation(s) in RCA: 33] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022]
Abstract
Jingsong Xu and colleagues investigate how neutrophils initiate polarized migration toward bacteria or chemoattractants. They find that attractant-induced activation of myosin phosphatase results in the deactivation of moesin at the prospective leading edge and its redistribution to the trailing edge, establishing polarity and directional pseudopod formation. Neutrophils respond to invading bacteria by adopting a polarized morphology, migrating in the correct direction, and engulfing the bacteria. How neutrophils establish and precisely orient this polarity toward pathogens remains unclear. Here we report that in resting neutrophils, the ERM (ezrin, radixin, and moesin) protein moesin in its active form (phosphorylated and membrane bound) prevented cell polarization by inhibiting the small GTPases Rac, Rho, and Cdc42. Attractant-induced activation of myosin phosphatase deactivated moesin at the prospective leading edge to break symmetry and establish polarity. Subsequent translocation of moesin to the trailing edge confined the formation of a prominent pseudopod directed toward pathogens and prevented secondary pseudopod formation in other directions. Therefore, both moesin-mediated inhibition and its localized deactivation by myosin phosphatase are essential for neutrophil polarization and effective neutrophil tracking of pathogens.
Collapse
Affiliation(s)
- Xiaowen Liu
- Department of Pharmacology and Department of Medicine, University of Illinois, Chicago, IL 60612
| | - Tao Yang
- State Key Laboratory of Biomembrane and Membrane Biotechnology, Institute of Zoology, Chinese Academy of Sciences, Beijing 100101, China
| | - Koya Suzuki
- Laboratory of Biological Science and Laboratory of Cellular Dynamics, Graduate School of Medicine, Osaka University, Suita 565-0871, Osaka, Japan
| | - Sachiko Tsukita
- Laboratory of Biological Science and Laboratory of Cellular Dynamics, Graduate School of Medicine, Osaka University, Suita 565-0871, Osaka, Japan
| | - Masaru Ishii
- Laboratory of Biological Science and Laboratory of Cellular Dynamics, Graduate School of Medicine, Osaka University, Suita 565-0871, Osaka, Japan
| | - Shuping Zhou
- State Key Laboratory of Biomembrane and Membrane Biotechnology, Institute of Zoology, Chinese Academy of Sciences, Beijing 100101, China
| | - Gang Wang
- Department of Pharmacology and Department of Medicine, University of Illinois, Chicago, IL 60612
| | - Luyang Cao
- Department of Pharmacology and Department of Medicine, University of Illinois, Chicago, IL 60612
| | - Feng Qian
- Department of Pharmacology and Department of Medicine, University of Illinois, Chicago, IL 60612
| | - Shalina Taylor
- Department of Pharmacology and Department of Medicine, University of Illinois, Chicago, IL 60612
| | - Myung-Jin Oh
- Department of Pharmacology and Department of Medicine, University of Illinois, Chicago, IL 60612
| | - Irena Levitan
- Department of Pharmacology and Department of Medicine, University of Illinois, Chicago, IL 60612
| | - Richard D Ye
- Department of Pharmacology and Department of Medicine, University of Illinois, Chicago, IL 60612
| | - Graeme K Carnegie
- Department of Pharmacology and Department of Medicine, University of Illinois, Chicago, IL 60612
| | - Yong Zhao
- State Key Laboratory of Biomembrane and Membrane Biotechnology, Institute of Zoology, Chinese Academy of Sciences, Beijing 100101, China
| | - Asrar B Malik
- Department of Pharmacology and Department of Medicine, University of Illinois, Chicago, IL 60612
| | - Jingsong Xu
- Department of Pharmacology and Department of Medicine, University of Illinois, Chicago, IL 60612
| |
Collapse
|
22
|
Feng S, Zhu W. Bidirectional molecular transport shapes cell polarization in a two-dimensional model of eukaryotic chemotaxis. J Theor Biol 2014; 363:235-46. [DOI: 10.1016/j.jtbi.2014.08.033] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/17/2014] [Revised: 08/17/2014] [Accepted: 08/18/2014] [Indexed: 12/17/2022]
|
23
|
O'Neill PR, Giri L, Karunarathne WKA, Patel AK, Venkatesh KV, Gautam N. The structure of dynamic GPCR signaling networks. WILEY INTERDISCIPLINARY REVIEWS-SYSTEMS BIOLOGY AND MEDICINE 2014; 6:115-23. [PMID: 24741711 DOI: 10.1002/wsbm.1249] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
G-protein-coupled receptors (GPCRs) stimulate signaling networks that control a variety of critical physiological processes. Static information on the map of interacting signaling molecules at the basis of many cellular processes exists, but little is known about the dynamic operation of these networks. Here we focus on two questions. First, Is the network architecture underlying GPCR-activated cellular processes unique in comparison with others such as transcriptional networks? We discuss how spatially localized GPCR signaling requires uniquely organized networks to execute polarized cell responses. Second, What approaches overcome challenges in deciphering spatiotemporally dynamic networks that govern cell behavior? We focus on recently developed microfluidic and optical approaches that allow GPCR signaling pathways to be triggered and perturbed with spatially and temporally variant input while simultaneously visualizing molecular and cellular responses. When integrated with mathematical modeling, these approaches can help identify design principles that govern cell responses to extracellular signals. We outline why optical approaches that allow the behavior of a selected cell to be orchestrated continually are particularly well suited for probing network organization in single cells.
Collapse
|
24
|
Evolutionarily conserved coupling of adaptive and excitable networks mediates eukaryotic chemotaxis. Nat Commun 2014; 5:5175. [PMID: 25346418 PMCID: PMC4211273 DOI: 10.1038/ncomms6175] [Citation(s) in RCA: 58] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/30/2014] [Accepted: 09/06/2014] [Indexed: 01/22/2023] Open
Abstract
Numerous models explain how cells sense and migrate towards shallow chemoattractant gradients. Studies show that an excitable signal transduction network acts as a pacemaker that controls the cytoskeleton to drive motility. Here we show that this network is required to link stimuli to actin polymerization and chemotactic motility and we distinguish the various models of chemotaxis. First, signalling activity is suppressed towards the low side in a gradient or following removal of uniform chemoattractant. Second, signalling activities display a rapid shut off and a slower adaptation during which responsiveness to subsequent test stimuli decline. Simulations of various models indicate that these properties require coupled adaptive and excitable networks. Adaptation involves a G-protein-independent inhibitor, as stimulation of cells lacking G-protein function suppresses basal activities. The salient features of the coupled networks were observed for different chemoattractants in Dictyostelium and in human neutrophils, suggesting an evolutionarily conserved mechanism for eukaryotic chemotaxis.
Collapse
|
25
|
Byrne MB, Kimura Y, Kapoor A, He Y, Mattam KS, Hasan KM, Olson LN, Wang F, Kenis PJA, Rao CV. Oscillatory behavior of neutrophils under opposing chemoattractant gradients supports a winner-take-all mechanism. PLoS One 2014; 9:e85726. [PMID: 24465668 PMCID: PMC3897492 DOI: 10.1371/journal.pone.0085726] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/27/2013] [Accepted: 11/29/2013] [Indexed: 12/19/2022] Open
Abstract
Neutrophils constitute the largest class of white blood cells and are the first responders in the innate immune response. They are able to sense and migrate up concentration gradients of chemoattractants in search of primary sites of infection and inflammation through a process known as chemotaxis. These chemoattractants include formylated peptides and various chemokines. While much is known about chemotaxis to individual chemoattractants, far less is known about chemotaxis towards many. Previous studies have shown that in opposing gradients of intermediate chemoattractants (interleukin-8 and leukotriene B4), neutrophils preferentially migrate toward the more distant source. In this work, we investigated neutrophil chemotaxis in opposing gradients of chemoattractants using a microfluidic platform. We found that primary neutrophils exhibit oscillatory motion in opposing gradients of intermediate chemoattractants. To understand this behavior, we constructed a mathematical model of neutrophil chemotaxis. Our results suggest that sensory adaptation alone cannot explain the observed oscillatory motion. Rather, our model suggests that neutrophils employ a winner-take-all mechanism that enables them to transiently lock onto sensed targets and continuously switch between the intermediate attractant sources as they are encountered. These findings uncover a previously unseen behavior of neutrophils in opposing gradients of chemoattractants that will further aid in our understanding of neutrophil chemotaxis and the innate immune response. In addition, we propose a winner-take-all mechanism allows the cells to avoid stagnation near local chemical maxima when migrating through a network of chemoattractant sources.
Collapse
Affiliation(s)
- Matthew B. Byrne
- Department of Chemical and Biomolecular Engineering, University of Illinois at Urbana-Champaign, Urbana, Illinois, United States of America
| | - Yuki Kimura
- Department of Chemical and Biomolecular Engineering, University of Illinois at Urbana-Champaign, Urbana, Illinois, United States of America
| | - Ashish Kapoor
- Department of Chemical and Biomolecular Engineering, University of Illinois at Urbana-Champaign, Urbana, Illinois, United States of America
| | - Yuan He
- Department of Pathology, University of Michigan, Ann Arbor, Michigan, United States of America
| | - Kewin S. Mattam
- Department of Chemical and Biomolecular Engineering, University of Illinois at Urbana-Champaign, Urbana, Illinois, United States of America
| | - Katherine M. Hasan
- Department of Chemical and Biomolecular Engineering, University of Illinois at Urbana-Champaign, Urbana, Illinois, United States of America
| | - Luke N. Olson
- Department of Computer Science, University of Illinois at Urbana-Champaign, Urbana, Illinois, United States of America
| | - Fei Wang
- Department of Cell and Developmental Biology, University of Illinois at Urbana-Champaign, Urbana, Illinois, United States of America
| | - Paul J. A. Kenis
- Department of Chemical and Biomolecular Engineering, University of Illinois at Urbana-Champaign, Urbana, Illinois, United States of America
| | - Christopher V. Rao
- Department of Chemical and Biomolecular Engineering, University of Illinois at Urbana-Champaign, Urbana, Illinois, United States of America
- * E-mail:
| |
Collapse
|
26
|
Ma S, Cai C, Ma Y, Bai Z, Meng X, Yang X, Zou F, Ge R. Store-operated Ca²⁺ entry mediated regulation of polarization in differentiated human neutrophil-like HL-60 cells under hypoxia. Mol Med Rep 2014; 9:819-24. [PMID: 24425141 DOI: 10.3892/mmr.2014.1894] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/29/2013] [Accepted: 11/06/2013] [Indexed: 11/05/2022] Open
Abstract
The regulation of neutrophil polarization by calcium entry is critical for maintaining an effective host response. Hypoxia has a major effect on the apoptosis of neutrophils, however the role of store-operated Ca2+ entry (SOCE) in neutrophil polarization under hypoxia remains to be elucidated. In the present study, we examined the polarization of differentiated human neutrophil-like HL-60 (dHL-60) cells exposed to hypoxia (3% O2) and the results demonstrated that the percentage of polarized cells following exposure to an N-formyl-Met-Leu-Phe (fMLP) gradient in the Zigmond chamber was increased. We examined stromal interaction molecule 1 (STIM1) and Orai1 expression in dHL-60 cells during hypoxia, and it was observed that the expression of STIM1 and Orai1 was significantly reduced at day 2. However, no apparent change was observed on the first day, indicating that this effect is dependent on stimulation time. Fluo-4/acetoxymethyl (AM) ester imaging also demonstrated that SOCE was decreased in dHL-60 cells. The plasmid overexpression assay demonstrated that the response of polarization was returned to the control level. We demonstrated the inhibitory role of SOCE on the polarization of dHL-60 cells under hypoxic conditions, which may be the mechanism for the adaptation of neutrophils to hypoxia. SOCE is also suggested to be a key modulator of immune deficiency under hypoxic conditions and is potentially a therapeutic target.
Collapse
Affiliation(s)
- Shuang Ma
- Research Centre for High Altitude Medicine, Qinghai University Medical College, Qinghai University, Xining, Qinghai 810000, P.R. China
| | - Chunqing Cai
- Department of Occupational Health and Occupational Medicine, School of Public Health and Tropical Medicine, Southern Medical University, Guangzhou, Guangdong 510515, P.R. China
| | - Yan Ma
- Research Centre for High Altitude Medicine, Qinghai University Medical College, Qinghai University, Xining, Qinghai 810000, P.R. China
| | - Zhengzhong Bai
- Research Centre for High Altitude Medicine, Qinghai University Medical College, Qinghai University, Xining, Qinghai 810000, P.R. China
| | - Xiaojing Meng
- Department of Occupational Health and Occupational Medicine, School of Public Health and Tropical Medicine, Southern Medical University, Guangzhou, Guangdong 510515, P.R. China
| | - Xinyi Yang
- Department of Occupational Health and Occupational Medicine, School of Public Health and Tropical Medicine, Southern Medical University, Guangzhou, Guangdong 510515, P.R. China
| | - Fei Zou
- Department of Occupational Health and Occupational Medicine, School of Public Health and Tropical Medicine, Southern Medical University, Guangzhou, Guangdong 510515, P.R. China
| | - Rili Ge
- Research Centre for High Altitude Medicine, Qinghai University Medical College, Qinghai University, Xining, Qinghai 810000, P.R. China
| |
Collapse
|
27
|
Shi C, Huang CH, Devreotes PN, Iglesias PA. Interaction of motility, directional sensing, and polarity modules recreates the behaviors of chemotaxing cells. PLoS Comput Biol 2013; 9:e1003122. [PMID: 23861660 PMCID: PMC3701696 DOI: 10.1371/journal.pcbi.1003122] [Citation(s) in RCA: 77] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/07/2013] [Accepted: 05/16/2013] [Indexed: 02/04/2023] Open
Abstract
Chemotaxis involves the coordinated action of separable but interrelated processes: motility, gradient sensing, and polarization. We have hypothesized that these are mediated by separate modules that account for these processes individually and that, when combined, recreate most of the behaviors of chemotactic cells. Here, we describe a mathematical model where the modules are implemented in terms of reaction-diffusion equations. Migration and the accompanying changes in cellular morphology are demonstrated in simulations using a mechanical model of the cell cortex implemented in the level set framework. The central module is an excitable network that accounts for random migration. The response to combinations of uniform stimuli and gradients is mediated by a local excitation, global inhibition module that biases the direction in which excitability is directed. A polarization module linked to the excitable network through the cytoskeleton allows unstimulated cells to move persistently and, for cells in gradients, to gradually acquire distinct sensitivity between front and back. Finally, by varying the strengths of various feedback loops in the model we obtain cellular behaviors that mirror those of genetically altered cell lines. Chemotaxis is the movement of cells in response to spatial gradients of chemical cues. While single-celled organisms rely on sensing and responding to chemical gradients to search for nutrients, chemotaxis is also an essential component of the mammalian immune system. However, chemotaxis can also be deleterious, since chemotactic tumor cells can lead to metastasis. Due to its importance, understanding the process by which cells sense and respond to chemical gradients has attracted considerable interest. Moreover, because of the complexity of chemotactic signaling, which includes multiple feedback loops and redundant pathways, this has been a research area in which computational models have had a significant impact in understanding experimental findings. Here, we propose a modular description of the signaling network that regulates chemotaxis. The modules describe different processes that are observed in chemotactic cells. In addition to accounting for these behaviors individually, we show that the overall system recreates many features of the directed motion of migrating cells. The signaling described by our modules is implemented as a series of equations, whereas movement and the accompanying cellular deformations are simulated using a mechanical model of the cell and implemented using level set methods, a method that allows simulations of cells as they change morphology.
Collapse
Affiliation(s)
- Changji Shi
- Department of Electrical and Computer Engineering, Whiting School of Engineering, Johns Hopkins University, Baltimore, Maryland, United States of America
| | - Chuan-Hsiang Huang
- Department of Cell Biology, School of Medicine, Johns Hopkins University, Baltimore, Maryland, United States of America
| | - Peter N. Devreotes
- Department of Cell Biology, School of Medicine, Johns Hopkins University, Baltimore, Maryland, United States of America
| | - Pablo A. Iglesias
- Department of Electrical and Computer Engineering, Whiting School of Engineering, Johns Hopkins University, Baltimore, Maryland, United States of America
- Biological Physics, Max Planck Institute for the Physics of Complex Systems, Dresden, Germany
- * E-mail:
| |
Collapse
|
28
|
|
29
|
Wang Y, Ku CJ, Zhang ER, Artyukhin AB, Weiner OD, Wu LF, Altschuler SJ. Identifying network motifs that buffer front-to-back signaling in polarized neutrophils. Cell Rep 2013; 3:1607-16. [PMID: 23665220 DOI: 10.1016/j.celrep.2013.04.009] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/10/2012] [Revised: 03/05/2013] [Accepted: 04/09/2013] [Indexed: 01/02/2023] Open
Abstract
Neutrophil polarity relies on local, mutual inhibition to segregate incompatible signaling circuits to the leading and trailing edges. Mutual inhibition alone should lead to cells having strong fronts and weak backs or vice versa. However, analysis of cell-to-cell variation in human neutrophils revealed that back polarity remains consistent despite changes in front strength. How is this buffering achieved? Pharmacological perturbations and mathematical modeling revealed a functional role for microtubules in buffering back polarity by mediating positive, long-range crosstalk from front to back; loss of microtubules inhibits buffering and results in anticorrelation between front and back signaling. Furthermore, a systematic, computational search of network topologies found that a long-range, positive front-to-back link is necessary for back buffering. Our studies suggest a design principle that can be employed by polarity networks: short-range mutual inhibition establishes distinct signaling regions, after which directed long-range activation insulates one region from variations in the other.
Collapse
Affiliation(s)
- Yanqin Wang
- Green Center for Systems Biology, Department of Pharmacology, Simmons Cancer Center, University of Texas Southwestern Medical Center, Dallas, TX 75390, USA
| | | | | | | | | | | | | |
Collapse
|
30
|
Wendelsdorf KV, Alam M, Bassaganya-Riera J, Bisset K, Eubank S, Hontecillas R, Hoops S, Marathe M. ENteric Immunity SImulator: a tool for in silico study of gastroenteric infections. IEEE Trans Nanobioscience 2013; 11:273-88. [PMID: 22987134 DOI: 10.1109/tnb.2012.2211891] [Citation(s) in RCA: 33] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023]
Abstract
Clinical symptoms of microbial infection of the gastrointestinal (GI) tract are often exacerbated by inflammation induced pathology. Identifying novel avenues for treating and preventing such pathologies is necessary and complicated by the complexity of interacting immune pathways in the gut, where effector and inflammatory immune cells are regulated by anti-inflammatory or regulatory cells. Here we present new advances in the development of the ENteric Immunity SImulator (ENISI), a simulator of GI immune mechanisms in response to resident commensal bacteria as well as invading pathogens and the effect on the development of intestinal lesions. ENISI is a tool for identifying potential treatment strategies that reduce inflammation-induced damage and, at the same time, ensure pathogen removal by allowing one to test plausibility of in vitro observed behavior as explanations for observations in vivo, propose behaviors not yet tested in vitro that could explain these tissue-level observations, and conduct low-cost, preliminary experiments of proposed interventions/treatments. An example of such application is shown in which we simulate dysentery resulting from Brachyispira hyodysenteriae infection and identify aspects of the host immune pathways that lead to continued inflammation-induced tissue damage even after pathogen elimination.
Collapse
Affiliation(s)
- Katherine V Wendelsdorf
- Network Dynamics and Simulation Science Laboratory, and Center for Modeling Immunity to Enteric Pathogens, Virginia Bioinformatics Institute, Virginia Polytechnic Institute and State University, Blacksburg, VA 24061, USA.
| | | | | | | | | | | | | | | |
Collapse
|
31
|
Chau AH, Walter JM, Gerardin J, Tang C, Lim WA. Designing synthetic regulatory networks capable of self-organizing cell polarization. Cell 2012; 151:320-32. [PMID: 23039994 DOI: 10.1016/j.cell.2012.08.040] [Citation(s) in RCA: 128] [Impact Index Per Article: 10.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/08/2012] [Revised: 07/07/2012] [Accepted: 08/16/2012] [Indexed: 11/26/2022]
Abstract
How cells form global, self-organized structures using genetically encoded molecular rules remains elusive. Here, we take a synthetic biology approach to investigate the design principles governing cell polarization. First, using a coarse-grained computational model, we searched for all possible simple networks that can achieve polarization. All solutions contained one of three minimal motifs: positive feedback, mutual inhibition, or inhibitor with positive feedback. These minimal motifs alone could achieve polarization under limited conditions; circuits that combined two or more of these motifs were significantly more robust. With these design principles as a blueprint, we experimentally constructed artificial polarization networks in yeast, using a toolkit of chimeric signaling proteins that spatially direct the synthesis and degradation of phosphatidylinositol (3,4,5)-trisphosphate (PIP(3)). Circuits with combinatorial motifs yielded clear foci of synthetic PIP(3) that can persist for nearly an hour. Thus, by harnessing localization-regulated signaling molecules, we can engineer simple molecular circuits that reliably execute spatial self-organized programs.
Collapse
Affiliation(s)
- Angela H Chau
- Howard Hughes Medical Institute and Department of Cellular and Molecular Pharmacology, University of California, San Francisco, San Francisco, CA 94158, USA
| | | | | | | | | |
Collapse
|
32
|
Holmes GR, Anderson SR, Dixon G, Robertson AL, Reyes-Aldasoro CC, Billings SA, Renshaw SA, Kadirkamanathan V. Repelled from the wound, or randomly dispersed? Reverse migration behaviour of neutrophils characterized by dynamic modelling. J R Soc Interface 2012; 9:3229-39. [PMID: 22951343 DOI: 10.1098/rsif.2012.0542] [Citation(s) in RCA: 32] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023] Open
Abstract
Following neutralization of infectious threats, neutrophils must be removed from inflammatory sites for normal tissue function to be restored. Recently, a new paradigm has emerged, in which viable neutrophils migrate away from inflammatory sites by a process best described as reverse migration. It has generally been assumed that this process is the mirror image of chemotaxis, where neutrophils are drawn into the areas of infection or tissue damage by gradients of chemotactic cues. Indeed, efforts are underway to identify cues that drive neutrophils away by the reverse process, fugetaxis. By using photoconvertible pigments expressed in neutrophils in transparent zebrafish larvae, we were able to image the position of each neutrophil during inflammation resolution in vivo. These neutrophil coordinates were analysed within a dynamic modelling framework, using different forms of the drift-diffusion equation with model selection and parameter estimation based on approximate Bayesian computation. This analysis revealed the experimental data were best fitted by a model incorporating a diffusion term but no drift term-where the presence of drift would indicate fugetaxis. This result, for the first time, provides rigorous data-driven evidence that reverse migration of neutrophils in vivo is not a form of fugetaxis, but rather a stochastic redistribution.
Collapse
Affiliation(s)
- Geoffrey R Holmes
- Complex Systems and Signal Processing Group, Department of Automatic Control and Systems Engineering, University of Sheffield, Sheffield, UK
| | | | | | | | | | | | | | | |
Collapse
|
33
|
Nishimura SI, Ueda M, Sasai M. Non-Brownian dynamics and strategy of amoeboid cell locomotion. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2012; 85:041909. [PMID: 22680500 DOI: 10.1103/physreve.85.041909] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/26/2011] [Indexed: 06/01/2023]
Abstract
Amoeboid cells such as Dictyostelium discoideum and Madin-Darby canine kidney cells show the non-Brownian dynamics of migration characterized by the superdiffusive increase of mean-squared displacement. In order to elucidate the physical mechanism of this non-Brownian dynamics, a computational model is developed which highlights a group of inhibitory molecules for actin polymerization. Based on this model, we propose a hypothesis that inhibitory molecules are sent backward in the moving cell to accumulate at the rear of cell. The accumulated inhibitory molecules at the rear further promote cell locomotion to form a slow positive feedback loop of the whole-cell scale. The persistent straightforward migration is stabilized with this feedback mechanism, but the fluctuation in the distribution of inhibitory molecules and the cell shape deformation concurrently interrupt the persistent motion to turn the cell into a new direction. A sequence of switching behaviors between persistent motions and random turns gives rise to the superdiffusive migration in the absence of the external guidance signal. In the complex environment with obstacles, this combined process of persistent motions and random turns drives the simulated amoebae to solve the maze problem in a highly efficient way, which suggests the biological advantage for cells to bear the non-Brownian dynamics.
Collapse
Affiliation(s)
- Shin I Nishimura
- Department of Mathematical and Life Sciences, Hiroshima University, Hiroshima 739-8526, Japan.
| | | | | |
Collapse
|
34
|
Liepe J, Taylor H, Barnes CP, Huvet M, Bugeon L, Thorne T, Lamb JR, Dallman MJ, Stumpf MPH. Calibrating spatio-temporal models of leukocyte dynamics against in vivo live-imaging data using approximate Bayesian computation. Integr Biol (Camb) 2012; 4:335-345. [PMID: 22327539 PMCID: PMC5058438 DOI: 10.1039/c2ib00175f] [Citation(s) in RCA: 30] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
In vivo studies allow us to investigate biological processes at the level of the organism. But not all aspects of in vivo systems are amenable to direct experimental measurements. In order to make the most of such data we therefore require statistical tools that allow us to obtain reliable estimates for e.g. kinetic in vivo parameters. Here we show how we can use approximate Bayesian computation approaches in order to analyse leukocyte migration in zebrafish embryos in response to injuries. We track individual leukocytes using live imaging following surgical injury to the embryos' tail-fins. The signalling gradient that leukocytes follow towards the site of the injury cannot be directly measured but we can estimate its shape and how it changes with time from the directly observed patterns of leukocyte migration. By coupling simple models of immune signalling and leukocyte migration with the unknown gradient shape into a single statistical framework we can gain detailed insights into the tissue-wide processes that are involved in the innate immune response to wound injury. In particular we find conclusive evidence for a temporally and spatially changing signalling gradient that modulates the changing activity of the leukocyte population in the embryos. We conclude with a robustness analysis which highlights the most important factors determining the leukocyte dynamics. Our approach relies only on the ability to simulate numerically the process under investigation and is therefore also applicable in other in vivo contexts and studies.
Collapse
Affiliation(s)
- Juliane Liepe
- Centre for Integrative Systems Biology and Bioinformatics, Division of Molecular Biosciences, Department of Life Sciences, Imperial College London, London, UK
| | - Harriet Taylor
- Division of Cell and Molecular Biology, Department of Life Sciences, Imperial College London, London, UK
- MRC Centre for Inflammation Research, Queens Medical Research Institute, College of Medicine and Veterinary Medicine, University of Edinburgh, Edinburgh, UK
| | - Chris P. Barnes
- Centre for Integrative Systems Biology and Bioinformatics, Division of Molecular Biosciences, Department of Life Sciences, Imperial College London, London, UK
| | - Maxime Huvet
- Centre for Integrative Systems Biology and Bioinformatics, Division of Molecular Biosciences, Department of Life Sciences, Imperial College London, London, UK
| | - Laurence Bugeon
- Division of Cell and Molecular Biology, Department of Life Sciences, Imperial College London, London, UK
| | - Thomas Thorne
- Centre for Integrative Systems Biology and Bioinformatics, Division of Molecular Biosciences, Department of Life Sciences, Imperial College London, London, UK
| | - Jonathan R. Lamb
- Division of Cell and Molecular Biology, Department of Life Sciences, Imperial College London, London, UK
| | - Margaret J. Dallman
- Division of Cell and Molecular Biology, Department of Life Sciences, Imperial College London, London, UK
- Centre for Integrative Systems Biology, Department of Life Sciences, Imperial College London, London, UK
| | - Michael P. H. Stumpf
- Centre for Integrative Systems Biology and Bioinformatics, Division of Molecular Biosciences, Department of Life Sciences, Imperial College London, London, UK
- Centre for Integrative Systems Biology, Department of Life Sciences, Imperial College London, London, UK
- Institute of Mathematical Sciences, Imperial College London, London, UK
| |
Collapse
|
35
|
Semplice M, Veglio A, Naldi G, Serini G, Gamba A. A bistable model of cell polarity. PLoS One 2012; 7:e30977. [PMID: 22383986 PMCID: PMC3285628 DOI: 10.1371/journal.pone.0030977] [Citation(s) in RCA: 32] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/12/2011] [Accepted: 12/29/2011] [Indexed: 12/16/2022] Open
Abstract
Ultrasensitivity, as described by Goldbeter and Koshland, has been considered for a long time as a way to realize bistable switches in biological systems. It is not as well recognized that when ultrasensitivity and reinforcing feedback loops are present in a spatially distributed system such as the cell plasmamembrane, they may induce bistability and spatial separation of the system into distinct signaling phases. Here we suggest that bistability of ultrasensitive signaling pathways in a diffusive environment provides a basic mechanism to realize cell membrane polarity. Cell membrane polarization is a fundamental process implicated in several basic biological phenomena, such as differentiation, proliferation, migration and morphogenesis of unicellular and multicellular organisms. We describe a simple, solvable model of cell membrane polarization based on the coupling of membrane diffusion with bistable enzymatic dynamics. The model can reproduce a broad range of symmetry-breaking events, such as those observed in eukaryotic directional sensing, the apico-basal polarization of epithelium cells, the polarization of budding and mating yeast, and the formation of Ras nanoclusters in several cell types.
Collapse
Affiliation(s)
- Matteo Semplice
- Department of Physics and Mathematics, Università dell'Insubria, Como, Italy
| | - Andrea Veglio
- Genomes and Genetics Department, Unit Physics of Biological Systems, Institut Pasteur, Paris, France
| | - Giovanni Naldi
- Department of Mathematics “F. Enriques”, Università degli studi di Milano, Milano, Italy
- * E-mail:
| | - Guido Serini
- Laboratory of Cell Adhesion Dynamics, Institute for Cancer Research and Treatment and Department of Oncological Sciences, School of Medicine, Università degli studi di Torino, Candiolo, Italy
| | - Andrea Gamba
- Department of Mathematics, Politecnico di Torino, Torino, Italy
- Laboratory of Systems Biology, Institute for Cancer Research and Treatment, Candiolo, Italy
- INFN, Torino, Italy
| |
Collapse
|
36
|
|
37
|
Abstract
Phosphatidylinositol lipids generated through the action of phosphinositide 3-kinase (PI3K) are key mediators of a wide array of biological responses. In particular, their role in the regulation of cell migration has been extensively studied and extends to amoeboid as well as mesenchymal migration. Through the emergence of fluorescent probes that target PI3K products as well as the use of specific inhibitors and knockout technologies, the spatio-temporal distribution of PI3K products in chemotaxing cells has been shown to represent a key anterior polarity signal that targets downstream effectors to actin polymerization. In addition, through intricate cross-talk networks PI3K products have been shown to regulate signals that control posterior effectors. Yet, in more complex environments or in conditions where chemoattractant gradients are steep, a variety of cell types can still chemotax in the absence of PI3K signals. Indeed, parallel signal transduction pathways have been shown to coordinately regulate cell polarity and directed movement. In this chapter, we will review the current role PI3K products play in the regulation of directed cell migration in various cell types, highlight the importance of mathematical modeling in the study of chemotaxis, and end with a brief overview of other signaling cascades known to also regulate chemotaxis.
Collapse
Affiliation(s)
- Michael C Weiger
- Laboratory of Cellular and Molecular Biology, National Cancer Institute, National Institutes of Health, 37 Convent Drive, Bldg.37/Rm2066, 20892-4256, Bethesda, MD, USA
| | | |
Collapse
|
38
|
Zheng Z, Chou CS, Yi TM, Nie Q. Mathematical analysis of steady-state solutions in compartment and continuum models of cell polarization. MATHEMATICAL BIOSCIENCES AND ENGINEERING : MBE 2011; 8:1135-1168. [PMID: 21936604 PMCID: PMC3806509 DOI: 10.3934/mbe.2011.8.1135] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/31/2023]
Abstract
Cell polarization, in which substances previously uniformly distributed become asymmetric due to external or/and internal stimulation, is a fundamental process underlying cell mobility, cell division, and other polarized functions. The yeast cell S. cerevisiae has been a model system to study cell polarization. During mating, yeast cells sense shallow external spatial gradients and respond by creating steeper internal gradients of protein aligned with the external cue. The complex spatial dynamics during yeast mating polarization consists of positive feedback, degradation, global negative feedback control, and cooperative effects in protein synthesis. Understanding such complex regulations and interactions is critical to studying many important characteristics in cell polarization including signal amplification, tracking dynamic signals, and potential trade-off between achieving both objectives in a robust fashion. In this paper, we study some of these questions by analyzing several models with different spatial complexity: two compartments, three compartments, and continuum in space. The step-wise approach allows detailed characterization of properties of the steady state of the system, providing more insights for biological regulations during cell polarization. For cases without membrane diffusion, our study reveals that increasing the number of spatial compartments results in an increase in the number of steady-state solutions, in particular, the number of stable steady-state solutions, with the continuum models possessing infinitely many steady-state solutions. Through both analysis and simulations, we find that stronger positive feedback, reduced diffusion, and a shallower ligand gradient all result in more steady-state solutions, although most of these are not optimally aligned with the gradient. We explore in the different settings the relationship between the number of steady-state solutions and the extent and accuracy of the polarization. Taken together these results furnish a detailed description of the factors that influence the tradeoff between a single correctly aligned but poorly polarized stable steady-state solution versus multiple more highly polarized stable steady-state solutions that may be incorrectly aligned with the external gradient.
Collapse
Affiliation(s)
- Zhenzhen Zheng
- Department of Mathematics, Center for Complex Biological Systems and Center for Mathematical and Computational Biology, University of California-Irvine, CA 92697, United States.
| | | | | | | |
Collapse
|
39
|
Dawes AT, Munro EM. PAR-3 oligomerization may provide an actin-independent mechanism to maintain distinct par protein domains in the early Caenorhabditis elegans embryo. Biophys J 2011; 101:1412-22. [PMID: 21943422 DOI: 10.1016/j.bpj.2011.07.030] [Citation(s) in RCA: 48] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/11/2011] [Revised: 06/11/2011] [Accepted: 07/06/2011] [Indexed: 02/06/2023] Open
Abstract
Par proteins establish discrete intracellular spatial domains to polarize many different cell types. In the single-cell embryo of the nematode worm Caenorhabditis elegans, the segregation of Par proteins is crucial for proper division and cell fate specification. Actomyosin-based cortical flows drive the initial formation of anterior and posterior Par domains, but cortical actin is not required for the maintenance of these domains. Here we develop a model of interactions between the Par proteins that includes both mutual inhibition and PAR-3 oligomerization. We show that this model gives rise to a bistable switch mechanism, allowing the Par proteins to occupy distinct anterior and posterior domains seen in the early C. elegans embryo, independent of dynamics or asymmetries in the actin cortex. The model predicts a sharp loss of cortical Par protein asymmetries during gradual depletion of the Par protein PAR-6, and we confirm this prediction experimentally. Together, these results suggest both mutual inhibition and PAR-3 oligomerization are sufficient to maintain distinct Par protein domains in the early C. elegans embryo.
Collapse
Affiliation(s)
- Adriana T Dawes
- Department of Mathematical and Statistical Statistical Sciences, University of Alberta, Edmonton, Alberta, Canada.
| | | |
Collapse
|
40
|
Welf ES, Haugh JM. Signaling pathways that control cell migration: models and analysis. WILEY INTERDISCIPLINARY REVIEWS-SYSTEMS BIOLOGY AND MEDICINE 2011; 3:231-40. [PMID: 21305705 DOI: 10.1002/wsbm.110] [Citation(s) in RCA: 51] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/16/2023]
Abstract
Dissecting the intracellular signaling mechanisms that govern the movement of eukaryotic cells presents a major challenge, not only because of the large number of molecular players involved, but even more so because of the dynamic nature of their regulation by both biochemical and mechanical interactions. Computational modeling and analysis have emerged as useful tools for understanding how the physical properties of cells and their microenvironment are coupled with certain biochemical pathways to actuate and control cell motility. In this focused review, we highlight some of the more recent applications of quantitative modeling and analysis in the field of cell migration. Both in modeling and experiment, it has been prudent to follow a reductionist approach in order to characterize what are arguably the principal modules: spatial polarization of signaling pathways, regulation of the actin cytoskeleton, and dynamics of focal adhesions. While it is important that we 'cut our teeth' on these subsystems, focusing on the details of certain aspects while ignoring or coarse-graining others, it is clear that the challenge ahead will be to characterize the couplings between them in an integrated framework.
Collapse
Affiliation(s)
- Erik S Welf
- Department of Chemical and Biomolecular Engineering, North Carolina State University, Raleigh, NC, USA
| | | |
Collapse
|
41
|
Jilkine A, Edelstein-Keshet L. A comparison of mathematical models for polarization of single eukaryotic cells in response to guided cues. PLoS Comput Biol 2011; 7:e1001121. [PMID: 21552548 PMCID: PMC3084230 DOI: 10.1371/journal.pcbi.1001121] [Citation(s) in RCA: 171] [Impact Index Per Article: 13.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022] Open
Abstract
Polarization, a primary step in the response of an individual eukaryotic cell to a spatial stimulus, has attracted numerous theoretical treatments complementing experimental studies in a variety of cell types. While the phenomenon itself is universal, details differ across cell types, and across classes of models that have been proposed. Most models address how symmetry breaking leads to polarization, some in abstract settings, others based on specific biochemistry. Here, we compare polarization in response to a stimulus (e.g., a chemoattractant) in cells typically used in experiments (yeast, amoebae, leukocytes, keratocytes, fibroblasts, and neurons), and, in parallel, responses of several prototypical models to typical stimulation protocols. We find that the diversity of cell behaviors is reflected by a diversity of models, and that some, but not all models, can account for amplification of stimulus, maintenance of polarity, adaptation, sensitivity to new signals, and robustness.
Collapse
Affiliation(s)
- Alexandra Jilkine
- Green Comprehensive Center for Computational and Systems Biology, Department of Pharmacology, University of Texas Southwestern Medical Center, Dallas, Texas, United States of America.
| | | |
Collapse
|
42
|
Iglesias PA. Spatial regulation of PI3K signaling during chemotaxis. WILEY INTERDISCIPLINARY REVIEWS-SYSTEMS BIOLOGY AND MEDICINE 2011; 1:247-253. [PMID: 20835994 DOI: 10.1002/wsbm.13] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/05/2022]
Abstract
Phosphoinositide 3-kinases (PI3Ks) are a family of lipid kinases that phosphorylate the 3' OH position of the inositol ring of phosphoinositides on the inner leaf of the plasma membrane. Receptor-mediated activation of the PI3K pathway plays a crucial role in numerous signaling pathways and regulates a number of critical cellular processes, including growth, differentiation, survival and directed migration. In this focus article, we review the temporal and spatial regulation of PI3K in chemotaxing cells with particular emphasis on the amoeba Dictyostelium as well as neutrophils. We also briefly discuss one model used to elucidate the PI3K pathway.
Collapse
Affiliation(s)
- Pablo A Iglesias
- Department of Electrical and Computer Engineering, The Johns Hopkins University, 3400 N. Charles St., Baltimore, MD 21218, USA
| |
Collapse
|
43
|
Cheong R, Paliwal S, Levchenko A. Models at the single cell level. WILEY INTERDISCIPLINARY REVIEWS-SYSTEMS BIOLOGY AND MEDICINE 2011; 2:34-48. [PMID: 20836009 DOI: 10.1002/wsbm.49] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
Many cellular behaviors cannot be completely captured or appropriately described at the cell population level. Noise induced by stochastic chemical reactions, spatially polarized signaling networks, and heterogeneous cell-cell communication are among the many phenomena that require fine-grained analysis. Accordingly, the mathematical models used to describe such systems must be capable of single cell or subcellular resolution. Here, we review techniques for modeling single cells, including models of stochastic chemical kinetics, spatially heterogeneous intracellular signaling, and spatial stochastic systems. We also briefly discuss applications of each type of model.
Collapse
Affiliation(s)
- Raymond Cheong
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD 21218, USA
| | - Saurabh Paliwal
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD 21218, USA
| | - Andre Levchenko
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD 21218, USA.,Whitaker Institute of Biomedical Engineering and Institute for Cell Engineering, Johns Hopkins University, Baltimore, MD, USA
| |
Collapse
|
44
|
Otsuji M, Terashima Y, Ishihara S, Kuroda S, Matsushima K. A conceptual molecular network for chemotactic behaviors characterized by feedback of molecules cycling between the membrane and the cytosol. Sci Signal 2010; 3:ra89. [PMID: 21156936 DOI: 10.1126/scisignal.2001056] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/02/2022]
Abstract
Cell chemotaxis has been characterized as the formation of a front-back axis that is triggered by a gradient of chemoattractant; however, chemotaxis is accompanied by more complicated behaviors. These include migration in a straight line with a stable axis [the stable single-axis (SSA) pattern] and repeated splitting of the leading edge of the cell into two regions, followed by the "choice" of one of these as the new leading edge [the split and choice (S&C) pattern]. Indeed, transition between these two behaviors can be observed in individual cells. However, the conceptual framework of the network of signaling molecules that generates these patterns remains to be clarified. We confirmed theoretically that a system that has positive and negative feedback loops involving the reciprocal cycling between the membrane and the cytosol of molecules that promote membrane protrusion or retraction generates SSA and S&C patterns of migratory behavior under similar conditions. We also predicted properties of the instabilities of such a system, which are essential for the generation of these behaviors, and we verified their existence in chemotaxing cells. Our research provides a simple model of network structure for chemotactic behaviors, including cell polarization.
Collapse
Affiliation(s)
- Mikiya Otsuji
- Department of Molecular Preventive Medicine, Graduate School of Medicine, University of Tokyo, Tokyo 113-0033, Japan.
| | | | | | | | | |
Collapse
|
45
|
Abstract
Chemotaxis plays a crucial role in many biological processes, including nervous system development. However, fundamental physical constraints limit the ability of a small sensing device such as a cell or growth cone to detect an external chemical gradient. One of these is the stochastic nature of receptor binding, leading to a constantly fluctuating binding pattern across the cell's array of receptors. This is analogous to the uncertainty in sensory information often encountered by the brain at the systems level. Here we derive analytically the Bayes-optimal strategy for combining information from a spatial array of receptors in both one and two dimensions to determine gradient direction. We also show how information from more than one receptor species can be optimally integrated, derive the gradient shapes that are optimal for guiding cells or growth cones over the longest possible distances, and illustrate that polarized cell behavior might arise as an adaptation to slowly varying environments. Together our results provide closed-form predictions for variations in chemotactic performance over a wide range of gradient conditions.
Collapse
Affiliation(s)
- Duncan Mortimer
- Queensland Brain Institute, University of Queensland, St. Lucia QLD 4072, Australia.
| | | | | | | |
Collapse
|
46
|
Dehmelt L, Bastiaens PIH. Spatial organization of intracellular communication: insights from imaging. Nat Rev Mol Cell Biol 2010; 11:440-52. [DOI: 10.1038/nrm2903] [Citation(s) in RCA: 111] [Impact Index Per Article: 7.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023]
|
47
|
Weiger MC, Ahmed S, Welf ES, Haugh JM. Directional persistence of cell migration coincides with stability of asymmetric intracellular signaling. Biophys J 2010; 98:67-75. [PMID: 20085720 DOI: 10.1016/j.bpj.2009.09.051] [Citation(s) in RCA: 41] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/27/2009] [Revised: 09/05/2009] [Accepted: 09/25/2009] [Indexed: 12/23/2022] Open
Abstract
It has long been appreciated that spatiotemporal dynamics of cell migration are under the control of intracellular signaling pathways, which are mediated by adhesion receptors and other transducers of extracellular cues. Further, there is ample evidence that aspects of cell migration are stochastic: how else could it exhibit directional persistence over timescales much longer than typical signal transduction processes, punctuated by abrupt changes in direction? Yet the mechanisms by which signaling processes affect those behaviors remain unclear. We have developed analytical methods for relating parallel live-cell microscopy measurements of cell migration dynamics to the intracellular signaling processes that govern them. In this analysis of phosphoinositide 3-kinase signaling in randomly migrating fibroblasts, we observe that hot spots of intense signaling coincide with localized cell protrusion and endure with characteristic lifetimes that correspond to those of cell migration persistence. We further show that distant hot spots are dynamically and stochastically coupled. These results are indicative of a mechanism by which changes in a cell's direction of migration are determined by a fragile balance of relatively rapid intracellular signaling processes.
Collapse
Affiliation(s)
- Michael C Weiger
- Department of Chemical and Biomolecular Engineering, North Carolina State University, Raleigh, North Carolina, USA
| | | | | | | |
Collapse
|
48
|
Hawkins RJ, Bénichou O, Piel M, Voituriez R. Rebuilding cytoskeleton roads: active-transport-induced polarization of cells. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2009; 80:040903. [PMID: 19905266 DOI: 10.1103/physreve.80.040903] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/09/2009] [Indexed: 05/28/2023]
Abstract
Many cellular processes require a polarization axis which generally initially emerges as an inhomogeneous distribution of molecular markers in the cell. We present a simple analytical model of a general mechanism of cell polarization taking into account the positive feedback due to the coupled dynamics of molecular markers and cytoskeleton filaments. We find that the geometry of the organization of cytoskeleton filaments, nucleated on the membrane (e.g., cortical actin) or from a center in the cytoplasm (e.g., microtubule asters), dictates whether the system is capable of spontaneous polarization or polarizes only in response to external asymmetric signals. Our model also captures the main features of recent experiments of cell polarization in two considerably different biological systems, namely, mating budding yeast and neuron growth cones.
Collapse
Affiliation(s)
- R J Hawkins
- UMR 7600, Université Pierre et Marie Curie/CNRS, 4 Place Jussieu, 75255 Paris Cedex 05, France
| | | | | | | |
Collapse
|
49
|
Endres RG, Wingreen NS. Accuracy of direct gradient sensing by cell-surface receptors. PROGRESS IN BIOPHYSICS AND MOLECULAR BIOLOGY 2009; 100:33-9. [DOI: 10.1016/j.pbiomolbio.2009.06.002] [Citation(s) in RCA: 40] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/01/2022]
|
50
|
Causin P, Facchetti G. Autocatalytic loop, amplification and diffusion: a mathematical and computational model of cell polarization in neural chemotaxis. PLoS Comput Biol 2009; 5:e1000479. [PMID: 19714204 PMCID: PMC2722090 DOI: 10.1371/journal.pcbi.1000479] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/07/2008] [Accepted: 07/21/2009] [Indexed: 12/11/2022] Open
Abstract
The chemotactic response of cells to graded fields of chemical cues is a complex process that requires the coordination of several intracellular activities. Fundamental steps to obtain a front vs. back differentiation in the cell are the localized distribution of internal molecules and the amplification of the external signal. The goal of this work is to develop a mathematical and computational model for the quantitative study of such phenomena in the context of axon chemotactic pathfinding in neural development. In order to perform turning decisions, axons develop front-back polarization in their distal structure, the growth cone. Starting from the recent experimental findings of the biased redistribution of receptors on the growth cone membrane, driven by the interaction with the cytoskeleton, we propose a model to investigate the significance of this process. Our main contribution is to quantitatively demonstrate that the autocatalytic loop involving receptors, cytoplasmic species and cytoskeleton is adequate to give rise to the chemotactic behavior of neural cells. We assess the fact that spatial bias in receptors is a precursory key event for chemotactic response, establishing the necessity of a tight link between upstream gradient sensing and downstream cytoskeleton dynamics. We analyze further crosslinked effects and, among others, the contribution to polarization of internal enzymatic reactions, which entail the production of molecules with a one-to-more factor. The model shows that the enzymatic efficiency of such reactions must overcome a threshold in order to give rise to a sufficient amplification, another fundamental precursory step for obtaining polarization. Eventually, we address the characteristic behavior of the attraction/repulsion of axons subjected to the same cue, providing a quantitative indicator of the parameters which more critically determine this nontrivial chemotactic response.
Collapse
Affiliation(s)
- Paola Causin
- Department of Mathematics F Enriques, Università degli Studi di Milano, Milano, Italy.
| | | |
Collapse
|