A generalizable data-driven multicellular model of pancreatic ductal adenocarcinoma.

BACKGROUND: Mechanistic models, when combined with pertinent data, can improve our knowledge regarding important molecular and cellular mechanisms found in cancer. These models make the prediction of tissue-level response to drug treatment possible, which can lead to new therapies and improved patient outcomes. Here we present a data-driven multiscale modeling framework to study molecular interactions between cancer, stromal, and immune cells found in the tumor microenvironment. We also develop methods to use molecular data available in The Cancer Genome Atlas to generate sample-specific models of cancer. RESULTS: By combining published models of different cells relevant to pancreatic ductal adenocarcinoma (PDAC), we built an agent-based model of the multicellular pancreatic tumor microenvironment, formally describing cell type-specific molecular interactions and cytokine-mediated cell-cell communications. We used an ensemble-based modeling approach to systematically explore how variations in the tumor microenvironment affect the viability of cancer cells. The results suggest that the autocrine loop involving EGF signaling is a key interaction modulator between pancreatic cancer and stellate cells. EGF is also found to be associated with previously described subtypes of PDAC. Moreover, the model allows a systematic exploration of the effect of possible therapeutic perturbations; our simulations suggest that reducing bFGF secretion by stellate cells will have, on average, a positive impact on cancer apoptosis. CONCLUSIONS: The developed framework allows model-driven hypotheses to be generated regarding therapeutically relevant PDAC states with potential molecular and cellular drivers indicating specific intervention strategies

Trade-off between Responsiveness and Noise Suppression in Biomolecular System Responses to Environmental Cues

When living systems detect changes in their external environment their response must be measured to balance the need to react appropriately with the need to remain stable, ignoring insignificant signals. Because this is a fundamental challenge of all biological systems that execute programs in response to stimuli, we developed a generalized time-frequency analysis (TFA) framework to systematically explore the dynamical properties of biomolecular networks. Using TFA, we focused on two well-characterized yeast gene regulatory networks responsive to carbon-source shifts and a mammalian innate immune regulatory network responsive to lipopolysaccharides (LPS). The networks are comprised of two different basic architectures. Dual positive and negative feedback loops make up the yeast galactose network; whereas overlapping positive and negative feed-forward loops are common to the yeast fatty-acid response network and the LPS-induced network of macrophages. TFA revealed remarkably distinct network behaviors in terms of trade-offs in responsiveness and noise suppression that are appropriately tuned to each biological response. The wild type galactose network was found to be highly responsive while the oleate network has greater noise suppression ability. The LPS network appeared more balanced, exhibiting less bias toward noise suppression or responsiveness. Exploration of the network parameter space exposed dramatic differences in system behaviors for each network. These studies highlight fundamental structural and dynamical principles that underlie each network, reveal constrained parameters of positive and negative feedback and feed-forward strengths that tune the networks appropriately for their respective biological roles, and demonstrate the general utility of the TFA approach for systems and synthetic biology

Co-dependence between trypanosome nuclear lamina components in nuclear stability and control of gene expression

The nuclear lamina is a filamentous structure subtending the nuclear envelope and required for chromatin organization, transcriptional regulation and maintaining nuclear structure. The trypanosomatid coiled-coil NUP-1 protein is a lamina component functionally analogous to lamins, the major lamina proteins of metazoa. There is little evidence for shared ancestry, suggesting the presence of a distinct lamina system in trypanosomes. To find additional trypanosomatid lamina components we identified NUP-1 interacting proteins by affinity capture and mass-spectrometry. Multiple components of the nuclear pore complex (NPC) and a second coiled-coil protein, which we termed NUP-2, were found. NUP-2 has a punctate distribution at the nuclear periphery throughout the cell cycle and is in close proximity to NUP-1, the NPCs and telomeric chromosomal regions. RNAi-mediated silencing of NUP-2 leads to severe proliferation defects, gross alterations to nuclear structure, chromosomal organization and nuclear envelope architecture. Further, transcription is altered at telomere-proximal variant surface glycoprotein (VSG) expression sites (ESs), suggesting a role in controlling ES expression, although NUP-2 silencing does not increase VSG switching. Transcriptome analysis suggests specific alterations to Pol I-dependent transcription. NUP-1 is mislocalized in NUP-2 knockdown cells and vice versa, implying that NUP-1 and NUP-2 form a co-dependent network and identifying NUP-2 as a second trypanosomatid nuclear lamina component

Involvement in surface antigen expression by a moonlighting FG-repeat nucleoporin in trypanosomes

Components of the nuclear periphery coordinate a multitude of activities, including macromolecular transport, cell-cycle progression, and chromatin organization. Nuclear pore complexes (NPCs) mediate nucleocytoplasmic transport, mRNA processing, and transcriptional regulation, and NPC components can define regions of high transcriptional activity in some organisms at the nuclear periphery and nucleoplasm. Lineage-specific features underpin several core nuclear functions and in trypanosomatids, which branched very early from other eukaryotes, unique protein components constitute the lamina, kinetochores, and parts of the NPCs. Here we describe a phenylalanine-glycine (FG)-repeat nucleoporin, TbNup53b, that has dual localizations within the nucleoplasm and NPC. In addition to association with nucleoporins, TbNup53b interacts with a known trans-splicing component, TSR1, and has a role in controlling expression of surface proteins including the nucleolar periphery-located, procyclin genes. Significantly, while several nucleoporins are implicated in intranuclear transcriptional regulation in metazoa, TbNup53b appears orthologous to components of the yeast/human Nup49/Nup58 complex, for which no transcriptional functions are known. These data suggest that FG-Nups are frequently co-opted to transcriptional functions during evolution and extend the presence of FG-repeat nucleoporin control of gene expression to trypanosomes, suggesting that this is a widespread and ancient eukaryotic feature, as well as underscoring once more flexibility within nucleoporin function

Patient-Specific Cell Communication Networks Associate With Disease Progression in Cancer.

The maintenance and function of tissues in health and disease depends on cell-cell communication. This work shows how high-level features, representing cell-cell communication, can be defined and used to associate certain signaling axes with clinical outcomes. We generated a scaffold of cell-cell interactions and defined a probabilistic method for creating per-patient weighted graphs based on gene expression and cell deconvolution results. With this method, we generated over 9,000 graphs for The Cancer Genome Atlas (TCGA) patient samples, each representing likely channels of intercellular communication in the tumor microenvironment (TME). It was shown that cell-cell edges were strongly associated with disease severity and progression, in terms of survival time and tumor stage. Within individual tumor types, there are predominant cell types, and the collection of associated edges were found to be predictive of clinical phenotypes. Additionally, genes associated with differentially weighted edges were enriched in Gene Ontology terms associated with tissue structure and immune response. Code, data, and notebooks are provided to enable the application of this method to any expression dataset (https://github.com/IlyaLab/Pan-Cancer-Cell-Cell-Comm-Net)

One-Cell Doubling Evaluation by Living Arrays of Yeast, ODELAY!

Cell growth is a complex phenotype widely used in systems biology to gauge the impact of genetic and environmental perturbations. Due to the magnitude of genome-wide studies, resolution is often sacrificed in favor of throughput, creating a demand for scalable, time-resolved, quantitative methods of growth assessment. We present ODELAY (One-cell Doubling Evaluation by Living Arrays of Yeast), an automated and scalable growth analysis platform. High measurement density and single-cell resolution provide a powerful tool for large-scale multiparameter growth analysis based on the modeling of microcolony expansion on solid media. Pioneered in yeast but applicable to other colony forming organisms, ODELAY extracts the three key growth parameters (lag time, doubling time, and carrying capacity) that define microcolony expansion from single cells, simultaneously permitting the assessment of population heterogeneity. The utility of ODELAY is illustrated using yeast mutants, revealing a spectrum of phenotypes arising from single and combinatorial growth parameter perturbations

ADAPTS: Automated deconvolution augmentation of profiles for tissue specific cells.

Immune cell infiltration of tumors and the tumor microenvironment can be an important component for determining patient outcomes. For example, immune and stromal cell presence inferred by deconvolving patient gene expression data may help identify high risk patients or suggest a course of treatment. One particularly powerful family of deconvolution techniques uses signature matrices of genes that uniquely identify each cell type as determined from single cell type purified gene expression data. Many methods from this family have been recently published, often including new signature matrices appropriate for a single purpose, such as investigating a specific type of tumor. The package ADAPTS helps users make the most of this expanding knowledge base by introducing a framework for cell type deconvolution. ADAPTS implements modular tools for customizing signature matrices for new tissue types by adding custom cell types or building new matrices de novo, including from single cell RNAseq data. It includes a common interface to several popular deconvolution algorithms that use a signature matrix to estimate the proportion of cell types present in heterogenous samples. ADAPTS also implements a novel method for clustering cell types into groups that are difficult to distinguish by deconvolution and then re-splitting those clusters using hierarchical deconvolution. We demonstrate that the techniques implemented in ADAPTS improve the ability to reconstruct the cell types present in a single cell RNAseq data set in a blind predictive analysis. ADAPTS is currently available for use in R on CRAN and GitHub