Network Analyses of Brain Tumor Patients' Multiomic Data Reveals Pharmacological Opportunities to Alter Cell State Transitions.

Glioblastoma Multiforme (GBM) remains a particularly difficult cancer to treat, and survival outcomes remain poor. In addition to the lack of dedicated drug discovery programs for GBM, extensive intratumor heterogeneity and epigenetic plasticity related to cell-state transitions are major roadblocks to successful drug therapy in GBM. To study these phenomenon, publicly available snRNAseq and bulk RNAseq data from patient samples were used to categorize cells from patients into four cell states (i.e. phenotypes), namely: (i) neural progenitor-like (NPC-like), (ii) oligodendrocyte progenitor-like (OPC-like), (iii) astrocyte- like (AC-like), and (iv) mesenchymal-like (MES-like). Patients were subsequently grouped into subpopulations based on which cell-state was the most dominant in their respective tumor. By incorporating phosphoproteomic measurements from the same patients, a protein-protein interaction network (PPIN) was constructed for each cell state. These four-cell state PPINs were pooled to form a single Boolean network that was used for in silico protein knockout simulations to investigate mechanisms that either promote or prevent cell state transitions. Simulation results were input into a boosted tree machine learning model which predicted the cell states or phenotypes of GBM patients from an independent public data source, the Glioma Longitudinal Analysis (GLASS) Consortium. Combining the simulation results and the machine learning predictions, we generated hypotheses for clinically relevant causal mechanisms of cell state transitions. For example, the transcription factor TFAP2A can be seen to promote a transition from the NPC-like to the MES-like state. Such protein nodes and the associated signaling pathways provide potential drug targets that can be further tested in vitro and support cell state-directed (CSD) therapy.

Low-frequency ERK and Akt activity dynamics are predictive of stochastic cell division events.

