Search

IntroductionThe 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. MethodsExtensive 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.ResultsResults 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. DiscussionThe 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.

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.

Abstract 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.

Abstract 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.

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.

Abstract CD4+ T cells provide adaptive immunity against pathogens and abnormal cells, and they are also associated with various immune-related diseases. CD4+ T cells’ metabolism is dysregulated in these pathologies and represents an opportunity for drug discovery and development. Genome-scale metabolic modeling offers an opportunity to accelerate drug discovery by providing high-quality information about possible target space in the context of a modeled disease. Here, we develop genome-scale models of naïve, Th1, Th2, and Th17 CD4+ T-cell subtypes to map metabolic perturbations in rheumatoid arthritis, multiple sclerosis, and primary biliary cholangitis. We subjected these models to in silico simulations for drug response analysis of existing FDA-approved drugs and compounds. Integration of disease-specific differentially expressed genes with altered reactions in response to metabolic perturbations identified 68 drug targets for the three autoimmune diseases. In vitro experimental validation, together with literature-based evidence, showed that modulation of fifty percent of identified drug targets suppressed CD4+ T cells, further increasing their potential impact as therapeutic interventions. Our approach can be generalized in the context of other diseases, and the metabolic models can be further used to dissect CD4+ T-cell metabolism.

Abstract Systems biology has experienced dramatic growth in the number, size, and complexity of computational models. To reproduce simulation results and reuse models, researchers must exchange unambiguous model descriptions. We review the latest edition of the Systems Biology Markup Language (SBML), a format designed for this purpose. A community of modelers and software authors developed SBML Level 3 over the past decade. Its modular form consists of a core suited to representing reaction‐based models and packages that extend the core with features suited to other model types including constraint‐based models, reaction‐diffusion models, logical network models, and rule‐based models. The format leverages two decades of SBML and a rich software ecosystem that transformed how systems biologists build and interact with models. More recently, the rise of multiscale models of whole cells and organs, and new data sources such as single‐cell measurements and live imaging, has precipitated new ways of integrating data with models. We provide our perspectives on the challenges presented by these developments and how SBML Level 3 provides the foundation needed to support this evolution.

Abstract Background Computational modeling is an increasingly common practice for disciplinary experts and therefore necessitates integration into science curricula. Computational models afford an opportunity for students to investigate the dynamics of biological systems, but there is significant gap in our knowledge of how these activities impact student knowledge of the structures, relationships, and dynamics of the system. We investigated how a computational modeling activity affected introductory biology students’ mental models of a prokaryotic gene regulatory system (lac operon) by analyzing conceptual models created before and after the activity. Results Students’ pre-lesson conceptual models consisted of provided, system-general structures (e.g., activator, repressor) connected with predominantly incorrect relationships, representing an incomplete mental model of gene regulation. Students’ post-lesson conceptual models included more context-specific structures (e.g., cAMP, lac repressor) and increased in total number of structures and relationships. Student conceptual models also included higher quality relationships among structures, indicating they learned about these context-specific structures through integration with their expanding mental model rather than in isolation. Conclusions Student mental models meshed structures in a manner indicative of knowledge accretion while they were productively re-constructing their understanding of gene regulation. Conceptual models can inform instructors about how students are relating system structures and whether students are developing more sophisticated models of system-general and system-specific dynamics.

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.

Graphical Abstract

