Dynamical modeling of the cholesterol regulatory pathway with Boolean networks

Mathematical Models for Cholesterol Metabolism and Transport

Cholesterol is an essential component of eukaryotic cellular membranes. It is also an important precursor for making other molecules needed by the body. Cholesterol homeostasis plays an essential role in human health. Having high cholesterol can increase the chances of getting heart disease. As a result of the risks associated with high cholesterol, it is imperative that studies are conducted to determine the best course of action to reduce whole body cholesterol levels. Mathematical models can provide direction on this. By examining existing models, the suitable reactions or processes for drug targeting to lower whole-body cholesterol can be determined. This paper examines existing models in the literature that, in total, cover most of the processes involving cholesterol metabolism and transport, including: the absorption of cholesterol in the intestine; the cholesterol biosynthesis in the liver; the storage and transport of cholesterol between the intestine, the liver, blood vessels, and peripheral cells. The findings presented in these models will be discussed for potential combination to form a comprehensive model of cholesterol within the entire body, which is then taken as an in-silico patient for identifying drug targets, screening drugs, and designing intervention strategies to regulate cholesterol levels in the human body.

Modeling cholesterol metabolism and atherosclerosis

Atherosclerotic cardiovascular disease (ASCVD) is the leading cause of morbidity and mortality among Western populations. Many risk factors have been identified for ASCVD; however, elevated low-density lipoprotein cholesterol (LDL-C) remains the gold standard. Cholesterol metabolism at the cellular and whole-body level is maintained by an array of interacting components. These regulatory mechanisms have complex behavior. Likewise, the mechanisms which underpin atherogenesis are nontrivial and multifaceted. To help overcome the challenge of investigating these processes mathematical modeling, which is a core constituent of the systems biology paradigm has played a pivotal role in deciphering their dynamics. In so doing models have revealed new insights about the key drivers of ASCVD. The aim of this review is fourfold; to provide an overview of cholesterol metabolism and atherosclerosis, to briefly introduce mathematical approaches used in this field, to critically discuss models of cholesterol metabolism and atherosclerosis, and to highlight areas where mathematical modeling could help to investigate in the future. This article is categorized under: Cardiovascular Diseases > Computational Models.

Insights on TAM Formation from a Boolean Model of Macrophage Polarization Based on In Vitro Studies

Simple Summary

The recent success of immunotherapy treatments against cancer relies on helping our own body’s defenses in the fight against tumours, namely reinvigorating the cancer killing action of T cells. Unfortunately, in a large proportion of patients these therapies are ineffective, in part due to the presence of other immune cells, macrophages, which are mis-educated by the cancer cells into promoting tumour growth. Here we start from an existing model of macrophage polarization and extend it to the specific conditions encountered inside a tumour by adding signals, receptors, transcription factors and cytokines that are known to be the key components in establishing the cancer cell-macrophage interaction. Then we use a mathematical Boolean model applied to a gene regulatory network of this biological process to simulate its temporal behaviour and explore scenarios that have not been experimentally tested so far. Additionally, the KO and overexpression simulations successfully reproduce the known experimental results while predicting the potential role of regulators (such as STAT1 and EGF) in preventing the formation of pro-tumoural macrophages, which can be tested experimentally.

Abstract

The tumour microenvironment is the surrounding of a tumour, including blood vessels, fibroblasts, signaling molecules, the extracellular matrix and immune cells, especially neutrophils and monocyte-derived macrophages. In a tumour setting, macrophages encompass a spectrum between a tumour-suppressive (M1) or tumour-promoting (M2) state. The biology of macrophages found in tumours (Tumour Associated Macrophages) remains unclear, but understanding their impact on tumour progression is highly important. In this paper, we perform a comprehensive analysis of a macrophage polarization network, following two lines of enquiry: (i) we reconstruct the macrophage polarization network based on literature, extending it to include important stimuli in a tumour setting, and (ii) we build a dynamical model able to reproduce macrophage polarization in the presence of different stimuli, including the contact with cancer cells. Our simulations recapitulate the documented macrophage phenotypes and their dependencies on specific receptors and transcription factors, while also unravelling the formation of a special type of tumour associated macrophages in an in vitro model of chronic lymphocytic leukaemia. This model constitutes the first step towards elucidating the cross-talk between immune and cancer cells inside tumours, with the ultimate goal of identifying new therapeutic targets that could control the formation of tumour associated macrophages in patients.

