Please wait a minute...
Quantitative Biology

ISSN 2095-4689

ISSN 2095-4697(Online)

CN 10-1028/TM

Postal Subscription Code 80-971

Quant. Biol.    2015, Vol. 3 Issue (2) : 55-68    https://doi.org/10.1007/s40484-015-0043-0
RESEARCH ARTICLE
Mathematical modeling reveals the mechanisms of feedforward regulation in cell fate decisions in budding yeast
Wenlong Li1,Ming Yi2,3,*(),Xiufen Zou1,*()
1. School of Mathematics and Statistics, Wuhan University, Wuhan 430072, China
2. Key Laboratory of Magnetic Resonance in Biological Systems, Wuhan Institute of Physics and Mathematics, Chinese Academy of Sciences, Wuhan 430071, China
3. Department of Physics, College of Science, Huazhong Agricultural University, Wuhan 430070, China
 Download: PDF(1870 KB)   HTML
 Export: BibTeX | EndNote | Reference Manager | ProCite | RefWorks
Abstract

The determination of cell fate is one of the key questions of developmental biology. Recent experiments showed that feedforward regulation is a novel feature of regulatory networks that controls reversible cellular transitions. However, the underlying mechanism of feedforward regulation-mediated cell fate decision is still unclear. Therefore, using experimental data, we develop a full mathematical model of the molecular network responsible for cell fate selection in budding yeast. To validate our theoretical model, we first investigate the dynamical behaviors of key proteins at the Start transition point and the G1/S transition point; a crucial three-node motif consisting of cyclin (Cln1/2), Substrate/Subunit Inhibitor of cyclin-dependent protein kinase (Sic1) and cyclin B (Clb5/6) is considered at these points. The rapid switches of these important components between high and low levels at two transition check points are demonstrated reasonably by our model. Many experimental observations about cell fate decision and cell size control are also theoretically reproduced. Interestingly, the feedforward regulation provides a reliable separation between different cell fates. Next, our model reveals that the threshold for the amount of WHIskey (Whi5) removed from the nucleus is higher at the Reentry point in pheromone-arrested cells compared with that at the Start point in cycling cells. Furthermore, we analyze the hysteresis in the cell cycle kinetics in response to changes in pheromone concentration, showing that Cln3 is the primary driver of reentry and Cln1/2 is the secondary driver of reentry. In particular, we demonstrate that the inhibition of Cln1/2 due to the accumulation of Factor ARrest (Far1) directly reinforces arrest. Finally, theoretical work verifies that the three-node coherent feedforward motif created by cell FUSion (Fus3), Far1 and STErile (Ste12) ensures the rapid arrest and reversibility of a cellular state. The combination of our theoretical model and the previous experimental data contributes to the understanding of the molecular mechanisms of the cell fate decision at the G1 phase in budding yeast and will stimulate further biological experiments in future.

Keywords cell fate decision      feedforward mechanism      mathematical modeling      hysteresis      reversibility     
Corresponding Author(s): Ming Yi,Xiufen Zou   
Just Accepted Date: 01 April 2015   Online First Date: 28 April 2015    Issue Date: 21 August 2015
 Cite this article:   
Wenlong Li,Ming Yi,Xiufen Zou. Mathematical modeling reveals the mechanisms of feedforward regulation in cell fate decisions in budding yeast[J]. Quant. Biol., 2015, 3(2): 55-68.
 URL:  
