Phenotypic Plasticity, Bet-Hedging, and Androgen Independence in Prostate Cancer: Role of Non-Genetic Heterogeneity

Phenotypic Plasticity and Cancer: A System Biology Perspective

Epithelial-to-mesenchymal transition (EMT) is a major axis of phenotypic plasticity not only in diseased conditions such as cancer metastasis and fibrosis but also during normal development and wound healing. Yet-another important axis of plasticity with metastatic implications includes the cancer stem cell (CSCs) and non-CSC transitions. However, in both processes, epithelial (E) and mesenchymal (M) phenotypes are not merely binary states. Cancer cells acquire a spectrum of phenotypes with traits, properties, and markers of both E and M phenotypes, giving rise to intermediary hybrid (E/M) phenotypes. E/M cells play an important role in tumor initiation, metastasis, and disease progression in multiple cancers. Furthermore, the hybrid phenotypes also play a major role in causing therapeutic resistance in cancer. Here, we discuss how a systems biology perspective on the problem, which is implicit in the 'Team Medicine' approach outlined in the theme of this Special Issue of The Journal of Clinical Medicine and includes an interdisciplinary team of experts, is more likely to shed new light on EMT in cancer and help us to identify novel therapeutics and strategies to target phenotypic plasticity in cancer.

Mediator kinase inhibition reverses castration resistance of advanced prostate cancer

Mediator kinases CDK19 and CDK8, pleiotropic regulators of transcriptional reprogramming, are differentially regulated by androgen signaling, but both kinases are upregulated in castration-resistant prostate cancer (CRPC). Genetic or pharmacological inhibition of CDK8 and CDK19 reverses the castration-resistant phenotype and restores the sensitivity of CRPC xenografts to androgen deprivation in vivo. Prolonged CDK8/19 inhibitor treatment combined with castration not only suppressed the growth of CRPC xenografts but also induced tumor regression and cures. Transcriptomic analysis revealed that Mediator kinase inhibition amplified and modulated the effects of castration on gene expression, disrupting CRPC adaptation to androgen deprivation. Mediator kinase inactivation in tumor cells also affected stromal gene expression, indicating that Mediator kinase activity in CRPC molded the tumor microenvironment. The combination of castration and Mediator kinase inhibition downregulated the MYC pathway, and Mediator kinase inhibition suppressed a MYC-driven CRPC tumor model even without castration. CDK8/19 inhibitors showed efficacy in patient-derived xenograft models of CRPC, and a gene signature of Mediator kinase activity correlated with tumor progression and overall survival in clinical samples of metastatic CRPC. These results indicate that Mediator kinases mediated androgen-independent in vivo growth of CRPC, supporting the development of CDK8/19 inhibitors for the treatment of this presently incurable disease.

Gompertz models with periodical treatment and applications to prostate cancer

In this paper, Gompertz type models are proposed to understand the temporal tumor volume behavior of prostate cancer when a periodical treatment is provided. Existence, uniqueness, and stability of periodic solutions are established. The models are used to fit the data and to forecast the tumor growth behavior based on prostate cancer treatments using capsaicin and docetaxel anticancer drugs. Numerical simulations show that the combination of capsaicin and docetaxel is the most efficient treatment of prostate cancer.

m6A epitranscriptome analysis reveals differentially methylated transcripts that drive early chemoresistance in bladder cancer

N 6-Methyladenosine (m6A) RNA modifications dynamically regulate messenger RNA processing, differentiation and cell fate. Given these functions, we hypothesized that m6A modifications play a role in the transition to chemoresistance. To test this, we took an agnostic discovery approach anchored directly to chemoresistance rather than to any particular m6A effector protein. Specifically, we used methyl-RNA immunoprecipitation followed by sequencing (MeRIP-seq) in parallel with RNA sequencing to identify gene transcripts that were both differentially methylated and differentially expressed between cisplatin-sensitive and cisplatin-resistant bladder cancer (BC) cells. We filtered and prioritized these genes using clinical and functional database tools, and then validated several of the top candidates via targeted quantitative polymerase chain reaction (qPCR) and MeRIP-PCR. In cisplatin-resistant cells, SLC7A11 transcripts had decreased methylation associated with decreased m6A reader YTHDF3 binding, prolonged RNA stability, and increased RNA and protein levels, leading to reduced ferroptosis and increased survival. Consistent with this, cisplatin-sensitive BC cell lines and patient-derived organoids exposed to cisplatin for as little as 48 h exhibited similar mechanisms of SLC7A11 upregulation and chemoresistance, trends that were also reflected in public cancer survival databases. Collectively, these findings highlight epitranscriptomic plasticity as a mechanism of rapid chemoresistance and a potential therapeutic target.

An integrative analysis of GEO data to identify possible therapeutic biomarkers of prostate cancer and targeting potential protein through Zea mays phytochemicals by virtual screening approaches

Prostate cancer (PC) is a prevalent type of cancer among men. Delaying the treatment of patients with upgraded or upstaged cancer may lead to unmanageable circumstances. The aim of this study is to contribute to the finding of biomarkers that are specific to PC and identify drug candidates derived from plants. The information about cancer is critical for clinicians to make decisions about patient treatment in the era of precision medicine. Advances in genomics technology have opened up new possibilities for identifying genes that are associated with cancer, including PC. This study identifies novel differentially expressed genes for PC. The seven PC microarray datasets were selected from the National Center for Biotechnology Information (NCBI)/Gene Expression Omnibus (GEO). The differentially expressed genes (DEGs) were found based on a fold change of |logFC| ≥ 1 and an adjusted p-value of <0.05. The DEGs were further studied using several bioinformatics tools, including STRING, CytoHubba, SRplot, Coremine Medical database, FunRich and GeneMANIA, cBioPortal. The six new potential biomarkers, GAGE2A, GAGE12G, GAGE2E, GAGE13, GAGE12F and CSAG1 were identified. These biomarkers are associated with biological processes (BPs) such as cell division, and gene expression regulation, so these genes may have a crucial role in PC progression and may serve as potential biomarkers for PC. A total of 497 phytochemicals from corn plants have been screened against the target protein and found LTS0176591 as the best lead molecule with docking score of -6.31 kcal/mol. Further, molecular mechanics-generalized born surface area (MM-GBSA), molecular dynamics simulation, principal component analysis (PCA), free energy landscape (FEL) and molecular mechanics-Poisson-Boltzmann surface area (MM-PBSA) were carried out to validate the findings.Communicated by Ramaswamy H. Sarma.

Dynamical hallmarks of cancer: Phenotypic switching in melanoma and epithelial-mesenchymal plasticity

