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.

A practical guide to mechanistic systems modeling in biology using a logic-based approach.

Mechanistic computational models enable the study of regulatory mechanisms implicated in various biological processes. These models provide a means to analyze the dynamics of the systems they describe, and to study and interrogate their properties, and provide insights about the emerging behavior of the system in the presence of single or combined perturbations. Aimed at those who are new to computational modeling, we present here a practical hands-on protocol breaking down the process of mechanistic modeling of biological systems in a succession of precise steps. The protocol provides a framework that includes defining the model scope, choosing validation criteria, selecting the appropriate modeling approach, constructing a model and simulating the model. To ensure broad accessibility of the protocol, we use a logical modeling framework, which presents a lower mathematical barrier of entry, and two easy-to-use and popular modeling software tools: Cell Collective and GINsim. The complete modeling workflow is applied to a well-studied and familiar biological process-the lac operon regulatory system. The protocol can be completed by users with little to no prior computational modeling experience approximately within 3 h.

The status of causality in biological databases: data resources and data retrieval possibilities to support logical modeling.

Causal molecular interactions represent key building blocks used in computational modeling, where they facilitate the assembly of regulatory networks. Logical regulatory networks can be used to predict biological and cellular behaviors by system perturbations and in silico simulations. Today, broad sets of causal interactions are available in a variety of biological knowledge resources. However, different visions, based on distinct biological interests, have led to the development of multiple ways to describe and annotate causal molecular interactions. It can therefore be challenging to efficiently explore various resources of causal interaction and maintain an overview of recorded contextual information that ensures valid use of the data. This review lists the different types of public resources with causal interactions, the different views on biological processes that they represent, the various data formats they use for data representation and storage, and the data exchange and conversion procedures that are available to extract and download these interactions. This may further raise awareness among the targeted audience, i.e. logical modelers and other scientists interested in molecular causal interactions, but also database managers and curators, about the abundance and variety of causal molecular interaction data, and the variety of tools and approaches to convert them into one interoperable resource.

Setting the basis of best practices and standards for curation and annotation of logical models in biology-highlights of the [BC]2 2019 CoLoMoTo/SysMod Workshop.

The fast accumulation of biological data calls for their integration, analysis and exploitation through more systematic approaches. The generation of novel, relevant hypotheses from this enormous quantity of data remains challenging. Logical models have long been used to answer a variety of questions regarding the dynamical behaviours of regulatory networks. As the number of published logical models increases, there is a pressing need for systematic model annotation, referencing and curation in community-supported and standardised formats. This article summarises the key topics and future directions of a meeting entitled 'Annotation and curation of computational models in biology', organised as part of the 2019 [BC]2 conference. The purpose of the meeting was to develop and drive forward a plan towards the standardised annotation of logical models, review and connect various ongoing projects of experts from different communities involved in the modelling and annotation of molecular biological entities, interactions, pathways and models. This article defines a roadmap towards the annotation and curation of logical models, including milestones for best practices and minimum standard requirements.

Metabolic reprogramming in Rheumatoid Arthritis Synovial Fibroblasts: A hybrid modeling approach.

Rheumatoid Arthritis (RA) is an autoimmune disease characterized by a highly invasive pannus formation consisting mainly of Synovial Fibroblasts (RASFs). This pannus leads to cartilage, bone, and soft tissue destruction in the affected joint. RASFs' activation is associated with metabolic alterations resulting from dysregulation of extracellular signals' transduction and gene regulation. Deciphering the intricate mechanisms at the origin of this metabolic reprogramming may provide significant insight into RASFs' involvement in RA's pathogenesis and offer new therapeutic strategies. Qualitative and quantitative dynamic modeling can address some of these features, but hybrid models represent a real asset in their ability to span multiple layers of biological machinery. This work presents the first hybrid RASF model: the combination of a cell-specific qualitative regulatory network with a global metabolic network. The automated framework for hybrid modeling exploits the regulatory network's trap-spaces as additional constraints on the metabolic network. Subsequent flux balance analysis allows assessment of RASFs' regulatory outcomes' impact on their metabolic flux distribution. The hybrid RASF model reproduces the experimentally observed metabolic reprogramming induced by signaling and gene regulation in RASFs. Simulations also enable further hypotheses on the potential reverse Warburg effect in RA. RASFs may undergo metabolic reprogramming to turn into "metabolic factories", producing high levels of energy-rich fuels and nutrients for neighboring demanding cells through the crucial role of HIF1.

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.

