Next Article in Journal
The Oncolytic Adenovirus XVir-N-31, in Combination with the Blockade of the PD-1/PD-L1 Axis, Conveys Abscopal Effects in a Humanized Glioblastoma Mouse Model
Next Article in Special Issue
Corydalis saxicola Bunting: A Review of Its Traditional Uses, Phytochemistry, Pharmacology, and Clinical Applications
Previous Article in Journal
Recent Advances in Molecular Improvement for Potato Tuber Traits
Previous Article in Special Issue
Investigation of Chemical Constituents of Eranthis longistipitata (Ranunculaceae): Coumarins and Furochromones
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 23
Issue 17
10.3390/ijms23179986
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Silibinin Suppresses the Hyperlipidemic Effects of the ALK-Tyrosine Kinase Inhibitor Lorlatinib in Hepatic Cells

by
Sara Verdura
1,2,
José Antonio Encinar
3,
Salvador Fernández-Arroyo
4,5,
Jorge Joven
4,5,
Elisabet Cuyàs
1,2,*,
Joaquim Bosch-Barrera
2,6,7 and
Javier A. Menendez
1,2,*
1
Metabolism and Cancer Group, Program Against Cancer Therapeutic Resistance (ProCURE), Catalan Institute of Oncology, 17007 Girona, Spain
2
Girona Biomedical Research Institute (IDIBGI), Salt, 17190 Girona, Spain
3
Institute of Research, Development and Innovation in Biotechnology of Elche (IDiBE) and Molecular and Cell Biology Institute (IBMC), Miguel Hernández University (UMH), 03207 Elche, Spain
4
Department of Medicine and Surgery, Universitat Rovira i Virgili, 43204 Reus, Spain
5
Unitat de Recerca Biomèdica (URB-CRB), Hospital Universitari de Sant Joan, Institut d’Investigació Sanitaria Pere Virgili, Universitat Rovira i Virgili, 43204 Reus, Spain
6
Medical Oncology, Catalan Institute of Oncology, 17007 Girona, Spain
7
Department of Medical Sciences, Medical School, University of Girona, 17071 Girona, Spain
*
Authors to whom correspondence should be addressed.
Int. J. Mol. Sci. 2022, 23(17), 9986; https://doi.org/10.3390/ijms23179986
Submission received: 5 August 2022 / Revised: 27 August 2022 / Accepted: 30 August 2022 / Published: 1 September 2022
(This article belongs to the Special Issue Frontiers in Plant Pharmacology: The Role of Medicinal Plants in Health Systems)
Browse Figures
Review Reports Versions Notes

Abstract

:
The third-generation anaplastic lymphoma tyrosine kinase inhibitor (ALK-TKI) lorlatinib has a unique side effect profile that includes hypercholesteremia and hypertriglyceridemia in >80% of lung cancer patients. Here, we tested the hypothesis that lorlatinib might directly promote the accumulation of cholesterol and/or triglycerides in human hepatic cells. We investigated the capacity of the hepatoprotectant silibinin to modify the lipid-modifying activity of lorlatinib. To predict clinically relevant drug–drug interactions if silibinin were used to clinically manage lorlatinib-induced hyperlipidemic effects in hepatic cells, we also explored the capacity of silibinin to interact with and block CYP3A4 activity using in silico computational descriptions and in vitro biochemical assays. A semi-targeted ultrahigh pressure liquid chromatography accurate mass quadrupole time-of-flight mass spectrometry with electrospray ionization (UHPLC-ESI-QTOF-MS/MS)-based lipidomic approach revealed that short-term treatment of hepatic cells with lorlatinib promotes the accumulation of numerous molecular species of cholesteryl esters and triglycerides. Silibinin treatment significantly protected the steady-state lipidome of hepatocytes against the hyperlipidemic actions of lorlatinib. Lipid staining confirmed the ability of lorlatinib to promote neutral lipid overload in hepatocytes upon long-term exposure, which was prevented by co-treatment with silibinin. Computational analyses and cell-free biochemical assays predicted a weak to moderate inhibitory activity of clinically relevant concentrations of silibinin against CYP3A4 when compared with recommended (rosuvastatin) and non-recommended (simvastatin) statins for lorlatinib-associated dyslipidemia. The elevated plasma cholesterol and triglyceride levels in lorlatinib-treated lung cancer patients might involve primary alterations in the hepatic accumulation of lipid intermediates. Silibinin could be clinically explored to reduce the undesirable hyperlipidemic activity of lorlatinib in lung cancer patients.

1. Introduction

The small macrocyclic compound lorlatinib (product name PF-06463922) is a third-generation anaplastic lymphoma kinase (ALK) tyrosine kinase inhibitor (TKI) originally developed to inhibit ALK mutant forms causing resistance to first- and second-generation ALK-TKIs (e.g., crizotinib, ceritinib, alectinib or brigatinib) in patients with advanced non-small-cell lung cancer (NSCLC) [1,2,3,4,5,6]. Consistent with its broad ALK mutational coverage and optimized central nervous system penetration through the blood–brain barrier, lorlatinib has shown substantial systemic and intracranial activity both in treatment-naïve patients and in patients with relapse after first- and second-generation ALK TKIs [7,8,9,10,11,12]. It has been approved for the treatment of patients with ALK-positive metastatic NSCLC whose disease has progressed after alectinib or ceritinib as the first ALK-TKI, or after crizotinib and at least one other ALK-TKI [13].
Lorlatinib has a unique treatment-related adverse event profile compared with other ALK-TKIs that is characterized by the frequent incidence of hypercholesteremia and hypertriglyceridemia, which have been found to occur in 61% and 82% of lorlatinib-treated patients included in a pooled safety analysis (n = 295) [14,15]. Results from the global phase 2 study B7461001 (NCT01970865) revealed that grade 3 or 4 elevations in total cholesterol and triglycerides both occurred at a frequency of 16% [2]. Hyperlipidemia—comprising the cluster terms hypercholesterolemia and hypertriglyceridemia—is the most commonly reported adverse reaction to lorlatinib, which usually arises within the first few weeks of treatment (median time to onset, 15 days). Lorlatinib-induced hyperlipidemia is a common reason for temporary (but not permanent) discontinuation and dose-reduction of lorlatinib in the clinical setting, which necessitates the rapid initiation of lipid-lowering medications. Accordingly, the majority (>80%) of patients receive at least one lipid-lowering agent within 3 weeks of the first lorlatinib dose (median time to onset, 20 days), and ~25% of patients who received a lipid-lowering agent for hypercholesterolemia or hypertriglyceridemia require the addition of a second agent. While statins are the most commonly prescribed lipid-lowering agents to clinically manage lorlatinib-related hyperlipidemia, the choice and dosing of statins should be guided by the differential metabolism of the CYP450 pathway [16]. As lorlatinib is metabolized primarily by CYP3A4, the so-called 3A4 statin substrates atorvastatin, lovastatin, and simvastatin should not be co-administered as lipid-lowering agents in lorlatinib-treated patients. The non-3A4 statin substrates rosuvastatin, pravastatin, and pitavastatin are not significantly metabolized by CYP3A4 enzymes, and so the potential for adverse statin-lorlatinib interactions is notably reduced. Rosuvastatin is currently the only lipid-lowering agent recommended for lorlatinib-associated hyperlipidemia based on its low involvement with CYP450 enzymes that are able to interact with lorlatinib [17].
Given the improvements in prognosis and survival of patients with ALK translocations, a better understanding of the primary causes underlying lorlatinib-triggered hypercholesterolemia/hypertriglyceridemia might inform new therapeutic strategies to prevent or manage the undesirable lipid-modifying activity of lorlatinib in patients with ALK-positive NSCLC. Here, we tested the hypothesis that lorlatinib might directly promote hypercholesteremic and hypertriglyceridemic effects in human hepatocytes. First, using non-targeted lipidomics comprising ultrahigh pressure liquid chromatography accurate-mass quadrupole time-of-flight mass spectrometry with electrospray ionization (UHPLC-ESI-QTOF-MS/MS) [18], and imaging of neutral lipids, we explored the ability of lorlatinib to alter the lipidome of hepatoma tissue-derived Huh-7 and HepG2 cells [19,20,21,22,23,24], which were employed as substitutes for primary hepatocytes. Second, we investigated whether silibinin––a flavonolignan that functions as a hepatoprotectant in patients with acute and chronic liver injury [25,26,27,28,29,30]––might prevent the lipid-modifying activity of lorlatinib in hepatocytes. To predict potentially relevant drug–drug interactions if silibinin were used to clinically manage lorlatinib-associated hyperlipidemia, we finally explored the capacity of silibinin to interact with and block CYP3A4 activity in comparison with currently employed statins. We now provide first-in-class evidence that the hyperlipidemic effects of lorlatinib might involve, at least in part, the induction of an increased content of cholesterol and triglycerides in hepatic cells, which can be prevented or reversed by silibinin.

2. Results

2.1. Hepatic Cells Treated with Lorlatinib Accumulate Cholesteryl Esters and Triglycerides

We first carried out MTT-based viability tests to evaluate the sensitivity of human hepatoma-derived Huh-7 cells, which closely mimic primary hepatocytes for drug metabolism and toxicity studies, to lorlatinib. As shown in Figure 1A, the IC50 of lorlatinib for Huh-7 cells (22 ± 5 μmol/L) was similar to that for ALK-negative A549 NSCLC cells (19 ± 5 μmol/L), but was tens or even hundreds of thousands higher than that obtained for H2228 and H3122 NSCLC cell lines (0.9 ± 0.2 and 0.003 ± 0.0002 μmol/L, respectively) harboring ALK rearrangements (variants 3a/b and 1, respectively) [31,32].
We then questioned whether lorlatinib used at concentrations efficacious against ALK-rearranged lung cancer cells but non-toxic to ALK-negative cells could promote substantial lipidome changes in Huh-7 cells (Figure 1B; Table 1). Accordingly, Huh-7 cells were treated with 1 μmol/L lorlatinib for 48 h, and the lipidome was characterized in a semi-targeted approach using a UHPLC-ESI-QTOF-MS/MS method, in which we matched mass spectra with commercial standards (SPLAH mix), databases such as METLIN-PCDL and Lipid MAPS (accessed on 1 July 2022), as well as corroboration of tentative lipid species upon obtaining their MS/MS spectrum. Representative chromatograms of lipidomic analyses in Huh-7 cells can be found in Figure S1.
We identified and quantified more than 100 lipid species belonging to six different families: diglycerides (DG), triglycerides (TG), phosphatidylcholines (PC), lysosophosphatidylcholines (LPC), phosphatidylethanolamines (PE), sphingomyelins (SM), and cholesterol esters (ChoE). Lorlatinib treatment was found to significantly increase several species of ChoE, numerous molecular species of long-chain TG, and several DG, LPC, and PC species (Figure 1C; Figure S2; Tables S1–S7).

