Decomposition of energy and free energy changes by following the flow of work along reaction path

被引:18
|
作者
Haas, Kevin [1 ]
Chu, Jhih-Wei [1 ]
机构
[1] Univ Calif Berkeley, Dept Chem Engn, Berkeley, CA 94720 USA
来源
JOURNAL OF CHEMICAL PHYSICS | 2009年 / 131卷 / 14期
基金
美国国家科学基金会;
关键词
MOLECULAR-DYNAMICS; RANDOM-WALK; COMPLEX; KINETICS; BINDING; THERMODYNAMICS; POLARIZATION; OPTIMIZATION; PERTURBATION; METHOTREXATE;
D O I
10.1063/1.3243080
中图分类号
O64 [物理化学(理论化学)、化学物理学];
学科分类号
070304 ; 081704 ;
摘要
To extract mechanistic information of activated processes, we propose to decompose potential energy and free energy differences between configurations into contributions from individual atoms, functional groups, or residues. Decomposition is achieved by calculating the mechanical work associated with the displacements and forces of each atom along a path that connects two states, i.e., following the flow of work. Specifically, we focus on decomposing energy or free energy differences along representative pathways such as minimum energy paths (MEPs) and minimum free energy paths (MFEPs), and a numerical metric is developed to quantify the required accuracy of the reaction path. A statistical mechanical analysis of energy decomposition is also presented to illustrate the generality of this approach. Decomposition along MEP and MFEP is demonstrated on two test cases to illustrate the ability to derive quantitative mechanistic information for different types of activated processes. First, the MEP of alanine dipeptide isomerization in vacuum and the MFEP of isomerization in explicit water is studied. Our analysis shows that carbonyl oxygen and amide hydrogen contribute to most of the energetic cost for isomerization and that explicit water solvation modulates the free energy landscape primarily through hydrogen bonding with these atoms. The second test case concerns the formation of tetrahedral intermediate during a transesterification reaction. Decomposition analysis shows that water molecules not only have strong stabilization effects on the tetrahedral intermediate but also constitute a sizable potential energy barrier due to their significant structural rearrangement during the reaction. We expect that the proposed method can be generally applied to develop mechanistic understanding of catalytic and biocatalytic processes and provide useful insight for strategies of molecular engineering. (C) 2009 American Institute of Physics. [doi:10.1063/1.3243080]
引用
收藏
页数:11
相关论文
共 50 条