Technologies for whole-cell modeling: Genome-wide reconstruction of a cell in silico

With advances in high-throughput, large-scale in vivo measurement and genome modification techniques at the single-nucleotide level, there is an increasing demand for the development of new technologies for the flexible design and control of cellular systems. Computer-aided design is a powerful tool to design new cells. Whole-cell modeling aims to integrate various cellular subsystems, determine their interactions and cooperative mechanisms, and predict comprehensive cellular behaviors by computational simulations on a genome-wide scale. It has been applied to prokaryotes, yeasts, and higher eukaryotic cells, and utilized in a wide range of applications, including production of valuable substances, drug discovery, and controlled differentiation. Whole-cell modeling, consisting of several thousand elements with diverse scales and properties, requires innovative model construction, simulation, and analysis techniques. Furthermore, whole-cell modeling has been extended to multiple scales, including high-resolution modeling at the single-nucleotide and single-amino acid levels and multicellular modeling of tissues and organs. This review presents an overview of the current state of whole-cell modeling, discusses the novel computational and experimental technologies driving it, and introduces further developments toward multihierarchical modeling on a whole-genome scale.

Development and Analysis of Multiscale Models for Tuberculosis: From Molecules to Populations

Although infectious disease dynamics are often analyzed at the macro-scale, increasing numbers of drug-resistant infections highlight the importance of within-host modeling that simultaneously solves across multiple scales to effectively respond to epidemics. We review multiscale modeling approaches for complex, interconnected biological systems and discuss critical steps involved in building, analyzing, and applying such models within the discipline of model credibility. We also present our two tools: CaliPro, for calibrating multiscale models (MSMs) to datasets, and tunable resolution, for fine- and coarse-graining sub-models while retaining insights. We include as an example our work simulating infection with Mycobacterium tuberculosis to demonstrate modeling choices and how predictions are made to generate new insights and test interventions. We discuss some of the current challenges of incorporating novel datasets, rigorously training computational biologists, and increasing the reach of MSMs. We also offer several promising future research directions of incorporating within-host dynamics into applications ranging from combinatorial treatment to epidemic response.

Generative models of morphogenesis in developmental biology

Understanding the mechanism by which cells coordinate their differentiation and migration is critical to our understanding of many fundamental processes such as wound healing, disease progression, and developmental biology. Mathematical models have been an essential tool for testing and developing our understanding, such as models of cells as soft spherical particles, reaction-diffusion systems that couple cell movement to environmental factors, and multi-scale multi-physics simulations that combine bottom-up rule-based models with continuum laws. However, mathematical models can often be loosely related to data or have so many parameters that model behaviour is weakly constrained. Recent methods in machine learning introduce new means by which models can be derived and deployed. In this review, we discuss examples of mathematical models of aspects of developmental biology, such as cell migration, and how these models can be combined with these recent machine learning methods.

Hybrid computational models of multicellular tumour growth considering glucose metabolism

Cancer cells metabolize glucose through metabolic pathways that differ from those used by healthy and differentiated cells. In particular, tumours have been shown to consume more glucose than their healthy counterparts and to use anaerobic metabolic pathways, even under aerobic conditions. Nevertheless, scientists have still not been able to explain why cancer cells evolved to present an altered metabolism and what evolutionary advantage this might provide them. Experimental and computational models have been increasingly used in recent years to understand some of these biological questions. Multicellular tumour spheroids are effective experimental models as they replicate the initial stages of avascular solid tumour growth. Furthermore, these experiments generate data which can be used to calibrate and validate computational studies that aim to simulate tumour growth. Hybrid models are of particular relevance in this field of research because they model cells as individual agents while also incorporating continuum representations of the substances present in the surrounding microenvironment that may participate in intracellular metabolic networks as concentration or density distributions. Henceforth, in this review, we explore the potential of computational modelling to reveal the role of metabolic reprogramming in tumour growth.

Alternative model systems for cognitive variation: eusocial-insect colonies

Understanding the origins and maintenance of cognitive variation in animal populations is central to the study of the evolution of cognition. However, the brain is itself a complex, hierarchical network of heterogeneous components, from diverse cell types to diverse neuropils, each of which may be of limited use to study in isolation or prohibitively challenging to manipulate in situ. Consequently, highly tractable alternative model systems may be valuable tools. Eusocial-insect colonies display emergent cognitive-like properties from relatively simple social interactions between diverse subunits that can be observed and manipulated while operating collectively. Here, we review the individual-scale mechanisms that cause group-level variation in how colonies solve problems analogous to cognitive challenges faced by brains, like decision-making, attention, and search.

The Energy Homeostasis Principle: A Naturalistic Approach to Explain the Emergence of Behavior