2.2. Silibinin Fully Protects the Steady-State Lipidome of Hepatic Cells against the Hyperlipidemic Effects of Lorlatinib

Volcano plots and Venn diagrams of differentially accumulated lipid species revealed that, when used as a single agent, silibinin (100 μmol/L, 48 h) treatment significantly lowered the baseline abundance of many ChoE species in Huh-7 cells. Similarly, silibinin co-treatment completely blunted the significant lorlatinib-induced rise in numerous ChoE and TG/DG species in Huh-7 cells. Indeed, cells co-treated with silibinin and lorlatinib showed notable decreases in the abundance of various ChoE species, and this was accompanied by the compensatory generation of PC species (Figure 1C; Figure S2). Of note, no significant changes in any lipid species were observed when comparing silibinin-treated hepatocytes with those co-treated with silibinin and lorlatinib (Figure 1C; Figure S2).
Colorimetric enzymatic assays confirmed the ability of silibinin to inhibit the accumulation of cholesterol esters and triglycerides in lorlatinib-treated Huh-7 cells (Figure 1D).

2.3. Silibinin Prevents the Lorlatinib-Induced Chronic Accumulation of Neutral Lipids in Hepatic Cell

To evaluate the potential clinical relevance of silibinin’s ability to acutely suppress lorlatinib-induced hyperlipidemia in hepatic cells, we designed a protocol of repeated daily exposure to lorlatinib and/or silibinin using a silibinin concentration (10 μmol/L) that can be practically achieved in a clinical setting. We first employed a conventional index of fatty liver disease––the intracellular accumulation of neutral lipids (steatosis)––to confirm the ability of lorlatinib to trigger hypertriglyceridemia and hypercholesteremia. The results of Oil red O- and LipidTOX™-based staining of neutral lipids (cholesteryl esters and triglycerides) in lorlatinib-treated Huh-7 cells cultured with or without silibinin are shown in Figure 2A. As expected, in the lorlatinib only-treated arm, hepatic cells accumulated large amounts of lipid droplet-like neutral lipids, whereas this pattern was largely lost in the presence of clinically relevant concentrations of silibinin, which notably decreased the quantity and size of the lipid droplets formed by exposure to lorlatinib.
To confirm these findings, we employed an independent in vitro model of fatty liver disease involving non-tumorigenic HepG2 cells, which are routinely employed as an inducible fatty liver cell model to screen for therapeutic drugs to eliminate lipid accumulation. The results of Oil red O staining of lorlatinib-treated HepG2 cells cultured with or without silibinin are shown in Figure 2B. Lorlatinib-treated HepG2 cells contained a greater number of larger lipid droplets than untreated cells. By contrast, there were markedly fewer and smaller neutral lipid depots in the presence of silibinin even at the highest tested concentration of lorlatinib.

2.4. Silibinin Does Not Overlap the Binding Mode of Lorlatinib to Cytochrome P450 3A4 (CYP3A4) In Silico

To predict clinically relevant drug–drug interactions if silibinin were used to clinically manage lorlatinib-induced hyperlipidemic effects in hepatic cells, we explored the capacity of silibinin to interact with and block the lorlatinib-metabolizing CYP3A4 isoenzyme using in silico computational descriptions and in vitro biochemical assays. First, to shed some light on the capacity of silibinin to interact with CYP3A4, we performed a computational study comparing the predicted molecular interactions between silibinin and CYP3A4 with those of lorlatinib and ritonavir, which is the most potent CYP3A4A inhibitor in clinical use. As silibinin is almost a 1:1 mixture of the diastereomers A and B, we performed classical molecular dockings of silibinin A and B against the 3D crystal structure of human CYP3A4 bound to an inhibitor (PDB:7KVS) [33] (Figure 3).
To characterize the potential binding of silibinin A and B to the large catalytic cavity of CYP3A4, which allows CYP3A4 to bind a broad range of substrates as compared with other CYP isoenzymes, a total of 75 flexible docking runs were set around the putative binding sites of lorlatinib and ritonavir using AutoDock Vina implemented by YASARA. The resulting binding energies with the docking simulations of silibinin A (−10.323 kcal/mol) and B (−9.561 kcal/mol) were similar to those obtained with lorlatinib (−9.876 kcal/mol) and ritonavir (−9.962 kcal/mol). However, closer inspection of the different conformations revealed that the predicted interaction of silibinin A and B with the active site cavity of CYP3A4 shared only 39–46% of the 12 amino acid residues involved in the lorlatinib binding mode, and 45–59% of the 22 participating amino acids involved in the ritonavir binding mode (Table 2).
Unlike lorlatinib and ritonavir, neither silibinin A nor silibinin B were predicted to interact with PHE304, which is critically important for CYP3A4 substrate specificity and catalytic capacity [34,35,36]. Moreover, neither silibinin A nor silibinin B were predicted to interact with a majority of the hydrophobic pocket near the I-helix (the so-called Phe-1 site) including PHE108, LEU210, LEU211, PHE241, ILE301, and PHE304 [37,38] (Figure 3). Silibinin is, therefore, predicted to interact with the catalytic site of CYP3A4 through a unique binding mode lacking the involvement of the gatekeeper residues regulating the access of CYP3A4 substrates to the heme group, suggesting a highly improbable interaction between silibinin and lorlatinib at the CYP3A4 active site.

2.5. Silibinin Is a Weak Inhibitor of the Lorlatinib-Metabolizing Cytochrome P450 3A4 (CYP3A4) Isoenzyme

The CYP3A4 inhibition propensity of silibinin in comparison with that of recommended (non-3A4 substrates such as rosuvastatin) and non-recommended (3A4 substrates such as simvastatin) statins to clinically manage lorlatinib-associated dyslipidemia was initially predicted using the DL-CYP prediction server [39], a free web tool to evaluate the tendency of small molecules to inhibit five major CYP450 isoforms (1A2, 3C9, 2C19, 2D6, and 3A4). Using this deep autoencoder multitask neural network trained on more than 13,000 compounds, simvastatin showed the highest predicted value (0.93), which ranked very close to that obtained with the positive control ritonavir (0.97) (Table 3). The tendency of CYP3A4 inhibition by rosuvastatin was predicted to be as low as 0.2, whereas silibinin had no predicted inhibitory effect in silico (0.06) (Table 3).
To facilitate the understanding of the predicted tendencies, we performed MD simulations for each of the CYP3A4-drug complexes. This approach considers the protein flexibility at the target-binding site during the molecular recognition process, allowing the user to confirm the kinetic stability and validate the binding poses obtained by docking. The CYP3A4 protein backbone root mean square deviation (RMSD) plots of the drugs’ heavy atoms, measured after superimposing CYP3A4 onto its reference structure during the MD simulation, were prepared in parallel. The approach is summarized in Figure 4 and shows as follows: the best poses of ritonavir, silibinin A, silibinin B, simvastatin, and rosuvastatin coupled to the catalytic cavity of CYP3A4 before (0 ns) and after (100 ns) the MD simulation, the time evolution of RMSD relative to the initial structure of CYP3A4 in the absence and presence of drugs, the binding free energy calculations under the MM/PBSA approximation from the entire MD simulation trajectory of 100 ns (or last 30 ns), and the identification of amino acid residues participating in the drug-CYP3A4 binding pocket. The MM/PBSA parameters, which estimate the free energy of the binding of small ligands to biological macromolecules, allows the correct ranking of drug candidates based on the fact that more positive energies indicate the occurrence of stronger binders (e.g., “strong” > 25 kcal/mol versus “weak” < 25 kcal/mol). Likewise, the results predicted a high affinity of ritonavir and simvastatin for the catalytic site of CYP3A4, reaching MM/PBSA values of ~70 and 40 kcal/mol, respectively (Figure 4). Silibinin A mimicked rosuvastatin as a predicted weak inhibitor of CYP3A4, showing MM/PBSA values of <20 kcal/mol. Silibinin B was predicted to encounter steric hindrance at the catalytic activity of CYP3A4 given its notable −60 kcal/mol MM/PBSA value (Figure 4).
The inhibitory effects of graded concentrations of simvastatin, rosuvastatin, silibinin, and Eurosil 85®/Euromed –a standardized pharmaceutical preparation containing 60% silibinin with an enhanced bioavailability (>80%) that is commonly used in clinical research [40,41]—on CYP3A4 activity were finally evaluated using the Vivid® CYP450 Red screening assay, which employs microsomes prepared from insect cells (baculosomes) expressing a human CYP450 isoenzyme (CYP3A4 in this case) and human cytochrome P450 reductase, and utilizing the benzyloxy–methyl–resofurin (BOMR) red substrate. The Vivid® BOMR substrate is a blocked dye that yields a minimal fluorescence signal until metabolized by CYP3A4 into products that are highly fluorescent in aqueous solutions (Figure 5, left panel). Simvastatin strongly inhibited the enzymatic activity of CYP3A4, exhibiting an IC50 value as low as 1.8 ± 0.9 μmol/L. Up to 10-fold higher concentrations of silibinin (IC50 = 16.7 ± 8.9 μmol/L) and Eurosil 85®/Euromed (12.1 ± 1.2 μmol/L) were necessary to achieve the half-inhibitory concentration of CYP3A4 activity (Figure 5, right panel). Rosuvastatin likewise failed to show any significant inhibitory effect on CYP3A4 activity even when tested at very high concentrations (>300 μmol/L).

3. Discussion

