Hybrid modelling of biological systems: current progress and future prospects

Systematic simulation of tumor cell invasion and migration in response to time-varying rotating magnetic field

Cancer invasion and migration play a pivotal role in tumor malignancy, which is a major cause of most cancer deaths. Rotating magnetic field (RMF), one of the typical dynamic magnetic fields, can exert substantial mechanical influence on cells. However, studying the effects of RMF on cell is challenging due to its complex parameters, such as variation of magnetic field intensity and direction. Here, we developed a systematic simulation method to explore the influence of RMF on tumor invasion and migration, including a finite element method (FEM) model and a cell-based hybrid numerical model. Coupling with the data of magnetic field from FEM, the cell-based hybrid numerical model was established to simulate the tumor cell invasion and migration. This model employed partial differential equations (PDEs) and finite difference method to depict cellular activities and solve these equations in a discrete system. PDEs were used to depict cell activities, and finite difference method was used to solve the equations in discrete system. As a result, this study provides valuable insights into the potential applications of RMF in tumor treatment, and a series of in vitro experiments were performed to verify the simulation results, demonstrating the model's reliability and its capacity to predict experimental outcomes and identify pertinent factors. Furthermore, these findings shed new light on the mechanical and chemical interplay between cells and the ECM, offering new insights and providing a novel foundation for both experimental and theoretical advancements in tumor treatment by using RMF.

Mathematical Modeling and Inference of Epidermal Growth Factor-Induced Mitogen-Activated Protein Kinase Cell Signaling Pathways

The mitogen-activated protein kinase (MAPK) pathway is an important intracellular signaling cascade that plays a key role in various cellular processes. Understanding the regulatory mechanisms of this pathway is essential for developing effective interventions and targeted therapies for related diseases. Recent advances in single-cell proteomic technologies have provided unprecedented opportunities to investigate the heterogeneity and noise within complex, multi-signaling networks across diverse cells and cell types. Mathematical modeling has become a powerful interdisciplinary tool that bridges mathematics and experimental biology, providing valuable insights into these intricate cellular processes. In addition, statistical methods have been developed to infer pathway topologies and estimate unknown parameters within dynamic models. This review presents a comprehensive analysis of how mathematical modeling of the MAPK pathway deepens our understanding of its regulatory mechanisms, enhances the prediction of system behavior, and informs experimental research, with a particular focus on recent advances in modeling and inference using single-cell proteomic data.

Network Modeling and Control of Dynamic Disease Pathways, Review and Perspectives

Dynamic disease pathways are a combination of complex dynamical processes among bio-molecules in a cell that leads to diseases. Network modeling of disease pathways considers disease-related bio-molecules (e.g. DNA, RNA, transcription factors, enzymes, proteins, and metabolites) and their interaction (e.g. DNA methylation, histone modification, alternative splicing, and protein modification) to study disease progression and predict therapeutic responses. These bio-molecules and their interactions are the basic elements in the study of the misregulation in the disease-related gene expression that lead to abnormal cellular responses. Gene regulatory networks, cell signaling networks, and metabolic networks are the three major types of intracellular networks for the study of the cellular responses elicited from extracellular signals. The disease-related cellular responses can be prevented or regulated by designing control strategies to manipulate these extracellular or other intracellular signals. The paper reviews the regulatory mechanisms, the dynamic models, and the control strategies for each intracellular network. The applications, limitations and the prospective for modeling and control are also discussed.

Bridging systems biology and tissue engineering: Unleashing the full potential of complex 3D in vitro tissue models of disease

Rapid advances in tissue engineering have resulted in more complex and physiologically relevant 3D in vitro tissue models with applications in fundamental biology and therapeutic development. However, the complexity provided by these models is often not leveraged fully due to the reductionist methods used to analyze them. Computational and mathematical models developed in the field of systems biology can address this issue. Yet, traditional systems biology has been mostly applied to simpler in vitro models with little physiological relevance and limited cellular complexity. Therefore, integrating these two inherently interdisciplinary fields can result in new insights and move both disciplines forward. In this review, we provide a systematic overview of how systems biology has been integrated with 3D in vitro tissue models and discuss key application areas where the synergies between both fields have led to important advances with potential translational impact. We then outline key directions for future research and discuss a framework for further integration between fields.

Design patterns for the construction of computational biological models

Computational biological models have proven to be an invaluable tool for understanding and predicting the behaviour of many biological systems. While it may not be too challenging for experienced researchers to construct such models from scratch, it is not a straightforward task for early stage researchers. Design patterns are well-known techniques widely applied in software engineering as they provide a set of typical solutions to common problems in software design. In this paper, we collect and discuss common patterns that are usually used during the construction and execution of computational biological models. We adopt Petri nets as a modelling language to provide a visual illustration of each pattern; however, the ideas presented in this paper can also be implemented using other modelling formalisms. We provide two case studies for illustration purposes and show how these models can be built up from the presented smaller modules. We hope that the ideas discussed in this paper will help many researchers in building their own future models.

MCell4 with BioNetGen: A Monte Carlo simulator of rule-based reaction-diffusion systems with Python interface