Cholesterol Homeostasis: An In Silico Investigation into How Aging Disrupts Its Key Hepatic Regulatory Mechanisms

The dysregulation of intracellular cholesterol homeostasis is associated with several age-related diseases, most notably cardiovascular disease (CVD). Research in this area has benefitted from using computational modelling to study the inherent complexity associated with the regulation of this system. In addition to facilitating hypothesis exploration, the utility of modelling lies in its ability to represent an array of rate limiting enzymatic reactions, together with multiple feedback loops, which collectively define the dynamics of cholesterol homeostasis. However, to date no model has specifically investigated the effects aging has on this system. This work addresses this shortcoming by explicitly focusing on the impact of aging on hepatic intracellular cholesterol homeostasis. The model was used to investigate the experimental findings that reactive oxygen species induce the total activation of 3-hydroxy-3-methylglutaryl-coenzyme A (HMG-CoA) reductase (HMGCR). Moreover, the model explored the impact of an age-related decrease in hepatic acetyl-CoA acetyltransferase 2 (ACAT2). The model suggested that an increase in the activity of HMGCR does not have as significant an impact on cholesterol homeostasis as a decrease in hepatic ACAT2 activity. According to the model, a decrease in the activity of hepatic ACAT2 raises free cholesterol (FC) and decreases low-density lipoprotein cholesterol (LDL-C) levels. Increased acetyl CoA synthesis resulted in a reduction in the number of hepatic low-density lipoprotein receptors, and increased LDL-C, FC, and cholesterol esters. The rise in LDL-C was restricted by elevated hepatic FC accumulation. Taken together these findings have important implications for healthspan. This is because emerging clinical data suggest hepatic FC accumulation is relevant to the pathogenesis of non-alcoholic fatty liver disease (NAFLD), which is associated with an increased risk of CVD. These pathophysiological changes could, in part, help to explain the phenomenon of increased mortality associated with low levels of LDL-C which have been observed in certain studies involving the oldest old (≥85 years).

Identification of Biologically Essential Nodes via Determinative Power in Logical Models of Cellular Processes

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.

A Boolean network of the crosstalk between IGF and Wnt signaling in aging satellite cells

Aging is a complex biological process, which determines the life span of an organism. Insulin-like growth factor (IGF) and Wnt signaling pathways govern the process of aging. Both pathways share common downstream targets that allow competitive crosstalk between these branches. Of note, a shift from IGF to Wnt signaling has been observed during aging of satellite cells. Biological regulatory networks necessary to recreate aging have not yet been discovered. Here, we established a mathematical in silico model that robustly recapitulates the crosstalk between IGF and Wnt signaling. Strikingly, it predicts critical nodes following a shift from IGF to Wnt signaling. These findings indicate that this shift might cause age-related diseases.

A new discrete dynamic model of ABA-induced stomatal closure predicts key feedback loops

