Probabilistic and maximum entropy modeling of chemical reaction systems: Characteristics and comparisons to mass action kinetic models.

We demonstrate and characterize a first-principles approach to modeling the mass action dynamics of metabolism. Starting from a basic definition of entropy expressed as a multinomial probability density using Boltzmann probabilities with standard chemical potentials, we derive and compare the free energy dissipation and the entropy production rates. We express the relation between entropy production and the chemical master equation for modeling metabolism, which unifies chemical kinetics and chemical thermodynamics. Because prediction uncertainty with respect to parameter variability is frequently a concern with mass action models utilizing rate constants, we compare and contrast the maximum entropy model, which has its own set of rate parameters, to a population of standard mass action models in which the rate constants are randomly chosen. We show that a maximum entropy model is characterized by a high probability of free energy dissipation rate and likewise entropy production rate, relative to other models. We then characterize the variability of the maximum entropy model predictions with respect to uncertainties in parameters (standard free energies of formation) and with respect to ionic strengths typically found in a cell.

Role of turgor-pressure induced boundary tension in the maintenance of the shoot apical meristem of Arabidopsis thaliana.

In plants, the robust maintenance of tissue structure is crucial to supporting its functionality. The multi-layered shoot apical meristem (SAM) of Arabidopsis, containing stem cells, is an approximately radially symmetric tissue whose shape and structure is maintained throughout the life of the plant. In this paper, a new biologically calibrated pseudo-three-dimensional (P3D) computational model of a longitudinal section of the SAM is developed. It includes anisotropic expansion and division of cells out of the cross-section plane, as well as representation of tension experienced by the SAM epidermis. Results from the experimentally calibrated P3D model provide new insights into maintenance of the structure of the SAM epidermal cell monolayer under tension and quantify dependence of epidermal and subepidermal cell anisotropy on the amount of tension. Moreover, the model simulations revealed that out-of-plane cell growth is important in offsetting cell crowding and regulating mechanical stresses experienced by tunica cells. Predictive model simulations show that tension-determined cell division plane orientation in the apical corpus may be regulating cell and tissue shape distributions needed for maintaining structure of the wild-type SAM. This suggests that cells' responses to local mechanical cues may serve as a mechanism to regulate cell- and tissue-scale patterning.

Combined computational modeling and experimental analysis integrating chemical and mechanical signals suggests possible mechanism of shoot meristem maintenance.

Stem cell maintenance in multilayered shoot apical meristems (SAMs) of plants requires strict regulation of cell growth and division. Exactly how the complex milieu of chemical and mechanical signals interact in the central region of the SAM to regulate cell division plane orientation is not well understood. In this paper, simulations using a newly developed multiscale computational model are combined with experimental studies to suggest and test three hypothesized mechanisms for the regulation of cell division plane orientation and the direction of anisotropic cell expansion in the corpus. Simulations predict that in the Apical corpus, WUSCHEL and cytokinin regulate the direction of anisotropic cell expansion, and cells divide according to tensile stress on the cell wall. In the Basal corpus, model simulations suggest dual roles for WUSCHEL and cytokinin in regulating both the direction of anisotropic cell expansion and cell division plane orientation. Simulation results are followed by a detailed analysis of changes in cell characteristics upon manipulation of WUSCHEL and cytokinin in experiments that support model predictions. Moreover, simulations predict that this layer-specific mechanism maintains both the experimentally observed shape and structure of the SAM as well as the distribution of WUSCHEL in the tissue. This provides an additional link between the roles of WUSCHEL, cytokinin, and mechanical stress in regulating SAM growth and proper stem cell maintenance in the SAM.

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.

Cell-Based Model of the Generation and Maintenance of the Shape and Structure of the Multilayered Shoot Apical Meristem of Arabidopsis thaliana.

One of the central problems in animal and plant developmental biology is deciphering how chemical and mechanical signals interact within a tissue to produce organs of defined size, shape, and function. Cell walls in plants impose a unique constraint on cell expansion since cells are under turgor pressure and do not move relative to one another. Cell wall extensibility and constantly changing distribution of stress on the wall are mechanical properties that vary between individual cells and contribute to rates of expansion and orientation of cell division. How exactly cell wall mechanical properties influence cell behavior is still largely unknown. To address this problem, a novel, subcellular element computational model of growth of stem cells within the multilayered shoot apical meristem (SAM) of Arabidopsis thaliana is developed and calibrated using experimental data. Novel features of the model include separate, detailed descriptions of cell wall extensibility and mechanical stiffness, deformation of the middle lamella, and increase in cytoplasmic pressure generating internal turgor pressure. The model is used to test novel hypothesized mechanisms of formation of the shape and structure of the growing, multilayered SAM based on WUS concentration of individual cells controlling cell growth rates and layer-dependent anisotropic mechanical properties of subcellular components of individual cells determining anisotropic cell expansion directions. Model simulations also provide a detailed prediction of distribution of stresses in the growing tissue which can be tested in future experiments.