Hierarchical semantic composition of biosimulation models using bond graphs

Computational fluid dynamic modeling of the lymphatic system: a review of existing models and future directions

Historically, research into the lymphatic system has been overlooked due to both a lack of knowledge and limited recognition of its importance. In the last decade however, lymphatic research has gained substantial momentum and has included the development of a variety of computational models to aid understanding of this complex system. This article reviews existing computational fluid dynamic models of the lymphatics covering each structural component including the initial lymphatics, pre-collecting and collecting vessels, and lymph nodes. This is followed by a summary of limitations and gaps in existing computational models and reasons that development in this field has been hindered to date. Over the next decade, efforts to further characterize lymphatic anatomy and physiology are anticipated to provide key data to further inform and validate lymphatic fluid dynamic models. Development of more comprehensive multiscale- and multi-physics computational models has the potential to significantly enhance the understanding of lymphatic function in both health and disease.

Building a search tool for compositely annotated entities using Transformer-based approach: Case study in Biosimulation Model Search Engine (BMSE)

The Transformer-based approaches to solving natural language processing (NLP) tasks such as BERT and GPT are gaining popularity due to their ability to achieve high performance. These approaches benefit from using enormous data sizes to create pre-trained models and the ability to understand the context of words in a sentence. Their use in the information retrieval domain is thought to increase effectiveness and efficiency. This paper demonstrates a BERT-based method (CASBERT) implementation to build a search tool over data annotated compositely using ontologies. The data was a collection of biosimulation models written using the CellML standard in the Physiome Model Repository (PMR). A biosimulation model structurally consists of basic entities of constants and variables that construct higher-level entities such as components, reactions, and the model. Finding these entities specific to their level is beneficial for various purposes regarding variable reuse, experiment setup, and model audit. Initially, we created embeddings representing compositely-annotated entities for constant and variable search (lowest level entity). Then, these low-level entity embeddings were vertically and efficiently combined to create higher-level entity embeddings to search components, models, images, and simulation setups. Our approach was general, so it can be used to create search tools with other data semantically annotated with ontologies - biosimulation models encoded in the SBML format, for example. Our tool is named Biosimulation Model Search Engine (BMSE).

Open problems in mathematical biology

Biology is data-rich, and it is equally rich in concepts and hypotheses. Part of trying to understand biological processes and systems is therefore to confront our ideas and hypotheses with data using statistical methods to determine the extent to which our hypotheses agree with reality. But doing so in a systematic way is becoming increasingly challenging as our hypotheses become more detailed, and our data becomes more complex. Mathematical methods are therefore gaining in importance across the life- and biomedical sciences. Mathematical models allow us to test our understanding, make testable predictions about future behaviour, and gain insights into how we can control the behaviour of biological systems. It has been argued that mathematical methods can be of great benefit to biologists to make sense of data. But mathematics and mathematicians are set to benefit equally from considering the often bewildering complexity inherent to living systems. Here we present a small selection of open problems and challenges in mathematical biology. We have chosen these open problems because they are of both biological and mathematical interest.

Energy-based advection modelling using bond graphs

Advection, the transport of a substance by the flow of a fluid, is a key process in biological systems. The energy-based bond graph approach to modelling chemical transformation within reaction networks is extended to include transport and thus advection. The approach is illustrated using a simple model of advection via circulating flow and by a simple pharmacokinetic model of anaesthetic gas uptake. This extension provides a physically consistent framework for linking advective flows with the fluxes associated with chemical reactions within the context of physiological systems in general and the human physiome in particular.

SBML to bond graphs: From conversion to composition

The Systems Biology Markup Language (SBML) is a popular software-independent XML-based format for describing models of biological phenomena. The BioModels Database is the largest online repository of SBML models. Several tools and platforms are available to support the reuse and composition of SBML models. However, these tools do not explicitly assess whether models are physically plausible or thermodynamically consistent. This often leads to ill-posed models that are physically impossible, impeding the development of realistic complex models in biology. Here, we present a framework that can automatically convert SBML models into bond graphs, which imposes energy conservation laws on these models. The new bond graph models are easily mergeable, resulting in physically plausible coupled models. We illustrate this by automatically converting and coupling a model of pyruvate distribution to a model of the pentose phosphate pathway.

Network thermodynamics of biological systems: A bond graph approach

