Identifying metabolic adaptations characteristic of cardiotoxicity using paired transcriptomics and metabolomics data integrated with a computational model of heart metabolism.

Improvements in the diagnosis and treatment of cancer have revealed long-term side effects of chemotherapeutics, particularly cardiotoxicity. Here, we present paired transcriptomics and metabolomics data characterizing in vitro cardiotoxicity to three compounds: 5-fluorouracil, acetaminophen, and doxorubicin. Standard gene enrichment and metabolomics approaches identify some commonly affected pathways and metabolites but are not able to readily identify metabolic adaptations in response to cardiotoxicity. The paired data was integrated with a genome-scale metabolic network reconstruction of the heart to identify shifted metabolic functions, unique metabolic reactions, and changes in flux in metabolic reactions in response to these compounds. Using this approach, we confirm previously seen changes in the p53 pathway by doxorubicin and RNA synthesis by 5-fluorouracil, we find evidence for an increase in phospholipid metabolism in response to acetaminophen, and we see a shift in central carbon metabolism suggesting an increase in metabolic demand after treatment with doxorubicin and 5-fluorouracil.

Computational model of cardiomyocyte apoptosis identifies mechanisms of tyrosine kinase inhibitor-induced cardiotoxicity.

Despite clinical observations of cardiotoxicity among cancer patients treated with tyrosine kinase inhibitors (TKIs), the molecular mechanisms by which these drugs affect the heart remain largely unknown. Mechanistic understanding of TKI-induced cardiotoxicity has been limited in part due to the complexity of tyrosine kinase signaling pathways and the multi-targeted nature of many of these drugs. TKI treatment has been associated with reactive oxygen species generation, mitochondrial dysfunction, and apoptosis in cardiomyocytes. To gain insight into the mechanisms mediating TKI-induced cardiotoxicity, this study constructs and validates a computational model of cardiomyocyte apoptosis, integrating intrinsic apoptotic and tyrosine kinase signaling pathways. The model predicts high levels of apoptosis in response to sorafenib, sunitinib, ponatinib, trastuzumab, and gefitinib, and lower levels of apoptosis in response to nilotinib and erlotinib, with the highest level of apoptosis induced by sorafenib. Knockdown simulations identified AP1, ASK1, JNK, MEK47, p53, and ROS as positive functional regulators of sorafenib-induced apoptosis of cardiomyocytes. Overexpression simulations identified Akt, IGF1, PDK1, and PI3K among the negative functional regulators of sorafenib-induced cardiomyocyte apoptosis. A combinatorial screen of the positive and negative regulators of sorafenib-induced apoptosis revealed ROS knockdown coupled with overexpression of FLT3, FGFR, PDGFR, VEGFR, or KIT as a particularly potent combination in reducing sorafenib-induced apoptosis. Network simulations of combinatorial treatment with sorafenib and the antioxidant N-acetyl cysteine (NAC) suggest that NAC may protect cardiomyocytes from sorafenib-induced apoptosis.

Mechanistic Systems Modeling to Improve Understanding and Prediction of Cardiotoxicity Caused by Targeted Cancer Therapeutics.

Tyrosine kinase inhibitors (TKIs) are highly potent cancer therapeutics that have been linked with serious cardiotoxicity, including left ventricular dysfunction, heart failure, and QT prolongation. TKI-induced cardiotoxicity is thought to result from interference with tyrosine kinase activity in cardiomyocytes, where these signaling pathways help to control critical processes such as survival signaling, energy homeostasis, and excitation-contraction coupling. However, mechanistic understanding is limited at present due to the complexities of tyrosine kinase signaling, and the wide range of targets inhibited by TKIs. Here, we review the use of TKIs in cancer and the cardiotoxicities that have been reported, discuss potential mechanisms underlying cardiotoxicity, and describe recent progress in achieving a more systematic understanding of cardiotoxicity via the use of mechanistic models. In particular, we argue that future advances are likely to be enabled by studies that combine large-scale experimental measurements with Quantitative Systems Pharmacology (QSP) models describing biological mechanisms and dynamics. As such approaches have proven extremely valuable for understanding and predicting other drug toxicities, it is likely that QSP modeling can be successfully applied to cardiotoxicity induced by TKIs. We conclude by discussing a potential strategy for integrating genome-wide expression measurements with models, illustrate initial advances in applying this approach to cardiotoxicity, and describe challenges that must be overcome to truly develop a mechanistic and systematic understanding of cardiotoxicity caused by TKIs.