将AMBER 14SB力场添加到随机定位CpHMD方法中
João G N Sequeira1, Adrian E Roitberg2, Miguel Machuqueiro1
1BioISI─Instituto de Biossistemas e Ciências Integrativas, Departamento de Química e Bioquímica, Faculdade de Ciências, Universidade de Lisboa, Lisboa 1749-016, Portugal.
Journal of chemical theory and computation
|June 16, 2025
概括
恒定pH分子动力学 (CpHMD) 模拟现在支持AMBER 14SB力场,增强质子化景观研究. 这一进步使主要蛋白质力场的直接比较成为可能,以改进生物分子模拟.
科学领域:
- 计算化学是一种计算化学.
- 生物分子模拟的模拟.
- 蛋白质动力学 蛋白质动力学
背景情况:
- 精确的分子动力学 (MD) 模拟需要结合pH效应来捕捉构造性,能量和质子化状态.
- 恒定pH分子动力学 (CpHMD) 方法,特别是随机定位CpHMD (st-CpHMD),是用于pH依赖模拟的先进技术.
- st-CpHMD目前支持GROMOS 54A7和CHARMM 36m的力场.
研究的目的:
- 在GROMACS软件包中将随机定位CpHMD (st-CpHMD) 兼容性扩展到AMBER 14SB力场.
- 为了验证 st-CpHMD 模拟的 AMBER 14SB 修改后的力场.
- 为了能够在CpHMD中直接比较GROMOS 54A7,CHARMM 36m和AMBER 14SB力场之间的性能.
主要方法:
- 实施了st-CpHMD对GROMACS.的AMBER 14SB力场的支持.
- 引入了对AMBER 14SB原子部分电荷进行轻微修改,以实现主链中和.
- 使用lyszyme和Staphylococcal核酶蛋白进行基准测试,比较pKa预测和计算成本.
主要成果:
- 珀14SB实现实现了与其他力场相比的pKa预测准确度,并且具有较低的根平均平方误差 (RMSE).
- 珀14SB模拟显示的计算成本低于CHARMM 36m,但高于GROMOS 54A7.7的计算成本.
- 确定了需要进一步改进方法的具有挑战性的案例,并建议加快PB/MC步骤以提高计算速度.
结论:
- 开发了第一个与GROMOS 54A7,CHARMM 36m和AMBER 14SB力场兼容的CpHMD方法.
- 实现了主要蛋白质力场的直接,定量性能比较,用于CpHMD模拟.
- 强调需要进一步优化,特别是在PB/MC步骤中,以提高计算效率.
相关概念视频
Calculating pH Changes in a Buffer Solution
54.1K
A buffer can prevent a sudden drop or increase in the pH of a solution after the addition of a strong acid or base up to its buffering capacity; however, such addition of a strong acid or base does result in the slight pH change of the solution. The small pH change can be calculated by determining the resulting change in the concentration of buffer components, i.e., a weak acid and its conjugate base or vice versa. The concentrations obtained using these stoichiometric calculations can be used...
54.1K
Titration of a Weak Base with a Strong Acid
5.8K
The titration curve of a weak base like ammonia with a strong acid like hydrochloric acid is the mirror image of the titration curve of a weak acid with a strong base.
Using the ICE table and substituting the Kb value, we calculate the initial pH of 50 mL of 0.1 M ammonia to be 11.11. Addition of 25 mL of 0.1 M hydrochloric acid to this solution of ammonia results in a buffer with an equal concentration of ammonia and ammonium ions. The pH of this buffer can be calculated by substituting these...
Using the ICE table and substituting the Kb value, we calculate the initial pH of 50 mL of 0.1 M ammonia to be 11.11. Addition of 25 mL of 0.1 M hydrochloric acid to this solution of ammonia results in a buffer with an equal concentration of ammonia and ammonium ions. The pH of this buffer can be calculated by substituting these...
5.8K
Titration of a Weak Acid with a Strong Base
2.4K
In titrating a weak acid with a strong base, different calculation methods are applied at various stages. Initially, the pH of a weak acid like acetic acid is calculated using its dissociation constant (Ka) and an ICE table. Upon addition of a strong base such as sodium hydroxide, a buffer forms, and its pH is determined using the Henderson-Hasselbalch equation. As more base is added and the titration reaches the halfway point, the pH becomes equal to the pKa of the acid, indicating equal...
2.4K
Acid-Base Titration Curves
129.7K
A titration curve is a plot of some solution property versus the amount of added titrant. For acid-base titrations, solution pH is a useful property to monitor because it varies predictably with the solution composition and, therefore, may be used to monitor the titration’s progress and detect its endpoint. Acid-base titration can be performed with a strong acid and a strong base, a strong acid and a weak base, or a strong base and a weak acid.
For a titration carried out for 25.00 mL of...
For a titration carried out for 25.00 mL of...
129.7K
Titration Calculations: Weak Acid - Strong Base
45.3K
Calculating pH for Titration Solutions: Weak Acid/Strong Base
For the titration of 25.00 mL of 0.100 M CH3CO2H with 0.100 M NaOH, the reaction can be represented as:
For the titration of 25.00 mL of 0.100 M CH3CO2H with 0.100 M NaOH, the reaction can be represented as:
45.3K
Titration of Polyprotic Base with a Strong Acid
884
The titration of a polyprotic base such as sodium carbonate with a strong acid such as hydrochloric acid results in two equivalence points on the titration curve. At the first equivalence point, the carbonate ions in the base are completely converted to bicarbonate ions. The second equivalence point corresponds to the complete conversion of bicarbonate ions to carbonic acid, which dissociates into carbon dioxide and water. The region before the first equivalence point corresponds to the...
884