Teachers’ integration of the Next Generation Science Standards and corresponding Science and Engineering Practices (SEPs) illustrate current science education reform in the United States. Effective teacher professional development (PD) on SEPs is essential for reform success. In this study, we evaluated the Nebraska STEM Education Conference, a PD program for middle school, high school, and first- and second-year post-secondary STEM teachers. This SEP-oriented PD program focused predominantly on the SEPs ‘developing and using models’ and ‘using mathematics and computational thinking.’ An electronic survey was used to measure participants’ (n = 45) prior integration of SEPs, influential factors and barriers to using SEPs, and changes to interest and confidence in using SEPs as a result of attending the PD program. Our results showed that teachers had limited prior use of SEPs in their teaching. Student interest and learning outcomes were the factors found to be most influential to teachers’ use of SEPs, while limited knowledge, confidence, and resources were the most commonly identified barriers. As a result of attending the PD program, participants significantly improved their confidence and interest to incorporate SEPs. We recommend continued SEP-oriented PD to foster successful NGSS integration and to advance reforms in science education.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Francisella tularensis can cause the zoonotic disease tularemia and is partitioned into subspecies due to differences in chromosomal organization and virulence. The subspecies holarctica (type B) is generally considered more clonal than the other subpopulations with moderate virulence compared to the hypervirulent A.I clade. We performed whole genome sequencing (WGS) on six type B strains isolated from the blood of patients with tularemia within a one-year period from the same United States region, to better understand the associated pathogenicity. The WGS data were compared to the prototype strain for this subspecies, specifically FSC200, which was isolated from a patient with tularemia in Europe. These findings revealed 520–528 single nucleotide polymorphisms (SNPs) between the six United States type B strains compared to FSC200, with slightly higher A+T content in the latter strain. In contrast, comparisons between the six type B isolates showed that five of the six type B isolates had only 4–22 SNPs, while one of the strains had 47–53 SNPs. Analysis of SNPs in the core genome for the six United States type B isolates and the FSC200 strain gave similar results, suggesting that some of these mutations may have been nonsynonymous, resulting in altered protein function and pathogenicity.

Triggering an appropriate protective response against invading agents is crucial to the effectiveness of human innate and adaptive immunity. Pathogen recognition and elimination requires integration of a myriad of signals from many different immune cells. For example, T cell functioning is not qualitatively, but quantitatively determined by cellular and humoral signals. Tipping the balance of signals, such that one of these is favored or gains advantage on another one, may impact the plasticity of T cells. This may lead to switching their phenotypes and, ultimately, modulating the balance between proliferating and memory T cells to sustain an appropriate immune response. We hypothesize that, similar to other intracellular processes such as the cell cycle, the process of T cell differentiation is the result of: (i) pleiotropy (pattern) and (ii) magnitude (dosage/concentration) of input signals, as well as (iii) their timing and duration. That is, a flexible, yet robust immune response upon recognition of the pathogen may result from the integration of signals at the right dosage and timing. To investigate and understand how system’s properties such as T cell plasticity and T cell-mediated robust response arise from the interplay between these signals, the use of experimental toolboxes that modulate immune proteins may be explored. Currently available methodologies to engineer T cells and a recently devised strategy to measure protein dosage may be employed to precisely determine, for example, the expression of transcription factors responsible for T cell differentiation into various subtypes. Thus, the immune response may be systematically investigated quantitatively. Here, we provide a perspective of how pattern, dosage and timing of specific signals, called interleukins, may influence T cell activation and differentiation during the course of the immune response. We further propose that interleukins alone cannot explain the phenotype variability observed in T cells. Specifically, we provide evidence that the dosage of intercellular components of both the immune system and the cell cycle regulating cell proliferation may contribute to T cell activation, differentiation, as well as T cell memory formation and maintenance. Altogether, we envision that a qualitative (pattern) and quantitative (dosage) crosstalk between the extracellular milieu and intracellular proteins leads to T cell plasticity and robustness. The understanding of this complex interplay is crucial to predict and prevent scenarios where tipping the balance of signals may be compromised, such as in autoimmunity.