Stomata, microscopic pores in leaf surfaces through which water loss and carbon dioxide uptake occur, are closed in response to drought by the phytohormone abscisic acid (ABA). This process is vital for drought tolerance and has been the topic of extensive experimental investigation in the last decades. Although a core signaling chain has been elucidated consisting of ABA binding to receptors, which alleviates negative regulation by protein phosphatases 2C (PP2Cs) of the protein kinase OPEN STOMATA 1 (OST1) and ultimately results in activation of anion channels, osmotic water loss, and stomatal closure, over 70 additional components have been identified, yet their relationships with each other and the core components are poorly elucidated. We integrated and processed hundreds of disparate observations regarding ABA signal transduction responses underlying stomatal closure into a network of 84 nodes and 156 edges and, as a result, established those relationships, including identification of a 36-node, strongly connected (feedback-rich) component as well as its in- and out-components. The network's domination by a feedback-rich component may reflect a general feature of rapid signaling events. We developed a discrete dynamic model of this network and elucidated the effects of ABA plus knockout or constitutive activity of 79 nodes on both the outcome of the system (closure) and the status of all internal nodes. The model, with more than 1024 system states, is far from fully determined by the available data, yet model results agree with existing experiments in 82 cases and disagree in only 17 cases, a validation rate of 75%. Our results reveal nodes that could be engineered to impact stomatal closure in a controlled fashion and also provide over 140 novel predictions for which experimental data are currently lacking. Noting the paucity of wet-bench data regarding combinatorial effects of ABA and internal node activation, we experimentally confirmed several predictions of the model with regard to reactive oxygen species, cytosolic Ca2+ (Ca2+c), and heterotrimeric G-protein signaling. We analyzed dynamics-determining positive and negative feedback loops, thereby elucidating the attractor (dynamic behavior) repertoire of the system and the groups of nodes that determine each attractor. Based on this analysis, we predict the likely presence of a previously unrecognized feedback mechanism dependent on Ca2+c. This mechanism would provide model agreement with 10 additional experimental observations, for a validation rate of 85%. Our research underscores the importance of feedback regulation in generating robust and adaptable biological responses. The high validation rate of our model illustrates the advantages of discrete dynamic modeling for complex, nonlinear systems common in biology.

Cytochrome P450 metabolism of the post-lanosterol intermediates explains enigmas of cholesterol synthesis

Cholesterol synthesis is among the oldest metabolic pathways, consisting of the Bloch and Kandutch-Russell branches. Following lanosterol, sterols of both branches are proposed to be dedicated to cholesterol. We challenge this dogma by mathematical modeling and with experimental evidence. It was not possible to explain the sterol profile of testis in cAMP responsive element modulator tau (Crem τ) knockout mice with mathematical models based on textbook pathways of cholesterol synthesis. Our model differs in the inclusion of virtual sterol metabolizing enzymes branching from the pathway. We tested the hypothesis that enzymes from the cytochrome P450 (CYP) superfamily can participate in the catalysis of non-classical reactions. We show that CYP enzymes can metabolize multiple sterols in vitro, establishing novel branching points of cholesterol synthesis. In conclusion, sterols of cholesterol synthesis can be oxidized further to metabolites not dedicated to production of cholesterol. Additionally, CYP7A1, CYP11A1, CYP27A1, and CYP46A1 are parts of a broader cholesterol synthesis network.

Cholesterol metabolism: A review of how ageing disrupts the biological mechanisms responsible for its regulation

Cholesterol plays a vital role in the human body as a precursor of steroid hormones and bile acids, in addition to providing structure to cell membranes. Whole body cholesterol metabolism is maintained by a highly coordinated balancing act between cholesterol ingestion, synthesis, absorption, and excretion. The aim of this review is to discuss how ageing interacts with these processes. Firstly, we will present an overview of cholesterol metabolism. Following this, we discuss how the biological mechanisms which underpin cholesterol metabolism are effected by ageing. Included in this discussion are lipoprotein dynamics, cholesterol absorption/synthesis and the enterohepatic circulation/synthesis of bile acids. Moreover, we discuss the role of oxidative stress in the pathological progression of atherosclerosis and also discuss how cholesterol biosynthesis is effected by both the mammalian target of rapamycin and sirtuin pathways. Next, we examine how diet and alterations to the gut microbiome can be used to mitigate the impact ageing has on cholesterol metabolism. We conclude by discussing how mathematical models of cholesterol metabolism can be used to identify therapeutic interventions.

An algebra-based method for inferring gene regulatory networks

BACKGROUND

