The establishment and validation of the dual in vitro mitochondrial toxicity assay system
In this study, we established a cell-based in vitro mitochondrial toxicity assay system coupling a first-step screening and a second-step reverse validation. The “glucose/galactose” assay was adopted in the first-step screening, in which cells switch their energy source from glycolysis to OXPHOS. So, cells are more sensitive to mitochondrial toxicants when glucose is replaced with galactose (Kamalian et al. 2015; Orlicka-Płocka et al. 2020; Swiss and Will 2011). In this step, the human hepatocarcinoma HepG2 cell line was selected because of its super sensitivity for mitochondrial toxicity assay than other cell lines (Bugelski et al. 2000; Ma et al. 2001; Schoonen et al. 2005). In the second step, we took advantage of the human colorectal adenocarcinoma cell line HCT116 SCO2 KO, in which the mitochondrial complex IV assembly SCO2 gene is depleted. Compared to parental wild-type cells (HCT116 WT), HCT116 SCO2 KO cells are deficient in respiration and resistant to cell proliferative inhibition by the mitochondrial toxicant metformin (Sung et al. 2010; Wang et al. 2017).
Mitochondrial toxicants will decrease intracellular ATP levels in galactose medium in the first assay, while are less cytotoxic in HCT116 SCO2 KO cells in the second assay. To validate this system, we chose a well-defined set of mitochondrial toxicants, including rotenone (complex I), metformin (complex I), antimycin A (complex III), and oligomycin A (complex V), as well as cytotoxic digoxin and tamoxifen, as positive and negative control compounds, respectively. In the first assay, the IC50-ATP value, the concentration that results in a 50% reduction in ATP content, was calculated from dose-response curves for HepG2 cells cultured in glucose or galactose medium. As shown in Fig. 1A and Supplementary Table S3, the IC50-ATP values of rotenone, metformin, oligomycin A, and antimycin A were significantly reduced in galactose medium compared with that in glucose medium, while digoxin and tamoxifen did not differ significantly. According to previous publications, an IC50-ATP-Glu/IC50-ATP-Gal ratio ≥ 2 was defined as mitochondrial toxicity (Kamalian et al. 2015; Orlicka-Płocka et al. 2020; Swiss and Will 2011).
In the second assay, cytotoxicity was assessed by the SRB assay. Half-maximum proliferative concentrations, IC50 values, were calculated from dose-response curves. As expected, the positive compounds rotenone and metformin significantly inhibited cell proliferation in HCT116 WT cells compared to HCT116 SCO2 KO cells (Fig. 1B; Supplementary Table S4). The dose-response curves and the IC50 values of the cytotoxic compound digoxin were not significantly different between the two cell lines, excluding its effect on mitochondrial function (Fig. 1; Supplementary Table S3, S4). Similarly, we defined the IC50-SCO2 KO/ IC50-WT ratio ≥ 2 as mitochondrial toxicity.
Systematic evaluation of mitochondrial toxicity of TCM injections
Using the dual in vitro assay approach, we systematically evaluated the mitochondrial toxicity of TCM injections we could collect. Among the 35 TCM injections tested, Xiyanping (XYP), Dengzhanhuasu (DZHS), Shuanghuanglian (SHL) and Yinzhihuang (YZH) injections exhibited the IC50-ATP-Glu/IC50-ATP-Gal ratio ≥ 2, ranging from 2.1061 to 5.5678, implying the potential mitochondrial toxicity (Fig. 2A; Table 1). This was validated by the secondary cytotoxic assay, in which the 4 TCM injections showed greater cell proliferation inhibition in HCT116 WT cells than in HCT116 SCO2 KO cells, and all the IC50-SCO2 KO/IC50-WT ratios were ≥ 2 (Fig. 2B; Table 2). In addition, three TCM injections, including Salvianolate, Sodium aescinate, and Shugannin, were moderately toxic to the mitochondrion in the first assay, with IC50-ATP-Glu/IC50-ATP-Gal ratios between 1–2 (Table 1). Two out of the three TCM injections, except Sodium aescinate, were also defined as weak mitochondrial toxicants in the second assay (Table 2). Therefore, the dual in vitro assays exhibited high consistency.
The safety issues of the above 4 TCM injections have been widely reported. As summarized in Table 3, body as a whole-general disorders as well as skin and appendages disorders, both have allergic pathology, are the most common ADRs. The Adverse Drug Reaction Information Bulletin, issued by the National Medical Products Administration (NMPA), is the major source of information on drug safety issues in China. Due to severe allergic reactions, SHL injection was notified twice in 2005 (National Center for ADR Monitoring 2005) and 2009 (National Center for ADR Monitoring 2009), and XYP injection was notified once in 2012 (National Center for ADR Monitoring 2012). In addition, DZHS injection caused liver and kidney dysfunction (Chengfeng et al. 2020; Zhang et al. 2021). YZH and XYP injections affected the cardiovascular and nervous systems (Administration 2016; Chen et al. 2012; Chen et al. 2016). As drug-induced mitochondrial toxicity can affect multiple organs, such as liver, heart, kidney, skeletal muscle, and brain (Will and Dykens 2014), we speculate that these ADRs may be related to their effects on mitochondria.
Identification of mitochondrial toxic ingredients by molecular docking
It is of great significance to identify mitochondrial toxic ingredients from the 4 TCM injections. On the one hand, it can be inferred that other TCM products containing the same ingredients may have similar concerns. On the other hand, by minimizing the content of the harmful components, the safety of related products can be improved. To this end, we first retrieved all known ingredients of the 4 TCM injections from the TCMSP database (https://tcmsp-e.com/). A total of 622 ingredients from the 7 common herbs were found after eliminating the overlaps, including 143 in Scutellariae Radix, 53 in Artemisiae Scopariae Herba, 49 in Erigeron Breviscapus, 236 in Lonicerae Japonicae Flos, 150 in Forsythiae Fructus, 98 in Gardeniae Fructus, and 1 in Andrographis Herba (Fig. 3A). Next, we employed molecular docking to identify potential mitochondrial toxic compounds. Complex I (NADH: ubiquinone oxidoreductase), the largest enzyme of the mitochondrial respiratory chain, is the most common binding target of mitochondrial toxic drugs (Imaizumi et al. 2015; Murai and Miyoshi 2016). NADH dehydrogenase [ubiquinone] flavoprotein 1, mitochondrial (NDUFV1) is the core subunit of complex I, and its binding with compounds results in defects in the electron transfer chain (Zhu et al. 2016). In addition, by empirical judgment and Sitemap (the module of Schrödinger (Maestro Version 11.1.011) software) analysis, NDUFV1 is the best protein-ligand binding pocket in complex I. Therefore, we took NDUFV1 as an example to perform molecular docking with 622 compounds from 4 TCM injections.
After eliminating the overlaps, a total of 87 compounds showed a strong binding ability to NDUFV1 with a docking score lower than − 6.804, the value for the natural ligand flavin mononucleotide. Among the 87 compounds, there are 9 in Scutellariae Radix, 4 in Artemisiae Scopariae Herba, 6 in Erigeron Breviscapus, 31 in Lonicerae Japonicae Flos, 37 in Forsythiae Fructus, and 13 in Gardeniae Fructus (Fig. 3B; Supplementary Table S5). Baicalin and phillyrin in SHL and YZH and scutellarin in DZHS are quality control ingredients in the “Chinese Pharmacopeia 2020 Edition” (Commission 2020). In addition, baicalin, rutin, and phillyrin have been reported to be associated with various ADRs (Table 4). XYP injection is a mixture of andrographolide sulfate, which is a water-soluble medicament prepared from andrographolide extracted from Andrographis Herba through sulfonating reaction (Chong et al. 2013; Zheng et al. 2016). Therefore, we further analyzed the binding modes of interaction between these compounds and the target protein. From the 3D crystal structure of the small molecule-NDUFV1 complex, it can be seen that the distance between the amino acid residue and the small molecule was less than 5 Å (Fig. 3C). The docking scores of scutellarin, rutin, phillyrin, and baicalin were − 7.456 kcal/mol, -7.953 kcal/mol, -7.471 kcal/mol, and − 8.459 kcal/mol, respectively (Fig. 3C). So, these compounds can bind strongly to NDUFV1, suggesting their potential mitochondrial toxicity. In line with our findings, it has been recently reported that anthracene-9,10-dione and phthalaldehyde, a series of flavonoid derivative substructures, including baicalin, scutellarin, and rutin, were mitochondrial toxic identified by machine learning (Zhao et al. 2021).
In vitro validation of mitochondrial toxicity of identified TCM ingredients
To validate the findings of molecular docking, these compounds were subjected to the dual in vitro mitochondrial toxicity assay. As expected, all compounds showed an IC50-ATP-Glu/IC50-ATP-Gal ratio ≥ 2 (Fig. 4A; Table 5) and an IC50-SCO2 KO/IC50-WT ratio ≥ 2 (Fig. 4B; Table 6), which demonstrated to be mitochondrial toxic. As summarized in Table 4, baicalin and rutin are the main components in SHL injection causing the anaphylactoid reaction (Wang et al. 2020; Zhang et al. 2017). Baicalin, also present in YZH injection, has been reported to induce kidney damage and renal fibrosis (Cai et al. 2017). Rutin, a plant pigment contained in SHL and YZH injections, is generally safe, but overdose can cause cardiovascular and neurological disorders (Wilson 2017). Although phillyrin in SHL injection and scutellarin in DZHS injection were identified as potential mitochondrial toxicants by molecular docking analysis (Fig. 3C) and our dual in vitro approach (Fig. 4A, B), there are few reports of their ADRs (Tabel 4).
XYP injection is a mixture of andrographolide sulfate made from andrographolide extracted from Andrographis Herba through sulfonation reaction to increase water solubility, bioavailability, and stability (Loureiro Damasceno et al. 2022; Zou et al. 2022). Andrographolide is the main component of XYP injection. However, its docking score with NDUFV1 was not significant (-6.101 kcal/mol) (Fig. 3C), which led us to speculate that its mitochondrial toxicity may not be due to its binding to complex I. To our surprise, neither the IC50-ATP-Glu/IC50-ATP-Gal ratio (Fig. 4A; Table 5) nor the IC50-SCO2 KO/IC50-WT ratio (Fig. 4B; Table 6) of andrographolide was ≥ 2, indicating lower mitochondrial toxicity compared with XYP injection (Fig. 2; Table 1, 2). This discrepancy indicates that the mitochondrial toxicity of XYP injection may come from uncharacterized components during the extraction process or during the sulfonation reaction, which needs further investigation.
Effects of mitochondrial toxic TCM injections and ingredients on the respiration of HepG2 cells
Measurement of oxygen consumption rate (OCR) has long been used as the gold standard for evaluating drug-induced mitochondrial toxicity (Meyer et al. 2018). To validate the identified mitochondrial toxic TCM injections and ingredients, we performed OCR assays in HepG2 cells using a fluorescence lifetime micro-oxygen monitoring system. As expected, positive mitochondrial toxicants, such as rotenone, oligomycin A and metformin, significantly inhibited OCR levels (Fig. 5A). Similar effects were observed for the identified TCM injections and their ingredients (Fig. 5B, C). Taken together, these results suggest that XYP, DZHS, SHL, and YZH injections have potential mitochondrial toxicity, which may be related to their ADRs. Their adverse effects on mitochondria are due at least in part to the toxic components baicalin, rutin, scutellarin, and phillyrin.