靶向铜死亡相关基因治疗类风湿关节炎生物信息学分析及干预中药的预测

杨会军，郭 响，李伟青，苏小军，田雪梅，王海东，金芳梅*

·数据挖掘与循证医学·

靶向铜死亡相关基因治疗类风湿关节炎生物信息学分析及干预中药的预测

杨会军1，郭 响2，李伟青1，苏小军1，田雪梅1，王海东1，金芳梅1*

1. 甘肃省中医院 风湿骨病中心，甘肃 兰州 730050 2. 浙江大学医学院附属第一医院，浙江 杭州 310003

探索铜死亡相关基因在类风湿关节炎（rheumatoid arthritis，RA）中的表达及其作用机制，并筛选通过靶向铜死亡基因治疗RA的潜在中药。通过检索GEO数据库获得RA芯片数据，分析铜死亡基因在RA中表达水平；依据铜死亡基因表达情况构建高表达亚型和低表达亚型，并进行差异分析，筛选出与RA相关的铜死亡基因。按不同分组对基因表达矩阵进行免疫浸润分析，并分析免疫浸润细胞与铜死亡基因之间的相关性。利用Kobas在线工具对差异基因进行京都基因与基因组百科全书（Kyoto encyclopedia of genes and genomes，KEGG）通路富集分析。通过Herb数据库检索靶向铜死亡基因治疗RA的潜在中药及其有效成分，使用Autodock vina软件模拟药物与铜死亡靶蛋白结合活性。通过收集临床样本验证铜死亡相关基因在RA中表达水平。、、、、、和7个铜死亡相关基因在RA中差异表达。免疫浸润分析结果显示，在RA中活化的记忆CD4+T细胞、静息自然杀伤（natural killer，NK）细胞和中性粒细胞与铜死亡相关。富集分析结果显示，铜死亡主要与趋化因子信号通路、Toll样受体信号通路、白细胞介素-17（interleukin-17，IL-17）信号通路、肿瘤坏死因子信号通路、核因子-κB信号通路、p53信号通路和过氧化物酶体增殖物激活受体信号通路相关。虎杖叶可能是靶向铜死亡基因治疗RA的潜在中药。分子对接结果显示，虎杖叶有效成分均有与铜死亡靶蛋白结合的可能。、、、在RA患者外周血单个核细胞中mRNA表达水平升高。铜死亡相关基因表达水平与RA病程进展相关，并预测虎杖叶是潜在靶向铜死亡基因治疗RA的药物，为RA的临床诊疗及药物开发提供了新的方向和理论基础。

类风湿关节炎；铜死亡相关基因；免疫浸润；生物信息学；虎杖叶；二氢硫辛酰胺脱氢酶；谷氨酰胺酶；铁氧还原蛋白1；二氢硫辛酸转乙酰化酶

类风湿关节炎（rheumatoid arthritis，RA）是临床上常见的一种炎症性疾病，症状主要表现为持续性疼痛、压痛和关节破坏[1]。祖国医学认为RA属于“痹证”范畴，主要因气血痹阻，经脉不通而痛；气行不畅而痹阻经络；血瘀气阻而化热[2]。目前我国RA患病率约0.42%，全球患病率约0.27%[3-4]。RA作为慢性疾病已成为全世界关注的公共卫生问题之一。RA的标准治疗手段主要用于控制炎症和关节疼痛，但常用的非甾体抗炎药、糖皮质激素、改善病情的抗风湿药物和生物制剂长期服用会引起不良反应[4-5]。因此，寻找安全、有效的治疗药物或治疗策略对RA临床治疗具有重要意义。RA的病理过程受T细胞、B细胞、单核细胞、巨噬细胞、中性粒细胞和肥大细胞等免疫细胞与成纤维样滑膜细胞之间的串扰影响[6]。免疫细胞的激活会分泌多种炎症因子，进而影响软骨细胞、破骨细胞和成纤维细胞导致软骨损伤和骨骼破坏[5-6]。干预这些免疫细胞或免疫因子作为RA治疗相关研究的新方法和新靶点具有很大的潜力，相关的治疗方案仍在研究中[6]。

