The statistical mechanics of complex signaling networks : nerve growth factor signaling

It is becoming increasingly appreciated that the signal transduction systems used by eukaryotic cells to achieve a variety of essential responses represent highly complex networks rather than simple linear pathways. While significant effort is being made to experimentally measure the rate constants for individual steps in these signaling networks, many of the parameters required to describe the behavior of these systems remain unknown, or at best, estimates. With these goals and caveats in mind, we use methods of statistical mechanics to extract useful predictions for complex cellular signaling networks. To establish the usefulness of our approach, we have applied our methods towards modeling the nerve growth factor (NGF)-induced differentiation of neuronal cells. Using our approach, we are able to extract predictions that are highly specific and accurate, thereby enabling us to predict the influence of specific signaling modules in determining the integrated cellular response to the two growth factors. We show that extracting biologically relevant predictions from complex signaling models appears to be possible even in the absence of measurements of all the individual rate constants. Our methods also raise some interesting insights into the design and possible evolution of cellular systems, highlighting an inherent property of these systems wherein particular ''soft'' combinations of parameters can be varied over wide ranges without impacting the final output and demonstrating that a few ''stiff'' parameter combinations center around the paramount regulatory steps of the network. We refer to this property -- which is distinct from robustness -- as ''sloppiness.''Comment: 24 pages, 10 EPS figures, 1 GIF (makes 5 multi-panel figs + caption for GIF), IOP style; supp. info/figs. included as brown_supp.pd

How to understand the cell by breaking it: network analysis of gene perturbation screens

Modern high-throughput gene perturbation screens are key technologies at the forefront of genetic research. Combined with rich phenotypic descriptors they enable researchers to observe detailed cellular reactions to experimental perturbations on a genome-wide scale. This review surveys the current state-of-the-art in analyzing perturbation screens from a network point of view. We describe approaches to make the step from the parts list to the wiring diagram by using phenotypes for network inference and integrating them with complementary data sources. The first part of the review describes methods to analyze one- or low-dimensional phenotypes like viability or reporter activity; the second part concentrates on high-dimensional phenotypes showing global changes in cell morphology, transcriptome or proteome.Comment: Review based on ISMB 2009 tutorial; after two rounds of revisio

A predictive computational model reveals that GIV/girdin serves as a tunable valve for EGFR-stimulated cyclic AMP signals.

Cellular levels of the versatile second messenger cyclic (c)AMP are regulated by the antagonistic actions of the canonical G protein → adenylyl cyclase pathway that is initiated by G-protein-coupled receptors (GPCRs) and attenuated by phosphodiesterases (PDEs). Dysregulated cAMP signaling drives many diseases; for example, its low levels facilitate numerous sinister properties of cancer cells. Recently, an alternative paradigm for cAMP signaling has emerged in which growth factor-receptor tyrosine kinases (RTKs; e.g., EGFR) access and modulate G proteins via a cytosolic guanine-nucleotide exchange modulator (GEM), GIV/girdin; dysregulation of this pathway is frequently encountered in cancers. In this study, we present a network-based compartmental model for the paradigm of GEM-facilitated cross-talk between RTKs and G proteins and how that impacts cellular cAMP. Our model predicts that cross-talk between GIV, Gαs, and Gαi proteins dampens ligand-stimulated cAMP dynamics. This prediction was experimentally verified by measuring cAMP levels in cells under different conditions. We further predict that the direct proportionality of cAMP concentration as a function of receptor number and the inverse proportionality of cAMP concentration as a function of PDE concentration are both altered by GIV levels. Taking these results together, our model reveals that GIV acts as a tunable control valve that regulates cAMP flux after growth factor stimulation. For a given stimulus, when GIV levels are high, cAMP levels are low, and vice versa. In doing so, GIV modulates cAMP via mechanisms distinct from the two most often targeted classes of cAMP modulators, GPCRs and PDEs

Alteration of gene expression in mammary gland tissue of dairy cows in response to dietary unsaturated fatty acids