Edmund Crampin (1973-2021) was at the forefront of Systems Biology research and his work will influence the field for years to come. This paper brings together and summarises the seminal work of his group in applying energy-based bond graph methods to biological systems. In particular, this paper: (a) motivates the need to consider energy in modelling biology; (b) introduces bond graphs as a methodology for achieving this; (c) describes extensions to modelling electrochemical transduction; (d) outlines how bond graph models can be constructed in a modular manner and (e) describes stoichiometric approaches to deriving fundamental properties of reaction networks. These concepts are illustrated using a new bond graph model of photosynthesis in chloroplasts.

A semantics, energy-based approach to automate biomodel composition

Hierarchical modelling is essential to achieving complex, large-scale models. However, not all modelling schemes support hierarchical composition, and correctly mapping points of connection between models requires comprehensive knowledge of each model's components and assumptions. To address these challenges in integrating biosimulation models, we propose an approach to automatically and confidently compose biosimulation models. The approach uses bond graphs to combine aspects of physical and thermodynamics-based modelling with biological semantics. We improved on existing approaches by using semantic annotations to automate the recognition of common components. The approach is illustrated by coupling a model of the Ras-MAPK cascade to a model of the upstream activation of EGFR. Through this methodology, we aim to assist researchers and modellers in readily having access to more comprehensive biological systems models.

The Cell Physiome: What Do We Need in a Computational Physiology Framework for Predicting Single-Cell Biology?

Modern biology and biomedicine are undergoing a big data explosion, needing advanced computational algorithms to extract mechanistic insights on the physiological state of living cells. We present the motivation for the Cell Physiome project: a framework and approach for creating, sharing, and using biophysics-based computational models of single-cell physiology. Using examples in calcium signaling, bioenergetics, and endosomal trafficking, we highlight the need for spatially detailed, biophysics-based computational models to uncover new mechanisms underlying cell biology. We review progress and challenges to date toward creating cell physiome models. We then introduce bond graphs as an efficient way to create cell physiome models that integrate chemical, mechanical, electromagnetic, and thermal processes while maintaining mass and energy balance. Bond graphs enhance modularization and reusability of computational models of cells at scale. We conclude with a look forward at steps that will help fully realize this exciting new field of mechanistic biomedical data science. Expected final online publication date for the Annual Review of Biomedical Data Science, Volume 5 is August 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Simplifying the Process of Going From Cells to Tissues Using Statistical Mechanics

The value of digital twins for prototyping controllers or interventions in a sandbox environment are well-established in engineering and physics. However, this is challenging for biophysics trying to seamlessly compose models of multiple spatial and temporal scale behavior into the digital twin. Two challenges stand out as constraining progress: (i) ensuring physical consistency of conservation laws across composite models and (ii) drawing useful and timely clinical and scientific information from conceptually and computationally complex models. Challenge (i) can be robustly addressed with bondgraphs. However, challenge (ii) is exacerbated using this approach. The complexity question can be looked at from multiple angles. First from the perspective of discretizations that reflect underlying biophysics (functional tissue units) and secondly by exploring maximum entropy as the principle guiding multicellular biophysics. Statistical mechanics, long applied to understanding emergent phenomena from atomic physics, coupled with the observation that cellular architecture in tissue is orchestrated by biophysical constraints on metabolism and communication, shows conceptual promise. This architecture along with cell specific properties can be used to define tissue specific network motifs associated with energetic contributions. Complexity can be addressed based on energy considerations and finding mean measures of dependent variables. A probability distribution of the tissue's network motif can be approximated with exponential random graph models. A prototype problem shows how these approaches could be implemented in practice and the type of information that could be extracted.

Modular assembly of dynamic models in systems biology

It is widely acknowledged that the construction of large-scale dynamic models in systems biology requires complex modelling problems to be broken up into more manageable pieces. To this end, both modelling and software frameworks are required to enable modular modelling. While there has been consistent progress in the development of software tools to enhance model reusability, there has been a relative lack of consideration for how underlying biophysical principles can be applied to this space. Bond graphs combine the aspects of both modularity and physics-based modelling. In this paper, we argue that bond graphs are compatible with recent developments in modularity and abstraction in systems biology, and are thus a desirable framework for constructing large-scale models. We use two examples to illustrate the utility of bond graphs in this context: a model of a mitogen-activated protein kinase (MAPK) cascade to illustrate the reusability of modules and a model of glycolysis to illustrate the ability to modify the model granularity.