The silicon trypanosome

African trypanosomes have emerged as promising unicellular model organisms for the next generation of systems biology. They offer unique advantages, due to their relative simplicity, the availability of all standard genomics techniques and a long history of quantitative research. Reproducible cultivation methods exist for morphologically and physiologically distinct life-cycle stages. The genome has been sequenced, and microarrays, RNA-interference and high-accuracy metabolomics are available. Furthermore, the availability of extensive kinetic data on all glycolytic enzymes has led to the early development of a complete, experiment-based dynamic model of an important biochemical pathway. Here we describe the achievements of trypanosome systems biology so far and outline the necessary steps towards the ambitious aim of creating a , a comprehensive, experiment-based, multi-scale mathematical model of trypanosome physiology. We expect that, in the long run, the quantitative modelling enabled by the Silicon Trypanosome will play a key role in selecting the most suitable targets for developing new anti-parasite drugs

Local and global modes of drug action in biochemical networks

It becomes increasingly accepted that a shift is needed from the traditional target-based approach of drug development to an integrated perspective of drug action in biochemical systems. We here present an integrative analysis of the interactions between drugs and metabolism based on the concept of drug scope. The drug scope represents the set of metabolic compounds and reactions that are potentially affected by a drug. We constructed and analyzed the scopes of all US approved drugs having metabolic targets. Our analysis shows that the distribution of drug scopes is highly uneven, and that drugs can be classified into several categories based on their scopes. Some of them have small scopes corresponding to localized action, while others have large scopes corresponding to potential large-scale systemic action. These groups are well conserved throughout different topologies of the underlying metabolic network. They can furthermore be associated to specific drug therapeutic properties

Increased entropy of signal transduction in the cancer metastasis phenotype

Studies into the statistical properties of biological networks have led to important biological insights, such as the presence of hubs and hierarchical modularity. There is also a growing interest in studying the statistical properties of networks in the context of cancer genomics. However, relatively little is known as to what network features differ between the cancer and normal cell physiologies, or between different cancer cell phenotypes. Based on the observation that frequent genomic alterations underlie a more aggressive cancer phenotype, we asked if such an effect could be detectable as an increase in the randomness of local gene expression patterns. Using a breast cancer gene expression data set and a model network of protein interactions we derive constrained weighted networks defined by a stochastic information flux matrix reflecting expression correlations between interacting proteins. Based on this stochastic matrix we propose and compute an entropy measure that quantifies the degree of randomness in the local pattern of information flux around single genes. By comparing the local entropies in the non-metastatic versus metastatic breast cancer networks, we here show that breast cancers that metastasize are characterised by a small yet significant increase in the degree of randomness of local expression patterns. We validate this result in three additional breast cancer expression data sets and demonstrate that local entropy better characterises the metastatic phenotype than other non-entropy based measures. We show that increases in entropy can be used to identify genes and signalling pathways implicated in breast cancer metastasis. Further exploration of such integrated cancer expression and protein interaction networks will therefore be a fruitful endeavour.Comment: 5 figures, 2 Supplementary Figures and Table

Computational and Mathematical Modelling of the EGF Receptor System

This chapter gives an overview of computational and mathematical modelling of the EGF receptor system. It begins with a survey of motivations for producing such models, then describes the main approaches that are taken to carrying out such modelling, viz. differential equations and individual-based modelling. Finally, a number of projects that applying modelling and simulation techniques to various aspects of the EGF receptor system are described

The pseudophosphatase MK-STYX interacts with G3BP and decreases stress granule formation

MK-STYX [MAPK (mitogen-activated protein kinase) phospho-serine/threonine/tyrosine-binding protein] is a pseudophosphatase member of the dual-specificity phosphatase subfamily of the PTPs (protein tyrosine phosphatases). MK-STYX is catalytically inactive due to the absence of two amino acids from the signature motif that are essential for phosphatase activity. The nucleophilic cysteine residue and the adjacent histidine residue, which are conserved in all active dual-specificity phosphatases, are replaced by serine and phenylalanine residues respectively in MK-STYX. Mutations to introduce histidine and cysteine residues into the active site of MK-STYX generated an active phosphatase. Using MS, we identified G3BP1 [Ras-GAP (GTPase-activating protein) SH3 (Src homology 3) domain-binding protein-1], a regulator of Ras signalling, as a binding partner of MK-STYX. We observed that G3BP1 bound to native MK-STYX; however, binding to the mutant catalytically active form of MK-STYX was dramatically reduced. G3BP1 is also an RNA-binding protein with endoribonuclease activity that is recruited to ‘stress granules’ after stress stimuli. Stress granules are large subcellular structures that serve as sites of mRNA sorting, in which untranslated mRNAs accumulate. We have shown that expression of MK-STYX inhibited stress granule formation induced either by aresenite or expression of G3BP itself; however, the catalytically active mutant MK-STYX was impaired in its ability to inhibit G3BP-induced stress granule assembly. These results reveal a novel facet of the function of a member of the PTP family, illustrating a role for MK-STYX in regulating the ability of G3BP1 to integrate changes in growth-factor stimulation and environmental stress with the regulation of protein synthesis

Signal duration and the time scale dependence of signal integration in biochemical pathways