https://academic.hep.com.cn/qb/EN/10.1007/s40484-015-0043-0
https://academic.hep.com.cn/qb/EN/Y2015/V3/I2/55
Fig.1  Schematic illustration of two cell network subsystems. All components and reactions considered in our mathematical model are included. All arrows for individual reactions are marked with the corresponding number of this reaction in the model. (A) Cell cycle subsystem of budding yeast in G1 phase. (B) Pheromone-induced MAPK pathway subsystem of budding yeast.
Fig.2  Numerical experiments for cell size control in budding yeast. (A) The definition of the Start point and the G1/S transition point. The Start point occurs when 50% of Whi5 is removed from the nucleus, and the crossing point of the Clb5/6free time course and the Sic1 time course is considered the G1/S transition point. (B) The relationship between the duration of cell cycle and the birth size.
Fig.3  Relationship between the time of pheromone addition and Whi5P activation for pheromone-arrested cells. At the Reentry point, the maximum slope is achieved in this curve of cell fate decision. In pheromone-arrested cells, the Reentry point is defined by the export of approximately 58% of the nuclear Whi5.
Fig.4  Hysteresis in the cell cycle kinetics in response to changes in pheromone concentration. (A) Schematic of the hysteresis experiment. Experimental group cells are exposed to 240 nM α-factor for 30 min (pheromone-arrested cells), and control group cells are not treated with anything (cycling cells). Then, the two groups of cells are exposed to different α-factor concentrations (different conditions). (B) Duration of arrest in pheromone-arrested cells and cycling cells in different conditions. (C) Hysteresis experiment for cells lacking the gene cln1/2. (D) Arrest duration of pheromone-arrested cells lacking the gene cln1/2, pheromone-arrested cells lacking the gene cln3 and pheromone-arrested WT cells in different conditions.
Fig.5  Dynamical analysis for the feedforward motif with slow transcription and fast phosphorylation. (A) Time courses of Fus3PP, synthesized Far1 (i.e., Far1T) and phosphorylated Far1 (i.e., Far1PP) in WT cells pre-exposed to 240 nM pheromone. (B) Phase diagram of the arrest duration and the Far1T abundance at the Reentry point for WT cells, cln1/2Δ cells, and cln3Δ cells for changes in the final pheromone concentration in the hysteresis experiment.
Fig.6  Whi5-threshold and arrest duration are enhanced by increases in the pre-exposed duration. (A) The relationship between the pre-exposed duration and Whi5-threshold. (B) The relationship of the arrest duration and the critical addition time with pre-exposed duration.
Fig.7  . Comparison between the feedforward motif and other structures. (A) Here, we consider four categories of networks that control the level of active Far1. (B) A general response for the networks to a square pulse of the mating pheromone. (C) The steady state relationship between MAPK activity and the mating pheromone. (D). Relationship between the time of pheromone addition and Whi5P activation. (E) Hysteresis experiment for different Far1 regulatory mechanisms. Only the feedforward regulation fulfills all requirements (rapid activation and reentry, ability to continually sense extracellular information, reliable separation of cell fates, and hysteresis effect).
1 Furlong, E. E. (2010) The importance of being specified: cell fate decisions and their role in cell biology. Mol. Biol. Cell, 21, 3797–3798
https://doi.org/10.1091/mbc.E10-05-0436. pmid: 21079016
2 Morgan, D. O. (2007) The cell cycle. London: New Science Press
3 Yang, L., MacLellan, W. R., Han, Z., Weiss, J. N. and Qu, Z. (2004) Multisite phosphorylation and network dynamics of cyclin-dependent kinase signaling in the eukaryotic cell cycle. Biophys. J., 86, 3432–3443
https://doi.org/10.1529/biophysj.103.036558. pmid: 15189845
4 Bardwell, L., Zou, X., Nie, Q. and Komarova, N. L. (2007) Mathematical models of specificity in cell signaling. Biophys. J., 92, 3425–3441
https://doi.org/10.1529/biophysj.106.090084. pmid: 17325015
5 Zou, X., Peng, T. and Pan, Z. (2008) Modeling specificity in the yeast MAPK signaling networks. J. Theor. Biol., 250, 139–155
https://doi.org/10.1016/j.jtbi.2007.09.024. pmid: 17977559
6 Kofahl, B. and Klipp, E. (2004) Modelling the dynamics of the yeast pheromone pathway. Yeast, 21, 831–850
https://doi.org/10.1002/yea.1122 pmid: 15300679
7 Li, Y., Yi, M. and Zou, X. (2013) Identification of the molecular mechanisms for cell-fate selection in budding yeast through mathematical modeling. Biophys. J., 104, 2282–2294
https://doi.org/10.1016/j.bpj.2013.03.057. pmid: 23708368
8 Chang, F. and Herskowitz, I. (1990) Identification of a gene necessary for cell cycle arrest by a negative growth factor of yeast: FAR1 is an inhibitor of a G1 cyclin, CLN2. Cell, 63, 999–1011
https://doi.org/10.1016/0092-8674(90)90503-7. pmid: 2147873
9 McKinney, J. D., Chang, F., Heintz, N. and Cross, F. R. (1993) Negative regulation of FAR1 at the Start of the yeast cell cycle. Genes Dev., 7, 833–843
https://doi.org/10.1101/gad.7.5.833. pmid: 8491380
10 Yang, X., Lau, K. Y., Sevim, V. and Tang, C. (2013) Design principles of the yeast G1/S switch. PLoS Biol., 11, e1001673
https://doi.org/10.1371/journal.pbio.1001673. pmid: 24130459
11 Mangan, S. and Alon, U. (2003) Structure and function of the feed-forward loop network motif. Proc. Natl. Acad. Sci. USA, 100, 11980–11985
https://doi.org/10.1073/pnas.2133841100 pmid: 14530388
12 Doncic, A. and Skotheim, J. M. (2013) Feedforward regulation ensures stability and rapid reversibility of a cellular state. Mol. Cell, 50, 856–868
https://doi.org/10.1016/j.molcel.2013.04.014. pmid: 23685071
13 Doncic, A., Falleur-Fettig, M. and Skotheim, J. M. (2011) Distinct interactions select and maintain a specific cell fate. Mol. Cell, 43, 528–539
https://doi.org/10.1016/j.molcel.2011.06.025 pmid: 21855793
14 Novak, B., Tyson, J. J., Gyorffy, B. and Csikasz-Nagy, A. (2007) Irreversible cell-cycle transitions are due to systems-level feedback. Nat. Cell Biol., 9, 724–728
https://doi.org/10.1038/ncb0707-724 pmid: 17603504
15 Skotheim, J. M., Di Talia, S., Siggia, E. D. and Cross, F. R. (2008) Positive feedback of G1 cyclins ensures coherent cell cycle entry. Nature, 454, 291–296
https://doi.org/10.1038/nature07118. pmid: 18633409
16 Charvin, G., Oikonomou, C., Siggia, E. D. and Cross, F. R. (2010) Origin of irreversibility of cell cycle start in budding yeast. PLoS Biol., 8, e1000284
https://doi.org/10.1371/journal.pbio.1000284 pmid: 20087409
17 Hartwell, L. H., Culotti, J., Pringle, J. R. and Reid, B. J. (1974) Genetic control of the cell division cycle in yeast. Science, 183, 46–51
https://doi.org/10.1126/science.183.4120.46. pmid: 4587263
18 Cross, F. R. (1995) Starting the cell cycle: what’s the point? Curr. Opin. Cell Biol., 7, 790–797
https://doi.org/10.1016/0955-0674(95)80062-X. pmid: 8608009
19 Talia, S. D., Skotheim, J. M., Bean, J. M., Siggia, E. D. and Cross, F. R. (2007) The effects of molecular noise and size control on variability in the budding yeast cell cycle. Nature, 448, 947–951
https://doi.org/10.1038/nature06072. pmid: 17713537
20 Di Talia, S., Wang, H., Skotheim, J. M., Rosebrock, A. P., Futcher, B. and Cross, F. R. (2009) Daughter-specific transcription factors regulate cell size control in budding yeast. PLoS Biol., 7, e1000221
https://doi.org/10.1371/journal.pbio.1000221. pmid: 19841732
21 Turner, J. J., Ewald, J. C. and Skotheim, J. M. (2012) Cell size control in yeast. Curr. Biol., 22, R350–R359
https://doi.org/10.1016/j.cub.2012.02.041 pmid: 22575477
22 Hao, N., Nayak, S., Behar, M., Shanks, R. H., Nagiec, M. J., Errede, B., Hasty, J., Elston, T. C. and Dohlman, H. G. (2008) Regulation of cell signaling dynamics by the protein kinase-scaffold Ste5. Mol. Cell, 30, 649–656
https://doi.org/10.1016/j.molcel.2008.04.016. pmid: 18538663
23 Tyers, M. and Futcher, B. (1993) Far1 and Fus3 link the mating pheromone signal transduction pathway to three G1-phase Cdc28 kinase complexes. Mol. Cell. Biol., 13, 5659–5669
pmid: 8395009
24 Takahashi, S. and Pryciak, P. M. (2008) Membrane localization of scaffold proteins promotes graded signaling in the yeast MAP kinase cascade. Curr. Biol., 18, 1184–1191
https://doi.org/10.1016/j.cub.2008.07.050. pmid: 18722124
25 Hartwell, L. H. and Unger, M. W. (1977) Unequal division in Saccharomyces cerevisiae and its implications for the control of cell division. J. Cell Biol., 75, 422–435
https://doi.org/10.1083/jcb.75.2.422. pmid: 400873
26 Qu, Z., MacLellan, W. R. and Weiss, J. N. (2003) Dynamics of the cell cycle: checkpoints, sizers, and timers. Biophys. J., 85, 3600–3611
https://doi.org/10.1016/S0006-3495(03)74778-X. pmid: 14645053
27 Colman-Lerner, A., Gordon, A., Serra, E., Chin, T., Orna, R., Endy, D., Pesce, C. G. and Brent, R. (2005) Regulated cell-to-cell variation in a cell-fate decision system. Nature, 437, 699–706
https://doi.org/10.1038/nature03998. pmid: 16170311
28 Paliwal, S., Iglesias, P. A., Campbell, K., Hilioti, Z., Groisman, A. and Levchenko, A. (2007) MAPK-mediated bimodal gene expression and adaptive gradient sensing in yeast. Nature, 446, 46–51
https://doi.org/10.1038/nature05561. pmid: 17310144
29 Yu, R. C., Pesce, C. G., Colman-Lerner, A., Lok, L., Pincus, D., Serra, E., Holl, M., Benjamin, K., Gordon, A. and Brent, R. (2008) Negative feedback that improves information transmission in yeast signalling. Nature, 456, 755–761
https://doi.org/10.1038/nature07513. pmid: 19079053
30 Ubersax, J. A. and Ferrell, J. E. (2006) A noisy ‘Start’ to the cell cycle. Mol. Syst. Biol., 2, 0014
https://doi.org/10.1038/msb4100056. pmid: 16738559
31 Kar, S., Baumann, W. T., Paul, M. R. and Tyson, J. J. (2009) Exploring the roles of noise in the eukaryotic cell cycle. Proc. Natl. Acad. Sci. USA, 106, 6471–6476
https://doi.org/10.1073/pnas.0810034106. pmid: 19246388
32 Li Y. K., Yi M., Zou X. F. (2014) The linear interplay of intrinsic and extrinsic noises ensures a high accuracy of cell fate selection in budding yeast. Sci. Rep, 4, 5764
33 Zou, X., Xiang, X., Chen, Y., Peng, T., Luo, X. and Pan, Z. (2010) Understanding inhibition of viral proteins on type I IFN signaling pathways with modeling and optimization. J. Theor. Biol., 265, 691–703
https://doi.org/10.1016/j.jtbi.2010.05.001. pmid: 20553733
34 Sun, J., Yi, M., Yang, L., Wei, W., Ding, Y. and Jia, Y. (2014) Enhancement of tunability of MAPK cascade due to coexistence of processive and distributive phosphorylation mechanisms. Biophys. J., 106, 1215–1226
https://doi.org/10.1016/j.bpj.2014.01.036. pmid: 24606945
[1] Supplementary Material Download
[1] Naoki Matsuda, Ken-ichi Hironaka, Masashi Fujii, Takumi Wada, Katsuyuki Kunida, Haruki Inoue, Miki Eto, Daisuke Hoshino, Yasuro Furuichi, Yasuko Manabe, Nobuharu L. Fujii, Hiroyuki Noji, Hiromi Imamura, Shinya Kuroda. Monitoring and mathematical modeling of mitochondrial ATP in myotubes at single-cell level reveals two distinct population with different kinetics[J]. Quant. Biol., 2020, 8(3): 228-237.
[2] Petr Kloucek, Armin von Gunten. On the possibility of identifying human subjects using behavioural complexity analyses[J]. Quant. Biol., 2016, 4(4): 261-269.
[3] Derek Eidum,Kanishk Asthana,Samir Unni,Michael Deng,Lingchong You. Construction, visualization, and analysis of biological network models in Dynetica[J]. Quant. Biol., 2014, 2(4): 142-150.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed