The simulation experiment description markup language (SED-ML): language specification for level 1 version 5.

Modern biological research is increasingly informed by computational simulation experiments, which necessitate the development of methods for annotating, archiving, sharing, and reproducing the conducted experiments. These simulations increasingly require extensive collaboration among modelers, experimentalists, and engineers. The Minimum Information About a Simulation Experiment (MIASE) guidelines outline the information needed to share simulation experiments. SED-ML is a computer-readable format for the information outlined by MIASE, created as a community project and supported by many investigators and software tools. Level 1 Version 5 of SED-ML expands the ability of modelers to define simulations in SED-ML using the Kinetic Simulation Algorithm Onotoloy (KiSAO). While it was possible in Version 4 to define a simulation entirely using KiSAO, Version 5 now allows users to define tasks, model changes, ranges, and outputs using the ontology as well. SED-ML is supported by a growing ecosystem of investigators, model languages, and software tools, including various languages for constraint-based, kinetic, qualitative, rule-based, and spatial models, and many simulation tools, visual editors, model repositories, and validators. Additional information about SED-ML is available at https://sed-ml.org/.

Toward mechanistic medical digital twins: some use cases in immunology.

A fundamental challenge for personalized medicine is to capture enough of the complexity of an individual patient to determine an optimal way to keep them healthy or restore their health. This will require personalized computational models of sufficient resolution and with enough mechanistic information to provide actionable information to the clinician. Such personalized models are increasingly referred to as medical digital twins. Digital twin technology for health applications is still in its infancy, and extensive research and development is required. This article focuses on several projects in different stages of development that can lead to specific-and practical-medical digital twins or digital twin modeling platforms. It emerged from a two-day forum on problems related to medical digital twins, particularly those involving an immune system component. Open access video recordings of the forum discussions are available.

Forum on immune digital twins: a meeting report.

Medical digital twins are computational models of human biology relevant to a given medical condition, which are tailored to an individual patient, thereby predicting the course of disease and individualized treatments, an important goal of personalized medicine. The immune system, which has a central role in many diseases, is highly heterogeneous between individuals, and thus poses a major challenge for this technology. In February 2023, an international group of experts convened for two days to discuss these challenges related to immune digital twins. The group consisted of clinicians, immunologists, biologists, and mathematical modelers, representative of the interdisciplinary nature of medical digital twin development. A video recording of the entire event is available. This paper presents a synopsis of the discussions, brief descriptions of ongoing digital twin projects at different stages of progress. It also proposes a 5-year action plan for further developing this technology. The main recommendations are to identify and pursue a small number of promising use cases, to develop stimulation-specific assays of immune function in a clinical setting, and to develop a database of existing computational immune models, as well as advanced modeling technology and infrastructure.

Augusta: From RNA-Seq to gene regulatory networks and Boolean models.

Computational models of gene regulations help to understand regulatory mechanisms and are extensively used in a wide range of areas, e.g., biotechnology or medicine, with significant benefits. Unfortunately, there are only a few computational gene regulatory models of whole genomes allowing static and dynamic analysis due to the lack of sophisticated tools for their reconstruction. Here, we describe Augusta, an open-source Python package for Gene Regulatory Network (GRN) and Boolean Network (BN) inference from the high-throughput gene expression data. Augusta can reconstruct genome-wide models suitable for static and dynamic analyses. Augusta uses a unique approach where the first estimation of a GRN inferred from expression data is further refined by predicting transcription factor binding motifs in promoters of regulated genes and by incorporating verified interactions obtained from databases. Moreover, a refined GRN is transformed into a draft BN by searching in the curated model database and setting logical rules to incoming edges of target genes, which can be further manually edited as the model is provided in the SBML file format. The approach is applicable even if information about the organism under study is not available in the databases, which is typically the case for non-model organisms including most microbes. Augusta can be operated from the command line and, thus, is easy to use for automated prediction of models for various genomes. The Augusta package is freely available at github.com/JanaMus/Augusta. Documentation and tutorials are available at augusta.readthedocs.io.

COMO: a pipeline for multi-omics data integration in metabolic modeling and drug discovery.

