Use of logic theory in understanding regulatory pathway signaling in response to infection

BDEtools: A MATLAB Package for Boolean Delay Equation Modeling

Boolean Delay Equations (BDEs) can simulate surprisingly complex behavior, despite their relative simplicity. In addition to steady-state dynamics, BDEs can also generate periodic and quasiperiodic oscillations, m:n frequency locking, and even chaos. Further, the enumerability of Boolean update functions and their compact parametrization means that BDEs can be leveraged to generate low-level descriptions of biological networks, from which more detailed formulations (e.g., differential equation models) can be constructed. However, although several studies have demonstrated the utility of BDE modeling in computational biology, a current barrier to the wider adoption of the BDE approach is the absence of freely available simulation software. In this work, we present BDEtools-an open-source MATLAB package for numerically solving BDE models. After giving a brief introduction to BDE modeling, we describe the package's solver algorithms, together with several utility functions that can be used to provide solver inputs and to process solver outputs. We also demonstrate the functionality of BDEtools by illustrating its application to an established model of a gene regulatory network that controls circadian rhythms. BDEtools makes it straightforward for researchers to quickly build reliable BDE models of biological networks. We hope that its ease of use and free availability will encourage more researchers to explore BDE formulations of their systems of interest. Through the continued use of BDEs by the computational biology community, we will, no doubt, discover their potential applicability to a broader class of biological networks.

nSeP: immune and metabolic biomarkers for early detection of neonatal sepsis-protocol for a prospective multicohort study

INTRODUCTION

Diagnosing neonatal sepsis is heavily dependent on clinical phenotyping as culture-positive body fluid has poor sensitivity, and existing blood biomarkers have poor specificity.A combination of machine learning, statistical and deep pathway biology analyses led to the identification of a tripartite panel of biologically connected immune and metabolic markers that showed greater than 99% accuracy for detecting bacterial infection with 100% sensitivity. The cohort study described here is designed as a large-scale clinical validation of this previous work.

METHODS AND ANALYSIS

This multicentre observational study will prospectively recruit a total of 1445 newborn infants (all gestations)-1084 with suspected early-or late-onset sepsis, and 361 controls-over 4 years. A small volume of whole blood will be collected from infants with suspected sepsis at the time of presentation. This sample will be used for integrated transcriptomic, lipidomic and targeted proteomics profiling. In addition, a subset of samples will be subjected to cellular phenotype and proteomic analyses. A second sample from the same patient will be collected at 24 hours, with an opportunistic sampling for stool culture. For control infants, only one set of blood and stool sample will be collected to coincide with clinical blood sampling. Along with detailed clinical information, blood and stool samples will be analysed and the information will be used to identify and validate the efficacy of immune-metabolic networks in the diagnosis of bacterial neonatal sepsis and to identify new host biomarkers for viral sepsis.

ETHICS AND DISSEMINATION

The study has received research ethics committee approval from the Wales Research Ethics Committee 2 (reference 19/WA/0008) and operational approval from Health and Care Research Wales. Submission of study results for publication will involve making available all anonymised primary and processed data on public repository sites.

TRIAL REGISTRATION NUMBER

NCT03777670.

The Interdependency and Co-Regulation of the Vitamin D and Cholesterol Metabolism

Vitamin D and cholesterol metabolism overlap significantly in the pathways that contribute to their biosynthesis. However, our understanding of their independent and co-regulation is limited. Cardiovascular disease is the leading cause of death globally and atherosclerosis, the pathology associated with elevated cholesterol, is the leading cause of cardiovascular disease. It is therefore important to understand vitamin D metabolism as a contributory factor. From the literature, we compile evidence of how these systems interact, relating the understanding of the molecular mechanisms involved to the results from observational studies. We also present the first systems biology pathway map of the joint cholesterol and vitamin D metabolisms made available using the Systems Biology Graphical Notation (SBGN) Markup Language (SBGNML). It is shown that the relationship between vitamin D supplementation, total cholesterol, and LDL-C status, and between latitude, vitamin D, and cholesterol status are consistent with our knowledge of molecular mechanisms. We also highlight the results that cannot be explained with our current knowledge of molecular mechanisms: (i) vitamin D supplementation mitigates the side-effects of statin therapy; (ii) statin therapy does not impact upon vitamin D status; and critically (iii) vitamin D supplementation does not improve cardiovascular outcomes, despite improving cardiovascular risk factors. For (iii), we present a hypothesis, based on observations in the literature, that describes how vitamin D regulates the balance between cellular and plasma cholesterol. Answering these questions will create significant opportunities for advancement in our understanding of cardiovascular health.

Opportunities for multiscale computational modelling of serotonergic drug effects in Alzheimer's disease

