J Ethnopharmacol I 他克莫司腎毒性如何破解？

2026年1月13日，海軍軍醫大學（第二軍醫大學）第二附屬醫院（上海長征醫院）藥學部陳萬生、陶霞、張鳳教授團隊在Journal of Ethnopharmacology (IF=5.4，中科院2區，Q1)在線發表題為“Integrated metabolomics and proteomics analysis elucidated the therapeutic effect of Huangkui Capsule on tacrolimus-induced chronic nephrotoxicity in rats”（整合代謝組學和蛋白質組學分析闡明了黃葵膠囊對大鼠他克莫司誘導的慢性腎毒性的治療作用）的研究論文。

海軍軍醫大學（第二軍醫大學）長征醫院藥學部龐濤主管藥師、博士研究生唐婷婷、碩士研究生王業健為本論文共同第一作者。

亮點：

?HK 可減輕體內他克莫司引起的慢性腎毒性和體外他克莫司引起的細胞毒性。

?HK 通過激活 xCT/Gpx4 通路抑制鐵死亡，從而對抗 TCN。

?來自HK的槲皮素、槲皮素-7-O-β-葡糖醛酸苷和檉柳素可通過上調Slc7a11/Slc3a2系統來增強GSH合成并抑制脂質過氧化。

民族藥理學相關性

黃葵膠囊（HK）在治療慢性腎臟病方面已顯示出顯著療效。既往研究表明，HK可以減輕腎臟損傷；然而，HK對他克莫司誘導的慢性腎毒性（TCN）的治療作用及其潛在的分子機制尚不明確。

研究目的

評價HK在緩解TCN方面的有效性，并探索其潛在的作用機制。

材料與方法

本研究采用生化指標和組織病理學檢查，在TCN大鼠模型上評價了HK的治療效果。通過整合腎臟代謝組學和蛋白質組學分析，揭示了TCN相關的獨特失調通路和潛在靶點，并深入探討了靶點與HK活性成分之間的關系。此外，還利用分子對接和Western blot分析對結果進行了驗證。

結果

TCN大鼠表現出蛋白尿增加和腎功能惡化，并伴有炎癥和氧化應激增強，而HK則以劑量依賴的方式改善腎損傷。整合代謝組學和蛋白質組學分析表明，HK通過恢復谷胱甘肽代謝和上調Gclc、Gclm、Slc7a11、Slc3a2和Gpx4來抑制TCN，并通過構建的“蛋白質-通路-代謝物”網絡分析揭示了鐵死亡的特異性調控機制。機制上，他克莫司干預增加了細胞內Fe2?水平，通過下調Slc7a11和Slc3a2抑制了xCT系統活性，并降低了Gpx4水平，從而增加了鐵死亡的易感性。來自 HK 的 Tamarixetin 和槲皮素-7-O-β-葡糖醛酸苷可以逆轉這些影響，恢復 Slc7a11、Slc3a2 和 Gpx4 的表達水平，表現出與 Fer-1 相當的效果。

結論

代謝組學和蛋白質組學綜合分析表明，HK通過抑制鐵死亡來減輕TCN。這種保護作用是由其生物活性成分介導的，特別是槲皮素、檉柳素和槲皮素-7-O-β-葡糖醛酸苷，它們通過上調Slc7a11/Slc3a2系統來增強GSH合成并抑制脂質過氧化。

01

研究背景及科學問題

他克莫司是一種具有強效免疫抑制作用的鈣調磷酸酶抑制劑，已被廣泛確立為腎移植受者一線治療藥物，用于減輕移植物排斥反應（Ong和Gaston，2021）。其長期腎毒性（TCN）是一個重要的臨床挑戰，即使他克莫司谷濃度維持在治療范圍內（5-8 ng/mL），仍有約20%的患者會出現TCN（Wojciechowski和Wiseman，2021；van Gelder等，2014；Kosoku等，2022）。TCN的早期表現通常包括腎功能受損和蛋白尿，這兩者都被認為是導致慢性移植物功能障礙的主要因素（Bentata，2020）。臨床證據表明，對血清和/或尿液樣本中特定代謝組學和蛋白質組學特征進行定量分析，能夠鑒別診斷腎功能損害，從而將這些分子特征確定為早期腎病檢測的潛在生物標志物（Guo et al., 2019;Medyńska et al., 2024）。目前，TCN的治療選擇仍然有限，主要干預措施是降低劑量或停藥。由此產生的臨床困境凸顯了開發安全有效的TCN治療策略并保護患者長期健康的迫切性。