Identifying potential drug targets using metabolic modeling requires integrating multiple modeling methods and heterogeneous biological datasets, which can be challenging without efficient tools. We developed Constraint-based Optimization of Metabolic Objectives (COMO), a user-friendly pipeline that integrates multi-omics data processing, context-specific metabolic model development, simulations, drug databases and disease data to aid drug discovery. COMO can be installed as a Docker Image or with Conda and includes intuitive instructions within a Jupyter Lab environment. It provides a comprehensive solution for the integration of bulk and single-cell RNA-seq, microarrays and proteomics outputs to develop context-specific metabolic models. Using public databases, open-source solutions for model construction and a streamlined approach for predicting repurposable drugs, COMO enables researchers to investigate low-cost alternatives and novel disease treatments. As a case study, we used the pipeline to construct metabolic models of B cells, which simulate and analyze them to predict metabolic drug targets for rheumatoid arthritis and systemic lupus erythematosus, respectively. COMO can be used to construct models for any cell or tissue type and identify drugs for any human disease where metabolic inhibition is relevant. The pipeline has the potential to improve the health of the global community cost-effectively by providing high-confidence targets to pursue in preclinical and clinical studies. The source code of the COMO pipeline is available at https://github.com/HelikarLab/COMO. The Docker image can be pulled at https://github.com/HelikarLab/COMO/pkgs/container/como.

Modeling and Remodeling the Cell: How Digital Twins and HCMV Can Elucidate the Complex Interactions of Viral Latency, Epigenetic Regulation, and Immune Responses.

Purpose of Review: Human cytomegalovirus (HCMV), while asymptomatic in most, causes significant complications during fetal development, following transplant or in immunosuppressed individuals. The host-virus interactions regulating viral latency and reactivation and viral control of the cellular environment (immune regulation, differentiation, epigenetics) are highly complex. Understanding these processes is essential to controlling infection and can be leveraged as a novel approach for understanding basic cell biology. Recent Findings: Immune digital twins (IDTs) are digital simulations integrating knowledge of human immunology, physiology, and patient-specific clinical data to predict individualized immune responses and targeted treatments. Recent studies used IDTs to elucidate mechanisms of T cells, dendritic cells, and epigenetic control-all key to HCMV biology. Summary: Here, we discuss how leveraging the unique biology of HCMV and IDTs will clarify immune response dynamics, host-virus interactions, and viral latency and reactivation and serve as a powerful IDT-validation platform for individualized and holistic health management.

Factors Influencing Instructors' Adoption and Continued Use of Computing Science Technologies: A Case Study in the Context of Cell Collective.

Acquiring computational modeling and simulation skills has become ever more critical for students in life sciences courses at the secondary and tertiary levels. Many modeling and simulation tools have been created to help instructors nurture those skills in their classrooms. Understanding the factors that may motivate instructors to use such tools is crucial to improve students' learning, especially for having authentic modeling and simulation learning experiences. This study designed and tested a decomposed technology acceptance model in which the perceived usefulness and perceived ease of use constructs are split between the teaching and learning sides of the technology to examine their relative weight in a single model. Using data from instructors using the Cell Collective modeling and simulation software, this study found that the relationship between perceived usefulness-teaching and attitude toward behavior was insignificant. Similarly, all relationships between perceived ease of use-teaching and the other variables (i.e., perceived usefulness-teaching and attitude toward behavior) became insignificant. In contrast, we found the relationships between perceived ease of use-learning and the other variables (i.e., perceived usefulness-teaching, perceived usefulness-learning, and attitude toward behavior) significant. These results suggest that priority should be given to the development of features improving learning over features facilitating teaching.

A multiscale mechanistic model of human dendritic cells for in-silico investigation of immune responses and novel therapeutics discovery.

Dendritic cells (DCs) are professional antigen-presenting cells (APCs) with the unique ability to mediate inflammatory responses of the immune system. Given the critical role of DCs in shaping immunity, they present an attractive avenue as a therapeutic target to program the immune system and reverse immune disease disorders. To ensure appropriate immune response, DCs utilize intricate and complex molecular and cellular interactions that converge into a seamless phenotype. Computational models open novel frontiers in research by integrating large-scale interaction to interrogate the influence of complex biological behavior across scales. The ability to model large biological networks will likely pave the way to understanding any complex system in more approachable ways. We developed a logical and predictive model of DC function that integrates the heterogeneity of DCs population, APC function, and cell-cell interaction, spanning molecular to population levels. Our logical model consists of 281 components that connect environmental stimuli with various layers of the cell compartments, including the plasma membrane, cytoplasm, and nucleus to represent the dynamic processes within and outside the DC, such as signaling pathways and cell-cell interactions. We also provided three sample use cases to apply the model in the context of studying cell dynamics and disease environments. First, we characterized the DC response to Sars-CoV-2 and influenza co-infection by in-silico experiments and analyzed the activity level of 107 molecules that play a role in this co-infection. The second example presents simulations to predict the crosstalk between DCs and T cells in a cancer microenvironment. Finally, for the third example, we used the Kyoto Encyclopedia of Genes and Genomes enrichment analysis against the model's components to identify 45 diseases and 24 molecular pathways that the DC model can address. This study presents a resource to decode the complex dynamics underlying DC-derived APC communication and provides a platform for researchers to perform in-silico experiments on human DC for vaccine design, drug discovery, and immunotherapies.

