Single-cell transcriptomics reveals stage- and side-specificity of gene modules in colorectal cancer.

BACKGROUND: An understanding of mechanisms underlying colorectal cancer (CRC) development and progression is yet to be fully elucidated. This study aims to employ network theoretic approaches to analyse single cell transcriptomic data from CRC to better characterize its progression and sided-ness. METHODS: We utilized a recently published single-cell RNA sequencing data (GEO-GSE178341) and parsed the cell X gene data by stage and side (right and left colon). Using Weighted Gene Co-expression Network Analysis (WGCNA), we identified gene modules with varying preservation levels (weak or strong) of network topology between early (pT1) and late stages (pT234), and between right and left colons. Spearman's rank correlation (ρ) was used to assess the similarity or dissimilarity in gene connectivity. RESULTS: Equalizing cell counts across different stages, we detected 13 modules for the early stage, two of which were non-preserved in late stages. Both non-preserved modules displayed distinct gene connectivity patterns between the early and late stages, characterized by low ρ values. One module predominately dealt with myeloid cells, with genes mostly enriched for cytokine-cytokine receptor interaction potentiallystimulating myeloid cells to participate in angiogenesis. The second module, representing a subset of epithelial cells, was mainly enriched for carbohydrate digestion and absorption, influencing the gut microenvironment through the breakdown of carbohydrates. In the comparison of left vs. right colons, two of 12 modules identified in the right colon were non-preserved in the left colon. One captured a small fraction of epithelial cells and was enriched for transcriptional misregulation in cancer, potentially impacting communication between epithelial cells and the tumor microenvironment. The other predominantly contained B cells with a crucial role in maintaining human gastrointestinal health and was enriched for B-cell receptor signalling pathway. CONCLUSIONS: We identified modules with topological and functional differences specific to cell types between the early and late stages, and between the right and left colons. This study enhances the understanding of roles played by different cell types at different stages and sides, providing valuable insights for future studies focused on the diagnosis and treatment of CRC.

Stochastic operator-splitting method for reaction-diffusion systems.

Many biochemical processes at the sub-cellular level involve a small number of molecules. The local numbers of these molecules vary in space and time, and exhibit random fluctuations that can only be captured with stochastic simulations. We present a novel stochastic operator-splitting algorithm to model such reaction-diffusion phenomena. The reaction and diffusion steps employ stochastic simulation algorithms and Brownian dynamics, respectively. Through theoretical analysis, we have developed an algorithm to identify if the system is reaction-controlled, diffusion-controlled or is in an intermediate regime. The time-step size is chosen accordingly at each step of the simulation. We have used three examples to demonstrate the accuracy and robustness of the proposed algorithm. The first example deals with diffusion of two chemical species undergoing an irreversible bimolecular reaction. It is used to validate our algorithm by comparing its results with the solution obtained from a corresponding deterministic partial differential equation at low and high number of molecules. In this example, we also compare the results from our method to those obtained using a Gillespie multi-particle (GMP) method. The second example, which models simplified RNA synthesis, is used to study the performance of our algorithm in reaction- and diffusion-controlled regimes and to investigate the effects of local inhomogeneity. The third example models reaction-diffusion of CheY molecules through the cytoplasm of Escherichia coli during chemotaxis. It is used to compare the algorithm's performance against the GMP method. Our analysis demonstrates that the proposed algorithm enables accurate simulation of the kinetics of complex and spatially heterogeneous systems. It is also computationally more efficient than commonly used alternatives, such as the GMP method.

Identification of crosstalk between phosphoprotein signaling pathways in RAW 264.7 macrophage cells.

Signaling pathways mediate the effect of external stimuli on gene expression in cells. The signaling proteins in these pathways interact with each other and their phosphorylation levels often serve as indicators for the activity of signaling pathways. Several signaling pathways have been identified in mammalian cells but the crosstalk between them is not well understood. Alliance for Cellular Signaling (AfCS) has measured time-course data in RAW 264.7 macrophage cells on important phosphoproteins, such as the mitogen-activated protein kinases (MAPKs) and signal transducer and activator of transcription (STATs), in single- and double-ligand stimulation experiments for 22 ligands. In the present work, we have used a data-driven approach to analyze the AfCS data to decipher the interactions and crosstalk between signaling pathways in stimulated macrophage cells. We have used dynamic mapping to develop a predictive model using a partial least squares approach. Significant interactions were selected through statistical hypothesis testing and were used to reconstruct the phosphoprotein signaling network. The proposed data-driven approach is able to identify most of the known signaling interactions such as protein kinase B (Akt) --> glycogen synthase kinase 3alpha/beta (GSKalpha/beta) etc., and predicts potential novel interactions such as P38 --> RSK and GSK --> ezrin/radixin/moesin. We have also shown that the model has good predictive power for extrapolation. Our novel approach captures the temporal causality and directionality in intracellular signaling pathways. Further, case specific analysis of the phosphoproteins in the network has led us to propose hypothesis about inhibition (phosphorylation) of GSKalpha/beta via P38.

