Research in the lab is built on our multi-disciplinary experience in bioengineering, quantitative pharmacology, and single-cell measurement and analysis, and benefits from active collaboration with leaders in clinical oncology, epigenetics and proteomics research. The lab is currently focused on the following areas of research:

Multi-parametric analysis of the heterogeneous cellular drug response
Analyses of cellular response to cancer drugs have been historically focused on variation in drug potency (e.g. IC50) derived from conventional dose-response assays. Assuming that this parameter is the most important difference between effective and ineffective drugs or between drug-sensitive and resistant cells, these assays have played a key role in cancer research and drug discovery. However, taking a multi-parametric information theory-based approach, our recent work has shown that variations in maximal effect (Emax) and the slope of the dose-response curve are as important in limiting the effectiveness of many of the cancer drugs, including targeted therapeutics (e.g. inhibitors of PI3K/AKT/mTOR or RAF/MEK signaling) and cytotoxic compounds (e.g. DNA-damaging agents). Our recent work has addressed the causes of such variations and ascribed them to cell-to-cell heterogeneity in the biochemical and phenotypic nature, extent and timing of drug response, even in genetically homogeneous populations of cancer cells. An important goal of our research is to develop novel and simple, but efficient experimental and analytical methods to precisely score drug and drug-combination effectiveness in a way that captures complexities overlooked in bulk assays. These methods can provide a basis for subsequent studies to elucidate the mechanisms and biomarkers of drug sensitivity, adaptation or resistance that will ultimately guide novel therapies.

  • Fallahi-Sichani M, Honarnejad S, Heiser LM, Gray JW, Sorger PK. Metrics other than potency reveal systematic variation in responses to cancer drugs, Nature Chemical Biology (2013). (PubMed Link) (PDF)
  • Lin JR, Fallahi-Sichani M, Sorger PK. Highly multiplexed imaging of single cells using a high-throughput cyclic immunofluorescence method, Nature Communications (2015). (PubMed Link) (PDF)
  • Lin JR, Fallahi-Sichani M, Chen JY, Sorger PK. Cyclic immunofluorescence (CycIF), a highly multiplexed method for single-cell imaging, Current Protocols in Chemical Biology (2016). (PubMed Link) (PDF)

Studies on adaptive regulation of cancer cell fate following targeted inhibition of oncogenic pathways
An increased understanding of the molecular biology of oncogenesis provides new opportunities for the design of molecularly targeted cancer therapies. However, therapeutic effects are often temporary and limited by drug adaptation, partial response and ultimately resistance. Understanding and ultimately preventing adaptation is likely to be key to durable therapy. Despite a wealth of data on signaling networks involved in adaptation, current understanding of these networks is incomplete and many questions remain unanswered: Is adaptation fundamentally similar across different cell types or individual cells that experience the same oncogenic pathway? How informative is the initial (pre-treatment) state for drug adaptation or resistance? How do adaptive changes in network state depend on the biochemical node inhibited and how do they translate to subsequent phenotypic states (i.e. cell death, proliferation, quiescence or differentiation)? What is the relative contribution of cell-autonomous versus non-cell-autonomous factors? Can we identify single or multi-component biomarkers of adaptation to guide precision treatment strategies that mitigate drug resistance? In addition, most of our knowledge about adaptive resistance comes from studying bulk tumor cell populations. It therefore remains unclear whether individual cells in a population adapt to drug effect (i.e. inhibition of oncogenic signaling) in the same way or in different ways. To address these important questions, we use a combination of long-term live-cell microscopy, high-throughput biochemical measurement, multiplex single-cell analysis and computational modeling to quantify signaling mechanisms involved in adaptive drug responses and their consequences for cancer cell fate and to use that knowledge to propose novel approaches to enhance drug maximal effect and to prevent or overcome drug resistance. Our work is currently focused on BRAF-mutant cancers and drugs targeting the BRAF oncoprotein or other components of ERK signaling, but the approach is applicable to other cancers and other oncogenic pathways.

  • Fallahi-Sichani M#, Becker V, Izar B, Baker GJ, Lin JR, Boswell SA, Shah P, Rotem A, Garraway LA, Sorger PK#. Adaptive resistance of melanoma cells to RAF inhibition via reversible induction of a slowly dividing de-differentiated state, Molecular Systems Biology (2017). (PubMed Link) (PDF (#Co-corresponding authors)
  • Fallahi-Sichani M, Moerke NJ, Niepel M, Zhang T, Gray NS, Sorger PK. Systematic analysis of BRAF(V600E) melanomas reveals a role for JNK/c-Jun pathway in adaptive resistance to drug-induced apoptosis, Molecular Systems Biology (2015). (PubMed Link) (PDF)
  • Moerke N, Fallahi-Sichani M. Reverse phase protein arrays for compound profiling, Current Protocols in Chemical Biology (2016). (PubMed Link) (PDF)
  • Tirosh I, Izar B, Prakadan SM, Wadsworth MH, Treacy D, Trombetta JJ, Rotem A, Rodman C, Lian C, Murphy G, Fallahi-Sichani M, Dutton-Regester K, Lin JR, Cohen O, Shah P, Lu D, Genshaft AS, Hughes TK, Ziegler CGK, Kazer SW, Gaillard A, Kolb KE, Villani AC, Johannessen CM, Andreev AY, Van Allen EM, Bertagnolli M, Sorger PK, Sullivan RJ, Flaherty KT, Frederick DT, Jané-Valbuena J, Yoon C, Rozenblatt-Rosen O, Shalek AK, Regev A, Garraway LA. Dissecting the multicellular ecosystem of metastatic melanoma by single-cell RNA-seq, Science (2016). (PubMed Link) (PDF)

