硫酸根离子促进Cnm–Cnm复合体稳定的多尺度分子动力学机制：口腔–全身关联的分子线索

《Frontiers in Microbiology》：Sulfate ions promote Cnm-mediated aggregation in Streptococcus mutans: molecular dynamics insights into the oral–systemic connection

　　

来自变形链球菌（Streptococcus mutans）的胶原结合蛋白Cnm与脑微出血相关，该病理状态常与缺血性卒中相关。近期实验证据提示，Cnm介导的相互作用可能驱动蛋白质–蛋白质聚集，并可能受离子环境影响。本研究试图通过多尺度分子动力学（molecular dynamics, MD）方法研究Cnm–Cnm缔合的分子起源决定因素，以阐明该问题。研究进行了全原子（all-atom, AA）模拟，以评估硫酸根离子（SO₄^2?）对复合物稳定性的影响；该离子选择基于在Cnm阳性S. mutans中实验观察到的含硫物种富集。为将分析扩展至更大尺寸和更长时间尺度，研究人员首先实施基于MARTINI力场的额外粗粒化（coarse-grained, CG）模拟，并针对原子水平参考进行验证，随后用于探索蛋白质–蛋白质缔合的热力学。结果表明，硫酸根离子显著增强Cnm–Cnm复合物的稳定性，表现为蛋白质间距离波动减小、分子间接触数增加以及结合自由能曲线中更深的极小值。这些发现提示，离子屏蔽可间接促进Cnm聚集，为可能与S. mutans致病性相关的聚集过程提供物理机制见解。

研究背景：口腔微生物组与宿主组织及免疫系统持续相互作用，口腔–全身连接逐渐受到关注。变形链球菌（Streptococcus mutans）经典上与龋病相关，但表达胶原结合蛋白Cnm的Cnm阳性S. mutans与系统性血管病理相关，并可能与脑微出血、缺血性卒中等脑血管/心血管疾病相关。Cnm可促进细菌黏附于富含胶原的组织，菌血症后可能参与血管表面相互作用；然而，Cnm相关表型与系统性血管效应之间的分子机制尚不完整。已有研究提示Cnm阳性菌株聚集表型增强，拉曼光谱显示含硫物种富集，但Cnm–Cnm自我关联及离子环境如何调控该复合物稳定性仍不完全清楚。因此，研究人员聚焦Cnm–Cnm缔合，整合MD模拟与既往显微和拉曼光谱观察，考察离子环境对蛋白质–蛋白质复合物稳定性的潜在调节作用。该研究不直接模拟细菌黏附血管组织或卒中发病机制，而是提供Cnm–Cnm相互作用的分子层面表征，以帮助理解更广泛的生物学表型。

关键技术方法概括：研究人员以Cnm二聚体晶体结构（RCSB PDB，PDB ID: pdb_00007lgr）为起点，构建含与不含SO₄^2?的水溶液体系。采用GROMACS 2022.3进行全原子MD（Amber03力场、SPC/E水、100 ns），并采用MARTINI 3力场进行粗粒化MD（1 μs）；以全原子轨迹为参考验证CG模型。轨迹分析包括均方根偏差（root-mean-square deviation, RMSD）、质心（center of mass, COM）距离、分子间接触、溶剂可及表面积（solvent-accessible surface area, SASA）与埋藏表面积（buried surface area, BSA）、平均力势（potential of mean force, PMF）及解离自由能（ΔG）估计。

研究结果：
3.1 显微观察与拉曼光谱证据：通过无蔗糖条件下Cnm阴性与Cnm阳性菌株的激光显微观察，发现Cnm阳性细胞聚集倾向更强、聚集体更大且更多；通过200–700 cm^?1区域拉曼光谱，发现Cnm阳性菌株中39条谱带中有25条可检测，包含与二硫键S–S伸缩相关的约527、511、547 cm^?1信号，以及多种氧化硫物种相关谱带（如253、635、651 cm^?1）；Cnm阴性菌株则显示糖原及多糖相关谱带（如356、440、485、574、600 cm^?1）等。结果表明两菌株在生化组成和适应策略上存在明显差异，Cnm阳性菌株表现为更强聚集倾向和含硫表面组分富集。

3.2 分子动力学模拟：由于完整体系过于复杂，研究人员选取硫酸根离子对Cnm–Cnm复合物稳定性的影响进行模拟分析。