COVID19 Disease Map, a computational knowledge repository of virus-host interaction mechanisms.

We need to effectively combine the knowledge from surging literature with complex datasets to propose mechanistic models of SARS-CoV-2 infection, improving data interpretation and predicting key targets of intervention. Here, we describe a large-scale community effort to build an open access, interoperable and computable repository of COVID-19 molecular mechanisms. The COVID-19 Disease Map (C19DMap) is a graphical, interactive representation of disease-relevant molecular mechanisms linking many knowledge sources. Notably, it is a computational resource for graph-based analyses and disease modelling. To this end, we established a framework of tools, platforms and guidelines necessary for a multifaceted community of biocurators, domain experts, bioinformaticians and computational biologists. The diagrams of the C19DMap, curated from the literature, are integrated with relevant interaction and text mining databases. We demonstrate the application of network analysis and modelling approaches by concrete examples to highlight new testable hypotheses. This framework helps to find signatures of SARS-CoV-2 predisposition, treatment response or prioritisation of drug candidates. Such an approach may help deal with new waves of COVID-19 or similar pandemics in the long-term perspective.

Community-driven roadmap for integrated disease maps.

The Disease Maps Project builds on a network of scientific and clinical groups that exchange best practices, share information and develop systems biomedicine tools. The project aims for an integrated, highly curated and user-friendly platform for disease-related knowledge. The primary focus of disease maps is on interconnected signaling, metabolic and gene regulatory network pathways represented in standard formats. The involvement of domain experts ensures that the key disease hallmarks are covered and relevant, up-to-date knowledge is adequately represented. Expert-curated and computer readable, disease maps may serve as a compendium of knowledge, allow for data-supported hypothesis generation or serve as a scaffold for the generation of predictive mathematical models. This article summarizes the 2nd Disease Maps Community meeting, highlighting its important topics and outcomes. We outline milestones on the roadmap for the future development of disease maps, including creating and maintaining standardized disease maps; sharing parts of maps that encode common human disease mechanisms; providing technical solutions for complexity management of maps; and Web tools for in-depth exploration of such maps. A dedicated discussion was focused on mathematical modeling approaches, as one of the main goals of disease map development is the generation of mathematically interpretable representations to predict disease comorbidity or drug response and to suggest drug repositioning, altogether supporting clinical decisions.

Automated inference of Boolean models from molecular interaction maps using CaSQ.

MOTIVATION: Molecular interaction maps have emerged as a meaningful way of representing biological mechanisms in a comprehensive and systematic manner. However, their static nature provides limited insights to the emerging behaviour of the described biological system under different conditions. Computational modelling provides the means to study dynamic properties through in silico simulations and perturbations. We aim to bridge the gap between static and dynamic representations of biological systems with CaSQ, a software tool that infers Boolean rules based on the topology and semantics of molecular interaction maps built with CellDesigner. RESULTS: We developed CaSQ by defining conversion rules and logical formulas for inferred Boolean models according to the topology and the annotations of the starting molecular interaction maps. We used CaSQ to produce executable files of existing molecular maps that differ in size, complexity and the use of Systems Biology Graphical Notation (SBGN) standards. We also compared, where possible, the manually built logical models corresponding to a molecular map to the ones inferred by CaSQ. The tool is able to process large and complex maps built with CellDesigner (either following SBGN standards or not) and produce Boolean models in a standard output format, Systems Biology Marked Up Language-qualitative (SBML-qual), that can be further analyzed using popular modelling tools. References, annotations and layout of the CellDesigner molecular map are retained in the obtained model, facilitating interoperability and model reusability. AVAILABILITY AND IMPLEMENTATION: The present tool is available online: https://lifeware.inria.fr/∼soliman/post/casq/ and distributed as a Python package under the GNU GPLv3 license. The code can be accessed here: https://gitlab.inria.fr/soliman/casq. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Effects of proteasome inhibitors on cytokines, metalloproteinases and their inhibitors and collagen type-I expression in periprosthetic tissues and fibroblasts from loose arthroplasty endoprostheses.