Elevated plasma levels of cholesterol and triglycerides are a common adverse effect in the majority of patients with ALK-positive NSCLC treated with lorlatinib [14,15,17]. The molecular mechanisms of lorlatinib-related lipid disorders are likely multifactorial and are currently unclear. Because hyperlipidemia is a classic feature of nephrotic syndrome, the fact that lorlatinib is apparently the only ALK-TKI reported to induce acute kidney injury involving renal cyst formation or progression suggests that minimal change disease (nephrotic syndrome) should be viewed as a possible secondary cause of lorlatinib-induced hyperlipidemia [42,43,44,45]. Here, we provide pre-clinical evidence that the elevated plasma cholesterol and triglyceride levels in lorlatinib-treated lung cancer patients might be due, at least in part, to direct alterations in the hepatic accumulation of lipid intermediates. Moreover, we reveal that the flavonolignan silibinin has the capacity to protect hepatic cells against the lipid-modifying activity of lorlatinib (Figure 6).
Lorlatinib is rapidly absorbed to achieve peak plasma concentrations of >2 μmol/L 0.5–4 h after a single 100 mg oral dose and is widely distributed in different tissues [46,47], with the highest concentration found in the liver [48]. The hepatic accumulation of lorlatinib can be explained in terms of its low metabolism in this organ, with only 12% of lorlatinib metabolized during first-pass hepatic metabolism or pre-systemic metabolism. We show that acute exposure to a clinically relevant concentration of lorlatinib (1 μmol/L) promotes a rapid and significant buildup of cholesteryl esters and tri-/di-glycerides in hepatic cells. We also found that longer exposure to lorlatinib drives the hallmarks of hepatic steatosis, such as the abnormal and excessive accumulation of neutral lipids in liver cells. While our study was not designed to investigate the precise molecular mechanisms involved in the hepatic alterations in cholesterol and triglycerides homeostasis in response to lorlatinib, it is reasonable to suggest that there is overlap in the targeted pathways for ALK-rearranged NSCLC growth and those that regulate hepatic function, leading to on-target or closely related off-target hepatic toxicity. It will be critical to understand the mechanistic pathways by which lorlatinib induces a dysregulation of lipid and cholesterol metabolism in hepatic cells. Because one of the major mechanisms of drug-induced liver injury (including TKIs such as lorlatinib) centrally involves mitochondrial perturbation and dysfunction [49,50], mitochondria might be viewed as a probable target through which lorlatinib alters lipid metabolism in human hepatocytes. Particularly, drugs inducing impairment of fatty acids oxidation are known to lead to lipid accumulation (steatosis) and steatohepatitis [50], thus suggesting that a plausible hypothesis for the hyperlipidemic mechanism of action of lorlatinib might involve a decrease in lipid catabolism in response to TKI-related mitochondrial dysfunction in hepatic cells, leading to a disbalance between fatty acid β-oxidation activity and cholesteryl esters/triglycerides synthesis [51]. Future studies will be needed to address this possibility. We acknowledge, however, that our pre-clinical findings in cell culture do not preclude the possibility that other mechanisms involving changes in peripheral tissues (e.g., adipose tissue, muscle), lipid uptake by the liver, or intestinal lipid absorption can contribute to the lorlatinib-induced lipid metabolic effects.
The flavonolignan silibinin has been shown to significantly temper the accumulation of hepatic and biliary triglycerides and cholesterol in cultured cells and animal models, and in patients with type 2 diabetes and chronic liver disease. We show here that silibinin suppresses the lipid-modifying activity of lorlatinib in human hepatic cells. Regardless of the mechanism through which silibinin protects the lipidome of hepatic cells against lorlatinib (e.g., by increasing the shift of fatty acids from triglycerides towards phospholipids and/or increasing the endogenous cholesterol conversion to bile acids), silibinin might represent an idoneous lipid-lowering agent in new treatment regimens for patients with NSCLC at the highest risk of developing brain metastases, such as those continuously exposed to ALK-TKIs [52]. As lorlatinib is effective against central nervous system metastasis even in patients pretreated with first- and second-generation ALK-TKIs [10], the known capacity of silibinin to provide better control of brain disease and allow a higher proportion of patients with NSCLC to receive additional lines of treatment [53,54] could be exploited by utilizing its ability to ameliorate lorlatinib-induced hyperlipidemia.
The low hepatic extraction ratio (<30%) of “low” clearance drugs such as lorlatinib can drive a change in the plasma counterpart when co-administered with another drug that is a metabolic inhibitor or inducer. Because lorlatinib is metabolized primarily by the CYP3A4 isoenzyme [55,56], lipid-lowering agents that do not interact with CYP3A4, such as rosuvastatin, should be co-administered in lorlatinib-treated patients. Our multi-faceted approach involving computational modeling, neural network-based prediction models to classify CYP inhibition, and in vitro validation using recombinant human CYP3A4, confirms and further extends the notion that exposure to standard doses of silibinin-containing milk thistle preparations leading to relatively low peak systemic concentrations (1–10 μmol/L) are expected to be accompanied by low-risk for an inhibitory interaction of the first-pass clearance of lorlatinib at the hepatic level [57]. Nonetheless, because notably higher intestinal concentrations of silibinin can be achieved when using gram doses of certain silibinin formulations that have been tested in patient populations [57,58], it cannot be excluded that clinically relevant inhibitory interactions with the lorlatinib-metabolizing CYP3A4 isoenzyme might occur at the intestinal level. Indeed, a limitation of our study is that the capacity of silibinin to protect hepatocytes from the hyperlipidemic activity of lorlatinib was observed in two hepatoma-derived cell lines (Huh-7 and HepG2) that naturally express low expression of drug metabolizing enzymes such as CYP3A4 [59,60,61]. We therefore acknowledge that our cell-free evaluations of the in silico and in vitro metabolic interactions between lorlatinib and silibinin (versus statins) using computational descriptions (with crystal structures of CYP3A4) and biochemical assays (with recombinant CYP3A4) cannot be mechanistically extrapolated to the cell-based experiments showing the ability of silibinin to suppress the hypercholesterolemic and hypertriglyceridemic effects of lorlatinib in Huh-7 and HepG2 cells. We failed to observed significant toxic effects when combining lorlatinib with silibinin concentrations as high as 100 μmol/L to carry out lipidomic studies with Huh-7 cells growing confluent for more than 1 month, thereby ensuring an increased expression and catalytic activity of CYP3A4 [62,63]. However, the extrapolation of the in vitro results to accurately predict the toxicity as it would occur in vivo might require the usage of the newly developed human hepatoma cell lines such as HepaRG [64,65,66], which shows stable expression of liver-specific functions over a long period in culture, including high levels of CYP3A4.

4. Materials and Methods

4.1. Materials

