Bisbiguanide analogs induce mitochondrial stress to inhibit lung cancer cell invasion.

Targeting cancer metabolism to limit cellular energy and metabolite production is an attractive therapeutic approach. Here, we developed analogs of the bisbiguanide, alexidine, to target lung cancer cell metabolism and assess a structure-activity relationship (SAR). The SAR led to the identification of two analogs, AX-4 and AX-7, that limit cell growth via G1/G0 cell-cycle arrest and are tolerated in vivo with favorable pharmacokinetics. Mechanistic evaluation revealed that AX-4 and AX-7 induce potent mitochondrial defects; mitochondrial cristae were deformed and the mitochondrial membrane potential was depolarized. Additionally, cell metabolism was rewired, as indicated by reduced oxygen consumption and mitochondrial ATP production, with an increase in extracellular lactate. Importantly, AX-4 and AX-7 impacted overall cell behavior, as these compounds reduced collective cell invasion. Taken together, our study establishes a class of bisbiguanides as effective mitochondria and cell invasion disrupters, and proposes bisbiguanides as promising approaches to limiting cancer metastasis.

Subpopulation commensalism promotes Rac1-dependent invasion of single cells via laminin-332.

Phenotypic heterogeneity poses a significant hurdle for cancer treatment but is under-characterized in the context of tumor invasion. Amidst the range of phenotypic heterogeneity across solid tumor types, collectively invading cells and single cells have been extensively characterized as independent modes of invasion, but their intercellular interactions have rarely been explored. Here, we isolated collectively invading cells and single cells from the heterogeneous 4T1 cell line and observed extensive transcriptional and epigenetic diversity across these subpopulations. By integrating these datasets, we identified laminin-332 as a protein complex exclusively secreted by collectively invading cells. Live-cell imaging revealed that laminin-332 derived from collectively invading cells increased the velocity and directionality of single cells. Despite collectively invading and single cells having similar expression of the integrin α6ß4 dimer, single cells demonstrated higher Rac1 activation upon laminin-332 binding to integrin α6ß4. This mechanism suggests a novel commensal relationship between collectively invading and single cells, wherein collectively invading cells promote the invasive potential of single cells through a laminin-332/Rac1 axis.

A live-cell platform to isolate phenotypically defined subpopulations for spatial multi-omic profiling.

Numerous techniques have been employed to deconstruct the heterogeneity observed in normal and diseased cellular populations, including single cell RNA sequencing, in situ hybridization, and flow cytometry. While these approaches have revolutionized our understanding of heterogeneity, in isolation they cannot correlate phenotypic information within a physiologically relevant live-cell state with molecular profiles. This inability to integrate a live-cell phenotype-such as invasiveness, cell:cell interactions, and changes in spatial positioning-with multi-omic data creates a gap in understanding cellular heterogeneity. We sought to address this gap by employing lab technologies to design a detailed protocol, termed Spatiotemporal Genomic and Cellular Analysis (SaGA), for the precise imaging-based selection, isolation, and expansion of phenotypically distinct live cells. This protocol requires cells expressing a photoconvertible fluorescent protein and employs live cell confocal microscopy to photoconvert a user-defined single cell or set of cells displaying a phenotype of interest. The total population is then extracted from its microenvironment, and the optically highlighted cells are isolated using fluorescence activated cell sorting. SaGA-isolated cells can then be subjected to multi-omics analysis or cellular propagation for in vitro or in vivo studies. This protocol can be applied to a variety of conditions, creating protocol flexibility for user-specific research interests. The SaGA technique can be accomplished in one workday by non-specialists and results in a phenotypically defined cellular subpopulations for integration with multi-omics techniques. We envision this approach providing multi-dimensional datasets exploring the relationship between live cell phenotypes and multi-omic heterogeneity within normal and diseased cellular populations.

Live-cell invasive phenotyping uncovers the ALK2/BMP6 iron homeostasis pathway as a therapeutic vulnerability in LKB1-mutant lung cancer.