Objective: Aseptic loosening is a major problem in total joint replacement. Implant wear debris provokes a foreign body host response and activates cells to produce a variety of mediators and ROS, leading to periprosthetic osteolysis. Elevated ROS levels can harm proteasome function. Proteasome inhibitors have been reported to alter the secretory profile of cells involved in inflammation and also to induce ROS production. In this work, we aimed to document the effects of proteasome inhibitors MG-132 and Epoxomicin, on the production of factors involved in aseptic loosening, in periprosthetic tissues and fibroblasts, and investigate the role of proteasome impairment in periprosthetic osteolysis. Materials and methods: IL-6 levels in tissue cultures were determined by sandwich ELISA. MMP-1, -3, -13, -14 and TIMP-1 levels in tissue or cell cultures were determined by indirect ELISA. Results for MMP-1 and TIMP-1 in tissue cultures were confirmed by Western blotting. MMP-2 and MMP-9 levels were determined by gelatin zymography. Gene expression of IL-6, MMP-1,-3,-14, TIMP-1 and collagen type-I was determined by RT-PCR. Results: Results show that proteasome inhibition induces the expression of ΜΜΡ-1, -2, -3, -9 and suppresses that of IL-6, MMP-14, -13, TIMP-1 and collagen type I, enhancing the collagenolytic and gelatinolytic activity already present in periprosthetic tissues, as documented in various studies. Conclusions: These findings suggest that proteasome impairment could be a contributing factor to aseptic loosening. Protection and enhancement of proteasome efficacy could thus be considered as an alternative strategy toward disease treatment.

Computational modeling of the main signaling pathways involved in mast cell activation.

A global and rigorous understanding of the signaling pathways and cross-regulatory processes involved in mast cell activation requires the integration of published information with novel functional datasets into a comprehensive computational model. Based on an exhaustive curation of the existing literature and using the software CellDesigner, we have built and annotated a comprehensive molecular map for the FcεRI signaling network. This map can be used to visualize and interpret high-throughput expression data. Furthermore, leaning on this map and using the logical modeling software GINsim, we have derived a qualitative dynamical model, which recapitulates the most salient features of mast cell activation. The resulting logical model can be used to explore the dynamical properties of the system and its responses to different stimuli, in normal or mutant conditions.

Proteomic analysis of the SH2 domain-containing leukocyte protein of 76 kDa (SLP76) interactome in resting and activated primary mast cells [corrected].

We report the first proteomic analysis of the SLP76 interactome in resting and activated primary mouse mast cells. This was made possible by a novel genetic approach used for the first time here. It consists in generating knock-in mice that express signaling molecules bearing a C-terminal tag that has a high affinity for a streptavidin analog. Tagged molecules can be used as molecular baits to affinity-purify the molecular complex in which they are engaged, which can then be studied by mass spectrometry. We examined first SLP76 because, although this cytosolic adapter is critical for both T cell and mast cell activation, its role is well known in T cells but not in mast cells. Tagged SLP76 was expressed in physiological amounts and fully functional in mast cells. We unexpectedly found that SLP76 is exquisitely sensitive to mast cell granular proteases, that Zn(2+)-dependent metalloproteases are especially abundant in mast cells and that they were responsible for SLP76 degradation. Adding a Zn(2+) chelator fully protected SLP76 in mast cell lysates, thereby enabling an efficient affinity-purification of this adapter with its partners. Label-free quantitative mass spectrometry analysis of affinity-purified SLP76 interactomes uncovered both partners already described in T cells and novel partners seen in mast cells only. Noticeably, molecules inducibly recruited in both cell types primarily concur to activation signals, whereas molecules recruited in activated mast cells only are mostly associated with inhibition signals. The transmembrane adapter LAT2, and the serine/threonine kinase with an exchange factor activity Bcr were the most recruited molecules. Biochemical and functional validations established the unexpected finding that Bcr is recruited by SLP76 and positively regulates antigen-induced mast cell activation. Knock-in mice expressing tagged molecules with a normal tissue distribution and expression therefore provide potent novel tools to investigate signalosomes and to uncover novel signaling molecules in mast cells.

