Open problems in mathematical biology

Design Principles for Biological Adaptation: A Systems and Control-Theoretic Treatment

Establishing a mapping between (from and to) the functionality of interest and the underlying network structure (design principles) remains a crucial step toward understanding and design of bio-systems. Perfect adaptation is one such crucial functionality that enables every living organism to regulate its essential activities in the presence of external disturbances. Previous approaches to deducing the design principles for adaptation have either relied on computationally burdensome brute-force methods or rule-based design strategies detecting only a subset of all possible adaptive network structures. This chapter outlines a scalable and generalizable method inspired by systems theory that unravels an exhaustive set of adaptation-capable structures. We first use the well-known performance parameters to characterize perfect adaptation. These performance parameters are then mapped back to a few parameters (poles, zeros, gain) characteristic of the underlying dynamical system constituted by the rate equations. Therefore, the performance parameters evaluated for the scenario of perfect adaptation can be expressed as a set of precise mathematical conditions involving the system parameters. Finally, we use algebraic graph theory to translate these abstract mathematical conditions to certain structural requirements for adaptation. The proposed algorithm does not assume any particular dynamics and is applicable to networks of any size. Moreover, the results offer a significant advancement in the realm of understanding and designing complex biochemical networks.

Review: The prevailing mathematical modeling classifications and paradigms to support the advancement of sustainable animal production

Mathematical modeling is typically framed as the art of reductionism of scientific knowledge into an arithmetical layout. However, most untrained people get the art of modeling wrong and end up neglecting it because modeling is not simply about writing equations and generating numbers through simulations. Models tell not only about a story; they are spoken to by the circumstances under which they are envisioned. They guide apprentice and experienced modelers to build better models by preventing known pitfalls and invalid assumptions in the virtual world and, most importantly, learn from them through simulation and identify gaps in pushing scientific knowledge further. The power of the human mind is well-documented for idealizing concepts and creating virtual reality models, and as our hypotheses grow more complicated and more complex data become available, modeling earns more noticeable footing in biological sciences. The fundamental modeling paradigms include discrete-events, dynamic systems, agent-based (AB), and system dynamics (SD). The source of knowledge is the most critical step in the model-building process regardless of the paradigm, and the necessary expertise includes (a) clear and concise mental concepts acquired through different ways that provide the fundamental structure and expected behaviors of the model and (b) numerical data necessary for statistical analysis, not for building the model. The unreasonable effectiveness of models to grow scientific learning and knowledge in sciences arise because different researchers would model the same problem differently, given their knowledge and experiential background, leading to choosing different variables and model structures. Secondly, different researchers might use different paradigms and even unalike mathematics to resolve the same problem; thus, model needs are intrinsic to their perceived assumptions and structures. Thirdly, models evolve as the scientific community knowledge accumulates and matures over time, hopefully resulting in improved modeling efforts; thus, the perfect model is fictional. Some paradigms are most appropriate for macro, high abstraction with less detailed-oriented scenarios, while others are most suitable for micro, low abstraction with higher detailed-oriented strategies. Modern hybridization aggregating artificial intelligence (AI) to mathematical models can become the next technological wave in modeling. AI can be an integral part of the SD/AB models and, before long, write the model code by itself. Success and failures in model building are more related to the ability of the researcher to interpret the data and understand the underlying principles and mechanisms to formulate the correct relationship among variables rather than profound mathematical knowledge.

Distilling identifiable and interpretable dynamic models from biological data

Mechanistic dynamical models allow us to study the behavior of complex biological systems. They can provide an objective and quantitative understanding that would be difficult to achieve through other means. However, the systematic development of these models is a non-trivial exercise and an open problem in computational biology. Currently, many research efforts are focused on model discovery, i.e. automating the development of interpretable models from data. One of the main frameworks is sparse regression, where the sparse identification of nonlinear dynamics (SINDy) algorithm and its variants have enjoyed great success. SINDy-PI is an extension which allows the discovery of rational nonlinear terms, thus enabling the identification of kinetic functions common in biochemical networks, such as Michaelis-Menten. SINDy-PI also pays special attention to the recovery of parsimonious models (Occam's razor). Here we focus on biological models composed of sets of deterministic nonlinear ordinary differential equations. We present a methodology that, combined with SINDy-PI, allows the automatic discovery of structurally identifiable and observable models which are also mechanistically interpretable. The lack of structural identifiability and observability makes it impossible to uniquely infer parameter and state variables, which can compromise the usefulness of a model by distorting its mechanistic significance and hampering its ability to produce biological insights. We illustrate the performance of our method with six case studies. We find that, despite enforcing sparsity, SINDy-PI sometimes yields models that are unidentifiable. In these cases we show how our method transforms their equations in order to obtain a structurally identifiable and observable model which is also interpretable.

On the Analytical Solution of Fractional SIR Epidemic Model

This article presents the solution of the fractional SIR epidemic model using the Laplace residual power series method. We introduce the fractional SIR model in the sense of Caputo’s derivative; it is presented by three fractional differential equations, in which the third one depends on the first coupled equations. The Laplace residual power series method (LRPSM) is implemented in this research to solve the proposed model, in which we present the solution in a form of convergent series expansion that converges rapidly to the exact one. We analyze the results and compare the obtained approximate solutions to those obtained from other methods. Figures and tables are illustrated to show the efficiency of the LRPSM in handling the proposed SIR model.

Bayesian parameter estimation for dynamical models in systems biology

Dynamical systems modeling, particularly via systems of ordinary differential equations, has been used to effectively capture the temporal behavior of different biochemical components in signal transduction networks. Despite the recent advances in experimental measurements, including sensor development and '-omics' studies that have helped populate protein-protein interaction networks in great detail, modeling in systems biology lacks systematic methods to estimate kinetic parameters and quantify associated uncertainties. This is because of multiple reasons, including sparse and noisy experimental measurements, lack of detailed molecular mechanisms underlying the reactions, and missing biochemical interactions. Additionally, the inherent nonlinearities with respect to the states and parameters associated with the system of differential equations further compound the challenges of parameter estimation. In this study, we propose a comprehensive framework for Bayesian parameter estimation and complete quantification of the effects of uncertainties in the data and models. We apply these methods to a series of signaling models of increasing mathematical complexity. Systematic analysis of these dynamical systems showed that parameter estimation depends on data sparsity, noise level, and model structure, including the existence of multiple steady states. These results highlight how focused uncertainty quantification can enrich systems biology modeling and enable additional quantitative analyses for parameter estimation.