Biochemical signaling pathways in living cells are often highly organized into spatially segregated volumes, membranes, scaffolds, subcellular compartments, and organelles comprising small numbers of interacting molecules. At this level of granularity stochastic behavior dominates, well-mixed continuum approximations based on concentrations break down and a particle-based approach is more accurate and more efficient. We describe and validate a new version of the open-source MCell simulation program (MCell4), which supports generalized 3D Monte Carlo modeling of diffusion and chemical reaction of discrete molecules and macromolecular complexes in solution, on surfaces representing membranes, and combinations thereof. The main improvements in MCell4 compared to the previous versions, MCell3 and MCell3-R, include a Python interface and native BioNetGen reaction language (BNGL) support. MCell4's Python interface opens up completely new possibilities for interfacing with external simulators to allow creation of sophisticated event-driven multiscale/multiphysics simulations. The native BNGL support, implemented through a new open-source library libBNG (also introduced in this paper), provides the capability to run a given BNGL model spatially resolved in MCell4 and, with appropriate simplifying assumptions, also in the BioNetGen simulation environment, greatly accelerating and simplifying model validation and comparison.

Community-scale models of microbiomes: Articulating metabolic modelling and metagenome sequencing

Building models is essential for understanding the functions and dynamics of microbial communities. Metabolic models built on genome-scale metabolic network reconstructions (GENREs) are especially relevant as a means to decipher the complex interactions occurring among species. Model reconstruction increasingly relies on metagenomics, which permits direct characterisation of naturally occurring communities that may contain organisms that cannot be isolated or cultured. In this review, we provide an overview of the field of metabolic modelling and its increasing reliance on and synergy with metagenomics and bioinformatics. We survey the means of assigning functions and reconstructing metabolic networks from (meta-)genomes, and present the variety and mathematical fundamentals of metabolic models that foster the understanding of microbial dynamics. We emphasise the characterisation of interactions and the scaling of model construction to large communities, two important bottlenecks in the applicability of these models. We give an overview of the current state of the art in metagenome sequencing and bioinformatics analysis, focusing on the reconstruction of genomes in microbial communities. Metagenomics benefits tremendously from third-generation sequencing, and we discuss the opportunities of long-read sequencing, strain-level characterisation and eukaryotic metagenomics. We aim at providing algorithmic and mathematical support, together with tool and application resources, that permit bridging the gap between metagenomics and metabolic modelling.

Artificial intelligence and bioinformatics: a journey from traditional techniques to smart approaches

The incorporation of AI models into bioinformatics has brought about a revolutionary era in the analysis and interpretation of biological data. This mini-review offers a succinct overview of the indispensable role AI plays in the convergence of computational techniques and biological research. The search strategy followed PRISMA guidelines, encompassing databases such as PubMed, Embase, and Google Scholar to include studies published between 2018 and 2024, utilizing specific keywords. We explored the diverse applications of AI methodologies, including machine learning (ML), deep learning (DL), and natural language processing (NLP), across various domains of bioinformatics. These domains encompass genome sequencing, protein structure prediction, drug discovery, systems biology, personalized medicine, imaging, signal processing, and text mining. AI algorithms have exhibited remarkable efficacy in tackling intricate biological challenges, spanning from genome sequencing to protein structure prediction, and from drug discovery to personalized medicine. In conclusion, this study scrutinizes the evolving landscape of AI-driven tools and algorithms, emphasizing their pivotal role in expediting research, facilitating data interpretation, and catalyzing innovations in biomedical sciences.

Early acquired resistance to EGFR-TKIs in lung adenocarcinomas before radiographic advanced identified by CT radiomic delta model based on two central studies

Early acquired resistance (EAR) to epidermal growth factor receptor tyrosine kinase inhibitors (EGFR-TKIs) in lung adenocarcinomas before radiographic advance cannot be perceived by the naked eye. This study aimed to discover and validate a CT radiomic model to precisely identify the EAR. Training cohort (n = 67) and internal test cohort (n = 29) were from the First Affiliated Hospital of Fujian Medical University, and external test cohort (n = 29) was from the Second Affiliated Hospital of Xiamen Medical College. Follow-up CT images at three different times of each patient were collected: (1) baseline images before EGFR-TKIs therapy; (2) first follow-up images after EGFR-TKIs therapy (FFT); (3) EAR images, which were the last follow-up images before radiographic advance. The features extracted from FFT and EAR were used to construct the classic radiomic model. The delta features which were calculated by subtracting the baseline from either FFT or EAR were used to construct the delta radiomic model. The classic radiomic model achieved AUC 0.682 and 0.641 in training and internal test cohorts, respectively. The delta radiomic model achieved AUC 0.730 and 0.704 in training and internal test cohorts, respectively. Over the external test cohort, the delta radiomic model achieved AUC 0.661. The decision curve analysis showed that when threshold of the probability of the EAR to the EGFR-TKIs was between 0.3 and 0.82, the proposed model was more benefit than treating all patients. Based on two central studies, the delta radiomic model derived from the follow-up non-enhanced CT images can help clinicians to identify the EAR to EGFR-TKIs in lung adenocarcinomas before radiographic advance and optimize clinical outcomes.