Stochastic hybrid modeling of intracellular calcium dynamics.

Deterministic models of biochemical processes at the subcellular level might become inadequate when a cascade of chemical reactions is induced by a few molecules. Inherent randomness of such phenomena calls for the use of stochastic simulations. However, being computationally intensive, such simulations become infeasible for large and complex reaction networks. To improve their computational efficiency in handling these networks, we present a hybrid approach, in which slow reactions and fluxes are handled through exact stochastic simulation and their fast counterparts are treated partially deterministically through chemical Langevin equation. The classification of reactions as fast or slow is accompanied by the assumption that in the time-scale of fast reactions, slow reactions do not occur and hence do not affect the probability of the state. Our new approach also handles reactions with complex rate expressions such as Michaelis-Menten kinetics. Fluxes which cannot be modeled explicitly through reactions, such as flux of Ca(2+) from endoplasmic reticulum to the cytosol through inositol 1,4,5-trisphosphate receptor channels, are handled deterministically. The proposed hybrid algorithm is used to model the regulation of the dynamics of cytosolic calcium ions in mouse macrophage RAW 264.7 cells. At relatively large number of molecules, the response characteristics obtained with the stochastic and deterministic simulations coincide, which validates our approach in the limit of large numbers. At low doses, the response characteristics of some key chemical species, such as levels of cytosolic calcium, predicted with stochastic simulations, differ quantitatively from their deterministic counterparts. These observations are ubiquitous throughout dose response, sensitivity, and gene-knockdown response analyses. While the relative differences between the peak-heights of the cytosolic [Ca(2+)] time-courses obtained from stochastic (mean of 16 realizations) and deterministic simulations are merely 1%-4% for most perturbations, it is specially sensitive to levels of G(ßγ) (relative difference as large as 90% at very low G(ßγ)).

An integrated model of eicosanoid metabolism and signaling based on lipidomics flux analysis.

There is increasing evidence for a major and critical involvement of lipids in signal transduction and cellular trafficking, and this has motivated large-scale studies on lipid pathways. The Lipid Metabolites and Pathways Strategy consortium is actively investigating lipid metabolism in mammalian cells and has made available time-course data on various lipids in response to treatment with KDO(2)-lipid A (a lipopolysaccharide analog) of macrophage RAW 264.7 cells. The lipids known as eicosanoids play an important role in inflammation. We have reconstructed an integrated network of eicosanoid metabolism and signaling based on the KEGG pathway database and the literature and have developed a kinetic model. A matrix-based approach was used to estimate the rate constants from experimental data and these were further refined using generalized constrained nonlinear optimization. The resulting model fits the experimental data well for all species, and simulated enzyme activities were similar to their literature values. The quantitative model for eicosanoid metabolism that we have developed can be used to design experimental studies utilizing genetic and pharmacological perturbations to probe fluxes in lipid pathways.

Systems biology of macrophages.

Cells and tissues function in context. Under a given growth or survival medium they perform tasks, replicate and die. Given a stimulus they respond by invoking myriad biomolecular networks that result in a specified cellular outcome. At any given instant it can be argued that the cell is in a "state" defined by its components--their concentrations and locations, the interactions between components--that are modulated in space and time, and the complex circuitry--that involves a large number of interacting networks and a snapshot of the dynamical processes--such as gene expression, cell cycle, transport of components, etc. At present, we can measure, using high and low throughput methods, several cellular components in a context-dependent manner and obtain a partial picture of cellular networks and dynamical processes. Are these measurements sufficient to answer important biological questions and help reconstruct a systems-level understanding of a mammalian cell? This chapter will address systems biology strategies developed to address this question and demonstrate the power of integration of diverse cellular data for answering interesting biological questions in macrophages. We will use this systems biology approach to address the following questions: (1) How good are macrophage cell lines in addressing phenotypic biology of primary macrophages? (2) How do signals associated with inflammatory molecules regulate gene transcription in macrophages? (3) How can we combine proteomic and other cellular measurements to characterize the repertoire of upstream signaling networks invoked by macrophages? (4) How do designed knockdowns of proteins influence cellular phenotypes?

Time-varying causal inference from phosphoproteomic measurements in macrophage cells.