细胞死亡与RA的疾病发展密切相关。如成纤维样滑膜细胞死亡减少导致滑膜细胞过度增殖，从而引发滑膜炎[7]。自身免疫性T细胞和B细胞的细胞死亡失衡也会增加自身免疫反应和炎症反应[8]。随着细胞死亡研究领域的快速发展，针对RA细胞死亡的相关临床治疗策略也取得了一定的进展[7-8]。铜死亡是一种新发现的细胞死亡方式，其机制与线粒体呼吸相关[7]。铜离子通过与三羧酸循环过程中的脂酰化蛋白质结合，使脂酰化蛋白聚集导致铁硫簇蛋白丢失，进而诱发蛋白质毒性应激并造成细胞死亡[7,9]。相关研究结果显示RA患者血液中铜含量升高，且与炎症标志物呈正相关[10]。铜死亡相关基因丙酮酸脱氢酶E1 α1亚基（pyruvate dehydrogenase E1 subunit alpha 1，）可能与转录激活因子3（signal transducer and activator of transcription 3，STAT3）协同调节RA中的巨噬细胞炎症反应，且铜死亡另一相关基因丙酮酸脱氢酶E1亚基β（pyruvate dehydrogenase E1 subunit beta，）可能与成纤维样滑膜细胞的异常增殖相关[11]。由此可见，阐明铜死亡在RA发病过程中的作用机制并探寻相关治疗药物具有重要的研究意义。

传统中医药是中华文化的瑰宝，治疗RA已有近千年历史。临床上中医主要通过祛风除湿、散瘀活血、通络止痛来消除痹症，代表性的复方有蠲痹汤、独活寄生汤、乌头汤和桂枝芍药知母汤[2,12]。有实验表明传统中药及其有效成分可以通过促进细胞凋亡、抑制成纤维样滑膜细胞异常增殖和抗炎等作用缓解RA症状[2,13]。然而中药及其活性成分能否通过靶向铜死亡治疗RA尚未见相关研究报道。本研究旨在通过生物信息学方法阐明铜死亡相关基因在RA中的表达及作用机制，筛选靶向铜死亡相关基因的中药及其活性成分。为明确铜死亡在RA中的调控机制并开发用于防治RA的靶向调节药物提供理论基础。

1 方法

1.1 铜死亡相关基因集获取

以“cuproptosis”和“铜死亡”为关键词分别检索PubMed数据库（https://pubmed.ncbi.nlm.nih.gov/）和中国知网数据库（https://www.cnki.net/）截至2022年9月30日收录的与铜死亡相关文献，构建铜死亡相关基因集。共获得10个与铜死亡密切相关的基因：细胞周期依赖性激酶抑制基因（cell cycle-dependent kinase inhibitory gene，）、铁氧还原蛋白1（ferridoxin 1，）、二氢硫辛酰胺脱氢酶（dihydrolipoamide dehydrogenase，）、二氢硫辛酸转乙酰化酶（dihydrolipoamide acetyltransferase，）、硫辛酸合成酶（lipoic acid synthetase，）、谷氨酰胺酶（glutaminase，）、脂酰转移酶1（lipoyltransferase 1，）、金属调节转录因子1（metal regulatory transcription factor 1，）用于后续分析[9]。

1.2 RA滑膜基因表达数据集获取

以“rheumatoid arthritis”为关键词检索由美国国立生物技术信息中心（National Center of Biotechnology Information，NCBI）创建并维护的基因表达数据库（Gene Expression Omnibus，GEO，https://www.ncbi.nlm.nih.gov/geo/），获取RA基因表达数据集。数据集筛选要求：（1）数据集包含RA组和正常对照组，且每组样本量不少于3；（2）RA组未经药物或其他手段干预治疗；（3）数据集信息完整无缺失。

1.3 RA样本分型

基于铜死亡相关基因在GEO数据集RA样本中的表达水平，使用Consensus Cluster Plus R语言包对所有RA样本进行表达一致性聚类分析，将样本分为铜死亡基因表达模式不同的亚型。

1.4 基因表达差异分析

通过RStudio 2022.02.1软件的“Limma”R语言包，对检索得到的基因表达数据集根据不同分组进行基因表达差异分析。以差异倍数（fold change，FC）绝对值大于2及矫正后＜0.05为差异表达基因（differentially expressed gene，DEG）筛选标准。差异分析结果通过微生信（https://www. bioinformatics.com.cn/）在线网页工具对结果进行可视化处理。

1.5 免疫浸润分析

基于基因表达数据矩阵使用CIBERSORT去卷积法评估每个样本中22种免疫细胞占比。比较不同分组间（正常对照组RA组、RA1组RA2组）22种免疫细胞丰度差异。

1.6 富集分析

将需要分析的基因简称导入在线网页工具Kobas（http://bioinfo.org/kobas）进行京都基因与基因组百科全书（Kyoto Encyclopedia of Genes and Genomes，KEGG）通路富集分析。分析物种设置为“智人（homo sapiens）”，以＜0.05为筛选标准。

1.7 潜在治疗中药预测