Phenotypic plasticity was recently incorporated as a hallmark of cancer. This plasticity can manifest along many interconnected axes, such as stemness and differentiation, drug-sensitive and drug-resistant states, and between epithelial and mesenchymal cell-states. Despite growing acceptance for phenotypic plasticity as a hallmark of cancer, the dynamics of this process remains poorly understood. In particular, the knowledge necessary for a predictive understanding of how individual cancer cells and populations of cells dynamically switch their phenotypes in response to the intensity and/or duration of their current and past environmental stimuli remains far from complete. Here, we present recent investigations of phenotypic plasticity from a systems-level perspective using two exemplars: epithelial-mesenchymal plasticity in carcinomas and phenotypic switching in melanoma. We highlight how an integrated computational-experimental approach has helped unravel insights into specific dynamical hallmarks of phenotypic plasticity in different cancers to address the following questions: a) how many distinct cell-states or phenotypes exist?; b) how reversible are transitions among these cell-states, and what factors control the extent of reversibility?; and c) how might cell-cell communication be able to alter rates of cell-state switching and enable diverse patterns of phenotypic heterogeneity? Understanding these dynamic features of phenotypic plasticity may be a key component in shifting the paradigm of cancer treatment from reactionary to a more predictive, proactive approach.

A Nexus between Genetic and Non-Genetic Mechanisms Guides KRAS Inhibitor Resistance in Lung Cancer

Several studies in the last few years have determined that, in contrast to the prevailing dogma that drug resistance is simply due to Darwinian evolution-the selection of mutant clones in response to drug treatment-non-genetic changes can also lead to drug resistance whereby tolerant, reversible phenotypes are eventually relinquished by resistant, irreversible phenotypes. Here, using KRAS as a paradigm, we illustrate how this nexus between genetic and non-genetic mechanisms enables cancer cells to evade the harmful effects of drug treatment. We discuss how the conformational dynamics of the KRAS molecule, that includes intrinsically disordered regions, is influenced by the binding of the targeted therapies contributing to conformational noise and how this noise impacts the interaction of KRAS with partner proteins to rewire the protein interaction network. Thus, in response to drug treatment, reversible drug-tolerant phenotypes emerge via non-genetic mechanisms that eventually enable the emergence of irreversible resistant clones via genetic mutations. Furthermore, we also discuss the recent data demonstrating how combination therapy can help alleviate KRAS drug resistance in lung cancer, and how new treatment strategies based on evolutionary principles may help minimize or even preclude the emergence of drug resistance.

Acquired resistance to KRAS G12C small-molecule inhibitors via genetic/nongenetic mechanisms in lung cancer

Inherent or acquired resistance to sotorasib poses a substantialt challenge for NSCLC treatment. Here, we demonstrate that acquired resistance to sotorasib in isogenic cells correlated with increased expression of integrin β4 (ITGB4), a component of the focal adhesion complex. Silencing ITGB4 in tolerant cells improved sotorasib sensitivity, while overexpressing ITGB4 enhanced tolerance to sotorasib by supporting AKT-mTOR bypass signaling. Chronic treatment with sotorasib induced WNT expression and activated the WNT/β-catenin signaling pathway. Thus, silencing both ITGB4 and β-catenin significantly improved sotorasib sensitivity in tolerant, acquired, and inherently resistant cells. In addition, the proteasome inhibitor carfilzomib (CFZ) exhibited synergism with sotorasib by down-regulating ITGB4 and β-catenin expression. Furthermore, adagrasib phenocopies the combination effect of sotorasib and CFZ by suppressing KRAS activity and inhibiting cell cycle progression in inherently resistant cells. Overall, our findings unveil previously unrecognized nongenetic mechanisms underlying resistance to sotorasib and propose a promising treatment strategy to overcome resistance.

"Bet hedging" against climate change in developing and adult animals: roles for stochastic gene expression, phenotypic plasticity, epigenetic inheritance and adaptation

Animals from embryos to adults experiencing stress from climate change have numerous mechanisms available for enhancing their long-term survival. In this review we consider these options, and how viable they are in a world increasingly experiencing extreme weather associated with climate change. A deeply understood mechanism involves natural selection, leading to evolution of new adaptations that help cope with extreme and stochastic weather events associated with climate change. While potentially effective at staving off environmental challenges, such adaptations typically occur very slowly and incrementally over evolutionary time. Consequently, adaptation through natural selection is in most instances regarded as too slow to aid survival in rapidly changing environments, especially when considering the stochastic nature of extreme weather events associated with climate change. Alternative mechanisms operating in a much shorter time frame than adaptation involve the rapid creation of alternate phenotypes within a life cycle or a few generations. Stochastic gene expression creates multiple phenotypes from the same genotype even in the absence of environmental cues. In contrast, other mechanisms for phenotype change that are externally driven by environmental clues include well-understood developmental phenotypic plasticity (variation, flexibility), which can enable rapid, within-generation changes. Increasingly appreciated are epigenetic influences during development leading to rapid phenotypic changes that can also immediately be very widespread throughout a population, rather than confined to a few individuals as in the case of favorable gene mutations. Such epigenetically-induced phenotypic plasticity can arise rapidly in response to stressors within a generation or across a few generations and just as rapidly be "sunsetted" when the stressor dissipates, providing some capability to withstand environmental stressors emerging from climate change. Importantly, survival mechanisms resulting from adaptations and developmental phenotypic plasticity are not necessarily mutually exclusive, allowing for classic "bet hedging". Thus, the appearance of multiple phenotypes within a single population provides for a phenotype potentially optimal for some future environment. This enhances survival during stochastic extreme weather events associated with climate change. Finally, we end with recommendations for future physiological experiments, recommending in particular that experiments investigating phenotypic flexibility adopt more realistic protocols that reflect the stochastic nature of weather.

Quantifying cancer cell plasticity with gene regulatory networks and single-cell dynamics

Phenotypic plasticity of cancer cells can lead to complex cell state dynamics during tumor progression and acquired resistance. Highly plastic stem-like states may be inherently drug-resistant. Moreover, cell state dynamics in response to therapy allow a tumor to evade treatment. In both scenarios, quantifying plasticity is essential for identifying high-plasticity states or elucidating transition paths between states. Currently, methods to quantify plasticity tend to focus on 1) quantification of quasi-potential based on the underlying gene regulatory network dynamics of the system; or 2) inference of cell potency based on trajectory inference or lineage tracing in single-cell dynamics. Here, we explore both of these approaches and associated computational tools. We then discuss implications of each approach to plasticity metrics, and relevance to cancer treatment strategies.

A Continuation Technique for Maximum Likelihood Estimators in Biological Models

Estimating model parameters is a crucial step in mathematical modelling and typically involves minimizing the disagreement between model predictions and experimental data. This calibration data can change throughout a study, particularly if modelling is performed simultaneously with the calibration experiments, or during an on-going public health crisis as in the case of the COVID-19 pandemic. Consequently, the optimal parameter set, or maximal likelihood estimator (MLE), is a function of the experimental data set. Here, we develop a numerical technique to predict the evolution of the MLE as a function of the experimental data. We show that, when considering perturbations from an initial data set, our approach is significantly more computationally efficient that re-fitting model parameters while producing acceptable model fits to the updated data. We use the continuation technique to develop an explicit functional relationship between fit model parameters and experimental data that can be used to measure the sensitivity of the MLE to experimental data. We then leverage this technique to select between model fits with similar information criteria, a priori determine the experimental measurements to which the MLE is most sensitive, and suggest additional experiment measurements that can resolve parameter uncertainty.