MetaLo: metabolic analysis of Logical models extracted from molecular interaction maps.

Molecular interaction maps (MIMs) are static graphical representations depicting complex biochemical networks that can be formalized using one of the Systems Biology Graphical Notation languages. Regardless of their extensive coverage of various biological processes, they are limited in terms of dynamic insights. However, MIMs can serve as templates for developing dynamic computational models. We present MetaLo, an open-source Python package that enables the coupling of Boolean models inferred from process description MIMs with generic core metabolic networks. MetaLo provides a framework to study the impact of signaling cascades, gene regulation processes, and metabolic flux distribution of central energy production pathways. MetaLo computes the Boolean model's asynchronous asymptotic behavior, through the identification of trap-spaces, and extracts metabolic constraints to contextualize the generic metabolic network. MetaLo is able to handle large-scale Boolean models and genome-scale metabolic models without requiring kinetic information or manual tuning. The framework behind MetaLo enables in depth analysis of the regulatory model, and may allow tackling a lack of omics data in poorly addressed biological fields to contextualize generic metabolic networks along with improper automatic reconstructions of cell- and/or disease-specific metabolic networks. MetaLo is available at https://pypi.org/project/metalo/ under the terms of the GNU General Public License v3.

Large-scale computational modelling of the M1 and M2 synovial macrophages in rheumatoid arthritis.

Macrophages play an essential role in rheumatoid arthritis. Depending on their phenotype (M1 or M2), they can play a role in the initiation or resolution of inflammation. The M1/M2 ratio in rheumatoid arthritis is higher than in healthy controls. Despite this, no treatment targeting specifically macrophages is currently used in clinics. Thus, devising strategies to selectively deplete proinflammatory macrophages and promote anti-inflammatory macrophages could be a promising therapeutic approach. State-of-the-art molecular interaction maps of M1 and M2 macrophages in rheumatoid arthritis are available and represent a dense source of knowledge; however, these maps remain limited by their static nature. Discrete dynamic modelling can be employed to study the emergent behaviours of these systems. Nevertheless, handling such large-scale models is challenging. Due to their massive size, it is computationally demanding to identify biologically relevant states in a cell- and disease-specific context. In this work, we developed an efficient computational framework that converts molecular interaction maps into Boolean models using the CaSQ tool. Next, we used a newly developed version of the BMA tool deployed to a high-performance computing cluster to identify the models' steady states. The identified attractors are then validated using gene expression data sets and prior knowledge. We successfully applied our framework to generate and calibrate the M1 and M2 macrophage Boolean models for rheumatoid arthritis. Using KO simulations, we identified NFkB, JAK1/JAK2, and ERK1/Notch1 as potential targets that could selectively suppress proinflammatory macrophages and GSK3B as a promising target that could promote anti-inflammatory macrophages in rheumatoid arthritis.

A large-scale Boolean model of the rheumatoid arthritis fibroblast-like synoviocytes predicts drug synergies in the arthritic joint.

