# A unifying model for the propagation of prion proteins in yeast brings insight into the [PSI+] prion

Paul Lemarre, Laurent Pujo-Menjouet, Suzanne S Sindi

► To cite this version:

Paul Lemarre, Laurent Pujo-Menjouet, Suzanne S Sindi. A unifying model for the propagation of prion proteins in yeast brings insight into the [PSI+] prion. PLoS Computational Biology, 2020, 16 (5), pp.e1007647. 10.1371/journal.pcbi.1007647. hal-02732998

HAL Id: hal-02732998

https://hal.science/hal-02732998v1

Submitted on 14 Mar 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

RESEARCH ARTICLE

# A unifying model for the propagation of prion proteins in yeast brings insight into the [$PSI^+$] prion

Paul Lemanre¹,²,³, Laurent Pujo-Menjouet¹,²✉, Suzanne S. Sindi³✉*

1 Univ Lyon, Université Claude Bernard Lyon 1, CNRS UMR 5208, Institut Camille Jordan, 43 blvd. du 11 novembre 1918, F-69622 Villeurbanne cedex, France, 2 INRIA Rhônes-Alpes, INRIA, Villeurbanne, France,

3 Department of Applied Mathematics, University of California Merced, Merced, California, United States of America

✉ These authors contributed equally to this work.

* ssindi@ucmerced.edu

## Abstract

The use of yeast systems to study the propagation of prions and amyloids has emerged as a crucial aspect of the global endeavor to understand those mechanisms. Yeast prion systems are intrinsically multi-scale: the molecular chemical processes are indeed coupled to the cellular processes of cell growth and division to influence phenotypical traits, observable at the scale of colonies. We introduce a novel modeling framework to tackle this difficulty using impulsive differential equations. We apply this approach to the [$PSI^+$] yeast prion, which is associated with the misconformation and aggregation of Sup35. We build a model that reproduces and unifies previously conflicting experimental observations on [$PSI^+$] and thus sheds light onto characteristics of the intracellular molecular processes driving aggregate replication. In particular our model uncovers a kinetic barrier for aggregate replication at low densities, meaning the change between prion or prion-free phenotype is a bi-stable transition. This result is based on the study of prion curing experiments, as well as the phenomenon of colony sectoring, a phenotype which is often ignored in experimental assays and has never been modeled. Furthermore, our results provide further insight into the effect of guanidine hydrochloride (GdnHCl) on Sup35 aggregates. To qualitatively reproduce the GdnHCl curing experiment, aggregate replication must not be completely inhibited, which suggests the existence of a mechanism different than Hsp104-mediated fragmentation. Those results are promising for further development of the [$PSI^+$] model, but also for extending the use of this novel framework to other yeast prion or amyloid systems.

## Author summary

In the study of yeast prions, mathematical modeling is a powerful tool, in particular when it comes to facing the difficulties of multi-scale systems. In this study, we introduce a mathematical framework for investigating this problem in a unifying way. We focus on the yeast prion [$PSI^+$] and present a simple molecular scheme for prion replication and a

<figure><img src="image_3.png" /></figure>

## OPEN ACCESS

**Citation:** Lemanre P, Pujo-Menjouet L, Sindi SS (2020) A unifying model for the propagation of prion proteins in yeast brings insight into the [$PSI^+$] prion. PLoS Comput Biol 16(5): e1007647. https://doi.org/10.1371/journal.pcbi.1007647

**Editor:** Philip K Maini, Oxford, UNITED KINGDOM

**Received:** January 7, 2020

**Accepted:** March 6, 2020

**Published:** May 26, 2020

**Copyright:** © 2020 Lemanre et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** All relevant data are within the manuscript and its Supporting Information files.

**Funding:** This work is supported by the Joint DMS/NIGMS Initiative to Support Research at the Interface of the Biological and Mathematical Sciences between the National Science Foundation (NSF) and the National Institutes of Health NIH (R01-GM126548, https://www.nsf.gov/funding/pgm_summ.jsp?pims_id=5300), awarded to SS. This work was supported by the LABEX MILYON (ANR-10-LABX-0070) of Université de Lyon, within the program "Investissements d'Avenir" (ANR-11-
IDEX-0007) operated by the French National Research Agency (ANR), awarded to LPM and PL. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

**Competing interests:** The authors have declared that no competing interests exist.

model of yeast budding. In order to qualitatively reproduce experiments, we need to introduce a non-linear mechanism in the molecular rates. This transforms the intracellular system into a bi-stable switch and allows for curing to occur, which is a crucial phenomenon for the study of yeast prions. To the best of our knowledge, no model in the literature includes such a mechanism, at least not explicitly. We also describe the GdnHCl curing experiment, and the propagon counting procedure. Reproducing this result requires challenging hypotheses that are commonly accepted, and our interpretation gives a new perspective on the concept of propagon. This study may be considered as a good example of how mathematical modeling can bring valuable insight into biological concepts and observations.

## Introduction

Amyloids are self-perpetuating protein aggregates, involved in many neurodegenerative diseases in mammals such as Alzheimer's Disease, Parkinson's Disease, Creutzfeldt-Jakob Disease, Huntington's Disease [1]. Prions are a particular case of amyloids, where a conformational change is infectious and transmissible through self-catalyzed aggregation [2]. Understanding the molecular mechanisms associated with amyloid growth and prion replication is crucial in the effort against neurodegenerative diseases. The yeast *Saccharomyces cerevisiae* has emerged as a powerful model system in the study of these processes [3]. Native yeast prions share many specificities with mammalian prions [4], and the protein quality control machinery involved in their propagation is highly conserved between mammals and yeast. Furthermore, yeast models are used to screen for anti-amyloid drugs through the artificial expression of mammalian proteins [5]. Many kinetic models of aggregate growth and nucleation have been proposed [6], but their validation using data from yeast colonies is challenging. Indeed, the propagation of yeast prions is a multi-scale system where molecular processes are coupled to cellular processes such as cell growth and division to produce effects that are observable at the phenotypical scale, i.e. at the scale of colonies. An illustration is proposed in Fig 1. Taking into account and relating those scales calls for specific modeling approaches, but to our knowledge no mathematical modeling work describes the whole system rigorously and in a controlled way. We propose a novel approach based on the use of impulsive differential equations. It is important to emphasize that this framework is versatile and can be used to study the evolution of any chemically active intracellular components through mass-action kinetics inside growing yeast colonies.