Intrinsically disordered proteins and conformational noise: The hypothesis a decade later

Phenotypic plasticity is the ability of individual genotypes to produce different phenotypes in response to environmental perturbations. We previously postulated how conformational noise emanating from conformational dynamics of intrinsically disordered proteins (IDPs) which is distinct from transcriptional noise, can contribute to phenotypic switching by rewiring the cellular protein interaction network. Since most transcription factors are IDPs, we posited that conformational noise is an integral component of transcriptional noise implying that IDPs may amplify total noise in the system either stochastically or in response to environmental changes. Here, we review progress in elucidating the details of the hypothesis. We highlight empirical evidence supporting the hypothesis, discuss conceptual advances that underscore its fundamental importance and implications, and identify areas for future investigations.

Understanding human aging and the fundamental cell signaling link in age-related diseases: the middle-aging hypovascularity hypoxia hypothesis

Aging-related hypoxia, oxidative stress, and inflammation pathophysiology are closely associated with human age-related carcinogenesis and chronic diseases. However, the connection between hypoxia and hormonal cell signaling pathways is unclear, but such human age-related comorbid diseases do coincide with the middle-aging period of declining sex hormonal signaling. This scoping review evaluates the relevant interdisciplinary evidence to assess the systems biology of function, regulation, and homeostasis in order to discern and decipher the etiology of the connection between hypoxia and hormonal signaling in human age-related comorbid diseases. The hypothesis charts the accumulating evidence to support the development of a hypoxic milieu and oxidative stress-inflammation pathophysiology in middle-aged individuals, as well as the induction of amyloidosis, autophagy, and epithelial-to-mesenchymal transition in aging-related degeneration. Taken together, this new approach and strategy can provide the clarity of concepts and patterns to determine the causes of declining vascularity hemodynamics (blood flow) and physiological oxygenation perfusion (oxygen bioavailability) in relation to oxygen homeostasis and vascularity that cause hypoxia (hypovascularity hypoxia). The middle-aging hypovascularity hypoxia hypothesis could provide the mechanistic interface connecting the endocrine, nitric oxide, and oxygen homeostasis signaling that is closely linked to the progressive conditions of degenerative hypertrophy, atrophy, fibrosis, and neoplasm. An in-depth understanding of these intrinsic biological processes of the developing middle-aged hypoxia could provide potential new strategies for time-dependent therapies in maintaining healthspan for healthy lifestyle aging, medical cost savings, and health system sustainability.

Single-cell transcriptomes in the heart: when every epigenome counts

The response of an organ to stimuli emerges from the actions of individual cells. Recent cardiac single-cell RNA-sequencing studies of development, injury, and reprogramming have uncovered heterogeneous populations even among previously well-defined cell types, raising questions about what level of experimental resolution corresponds to disease-relevant, tissue-level phenotypes. In this review, we explore the biological meaning behind this cellular heterogeneity by undertaking an exhaustive analysis of single-cell transcriptomics in the heart (including a comprehensive, annotated compendium of studies published to date) and evaluating new models for the cardiac function that have emerged from these studies (including discussion and schematics that depict new hypotheses in the field). We evaluate the evidence to support the biological actions of newly identified cell populations and debate questions related to the role of cell-to-cell variability in development and disease. Finally, we present emerging epigenomic approaches that, when combined with single-cell RNA-sequencing, can resolve basic mechanisms of gene regulation and variability in cell phenotype.

Illuminating Intrinsically Disordered Proteins with Integrative Structural Biology

Intense study of intrinsically disordered proteins (IDPs) did not begin in earnest until the late 1990s when a few groups, working independently, convinced the community that these 'weird' proteins could have important functions. Over the past two decades, it has become clear that IDPs play critical roles in a multitude of biological phenomena with prominent examples including coordination in signaling hubs, enabling gene regulation, and regulating ion channels, just to name a few. One contributing factor that delayed appreciation of IDP functional significance is the experimental difficulty in characterizing their dynamic conformations. The combined application of multiple methods, termed integrative structural biology, has emerged as an essential approach to understanding IDP phenomena. Here, we review some of the recent applications of the integrative structural biology philosophy to study IDPs.

Cancer Stem Cells-The Insight into Non-Coding RNAs

Since their initial identification three decades ago, there has been extensive research regarding cancer stem cells (CSCs). It is important to consider the biology of cancer stem cells with a particular focus on their phenotypic and metabolic plasticity, the most important signaling pathways, and non-coding RNAs (ncRNAs) regulating these cellular entities. Furthermore, the current status of therapeutic approaches against CSCs is an important consideration regarding employing the technology to improve human health. Cancer stem cells have claimed to be one of the most important group of cells for the development of several common cancers as they dictate features, such as resistance to radio- and chemotherapy, metastasis, and secondary tumor formation. Therapies which could target these cells may develop into an effective strategy for tumor eradication and a hope for patients for whom this disease remains uncurable.

Microfluidics for long-term single-cell time-lapse microscopy: Advances and applications

Cells are inherently dynamic, whether they are responding to environmental conditions or simply at equilibrium, with biomolecules constantly being made and destroyed. Due to their small volumes, the chemical reactions inside cells are stochastic, such that genetically identical cells display heterogeneous behaviors and gene expression profiles. Studying these dynamic processes is challenging, but the development of microfluidic methods enabling the tracking of individual prokaryotic cells with microscopy over long time periods under controlled growth conditions has led to many discoveries. This review focuses on the recent developments of one such microfluidic device nicknamed the mother machine. We overview the original device design, experimental setup, and challenges associated with this platform. We then describe recent methods for analyzing experiments using automated image segmentation and tracking. We further discuss modifications to the experimental setup that allow for time-varying environmental control, replicating batch culture conditions, cell screening based on their dynamic behaviors, and to accommodate a variety of microbial species. Finally, this review highlights the discoveries enabled by this technology in diverse fields, such as cell-size control, genetic mutations, cellular aging, and synthetic biology.

Addressing Drug Resistance in Cancer: A Team Medicine Approach

Drug resistance remains one of the major impediments to treating cancer. Although many patients respond well initially, resistance to therapy typically ensues. Several confounding factors appear to contribute to this challenge. Here, we first discuss some of the challenges associated with drug resistance. We then discuss how a 'Team Medicine' approach, involving an interdisciplinary team of basic scientists working together with clinicians, has uncovered new therapeutic strategies. These strategies, referred to as intermittent or 'adaptive' therapy, which are based on eco-evolutionary principles, have met with remarkable success in potentially precluding or delaying the emergence of drug resistance in several cancers. Incorporating such treatment strategies into clinical protocols could potentially enhance the precision of delivering personalized medicine to patients. Furthermore, reaching out to patients in the network of hospitals affiliated with leading academic centers could help them benefit from such innovative treatment options. Finally, lowering the dose of the drug and its frequency (because of intermittent rather than continuous therapy) can also have a significant impact on lowering the toxicity and undesirable side effects of the drugs while lowering the financial burden carried by the patient and insurance providers.

