Drug GRADE: An Integrated Analysis of Population Growth and Cell Death Reveals Drug-Specific and Cancer Subtype-Specific Response Profiles

High content imaging of relative ATP levels for mitochondrial toxicity prediction in human induced pluripotent stem cell derived cardiomyocytes

Human induced pluripotent stem cell derived cardiomyocytes (hiPSC-CMs) are increasingly being evaluated in assays aimed at better understanding potential cardiotoxic liability of newly developed therapeutic compounds. Disruption of mitochondria has been implicated in the mechanism of drug-induced cardiotoxicity of some compounds. Therefore, we have developed a high content imaging assay for the investigation of mitochondrial toxicity in hiPSC-CMs using ATP-Red, a fluorescent dye capable of detecting subcellular localization of relative ATP levels in living cells. We demonstrated time-dependent decreases in ATP-Red signal over 6h treatment with known mitochondrial toxicants antimycin (0.03, 0.1µM) or oligomycin (3, 10µM). Concentration-dependent decreases in ATP-Red signal with antimycin (0.001-0.3µM) and oligomycin (0.003-1µM) were rescued by restoring glycolysis through glucose supplementation. Decreased ATP levels were also identified in a selection of clinically available drugs with reported association with mitochondrial toxicity but absent in compounds with no known association with mitochondrial dysfunction. ATP measurements using the ATP-Red imaging assay were consistent with orthogonal measurements of whole cell ATP levels in lysed hiPSC-CMs following compound treatment. Similar findings were also obtained with measurement of mitochondrial membrane potential, except amiodarone which had no change despite decreased ATP levels. The developed high throughput imaging assay, assessing subcellular ATP dynamics, could provide mechanistic insights for tested compounds.

Pol II degradation activates cell death independently from the loss of transcription

Pol II-mediated transcription is essential for eukaryotic life. While loss of transcription is thought to be universally lethal, the associated mechanisms promoting cell death are not yet known. Here, we show that death following loss of Pol II is not caused by dysregulated gene expression. Instead, death occurs in response to the loss of Pol II protein itself, specifically loss of the enzymatic subunit, Rbp1. Loss of Pol II exclusively activates apoptosis, and expression of a transcriptionally inactive version of Rpb1 rescues cell viability. Using functional genomics, we identify a previously uncharacterized mechanism that regulates lethality following loss of Pol II, which we call the Pol II Degradation-dependent Apoptotic Response (PDAR). Using the genetic dependencies of PDAR, we identify clinically used drugs that owe their efficacy to a PDAR-dependent mechanism. Our findings unveil a novel apoptotic signaling response that contributes to the efficacy of a wide array of anti-cancer therapies.

Inhibitory effects of flavonoids on organic cation transporter 1: Implications for food/herb-drug interactions and hepatoprotective effects

Organic cation transporter 1 (OCT1, gene symbol: SLC22A1) is mainly responsible for the hepatic uptake of various cationic drugs, closely associated with drug-induced liver injury (DILI). Screening and identifying potent OCT1 inhibitors with little toxicity in natural products is of great value in alleviating OCT1-mediated liver injury. Flavonoids, a group of polyphenols commonly found in foodstuffs and herbal products, have been reported to cause transporter-mediated food/herb-drug interactions (FDIs). Our objective was to investigate potential inhibitors of OCT1 from 96 flavonoids, evaluate the hepatoprotective effects on retrorsine-induced liver injury, and clarify the structure-activity relationships of flavonoids with OCT1. Thirteen flavonoids exhibited significant inhibition (>50%) on OCT1 in OCT1-HEK293 cells. Among them, the five strongest flavonoid inhibitors (IC50 < 10 μM), including α-naphthoflavone, apigenin, 6-hydroxyflavone, luteolin, and isosilybin markedly decreased oxaliplatin-induced cytotoxicity. In retrorsine-induced liver injury models, they also reduced alanine aminotransferase (ALT) and aspartate aminotransferase (AST) to different levels, the best of which was 6-hydroxyflavone. The pharmacophore model clarified that hydrogen bond acceptors at the 4,8,5' position might play a vital role in the inhibitory effect of flavonoids on OCT1. Taken together, our findings would pave the way to predicting the potential risks of flavonoid-related FDIs in humans and optimizing flavonoid structure to alleviate OCT1-mediated liver injury.

Functional genomics reveals an off-target dependency of drug synergy in gastric cancer therapy

BACKGROUND