黃葵（Abelmoschus manihot(L.) Medic.）是錦葵科的一種植物，因其富含黃酮類化合物而具有藥理學價值。其主要生物活性成分包括槲皮素、檉柳素、山奈酚及其衍生物（Wei et al., 2023）。在中國，黃秋葵在慢性腎臟?。–KD）的治療中有著悠久的傳統應用歷史（Luan et al., 2020;Li et al., 2021b;Chen et al., 2025）。黃秋葵的乙醇提取物制成黃葵膠囊（HK），具有顯著的腎臟保護作用（Zhang et al., 2014;Yan et al., 2025）。其治療CKD的作用機制涉及多種途徑，包括抗炎和抗氧化途徑。己酮可可堿（HK）因其抗炎作用而聞名，其作用機制是通過阻斷炎癥通路和降低炎癥標志物（Yan et al., 2025）。同時，HK還能通過增強抗氧化酶（如超氧化物歧化酶（SOD））的活性，并降低丙二醛（MDA）等毒性代謝物的水平來減輕氧化應激，從而減輕腎小管損傷（Li et al., 2021b）。此外，HK還能通過抑制免疫復合物的積累來重塑腎臟免疫微環境。這種重塑作用可以保護腎小球濾過屏障并減少尿蛋白漏（Wei et al., 2023）。臨床上，HK單藥治療及其與標準腎素-血管緊張素-醛固酮系統抑制劑聯合治療均已證實對慢性腎臟?。–KD）具有良好的療效。研究證實，這些治療方法能有效維持腎功能的長期穩定，延緩終末期腎病的進展，且不良事件發生率始終較低（Zhang et al, 2014;Yan et al, 2025）。此外，HK還能有效改善腎移植后患者的蛋白尿并降低炎癥因子水平（Yao et al, 2025;Xing et al., 2019;Wu et al., 2024）。盡管近年來HK的臨床和基礎研究不斷增加，但其在TCN中的具體腎臟保護作用仍有待明確。此外，HK發揮抗TCN保護作用的完整分子機制和信號通路網絡也尚未完全闡明。

代謝組學是一種通過分析小分子代謝物來表征病理狀態的寶貴工具。它已被廣泛用于闡明藥物不良反應通路，并日益融入毒理學研究和安全性藥理學評價中（Guo et al., 2022;Lim et al., 2023）。同樣，蛋白質組學也已成為另一種強大的系統生物學工具，用于識別與毒性相關的蛋白質組學改變，同時闡明其潛在的機制網絡（Deng et al., 2024;Fartade et al., 2025）。通過整合代謝組學和蛋白質組學數據，我們可以系統地分析細胞網絡和疾病機制，從而顯著增強我們對生物通路如何響應擾動的理解（Canzler et al., 2020;Dubin and Rhee, 2020;Hussain and Kingsley, 2025）。這種聯合組學方法已在我們之前的研究中得到有效應用，用于闡明化療引起的副作用的潛在機制并確定治療策略（Li et al., 2021a;Yao et al., 2022）。

迄今為止，代謝組學已成為闡明HK對糖尿病腎病和順鉑誘導的腎損傷的腎臟保護機制的關鍵工具（Yan et al., 2025;Han et al., 2025;Liao et al., 2025）。代謝組學分析表明，HK可通過調節鞘脂和視黃醇代謝等代謝物來改善腎功能障礙。借助質譜成像分析進行化合物鑒定或16S rRNA測序，研究人員還鑒定了核心治療靶點和藥物活性成分。這些研究為HK治療腎臟疾病的機制驗證提供了寶貴的思路和整體框架。

在本研究中，我們建立了大鼠TCN模型，并采用腎臟代謝組學分析系統地定量分析了HK治療后內源性代謝物的變化，從而闡明了與TCN病理生理相關的生化通路紊亂。同時，我們通過蛋白質組學分析鑒定了與TCN發病機制相關的潛在治療靶點。隨后，我們利用他克莫司誘導的腎損傷體外模型和分子對接技術驗證了HK類黃酮與這些靶點的相互作用。綜上所述，這一綜合策略闡明了HK對TCN的保護機制，為其治療潛力提供了新的分子層面的見解，并為未來HK介導的腎損傷代謝調控研究奠定了基礎。