A Review of Radiation-Induced Alterations of Multi-Omic Profiles, Radiation Injury Biomarkers, and Countermeasures.

Increasing utilization of nuclear power enhances the risks associated with industrial accidents, occupational hazards, and the threat of nuclear terrorism. Exposure to ionizing radiation interferes with genomic stability and gene expression resulting in the disruption of normal metabolic processes in cells and organs by inducing complex biological responses. Exposure to high-dose radiation causes acute radiation syndrome, which leads to hematopoietic, gastrointestinal, cerebrovascular, and many other organ-specific injuries. Altered genomic variations, gene expression, metabolite concentrations, and microbiota profiles in blood plasma or tissue samples reflect the whole-body radiation injuries. Hence, multi-omic profiles obtained from high-resolution omics platforms offer a holistic approach for identifying reliable biomarkers to predict the radiation injury of organs and tissues resulting from radiation exposures. In this review, we performed a literature search to systematically catalog the radiation-induced alterations from multi-omic studies and radiation countermeasures. We covered radiation-induced changes in the genomic, transcriptomic, proteomic, metabolomic, lipidomic, and microbiome profiles. Furthermore, we have covered promising multi-omic biomarkers, FDA-approved countermeasure drugs, and other radiation countermeasures that include radioprotectors and radiomitigators. This review presents an overview of radiation-induced alterations of multi-omics profiles and biomarkers, and associated radiation countermeasures.

Drug-target identification in COVID-19 disease mechanisms using computational systems biology approaches.

Introduction: The COVID-19 Disease Map project is a large-scale community effort uniting 277 scientists from 130 Institutions around the globe. We use high-quality, mechanistic content describing SARS-CoV-2-host interactions and develop interoperable bioinformatic pipelines for novel target identification and drug repurposing. Methods: Extensive community work allowed an impressive step forward in building interfaces between Systems Biology tools and platforms. Our framework can link biomolecules from omics data analysis and computational modelling to dysregulated pathways in a cell-, tissue- or patient-specific manner. Drug repurposing using text mining and AI-assisted analysis identified potential drugs, chemicals and microRNAs that could target the identified key factors. Results: Results revealed drugs already tested for anti-COVID-19 efficacy, providing a mechanistic context for their mode of action, and drugs already in clinical trials for treating other diseases, never tested against COVID-19. Discussion: The key advance is that the proposed framework is versatile and expandable, offering a significant upgrade in the arsenal for virus-host interactions and other complex pathologies.

Arginine Catabolism and Polyamine Biosynthesis Pathway Disparities Within Francisella tularensis Subpopulations.

Francisella tularensis is a highly infectious zoonotic pathogen with as few as 10 organisms causing tularemia, a disease that is fatal if untreated. Although F. tularensis subspecies tularensis (type A) and subspecies holarctica (type B) share over 99.5% average nucleotide identity, notable differences exist in genomic organization and pathogenicity. The type A clade has been further divided into subtypes A.I and A.II, with A.I strains being recognized as some of the most virulent bacterial pathogens known. In this study, we report on major disparities that exist between the F. tularensis subpopulations in arginine catabolism and subsequent polyamine biosynthesis. The genes involved in these pathways include the speHEA and aguAB operons, along with metK. In the hypervirulent F. tularensis A.I clade, such as the A.I prototype strain SCHU S4, these genes were found to be intact and highly transcribed. In contrast, both subtype A.II and type B strains have a truncated speA gene, while the type B clade also has a disrupted aguA and truncated aguB. Ablation of the chromosomal speE gene that encodes a spermidine synthase reduced subtype A.I SCHU S4 growth rate, whereas the growth rate of type B LVS was enhanced. These results demonstrate that spermine synthase SpeE promotes faster replication in the F. tularensis A.I clade, whereas type B strains do not rely on this enzyme for in vitro fitness. Our ongoing studies on amino acid and polyamine flux within hypervirulent A.I strains should provide a better understanding of the factors that contribute to F. tularensis pathogenicity.