The inference of gene regulatory networks (GRNs) from experimental observations is at the heart of systems biology. This includes the inference of both the network topology and its dynamics. While there are many algorithms available to infer the network topology from experimental data, less emphasis has been placed on methods that infer network dynamics. Furthermore, since the network inference problem is typically underdetermined, it is essential to have the option of incorporating into the inference process, prior knowledge about the network, along with an effective description of the search space of dynamic models. Finally, it is also important to have an understanding of how a given inference method is affected by experimental and other noise in the data used.

RESULTS

This paper contains a novel inference algorithm using the algebraic framework of Boolean polynomial dynamical systems (BPDS), meeting all these requirements. The algorithm takes as input time series data, including those from network perturbations, such as knock-out mutant strains and RNAi experiments. It allows for the incorporation of prior biological knowledge while being robust to significant levels of noise in the data used for inference. It uses an evolutionary algorithm for local optimization with an encoding of the mathematical models as BPDS. The BPDS framework allows an effective representation of the search space for algebraic dynamic models that improves computational performance. The algorithm is validated with both simulated and experimental microarray expression profile data. Robustness to noise is tested using a published mathematical model of the segment polarity gene network in Drosophila melanogaster. Benchmarking of the algorithm is done by comparison with a spectrum of state-of-the-art network inference methods on data from the synthetic IRMA network to demonstrate that our method has good precision and recall for the network reconstruction task, while also predicting several of the dynamic patterns present in the network.

CONCLUSIONS

Boolean polynomial dynamical systems provide a powerful modeling framework for the reverse engineering of gene regulatory networks, that enables a rich mathematical structure on the model search space. A C++ implementation of the method, distributed under LPGL license, is available, together with the source code, at http://www.paola-vera-licona.net/Software/EARevEng/REACT.html.

Mathematical models of mitochondrial aging and dynamics

Research on the role of mitochondria in aging and disease is rapidly growing. Furthermore, in recent years, it also became clear that mitochondria are dynamic structures undergoing constant and rapid cycles of fusion and fission. The involvement of mitochondria in multiple complex processes makes them a prime target for mathematical and computational modeling. This review consists of two parts. In the first (Section 2), we provide a detailed introduction to the underlying concepts of mathematical modeling to help the reader who is not so familiar with these techniques to judge the requirements and results that can be obtained through modeling. In the second part (Section 3), we review existing mathematical and computational models that investigate mitochondrial dynamics and the role of mitochondria for the aging process.

Learning contextual gene set interaction networks of cancer with condition specificity

Background

Identifying similarities and differences in the molecular constitutions of various types of cancer is one of the key challenges in cancer research. The appearances of a cancer depend on complex molecular interactions, including gene regulatory networks and gene-environment interactions. This complexity makes it challenging to decipher the molecular origin of the cancer. In recent years, many studies reported methods to uncover heterogeneous depictions of complex cancers, which are often categorized into different subtypes. The challenge is to identify diverse molecular contexts within a cancer, to relate them to different subtypes, and to learn underlying molecular interactions specific to molecular contexts so that we can recommend context-specific treatment to patients.

Results

In this study, we describe a novel method to discern molecular interactions specific to certain molecular contexts. Unlike conventional approaches to build modular networks of individual genes, our focus is to identify cancer-generic and subtype-specific interactions between contextual gene sets, of which each gene set share coherent transcriptional patterns across a subset of samples, termed contextual gene set. We then apply a novel formulation for quantitating the effect of the samples from each subtype on the calculated strength of interactions observed. Two cancer data sets were analyzed to support the validity of condition-specificity of identified interactions. When compared to an existing approach, the proposed method was much more sensitive in identifying condition-specific interactions even in heterogeneous data set. The results also revealed that network components specific to different types of cancer are related to different biological functions than cancer-generic network components. We found not only the results that are consistent with previous studies, but also new hypotheses on the biological mechanisms specific to certain cancer types that warrant further investigations.

Conclusions