Model-informed experimental design recommendations for distinguishing intrinsic and acquired targeted therapeutic resistance in head and neck cancer

The promise of precision medicine has been limited by the pervasive resistance to many targeted therapies for cancer. Inferring the timing (i.e., pre-existing or acquired) and mechanism (i.e., drug-induced) of such resistance is crucial for designing effective new therapeutics. This paper studies cetuximab resistance in head and neck squamous cell carcinoma (HNSCC) using tumor volume data obtained from patient-derived tumor xenografts. We ask if resistance mechanisms can be determined from this data alone, and if not, what data would be needed to deduce the underlying mode(s) of resistance. To answer these questions, we propose a family of mathematical models, with each member of the family assuming a different timing and mechanism of resistance. We present a method for fitting these models to individual volumetric data, and utilize model selection and parameter sensitivity analyses to ask: which member(s) of the family of models best describes HNSCC response to cetuximab, and what does that tell us about the timing and mechanisms driving resistance? We find that along with time-course volumetric data to a single dose of cetuximab, the initial resistance fraction and, in some instances, dose escalation volumetric data are required to distinguish among the family of models and thereby infer the mechanisms of resistance. These findings can inform future experimental design so that we can best leverage the synergy of wet laboratory experimentation and mathematical modeling in the study of novel targeted cancer therapeutics.

Stochastic population dynamics of cancer stemness and adaptive response to therapies

Intratumoral heterogeneity can exist along multiple axes: Cancer stem cells (CSCs)/non-CSCs, drug-sensitive/drug-tolerant states, and a spectrum of epithelial-hybrid-mesenchymal phenotypes. Further, these diverse cell-states can switch reversibly among one another, thereby posing a major challenge to therapeutic efficacy. Therefore, understanding the origins of phenotypic plasticity and heterogeneity remains an active area of investigation. While genomic components (mutations, chromosomal instability) driving heterogeneity have been well-studied, recent reports highlight the role of non-genetic mechanisms in enabling both phenotypic plasticity and heterogeneity. Here, we discuss various processes underlying phenotypic plasticity such as stochastic gene expression, chromatin reprogramming, asymmetric cell division and the presence of multiple stable gene expression patterns ('attractors'). These processes can facilitate a dynamically evolving cell population such that a subpopulation of (drug-tolerant) cells can survive lethal drug exposure and recapitulate population heterogeneity on drug withdrawal, leading to relapse. These drug-tolerant cells can be both pre-existing and also induced by the drug itself through cell-state reprogramming. The dynamics of cell-state transitions both in absence and presence of the drug can be quantified through mathematical models. Such a dynamical systems approach to elucidating patterns of intratumoral heterogeneity by integrating longitudinal experimental data with mathematical models can help design effective combinatorial and/or sequential therapies for better clinical outcomes.

Neuropilin-2 axis in regulating secretory phenotype of neuroendocrine-like prostate cancer cells and its implication in therapy resistance

Neuroendocrine (NE)-like tumors secrete various signaling molecules to establish paracrine communication within the tumor milieu and to create a therapy-resistant environment. It is important to identify molecular mediators that regulate this secretory phenotype in NE-like cancer. The current study highlights the importance of a cell surface molecule, Neuropilin-2 (NRP2), for the secretory function of NE-like prostate cancer (PCa). Our analysis on different patient cohorts suggests that NRP2 is high in NE-like PCa. We have developed cell line models to investigate NRP2's role in NE-like PCa. Our bioinformatics, mass spectrometry, cytokine array, and other supporting experiments reveal that NRP2 regulates robust secretory phenotype in NE-like PCa and controls the secretion of factors promoting cancer cell survival. Depletion of NRP2 reduces the secretion of these factors and makes resistant cancer cells sensitive to chemotherapy in vitro and in vivo. Therefore, targeting NRP2 can revert cellular secretion and sensitize PCa cells toward therapy.

From developmental to atavistic bet-hedging: How cancer cells pervert the exploitation of random single-cell phenotypic fluctuations

Stochastic gene expression plays a leading developmental role through its contribution to cell differentiation. It is also proposed to promote phenotypic diversification in malignant cells. However, it remains unclear if these two forms of cellular bet-hedging are identical or rather display distinct features. Here we argue that bet-hedging phenomena in cancer cells are more similar to those occurring in unicellular organisms than to those of normal metazoan cells. We further propose that the atavistic bet-hedging strategies in cancer originate from a hijacking of the normal developmental bet-hedging of metazoans. Finally, we discuss the constraints that may shape the atavistic bet-hedging strategies of cancer cells.

Emergence of hybrid states of stem-like cancer cells correlates with poor prognosis in oral cancer

Cancer cell state transitions emerged as powerful mechanisms responsible for drug tolerance and overall poor prognosis; however, evidences were largely missing in oral cancer. Here, by multiplexing phenotypic markers of stem-like cancer cells (SLCCs); CD44, CD24 and aldehyde dehydrogenase (ALDH), we characterized diversity among multiple oral tumor tissues and cell lines. Two distinct patterns of spontaneous transitions with stochastic bidirectional interconversions on ‘ALDH-axis’, and unidirectional non-interconvertible transitions on ‘CD24-axis’ were observed. Interestingly, plastic ‘ALDH-axis’ was harnessed by cells to adapt to a Cisplatin tolerant state. Furthermore, phenotype-specific RNA sequencing suggested the possible maintenance of intermediate hybrid cell states maintaining stemness within the differentiating subpopulations. Importantly, survival analysis with subpopulation-specific gene sets strongly suggested that cell-state transitions may drive non-genetic heterogeneity, resulting in poor prognosis. Therefore, we have described the phenotypic-composition of heterogeneous subpopulations critical for global tumor behavior in oral cancer; which may provide prerequisite knowledge for treatment strategies.

•

Demonstrated population trajectory driven non-genetic heterogeneity in oral cancer

•

Created transition maps for subpopulations using discrete time Markov chain model

•

Demonstrated maintenance of stemness in cells undergoing differentiation

•

Uniquely expressed genes of these subpopulations associated with disease prognosis

Tumour Stem Cells in Breast Cancer

Tumour stem cells (CSCs) are a self-renewing population that plays important roles in tumour initiation, recurrence, and metastasis. Although the medical literature is extensive, problems with CSC identification and cancer therapy remain. This review provides the main mechanisms of CSC action in breast cancer (BC): CSC markers and signalling pathways, heterogeneity, plasticity, and ecological behaviour. The dynamic heterogeneity of CSCs and the dynamic transitions of CSC⁻ non-CSCs and their significance for metastasis are considered.

