摘要
葫芦巴碱作为咖啡来源的生物碱，具有抗氧化和线粒体保护作用，但其在卵巢衰老中的具体机制尚不明确。本研究通过整合网络药理学、蛋白互作分析、分子对接、100纳秒分子动力学模拟及单细胞转录组数据，筛选出57个潜在靶点，其中MMP9、JAK2、PARP1、HDAC1和CYP3A4被确定为核心靶点，涉及细胞外基质重塑、基因组与表观基因组维护、内分泌/外源物代谢等关键通路。结构模拟与单细胞表达谱共同支持PARP1和MMP9为葫芦巴碱调控卵巢衰老的核心介质，JAK2和HDAC1为中间候选靶点，CYP3A4可能介导系统代谢效应。本研究为阐明咖啡衍生葫芦巴碱在卵巢衰老中的作用提供了理论框架和实验验证方向。
1 引言
卵巢功能维持对生育、内分泌健康及长期代谢风险至关重要。卵巢衰老涉及细胞外基质重塑、慢性炎症、氧化应激、基因组不稳定及激素代谢紊乱等多重过程，但膳食生物活性分子对此过程的影响尚未明确。咖啡作为广泛摄入的生物活性物质来源，其成分葫芦巴碱近年被发现具有线粒体保护和氧化还原调节作用，且在人体内吸收迅速、半衰期较长（约13–14小时），可能通过多靶点途径影响卵巢衰老进程。
2 材料与方法
2.1 靶点预测
通过SwissTargetPrediction数据库基于葫芦巴碱的SMILES结构预测人源蛋白靶点。
2.2 卵巢衰老相关基因筛选
从OMIM和GeneCards数据库以“卵巢衰老”为关键词检索相关基因，筛选评分≥5的基因。
2.3 蛋白互作网络构建
利用STRING数据库构建蛋白互作网络，通过Cytoscape的CytoHubba插件计算节点中心性，综合度值、介数、最大团中心性（MCC）和最大邻域成分（MNC）四种指标筛选核心靶点。
2.4 通路富集分析
使用DAVID和Reactome数据库进行基因本体（GO）和京都基因与基因组百科全书（KEGG）通路富集分析，显著性阈值设为p<0.05。
2.5 单细胞转录组分析
基于GEO数据集GSE255690中青年（18–28岁）与中年（36–39岁）卵巢样本的单细胞RNA测序数据，通过Seurat流程分析靶点基因的细胞类型特异性表达及年龄相关变化。
2.6 分子对接与动力学模拟
采用AutoDock Vina进行分子对接，结合GROMACS进行100纳秒分子动力学模拟，评估结合自由能（MM/PBSA ΔGbind）、均方根偏差（RMSD）、回转半径（Rg）等参数。
3 结果
3.1 核心靶点筛选
从84个葫芦巴碱预测靶点与4186个卵巢衰老相关基因中交集得到57个共同靶点。PPI网络分析显示MMP9、JAK2、PARP1、HDAC1和CYP3A4在四种中心性指标中均位列前十，被定义为核心枢纽基因。
3.2 功能富集分析
KEGG富集通路包括类固醇激素生物合成、PI3K-Akt信号通路、松弛素信号通路等；GO分析显示靶点主要参与UV-A响应、蛋白去乙酰化、染色质重塑等生物过程。Reactome模块进一步提示靶点聚集于免疫信号、ECM组织、DNA修复及线粒体代谢等通路。
3.3 单细胞表达谱
中年卵巢组织中，HDAC1在内皮、颗粒细胞等多类细胞中表达显著下调；PARP1在血管基质细胞中下调而在颗粒细胞中轻微上调；MMP9主要在单核细胞中表达升高；CYP3A4在卵巢组织中几乎不表达，提示其可能通过系统代谢间接作用。
3.4 分子对接与动力学验证
葫芦巴碱与各靶点的对接能介于-4.6至-6.5 kcal/mol，其中与CYP3A4（-6.5 kcal/mol）、PARP1（-6.2 kcal/mol）和MMP9（-6.1 kcal/mol）结合较强。分子动力学模拟显示所有复合物结构稳定，结合自由能计算进一步支持PARP1（-10.38 kcal/mol）和MMP9（-12.28 kcal/mol）为高亲和力靶点。
4 讨论
本研究首次通过计算生物学手段系统阐释葫芦巴碱可能通过多靶点调控卵巢衰老的分子网络。其中PARP1和MMP9在结构模拟和转录组数据中均表现出核心地位，分别关联DNA修复效能和ECM动态平衡，而JAK2和HDAC1可能通过炎症信号和表观遗传修饰参与调节。未来需通过体外实验验证葫芦巴碱对这些靶点的直接作用，并探索其与咖啡其他成分（如绿原酸）的协同效应。本研究为开发以膳食干预延缓生殖衰老的策略提供了理论基础。
