含铁硫簇蛋白结构优化问题咨询

模小助

整理一个 ITP/TOP 相关问题，重点看参数和连接方式是否合理。

我研究的蛋白为活性SAM酶，含4Fe-4S 和2Fe-2S，做分子模拟时一直不知道怎么获得拓扑结构，看到 moleculetype ]，将之和铁硫簇部分的[ moleculetype ]都引入主top里，把铁硫簇和配位残基的原子电荷替换为前面计算的，把铁硫簇部分的参数替换为簇模型中铁硫簇部分的

主要想确认：这种拓扑/力场处理是否合理，有没有更稳妥的设置方式？

admin

金属中心是这类体系最需要谨慎的部分。不能只根据“有一个 Zn/Cu/Fe 就加一个普通离子类型”，要先判断金属是游离离子、配位中心还是共价/半共价结构；不同情况适用的模型完全不同。

建议先确认实验配位残基、配位数和几何构型，再决定非键模型、bonded model 或专门参数化方案。模拟过程中最好监控金属-配位原子距离和配位数，避免“跑得稳定但化学结构已经错了”。