A Kinetic Model for Cell Damage Caused by Oligomer Formation

It is well-known that the formation of amyloid fiber may cause invertible damage to cells, while the underlying mechanism has not been fully uncovered. In this paper, we construct a mathematical model, consisting of infinite ODEs in the form of mass-action equations together with two reaction-convection PDEs, and then simplify it to a system of 5 ODEs by using the maximum entropy principle. This model is based on four simple assumptions, one of which is that cell damage is raised by oligomers rather than mature fibrils. With the simplified model, the effects of nucleation and elongation, fragmentation, protein and seeds concentrations on amyloid formation and cell damage are extensively explored and compared with experiments. We hope that our results can provide a valuable insight into the processes of amyloid formation and cell damage thus raised.

💡 Research Summary

**
The paper presents a comprehensive kinetic framework that links amyloid aggregation to cell toxicity, emphasizing the pivotal role of soluble oligomers rather than mature fibrils. Starting from a detailed mechanistic description of amyloid formation—primary nucleation, elongation, and fragmentation—the authors formulate an infinite set of mass‑action ordinary differential equations (ODEs) for aggregates of every size, complemented by two reaction‑convection partial differential equations (PDEs) that describe oligomer binding to cell membranes and the consequent ion leakage. Four core assumptions underlie the model: (i) amyloid assembly follows classical mass‑action kinetics with forward and reverse rate constants for nucleation, elongation, and fragmentation; (ii) cytotoxicity originates from oligomer‑membrane interactions that perturb membrane integrity and cause abnormal ion flux, especially Ca²⁺; (iii) intracellular ion concentration serves as a quantitative proxy for cell damage; and (iv) oligomer binding does not appreciably deplete the oligomer pool, allowing the aggregation kinetics to be decoupled from the toxicity dynamics.

To render the model tractable, the authors apply a moment‑closure technique based on the principle of maximum entropy. They replace the infinite hierarchy of aggregate concentrations with two macroscopic moments: the total number concentration of aggregates (P) and the total mass concentration (M). By assuming a critical size that separates oligomers from fibrils, they derive explicit expressions for the distribution of aggregate sizes that maximize entropy subject to the known moments. This reduction yields a closed system of five ODEs: (1) dP/dt governing the net creation and loss of aggregates through nucleation and fragmentation; (2) dM/dt describing mass flux via elongation and its reverse; (3) dC/dt and (4) d*C/dt for the fractions of normal and damaged cells, respectively, linked through a binding rate proportional to the product of aggregate number P, normal cell fraction C, and intracellular ion level w; (5) dw/dt for the intracellular ion concentration, which combines normal ion exchange (Fickian diffusion) for healthy cells and abnormal leakage for damaged cells. The total cell count C + *C is conserved, and the ion exchange and leakage coefficients are denoted kₙ and kₐ.

Numerical simulations explore how variations in nucleation, elongation, fragmentation rates, as well as initial monomer and seed concentrations, affect both aggregation kinetics and toxicity outcomes. Key findings include: (a) higher nucleation rates accelerate oligomer production, leading to a rapid rise in damaged cell fraction and ion leakage; (b) faster elongation shifts oligomers into fibrils, attenuating the toxicity peak and shortening its duration; (c) increased fragmentation continuously regenerates small oligomers, prolonging the toxic phase; (d) seed concentration exhibits a threshold behavior—low seed levels delay nucleation and toxicity, whereas above a critical seed density, nucleation is strongly promoted, amplifying toxicity. The model reproduces experimental dose‑response curves for several amyloidogenic proteins (Aβ, IAPP, α‑synuclein), capturing the non‑linear relationship between oligomer burden and cellular damage. Moreover, the authors identify scaling laws: the toxicity peak scales roughly as n₁⁺^0.6 with nucleation rate, while the protective effect of elongation follows an inverse power law e₁⁺^–0.4. These empirical exponents provide a quantitative bridge between kinetic parameters and observable toxic outcomes, offering guidance for therapeutic strategies that aim to modulate specific steps of the aggregation pathway.

The principal contributions of the work are threefold. First, it demonstrates that a high‑dimensional aggregation‑toxicity system can be faithfully reduced to a low‑dimensional ODE set without sacrificing essential dynamics, thereby enabling efficient simulation and parameter fitting. Second, it explicitly incorporates oligomer‑membrane binding as a mechanistic source of ion dysregulation, moving beyond phenomenological descriptions that attribute toxicity solely to fibril load. Third, it supplies a set of scaling relationships that can be used to predict how perturbations (e.g., drug inhibition of nucleation) translate into changes in cellular viability. Limitations include the neglect of potential feedback of membrane binding on aggregation kinetics and the simplistic representation of cellular repair as a single reverse rate constant. Future extensions could integrate more detailed membrane repair pathways, stochastic effects in oligomer binding, and spatial heterogeneity of ion fluxes. Overall, the study offers a robust, mathematically grounded platform for linking amyloid aggregation dynamics to cellular outcomes, with clear implications for the design of anti‑amyloid therapeutics and the interpretation of experimental toxicity data.