The acquisition of invasive properties is a prerequisite for tumor progression and metastasis. Molecular subtypes of KRAS-driven lung cancer exhibit distinct modes of invasion that likely contribute to unique growth properties and therapeutic susceptibilities. Despite this, pre-clinical discovery strategies designed to exploit invasive phenotypes are lacking. To address this, we designed an experimental system to screen for targetable signaling pathways linked to active early invasion phenotypes in the two most prominent molecular subtypes, TP53 and LKB1, of KRAS-driven lung adenocarcinoma (LUAD). By combining live-cell imaging of human bronchial epithelial cells in a 3D invasion matrix with RNA transcriptome profiling, we identified the LKB1-specific upregulation of bone morphogenetic protein 6 (BMP6). Examination of early-stage lung cancer patients confirmed upregulation of BMP6 in LKB1-mutant lung tumors. At the molecular level, we find that the canonical iron regulatory hormone Hepcidin is induced via BMP6 signaling upon LKB1 loss, where intact LKB1 kinase activity is necessary to maintain signaling homeostasis. Furthermore, pre-clinical studies in a novel Kras/Lkb1-mutant syngeneic mouse model show that potent growth suppression was achieved by inhibiting the ALK2/BMP6 signaling axis with single agents that are currently in clinical trials. We show that alterations in the iron homeostasis pathway are accompanied by simultaneous upregulation of ferroptosis protection proteins. Thus, LKB1 is sufficient to regulate both the 'gas' and 'breaks' to finely tune iron-regulated tumor progression.

A live-cell platform to isolate phenotypically defined subpopulations for spatial multi-omic profiling.

Numerous techniques have been employed to deconstruct the heterogeneity observed in normal and diseased cellular populations, including single cell RNA sequencing, in situ hybridization, and flow cytometry. While these approaches have revolutionized our understanding of heterogeneity, in isolation they cannot correlate phenotypic information within a physiologically relevant live-cell state, with molecular profiles. This inability to integrate a historical live-cell phenotype, such as invasiveness, cell:cell interactions, and changes in spatial positioning, with multi-omic data, creates a gap in understanding cellular heterogeneity. We sought to address this gap by employing lab technologies to design a detailed protocol, termed Spatiotemporal Genomics and Cellular Analysis (SaGA), for the precise imaging-based selection, isolation, and expansion of phenotypically distinct live-cells. We begin with cells stably expressing a photoconvertible fluorescent protein and employ live cell confocal microscopy to photoconvert a user-defined single cell or set of cells displaying a phenotype of interest. The total population is then extracted from its microenvironment, and the optically highlighted cells are isolated using fluorescence activated cell sorting. SaGA-isolated cells can then be subjected to multi-omics analysis or cellular propagation for in vitro or in vivo studies. This protocol can be applied to a variety of conditions, creating protocol flexibility for user-specific research interests. The SaGA technique can be accomplished in one workday by non-specialists and results in a phenotypically defined cellular subpopulation for integration with multi-omics techniques. We envision this approach providing multi-dimensional datasets exploring the relationship between live-cell phenotype and multi-omic heterogeneity within normal and diseased cellular populations.

Loss of the endocytic tumor suppressor HD-PTP phenocopies LKB1 and promotes RAS-driven oncogenesis.

Oncogenic RAS mutations drive aggressive cancers that are difficult to treat in the clinic, and while direct inhibition of the most common KRAS variant in lung adenocarcinoma (G12C) is undergoing clinical evaluation, a wide spectrum of oncogenic RAS variants together make up a large percentage of untargetable lung and GI cancers. Here we report that loss-of-function alterations (mutations and deep deletions) in the gene that encodes HD-PTP (PTPN23) occur in up to 14% of lung cancers in the ORIEN Avatar lung cancer cohort, associate with adenosquamous histology, and occur alongside an altered spectrum of KRAS alleles. Furthermore, we show that in publicly available early-stage NSCLC studies loss of HD-PTP is mutually exclusive with loss of LKB1, which suggests they restrict a common oncogenic pathway in early lung tumorigenesis. In support of this, knockdown of HD-PTP in RAS-transformed lung cancer cells is sufficient to promote FAK-dependent invasion. Lastly, knockdown of the Drosophila homolog of HD-PTP (dHD-PTP/Myopic) synergizes to promote RAS-dependent neoplastic progression. Our findings highlight a novel tumor suppressor that can restrict RAS-driven lung cancer oncogenesis and identify a targetable pathway for personalized therapeutic approaches for adenosquamous lung cancer.