Integrating molecular-targeted agents into combination chemotherapy is transformative for enhancing treatment outcomes in cancer. However, realizing the full potential of this approach requires a clear comprehension of the genetic dependencies underlying drug synergy. While the interactions between conventional chemotherapeutics are well-explored, the interplay of molecular-targeted agents with conventional chemotherapeutics remains a frontier in cancer treatment. Hence, we leveraged a powerful functional genomics approach to decode genomic dependencies that drive synergy in molecular-targeted agent/chemotherapeutic combinations in gastric adenocarcinoma, addressing a critical need in gastric cancer therapy.

METHODS

We screened pharmacological interactions between fifteen molecular-targeted agent/conventional chemotherapeutic pairs in gastric adenocarcinoma cells, and examined the genome-scale genetic dependencies of synergy integrating genome-wide CRISPR screening with the shRNA-based signature assay. We validated the synergy in cell death using fluorescence-based and lysis-dependent inference of cell death kinetics assay, and validated the genetic dependencies by single-gene knockout experiments.

RESULTS

Our combination screen identified SN-38/erlotinib as the drug pair with the strongest synergism. Functional genomics assays unveiled a genetic dependency signature of SN-38/erlotinib identical to SN-38. Remarkably, the enhanced cell death with improved kinetics induced by SN-38/erlotinib was attributed to erlotinib's off-target effect, inhibiting ABCG2, rather than its on-target effect on EGFR.

CONCLUSION

In the era of precision medicine, where emphasis on primary drug targets prevails, our research challenges this paradigm by showcasing a robust synergy underpinned by an off-target dependency. Further dissection of the intricate genetic dependencies that underlie synergy can pave the way to developing more effective combination strategies in gastric cancer therapy.

Functional genomic screens with death rate analyses reveal mechanisms of drug action

A common approach for understanding how drugs induce their therapeutic effects is to identify the genetic determinants of drug sensitivity. Because 'chemo-genetic profiles' are performed in a pooled format, inference of gene function is subject to several confounding influences related to variation in growth rates between clones. In this study, we developed Method for Evaluating Death Using a Simulation-assisted Approach (MEDUSA), which uses time-resolved measurements, along with model-driven constraints, to reveal the combination of growth and death rates that generated the observed drug response. MEDUSA is uniquely effective at identifying death regulatory genes. We apply MEDUSA to characterize DNA damage-induced lethality in the presence and absence of p53. Loss of p53 switches the mechanism of DNA damage-induced death from apoptosis to a non-apoptotic death that requires high respiration. These findings demonstrate the utility of MEDUSA both for determining the genetic dependencies of lethality and for revealing opportunities to potentiate chemo-efficacy in a cancer-specific manner.

Calcineurin promotes adaptation to chronic stress through two distinct mechanisms

Adaptation to environmental stress requires coordination between stress-defense programs and cell cycle progression. The immediate response to many stressors has been well characterized, but how cells survive in challenging environments long term is unknown. Here, we investigate the role of the stress-activated phosphatase calcineurin (CN) in adaptation to chronic CaCl2 stress in Saccharomyces cerevisiae. We find that prolonged exposure to CaCl2 impairs mitochondrial function and demonstrate that cells respond to this stressor using two CN-dependent mechanisms-one that requires the downstream transcription factor Crz1 and another that is Crz1 independent. Our data indicate that CN maintains cellular fitness by promoting cell cycle progression and preventing CaCl2-induced cell death. When Crz1 is present, transient CN activation suppresses cell death and promotes adaptation despite high levels of mitochondrial loss. However, in the absence of Crz1, prolonged activation of CN prevents mitochondrial loss and further cell death by upregulating glutathione biosynthesis genes thereby mitigating damage from reactive oxygen species. These findings illustrate how cells maintain long-term fitness during chronic stress and suggest that CN promotes adaptation in challenging environments by multiple mechanisms.

Calcineurin promotes adaptation to chronic stress through two distinct mechanisms

Adaptation to environmental stress requires coordination between stress-defense programs and cell cycle progression. The immediate response to many stressors has been well characterized, but how cells survive in challenging environments long-term is unknown. Here, we investigate the role of the stress-activated phosphatase calcineurin (CN) in adaptation to chronic CaCl2 stress in Saccharomyces cerevisiae. We find that prolonged exposure to CaCl2 impairs mitochondrial function and demonstrate that cells respond to this stressor using two CN-dependent mechanisms - one that requires the downstream transcription factor Crz1 and another that is Crz1-independent. Our data indicate that CN maintains cellular fitness by promoting cell cycle progression and preventing CaCl2-induced cell death. When Crz1 is present, transient CN activation suppresses cell death and promotes adaptation despite high levels of mitochondrial loss. However, in the absence of Crz1, prolonged activation of CN prevents mitochondrial loss and further cell death by upregulating glutathione (GSH) biosynthesis genes thereby mitigating damage from reactive oxygen species. These findings illustrate how cells maintain long-term fitness during chronic stress and suggest that CN promotes adaptation in challenging environments by multiple mechanisms.