02

重要發現及亮點

3.1 . HK治療顯著緩解了大鼠的TCN癥狀并改善了腎功能。

與正常組（NOR組）相比，中毒組（MOD組）大鼠的體重增加顯著降低（57.84 ± 16.27 g vs. 110.40 ± 3.71 g）（p< 0.001，圖S1）。HK治療逆轉了這一趨勢，其中HK-H組（2.8 g/kg/天）的恢復最為顯著。同時，與正常組相比，中毒組的血清肌酐（Scr）、尿素氮（BUN）、尿肌酐（Ucr）和24小時尿蛋白（24-h UProt）水平顯著升高（p< 0.001）。低劑量、中劑量和高劑量HK（HK-L、HK-M和HK-H）治療均顯著降低了這些指標，且呈明顯的劑量依賴性。值得注意的是，HK-H組的血清肌酐、尿素氮、尿素氮和24小時尿蛋白水平顯著低于正常組（圖1A-D）。組織學分析（H&E染色和Masson染色）顯示，MOD組腎小管管腔擴張、腎小球萎縮、炎性細胞浸潤、大量藍染膠原纖維以及明顯的間質纖維增生，表明他克莫司誘導了明顯的腎損傷。相比之下，所有HK治療組的腎損傷均有不同程度的減輕，證實了HK對他克莫司腎病（TCN）的保護作用（圖1E）。

圖 1.HK可顯著改善由他克莫司引起的慢性腎損傷導致的腎功能退化和腎組織病理損傷。

（AD）血清肌酐（Scr）、尿肌酐（Ucr）、血尿素氮（BUN）和24小時尿蛋白水平（n=8）。（E）蘇木精-伊紅（H&E）染色和Masson染色（×200）。（F）TUNEL染色（×40）。（GI）腎小球髓鞘形成蛋白-1（KIM-1）、中性粒細胞明膠酶相關脂質運載蛋白（NGAL）和8-羥基脫氧鳥苷（8-OHdG）水平（n=8）。（JN）腎臟中轉化生長因子-β1（TGF-β1）、I型膠原、α-平滑肌肌動蛋白（α-SMA）和caspase-3的蛋白表達（n=3）。###與正常組（NOR ）相比，p < 0.001；* 與中度偏頭痛組（MOD）相比，p< 0.05，** 與中度偏頭痛組（MOD）相比，p< 0.01，***與中度偏頭痛組（MOD）相比，p < 0.001。

3.2 . HK減輕了TCN大鼠腎組織細胞的凋亡，降低了腎損傷標志物水平和腎纖維化程度。

TUNEL染色結果顯示，MOD組細胞凋亡顯著增加，而HK-L、HK-M和HK-H治療可顯著降低細胞凋亡（圖1F）。與此一致，HK給藥后尿液中KIM-1、NGAL和8-OHdG的水平顯著降低（圖1G-I）（p< 0.01，p< 0.001）。高劑量HK使KIM-1從22.55 ± 3.2 pg/mL降至13.42 ± 2.10 pg/mL，NGAL從33.59 ± 1.09 ng/mL降至20.35 ± 1.41 ng/mL，8-OHdG從29.16 ± 2.25 ng/mL降至16.95 ± 2.31 ng/mL（p< 0.001）。 Western blotting進一步證實，與NOR組相比，他克莫司上調了MOD組中纖維化標志物（TGF-β1、膠原蛋白1和α-SMA）和凋亡標志物caspase-3的表達（p< 0.001）。HK治療表現出這些蛋白的劑量依賴性降低，有效改善了腎纖維化（圖1J-N）。

3.3 . HK改善了TCN大鼠的腎臟代謝紊亂