Rheumatoid arthritis (RA) is a complex autoimmune disease with an unknown aetiology. However, rheumatoid arthritis fibroblast-like synoviocytes (RA-FLS) play a significant role in initiating and perpetuating destructive joint inflammation by expressing immuno-modulating cytokines, adhesion molecules, and matrix remodelling enzymes. In addition, RA-FLS are primary drivers of inflammation, displaying high proliferative rates and an apoptosis-resistant phenotype. Thus, RA-FLS-directed therapies could become a complementary approach to immune-directed therapies by predicting the optimal conditions that would favour RA-FLS apoptosis, limit inflammation, slow the proliferation rate and minimise bone erosion and cartilage destruction. In this paper, we present a large-scale Boolean model for RA-FLS that consists of five submodels focusing on apoptosis, cell proliferation, matrix degradation, bone erosion and inflammation. The five-phenotype-specific submodels can be simulated independently or as a global model. In silico simulations and perturbations reproduced the expected biological behaviour of the system under defined initial conditions and input values. The model was then used to mimic the effect of mono or combined therapeutic treatments and predict novel targets and drug candidates through drug repurposing analysis.

Hybrid computational modeling highlights reverse warburg effect in breast cancer-associated fibroblasts.

Cancer-associated fibroblasts (CAFs) are amongst the key players of the tumor microenvironment (TME) and are involved in cancer initiation, progression, and resistance to therapy. They exhibit aggressive phenotypes affecting extracellular matrix remodeling, angiogenesis, immune system modulation, tumor growth, and proliferation. CAFs phenotypic changes appear to be associated with metabolic alterations, notably a reverse Warburg effect that may drive fibroblasts transformation. However, its precise molecular mechanisms and regulatory drivers are still under investigation. Deciphering the reverse Warburg effect in breast CAFs may contribute to a better understanding of the interplay between TME and tumor cells, leading to new treatment strategies. In this regard, dynamic modeling approaches able to span multiple biological layers are essential to capture the emergent properties of various biological entities when complex and intertwined pathways are involved. This work presents the first hybrid large-scale computational model for breast CAFs covering major cellular signaling, gene regulation, and metabolic processes. It was generated by combining a cell- and disease-specific asynchronous Boolean model with a generic core metabolic network leveraging both data-driven and manual curation approaches. This model reproduces the experimentally observed reverse Warburg effect in breast CAFs and further identifies Hypoxia-Inducible Factor 1 (HIF-1) as its key molecular driver. Targeting HIF-1 as part of a TME-centered therapeutic strategy may prove beneficial in the treatment of breast cancer by addressing the reverse Warburg effect. Such findings in CAFs, in light of our previously published results in rheumatoid arthritis synovial fibroblasts, point to a common HIF-1-driven metabolic reprogramming of fibroblasts in breast cancer and rheumatoid arthritis.

A guide for developing comprehensive systems biology maps of disease mechanisms: planning, construction and maintenance.

As a conceptual model of disease mechanisms, a disease map integrates available knowledge and is applied for data interpretation, predictions and hypothesis generation. It is possible to model disease mechanisms on different levels of granularity and adjust the approach to the goals of a particular project. This rich environment together with requirements for high-quality network reconstruction makes it challenging for new curators and groups to be quickly introduced to the development methods. In this review, we offer a step-by-step guide for developing a disease map within its mainstream pipeline that involves using the CellDesigner tool for creating and editing diagrams and the MINERVA Platform for online visualisation and exploration. We also describe how the Neo4j graph database environment can be used for managing and querying efficiently such a resource. For assessing the interoperability and reproducibility we apply FAIR principles.

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.

Boolean modelling as a logic-based dynamic approach in systems medicine.

Molecular mechanisms of health and disease are often represented as systems biology diagrams, and the coverage of such representation constantly increases. These static diagrams can be transformed into dynamic models, allowing for in silico simulations and predictions. Boolean modelling is an approach based on an abstract representation of the system. It emphasises the qualitative modelling of biological systems in which each biomolecule can take two possible values: zero for absent or inactive, one for present or active. Because of this approximation, Boolean modelling is applicable to large diagrams, allowing to capture their dynamic properties. We review Boolean models of disease mechanisms and compare a range of methods and tools used for analysis processes. We explain the methodology of Boolean analysis focusing on its application in disease modelling. Finally, we discuss its practical application in analysing signal transduction and gene regulatory pathways in health and disease.