Silibinin was purchased from Sigma-Aldrich (St. Louis, MO, USA; Cat. #S0417). Simvastatin (MK 733) and rosuvastatin (ZD4522) were obtained from Selleck Chemicals LLC (Houston, TX, USA; Cat. #S1792 and S2169, respectively). Eurosil85/Euromed was kindly provided by Meda Pharma S.L. (Barcelona, Spain). All reagents were dissolved in sterile dimethylsulfoxide to prepare 10 mmol/L stock solutions, which were stored in aliquots at −20°C until use. Working concentrations were diluted in culture medium prior each experiment.
Grade MS methanol (MeOH) and 2-propanol, ammonium formate, formic acid, acetic acid, methyl tert-buthyl ether (MTBE) and standards for calibration curves were purchased from Sigma-Aldrich. Water (Milli-Q grade) was obtained from a Milli-Q integral water purification system (Millipore Corp., Burlington, MA, USA). For internal standards, SPLASH Lipidomix was purchased from Avanti Polar lipids (Alabaster, AL, USA; Cat. #330707).

4.2. Cell Lines

Huh-7 and HepG2 cells were kindly provided by Dr. Jose Manuel Fernández-Real in our institution (IDIBGI), and were cultured in complete Dulbecco′s modified Eagle′s medium (Gibco/Invitrogen, Carlsbad, CA, USA) supplemented with 10% fetal bovine serum, 2 mmol/L L-glutamine, and 100 IU/mL penicillin-streptomycin (all from Gibco/Invitrogen, Carlsbad, CA, USA).

4.3. Cell Viability Assays

Cell viability effects of lorlatinib were determined using the colorimetric MTT (3-4,5-dimethylthiazol-2-yl-2,5-diphenyl-tetrazolium bromide) reduction assay. Dose–response curves to graded concentrations of lorlatinib were plotted as a percentage of the control cell absorbance, which was obtained from control cells containing the vehicle (DMSO v/v) processed simultaneously. For each treatment, cell viability was evaluated using the following equation: (OD570 of the treated sample/OD570 of the untreated sample) × 100. Sensitivity to agents was expressed in terms of the concentrations required for a 50% (IC50) reduction in cell viability. Since the percentage of control absorbance was considered to be the surviving fraction of cells, the IC50 values were defined as the concentration of drug that produced 50% reduction in control absorbance (by interpolation).

4.4. Non-Targeted Lipidomics

Lipid extraction, preparation of external calibration curves, and quality control approaches were carried out as described [18]. Samples (2 μL) were injected directly into a 1290 Infinity UHPLC system coupled with a dual Agilent jet stream ESI source to a 6550 QTOF-MS instrument (Agilent Technologies, Santa Clara, CA, USA). UHPLC-ESI-QTOF-MS conditions, data analysis and indirect quantification of selected lipid species using MassHunter Qualitative Analysis B.07.00 software (Agilent Technologies), and corroboration of the tentative characterization of compound match entities using the Lipid Maps database (www.lipidmaps.org, accessed on 1 July 2022) and targeted MS/MS acquisition on LC-QTOF-MS instrument, have been previously described in detail [18].
Differentially expressed lipid species from untargeted lipidomic data were analyzed by generating volcano plots indicating statistical significance versus magnitude of fold-change for each lipid species (cutoff of Benjamini–Hochberg corrected p-value < 0.05 and log2 fold-change > 1 or <1).

4.5. Triglyceride and Cholesterol Quantification

Quantification of cholesteryl esters and triglycerides in Huh-7 cells cultured in the absence/presence of lorlatinib and/or silibinin was performed using cholesterol (Cat. #MAK043) and triglyceride (Cat. #MAK266) quantification kits (Sigma-Aldrich).

4.6. Accumulation of Neutral Lipids

For assessment of neutral lipid formation, cells were washed in PBS, fixed with 4% paraformaldehyde for 30 min, washed in distilled water, and then briefly incubated in 60% isopropanol. Following equilibration, cells were stained with a working solution (60%) of Oil Red O for 20 min. An Oil Red O stock solution was prepared by dissolving 0.35% w/v Oil Red O in 100 mL 100% isopropanol. HCS LipidTOX™ Green Neutral Lipid Stain (Invitrogen; Cat. #H34475) was employed as a second staining cocktail for visualization of intracellular neutral lipids. Formaldehyde fixation and preparing/using the labeling solution was carried out as per the manufacturer’s instructions.

4.7. Prediction of CYP450 Inhibition In Silico

The CYP450 inhibition propensity of silibinin, simvastatin, and rosuvastatin was predicted with the DL-CYP Prediction Server (http://www.pkumdl.cn:8000/deepcyp/home.php, accessed on 1 July 2022), a free web tool to evaluate the tendency of small molecules to inhibit five CYP450 major isoforms, namely 1A2, 2C9, 2C19, 2D6, and 3A4 [39]. All structure-data files (sdf) for tested compounds were downloaded from PubChem (National Center for Biotechnology Information) and used as the inputs to predict their level of inhibition against human CYP450 isoforms. The CYP3A4 inhibitor ritonavir was used as a positive control to predict the inhibition of CYP3A4 activity. The results were expressed as values between 0 and 1.

4.8. Docking Calculations, Molecular Dynamics Simulations, and Binding Free Energy Analysis

Docking calculations, molecular dynamics (MD) simulations, and molecular mechanics Poisson–Boltzmann surface area (MM/PBSA) calculations to determine the alchemical binding free energy of silibinin, simvastatin, and rosuvastatin against CYP3A4 7KVS were performed using procedures described in previous works from our group [67,68,69,70,71,72,73]. To perform the docking studies with AutoDockVina (v1.1.2, San Diego, CA, USA), the CYP3A4 7KVS crystal structure was transformed to the PDBQT format, including the atomic charges and atom-type definitions. These preparations were performed using the AutoDock/Vina plugin with scripts from the AutoDock Tools package [74]. YASARA dynamics v19.9.17 (Vienna, Austria) was employed for all MD simulations with the AMBER14 force field. All simulation steps were run using a pre-installed macro (md_run.mcr) within the YASARA suite. Data were collected every 100 ps during 100 ns. The MM/PBSA calculations of solvation binding energy were calculated using the YASARA macro md_analyzebindenergy.mcr, with more negative values indicating instability. MM/PBSA was implemented with the YASARA macro md_analyzebindenergy.mcr to calculate the binding free energy with solvation of the ligand, complex, and free protein, as previously described [67,68,69,70,71,72,73]. All of the figures were prepared using PyMol 2.0 software and all interactions were detected using the protein-ligand interaction profiler (PLIP) algorithm [75].

4.9. CYP3A4 Inhibition

The characterization of the CYP3A4 inhibitory potency of silibinin, simvastatin, and rosuvastatin was outsourced to the SelectScreen P450 Profiling Service (ThermoFisher Scientific, Waltham, MA, USA) using CYP3A4-transfected baculosomes from the Vivid® CYP450 Red assay platform (Cat. No. P2856). The stringent validation process ensures the highest quality data possible for the 10-point titration, which was performed in duplicate (n = 3). An additional two control wells for each compound concentration tested are conducted to detect autofluorescent compound interference. Strict quality control protocols ensure that any assay results not meeting set specifications will be automatically repeated, this includes performing IC50 assays with a known inhibitor for each P450 on a per plate basis to ensure the validity of the data. To calculate percent inhibition due to presence of test compounds (i.e., silibinin, simvastatin, and rosuvastatin) or positive inhibition (i.e., ketoconazole for CYP3A4, which enables the subtraction of background during data analysis) we used the following equation:
%   inhibition = ( 1   –   X – B A – B ) × 100 %
where X is the fluorescence intensity of test compound, A is the fluorescence intensity observed in the absence of inhibitor (DMSO-only control that accounts for possible inhibition caused by introduction of test compounds originally dissolved in organic solvents), and B is the fluorescence intensity observed in the presence of positive inhibition control.

4.10. Statistical Analysis

All cell-based observations were confirmed by at least three independent experiments performed in triplicate for each cell line and for each condition. Data are presented as mean ± SD. Two-group comparisons were performed using Student′s t-test for paired and unpaired values. Comparisons of means of ≥3 groups were performed by ANOVA, and the existence of individual differences, in case of significant F values with ANOVA, was tested by Scheffé′s multiple contrasts: p-values < 0.01 and <0.001 were considered to be statistically significant (denoted as * and **, respectively). All statistical tests were two-sided.

5. Conclusions

Hyperlipidemia, which encompasses hypercholesteremia and hypertriglyceridemia, is a unique adverse side effect of lorlatinib that is mainly controlled with dose interruptions/modification and lipid-lowering therapies. We report here that clinically relevant concentrations of lorlatinib directly modify the lipid profile of human hepatocytes. Silibinin, a flavonolignan with remarkable clinical efficacy against brain metastases in patients with ALK-rearranged NSCLC and capable of overcoming resistance to ALK-TKIs in vitro [76,77], protects the native lipidome of hepatic cells against the hyperlipidemic effects of lorlatinib and prevents lipid accumulation at therapeutically relevant concentrations. Although silibinin might become a new candidate to clinically manage the undesirable lorlatinib-induced hypercholesterolemic and hyperlipidemic effects in human hepatic cells, further studies will be required to fully define the inhibitory potential of silibinin against the lorlatinib-metabolizing CYP3A4 isoenzyme when used as specific orally bioavailable formulations (e.g., Eurosil 85®) before suggesting its therapeutic combination with lorlatinib in a clinical setting with ALK-rearranged lung cancer patients.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms23179986/s1.

Author Contributions

Conceptualization, S.V., J.A.E., J.B-B., and J.A.M.; methodology, S.V., J.A.E., S.F.-A., J.J., E.C., and J.A.M.; formal analysis, S.V., J.A.E., S.F.-A., J.J., E.C., and J.A.M.; investigation, S.V., J.A.E., S.F.-A., J.J., E.C., and J.A.M.; validation, S.V., E.C., and J.A.M.; data curation, S.V., J.A.E., S.F.-A., J.J., E.C., J.B.-B., and J.A.M.; writing—original draft preparation, E.C., J.B.-B., and J.A.M..; writing—review and editing, J.A.M.; visualization, E.C. and J.A.M.; supervision, S.V., E.C., and J.A.M.; project administration, E.C.; funding acquisition, E.C., J.B.-B., and J.A.M. All authors have read and agreed to the published version of the manuscript.

Funding

Work in the Menendez laboratory is supported by the Spanish Ministry of Science and Innovation (Grant PID2019-10455GB-I00, Plan Nacional de l + D + I, founded by the European Regional Development Fund, Spain) and by an unrestricted research grant from the Fundació Oncolliga Girona (Lliga catalana d’ajuda al malalt de càncer, Girona). Work in the Jose A. Encinar laboratory is supported by the Spanish Ministry of Economy and Competitiveness (MINECO, Grant RTI2019-096724-B-C21) and the Generalitat Valenciana (PROMETEO/2021/059). Joaquim Bosch-Barrera is the recipient of Research Grants from Grupo Español de Cáncer de Pulmón (GECP), La Marató de TV3 foundation (201906), and the Health Research and Innovation Strategic Plan (SLT006/17/114; PERIS 2016–2020); Pla stratègic de recerca i innovació en salut; Departament de Salut, Generalitat de Catalunya. Elisabet Cuyàs holds a “Miguel Servet” research contract (CP20/00003) from the Instituto de Salud Carlos III (Spain) and is supported by the Spanish Ministry of Science and Innovation (Grant PI22/00297, Proyectos de I + D + I en Salud, Acción Estratégica en Salud 2021–2023, founded by the European Regional Development Fund, Spain).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All data generated or analyzed during this study are included in this published article (and its supplementary information file).

Acknowledgments

The authors would like to heartily thank König (www.konig.cat) for their generous donation to fund our cancer research on silibinin in Girona. We are grateful to the Servicio de Supercomputación of the Granada University for letting us take advantage of the ALBAICIN computer cluster (https://supercomputacion.ugr.es/arquitecturas/albaicin/), which has provided ~80% of the calculation time necessary to prepare this manuscript. We wish to especially thank Santiago Melchor for his kind assistance and cooperation. We are also grateful to the Cluster of Scientific Computing (http://ccc.umh.es) of the Miguel Hernández University (UMH) for providing computer facilities. The authors would like to thank Kenneth McCreath for detailed editing of this manuscript.

Conflicts of Interest

J.B.-B. reports grants and personal fees from Roche-Genentech, grants from Pfizer and Pierre Fabre, and personal fees from MSD, BMS, AstraZeneca, Boehringer-Ingelheim, and Novartis, outside the submitted work. This study was partially supported by a research grant from Pfizer Global Research and Development to J.B.-B and J.A.M. The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

  1. Shaw, A.T.; Felip, E.; Bauer, T.M.; Besse, B.; Navarro, A.; Postel-Vinay, S.; Gainor, J.F.; Johnson, M.; Dietrich, J.; James, L.P.; et al. Lorlatinib in non-small-cell lung cancer with ALK or ROS1 rearrangement: An international, multicentre, open-label, single-arm first-in-man phase 1 trial. Lancet Oncol. 2017, 18, 1590–1599. [Google Scholar] [CrossRef]
  2. Solomon, B.J.; Besse, B.; Bauer, T.M.; Felip, E.; Soo, R.A.; Camidge, D.R.; Chiari, R.; Bearz, A.; Lin, C.C.; Gadgeel, S.M.; et al. Lorlatinib in patients with ALK-positive non-small-cell lung cancer: Results from a global phase 2 study. Lancet Oncol. 2018, 19, 1654–1667. [Google Scholar] [CrossRef]
  3. Shaw, A.T.; Bauer, T.M.; de Marinis, F.; Felip, E.; Goto, Y.; Liu, G.; Mazieres, J.; Kim, D.W.; Mok, T.; Polli, A.; et al. First-Line Lorlatinib or Crizotinib in Advanced ALK-Positive Lung Cancer. N. Engl. J. Med. 2020, 383, 2018–2029. [Google Scholar] [CrossRef] [PubMed]
  4. Peng, L.; Zhu, L.; Sun, Y.; Stebbing, J.; Selvaggi, G.; Zhang, Y.; Yu, Z. Targeting ALK Rearrangements in NSCLC: Current State of the Art. Front. Oncol. 2022, 12, 863461. [Google Scholar] [CrossRef] [PubMed]
  5. Yang, J.; Gong, W. Lorlatinib for the treatment of anaplastic lymphoma kinase-positive non-small cell lung cancer. Expert Rev. Clin. Pharmacol. 2019, 12, 173–178. [Google Scholar] [CrossRef]
  6. Choo, J.R.; Soo, R.A. Lorlatinib for the treatment of ALK-positive metastatic non-small cell lung cancer. Expert Rev. Anticancer Ther. 2020, 20, 233–240. [Google Scholar] [CrossRef]
  7. Johnson, T.W.; Richardson, P.F.; Bailey, S.; Brooun, A.; Burke, B.J.; Collins, M.R.; Cui, J.J.; Deal, J.G.; Deng, Y.L.; Dinh, D.; et al. Discovery of (10R)-7-amino-12-fluoro-2,10,16-trimethyl-15-oxo-10,15,16,17-tetrahydro-2H-8,4-(metheno)pyrazolo [4,3-h][2,5,11]-benzoxadiazacyclotetradecine-3-carbonitrile (PF-06463922), a macrocyclic inhibitor of anaplastic lymphoma kinase (ALK) and c-ros oncogene 1 (ROS1) with preclinical brain exposure and broad-spectrum potency against ALK-resistant mutations. J. Med. Chem. 2014, 57, 4720–4744. [Google Scholar]
  8. Collier, T.L.; Normandin, M.D.; Stephenson, N.A.; Livni, E.; Liang, S.H.; Wooten, D.W.; Esfahani, S.A.; Stabin, M.G.; Mahmood, U.; Chen, J.; et al. Synthesis and preliminary PET imaging of 11C and 18F isotopologues of the ROS1/ALK inhibitor lorlatinib. Nat. Commun. 2017, 8, 15761. [Google Scholar] [CrossRef]
  9. Costa, D.B.; Shaw, A.T.; Ou, S.H.; Solomon, B.J.; Riely, G.J.; Ahn, M.J.; Zhou, C.; Shreeve, S.M.; Selaru, P.; Polli, A.; et al. Clinical Experience with Crizotinib in Patients with Advanced ALK-Rearranged Non-Small-Cell Lung Cancer and Brain Metastases. J. Clin. Oncol. 2015, 33, 1881–1888. [Google Scholar] [CrossRef]
  10. Bauer, T.M.; Shaw, A.T.; Johnson, M.L.; Navarro, A.; Gainor, J.F.; Thurm, H.; Pithavala, Y.K.; Abbattista, A.; Peltz, G.; Felip, E. Brain Penetration of Lorlatinib: Cumulative Incidences of CNS and Non-CNS Progression with Lorlatinib in Patients with Previously Treated ALK-Positive Non-Small-Cell Lung Cancer. Target Oncol. 2020, 15, 55–65. [Google Scholar] [CrossRef]
  11. Felip, E.; Shaw, A.T.; Bearz, A.; Camidge, D.R.; Solomon, B.J.; Bauman, J.R.; Bauer, T.M.; Peters, S.; Toffalorio, F.; Abbattista, A.; et al. Intracranial and extracranial efficacy of lorlatinib in patients with ALK-positive non-small-cell lung cancer previously treated with second-generation ALK TKIs. Ann. Oncol. 2021, 32, 620–630. [Google Scholar] [CrossRef] [PubMed]
  12. Yun, K.M.; Bazhenova, L.A. Update on Lorlatinib: Role in Reducing the Risk of Disease Progression in ALK-Positive NSCLC. Cancer Manag. Res. 2022, 14, 843–850. [Google Scholar] [CrossRef] [PubMed]
  13. Syed, Y.Y. Lorlatinib: First Global Approval. Drugs 2019, 79, 93–98. [Google Scholar] [CrossRef]
  14. Bauer, T.M.; Felip, E.; Solomon, B.J.; Thurm, H.; Peltz, G.; Chioda, M.D.; Shaw, A.T. Clinical Management of Adverse Events Associated with Lorlatinib. Oncologist 2019, 24, 1103–1110. [Google Scholar] [CrossRef] [PubMed]
  15. Blais, N.; Adam, J.-P.; Nguyen, J.; Grégoire, J.C. Evaluation and Management of Dyslipidemia in Patients Treated with Lorlatinib. Curr. Oncol. 2021, 28, 265–272. [Google Scholar] [CrossRef]
  16. Neuvonen, P.J.; Niemi, M.; Backman, J.T. Drug interactions with lipid-lowering drugs: Mechanisms and clinical relevance. Clin. Pharmacol. Ther. 2006, 80, 565–581. [Google Scholar] [CrossRef] [PubMed]
  17. Reed, M.; Rosales, A.S.; Chioda, M.D.; Parker, L.; Devgan, G.; Kettle, J. Consensus Recommendations for Management and Counseling of Adverse Events Associated With Lorlatinib: A Guide for Healthcare Practitioners. Adv. Ther. 2020, 37, 3019–3030. [Google Scholar] [CrossRef]
  18. Fernández-Arroyo, S.; Hernández-Aguilera, A.; de Vries, M.A.; Burggraaf, B.; van der Zwan, E.; Pouw, N.; Joven, J.; Cabezas, M.C. Effect of Vitamin D3 on the Postprandial Lipid Profile in Obese Patients: A Non-Targeted Lipidomics Study. Nutrients 2019, 11, 1194. [Google Scholar] [CrossRef]
  19. Nakabayashi, H.; Taketa, K.; Miyano, K.; Yamane, T.; Sato, J. Growth of human hepatoma cells lines with differentiated functions in chemically defined medium. Cancer Res. 1982, 42, 3858–3863. [Google Scholar]
  20. Kawamoto, M.; Yamaji, T.; Saito, K.; Shirasago, Y.; Satomura, K.; Endo, T.; Fukasawa, M.; Hanada, K.; Osada, N. Identification of Characteristic Genomic Markers in Human Hepatoma HuH-7 and Huh7.5.1-8 Cell Lines. Front. Genet. 2020, 11, 546106. [Google Scholar] [CrossRef]
  21. Ramboer, E.; Vanhaecke, T.; Rogiers, V.; Vinken, M. Immortalized Human Hepatic Cell Lines for In Vitro Testing and Research Purposes. Methods Mol. Biol. 2015, 1250, 53–76. [Google Scholar] [PubMed]
  22. Berger, E.; Vega, N.; Weiss-Gayet, M.; Géloën, A. Gene Network Analysis of Glucose Linked Signaling Pathways and Their Role in Human Hepatocellular Carcinoma Cell Growth and Survival in HuH7 and HepG2 Cell Lines. Biomed. Res. Int. 2015, 2015, 821761. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Bulutoglu, B.; Rey-Bedón, C.; Mert, S.; Tian, L.; Jang, Y.Y.; Yarmush, M.L.; Usta, O.B. A comparison of hepato-cellular in vitro platforms to study CYP3A4 induction. PLoS ONE 2020, 15, e0229106. [Google Scholar] [CrossRef] [PubMed]
  24. Arzumanian, V.A.; Kiseleva, O.I.; Poverennaya, E.V. The Curious Case of the HepG2 Cell Line: 40 Years of Expertise. Int. J. Mol. Sci. 2021, 22, 13135. [Google Scholar] [CrossRef] [PubMed]
  25. Bijak, M. Silybin, a Major Bioactive Component of Milk Thistle (Silybum marianum L. Gaernt.)-Chemistry, Bioavailability, and Metabolism. Molecules 2017, 22, 1942. [Google Scholar] [CrossRef]
  26. Křen, V.; Valentová, K. Silybin and its congeners: From traditional medicine to molecular effects. Nat. Prod. Rep. 2022, 39, 1264–1281. [Google Scholar] [CrossRef]
  27. Wang, X.; Zhang, Z.; Wu, S.C. Health Benefits of Silybum marianum: Phytochemistry, Pharmacology, and Applications. J. Agric. Food Chem. 2020, 68, 11644–11664. [Google Scholar] [CrossRef]
  28. Federico, A.; Dallio, M.; Loguercio, C. Silymarin/Silybin and Chronic Liver Disease: A Marriage of Many Years. Molecules 2017, 22, 191. [Google Scholar] [CrossRef]
  29. Hackett, E.S.; Twedt, D.C.; Gustafson, D.L. Milk thistle and its derivative compounds: A review of opportunities for treatment of liver disease. J. Vet. Intern. Med. 2013, 27, 10–16. [Google Scholar] [CrossRef]
  30. Abenavoli, L.; Capasso, R.; Milic, N.; Capasso, F. Milk thistle in liver diseases: Past, present, future. Phytother. Res. 2010, 24, 1423–1432. [Google Scholar] [CrossRef]
  31. Rikova, K.; Guo, A.; Zeng, Q.; Possemato, A.; Yu, J.; Haack, H.; Nardone, J.; Lee, K.; Reeves, C.; Li, Y.; et al. Global survey of phosphotyrosine signaling identifies oncogenic kinases in lung cancer. Cell 2007, 131, 1190–1203. [Google Scholar] [CrossRef] [PubMed]
  32. Li, Y.; Ye, X.; Liu, J.; Zha, J.; Pei, L. Evaluation of EML4-ALK fusion proteins in non-small cell lung cancer using small molecule inhibitors. Neoplasia 2011, 13, 1–11. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Samuels, E.R.; Sevrioukova, I.F. Rational Design of CYP3A4 Inhibitors: A One-Atom Linker Elongation in Ritonavir-Like Compounds Leads to a Marked Improvement in the Binding Strength. Int. J. Mol. Sci. 2021, 22, 852. [Google Scholar] [CrossRef] [PubMed]
  34. Domanski, T.L.; He, Y.A.; Harlow, G.R.; Halpert, J.R. Dual role of human cytochrome P450 3A4 residue Phe-304 in substrate specificity and cooperativity. J. Pharmacol. Exp. Ther. 2000, 293, 585–591. [Google Scholar]
  35. Hackett, J.C. Membrane-embedded substrate recognition by cytochrome P450 3A4. J. Biol. Chem. 2018, 293, 4037–4046. [Google Scholar] [CrossRef]
  36. Sevrioukova, I.F.; Poulos, T.L. Understanding the mechanism of cytochrome P450 3A4: Recent advances and remaining problems. Dalton Trans. 2013, 42, 3116–3126. [Google Scholar] [CrossRef]
  37. Sevrioukova, I.F.; Poulos, T.L. Structure and mechanism of the complex between cytochrome P4503A4 and ritonavir. Proc. Natl. Acad. Sci. USA 2010, 107, 18422–18427. [Google Scholar] [CrossRef]
  38. Sevrioukova, I.F.; Poulos, T.L. Ritonavir analogues as a probe for deciphering the cytochrome P450 3A4 inhibitory mechanism. Curr. Top. Med. Chem. 2014, 14, 1348–1355. [Google Scholar] [CrossRef]
  39. Li, X.; Xu, Y.; Lai, L.; Pei, J. Prediction of Human Cytochrome P450 Inhibition Using a Multitask Deep Autoencoder Neural Network. Mol. Pharm. 2018, 15, 4336–4345. [Google Scholar] [CrossRef]
  40. Gillessen, A.; Schmidt, H.H. Silymarin as Supportive Treatment in Liver Diseases: A Narrative Review. Adv. Ther. 2020, 37, 1279–1301. [Google Scholar] [CrossRef]
  41. Pérez-Sánchez, A.; Cuyàs, E.; Ruiz-Torres, V.; Agulló-Chazarra, L.; Verdura, S.; González-Álvarez, I.; Bermejo, M.; Joven, J.; Micol, V.; Bosch-Barrera, J.; et al. Intestinal Permeability Study of Clinically Relevant Formulations of Silibinin in Caco-2 Cell Monolayers. Int. J. Mol. Sci. 2019, 20, 1606. [Google Scholar] [CrossRef] [PubMed]
  42. Schnell, P.; Bartlett, C.H.; Solomon, B.J.; Tassell, V.; Shaw, A.T.; de Pas, T.; Lee, S.H.; Lee, G.K.; Tanaka, K.; Tan, W.; et al. Complex renal cysts associated with crizotinib treatment. Cancer Med. 2015, 4, 887–896. [Google Scholar] [CrossRef] [PubMed]
  43. Camidge, D.R.; Kim, E.E.; Usari, T.; Polli, A.; Lewis, I.; Wilner, K.D. Renal effects of crizotinib in patients with ALK-positive advanced NSCLC. J. Thorac. Oncol. 2019, 14, 1077–1085. [Google Scholar] [CrossRef] [PubMed]
  44. Hong, M.H. Nephrotoxicity of cancer therapeutic drugs: Focusing on novel agents. Kidney Res. Clin. Pract. 2021, 40, 344–354. [Google Scholar] [CrossRef] [PubMed]
  45. McGee, K.; Stone, N.J.; Wadhwani, S.; Kanwar, Y.S.; Villaflor, V.; Akhter, N. A possible mechanism of hyperlipidemia in a patient with metastatic non-small cell lung cancer on lorlatinib therapy. J. Oncol. Pharm. Pract. 2021, 27, 2010–2013. [Google Scholar] [CrossRef]
  46. Hibma, J.E.; O’Gorman, M.; Nepal, S.; Pawlak, S.; Ginman, K.; Pithavala, Y.K. Evaluation of the absolute oral bioavailability of the anaplastic lymphoma kinase/c-ROS oncogene 1 kinase inhibitor lorlatinib in healthy participants. Cancer Chemother. Pharmacol. 2022, 89, 71–81. [Google Scholar] [CrossRef] [PubMed]
  47. Chen, J.; Houk, B.; Pithavala, Y.K.; Ruiz-Garcia, A. Population pharmacokinetic model with time-varying clearance for lorlatinib using pooled data from patients with non-small cell lung cancer and healthy participants. CPT Pharmacomet. Syst. Pharmacol. 2021, 10, 148–160. [Google Scholar] [CrossRef]
  48. Chen, W.; Shi, Y.; Qi, S.; Zhou, H.; Li, C.; Jin, D.; Li, G. Pharmacokinetic Study and Tissue Distribution of Lorlatinib in Mouse Serum and Tissue Samples by Liquid Chromatography-Mass Spectrometry. J. Anal. Methods Chem. 2019, 2019, 7574369. [Google Scholar] [CrossRef]
  49. Rodríguez-Hernández, M.A.; de la Cruz-Ojeda, P.; López-Grueso, M.J.; Navarro-Villarán, E.; Requejo-Aguilar, R.; Castejón-Vega, B.; Negrete, M.; Gallego, P.; Vega-Ochoa, Á.; Victor, V.M.; et al. Integrated molecular signaling involving mitochondrial dysfunction and alteration of cell metabolism induced by tyrosine kinase inhibitors in cancer. Redox Biol. 2020, 36, 101510. [Google Scholar] [CrossRef]
  50. Mihajlovic, M.; Vinken, M. Mitochondria as the Target of Hepatotoxicity and Drug-Induced Liver Injury: Molecular Mechanisms and Detection Methods. Int. J. Mol. Sci. 2022, 23, 3315. [Google Scholar] [CrossRef]
  51. Yin, F.; Gupta, R.; Vergnes, L.; Driscoll, W.S.; Ricks, J.; Ramanathan, G.; Stewart, J.A.; Shih, D.M.; Faull, K.F.; Beaven, S.W.; et al. Diesel Exhaust Induces Mitochondrial Dysfunction, Hyperlipidemia, and Liver Steatosis. Arterioscler. Thromb. Vasc. Biol. 2019, 39, 1776–1786. [Google Scholar] [CrossRef]
  52. Remon, J.; Besse, B. Brain Metastases in Oncogene-Addicted Non-Small Cell Lung Cancer Patients: Incidence and Treatment. Front Oncol. 2018, 8, 88. [Google Scholar] [CrossRef] [PubMed]
  53. Bosch-Barrera, J.; Sais, E.; Cañete, N.; Marruecos, J.; Cuyàs, E.; Izquierdo, A.; Porta, R.; Haro, M.; Brunet, J.; Pedraza, S.; et al. Response of brain metastasis from lung cancer patients to an oral nutraceutical product containing silibinin. Oncotarget 2016, 7, 32006–32014. [Google Scholar] [CrossRef] [Green Version]
  54. Priego, N.; Zhu, L.; Monteiro, C.; Mulders, M.; Wasilewski, D.; Bindeman, W.; Doglio, L.; Martínez, L.; Martínez-Saez, E.; Ramón, Y.; et al. STAT3 labels a subpopulation of reactive astrocytes required for brain metastasis. Nat. Med. 2018, 24, 1024–1035. [Google Scholar] [CrossRef] [PubMed]
  55. Li, W.; Sparidans, R.W.; Wang, Y.; Lebre, M.C.; Wagenaar, E.; Beijnen, J.H.; Schinkel, A.H. P-glycoprotein (MDR1/ABCB1) restricts brain accumulation and cytochrome P450-3A (CYP3A) limits oral availability of the novel ALK/ROS1 inhibitor lorlatinib. Int. J. Cancer 2018, 143, 2029–2038. [Google Scholar] [CrossRef] [PubMed]
  56. Hu, W.; Lettiere, D.; Tse, S.; Johnson, T.R.; Biddle, K.E.; Thibault, S.; Palazzi, X.; Chen, J.; Pithavala, Y.K.; Finkelstein, M. Liver Toxicity Observed With Lorlatinib When Combined With Strong CYP3A Inducers: Evaluation of Cynomolgus Monkey as a Nonclinical Model for Assessing the Mechanism of Combinational Toxicity. Toxicol. Sci. 2021, 182, 183–194. [Google Scholar] [CrossRef]
  57. Brantley, S.J.; Graf, T.N.; Oberlies, N.H.; Paine, M.F. A systematic approach to evaluate herb-drug interaction mechanisms: Investigation of milk thistle extracts and eight isolated constituents as CYP3A inhibitors. Drug Metab. Dispos. 2013, 41, 1662–1670. [Google Scholar] [CrossRef] [PubMed]
  58. Beckmann-Knopp, S.; Rietbrock, S.; Weyhenmeyer, R.; Böcker, R.H.; Beckurts, K.T.; Lang, W.; Hunz, M.; Fuhr, U. Inhibitory effects of silibinin on cytochrome P-450 enzymes in human liver microsomes. Pharmacol. Toxicol. 2000, 86, 250–256. [Google Scholar] [CrossRef]
  59. Lin, J.; Schyschka, L.; Mühl-Benninghaus, R.; Neumann, J.; Hao, L.; Nussler, N.; Dooley, S.; Liu, L.; Stöckle, U.; Nussler, A.K.; et al. Comparative analysis of phase I and II enzyme activities in 5 hepatic cell lines identifies Huh-7 and HCC-T cells with the highest potential to study drug metabolism. Arch. Toxicol. 2012, 86, 87–95. [Google Scholar] [CrossRef]
  60. Liu, Y.; Flynn, T.J.; Xia, M.; Wiesenfeld, P.L.; Ferguson, M.S. Evaluation of CYP3A4 inhibition and hepatotoxicity using DMSO-treated human hepatoma HuH-7 cells. Cell. Biol. Toxicol. 2015, 31, 221–230. [Google Scholar] [CrossRef]
  61. Guo, L.; Dial, S.; Branham, W.; Liu, J.; Fang, J.L.; Green, B.; Deng, H.; Kaput, J.; Ning, B. Similarities and differences in the expression of drug-metabolizing enzymes between human hepatic cell lines and primary human hepatocytes. Drug Metab. Dispos. 2011, 39, 528–538. [Google Scholar] [CrossRef] [PubMed]
  62. Sivertsson, L.; Ek, M.; Darnell, M.; Edebert, I.; Ingelman-Sundberg, M.; Neve, E.P. CYP3A4 catalytic activity is induced in confluent Huh7 hepatoma cells. Drug Metab. Dispos. 2010, 38, 995–1002. [Google Scholar] [CrossRef] [PubMed]
  63. Sivertsson, L.; Edebert, I.; Palmertz, M.P.; Ingelman-Sundberg, M.; Neve, E.P. Induced CYP3A4 expression in confluent Huh7 hepatoma cells as a result of decreased cell proliferation and subsequent pregnane X receptor activation. Mol. Pharmacol. 2013, 83, 659–670. [Google Scholar] [CrossRef]
  64. Kanebratt, K.P.; Andersson, T.B. HepaRG cells as an in vitro model for evaluation of cytochrome P450 induction in humans. Drug Metab. Dispos. 2008, 36, 137–145. [Google Scholar] [CrossRef] [PubMed]
  65. Kanebratt, K.P.; Andersson, T.B. Evaluation of HepaRG cells as an in vitro model for human drug metabolism studies. Drug Metab. Dispos. 2008, 36, 1444–1452. [Google Scholar] [CrossRef] [PubMed]
  66. Andersson, T.B.; Kanebratt, K.P.; Kenna, J.G. The HepaRG cell line: A unique in vitro tool for understanding drug metabolism and toxicology in human. Expert Opin. Drug Metab. Toxicol. 2012, 8, 909–920. [Google Scholar] [CrossRef] [PubMed]
  67. Corominas-Faja, B.; Cuyàs, E.; Lozano-Sánchez, J.; Cufí, S.; Verdura, S.; Fernández-Arroyo, S.; Borrás-Linares, I.; Martin-Castillo, B.; Martin, Á.G.; Lupu, R.; et al. Extra-virgin olive oil contains a metabolo-epigenetic inhibitor of cancer stem cells. Carcinogenesis 2018, 39, 601–613. [Google Scholar] [CrossRef]
  68. Verdura, S.; Cuyàs, E.; Lozano-Sánchez, J.; Bastidas-Velez, C.; Llorach-Parés, L.; Fernández-Arroyo, S.; Hernández-Aguilera, A.; Joven, J.; Nonell-Canals, A.; Bosch-Barrera, J.; et al. An olive oil phenolic is a new chemotype of mutant isocitrate dehydrogenase 1 (IDH1) inhibitors. Carcinogenesis 2019, 40, 27–40. [Google Scholar] [CrossRef]
  69. Cuyàs, E.; Castillo, D.; Llorach-Parés, L.; Lozano-Sánchez, J.; Verdura, S.; Nonell-Canals, A.; Brunet, J.; Bosch-Barrera, J.; Joven, J.; Valdés, R.; et al. Computational de-orphanization of the olive oil biophenol oleacein: Discovery of new metabolic and epigenetic targets. Food Chem. Toxicol. 2019, 131, 110529. [Google Scholar] [CrossRef]
  70. Encinar, J.A.; Menendez, J.A. Potential Drugs Targeting Early Innate Immune Evasion of SARS-Coronavirus 2 via 2′-O-Methylation of Viral RNA. Viruses 2020, 12, 525. [Google Scholar] [CrossRef]
  71. Tramonti, A.; Cuyàs, E.; Encinar, J.A.; Pietzke, M.; Paone, A.; Verdura, S.; Arbusà, A.; Martin-Castillo, B.; Giardina, G.; Joven, J.; et al. Metformin Is a Pyridoxal-5′-phosphate (PLP)-Competitive Inhibitor of SHMT2. Cancers 2021, 13, 4009. [Google Scholar] [CrossRef] [PubMed]
  72. Verdura, S.; Cuyàs, E.; Cortada, E.; Brunet, J.; Lopez-Bonet, E.; Martin-Castillo, B.; Bosch-Barrera, J.; Encinar, J.A.; Menendez, J.A. Resveratrol targets PD-L1 glycosylation and dimerization to enhance antitumor T-cell immunity. Aging 2020, 12, 8–34. [Google Scholar] [CrossRef] [PubMed]
  73. Cuyàs, E.; Verdura, S.; Lozano-Sánchez, J.; Viciano, I.; Llorach-Parés, L.; Nonell-Canals, A.; Bosch-Barrera, J.; Brunet, J.; Segura-Carretero, A.; Sanchez-Martinez, M.; et al. The extra virgin olive oil phenolic oleacein is a dual substrate-inhibitor of catechol-O-methyltransferase. Food Chem. Toxicol. 2019, 128, 35–45. [Google Scholar] [CrossRef]
  74. Seeliger, D.; de Groot, B.L. Ligand docking and binding site analysis with PyMOL and Autodock/Vina. J. Comput. Aided Mol. Des. 2010, 24, 417–422. [Google Scholar] [CrossRef]
  75. Salentin, S.; Schreiber, S.; Haupt, V.J.; Adasme, M.F.; Schroeder, M. PLIP: Fully automated protein-ligand interaction profiler. Nucleic Acids Res. 2015, 43, W443–W447. [Google Scholar] [CrossRef] [PubMed]
  76. Cuyàs, E.; Pérez-Sánchez, A.; Micol, V.; Menendez, J.A.; Bosch-Barrera, J. STAT3-targeted treatment with silibinin overcomes the acquired resistance to crizotinib in ALK-rearranged lung cancer. Cell Cycle 2016, 15, 3413–3418. [Google Scholar] [CrossRef] [PubMed]
  77. Verdura, S.; Cuyàs, E.; Ruiz-Torres, V.; Micol, V.; Joven, J.; Bosch-Barrera, J.; Menendez, J.A. Lung Cancer Management with Silibinin: A Historical and Translational Perspective. Pharmaceuticals 2021, 14, 559. [Google Scholar] [CrossRef]
Ijms 23 09986 g001 550
Figure 1. Silibinin prevents lorlatinib-driven lipidome alterations in human hepatic cells: (A) bar graphs of the IC50 values for each cell line calculated from the MTT assays as described in “Materials and Methods”; the results are presented as the means (columns) ± S.D. (bars) (n ≥ 3 in technical replicates); (B) to test the hypothesis that lorlatinib could directly promote the accumulation of triglycerides and/or cholesterol in human hepatocytes, we carried out an untargeted UHPLC-ESI-QTOF-MS/MS-based lipidomic analysis of > 100 molecular lipid species in Huh-7 cells cultured in the absence/presence of lorlatinib and/or silibinin; (C) Volcano plots and Venn diagrams of the results from lipidomic analyses (B) in human Huh-7 cells treated with lorlatinib (1 μmol/L), silibinin (100 μmol/L) or their combination for 48 h. Each dot represents a lipid species with its corresponding mean Log2 fold-change (FC) (x axis) and Benjamini–Hochberg corrected p-value (-log10, y axis). Colored dots illustrate differential lipid species, using a cutoff of p < 0.05 and log2FC > 1 or <1. (D). ELISA-based quantification of cholesteryl esters (left) and triglycerides (right) in Huh-7 cells treated with lorlatinib (1 μmol/L), silibinin (100 μmol/L) or their combination for 48 h. p-values < 0.01 (*) and <0.001 (**).
Figure 1. Silibinin prevents lorlatinib-driven lipidome alterations in human hepatic cells: (A) bar graphs of the IC50 values for each cell line calculated from the MTT assays as described in “Materials and Methods”; the results are presented as the means (columns) ± S.D. (bars) (n ≥ 3 in technical replicates); (B) to test the hypothesis that lorlatinib could directly promote the accumulation of triglycerides and/or cholesterol in human hepatocytes, we carried out an untargeted UHPLC-ESI-QTOF-MS/MS-based lipidomic analysis of > 100 molecular lipid species in Huh-7 cells cultured in the absence/presence of lorlatinib and/or silibinin; (C) Volcano plots and Venn diagrams of the results from lipidomic analyses (B) in human Huh-7 cells treated with lorlatinib (1 μmol/L), silibinin (100 μmol/L) or their combination for 48 h. Each dot represents a lipid species with its corresponding mean Log2 fold-change (FC) (x axis) and Benjamini–Hochberg corrected p-value (-log10, y axis). Colored dots illustrate differential lipid species, using a cutoff of p < 0.05 and log2FC > 1 or <1. (D). ELISA-based quantification of cholesteryl esters (left) and triglycerides (right) in Huh-7 cells treated with lorlatinib (1 μmol/L), silibinin (100 μmol/L) or their combination for 48 h. p-values < 0.01 (*) and <0.001 (**).
Ijms 23 09986 g001
Ijms 23 09986 g002 550
Figure 2. Silibinin inhibits the chronic accumulation of neutral lipids in lorlatinib-treated hepatic cells (A) Conventional indices of fatty acid liver disease ––lipid droplet content in hepatic cells––were assessed using Oil Red O- and LipidTOX green-based staining of neutral lipids in Huh-7 (A) and HepG2 (B) cells chronically (A) or acutely (B) exposed to lorlatinib in the absence or presence of silibinin. (Chronic exposure: 7 days at 10 μmol/L silibinin; acute exposure: 48 h at 100 μmol/L silibinin). Scale bar = 100 μm.
Figure 2. Silibinin inhibits the chronic accumulation of neutral lipids in lorlatinib-treated hepatic cells (A) Conventional indices of fatty acid liver disease ––lipid droplet content in hepatic cells––were assessed using Oil Red O- and LipidTOX green-based staining of neutral lipids in Huh-7 (A) and HepG2 (B) cells chronically (A) or acutely (B) exposed to lorlatinib in the absence or presence of silibinin. (Chronic exposure: 7 days at 10 μmol/L silibinin; acute exposure: 48 h at 100 μmol/L silibinin). Scale bar = 100 μm.
Ijms 23 09986 g002
Ijms 23 09986 g003 550
Figure 3. Silibinin does not share the binding mode of lorlatinib to CYP3A4. Figure depicts the backbone of the overall crystal structure of CYP3A4 (7KVS) with rainbow colors showing the best docked poses of silibinin A, silibinin B, ritonavir, and lorlatinib at the catalytic site. Each inset shows the detailed interactions of silibinin docked to CYP3A4, indicating the participating amino acids involved in the interaction and the type of interaction (hydrogen bonds, hydrophilic interactions, salt bridges, II-stacking, etc).
Figure 3. Silibinin does not share the binding mode of lorlatinib to CYP3A4. Figure depicts the backbone of the overall crystal structure of CYP3A4 (7KVS) with rainbow colors showing the best docked poses of silibinin A, silibinin B, ritonavir, and lorlatinib at the catalytic site. Each inset shows the detailed interactions of silibinin docked to CYP3A4, indicating the participating amino acids involved in the interaction and the type of interaction (hydrogen bonds, hydrophilic interactions, salt bridges, II-stacking, etc).
Ijms 23 09986 g003
Ijms 23 09986 g004 550
Figure 4. Incorporation models of CYP3A4-targeting drugs: ritonavir and statins vs silibinin. The root mean square deviation (RMSD, Å) of each drug’s heavy atoms over the simulation time, measured after superposing the protein onto its reference structure, and the molecular mechanics Poisson–Boltzmann surface area (MM/PBSA) binding energy analyses calculated from the entire trajectory of the 100 ns (or last 30 ns) MD simulation, are shown. The best poses of ritonavir, silibinin A, silibinin B, simvastatin, and rosuvastatin coupled to the catalytic site of CYP3A4 before (0 ns) and after (100 ns) the molecular dynamics (MD) simulation are shown. The protein is represented as a function of the hydrophobicity of its surface amino acids and the Na+ and Cl− ions have been eliminated to facilitate visualization. Each inset shows the detailed interactions of the participating amino acids involved and the type of interaction (hydrogen bonds, hydrophilic interactions, salt bridges, Π-stacking, etc.).
Figure 4. Incorporation models of CYP3A4-targeting drugs: ritonavir and statins vs silibinin. The root mean square deviation (RMSD, Å) of each drug’s heavy atoms over the simulation time, measured after superposing the protein onto its reference structure, and the molecular mechanics Poisson–Boltzmann surface area (MM/PBSA) binding energy analyses calculated from the entire trajectory of the 100 ns (or last 30 ns) MD simulation, are shown. The best poses of ritonavir, silibinin A, silibinin B, simvastatin, and rosuvastatin coupled to the catalytic site of CYP3A4 before (0 ns) and after (100 ns) the molecular dynamics (MD) simulation are shown. The protein is represented as a function of the hydrophobicity of its surface amino acids and the Na+ and Cl− ions have been eliminated to facilitate visualization. Each inset shows the detailed interactions of the participating amino acids involved and the type of interaction (hydrogen bonds, hydrophilic interactions, salt bridges, Π-stacking, etc.).
Ijms 23 09986 g004
Ijms 23 09986 g005 550
Figure 5. Potency of human CYP3A4 inhibition: statins vs silibinin. Percent inhibition and IC50 values of simvastatin, rosuvastatin, silibinin, and Eurosil 85® against recombinant human CYP3A4 expressed in baculosomes using the Vivid BOMR Red substrate. The results are presented as the means (columns) ± S.D. (bars) (n = 3 experimental curves in technical duplicates).
Figure 5. Potency of human CYP3A4 inhibition: statins vs silibinin. Percent inhibition and IC50 values of simvastatin, rosuvastatin, silibinin, and Eurosil 85® against recombinant human CYP3A4 expressed in baculosomes using the Vivid BOMR Red substrate. The results are presented as the means (columns) ± S.D. (bars) (n = 3 experimental curves in technical duplicates).
Ijms 23 09986 g005
Ijms 23 09986 g006 550
Figure 6. Silibinin suppresses the hyperlipidemic effects of lorlatinib in hepatic cells. The third generation ALK-TKI lorlatinib promotes the accumulation of cholesterol and triglycerides in human hepatic cells. The flavonolignan silibinin protects hepatic cells against the hypertriglyceridemic and hypercholesteremic effects of lorlatinib. Silibinin is a new candidate to clinically manage the undesirable hyperlipidemic activity of lorlatinib in patients with ALK-rearranged lung cancer.
Figure 6. Silibinin suppresses the hyperlipidemic effects of lorlatinib in hepatic cells. The third generation ALK-TKI lorlatinib promotes the accumulation of cholesterol and triglycerides in human hepatic cells. The flavonolignan silibinin protects hepatic cells against the hypertriglyceridemic and hypercholesteremic effects of lorlatinib. Silibinin is a new candidate to clinically manage the undesirable hyperlipidemic activity of lorlatinib in patients with ALK-rearranged lung cancer.
Ijms 23 09986 g006
Table
Table 1. Lipid species (n = 124) identified using a UHPLC-ESI-QTOF-MS/MS method (e:ester).
Table 1. Lipid species (n = 124) identified using a UHPLC-ESI-QTOF-MS/MS method (e:ester).
DG (n = 9)TG (n = 25)PC (n = 40)LPC (n = 9)PE (n = 4)SM (n = 24)ChoE (n = 13)
DG 34:1TG 46:0PC 30:0LPC 15:0PE 32:0SM 32:0ChoE 16:0
DG 34:2TG 46:1PC 31:0LPC 16:0PE 36:4SM 32:1ChoE 16:1
DG 34:3TG 46:2PC 32:0LPC 16:0ePE 36:5eSM 32:2ChoE 17:0
DG 36:0TG 48:0PC 32:1LPC 18:0PE 38:5eSM 33:1ChoE 17:1
DG 36:1TG 48:1PC 32:1eLPC 18:0e SM 34:1ChoE 18:0
DG 36:2TG 48:2PC 32:2LPC 18:1 SM 34:2ChoE 18:2
DG 36:3TG 48:3PC 33:0LPC 18:2 SM 35:0ChoE 18:3
DG 36:4TG 50:0PC 33:1LPC 20:0 SM 35:1ChoE 20:2
DG 40:4TG 50:1PC 33:2LPC 20:3 SM 36:0ChoE 20:3
TG 50:2PC 34:0 SM 36:1ChoE 20:4
TG 50:3PC 34:1 SM 36:2ChoE 20:5
TG 50:4PC 34:1e SM 38:1ChoE 22:4
TG 51:2PC 34:2 SM 38:2ChoE 22:5
TG 52:1PC 34:2e SM 39:1
TG 52:2PC 34:3 SM 40:0
TG 52:3PC 34:4 SM 40:1
TG 52:4PC 35:1 SM 40:2
TG 52:5PC 35:2 SM 41:1
TG 52:6PC 35:4 SM 41:2
TG 54:2PC 36:0 SM 42:1
TG 54:3PC 36:1 SM 42:2
TG 54:4PC 36:2 SM 42:3
TG 54:5PC 36:2e SM 43:1
TG 54:6PC 36:3 SM 43:2
TG 54:7PC 36:4
PC 36:5
PC 36:5e
PC 38:2
PC 38:3
PC 38:4
PC 38:5
PC 38:5e
PC 38:6
PC 38:6e
PC 40:4
PC 40:4e
PC 40:5
PC 40:6
PC 42:4e
PC 42:5e
Table
Table 2. Details of the interaction between silibinin A/B, ritonavir and lorlatinib to CYP3A4.
Table 2. Details of the interaction between silibinin A/B, ritonavir and lorlatinib to CYP3A4.
ΔG (kcal/mol)Kd [nM]DrugResidues Involved in the Interaction (7KVS.pdb)
−10.32327.1silibinin ATYR 53, PHE 57, ASP 76, ARG 105, ARG 106, PHE 108,
GLY 109, PHE 213, PHE 215, PHE 220, ILE 223, THR 224,
PRO 227, ILE 230, VAL 240, ALA 370, MET 371,
ARG 372, LEU 373, GLU 374, ARG 375, HEME 601
−9.96249.9ritonavirPHE 57, ARG 105, ARG 106, PHE 108, MET 114, SER 119,
ILE 120, LEU 210, LEU 211, PHE 213, PHE 241, ILE 300,
ILE 301, PHE 304, ALA 305, THR 309, ILE 369, ALA 370,
MET 371, ARG 372, LEU 373, GLU 374, HEME 601
−9.87657.6lorlatinibPHE 57, ARG 105, SER 119, LEU 211, PHE 304,
GLU 308, THR 309, SER 312, ILE 369, ALA 370,
MET 371, LEU 373, LEU 482, HEME 601
−9.65184.3silibinin BPHE 57, ARG 105, ARG 106, PRO 107, PHE 108, SER 119,
ILE 301, ALA 305, THR 309, ALA 370, MET 371,
ARG 372, LEU 373, GLU 374, HEME 601
Table
Table 3. Prediction of human cytochrome P450 inhibition.
Table 3. Prediction of human cytochrome P450 inhibition.
Cytochrome P450 Isoforms
Drug1A22C92C192D63A4
Ritonavir0.000.340.360.010.97
Simvastatin0.000.020.040.000.93
Rosuvastatin0.000.450.180.000.2
Silibinin0.000.020.040.000.06
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Verdura, S.; Encinar, J.A.; Fernández-Arroyo, S.; Joven, J.; Cuyàs, E.; Bosch-Barrera, J.; Menendez, J.A. Silibinin Suppresses the Hyperlipidemic Effects of the ALK-Tyrosine Kinase Inhibitor Lorlatinib in Hepatic Cells. Int. J. Mol. Sci. 2022, 23, 9986. https://doi.org/10.3390/ijms23179986

AMA Style

Verdura S, Encinar JA, Fernández-Arroyo S, Joven J, Cuyàs E, Bosch-Barrera J, Menendez JA. Silibinin Suppresses the Hyperlipidemic Effects of the ALK-Tyrosine Kinase Inhibitor Lorlatinib in Hepatic Cells. International Journal of Molecular Sciences. 2022; 23(17):9986. https://doi.org/10.3390/ijms23179986

Chicago/Turabian Style

Verdura, Sara, José Antonio Encinar, Salvador Fernández-Arroyo, Jorge Joven, Elisabet Cuyàs, Joaquim Bosch-Barrera, and Javier A. Menendez. 2022. "Silibinin Suppresses the Hyperlipidemic Effects of the ALK-Tyrosine Kinase Inhibitor Lorlatinib in Hepatic Cells" International Journal of Molecular Sciences 23, no. 17: 9986. https://doi.org/10.3390/ijms23179986

APA Style

Verdura, S., Encinar, J. A., Fernández-Arroyo, S., Joven, J., Cuyàs, E., Bosch-Barrera, J., & Menendez, J. A. (2022). Silibinin Suppresses the Hyperlipidemic Effects of the ALK-Tyrosine Kinase Inhibitor Lorlatinib in Hepatic Cells. International Journal of Molecular Sciences, 23(17), 9986. https://doi.org/10.3390/ijms23179986

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop