Spatial Stochastic Model of the Pre-B Cell Receptor.

Survival and proliferation of immature B lymphocytes requires expression and tonic signaling of the pre-B cell receptor (pre-BCR). This low level, ligand-independent signaling is likely achieved through frequent, but short-lived, homo interactions. Tonic signaling is also central in the pathology of precursor B acute lymphoblastic leukemia (B-ALL). In order to understand how repeated, transient events can lead to sustained signaling and to assess the impact of receptor accumulation induced by the membrane landscape, we developed a spatial stochastic model of receptor aggregation and downstream signaling events. Our rule- and agent-based model builds on previous mature BCR signaling models and incorporates novel parameters derived from single particle tracking of pre-BCR on surfaces of two different B-ALL cell lines, 697 and Nalm6. Live cell tracking of receptors on the two cell lines revealed characteristic differences in their dimer dissociation rates and diffusion coefficients. We report here that these differences affect pre-BCR aggregation and consequent signal initiation events. Receptors on Nalm6 cells, which have a lower off-rate and lower diffusion coefficient, more frequently form higher order oligomers than pre-BCR on 697 cells, resulting in higher levels of downstream phosphorylation in the Nalm6 cell line.

KrasG12D induces changes in chromatin territories that differentially impact early nuclear reprogramming in pancreatic cells.

BACKGROUND: Pancreatic ductal adenocarcinoma initiation is most frequently caused by Kras mutations. RESULTS: Here, we apply biological, biochemical, and network biology methods to validate GEMM-derived cell models using inducible KrasG12D expression. We describe the time-dependent, chromatin remodeling program that impacts function during early oncogenic signaling. We find that the KrasG12D-induced transcriptional response is dominated by downregulated expression concordant with layers of epigenetic events. More open chromatin characterizes the ATAC-seq profile associated with a smaller group of upregulated genes and epigenetic marks. RRBS demonstrates that promoter hypermethylation does not account for the silencing of the extensive gene promoter network. Moreover, ChIP-Seq reveals that heterochromatin reorganization plays little role in this early transcriptional program. Notably, both gene activation and silencing primarily depend on the marking of genes with a combination of H3K27ac, H3K4me3, and H3K36me3. Indeed, integrated modeling of all these datasets shows that KrasG12D regulates its transcriptional program primarily through unique super-enhancers and enhancers, and marking specific gene promoters and bodies. We also report chromatin remodeling across genomic areas that, although not contributing directly to cis-gene transcription, are likely important for KrasG12D functions. CONCLUSIONS: In summary, we report a comprehensive, time-dependent, and coordinated early epigenomic program for KrasG12D in pancreatic cells, which is mechanistically relevant to understanding chromatin remodeling events underlying transcriptional outcomes needed for the function of this oncogene.

Inactivation of the Euchromatic Histone-Lysine N-Methyltransferase 2 Pathway in Pancreatic Epithelial Cells Antagonizes Cancer Initiation and Pancreatitis-Associated Promotion by Altering Growth and Immune Gene Expression Networks.

Pancreatic ductal adenocarcinoma (PDAC) is an aggressive, painful disease with a 5-year survival rate of only 9%. Recent evidence indicates that distinct epigenomic landscapes underlie PDAC progression, identifying the H3K9me pathway as important to its pathobiology. Here, we delineate the role of Euchromatic Histone-lysine N-Methyltransferase 2 (EHMT2), the enzyme that generates H3K9me, as a downstream effector of oncogenic KRAS during PDAC initiation and pancreatitis-associated promotion. EHMT2 inactivation in pancreatic cells reduces H3K9me2 and antagonizes Kras G12D -mediated acinar-to-ductal metaplasia (ADM) and Pancreatic Intraepithelial Neoplasia (PanIN) formation in both the Pdx1-Cre and P48 Cre/+ Kras G12D mouse models. Ex vivo acinar explants also show impaired EGFR-KRAS-MAPK pathway-mediated ADM upon EHMT2 deletion. Notably, Kras G12D increases EHMT2 protein levels and EHMT2-EHMT1-WIZ complex formation. Transcriptome analysis reveals that EHMT2 inactivation upregulates a cell cycle inhibitory gene expression network that converges on the Cdkn1a/p21-Chek2 pathway. Congruently, pancreas tissue from Kras G12D animals with EHMT2 inactivation have increased P21 protein levels and enhanced senescence. Furthermore, loss of EHMT2 reduces inflammatory cell infiltration typically induced during Kras G12D -mediated initiation. The inhibitory effect on Kras G12D -induced growth is maintained in the pancreatitis-accelerated model, while simultaneously modifying immunoregulatory gene networks that also contribute to carcinogenesis. This study outlines the existence of a novel KRAS-EHMT2 pathway that is critical for mediating the growth-promoting and immunoregulatory effects of this oncogene in vivo, extending human observations to support a pathophysiological role for the H3K9me pathway in PDAC.

Inorganic arsenic promotes luminal to basal transition and metastasis of breast cancer.

Inorganic arsenic (iAs/As2 O32- ) is an environmental toxicant found in watersheds around the world including in densely populated areas. iAs is a class I carcinogen known to target the skin, lungs, bladder, and digestive organs, but its role as a primary breast carcinogen remains controversial. Here, we examined a different possibility: that exposure to iAs promotes the transition of well-differentiated epithelial breast cancer cells characterized by estrogen and progesterone receptor expression (ER+/PR+), to more basal phenotypes characterized by active proliferation, and propensity to metastasis in vivo. Our results indicate two clear phenotypic responses to low-level iAs that depend on the duration of the exposure. Short-term pulses of iAs activate ER signaling, consistent with its reported pseudo-estrogen activity, but longer-term, chronic treatments for over 6 months suppresses both ER and PR expression and signaling. In fact, washout of these chronically exposed cells for up to 1 month failed to fully reverse the transcriptional and phenotypic effects of prolonged treatments, indicating durable changes in cellular physiologic identity. RNA-seq studies found that chronic iAs drives the transition toward more basal phenotypes characterized by impaired hormone receptor signaling despite the conservation of estrogen receptor expression. Because treatments for breast cancer patients are largely designed based on the detection of hormone receptor expression, our results suggest greater scrutiny of ER+ cancers in patients exposed to iAs, because these tumors may spawn more aggressive phenotypes than unexposed ER+ tumors, in particular, basal subtypes that tend to develop therapy resistance and metastasis.