Addressing barriers in comprehensiveness, accessibility, reusability, interoperability and reproducibility of computational models in systems biology.

Computational models are often employed in systems biology to study the dynamic behaviours of complex systems. With the rise in the number of computational models, finding ways to improve the reusability of these models and their ability to reproduce virtual experiments becomes critical. Correct and effective model annotation in community-supported and standardised formats is necessary for this improvement. Here, we present recent efforts toward a common framework for annotated, accessible, reproducible and interoperable computational models in biology, and discuss key challenges of the field.

Two heads are better than one: current landscape of integrating QSP and machine learning : An ISoP QSP SIG white paper by the working group on the integration of quantitative systems pharmacology and machine learning.

Quantitative systems pharmacology (QSP) modeling is applied to address essential questions in drug development, such as the mechanism of action of a therapeutic agent and the progression of disease. Meanwhile, machine learning (ML) approaches also contribute to answering these questions via the analysis of multi-layer 'omics' data such as gene expression, proteomics, metabolomics, and high-throughput imaging. Furthermore, ML approaches can also be applied to aspects of QSP modeling. Both approaches are powerful tools and there is considerable interest in integrating QSP modeling and ML. So far, a few successful implementations have been carried out from which we have learned about how each approach can overcome unique limitations of the other. The QSP + ML working group of the International Society of Pharmacometrics QSP Special Interest Group was convened in September, 2019 to identify and begin realizing new opportunities in QSP and ML integration. The working group, which comprises 21 members representing 18 academic and industry organizations, has identified four categories of current research activity which will be described herein together with case studies of applications to drug development decision making. The working group also concluded that the integration of QSP and ML is still in its early stages of moving from evaluating available technical tools to building case studies. This paper reports on this fast-moving field and serves as a foundation for future codification of best practices.

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.

Essential role of systemic iron mobilization and redistribution for adaptive thermogenesis through HIF2-α/hepcidin axis.

Iron is an essential biometal, but is toxic if it exists in excess. Therefore, iron content is tightly regulated at cellular and systemic levels to meet metabolic demands but to avoid toxicity. We have recently reported that adaptive thermogenesis, a critical metabolic pathway to maintain whole-body energy homeostasis, is an iron-demanding process for rapid biogenesis of mitochondria. However, little information is available on iron mobilization from storage sites to thermogenic fat. This study aimed to determine the iron-regulatory network that underlies beige adipogenesis. We hypothesized that thermogenic stimulus initiates the signaling interplay between adipocyte iron demands and systemic iron liberation, resulting in iron redistribution into beige fat. To test this hypothesis, we induced reversible activation of beige adipogenesis in C57BL/6 mice by administering a ß3-adrenoreceptor agonist CL 316,243 (CL). Our results revealed that CL stimulation induced the iron-regulatory protein-mediated iron import into adipocytes, suppressed hepcidin transcription, and mobilized iron from the spleen. Mechanistically, CL stimulation induced an acute activation of hypoxia-inducible factor 2-α (HIF2-α), erythropoietin production, and splenic erythroid maturation, leading to hepcidin suppression. Disruption of systemic iron homeostasis by pharmacological HIF2-α inhibitor PT2385 or exogenous administration of hepcidin-25 significantly impaired beige fat development. Our findings suggest that securing iron availability via coordinated interplay between renal hypoxia and hepcidin down-regulation is a fundamental mechanism to activate adaptive thermogenesis. It also provides an insight into the effects of adaptive thermogenesis on systemic iron mobilization and redistribution.

A multi-approach and multi-scale platform to model CD4+ T cells responding to infections.