Cellular signaling circuitry in eukaryotes can be studied by analyzing the regulation of protein phosphorylation and its impact on downstream mechanisms leading to a phenotype. A primary role of phosphorylation is to act as a switch to turn "on" or "off" a protein activity or a cellular pathway. Specifically, protein phosphorylation is a major leit motif for transducing molecular signals inside the cell. Errors in transferring cellular information can alter the normal function and may lead to diseases such as cancer; an accurate reconstruction of the "true" signaling network is essential for understanding the molecular machinery involved in normal and pathological function. In this study, we have developed a novel framework for time-dependent reconstruction of signaling networks involved in the activation of macrophage cells leading to an inflammatory response. Several signaling pathways have been identified in macrophage cells, but the time-varying causal relationship that can produce a dynamic directed graph of these molecules has not been explored in detail. Here, we use the notion of Granger causality, and apply a vector autoregressive model to phosphoprotein time-course data in RAW 264.7 macrophage cells. Through the reconstruction of the phosphoprotein network, we were able to estimate the directionality and the dynamics of information flow. Significant interactions were selected through statistical hypothesis testing ( t-test) of the coefficients of a linear model and were used to reconstruct the phosphoprotein signaling network. Our approach results in a three-stage phosphoprotein network that represents the evolution of the causal interactions in the intracellular signaling pathways.

Information theoretic approach to complex biological network reconstruction: application to cytokine release in RAW 264.7 macrophages.

BACKGROUND: High-throughput methods for biological measurements generate vast amounts of quantitative data, which necessitate the development of advanced approaches to data analysis to help understand the underlying mechanisms and networks. Reconstruction of biological networks from measured data of different components is a significant challenge in systems biology. RESULTS: We use an information theoretic approach to reconstruct phosphoprotein-cytokine networks in RAW 264.7 macrophage cells. Cytokines are secreted upon activation of a wide range of regulatory signals transduced by the phosphoprotein network. Identifying these components can help identify regulatory modules responsible for the inflammatory phenotype. The information theoretic approach is based on estimation of mutual information of interactions by using kernel density estimators. Mutual information provides a measure of statistical dependencies between interacting components. Using the topology of the network derived, we develop a data-driven parsimonious input-output model of the phosphoprotein-cytokine network. CONCLUSIONS: We demonstrate the applicability of our information theoretic approach to reconstruction of biological networks. For the phosphoprotein-cytokine network, this approach not only captures most of the known signaling components involved in cytokine release but also predicts new signaling components involved in the release of cytokines. The results of this study are important for gaining a clear understanding of macrophage activation during the inflammation process.

A kinetic model for calcium dynamics in RAW 264.7 cells: 1. Mechanisms, parameters, and subpopulational variability.

Calcium (Ca(2+)) is an important second messenger and has been the subject of numerous experimental measurements and mechanistic studies in intracellular signaling. Calcium profile can also serve as a useful cellular phenotype. Kinetic models of calcium dynamics provide quantitative insights into the calcium signaling networks. We report here the development of a complex kinetic model for calcium dynamics in RAW 264.7 cells stimulated by the C5a ligand. The model is developed using the vast number of measurements of in vivo calcium dynamics carried out in the Alliance for Cellular Signaling (AfCS) Laboratories. Ligand binding, phospholipase C-beta (PLC-beta) activation, inositol 1,4,5-trisphosphate (IP(3)) receptor (IP(3)R) dynamics, and calcium exchange with mitochondria and extracellular matrix have all been incorporated into the model. The experimental data include data from both native and knockdown cell lines. Subpopulational variability in measurements is addressed by allowing nonkinetic parameters to vary across datasets. The model predicts temporal response of Ca(2+) concentration for various doses of C5a under different initial conditions. The optimized parameters for IP(3)R dynamics are in agreement with the legacy data. Further, the half-maximal effect concentration of C5a and the predicted dose response are comparable to those seen in AfCS measurements. Sensitivity analysis shows that the model is robust to parametric perturbations.

A kinetic model for calcium dynamics in RAW 264.7 cells: 2. Knockdown response and long-term response.

This article addresses how quantitative models such as the one proposed in the companion article can be used to study cellular network perturbations such as knockdowns and pharmacological perturbations in a predictive manner. Using the kinetic model for cytosolic calcium dynamics in RAW 264.7 cells developed in the companion article, the calcium response to complement 5a (C5a) for the knockdown of seven proteins (C5a receptor; G-beta-2; G-alpha,i-2,3; regulator of G-protein signaling-10; G-protein coupled receptor kinase-2; phospholipase C beta-3; arrestin) is predicted and validated against the data from the Alliance for Cellular Signaling. The knockdown responses provide insights into how altered expressions of important proteins in disease states result in intermediate measurable phenotypes. Long-term response and long-term dose response have also been predicted, providing insights into how the receptor desensitization, internalization, and recycle result in tolerance. Sensitivity analysis of long-term response shows that the mechanisms and parameters in the receptor recycle path are important for long-term calcium dynamics.