As a first application of this approach, we focus here on the yeast prion [$PSI^+$. It has been studied in detail [7], and experimental results on [$PSI^+$] have become essential groundwork for the development of prion biology [8, 9]. The [$PSI^+$] phenotype is associated with the oligomerization and loss-of-function of the Sup35 protein, an enzyme release factor that promotes stop-codon recognition. When aggregates are formed, a change of color of the yeast colonies is observed (from red to white) if they encode a stop-codon in either the ADE1 or ADE2 genes [10]. See Fig 1C for an illustration of the change in phenotype. As is the case with most native yeast prions, the [$PSI^+$] phenotype is reversible. It can be cured using different physico-chemical treatments [11, 12] but also by introducing point mutations in Sup35 or its co-chaperones [13, 14]. Curing experiments are widely used to derive information on the molecular properties of Sup35 aggregates, depending on the curing treatment applied or the strain studied. The most exploited one is Guanidine Hydrochloride (GdnHCl) treatment, because it allows biologists to infer prion seed numbers in [$PSI^+$] cells [15, 16]. We reproduce some of those curing
<figure><img src="image_4.png" /><figcaption>**Fig 1. Different scales of yeast prion propagation.** (A) Molecular processes of protein conversion and aggregation. (B) Growth and budding, along with aggregate transmission from mother to daughter. (C) Typical experimental observation of yeast colonies (circles) in a culture plate (courtesy of T. Serio). Each colony was founded by an isolated single cell.</figcaption></figure>

https://doi.org/10.1371/journal.pcbi.1007647.g001

experiments using a model for aggregate replication at the molecular scale combined with a model for cell growth and division, as described in Fig 2. By coupling those aspects, we build a full model of prion propagation inside yeast colonies, which is detailed in Methods. Starting as simple as possible, our model sheds light onto two specific characteristics of the [PSI⁺] prion that were not uncovered in the past, as we detail in Results. Firstly, the replication of aggregates at low densities follows non-linear effects, with a kinetic barrier preventing the expansion of the population of aggregates below a certain threshold. A second implication of our modeling work is that the effect of GdnHCl is not a complete block of aggregate fragmentation, as was suggested before [11], rather a significant reduction of the replication efficiency.

## Results

We develop a multi-scale model of prion replication and propagation in yeast to investigate the relationship between molecular aggregation mechanisms and colony level phenotypical properties. The model is presented in Fig 2 and Table 1. A detailed mathematical introduction is led in Methods and an analytical study in S1 Text. In order to understand the results one should keep in mind that, in our model, newly born (daughter) yeast cells have a longer cell-cycle than their mothers. Furthermore, aggregate transmission at the moment of cell division is biased towards retention in the mother. Our goal was to develop the minimal molecular model required to reproduce experimental observations. In this work, we compare to observations for the [PSI⁺] prion, but it must be stressed that it could be adapted to any amyloid or prion propagating in growing yeast colonies. Using experimental results from the literature, we build a simple yet satisfying model and derive conclusions regarding the molecular mechanism of aggregate replication.

## Aggregate replication is limited by a kinetic barrier

**The [PSI⁺] phenotype is reversible.** We investigate the general concept of [PSI⁺] curing. A cell is considered free of [PSI⁺] if, when plated onto regular medium (free of any treatment) it grows into a [psi⁻] colony, as determined by a color assay (i.e. it grows into a dark red colony). In the same exact environmental and genetic conditions cells can grow either into [PSI⁺]
A

<figure><img src="image_5.png" /></figure>

B

<figure><img src="image_6.png" /></figure>

C

$$ 
\begin{aligned}
\frac{dV}{dt} &= \lambda - \gamma_k V - \rho V f(S), \\
\frac{dS}{dt} &= \rho V f(S) - \gamma_k S, \\
V(t_k^+) &= V(t_k^-), \\
S(t_k^+) &= (1 + \alpha_k) S(t_k^-).
\end{aligned}
 $$

<caption>**Fig 2.** Multi-scale model of yeast prion propagation. (A) Cell growth and division model used for yeast budding. Cells are growing exponentially, and divide when they reach the volume $V_0$. Division is asymmetrical and daughter cells have a smaller volume than mother cells. (B) Chemical model of aggregate replication, where V is the concentration of normal soluble Sup35 and S is the concentration of Sup35 in the prion conformation. (C) Equation system used in the model. The subscript k denotes the impulsion number. The function f is a sigmoid that accounts for non-linear cooperative effects of aggregate replication. See Methods for a description of each term.</caption>

https://doi.org/10.1371/journal.pcbi.1007647.g002

or [psi⁻] colonies, depending on the history of treatments they have sustained, as is illustrated in Fig 1C. This property is exploited to track the dynamics of curing during continuous treatment with a de-stabilizing agent [11, 12, 15]. In order to reproduce experimental results involving [PSI⁺] curing, this is an essential aspect that must be captured by the model.

**A multi-stable molecular model.** In terms of mathematical modeling, the fact that different outcomes are possible in the same conditions is interpreted as multi-stability. Different stable solutions are able to attract the trajectories of the system, and the outcome depends on the initial condition. Our model introduces this property by including a non-linearity in the replication rate of aggregates, with a sigmoidal dependency on the concentration of aggregates portrayed by the function f in the model equations (see Fig 2). This shape of function is frequently used to model cooperative reactions [20], in particular it appears naturally from ligand-receptor interaction equations. Without knowing the specific kinetic scheme, we empirically choose a Hill function. With this non-linearity, the colony grown from a cell will exhibit one of three possible behaviors depending only on its initial condition.
<table><caption>Table 1. Parameter values used by default in the simulations.</caption><thead><tr><th>Parameter</th><th>Value</th><th>Unit</th><th>Description</th></tr></thead><tbody><tr><td>T<sub>M</sub></td><td>2*</td><td>hr</td><td>Mother doubling time</td></tr><tr><td>T<sub>D</sub></td><td>3*</td><td>hr</td><td>Daughter doubling time</td></tr><tr><td>π</td><td>0.6<sup>†</sup></td><td>-</td><td>Mother/daughter volume ratio</td></tr><tr><td>ε</td><td>0.1</td><td>-</td><td>Mass transmission bias</td></tr><tr><td>α<sub>M</sub></td><td>ε/π</td><td>-</td><td>Mother concentration bias</td></tr><tr><td>α<sub>D</sub></td><td>-ε/(1-π)</td><td>-</td><td>Daughter concentration bias</td></tr><tr><td>γ<sub>M</sub></td><td>-1/T<sub>M</sub> ln(π) = 0.26</td><td>hr<sup>-1</sup></td><td>Mother growth rate</td></tr><tr><td>γ<sub>D</sub></td><td>-1/T<sub>D</sub> ln(1-π) = 0.31</td><td>hr<sup>-1</sup></td><td>Daughter growth rate</td></tr><tr><td>λ</td><td>0.7<sup>‡</sup></td><td>μM.hr<sup>-1</sup></td><td>Sup35 monomer basal production rate</td></tr><tr><td>ρ</td><td>10</td><td>hr<sup>-1</sup></td><td>Maximal aggregate replication rate</td></tr><tr><td>K</td><td>1</td><td>μM</td><td>Replication efficiency threshold</td></tr><tr><td>n</td><td>5</td><td>-</td><td>Replication efficiency order</td></tr></tbody></table>

Simulations presented in the text use those values unless specified otherwise. See Fig 2 and Methods for a description of the equations.

https://doi.org/10.1371/journal.pcbi.1007647.t001

* A full [PSI<sup>+</sup>] (totally white) colony, when all cells have a positive concentration of aggregates.

* A full [psi<sup>-</sup>] (dark red) colony, when all cells have a zero concentration of aggregates.

* A sectored colony (part white, part red), when parts of the colony are [PSI<sup>+</sup>] and others are [psi<sup>-</sup>].

Numerically, the aggregate concentration never reaches exactly zero even when the solution is converging towards the prion-free solution. In practice, when needed, we differentiate between a [PSI<sup>+</sup>] cell and a [psi<sup>-</sup>] cell by comparing the aggregate concentration to a critical value, that was set to 0.5μM. Fig 3 shows an example of simulation output for each of those cases, and Fig 1C shows a typical experiment exhibiting all three possible cases of phenotypes. It must be emphasized that without the bi-stability introduced in the molecular kinetic scheme, those three outcomes would not be possible simultaneously (i.e. without changing the model parameters). This result is justified mathematically in Methods. These only represent the fate of three different cells growing in the same exact conditions, but starting with a different aggregated state of Sup35. Those three types of behavior are observed in every curing experiment monitored by color assays [7]. By numerically investigating the model, we derive a map linking the aggregated state of Sup35 in the founder cell to the phenotype of the colony grown from that cell, as shown in Fig 4. This map is numerically derived from the study of the two extreme lineages, mother-only and daughter-only. Starting from an initial condition in the (V, S) plane, if the mother system converges to its prion-free solution then this initial condition is scored as [psi<sup>-</sup>]. If the daughter system converges to its prion solution, this initial condition is scored as [PSI<sup>+</sup>]. For cases in-between when the mother system converges to the prion solution and the daughter system converges to the prion-free solution, the initial condition is scored as sectored. Those two lineages correspond to periodic impulsive systems, and their trajectories are attracted to periodic solutions. In each case (mother system or daughter system), two different periodic solutions are locally stable, a prion-free periodic solution and a

<footnote>*:[17].</footnote>

<footnote>†:[18].</footnote>

<footnote>‡:[19]</footnote>
<figure><img src="image_7.png" /><figcaption>**Fig 3.** Representative simulation outputs of the bi-stable impulsive model. Evolution of the aggregate S and monomer V concentration for three different initial conditions. Each trajectory shown corresponds to a different lineage. Because cells are dividing the number of trajectories plotted increases in time, as does the number of cells. The trajectories for S and V are shown in the same color for the same cell. The parameters used are described in Table 1. The initial conditions used are to be related with the attraction basin showed in Fig 4. (A) Completely cured ($[\text{psi}^-]$) colony ($V(0) = 0.07\mu M, S(0) = 0.3\mu M$) (B) Full $[\text{PSI}^+]$ colony ($V(0) = 0.07\mu M, S(0) = 0.5\mu M$) (C) Sectored colony ($V(0) = 0.07\mu M, S(0) = 0.4\mu M$). The simulation codes used to create those figures are provided in S1 Codes.</figcaption></figure>

https://doi.org/10.1371/journal.pcbi.1007647.g003

$[\text{PSI}^+]$ periodic solution. A similar behavior was proven with a one-dimensional bi-stable equation in [21], but our results are only numerical so far. Those periodic solutions are denoted on the phenotype map in Fig 4. Each panel contains a total of four periodic solutions, as there are two for each system (mother-only or daughter-only). The prion-free periodic solution appears on the horizontal axis (since $S(t) = 0$), and the prion periodic solution is well above the axis. They provide information regarding the whole colony, because they represent the most extreme behavior imposed by the transmission bias. After a sufficient amount of time, any lineage in the colony is predicted to be in-between those two lineages.

**The particular case of colony sectoring.** The phenomenon of colony sectoring is the consequence of both the kinetic barrier and asymmetrical division. Indeed, our model includes a transmission bias that favors retention of aggregates by mother cells. (This bias is supported by numerous studies [7, 9, 22]). This means that under conditions of $[\text{PSI}^+]$ de-stabilization, daughter cells are more likely to fall under the threshold of efficient aggregate sustaining, and become $[\text{psi}^-]$. This provides an explanation for sectoring that is relevant with many aspects of the biology. First, sectoring occurs more often when cells are treated with de-stabilizing agents. In Fig 4, we see how reducing the proportion of aggregated Sup35 will bring cells into the region corresponding to a sectored phenotype. Furthermore, we also know that “weak” strains are more likely to exhibit sectors [9]. Fig 4(B) shows the phenotype map for a strain that has a slightly reduced aggregate replication rate. We see immediately that the region leading to sectored phenotypes is enlarged, and also closer to the periodic solutions. This means that such a
<figure><img src="image_8.png" /><figcaption>**Fig 4. Predicting colony color phenotype.** Numerical prediction of colony phenotype based on the state of the founder cell, for two different strains. A cell is predicted to be [*psi*⁻] if the mother-only lineage converges to the prion-free solution. If the daughter-only lineage converges to the prion solution, the cell is scored as [*PSI*⁺]. If the daughter-only lineage is cured of prion and the mother-only lineage is not, the cell is scored as sectored. The stable periodic solutions for the mother-only and daughter-only solutions are also shown. Each periodic solution is shown as two dots linked together by a trajectory of the differential equation system as well as a vertical line representing the effect of cell division. Those elements may overlap due to the scaling. For each system (mother or daughter) a prion-free solution is visible on the horizontal axis, as well as a prion solution in the white region of the diagram. The parameters used for this diagram are described in Table 1. Panel (A) shows a strong strain with a maximal replication rate ρ = 10hr⁻¹ and panel (B) is a weak strain with ρ = 0.87hr⁻¹. The simulation codes used to create those figures are provided in S1 Codes.</figcaption></figure>

https://doi.org/10.1371/journal.pcbi.1007647.g004

strain is more likely to be de-stabilized by random perturbations or by chemical treatment. In the literature, sectoring is often dismissed by counting mosaic colonies as either full [*PSI*⁺] or half and half [3, 11]. Perhaps some valuable information could be obtained by studying the phenomenon in more detail.

**GdnHCl does not completely block aggregate replication**

**The propagon count experiment.** Among the various agents that de-stabilize [*PSI*⁺], GdnHCl is one of the most studied. It acts by impairing the chaperone Hsp104, which fragments Sup35 aggregates, thus preventing the creation of new aggregates [11]. [*PSI*⁺] colonies growing under GdnHCl treatment progressively lose the prion phenotype, by dilution of the aggregates between dividing cells. This property was used to infer the number of prion “seeds” in individual cells [15] by fitting an exponential model to the curing curve. This approach was supported by [3], with the discovery that during GdnHCl treatment, the number of [*PSI*⁺] revertants (cells that will grow into a [*PSI*⁺] colony when transferred back onto regular medium) reaches a plateau. This number was interpreted as the number of propagons in the initial cell, where propagon refers to the hypothetical prion elementary seed. The usual interpretation of the GdnHCl curing experiment is illustrated by Fig 5, along with the typical curing curves obtained from this experiment. Three phases are identified. The fraction of [*PSI*⁺] cells in the colony starts at 1. After a lag time, [*psi*⁻] cells are created so the number of [*PSI*⁺] cells increases more slowly than the exponential growth of the whole colony. Finally the number of [*PSI*⁺] cells reaches a plateau which means the fraction of [*PSI*⁺] cells decreases exponentially as the colony keeps growing. Propagon counts and statistics on those counts are often used for instance to characterize strains [23], to quantify the bias of transmission between mother and
<figure><img src="image_9.png" /><figcaption>**Fig 5. The GdnHCl curing experiment and propagon counting.** (A) Two phase propagon counting experiment. Propagons are distributed in dividing cells under GdnHCl treatment. Cells are then plated onto regular medium and the color phenotype ([$PSI^+$] or [$psi^−$]) of each subsequent colony is determined by color assay. The number of [$PSI^+$] colonies is counted and gives the propagon count. (B) Qualitative representation of the curves obtained during the GdnHCl curing experiment. In black is the percentage of [$PSI^+$] cells in the colony (logarithmic scale), in red is the number of [$PSI^+$] cells in the colony (logarithmic scale). The three phases are separated by solid vertical lines. First the number of [$PSI^+$] cells grows as fast as the whole colony, which is growing exponentially all through the experiment. Then, [$psi^−$] cells appear and the fraction of [$PSI^+$] cells starts decreasing. Finally, the number of [$PSI^+$] cells reaches a plateau but the whole colony keeps growing exponentially, hence the exponential decrease of the fraction of [$PSI^+$] cells. See [9] for an experimental example of such curves.</figcaption></figure>

https://doi.org/10.1371/journal.pcbi.1007647.g005

daughter cells [17] and to study the effect of different point mutations [13]. It needs to be emphasized that this interpretation is based on two hypotheses. The first hypothesis is that the threshold between a [$PSI^+$] cell and a [$psi^−$] cell is the presence of aggregates or not. In our modeling framework, this hypothesis translates to the presence of a concentration threshold under which aggregate replication is inefficient. The second hypothesis is that GdnHCl instantaneously and totally blocks the creation of new aggregates. Our results contradict this second hypothesis as explained below.

**The propagon plateau.** During GdnHCl treatment, the number of [$PSI^+$] cells in a colony reaches a plateau [9]. Our model reproduces this property under specific parameter choices, as shown in Fig 6. In those simulations, cells are scored as [$PSI^+$] as soon as the concentration of aggregated Sup35 is above a given threshold (we used 0.5µM). However the qualitative result

<figure><img src="image_10.png" /><figcaption>**Fig 6. Reproducing a propagon count experiment.** Evolution of the [$PSI^+$] cells count and [$PSI^+$] proportion simulated by our model, in the case of $\rho = 0.21\text{hr}^{-1}$ (all other parameters as in Table 1). Cells are scored as [$PSI^+$] as soon as they contain a concentration of aggregated Sup35 higher than 0.5µM. The simulation codes used to create those figures are provided in S1 Codes.</figcaption></figure>

https://doi.org/10.1371/journal.pcbi.1007647.g006
does not depend on the scoring method, because we know that a finite number of lineages in the colony are attracted by a solution with a positive concentration of aggregates. All the other lineages are attracted by the prion-free solution. This behavior is possible with our model but only if aggregates continue to replicate in GdnHCl conditions, which strongly contradicts former experimental studies [11, 15]. The hypothesis that GdnHCl interrupts all chemical activity of aggregates was nuanced in previous experimental work. Indeed it was shown GdnHCl does not stop aggregate growth [11, 17, 24]. Furthermore, [25] used fluorescent tagging to track aggregates during GdnHCl treatment, and report a decrease in cells with foci slower than the halving predicted by the segregation hypothesis.

Reproducing the propagon count experiment with our model requires very specific kinetic parameters for the replication reaction. Those conditions are found numerically by studying the mother-only lineage in the model. Indeed when this lineage is attracted to a [PSI⁺] solution but each of its daughters becomes [psi⁻], we are assured that the total number of [PSI⁺] cells in the colony reaches a plateau. Once again this is the consequence of bi-stability in the molecular model as well as the asymmetric division favoring aggregate retention by mother cells. The fact that those conditions correspond to a very narrow parameter range is worth emphasizing and discussing further.

# Discussion

## How to explain the concentration threshold for aggregate replication

In previous modeling work, one approach was to study aggregates as discrete entities [15, 26, 27]. When doing so, the kinetic barrier hypothesis is implicit and the threshold lies at the difference between 0 and 1 aggregate. However, it is not a proper modeling approach when it comes to studying chemical reactions, because mass-action kinetics cannot be applied to numbers in isolation but rather to densities. Our model implements chemical modeling inside a multi-scale phenomenon in a rigorous way and provides a powerful tool for generating insights into existing experimental outcomes: colony sectoring and propagon counting assays.

Our work suggests that the multi-stability in the intracellular scale is essential to reproduce the experimental results for yeast prions. That is, the prion-free state and the prion state both need to be locally stable. But what is the mechanistic origin of this multi-stability? The non-linearity we introduced could be interpreted by a cooperativity phenomenon, where the replication of aggregates (for instance fragmentation) is chemical reaction that is catalyzed by the presence of other aggregates. Due to the lack of clear biological evidence, there is so far no real hypothesis to explain such a mechanism. One idea is to take into account the action of chaperones, in particular Hsp104. However, even if fragmentation is limited by the presence of Hsp104 the model does not exhibit bi-stability [28]. We would like to emphasize that this type of modeling (global models) is universally used in the prion biology, in particular with the widespread use of the Nucleated Polymerization Model (NPM) [29, 30]. The introduction of bi-stable systems was suggested in earlier work [31, 32], in the case of mammals and more recently to study interaction between prion strains [33]. The benefit of the framework we suggest here is that it could allow to use data to validate a precise kinetic scheme of aggregate replication. In particular, the sectoring phenomenon should be studied with a new perspective, especially during curing experiments.

## What is the effect of GdnHCl?

Our results question all previous assumptions made about the effect of GdnHCl. With our model, the propagon experiment is reproduced only if the chemical replication rates are chosen very carefully. This is potentially explained by two reasons. The first is that our model
might be too simple to capture the possibilities of the full biological system. This would also
explain why sectoring is reduced to such a narrow region in the phenotype map Fig 4.

A second explanation is that our model may be missing an essential chemical process,
which remains unaffected when GdnHCl is present. There is precedent for this assumption,
because there is evidence that aggregates are still chemically active under GdnHCl treatment.
First GdnHCl does not prevent aggregates from growing by polymerizing newly synthesized
Sup35 [17, 34]. Furthermore, there is evidence for an action of Hsp104 that is not affected by
GdnHCl [35, 36]. This brings up the controversial question of the roles of Hsp104 in the prop-
agation of [PSI<sup>+</sup>] and other yeast prions. The only conclusion our results bring is that GdnHCl
curing is not explained by an exponential dilution model as suggested previously [15, 18],
because of the plateau of [PSI<sup>+</sup>] cells. Our study is a first step in the design of a more elaborate
kinetic model that includes the effect of Hsp104.

Conclusion and perspectives

We introduced a novel modeling tool to the field of yeast prions, with the major benefit of relating different scales in a controlled and rigorous way. From the molecular mechanisms, with a kinetic scheme built using mass-action kinetics, to the phenotypical traits at the col-ony level, this framework has the potential of taking into account every aspect of the system. In the case of the [PSI<sup>+</sup>] prion, we build a simple model with the primary goal to qualitatively reproduce experimental observations from the literature. We focus on curing experiments, and reproducing any curing experiment with our framework requires introducing a very particular characteristic into the molecular model. The kinetic scheme needs to be bi-stable, where the prion-free state and the prion state are both simultaneously stable, and the transi-tion between the two of them is a bi-stable switch. This allows the possibility of curing, and it concomitantly explains the phenomenon of sectoring in a deterministic way. This phenome-non is often dismissed but is in fact instructive with regards to the molecular processes. By investigating in more detail the case of GdnHCl curing, we have reason to question the suggested effect of this agent. Indeed this experiment is reproduced by our model, but the fragmentation of aggregates must not be completely inhibited contrary to the commonly accepted effect of GdnHCl.

Overall, the framework of impulsive differential equations is versatile and could be adapted
to many different cases. Studying the [PSI<sup>+</sup>] prion already revealed instructive, even though
work is still in progress. In the future, we aim to use this framework to build and validate a
complete model of aggregate replication including the role of Hsp104 and possibly its co-chap-
erones, a size-distribution of aggregates, stochasticity in the cell division events. This would be
done through hypothesis testing and close collaboration with biologists. Inferring parameters
is a long-term goal, that first requires understanding the very structure of the molecular pro-
cesses. In particular, it needs to be clear what is the mechanistic origin of the cooperativity
in the replication of aggregates. Another use of the model is to extend it to other yeast prions
and amyloid models, as yeast models are used to screen for anti-amyloid drugs, and a specific
modeling framework would help interpreting those experiments.

Methods

Using impulsive differential equations to model cell division

We first introduce impulsive differential equations (IDEs) and use them to model a population
of dividing yeast cells. Given a system of ordinary differential equations dZ/dt = f(Z, t), where
Z ∈ ℝⁿ represents the concentrations of the biochemical species we are tracking, and given
$Z_0 \in \mathbb{R}^n, Z_0 \ge 0$ an initial condition, we define the impulsive system

$$ 
\left\{
\begin{array}{l}
\displaystyle \frac{dZ}{dt} = f(Z, t), \quad t \in \mathbb{R}_+ \setminus \{t_k, k \in \mathbb{N}\}, \\
Z(t_k^+) - Z(t_k^-) = \alpha_k \cdot Z(t_k^-), \quad k \in \mathbb{N}, \\
Z(0) = Z_0.
\end{array}
\right.
$$

The sequence $(t_k)_{k \in \mathbb{N}^*}$ represents the impulsion times, which we choose to be fixed. We require $0 \le t_0 < \dots < t_k < \dots$ and $\lim_{k \to \infty} t_k = \infty$. Each impulsion corresponds to a cell division, and the vector $\alpha_k \in \mathbb{R}^n$ represents the partition rate at the moment of division $k$. A solution to this system is a piece-wise continuous function $Z$ on $(t_k, t_{k+1})$, for $k \in \mathbb{N}^*$. The existence and the properties of such solutions are studied theoretically [37, 38]. IDEs are used in various modeling fields, including epidemiology [39], population dynamics [40] and more recently T-cell differentiation [21].

In order to model cell division with IDEs, we need to clarify how species are distributed between mother and daughter cells. Consider a yeast cell about to bud and produce a daughter. The volume of the mother cell is $V_0$, which will split into two parts right after division. The mother's volume becomes $V_M = \pi V_0$ and the daughter gets $V_D = (1 - \pi)V_0$, where $\pi$ is the volume asymmetry. For yeast cells it has been observed that $\pi = 0.6$ [18]. Next consider a chemical component in the mother of mass $M_0$ and concentration $C_0 = \frac{M_0}{V_0}$. During division this mass is distributed between the mother and the daughter cells. With no bias (fast diffusion) it should distribute with the same ratio as the volume. If we introduce a bias, the mother retains the mass $M_M = (\pi + \varepsilon)M_0$ and the daughter receives $M_D = (1 - \pi - \varepsilon)M_0$. The parameter $\varepsilon$ represents the bias, with a positive value corresponding to a retention of the mass by the mother. In order to ensure positive mass we require

$$ -\pi < \varepsilon < 1 - \pi. $$

After division, the mother concentration ($C_M$) and the daughter concentration ($C_D$) are given by

$$ C_M = (1 + \alpha_M)C_0, $$

$$ C_D = (1 + \alpha_D)C_0, $$

$$ \text{where } \alpha_M = \frac{\varepsilon}{\pi}, \alpha_D = -\frac{\varepsilon}{1 - \pi}. $$

We do not consider degenerate cases where $\varepsilon = -\pi$ or $\varepsilon = 1 - \pi$, corresponding to the case where the transmission from mother to daughter is (respectively) full or null. This hypothesis is biologically relevant in the sense that it seems highly unlikely that exactly no misconformed protein be transmitted at the moment of division.

We assume that cells are exponentially growing, and that they divide when they reach the volume $V_0$. In other words, the intracellular dynamics of prion aggregation have no impact on the division cycle. This is consistent with [17]. This assumption imposes a relation on the growth rate $\gamma$ and the division time $T$. As before, if we denote the mother parameters with a subscript $M$ and the daughter parameters with a subscript $D$, we have the following two relations

$$ e^{-\gamma_M T_M} = \pi, $$

$$ e^{-\gamma_D T_D} = 1 - \pi. $$

As a consequence when $\pi$, $T_M$ and $T_D$ are fixed (as measured experimentally [9, 17]), $\gamma_M$ and
$\gamma_D$ are determined by these relations. Fig 2A shows a visual representation of the cell division model we are using. (The parameters we used in results are given in Table 1).

An impulsive differential equation model allows one to follow a cell lineage and its intracellular components in time. The cell history is described by the sequence of impulsions it undergoes, either mother or daughter impulsions. Analytically, we study in detail the two extreme lineages, mother-only and daughter-only. Indeed those lineages are subject to periodic and identical impulsions, which facilitates the analysis [37]. If we use the IDEs to track all lineages starting from a single cell, we can study the complete yeast colony. Having formalized the cell division dynamics model, we next introduce the model for intracellular dynamics.

## A bi-stable prion replication process

In this section, we introduce a two-dimensional bi-stable model of prion dynamics, at first without considering the effect of impulsions. In each yeast cell, we track the concentration of soluble Sup35 (V), and of Sup35 in the prion conformation (S). The model is illustrated in Fig 2B, and described by the following system of ordinary differential equations

$$ \frac{dV}{dt} = \lambda - \gamma V - \rho Vf(S), $$

$$ \frac{dS}{dt} = \rho Vf(S) - \gamma S. $$

The parameters $\lambda$ and $\gamma$ are respectively the Sup35 monomer production rate and the cell growth rate. For simplicity we do not include the aggregate size dynamics but model the recruitment of monomers into prion aggregates. It occurs with a cooperative reaction, at maximal rate $\rho$ and non-linear efficiency $f(S)$. The function $f$ is chosen to be a Hill function of order $n(>1)$ and threshold $K$

$$ f(S) = \frac{S^n}{K^n + S^n}. $$

This function is chosen for one distinct feature: it makes the system multi-stable as soon as $n > 1$. The prion-free equilibrium (with $S = 0$) exists under any choice of parameters and remains locally stable. Two other equilibria appear from a saddle-node bifurcation, with one of them unstable and the other locally stable. In conditions when the three equilibria exist (one prion-free and two prion equilibria), the asymptotic outcome of the model depends on the initial condition. Starting with a concentration of aggregates too low will cause the solution to converge to the prion-free equilibrium, whereas if the initial concentration of aggregates is sufficient their population will be stably maintained. See S1 Text for the mathematical description and proof of this property, as well as some numerical illustration of the bi-stability. As detailed in Discussion, it is essential that we use a model with multi-stability because we are investigating a curing experiment, where cells can lose the prion phenotype.

## Complete bi-stable model with impulsions

When combining the bi-stable model and the impulsive differential equation framework, we obtain a complete model of aggregate replication and transmission in growing yeast colonies. The choice of parameters used by default is detailed in Table 1. The cell division parameters (doubling time and mother-daughter volume ratio) are measured experimentally and we choose values within the typical observed range [16, 17]. The mass transmission bias is a parameter that will require a thorough investigation in further work but we choose a reasonable value of $\varepsilon = 0.1$ for the time being. The Sup35 monomer production rate is such that the
prion-free steady-state concentration of Sup35 in cells is $\frac{\lambda}{\gamma} \approx 2.5\mu M$ [19]. The prion replication parameters $\rho$, $K$ and $n$ are used as adjustment parameters to investigate the behavior of the model.

**Supporting information**

**S1 Text. Analytical study of the bi-stable impulsive system for yeast prion propagation.**
(PDF)

**S1 Codes. Matlab simulation codes used to produce Figs 3, 4 and 6.**
(ZIP)

**Acknowledgments**

The authors would like to thank Dr. Tricia R. Serio for contributing a figure of yeast colonies and offering valuable comments.

**Author Contributions**

**Conceptualization:** Paul Lemarre.

**Formal analysis:** Paul Lemarre.

**Funding acquisition:** Laurent Pujo-Menjouet, Suzanne S. Sindi.

**Investigation:** Paul Lemarre, Laurent Pujo-Menjouet, Suzanne S. Sindi.

**Software:** Paul Lemarre.

**Supervision:** Laurent Pujo-Menjouet, Suzanne S. Sindi.

**Validation:** Laurent Pujo-Menjouet, Suzanne S. Sindi.

**Writing – original draft:** Paul Lemarre.

**Writing – review & editing:** Paul Lemarre, Laurent Pujo-Menjouet, Suzanne S. Sindi.

**References**

1. Soto C. Unfolding the role of protein misfolding in neurodegenerative diseases. Nature Reviews Neuroscience. 2003; 4(1):49. https://doi.org/10.1038/nrn1007 PMID: 12511861

2. Weissman JS, Chien P, DePace AH. Emerging Principles of Conformation-Based Prion Inheritance. Annual Reviews Biochem. 2004; 73:617–656. https://doi.org/10.1146/annurev.biochem.72.121801.161837

3. Tuite MF, Cox BS. Propagation of yeast prions. Nature Reviews Molecular Cell Biology. 2003; 4:878–889. https://doi.org/10.1038/nrm1247 PMID: 14625537

4. Tuite MF, Serio TR. The prion hypothesis: from biological anomaly to regulatory mechanism. Nature Reviews Molecular Cell Biology. 2010; 11(12):823–833. https://doi.org/10.1038/nrm3007 PMID: 21081963

5. Chernova TA, Chernoff YO, Wilkinson KD. Yeast Models for Amyloids and Prions: Environmental Modulation and Drug Discovery. Molecules. 2019; 24(18):3388. https://doi.org/10.3390/molecules24183388

6. Sindi SS. Mathematical Modeling of Prion Disease. In: Prion-An Overview. InTech; 2017.

7. Cox B, Tuite M. Life of [PSI]. Curr Genet. 2018; 64:1–8. https://doi.org/10.1007/s00294-017-0714-7 PMID: 28653109

8. Sindi SS, Serio TR. Prion Dynamics and the Quest for the Genetic Determinant in Protein-Only Inheritance. Current Opinion in Microbiology. 2009; 12(6):623–630. https://doi.org/10.1016/j.mib.2009.09.003 PMID: 19864176
9. Cox BS, Byrne L, Tuite MF. Prion stability. Prion. 2007; 1(3):170-178. https://doi.org/10.4161/pri.1.3.4839 PMID: 19164897

10. Liebman SW, Chernoff YO. Prions in yeast. Genetics. 2012; 191(4):1041-1072. https://doi.org/10.1534/genetics.111.137760 PMID: 22879407

11. Wegrzyn RD, Bapat K, Newnam GP, Zink AD, Chernoff YO. Mechanism of Prion Loss after Hsp104 Inactivation in Yeast. Molecular and Cellular Biology. 2001; 21(14):4656-4669. https://doi.org/10.1128/MCB.21.14.4656-4669.2001 PMID: 11416143

12. Zhao X, Rodriguez R, Silberman RE, Ahearn JM, Saidha S, Cummins KC, et al. Heat shock protein 104 (Hsp104)-mediated curing of [PSI'] yeast prions depends on both [PSI'] conformation and the properties of the Hsp104 homologs. J Biol Chem. 2017; 292(21):8630-8641. https://doi.org/10.1074/jbc.M116.770719 PMID: 28373280

13. Langlois CR, Pei F, Sindi SS, Serio TR. Distinct prion domain sequences ensure efficient amyloid propagation by promoting chaperone binding or processing in vivo. PLoS genetics. 2016; 12(11):e1006417. https://doi.org/10.1371/journal.pgen.1006417 PMID: 27814358

14. Pei F, DiSalvo S, Sindi SS, Serio TR. A dominant-negative mutant inhibits multiple prion variants through a common mechanism. PLoS genetics. 2017; 13(10):e1007085. https://doi.org/10.1371/journal.pgen.1007085 PMID: 29084237

15. Eaglestone SS, Ruddock LW, Cox BS, Tuite MF. Guanidine hydrochloride blocks a critical step in the propagation of the prion-like determinant [PSI'] of Saccharomyces cerevisiae. PNAS. 2000; 97(1). https://doi.org/10.1073/pnas.97.1.240 PMID: 10618402

16. Cox B, Ness F, Tuite M. Analysis of the generation and segregation of propagons: entities that propagate the [PSI+] prion in yeast. Genetics. 2003; 165(1):23-33. PMID: 14504215

17. Satpute-Krishnan P, Langseth SX, Serio TR. Hsp104-Dependent Remodeling of Prion Complexes Mediates Protein-Only Inheritance. PLoS Biology. 2007; 5(2). https://doi.org/10.1371/journal.pbio.0050024 PMID: 17253904

18. Byrne LJ, Cole DJ, Cox BS, Ridout MS, Morgan BJT, Tuite MF. The Number and Transmission of [PSI'] Prion Seeds (Propagons) in the Yeast Saccharomyces cerevisiae. PLoS ONE. 2009; 3(4).

19. Ho B, Baryshnikova A, Brown GW. Unification of protein abundance datasets yields a quantitative Saccharomyces cerevisiae proteome. Cell systems. 2018; 6(2):192-205. https://doi.org/10.1016/j.cels.2017.12.004 PMID: 29361465

20. Weiss JN. The Hill equation revisited: uses and misuses. FASEB J. 1997; 11. https://doi.org/10.1096/fasebj.11.11.9285481 PMID: 9285481

21. Girel S, Crauste F. Existence and stability of periodic solutions of an impulsive differential equation and application to CD8 T-cell differentiation. Journal of Mathematical Biology. 2018; https://doi.org/10.1007/s00285-018-1220-3 PMID: 29500513

22. Derdowski A, Sindi SS, Klaips CL, DiSalvo S, Serio TR. A size threshold limits prion transmission and establishes phenotypic diversity. Science. 2010; 330(6004):680-683. https://doi.org/10.1126/science.1197785 PMID: 21030659

23. Bradley ME, Liebman SW. Destabilizing interactions among [PSI+] and [PIN+] yeast prion variants. Genetics. 2003; 165(4):1675-1685. PMID: 14704158

24. Ness F, Ferreira P, Cox BS, Tuite MF. Guanidine Hydrochloride Inhibits the Generation of Prion "Seeds" but Not Prion Protein Aggregation in Yeast. Molecular and Cellular Biology. 2002; 22 (15):5593-5605. https://doi.org/10.1128/MCB.22.15.5593-5605.2002 PMID: 12101251

25. Park YN, Morales D, Rubinson EH, Maison D, Eisenberg E, Greene LE. Differences in the Curing of [PSI'] Prion by Various Methods of Hsp104 Inactivation. PLoS ONE. 2012; 7(6). https://doi.org/10.1371/journal.pone.0037692

26. Sindi SS, Olofsson P. A discrete-time branching process model of yeast prion curing curves. Mathematical Population Studies. 2013; 20(1). https://doi.org/10.1080/08898480.2013.748566 PMID: 30287982

27. Banks HT, Flores KB, Langlois CR, Serio TR, Sindi SS. Estimating the rate of prion aggregate amplification in yeast with a generation and structured population model. Inverse Problems in Science and Engineering. 2017;.

28. Davis JK, Sindi SS. A mathematical model of the dynamics of prion aggregates with chaperone-mediated fragmentation. Journal of Mathematical Biology. 2016; 72(6):1555-1578. https://doi.org/10.1007/s00285-015-0921-0 PMID: 26297259

29. Greer ML, Pujo-Menjouet L, Webb GF. A mathematical analysis of the dynamics of prion proliferation. Journal of Theoretical Biology. 2006; 242(3):598-606. http://dx.doi.org/10.1016/j.jtbi.2006.04.010 PMID: 16753184

30. Tanaka M, Collins SR, Toyama BH, Weissman JS. The physical basis of how prions conformations determine strain phenotypes. Nature. 2006; 442:585-589. https://doi.org/10.1038/nature04922 PMID: 16810177
31. Laurent M. Bistability and the species barrier in prion diseases: stepping across the threshold or not. Biophysical chemistry. 1998; 72(1-2):211-222. https://doi.org/10.1016/S0301-4622(98)00135-5 PMID: 17029708

32. Kellershohn N, Laurent M. Prion diseases: dynamics of the infection and properties of the bistable transition. Biophysical journal. 2001; 81(5):2517-2529. https://doi.org/10.1016/S0006-3495(01)75897-3 PMID: 11606267

33. Lemanre P, Pujo-Menjouet L, Sindi SS. Generalizing a mathematical model of prion aggregation allows strain coexistence and co-stability by including a novel misfolded species. Journal of mathematical biology. 2018; p. 1-31.

34. Kryndushkin DS, Alexandrov IM, Ter-Avanesyan MD, Kushnirov VV. Yeast [PSI*] Prion Aggregates Are Formed by Small Sup35 Polymers Fragmented by Hsp104. Journal of Biological Chemistry. 2003; 278(49). https://doi.org/10.1074/jbc.M307996200 PMID: 14507919

35. Park YN, Zhao X, Yim YI, Todor H, Ellerbrock R, Reidy M, et al. Hsp104 Overexpression Cures *Saccharomyces cerevisiae* [*PSI*] by Causing Dissolution of the Prion Seeds. Eukaryotic Cell. 2014; 13(5):635-647. https://doi.org/10.1128/EC.00300-13 PMID: 24632242

36. Ness F, Cox BS, Wongwigkarn J, Naeimi WR, Tuite MF. Over-expression of the molecular chaperone Hsp104 in Saccharomyces cerevisiae results in the malpartition of [*PSI*] propagons. Molecular Biology. 2017; 104(1):125-143.

37. Bainov D, Simeonov P. Impulsive differential equations: periodic solutions and applications. Pitmans Monographs and Surveys in Pure and Applied Mathematics. Longmans Scientific and Technical; 1993.

38. Lakshmikanthan V, Bainov DD, Simeonov PS. Theory of impulsive differential equations. vol. 6 of Series in Modern Applied Mathematics. World Scientific; 1989.

39. Gao S, Teng Z, Nieto JJ, Torres A. Analysis of an SIR Epidemic Model with Pulse Vaccination and Distributed Time Delay. Journal of Biomedicine and Biotechnology. 2007; https://doi.org/10.1155/2007/64870

40. Liu X, Takeuchi Y. Periodicity and global dynamics of an impulsive delay Lasota-Wazewska model. Journal of mathematical analysis and applications. 2007; 327(1):326-341. https://doi.org/10.1016/j.jmaa.2006.04.026