Epigenetically heterogeneous tumor cells direct collective invasion through filopodia-driven fibronectin micropatterning.

Tumor heterogeneity drives disease progression, treatment resistance, and patient relapse, yet remains largely underexplored in invasion and metastasis. Here, we investigated heterogeneity within collective cancer invasion by integrating DNA methylation and gene expression analysis in rare purified lung cancer leader and follower cells. Our results showed global DNA methylation rewiring in leader cells and revealed the filopodial motor MYO10 as a critical gene at the intersection of epigenetic heterogeneity and three-dimensional (3D) collective invasion. We further identified JAG1 signaling as a previously unknown upstream activator of MYO10 expression in leader cells. Using live-cell imaging, we found that MYO10 drives filopodial persistence necessary for micropatterning extracellular fibronectin into linear tracks at the edge of 3D collective invasion exclusively in leaders. Our data fit a model where epigenetic heterogeneity and JAG1 signaling jointly drive collective cancer invasion through MYO10 up-regulation in epigenetically permissive leader cells, which induces filopodia dynamics necessary for linearized fibronectin micropatterning.

Prognostic significance of an invasive leader cell-derived mutation cluster on chromosome 16q.

BACKGROUND: Intratumoral heterogeneity is defined by subpopulations with varying genotypes and phenotypes. Specialized, highly invasive leader cells and less invasive follower cells are phenotypically distinct subpopulations that cooperate during collective cancer invasion. Because leader cells are a rare subpopulation that would be missed by bulk sequencing, a novel image-guided genomics platform was used to precisely select this subpopulation. This study identified a novel leader cell mutation signature and tested its ability to predict prognosis in non-small cell lung cancer (NSCLC) patient cohorts. METHODS: Spatiotemporal genomic and cellular analysis was used to isolate and perform RNA sequencing on leader and follower populations from the H1299 NSCLC cell line, and it revealed a leader-specific mutation cluster on chromosome 16q. Genomic data from patients with lung squamous cell carcinoma (LUSC; n = 475) and lung adenocarcinoma (LUAD; n = 501) from The Cancer Genome Atlas were stratified by 16q mutation cluster (16qMC) status (16qMC+ vs 16qMC-) and compared for overall survival (OS), progression-free survival (PFS), and gene set enrichment analysis (GSEA). RESULTS: Poorer OS, poorer PFS, or both were found across all stages and among early-stage patients with 16qMC+ tumors within the LUSC and LUAD cohorts. GSEA revealed 16qMC+ tumors to be enriched for the expression of metastasis- and survival-associated gene sets. CONCLUSIONS: This represents the first leader cell mutation signature identified in patients and has the potential to better stratify high-risk NSCLC and ultimately improve patient outcomes.

Genetic heterogeneity within collective invasion packs drives leader and follower cell phenotypes.

Collective invasion, the coordinated movement of cohesive packs of cells, has become recognized as a major mode of metastasis for solid tumors. These packs are phenotypically heterogeneous and include specialized cells that lead the invasive pack and others that follow behind. To better understand how these unique cell types cooperate to facilitate collective invasion, we analyzed transcriptomic sequence variation between leader and follower populations isolated from the H1299 non-small cell lung cancer cell line using an image-guided selection technique. We now identify 14 expressed mutations that are selectively enriched in leader or follower cells, suggesting a novel link between genomic and phenotypic heterogeneity within a collectively invading tumor cell population. Functional characterization of two phenotype-specific candidate mutations showed that ARP3 enhances collective invasion by promoting the leader cell phenotype and that wild-type KDM5B suppresses chain-like cooperative behavior. These results demonstrate an important role for distinct genetic variants in establishing leader and follower phenotypes and highlight the necessity of maintaining a capacity for phenotypic plasticity during collective cancer invasion.

Analyzing Mechanisms of Metastatic Cancer Cell Adhesive Phenotype Leveraging Preparative Adhesion Chromatography Microfluidic.