Data-driven and multi-scale modeling of complex biological processes and bio-molecular networks
Advance in biomedical research is increasingly dependent on quantitative experimental methods that use emerging technologies to profile complex and multi-scale biological systems. At the same time sophisticated analytical methods are required to dissect the complex interactions between the measured components and uncover how these interactions at a range of spatial and temporal scales give rise to the function and behavior of the system. We have contributed to the development of both experimental and computational procedures for analyzing cellular phenotypes and biochemical networks and modeling multi-scale molecular, cellular and tissue-level mechanisms. In the past, we have developed multi-scale computational modeling approaches that combine stochastic (Monte Carlo) and deterministic (ordinary or partial differential equation) modeling to investigate the dynamics of G-protein coupled receptor signaling, and to study the mechanisms by which TNF signaling determines immunity to M. tuberculosis infection. We are working to apply similar and novel approaches to the studies of other complex diseases, particularly cancer, where there is a critical need for bridging studies of therapeutic targets at the molecular scale to the effects of modulation of these targets at the cellular and tissue (tumor) levels. More recently, we have used statistical modeling approaches such as partial least squares regression and mutual information analysis with the goal of uncovering important determinants of regulation from complex high-dimensional data spaces and to identify measurements of network activity that are most predictive of a relevant phenotype. These and other mathematical approaches that capture dynamic, multi-factorial differences in signaling and biochemical networks across cell types and drugs are also essential components of our research.

  • Kirschner DE, Hunt CA, Marino S, Fallahi-Sichani M, Linderman JJ. Tuneable resolution as a systems biology approach for multi-scale, multi-compartment computational models, WIREs Systems Biology and Medicine (2014). (PubMed Link) (PDF)
  • Fallahi-Sichani M, Moerke NJ, Niepel M, Zhang T, Gray NS, Sorger PK. Systematic analysis of BRAF(V600E) melanomas reveals a role for JNK/c-Jun pathway in adaptive resistance to drug-induced apoptosis, Molecular Systems Biology (2015). (PubMed Link) (PDF)
  • Fallahi-Sichani M, Kirschner DE, Linderman JJ. NF-kappaB signaling dynamics play a key role in infection control in tuberculosis, Frontiers in Physiology 3:170 (2012). (PubMed Link) (PDF)
  • Fallahi-Sichani M, Flynn JL, Linderman JJ, Kirschner DE. Differential risk of tuberculosis reactivation among anti-TNF therapies is due to drug binding kinetics and permeability, Journal of Immunology (2012). (PubMed Link) (PDF)
  • Fallahi-Sichani M, El-Kebir M, Marino S, Kirschner DE, Linderman JJ. Multiscale computational modeling reveals a critical role for TNF-α receptor 1 dynamics in tuberculosis granuloma formation, Journal of Immunology (2011). (PubMed Link) (PDF)
  • Fallahi-Sichani M, Schaller MA, Kirschner DE, Kunkel SL, Linderman JJ. Identification of key processes that control tumor necrosis factor availability in a tuberculosis granuloma, PLoS Computational Biology (2010). (PubMed Link) (PDF)
  • Fallahi-Sichani M, Linderman JJ, Lipid raft-mediated regulation of G-protein coupled receptor signaling by ligands which influence receptor dimerization: A computational study, PLoS ONE (2009). (PubMed Link) (PDF)
  • Fallahi-Sichani M, Marino S, Flynn JL, Linderman JJ, Kirschner DE. A systems biology approach for understanding granuloma formation and function in tuberculosis, In McFadden J, Beste D, Kierzek A (Ed.), Systems biology of tuberculosis. Springer (2013). (Publisher Link) (PDF)
  • Marino S, Fallahi-Sichani M, Linderman JJ, Kirschner DE. Mathematical models of anti-TNF therapies and their correlation with tuberculosis, In Pathak Y, Benita S (Ed.), Antibody-mediated drug delivery systems: Concepts, Technology and Applications. John Wiley and Sons (2012). (Publisher Link) (PDF)