The analysis on the contextual gene sets and characterization of networks of interaction composed of these sets discovered distinct functional differences underlying various types of cancer. The results show that our method successfully reveals many subtype-specific regions in the identified maps of biological contexts, which well represent biological functions that can be connected to specific subtypes.

A whole-body mathematical model of cholesterol metabolism and its age-associated dysregulation

BACKGROUND

Global demographic changes have stimulated marked interest in the process of aging. There has been, and will continue to be, an unrelenting rise in the number of the oldest old ( >85 years of age). Together with an ageing population there comes an increase in the prevalence of age related disease. Of the diseases of ageing, cardiovascular disease (CVD) has by far the highest prevalence. It is regarded that a finely tuned lipid profile may help to prevent CVD as there is a long established relationship between alterations to lipid metabolism and CVD risk. In fact elevated plasma cholesterol, particularly Low Density Lipoprotein Cholesterol (LDL-C) has consistently stood out as a risk factor for having a cardiovascular event. Moreover it is widely acknowledged that LDL-C may rise with age in both sexes in a wide variety of groups. The aim of this work was to use a whole-body mathematical model to investigate why LDL-C rises with age, and to test the hypothesis that mechanistic changes to cholesterol absorption and LDL-C removal from the plasma are responsible for the rise. The whole-body mechanistic nature of the model differs from previous models of cholesterol metabolism which have either focused on intracellular cholesterol homeostasis or have concentrated on an isolated area of lipoprotein dynamics. The model integrates both current and previously published data relating to molecular biology, physiology, ageing and nutrition in an integrated fashion.

RESULTS

The model was used to test the hypothesis that alterations to the rate of cholesterol absorption and changes to the rate of removal of LDL-C from the plasma are integral to understanding why LDL-C rises with age. The model demonstrates that increasing the rate of intestinal cholesterol absorption from 50% to 80% by age 65 years can result in an increase of LDL-C by as much as 34 mg/dL in a hypothetical male subject. The model also shows that decreasing the rate of hepatic clearance of LDL-C gradually to 50% by age 65 years can result in an increase of LDL-C by as much as 116 mg/dL.

CONCLUSIONS

Our model clearly demonstrates that of the two putative mechanisms that have been implicated in the dysregulation of cholesterol metabolism with age, alterations to the removal rate of plasma LDL-C has the most significant impact on cholesterol metabolism and small changes to the number of hepatic LDL receptors can result in a significant rise in LDL-C. This first whole-body systems based model of cholesterol balance could potentially be used as a tool to further improve our understanding of whole-body cholesterol metabolism and its dysregulation with age. Furthermore, given further fine tuning the model may help to investigate potential dietary and lifestyle regimes that have the potential to mitigate the effects aging has on cholesterol metabolism.

A physiologically based in silico kinetic model predicting plasma cholesterol concentrations in humans

Increased plasma cholesterol concentration is associated with increased risk of cardiovascular disease. This study describes the development, validation, and analysis of a physiologically based kinetic (PBK) model for the prediction of plasma cholesterol concentrations in humans. This model was directly adapted from a PBK model for mice by incorporation of the reaction catalyzed by cholesterol ester transfer protein and contained 21 biochemical reactions and eight different cholesterol pools. The model was calibrated using published data for humans and validated by comparing model predictions on plasma cholesterol levels of subjects with 10 different genetic mutations (including familial hypercholesterolemia and Smith-Lemli-Opitz syndrome) with experimental data. Average model predictions on total cholesterol were accurate within 36% of the experimental data, which was within the experimental margin. Sensitivity analysis of the model indicated that the HDL cholesterol (HDL-C) concentration was mainly dependent on hepatic transport of cholesterol to HDL, cholesterol ester transfer from HDL to non-HDL, and hepatic uptake of cholesterol from non-HDL-C. Thus, the presented PBK model is a valid tool to predict the effect of genetic mutations on cholesterol concentrations, opening the way for future studies on the effect of different drugs on cholesterol levels in various subpopulations in silico.