It is still elusive to explain the emergence of behavior and understanding based on its neural mechanisms. One renowned proposal is the Free Energy Principle (FEP), which uses an information-theoretic framework derived from thermodynamic considerations to describe how behavior and understanding emerge. FEP starts from a whole-organism approach, based on mental states and phenomena, mapping them into the neuronal substrate. An alternative approach, the Energy Homeostasis Principle (EHP), initiates a similar explanatory effort but starts from single-neuron phenomena and builds up to whole-organism behavior and understanding. In this work, we further develop the EHP as a distinct but complementary vision to FEP and try to explain how behavior and understanding would emerge from the local requirements of the neurons. Based on EHP and a strict naturalist approach that sees living beings as physical and deterministic systems, we explain scenarios where learning would emerge without the need for volition or goals. Given these starting points, we state several considerations of how we see the nervous system, particularly the role of the function, purpose, and conception of goal-oriented behavior. We problematize these conceptions, giving an alternative teleology-free framework in which behavior and, ultimately, understanding would still emerge. We reinterpret neural processing by explaining basic learning scenarios up to simple anticipatory behavior. Finally, we end the article with an evolutionary perspective of how this non-goal-oriented behavior appeared. We acknowledge that our proposal, in its current form, is still far from explaining the emergence of understanding. Nonetheless, we set the ground for an alternative neuron-based framework to ultimately explain understanding.

Towards the First Principles in Biology and Cancer: New Vistas in Computational Systems Biology of Cancer

These days many leading scientists argue for a new paradigm for cancer research and propose a complex systems-view of cancer supported by empirical evidence. As an example, Thea Newman (2021) has applied "the lessons learned from physical systems to a critique of reductionism in medical research, with an emphasis on cancer". It is the understanding of this author that the mesoscale constructs that combine the bottom-up as well as top-down approaches, are very close to the concept of emergence. The mesoscale constructs can be said to be those effective components through which the system allows itself to be understood. A short list of basic concepts related to life/biology fundamentals are first introduced to demonstrate a lack of emphasis on these matters in literature. It is imperative that physical and chemical approaches are introduced and incorporated in biology to make it more conceptually sound, quantitative, and based on the first principles. Non-equilibrium thermodynamics is the only tool currently available for making progress in this direction. A brief outline of systems biology, the discovery of emergent properties, and metabolic modeling are introduced in the second part. Then, different cancer initiation concepts are reviewed, followed by application of non-equilibrium thermodynamics in the metabolic and genomic analysis of initiation and development of cancer, stressing the endogenous network hypothesis (ENH). Finally, extension of the ENH is suggested to include a cancer niche (exogenous network hypothesis). It is expected that this will lead to a unifying systems-biology approach for a future combination of the analytical and synthetic arms of two major hypotheses of cancer models (SMT and TOFT).

Recent applications of quantitative systems pharmacology and machine learning models across diseases

Quantitative systems pharmacology (QSP) is a quantitative and mechanistic platform describing the phenotypic interaction between drugs, biological networks, and disease conditions to predict optimal therapeutic response. In this meta-analysis study, we review the utility of the QSP platform in drug development and therapeutic strategies based on recent publications (2019-2021). We gathered recent original QSP models and described the diversity of their applications based on therapeutic areas, methodologies, software platforms, and functionalities. The collection and investigation of these publications can assist in providing a repository of recent QSP studies to facilitate the discovery and further reusability of QSP models. Our review shows that the largest number of QSP efforts in recent years is in Immuno-Oncology. We also addressed the benefits of integrative approaches in this field by presenting the applications of Machine Learning methods for drug discovery and QSP models. Based on this meta-analysis, we discuss the advantages and limitations of QSP models and propose fields where the QSP approach constitutes a valuable interface for more investigations to tackle complex diseases and improve drug development.

Characteristics of mathematical modeling languages that facilitate model reuse in systems biology: a software engineering perspective

Reuse of mathematical models becomes increasingly important in systems biology as research moves toward large, multi-scale models composed of heterogeneous subcomponents. Currently, many models are not easily reusable due to inflexible or confusing code, inappropriate languages, or insufficient documentation. Best practice suggestions rarely cover such low-level design aspects. This gap could be filled by software engineering, which addresses those same issues for software reuse. We show that languages can facilitate reusability by being modular, human-readable, hybrid (i.e., supporting multiple formalisms), open, declarative, and by supporting the graphical representation of models. Modelers should not only use such a language, but be aware of the features that make it desirable and know how to apply them effectively. For this reason, we compare existing suitable languages in detail and demonstrate their benefits for a modular model of the human cardiac conduction system written in Modelica.

Advancing crop genomics from lab to field

Crop genomics remains a key element in ensuring scientific progress to secure global food security. It has been two decades since the sequence of the first plant genome, that of Arabidopsis thaliana, was released, and soon after that the draft sequencing of the rice genome was completed. Since then, the genomes of more than 100 crops have been sequenced, plant genome research has expanded across multiple fronts and the next few years promise to bring further advances spurred by the advent of new technologies and approaches. We are likely to see continued innovations in crop genome sequencing, genetic mapping and the acquisition of multiple levels of biological data. There will be exciting opportunities to integrate genome-scale information across multiple scales of biological organization, leading to advances in our mechanistic understanding of crop biological processes, which will, in turn, provide greater impetus for translation of laboratory results to the field.

