ObjectiveTo investigate the mechanism of liver injury induced by tripterygium glycosides tablets （TG） and the molecular mechanism of compound glycyrrhizin tablets （CG） in alleviating the abnormalities of cholesterol metabolism caused by TG via cholesterol metabolism. MethodsAccording to the body weights， male Sprague-Dawley （SD） rats were randomly grouped as follows： control （pure water）， low-dose TG （TG-L， 189.0 mg·kg^-1·d^-1）， high-dose TG （TG-H， 472.5 mg·kg^-1·d^-1）， TG-L+CG （189.0 mg·kg^-1·d^-1 TG + 20.25 mg·kg^-1·d^-1 CG）， and TG-H+CG （472.5 mg·kg^-1·d^-1 TG + 20.25 mg·kg^-1·d^-1 CG）， with 6 rats in each group. Rats were administrated with corresponding drugs once daily for 3 weeks. At the end of the last administration， the mRNA and protein levels of liver X receptor-alpha （LXR-α）， low-density lipoprotein receptor （LDLR）， adenosine triphosphate-binding cassette transporter A1 （ABCA1）， adenosine triphosphate-binding cassette transporter G1 （ABCG1）， 3-hydroxy-3-methylglutaryl coenzyme A reductase （HMGCR）， cholesterol 7α-hydroxylase （CYP7A1）， cholesterol 12α-hydroxylase （CYP8B1）， and sterol 27-hydroxylase （CYP27A1） in the liver tissue were determined by Real-time PCR and Western blotting， respectively. The level of 3-hydroxy-3-methylglutaryl coenzyme A reductase （HMG-CoAR）， a regulatory enzyme of cholesterol synthesis， was measured by enzyme-linked immunosorbent assay （ELISA）. HepG2 cells were used to observe the effect of TG on the cell proliferation in vitro. Specifically， HepG2 cells were grouped as follows： Low-dose TG （TG-l， 15 mg·L^-1）， medium-dose TG （TG-m， 45 mg·L^-1）， high-dose TG （TG-h， 135 mg·L^-1）， fenofibrate （FB， 10 μmol·L^-1）， CG extract， TG-h+FB （135 mg·L^-1 TG + 10 μmol·L^-1 FB）， TG-m+FB （45 mg·L^-1 TG + 10 μmol·L^-1 FB）， TG-l+FB （15 mg·L^-1 TG + 10 μmol·L^-1 FB）， TG-h+CG （135 mg·L^-1 TG + 60 μmol·L^-1 CG）， TG-m+CG （45 mg·L^-1 TG + 60 μmol·L^-1 CG）， and TG-l+CG （15 mg·L^-1 TG + 60 μmol·L^-1 CG）. The mRNA and protein levels of LXR-α， ABCG1， LDLR， CYP7A1， CYP8B1， and CYP27A1 in HepG2 cells were determined by Real-time PCR and Western blotting， respectively. ResultsThe rat experiment showed that compared with the control group， the TG-H group showed down-regulated mRNA levels of CYP7A1， CYP8B1， and CYP27A1 in the liver tissue （P<0.05， P<0.01）， which were up-regulated by the application of CG （P<0.05， P<0.01）， and the TG-H+CG group showed up-regulated mRNA level of LDLR （P<0.01）. Compared with the control group， the TG-L and TG-H groups showed down-regulated protein levels of LDLR， CYP7A1， and CYP8B1 in the liver tissue （P<0.05， P<0.01）. In addition， the protein levels of ABCG1 and LXR-α were down-regulated in the TG-H and TG-L groups， respectively （P<0.05）. Compared with TG alone， TG+CG up-regulated the protein levels of ABCG1 and LDLR （P<0.05， P<0.01）， and the protein levels of CYP7A1 and CYP8B1 in the TG-H+CG group were up-regulated （P<0.05， P<0.01）. The cell experiment showed that compared with the control group， the TG-h group presented up-regulated mRNA level of LXR-α （P<0.01）， and the TG-m and TG-h groups showcased down-regulated mRNA levels of LDLR and CYP7A1 （P<0.01） and up-regulated mRNA level of CYP27A1 （P<0.01） in HepG2 cells. The combination of CG with TG restored the above changes （P<0.01）. Western blotting results showed that compared with the control group， the TG-m and TG-h groups showed down-regulated protein levels of LXR-α， ABCG1， LDLR， CYP7A1， CYP8B1， and CYP27A1 in HepG2 cells （P<0.01）. Compared with the TG-h group， the TG-h+CG group showed up-regulated protein level of LDLR （P<0.05）. Compared with the TG-m group， the TG-m+CG group showcased up-regulated protein levels of LDLR， ABCG1， CYP7A1， and CYP27A1 （P<0.05， P<0.01）. ConclusionThe administration of TG at 189.0， 472.5 mg·kg^-1 for 3 weeks could modulate the signaling pathways associated with cholesterol efflux， endocytosis， and cholesterol biotransformation in hepatocytes， leading to the accumulation of cholesterol and subsequent liver injury in rats. CG could ameliorate the liver injury induced by lipid metabolism disorders caused by TG by up-regulating the expression of LXR-α， LDLR， ABCG1， CYP7A1， CYP8B1， and CYP27A1 to promote cholesterol biotransformation.