代謝組學數據的PCA分析表明，在正負離子模式下，質控樣品均呈現緊密聚類，證實了分析流程的可靠性（圖S2A-H）。PCA模型（圖2A和2B）顯示出清晰的組間分離，HK處理組的代謝譜向NOR組偏移。這一趨勢表明，HK處理后他克莫司誘導的代謝功能障礙得到逆轉。OPLS-DA進一步驗證了這種分離，在NOR組和MOD組之間以及MOD組和HK-H組之間觀察到了明顯的聚類（圖2C和2D）。代謝物鑒定后，分別從NOR組與MOD組以及MOD組與HK-H組的比較中鑒定出106種和125種差異代謝物（圖2E）。隨后的維恩圖分析鑒定出80種在兩組比較中均存在的差異代謝物（圖2F）。表S1總結了NOR組、MOD組和HK-H組中80種差異代謝物的變化。對這些代謝物進行KEGG通路富集分析，結果顯示33條通路顯著富集。關鍵的顯著富集通路（Impact > 0.1，p< 0.05）包括谷胱甘肽代謝、精氨酸和脯氨酸代謝、?；撬岷蛠喤；撬岽x以及精氨酸生物合成（圖2H）。在這些通路中鑒定出9種關鍵差異代謝物：L-脯氨酸、L-半胱氨酸、L-谷氨酸、L-精氨酸、γ-谷氨酰半胱氨酸、焦谷氨酸、肌酸、谷胱甘肽和?；撬幔▓D2I）。這些代謝物在NOR組、MOD組和HK-H組之間均表現出顯著變化。同樣，我們也分析了NOR組、MOD組和HK-L組或HK-M組之間的代謝譜變化。 OPLS模型清晰地展示了HK-L組（圖S3A和B）和HK-M組（圖S3C和D）與MOD組相比的代謝譜差異。差異代謝物篩選后，分別在HK-L組和HK-M組中鑒定出60種和62種差異代謝物（表S2和S3）。后續的KEGG通路富集分析表明，HK-L組和HK-M組中的差異表達代謝物顯著富集于四個共同通路：谷胱甘肽代謝、精氨酸和脯氨酸代謝、?；撬岷蛠喤；撬岽x以及精氨酸生物合成（影響值>0.1，p<0.05，圖S3B和D）。值得注意的是，這種富集模式與HK-H組的觀察結果高度一致，表明HK治療范圍內代謝反應具有一致性。

圖 2.非靶向腎臟代謝組學分析。

（A）正離子模式和（B）負離子模式下所有組的PCA二維得分圖；（C，D）正負離子模式下兩組之間的OPLS-DA二維得分圖和置換檢驗結果；（E）兩組之間上調和下調的代謝物；（F）差異代謝物的維恩圖；（G）80個差異代謝物的熱圖；（H）差異代謝物通路圖；（I）關鍵通路中9個差異代謝物的熱圖。

3.4 . HK對TCN大鼠腎臟緩解作用的蛋白質組學特征分析

質量控制分析證實了蛋白質組學數據的可靠性（圖S4A-D）。在NOR組和MOD組之間共鑒定出364個差異表達蛋白（DEPs）（|log2 FC| > 1，p < 0.05），其中154個上調，210個下調（圖3A）。同樣，在MOD組和HK-H組之間鑒定出145個DEPs，其中64個上調，81個下調（圖3B）。熱圖可視化了DEPs的表達模式，而維恩圖則展示了兩個比較中共有的112個蛋白（圖3C和D）?；虮倔w（GO）富集分析揭示了各組之間不同的分子特征（圖3E）。具體而言，DEPs主要參與細胞氨基酸分解代謝和谷胱甘肽代謝等生物學過程。這些分子功能包括抗氧化/氧化還原酶、谷氨酸-半胱氨酸連接酶和跨膜轉運蛋白活性。細胞成分包括細胞外囊泡和線粒體基質。KEGG通路分析顯示，差異表達蛋白富集于十條通路（p< 0.05），包括鐵死亡、谷胱甘肽代謝以及蛋白質消化吸收（圖3F，表S4）。隨后，將這些通路中富集的24個蛋白組成一個關鍵蛋白簇，用于后續的整合分析。

圖 3.MOD組的腎臟蛋白質組學特征與 NOR 組或 HK-H 組不同。

火山圖：（A）NOR vs. MOD 和（B）MOD vs. HK-H；（C）差異蛋白的維恩圖；（D）GO 功能分析結果；（E）差異蛋白的熱圖；（F）KEGG 通路的?；鶊D；（G）顯著差異代謝物與蛋白相關性分析的熱圖；（H）蛋白-KEGG 通路-代謝物關聯圖，篩選出 7 條主要通路；（I）5 個關鍵蛋白和（J）3 個關鍵代謝物的表達水平變化；#p<0.001，###p<0.001 vs NOR 組；*p<0.05，**p<0.01，***p<0.001 vs MOD 組。

