队列关联分析
探讨MEHP暴露与骨质疏松风险之间的关联。
基于NHANES 2013-2018数据，纳入2846名参与者，通过多因素logistic回归和限制性立方样条分析MEHP与OP的关联。
关于丁香通
公司信息
个人用户
企业机构
无忧采购 轻松科研
- 移动端
探索增塑剂与骨质疏松症的关联及机制：一项队列和网络毒理学研究。
Exploring plasticisers-osteoporosis links and mechanisms: a cohort and network toxicology study.
实验完整度
高包含NHANES队列关联分析、网络毒理学与分子模拟，以及BMSC细胞功能验证等多个独立证据层级。
主要模型
NHANES 2013-2018队列（2846名参与者） 小鼠骨髓基质细胞（BMSCs）
重点核对
MEHP浓度检测：HPLC-ESI-MS/MS，LLOD为0.8 ng/mL，低于LLOD的值按LLOD/√2估算 骨质疏松定义：BMD低于年轻成人参考均值2.5个标准差以上 MEHP处理浓度：500、600、700 ng/mL，处理时间1、3、5、7天 siRNA敲低CTSD和VCP，转染浓度50 nM 分子动力学模拟：100 ns，温度310 K，压力1.0 bar
摘要
背景：增塑剂广泛存在于日常生活中，已被认为与骨质疏松症（OP）相关，但确切机制仍不清楚。方法：本研究基于NHANES数据，采用多因素logistic回归分析邻苯二甲酸单（2-乙基己基）酯（MEHP）与OP的关联。网络毒理学确定了MEHP诱导OP的关键靶点和通路。分子对接和动力学模拟验证了MEHP-靶点相互作用的稳定性。进一步在小鼠骨髓基质细胞（BMSCs）中评估了MEHP对成骨分化的影响。结果：所有logistic回归模型均证实MEHP水平与OP呈显著正相关。网络毒理学分析确定CTSD、SOAT1和VCP为关键靶点，凋亡通路为MEHP诱导OP的关键机制。分子模拟表明MEHP与这些靶点的结合稳定。细胞实验显示，MEHP通过下调CTSD和VCP显著抑制BMSC成骨，而SOAT1相关性较弱。结论：MEHP暴露与OP风险正相关，CTSD、VCP和凋亡通路可能在损害BMSC成骨中发挥关键作用。
实验结论
提炼研究问题、关键发现与证据，快速把握文章的核心贡献。
- 研究问题
- 增塑剂代谢物MEHP暴露是否与骨质疏松风险相关，其潜在机制如何？
- 核心机制
- MEHP可能通过下调CTSD和VCP表达，激活凋亡通路，从而抑制BMSC成骨分化，参与骨质疏松的发生。
- 主要证据
- NHANES队列分析显示MEHP水平与OP风险正相关（最高四分位OR=4.228, 95% CI 2.112-8.461）；分子模拟显示MEHP与CTSD、VCP、SOAT1稳定结合；细胞实验证实MEHP抑制BMSC成骨，siRNA敲低CTSD和VCP后成骨能力下降。
- 研究意义
- 为理解环境暴露因素MEHP对骨代谢的影响提供新证据，可能为骨质疏松的预防和治疗提供潜在靶点。
研究路径
按研究推进顺序梳理实验设计、验证步骤与关键观察。
网络毒理学分析
识别MEHP诱导OP的潜在靶点和通路。
从GEO下载GSE56815数据集，鉴定差异表达基因；从PubChem获取MEHP结构，利用SwissTargetPrediction、SEA和SuperPred预测靶点；取交集后进行GO和KEGG富集分析。
关键靶点筛选
利用机器学习和PPI网络筛选MEHP诱导OP的关键靶点。
构建PPI网络，使用cytoHubba四种算法筛选hub基因；采用LASSO、SVM-RFE和Boruta三种机器学习算法进一步筛选，取交集得到关键靶点CTSD、VCP和SOAT1。
分子对接与动力学模拟
验证MEHP与关键靶点蛋白的结合稳定性和亲和力。
从PDB获取靶点蛋白结构，使用AutoDock Vina进行分子对接；利用GROMACS进行100 ns分子动力学模拟，分析RMSD、Rg、SASA和结合自由能。
细胞实验验证
验证MEHP对BMSC成骨分化的影响及关键靶点的作用。
分离培养小鼠BMSCs，用MEHP处理，通过CCK8检测增殖，ALP染色、qRT-PCR和Western blot检测成骨标志物及CTSD、VCP表达；siRNA敲低CTSD和VCP验证其对成骨的影响。
研究方法
按研究目的归类文中使用的方法，便于定位所需技术。
产品清单
| 实验环节 | 名称 | 品牌 | 货号 |
|---|---|---|---|
| 分子对接 | AutoDock Vina | -- | -- |
| 分子动力学模拟 | GROMACS 2020.6 | -- | -- |
| Sobtop | -- | -- | |
| gmx_MMPBSA | -- | -- | |
| 分子对接 | Open Babel GUI | -- | -- |
| PyMOL | -- | -- | |
| 细胞培养 | 高糖DMEM培养基 | -- | -- |
| 胎牛血清 | -- | -- | |
| 青霉素-链霉素 | -- | -- | |
| I型胶原酶 | -- | -- | |
| 胰蛋白酶 | -- | -- | |
| 磷酸盐缓冲液 | -- | -- | |
| RNA提取与PCR | TRIzol试剂 | -- | -- |
| ReverTra Ace逆转录酶 | Toyobo | -- | |
| SYBR FAST qPCR预混液 | Tsingke | -- | |
| ALP染色 | 碱性磷酸酶染色试剂盒 | Sigma-Aldrich | SCR004 |
| Western blot | RIPA裂解液 | -- | -- |
| PVDF膜 | -- | -- | |
| BMP2一抗 | Bioss | bs-0514R | |
| RUNX2一抗 | Bioss | BSM-52672R | |
| GAPDH一抗 | Proteintech | 10494-1-AP | |
| COL1一抗 | Proteintech | 14695-1-AP | |
| VCP一抗 | Proteintech | 82463-1-RR | |
| CTSD一抗 | Proteintech | 21327-1-AP | |
| HRP标记的二抗 | -- | -- | |
| Super ECL Plus化学发光HRP底物 | Solarbio | -- | |
| siRNA转染 | siRNA | Sangon | -- |
| Lipofectamine 3000转染试剂 | Invitrogen | -- | |
| 细胞处理 | 邻苯二甲酸单(2-乙基己基)酯 | -- | -- |
关键环节
汇总复现实验时建议重点确认的条件及原文阅读提示。
| 环节 | 核对要点 |
|---|---|
| 队列分析 | NHANES周期（2013-2018）、样本筛选标准（年龄≥20岁，排除MEHP或OP数据缺失）、MEHP测量方法（HPLC-ESI-MS/MS，LLOD=0.8 ng/mL）、OP定义（WHO标准）、协变量调整（年龄、性别、种族等） 阅读提示：参见Methods 2.1.1和2.1.2 |
| 网络毒理学 | GSE56815数据集、差异表达基因筛选标准（P<0.05且|log2FC|>0.1）、MEHP靶点预测数据库（SwissTargetPrediction、SEA、SuperPred）、GO和KEGG富集分析方法 阅读提示：参见Methods 2.2 |
| 关键靶点筛选 | PPI网络构建参数（最小交互分数0.15）、cytoHubba算法（Betweenness、Degree、MCC、Stress）、机器学习算法（LASSO、SVM-RFE、Boruta）及交叉验证设置 阅读提示：参见Methods 2.3 |
| 分子对接与动力学模拟 | 对接软件（AutoDock Vina）、蛋白质结构来源（UniProt/PDB）、对接格子参数、MD模拟力场（AMBER99SB-ILDN/GAFF）、水模型（TIP3P）、温度和压力、模拟时长（100 ns） 阅读提示：参见Methods 2.5和2.6 |
| 细胞实验 | BMSCs来源（2日龄新生小鼠）、MEHP处理浓度（500、600、700 ng/mL）及时间（1-7天）、成骨诱导时间（7天）、siRNA序列和转染浓度（50 nM）、检测指标（ALP、COL1、RUNX2、BMP2、CTSD、VCP） 阅读提示：参见Methods 2.7-2.9 |