EMT, Stemness, and Drug Resistance in Biological Context: A 3D Tumor Tissue/In Silico Platform for Analysis of Combinatorial Treatment in NSCLC with Aggressive KRAS-Biomarker Signatures

Simple Summary

The phenotypic transition of tumor cells from epithelial to mesenchymal characteristics is called EMT and is widely discussed in the scientific community as a game changer in drug resistance and metastasis formation. However, clinical studies could not prove the efficacy of EMT-interfering treatments, and in clinical routine, EMT is not investigated to assess invasion. To fill this gap between bench and bedside, we use in this study a lung tumor tissue model with a preserved basement membrane for investigation of EMT functions with respect to invasion across this membrane and drug resistance. Our results suggest EMT is more a marker of drug resistance than a maker. Invasion is enhanced by EMT but more dependent on intrinsic factors, and EMT is not detected in the center of invasive tumor nodules. An in silico signaling network model is used to integrate these in vitro results and to reveal determinants for drug response.

Abstract

Epithelial-to-mesenchymal transition (EMT) is discussed to be centrally involved in invasion, stemness, and drug resistance. Experimental models to evaluate this process in its biological complexity are limited. To shed light on EMT impact and test drug response more reliably, we use a lung tumor test system based on a decellularized intestinal matrix showing more in vivo-like proliferation levels and enhanced expression of clinical markers and carcinogenesis-related genes. In our models, we found evidence for a correlation of EMT with drug resistance in primary and secondary resistant cells harboring KRAS^G12C or EGFR mutations, which was simulated in silico based on an optimized signaling network topology. Notably, drug resistance did not correlate with EMT status in KRAS-mutated patient-derived xenograft (PDX) cell lines, and drug efficacy was not affected by EMT induction via TGF-β. To investigate further determinants of drug response, we tested several drugs in combination with a KRAS^G12C inhibitor in KRAS^G12C mutant HCC44 models, which, besides EMT, display mutations in P53, LKB1, KEAP1, and high c-MYC expression. We identified an aurora-kinase A (AURKA) inhibitor as the most promising candidate. In our network, AURKA is a centrally linked hub to EMT, proliferation, apoptosis, LKB1, and c-MYC. This exemplifies our systemic analysis approach for clinical translation of biomarker signatures.

Bladder cancer cells shift rapidly and spontaneously to cisplatin-resistant oxidative phosphorylation that is trackable in real time

Genetic mutations have long been recognized as drivers of cancer drug resistance, but recent work has defined additional non-genetic mechanisms of plasticity, wherein cancer cells assume a drug resistant phenotype marked by altered epigenetic and transcriptional states. Currently, little is known about the real-time, dynamic nature of this phenotypic shift. Using a bladder cancer model of nongenetic plasticity, we discovered that rapid transition to drug resistance entails upregulation of mitochondrial gene expression and a corresponding metabolic shift towards the tricarboxylic acid cycle and oxidative phosphorylation. Based on this distinction, we were able to track cancer cell metabolic profiles in real time using fluorescence lifetime microscopy (FLIM). We observed single cells transitioning spontaneously to an oxidative phosphorylation state over hours to days, a trend that intensified with exposure to cisplatin chemotherapy. Conversely, pharmacological inhibition of oxidative phosphorylation significantly reversed the FLIM metabolic signature and reduced cisplatin resistance. These rapid, spontaneous metabolic shifts offer a new means of tracking nongenetic cancer plasticity and forestalling the emergence of drug resistance.

Quantification of long-term doxorubicin response dynamics in breast cancer cell lines to direct treatment schedules

While acquired chemoresistance is recognized as a key challenge to treating many types of cancer, the dynamics with which drug sensitivity changes after exposure are poorly characterized. Most chemotherapeutic regimens call for repeated dosing at regular intervals, and if drug sensitivity changes on a similar time scale then the treatment interval could be optimized to improve treatment performance. Theoretical work suggests that such optimal schedules exist, but experimental confirmation has been obstructed by the difficulty of deconvolving the simultaneous processes of death, adaptation, and regrowth taking place in cancer cell populations. Here we present a method of optimizing drug schedules in vitro through iterative application of experimentally calibrated models, and demonstrate its ability to characterize dynamic changes in sensitivity to the chemotherapeutic doxorubicin in three breast cancer cell lines subjected to treatment schedules varying in concentration, interval between pulse treatments, and number of sequential pulse treatments. Cell populations are monitored longitudinally through automated imaging for 600–800 hours, and this data is used to calibrate a family of cancer growth models, each consisting of a system of ordinary differential equations, derived from the bi-exponential model which characterizes resistant and sensitive subpopulations. We identify a model incorporating both a period of growth arrest in surviving cells and a delay in the death of chemosensitive cells which outperforms the original bi-exponential growth model in Akaike Information Criterion based model selection, and use the calibrated model to quantify the performance of each drug schedule. We find that the inter-treatment interval is a key variable in determining the performance of sequential dosing schedules and identify an optimal retreatment time for each cell line which extends regrowth time by 40%-239%, demonstrating that the time scale of changes in chemosensitivity following doxorubicin exposure allows optimization of drug scheduling by varying this inter-treatment interval.

Acquired chemoresistance is a common cause of treatment failure in cancer. The scheduling of a multi-dose course of chemotherapeutic treatment may influence the dynamics of acquired chemoresistance, and drug schedule optimization may increase the duration of effectiveness of a particular chemotherapeutic agent for a particular patient. Here we present a method for experimentally optimizing an in vitro drug schedule through iterative rounds of experimentation and computational analysis, and demonstrate the method’s ability to improve the performance of doxorubicin treatment in three breast carcinoma cell lines. Specifically, we find that the interval between drug exposures can be optimized while holding drug concentration and number of treatments constant, suggesting that this may be a key variable to explore in future drug schedule optimization efforts. We further use this method’s model calibration and selection process to extract information about the underlying biology of the doxorubicin response, and find that the incorporation of delays on both cell death and regrowth are necessary for accurate parameterization of cell growth data.

Intrinsically disordered proteins: Ensembles at the limits of Anfinsen's dogma