Scientific novelty beyond the experiment

Practical experiments drive important scientific discoveries in biology, but theory-based research studies also contribute novel-sometimes paradigm-changing-findings. Here, we appraise the roles of theory-based approaches focusing on the experiment-dominated wet-biology research areas of microbial growth and survival, cell physiology, host-pathogen interactions, and competitive or symbiotic interactions. Additional examples relate to analyses of genome-sequence data, climate change and planetary health, habitability, and astrobiology. We assess the importance of thought at each step of the research process; the roles of natural philosophy, and inconsistencies in logic and language, as drivers of scientific progress; the value of thought experiments; the use and limitations of artificial intelligence technologies, including their potential for interdisciplinary and transdisciplinary research; and other instances when theory is the most-direct and most-scientifically robust route to scientific novelty including the development of techniques for practical experimentation or fieldwork. We highlight the intrinsic need for human engagement in scientific innovation, an issue pertinent to the ongoing controversy over papers authored using/authored by artificial intelligence (such as the large language model/chatbot ChatGPT). Other issues discussed are the way in which aspects of language can bias thinking towards the spatial rather than the temporal (and how this biased thinking can lead to skewed scientific terminology); receptivity to research that is non-mainstream; and the importance of theory-based science in education and epistemology. Whereas we briefly highlight classic works (those by Oakes Ames, Francis H.C. Crick and James D. Watson, Charles R. Darwin, Albert Einstein, James E. Lovelock, Lynn Margulis, Gilbert Ryle, Erwin R.J.A. Schrödinger, Alan M. Turing, and others), the focus is on microbiology studies that are more-recent, discussing these in the context of the scientific process and the types of scientific novelty that they represent. These include several studies carried out during the 2020 to 2022 lockdowns of the COVID-19 pandemic when access to research laboratories was disallowed (or limited). We interviewed the authors of some of the featured microbiology-related papers and-although we ourselves are involved in laboratory experiments and practical fieldwork-also drew from our own research experiences showing that such studies can not only produce new scientific findings but can also transcend barriers between disciplines, act counter to scientific reductionism, integrate biological data across different timescales and levels of complexity, and circumvent constraints imposed by practical techniques. In relation to urgent research needs, we believe that climate change and other global challenges may require approaches beyond the experiment.

Targeted Quantification of Protein Phosphorylation and Its Contributions towards Mathematical Modeling of Signaling Pathways

Post-translational modifications (PTMs) are key regulatory mechanisms that can control protein function. Of these, phosphorylation is the most common and widely studied. Because of its importance in regulating cell signaling, precise and accurate measurements of protein phosphorylation across wide dynamic ranges are crucial to understanding how signaling pathways function. Although immunological assays are commonly used to detect phosphoproteins, their lack of sensitivity, specificity, and selectivity often make them unreliable for quantitative measurements of complex biological samples. Recent advances in Mass Spectrometry (MS)-based targeted proteomics have made it a more useful approach than immunoassays for studying the dynamics of protein phosphorylation. Selected reaction monitoring (SRM)-also known as multiple reaction monitoring (MRM)-and parallel reaction monitoring (PRM) can quantify relative and absolute abundances of protein phosphorylation in multiplexed fashions targeting specific pathways. In addition, the refinement of these tools by enrichment and fractionation strategies has improved measurement of phosphorylation of low-abundance proteins. The quantitative data generated are particularly useful for building and parameterizing mathematical models of complex phospho-signaling pathways. Potentially, these models can provide a framework for linking analytical measurements of clinical samples to better diagnosis and treatment of disease.

Arduino Soft Sensor for Monitoring Schizochytrium sp. Fermentation, a Proof of Concept for the Industrial Application of Genome-Scale Metabolic Models in the Context of Pharma 4.0

Schizochytrium sp. is a microorganism cultured for producing docosahexaenoic acid (DHA). Genome-scale metabolic modeling (GEM) is a promising technique for describing gen-protein-reactions in cells, but with still limited industrial application due to its complexity and high computation requirements. In this work, we simplified GEM results regarding the relationship between the specific oxygen uptake rate (−rO2), the specific growth rate (µ), and the rate of lipid synthesis (rL) using an evolutionary algorithm for developing a model that can be used by a soft sensor for fermentation monitoring. The soft sensor estimated the concentration of active biomass (X), glutamate (N), lipids (L), and DHA in a Schizochytrium sp. fermentation using the dissolved oxygen tension (DO) and the oxygen mass transfer coefficient (kLa) as online input variables. The soft sensor model described the biomass concentration response of four reported experiments characterized by different kLa values. The average range normalized root-mean-square error for X, N, L, and DHA were equal to 1.1, 1.3, 1.1, and 3.2%, respectively, suggesting an acceptable generalization capacity. The feasibility of implementing the soft sensor over a low-cost electronic board was successfully tested using an Arduino UNO, showing a novel path for applying GEM-based soft sensors in the context of Pharma 4.0.