Measuring Drug Response with Single-Cell Growth Rate Quantification

Intratumoral heterogeneity is a substantial cause of drug resistance development during chemotherapy or other drug treatments for cancer. Therefore, monitoring and measuring cell exposure and response to drugs at the single-cell level are crucial. Previous research suggested that the single-cell growth rate can be used to investigate drug-cell interactions. However, currently established methods for quantifying single-cell growth are limited to isolated or monolayer cells. Here, we introduce a technique that accurately measures both 2D and 3D cell growth rates using label-free ratiometric stimulated Raman scattering (SRS) microscopy. We use deuterated amino acids, leucine, isoleucine, and valine, as tracers and measure the C-D SRS signal from deuterium-labeled proteins and the C-H SRS signal from unlabeled proteins simultaneously to determine the cell growth rate at the single-cell level. The technique offers single-cell level drug sensitivity measurement with a shorter turnaround time (within 12 h) than most traditional assays. The submicrometer resolution of the imaging technique allows us to examine the effects of chemotherapeutic drugs, including kinase inhibitors, mitotic inhibitors, and topoisomerase II inhibitors, on both the cell growth rate and morphology. The capability of quantifying 3D cell growth rates provides insight into a deeper understanding of the cell-drug interaction in the actual tumor environment.

Functional genomics reveals an off-target dependency of drug synergy in gastric cancer therapy

The rational combination of anticancer agents is critical to improving patient outcomes in cancer. Nonetheless, most combination regimens in the clinic result from empirical methodologies disregarding insight into the mechanism of action and missing the opportunity to improve therapy outcomes incrementally. Deciphering the genetic dependencies and vulnerabilities responsible for synergistic interactions is crucial for rationally developing effective anticancer drug combinations. Hence, we screened pairwise pharmacological interactions between molecular-targeted agents and conventional chemotherapeutics and examined the genome-scale genetic dependencies in gastric adenocarcinoma cell models. Since this type of cancer is mainly chemoresistant and incurable, clinical situations demand effective combination strategies. Our pairwise combination screen revealed SN38/erlotinib as the drug pair with the most robust synergism. Genome-wide CRISPR screening and a shRNA-based signature assay indicated that the genetic dependency/vulnerability signature of SN38/erlotinib is the same as SN38 alone. Additional investigation revealed that the enhanced cell death with improved death kinetics caused by the SN38/erlotinib combination is surprisingly due to erlotinib's off-target effect that inhibits ABCG2 but not its on-target effect on EGFR. Our results confirm that a genetic dependency signature different from the single-drug application may not be necessary for the synergistic interaction of molecular-targeted agents with conventional chemotherapeutics in gastric adenocarcinoma. The findings also demonstrated the efficacy of functional genomics approaches in unveiling biologically validated mechanisms of pharmacological interactions.

p53 controls choice between apoptotic and non-apoptotic death following DNA damage

DNA damage can activate apoptotic and non-apoptotic forms of cell death; however, it remains unclear what features dictate which type of cell death is activated. We report that p53 controls the choice between apoptotic and non-apoptotic death following exposure to DNA damage. In contrast to the conventional model, which suggests that p53-deficient cells should be resistant to DNA damage-induced cell death, we find that p53-deficient cells die at high rates following DNA damage, but exclusively using non-apoptotic mechanisms. Our experimental data and computational modeling reveal that non-apoptotic death in p53-deficient cells has not been observed due to use of assays that are either insensitive to cell death, or that specifically score apoptotic cells. Using functional genetic screening - with an analysis that enables computational inference of the drug-induced death rate - we find in p53-deficient cells that DNA damage activates a mitochondrial respiration-dependent form of cell death, called MPT-driven necrosis. Cells deficient for p53 have high basal respiration, which primes MPT-driven necrosis. Finally, using metabolite profiling, we identified mitochondrial activity-dependent metabolic vulnerabilities that can be targeted to potentiate the lethality of DNA damage specifically in p53-deficient cells. Our findings reveal how the dual functions of p53 in regulating mitochondrial activity and the DNA damage response combine to facilitate the choice between apoptotic and non-apoptotic death.

Inferring density-dependent population dynamics mechanisms through rate disambiguation for logistic birth-death processes