A Holistic Approach to Study Photosynthetic Acclimation Responses of Plants to Fluctuating Light

Plants in natural environments receive light through sunflecks, the duration and distribution of these being highly variable across the day. Consequently, plants need to adjust their photosynthetic processes to avoid photoinhibition and maximize yield. Changes in the composition of the photosynthetic apparatus in response to sustained changes in the environment are referred to as photosynthetic acclimation, a process that involves changes in protein content and composition. Considering this definition, acclimation differs from regulation, which involves processes that alter the activity of individual proteins over short-time periods, without changing the abundance of those proteins. The interconnection and overlapping of the short- and long-term photosynthetic responses, which can occur simultaneously or/and sequentially over time, make the study of long-term acclimation to fluctuating light in plants challenging. In this review we identify short-term responses of plants to fluctuating light that could act as sensors and signals for acclimation responses, with the aim of understanding how plants integrate environmental fluctuations over time and tailor their responses accordingly. Mathematical modeling has the potential to integrate physiological processes over different timescales and to help disentangle short-term regulatory responses from long-term acclimation responses. We review existing mathematical modeling techniques for studying photosynthetic responses to fluctuating light and propose new methods for addressing the topic from a holistic point of view.

High-speed parameter search of dynamic biological pathways from time-course transcriptomic profiles using high-level Petri net

Dynamic simulation promises a deeper understanding of complex molecular mechanisms of biological pathways. How to determine the reaction kinetic parameters which govern the simulation results is still an open question in the field of systems biology. (1) Background: To execute simulation experiments, it is an essential first step to search effective values of model parameters. The complexity of biological systems and the experimental measurement technology severely limit the acquirement of accurate kinetic parameters. Previously proposed genomic data assimilation (GDA) approach enables users to handle parameter estimation using time-course information. However, it highly depends on successive time points and costs massive computational resource; (2) Methods: To address this problem, we present a new high-speed parameter search method for estimating the kinetic parameters of quantitative biological pathways using time-course transcriptomic profiles. The key idea of our method is to interactively prune the search space by introducing Probabilistic Linear-time Temporal Logic (PLTL) based model checking into GDA. (3) Results and conclusion: We demonstrated the effectiveness of our method by comparing with GDA on Mus musculus transcription circuits modelled by hybrid functional Petri net with extension. As a result, our method works faster and more accurate than GDA for both time-course datasets with dense and sparse observed values.

The potential for complex computational models of aging

The gradual accumulation of damage and dysregulation during the aging of living organisms can be quantified. Even so, the aging process is complex and has multiple interacting physiological scales - from the molecular to cellular to whole tissues. In the face of this complexity, we can significantly advance our understanding of aging with the use of computational models that simulate realistic individual trajectories of health as well as mortality. To do so, they must be systems-level models that incorporate interactions between measurable aspects of age-associated changes. To incorporate individual variability in the aging process, models must be stochastic. To be useful they should also be predictive, and so must be fit or parameterized by data from large populations of aging individuals. In this perspective, we outline where we have been, where we are, and where we hope to go with such computational models of aging. Our focus is on data-driven systems-level models, and on their great potential in aging research.

V-Mango: a functional-structural model of mango tree growth, development and fruit production

BACKGROUND AND AIMS

Mango (Mangifera indica L.) is the fifth most widely produced fruit in the world. Its cultivation, mainly in tropical and sub-tropical regions, raises a number of issues such as the irregular fruit production across years, phenological asynchronisms that lead to long periods of pest and disease susceptibility, and the heterogeneity of fruit quality and maturity at harvest. To address these issues, we developed an integrative functional-structural plant model that synthesizes knowledge about the vegetative and reproductive development of the mango tree and opens up the possible simulation of cultivation practices.

METHODS

We designed a model of architectural development in order to precisely characterize the intricate developmental processes of the mango tree. The appearance of botanical entities was decomposed into elementary stochastic events describing occurrence, intensity and timing of development. These events were determined by structural (position and fate of botanical entities) and temporal (appearance dates) factors. Daily growth and development of growth units and inflorescences were modelled using empirical distributions and thermal time. Fruit growth was determined using an ecophysiological model that simulated carbon- and water-related processes at the fruiting branch scale.

KEY RESULTS

The model simulates the dynamics of the population of growth units, inflorescences and fruits at the tree scale during a growing cycle. Modelling the effects of structural and temporal factors makes it possible to simulate satisfactorily the complex interplays between vegetative and reproductive development. The model allowed the characterization of the susceptibility of mango tree to pests and the investigatation of the influence of tree architecture on fruit growth.

CONCLUSIONS