对在RA病理过程中发挥重要作用的铜死亡相关基因，通过本草图鉴数据库（HERB，http://herb.ac.cn/）检索其相应中药。针对涵盖靶点较多的中药，通过检索相关文献，确定其有效成分。通过SwissADME在线工具（http://www.swissadme. ch）依据里宾斯基五规则（Lipinski’s rule of five，RO5）筛选出成药性较好的有效成分，并使用ProTox-II在线工具对有效成分进行安全性评估[14]。借助AutoDockTools 1.5.6软件模拟中药有效成分结合铜死亡相关基因，预测其治疗RA的潜在可能。

1.8 临床样本收集

本研究获得甘肃省中医院伦理委员会批准（批准文号2020-031-01），所有研究对象均签署知情同意书。本研究于甘肃省中医院共收集5例健康志愿者外周血样本（健康对照组）及5例临床诊断为RA患者（RA组）外周血样本，所有外周血样本通过提取试剂分离出单核细胞，放入−80 ℃冰箱保存。

纳入标准：（1）RA组患者经临床医生诊断符合2010年美国风湿病学会联合欧洲抗风湿病联盟提出的RA诊断标准；（2）RA组患者处于RA疾病活动期；（3）年龄18～80岁；（4）性别不限；（5）签署知情同意书者。排除标准：（1）患有精神疾病或其他原因无法配合者；（2）合并其他急慢性疾病者；（3）近期有生育要求或正处于妊娠期或哺乳期的女性；（4）入组前4周内使用过药物干预治疗者。

1.9 实时荧光定量-PCR（real-time fluorescence quantification-PCR，RT-qPCR）分析

通过Trizol法提取外周血样本单个核细胞总RNA，检测浓度后取1000 ng RNA为模板逆转录为cDNA（Takara，RR047A）。引物序列如表1所示，使用SYBR green染料（诺唯赞，批号Q141-02）按20 μL体系进行40个循环的扩增。t值以磷酸甘油醛脱氢酶（reduced glyceraldehyde-phosphate dehydrogenase，GAPDH）为内参按2−∆∆Ct法分析。

表1 引物序列

1.10 统计分析

2 结果

2.1 铜死亡相关基因在RA中表达情况

通过检索GEO数据库得到物种为人的RA滑膜基因表达数据集GSE89408。数据集基于芯片GPL11154，包括对照组样本28例，RA样本152例。基于GSE89408数据集，10个铜死亡相关基因表达差异分析结果显示（图1-A），与正常对照组相比，RA组关节滑膜、、、、、、、和表达上调且差异具有统计学意义（＜0.05）。

2.2 铜死亡相关基因表达分型

基于10个铜死亡相关基因在RA样本关节滑膜表达情况，将152例RA滑膜样本进行一致性聚类分析（图1-B），分为RA1和RA2 2种亚型。主成分分析（principal component analysis，PCA）结果显示（图1-C），RA1与RA2 2种亚型间存在差异，可用于进一步分析。铜死亡相关基因在RA1和RA2 2种亚型表达差异分析结果显示（图1-D），与RA1组相比，RA2组、、、、、、、和表达下调且差异具有统计学意义（＜0.05）。

2.3 差异分析结果

基于数据集GSE89408将RA组与正常对照组相比，根据DEGs筛选条件，共得到10 558个RA滑膜DEGs（图2-A）。铜死亡2个亚型RA1与RA2分组间差异分析共得到652个DEGs（图2-B）。对2次差异分析结果进行交集分析，得到501个共同DEGs。

2.4 通路富集分析结果

501个共同DEGs通路分析结果（图2-C）显示，趋化因子（chemokine）信号通路、Toll样受体（Toll-like receptor）信号通路、白细胞介素-17（interleukin-17，IL-17）信号通路、肿瘤坏死因子（tumor necrosis factor，TNF）信号通路、核因子-κB（nuclear factor-κB，NF-κB）信号通路、p53信号通路和过氧化物酶体增殖物激活受体（peroxisomeproliferators-activated receptor，PPAR）信号通路与铜死亡相关，在RA病理过程中发挥重要作用。

2.5 免疫浸润分析结果

免疫浸润分析结果显示，与正常对照组相比（图3-A），RA组关节滑膜样本中记忆B细胞（B cell memory，＜0.001）、浆细胞（plasma cells，＜0.05）、未活化的记忆CD4+T细胞（T cells CD4 memory resting，＜0.000 1）、活化的记忆CD4+T细胞（T cells CD4 memory activated，＜0.000 1）、滤泡辅助性T细胞（T cells follicular helper，＜0.05）、调节性T细胞（T cells regulatory，＜0.000 1）、γδT细胞（T cells gamma delta，＜0.05）、静息自然杀伤（natural killer，NK）细胞（NK cells resting，＜0.01）、活化NK细胞（NK cells activated，＜0.000 1）、单核细胞（monocytes，＜0.05）、M1型巨噬细胞（macrophages M1，＜0.000 1）、静息树突细胞（dendritic cells resting，＜0.05）、静息肥大细胞（mast cells resting，＜0.01）和中性粒细胞（neutrophils，＜0.000 1）丰度存在差异。