Signal duration (e.g. the time scales over which an active signaling intermediate persists) is a key regulator of biological decisions in myriad contexts such as cell growth, proliferation, and developmental lineage commitments. Accompanying differences in signal duration are numerous downstream biological processes that require multiple steps of biochemical regulation. Here, we present an analysis that investigates how simple biochemical motifs that involve multiple stages of regulation can be constructed to differentially process signals that persist at different time scales. We compute the dynamic gain within these networks and resulting power spectra to better understand how biochemical networks can integrate signals at different time scales. We identify topological features of these networks that allow for different frequency dependent signal processing properties. Our studies suggest design principles for why signal duration in connection with multiple steps of downstream regulation is a ubiquitous control motif in biochemical systems.Comment: 27 pages, 4 figure

Linking Proteins to Signaling Pathways for Experiment Design and Evaluation

Biomedical experimental work often focuses on altering the functions of selected proteins. These changes can hit signaling pathways, and can therefore unexpectedly and non-specifically affect cellular processes. We propose PathwayLinker, an online tool that can provide a first estimate of the possible signaling effects of such changes, e.g., drug or microRNA treatments. PathwayLinker minimizes the users' efforts by integrating protein-protein interaction and signaling pathway data from several sources with statistical significance tests and clear visualization. We demonstrate through three case studies that the developed tool can point out unexpected signaling bias in normal laboratory experiments and identify likely novel signaling proteins among the interactors of known drug targets. In our first case study we show that knockdown of the Caenorhabditis elegans gene cdc-25.1 (meant to avoid progeny) may globally affect the signaling system and unexpectedly bias experiments. In the second case study we evaluate the loss-of-function phenotypes of a less known C. elegans gene to predict its function. In the third case study we analyze GJA1, an anti-cancer drug target protein in human, and predict for this protein novel signaling pathway memberships, which may be sources of side effects. Compared to similar services, a major advantage of PathwayLinker is that it drastically reduces the necessary amount of manual literature searches and can be used without a computational background. PathwayLinker is available at http://PathwayLinker.org. Detailed documentation and source code are available at the website

A Hidden Feedback in Signaling Cascades Is Revealed

Cycles involving covalent modification of proteins are key components of the intracellular signaling machinery. Each cycle is comprised of two interconvertable forms of a particular protein. A classic signaling pathway is structured by a chain or cascade of basic cycle units in such a way that the activated protein in one cycle promotes the activation of the next protein in the chain, and so on. Starting from a mechanistic kinetic description and using a careful perturbation analysis, we have derived, to our knowledge for the first time, a consistent approximation of the chain with one variable per cycle. The model we derive is distinct from the one that has been in use in the literature for several years, which is a phenomenological extension of the Goldbeter-Koshland biochemical switch. Even though much has been done regarding the mathematical modeling of these systems, our contribution fills a gap between existing models and, in doing so, we have unveiled critical new properties of this type of signaling cascades. A key feature of our new model is that a negative feedback emerges naturally, exerted between each cycle and its predecessor. Due to this negative feedback, the system displays damped temporal oscillations under constant stimulation and, most important, propagates perturbations both forwards and backwards. This last attribute challenges the widespread notion of unidirectionality in signaling cascades. Concrete examples of applications to MAPK cascades are discussed. All these properties are shared by the complete mechanistic description and our simplified model, but not by previously derived phenomenological models of signaling cascades

Immunohistochemical analysis of changes in signaling pathway activation downstream of growth factor receptors in pancreatic duct cell carcinogenesis

<p>Abstract</p> <p>Background</p> <p>The pathogenesis of pancreatic ductal adenocarcinoma (PDAC) involves multi-stage development of molecular aberrations affecting signaling pathways that regulate cancer growth and progression. This study was performed to gain a better understanding of the abnormal signaling that occurs in PDAC compared with normal duct epithelia.</p> <p>Methods</p> <p>We performed immunohistochemistry on a tissue microarray of 26 PDAC, 13 normal appearing adjacent pancreatic ductal epithelia, and 12 normal non-PDAC ducts. We compared the levels of 18 signaling proteins including growth factor receptors, tumor suppressors and 13 of their putative downstream phosphorylated (p-) signal transducers in PDAC to those in normal ductal epithelia.</p> <p>Results</p> <p>The overall profiles of signaling protein expression levels, activation states and sub-cellular distribution in PDAC cells were distinguishable from non-neoplastic ductal epithelia. The ERK pathway activation was correlated with high levels of ^S2448p-mTOR (100%, p = 0.05), ^T389p-S6K (100%, p = 0.02 and ^S235/236p-S6 (86%, p = 0.005). Additionally, ^T389p-S6K correlated with ^S727p-STAT3 (86%, p = 0.005). Advanced tumors with lymph node metastasis were characterized by high levels of ^S276p-NFκB (100%, p = 0.05) and ^S9p-GSK3β (100%, p = 0.05). High levels of PKBβ/AKT2, EGFR, as well as nuclear ^T202/Y204p-ERK and ^T180/Y182p-p38 were observed in normal ducts adjacent to PDAC compared with non-cancerous pancreas.</p> <p>Conclusion</p> <p>Multiple signaling proteins are activated in pancreatic duct cell carcinogenesis including those associated with the ERK, PKB/AKT, mTOR and STAT3 pathways. The ERK pathway activation appears also increased in duct epithelia adjacent to carcinoma, suggesting tumor micro-environmental effects.</p