OBJECTIVE: To analyze the acupoint selection patterns and core prescriptions of acupuncture for polycystic ovary syndrome (PCOS) using data mining, and to explore the molecular mechanisms of core acupoints through network acupuncture medicine. METHODS: The randomized controlled trials (RCTs) on acupuncture for PCOS published from January 1, 2004 to July 21, 2024 were retrieved from CNKI, VIP, Wanfang, PubMed, and Web of Science databases. R software (version 4.4.0) was used for acupoint frequency and association rule analysis to identify core acupoint prescriptions. Potential targets were predicted via the STITCH and Swiss Target Prediction databases, and a "core prescription-active compounds-targets- PCOS" network was constructed. Cytoscape 3.7.1 was applied to build protein-protein interaction (PPI) networks of potential targets of core acupoint prescriptions. Key therapeutic targets were subjected to gene ontology (GO) and Kyoto encyclopedia of genes and genomes (KEGG) enrichment analyses using the DAVID and Microbioinformatics platforms. RESULTS: A total of 176 RCTs were included, covering 208 prescriptions and 89 acupoints. The five most frequently used acupoints were Guanyuan (CV4), Sanyinjiao (SP6), Zigong (EX-CA1), Zusanli (ST36) and Zhongji (CV3). Association rule analysis yielded 13 core acupoint combinations, with Guanyuan (CV4), Sanyinjiao (SP6), Zigong (EX-CA1) and Zusanli (ST36) as the core prescription. Twenty-seven active compounds were involved, with 852 potential therapeutic targets, among which 208 targets overlapped with PCOS-related targets. Network acupuncture medicine analysis suggested that the core prescription may act through targets such as estrogen receptor 1 (ESR1), proto-oncogene tyrosine-protein kinase Src (SRC), signal transducer and activator of transcription 3 (STAT3), peroxisome proliferator-activated receptor gamma (PPARG), and RAC-alpha serine/threonine-protein kinase (AKT1). GO and KEGG analyses indicated that the main pathways included the hypoxia-inducible factor 1 (HIF-1) signaling pathway, phosphatidylinositol 3-kinase-protein kinase B (PI3K-AKT) signaling pathway, and advanced glycation end products-receptor for advanced glycation end products (AGE-RAGE) signaling pathway, involving processes such as signal transduction, receptor complex formation, and cytokine activity. CONCLUSION The core acupoint prescription for PCOS might exert therapeutic effects through multiple targets and pathways, providing a theoretical basis for mechanistic research on acupoint prescriptions.
Humans ; Acupuncture Therapy ; Data Mining ; Acupuncture Points ; Polycystic Ovary Syndrome/metabolism* ; Female ; Protein Interaction Maps ; Randomized Controlled Trials as Topic

Proteolysis-targeting chimeras (PROTACs) represent a promising class of drugs that can target disease-causing proteins more effectively than traditional small molecule inhibitors can, potentially revolutionizing drug discovery and treatment strategies. However, the links between in vitro and in vivo data are poorly understood, hindering a comprehensive understanding of the absorption, distribution, metabolism, and excretion (ADME) of PROTACs. In this work, ¹⁴C-labeled vepdegestrant (ARV-471), which is currently in phase III clinical trials for breast cancer, was synthesized as a model PROTAC to characterize its preclinical ADME properties and simulate its clinical pharmacokinetics (PK) by establishing a physiologically based pharmacokinetics (PBPK) model. For in vitro-in vivo extrapolation (IVIVE), hepatocyte clearance correlated more closely with in vivo rat PK data than liver microsomal clearance did. PBPK models, which were initially developed and validated in rats, accurately simulate ARV-471's PK across fed and fasted states, with parameters within 1.75-fold of the observed values. Human models, informed by in vitro ADME data, closely mirrored postoral dose plasma profiles at 30 mg. Furthermore, no human-specific metabolites were identified in vitro and the metabolic profile of rats could overlap that of humans. This work presents a roadmap for developing future PROTAC medications by elucidating the correlation between in vitro and in vivo characteristics.