RA1组与RA2组免疫浸润分析结果（图3-B）显示，两组间初始CD4+T细胞（T cells CD4 naive，＜0.01）、活化的记忆CD4+T细胞（＜0.05）、静息NK细胞（＜0.05）、M2型巨噬细胞（macrophages M2，＜0.001）、嗜酸性粒细胞（eosinophils，＜0.05）和中性粒细胞（＜0.01）丰度存在差异。

2.6 铜死亡相关基因诊断价值分析

结合在正常对照组与RA组及RA1与RA2组间差异表达的铜死亡相关基因，共有8个共同铜死亡相关基因（、、、、、、和）表达存在差异。借助微生信在线网页工具构建基因受试者工作特征（receiver operator characteristic，ROC）曲线，获得基因对RA诊断的灵敏度和特异度，评估基因对RA诊断价值。结果（图4）显示，（AUC＝0.862）、（AUC＝0.826）、（AUC＝0.880）、（AUC＝0.901）、LIAS（AUC＝0.852）、LIPT1（AUC＝0.851）和PDHB（AUC＝0.857）特异性较高。此外，PDHA1（AUC＝0.658）特异性较低。

A-铜死亡相关基因在正常对照组和RA组表达水平，与正常对照组比较：*P＜0.05 ***P＜0.001 ****P＜0.000 1；B-基于铜死亡相关基因表达水平对RA样本聚类分析；C-RA不同亚型组PCA分析结果；D-RA不同亚型组铜死亡相关基因表达水平，与RA1组比较：****P＜0.000 1

2.7 相关性分析结果

7个共同铜死亡相关基因与3个共同差异免疫细胞（活化的记忆CD4+T细胞、静息NK细胞和中性粒细胞）相关性分析结果如图5所示，、、、、、和之间表达相关性均为正相关（＞0.05）。活化的记忆CD4+T细胞与7个铜死亡相关基因表达均呈正相关。与之相反，静息NK细胞与7个铜死亡相关基因表达的相关性系数均小于0。中性粒细胞与7个铜死亡相关基因表达相关性较弱。

2.8 靶向铜死亡基因治疗RA的潜在中药及其有效成分确定

基于Herb数据库检索能靶向7个共同铜死亡相关基因的潜在中药，其中未能找到能靶向、和基因的中药，潜在靶向、、和基因的治疗RA的中药如表2所示。通过检索相关文献确定中药的有效成分，无相关文献则通过Herb数据库补充。共获得7种潜在中药（草莓车轴草、虎杖叶、白罂粟、丽春花果实、水团花、粟壳和罂粟嫩苗），并确定41种有效成分。

A-正常对照组与RA组差异分析火山图；B-RA1组与RA2组差异分析火山图，蓝色代表表达下调基因，红色代表表达上调基因，灰色代表无统计学差异基因；C-RA不同亚型组差异基因与RA相关基因韦恩图；D-交集基因相关通路

2.9 分子对接模拟药物分子与蛋白结合

通过SwissADME在线工具共筛选出符合RO5原则的药物成分13种（表3），安全性评估结果显示咖啡酸、库明吉啶、大黄素、考迈斯托醇、阿魏酸、没食子酸和白藜芦醇毒性较低（毒性等级＞3，出现毒性反应数≤1）。以甲氨蝶呤为阳性对照。将以上7个安全性较高的成分与7个共同铜死亡相关基因进行分子对接，药物分子配体与蛋白受体结合能如图6-A所示。结合能均小于0，表明药物分子与蛋白受体有结合可能。尽管吗啡与DLAT结合能最低，但鉴于吗啡属于危险药物，成为RA治疗药物可能性较低。本研究展示其他结合能较低的构象。DLD与大黄素结合构象、库明吉啶与FDX1结合构象如图6-B、C所示。

2.10 两组受试者基线资料

最终正常对照组和RA组各纳入5例受试者。两组受试者基线资料如表4所示。

2.11 铜死亡相关基因在RA患者中mRNA表达水平

铜死亡相关基因、、、在健康受试者和RA患者外周血单个核细胞中mRNA表达水平如图7所示。与健康对照组相比，RA患者外周血单个核细胞中、、、的mRNA表达水平明显升高且差异具有统计差异（＜0.05、0.01）。

3 讨论