Intrinsically disordered proteins (IDPs) are proteins that lack rigid 3D structure. Hence, they are often misconceived to present a challenge to Anfinsen's dogma. However, IDPs exist as ensembles that sample a quasi-continuum of rapidly interconverting conformations and, as such, may represent proteins at the extreme limit of the Anfinsen postulate. IDPs play important biological roles and are key components of the cellular protein interaction network (PIN). Many IDPs can interconvert between disordered and ordered states as they bind to appropriate partners. Conformational dynamics of IDPs contribute to conformational noise in the cell. Thus, the dysregulation of IDPs contributes to increased noise and "promiscuous" interactions. This leads to PIN rewiring to output an appropriate response underscoring the critical role of IDPs in cellular decision making. Nonetheless, IDPs are not easily tractable experimentally. Furthermore, in the absence of a reference conformation, discerning the energy landscape representation of the weakly funneled IDPs in terms of reaction coordinates is challenging. To understand conformational dynamics in real time and decipher how IDPs recognize multiple binding partners with high specificity, several sophisticated knowledge-based and physics-based in silico sampling techniques have been developed. Here, using specific examples, we highlight recent advances in energy landscape visualization and molecular dynamics simulations to discern conformational dynamics and discuss how the conformational preferences of IDPs modulate their function, especially in phenotypic switching. Finally, we discuss recent progress in identifying small molecules targeting IDPs underscoring the potential therapeutic value of IDPs. Understanding structure and function of IDPs can not only provide new insight on cellular decision making but may also help to refine and extend Anfinsen's structure/function paradigm.

Population Dynamics of Epithelial-Mesenchymal Heterogeneity in Cancer Cells

Phenotypic heterogeneity is a hallmark of aggressive cancer behaviour and a clinical challenge. Despite much characterisation of this heterogeneity at a multi-omics level in many cancers, we have a limited understanding of how this heterogeneity emerges spontaneously in an isogenic cell population. Some longitudinal observations of dynamics in epithelial-mesenchymal heterogeneity, a canonical example of phenotypic heterogeneity, have offered us opportunities to quantify the rates of phenotypic switching that may drive such heterogeneity. Here, we offer a mathematical modeling framework that explains the salient features of population dynamics noted in PMC42-LA cells: (a) predominance of EpCAM^high subpopulation, (b) re-establishment of parental distributions from the EpCAM^high and EpCAM^low subpopulations, and (c) enhanced heterogeneity in clonal populations established from individual cells. Our framework proposes that fluctuations or noise in content duplication and partitioning of SNAIL—an EMT-inducing transcription factor—during cell division can explain spontaneous phenotypic switching and consequent dynamic heterogeneity in PMC42-LA cells observed experimentally at both single-cell and bulk level analysis. Together, we propose that asymmetric cell division can be a potential mechanism for phenotypic heterogeneity.

Roadmap on plasticity and epigenetics in cancer

The role of plasticity and epigenetics in shaping cancer evolution and response to therapy has taken center stage with recent technological advances including single cell sequencing. This roadmap article is focused on state-of-the-art mathematical and experimental approaches to interrogate plasticity in cancer, and addresses the following themes and questions: is there a formal overarching framework that encompasses both non-genetic plasticity and mutation-driven somatic evolution? How do we measure and model the role of the microenvironment in influencing/controlling non-genetic plasticity? How can we experimentally study non-genetic plasticity? Which mathematical techniques are required or best suited? What are the clinical and practical applications and implications of these concepts?

Find the path of least resistance: Adaptive therapy to delay treatment failure and improve outcomes

Cytotoxic chemotherapy and targeted therapies help people with advanced cancers, but for most, treatment fails. Cancer heterogeneity is one cause of treatment failure, but also suggests an opportunity to improve outcomes; reconceptualising cancer therapy as an ecological problem offers the strategy of adaptive therapy. If an agent is active against a patient's cancer, instead of traditional continuous dosing at the maximum tolerated dose until treatment failure, the patient and their oncologist may instead choose to pause treatment as soon as the cancer responds. When tumour burden increases, the cancer is rechallenged with the same agent in hope of delivering another response, ideally before symptoms occur or quality-of-life is impacted. These 'loops' of 'pause/restart' allows an active treatment to be used strategically, to delay the development of evolutionary selection within the cancer, delaying the onset of treatment resistance, controlling the cancer for longer. Modelling predicts patients can navigate several 'loops', potentially increasing the utility of an active treatment by multiples, and early trials suggest at least doubling of progression-free survival. In this narrative review we confront how cancer heterogeneity limits treatment effectiveness, re-examine cancer as an ecological problem, review the data supporting adaptive therapy and outline the challenges and opportunities faced in clinical practice to implement this evolutionary concept. In an era where multiple novel active anti-neoplastic agents are being used with ancient inflexibile maximum tolerated dose for maximum duration approaches, adaptive dosing offers a personalised, n = 1 approach to cancer therapy selection.

Synthesis and preliminary evaluation of octreotate conjugates of bioactive synthetic amatoxins for targeting somatostatin receptor (sstr2) expressing cells

Targeted cancer therapy represents a paradigm-shifting approach that aims to deliver a toxic payload selectively to target-expressing cells thereby sparing normal tissues the off-target effects associated with traditional chemotherapeutics. Since most targeted constructs rely on standard microtubule inhibitors or DNA-reactive molecules as payloads, new toxins that inhibit other intracellular targets are needed to realize the full potential of targeted therapy. Among these new payloads, α-amanitin has gained attraction as a payload in targeted therapy. Here, we conjugate two synthetic amanitins at different sites to demonstrate their utility as payloads in peptide drug conjugates (PDCs). As an exemplary targeting agent, we chose octreotate, a well-studied somatostatin receptor (sstr2) peptide agonist for the conjugation to synthetic amatoxins via three tailor-built linkers. The linker chemistry permitted the evaluation of one non-cleavable and two cleavable self-immolative conjugates. The immolating linkers were chosen to take advantage of either the reducing potential of the intracellular environment or the high levels of lysosomal proteases in tumor cells to trigger toxin release. Cell-based assays on target-positive Ar42J cells revealed target-specific reduction in viability with up to 1000-fold enhancement in bioactivity compared to the untargeted amatoxins. Altogether, this preliminary study enabled the development of a highly modular synthetic platform for the construction of amanitin-based conjugates that can be readily extended to various targeting moieties.

Transition Therapy: Tackling the Ecology of Tumor Phenotypic Plasticity

Phenotypic switching in cancer cells has been found to be present across tumor types. Recent studies on Glioblastoma report a remarkably common architecture of four well-defined phenotypes coexisting within high levels of intra-tumor genetic heterogeneity. Similar dynamics have been shown to occur in breast cancer and melanoma and are likely to be found across cancer types. Given the adaptive potential of phenotypic switching (PHS) strategies, understanding how it drives tumor evolution and therapy resistance is a major priority. Here we present a mathematical framework uncovering the ecological dynamics behind PHS. The model is able to reproduce experimental results, and mathematical conditions for cancer progression reveal PHS-specific features of tumors with direct consequences on therapy resistance. In particular, our model reveals a threshold for the resistant-to-sensitive phenotype transition rate, below which any cytotoxic or switch-inhibition therapy is likely to fail. The model is able to capture therapeutic success thresholds for cancers where nonlinear growth dynamics or larger PHS architectures are in place, such as glioblastoma or melanoma. By doing so, the model presents a novel set of conditions for the success of combination therapies able to target replication and phenotypic transitions at once. Following our results, we discuss transition therapy as a novel scheme to target not only combined cytotoxicity but also the rates of phenotypic switching.

