1
|
Hladyshau S, Guan K, Nivedita N, Errede B, Tsygankov D, Elston TC. Multiscale Modeling of Bistability in the Yeast Polarity Circuit. Cells 2024; 13:1358. [PMID: 39195248 DOI: 10.3390/cells13161358] [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/03/2024] [Revised: 08/05/2024] [Accepted: 08/13/2024] [Indexed: 08/29/2024] Open
Abstract
Cell polarity refers to the asymmetric distribution of proteins and other molecules along a specified axis within a cell. Polarity establishment is the first step in many cellular processes. For example, directed growth or migration requires the formation of a cell front and back. In many cases, polarity occurs in the absence of spatial cues. That is, the cell undergoes symmetry breaking. Understanding the molecular mechanisms that allow cells to break symmetry and polarize requires computational models that span multiple spatial and temporal scales. Here, we apply a multiscale modeling approach to examine the polarity circuit of yeast. In addition to symmetry breaking, experiments revealed two key features of the yeast polarity circuit: bistability and rapid dismantling of the polarity site following a loss of signal. We used modeling based on ordinary differential equations (ODEs) to investigate mechanisms that generate these behaviors. Our analysis revealed that a model involving positive and negative feedback acting on different time scales captured both features. We then extend our ODE model into a coarse-grained reaction-diffusion equation (RDE) model to capture the spatial profiles of polarity factors. After establishing that the coarse-grained RDE model qualitatively captures key features of the polarity circuit, we expand it to more accurately capture the biochemical reactions involved in the system. We convert the expanded model to a particle-based model that resolves individual molecules and captures fluctuations that arise from the stochastic nature of biochemical reactions. Our models assume that negative regulation results from negative feedback. However, experimental observations do not rule out the possibility that negative regulation occurs through an incoherent feedforward loop. Therefore, we conclude by using our RDE model to suggest how negative feedback might be distinguished from incoherent feedforward regulation.
Collapse
Affiliation(s)
- Siarhei Hladyshau
- School of Biological Sciences, Georgia Institute of Technology, Atlanta, GA 30332, USA
- Wallace H. Coulter Department of Biomedical Engineering, Georgia Institute of Technology and Emory University, Atlanta, GA 30332, USA
| | - Kaiyun Guan
- Curriculum in Bioinformatics and Computational Biology, University of North Carolina at Chapel Hill, Chapel Hill, NC 27599, USA
| | - Nivedita Nivedita
- Department of Pharmacology, University of North Carolina at Chapel Hill, Chapel Hill, NC 27599, USA
| | - Beverly Errede
- Department of Biochemistry and Biophysics, University of North Carolina at Chapel Hill, Chapel Hill, NC 27599, USA
| | - Denis Tsygankov
- Wallace H. Coulter Department of Biomedical Engineering, Georgia Institute of Technology and Emory University, Atlanta, GA 30332, USA
| | - Timothy C Elston
- Department of Pharmacology, University of North Carolina at Chapel Hill, Chapel Hill, NC 27599, USA
- Computational Medicine Program, University of North Carolina at Chapel Hill, Chapel Hill, NC 27599, USA
| |
Collapse
|
2
|
Raja Venkatesh A, Le KH, Weld DM, Brandman O. Diffusive lensing as a mechanism of intracellular transport and compartmentalization. eLife 2024; 12:RP89794. [PMID: 38896469 PMCID: PMC11186627 DOI: 10.7554/elife.89794] [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] [Indexed: 06/21/2024] Open
Abstract
While inhomogeneous diffusivity has been identified as a ubiquitous feature of the cellular interior, its implications for particle mobility and concentration at different length scales remain largely unexplored. In this work, we use agent-based simulations of diffusion to investigate how heterogeneous diffusivity affects the movement and concentration of diffusing particles. We propose that a nonequilibrium mode of membrane-less compartmentalization arising from the convergence of diffusive trajectories into low-diffusive sinks, which we call 'diffusive lensing,' is relevant for living systems. Our work highlights the phenomenon of diffusive lensing as a potentially key driver of mesoscale dynamics in the cytoplasm, with possible far-reaching implications for biochemical processes.
Collapse
Affiliation(s)
- Achuthan Raja Venkatesh
- Department of Biochemistry, Stanford UniversityStanfordUnited States
- Department of Biological Sciences, Indian Institute of Science Education and Research (IISER) MohaliMohaliIndia
| | - Kathy H Le
- Department of Biochemistry, Stanford UniversityStanfordUnited States
| | - David M Weld
- Department of Physics, University of California, Santa BarbaraSanta BarbaraUnited States
| | - Onn Brandman
- Department of Biochemistry, Stanford UniversityStanfordUnited States
| |
Collapse
|
3
|
Kubečka J, Besel V, Neefjes I, Knattrup Y, Kurtén T, Vehkamäki H, Elm J. Computational Tools for Handling Molecular Clusters: Configurational Sampling, Storage, Analysis, and Machine Learning. ACS OMEGA 2023; 8:45115-45128. [PMID: 38046354 PMCID: PMC10688175 DOI: 10.1021/acsomega.3c07412] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 09/25/2023] [Revised: 10/25/2023] [Accepted: 10/26/2023] [Indexed: 12/05/2023]
Abstract
Computational modeling of atmospheric molecular clusters requires a comprehensive understanding of their complex configurational spaces, interaction patterns, stabilities against fragmentation, and even dynamic behaviors. To address these needs, we introduce the Jammy Key framework, a collection of automated scripts that facilitate and streamline molecular cluster modeling workflows. Jammy Key handles file manipulations between varieties of integrated third-party programs. The framework is divided into three main functionalities: (1) Jammy Key for configurational sampling (JKCS) to perform systematic configurational sampling of molecular clusters, (2) Jammy Key for quantum chemistry (JKQC) to analyze commonly used quantum chemistry output files and facilitate database construction, handling, and analysis, and (3) Jammy Key for machine learning (JKML) to manage machine learning methods in optimizing molecular cluster modeling. This automation and machine learning utilization significantly reduces manual labor, greatly speeds up the search for molecular cluster configurations, and thus increases the number of systems that can be studied. Following the example of the Atmospheric Cluster Database (ACDB) of Elm (ACS Omega, 4, 10965-10984, 2019), the molecular clusters modeled in our group using the Jammy Key framework have been stored in an improved online GitHub repository named ACDB 2.0. In this work, we present the Jammy Key package alongside its assorted applications, which underline its versatility. Using several illustrative examples, we discuss how to choose appropriate combinations of methodologies for treating particular cluster types, including reactive, multicomponent, charged, or radical clusters, as well as clusters containing flexible or multiconformer monomers or heavy atoms. Finally, we present a detailed example of using the tools for atmospheric acid-base clusters.
Collapse
Affiliation(s)
- Jakub Kubečka
- Aarhus
University, Department of Chemistry, Langelandsgade 140, Aarhus 8000, Denmark
| | - Vitus Besel
- University
of Helsinki, Institute for Atmospheric and
Earth System Research/Physics, Faculty of Science, P.O. Box 64, Helsinki 00140, Finland
| | - Ivo Neefjes
- University
of Helsinki, Institute for Atmospheric and
Earth System Research/Physics, Faculty of Science, P.O. Box 64, Helsinki 00140, Finland
| | - Yosef Knattrup
- Aarhus
University, Department of Chemistry, Langelandsgade 140, Aarhus 8000, Denmark
| | - Theo Kurtén
- University
of Helsinki, Institute for Atmospheric and
Earth System Research/Chemistry, Faculty of Science, P.O. Box 64, Helsinki 00140, Finland
| | - Hanna Vehkamäki
- University
of Helsinki, Institute for Atmospheric and
Earth System Research/Physics, Faculty of Science, P.O. Box 64, Helsinki 00140, Finland
| | - Jonas Elm
- Aarhus
University, Department of Chemistry, Langelandsgade 140, Aarhus 8000, Denmark
| |
Collapse
|
4
|
Waters CS, Angenent SB, Altschuler SJ, Wu LF. A PINK1 input threshold arises from positive feedback in the PINK1/Parkin mitophagy decision circuit. Cell Rep 2023; 42:113260. [PMID: 37851575 PMCID: PMC10668033 DOI: 10.1016/j.celrep.2023.113260] [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: 11/11/2022] [Revised: 08/25/2023] [Accepted: 09/28/2023] [Indexed: 10/20/2023] Open
Abstract
Mechanisms that prevent accidental activation of the PINK1/Parkin mitophagy circuit on healthy mitochondria are poorly understood. On the surface of damaged mitochondria, PINK1 accumulates and acts as the input signal to a positive feedback loop of Parkin recruitment, which in turn promotes mitochondrial degradation via mitophagy. However, PINK1 is also present on healthy mitochondria, where it could errantly recruit Parkin and thereby activate this positive feedback loop. Here, we explore emergent properties of the PINK1/Parkin circuit by quantifying the relationship between mitochondrial PINK1 concentrations and Parkin recruitment dynamics. We find that Parkin is recruited to mitochondria only if PINK1 levels exceed a threshold and then only after a delay that is inversely proportional to PINK1 levels. Furthermore, these two regulatory properties arise from the input-coupled positive feedback topology of the PINK1/Parkin circuit. These results outline an intrinsic mechanism by which the PINK1/Parkin circuit can avoid errant activation on healthy mitochondria.
Collapse
Affiliation(s)
- Christopher S Waters
- Department of Pharmaceutical Chemistry, University of California, San Francisco, San Francisco, CA 94158, USA
| | - Sigurd B Angenent
- Mathematics Department, University of Wisconsin Madison, Madison, WI 53706, USA
| | - Steven J Altschuler
- Department of Pharmaceutical Chemistry, University of California, San Francisco, San Francisco, CA 94158, USA.
| | - Lani F Wu
- Department of Pharmaceutical Chemistry, University of California, San Francisco, San Francisco, CA 94158, USA.
| |
Collapse
|
5
|
Guan K, Curtis ER, Lew DJ, Elston TC. Particle-based simulations reveal two positive feedback loops allow relocation and stabilization of the polarity site during yeast mating. PLoS Comput Biol 2023; 19:e1011523. [PMID: 37782676 PMCID: PMC10569529 DOI: 10.1371/journal.pcbi.1011523] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/04/2023] [Revised: 10/12/2023] [Accepted: 09/17/2023] [Indexed: 10/04/2023] Open
Abstract
Many cells adjust the direction of polarized growth or migration in response to external directional cues. The yeast Saccharomyces cerevisiae orient their cell fronts (also called polarity sites) up pheromone gradients in the course of mating. However, the initial polarity site is often not oriented towards the eventual mating partner, and cells relocate the polarity site in an indecisive manner before developing a stable orientation. During this reorientation phase, the polarity site displays erratic assembly-disassembly behavior and moves around the cell cortex. The mechanisms underlying this dynamic behavior remain poorly understood. Particle-based simulations of the core polarity circuit revealed that molecular-level fluctuations are unlikely to overcome the strong positive feedback required for polarization and generate relocating polarity sites. Surprisingly, inclusion of a second pathway that promotes polarity site orientation generated relocating polarity sites with properties similar to those observed experimentally. This pathway forms a second positive feedback loop involving the recruitment of receptors to the cell membrane and couples polarity establishment to gradient sensing. This second positive feedback loop also allows cells to stabilize their polarity site once the site is aligned with the pheromone gradient.
Collapse
Affiliation(s)
- Kaiyun Guan
- Curriculum in Bioinformatics and Computational Biology, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina, United States of America
| | - Erin R. Curtis
- Department of Biology, Massachusetts Institute of Technology, Cambridge, Massachusetts, United States of America
| | - Daniel J. Lew
- Department of Biology, Massachusetts Institute of Technology, Cambridge, Massachusetts, United States of America
| | - Timothy C. Elston
- Department of Pharmacology and Computational Medicine Program, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina, United States of America
| |
Collapse
|
6
|
Bernoff AJ, Jilkine A, Navarro Hernández A, Lindsay AE. Single-cell directional sensing from just a few receptor binding events. Biophys J 2023; 122:3108-3116. [PMID: 37355773 PMCID: PMC10432224 DOI: 10.1016/j.bpj.2023.06.015] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/11/2023] [Revised: 06/08/2023] [Accepted: 06/20/2023] [Indexed: 06/26/2023] Open
Abstract
Identifying the directionality of signaling sources from noisy input to membrane receptors is an essential task performed by many cell types. A variety of models have been proposed to explain directional sensing in cells. However, many of these require significant computational and memory capacities for the cell. We propose and analyze a simple mechanism in which a cell adopts the direction associated with the first few membrane binding events. This model yields an accurate angular estimate to the source long before steady state is reached in biologically relevant scenarios. Our proposed mechanism allows for reliable estimates of the directionality of external signals using temporal information and assumes minimal computational capacities of the cell.
Collapse
Affiliation(s)
- Andrew J Bernoff
- Department of Mathematics, Harvey Mudd College, Claremont, California
| | - Alexandra Jilkine
- Department of Applied & Computational Mathematics & Statistics, University of Notre Dame, South Bend, Indiana
| | - Adrián Navarro Hernández
- Department of Applied & Computational Mathematics & Statistics, University of Notre Dame, South Bend, Indiana
| | - Alan E Lindsay
- Department of Applied & Computational Mathematics & Statistics, University of Notre Dame, South Bend, Indiana.
| |
Collapse
|
7
|
Lang CF, Munro EM. Oligomerization of peripheral membrane proteins provides tunable control of cell surface polarity. Biophys J 2022; 121:4543-4559. [PMID: 36815706 PMCID: PMC9750853 DOI: 10.1016/j.bpj.2022.10.035] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/18/2022] [Revised: 08/31/2022] [Accepted: 10/24/2022] [Indexed: 11/02/2022] Open
Abstract
Asymmetric distributions of peripheral membrane proteins define cell polarity across all kingdoms of life. Non-linear positive feedback on membrane binding is essential to amplify and stabilize these asymmetries, but how specific molecular sources of non-linearity shape polarization dynamics remains poorly understood. Here we show that the ability to oligomerize, which is common to many peripheral membrane proteins, can play a profound role in shaping polarization dynamics in simple feedback circuits. We show that size-dependent binding avidity and mobility of membrane-bound oligomers endow polarity circuits with several key properties. Size-dependent membrane binding avidity confers a form of positive feedback on the accumulation of oligomer subunits. Although insufficient by itself, this sharply reduces the amount of additional feedback required for spontaneous emergence and stable maintenance of polarized states. Size-dependent oligomer mobility makes symmetry breaking and stable polarity more robust with respect to variation in subunit diffusivities and cell sizes, and slows the approach to a final stable spatial distribution, allowing cells to "remember" polarity boundaries imposed by transient external cues. Together, these findings reveal how oligomerization of peripheral membrane proteins can provide powerful and highly tunable sources of non-linear feedback in biochemical circuits that govern cell surface polarity. Given its prevalence and widespread involvement in cell polarity, we speculate that self-oligomerization may have provided an accessible path to evolving simple polarity circuits.
Collapse
Affiliation(s)
- Charles F Lang
- Department of Molecular Genetics and Cell Biology, University of Chicago, Chicago, Illinois; Committee on Genetics, Genomics, and Systems Biology, University of Chicago, Chicago, Illinois
| | - Edwin M Munro
- Department of Molecular Genetics and Cell Biology, University of Chicago, Chicago, Illinois; Committee on Genetics, Genomics, and Systems Biology, University of Chicago, Chicago, Illinois.
| |
Collapse
|
8
|
García-Tejera R, Schumacher L, Grima R. Regulation of stem cell dynamics through volume exclusion. Proc Math Phys Eng Sci 2022. [DOI: 10.1098/rspa.2022.0376] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022] Open
Abstract
The maintenance and regeneration of adult tissues rely on the self-renewal of stem cells. Regeneration without over-proliferation requires precise regulation of the stem cell proliferation and differentiation rates. The nature of such regulatory mechanisms in different tissues, and how to incorporate them in models of stem cell population dynamics, is incompletely understood. The critical birth-death (CBD) process is widely used to model stem cell populations, capturing key phenomena, such as scaling laws in clone size distributions. However, the CBD process neglects regulatory mechanisms. Here, we propose the birth-death process with volume exclusion (vBD), a variation of the birth-death process that considers crowding effects, such as may arise due to limited space in a stem cell niche. While the deterministic rate equations predict a single non-trivial attracting steady state, the master equation predicts extinction and transient distributions of stem cell numbers with three possible behaviours: long-lived quasi-steady state (QSS), and short-lived bimodal or unimodal distributions. In all cases, we approximate solutions to the vBD master equation using a renormalized system-size expansion, QSS approximation and the Wentzel–Kramers–Brillouin method. Our study suggests that the size distribution of a stem cell population bears signatures that are useful to detect negative feedback mediated via volume exclusion.
Collapse
Affiliation(s)
- Rodrigo García-Tejera
- Centre for Regenerative Medicine, University of Edinburgh, 5 Little France Dr, Edinburgh EH16 4UU, UK
- School of Biological Sciences, University of Edinburgh, Kings Buildings, Mayfield Road, Edinburgh EH9 3JF, UK
| | - Linus Schumacher
- Centre for Regenerative Medicine, University of Edinburgh, 5 Little France Dr, Edinburgh EH16 4UU, UK
- School of Biological Sciences, University of Edinburgh, Kings Buildings, Mayfield Road, Edinburgh EH9 3JF, UK
| | - Ramon Grima
- School of Biological Sciences, University of Edinburgh, Kings Buildings, Mayfield Road, Edinburgh EH9 3JF, UK
| |
Collapse
|
9
|
Lawson MJ, Drawert B, Petzold L, Yi TM. A positive feedback loop involving the Spa2 SHD domain contributes to focal polarization. PLoS One 2022; 17:e0263347. [PMID: 35134079 PMCID: PMC8824340 DOI: 10.1371/journal.pone.0263347] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/03/2021] [Accepted: 01/16/2022] [Indexed: 11/18/2022] Open
Abstract
Focal polarization is necessary for finely arranged cell-cell interactions. The yeast mating projection, with its punctate polarisome, is a good model system for this process. We explored the critical role of the polarisome scaffold protein Spa2 during yeast mating with a hypothesis motivated by mathematical modeling and tested by in vivo experiments. Our simulations predicted that two positive feedback loops generate focal polarization, including a novel feedback pathway involving the N-terminal domain of Spa2. We characterized the latter using loss-of-function and gain-of-function mutants. The N-terminal region contains a Spa2 Homology Domain (SHD) which is conserved from yeast to humans, and when mutated largely reproduced the spa2Δ phenotype. Our work together with published data show that the SHD domain recruits Msb3/4 that stimulates Sec4-mediated transport of Bud6 to the polarisome. There, Bud6 activates Bni1-catalyzed actin cable formation, recruiting more Spa2 and completing the positive feedback loop. We demonstrate that disrupting this loop at any point results in morphological defects. Gain-of-function perturbations partially restored focal polarization in a spa2 loss-of-function mutant without restoring localization of upstream components, thus supporting the pathway order. Thus, we have collected data consistent with a novel positive feedback loop that contributes to focal polarization during pheromone-induced polarization in yeast.
Collapse
Affiliation(s)
- Michael J. Lawson
- Division of Biology and Biological Engineering, California Institute of Technology, Pasadena, CA, United States of America
| | - Brian Drawert
- Department of Computer Science, University of North Carolina Asheville, Asheville, NC, United States of America
| | - Linda Petzold
- Department of Computer Science, University of California, Santa Barbara, Santa Barbara, CA, United States of America
| | - Tau-Mu Yi
- Molecular, Cellular, and Developmental Biology, 3131 Biological Sciences II, University of California, Santa Barbara, Santa Barbara, CA, United States of America
- * E-mail:
| |
Collapse
|
10
|
Abstract
Accurate decoding of spatial chemical landscapes is critical for many cell functions. Eukaryotic cells decode local chemical gradients to orient growth or movement in productive directions. Recent work on yeast model systems, whose gradient sensing pathways display much less complexity than those in animal cells, has suggested new paradigms for how these very small cells successfully exploit information in noisy and dynamic pheromone gradients to identify their mates. Pheromone receptors regulate a polarity circuit centered on the conserved Rho-family GTPase, Cdc42. The polarity circuit contains both positive and negative feedback pathways, allowing spontaneous symmetry breaking and also polarity site disassembly and relocation. Cdc42 orients the actin cytoskeleton, leading to focused vesicle traffic that promotes movement of the polarity site and also reshapes the cortical distribution of receptors at the cell surface. In this article, we review the advances from work on yeasts and compare them with the excitable signaling pathways that have been revealed in chemotactic animal cells. Expected final online publication date for the Annual Review of Biophysics, Volume 51 is May 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.
Collapse
Affiliation(s)
- Debraj Ghose
- Department of Pharmacology and Cancer Biology, Duke University Medical Center, Durham, North Carolina, USA;
| | - Timothy Elston
- Department of Pharmacology, University of North Carolina at Chapel Hill, North Carolina, USA
| | - Daniel Lew
- Department of Pharmacology and Cancer Biology, Duke University Medical Center, Durham, North Carolina, USA;
| |
Collapse
|
11
|
Ramirez SA, Pablo M, Burk S, Lew DJ, Elston TC. A novel stochastic simulation approach enables exploration of mechanisms for regulating polarity site movement. PLoS Comput Biol 2021; 17:e1008525. [PMID: 34264926 PMCID: PMC8315557 DOI: 10.1371/journal.pcbi.1008525] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/03/2020] [Revised: 07/27/2021] [Accepted: 06/24/2021] [Indexed: 12/23/2022] Open
Abstract
Cells polarize their movement or growth toward external directional cues in many different contexts. For example, budding yeast cells grow toward potential mating partners in response to pheromone gradients. Directed growth is controlled by polarity factors that assemble into clusters at the cell membrane. The clusters assemble, disassemble, and move between different regions of the membrane before eventually forming a stable polarity site directed toward the pheromone source. Pathways that regulate clustering have been identified but the molecular mechanisms that regulate cluster mobility are not well understood. To gain insight into the contribution of chemical noise to cluster behavior we simulated clustering using the reaction-diffusion master equation (RDME) framework to account for molecular-level fluctuations. RDME simulations are a computationally efficient approximation, but their results can diverge from the underlying microscopic dynamics. We implemented novel concentration-dependent rate constants that improved the accuracy of RDME-based simulations, allowing us to efficiently investigate how cluster dynamics might be regulated. Molecular noise was effective in relocating clusters when the clusters contained low numbers of limiting polarity factors, and when Cdc42, the central polarity regulator, exhibited short dwell times at the polarity site. Cluster stabilization occurred when abundances or binding rates were altered to either lengthen dwell times or increase the number of polarity molecules in the cluster. We validated key results using full 3D particle-based simulations. Understanding the mechanisms cells use to regulate the dynamics of polarity clusters should provide insights into how cells dynamically track external directional cues. Cells localize polarity molecules in a small region of the plasma membrane forming a polarity cluster that directs functions such as migration, reproduction, and growth. Guided by external signals, these clusters move across the membrane allowing cells to reorient growth or motion. The polarity molecules continuously and randomly shuttle between the cluster and the cell cytosol and, as a result, the number and distribution of molecules at the cluster constantly changes. Here we present an improved stochastic simulation algorithm to investigate how such molecular-scale fluctuations induce cluster movement across the cell membrane. Unexpectedly, cluster mobility does not correlate with variations in total molecule abundance within the cluster, but rather with changes in the spatial distribution of molecules that form the cluster. Cluster motion is faster when polarity molecules are scarce and when they shuttle rapidly between the cluster and the cytosol. Our results suggest that cells control cluster mobility by regulating the abundance of polarity molecules and biochemical reactions that affect the time molecules spend at the cluster. We provide insights into how cells harness random molecular behavior to perform functions important for survival, such as detecting the direction of external signals.
Collapse
Affiliation(s)
- Samuel A. Ramirez
- Department of Pharmacology and Computational Medicine Program, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina, United States of America
- * E-mail: (SAR); (TCE)
| | - Michael Pablo
- Department of Chemistry, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina, United States of America
- Program in Molecular and Cellular Biophysics, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina, United States of America
| | - Sean Burk
- Department of Pharmacology and Computational Medicine Program, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina, United States of America
| | - Daniel J. Lew
- Department of Pharmacology and Cancer Biology, Duke University, Durham, North Carolina, United States of America
| | - Timothy C. Elston
- Department of Pharmacology and Computational Medicine Program, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina, United States of America
- * E-mail: (SAR); (TCE)
| |
Collapse
|
12
|
Banavar SP, Trogdon M, Drawert B, Yi TM, Petzold LR, Campàs O. Coordinating cell polarization and morphogenesis through mechanical feedback. PLoS Comput Biol 2021; 17:e1007971. [PMID: 33507956 PMCID: PMC7872284 DOI: 10.1371/journal.pcbi.1007971] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/15/2020] [Revised: 02/09/2021] [Accepted: 12/21/2020] [Indexed: 12/30/2022] Open
Abstract
Many cellular processes require cell polarization to be maintained as the cell changes shape, grows or moves. Without feedback mechanisms relaying information about cell shape to the polarity molecular machinery, the coordination between cell polarization and morphogenesis, movement or growth would not be possible. Here we theoretically and computationally study the role of a genetically-encoded mechanical feedback (in the Cell Wall Integrity pathway) as a potential coordination mechanism between cell morphogenesis and polarity during budding yeast mating projection growth. We developed a coarse-grained continuum description of the coupled dynamics of cell polarization and morphogenesis as well as 3D stochastic simulations of the molecular polarization machinery in the evolving cell shape. Both theoretical approaches show that in the absence of mechanical feedback (or in the presence of weak feedback), cell polarity cannot be maintained at the projection tip during growth, with the polarization cap wandering off the projection tip, arresting morphogenesis. In contrast, for mechanical feedback strengths above a threshold, cells can robustly maintain cell polarization at the tip and simultaneously sustain mating projection growth. These results indicate that the mechanical feedback encoded in the Cell Wall Integrity pathway can provide important positional information to the molecular machinery in the cell, thereby enabling the coordination of cell polarization and morphogenesis.
Collapse
Affiliation(s)
- Samhita P. Banavar
- Department of Physics, University of California, University of California, Santa Barbara, California, United States of America
- California NanoSystems Institute, University of California, Santa Barbara, California, United States of America
| | - Michael Trogdon
- Department of Mechanical Engineering, University of California, Santa Barbara, California, United States of America
| | - Brian Drawert
- Department of Computer Science, University of North Carolina, Asheville, North Carolina, United States of America
| | - Tau-Mu Yi
- Department of Molecular, Cell and Developmental Biology, University of California, Santa Barbara, California, United States of America
| | - Linda R. Petzold
- Department of Mechanical Engineering, University of California, Santa Barbara, California, United States of America
- Center for Bioengineering, University of California, Santa Barbara, California, United States of America
| | - Otger Campàs
- California NanoSystems Institute, University of California, Santa Barbara, California, United States of America
- Department of Mechanical Engineering, University of California, Santa Barbara, California, United States of America
- Department of Molecular, Cell and Developmental Biology, University of California, Santa Barbara, California, United States of America
- Center for Bioengineering, University of California, Santa Barbara, California, United States of America
- Cluster of Excellence Physics of Life, TU Dresden, Dresden, Germany
| |
Collapse
|
13
|
Ghose D, Lew D. Mechanistic insights into actin-driven polarity site movement in yeast. Mol Biol Cell 2020; 31:1085-1102. [PMID: 32186970 PMCID: PMC7346724 DOI: 10.1091/mbc.e20-01-0040] [Citation(s) in RCA: 18] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/15/2020] [Revised: 03/02/2020] [Accepted: 03/10/2020] [Indexed: 11/11/2022] Open
Abstract
Directed cell growth or migration are critical for the development and function of many eukaryotic cells. These cells develop a dynamic "front" (also called "polarity site") that can change direction. Polarity establishment involves autocatalytic accumulation of polarity regulators, including the conserved Rho-family GTPase Cdc42, but the mechanisms underlying polarity reorientation remain poorly understood. The tractable model yeast, Saccharomyces cerevisiae, relocates its polarity site when searching for mating partners. Relocation requires polymerized actin, and is thought to involve actin-mediated vesicle traffic to the polarity site. In this study, we provide a quantitative characterization of spontaneous polarity site movement as a search process and use a mechanistic computational model that combines polarity protein biochemical interactions with vesicle trafficking to probe how various processes might affect polarity site movement. Our findings identify two previously documented features of yeast vesicle traffic as being particularly relevant to such movement: tight spatial focusing of exocytosis enhances the directional persistence of movement, and association of Cdc42-directed GTPase-Activating Proteins with secretory vesicles increases the distance moved. Furthermore, we suggest that variation in the rate of exocytosis beyond simple Poisson dynamics may be needed to fully account for the characteristics of polarity site movement in vivo.
Collapse
Affiliation(s)
- Debraj Ghose
- Computational Biology and Bioinformatics, Duke University Medical Center, Durham, NC 27710
- Department of Pharmacology and Cancer Biology, Duke University Medical Center, Durham, NC 27710
| | - Daniel Lew
- Department of Pharmacology and Cancer Biology, Duke University Medical Center, Durham, NC 27710
| |
Collapse
|
14
|
Bezeljak U, Loya H, Kaczmarek B, Saunders TE, Loose M. Stochastic activation and bistability in a Rab GTPase regulatory network. Proc Natl Acad Sci U S A 2020; 117:6540-6549. [PMID: 32161136 PMCID: PMC7104049 DOI: 10.1073/pnas.1921027117] [Citation(s) in RCA: 18] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022] Open
Abstract
The eukaryotic endomembrane system is controlled by small GTPases of the Rab family, which are activated at defined times and locations in a switch-like manner. While this switch is well understood for an individual protein, how regulatory networks produce intracellular activity patterns is currently not known. Here, we combine in vitro reconstitution experiments with computational modeling to study a minimal Rab5 activation network. We find that the molecular interactions in this system give rise to a positive feedback and bistable collective switching of Rab5. Furthermore, we find that switching near the critical point is intrinsically stochastic and provide evidence that controlling the inactive population of Rab5 on the membrane can shape the network response. Notably, we demonstrate that collective switching can spread on the membrane surface as a traveling wave of Rab5 activation. Together, our findings reveal how biochemical signaling networks control vesicle trafficking pathways and how their nonequilibrium properties define the spatiotemporal organization of the cell.
Collapse
Affiliation(s)
- Urban Bezeljak
- Institute of Science and Technology Austria, 3400 Klosterneuburg, Austria
| | - Hrushikesh Loya
- Department of Electrical Engineering, Indian Institute of Technology Bombay, Mumbai, Maharashtra 400076, India
| | - Beata Kaczmarek
- Institute of Science and Technology Austria, 3400 Klosterneuburg, Austria
| | - Timothy E Saunders
- Mechanobiology Institute, National University of Singapore, Singapore 119077;
- Department of Biological Sciences, National University of Singapore, Singapore 119077
| | - Martin Loose
- Institute of Science and Technology Austria, 3400 Klosterneuburg, Austria;
| |
Collapse
|
15
|
Paquin-Lefebvre F, Xu B, DiPietro KL, Lindsay AE, Jilkine A. Pattern formation in a coupled membrane-bulk reaction-diffusion model for intracellular polarization and oscillations. J Theor Biol 2020; 497:110242. [PMID: 32179107 DOI: 10.1016/j.jtbi.2020.110242] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/26/2019] [Revised: 03/02/2020] [Accepted: 03/05/2020] [Indexed: 01/19/2023]
Abstract
Reaction-diffusion systems have been widely used to study spatio-temporal phenomena in cell biology, such as cell polarization. Coupled bulk-surface models naturally include compartmentalization of cytosolic and membrane-bound polarity molecules. Here we study the distribution of the polarity protein Cdc42 in a mass-conserved membrane-bulk model, and explore the effects of diffusion and spatial dimensionality on spatio-temporal pattern formation. We first analyze a one-dimensional (1-D) model for Cdc42 oscillations in fission yeast, consisting of two diffusion equations in the bulk domain coupled to nonlinear ODEs for binding kinetics at each end of the cell. In 1-D, our analysis reveals the existence of symmetric and asymmetric steady states, as well as anti-phase relaxation oscillations typical of slow-fast systems. We then extend our analysis to a two-dimensional (2-D) model with circular bulk geometry, for which species can either diffuse inside the cell or become bound to the membrane and undergo a nonlinear reaction-diffusion process. We also consider a nonlocal system of PDEs approximating the dynamics of the 2-D membrane-bulk model in the limit of fast bulk diffusion. In all three model variants we find that mass conservation selects perturbations of spatial modes that simply redistribute mass. In 1-D, only anti-phase oscillations between the two ends of the cell can occur, and in-phase oscillations are excluded. In higher dimensions, no radially symmetric oscillations are observed. Instead, the only instabilities are symmetry-breaking, either corresponding to stationary Turing instabilities, leading to the formation of stationary patterns, or to oscillatory Turing instabilities, leading to traveling and standing waves. Codimension-two Bogdanov-Takens bifurcations occur when the two distinct instabilities coincide, causing traveling waves to slow down and to eventually become stationary patterns. Our work clarifies the effect of geometry and dimensionality on behaviors observed in mass-conserved cell polarity models.
Collapse
Affiliation(s)
- Frédéric Paquin-Lefebvre
- Department of Mathematics and Institute of Applied Mathematics, University of British Columbia, Vancouver, Canada
| | - Bin Xu
- Department of Applied and Computational Mathematics and Statistics, University of Notre Dame, Notre Dame, IN, 46556, USA
| | - Kelsey L DiPietro
- Department of Applied and Computational Mathematics and Statistics, University of Notre Dame, Notre Dame, IN, 46556, USA; Sandia National Laboratories, NM, 46556, USA
| | - Alan E Lindsay
- Department of Applied and Computational Mathematics and Statistics, University of Notre Dame, Notre Dame, IN, 46556, USA
| | - Alexandra Jilkine
- Department of Applied and Computational Mathematics and Statistics, University of Notre Dame, Notre Dame, IN, 46556, USA.
| |
Collapse
|
16
|
Comparison of Deterministic and Stochastic Regime in a Model for Cdc42 Oscillations in Fission Yeast. Bull Math Biol 2019; 81:1268-1302. [PMID: 30756233 DOI: 10.1007/s11538-019-00573-5] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/27/2018] [Accepted: 01/29/2019] [Indexed: 01/13/2023]
Abstract
Oscillations occur in a wide variety of essential cellular processes, such as cell cycle progression, circadian clocks and calcium signaling in response to stimuli. It remains unclear how intrinsic stochasticity can influence these oscillatory systems. Here, we focus on oscillations of Cdc42 GTPase in fission yeast. We extend our previous deterministic model by Xu and Jilkine to construct a stochastic model, focusing on the fast diffusion case. We use SSA (Gillespie's algorithm) to numerically explore the low copy number regime in this model, and use analytical techniques to study the long-time behavior of the stochastic model and compare it to the equilibria of its deterministic counterpart. Numerical solutions suggest noisy limit cycles exist in the parameter regime in which the deterministic system converges to a stable limit cycle, and quasi-cycles exist in the parameter regime where the deterministic model has a damped oscillation. Near an infinite period bifurcation point, the deterministic model has a sustained oscillation, while stochastic trajectories start with an oscillatory mode and tend to approach deterministic steady states. In the low copy number regime, metastable transitions from oscillatory to steady behavior occur in the stochastic model. Our work contributes to the understanding of how stochastic chemical kinetics can affect a finite-dimensional dynamical system, and destabilize a deterministic steady state leading to oscillations.
Collapse
|
17
|
The Interplay of Structural and Cellular Biophysics Controls Clustering of Multivalent Molecules. Biophys J 2019; 116:560-572. [PMID: 30661665 DOI: 10.1016/j.bpj.2019.01.001] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/27/2018] [Revised: 12/24/2018] [Accepted: 01/02/2019] [Indexed: 12/12/2022] Open
Abstract
Dynamic molecular clusters are assembled through weak multivalent interactions and are platforms for cellular functions, especially receptor-mediated signaling. Clustering is also a prerequisite for liquid-liquid phase separation. It is not well understood, however, how molecular structure and cellular organization control clustering. Using coarse-grained kinetic Langevin dynamics, we performed computational experiments on a prototypical ternary system modeled after membrane-bound nephrin, the adaptor Nck1, and the actin nucleation promoting factor NWASP. Steady-state cluster size distributions favored stoichiometries that optimized binding (stoichiometry matching) but still were quite broad. At high concentrations, the system can be driven beyond the saturation boundary such that cluster size is limited only by the number of available molecules. This behavior would be predictive of phase separation. Domains close to binding sites sterically inhibited clustering much less than terminal domains because the latter effectively restrict access to the cluster interior. Increased flexibility of interacting molecules diminished clustering by shielding binding sites within compact conformations. Membrane association of nephrin increased the cluster size distribution in a density-dependent manner. These properties provide insights into how molecular ensembles function to localize and amplify cell signaling.
Collapse
|
18
|
Johnson ME. Modeling the Self-Assembly of Protein Complexes through a Rigid-Body Rotational Reaction-Diffusion Algorithm. J Phys Chem B 2018; 122:11771-11783. [PMID: 30256109 DOI: 10.1021/acs.jpcb.8b08339] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
The reaction-diffusion equations provide a powerful framework for modeling nonequilibrium, cell-scale dynamics over the long time scales that are inaccessible by traditional molecular modeling approaches. Single-particle reaction-diffusion offers the highest resolution technique for tracking such dynamics, but it has not been applied to the study of protein self-assembly due to its treatment of reactive species as single-point particles. Here, we develop a relatively simple but accurate approach for building rigid structure and rotation into single-particle reaction-diffusion methods, providing a rate-based method for studying protein self-assembly. Our simplifying assumption is that reactive collisions can be evaluated purely on the basis of the separations between the sites, and not their orientations. The challenge of evaluating reaction probabilities can then be performed using well-known equations based on translational diffusion in both 3D and 2D, by employing an effective diffusion constant we derive here. We show how our approach reproduces both the kinetics of association, which is altered by rotational diffusion, and the equilibrium of reversible association, which is not. Importantly, the macroscopic kinetics of association can be predicted on the basis of the microscopic parameters of our structurally resolved model, allowing for critical comparisons with theory and other rate-based simulations. We demonstrate this method for efficient, rate-based simulations of self-assembly of clathrin trimers, highlighting how formation of regular lattices impacts the kinetics of association.
Collapse
Affiliation(s)
- Margaret E Johnson
- TC Jenkins Department of Biophysics , The Johns Hopkins University , 3400 North Charles Street , Baltimore , Maryland 21218 , United States
| |
Collapse
|
19
|
Khalili B, Merlini L, Vincenzetti V, Martin SG, Vavylonis D. Exploration and stabilization of Ras1 mating zone: A mechanism with positive and negative feedbacks. PLoS Comput Biol 2018; 14:e1006317. [PMID: 30028833 PMCID: PMC6070293 DOI: 10.1371/journal.pcbi.1006317] [Citation(s) in RCA: 14] [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: 12/12/2017] [Revised: 08/01/2018] [Accepted: 06/21/2018] [Indexed: 02/07/2023] Open
Abstract
In mating fission yeast cells, sensing and response to extracellular pheromone concentrations occurs through an exploratory Cdc42 patch that stochastically samples the cell cortex before stabilizing towards a mating partner. Active Ras1 (Ras1-GTP), an upstream regulator of Cdc42, and Gap1, the GTPase-activating protein for Ras1, localize at the patch. We developed a reaction-diffusion model of Ras1 patch appearance and disappearance with a positive feedback by a Guanine nucleotide Exchange Factor (GEF) and Gap1 inhibition. The model is based on new estimates of Ras1-GDP, Ras1-GTP and Gap1 diffusion coefficients and rates of cytoplasmic exchange studied by FRAP. The model reproduces exploratory patch behavior and lack of Ras1 patch in cells lacking Gap1. Transition to a stable patch can occur by change of Gap1 rates constants or local increase of the positive feedback rate constants. The model predicts that the patch size and number of patches depend on the strength of positive and negative feedbacks. Measurements of Ras1 patch size and number in cells overexpressing the Ras1 GEF or Gap1 are consistent with the model. Unicellular fission yeasts mate by fusing with partners of the opposite mating type. Each pair member grows towards its selected partner that signals its presence through secreted pheromone. The process of partner selection occurs through an exploratory patch (containing activated signaling protein Cdc42 and upstream regulator Ras1) that assembles and disassembles on the cell cortex, stabilizing in regions of higher opposite pheromone concentration. We present a computational model of the molecular mechanisms driving the dynamical pattern of patch exploration and stabilization. The model is based on reaction and diffusion along the curved cell membrane, with diffusion coefficients measured experimentally. In the model, a positive Ras1 activation feedback loop generates a patch containing most of the activating protein (Ras1 GEF). The fast diffusing inhibitor Gap1 that is recruited locally from the cytoplasm spreads on the cell membrane, limiting patch size and causing its decay. Spontaneous reinitiation of Ras1 activation elsewhere on the cortex provides a mechanism for exploration. Transition of the system’s behavior to that of a single stable patch is possible upon simulated pheromone sensing. The computational model provides predictions for the number of patches and patch size dependence on parameters that we tested experimentally.
Collapse
Affiliation(s)
- Bita Khalili
- Department of Physics, Lehigh University, Bethlehem, Pennsylvania, United States of America
- Department of Fundamental Microbiology, University of Lausanne, Lausanne, Switzerland
| | - Laura Merlini
- Department of Fundamental Microbiology, University of Lausanne, Lausanne, Switzerland
| | - Vincent Vincenzetti
- Department of Fundamental Microbiology, University of Lausanne, Lausanne, Switzerland
| | - Sophie G. Martin
- Department of Fundamental Microbiology, University of Lausanne, Lausanne, Switzerland
| | - Dimitrios Vavylonis
- Department of Physics, Lehigh University, Bethlehem, Pennsylvania, United States of America
- * E-mail:
| |
Collapse
|
20
|
Trogdon M, Drawert B, Gomez C, Banavar SP, Yi TM, Campàs O, Petzold LR. The effect of cell geometry on polarization in budding yeast. PLoS Comput Biol 2018; 14:e1006241. [PMID: 29889845 PMCID: PMC6013239 DOI: 10.1371/journal.pcbi.1006241] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/13/2017] [Revised: 06/21/2018] [Accepted: 05/29/2018] [Indexed: 11/19/2022] Open
Abstract
The localization (or polarization) of proteins on the membrane during the mating of budding yeast (Saccharomyces cerevisiae) is an important model system for understanding simple pattern formation within cells. While there are many existing mathematical models of polarization, for both budding and mating, there are still many aspects of this process that are not well understood. In this paper we set out to elucidate the effect that the geometry of the cell can have on the dynamics of certain models of polarization. Specifically, we look at several spatial stochastic models of Cdc42 polarization that have been adapted from published models, on a variety of tip-shaped geometries, to replicate the shape change that occurs during the growth of the mating projection. We show here that there is a complex interplay between the dynamics of polarization and the shape of the cell. Our results show that while models of polarization can generate a stable polarization cap, its localization at the tip of mating projections is unstable, with the polarization cap drifting away from the tip of the projection in a geometry dependent manner. We also compare predictions from our computational results to experiments that observe cells with projections of varying lengths, and track the stability of the polarization cap. Lastly, we examine one model of actin polarization and show that it is unlikely, at least for the models studied here, that actin dynamics and vesicle traffic are able to overcome this effect of geometry.
Collapse
Affiliation(s)
- Michael Trogdon
- Department of Mechanical Engineering, University of California, Santa Barbara, Santa Barbara, California, United States of America
- * E-mail:
| | - Brian Drawert
- Department of Computer Science, University of North Carolina, Asheville, Asheville, North Carolina, United States of America
| | - Carlos Gomez
- Department of Molecular, Cell and Developmental Biology, University of California, Santa Barbara, Santa Barbara, California, United States of America
- California NanoSystems Institute, University of California, Santa Barbara, Santa Barbara, California, United States of America
| | - Samhita P. Banavar
- California NanoSystems Institute, University of California, Santa Barbara, Santa Barbara, California, United States of America
- Department of Physics, University of California, Santa Barbara, Santa Barbara, California, United States of America
| | - Tau-Mu Yi
- Department of Molecular, Cell and Developmental Biology, University of California, Santa Barbara, Santa Barbara, California, United States of America
| | - Otger Campàs
- Department of Mechanical Engineering, University of California, Santa Barbara, Santa Barbara, California, United States of America
- Department of Molecular, Cell and Developmental Biology, University of California, Santa Barbara, Santa Barbara, California, United States of America
- California NanoSystems Institute, University of California, Santa Barbara, Santa Barbara, California, United States of America
- Center for Bioengineering, University of California, Santa Barbara, Santa Barbara, California, United States of America
| | - Linda R. Petzold
- Department of Mechanical Engineering, University of California, Santa Barbara, Santa Barbara, California, United States of America
- Center for Bioengineering, University of California, Santa Barbara, Santa Barbara, California, United States of America
- Department of Computer Science, University of California, Santa Barbara, Santa Barbara, California, United States of America
| |
Collapse
|
21
|
Pablo M, Ramirez SA, Elston TC. Particle-based simulations of polarity establishment reveal stochastic promotion of Turing pattern formation. PLoS Comput Biol 2018. [PMID: 29529021 PMCID: PMC5864077 DOI: 10.1371/journal.pcbi.1006016] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
Polarity establishment, the spontaneous generation of asymmetric molecular distributions, is a crucial component of many cellular functions. Saccharomyces cerevisiae (yeast) undergoes directed growth during budding and mating, and is an ideal model organism for studying polarization. In yeast and many other cell types, the Rho GTPase Cdc42 is the key molecular player in polarity establishment. During yeast polarization, multiple patches of Cdc42 initially form, then resolve into a single front. Because polarization relies on strong positive feedback, it is likely that the amplification of molecular-level fluctuations underlies the generation of multiple nascent patches. In the absence of spatial cues, these fluctuations may be key to driving polarization. Here we used particle-based simulations to investigate the role of stochastic effects in a Turing-type model of yeast polarity establishment. In the model, reactions take place either between two molecules on the membrane, or between a cytosolic and a membrane-bound molecule. Thus, we developed a computational platform that explicitly simulates molecules at and near the cell membrane, and implicitly handles molecules away from the membrane. To evaluate stochastic effects, we compared particle simulations to deterministic reaction-diffusion equation simulations. Defining macroscopic rate constants that are consistent with the microscopic parameters for this system is challenging, because diffusion occurs in two dimensions and particles exchange between the membrane and cytoplasm. We address this problem by empirically estimating macroscopic rate constants from appropriately designed particle-based simulations. Ultimately, we find that stochastic fluctuations speed polarity establishment and permit polarization in parameter regions predicted to be Turing stable. These effects can operate at Cdc42 abundances expected of yeast cells, and promote polarization on timescales consistent with experimental results. To our knowledge, our work represents the first particle-based simulations of a model for yeast polarization that is based on a Turing mechanism. Many cells need to generate and maintain biochemical signals in specific subcellular regions. This phenomenon is broadly called polarity establishment, and is important in fundamental processes such as cell migration and differentiation. A key polarity factor found in diverse organisms, including yeast and humans, is the protein Cdc42. In yeast, Cdc42-dependent polarization occurs through a self-reinforcing biochemical signaling loop. Directional cues can guide polarity establishment, but interestingly, yeast can polarize in the absence of such a cue. The mechanism thought to underlie this symmetry breaking involves the amplification of inhomogeneities in molecular distributions that arise from molecular-level fluctuations. We investigated the effects of random fluctuations on polarization by performing particle-based simulations of the Cdc42 signaling network. We found that fluctuations can facilitate polarization, allowing faster polarization, and polarization over a broader range of concentrations. Our observations may help understand how polarity works in other systems.
Collapse
Affiliation(s)
- Michael Pablo
- Department of Chemistry, The University of North Carolina, Chapel Hill, NC, United States of America
- Program in Molecular and Cellular Biophysics, The University of North Carolina, Chapel Hill, NC, United States of America
| | - Samuel A. Ramirez
- Department of Pharmacology, The University of North Carolina, Chapel Hill, NC, United States of America
| | - Timothy C. Elston
- Department of Pharmacology, The University of North Carolina, Chapel Hill, NC, United States of America
- * E-mail:
| |
Collapse
|
22
|
Andrews SS. Smoldyn: particle-based simulation with rule-based modeling, improved molecular interaction and a library interface. Bioinformatics 2017; 33:710-717. [PMID: 28365760 DOI: 10.1093/bioinformatics/btw700] [Citation(s) in RCA: 53] [Impact Index Per Article: 7.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/16/2016] [Accepted: 11/03/2016] [Indexed: 12/17/2022] Open
Abstract
Motivation Smoldyn is a spatial and stochastic biochemical simulator. It treats each molecule of interest as an individual particle in continuous space, simulating molecular diffusion, molecule-membrane interactions and chemical reactions, all with good accuracy. This article presents several new features. Results Smoldyn now supports two types of rule-based modeling. These are a wildcard method, which is very convenient, and the BioNetGen package with extensions for spatial simulation, which is better for complicated models. Smoldyn also includes new algorithms for simulating the diffusion of surface-bound molecules and molecules with excluded volume. Both are exact in the limit of short time steps and reasonably good with longer steps. In addition, Smoldyn supports single-molecule tracking simulations. Finally, the Smoldyn source code can be accessed through a C/C ++ language library interface. Availability and Implementation Smoldyn software, documentation, code, and examples are at http://www.smoldyn.org . Contact steven.s.andrews@gmail.com.
Collapse
Affiliation(s)
- Steven S Andrews
- Basic Sciences Division, Fred Hutchinson Cancer Research Center, Seattle, WA 98109, USA.,Isaac Newton Institute for Mathematical Sciences, Cambridge CB3 0EH, UK
| |
Collapse
|
23
|
Foteinopoulos P, Mulder BM. A microtubule-based minimal model for spontaneous and persistent spherical cell polarity. PLoS One 2017; 12:e0184706. [PMID: 28931032 PMCID: PMC5607169 DOI: 10.1371/journal.pone.0184706] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/20/2016] [Accepted: 08/22/2017] [Indexed: 11/18/2022] Open
Abstract
We propose a minimal model for the spontaneous and persistent generation of polarity in a spherical cell based on dynamic microtubules and a single mobile molecular component. This component, dubbed the polarity factor, binds to microtubules nucleated from a centrosome located in the center of the cell, is subsequently delivered to the cell membrane, where it diffuses until it unbinds. The only feedback mechanism we impose is that the residence time of the microtubules at the membrane increases with the local density of the polarity factor. We show analytically that this system supports a stable unipolar symmetry-broken state for a wide range of parameters. We validate the predictions of the model by 2D particle-based simulations. Our model provides a route towards the creation of polarity in a minimal cell-like environment using a biochemical reconstitution approach.
Collapse
Affiliation(s)
| | - Bela M. Mulder
- Systems Biophysics Department, Institute AMOLF, Amsterdam, the Netherlands
- * E-mail:
| |
Collapse
|
24
|
Hegemann B, Peter M. Local sampling paints a global picture: Local concentration measurements sense direction in complex chemical gradients. Bioessays 2017; 39. [PMID: 28556309 DOI: 10.1002/bies.201600134] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Abstract
Detecting and interpreting extracellular spatial signals is essential for cellular orientation within complex environments, such as during directed cell migration or growth in multicellular development. Although the molecular understanding of how cells read spatial signals like chemical gradients is still lacking, recent work has revealed that stochastic processes at different temporal and spatial scales are at the core of this gradient sensing process in a wide range of eukaryotes. Fast biochemical reactions like those underlying GTPase activity dynamics form a functional module together with slower cell morphological changes driven by membrane remodelling. This biochemical-morphological module explores the environment by stochastic local concentration sampling to determine the source of the gradient signal, enabling efficient signal detection and interpretation before polarised growth or migration towards the gradient source is initiated. Here we review recent data describing local sampling and propose a model of local fast and slow feedback counteracted by gradient-dependent substrate limitation to be at the core of gradient sensing by local sampling.
Collapse
Affiliation(s)
- Björn Hegemann
- Department of Biology, Institute of Biochemistry, ETH Zurich, Zürich, Switzerland
| | - Matthias Peter
- Department of Biology, Institute of Biochemistry, ETH Zurich, Zürich, Switzerland
| |
Collapse
|
25
|
Abstract
In this paper we construct and analyze a model of cell receptor aggregation. Experiments have shown that receptors in an aggregated state have greatly reduced mobility. We model the effects of this reduced mobility with a density dependent diffusion and study the impact of density dependent diffusion on aggregate formation in a one-dimensional domain. Critical values of receptor diffusivity and receptor activation are found and compared with numerical simulations. We find that the role of density dependant diffusion is quite limited in the formation of aggregate structures. In the case of receptor activation, the analytical results agree very well with the numerical calculations. Finally, we consider our model in higher dimensional domains. In this case our analysis is primarily numerical.
Collapse
|
26
|
Schaff JC, Gao F, Li Y, Novak IL, Slepchenko BM. Numerical Approach to Spatial Deterministic-Stochastic Models Arising in Cell Biology. PLoS Comput Biol 2016; 12:e1005236. [PMID: 27959915 PMCID: PMC5154471 DOI: 10.1371/journal.pcbi.1005236] [Citation(s) in RCA: 25] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/23/2016] [Accepted: 11/03/2016] [Indexed: 01/01/2023] Open
Abstract
Hybrid deterministic-stochastic methods provide an efficient alternative to a fully stochastic treatment of models which include components with disparate levels of stochasticity. However, general-purpose hybrid solvers for spatially resolved simulations of reaction-diffusion systems are not widely available. Here we describe fundamentals of a general-purpose spatial hybrid method. The method generates realizations of a spatially inhomogeneous hybrid system by appropriately integrating capabilities of a deterministic partial differential equation solver with a popular particle-based stochastic simulator, Smoldyn. Rigorous validation of the algorithm is detailed, using a simple model of calcium 'sparks' as a testbed. The solver is then applied to a deterministic-stochastic model of spontaneous emergence of cell polarity. The approach is general enough to be implemented within biologist-friendly software frameworks such as Virtual Cell.
Collapse
Affiliation(s)
- James C. Schaff
- Richard D. Berlin Center for Cell Analysis and Modeling, Department of Cell Biology, University of Connecticut Health Center, Farmington, Connecticut, United States of America
| | - Fei Gao
- Richard D. Berlin Center for Cell Analysis and Modeling, Department of Cell Biology, University of Connecticut Health Center, Farmington, Connecticut, United States of America
| | - Ye Li
- Richard D. Berlin Center for Cell Analysis and Modeling, Department of Cell Biology, University of Connecticut Health Center, Farmington, Connecticut, United States of America
| | - Igor L. Novak
- Richard D. Berlin Center for Cell Analysis and Modeling, Department of Cell Biology, University of Connecticut Health Center, Farmington, Connecticut, United States of America
| | - Boris M. Slepchenko
- Richard D. Berlin Center for Cell Analysis and Modeling, Department of Cell Biology, University of Connecticut Health Center, Farmington, Connecticut, United States of America
| |
Collapse
|
27
|
Drawert B, Hellander S, Trogdon M, Yi TM, Petzold L. A framework for discrete stochastic simulation on 3D moving boundary domains. J Chem Phys 2016; 145:184113. [PMID: 27846706 PMCID: PMC5116030 DOI: 10.1063/1.4967338] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/24/2016] [Accepted: 10/26/2016] [Indexed: 12/20/2022] Open
Abstract
We have developed a method for modeling spatial stochastic biochemical reactions in complex, three-dimensional, and time-dependent domains using the reaction-diffusion master equation formalism. In particular, we look to address the fully coupled problems that arise in systems biology where the shape and mechanical properties of a cell are determined by the state of the biochemistry and vice versa. To validate our method and characterize the error involved, we compare our results for a carefully constructed test problem to those of a microscale implementation. We demonstrate the effectiveness of our method by simulating a model of polarization and shmoo formation during the mating of yeast. The method is generally applicable to problems in systems biology where biochemistry and mechanics are coupled, and spatial stochastic effects are critical.
Collapse
Affiliation(s)
- Brian Drawert
- Department of Computer Science, University of California-Santa Barbara, Santa Barbara, California 93106, USA
| | - Stefan Hellander
- Department of Computer Science, University of California-Santa Barbara, Santa Barbara, California 93106, USA
| | - Michael Trogdon
- Department of Mechanical Engineering, University of California-Santa Barbara, Santa Barbara, California 93106, USA
| | - Tau-Mu Yi
- Department of Molecular, Cellular, and Developmental Biology, University of California-Santa Barbara, Santa Barbara, California 93106, USA
| | - Linda Petzold
- Department of Computer Science, University of California-Santa Barbara, Santa Barbara, California 93106, USA
| |
Collapse
|
28
|
Sharma R, Roberts E. Gradient sensing by a bistable regulatory motif enhances signal amplification but decreases accuracy in individual cells. Phys Biol 2016; 13:036003. [DOI: 10.1088/1478-3975/13/3/036003] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
|
29
|
Muller N, Piel M, Calvez V, Voituriez R, Gonçalves-Sá J, Guo CL, Jiang X, Murray A, Meunier N. A Predictive Model for Yeast Cell Polarization in Pheromone Gradients. PLoS Comput Biol 2016; 12:e1004795. [PMID: 27077831 PMCID: PMC4831791 DOI: 10.1371/journal.pcbi.1004795] [Citation(s) in RCA: 15] [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: 11/20/2015] [Accepted: 02/08/2016] [Indexed: 11/18/2022] Open
Abstract
Budding yeast cells exist in two mating types, a and α, which use peptide pheromones to communicate with each other during mating. Mating depends on the ability of cells to polarize up pheromone gradients, but cells also respond to spatially uniform fields of pheromone by polarizing along a single axis. We used quantitative measurements of the response of a cells to α-factor to produce a predictive model of yeast polarization towards a pheromone gradient. We found that cells make a sharp transition between budding cycles and mating induced polarization and that they detect pheromone gradients accurately only over a narrow range of pheromone concentrations corresponding to this transition. We fit all the parameters of the mathematical model by using quantitative data on spontaneous polarization in uniform pheromone concentration. Once these parameters have been computed, and without any further fit, our model quantitatively predicts the yeast cell response to pheromone gradient providing an important step toward understanding how cells communicate with each other.
Collapse
Affiliation(s)
- Nicolas Muller
- MAP5, CNRS UMR 8145, Université Paris Descartes, Paris, France
| | - Matthieu Piel
- Institut Curie, CNRS UMR 144, Paris, France
- * E-mail: (MP); (AM); (NM)
| | - Vincent Calvez
- Unité de Mathématiques Pures et Appliquées, CNRS UMR 5669 and équipe-projet INRIA NUMED, École Normale Supérieure de Lyon, Lyon, France
| | - Raphaël Voituriez
- Laboratoire Jean Perrin and Laboratoire de Physique Théorique de la Matière Condensée, UMR 7600 CNRS /UPMC, Paris, France
| | - Joana Gonçalves-Sá
- Molecular and Cell Biology and FAS Center for Systems Biology, Harvard University, Cambridge, Massachusetts, United States of America
| | - Chin-Lin Guo
- Molecular and Cell Biology and FAS Center for Systems Biology, Harvard University, Cambridge, Massachusetts, United States of America
- Institute of Physics, Academia Sinica, Taiwan
| | - Xingyu Jiang
- Department of Chemistry and Chemical Biology, Harvard University, Cambridge, Massachusetts, United States of America
- CAS Key Laboratory for Biological Effects of Nanomaterials and Nanosafety, National Center for Nanoscience and Technology, Beijing, People’s Republic of China
| | - Andrew Murray
- Molecular and Cell Biology and FAS Center for Systems Biology, Harvard University, Cambridge, Massachusetts, United States of America
- * E-mail: (MP); (AM); (NM)
| | - Nicolas Meunier
- MAP5, CNRS UMR 8145, Université Paris Descartes, Paris, France
- * E-mail: (MP); (AM); (NM)
| |
Collapse
|
30
|
Mlynarczyk PJ, Pullen RH, Abel SM. Confinement and diffusion modulate bistability and stochastic switching in a reaction network with positive feedback. J Chem Phys 2016; 144:015102. [DOI: 10.1063/1.4939219] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Paul J. Mlynarczyk
- Department of Chemical and Biomolecular Engineering, University of Tennessee, Knoxville, Tennessee 37996, USA
| | - Robert H. Pullen
- Department of Chemical and Biomolecular Engineering, University of Tennessee, Knoxville, Tennessee 37996, USA
| | - Steven M. Abel
- Department of Chemical and Biomolecular Engineering, University of Tennessee, Knoxville, Tennessee 37996, USA
| |
Collapse
|
31
|
Drawert B, Trogdon M, Toor S, Petzold L, Hellander A. MOLNs: A CLOUD PLATFORM FOR INTERACTIVE, REPRODUCIBLE, AND SCALABLE SPATIAL STOCHASTIC COMPUTATIONAL EXPERIMENTS IN SYSTEMS BIOLOGY USING PyURDME. SIAM JOURNAL ON SCIENTIFIC COMPUTING : A PUBLICATION OF THE SOCIETY FOR INDUSTRIAL AND APPLIED MATHEMATICS 2016; 38:C179-C202. [PMID: 28190948 PMCID: PMC5302863 DOI: 10.1137/15m1014784] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/15/2023]
Abstract
Computational experiments using spatial stochastic simulations have led to important new biological insights, but they require specialized tools and a complex software stack, as well as large and scalable compute and data analysis resources due to the large computational cost associated with Monte Carlo computational workflows. The complexity of setting up and managing a large-scale distributed computation environment to support productive and reproducible modeling can be prohibitive for practitioners in systems biology. This results in a barrier to the adoption of spatial stochastic simulation tools, effectively limiting the type of biological questions addressed by quantitative modeling. In this paper, we present PyURDME, a new, user-friendly spatial modeling and simulation package, and MOLNs, a cloud computing appliance for distributed simulation of stochastic reaction-diffusion models. MOLNs is based on IPython and provides an interactive programming platform for development of sharable and reproducible distributed parallel computational experiments.
Collapse
Affiliation(s)
- Brian Drawert
- Department of Computer Science, University of California, Santa Barbara, Santa Barbara, CA 93106
| | - Michael Trogdon
- Department of Mechanical Engineering, University of California, Santa Barbara, Santa Barbara, CA 93106
| | - Salman Toor
- Department of Computer Science, University of Helsinki, Helsinki FI-00014, Finland, and Department of Information Technology, Division of Scientific Computing, Uppsala University, Uppsala, 75105 Sweden
| | - Linda Petzold
- Departments of Computer Science and Mechanical Engineering, University of California, Santa Barbara, Santa Barbara, CA 93106
| | - Andreas Hellander
- Corresponding author. Department of Information Technology, Division of Scientific Computing, Uppsala University, Uppsala, 75105 Sweden ()
| |
Collapse
|
32
|
Hegemann B, Unger M, Lee SS, Stoffel-Studer I, van den Heuvel J, Pelet S, Koeppl H, Peter M. A Cellular System for Spatial Signal Decoding in Chemical Gradients. Dev Cell 2015; 35:458-70. [PMID: 26585298 DOI: 10.1016/j.devcel.2015.10.013] [Citation(s) in RCA: 33] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/25/2015] [Revised: 07/29/2015] [Accepted: 10/18/2015] [Indexed: 12/21/2022]
Abstract
Directional cell growth requires that cells read and interpret shallow chemical gradients, but how the gradient directional information is identified remains elusive. We use single-cell analysis and mathematical modeling to define the cellular gradient decoding network in yeast. Our results demonstrate that the spatial information of the gradient signal is read locally within the polarity site complex using double-positive feedback between the GTPase Cdc42 and trafficking of the receptor Ste2. Spatial decoding critically depends on low Cdc42 activity, which is maintained by the MAPK Fus3 through sequestration of the Cdc42 activator Cdc24. Deregulated Cdc42 or Ste2 trafficking prevents gradient decoding and leads to mis-oriented growth. Our work discovers how a conserved set of components assembles a network integrating signal intensity and directionality to decode the spatial information contained in chemical gradients.
Collapse
Affiliation(s)
- Björn Hegemann
- Department of Biology, Institute of Biochemistry, ETH Zurich, Otto-Stern-Weg 3, 8093 Zürich, Switzerland.
| | - Michael Unger
- Department of Biology, Institute of Biochemistry, ETH Zurich, Otto-Stern-Weg 3, 8093 Zürich, Switzerland; Automatic Control Laboratory, ETH Zurich, Physikstrasse 3, 8092 Zürich, Switzerland
| | - Sung Sik Lee
- Department of Biology, Institute of Biochemistry, ETH Zurich, Otto-Stern-Weg 3, 8093 Zürich, Switzerland
| | - Ingrid Stoffel-Studer
- Department of Biology, Institute of Biochemistry, ETH Zurich, Otto-Stern-Weg 3, 8093 Zürich, Switzerland
| | - Jasmin van den Heuvel
- Department of Biology, Institute of Biochemistry, ETH Zurich, Otto-Stern-Weg 3, 8093 Zürich, Switzerland
| | - Serge Pelet
- Department of Fundamental Microbiology, University of Lausanne, Biophore Building, 1015 Lausanne, Switzerland
| | - Heinz Koeppl
- Automatic Control Laboratory, ETH Zurich, Physikstrasse 3, 8092 Zürich, Switzerland; Department of Electrical Engineering and Information Technology, Technische Universität Darmstadt, Rundeturmstrasse 12, 64283 Darmstadt, Germany
| | - Matthias Peter
- Department of Biology, Institute of Biochemistry, ETH Zurich, Otto-Stern-Weg 3, 8093 Zürich, Switzerland.
| |
Collapse
|
33
|
Wehrens M, ten Wolde PR, Mugler A. Positive feedback can lead to dynamic nanometer-scale clustering on cell membranes. J Chem Phys 2015; 141:205102. [PMID: 25429963 DOI: 10.1063/1.4901888] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022] Open
Abstract
Clustering of molecules on biological membranes is a widely observed phenomenon. A key example is the clustering of the oncoprotein Ras, which is known to be important for signal transduction in mammalian cells. Yet, the mechanism by which Ras clusters form and are maintained remains unclear. Recently, it has been discovered that activated Ras promotes further Ras activation. Here we show using particle-based simulation that this positive feedback is sufficient to produce persistent clusters of active Ras molecules at the nanometer scale via a dynamic nucleation mechanism. Furthermore, we find that our cluster statistics are consistent with experimental observations of the Ras system. Interestingly, we show that our model does not support a Turing regime of macroscopic reaction-diffusion patterning, and therefore that the clustering we observe is a purely stochastic effect, arising from the coupling of positive feedback with the discrete nature of individual molecules. These results underscore the importance of stochastic and dynamic properties of reaction diffusion systems for biological behavior.
Collapse
Affiliation(s)
- Martijn Wehrens
- FOM Institute AMOLF, Science Park 104, 1098 XG Amsterdam, The Netherlands
| | | | - Andrew Mugler
- FOM Institute AMOLF, Science Park 104, 1098 XG Amsterdam, The Netherlands
| |
Collapse
|
34
|
Razooky BS, Pai A, Aull K, Rouzine IM, Weinberger LS. A hardwired HIV latency program. Cell 2015; 160:990-1001. [PMID: 25723172 DOI: 10.1016/j.cell.2015.02.009] [Citation(s) in RCA: 171] [Impact Index Per Article: 19.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/02/2014] [Revised: 12/09/2014] [Accepted: 02/05/2015] [Indexed: 12/11/2022]
Abstract
Biological circuits can be controlled by two general schemes: environmental sensing or autonomous programs. For viruses such as HIV, the prevailing hypothesis is that latent infection is controlled by cellular state (i.e., environment), with latency simply an epiphenomenon of infected cells transitioning from an activated to resting state. However, we find that HIV expression persists despite the activated-to-resting cellular transition. Mathematical modeling indicates that HIV's Tat positive-feedback circuitry enables this persistence and strongly controls latency. To overcome the inherent crosstalk between viral circuitry and cellular activation and to directly test this hypothesis, we synthetically decouple viral dependence on cellular environment from viral transcription. These circuits enable control of viral transcription without cellular activation and show that Tat feedback is sufficient to regulate latency independent of cellular activation. Overall, synthetic reconstruction demonstrates that a largely autonomous, viral-encoded program underlies HIV latency—potentially explaining why cell-targeted latency-reversing agents exhibit incomplete penetrance.
Collapse
Affiliation(s)
- Brandon S Razooky
- The Gladstone Institutes (Virology and Immunology), San Francisco, California Institute of Quantitative Biosciences, University of California, San Francisco, 94158; Biophysics Graduate Group, California Institute of Quantitative Biosciences, University of California, San Francisco, 94158
| | - Anand Pai
- The Gladstone Institutes (Virology and Immunology), San Francisco, California Institute of Quantitative Biosciences, University of California, San Francisco, 94158; Department of Biochemistry and Biophysics, California Institute of Quantitative Biosciences, University of California, San Francisco, 94158
| | - Katherine Aull
- Biophysics Graduate Group, California Institute of Quantitative Biosciences, University of California, San Francisco, 94158
| | - Igor M Rouzine
- The Gladstone Institutes (Virology and Immunology), San Francisco, California Institute of Quantitative Biosciences, University of California, San Francisco, 94158
| | - Leor S Weinberger
- The Gladstone Institutes (Virology and Immunology), San Francisco, California Institute of Quantitative Biosciences, University of California, San Francisco, 94158; Department of Biochemistry and Biophysics, California Institute of Quantitative Biosciences, University of California, San Francisco, 94158; QB3, California Institute of Quantitative Biosciences, University of California, San Francisco, 94158.
| |
Collapse
|
35
|
Saunders TE. Aggregation-fragmentation model of robust concentration gradient formation. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2015; 91:022704. [PMID: 25768528 DOI: 10.1103/physreve.91.022704] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/28/2014] [Indexed: 06/04/2023]
Abstract
Concentration gradients of signaling molecules are essential for patterning during development and they have been observed in both unicellular and multicellular systems. In subcellular systems, clustering of the signaling molecule has been observed. We develop a theoretical model of cluster-mediated concentration gradient formation based on the Becker-Döring equations of aggregation-fragmentation processes. We show that such a mechanism produces robust concentration gradients on realistic time and spatial scales so long as the process of clustering does not significantly stabilize the signaling molecule. Finally, we demonstrate that such a model is applicable to the pom1p subcellular gradient in fission yeast.
Collapse
Affiliation(s)
- Timothy E Saunders
- Mechanobiology Institute, National University of Singapore, Singapore; Department of Biological Sciences, National University of Singapore, Singapore; and Institute of Molecular and Cell Biology, Proteos, Singapore
| |
Collapse
|
36
|
Krishnan J, Mois K, Suwanmajo T. The behaviour of basic autocatalytic signalling modules in isolation and embedded in networks. J Chem Phys 2014; 141:175102. [DOI: 10.1063/1.4898370] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
37
|
Dodgson J, Chessel A, Yamamoto M, Vaggi F, Cox S, Rosten E, Albrecht D, Geymonat M, Csikasz-Nagy A, Sato M, Carazo-Salas RE. Spatial segregation of polarity factors into distinct cortical clusters is required for cell polarity control. Nat Commun 2013; 4:1834. [PMID: 23673619 PMCID: PMC3674234 DOI: 10.1038/ncomms2813] [Citation(s) in RCA: 43] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/03/2012] [Accepted: 03/26/2013] [Indexed: 01/06/2023] Open
Abstract
Cell polarity is regulated by evolutionarily conserved polarity factors whose precise higher-order organization at the cell cortex is largely unknown. Here we image frontally the cortex of live fission yeast cells using time-lapse and super-resolution microscopy. Interestingly, we find that polarity factors are organized in discrete cortical clusters resolvable to ~50–100 nm in size, which can form and become cortically enriched by oligomerization. We show that forced co-localization of the polarity factors Tea1 and Tea3 results in polarity defects, suggesting that the maintenance of both factors in distinct clusters is required for polarity. However, during mitosis, their co-localization increases, and Tea3 helps to retain the cortical localization of the Tea1 growth landmark in preparation for growth reactivation following mitosis. Thus, regulated spatial segregation of polarity factor clusters provides a means to spatio-temporally control cell polarity at the cell cortex. We observe similar clusters in Saccharomyces cerevisiae and Caenorhabditis elegans cells, indicating this could be a universal regulatory feature. Cell polarity is generated and maintained by the spatial accumulation of polarity factors. By imaging fission yeast cells ‘end-on’, the authors show that the polarity factors Tea1 and Tea3 segregate into distinct clusters, and that surprisingly, their segregation is critical for cell polarization.
Collapse
Affiliation(s)
- James Dodgson
- The Gurdon Institute, University of Cambridge, Cambridge CB2 1QN, UK
| | | | | | | | | | | | | | | | | | | | | |
Collapse
|
38
|
Freisinger T, Klünder B, Johnson J, Müller N, Pichler G, Beck G, Costanzo M, Boone C, Cerione RA, Frey E, Wedlich-Söldner R. Establishment of a robust single axis of cell polarity by coupling multiple positive feedback loops. Nat Commun 2013; 4:1807. [PMID: 23651995 PMCID: PMC3674238 DOI: 10.1038/ncomms2795] [Citation(s) in RCA: 75] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/13/2012] [Accepted: 03/22/2013] [Indexed: 01/06/2023] Open
Abstract
Establishment of cell polarity—or symmetry breaking—relies on local accumulation of polarity regulators. Although simple positive feedback is sufficient to drive symmetry breaking, it is highly sensitive to stochastic fluctuations typical for living cells. Here, by integrating mathematical modelling with quantitative experimental validations, we show that in the yeast Saccharomyces cerevisiae a combination of actin- and guanine nucleotide dissociation inhibitor-dependent recycling of the central polarity regulator Cdc42 is needed to establish robust cell polarity at a single site during yeast budding. The guanine nucleotide dissociation inhibitor pathway consistently generates a single-polarization site, but requires Cdc42 to cycle rapidly between its active and inactive form, and is therefore sensitive to perturbations of the GTPase cycle. Conversely, actin-mediated recycling of Cdc42 induces robust symmetry breaking but cannot restrict polarization to a single site. Our results demonstrate how cells optimize symmetry breaking through coupling between multiple feedback loops. A positive feedback loop which results in localized accumulation of the small GTPase Cdc42 generates cell polarity in budding yeast; however, such loops are inherently susceptible to noise. Here the authors demonstrate how two pathways that mediate Cdc42 recycling work together to ensure the robustness of symmetry breaking.
Collapse
Affiliation(s)
- Tina Freisinger
- Max-Planck-Institute of Biochemistry, Cellular Dynamics and Cell Patterning, Am Klopferspitz 18, D-82152 Martinsried, Germany
| | | | | | | | | | | | | | | | | | | | | |
Collapse
|
39
|
Klann M, Koeppl H. Reaction schemes, escape times and geminate recombinations in particle-based spatial simulations of biochemical reactions. Phys Biol 2013; 10:046005. [DOI: 10.1088/1478-3975/10/4/046005] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
|
40
|
Kochanczyk M, Jaruszewicz J, Lipniacki T. Stochastic transitions in a bistable reaction system on the membrane. J R Soc Interface 2013; 10:20130151. [PMID: 23635492 DOI: 10.1098/rsif.2013.0151] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022] Open
Abstract
Transitions between steady states of a multi-stable stochastic system in the perfectly mixed chemical reactor are possible only because of stochastic switching. In realistic cellular conditions, where diffusion is limited, transitions between steady states can also follow from the propagation of travelling waves. Here, we study the interplay between the two modes of transition for a prototype bistable system of kinase-phosphatase interactions on the plasma membrane. Within microscopic kinetic Monte Carlo simulations on the hexagonal lattice, we observed that for finite diffusion the behaviour of the spatially extended system differs qualitatively from the behaviour of the same system in the well-mixed regime. Even when a small isolated subcompartment remains mostly inactive, the chemical travelling wave may propagate, leading to the activation of a larger compartment. The activating wave can be induced after a small subdomain is activated as a result of a stochastic fluctuation. Such a spontaneous onset of activity is radically more probable in subdomains characterized by slower diffusion. Our results show that a local immobilization of substrates can lead to the global activation of membrane proteins by the mechanism that involves stochastic fluctuations followed by the propagation of a semi-deterministic travelling wave.
Collapse
Affiliation(s)
- Marek Kochanczyk
- Institute of Fundamental Technological Research, Polish Academy of Sciences, Warsaw, Poland
| | | | | |
Collapse
|
41
|
Klann M, Koeppl H. Spatial simulations in systems biology: from molecules to cells. Int J Mol Sci 2012; 13:7798-7827. [PMID: 22837728 PMCID: PMC3397560 DOI: 10.3390/ijms13067798] [Citation(s) in RCA: 36] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/02/2012] [Revised: 06/08/2012] [Accepted: 06/12/2012] [Indexed: 12/23/2022] Open
Abstract
Cells are highly organized objects containing millions of molecules. Each biomolecule has a specific shape in order to interact with others in the complex machinery. Spatial dynamics emerge in this system on length and time scales which can not yet be modeled with full atomic detail. This review gives an overview of methods which can be used to simulate the complete cell at least with molecular detail, especially Brownian dynamics simulations. Such simulations require correct implementation of the diffusion-controlled reaction scheme occurring on this level. Implementations and applications of spatial simulations are presented, and finally it is discussed how the atomic level can be included for instance in multi-scale simulation methods.
Collapse
Affiliation(s)
- Michael Klann
- Authors to whom correspondence should be addressed; E-Mails: (M.K.); (H.K.); Tel.: +41-44-632-4274 (M.K.); +41-44-632-7288 (H.K.); Fax: +41-44-632-1211 (M.K.; H.K.)
| | - Heinz Koeppl
- Authors to whom correspondence should be addressed; E-Mails: (M.K.); (H.K.); Tel.: +41-44-632-4274 (M.K.); +41-44-632-7288 (H.K.); Fax: +41-44-632-1211 (M.K.; H.K.)
| |
Collapse
|
42
|
Houk AR, Jilkine A, Mejean CO, Boltyanskiy R, Dufresne ER, Angenent SB, Altschuler SJ, Wu LF, Weiner OD. Membrane tension maintains cell polarity by confining signals to the leading edge during neutrophil migration. Cell 2012; 148:175-88. [PMID: 22265410 PMCID: PMC3308728 DOI: 10.1016/j.cell.2011.10.050] [Citation(s) in RCA: 398] [Impact Index Per Article: 33.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/17/2011] [Revised: 09/15/2011] [Accepted: 10/24/2011] [Indexed: 12/24/2022]
Abstract
Little is known about how neutrophils and other cells establish a single zone of actin assembly during migration. A widespread assumption is that the leading edge prevents formation of additional fronts by generating long-range diffusible inhibitors or by sequestering essential polarity components. We use morphological perturbations, cell-severing experiments, and computational simulations to show that diffusion-based mechanisms are not sufficient for long-range inhibition by the pseudopod. Instead, plasma membrane tension could serve as a long-range inhibitor in neutrophils. We find that membrane tension doubles during leading-edge protrusion, and increasing tension is sufficient for long-range inhibition of actin assembly and Rac activation. Furthermore, reducing membrane tension causes uniform actin assembly. We suggest that tension, rather than diffusible molecules generated or sequestered at the leading edge, is the dominant source of long-range inhibition that constrains the spread of the existing front and prevents the formation of secondary fronts.
Collapse
Affiliation(s)
- Andrew R. Houk
- Cardiovascular Research Institute and Department of Biochemistry, University of California San Francisco, San Francisco, California
| | - Alexandra Jilkine
- Green Comprehensive Center for Molecular, Computational, and Systems Biology, Southwestern Medical School, Dallas, Texas
| | | | | | | | | | - Steven J. Altschuler
- Green Comprehensive Center for Molecular, Computational, and Systems Biology, Southwestern Medical School, Dallas, Texas
| | - Lani F. Wu
- Green Comprehensive Center for Molecular, Computational, and Systems Biology, Southwestern Medical School, Dallas, Texas
| | - Orion D. Weiner
- Cardiovascular Research Institute and Department of Biochemistry, University of California San Francisco, San Francisco, California
| |
Collapse
|