Alzheimer's disease (AD) is an age-specific neurodegenerative disease that compromises cognitive functioning and impacts the quality of life of an individual. Pathologically, AD is characterised by abnormal accumulation of beta-amyloid (Aβ) and hyperphosphorylated tau protein. Despite research advances over the last few decades, there is currently still no cure for AD. Although, medications are available to control some behavioural symptoms and slow the disease's progression, most prescribed medications are based on cholinesterase inhibitors. Over the last decade, there has been increased attention towards novel drugs, targeting alternative neurotransmitter pathways, particularly those targeting serotonergic (5-HT) system. In this review, we focused on 5-HT receptor (5-HTR) mediated signalling and drugs that target these receptors. These pathways regulate key proteins and kinases such as GSK-3 that are associated with abnormal levels of Aβ and tau in AD. We then review computational studies related to 5-HT signalling pathways with the potential for providing deeper understanding of AD pathologies. In particular, we suggest that multiscale and multilevel modelling approaches could potentially provide new insights into AD mechanisms, and towards discovering novel 5-HTR based therapeutic targets.

In silico disease model: from simple networks to complex diseases

In this chapter we consider in silico modeling of diseases starting from some simple to some complex (and mathematical) concepts. Examples and applications of in silico modeling for some important categories of diseases (such as for cancers, infectious diseases, and neuronal diseases) are also given.

Is systems pharmacology ready to impact upon therapy development? A study on the cholesterol biosynthesis pathway

Background and Purpose

An ever‐growing wealth of information on current drugs and their pharmacological effects is available from online databases. As our understanding of systems biology increases, we have the opportunity to predict, model and quantify how drug combinations can be introduced that outperform conventional single‐drug therapies. Here, we explore the feasibility of such systems pharmacology approaches with an analysis of the mevalonate branch of the cholesterol biosynthesis pathway.

Experimental Approach

Using open online resources, we assembled a computational model of the mevalonate pathway and compiled a set of inhibitors directed against targets in this pathway. We used computational optimization to identify combination and dose options that show not only maximal efficacy of inhibition on the cholesterol producing branch but also minimal impact on the geranylation branch, known to mediate the side effects of pharmaceutical treatment.

Key Results

We describe serious impediments to systems pharmacology studies arising from limitations in the data, incomplete coverage and inconsistent reporting. By curating a more complete dataset, we demonstrate the utility of computational optimization for identifying multi‐drug treatments with high efficacy and minimal off‐target effects.

Conclusion and Implications

We suggest solutions that facilitate systems pharmacology studies, based on the introduction of standards for data capture that increase the power of experimental data. We propose a systems pharmacology workflow for the refinement of data and the generation of future therapeutic hypotheses.

Exhaustive search of linear information encoding protein-peptide recognition

High-throughput in vitro methods have been extensively applied to identify linear information that encodes peptide recognition. However, these methods are limited in number of peptides, sequence variation, and length of peptides that can be explored, and often produce solutions that are not found in the cell. Despite the large number of methods developed to attempt addressing these issues, the exhaustive search of linear information encoding protein-peptide recognition has been so far physically unfeasible. Here, we describe a strategy, called DALEL, for the exhaustive search of linear sequence information encoded in proteins that bind to a common partner. We applied DALEL to explore binding specificity of SH3 domains in the budding yeast Saccharomyces cerevisiae. Using only the polypeptide sequences of SH3 domain binding proteins, we succeeded in identifying the majority of known SH3 binding sites previously discovered either in vitro or in vivo. Moreover, we discovered a number of sites with both non-canonical sequences and distinct properties that may serve ancillary roles in peptide recognition. We compared DALEL to a variety of state-of-the-art algorithms in the blind identification of known binding sites of the human Grb2 SH3 domain. We also benchmarked DALEL on curated biological motifs derived from the ELM database to evaluate the effect of increasing/decreasing the enrichment of the motifs. Our strategy can be applied in conjunction with experimental data of proteins interacting with a common partner to identify binding sites among them. Yet, our strategy can also be applied to any group of proteins of interest to identify enriched linear motifs or to exhaustively explore the space of linear information encoded in a polypeptide sequence. Finally, we have developed a webserver located at http://michnick.bcm.umontreal.ca/dalel, offering user-friendly interface and providing different scenarios utilizing DALEL.

Here we describe the first strategy for the exhaustive search of the linear information encoding protein-peptide recognition; an approach that has previously been physically unfeasible because the combinatorial space of polypeptide sequences is too vast. The search covers the entire space of sequences with no restriction on motif length or composition, and includes all possible combinations of amino acids at distinct positions of each sequence, as well as positions with correlated preferences for amino acids.

Digital clocks: simple Boolean models can quantitatively describe circadian systems