This integrative functional-structural model simulates mango tree vegetative and reproductive development over successive growing cycles, allowing a precise characterization of tree phenology and fruit growth and production. The next step is to integrate the effects of cultivation practices, such as pruning, into the model.

Advances in Astrocyte Computational Models: From Metabolic Reconstructions to Multi-omic Approaches

The growing importance of astrocytes in the field of neuroscience has led to a greater number of computational models devoted to the study of astrocytic functions and their metabolic interactions with neurons. The modeling of these interactions demands a combined understanding of brain physiology and the development of computational frameworks based on genomic-scale reconstructions, system biology, and dynamic models. These computational approaches have helped to highlight the neuroprotective mechanisms triggered by astrocytes and other glial cells, both under normal conditions and during neurodegenerative processes. In the present review, we evaluate some of the most relevant models of astrocyte metabolism, including genome-scale reconstructions and astrocyte-neuron interactions developed in the last few years. Additionally, we discuss novel strategies from the multi-omics perspective and computational models of other glial cell types that will increase our knowledge in brain metabolism and its association with neurodegenerative diseases.

Systems biology and big data in asthma and allergy: recent discoveries and emerging challenges

Asthma is a common condition caused by immune and respiratory dysfunction, and it is often linked to allergy. A systems perspective may prove helpful in unravelling the complexity of asthma and allergy. Our aim is to give an overview of systems biology approaches used in allergy and asthma research. Specifically, we describe recent "omic"-level findings, and examine how these findings have been systematically integrated to generate further insight.Current research suggests that allergy is driven by genetic and epigenetic factors, in concert with environmental factors such as microbiome and diet, leading to early-life disturbance in immunological development and disruption of balance within key immuno-inflammatory pathways. Variation in inherited susceptibility and exposures causes heterogeneity in manifestations of asthma and other allergic diseases. Machine learning approaches are being used to explore this heterogeneity, and to probe the pathophysiological patterns or "endotypes" that correlate with subphenotypes of asthma and allergy. Mathematical models are being built based on genomic, transcriptomic and proteomic data to predict or discriminate disease phenotypes, and to describe the biomolecular networks behind asthma.The use of systems biology in allergy and asthma research is rapidly growing, and has so far yielded fruitful results. However, the scale and multidisciplinary nature of this research means that it is accompanied by new challenges. Ultimately, it is hoped that systems medicine, with its integration of omics data into clinical practice, can pave the way to more precise, personalised and effective management of asthma.

Exploring cardiac form and function: A length-scale computational biology approach

The ability to adequately pump blood throughout the body is the result of tightly regulated feedback mechanisms that exist across many spatial scales in the heart. Diseases which impede the function at any one of the spatial scales can cause detrimental cardiac remodeling and eventual heart failure. An overarching goal of cardiac research is to use engineered heart tissue in vitro to study the physiology of diseased heart tissue, develop cell replacement therapies, and explore drug testing applications. A commonality within the field is to manipulate the flow of mechanical signals across the various spatial scales to direct self-organization and build functional tissue. Doing so requires an understanding of how chemical, electrical, and mechanical cues can be used to alter the cellular microenvironment. We discuss how mathematical models have been used in conjunction with experimental techniques to explore various structure-function relations that exist across numerous spatial scales. We highlight how a systems biology approach can be employed to recapitulate in vivo characteristics in vitro at the tissue, cell, and subcellular scales. Specific focus is placed on the interplay between experimental and theoretical approaches. Various modeling methods are showcased to demonstrate the breadth and power afforded to the systems biology approach. An overview of modeling methodologies exemplifies how the strengths of different scientific disciplines can be used to supplement and/or inspire new avenues of experimental exploration. This article is categorized under: Models of Systems Properties and Processes > Mechanistic Models Models of Systems Properties and Processes > Cellular Models Models of Systems Properties and Processes > Organ, Tissue, and Physiological Models.

Reliability-based topology design for large-scale networks

Network topology has a critical effect on its reliability. Due to the complexity of networks, how to design a more reliable network topology is one of the fascinating and significant challenges. In this paper, inspired by the multi-scale model of biological systems which range from gene to individual cells, and up to the individual organism, we provide a novel method to design a reliable topology of large-scale networks. By means of fractal theory and its application in complex networks, a network topology can be described as consisting of the superposition of basic units, called fractal-cells, which are made up of several components and their connections. Based on this fractal-cell structure, we find that a fractal-cell network has structural similarity with its primitive structure, and then develop a method to build a reliable topology which is generated from a primitive ring structure. Numerical simulation compared with the existing methods and performing on real networks show that our proposed method is effective.

Global sensitivity analysis of biological multi-scale models