Objective:To investigate the metabolism-related proteins and their presence in the plasma of female androgenetic alopecia (FAGA) patients.Methods:From March 2021 to March 2023, FAGA patients aged 18-50 (FAGA group) and healthy women (HC group) were recruited from the Dermatology Outpatient Department of Huashan Hospital. 3 ml of peripheral venous blood was collected from each participant and centrifuged to obtain plasma. Olink proteomics analysis was performed on the collected plasma, differentially expressed proteins were screened with R language, the diagnostic accuracy of the differentially expressed proteins was assessed using receiver operating characteristic (ROC) curve. Gene ontology (GO) analysis was performed on differentially expressed proteins. Immunofluorescence analysis on hair follicles in the parietal region of the FAGA group and the occipital region of the HC group was performed to validate the differentially expressed proteins identified. SPSS 25.0 software was used to analyze the data, with normal distribution metric data represented by Mean±SD. Student’s t-test was used to compare the basic information of two groups of subjects and the relative fluorescence intensity of differentially expressed proteins in hair follicles. Pearson correlation analysis was performed on plasma metabolism-related proteins and the basic information of subjects. P<0.05 indicates a statistically significant difference. Results:Sixty-one cases were included in the FAGA group, with an average age of (33.8±7.4) years and an onset age of (29.5±7.8) years. Among them, 38 cases were mild FAGA, 14 cases were moderate, and 9 cases were severe. Twenty-seven cases were included in the HC group, with an average age of (32.0±7.7) years. There was no statistically significant difference in the basic information (age, body mass index, testosterone, 25-hydroxyvitamin D, uric acid, and ferritin levels) between the two groups of subjects ( P>0.05). Compared to the HC group, the plasma of the FAGA group showed 26 significantly upregulated differentially expressed proteins ( P<0.05), with AHCY and NECTIN2 exhibiting the most significant differences (all P=0.003). The ROC curve evaluation revealed that the area under the curve for AHCY and NECTIN2 was greater than 0.7, indicating good diagnostic accuracy. The GO analysis revealed that the differentially expressed proteins were primarily enriched in the BAT3 complex (cellular component), ubiquitin-dependent ERAD pathway, natural killer cell activation (biological process), as well as ubiquitin protein ligase binding and ubiquitin-specific protease binding (molecular function). Pearson correlation analysis revealed that AHCY ( r=-0.23, P=0.010) and NECTIN2 ( r=-0.31, P=0.033) were negatively correlated with the severity of hair loss in FAGA patients. The results of hair follicle immunofluorescence analysis showed that the relative fluorescence intensity of AHCY and NECTIN2 in the FAGA group was higher than that in the HC group ( P<0.05). In other words, both AHCY and NECTIN2 were upregulated in the FAGA group. Conclusion:Metabolism-related proteins play an important role in FAGA. AHCY and NECTIN2 may serve as early diagnostic biomarkers for FAGA.

