黑皮质素‑4受体（melanocortin‑4 receptor, MC4R）基因在脊椎动物的食欲调控、体重及能量平衡中发挥核心作用，并对绵羊的生长、脂肪沉积、饲料效率和繁殖具有重要影响。然而，该基因中非同义单核苷酸多态性（non‑synonymous single nucleotide polymorphisms, nsSNPs）的功能后果仍不明确。本研究旨在采用整合的计算分析方法，评估绵羊MC4R基因中nsSNPs对结构与功能的潜在影响。研究人员从Ensembl和UniProt数据库中筛选出12个绵羊MC4R nsSNPs，并利用多种工具对其功能、稳定性及进化保守性进行系统预测。功能影响由SIFT、PolyPhen‑2、SNPs&GO和PhD‑SNP进行评估，蛋白稳定性由I‑Mutant、MUpro、CUPSAT和DynaMut分析，保守性分析采用ConSurf，结构建模使用SOPMA、I‑TASSER和GalaxyWEB，模型验证借助PROCHECK和ProSA，并进行了100 ns的全原子分子动力学模拟（WEBGRO）。结果显示，R220G、G98R、E100K和Y187C四个nsSNPs在多项预测中一致被判定为有害，且这些有害变异多位于高度保守区域。结构建模与分子动力学模拟表明，这些nsSNPs可能改变蛋白的稳定性、柔性和构象动态，从而影响野生型蛋白的结构与功能。上述发现提示这些nsSNPs可能是影响MC4R结构和功能的候选变异，但其生物学和实际育种意义应谨慎解读，因为结果仅基于计算预测。所识别的变异可作为进一步研究的靶点，但其在选择性育种中的应用仍需通过群体研究、基因型‑表型关联分析及实验手段加以验证。总体而言，本研究为优先筛选可能影响功能的MC4R变异提供了计算框架，而非直接支持其用于育种实践的证据。
研究背景方面，黑皮质素‑4受体（MC4R）基因编码G蛋白偶联受体，在脊椎动物能量平衡、食欲及体重调节中具有关键作用。在畜牧业中，MC4R遗传变异与绵羊等重要经济性状密切相关，包括生长速率、肉质、脂肪沉积、生殖激素水平及饲料效率。非同义单核苷酸多态性（nsSNPs）因可导致氨基酸序列变化，可能直接影响蛋白质结构与功能，因此在分子标记辅助选择中备受关注。尽管MC4R在猪和牛中的研究较多，但在绵羊中的结构与功能变异仍缺乏系统性探索。现有绵羊MC4R相关研究多集中于表型关联分析，缺乏基于蛋白质水平的结构与功能机制解析。随着基因组数据资源的丰富，整合多种计算工具预测nsSNPs对蛋白功能、稳定性及动态特性的影响，成为弥补功能注释不足的重要手段。在此背景下，研究人员开展本项研究，旨在通过系统的计算分析，识别可能显著影响绵羊MC4R结构与功能的nsSNPs，为后续实验验证与育种应用提供依据。
在研究方法上，研究人员首先从Ensembl数据库获取绵羊MC4R基因的转录本与蛋白序列，并从UniProt获得参考蛋白编号B2MV59。从初始鉴定的593个变异中筛选出位于编码区的12个nsSNPs作为研究对象。功能预测采用SIFT、PolyPhen‑2、SNPs&GO和PhD‑SNP四种工具；蛋白稳定性分析使用I‑Mutant、MUpro、CUPSAT和DynaMut；进化保守性评估采用ConSurf；二级结构预测由SOPMA完成；三级结构建模通过I‑TASSER构建并用GalaxyWEB优化，模型质量经PROCHECK和ProSA‑web验证；蛋白‑蛋白互作分析依托STRING数据库；突变体构建与相互作用分析在UCSF Chimera软件中实现；分子动力学模拟在WEBGRO平台进行，总时长100 ns，并分析RMSD、RMSF、Rg、SASA及氢键等参数。
研究结果方面，3.1节的功能预测显示，12个nsSNPs中，R220G、G98R、E100K和Y187C在多个工具中一致被判定为有害或致病性变异，其余变异则多为中性或混合结果。3.2节的稳定性分析表明，多数nsSNPs对MC4R蛋白产生去稳定作用，其中M79T、R220G和G98R的去稳定效应最为显著，ΔΔG值低于‑1.0 kcal/mol。3.3节的ConSurf保守性分析发现，部分有害变异位于高度保守且埋藏的残基位置，提示其对维持蛋白结构完整性至关重要；Q307E位于高度保守且暴露的区域，可能参与功能活性。3.4节的二级结构预测显示，MC4R蛋白以α‑螺旋（43.67%）和无规卷曲（43.37%）为主，符合G蛋白偶联受体的典型结构特征。3.5节的三级结构建模与验证结果表明，优化后的模型立体化学质量良好，Ramachandran图中92.1%的残基位于最优选区，ProSA‑web的Z值为‑2.98，处于天然蛋白合理范围。3.6节的蛋白‑蛋白互作网络分析显示，MC4R与POMC（得分0.999）、NPY（0.973）、AGRP（0.951）等关键能量平衡调控因子存在强相互作用，构成高度互联的功能网络。3.7节的分子动力学模拟结果显示，有害nsSNPs如R220G和M79T引起更高的结构波动和柔性变化，导致蛋白稳定性下降，而野生型及其他部分变异则保持较稳定的构象。
讨论部分指出，虽然整合计算分析识别出多个潜在有害nsSNPs，但这些结果仍需实验验证才能应用于育种。MC4R在多种物种中与生产性状密切相关，本研究首次在绵羊中系统结合功能预测、结构建模及动力学模拟，提供了机制层面的见解。研究优势在于多工具交叉验证提高了预测可靠性，但仍存在纯计算分析的局限性。研究人员强调，这些nsSNPs应视为候选位点，未来需通过体外功能实验和体内表型关联研究确认其生物学意义。
结论部分翻译为：本研究对绵羊MC4R基因中12个非同义单核苷酸多态性进行了全面的计算分析。功能预测确定R220G、G98R、E100K和Y187C为最具破坏性的变异。稳定性分析显示，M79T、R220G和G98R的去稳定效应最强，ΔΔG值低于‑1.0 kcal/mol。保守性定位表明，多个有害变异位于高度保守区域，提示其潜在的结构与功能重要性。优化后的三维模型立体化学质量良好，蛋白‑蛋白互作分析确认MC4R在能量调控通路中的核心地位。分子动力学模拟进一步支持了关键nsSNPs的去稳定作用，表现为柔性增加与结构稳定性下降。这些结果突出了特定MC4R nsSNPs作为候选变异的潜力，值得在未来针对绵羊生长和代谢性状的选择性育种计划中深入研究。后续实验验证对于确认这些计算预测及其在畜牧遗传学中的实际应用至关重要。
