表面加速弦方法用于定位最小自由能量路径
Timothy J Giese1, Şölen Ekesan1, Erika McCarthy1
1Laboratory for Biomolecular Simulation Research, Institute for Quantitative Biomedicine and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey 08854, United States.
Journal of chemical theory and computation
|February 17, 2024
概括
我们开发了一种表面加速弦方法 (SASM),以更有效地优化反应路径. SASM使用来自多次代的聚合采样,融合路径比SMCV和MSMCV等现有方法快三倍.
科学领域:
- 计算化学计算化学
- 生物物理化学 生物物理化学
- 化学物理 化学物理
背景情况:
- 优化反应通路对于理解化学和生物过程至关重要.
- 量子力学/分子力学 (QM/MM) 方法对于路径采样是计算上昂贵的.
- 现有的方法,如集体变量中的字符串方法 (SMCV) 和修改的SMCV (MSMCV),在融合和路径表示方面存在局限性.
研究的目的:
- 引入和评估表面加速弦方法 (SASM),以优化有效的反应路径.
- 为了证明SASM能够加快融合并提高免费能源配置的准确性.
- 使用QM/MM应用程序,比较SASM与SMCV和MSMCV的性能.
主要方法:
- SASM利用来自当前和之前代的总量采样来加快路径收.
- 它将采样和路径表示图像数字脱.
- 雨潜在的位置被优化,以提高自由能量表面探索和准确性.
主要成果:
- SASM在平坦的自由能源地区进行了改进的勘探,并且在稀疏的离散化下显示出更好的配置质量.
- 对 ribozyme 甲基转移酶,Hammerhead ribozyme 和 B-DNA 分聚的比较研究表明,SASM 的收路径大约是 SMCV 和 MSMCV 的三倍.
- 所有方法 (SASM,SMCV,MSMCV) 都在免费可用的FE-ToolKit包中实现.
结论:
- 使用QM/MM方法,SASM可显著提高反应路径优化的效率.
- 该方法提高了免费能源配置计算的准确性和稳定性.
- SASM为复杂化学和生物反应的计算研究提供了有价值的工具.
相关概念视频
Energy Diagrams - II
4.6K
Energy diagrams are important to understand the dynamics of a system. The topology of an energy diagram helps illustrate the equilibrium points of the system.
The point in the energy diagram at which the system’s potential energy is the lowest is known as the local minima. The system tends to stay in this position indefinitely unless acted upon by a net force. The slope of the potential energy diagram at the local minima is zero, indicating that zero net force is acting on the system. The...
The point in the energy diagram at which the system’s potential energy is the lowest is known as the local minima. The system tends to stay in this position indefinitely unless acted upon by a net force. The slope of the potential energy diagram at the local minima is zero, indicating that zero net force is acting on the system. The...
4.6K
Arrhenius Plots
39.4K
The Arrhenius equation relates the activation energy and the rate constant, k, for chemical reactions. In the Arrhenius equation, k = Ae−Ea/RT, R is the ideal gas constant, which has a value of 8.314 J/mol·K, T is the temperature on the kelvin scale, Ea is the activation energy in J/mole, e is the constant 2.7183, and A is a constant called the frequency factor, which is related to the frequency of collisions and the orientation of the reacting molecules.
The Arrhenius equation can be used...
The Arrhenius equation can be used...
39.4K
Energy Diagrams, Transition States, and Intermediates
16.5K
Free-energy diagrams, or reaction coordinate diagrams, are graphs showing the energy changes that occur during a chemical reaction. The reaction coordinate represented on the horizontal axis shows how far the reaction has progressed structurally. Positions along the x-axis close to the reactants have structures resembling the reactants, while positions close to the products resemble the products. Peaks on the energy diagram represent stable structures with measurable lifetimes, while...
16.5K
Energy Diagrams - I
5.0K
The dynamics of a mechanical system can be easily understood by interpreting a potential energy diagram. Since energy is a scalar quantity, the interpretation of the dynamics of the system becomes even simpler.
Take the example of a skater on a parabolic ramp. The potential energy at different points along the ramp will be proportional to the height of the ramp, which varies quadratically with the horizontal position on the ramp. As the skater moves down the ramp from the highest position,...
Take the example of a skater on a parabolic ramp. The potential energy at different points along the ramp will be proportional to the height of the ramp, which varies quadratically with the horizontal position on the ramp. As the skater moves down the ramp from the highest position,...
5.0K
Maxam-Gilbert Sequencing
11.2K
In the same year as the discovery of the Sanger sequencing method, another group of scientists, Allan Maxam and Walter Gilbert, demonstrated their chemical-cleavage method for DNA sequencing. The Maxam-Gilbert method relies on using different chemicals that can cleave the DNA sequence at specific sites, the separation of resulting DNA fragments of variable size using electrophoresis, and deciphering the DNA sequence from the resulting gel bands.
Challenges of the Maxam-Gilbert Method
The...
Challenges of the Maxam-Gilbert Method
The...
11.2K
Energy to Drive Translocation
2.1K
Mitochondrial protein import is powered by two distinct energy sources: ATP hydrolysis and electrochemical potential across the inner membrane. Newly synthesized precursors are bound by cytosolic chaperones of the Hsp70 family, which guide them to the import receptors on the mitochondrial surface. Utilizing the energy of ATP hydrolysis, Hsp70 chaperones transfer these precursors to the TOM receptors on the mitochondrial outer membrane.
Generally, polypeptides are unfolded by two distinct...
Generally, polypeptides are unfolded by two distinct...
2.1K


