Uniqueness of Weak Solutions to a Prion Equation with Polymer Joining
Elena Leis, Christoph Walker

TL;DR
This paper proves the uniqueness of weak solutions in a mathematical model describing prion proliferation involving polymerization, splitting, and joining, extending previous work on existence of solutions.
Contribution
It establishes the first proof of uniqueness for weak solutions in a complex prion proliferation model with polymer joining.
Findings
Uniqueness of weak solutions is proven for the prion model.
The model includes polymerization, splitting, and joining processes.
Global weak solutions exist for unbounded reaction rates.
Abstract
We consider a model for prion proliferation that includes prion polymerization, polymer splitting, and polymer joining. The model consists of an ordinary differential equation for the prion monomers and a hyperbolic nonlinear differential equation with integral terms for the prion polymers and was shown to possess global weak solutions for unbounded reaction rates [11]. Here we prove the uniqueness of weak solutions.
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Taxonomy
TopicsPrion Diseases and Protein Misfolding · Mathematical and Theoretical Epidemiology and Ecology Models
Uniqueness of Weak Solutions to
a Prion Equation with Polymer Joining
Elena Leis
Leibniz Universität Hannover, Institut für Angewandte Mathematik, Welfengarten 1, D–30167 Hannover, Germany
and
Christoph Walker
Leibniz Universität Hannover, Institut für Angewandte Mathematik, Welfengarten 1, D–30167 Hannover, Germany
Abstract.
We consider a model for prion proliferation that includes prion polymerization, polymer splitting, and polymer joining. The model consists of an ordinary differential equation for the prion monomers and a hyperbolic nonlinear differential equation with integral terms for the prion polymers and was shown to possess global weak solutions for unbounded reaction rates [11]. Here we prove the uniqueness of weak solutions.
Key words and phrases:
Prions, polymer joining, integro-differential equation, weak solutions, uniqueness.
1. Introduction
In this article we consider a mathematical model for the dynamics of prions which are thought to be misfolded proteins and cause deadly neurodegenerative diseases including “mad cow disease” in mammals. The model describes the proliferation of prions and was introduced in [8] to which we refer for more information regarding the biological background. The infectious prions are treated as polymers and interact with the noninfectious monomer form. The model includes polymerization, polymer joining, and polymer splitting. These processes can mathematically be described by a coupled system consisting of an ordinary differential equation for the number of noninfectious monomers , given by
[TABLE]
and an integro-differential equation for the density distribution function of infectious polymers of size of the form
[TABLE]
for . These equations are supplemented with the boundary condition
[TABLE]
and the initial values
[TABLE]
Here, is a constant monomer background source while and are the metabolic degradation rates for monomers, respectively, -polymers. The function is the splitting rate for a polymer of size into two polymers of size and , where is the probability (density) for this event. Any daughter polymer with size less than the critical size is assumed to disintegrate instantaneously into monomers. In the polymerization process, infectious polymers of size attach noninfectious monomers at rate . If there is a saturation effect when the number of monomers forming the infectious polymers becomes large resulting in less lengthening overall. Two polymers of size and may join at rate . Note that equation (1.2) is reminiscent of the continuous coagulation-fragmentation equation known from physics (see e.g. [5] and the references therein).
For the case , that is, when the bilinear polymer joining terms are neglected, equations (1.1)-(1.4) were studied in [6, 10, 13, 14] with respect to existence and uniqueness and in [1, 2, 4, 6, 7, 13] with respect to qualitative aspects. The model with polymer joining was introduced in [8]. There it was assumed that the rates have the particular form
[TABLE]
In this case, equation (1.2) can be integrated and a closed system of ordinary differential equations for the unknowns , , and can be derived that can be globally solved. The equations (1.1)-(1.2) then decouple since is determined (see also [6, 9, 12]).
When polymer joining is taken into account, but (1.5) is not assumed, the existence of solutions to (1.1)-(1.4) was established in [11]. More precisely, it was shown that for bounded reaction rates , , , and a unique global classical solution exists. For unbounded (and thus biologically more relevant) reaction rates satisfying certain growth restrictions, the existence of global weak solutions was established. However, the uniqueness of weak solutions was left open and it is the purpose of this article to fill this gap. We shall show herein that under reasonable growth conditions on the reaction rates there is at most one weak solution. We thus extend the result of [10] to include polymer joining using the same techniques. We shall point out here that we also use ideas from [5] on the coagulation equation to handle the latter (see also [3] where similar techniques are used to investigate uniqueness for the coagulation-fragmenation equations). The uniqueness result from the present work complements the existence result of [11] to provide the well-posedness of (1.1)-(1.4) in the framework of weak solutions.
Before introducing the notation of a weak solution we remark that solutions are supposed to preserve the number of monomers. More precisely, let the splitting kernel be a measurable function defined on satisfying the symmetry condition
[TABLE]
and being normalized according to
[TABLE]
Then splitting conserves the number of monomers and (1.6), (1.7) imply
[TABLE]
Note that if is e.g. of the form
[TABLE]
with a non-negative integrable function defined on satisfying
[TABLE]
then conditions (1.6), (1.7) hold. In particular, for one obtains the rate
[TABLE]
from (1.5) as considered in [6, 8, 9]. We further assume that the polymer joining rate is non-negative and symmetric, that is,
[TABLE]
Throughout this article we assume that satisfies conditions (1.6), (1.7) while satisfies (1.11). We also assume that
[TABLE]
and that is a positive measurable function on growing at most linearly. It is then straightforward to check that (1.6), (1.7), and (1.11) imply that any solution to (1.1)-(1.4) satisfies (formally) the monomer balance law
[TABLE]
at time . That is, the overall number of monomers changes only due to natural production or metabolic degradation. To keep track of the biologically important quantities
[TABLE]
of all polymers respectively monomers forming those polymers, we shall thus consider solutions with belonging to the positive cone of , denoted by .
Definition 1.1**.**
Given and we call a pair a monomer preserving (global) weak solution to (1.1)-(1.4) provided
- (i)
* is a non-negative solution to (1.1),*
- (ii)
u\in L_{\infty,\mathrm{loc}}\big{(}\mathbb{R}^{+},L_{1}^{+}(Y,y\mathrm{d}y)\big{)}* is a weak solution to (1.2), that is, it satisfies for all *
[TABLE]
and
[TABLE]
for any test function ,
- (iii)
the balance law (1.12) holds.
Note that the weak formulation in (ii) above is obtained by testing (1.2) against the test function and using for the operator
[TABLE]
the identity
[TABLE]
which follows from the symmetry of .
As pointed out above the existence of a global weak solution in the sense of Definition 1.1 was obtained in [11] under fairly general conditions on the reaction rates. We next state conditions under which such solutions are unique.
Main Results
We shall first state simplified and hence more illustrative versions of our actual results. To this end we temporarily assume that the rates are of the particular form
[TABLE]
for with , , and . Moreover, let us also assume that
[TABLE]
We set
[TABLE]
and
[TABLE]
Then we have the following uniqueness result for weak solutions with finite higher moments.
Theorem 1.2**.**
Suppose (1.17), (1.18) with . Further suppose that there is a constant such that, if , then
[TABLE]
while if , then
[TABLE]
Then there is (large enough and depending on , , , ) such that (1.1)-(1.4) has for each initial value with and at most one monomer-preserving weak solution in the sense of Definition 1.1 with u\in L_{\infty,\mathrm{loc}}\big{(}\mathbb{R}^{+},L_{1}(Y,y^{\sigma}\mathrm{d}y)\big{)}.
Theorem 1.2 is a special case of a more general result stated in Theorem 2.1. The latter does actually not require structural assumptions on the reaction rates as in (1.17) but rather suitable growth conditions. Note, however, that it does not provide uniqueness of weak solutions in the natural phase space . For rates and with at most linear growth we can though improve Theorem 1.2 to obtain a uniqueness result in which in particular includes the rates from (1.5). The following theorem is a special case of a more general result, see Theorem 2.2.
Theorem 1.3**.**
Suppose (1.17) with and (1.18). Further let (1.19) hold with and suppose that for . Then, given any there exists at most one monomer-preserving weak solution to (1.1)-(1.4) in the sense of Definition 1.1.
As pointed out before, the above theorems extend the results from [10] for the case , i.e. without polymer joining. To include polymer joining herein we use ideas from [5] on the coagulation equation. Assumptions (1.19), (1.20) correspond to the assumptions therein.
Combining now the uniqueness statements above with the existence results of [11, Theorem 2.3, Proposition 2.4] we obtain the well-posedness of (1.1)-(1.4) within the framework of weak solutions.
Corollary 1.4**.**
(a)* Let the assumptions of Theorem 1.2 hold for . Then there is (large enough) such that (1.1)-(1.4) has for each initial value with and a unique global monomer-preserving weak solution in the sense of Definition 1.1 such that u\in L_{\infty,\mathrm{loc}}\big{(}\mathbb{R}^{+},L_{1}(Y,y^{\sigma}\mathrm{d}y)\big{)}.
(b) Let the assumptions of Theorem 1.3 hold. Then, given any , there exists a unique global monomer-preserving weak solution to (1.1)-(1.4) in the sense of Definition 1.1.*
The restriction to in Corollary 1.4 (a) ensures the existence of a solution with belonging to for . However, the existence of a solution without this additional constraint can be shown also for , see [11].
2. Sharper Statements of the Main Results
We now state more general results than in Theorem 1.2 and Theorem 1.3 that do not rely on structural conditions on the reaction rates as in (1.17)-(1.18), but rather on growth conditions.
Let us recall that we assume throughout that , that satisfies conditions (1.6), (1.7) while satisfies (1.11). Then, as in [10], we shall further assume that
[TABLE]
and
[TABLE]
Let there be a strictly positive function
[TABLE]
for some constant such that
[TABLE]
Introducing ’s primitive
[TABLE]
we shall further assume that
[TABLE]
where
[TABLE]
In addition, assume that
[TABLE]
As for the polymer joining rate we suppose that there is a constant such that
[TABLE]
along with
[TABLE]
We then assume that there is a constant such that
[TABLE]
For the first uniqueness result we also require that
[TABLE]
and
[TABLE]
for some constant . Then we have the following uniqueness result for solutions with sufficient integrability.
Theorem 2.1**.**
Let (2.1)-(2.16) be satisfied. Then, for any initial value with and u^{0}\in L_{1}^{+}(Y,y\mathrm{d}y)\cap L_{1}\big{(}Y,(\mu+\beta)(y)G(y)\mathrm{d}y\big{)} there is at most one monomer-preserving weak solution to (1.1)-(1.4) in the sense of Definition 1.1 such that
[TABLE]
Although Theorem 2.1 applies to a wide class of reaction rates it does not include the rates from (1.5) to yield uniqueness when belongs to the natural phase space L_{\infty,\mathrm{loc}}\big{(}\mathbb{R}^{+},L_{1}(Y,y\mathrm{d}y)\big{)} but rather for u\in L_{\infty,\mathrm{loc}}\big{(}\mathbb{R}^{+},L_{1}(Y,y^{2}\mathrm{d}y)\big{)}. To remedy this issue we shall consider (2.1)-(2.13) in the particular case (then (2.14) trivially holds) and further suppose that there are and such that
[TABLE]
and
[TABLE]
Moreover, we suppose that we can decompose and in the form
[TABLE]
and that there is a constant such that
[TABLE]
Then we obtain the following uniqueness result in the natural phase space .
Theorem 2.2**.**
Suppose that (2.1)-(2.13) hold with and let (2.17)-(2.20) be satisfied. Then, for any initial value with and there is at most one monomer-preserving weak solution to (1.1)-(1.4) in the sense of Definition 1.1.
In the next section we derive suitable estimates on the primitive of the difference of two solutions. This first proves Theorem 2.1 and Theorem 2.2 which then entail Theorem 1.2 and Theorem 1.3.
3. Proofs
A Priori Estimates
Throughout this section we suppose that (2.1)-(2.14) are satisfied. Let be given and consider two monomer-preserving solutions and to (1.1)-(1.4) in the sense of Definition 1.1 such that
[TABLE]
Let us point out that then, in particular,
[TABLE]
where for and that (2.5) implies
[TABLE]
We now define
[TABLE]
Clearly, to prove that and coincide it suffices show that . Let us fix in the following. Note that
[TABLE]
and . Using integration by parts it then follows as in [10, Lemma 2.2] that satisfies the evolution equation
[TABLE]
for and . Therefore, introducing for and
[TABLE]
and
[TABLE]
it follows from , , that
[TABLE]
Our aim is to estimate
[TABLE]
and then apply Gronwall’s inequality to obtain that . This is implied by the subsequent lemmata.
Lemma 3.1**.**
There is such that
[TABLE]
Proof.
Choose in Definition 1.1 and recall from (1.15) the definition of . Then integration by parts yields
[TABLE]
Owing to (1.8) and (3.2) (applied to , see (3.1)), the boundary term at vanishes and we thus deduce with the help of (2.2), (2.6), and (2.7) that
[TABLE]
To estimate the last term we use (1.16) and obtain
[TABLE]
Integration by parts gives
[TABLE]
where the boundary term at vanishes since the monotonicity of and (3.1) imply
[TABLE]
Hence, together with (2.10) and (3.1) we deduce for
[TABLE]
From this along with (3.7), (3.8), and the symmetry of put us in a position to apply Gronwall’s inequality to conclude. ∎
We next estimate the difference between and .
Lemma 3.2**.**
There is such that
[TABLE]
Proof.
The proof is the same as in [10, Lemma 3.4] except that we have to treat also the case . It readily follows from (1.1), (2.4), and (3.1) that
[TABLE]
and hence
[TABLE]
Note that the boundary terms at vanish owing to (1.8), (2.4), (3.1), and (3.2). Thus, from (2.4) and (2.7) we deduce
[TABLE]
for , where we recall (2.14) for the case that . Gronwall’s lemma together with Lemma 3.1 yield the claim. ∎
For the first term in (3.5) we obtain:
Lemma 3.3**.**
If and , then
[TABLE]
where
[TABLE]
Proof.
This can be shown exactly as estimate (37) in [10] and using Lemma 3.1. ∎
For the second term in (3.5) we note:
Lemma 3.4**.**
If , then
[TABLE]
Proof.
We adapt parts of the proof of [5, Proposition 3.3]. Let . We set
[TABLE]
and note that, for and ,
[TABLE]
Then Fubini’s theorem (along with (2.12), (3.1)) and imply
[TABLE]
where the last inequality follows from the symmetry of . Thus, introducing
[TABLE]
we obtain
[TABLE]
Then (2.12), (3.1), and (3.10) entail
[TABLE]
For technical reasons we introduce for fixed the truncation
[TABLE]
of . As above, (2.12), (3.1), and (3.10) entail that
[TABLE]
Moreover, recalling that we have
[TABLE]
It then follows from (3.10) and (2.13) that
[TABLE]
while (2.11) entails that
[TABLE]
Consequently, (3.12), (3.13), and (3.1) imply that for each and . Hence, we can rewrite (3.11) in the form
[TABLE]
It follows from (3.12) that
[TABLE]
and, since is non-negative and non-decreasing,
[TABLE]
hence, by (3.1),
[TABLE]
Therefore, (3.14) implies for each the equality
[TABLE]
From (3.10) and (2.12) we obtain
[TABLE]
so that (3.1) guarantees
[TABLE]
Next, we invoke (2.13), (3.1), (3.3), (3.10) and apply Lebesgue’s theorem to get
[TABLE]
Moreover, we have
[TABLE]
and thus, by (2.11),
[TABLE]
We then pass to the limit in (3.16) and deduce from (3.10), (3.13) and (3.15)-(3.19) that
[TABLE]
The inequalities (2.11) and (2.13) yield together with (3.1) and Lemma 3.1 the assertion. ∎
We now obtain from Fatou’s lemma and (3.5), (3.6) that
[TABLE]
for . Lemma 3.3 and Lemma 3.4 then yield
[TABLE]
Using
[TABLE]
along with (2.14) when , it follows from (3.1), Lemma 3.1 and Lemma 3.2 that
[TABLE]
for . It then remains to show that for which we consider the cases of Theorem 2.1 and Theorem 2.2 separately.
Proof of Theorem 2.1
Let (2.1)-(2.16) be satisfied and given an initial value with and u^{0}\in L_{1}^{+}(Y,y\mathrm{d}y)\cap L_{1}\big{(}Y,(\mu+\beta)(y)G(y)\mathrm{d}y\big{)} consider two corresponding monomer-preserving weak solutions and to (1.1)-(1.4) in the sense of Definition 1.1 such that
[TABLE]
One then shows exactly as in the proof of [10, Theorem 3.1] (see equation (39) therein) that (2.15), (2.16) along with (3.22) imply
[TABLE]
Applying Gronwall’s lemma to inequality (3.21) yields and Theorem 2.1 follows.
Proof of Theorem 2.2
To prove Theorem 2.2 suppose (2.1)-(2.13) with (then (2.14) trivially holds) and (2.17)-(2.20). Given an initial value with and consider two corresponding monomer-preserving weak solutions and to (1.1)-(1.4) in the sense of Definition 1.1. Note that now , hence (3.1) holds. By Definition 1.1 (see in particular (1.12)) we have
[TABLE]
One then shows as in [10, Lemma 2.1] (using (1.12), (2.18)) that
[TABLE]
the only difference being that when testing (1.2) by with , an additional term
[TABLE]
comes in. However, this term tends to zero as by Lebesgue’s theorem since (2.17) together with (3.1) and (3.24) entail that
[TABLE]
Finally, exactly as in the proof of [10, Theorem 3.2], (3.24) and (3.25) imply (3.23). We may then apply again Gronwall’s inequality to (3.21) (with ) and conclude Theorem 2.2.
Proof of Theorem 1.2 and Theorem 1.3
Now, Theorem 1.2 and Theorem 1.3 are consequences of Theorem 2.1 respectively Theorem 2.2. Indeed, it was shown in the proof of [10, Theorem 1.2] that (1.17), (1.18) imply (2.1)-(2.9) and (2.15)-(2.16) when taking (to satisfy (2.14) one may take ) and that (2.18)-(2.20) hold under the assumptions of Theorem 1.3 while (2.17) follows from (1.19) with in this case. In [5, Lemma 3.4.] it was shown that (1.19), (1.20) imply (2.10)-(2.13) for these . This proves Theorem 1.2 and Theorem 1.3.
Acknowledgments
We thank the referee for carefully reading the paper and providing helpful comments.
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1[1] V. Calvez, N. Lenuzza, M. Doumic, J.-P. Deslys, F. Mouthon, B. Perthame. Prion dynamics with size dependency-strain phenomena. J. Biol. Dyn. 4 (2010), no. 1, 28–42.
- 2[2] V. Calvez, N. Lenuzza, D. Oelz, J.-P. Deslys, P. Laurent, F. Mouthon, B. Perthame. Size distribution dependence of prion aggregates infectivity. Math. Biosci. 217 (2009), no. 1, 89–99.
- 3[3] E. Cepeda. Well-posedness for a coagulation multiple-fragmenation equation. Differential Integral Equations 27 (2014),105–136.
- 4[4] M. Doumic, P. Gabriel. Eigenelements of a general aggregation-fragmentation model. Math. Models Methods Appl. Sci. 20 (2010), 757–783.
- 5[5] N. Fournier, Ph. Laurençot. Well-posedness of Smoluchowski’s coagulation equation for a class of homogeneous kernels. J. Funct. Analysis 233 (2006), 351–379.
- 6[6] H. Engler, J. Prüss, G. F. Webb. Analysis of a model for the dynamics of prions II. J. Math. Anal. Appl. 324 (2006), 98–117.
- 7[7] P. Gabriel. Global stability for the prion equation with general incidence. Math. Biosci. Eng. 12 (2015), no. 4, 789–801.
- 8[8] M. L. Greer, P. van den Driessche, L. Wang, G. F. Webb. Effects of general incidence and polymer joining on nucleated polymerzation in a model of prion proliferation. SIAM J. Appl. Math. 68 (2007), no. 1, 154–170.