The gene networks that comprise the circadian clock modulate biological function across a range of scales, from gene expression to performance and adaptive behaviour. The clock functions by generating endogenous rhythms that can be entrained to the external 24-h day–night cycle, enabling organisms to optimally time biochemical processes relative to dawn and dusk. In recent years, computational models based on differential equations have become useful tools for dissecting and quantifying the complex regulatory relationships underlying the clock's oscillatory dynamics. However, optimizing the large parameter sets characteristic of these models places intense demands on both computational and experimental resources, limiting the scope of in silico studies. Here, we develop an approach based on Boolean logic that dramatically reduces the parametrization, making the state and parameter spaces finite and tractable. We introduce efficient methods for fitting Boolean models to molecular data, successfully demonstrating their application to synthetic time courses generated by a number of established clock models, as well as experimental expression levels measured using luciferase imaging. Our results indicate that despite their relative simplicity, logic models can (i) simulate circadian oscillations with the correct, experimentally observed phase relationships among genes and (ii) flexibly entrain to light stimuli, reproducing the complex responses to variations in daylength generated by more detailed differential equation formulations. Our work also demonstrates that logic models have sufficient predictive power to identify optimal regulatory structures from experimental data. By presenting the first Boolean models of circadian circuits together with general techniques for their optimization, we hope to establish a new framework for the systematic modelling of more complex clocks, as well as other circuits with different qualitative dynamics. In particular, we anticipate that the ability of logic models to provide a computationally efficient representation of system behaviour could greatly facilitate the reverse-engineering of large-scale biochemical networks.

Systems vaccinology: its promise and challenge for HIV vaccine development

PURPOSE OF REVIEW

The use of systems biology approaches to understand and predict vaccine-induced immunity promises to revolutionize vaccinology. For centuries vaccines were developed empirically, with very little understanding of the mechanisms by which they mediate protective immunity. The so-called systems vaccinology approach employs high-throughput technologies (e.g. microarrays, RNA-seq and mass spectrometry-based proteomics and metabolomics) and computational modeling to describe the complex interactions between all the parts of immune system, with a view to elucidating new biological rules capable of predicting the behavior of the system.

RECENT FINDINGS

Systems biology successfully applied to yellow-fever and influenza vaccines has led to the discovery of signatures that predict vaccine immunogenicity, and promises to advance basic immunology research by providing novel mechanistic insights about immune regulation. However a major challenge of systems vaccinology concerns the analyses and interpretation of the large and noisy data sets generated by high-throughput techniques. Overcoming these issues, we envision that systems vaccinology will have a potential impact on vaccine development, including HIV vaccines.

SUMMARY

High-throughput technologies allow the investigation of vaccine-induced immune responses at system and molecular levels. These are currently being used to unravel new molecular insights about the immune system, and are on the verge of being integrated into clinical trials to enable rational vaccine design and development.

Analysis of cell adhesion during early stages of colon cancer based on an extended multi-valued logic approach

Cell adhesion in the normal colon is typically associated with differentiated cells, whereas in cancerous colon it is associated with advanced tumors. For advanced tumors growing evidence supports the existence of stem-like cells that have originated from transdifferentiation. Because stem cells can also be transformed in their own niche, at the base of the Lieberkühn's crypts, we conjectured that cell adhesion can also be critical in early tumorigenesis. To assess this hypothesis we built an annotated, multi-valued logic model addressing cell adhesion of normal and tumorigenic stem cells in the human colon. The model accounts for (i) events involving intercellular adhesion structures, (ii) interactions involving cytoskeleton-related structures, (iii) compartmental distribution of α/β/γ/δ-catenins, and (iv) variations in critical cell adhesion regulators (e.g., ILK, FAK, IQGAP, SNAIL, Caveolin). We developed a method that can deal with graded multiple inhibitions, something which is not possible with conventional logical approaches. The model comprises 315 species (including 26 genes), interconnected by 269 reactions. Simulations of the model covered six scenarios, which considered two types of colonic cells (stem vs. differentiated cells), under three conditions (normal, stressed and tumor). Each condition results from the combination of 92 inputs. We compared our multi-valued logic approach with the conventional Boolean approach for one specific example and validated the predictions against published data. Our analysis suggests that stem cells in their niche synthesize high levels of cytoplasmatic E-cadherin and CdhEP(Ser684,686,692), even under normal-mitogenic stimulus or tumorigenic conditions. Under these conditions, E-cadherin would be incorporated into the plasmatic membrane, but only as a non-adhesive CdhE_β-catenin_IQGAP complex. Under stress conditions, however, this complex could be displaced, yielding adhesive CdhE_β-catenin((cis/trans)) complexes. In the three scenarios tested with stem cells, desmosomes or tight junctions were not assembled. Other model predictions include expected levels of the nuclear complex β-catenin_TCF4 and the anti-apoptotic protein Survivin for both normal and tumorigenic colonic stem cells.

The in silico macrophage: toward a better understanding of inflammatory disease

Macrophages function as sentinel, cell-regulatory 'hubs' capable of initiating, perpetuating and contributing to the resolution of an inflammatory response, following their activation from a resting state. Highly complex and varied gene expression programs within the macrophage enable such functional diversity. To investigate how programs of gene expression relate to the phenotypic attributes of the macrophage, the development of in silico modeling methods is needed. Such models need to cover multiple scales, from molecular pathways in cell-autonomous immunity and intercellular communication pathways in tissue inflammation to whole organism response pathways in systemic disease. Here, we highlight the potential of in silico macrophage modeling as an amenable and important yet under-exploited tool in aiding in our understanding of the immune inflammatory response. We also discuss how in silico macrophage modeling can help in future therapeutic strategies for modulating both the acute protective effects of inflammation (such as host defense and tissue repair) and the harmful chronic effects (such as autoimmune diseases).