Dynamic Phenotypic Switching and Group Behavior Help Non-Small Cell Lung Cancer Cells Evade Chemotherapy

Drug resistance, a major challenge in cancer therapy, is typically attributed to mutations and genetic heterogeneity. Emerging evidence suggests that dynamic cellular interactions and group behavior also contribute to drug resistance. However, the underlying mechanisms remain poorly understood. Here, we present a new mathematical approach with game theoretical underpinnings that we developed to model real-time growth data of non-small cell lung cancer (NSCLC) cells and discern patterns in response to treatment with cisplatin. We show that the cisplatin-sensitive and cisplatin-tolerant NSCLC cells, when co-cultured in the absence or presence of the drug, display dynamic group behavior strategies. Tolerant cells exhibit a 'persister-like' behavior and are attenuated by sensitive cells; they also appear to 'educate' sensitive cells to evade chemotherapy. Further, tolerant cells can switch phenotypes to become sensitive, especially at low cisplatin concentrations. Finally, switching treatment from continuous to an intermittent regimen can attenuate the emergence of tolerant cells, suggesting that intermittent chemotherapy may improve outcomes in lung cancer.

Single-Cell Technologies to Study Phenotypic Heterogeneity and Bacterial Persisters

Antibiotic persistence is a phenomenon in which rare cells of a clonal bacterial population can survive antibiotic doses that kill their kin, even though the entire population is genetically susceptible. With antibiotic treatment failure on the rise, there is growing interest in understanding the molecular mechanisms underlying bacterial phenotypic heterogeneity and antibiotic persistence. However, elucidating these rare cell states can be technically challenging. The advent of single-cell techniques has enabled us to observe and quantitatively investigate individual cells in complex, phenotypically heterogeneous populations. In this review, we will discuss current technologies for studying persister phenotypes, including fluorescent tags and biosensors used to elucidate cellular processes; advances in flow cytometry, mass spectrometry, Raman spectroscopy, and microfluidics that contribute high-throughput and high-content information; and next-generation sequencing for powerful insights into genetic and transcriptomic programs. We will further discuss existing knowledge gaps, cutting-edge technologies that can address them, and how advances in single-cell microbiology can potentially improve infectious disease treatment outcomes.

Bridging Tumorigenesis and Therapy Resistance With a Non-Darwinian and Non-Lamarckian Mechanism of Adaptive Evolution

Although neo-Darwinian (and less often Lamarckian) dynamics are regularly invoked to interpret cancer's multifarious molecular profiles, they shine little light on how tumorigenesis unfolds and often fail to fully capture the frequency and breadth of resistance mechanisms. This uncertainty frames one of the most problematic gaps between science and practice in modern times. Here, we offer a theory of adaptive cancer evolution, which builds on a molecular mechanism that lies outside neo-Darwinian and Lamarckian schemes. This mechanism coherently integrates non-genetic and genetic changes, ecological and evolutionary time scales, and shifts the spotlight away from positive selection towards purifying selection, genetic drift, and the creative-disruptive power of environmental change. The surprisingly simple use-it or lose-it rationale of the proposed theory can help predict molecular dynamics during tumorigenesis. It also provides simple rules of thumb that should help improve therapeutic approaches in cancer.

Collective metastasis: coordinating the multicellular voyage

The metastatic process is arduous. Cancer cells must escape the confines of the primary tumor, make their way into and travel through the circulation, then survive and proliferate in unfavorable microenvironments. A key question is how cancer cells overcome these multiple barriers to orchestrate distant organ colonization. Accumulating evidence in human patients and animal models supports the hypothesis that clusters of tumor cells can complete the entire metastatic journey in a process referred to as collective metastasis. Here we highlight recent studies unraveling how multicellular coordination, via both physical and biochemical coupling of cells, induces cooperative properties advantageous for the completion of metastasis. We discuss conceptual challenges and unique mechanisms arising from collective dissemination that are distinct from single cell-based metastasis. Finally, we consider how the dissection of molecular transitions regulating collective metastasis could offer potential insight into cancer therapy.

The role of memory in non-genetic inheritance and its impact on cancer treatment resistance

Intra-tumour heterogeneity is a leading cause of treatment failure and disease progression in cancer. While genetic mutations have long been accepted as a primary mechanism of generating this heterogeneity, the role of phenotypic plasticity is becoming increasingly apparent as a driver of intra-tumour heterogeneity. Consequently, understanding the role of this plasticity in treatment resistance and failure is a key component of improving cancer therapy. We develop a mathematical model of stochastic phenotype switching that tracks the evolution of drug-sensitive and drug-tolerant subpopulations to clarify the role of phenotype switching on population growth rates and tumour persistence. By including cytotoxic therapy in the model, we show that, depending on the strategy of the drug-tolerant subpopulation, stochastic phenotype switching can lead to either transient or permanent drug resistance. We study the role of phenotypic heterogeneity in a drug-resistant, genetically homogeneous population of non-small cell lung cancer cells to derive a rational treatment schedule that drives population extinction and avoids competitive release of the drug-tolerant sub-population. This model-informed therapeutic schedule results in increased treatment efficacy when compared against periodic therapy, and, most importantly, sustained tumour decay without the development of resistance.

Lineage Plasticity in Cancer: The Tale of a Skin-Walker

Lineage plasticity, the switching of cells from one lineage to another, has been recognized as a cardinal property essential for embryonic development, tissue repair and homeostasis. However, such a highly regulated process goes awry when cancer cells exploit this inherent ability to their advantage, resulting in tumorigenesis, relapse, metastasis and therapy resistance. In this review, we summarize our current understanding on the role of lineage plasticity in tumor progression and therapeutic resistance in multiple cancers. Lineage plasticity can be triggered by treatment itself and is reported across various solid as well as liquid tumors. Here, we focus on the importance of lineage switching in tumor progression and therapeutic resistance of solid tumors such as the prostate, lung, hepatocellular and colorectal carcinoma and the myeloid and lymphoid lineage switch observed in leukemias. Besides this, we also discuss the role of epithelial-mesenchymal transition (EMT) in facilitating the lineage switch in biphasic cancers such as aggressive carcinosarcomas. We also discuss the mechanisms involved, current therapeutic approaches and challenges that lie ahead in taming the scourge of lineage plasticity in cancer.

A mechanistic model captures the emergence and implications of non-genetic heterogeneity and reversible drug resistance in ER+ breast cancer cells