Mathematical models of biological systems need to both reflect and manage the inherent complexities of biological phenomena. Through their versatility and ability to capture behavior at multiple scales, multi-scale models offer a valuable approach. Due to the typically nonlinear and stochastic nature of multi-scale models as well as unknown parameter values, various types of uncertainty are present; thus, effective assessment and quantification of such uncertainty through sensitivity analysis is important. In this review, we discuss global sensitivity analysis in the context of multi-scale and multi-compartment models and highlight its value in model development and analysis. We present an overview of sensitivity analysis methods, approaches for extending such methods to a multi-scale setting, and examples of how sensitivity analysis can inform model reduction. Through schematics and references to past work, we aim to emphasize the advantages and usefulness of such techniques.

Modeling heterogeneous tumor growth dynamics and cell-cell interactions at single-cell and cell-population resolution

Cancer is a complex, dynamic disease that despite recent advances remains mostly incurable. Inter- and intratumoral heterogeneity are generally considered major drivers of therapy resistance, metastasis, and treatment failure. Recent advances in high-throughput experimentation have produced a wealth of data on tumor heterogeneity and researchers are increasingly turning to mathematical modeling to aid in the interpretation of these complex datasets. In this mini-review, we discuss three important classes of approaches for modeling cellular dynamics within heterogeneous tumors: agent-based models, population dynamics, and multiscale models. An important new focus, for which we provide an example, is the role of intratumoral cell-cell interactions.

Application of Modeling Approaches to Explore Vaccine Adjuvant Mode-of-Action

Novel adjuvant technologies have a key role in the development of next-generation vaccines, due to their capacity to modulate the duration, strength and quality of the immune response. The AS01 adjuvant is used in the malaria vaccine RTS,S/AS01 and in the licensed herpes-zoster vaccine (Shingrix) where the vaccine has proven its ability to generate protective responses with both robust humoral and T-cell responses. For many years, animal models have provided insights into adjuvant mode-of-action (MoA), generally through investigating individual genes or proteins. Furthermore, modeling and simulation techniques can be utilized to integrate a variety of different data types; ranging from serum biomarkers to large scale “omics” datasets. In this perspective we present a framework to create a holistic integration of pre-clinical datasets and immunological literature in order to develop an evidence-based hypothesis of AS01 adjuvant MoA, creating a unified view of multiple experiments. Furthermore, we highlight how holistic systems-knowledge can serve as a basis for the construction of models and simulations supporting exploration of key questions surrounding adjuvant MoA. Using the Systems-Biology-Graphical-Notation, a tool for graphical representation of biological processes, we have captured high-level cellular behaviors and interactions, and cytokine dynamics during the early immune response, which are substantiated by a series of diagrams detailing cellular dynamics. Through explicitly describing AS01 MoA we have built a consensus of understanding across multiple experiments, and so we present a framework to integrate modeling approaches into exploring adjuvant MoA, in order to guide experimental design, interpret results and inform rational design of vaccines.

Integrative Computational Modeling of the Lymph Node Stromal Cell Landscape

Adaptive immune responses develop in secondary lymphoid organs such as lymph nodes (LNs) in a well-coordinated series of interactions between migrating immune cells and resident stromal cells. Although many processes that occur in LNs are well understood from an immunological point of view, our understanding of the fundamental organization and mechanisms that drive these processes is still incomplete. The aim of systems biology approaches is to unravel the complexity of biological systems and describe emergent properties that arise from interactions between individual constituents of the system. The immune system is greater than the sum of its parts, as is the case with any sufficiently complex system. Here, we review recent work and developments of computational LN models with focus on the structure and organization of the stromal cells. We explore various mathematical studies of intranodal T cell motility and migration, their interactions with the LN-resident stromal cells, and computational models of functional chemokine gradient fields and lymph flow dynamics. Lastly, we discuss briefly the importance of hybrid and multi-scale modeling approaches in immunology and the technical challenges involved.

Modeling the multi-scale mechanisms of macromolecular resource allocation

As microbes face changing environments, they dynamically allocate macromolecular resources to produce a particular phenotypic state. Broad 'omics' data sets have revealed several interesting phenomena regarding how the proteome is allocated under differing conditions, but the functional consequences of these states and how they are achieved remain open questions. Various types of multi-scale mathematical models have been used to elucidate the genetic basis for systems-level adaptations. In this review, we outline several different strategies by which microbes accomplish resource allocation and detail how mathematical models have aided in our understanding of these processes. Ultimately, such modeling efforts have helped elucidate the principles of proteome allocation and hold promise for further discovery.

A cell-based mechanical model of coronary artery tunica media

A three-dimensional cell-based mechanical model of coronary artery tunica media is proposed. The model is composed of spherical cells forming a hexagonal close-packed lattice. Tissue anisotropy is taken into account by varying interaction forces with the direction of intercellular connection. Several cell-centre interaction potentials for repulsion and attraction are considered, including the Hertz contact model and its neo-Hookean extension, the Johnson-Kendall-Roberts model of adhesive contact, and a wormlike chain model. The model is validated against data from in vitro uni-axial tension tests performed on dissected strips of tunica media. The wormlike chain potential in combination with the neo-Hookean Hertz contact model produces stress-stretch curves which represent the experimental data very well.