铜死亡是新发现的一种细胞死亡方式，可能会影响RA中的各种细胞，从而涉及炎症、血管生成和骨破坏的发展过程[11]。本研究通过下载GEO芯片数据，对RA滑膜基因表达数据分析结果显示，大部分铜死亡相关基因在RA患者样本中表达上调，这表明铜死亡与RA的疾病进展相关。基于铜死亡基因在RA样本中的表达水平，将RA样本分为铜死亡基因高表达亚型组和低表达亚型组。对正常对照组与RA组、铜死亡基因高表达组与低表达组分别进行免疫浸润分析，确定了RA中3个与铜死亡相关的免疫细胞，分别为活化的记忆CD4+T细胞、静息NK细胞和中性粒细胞。大量记忆CD4+T细胞、NK细胞存在于RA患者外周血、关节滑膜与滑液中，这些细胞能够高效地产生促炎细胞因子进而导致关节炎症，同时促进骨破坏[17-18]。同样在RA环境中，中性粒细胞表达多种趋化因子受体，并释放大量炎症介质[19]。因此，靶向抑制这些免疫细胞，降低促炎细胞因子水平是一种有前景的RA治疗策略[5]。铜死亡基因与免疫细胞可能在RA免疫浸润方面发挥重要作用，能否通过调节免疫细胞死亡缓解RA症状，有待进一步研究。

A-正常对照组vs RA组 B-RA1组vs RA2组 *P＜0.05 **P＜0.01 ***P＜0.001 ****P＜0.000 1

基于差异分析共得到501个DEGs，KEGG通路富集分析结果显示，RA中铜死亡可能与趋化因子信号通路、Toll样受体信号通路、IL-17信号通路、肿瘤坏死因子信号通路、NF-κB信号通路、p53信号通路和PPAR信号通路相关。NF-κB信号通路是细胞因子与炎症反应之间的桥梁，铜元素可以通过激活NF-κB信号通路促进炎症和线粒体自噬障碍[20]。在长期铜暴露或铜诱导的条件下，p53表达增加且肿瘤坏死因子信号通路被激活，进一步促进细胞凋亡[21-22]。Liu等[23]研究显示，铜元素可能通过激活PPAR信号通路调节脂质代谢。PPAR也已被证实对CD4+T细胞的存活、激活和分化具有调节作用[24]。在RA中，各种趋化因子通过靶向中性粒细胞、成骨细胞、NK细胞、T细胞、软骨细胞和巨噬细胞等促进血管生成，加重关节滑膜炎症，并使关节长期处于慢性炎症环境中[25]。此外，IL-17通过与其受体结合刺激TNF-α、IL-1β、IL-6和趋化因子产生，以及NF-κB活化[25]。而NF-κB通路被认为是RA炎症反应的主要途径之一[26]。由此可见，通过干预铜死亡相关通路促进促炎细胞群铜死亡可以有效缓解RA炎症。

图4 铜死亡相关基因ROC曲线

图5 相关性分析结果

表2 靶向铜死亡基因治疗RA的潜在中药

表3 潜在治疗RA药物成分化学性质与毒性参数

Y-出现毒性反应 N-无毒性反应

Y-toxic reaction occurs N-no toxic reaction occurs

A-潜在治疗药物与靶点蛋白分子对接结合能 B-DLD与大黄素结合构象 C-库明吉啶与FDX1结合构象

表4 两组受试者基线资料

HC-健康对照组 RA-RA组 与HC组比较：*P＜0.05 **P＜0.01

铜死亡基因在RA芯片数据中表达分析结果显示，、、、、、和7个基因在2次差异分析中均为差异表达基因，提示其可能是靶向铜死亡干预RA的重要基因。目前尚无在RA中的研究报道，主要促进糖酵解和脂肪酸氧化，并改变氨基酸代谢参与炎症反应[27]。主要与线粒体功能、氧化应激和炎症有关，过表达可以增加调节性T细胞（regulatory T cell，Treg）的数量并减少了T细胞浸润，抑制NF-κB活性，减轻慢性炎症反应并保护线粒体功能[11,28-29]。在RA成纤维样滑膜细胞中表达上调，抑制表达水平可以抑制RA成纤维样滑膜细胞增殖并抑制关节炎小鼠的炎症反应[30]。主要调节线粒体呼吸、乳酸和丙酮酸的产生，并与RA患者的滑膜细胞增殖和血管生成相关[11,31]。也能通过影响线粒体功能参与RA的发展[11]。已被确定为RA的易感基因，可能与帕金森病蛋白7共同调节Treg细胞并保持其功能完整性；表达下调可能导致RA成纤维样滑膜细胞异常增殖[11]。ROC曲线分析7个重要基因，AUC值均高于0.5。以上结果表明这些基因可作为RA的诊断标志物和潜在治疗靶点。