Density dependence is important in the ecology and evolution of microbial and cancer cells. Typically, we can only measure net growth rates, but the underlying density-dependent mechanisms that give rise to the observed dynamics can manifest in birth processes, death processes, or both. Therefore, we utilize the mean and variance of cell number fluctuations to separately identify birth and death rates from time series that follow stochastic birth-death processes with logistic growth. Our nonparametric method provides a novel perspective on stochastic parameter identifiability, which we validate by analyzing the accuracy in terms of the discretization bin size. We apply our method to the scenario where a homogeneous cell population goes through three stages: (1) grows naturally to its carrying capacity, (2) is treated with a drug that reduces its carrying capacity, and (3) overcomes the drug effect to restore its original carrying capacity. In each stage, we disambiguate whether the dynamics occur through the birth process, death process, or some combination of the two, which contributes to understanding drug resistance mechanisms. In the case of limited sample sizes, we provide an alternative method based on maximum likelihood and solve a constrained nonlinear optimization problem to identify the most likely density dependence parameter for a given cell number time series. Our methods can be applied to other biological systems at different scales to disambiguate density-dependent mechanisms underlying the same net growth rate.

Challenges and Opportunities Modeling the Dynamic Tumor Matrisome

We need novel strategies to target the complexity of cancer and, particularly, of metastatic disease. As an example of this complexity, certain tissues are particularly hospitable environments for metastases, whereas others do not contain fertile microenvironments to support cancer cell growth. Continuing evidence that the extracellular matrix (ECM) of tissues is one of a host of factors necessary to support cancer cell growth at both primary and secondary tissue sites is emerging. Research on cancer metastasis has largely been focused on the molecular adaptations of tumor cells in various cytokine and growth factor environments on 2-dimensional tissue culture polystyrene plates. Intravital imaging, conversely, has transformed our ability to watch, in real time, tumor cell invasion, intravasation, extravasation, and growth. Because the interstitial ECM that supports all cells in the tumor microenvironment changes over time scales outside the possible window of typical intravital imaging, bioengineers are continuously developing both simple and sophisticated in vitro controlled environments to study tumor (and other) cell interactions with this matrix. In this perspective, we focus on the cellular unit responsible for upholding the pathologic homeostasis of tumor-bearing organs, cancer-associated fibroblasts (CAFs), and their self-generated ECM. The latter, together with tumoral and other cell secreted factors, constitute the "tumor matrisome". We share the challenges and opportunities for modeling this dynamic CAF/ECM unit, the tools and techniques available, and how the tumor matrisome is remodeled (e.g., via ECM proteases). We posit that increasing information on tumor matrisome dynamics may lead the field to alternative strategies for personalized medicine outside genomics.

Bacterial diet modulates tamoxifen-induced death via host fatty acid metabolism

Tamoxifen is a selective estrogen receptor (ER) modulator that is used to treat ER-positive breast cancer, but that at high doses kills both ER-positive and ER-negative breast cancer cells. We recapitulate this off-target effect in Caenorhabditis elegans, which does not have an ER ortholog. We find that different bacteria dramatically modulate tamoxifen toxicity in C. elegans, with a three-order of magnitude difference between animals fed Escherichia coli, Comamonas aquatica, and Bacillus subtilis. Remarkably, host fatty acid (FA) biosynthesis mitigates tamoxifen toxicity, and different bacteria provide the animal with different FAs, resulting in distinct FA profiles. Surprisingly these bacteria modulate tamoxifen toxicity by different death mechanisms, some of which are modulated by FA supplementation and others by antioxidants. Together, this work reveals a complex interplay between microbiota, FA metabolism and tamoxifen toxicity that may provide a blueprint for similar studies in more complex mammals.

Here, Diot et al. use the nematode Caenorhabditis elegans as a model to identify off-target toxicity mechanisms for tamoxifen, and find that these include fatty acid metabolism and cell death, which can be modulated by different bacterial species.

Evolution of 3-(4-hydroxyphenyl)indoline-2-one as a scaffold for potent and selective anticancer activity

Development of targeted anticancer modalities has prompted a new era in cancer treatment that is notably different from the age of radical surgery and highly toxic chemotherapy. Behind each effective compound is a rich and complex history from first identification of chemical matter, detailed optimization, and mechanistic investigations, ultimately leading to exciting molecules for drug development. Herein we review the history and on-going journey of one such anticancer scaffold, the 3-(4-hydroxyphenyl)indoline-2-ones. With humble beginnings in 19th century Bavaria, we review this scaffold's synthetic history and anticancer optimization, including its recent demonstration of tumor eradication of drug-resistant, estrogen receptor-positive breast cancer. Compounds containing the 3-(4-hydroxyphenyl)indoline-2-one pharmacophore are emerging as intriguing candidates for the treatment of cancer.

