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.

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.

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.

Genomic signatures defining responsiveness to allopurinol and combination therapy for lung cancer identified by systems therapeutics analyses.

The ability to predict responsiveness to drugs in individual patients is limited. We hypothesized that integrating molecular information from databases would yield predictions that could be experimentally tested to develop transcriptomic signatures for specific drugs. We analyzed lung adenocarcinoma patient data from The Cancer Genome Atlas and identified a subset of patients in which xanthine dehydrogenase (XDH) expression correlated with decreased survival. We tested allopurinol, an FDA-approved drug that inhibits XDH, on human non-small-cell lung cancer (NSCLC) cell lines obtained from the Broad Institute Cancer Cell Line Encyclopedia and identified sensitive and resistant cell lines. We utilized the transcriptomic profiles of these cell lines to identify six-gene signatures for allopurinol-sensitive and allopurinol-resistant cell lines. Transcriptomic networks identified JAK2 as an additional target in allopurinol-resistant lines. Treatment of resistant cell lines with allopurinol and CEP-33779 (a JAK2 inhibitor) resulted in cell death. The effectiveness of allopurinol alone or allopurinol and CEP-33779 was verified in vivo using tumor formation in NCR-nude mice. We utilized the six-gene signatures to predict five additional allopurinol-sensitive NSCLC cell lines and four allopurinol-resistant cell lines susceptible to combination therapy. We searched the transcriptomic data from a library of patient-derived NSCLC tumors from the Jackson Laboratory to identify tumors that would be predicted to be sensitive to allopurinol or allopurinol + CEP-33779 treatment. Patient-derived tumors showed the predicted drug sensitivity in vivo. These data indicate that we can use integrated molecular information from cancer databases to predict drug responsiveness in individual patients and thus enable precision medicine.

Disruption of podocyte cytoskeletal biomechanics by dasatinib leads to nephrotoxicity.

Nephrotoxicity is a critical adverse event that leads to discontinuation of kinase inhibitor (KI) treatment. Here we show, through meta-analyses of FDA Adverse Event Reporting System, that dasatinib is associated with high risk for glomerular toxicity that is uncoupled from hypertension, suggesting a direct link between dasatinib and podocytes. We further investigate the cellular effects of dasatinib and other comparable KIs with varying risks of nephrotoxicity. Dasatinib treated podocytes show significant changes in focal adhesions, actin cytoskeleton, and morphology that are not observed with other KIs. We use phosphoproteomics and kinome profiling to identify the molecular mechanisms of dasatinib-induced injury to the actin cytoskeleton, and atomic force microscopy to quantify impairment to cellular biomechanics. Furthermore, chronic administration of dasatinib in mice causes reversible glomerular dysfunction, loss of stress fibers, and foot process effacement. We conclude that dasatinib induces nephrotoxicity through altered podocyte actin cytoskeleton, leading to injurious cellular biomechanics.

Interconnected network motifs control podocyte morphology and kidney function.

Podocytes are kidney cells with specialized morphology that is required for glomerular filtration. Diseases, such as diabetes, or drug exposure that causes disruption of the podocyte foot process morphology results in kidney pathophysiology. Proteomic analysis of glomeruli isolated from rats with puromycin-induced kidney disease and control rats indicated that protein kinase A (PKA), which is activated by adenosine 3',5'-monophosphate (cAMP), is a key regulator of podocyte morphology and function. In podocytes, cAMP signaling activates cAMP response element-binding protein (CREB) to enhance expression of the gene encoding a differentiation marker, synaptopodin, a protein that associates with actin and promotes its bundling. We constructed and experimentally verified a ß-adrenergic receptor-driven network with multiple feedback and feedforward motifs that controls CREB activity. To determine how the motifs interacted to regulate gene expression, we mapped multicompartment dynamical models, including information about protein subcellular localization, onto the network topology using Petri net formalisms. These computational analyses indicated that the juxtaposition of multiple feedback and feedforward motifs enabled the prolonged CREB activation necessary for synaptopodin expression and actin bundling. Drug-induced modulation of these motifs in diseased rats led to recovery of normal morphology and physiological function in vivo. Thus, analysis of regulatory motifs using network dynamics can provide insights into pathophysiology that enable predictions for drug intervention strategies to treat kidney disease.

A potential peptide therapeutic derived from the juxtamembrane domain of the epidermal growth factor receptor.

The epidermal growth factor receptor (EGFR) is involved in many cancers and EGFR has been heavily pursued as a drug target. Drugs targeting EGFR have shown promising clinical results for several cancer types. However, resistance to EGFR inhibitors often occurs, such as with KRAS mutant cancers, therefore new methods of targeting EGFR are needed. The juxtamembrane (JXM) domain of EGFR is critical for receptor activation and targeting this region could potentially be a new method of inhibiting EGFR. We hypothesized that the structural role of the JXM region could be mimicked by peptides encoding a JXM amino acid sequence, which could interfere with EGFR signaling and consequently could have anti-cancer activity. A peptide encoding EGFR 645-662 conjugated to the Tat sequence (TE-64562) displayed anti-cancer activity in multiple human cancer cell types with diminished activity in non-EGFR expressing cells and non-cancerous cells. In nude mice, TE-64562 delayed MDA-MB-231 tumor growth and prolonged survival, without inducing toxicity. TE-64562 induced non-apoptotic cell death after several hours and caspase-3-mediated apoptotic cell death with longer treatment. Mechanistically, TE-64562 bound to EGFR, inhibited its dimerization and caused its down-regulation. TE-64562 reduced phosphorylated and total EGFR levels but did not inhibit kinase activity and instead prolonged it. Our analysis of patient data from The Cancer Genome Atlas supported the hypothesis that down-regulation of EGFR is a potential therapeutic strategy, since phospho- and total-EGFR levels were strongly correlated in a large majority of patient tumor samples, indicating that lower EGFR levels are associated with lower phospho-EGFR levels and presumably less proliferative signals in breast cancer. Akt and Erk were inhibited by TE-64562 and this inhibition was observed in vivo in tumor tissue upon treatment with TE-64562. These results are the first to indicate that the JXM domain of EGFR is a viable drug target for several cancer types.