通过检索相关数据库获得干预RA铜死亡基因的中药及其有效成分，得到7种中药及其41个有效成分。其中罂粟虽然具有敛肺涩肠、固肾镇痛等功效，但其主要成分吗啡和可待因等生物碱具有成瘾性，已被国家列入麻醉精神类药品目录严格管制[32]。本研究不探讨此类药物及相关成分。《本草拾遗》和《本草推陈》记载传统中药虎杖叶具有平肝潜阳、祛痰熄风等功效，临床上主要应用于治疗风湿关节疼痛、消炎和消肿等[15]。相关研究表明，虎杖中有效成分具有抗炎和抗痛风作用，在血管生成方面具有双重作用，且长期高剂量应用虎杖会产生毒性，造成肝脏损伤，但毒性反应具有可逆性[33]。因此，虎杖叶有成为靶向铜死亡缓解RA炎症症状治疗药物的潜在可能，但虎杖叶中发挥主要治疗作用的有效成分尚不明确，需要进一步筛选鉴定。根据RO5成药原则及毒性预测结果，筛选出咖啡酸、库明吉啶、大黄素、考迈斯托醇、阿魏酸、没食子酸和白藜芦醇7个药物成分。分子对接结果显示，这些药物成分均有与RA铜死亡靶点结合的可能。白藜芦醇已被应用于临床治疗RA，且治疗效果较好。在RA中，白藜芦醇对p53通路具有调控作用并促进成纤维样滑膜细胞、巨噬细胞和浸润的炎性细胞凋亡[34]。通过构建佐剂性关节炎大鼠模型，发现咖啡酸和阿魏酸可减轻关节炎大鼠炎症，缓解爪子水肿和炎症细胞浸润，保护关节组织免受血管翳形成以及软骨和骨破坏[35-36]。大黄素通过抑制NF-κB通路在胶原诱导关节炎小鼠中发挥抗炎作用[37]。没食子酸可降低促炎细胞因子水平并诱导成纤维样滑膜细胞凋亡[38]。在骨关节炎中，考迈斯托醇可以抑制IL-1β诱导的炎症[39]。由此可见，虎杖叶可以作为靶向铜死亡治疗RA的潜在治疗药物。

4 结论

靶向铜死亡基因促进免疫细胞死亡，从而缓解RA炎症，是一种很有前景的治疗策略。本研究分析了铜死亡基因在RA中表达水平差异，且可能与记忆CD4+T细胞、NK细胞和中性粒细胞免疫浸润相关。同时本研究预测虎杖叶有效成分可能靶向铜死亡基因缓解RA炎症。本研究的分析结果为靶向铜死亡治疗RA的相关研究提供了新的方向和理论基础，为临床RA治疗提供新的选择药物。但研究仍存在一些问题，如未能检索到草莓车轴草有效成分及其药理作用的相关文献，有待进一步深入研究。

利益冲突 所有作者均声明不存在利益冲突

[1] Padyukov L. Genetics of rheumatoid arthritis [J]., 2022, 44(1): 47-62.

[2] 叶敏, 万军, 王李俊, 等. 中药通过调节肠道菌群治疗类风湿关节炎的研究进展[J]. 中草药, 2022, 53(22): 7279-7284.

[3] Finckh A, Gilbert B, Hodkinson B,. Global epidemiology of rheumatoid arthritis [J]., 2022, 18(10): 591-602.

[4] 耿研, 谢希, 王昱, 等. 类风湿关节炎诊疗规范 [J]. 中华内科杂志, 2022, 61(1): 51-59.

[5] Mueller A L, Payandeh Z, Mohammadkhani N,. Recent advances in understanding the pathogenesis of rheumatoid arthritis: New treatment strategies [J]., 2021, 10(11): 3017.

[6] Zhao J N, Guo S C, Schrodi S J,. Molecular and cellular heterogeneity in rheumatoid arthritis: Mechanisms and clinical implications [J]., 2021, 12: 790122.

[7] Cobine P A, Brady D C. Cuproptosis: Cellular and molecular mechanisms underlying copper-induced cell death [J]., 2022, 82(10): 1786-1787.

[8] Zhao J N, Jiang P, Guo S C,. Apoptosis, autophagy, NETosis, necroptosis, and pyroptosis mediated programmed cell death as targets for innovative therapy in rheumatoid arthritis [J]., 2021, 12: 809806.

[9] Tsvetkov P, Coy S, Petrova B,. Copper induces cell death by targeting lipoylated TCA cycle proteins [J]., 2022, 375(6586): 1254-1261.

[10] Wang H J, Zhang R R, Shen J,. Circulating level of blood iron and copper associated with inflammation and disease activity of rheumatoid arthritis [J]., 2023, 201(1): 90-97.

[11] Zhao J N, Guo S C, Schrodi S J,. Cuproptosis and cuproptosis-related genes in rheumatoid arthritis: Implication, prospects, and perspectives [J]., 2022, 13: 930278.