A closed-loop multi-level model of glucose homeostasis

BACKGROUND

The pathophysiologic processes underlying the regulation of glucose homeostasis are considerably complex at both cellular and systemic level. A comprehensive and structured specification for the several layers of abstraction of glucose metabolism is often elusive, an issue currently solvable with the hierarchical description provided by multi-level models. In this study we propose a multi-level closed-loop model of whole-body glucose homeostasis, coupled with the molecular specifications of the insulin signaling cascade in adipocytes, under the experimental conditions of normal glucose regulation and type 2 diabetes.

METHODOLOGY/PRINCIPAL FINDINGS

The ordinary differential equations of the model, describing the dynamics of glucose and key regulatory hormones and their reciprocal interactions among gut, liver, muscle and adipose tissue, were designed for being embedded in a modular, hierarchical structure. The closed-loop model structure allowed self-sustained simulations to represent an ideal in silico subject that adjusts its own metabolism to the fasting and feeding states, depending on the hormonal context and invariant to circadian fluctuations. The cellular level of the model provided a seamless dynamic description of the molecular mechanisms downstream the insulin receptor in the adipocytes by accounting for variations in the surrounding metabolic context.

CONCLUSIONS/SIGNIFICANCE

The combination of a multi-level and closed-loop modeling approach provided a fair dynamic description of the core determinants of glucose homeostasis at both cellular and systemic scales. This model architecture is intrinsically open to incorporate supplementary layers of specifications describing further individual components influencing glucose metabolism.

Making ecological models adequate

Critical evaluation of the adequacy of ecological models is urgently needed to enhance their utility in developing theory and enabling environmental managers and policymakers to make informed decisions. Poorly supported management can have detrimental, costly or irreversible impacts on the environment and society. Here, we examine common issues in ecological modelling and suggest criteria for improving modelling frameworks. An appropriate level of process description is crucial to constructing the best possible model, given the available data and understanding of ecological structures. Model details unsupported by data typically lead to over parameterisation and poor model performance. Conversely, a lack of mechanistic details may limit a model's ability to predict ecological systems' responses to management. Ecological studies that employ models should follow a set of model adequacy assessment protocols that include: asking a series of critical questions regarding state and control variable selection, the determinacy of data, and the sensitivity and validity of analyses. We also need to improve model elaboration, refinement and coarse graining procedures to better understand the relevancy and adequacy of our models and the role they play in advancing theory, improving hind and forecasting, and enabling problem solving and management.

Modeling of the jasmonate signaling pathway in Arabidopsis thaliana with respect to pathophysiology of Alternaria blight in Brassica

The productivity of Oilseed Brassica, one of the economically important crops of India, is seriously affected by the disease, Alternaria blight. The disease is mainly caused by two major necrotrophic fungi, Alternaria brassicae and Alternaria brassicicola which are responsible for significant yield losses. Till date, no resistant source is available against Alternaria blight, hence plant breeding methods can not be used to develop disease resistant varieties. Jasmonate mediated signalling pathway, which is known to play crucial role during defense response against necrotrophs, could be strengthened in Brassica plants to combat the disease. Since scanty information is available in Brassica-Alternaria pathosystems at molecular level therefore, in the present study efforts have been made to model jasmonic acid pathway in Arabidopsis thaliana to simulate the dynamic behaviour of molecular species in the model. Besides, the developed model was also analyzed topologically for investigation of the hubs node. COI1 is identified as one of the promising candidate genes in response to Alternaria and other linked components of plant defense mechanisms against the pathogens. The findings from present study are therefore informative for understanding the molecular basis of pathophysiology and rational management of Alternaria blight for securing food and nutritional security.

Systems approaches in osteoarthritis: Identifying routes to novel diagnostic and therapeutic strategies

Systems orientated research offers the possibility of identifying novel therapeutic targets and relevant diagnostic markers for complex diseases such as osteoarthritis. This review demonstrates that the osteoarthritis research community has been slow to incorporate systems orientated approaches into research studies, although a number of key studies reveal novel insights into the regulatory mechanisms that contribute both to joint tissue homeostasis and its dysfunction. The review introduces both top-down and bottom-up approaches employed in the study of osteoarthritis. A holistic and multiscale approach, where clinical measurements may predict dysregulation and progression of joint degeneration, should be a key objective in future research. The review concludes with suggestions for further research and emerging trends not least of which is the coupled development of diagnostic tests and therapeutics as part of a concerted effort by the osteoarthritis research community to meet clinical needs. © 2017 The Authors. Journal of Orthopaedic Research Published by Wiley Periodicals, Inc. on behalf of Orthopaedic Research Society. J Orthop Res 35:1573-1588, 2017.

Systems modelling ageing: from single senescent cells to simple multi-cellular models