An integrated, parallel-plate microfluidic device is engineered to interrogate and fractionate cells based on their adhesivity to a substrate surface functionalized with adhesive ligand in a tightly controlled flow environment to elucidate associated cell-intrinsic pathways. Wall shear stress levels and endothelial presentation of E-selectin are modeled after the inflamed vasculature microenvironment in order to simulate in vitro conditions under which in vivo hematogenous metastasis occurs. Based on elution time from the flow channel, the collection of separate fractions of cells-noninteracting and interacting-at high yields and viabilities enables multiple postperfusion analyses, including flow cytometry, in vivo metastasis modeling, and transcriptomic analysis. This platform enables the interrogation of flow-regulated cell molecular profiles, such as (co)expression levels of natively expressed selectin ligands sLex , CD44, and carcinoembryonic antigen, and cancer stem cell marker CD24. This additionally reveals E-selectin adhesivity exhibited by metastatic human colon carcinoma cells to be a transient phenotype. Facile and rapid, this methodology for unbiased, label free sorting of large populations of cells based on their adhesion in flow represents a method of studying flow-regulated adhesion in vitro for the identification of molecular drug targets for development as antimetastatic cancer therapeutics.

Inhibition of IGF1R enhances 2-deoxyglucose in the treatment of non-small cell lung cancer.

OBJECTIVE: We previously postulated that 2-deoxyglucose (2-DG) activates multiple pro-survival pathways through IGF1R to negate its inhibitory effect on glycolysis. Here, we evaluated whether IGF1R inhibitor synergizes with 2-DG to impede the growth of non-small cell lung cancer (NSCLC). MATERIALS AND METHODS: The activation of IGF1R signaling was assessed by the phosphorylation of IGF1R and its downstream target AKT using immunoblot. Drug dose response and combination index analyses were carried out according to the method of Chou and Talalay. Flow cytometry was used to evaluate cell cycle progression. Apoptosis was monitored by caspase-3/PARP cleavages or Annexin V staining. A subcutaneous xenograft model was used to assess this combination in vivo. RESULTS: 2-DG induces the phosphorylation of IGF1R in its kinase domain, which can be abolished by the IGF1R inhibitor BMS-754807. Furthermore, the combination of 2-DG and BMS-754807 synergistically inhibited the survival of several non-small cell lung cancer (NSCLC) cell lines both in vitro and in vivo. The mechanistic basis of this synergy was cell line-dependent, and LKB1-inactivated EKVX cells underwent apoptosis following treatment with a subtoxic dose of 2-DG and BMS-754807. For these cells, the restoration of LKB1 kinase activity suppressed apoptosis induced by this combination but enhanced G1 arrest. In H460 cells, the addition of 2-DG did not enhance the low level of apoptosis induced by BMS-754807. However, treatment with 0.75 µM of BMS-754807 resulted in the accumulation of H460 cells with 8n-DNA content without affecting cell density increases. Hence, H460 cells may escape BMS-754807-induced G2/M cell cycle arrest through polyploidy. The inclusion of 2-DG blocked formation of the 8n-DNA cell population and restored G2/M phase cell cycle arrest. CONCLUSION: The combination of 2-DG and IGF1R inhibitor BMS-754807 may be used to suppress the proliferation of NSCLC tumors through different mechanisms.

TASI: A software tool for spatial-temporal quantification of tumor spheroid dynamics.

Spheroid cultures derived from explanted cancer specimens are an increasingly utilized resource for studying complex biological processes like tumor cell invasion and metastasis, representing an important bridge between the simplicity and practicality of 2-dimensional monolayer cultures and the complexity and realism of in vivo animal models. Temporal imaging of spheroids can capture the dynamics of cell behaviors and microenvironments, and when combined with quantitative image analysis methods, enables deep interrogation of biological mechanisms. This paper presents a comprehensive open-source software framework for Temporal Analysis of Spheroid Imaging (TASI) that allows investigators to objectively characterize spheroid growth and invasion dynamics. TASI performs spatiotemporal segmentation of spheroid cultures, extraction of features describing spheroid morpho-phenotypes, mathematical modeling of spheroid dynamics, and statistical comparisons of experimental conditions. We demonstrate the utility of this tool in an analysis of non-small cell lung cancer spheroids that exhibit variability in metastatic and proliferative behaviors.

The complex ecosystem in non small cell lung cancer invasion.