A variety of biological networks can be modeled as logical or Boolean networks. However, a simplification of the reality to binary states of the nodes does not ease the difficulty of analyzing the dynamics of large, complex networks, such as signal transduction networks, due to the exponential dependence of the state space on the number of nodes. This paper considers a recently introduced method for finding a fairly small subnetwork, representing a collection of nodes that determine the states of most other nodes with a reasonable level of entropy. The subnetwork contains the most determinative nodes that yield the highest information gain. One of the goals of this paper is to propose an algorithm for finding a suitable subnetwork size. The information gain is quantified by the so-called determinative power of the nodes, which is obtained via the mutual information, a concept originating in information theory. We find the most determinative nodes for 36 network models available in the online database Cell Collective (http://cellcollective.org). We provide statistical information that indicates a weak correlation between the subnetwork size and other variables, such as network size, or maximum and average determinative power of nodes. We observe that the proportion represented by the subnetwork in comparison to the whole network shows a weak tendency to decrease for larger networks. The determinative power of nodes is weakly correlated to the number of outputs of a node, and it appears to be independent of other topological measures such as closeness or betweenness centrality. Once the subnetwork of the most determinative nodes is identified, we generate a biological function analysis of its nodes for some of the 36 networks. The analysis shows that a large fraction of the most determinative nodes are essential and involved in crucial biological functions. The biological pathway analysis of the most determinative nodes shows that they are involved in important disease pathways.

CD4+ T cells provide cell-mediated immunity in response to various antigens. During an immune response, naïve CD4+ T cells differentiate into specialized effector T helper (Th1, Th2, and Th17) cells and induced regulatory (iTreg) cells based on a cytokine milieu. In recent studies, complex phenotypes resembling more than one classical T cell lineage have been experimentally observed. Herein, we sought to characterize the capacity of T cell differentiation in response to the complex extracellular environment. We constructed a comprehensive mechanistic (logical) computational model of the signal transduction that regulates T cell differentiation. The model’s dynamics were characterized and analyzed under 511 different environmental conditions. Under these conditions, the model predicted the classical as well as the novel complex (mixed) T cell phenotypes that can co-express transcription factors (TFs) related to multiple differentiated T cell lineages. Analyses of the model suggest that the lineage decision is regulated by both compositions and dosage of signals that constitute the extracellular environment. In this regard, we first characterized the specific patterns of extracellular environments that result in novel T cell phenotypes. Next, we predicted the inputs that can regulate the transition between the canonical and complex T cell phenotypes in a dose-dependent manner. Finally, we predicted the optimal levels of inputs that can simultaneously maximize the activity of multiple lineage-specifying TFs and that can drive a phenotype toward one of the co-expressed TFs. In conclusion, our study provides new insights into the plasticity of CD4+ T cell differentiation, and also acts as a tool to design testable hypotheses for the generation of complex T cell phenotypes by various input combinations and dosages.

Analysing models of biological networks typically relies on workflows in which different software tools with sensitive parameters are chained together, many times with additional manual steps. The accessibility and reproducibility of such workflows is challenging, as publications often overlook analysis details, and because some of these tools may be difficult to install, and/or have a steep learning curve. The CoLoMoTo Interactive Notebook provides a unified environment to edit, execute, share, and reproduce analyses of qualitative models of biological networks. This framework combines the power of different technologies to ensure repeatability and to reduce users' learning curve of these technologies. The framework is distributed as a Docker image with the tools ready to be run without any installation step besides Docker, and is available on Linux, macOS, and Microsoft Windows. The embedded computational workflows are edited with a Jupyter web interface, enabling the inclusion of textual annotations, along with the explicit code to execute, as well as the visualization of the results. The resulting notebook files can then be shared and re-executed in the same environment. To date, the CoLoMoTo Interactive Notebook provides access to the software tools GINsim, BioLQM, Pint, MaBoSS, and Cell Collective, for the modeling and analysis of Boolean and multi-valued networks. More tools will be included in the future. We developed a Python interface for each of these tools to offer a seamless integration in the Jupyter web interface and ease the chaining of complementary analyses.

3D printing represents an emerging technology with significant potential to advance life-science education by allowing students to directly explore the relationship between macromolecular structure and function. In this article and supplemental video guide, we describe our development of a model-based instructional module on DNA supercoiling and outline practical tips for implementing models in undergraduate classrooms. We also present a procedure to design and print 3D dynamic models for classroom use. Furthermore, we describe repositories of 3D biomolecule files to make using models accessible and cost-effective.

While the use of computer tools to simulate complex processes such as computer circuits is normal practice in fields like engineering, the majority of life sciences/biological sciences courses continue to rely on the traditional textbook and memorization approach. To address this issue, we explored the use of the Cell Collective platform as a novel, interactive, and evolving pedagogical tool to foster student engagement, creativity, and higher-level thinking. Cell Collective is a Web-based platform used to create and simulate dynamical models of various biological processes. Students can create models of cells, diseases, or pathways themselves or explore existing models. This technology was implemented in both undergraduate and graduate courses as a pilot study to determine the feasibility of such software at the university level. First, a new (In Silico Biology) class was developed to enable students to learn biology by "building and breaking it" via computer models and their simulations. This class and technology also provide a non-intimidating way to incorporate mathematical and computational concepts into a class with students who have a limited mathematical background. Second, we used the technology to mediate the use of simulations and modeling modules as a learning tool for traditional biological concepts, such as T cell differentiation or cell cycle regulation, in existing biology courses. Results of this pilot application suggest that there is promise in the use of computational modeling and software tools such as Cell Collective to provide new teaching methods in biology and contribute to the implementation of the "Vision and Change" call to action in undergraduate biology education by providing a hands-on approach to biology.

The non-receptor tyrosine kinase Src and receptor tyrosine kinase epidermal growth factor receptor (EGFR/ErbB1) have been established as collaborators in cellular signaling and their combined dysregulation plays key roles in human cancers, including breast cancer. In part due to the complexity of the biochemical network associated with the regulation of these proteins as well as their cellular functions, the role of Src in EGFR regulation remains unclear. Herein we present a new comprehensive, multi-scale dynamical model of ErbB receptor signal transduction in human mammary epithelial cells. This model, constructed manually from published biochemical literature, consists of 245 nodes representing proteins and their post-translational modifications sites, and over 1,000 biochemical interactions. Using computer simulations of the model, we find it is able to reproduce a number of cellular phenomena. Furthermore, the model predicts that overexpression of Src results in increased endocytosis of EGFR in the absence/low amount of the epidermal growth factor (EGF). Our subsequent laboratory experiments also suggest increased internalization of EGFR upon Src overexpression under EGF-deprived conditions, further supporting this model-generated hypothesis.

It has been suggested that irreducible sets of states in Probabilistic Boolean Networks correspond to cellular phenotype. In this study, we identify such sets of states for each phase of the budding yeast cell cycle. We find that these "ergodic sets" underly the cyclin activity levels during each phase of the cell cycle. Our results compare to the observations made in several laboratory experiments as well as the results of differential equation models. Dynamical studies of this model: (i) indicate that under stochastic external signals the continuous oscillating waves of cyclin activity and the opposing waves of CKIs emerge from the logic of a Boolean-based regulatory network without the need for specific biochemical/kinetic parameters; (ii) suggest that the yeast cell cycle network is robust to the varying behavior of cell size (e.g., cell division under nitrogen deprived conditions); (iii) suggest the irreversibility of the Start signal is a function of logic of the G1 regulon, and changing the structure of the regulatory network can render start reversible.

Computational modeling of biological processes is a promising tool in biomedical research. While a large part of its potential lies in the ability to integrate it with laboratory research, modeling currently generally requires a high degree of training in mathematics and/or computer science. To help address this issue, we have developed a web-based tool, Bio-Logic Builder, that enables laboratory scientists to define mathematical representations (based on a discrete formalism) of biological regulatory mechanisms in a modular and non-technical fashion. As part of the user interface, generalized "bio-logic" modules have been defined to provide users with the building blocks for many biological processes. To build/modify computational models, experimentalists provide purely qualitative information about a particular regulatory mechanisms as is generally found in the laboratory. The Bio-Logic Builder subsequently converts the provided information into a mathematical representation described with Boolean expressions/rules. We used this tool to build a number of dynamical models, including a 130-protein large-scale model of signal transduction with over 800 interactions, influenza A replication cycle with 127 species and 200+ interactions, and mammalian and budding yeast cell cycles. We also show that any and all qualitative regulatory mechanisms can be built using this tool.

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.