Initial condition of stochastic self-assembly

How and why to build a mathematical model: A case study using prion aggregation

Biological systems are inherently complex, and the increasing level of detail with which we are able to experimentally probe such systems continually reveals new complexity. Fortunately, mathematical models are uniquely positioned to provide a tool suitable for rigorous analysis, hypothesis generation, and connecting results from isolated in vitro experiments with results from in vivo and whole-organism studies. However, developing useful mathematical models is challenging because of the often different domains of knowledge required in both math and biology. In this work, we endeavor to provide a useful guide for researchers interested in incorporating mathematical modeling into their scientific process. We advocate for the use of conceptual diagrams as a starting place to anchor researchers from both domains. These diagrams are useful for simplifying the biological process in question and distinguishing the essential components. Not only do they serve as the basis for developing a variety of mathematical models, but they ensure that any mathematical formulation of the biological system is led primarily by scientific questions. We provide a specific example of this process from our own work in studying prion aggregation to show the power of mathematical models to synergistically interact with experiments and push forward biological understanding. Choosing the most suitable model also depends on many different factors, and we consider how to make these choices based on different scales of biological organization and available data. We close by discussing the many opportunities that abound for both experimentalists and modelers to take advantage of collaborative work in this field.

Evolution of Coagulation-Fragmentation Stochastic Processes Using Accurate Chemical Master Equation Approach

Coagulation and fragmentation (CF) is a fundamental process in which smaller particles attach to each other to form larger clusters while existing clusters break up into smaller particles . It is a ubiquitous process that plays important roles in many physical and biological phenomena. CF is typically a stochastic process that often occurs in confined spaces with a limited number of available particles . Here, we study the CF process formulated with the discrete Chemical Master Equation (dCME). Using the newly developed Accurate Chemical Master Equation (ACME) method, we examine the time-dependent behavior of the CF system. We investigate the effects of a number of important factors that influence the overall behavior of the system, including the dimensionality, the ratio of attachment to detachment rates among clusters, and the initial conditions. By comparing CF in one and three dimensions, we conclude that systems in three dimensions are more likely to form large clusters. We also demonstrate how the ratio of the attachment to detachment rates affects the dynamics and the steady-state of the system. Finally, we demonstrate the relationship between the formation of large clusters and the initial condition.

Generalizing a mathematical model of prion aggregation allows strain coexistence and co-stability by including a novel misfolded species

Prions are proteins capable of adopting misfolded conformations and transmitting these conformations to other normally folded proteins. Prions are most commonly known for causing fatal neurodegenerative diseases in mammals but are also associated with several harmless phenotypes in yeast. A distinct feature of prion propagation is the existence of different phenotypical variants, called strains. It is widely accepted that these strains correspond to different conformational states of the protein, but the mechanisms driving their interactions remain poorly understood. This study uses mathematical modeling to provide insight into this problem. We show that the classical model of prion dynamics allows at most one conformational strain to stably propagate. In order to conform to biological observations of strain coexistence and co-stability, we develop an extension of the classical model by introducing a novel prion species consistent with biological studies. Qualitative analysis of this model reveals a new variety of behavior. Indeed, it allows for stable coexistence of different strains in a wide parameter range, and it also introduces intricate initial condition dependency. These new behaviors are consistent with experimental observations of prions in both mammals and yeast. As such, our model provides a valuable tool for investigating the underlying mechanisms of prion propagation and the link between prion strains and strain specific phenotypes. The consideration of a novel prion species brings a change in perspective on prion biology and we use our model to generate hypotheses about prion infectivity.