3.2.1 粗粒化模型验证：以含硫酸根的全原子100 ns模拟为参考，比较两蛋白COM间径向分布函数及归一化分子间接触，发现CG与AA模型一致，平衡蛋白间距离约2.2–2.6 nm，界面相互作用保持稳定，支持CG模型用于后续长时间尺度分析。

3.2.2 硫酸根离子对复合物稳定性的影响：基于CG 1 μs模拟，含SO₄^2?体系中两蛋白平均COM距离为29.35 ?、最小18.35 ?，而无硫酸根体系平均距离为76.70 ?并出现部分解离；含硫酸根体系平均分子间接触数为79.52，无硫酸根体系为10.53；含硫酸根体系平均BSA为676.45 ?²，无硫酸根体系为396.75 ?²；含硫酸根体系RMSD在约200 ns后稳定于19.83 ?，无硫酸根体系则在30–60 ?间波动。PMF分析显示含硫酸根体系具有明确极小值，解离ΔG为10.34 kJ mol^?1，无硫酸根体系为5.73 kJ mol^?1；时间窗分析中，含硫酸根体系PMF主极小值约在23.5–32.1 ?，无硫酸根体系趋于66–126 ?。全原子轨迹重现了含硫酸根时更近最小距离和更多分子间接触的趋势。残基接触分析显示16对残基在至少70%轨迹帧中保持接触，其中6对超过90%：Tyr91–Glu178（98.13%）、Thr78–Asn177（96.63%）、Gln73–Thr206（94.57%）、Asn69–Phe172（91.66%）、Gln65–Thr158（90.54%）、Arg36–Glu178（90.36%）。结果表明，在所模拟条件下硫酸根离子与更短蛋白间距离、更多接触、更大埋藏界面、更小结构波动和更深自由能极小值相关；但仅比较含与不含硫酸根体系，不能确定该稳定化是否特异于SO₄^2?。

讨论部分总结：作者指出，模拟的硫酸根依赖稳定化应视为可能解释Cnm阳性菌株增强聚集表型的分子假设，需要受控离子条件下的直接实验验证；研究不提供硫酸根促进Cnm阳性S. mutans细胞聚集的直接实验证据，也不确立拉曼光谱检测的含硫物种来源于Cnm。模型局限于蛋白质水平，描述两个孤立Cnm分子在简化水环境中的缔合，不复现S. mutans细胞表面的空间组织和分子复杂性，也未涵盖其他表面组分、多物种微生物群落、宿主蛋白或口腔微环境复杂理化条件；且不直接模拟细菌黏附血管组织或脑血管病发病机制。未来工作将研究其他生理相关离子、Cnm变体序列和结构变异、拉曼光谱差异的具体分子组分、更大的多聚组装体，并整合胶原和内皮组分等宿主因素。

研究结论翻译：本研究整合显微观察和拉曼光谱实验观察与多尺度MD模拟，研究可能促进S. mutans中Cnm相关聚集的分子因素。实验结果显示Cnm阳性与Cnm阴性菌株之间的生化和结构差异，包括Cnm阳性表型更强的聚集倾向及与含硫物种相关的拉曼光谱特征差异。由于全细胞拉曼光谱反映细菌整体生化组成，这些光谱差异不能特异归因于Cnm蛋白，也可能反映所研究菌株之间更广泛的表型差异。在分子水平，模拟提供了可随这些实验观察一并考虑的机制解释。在所研究模拟条件下，全原子和粗粒化模拟一致显示含硫酸根时Cnm–Cnm复合物更稳定，表现为蛋白质间距离减小、分子间接触增加、埋藏表面积增大和自由能极小值更深。原子与粗粒化描述的一致性支持两种模拟分辨率下观察到硫酸根相关稳定化，但计算自由能曲线应解释为采样相互作用景观的定性描述，而非绝对热力学自由能。综上，实验和计算结果支持离子环境可能影响Cnm–Cnm缔合的情景。但本研究不提供硫酸根促进Cnm阳性S. mutans细胞聚集的直接实验证据，也不确立拉曼光谱检测的含硫物种来源于Cnm。因此，模拟中观察到的硫酸根依赖稳定化应视为可能解释增强聚集表型的分子假设，并需在受控离子条件下直接实验验证。此外，由于模拟仅比较含硫酸根与无硫酸根体系，不能确定观察到的稳定化是否特异于硫酸根，还是反映二价离子或局部静电环境改变的一般效应。