3.5 .綜合代謝和原生質分析

相關性分析鑒定出7種關鍵差異代謝物，包括肌酸、谷胱甘肽、L-精氨酸、L-半胱氨酸、L-谷氨酸、L-脯氨酸和?；撬?，以及19種關鍵差異蛋白，包括Aadat、C1qb、Cndp1、Ddo、Folh1、Gclc、Gclm、Gpx1、Gpx4、Hnmt、Mep1a、Mme、Mrc2、Slc3a2、Slc6a19、Slc7a11、Tap1、Tubb6和Xpnpep2，這些蛋白與TCN的病理生理學相關（圖3G，表1）。通路富集分析表明，多個關鍵過程存在顯著失調，包括鐵死亡、蛋白質消化吸收以及與谷胱甘肽、組氨酸和氨基酸相關的多種代謝通路。進一步構建的“蛋白質-通路-代謝物”網絡分析表明，HK 通過增強谷胱甘肽合成來抑制鐵死亡，從而減弱 TCN（圖 3H-J）。

3.6 . HK抑制TCN中氧化應激誘導的鐵死亡

與NOR組相比，MOD組的抗氧化標志物水平顯著降低，包括GSH、SOD和CAT。相反，脂質過氧化產物MDA、鐵死亡相關標志物Fe2?以及DHE熒光（ROS的指標）的水平顯著升高（p< 0.001）。HK治療呈劑量依賴性地減輕了這些改變，表現為GSH、SOD和CAT水平的恢復，以及MDA、Fe2?和DHE熒光水平的降低（圖4A-F），其中HK-H組的效果最為顯著。同時，Western blot分析顯示，MOD組大鼠腎組織中關鍵鐵死亡調節因子（Slc7a11、Gpx4、Slc3a2）的蛋白表達顯著下調。 HK 治療可劑量依賴性地逆轉這種下調（p< 0.05，p< 0.01，p< 0.001），HK-H 組可顯著地將 Slc7a11 和 Gpx4 的表達恢復到接近正常水平（p< 0.001）（圖 4G-J）。

圖 4.腎臟組織中鐵死亡和谷胱甘肽代謝相關標志物。

（AE）SOD、MDA、Fe2+、CAT 和 GSH 水平（n=8）；（F）ROS 表達水平（×40）（n=3）；（G）谷胱甘肽代謝和鐵死亡關鍵蛋白條帶的蛋白質印跡結果；（HJ）Gpx4、Slc3a2 和 Slc7a11 的表達結果；###p<0.001 vs NOR 組；*p<0.05，**p<0.01，***p<0.001 vs MOD 組。

3.7 . HK類黃酮抑制他克莫司誘導的HK-2細胞損傷中氧化應激誘導的鐵死亡

我們采用先前研究（ Zhang et al., 2025）中建立的方法鑒定了六種HK類黃酮，并據此進行了重新分析（圖S5和S6，表S5）。它們對HK-2細胞活力的影響如圖S7所示。他克莫司處理顯著降低了HK-2細胞的活力，而六種黃酮類化合物中的任何一種都能減輕這種降低，其程度與陽性對照Fer-1相當（圖5A，p< 0.001）。此外，他克莫司在HK-2細胞中誘導了顯著的氧化應激反應，表現為MDA水平升高、GSH水平降低以及SOD和CAT活性降低（p< 0.001）。Fer-1的聯合給藥證實了鐵死亡的參與，Fer-1能有效逆轉這些生化改變。對HK類黃酮的比較評估發現，槲皮素、檉柳素和槲皮素-7-O-β-葡糖醛酸苷是鐵死亡最有效的抑制劑，因為它們顯著地使所有四種氧化應激標志物的水平恢復正常。山奈酚-3-O-β-D-葡糖醛酸苷和米奎利寧也能降低MDA水平并提高GSH水平，但它們對SOD和CAT活性的影響較?。▓D5B-E）。此外，Fer-1、槲皮素、檉柳素、槲皮素-7-O-β-葡糖醛酸苷和山奈酚-3-O-β-D-葡糖醛酸苷均能顯著減弱他克莫司誘導的Fe2?積累（圖5F）。透射電鏡提供了鐵死亡的超微結構證據，表明他克莫司導致嚴重的線粒體損傷，包括線粒體斷裂、外膜破裂和嵴丟失。與Fer-1或黃酮類化合物聯合治療可顯著改善這些形態學缺陷，其中槲皮素、檉柳素和槲皮素-7-O-β-葡糖醛酸苷的保護作用最為顯著（圖5G）。