[12] 王帆帆, 刘健, 方妍妍, 等. 中医药改善类风湿关节炎血瘀证的研究进展 [J]. 风湿病与关节炎, 2022, 11(9): 74-77.

[13] Wang Y, Chen S, Du K,. Traditional herbal medicine: Therapeutic potential in rheumatoid arthritis [J]., 2021, 279: 114368.

[14] Lipinski C A, Lombardo F, Dominy B W,. Experimental and computational approaches to estimate solubility and permeability in drug discovery and development settings [J]., 1997, 23(1/2/3): 3-25.

[15] 王欣, 覃瑶, 孙建彬, 等. 虎杖叶的化学成分、药理活性、临床应用及质量控制研究进展 [J]. 亚太传统医药, 2019, 15(10): 196-200.

[16] 郭跃伟, 黄伟晖, 陈雯婷, 等. 水团花化学成分的研究 [J]. 中国现代中药, 2012, 14(3): 15-19.

[17] Yang Y Y, Day J, Souza-Fonseca Guimaraes F,. Natural killer cells in inflammatory autoimmune diseases [J]., 2021, 10(2): e1250.

[18] Luque-Campos N, Contreras-López R A, Jose Paredes-Martínez M,. Mesenchymal stem cells improve rheumatoid arthritis progression by controlling memory T cell response [J]., 2019, 10: 798.

[19] Cecchi I, Arias de La Rosa I, Menegatti E,. Neutrophils: Novel key players in rheumatoid arthritis. Current and future therapeutic targets [J]., 2018, 17(11): 1138-1149.

[20] Zhou Q, Zhang Y, Lu L,. Copper induces microglia-mediated neuroinflammation through ROS/NF-κB pathway and mitophagy disorder [J]., 2022, 168: 113369.

[21] Li Y L, Chen H M, Liao J Z,. Long-term copper exposure promotes apoptosis and autophagy by inducing oxidative stress in pig testis [J]., 2021, 28(39): 55140-55153.

[22] Liu H, Guo H R, Jian Z J,. Copper induces oxidative stress and apoptosis in the mouse liver [J]., 2020, 2020: 1-20.

[23] Liu L, Sui X Y, Li F C. Effect of dietary copper addition on lipid metabolism in rabbits [J]., 2017, 61(1): 1348866.

[24] Le Menn G, Neels J. Regulation of immune cell function by PPARs and the connection with metabolic and neurodegenerative diseases [J]., 2018, 19(6): 1575.

[25] Fang Q H, Zhou C, Nandakumar K S. Molecular and cellular pathways contributing to joint damage in rheumatoid arthritis [J]., 2020, 2020: 1-20.

[26] Ilchovska D D, Barrow D M. An Overview of the NF-κB mechanism of pathophysiology in rheumatoid arthritis, investigation of the NF-kB ligand RANKL and related nutritional interventions [J]., 2021, 20(2): 102741.

[27] Zhang Z Y, Ma Y R, Guo X L,. FDX1 can impact the prognosis and mediate the metabolism of lung adenocarcinoma [J]., 2021, 12: 749134.

[28] Tian S, Nakamura J, Hiller S,. New insights into immunomodulation via overexpressing lipoic acid synthase as a therapeutic potential to reduce atherosclerosis [J]., 2020, 133/134: 106777.

[29] Zhao Y, Xu G, Li H,. Overexpression of endogenous lipoic acid synthase attenuates pulmonary fibrosis induced by crystalline silica in mice [J]., 2020, 323: 57-66.

[30] Takahashi S, Saegusa J, Sendo S,. Glutaminase 1 plays a key role in the cell growth of fibroblast-like synoviocytes in rheumatoid arthritis [J]., 2017, 19(1): 1-10.

[31] Chang X T, Wei C. Glycolysis and rheumatoid arthritis [J]., 2011, 14(3): 217-222.

[32] 宋小兰, 王丽萍, 赵宁娟, 等. 基于气相色谱-质谱和网络药理学探讨紫斑罂粟壳挥发油镇咳祛痰平喘的作用机制 [J]. 分析测试技术与仪器, 2022, 28(3): 343-352.

[33] 梁春晓, 王珊珊, 陈淑静, 等. 虎杖化学成分及药理活性研究进展 [J]. 中草药, 2022, 53(4): 1264-1276.

[34] Sheng S Y, Wang X Y, Liu X,. The role of resveratrol on rheumatoid arthritis: From bench to bedside [J]., 2022, 13: 829677.

[35] Fikry E M, Gad A M, Eid A H,. Caffeic acid and ellagic acid ameliorate adjuvant-induced arthritis in rats via targeting inflammatory signals, chitinase-3-like protein-1 and angiogenesis [J]., 2019, 110: 878-886.