The aim of this study was to determine the effects of supplementing unprotected dietary unsaturated fatty acids (UFAs) from different plant oils on gene expression in the mammary gland of grazing dairy cows. A total of 28 Holstein–Friesian dairy cows in mid-lactation were blocked according to parity, days in milk, milk yield and fat percentage. The cows were then randomly assigned to four UFA sources based on rapeseed, soybean, linseed or a mixture of the three oils for 23 days, after which, all 28 cows were switched to a control diet for an additional 28 days. On the last day of both periods, mammary gland biopsies were taken to study genome-wide differences in gene expression on Affymetrix GeneChip® Bovine Genome Arrays (no. 900493) by ServiceXS (Leiden, The Netherlands). Supplementation with UFAs resulted in increased milk yield but decreased milk fat and protein percentages. Furthermore, the proportion of de novo fatty acids (FAs) in the milk was reduced, whereas that of long-chain FAs increased. Applying a statistical cut-off of false discovery rate of q-value

Engineering simulations for cancer systems biology

Computer simulation can be used to inform in vivo and in vitro experimentation, enabling rapid, low-cost hypothesis generation and directing experimental design in order to test those hypotheses. In this way, in silico models become a scientific instrument for investigation, and so should be developed to high standards, be carefully calibrated and their findings presented in such that they may be reproduced. Here, we outline a framework that supports developing simulations as scientific instruments, and we select cancer systems biology as an exemplar domain, with a particular focus on cellular signalling models. We consider the challenges of lack of data, incomplete knowledge and modelling in the context of a rapidly changing knowledge base. Our framework comprises a process to clearly separate scientific and engineering concerns in model and simulation development, and an argumentation approach to documenting models for rigorous way of recording assumptions and knowledge gaps. We propose interactive, dynamic visualisation tools to enable the biological community to interact with cellular signalling models directly for experimental design. There is a mismatch in scale between these cellular models and tissue structures that are affected by tumours, and bridging this gap requires substantial computational resource. We present concurrent programming as a technology to link scales without losing important details through model simplification. We discuss the value of combining this technology, interactive visualisation, argumentation and model separation to support development of multi-scale models that represent biologically plausible cells arranged in biologically plausible structures that model cell behaviour, interactions and response to therapeutic interventions

On the Interaction between TCP and the Wireless Channel in CDMA2000 Networks

In this work, we conducted extensive active measurements on a large nationwide CDMA2000 1xRTT network in order to characterize the impact of both the Radio Link Protocol and more importantly, the wireless scheduler, on TCP. Our measurements include standard TCP/UDP logs, as well as detailed RF layer statistics that allow observability into RF dynamics. With the help of a robust correlation measure, normalized mutual information, we were able to quantify the impact of these two RF factors on TCP performance metrics such as the round trip time, packet loss rate, instantaneous throughput etc. We show that the variable channel rate has the larger impact on TCP behavior when compared to the Radio Link Protocol. Furthermore, we expose and rank the factors that influence the assigned channel rate itself and in particular, demonstrate the sensitivity of the wireless scheduler to the data sending rate. Thus, TCP is adapting its rate to match the available network capacity, while the rate allocated by the wireless scheduler is influenced by the sender's behavior. Such a system is best described as a closed loop system with two feedback controllers, the TCP controller and the wireless scheduler, each one affecting the other's decisions. In this work, we take the first steps in characterizing such a system in a realistic environment

A Review of Mathematical Models for the Formation of\ud Vascular Networks

Mainly two mechanisms are involved in the formation of blood vasculature: vasculogenesis and angiogenesis. The former consists of the formation of a capillary-like network from either a dispersed or a monolayered population of endothelial cells, reproducible also in vitro by specific experimental assays. The latter consists of the sprouting of new vessels from an existing capillary or post-capillary venule. Similar phenomena are also involved in the formation of the lymphatic system through a process generally called lymphangiogenesis.\ud \ud A number of mathematical approaches have analysed these phenomena. This paper reviews the different modelling procedures, with a special emphasis on their ability to reproduce the biological system and to predict measured quantities which describe the overall processes. A comparison between the different methods is also made, highlighting their specific features