Objective:To investigate the metabolism-related proteins and their presence in the plasma of female androgenetic alopecia (FAGA) patients.Methods:From March 2021 to March 2023, FAGA patients aged 18-50 (FAGA group) and healthy women (HC group) were recruited from the Dermatology Outpatient Department of Huashan Hospital. 3 ml of peripheral venous blood was collected from each participant and centrifuged to obtain plasma. Olink proteomics analysis was performed on the collected plasma, differentially expressed proteins were screened with R language, the diagnostic accuracy of the differentially expressed proteins was assessed using receiver operating characteristic (ROC) curve. Gene ontology (GO) analysis was performed on differentially expressed proteins. Immunofluorescence analysis on hair follicles in the parietal region of the FAGA group and the occipital region of the HC group was performed to validate the differentially expressed proteins identified. SPSS 25.0 software was used to analyze the data, with normal distribution metric data represented by Mean±SD. Student’s t-test was used to compare the basic information of two groups of subjects and the relative fluorescence intensity of differentially expressed proteins in hair follicles. Pearson correlation analysis was performed on plasma metabolism-related proteins and the basic information of subjects. P<0.05 indicates a statistically significant difference. Results:Sixty-one cases were included in the FAGA group, with an average age of (33.8±7.4) years and an onset age of (29.5±7.8) years. Among them, 38 cases were mild FAGA, 14 cases were moderate, and 9 cases were severe. Twenty-seven cases were included in the HC group, with an average age of (32.0±7.7) years. There was no statistically significant difference in the basic information (age, body mass index, testosterone, 25-hydroxyvitamin D, uric acid, and ferritin levels) between the two groups of subjects ( P>0.05). Compared to the HC group, the plasma of the FAGA group showed 26 significantly upregulated differentially expressed proteins ( P<0.05), with AHCY and NECTIN2 exhibiting the most significant differences (all P=0.003). The ROC curve evaluation revealed that the area under the curve for AHCY and NECTIN2 was greater than 0.7, indicating good diagnostic accuracy. The GO analysis revealed that the differentially expressed proteins were primarily enriched in the BAT3 complex (cellular component), ubiquitin-dependent ERAD pathway, natural killer cell activation (biological process), as well as ubiquitin protein ligase binding and ubiquitin-specific protease binding (molecular function). Pearson correlation analysis revealed that AHCY ( r=-0.23, P=0.010) and NECTIN2 ( r=-0.31, P=0.033) were negatively correlated with the severity of hair loss in FAGA patients. The results of hair follicle immunofluorescence analysis showed that the relative fluorescence intensity of AHCY and NECTIN2 in the FAGA group was higher than that in the HC group ( P<0.05). In other words, both AHCY and NECTIN2 were upregulated in the FAGA group. Conclusion:Metabolism-related proteins play an important role in FAGA. AHCY and NECTIN2 may serve as early diagnostic biomarkers for FAGA.

Objective:To investigate the risk factors of postoperative recurrence in patients with primary brain glioma and to construct a prediction model.Methods:A total of 98 patients with primary brain glioma treated by radical surgery in Nanchong Central Hospital from January 2018 to January 2021 were retrospectively included, and were divided into recurrent group (40 cases) and non-recurrent group (58 cases) according to whether there was recurrence or not during the follow-up period. The independent influencing factors for postoperative recurrence in patients with primary brain glioma were evaluated by multivariate logistic regression. Logistic prediction model of postoperative recurrence risk of patients with primary brain glioma was established, and the predictive efficacy of each index was calculated by receiver operator characteristic (ROC) curve.Results:There were statistically significant differences between recurrent group and non-recurrent group in glioma World Health Organization (WHO) grade ( χ2=12.48, P<0.001), isocitrate dehydrogenase (IDH) 1/2 mutation ( χ2=13.24, P<0.001), mean platelet volume (MPV) ( t=5.34, P<0.001), and MPV/platelet count (PLT) ( t=9.73, P<0.001). Multivariate analysis showed that WHO grade Ⅲ-Ⅳ ( OR=8.54, 95% CI: 1.62-44.99, P=0.011), IDH1/2 wild type ( OR=9.08, 95% CI: 1.68-49.19, P=0.010), low MPV ( OR=0.46, 95% CI: 0.21-0.99, P=0.048) and low MPV/PLT ( OR=0.02, 95% CI: 0.01-0.03, P<0.001) were independent risk factors for postoperative recurrence in patients with primary brain glioma. The logistic prediction model based on the above indicators was logit ( P) =11.78+2.15×WHO grade+2.21×IDH1/2 mutation situation-0.78×MPV-200.70×MPV/PLT ( R2=0.785). The ROC curve analysis results showed that WHO grade, IDH1/2 mutation, MPV, MPV/PLT, and logistic prediction model P-value could all be used to predict the risk of postoperative recurrence in patients with primary brain glioma; The areas under the curve were 0.681, 0.684, 0.783, 0.920 and 0.964, respectively. In the area under ROC curve comparison of each indicator, the predictive performance of the logistic regression model P-value was significantly higher than that of other indicators (all P<0.05) . Conclusion:Postoperative recurrence in patients with primary brain glioma may be related to glioma WHO grade, IDH1/2 mutation and platelet related laboratory indexes. The model constructed based on the above indicators can be used to predict the recurrence risk in patients with primary brain glioma after surgery.