圖 5.他克莫司誘導的 HK-2 損傷細胞中鐵死亡相關標志物以及 HK 黃酮類化合物和 Fer-1 的發現。

(A) 細胞活力； (BF) SOD、MDA、CAT、GSH 和 Fe2+水平； (G)透射電子顯微鏡(×8000)； TAC，他克莫司；闕，槲皮素； Tam，檉柳西汀； Iso，異鼠李素； Que-7-Glu、槲皮素-7-O-β-葡萄糖苷酸； Kae-D-Glu、山奈酚-3-O-β-D-葡萄糖苷酸； Miq，米奎寧。###與 NOR 組相比，p<0.001；與 MOD 組相比， *p<0.05，***p<0.001。

3.8分子對接分析和分子動力學模擬

六種HK類黃酮與Slc7a11（圖6A）和Slc3a2（圖6B）表現出很強的結合親和力，結合能低于-7.0 kcal/mol。其中三種黃酮的親和力最高：槲皮素與Slc7a11的結合能為-9.1 kcal/mol，與Slc3a2的結合能為-7.9 kcal/mol；檉柳素與Slc7a11的結合能為-9.2 kcal/mol，與Slc3a2的結合能為-7.9 kcal/mol；槲皮素-7-O-β-葡糖醛酸苷與Slc7a11的結合能為-9.5 kcal/mol，與Slc3a2的結合能為-8.7 kcal/mol（表2）。分子動力學模擬（模擬時間超過100納秒）表明，三種黃酮類化合物與Slc7a11（圖6C）和Slc3a2（圖6D）均能穩定結合，這體現在其較低的RMSD和RMSF值上。在這些黃酮類化合物中，槲皮素-7-O-β-葡糖醛酸苷與這兩種蛋白形成的復合物最為穩定。

圖 6.分子對接模擬和分子動力學模擬。

槲皮素、檉柳黃素、異鼠李素、槲皮素-7-O-β-葡萄糖苷酸、山奈酚-3-O-β-D-葡萄糖苷酸和米奎寧與(A) Slc7a11和(B) Slc3a2結合； (C) Slc7a11 和 (D) Slc3a2 的 RMSD 圖表。 TAC，他克莫司；闕，槲皮素； Tam，檉柳西汀； Iso，異鼠李素； Que-7-Glu、槲皮素-7-O-β-葡萄糖苷酸； Kae-D-Glu、山奈酚-3-O-β-D-葡萄糖苷酸； Miq，米奎寧。

表2.化合物-蛋白質復合物的結合親和力

3.9 . HK通過xCT/GPX-4信號通路抑制鐵死亡

在HK-2細胞中，他克莫司暴露顯著抑制了Gpx4、Slc7a11和Slc3a2的表達水平（p< 0.001）。槲皮素、槲皮素-7-O-β-葡糖醛酸苷、檉柳素和山奈酚-3-O-β-D-葡糖醛酸苷均能有效逆轉這種抑制作用（圖7A-D）。這些結果證實，HK及其活性成分在體內外均能通過xCT/Gpx4信號通路拮抗他克莫司誘導的鐵死亡。

圖7.類黃酮影響HK-2 細胞中 Slc7a11、Gpx4 和 Slc3a2 的蛋白表達水平。(A) 蛋白質條帶； (B)Slc7a11； (C)GPX4； (D)Slc3a2。 TAC，他克莫司；闕，槲皮素； Tam，檉柳西??； Iso，異鼠李素； Que-7-Glu、槲皮素-7-O-β-葡萄糖苷酸； Kae-D-Glu、山奈酚-3-O-β-D-葡萄糖苷酸； Miq，米奎寧。###與 NOR 組相比，p<0.001；與 MOD 組相比，*p<0.05，**p<0.01，***p <0.001。

【Citation】Pang T, Tang T, Wang Y, Li J, Yang L, He Y, Yang H, Han J, Tao X, Chen W, Zhang F. Integrated metabolomics and proteomics analysis elucidated the therapeutic effect of Huangkuui Capsule on tacrolimus-induced chronic nephrotoxicity in rats.Journal of Ethnopharmacology.2026.1.13.