Feature Article Basic Ingredients of Free Energy Calculations: A Review

被引:308
作者
Christ, Clara D. [1 ]
Mark, Alan E. [2 ]
van Gunsteren, Wilfred F. [1 ]
机构
[1] ETH, Phys Chem Lab, Swiss Fed Inst Technol, CH-8093 Zurich, Switzerland
[2] Univ Queensland, Sch Mol & Microbial Sci, Inst Mol Biosci, St Lucia, Qld 4072, Australia
基金
瑞士国家科学基金会;
关键词
free energy; affinity; molecular simulation; sampling; MOLECULAR-DYNAMICS SIMULATIONS; REPLICA-EXCHANGE METHOD; MONTE-CARLO-SIMULATION; EQUATION-OF-STATE; FORCE-FIELD; CONFORMATIONAL DYNAMICS; COMPUTER-SIMULATION; PROTEIN-STRUCTURE; EFFICIENT; OPTIMIZATION;
D O I
10.1002/jcc.21450
中图分类号
O6 [化学];
学科分类号
0703 ;
摘要
Methods to compute free energy differences between different states of a molecular system are reviewed with the aim of identifying their basic ingredients and their utility when applied in practice to biomolecular systems. A free energy calculation is comprised of three basic components: (i) a suitable model or Hamiltonian, (ii) a sampling protocol with which one can generate a representative ensemble of molecular configurations, and (iii) an estimator of the free energy difference itself. Alternative sampling protocols can be distinguished according to whether one or more states are to be sampled. In cases where only a single state is considered, six alternative techniques could be distinguished: (i) changing the dynamics, (ii) deforming the energy surface, (iii) extending the dimensionality, (iv) perturbing the forces, (v) reducing the number of degrees of freedom, and (vi) multi-copy approaches. In cases where multiple states are to be sampled, the three primary techniques are staging, importance sampling, and adiabatic decoupling. Estimators of the free energy can be classified as global methods that either count the number of times a given state is sampled or use energy differences. Or. they can be classified as local methods that either make use of the force or are based on transition probabilities. Finally, this overview of the available techniques and how they can be best used in a practical context is aimed Lit helping the reader choose the most appropriate combination of approaches for the biomolecular system. Hamiltonian and free energy difference of interest. (C) 2009 Wiley Periodicals, Inc. J Comput Chem 31: 1569-1582, 2010
引用
收藏
页码:1569 / 1582
页数:14
相关论文
共 139 条
[31]   Adaptive biasing force method for scalar and vector free energy calculations [J].
Darve, Eric ;
Rodriguez-Gomez, David ;
Pohorille, Andrew .
JOURNAL OF CHEMICAL PHYSICS, 2008, 128 (14)
[32]  
Daura X, 1998, J COMPUT CHEM, V19, P535, DOI 10.1002/(SICI)1096-987X(19980415)19:5<535::AID-JCC6>3.0.CO
[33]  
2-N
[34]   HYBRID MONTE-CARLO [J].
DUANE, S ;
KENNEDY, AD ;
PENDLETON, BJ ;
ROWETH, D .
PHYSICS LETTERS B, 1987, 195 (02) :216-222
[35]  
Feenstra KA, 1999, J COMPUT CHEM, V20, P786, DOI 10.1002/(SICI)1096-987X(199906)20:8<786::AID-JCC5>3.0.CO
[36]  
2-B
[37]   OPTIMIZED MONTE-CARLO DATA-ANALYSIS [J].
FERRENBERG, AM ;
SWENDSEN, RH .
PHYSICAL REVIEW LETTERS, 1989, 63 (12) :1195-1198
[38]  
Foloppe N, 2000, J COMPUT CHEM, V21, P86, DOI 10.1002/(SICI)1096-987X(20000130)21:2<86::AID-JCC2>3.0.CO
[39]  
2-G
[40]  
Frenkel D., 1996, UNDERSTANDING MOL SI, DOI DOI 10.1016/B978-0-12-267351-1.X5000-7