OBJECTIVE To study the improvement effect of total flavonoids from Rosa multiflora root on vascular injury in rheumatoid arthritis （RA） model rats and its potential mechanism. METHODS Female Wistar rats were randomly divided into normal control group， model group， aspirin group （positive control， 30 mg/kg）， low-dose and high-dose groups of total flavonoids from R. multiflora root （4.15， 8.30 g/kg， by crude drug）， with 10 rats in each group. Except for the normal control group， the RA model was induced in other groups by collagen induction and high-fat diet. After 14 days of modeling， they were given corresponding drug solution/0.5% sodium carboxymethyl cellulose solution intragastrically， once a day， for 36 consecutive days. The total body score， arthritis index （AI） and swollen joint count （SJC） of the rats were evaluated regularly. After the last medication， serum levels of interleukin-6 （IL-6）， intercellular adhesion molecule-1 （ICAM-1） and vascular cell adhesion molecule- 1 （VCAM-1） were determined. The pathological morphological changes in the vascular tissue of thoracic aorta were observed； the protein expression of Toll-like receptor 4 （TLR4） and the protein phosphorylation levels of Janus kinase 2 （JAK2） and signal transduction and activator of transcription 3 （STAT3） in vascular tissue of thoracic aorta were measured. RESULTS Compared with the normal control group， serum levels of IL-6， ICAM-1 and VCAM-1， protein expression of TLR4， and the protein phosphorylation levels of JAK2 and STAT3 in vascular tissue of thoracic aorta were increased significantly in model group （P＜ 0.01）. The atherosclerotic plaque （atheroma）， cholesterol crystal， lymphocyte infiltration and a small number of unbroken foam cell aggregation could be seen in the vascular tissue of thoracic aorta. Compared with the model group， total body score （except for the low-dose group）， AI and SJC were decreased significantly in groups of total flavonoids from R. multiflora root on the 28th day （P＜0.05 or P＜0.01）； total body score，AI and SJC were decreased significantly in low-dose group of total flavonoids from R. multiflora root on the 49th day （P＜0.05 or P＜0.01）； the other quantitative indicators in serum and vascular tissue were significantly reversed in groups of total flavonoids from R. multiflora root （P＜0.05 or P＜0.01）， and pathological damage of vascular tissue was significantly relieved. CONCLUSIONS Total flavonoids from R. multiflora root can significantly improve vascular injury in RA model rats， and its mechanism may be related to reducing the protein expression of TLR4 in vascular tissue and inhibiting the activation of IL-6/JAK2/ STAT3 signaling pathway.

Background According to the Classification and Catalogue of Occupational Diseases, brucellosis is one of the notifiable occupational infectious diseases, which occurs from time to time in the occupational population. Objective To compare the work-related injury appraisal process and results of 13 cases of brucellosis at both provincial and municipal levels, analyze and summarize the bias in the practical work of labor capacity identification for occupational diseases such as brucellosis by appraisal management departments and experts, and propose suggestions for optimizing appraisal work. Methods A comparative study was conducted on the objective examination results and labor capacity appraisal conclusions based on the occupational contact history, clinical diagnosis, occupational disease diagnosis staging, and labor capacity appraisal of 13 patients with brucellosis. The reasons for the inconsistency between the initial appraisal conclusion by institutions at the municipal level and the final appraisal conclusion by institutions at the provincial level were compared and analyzed. Results All of the 13 patients with brucellosis applied for municipal-level labor capacity identification after being identified as work-related injuries, 11 of which did not receive a disability level, and 2 were rated as level 10 disability. Four of those who did not receive the disability rate applied for provincial-level labor capacity identification. As a result, 2 cases were maintained original appraisal conclusions, while the other 2 changed the conclusions to level 9 disability and level 10 disability respectively. It was the first time in Shijiangzhuang municipal-level primary labor capacity appraisal and Hebei provincial-level labor capacity re-appraisal that the work-related injury caused by occupational brucellosis was rated as level 9 disability or level 10 disability. Hence, the lessons learned from this comparative analysis are that the degree of target organ damage and (or) organ dysfunction are the direct basis for work-related injury appraisal; an objective and scientific labor capacity identification for occupational brucellosis should base on the each case of disability evaluation, identify the relevant target organ damage and the degree of dysfunction, and rate the disability level after a comprehensive appraisal. Conclusion This analysis would be a guidance to the identification of labor capacity for occupational brucellosis in Hebei Province and the whole country. There is a hysteresis issue in the occupational disease provisions in the national standard GB/T 16180—2014 Standard for identify work ability—Gradation of disability caused by work-related injuries and occupatiaonal diseases. In current situation, appraisal experts should not only search for clauses that directly correspond to the occupational diseases and injuries, but also target conditions not covered in the clauses and conduct assessment based on the characteristics of occupational diseases, with scientific, accurate, and flexible application of the clauses in the standard and appendix, so as to make fair, just, and professional appraisal conclusions.