Immune responses rely on a complex adaptive system in which the body and infections interact at multiple scales and in different compartments. We developed a modular model of CD4+ T cells, which uses four modeling approaches to integrate processes at three spatial scales in different tissues. In each cell, signal transduction and gene regulation are described by a logical model, metabolism by constraint-based models. Cell population dynamics are described by an agent-based model and systemic cytokine concentrations by ordinary differential equations. A Monte Carlo simulation algorithm allows information to flow efficiently between the four modules by separating the time scales. Such modularity improves computational performance and versatility and facilitates data integration. We validated our technology by reproducing known experimental results, including differentiation patterns of CD4+ T cells triggered by different combinations of cytokines, metabolic regulation by IL2 in these cells, and their response to influenza infection. In doing so, we added multi-scale insights to single-scale studies and demonstrated its predictive power by discovering switch-like and oscillatory behaviors of CD4+ T cells that arise from nonlinear dynamics interwoven across three scales. We identified the inflamed lymph node's ability to retain naive CD4+ T cells as a key mechanism in generating these emergent behaviors. We envision our model and the generic framework encompassing it to serve as a tool for understanding cellular and molecular immunological problems through the lens of systems immunology.

SysMod: the ISCB community for data-driven computational modelling and multi-scale analysis of biological systems.

Computational models of biological systems can exploit a broad range of rapidly developing approaches, including novel experimental approaches, bioinformatics data analysis, emerging modelling paradigms, data standards and algorithms. A discussion about the most recent advances among experts from various domains is crucial to foster data-driven computational modelling and its growing use in assessing and predicting the behaviour of biological systems. Intending to encourage the development of tools, approaches and predictive models, and to deepen our understanding of biological systems, the Community of Special Interest (COSI) was launched in Computational Modelling of Biological Systems (SysMod) in 2016. SysMod's main activity is an annual meeting at the Intelligent Systems for Molecular Biology (ISMB) conference, which brings together computer scientists, biologists, mathematicians, engineers, computational and systems biologists. In the five years since its inception, SysMod has evolved into a dynamic and expanding community, as the increasing number of contributions and participants illustrate. SysMod maintains several online resources to facilitate interaction among the community members, including an online forum, a calendar of relevant meetings and a YouTube channel with talks and lectures of interest for the modelling community. For more than half a decade, the growing interest in computational systems modelling and multi-scale data integration has inspired and supported the SysMod community. Its members get progressively more involved and actively contribute to the annual COSI meeting and several related community workshops and meetings, focusing on specific topics, including particular techniques for computational modelling or standardisation efforts.

Integrative network analyses of transcriptomics data reveal potential drug targets for acute radiation syndrome.

Recent political unrest has highlighted the importance of understanding the short- and long-term effects of gamma-radiation exposure on human health and survivability. In this regard, effective treatment for acute radiation syndrome (ARS) is a necessity in cases of nuclear disasters. Here, we propose 20 therapeutic targets for ARS identified using a systematic approach that integrates gene coexpression networks obtained under radiation treatment in humans and mice, drug databases, disease-gene association, radiation-induced differential gene expression, and literature mining. By selecting gene targets with existing drugs, we identified potential candidates for drug repurposing. Eight of these genes (BRD4, NFKBIA, CDKN1A, TFPI, MMP9, CBR1, ZAP70, IDH3B) were confirmed through literature to have shown radioprotective effect upon perturbation. This study provided a new perspective for the treatment of ARS using systems-level gene associations integrated with multiple biological information. The identified genes might provide high confidence drug target candidates for potential drug repurposing for ARS.

Teaching Metabolism in Upper-Division Undergraduate Biochemistry Courses using Online Computational Systems and Dynamical Models Improves Student Performance.

Understanding metabolic function requires knowledge of the dynamics, interdependence, and regulation of metabolic networks. However, multiple professional societies have recognized that most undergraduate biochemistry students acquire only a surface-level understanding of metabolism. We hypothesized that guiding students through interactive computer simulations of metabolic systems would increase their ability to recognize how individual interactions between components affect the behavior of a system under different conditions. The computer simulations were designed with an interactive activity (i.e., module) that used the predict-observe-explain model of instruction to guide students through a process in which they iteratively predict outcomes, test their predictions, modify the interactions of the system, and then retest the outcomes. We found that biochemistry students using modules performed better on metabolism questions compared with students who did not use the modules. The average learning gain was 8% with modules and 0% without modules, a small to medium effect size. We also confirmed that the modules did not create or reinforce a gender bias. Our modules provide instructors with a dynamic, systems-driven approach to help students learn about metabolic regulation and equip students with important cognitive skills, such as interpreting and analyzing simulation results, and technical skills, such as building and simulating computer-based models.