Understanding the dynamics of intracellular signaling pathways, such as ERK1/2 (ERK) and Akt1/2 (Akt), in the context of cell fate decisions is important for advancing our knowledge of cellular processes and diseases, particularly cancer. While previous studies have established associations between ERK and Akt activities and proliferative cell fate, the heterogeneity of single-cell responses adds complexity to this understanding. This study employed a data-driven approach to address this challenge, developing machine learning models trained on a dataset of growth factor-induced ERK and Akt activity time courses in single cells, to predict cell division events. The most predictive models were developed by applying discrete wavelet transforms (DWTs) to extract low-frequency features from the time courses, followed by using Ensemble Integration, a data integration and predictive modeling framework. The results demonstrated that these models effectively predicted cell division events in MCF10A cells (F-measure=0.524, AUC=0.726). ERK dynamics were found to be more predictive than Akt, but the combination of both measurements further enhanced predictive performance. The ERK model`s performance also generalized to predicting division events in RPE cells, indicating the potential applicability of these models and our data-driven methodology for predicting cell division across different biological contexts. Interpretation of these models suggested that ERK dynamics throughout the cell cycle, rather than immediately after growth factor stimulation, were associated with the likelihood of cell division. Overall, this work contributes insights into the predictive power of intra-cellular signaling dynamics for cell fate decisions, and highlights the potential of machine learning approaches in unraveling complex cellular behaviors.

Network Analyses of Brain Tumor Patients' Multiomic Data Reveals Pharmacological Opportunities to Alter Cell State Transitions.

Glioblastoma Multiforme (GBM) remains a particularly difficult cancer to treat, and survival outcomes remain poor. In addition to the lack of dedicated drug discovery programs for GBM, extensive intratumor heterogeneity and epigenetic plasticity related to cell-state transitions are major roadblocks to successful drug therapy in GBM. To study these phenomenon, publicly available snRNAseq and bulk RNAseq data from patient samples were used to categorize cells from patients into four cell states (i.e. phenotypes), namely: (i) neural progenitor-like (NPC-like), (ii) oligodendrocyte progenitor-like (OPC-like), (iii) astrocyte-like (AC-like), and (iv) mesenchymal-like (MES-like). Patients were subsequently grouped into subpopulations based on which cell-state was the most dominant in their respective tumor. By incorporating phosphoproteomic measurements from the same patients, a protein-protein interaction network (PPIN) was constructed for each cell state. These four-cell state PPINs were pooled to form a single Boolean network that was used for in silico protein knockout simulations to investigate mechanisms that either promote or prevent cell state transitions. Simulation results were input into a boosted tree machine learning model which predicted the cell states or phenotypes of GBM patients from an independent public data source, the Glioma Longitudinal Analysis (GLASS) Consortium. Combining the simulation results and the machine learning predictions, we generated hypotheses for clinically relevant causal mechanisms of cell state transitions. For example, the transcription factor TFAP2A can be seen to promote a transition from the NPC-like to the MES-like state. Such protein nodes and the associated signaling pathways provide potential drug targets that can be further tested in vitro and support cell state-directed (CSD) therapy.

MOBILE pipeline enables identification of context-specific networks and regulatory mechanisms.

Robust identification of context-specific network features that control cellular phenotypes remains a challenge. We here introduce MOBILE (Multi-Omics Binary Integration via Lasso Ensembles) to nominate molecular features associated with cellular phenotypes and pathways. First, we use MOBILE to nominate mechanisms of interferon-γ (IFNγ) regulated PD-L1 expression. Our analyses suggest that IFNγ-controlled PD-L1 expression involves BST2, CLIC2, FAM83D, ACSL5, and HIST2H2AA3 genes, which were supported by prior literature. We also compare networks activated by related family members transforming growth factor-beta 1 (TGFß1) and bone morphogenetic protein 2 (BMP2) and find that differences in ligand-induced changes in cell size and clustering properties are related to differences in laminin/collagen pathway activity. Finally, we demonstrate the broad applicability and adaptability of MOBILE by analyzing publicly available molecular datasets to investigate breast cancer subtype specific networks. Given the ever-growing availability of multi-omics datasets, we envision that MOBILE will be broadly useful for identification of context-specific molecular features and pathways.

Predicting anti-cancer drug combination responses with a temporal cell state network model.

Cancer chemotherapy combines multiple drugs, but predicting the effects of drug combinations on cancer cell proliferation remains challenging, even for simple in vitro systems. We hypothesized that by combining knowledge of single drug dose responses and cell state transition network dynamics, we could predict how a population of cancer cells will respond to drug combinations. We tested this hypothesis here using three targeted inhibitors of different cell cycle states in two different cell lines in vitro. We formulated a Markov model to capture temporal cell state transitions between different cell cycle phases, with single drug data constraining how drug doses affect transition rates. This model was able to predict the landscape of all three different pairwise drug combinations across all dose ranges for both cell lines with no additional data. While further application to different cell lines, more drugs, additional cell state networks, and more complex co-culture or in vivo systems remain, this work demonstrates how currently available or attainable information could be sufficient for prediction of drug combination response for single cell lines in vitro.

Predicting individual-specific cardiotoxicity responses induced by tyrosine kinase inhibitors.

Introduction: Tyrosine kinase inhibitor drugs (TKIs) are highly effective cancer drugs, yet many TKIs are associated with various forms of cardiotoxicity. The mechanisms underlying these drug-induced adverse events remain poorly understood. We studied mechanisms of TKI-induced cardiotoxicity by integrating several complementary approaches, including comprehensive transcriptomics, mechanistic mathematical modeling, and physiological assays in cultured human cardiac myocytes. Methods: Induced pluripotent stem cells (iPSCs) from two healthy donors were differentiated into cardiac myocytes (iPSC-CMs), and cells were treated with a panel of 26 FDA-approved TKIs. Drug-induced changes in gene expression were quantified using mRNA-seq, changes in gene expression were integrated into a mechanistic mathematical model of electrophysiology and contraction, and simulation results were used to predict physiological outcomes. Results: Experimental recordings of action potentials, intracellular calcium, and contraction in iPSC-CMs demonstrated that modeling predictions were accurate, with 81% of modeling predictions across the two cell lines confirmed experimentally. Surprisingly, simulations of how TKI-treated iPSC-CMs would respond to an additional arrhythmogenic insult, namely, hypokalemia, predicted dramatic differences between cell lines in how drugs affected arrhythmia susceptibility, and these predictions were confirmed experimentally. Computational analysis revealed that differences between cell lines in the upregulation or downregulation of particular ion channels could explain how TKI-treated cells responded differently to hypokalemia. Discussion: Overall, the study identifies transcriptional mechanisms underlying cardiotoxicity caused by TKIs, and illustrates a novel approach for integrating transcriptomics with mechanistic mathematical models to generate experimentally testable, individual-specific predictions of adverse event risk.

Relating individual cell division events to single-cell ERK and Akt activity time courses.

Biochemical correlates of stochastic single-cell fates have been elusive, even for the well-studied mammalian cell cycle. We monitored single-cell dynamics of the ERK and Akt pathways, critical cell cycle progression hubs and anti-cancer drug targets, and paired them to division events in the same single cells using the non-transformed MCF10A epithelial line. Following growth factor treatment, in cells that divide both ERK and Akt activities are significantly higher within the S-G2 time window (~ 8.5-40 h). Such differences were much smaller in the pre-S-phase, restriction point window which is traditionally associated with ERK and Akt activity dependence, suggesting unappreciated roles for ERK and Akt in S through G2. Simple metrics of central tendency in this time window are associated with subsequent cell division fates. ERK activity was more strongly associated with division fates than Akt activity, suggesting Akt activity dynamics may contribute less to the decision driving cell division in this context. We also find that ERK and Akt activities are less correlated with each other in cells that divide. Network reconstruction experiments demonstrated that this correlation behavior was likely not due to crosstalk, as ERK and Akt do not interact in this context, in contrast to other transformed cell types. Overall, our findings support roles for ERK and Akt activity throughout the cell cycle as opposed to just before the restriction point, and suggest ERK activity dynamics may be more important than Akt activity dynamics for driving cell division in this non-transformed context.

A scalable, open-source implementation of a large-scale mechanistic model for single cell proliferation and death signaling.

Mechanistic models of how single cells respond to different perturbations can help integrate disparate big data sets or predict response to varied drug combinations. However, the construction and simulation of such models have proved challenging. Here, we developed a python-based model creation and simulation pipeline that converts a few structured text files into an SBML standard and is high-performance- and cloud-computing ready. We applied this pipeline to our large-scale, mechanistic pan-cancer signaling model (named SPARCED) and demonstrate it by adding an IFNγ pathway submodel. We then investigated whether a putative crosstalk mechanism could be consistent with experimental observations from the LINCS MCF10A Data Cube that IFNγ acts as an anti-proliferative factor. The analyses suggested this observation can be explained by IFNγ-induced SOCS1 sequestering activated EGF receptors. This work forms a foundational recipe for increased mechanistic model-based data integration on a single-cell level, an important building block for clinically-predictive mechanistic models.

Leveraging modeling and simulation to optimize the therapeutic window for epigenetic modifier drugs.

Dysregulated epigenetic processes can lead to altered gene expression and give rise to malignant transformation and tumorigenesis. Epigenetic drugs aim to revert the phenotype of cancer cells to normally functioning cells, and are developed and applied to treat both hematological and solid cancers. Despite this promising therapeutic avenue, the successful development of epigenetic modulators has been challenging. We argue that besides identifying the right responder patient population, the selection of an optimized dosing regimen is equally important. For the majority of epigenetic modulators, hematological adverse effects such as thrombocytopenia, anemia or neutropenia are frequently observed and may limit their therapeutic potential. Therefore, one of the key challenges is to identify a dosing regimen that maximizes drug efficacy and minimizes toxicity. This requires a good understanding of the quantitative relationship between the administered dose, the drug exposure and the magnitude and duration of drug response related to safety and efficacy. With case examples, we highlight how modeling and simulation has been successfully applied to address those questions. As an outlook, we suggest the combination of efficacy and safety prediction models that capture the quantitative, mechanistic relationships governing the balance between their safety and efficacy dynamics. A stepwise approach for its implementation is presented. Utilizing in silico explorations, the impact of dosing regimen on the therapeutic window can be explored. This will serve as a basis to select the most promising dosing regimen that maximizes efficacy while minimizing adverse effects and to increase the probability of success for the given epigenetic drug.

Proteomic cellular signatures of kinase inhibitor-induced cardiotoxicity.

Drug Toxicity Signature Generation Center (DToxS) at the Icahn School of Medicine at Mount Sinai is one of the centers for the NIH Library of Integrated Network-Based Cellular Signatures (LINCS) program. Its key aim is to generate proteomic and transcriptomic signatures that can predict cardiotoxic adverse effects of kinase inhibitors approved by the Food and Drug Administration. Towards this goal, high throughput shotgun proteomics experiments (308 cell line/drug combinations +64 control lysates) have been conducted. Using computational network analyses, these proteomic data can be integrated with transcriptomic signatures, generated in tandem, to identify cellular signatures of cardiotoxicity that may predict kinase inhibitor-induced toxicity and enable possible mitigation. Both raw and processed proteomics data have passed several quality control steps and been made publicly available on the PRIDE database. This broad protein kinase inhibitor-stimulated human cardiomyocyte proteomic data and signature set is valuable for prediction of drug toxicities.

Anti-invasive efficacy and survival benefit of the YAP-TEAD inhibitor verteporfin in preclinical glioblastoma models.

BACKGROUND: Glioblastoma (GBM) remains a largely incurable disease as current therapy fails to target the invasive nature of glioma growth in disease progression and recurrence. Here, we use the FDA-approved drug and small molecule Hippo inhibitor Verteporfin (VP) to target YAP-TEAD activity, known to mediate convergent aspects of tumor invasion/metastasis, and assess the drug's efficacy and survival benefit in GBM models. METHODS: Up to 8 low-passage patient-derived GBM cell lines with distinct genomic drivers, including 3 primary/recurrent pairs, were treated with VP or vehicle (VEH) to assess in vitro effects on proliferation, migration, invasion, YAP-TEAD activity, and transcriptomics. Patient-derived orthotopic xenograft (PDX) models were used to assess VP's brain penetrance and effects on tumor burden and survival. RESULTS: VP treatment disturbed YAP/TAZ-TEAD activity; disrupted transcriptome signatures related to invasion, epithelial-to-mesenchymal, and proneural-to-mesenchymal transition, phenocopying TEAD1-knockout effects; and impaired tumor migration/invasion dynamics across primary and recurrent GBM lines. In an aggressive orthotopic PDX GBM model, short-term VP treatment consistently diminished core and infiltrative tumor burden, which was associated with decreased tumor expression of Ki67, nuclear YAP, TEAD1, and TEAD-associated targets EGFR, CDH2, and ITGB1. Finally, long-term VP treatment appeared nontoxic and conferred survival benefit compared to VEH in 2 PDX models: as monotherapy in primary (de novo) GBM and in combination with Temozolomide chemoradiation in recurrent GBM, where VP treatment associated with increased MGMT methylation. CONCLUSIONS: We demonstrate combined anti-invasive and anti-proliferative efficacy for VP with survival benefit in preclinical GBM models, indicating potential therapeutic value of this already FDA-approved drug if repurposed for GBM patients.

Transcriptomic profiling of human cardiac cells predicts protein kinase inhibitor-associated cardiotoxicity.

Kinase inhibitors (KIs) represent an important class of anti-cancer drugs. Although cardiotoxicity is a serious adverse event associated with several KIs, the reasons remain poorly understood, and its prediction remains challenging. We obtain transcriptional profiles of human heart-derived primary cardiomyocyte like cell lines treated with a panel of 26 FDA-approved KIs and classify their effects on subcellular pathways and processes. Individual cardiotoxicity patient reports for these KIs, obtained from the FDA Adverse Event Reporting System, are used to compute relative risk scores. These are then combined with the cell line-derived transcriptomic datasets through elastic net regression analysis to identify a gene signature that can predict risk of cardiotoxicity. We also identify relationships between cardiotoxicity risk and structural/binding profiles of individual KIs. We conclude that acute transcriptomic changes in cell-based assays combined with drug substructures are predictive of KI-induced cardiotoxicity risk, and that they can be informative for future drug discovery.

Predicting In Vivo Efficacy from In Vitro Data: Quantitative Systems Pharmacology Modeling for an Epigenetic Modifier Drug in Cancer.

Reliably predicting in vivo efficacy from in vitro data would facilitate drug development by reducing animal usage and guiding drug dosing in human clinical trials. However, such prediction remains challenging. Here, we built a quantitative pharmacokinetic/pharmacodynamic (PK/PD) mathematical model capable of predicting in vivo efficacy in animal xenograft models of tumor growth while trained almost exclusively on in vitro cell culture data sets. We studied a chemical inhibitor of LSD1 (ORY-1001), a lysine-specific histone demethylase enzyme with epigenetic function, and drug-induced regulation of target engagement, biomarker levels, and tumor cell growth across multiple doses administered in a pulsed and continuous fashion. A PK model of unbound plasma drug concentration was linked to the in vitro PD model, which enabled the prediction of in vivo tumor growth dynamics across a range of drug doses and regimens. Remarkably, only a change in a single parameter-the one controlling intrinsic cell/tumor growth in the absence of drug-was needed to scale the PD model from the in vitro to in vivo setting. These findings create a framework for using in vitro data to predict in vivo drug efficacy with clear benefits to reducing animal usage while enabling the collection of dense time course and dose response data in a highly controlled in vitro environment.

Wilm's tumor 1 promotes memory flexibility.

Under physiological conditions, strength and persistence of memory must be regulated in order to produce behavioral flexibility. In fact, impairments in memory flexibility are associated with pathologies such as post-traumatic stress disorder or autism; however, the underlying mechanisms that enable memory flexibility are still poorly understood. Here, we identify transcriptional repressor Wilm's Tumor 1 (WT1) as a critical synaptic plasticity regulator that decreases memory strength, promoting memory flexibility. WT1 is activated in the hippocampus following induction of long-term potentiation (LTP) or learning. WT1 knockdown enhances CA1 neuronal excitability, LTP and long-term memory whereas its overexpression weakens memory retention. Moreover, forebrain WT1-deficient mice show deficits in both reversal, sequential learning tasks and contextual fear extinction, exhibiting impaired memory flexibility. We conclude that WT1 limits memory strength or promotes memory weakening, thus enabling memory flexibility, a process that is critical for learning from new experiences.

A Multi-center Study on the Reproducibility of Drug-Response Assays in Mammalian Cell Lines.

Evidence that some high-impact biomedical results cannot be repeated has stimulated interest in practices that generate findable, accessible, interoperable, and reusable (FAIR) data. Multiple papers have identified specific examples of irreproducibility, but practical ways to make data more reproducible have not been widely studied. Here, five research centers in the NIH LINCS Program Consortium investigate the reproducibility of a prototypical perturbational assay: quantifying the responsiveness of cultured cells to anti-cancer drugs. Such assays are important for drug development, studying cellular networks, and patient stratification. While many experimental and computational factors impact intra- and inter-center reproducibility, the factors most difficult to identify and control are those with a strong dependency on biological context. These factors often vary in magnitude with the drug being analyzed and with growth conditions. We provide ways to identify such context-sensitive factors, thereby improving both the theory and practice of reproducible cell-based assays.

Mitochondrial origins of fractional control in regulated cell death.

Individual cells in clonal populations often respond differently to environmental changes; for binary phenotypes, such as cell death, this can be measured as a fractional response. These types of responses have been attributed to cell-intrinsic stochastic processes and variable abundances of biochemical constituents, such as proteins, but the influence of organelles is still under investigation. We use the response to TNF-related apoptosis inducing ligand (TRAIL) and a new statistical framework for determining parameter influence on cell-to-cell variability through the inference of variance explained, DEPICTIVE, to demonstrate that variable mitochondria abundance correlates with cell survival and determines the fractional cell death response. By quantitative data analysis and modeling we attribute this effect to variable effective concentrations at the mitochondria surface of the pro-apoptotic proteins Bax/Bak. Further, our study suggests that inhibitors of anti-apoptotic Bcl-2 family proteins, used in cancer treatment, may increase the diversity of cellular responses, enhancing resistance to treatment.

A mechanistic pan-cancer pathway model informed by multi-omics data interprets stochastic cell fate responses to drugs and mitogens.

Most cancer cells harbor multiple drivers whose epistasis and interactions with expression context clouds drug and drug combination sensitivity prediction. We constructed a mechanistic computational model that is context-tailored by omics data to capture regulation of stochastic proliferation and death by pan-cancer driver pathways. Simulations and experiments explore how the coordinated dynamics of RAF/MEK/ERK and PI-3K/AKT kinase activities in response to synergistic mitogen or drug combinations control cell fate in a specific cellular context. In this MCF10A cell context, simulations suggest that synergistic ERK and AKT inhibitor-induced death is likely mediated by BIM rather than BAD, which is supported by prior experimental studies. AKT dynamics explain S-phase entry synergy between EGF and insulin, but simulations suggest that stochastic ERK, and not AKT, dynamics seem to drive cell-to-cell proliferation variability, which in simulations is predictable from pre-stimulus fluctuations in C-Raf/B-Raf levels. Simulations suggest MEK alteration negligibly influences transformation, consistent with clinical data. Tailoring the model to an alternate cell expression and mutation context, a glioma cell line, allows prediction of increased sensitivity of cell death to AKT inhibition. Our model mechanistically interprets context-specific landscapes between driver pathways and cell fates, providing a framework for designing more rational cancer combination therapy.

Analysis of copy number loss of the ErbB4 receptor tyrosine kinase in glioblastoma.

Current treatments for glioblastoma multiforme (GBM)-an aggressive form of brain cancer-are minimally effective and yield a median survival of 14.6 months and a two-year survival rate of 30%. Given the severity of GBM and the limitations of its treatment, there is a need for the discovery of novel drug targets for GBM and more personalized treatment approaches based on the characteristics of an individual's tumor. Most receptor tyrosine kinases-such as EGFR-act as oncogenes, but publicly available data from the Cancer Cell Line Encyclopedia (CCLE) indicates copy number loss in the ERBB4 RTK gene across dozens of GBM cell lines, suggesting a potential tumor suppressor role. This loss is mutually exclusive with loss of its cognate ligand NRG1 in CCLE as well, more strongly suggesting a functional role. The availability of higher resolution copy number data from clinical GBM patients in The Cancer Genome Atlas (TCGA) revealed that a region in Intron 1 of the ERBB4 gene was deleted in 69.1% of tumor samples harboring ERBB4 copy number loss; however, it was also found to be deleted in the matched normal tissue samples from these GBM patients (n = 81). Using the DECIPHER Genome Browser, we also discovered that this mutation occurs at approximately the same frequency in the general population as it does in the disease population. We conclude from these results that this loss in Intron 1 of the ERBB4 gene is neither a de novo driver mutation nor a predisposing factor to GBM, despite the indications from CCLE. A biological role of this significantly occurring genetic alteration is still unknown. While this is a negative result, the broader conclusion is that while copy number data from large cell line-based data repositories may yield compelling hypotheses, careful follow up with higher resolution copy number assays, patient data, and general population analyses are essential to codify initial hypotheses prior to investing experimental resources.

Integrating Transcriptomic Data with Mechanistic Systems Pharmacology Models for Virtual Drug Combination Trials.

Monotherapy clinical trials with mutation-targeted kinase inhibitors, despite some success in other cancers, have yet to impact glioblastoma (GBM). Besides insufficient blood-brain barrier penetration, combinations are key to overcoming obstacles such as intratumoral heterogeneity, adaptive resistance, and the epistatic nature of tumor genomics that cause mutation-targeted therapies to fail. With now hundreds of potential drugs, exploring the combination space clinically and preclinically is daunting. We are building a simulation-based approach that integrates patient-specific data with a mechanistic computational model of pan-cancer driver pathways (receptor tyrosine kinases, RAS/RAF/ERK, PI3K/AKT/mTOR, cell cycle, apoptosis, and DNA damage) to prioritize drug combinations by their simulated effects on tumor cell proliferation and death. Here we illustrate a first step, tailoring the model to 14 GBM patients from The Cancer Genome Atlas defined by an mRNA-seq transcriptome, and then simulating responses to three promiscuous FDA-approved kinase inhibitors (bosutinib, ibrutinib, and cabozantinib) with evidence for blood-brain barrier penetration. The model captures binding of the drug to primary targets and off-targets based on published affinity data and simulates responses of 100 heterogeneous tumor cells within a patient. Single drugs are marginally effective or even counterproductive. Common copy number alterations (PTEN loss, EGFR amplification, and NF1 loss) have a negligible correlation with single-drug or combination efficacy, reinforcing the importance of postgenetic approaches that account for kinase inhibitor promiscuity to match drugs to patients. Drug combinations tend to be either cytostatic or cytotoxic, but seldom both, highlighting the need for considering targeted and nontargeted therapy. Although we focus on GBM, the approach is generally applicable.

The Library of Integrated Network-Based Cellular Signatures NIH Program: System-Level Cataloging of Human Cells Response to Perturbations.

The Library of Integrated Network-Based Cellular Signatures (LINCS) is an NIH Common Fund program that catalogs how human cells globally respond to chemical, genetic, and disease perturbations. Resources generated by LINCS include experimental and computational methods, visualization tools, molecular and imaging data, and signatures. By assembling an integrated picture of the range of responses of human cells exposed to many perturbations, the LINCS program aims to better understand human disease and to advance the development of new therapies. Perturbations under study include drugs, genetic perturbations, tissue micro-environments, antibodies, and disease-causing mutations. Responses to perturbations are measured by transcript profiling, mass spectrometry, cell imaging, and biochemical methods, among other assays. The LINCS program focuses on cellular physiology shared among tissues and cell types relevant to an array of diseases, including cancer, heart disease, and neurodegenerative disorders. This Perspective describes LINCS technologies, datasets, tools, and approaches to data accessibility and reusability.