Stochastic Boolean networks: an efficient approach to modeling gene regulatory networks

BACKGROUND

Various computational models have been of interest due to their use in the modelling of gene regulatory networks (GRNs). As a logical model, probabilistic Boolean networks (PBNs) consider molecular and genetic noise, so the study of PBNs provides significant insights into the understanding of the dynamics of GRNs. This will ultimately lead to advances in developing therapeutic methods that intervene in the process of disease development and progression. The applications of PBNs, however, are hindered by the complexities involved in the computation of the state transition matrix and the steady-state distribution of a PBN. For a PBN with n genes and N Boolean networks, the complexity to compute the state transition matrix is O(nN22n) or O(nN2n) for a sparse matrix.

RESULTS

This paper presents a novel implementation of PBNs based on the notions of stochastic logic and stochastic computation. This stochastic implementation of a PBN is referred to as a stochastic Boolean network (SBN). An SBN provides an accurate and efficient simulation of a PBN without and with random gene perturbation. The state transition matrix is computed in an SBN with a complexity of O(nL2n), where L is a factor related to the stochastic sequence length. Since the minimum sequence length required for obtaining an evaluation accuracy approximately increases in a polynomial order with the number of genes, n, and the number of Boolean networks, N, usually increases exponentially with n, L is typically smaller than N, especially in a network with a large number of genes. Hence, the computational efficiency of an SBN is primarily limited by the number of genes, but not directly by the total possible number of Boolean networks. Furthermore, a time-frame expanded SBN enables an efficient analysis of the steady-state distribution of a PBN. These findings are supported by the simulation results of a simplified p53 network, several randomly generated networks and a network inferred from a T cell immune response dataset. An SBN can also implement the function of an asynchronous PBN and is potentially useful in a hybrid approach in combination with a continuous or single-molecule level stochastic model.

CONCLUSIONS

Stochastic Boolean networks (SBNs) are proposed as an efficient approach to modelling gene regulatory networks (GRNs). The SBN approach is able to recover biologically-proven regulatory behaviours, such as the oscillatory dynamics of the p53-Mdm2 network and the dynamic attractors in a T cell immune response network. The proposed approach can further predict the network dynamics when the genes are under perturbation, thus providing biologically meaningful insights for a better understanding of the dynamics of GRNs. The algorithms and methods described in this paper have been implemented in Matlab packages, which are attached as Additional files.

Dynamical systems approach to endothelial heterogeneity

Endothelial cells display remarkable phenotypic heterogeneity. An important goal is to elucidate the scope and mechanisms of endothelial heterogeneity and to use this information to develop vascular bed-specific therapies. We reexamine our current understanding of the molecular basis of endothelial heterogeneity. We introduce multistability as a new explanatory framework in vascular biology. We draw on the field of nonlinear dynamics to propose a dynamical systems framework for modeling multistability and its derivative properties, including robustness, memory, and plasticity. Our perspective allows for both a conceptual and quantitative description of system-level features of endothelial regulation.

Discrete dynamic modeling of signal transduction networks

Newly available experimental data characterizing different processes involved in signaling pathways have provided the opportunity for network analysis and modeling of these interacting pathways. Current approaches in studying the dynamics of signaling networks fall into two major groups, namely, continuous and discrete models. The lack of kinetic information for biochemical interactions has limited the wide applicability of continuous models. To address this issue, discrete dynamic models, based on a qualitative description of a system's variables, have been applied for the analysis of biological systems with many unknown parameters. The purpose of this chapter is to give a detailed description of Boolean modeling, the simplest type of discrete dynamic modeling, and the ways in which it can be applied to analyze the dynamics of signaling networks. This is followed by practical examples of a Boolean dynamic framework applied to the modeling of the abscisic acid signal transduction network in plants as well as the T-cell survival signaling network in humans.

Dynamical and structural analysis of a T cell survival network identifies novel candidate therapeutic targets for large granular lymphocyte leukemia