[36] Zhu L J, Zhang Z S, Xia N N,. Anti-arthritic activity of ferulic acid in complete Freund’s adjuvant (CFA)-induced arthritis in rats: JAK2 inhibition [J]., 2020, 28(2): 463-473.

[37] Hwang J K, Noh E M, Moon S J,. Emodin suppresses inflammatory responses and joint destruction in collagen-induced arthritic mice [J]., 2013, 52(9): 1583-1591.

[38] Yoon C H, Chung S J, Lee S W,. Gallic acid, a natural polyphenolic acid, induces apoptosis and inhibits proinflammatory gene expressions in rheumatoid arthritis fibroblast-like synoviocytes [J]., 2013, 80(3): 274-279.

[39] You J S, Cho I A, Kang K R,. Coumestrol counteracts interleukin-1beta-Induced catabolic effects by suppressing inflammation in primary rat chondrocytes [J]., 2017, 40(1): 79-91.

Bioinformatics analysis of targeting cuproptosis-related genes therapy for rheumatoid arthritis and prediction of intervention traditional Chinese medicine

YANG Hui-jun1, GUO Xiang2, LI Wei-qing1, SU Xiao-jun1, TIAN Xue-mei1, WANG Hai-dong1, JIN Fang-mei1

1. Osteopathy and Rheumatology Center, Lanzhou 730050, China 2. The First Affiliated Hospital, Zhejiang University School of Medicine, Hangzhou 310003, China

To explore the expression of cuproptosis-related genes in rheumatoid arthritis (RA) and mechanism, and to screen the potential traditional Chinese medicine in the treatment RA by targeting cuproptosis-related genes.The RA microarray data were obtained by searching the GEO database, and the expression levels of cuproptosis-related genes in RA were analyzed. The high expression subtypes and low expression subtypes were constructed according to the expression of cuproptosis-related genes, and differential analysis was performed to screen out the cuproptosis-related genes related to RA. Immune infiltration analysis was performed on the gene expression matrix by different groups and the correlation between immune infiltrating cells and cuproptosis-related genes was analyzed. The Kyoto encyclopedia of genes and genomes (KEGG) pathway enrichment analysis was performed on differential expressed genes by the Kobas online tool. The Herb database was used to search for potential traditional Chinese medicines and their active ingredients in the treatment of RA by targeting cuproptosis-related gene, and the Autodock vina software was used to simulate the possible binding of drugs and cuproptosis-related target proteins. The expression levels of cuproptosis-related genes in RA were verified by clinical samples.A total of seven cuproptosis-related genes (,,,,,, and) were differentially expressed in RA. The result of immune infiltration analysis showed that activated memory CD4+T cells, resting NK cells, and neutrophils were associated with cuproptosis in RA. The results of pathway enrichment analysis showed that cuproptosis was mainly associated with the chemokine signaling pathway, Toll-like receptor signaling pathway, interleukin-17 (IL-17) signaling pathway, tumor necrosis factor signaling pathway, nuclear factor-κB signaling pathway, p53 signaling pathway, and peroxisome proliferators-activated receptors signaling pathways. Leaf of Huzhang () may be a potential drug in the treatment of RA by targeting cuproptosis-related gene, and the molecular docking results showed that the active components in leaf ofmay bind to cuproptosis-related proteins. The mRNA expression levels of,,andin peripheral blood mononuclear cells of RA patients increased.The expression level of cuproptosis -related genes is related to the progression of RA disease, and predicted thatleaf is a potential drug in treatment of RA by targeting cuproptosis-related genes, which provides a new direction and theoretical basis for the clinical diagnosis and treatment of RA and drug development.

rheumatoid arthritis; cuproptosis-related genes; immune infiltraton; bioinformatics; leaf of; DLD; GLS; FDX1; DLAT

R285

A

0253 – 2670(2023)13 - 4253 - 13

10.7501/j.issn.0253-2670.2023.13.018

2022-12-23

国家中医药管理局高水平中医药重点学科建设项目实施方案（国中医药人教函[2022]226号）；国家自然科学基金项目（82260880）；国家重点研发计划（2018YFC1705202，2018YFC1705203）；国家中医药管理局行业专项（20150700107）；甘肃省中医药项目（GZKP-2021-12）；兰州市科技指导性计划项目（2022-ZD-73）；王海东甘肃省名中医传承工作室建设项目（甘卫中医函[2022]50号）

杨会军（1985—），男，医学博士，主治医师，研究方向为风湿骨病发病机制及中医特色疗法和中药干预研究。E-mail: yhj_16999@163.com

通信作者：金芳梅，女，主治医师，研究方向为风湿骨病发病机制及中医特色疗法和中药干预研究。E-mail: 13909406533@163.com

[责任编辑 潘明佳]