Systems modelling has been successfully used to investigate several key molecular mechanisms of ageing. Modelling frameworks to allow integration of models and methods to enhance confidence in models are now well established. In this article, we discuss these issues and work through the process of building an integrated model for cellular senescence as a single cell and in a simple tissue context.

Facing the challenges of multiscale modelling of bacterial and fungal pathogen-host interactions

Recent and rapidly evolving progress on high-throughput measurement techniques and computational performance has led to the emergence of new disciplines, such as systems medicine and translational systems biology. At the core of these disciplines lies the desire to produce multiscale models: mathematical models that integrate multiple scales of biological organization, ranging from molecular, cellular and tissue models to organ, whole-organism and population scale models. Using such models, hypotheses can systematically be tested. In this review, we present state-of-the-art multiscale modelling of bacterial and fungal infections, considering both the pathogen and host as well as their interaction. Multiscale modelling of the interactions of bacteria, especially Mycobacterium tuberculosis, with the human host is quite advanced. In contrast, models for fungal infections are still in their infancy, in particular regarding infections with the most important human pathogenic fungi, Candida albicans and Aspergillus fumigatus. We reflect on the current availability of computational approaches for multiscale modelling of host-pathogen interactions and point out current challenges. Finally, we provide an outlook for future requirements of multiscale modelling.

Parallelization and High-Performance Computing Enables Automated Statistical Inference of Multi-scale Models

Mechanistic understanding of multi-scale biological processes, such as cell proliferation in a changing biological tissue, is readily facilitated by computational models. While tools exist to construct and simulate multi-scale models, the statistical inference of the unknown model parameters remains an open problem. Here, we present and benchmark a parallel approximate Bayesian computation sequential Monte Carlo (pABC SMC) algorithm, tailored for high-performance computing clusters. pABC SMC is fully automated and returns reliable parameter estimates and confidence intervals. By running the pABC SMC algorithm for ∼10⁶ hr, we parameterize multi-scale models that accurately describe quantitative growth curves and histological data obtained in vivo from individual tumor spheroid growth in media droplets. The models capture the hybrid deterministic-stochastic behaviors of 10⁵-10⁶ of cells growing in a 3D dynamically changing nutrient environment. The pABC SMC algorithm reliably converges to a consistent set of parameters. Our study demonstrates a proof of principle for robust, data-driven modeling of multi-scale biological systems and the feasibility of multi-scale model parameterization through statistical inference.

Agent-Based Modelling in Multicellular Systems Biology

This chapter aims at discussing the content of multi-agent based simulation (MABS) applied to computational biology i.e., to modelling and simulating biological systems by means of computational models, methodologies, and frameworks. In particular, the adoption of agent-based modelling (ABM) in the field of multicellular systems biology is explored, focussing on the challenging scenarios of developmental biology. After motivating why agent-based abstractions are critical in representing multicellular systems behaviour, MABS is discussed as the source of the most natural and appropriate mechanism for analysing the self-organising behaviour of systems of cells. As a case study, an application of MABS to the development of Drosophila Melanogaster is finally presented, which exploits the ALCHEMIST platform for agent-based simulation.

A Novel Method to Verify Multilevel Computational Models of Biological Systems Using Multiscale Spatio-Temporal Meta Model Checking