Many tumors are characterized by genetic instability, producing an assortment of genetic variants of tumor cells called subclones. These tumors and their surrounding environments form complex multi-cellular ecosystems, where subclones compete for resources and cooperate to perform multiple tasks, including cancer invasion. Our recent empirical studies revealed existence of such distinct phenotypes of cancer cells, leaders and followers, in lung cancer. These two cellular subclones exchange a complex array of extracellular signals demonstrating a symbiotic relationship at the cellular level. Here, we develop a computational model of the microenvironment of the lung cancer ecosystem to explore how the interactions between subclones can advance or inhibit invasion. We found that, due to the complexity of the ecosystem, invasion may have very different dynamics characterized by the different levels of aggressiveness. By altering the signaling environment, we could alter the ecological relationship between the cell types and the overall ecosystem development. Competition between leader and follower cell populations (defined by the limited amount of resources), positive feedback within the leader cell population (controlled by the focal adhesion kinase and fibronectin signaling), and impact of the follower cells to the leaders (represented by yet undetermined proliferation signal) all had major effects on the outcome of the collective dynamics. Specifically, our analysis revealed a class of tumors (defined by the strengths of fibronectin signaling and competition) that are particularly sensitive to manipulations of the signaling environment. These tumors can undergo irreversible changes to the tumor ecosystem that outlast the manipulations of feedbacks and have a profound impact on invasive potential. Our study predicts a complex division of labor between cancer cell subclones and suggests new treatment strategies targeting signaling within the tumor ecosystem.

Vimentin Is Required for Lung Adenocarcinoma Metastasis via Heterotypic Tumor Cell-Cancer-Associated Fibroblast Interactions during Collective Invasion.

Purpose: Vimentin is an epithelial-to-mesenchymal transition (EMT) biomarker and intermediate filament protein that functions during cell migration to maintain structure and motility. Despite the abundance of clinical data linking vimentin to poor patient outcome, it is unclear if vimentin is required for metastasis or is a correlative biomarker. We developed a novel genetically engineered mouse model (GEMM) to probe vimentin in lung adenocarcinoma metastasis.Experimental Design: We used the LSL-KrasG12D/Lkb1fl/fl/Vim-/- model (KLV-/-), which incorporates a whole-body knockout of vimentin and is derived from the Cre-dependent LSL-KrasG12D/Lkb1fl/fl model (KLV+/+). We compared the metastatic phenotypes of the GEMMs and analyzed primary tumors from the KLV models and lung adenocarcinoma patients to assess vimentin expression and function.Results: Characterization of KLV+/+ and KLV-/- mice shows that although vimentin is not required for primary lung tumor growth, vimentin is required for metastasis, and vimentin loss generates lower grade primary tumors. Interestingly, in the KLV+/+ mice, vimentin was not expressed in tumor cells but in cancer-associated fibroblasts (CAFs) surrounding collective invasion packs (CIPs) of epithelial tumor cells, with significantly less CIPs in KLV-/- mice. CIPs correlate with tumor grade and are vimentin-negative and E-cadherin-positive, indicating a lack of cancer cell EMT. A similar heterotypic staining pattern was observed in human lung adenocarcinoma samples. In vitro studies show that vimentin is required for CAF motility to lead tumor cell invasion, supporting a vimentin-dependent model of collective invasion.Conclusions: These data show that vimentin is required for lung adenocarcinoma metastasis by maintaining heterotypic tumor cell-CAF interactions during collective invasion. Clin Cancer Res; 24(2); 420-32. ©2017 AACR.

Modulation of Bax and mTOR for Cancer Therapeutics.

A rationale exists for pharmacologic manipulation of the serine (S)184 phosphorylation site of the proapoptotic Bcl2 family member Bax as an anticancer strategy. Here, we report the refinement of the Bax agonist SMBA1 to generate CYD-2-11, which has characteristics of a suitable clinical lead compound. CYD-2-11 targeted the structural pocket proximal to S184 in the C-terminal region of Bax, directly activating its proapoptotic activity by inducing a conformational change enabling formation of Bax homooligomers in mitochondrial membranes. In murine models of small-cell and non-small cell lung cancers, including patient-derived xenograft and the genetically engineered mutant KRAS-driven lung cancer models, CYD-2-11 suppressed malignant growth without evident significant toxicity to normal tissues. In lung cancer patients treated with mTOR inhibitor RAD001, we observed enhanced S184 Bax phosphorylation in lung cancer cells and tissues that inactivates the propaoptotic function of Bax, contributing to rapalog resistance. Combined treatment of CYD-2-11 and RAD001 in murine lung cancer models displayed strong synergistic activity and overcame rapalog resistance in vitro and in vivo Taken together, our findings provide preclinical evidence for a pharmacologic combination of Bax activation and mTOR inhibition as a rational strategy to improve lung cancer treatment. Cancer Res; 77(11); 3001-12. ©2017 AACR.