The blood cancer T cell large granular lymphocyte (T-LGL) leukemia is a chronic disease characterized by a clonal proliferation of cytotoxic T cells. As no curative therapy is yet known for this disease, identification of potential therapeutic targets is of immense importance. In this paper, we perform a comprehensive dynamical and structural analysis of a network model of this disease. By employing a network reduction technique, we identify the stationary states (fixed points) of the system, representing normal and diseased (T-LGL) behavior, and analyze their precursor states (basins of attraction) using an asynchronous Boolean dynamic framework. This analysis identifies the T-LGL states of 54 components of the network, out of which 36 (67%) are corroborated by previous experimental evidence and the rest are novel predictions. We further test and validate one of these newly identified states experimentally. Specifically, we verify the prediction that the node SMAD is over-active in leukemic T-LGL by demonstrating the predominant phosphorylation of the SMAD family members Smad2 and Smad3. Our systematic perturbation analysis using dynamical and structural methods leads to the identification of 19 potential therapeutic targets, 68% of which are corroborated by experimental evidence. The novel therapeutic targets provide valuable guidance for wet-bench experiments. In addition, we successfully identify two new candidates for engineering long-lived T cells necessary for the delivery of virus and cancer vaccines. Overall, this study provides a bird's-eye-view of the avenues available for identification of therapeutic targets for similar diseases through perturbation of the underlying signal transduction network.

T-LGL leukemia is a blood cancer characterized by an abnormal increase in the abundance of a type of white blood cell called T cell. Since there is no known curative therapy for this disease, identification of potential therapeutic targets is of utmost importance. Experimental identification of manipulations capable of reversing the disease condition is usually a long, arduous process. Mathematical modeling can aid this process by identifying potential therapeutic interventions. In this work, we carry out a systematic analysis of a network model of T cell survival in T-LGL leukemia to get a deeper insight into the unknown facets of the disease. We identify the T-LGL status of 54 components of the system, out of which 36 (67%) are corroborated by previous experimental evidence and the rest are novel predictions, one of which we validate by follow-up experiments. By deciphering the structure and dynamics of the underlying network, we identify component perturbations that lead to programmed cell death, thereby suggesting several novel candidate therapeutic targets for future experiments.

Exploring phospholipase C-coupled Ca(2+) signalling networks using Boolean modelling

In this study, the authors explored the utility of a descriptive and predictive bionetwork model for phospholipase C-coupled calcium signalling pathways, built with non-kinetic experimental information. Boolean models generated from these data yield oscillatory activity patterns for both the endoplasmic reticulum resident inositol-1,4,5-trisphosphate receptor (IP(3)R) and the plasma-membrane resident canonical transient receptor potential channel 3 (TRPC3). These results are specific as randomisation of the Boolean operators ablates oscillatory pattern formation. Furthermore, knock-out simulations of the IP(3)R, TRPC3 and multiple other proteins recapitulate experimentally derived results. The potential of this approach can be observed by its ability to predict previously undescribed cellular phenotypes using in vitro experimental data. Indeed, our cellular analysis of the developmental and calcium-regulatory protein, DANGER1a, confirms the counter-intuitive predictions from our Boolean models in two highly relevant cellular models. Based on these results, the authors theorise that with sufficient legacy knowledge and/or computational biology predictions, Boolean networks can provide a robust method for predictive modelling of any biological system. [Includes supplementary material].

The differences of sarcopenia-related phenotypes: effects of gender and population