Insights gained from multilevel computational models of biological systems can be translated into real-life applications only if the model correctness has been verified first. One of the most frequently employed in silico techniques for computational model verification is model checking. Traditional model checking approaches only consider the evolution of numeric values, such as concentrations, over time and are appropriate for computational models of small scale systems (e.g. intracellular networks). However for gaining a systems level understanding of how biological organisms function it is essential to consider more complex large scale biological systems (e.g. organs). Verifying computational models of such systems requires capturing both how numeric values and properties of (emergent) spatial structures (e.g. area of multicellular population) change over time and across multiple levels of organization, which are not considered by existing model checking approaches. To address this limitation we have developed a novel approximate probabilistic multiscale spatio-temporal meta model checking methodology for verifying multilevel computational models relative to specifications describing the desired/expected system behaviour. The methodology is generic and supports computational models encoded using various high-level modelling formalisms because it is defined relative to time series data and not the models used to generate it. In addition, the methodology can be automatically adapted to case study specific types of spatial structures and properties using the spatio-temporal meta model checking concept. To automate the computational model verification process we have implemented the model checking approach in the software tool Mule (http://mule.modelchecking.org). Its applicability is illustrated against four systems biology computational models previously published in the literature encoding the rat cardiovascular system dynamics, the uterine contractions of labour, the Xenopus laevis cell cycle and the acute inflammation of the gut and lung. Our methodology and software will enable computational biologists to efficiently develop reliable multilevel computational models of biological systems.

Multi-class and multi-scale models of complex biological phenomena

Computational modeling has significantly impacted our ability to analyze vast (and exponentially increasing) quantities of experimental data for a variety of applications, such as drug discovery and disease forecasting. Single-scale, single-class models persist as the most common group of models, but biological complexity often demands more sophisticated approaches. This review surveys modeling approaches that are multi-class (incorporating multiple model types) and/or multi-scale (accounting for multiple spatial or temporal scales) and describes how these models, and combinations thereof, should be used within the context of the problem statement. We end by highlighting agent-based models as an intuitive, modular, and flexible framework within which multi-scale and multi-class models can be implemented.

Model calibration and uncertainty analysis in signaling networks

For a long time the biggest challenges in modeling cellular signal transduction networks has been the inference of crucial pathway components and the qualitative description of their interactions. As a result of the emergence of powerful high-throughput experiments, it is now possible to measure data of high temporal and spatial resolution and to analyze signaling dynamics quantitatively. In addition, this increase of high-quality data is the basis for a better understanding of model limitations and their influence on the predictive power of models. We review established approaches in signal transduction network modeling with a focus on ordinary differential equation models as well as related developments in model calibration. As central aspects of the calibration process we discuss possibilities of model adaptation based on data-driven parameter optimization and the concomitant objective of reducing model uncertainties.

Anatomy and Physiology of Multiscale Modeling and Simulation in Systems Medicine

Systems medicine is the application of systems biology concepts, methods, and tools to medical research and practice. It aims to integrate data and knowledge from different disciplines into biomedical models and simulations for the understanding, prevention, cure, and management of complex diseases. Complex diseases arise from the interactions among disease-influencing factors across multiple levels of biological organization from the environment to molecules. To tackle the enormous challenges posed by complex diseases, we need a modeling and simulation framework capable of capturing and integrating information originating from multiple spatiotemporal and organizational scales. Multiscale modeling and simulation in systems medicine is an emerging methodology and discipline that has already demonstrated its potential in becoming this framework. The aim of this chapter is to present some of the main concepts, requirements, and challenges of multiscale modeling and simulation in systems medicine.

Agent-Based Modeling in Systems Pharmacology

Modeling and simulation (M&S) techniques provide a platform for knowledge integration and hypothesis testing to gain insights into biological systems that would not be possible a priori. Agent‐based modeling (ABM) is an M&S technique that focuses on describing individual components rather than homogenous populations. This tutorial introduces ABM to systems pharmacologists, using relevant case studies to highlight how ABM‐specific strengths have yielded success in the area of preclinical mechanistic modeling.

Advancing a distributed multi-scale computing framework for large-scale high-throughput discovery in materials science

We describe the development of a large-scale high-throughput application for discovery in materials science. Our point of departure is a computational framework for distributed multi-scale computation. We augment the original framework with a specialized module whose role is to route evaluation requests needed by the high-throughput application to a collection of available computational resources. We evaluate the feasibility and performance of the resulting high-throughput computational framework by carrying out a high-throughput study of battery solvents. Our results indicate that distributed multi-scale computing, by virtue of its adaptive nature, is particularly well-suited for building high-throughput applications.

Automatising the analysis of stochastic biochemical time-series

BACKGROUND

Mathematical and computational modelling of biochemical systems has seen a lot of effort devoted to the definition and implementation of high-performance mechanistic simulation frameworks. Within these frameworks it is possible to analyse complex models under a variety of configurations, eventually selecting the best setting of, e.g., parameters for a target system.

MOTIVATION

This operational pipeline relies on the ability to interpret the predictions of a model, often represented as simulation time-series. Thus, an efficient data analysis pipeline is crucial to automatise time-series analyses, bearing in mind that errors in this phase might mislead the modeller's conclusions.

RESULTS

For this reason we have developed an intuitive framework-independent Python tool to automate analyses common to a variety of modelling approaches. These include assessment of useful non-trivial statistics for simulation ensembles, e.g., estimation of master equations. Intuitive and domain-independent batch scripts will allow the researcher to automatically prepare reports, thus speeding up the usual model-definition, testing and refinement pipeline.

Agent-based spatiotemporal simulation of biomolecular systems within the open source MASON framework

Agent-based modelling is being used to represent biological systems with increasing frequency and success. This paper presents the implementation of a new tool for biomolecular reaction modelling in the open source Multiagent Simulator of Neighborhoods framework. The rationale behind this new tool is the necessity to describe interactions at the molecular level to be able to grasp emergent and meaningful biological behaviour. We are particularly interested in characterising and quantifying the various effects that facilitate biocatalysis. Enzymes may display high specificity for their substrates and this information is crucial to the engineering and optimisation of bioprocesses. Simulation results demonstrate that molecule distributions, reaction rate parameters, and structural parameters can be adjusted separately in the simulation allowing a comprehensive study of individual effects in the context of realistic cell environments. While higher percentage of collisions with occurrence of reaction increases the affinity of the enzyme to the substrate, a faster reaction (i.e., turnover number) leads to a smaller number of time steps. Slower diffusion rates and molecular crowding (physical hurdles) decrease the collision rate of reactants, hence reducing the reaction rate, as expected. Also, the random distribution of molecules affects the results significantly.