Resistance to anti-estrogen therapy is an unsolved clinical challenge in successfully treating ER+ breast cancer patients. Recent studies have demonstrated the role of non-genetic (i.e. phenotypic) adaptations in tolerating drug treatments; however, the mechanisms and dynamics of such non-genetic adaptation remain elusive. Here, we investigate coupled dynamics of epithelial–mesenchymal transition (EMT) in breast cancer cells and emergence of reversible drug resistance. Our mechanism-based model for underlying regulatory network reveals that these two axes can drive one another, thus enabling non-genetic heterogeneity in a cell population by allowing for six co-existing phenotypes: epithelial-sensitive, mesenchymal-resistant, hybrid E/M-sensitive, hybrid E/M-resistant, mesenchymal-sensitive and epithelial-resistant, with the first two ones being most dominant. Next, in a population dynamics framework, we exemplify the implications of phenotypic plasticity (both drug-induced and intrinsic stochastic switching) and/or non-genetic heterogeneity in promoting population survival in a mixture of sensitive and resistant cells, even in the absence of any cell–cell cooperation. Finally, we propose the potential therapeutic use of mesenchymal–epithelial transition inducers besides canonical anti-estrogen therapy to limit the emergence of reversible drug resistance. Our results offer mechanistic insights into empirical observations on EMT and drug resistance and illustrate how such dynamical insights can be exploited for better therapeutic designs.

Aspects of the Tumor Microenvironment Involved in Immune Resistance and Drug Resistance

The tumor microenvironment (TME) is a complex and ever-changing "rogue organ" composed of its own blood supply, lymphatic and nervous systems, stroma, immune cells and extracellular matrix (ECM). These complex components, utilizing both benign and malignant cells, nurture the harsh, immunosuppressive and nutrient-deficient environment necessary for tumor cell growth, proliferation and phenotypic flexibility and variation. An important aspect of the TME is cellular crosstalk and cell-to-ECM communication. This interaction induces the release of soluble factors responsible for immune evasion and ECM remodeling, which further contribute to therapy resistance. Other aspects are the presence of exosomes contributed by both malignant and benign cells, circulating deregulated microRNAs and TME-specific metabolic patterns which further potentiate the progression and/or resistance to therapy. In addition to biochemical signaling, specific TME characteristics such as the hypoxic environment, metabolic derangements, and abnormal mechanical forces have been implicated in the development of treatment resistance. In this review, we will provide an overview of tumor microenvironmental composition, structure, and features that influence immune suppression and contribute to treatment resistance.

Group Behavior and Emergence of Cancer Drug Resistance

Drug resistance is a major impediment in cancer. Although it is generally thought that acquired drug resistance is due to genetic mutations, emerging evidence indicates that nongenetic mechanisms also play an important role. Resistance emerges through a complex interplay of clonal groups within a heterogeneous tumor and the surrounding microenvironment. Traits such as phenotypic plasticity, intercellular communication, and adaptive stress response, act in concert to ensure survival of intermediate reversible phenotypes, until permanent, resistant clones can emerge. Understanding the role of group behavior, and the underlying nongenetic mechanisms, can lead to more efficacious treatment designs and minimize or delay emergence of resistance.

Topological signatures in regulatory network enable phenotypic heterogeneity in small cell lung cancer

Phenotypic (non-genetic) heterogeneity has significant implications for the development and evolution of organs, organisms, and populations. Recent observations in multiple cancers have unraveled the role of phenotypic heterogeneity in driving metastasis and therapy recalcitrance. However, the origins of such phenotypic heterogeneity are poorly understood in most cancers. Here, we investigate a regulatory network underlying phenotypic heterogeneity in small cell lung cancer, a devastating disease with no molecular targeted therapy. Discrete and continuous dynamical simulations of this network reveal its multistable behavior that can explain co-existence of four experimentally observed phenotypes. Analysis of the network topology uncovers that multistability emerges from two teams of players that mutually inhibit each other, but members of a team activate one another, forming a 'toggle switch' between the two teams. Deciphering these topological signatures in cancer-related regulatory networks can unravel their 'latent' design principles and offer a rational approach to characterize phenotypic heterogeneity in a tumor.

Coupled Feedback Loops Involving PAGE4, EMT and Notch Signaling Can Give Rise to Non-genetic Heterogeneity in Prostate Cancer Cells

Non-genetic heterogeneity is emerging as a crucial factor underlying therapy resistance in multiple cancers. However, the design principles of regulatory networks underlying non-genetic heterogeneity in cancer remain poorly understood. Here, we investigate the coupled dynamics of feedback loops involving (a) oscillations in androgen receptor (AR) signaling mediated through an intrinsically disordered protein PAGE4, (b) multistability in epithelial-mesenchymal transition (EMT), and c) Notch-Delta-Jagged signaling mediated cell-cell communication, each of which can generate non-genetic heterogeneity through multistability and/or oscillations. Our results show how different coupling strengths between AR and EMT signaling can lead to monostability, bistability, or oscillations in the levels of AR, as well as propagation of oscillations to EMT dynamics. These results reveal the emergent dynamics of coupled oscillatory and multi-stable systems and unravel mechanisms by which non-genetic heterogeneity in AR levels can be generated, which can act as a barrier to most existing therapies for prostate cancer patients.

Cell memory of epithelial-mesenchymal plasticity in cancer

Fundamental biological processes of cell identity and cell fate determination are controlled by complex regulatory networks. These processes require molecular mechanisms that confer cellular phenotypic memory and state persistence. In this minireview, we will summarize mechanisms of cell memory based on regulatory hysteretic feedback loops and explore epigenetic mechanisms widely represented in nature, with special focus on epithelial-to-mesenchymal plasticity. We will also discuss the functional consequences of cell memory and epithelial-to-mesenchymal plasticity dynamics during development and cancer metastasis.

Multi-Stability and Consequent Phenotypic Plasticity in AMPK-Akt Double Negative Feedback Loop in Cancer Cells

Adaptation and survival of cancer cells to various stress and growth factor conditions is crucial for successful metastasis. A double-negative feedback loop between two serine/threonine kinases AMPK (AMP-activated protein kinase) and Akt can regulate the adaptation of breast cancer cells to matrix-deprivation stress. This feedback loop can significantly generate two phenotypes or cell states: matrix detachment-triggered pAMPK^high/ pAkt^low state, and matrix (re)attachment-triggered pAkt^high/ pAMPK^low state. However, whether these two cell states can exhibit phenotypic plasticity and heterogeneity in a given cell population, i.e., whether they can co-exist and undergo spontaneous switching to generate the other subpopulation, remains unclear. Here, we develop a mechanism-based mathematical model that captures the set of experimentally reported interactions among AMPK and Akt. Our simulations suggest that the AMPK-Akt feedback loop can give rise to two co-existing phenotypes (pAkt^high/ pAMPK^low and pAMPK^high/pAkt^low) in specific parameter regimes. Next, to test the model predictions, we segregated these two subpopulations in MDA-MB-231 cells and observed that each of them was capable of switching to another in adherent conditions. Finally, the predicted trends are supported by clinical data analysis of The Cancer Genome Atlas (TCGA) breast cancer and pan-cancer cohorts that revealed negatively correlated pAMPK and pAkt protein levels. Overall, our integrated computational-experimental approach unravels that AMPK-Akt feedback loop can generate multi-stability and drive phenotypic switching and heterogeneity in a cancer cell population.