Targeting adhesion signaling in KRAS, LKB1 mutant lung adenocarcinoma.

Loss of LKB1 activity is prevalent in KRAS mutant lung adenocarcinoma and promotes aggressive and treatment-resistant tumors. Previous studies have shown that LKB1 is a negative regulator of the focal adhesion kinase (FAK), but in vivo studies testing the efficacy of FAK inhibition in LKB1 mutant cancers are lacking. Here, we took a pharmacologic approach to show that FAK inhibition is an effective early-treatment strategy for this high-risk molecular subtype. We established a lenti-Cre-induced Kras and Lkb1 mutant genetically engineered mouse model (KLLenti) that develops 100% lung adenocarcinoma and showed that high spatiotemporal FAK activation occurs in collective invasive cells that are surrounded by high levels of collagen. Modeling invasion in 3D, loss of Lkb1, but not p53, was sufficient to drive collective invasion and collagen alignment that was highly sensitive to FAK inhibition. Treatment of early, stage-matched KLLenti tumors with FAK inhibitor monotherapy resulted in a striking effect on tumor progression, invasion, and tumor-associated collagen. Chronic treatment extended survival and impeded local lymph node spread. Lastly, we identified focally upregulated FAK and collagen-associated collective invasion in KRAS and LKB1 comutated human lung adenocarcinoma patients. Our results suggest that patients with LKB1 mutant tumors should be stratified for early treatment with FAK inhibitors.

Local alignment vectors reveal cancer cell-induced ECM fiber remodeling dynamics.

Invasive cancer cells interact with the surrounding extracellular matrix (ECM), remodeling ECM fiber network structure by condensing, degrading, and aligning these fibers. We developed a novel local alignment vector analysis method to quantitatively measure collagen fiber alignment as a vector field using Circular Statistics. This method was applied to human non-small cell lung carcinoma (NSCLC) cell lines, embedded as spheroids in a collagen gel. Collagen remodeling was monitored using second harmonic generation imaging under normal conditions and when the LKB1-MARK1 pathway was disrupted through RNAi-based approaches. The results showed that inhibiting LKB1 or MARK1 in NSCLC increases the collagen fiber alignment and captures outward alignment vectors from the tumor spheroid, corresponding to high invasiveness of LKB1 mutant cancer cells. With time-lapse imaging of ECM micro-fiber morphology, the local alignment vector can measure the dynamic signature of invasive cancer cell activity and cell-migration-induced ECM and collagen remodeling and realigning dynamics.

A role for activated Cdc42 in glioblastoma multiforme invasion.

Cdc42 is a Rho-GTPase which plays a major role in regulating cell polarity and migration by specifying the localization of filopodia. However, the role of Cdc42 in GBM invasion has not been thoroughly investigated. We generated stable doxycycline-inducible clones expressing wild type (WT)-, constitutively active (CA)-, and dominant negative (DN)-Cdc42 in three different human glioma cell lines. Expression of CA-Cdc42 significantly increased the migration and invasive properties of malignant glioma cells compared to WT and DN-Cdc42 cell clones, and this was accompanied by a greater number of filopodia and focal adhesion structures which co-localize with phosphorylated focal adhesion kinase (FAK). By mass spectrometry and immunoprecipitation studies, we demonstrated that activated Cdc42 binds to IQGAP1. When implanted orthotopically in mice, the CA-Cdc42 expressing glioma cells exhibited enhanced local migration and invasion, and led to larger tumors, which significantly reduced survival. Using the Cancer Genome Atlas dataset, we determined that high Cdc42 expression is associated with poorer progression free survival, and that Cdc42 expression is highest in the proneural and neural subgroups of GBM. In summary, our studies demonstrate that activated Cdc42 is a critical determinant of the migratory and invasive phenotype of malignant gliomas, and that its effect may be mediated, at least in part, through its interaction with IQGAP1 and phosphorylated FAK.