Sarcopenia is a serious condition especially in the elderly population mainly characterized by the loss of skeletal muscle mass and strength with aging. Extremity skeletal muscle mass index (EMMI) (sum of skeletal muscle mass in arms and legs/height2) is gaining popularity in sarcopenia definition (less than two standard deviations below the mean of a young adult reference group), but little is known about the gender- and population-specific differences of EMMI. This study aimed at investigating the differences of EMMI, arm muscle mass index (AMMI), and leg muscle mass index (LMMI) between gender groups and populations (Chinese vs. Caucasians). The participants included 1,809 Chinese and 362 Caucasians with normal weight aged from 19 to 45 years old. Extremity muscle mass, arm muscle mass, and leg muscle mass were measured by using dual energy x-ray absorptiometry. Independent sample t tests were used to analyze the differences in muscle mass indexes between the studied groups. All the study parameters including EMMIs, AMMIs, and LMMIs were significantly higher (P ≤ 0.0003) in the Caucasian group than in the Chinese group and also higher in the male group than in the female group, and these significant differences (P ≤ 0.0005) remained after adjusting for age by simple regressions. The detected differences of muscle mass indexes between different gender and ethnic groups may provide important implications in their different risk of future sarcopenia.

Adapted Boolean network models for extracellular matrix formation

BACKGROUND

Due to the rapid data accumulation on pathogenesis and progression of chronic inflammation, there is an increasing demand for approaches to analyse the underlying regulatory networks. For example, rheumatoid arthritis (RA) is a chronic inflammatory disease, characterised by joint destruction and perpetuated by activated synovial fibroblasts (SFB). These abnormally express and/or secrete pro-inflammatory cytokines, collagens causing joint fibrosis, or tissue-degrading enzymes resulting in destruction of the extra-cellular matrix (ECM).We applied three methods to analyse ECM regulation: data discretisation to filter out noise and to reduce complexity, Boolean network construction to implement logic relationships, and formal concept analysis (FCA) for the formation of minimal, but complete rule sets from the data.

RESULTS

First, we extracted literature information to develop an interaction network containing 18 genes representing ECM formation and destruction. Subsequently, we constructed an asynchronous Boolean network with biologically plausible time intervals for mRNA and protein production, secretion, and inactivation. Experimental gene expression data was obtained from SFB stimulated by TGFbeta1 or by TNFalpha and discretised thereafter. The Boolean functions of the initial network were improved iteratively by the comparison of the simulation runs to the experimental data and by exploitation of expert knowledge. This resulted in adapted networks for both cytokine stimulation conditions. The simulations were further analysed by the attribute exploration algorithm of FCA, integrating the observed time series in a fine-tuned and automated manner. The resulting temporal rules yielded new contributions to controversially discussed aspects of fibroblast biology (e.g., considerable expression of TNF and MMP9 by fibroblasts stimulation) and corroborated previously known facts (e.g., co-expression of collagens and MMPs after TNFalpha stimulation), but also revealed some discrepancies to literature knowledge (e.g., MMP1 expression in the absence of FOS).

CONCLUSION

The newly developed method successfully and iteratively integrated expert knowledge at different steps, resulting in a promising solution for the in-depth understanding of regulatory pathways in disease dynamics. The knowledge base containing all the temporal rules may be queried to predict the functional consequences of observed or hypothetical gene expression disturbances. Furthermore, new hypotheses about gene relations were derived which await further experimental validation.

Discrete dynamic modeling of cellular signaling networks

Understanding signal transduction in cellular systems is a central issue in systems biology. Numerous experiments from different laboratories generate an abundance of individual components and causal interactions mediating environmental and developmental signals. However, for many signal transduction systems there is insufficient information on the overall structure and the molecular mechanisms involved in the signaling network. Moreover, lack of kinetic and temporal information makes it difficult to construct quantitative models of signal transduction pathways. Discrete dynamic modeling, combined with network analysis, provides an effective way to integrate fragmentary knowledge of regulatory interactions into a predictive mathematical model which is able to describe the time evolution of the system without the requirement for kinetic parameters. This chapter introduces the fundamental concepts of discrete dynamic modeling, particularly focusing on Boolean dynamic models. We describe this method step-by-step in the context of cellular signaling networks. Several variants of Boolean dynamic models including threshold Boolean networks and piecewise linear systems are also covered, followed by two examples of successful application of discrete dynamic modeling in cell biology.