Materials-driven approaches to understand extrinsic drug resistance in cancer

Metastatic cancer has a poor prognosis, because it is broadly disseminated and associated with both intrinsic and acquired drug resistance. Critical unmet needs in effectively killing drug resistant cancer cells include overcoming the drug desensitization characteristics of some metastatic cancers/lesions, and tailoring therapeutic regimens to both the tumor microenvironment and the genetic profiles of the resident cancer cells. Bioengineers and materials scientists are developing technologies to determine how metastatic sites exclude therapies, and how extracellular factors (including cells, proteins, metabolites, extracellular matrix, and abiotic factors) at metastatic sites significantly affect drug pharmacodynamics. Two looming challenges are determining which feature, or combination of features, from the tumor microenvironment drive drug resistance, and what the relative impact is of extracellular signals vs. intrinsic cell genetics in determining drug response. Sophisticated systems biology tools that can de-convolve a crowded network of signals and responses, as well as controllable microenvironments capable of providing discrete and tunable extracellular cues can help us begin to interrogate the high dimensional interactions governing drug resistance in patients.

FLICK: an optimized plate reader-based assay to infer cell death kinetics

Evaluating drug sensitivity is improved by directly quantifying death kinetics, rather than correlates of viability, such as metabolic activity. This is challenging, requiring time-lapse microscopy and genetically encoded labels to distinguish live and dead cells. Here, we describe fluorescence-based and lysis-dependent inference of cell death kinetics (FLICK). This method requires only a standard fluorescence plate reader, retaining the high-throughput nature and broad accessibility of common viability assays. However, FLICK specifically quantifies death, including an accurate inference of death kinetics.

For complete details on the use and execution of this protocol, please refer to Richards et al. (2020).

•

FLICK generates a death-specific measure of drug response not a correlate of viability

•

FLICK is high-throughput and easy to use like common plate reader-based assays

•

FLICK can be used to generate any commonly used drug response metric

•

FLICK can be used to infer cell death kinetics

Human ALS/FTD brain organoid slice cultures display distinct early astrocyte and targetable neuronal pathology

Amyotrophic lateral sclerosis overlapping with frontotemporal dementia (ALS/FTD) is a fatal and currently untreatable disease characterized by rapid cognitive decline and paralysis. Elucidating initial cellular pathologies is central to therapeutic target development, but obtaining samples from presymptomatic patients is not feasible. Here, we report the development of a cerebral organoid slice model derived from human induced pluripotent stem cells (iPSCs) that recapitulates mature cortical architecture and displays early molecular pathology of C9ORF72 ALS/FTD. Using a combination of single-cell RNA sequencing and biological assays, we reveal distinct transcriptional, proteostasis and DNA repair disturbances in astroglia and neurons. We show that astroglia display increased levels of the autophagy signaling protein P62 and that deep layer neurons accumulate dipeptide repeat protein poly(GA), DNA damage and undergo nuclear pyknosis that could be pharmacologically rescued by GSK2606414. Thus, patient-specific iPSC-derived cortical organoid slice cultures are a reproducible translational platform to investigate preclinical ALS/FTD mechanisms as well as novel therapeutic approaches.

ELP-dependent expression of MCL1 promotes resistance to EGFR inhibition in triple-negative breast cancer cells

Targeted therapeutics for cancer generally exploit "oncogene addiction," a phenomenon in which the growth and survival of tumor cells depend on the activity of a particular protein. However, the efficacy of oncogene-targeted therapies varies substantially. For instance, targeting epidermal growth factor receptor (EGFR) signaling is effective in some non-small cell lung cancer (NSCLC) but not in triple-negative breast cancer (TNBC), although these cancers show a similar degree of increase in EGFR activity. Using a genome-wide CRISPR-Cas9 genetic knockout screen, we found that the Elongator (ELP) complex mediates insensitivity to the EGFR inhibitor erlotinib in TNBC cells by promoting the synthesis of the antiapoptotic protein Mcl-1. Depleting ELP proteins promoted apoptotic cell death in an EGFR inhibition-dependent manner. Pharmacological inhibition of Mcl-1 synergized with EGFR inhibition in a panel of genetically diverse TNBC cells. The findings indicate that TNBC "addiction" to EGFR signaling is masked by the ELP complex and that resistance to EGFR inhibitors in TNBC might be overcome by cotargeting Mcl-1.