ObjectiveTo explore the characteristics and regulatory mechanisms of cholestatic liver injury （CLI） caused by traditional Chinese medicine based on data mining， network pharmacology， and molecular docking. MethodChina National Knowledge Infrastructure and PubMed were searched for the relevant literature on CLI caused by traditional Chinese medicine from inception to 2024， and the information was standardized， summarized， and analyzed by cluster analysis. The Traditional Chinese Medicine Systems Pharmacology Database and Analysis Platform （TCMSP）， Integrative Pharmacology-based Research Platform of Traditional Chinese Medicine （TCMIP）， GeneCards， Online Mendelian Inheritance in Man （OMIM）， and DisGeNET were used to retrieve the active ingredients and targets of core medicines. The Venn diagram was established to map the common targets shared by the core medicines and CLI. Cytoscape 3.10.2 was used to construct the protein-protein interaction （PPI） network of the common targets. DAVID was used for gene ontology （GO） and Kyoto Encyclopedia of Genes and Genomes （KEGG） enrichment of the core targets. Finally， molecular docking was performed by AutoDock Vina. ResultA total of 849 eligible articles were included in this study， from which 64 active ingredients of 39 herbal medicines that can cause cholestasis were counted and categorized according to the 2020 edition of the Pharmacopoeia of the People's Republic of China. The frequency of the medidicnes followed a descending order of heat-clearing medicines， exterior-releasing medicines， blood-activating and stasis-resolving medicines， purgative medicines， phlegm-resolving and cough- and dyspnea-relieving medicines， tonics， wind- and dampness-expelling medicines， interior-warming medicines， urination-promoting and dampness-draining medicines， Qi movement-regulating medicines， hemostatics， toxin-removing， worm-killing， and itch-relieving medicines， astringent medicines， dampness-eliminating medicines， and tranquiling medicines. The cluster analysis revealed that the reports of CLI caused by heat-clearing medicines accounted for the highest proportion of 39.69%. Among the heat-clearing medicines， Gardeniae Fructus （92 articles， accounting for 10.84%）， Scutellariae Radix （76 articles， 8.95%）， and Sophorae Flavescentis Radix （69 articles， 6.95%） were frequently reported. The core targets of cholestasis induced by Chinese herbal medicines reported to cause CLI mainly included tumor necrosis factor （TNF）， peroxisome proliferator activated receptor alpha （PPARA）， farnesoid X receptor （FXR）， glutamic-pyruvic transaminase 2 （GPT2）， superoxide dismutase 1 （SOD1）， interferon gamma （IFN-γ）， interleukin-6（IL-6）， CD36， Apolipoprotein A1（APOA1）， Angiotensin converting enzyme（ACE）， Cytochrome P450 3A4 enzyme（CYP3A4）， Protein kinase B1（Akt1）， APOB， albumin（ALB）， ATP binding cassette transporter A4（ABCB4）， SLC10A1， Estrogen Receptor alpha （ESR1）， signal transducer and activator of transcription1（STAT1）， β-actin（ACTB）， Endothelin 1（EDN1）， ABCG2， and peroxisome proliferator activated receptor gamma（PPARG）. The signaling pathways involved included bile secretion， ABC transporter， steroid biosynthesis， DNA adducts， drug metabolism， cytochrome P450， peroxisomes， primary bile acid biosynthesis， retinol metabolism， and Toll-like receptor. The molecular docking results showed that the active ingredients （e.g.， baicalin and berberine） of the heat-clearing medicines reported by high frequency to cause CLI had high binding affinity to the targets including ABCG2， IFN-γ， EDN1， IL-6 and SOD1， with the binding energy in the range of -13 kcal·mol^-1 to -9 kcal·mol^-1， and the regulatory pathways were highly correlated with transporters， microvascular function regulation， inflammation， and oxidative stress， which was consistent with the cluster analysis. ConclusionThe available reports about the Chinese herbal medicines causing CLI mainly focused on heat-clearing medicines， and the core targets included ABCG2， IFN-γ， EDN1， IL-6， and SOD1. The regulatory pathways were mainly related to transporters， microvascular function regulation， inflammation， and oxidative stress.