Dynamic Boolean modeling of molecular and cellular interactions in psoriasis predicts drug target candidates.

Psoriasis arises from complex interactions between keratinocytes and immune cells, leading to uncontrolled inflammation, immune hyperactivation, and a perturbed keratinocyte life cycle. Despite the availability of drugs for psoriasis management, the disease remains incurable. Treatment response variability calls for new tools and approaches to comprehend the mechanisms underlying disease development. We present a Boolean multiscale population model that captures the dynamics of cell-specific phenotypes in psoriasis, integrating discrete logical formalism and population dynamics simulations. Through simulations and network analysis, the model predictions suggest that targeting neutrophil activation in conjunction with inhibition of either prostaglandin E2 (PGE2) or STAT3 shows promise comparable to interleukin-17 (IL-17) inhibition, one of the most effective treatment options for moderate and severe cases. Our findings underscore the significance of considering complex intercellular interactions and intracellular signaling in psoriasis and highlight the importance of computational approaches in unraveling complex biological systems for drug target identification.

TFCheckpoint database update, a cross-referencing system for transcription factors from human, mouse and rat.

Prior knowledge about DNA-binding transcription factors (dbTFs), transcription co-regulators (coTFs) and general transcriptional factors (GTFs) is crucial for the study and understanding of the regulation of transcription. This is reflected by the many publications and database resources describing knowledge about TFs. We previously launched the TFCheckpoint database, an integrated resource focused on human, mouse and rat dbTFs, providing users access to a comprehensive overview of these proteins. Here, we describe TFCheckpoint 2.0 (https://www.tfcheckpoint.org/index.php), comprising 13 collections of dbTFs, coTFs and GTFs. TFCheckpoint 2.0 provides an easy and versatile cross-referencing system for users to view and download collections that may otherwise be cumbersome to find, compare and retrieve.

Characterisation of Colorectal Cancer Cell Lines through Proteomic Profiling of Their Extracellular Vesicles.

Colorectal cancer (CRC) is one of the most prevalent cancers, driven by several factors including deregulations in intracellular signalling pathways. Small extracellular vesicles (sEVs) are nanosized protein-packaged particles released from cells, which are present in liquid biopsies. Here, we characterised the proteome landscape of sEVs and their cells of origin in three CRC cell lines HCT116, HT29 and SW620 to explore molecular traits that could be exploited as cancer biomarker candidates and how intracellular signalling can be assessed by sEV analysis instead of directly obtaining the cell of origin itself. Our findings revealed that sEV cargo clearly reflects its cell of origin with proteins of the PI3K-AKT pathway highly represented in sEVs. Proteins known to be involved in CRC were detected in both cells and sEVs including KRAS, ARAF, mTOR, PDPK1 and MAPK1, while TGFB1 and TGFBR2, known to be key players in epithelial cancer carcinogenesis, were found to be enriched in sEVs. Furthermore, the phosphopeptide-enriched profiling of cell lysates demonstrated a distinct pattern between cell lines and highlighted potential phosphoproteomic targets to be investigated in sEVs. The total proteomic and phosphoproteomics profiles described in the current work can serve as a source to identify candidates for cancer biomarkers that can potentially be assessed from liquid biopsies.

Halogenated boroxine increases propensity to apoptosis in leukemia (UT-7) but not non-tumor cells in vitro.

A hallmark of the development of solid and hematological malignancies is the dysregulation of apoptosis, which leads to an imbalance between cell proliferation, cell survival and death. Halogenated boroxine [K2 (B3 O3 F4 OH)] (HB) is a derivative of cyclic anhydride of boronic acid, with reproducible anti-tumor and anti-proliferative effects in different cell models. Notably, these changes are observed to be more profound in tumor cells than in normal cells. Here, we investigated the underlying mechanisms through an extensive evaluation of (a) deregulated target genes and (b) their interactions and links with main apoptotic pathway genes upon treatment with an optimized concentration of HB. To provide deeper insights into the mechanism of action of HB, we performed identification, visualization, and pathway association of differentially expressed genes (DEGs) involved in regulation of apoptosis among tumor and non-tumor cells upon HB treatment. We report that HB at a concentration of 0.2 mg·mL-1 drives tumor cells to apoptosis, whereas non-tumor cells are not affected. Comparison of DEG profiles, gene interactions and pathway associations suggests that the HB effect and tumor-'selectivity' can be explained by Bax/Bak-independent mitochondrial depolarization by ROS generation and TRAIL-like activation, followed by permanent inhibition of NFκB signaling pathway specifically in tumor cells.

Boolean function metrics can assist modelers to check and choose logical rules.

Computational models of biological processes provide one of the most powerful methods for a detailed analysis of the mechanisms that drive the behavior of complex systems. Logic-based modeling has enhanced our understanding and interpretation of those systems. Defining rules that determine how the output activity of biological entities is regulated by their respective inputs has proven to be challenging. Partly this is because of the inherent noise in data that allows multiple model parameterizations to fit the experimental observations, but some of it is also due to the fact that models become increasingly larger, making the use of automated tools to assemble the underlying rules indispensable. We present several Boolean function metrics that provide modelers with the appropriate framework to analyze the impact of a particular model parameterization. We demonstrate the link between a semantic characterization of a Boolean function and its consistency with the model's underlying regulatory structure. We further define the properties that outline such consistency and show that several of the Boolean functions under study violate them, questioning their biological plausibility and subsequent use. We also illustrate that regulatory functions can have major differences with regard to their asymptotic output behavior, with some of them being biased towards specific Boolean outcomes when others are dependent on the ratio between activating and inhibitory regulators. Application results show that in a specific signaling cancer network, the function bias can be used to guide the choice of logical operators for a model that matches data observations. Moreover, graph analysis indicates that commonly used Boolean functions become more biased with increasing numbers of regulators, supporting the idea that rule specification can effectively determine regulatory outcome despite the complex dynamics of biological networks.

The gene regulation knowledge commons: the action area of GREEKC.

As computational modeling becomes more essential to analyze and understand biological regulatory mechanisms, governance of the many databases and knowledge bases that support this domain is crucial to guarantee reliability and interoperability of resources. To address this, the COST Action Gene Regulation Ensemble Effort for the Knowledge Commons (GREEKC, CA15205, www.greekc.org) organized nine workshops in a four-year period, starting September 2016. The workshops brought together a wide range of experts from all over the world working on various steps in the knowledge management process that focuses on understanding gene regulatory mechanisms. The discussions between ontologists, curators, text miners, biologists, bioinformaticians, philosophers and computational scientists spawned a host of activities aimed to standardize and update existing knowledge management workflows and involve end-users in the process of designing the Gene Regulation Knowledge Commons (GRKC). Here the GREEKC consortium describes its main achievements in improving this GRKC.

ExTRI: Extraction of transcription regulation interactions from literature.

The regulation of gene transcription by transcription factors is a fundamental biological process, yet the relations between transcription factors (TF) and their target genes (TG) are still only sparsely covered in databases. Text-mining tools can offer broad and complementary solutions to help locate and extract mentions of these biological relationships in articles. We have generated ExTRI, a knowledge graph of TF-TG relationships, by applying a high recall text-mining pipeline to MedLine abstracts identifying over 100,000 candidate sentences with TF-TG relations. Validation procedures indicated that about half of the candidate sentences contain true TF-TG relationships. Post-processing identified 53,000 high confidence sentences containing TF-TG relationships, with a cross-validation F1-score close to 75%. The resulting collection of TF-TG relationships covers 80% of the relations annotated in existing databases. It adds 11,000 other potential interactions, including relationships for ~100 TFs currently not in public TF-TG relation databases. The high confidence abstract sentences contribute 25,000 literature references not available from other resources and offer a wealth of direct pointers to functional aspects of the TF-TG interactions. Our compiled resource encompassing ExTRI together with publicly available resources delivers literature-derived TF-TG interactions for more than 900 of the 1500-1600 proteins considered to function as specific DNA binding TFs. The obtained result can be used by curators, for network analysis and modelling, for causal reasoning or knowledge graph mining approaches, or serve to benchmark text mining strategies.

Dealing with different conceptions of pollution in the Gene Regulation Knowledge Commons.

Current research of gene regulatory mechanisms is increasingly dependent on the availability of high-quality information from manually curated databases. Biocurators undertake the task of extracting knowledge claims from scholarly publications, organizing these claims in a meaningful format and making them computable. In doing so, they enhance the value of existing scientific knowledge by making it accessible to the users of their databases. In this capacity, biocurators are well positioned to identify and weed out information that is of insufficient quality. The criteria that define information quality are typically outlined in curation guidelines developed by biocurators. These guidelines have been prudently developed to reflect the needs of the user community the database caters to. The guidelines depict the standard evidence that this community recognizes as sufficient justification for trustworthy data. Additionally, these guidelines determine the process by which data should be organized and maintained to be valuable to users. Following these guidelines, biocurators assess the quality, reliability, and validity of the information they encounter. In this article we explore to what extent different use cases agree with the inclusion criteria that define positive and negative data, implemented by the database. What are the drawbacks to users who have queries that would be well served by results that fall just short of the criteria used by a database? Finally, how can databases (and biocurators) accommodate the needs of such more explorative use cases?

Logical and experimental modeling of cytokine and eicosanoid signaling in psoriatic keratinocytes.

Psoriasis is a chronic skin disease, in which immune cells and keratinocytes keep each other in a state of inflammation. It is believed that phospholipase A2 (PLA2)-dependent eicosanoid release plays a key role in this. T-helper (Th) 1-derived cytokines are established activators of phospholipases in keratinocytes, whereas Th17-derived cytokines have largely unknown effects. Logical model simulations describing the function of cytokine and eicosanoid signaling networks combined with experimental data suggest that Th17 cytokines stimulate proinflammatory cytokine expression in psoriatic keratinocytes via activation of cPLA2α-Prostaglandin E2-EP4 signaling, which could be suppressed using the anti-psoriatic calcipotriol. cPLA2α inhibition and calcipotriol distinctly regulate expression of key psoriatic genes, possibly offering therapeutic advantage when applied together. Model simulations additionally suggest EP4 and protein kinase cAMP-activated catalytic subunit alpha as drug targets that may restore a normal phenotype. Our work illustrates how the study of complex diseases can benefit from an integrated systems approach.

A GO catalogue of human DNA-binding transcription factors.

To control gene transcription, DNA-binding transcription factors recognise specific sequence motifs in gene regulatory regions. A complete and reliable GO annotation of all DNA-binding transcription factors is key to investigating the delicate balance of gene regulation in response to environmental and developmental stimuli. The need for such information is demonstrated by the many lists of transcription factors that have been produced over the past decade. The COST Action Gene Regulation Ensemble Effort for the Knowledge Commons (GREEKC) Consortium brought together experts in the field of transcription with the aim of providing high quality and interoperable gene regulatory data. The Gene Ontology (GO) Consortium provides strict definitions for gene product function, including factors that regulate transcription. The collaboration between the GREEKC and GO Consortia has enabled the application of those definitions to produce a new curated catalogue of over 1400 human DNA-binding transcription factors, that can be accessed at https://www.ebi.ac.uk/QuickGO/targetset/dbTF. This catalogue has facilitated an improvement in the GO annotation of human DNA-binding transcription factors and led to the GO annotation of almost sixty thousand DNA-binding transcription factors in over a hundred species. Thus, this work will aid researchers investigating the regulation of transcription in both biomedical and basic science.

Formalization of gene regulation knowledge using ontologies and gene ontology causal activity models.

Gene regulation computational research requires handling and integrating large amounts of heterogeneous data. The Gene Ontology has demonstrated that ontologies play a fundamental role in biological data interoperability and integration. Ontologies help to express data and knowledge in a machine processable way, which enables complex querying and advanced exploitation of distributed data. Contributing to improve data interoperability in gene regulation is a major objective of the GREEKC Consortium, which aims to develop a standardized gene regulation knowledge commons. GREEKC proposes the use of ontologies and semantic tools for developing interoperable gene regulation knowledge models, which should support data annotation. In this work, we study how such knowledge models can be generated from cartoons of gene regulation scenarios. The proposed method consists of generating descriptions in natural language of the cartoons; extracting the entities from the texts; finding those entities in existing ontologies to reuse as much content as possible, especially from well known and maintained ontologies such as the Gene Ontology, the Sequence Ontology, the Relations Ontology and ChEBI; and implementation of the knowledge models. The models have been implemented using Protégé, a general ontology editor, and Noctua, the tool developed by the Gene Ontology Consortium for the development of causal activity models to capture more comprehensive annotations of genes and link their activities in a causal framework for Gene Ontology Annotations. We applied the method to two gene regulation scenarios and illustrate how to apply the models generated to support the annotation of data from research articles.

Gene Ontology representation for transcription factor functions.

Transcription plays a central role in defining the identity and functionalities of cells, as well as in their responses to changes in the cellular environment. The Gene Ontology (GO) provides a rigorously defined set of concepts that describe the functions of gene products. A GO annotation is a statement about the function of a particular gene product, represented as an association between a gene product and the biological concept a GO term defines. Critically, each GO annotation is based on traceable scientific evidence. Here, we describe the different GO terms that are associated with proteins involved in transcription and its regulation, focusing on the standard of evidence required to support these associations. This article is intended to help users of GO annotations understand how to interpret the annotations and can contribute to the consistency of GO annotations. We distinguish between three classes of activities involved in transcription or directly regulating it - general transcription factors, DNA-binding transcription factors, and transcription co-regulators.

CausalBuilder: bringing the MI2CAST causal interaction annotation standard to the curator.

Molecular causal interactions are defined as regulatory connections between biological components. They are commonly retrieved from biological experiments and can be used for connecting biological molecules together to enable the building of regulatory computational models that represent biological systems. However, including a molecular causal interaction in a model requires assessing its relevance to that model, based on the detailed knowledge about the biomolecules, interaction type and biological context. In order to standardize the representation of this knowledge in 'causal statements', we recently developed the Minimum Information about a Molecular Interaction Causal Statement (MI2CAST) guidelines. Here, we introduce causalBuilder: an intuitive web-based curation interface for the annotation of molecular causal interactions that comply with the MI2CAST standard. The causalBuilder prototype essentially embeds the MI2CAST curation guidelines in its interface and makes its rules easy to follow by a curator. In addition, causalBuilder serves as an original application of the Visual Syntax Method general-purpose curation technology and provides both curators and tool developers with an interface that can be fully configured to allow focusing on selected MI2CAST concepts to annotate. After the information is entered, the causalBuilder prototype produces genuine causal statements that can be exported in different formats.

The Minimum Information about a Molecular Interaction CAusal STatement (MI2CAST).

MOTIVATION: A large variety of molecular interactions occurs between biomolecular components in cells. When a molecular interaction results in a regulatory effect, exerted by one component onto a downstream component, a so-called 'causal interaction' takes place. Causal interactions constitute the building blocks in our understanding of larger regulatory networks in cells. These causal interactions and the biological processes they enable (e.g. gene regulation) need to be described with a careful appreciation of the underlying molecular reactions. A proper description of this information enables archiving, sharing and reuse by humans and for automated computational processing. Various representations of causal relationships between biological components are currently used in a variety of resources. RESULTS: Here, we propose a checklist that accommodates current representations, called the Minimum Information about a Molecular Interaction CAusal STatement (MI2CAST). This checklist defines both the required core information, as well as a comprehensive set of other contextual details valuable to the end user and relevant for reusing and reproducing causal molecular interaction information. The MI2CAST checklist can be used as reporting guidelines when annotating and curating causal statements, while fostering uniformity and interoperability of the data across resources. AVAILABILITY AND IMPLEMENTATION: The checklist together with examples is accessible at https://github.com/MI2CAST/MI2CAST. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

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.

UniBioDicts: Unified access to Biological Dictionaries.

SUMMARY: We present a set of software packages that provide uniform access to diverse biological vocabulary resources that are instrumental for current biocuration efforts and tools. The Unified Biological Dictionaries (UniBioDicts or UBDs) provide a single query-interface for accessing the online API services of leading biological data providers. Given a search string, UBDs return a list of matching term, identifier and metadata units from databases (e.g. UniProt), controlled vocabularies (e.g. PSI-MI) and ontologies (e.g. GO, via BioPortal). This functionality can be connected to input fields (user-interface components) that offer autocomplete lookup for these dictionaries. UBDs create a unified gateway for accessing life science concepts, helping curators find annotation terms across resources (based on descriptive metadata and unambiguous identifiers), and helping data users search and retrieve the right query terms. AVAILABILITY AND IMPLEMENTATION: The UBDs are available through npm and the code is available in the GitHub organisation UniBioDicts (https://github.com/UniBioDicts) under the Affero GPL license. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

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.

A Middle-Out Modeling Strategy to Extend a Colon Cancer Logical Model Improves Drug Synergy Predictions in Epithelial-Derived Cancer Cell Lines.

Cancer is a heterogeneous and complex disease and one of the leading causes of death worldwide. The high tumor heterogeneity between individuals affected by the same cancer type is accompanied by distinct molecular and phenotypic tumor profiles and variation in drug treatment response. In silico modeling of cancer as an aberrantly regulated system of interacting signaling molecules provides a basis to enhance our biological understanding of disease progression, and it offers the means to use computer simulations to test and optimize drug therapy designs on particular cancer types and subtypes. This sets the stage for precision medicine: the design of treatments tailored to individuals or groups of patients based on their tumor-specific molecular cancer profiles. Here, we show how a relatively large manually curated logical model can be efficiently enhanced further by including components highlighted by a multi-omics data analysis of data from Consensus Molecular Subtypes covering colorectal cancer. The model expansion was performed in a pathway-centric manner, following a partitioning of the model into functional subsystems, named modules. The resulting approach constitutes a middle-out modeling strategy enabling a data-driven expansion of a model from a generic and intermediate level of molecular detail to a model better covering relevant processes that are affected in specific cancer subtypes, comprising 183 biological entities and 603 interactions between them, partitioned in 25 functional modules of varying size and structure. We tested this model for its ability to correctly predict drug combination synergies, against a dataset of experimentally determined cell growth responses with 18 drugs in all combinations, on eight cancer cell lines. The results indicate that the extended model had an improved accuracy for drug synergy prediction for the majority of the experimentally tested cancer cell lines, although significant improvements of the model's predictive performance are still needed. Our study demonstrates how a tumor-data driven middle-out approach toward refining a logical model of a biological system can further customize a computer model to represent specific cancer cell lines and provide a basis for identifying synergistic effects of drugs targeting specific regulatory proteins. This approach bridges between preclinical cancer model data and clinical patient data and may thereby ultimately be of help to develop patient-specific in silico models that can steer treatment decisions in the clinic.

A foundation for reference models for drug combinations with an application to Loewe's reference model.

BACKGROUND: Treating patients with combinations of drugs that have synergistic effects has become widespread practice in the clinic. Drugs work synergistically when the observed effect of a drug combination is larger than the effect predicted by the reference model. The reference model is a theoretical null model that returns the combined effect of given doses of drugs under the assumption that these drugs do not interact. There is ongoing debate on what it means for drugs to not interact. The controversy transcends mathematical punctuality, as different non-interaction principles result in different reference models. A famous reference model that has been in existence for already a long time is Loewe's reference model. Loewe's vision on non-interaction was purely intuitive: two drugs do not interact if all combinations of doses that result in a certain given effect lie on a straight line. RESULTS: We show that Loewe's reference model can be obtained from much more fundamental principles. First, we introduce the new notion of complementary dose. Secondly, we reformulate the existing concept of equivalent dose, whereby our formulation is more general than existing ones. Finally, a very general non-interaction principle is put forward. The proposed non-interaction principle represents a certain interplay between complementary and equivalent doses: drugs are non-interacting if complementarity is preserved under equivalence. It is then shown that Loewe's reference model naturally follows from these principles by an appropriate choice of complementarity. CONCLUSIONS: The presented work increases insight into Loewe's reference model for drug combinations, which is realized by the introduction of a very general non-interaction principle that does not refer to any specific dose-response curve, nor to any property of applicable dose-response curves.

Network Building with the Cytoscape BioGateway App Explained in Five Use Cases.

The BioGateway App is a plugin for the Cytoscape network editor, allowing users to interactively build biological networks by querying the Biogateway Resource Description Framework (RDF) triple store. BioGateway contains information from several curated resources including UniProtKB, IntAct, Gene Ontology Annotations, various datasets containing transcription-factor regulatory relations to specific target genes, and more. The BioGateway App facilitates the step-by-step creation of complex SPARQL queries through an intuitive Graphical User Interface, allowing users to build and explore biological interaction networks to assess, among other things, gene regulatory relationships, gene ontology annotations, and protein-protein interactions. As the BioGateway information content is most abundant for human proteins and genes, this article describes the utility of the tool through a series of use cases on these human data, starting from the most basic levels and then detailing applications that address some of the rich complexity of the integrated data. Network refinement and display can be further optimized via the selection and filtering possibilities that the Cytoscape framework provides. The use cases also provide examples to explore network information in other species, as they become supported by BioGateway. © 2020 The Authors. Basic Protocol 1: Introducing a node from the canvas Basic Protocol 2: Introducing a node from the query builder Basic Protocol 3: Exploring molecular relationships between diseases Basic Protocol 4: Find proteins with protein kinase activity involved in a disease and explore the context around them Basic Protocol 5: Exploring the potential downstream effects after targeted inhibition of proteins Support Protocol: Installation of the BioGateway plugin through the Cytoscape App Manager and from source.