Zhang Shi, Giese Timothy J, Lee Tai-Sung, York Darrin M
Laboratory for Biomolecular Simulation Research, Institute for Quantitative Biomedicine and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey 08854, United States.
J Chem Theory Comput. 2024 May 14;20(9):3935-3953. doi: 10.1021/acs.jctc.4c00251. Epub 2024 Apr 26.
An alchemical enhanced sampling (ACES) method has recently been introduced to facilitate importance sampling in free energy simulations. The method achieves enhanced sampling from Hamiltonian replica exchange within a dual topology framework while utilizing new smoothstep softcore potentials. A common sampling problem encountered in lead optimization is the functionalization of aromatic rings that exhibit distinct conformational preferences when interacting with the protein. It is difficult to converge the distribution of ring conformations due to the long time scale of ring flipping events; however, the ACES method addresses this issue by modeling the and ring conformations within a dual topology. ACES thereby samples the conformer distributions by alchemically tunneling between states, as opposed to traversing a physical pathway with a high rotational barrier. We demonstrate the use of ACES to overcome conformational sampling issues involving ring flipping in ML300-derived noncovalent inhibitors of SARS-CoV-2 Main Protease (M). The demonstrations explore how the use of replica exchange and the choice of softcore selection affects the convergence of the ring conformation distributions. Furthermore, we examine how the accuracy of the calculated free energies is affected by the degree of phase space overlap (PSO) between adjacent states (i.e., between neighboring λ-windows) and the Hamiltonian replica exchange acceptance ratios. Both of these factors are sensitive to the spacing between the intermediate states. We introduce a new method for choosing a schedule of λ values. The method analyzes short "burn-in" simulations to construct a 2D map of the nonlocal PSO. The schedule is obtained by optimizing an alchemical pathway on the 2D map that equalizes the PSO between the λ intervals. The optimized phase space overlap λ-spacing method (Opt-PSO) leads to more numerous end-to-end single passes and round trips due to the correlation between PSO and Hamiltonian replica exchange acceptance ratios. The improved exchange statistics enhance the efficiency of ACES method. The method has been implemented into the FE-ToolKit software package, which is freely available.
最近引入了一种炼金术增强采样(ACES)方法,以促进自由能模拟中的重要性采样。该方法在双拓扑框架内通过哈密顿量副本交换实现增强采样,同时利用新的平滑步长软核势。在先导优化中遇到的一个常见采样问题是芳香环的官能化,当与蛋白质相互作用时,芳香环表现出不同的构象偏好。由于环翻转事件的时间尺度较长,难以使环构象的分布收敛;然而,ACES方法通过在双拓扑内对环构象进行建模来解决这个问题。ACES通过在状态之间进行炼金术隧穿来采样构象异构体分布,而不是通过具有高旋转势垒的物理途径。我们展示了使用ACES来克服涉及ML300衍生的SARS-CoV-2主蛋白酶(M)非共价抑制剂中环翻转的构象采样问题。这些演示探讨了副本交换的使用和软核选择的选择如何影响环构象分布的收敛。此外,我们研究了计算的自由能的准确性如何受到相邻状态(即相邻λ窗口之间)之间的相空间重叠(PSO)程度和哈密顿量副本交换接受率的影响。这两个因素都对中间状态之间的间距敏感。我们引入了一种选择λ值调度的新方法。该方法分析短的“预烧”模拟以构建非局部PSO的二维图。通过在二维图上优化炼金术途径来获得调度,该途径使λ区间之间的PSO相等。由于PSO与哈密顿量副本交换接受率之间的相关性,优化的相空间重叠λ间距方法(Opt-PSO)导致更多的端到端单次通过和往返。改进的交换统计提高了ACES方法的效率。该方法已被实现到FE-ToolKit软件包中,该软件包可免费获得。