Effect of Spatial Inhomogeneities on the Membrane Surface on Receptor Dimerization and Signal Initiation.

Important signal transduction pathways originate on the plasma membrane, where microdomains may transiently entrap diffusing receptors. This results in a non-random distribution of receptors even in the resting state, which can be visualized as "clusters" by high resolution imaging methods. Here, we explore how spatial in-homogeneities in the plasma membrane might influence the dimerization and phosphorylation status of ErbB2 and ErbB3, two receptor tyrosine kinases that preferentially heterodimerize and are often co-expressed in cancer. This theoretical study is based upon spatial stochastic simulations of the two-dimensional membrane landscape, where variables include differential distributions and overlap of transient confinement zones ("domains") for the two receptor species. The in silico model is parameterized and validated using data from single particle tracking experiments. We report key differences in signaling output based on the degree of overlap between domains and the relative retention of receptors in such domains, expressed as escape probability. Results predict that a high overlap of domains, which favors transient co-confinement of both receptor species, will enhance the rate of hetero-interactions. Where domains do not overlap, simulations confirm expectations that homo-interactions are favored. Since ErbB3 is uniquely dependent on ErbB2 interactions for activation of its catalytic activity, variations in domain overlap or escape probability markedly alter the predicted patterns and time course of ErbB3 and ErbB2 phosphorylation. Taken together, these results implicate membrane domain organization as an important modulator of signal initiation, motivating the design of novel experimental approaches to measure these important parameters across a wider range of receptor systems.

Theory and Experimental Validation of a Spatio-temporal Model of Chemotherapy Transport to Enhance Tumor Cell Kill.

AUTHOR SUMMARY: Cancer treatment efficacy can be significantly enhanced through the elution of drug from nano-carriers that can temporarily stay in the tumor vasculature. Here we present a relatively simple yet powerful mathematical model that accounts for both spatial and temporal heterogeneities of drug dosing to help explain, examine, and prove this concept. We find that the delivery of systemic chemotherapy through a certain form of nano-carriers would have enhanced tumor kill by a factor of 2 to 4 over the standard therapy that the patients actually received. We also find that targeting blood volume fraction (a parameter of the model) through vascular normalization can achieve more effective drug delivery and tumor kill. More importantly, this model only requires a limited number of parameters which can all be readily assessed from standard clinical diagnostic measurements (e.g., histopathology and CT). This addresses an important challenge in current translational research and justifies further development of the model towards clinical translation.

Simulating cancer growth with multiscale agent-based modeling.

There have been many techniques developed in recent years to in silico model a variety of cancer behaviors. Agent-based modeling is a specific discrete-based hybrid modeling approach that allows simulating the role of diversity in cell populations as well as within each individual cell; it has therefore become a powerful modeling method widely used by computational cancer researchers. Many aspects of tumor morphology including phenotype-changing mutations, the adaptation to microenvironment, the process of angiogenesis, the influence of extracellular matrix, reactions to chemotherapy or surgical intervention, the effects of oxygen and nutrient availability, and metastasis and invasion of healthy tissues have been incorporated and investigated in agent-based models. In this review, we introduce some of the most recent agent-based models that have provided insight into the understanding of cancer growth and invasion, spanning multiple biological scales in time and space, and we further describe several experimentally testable hypotheses generated by those models. We also discuss some of the current challenges of multiscale agent-based cancer models.

Development of a diffusion-based mathematical model for predicting chemotherapy effects.

Mathematical modeling of drug transport can complement current experimental and clinical investigations to understand drug resistance mechanisms, which eventually will help to develop patient-specific chemotherapy treatments. In this paper, we present a general time- and space-dependent mathematical model based on diffusion theory for predicting chemotherapy outcome. This model has two important parameters: the blood volume fraction and radius of blood vessels divided by drug diffusion penetration length. Model analysis finds that a larger ratio of the radius of blood vessel to diffusion penetration length resulted in to a larger fraction of tumor killed, thereby leading to a better treatment outcome. Clinical translation of the model can help quantify and predict the optimal dosage size and frequency of chemotherapy for individual patients.

The signaling suppressor CIS controls proallergic T cell development and allergic airway inflammation.

Transcription factors of the STAT family are critical in the cytokine-mediated functional differentiation of CD4(+) helper T cells. Signaling inhibitors of the SOCS family negatively regulate the activation of STAT proteins; however, their roles in the differentiation and function of helper T cells are not well understood. Here we found that the SOCS protein CIS, which was substantially induced by interleukin 4 (IL-4), negatively regulated the activation of STAT3, STAT5 and STAT6 in T cells. CIS-deficient mice spontaneously developed airway inflammation, and CIS deficiency in T cells led to greater susceptibility to experimental allergic asthma. CIS-deficient T cells showed enhanced differentiation into the TH2 and TH9 subsets of helper T cells. STAT5 and STAT6 regulated IL-9 expression by directly binding to the Il9 promoter. Our data thus demonstrate a critical role for CIS in controlling the proallergic generation of helper T cells.