The SAAP force field. A simple approach to a new all-atom protein force field by using single amino acid potential (SAAP) functions in various solvents

被引:16
作者
Iwaoka, M [1 ]
Tomoda, S [1 ]
机构
[1] Univ Tokyo, Grad Sch Arts & Sci, Dept Life Sci, Meguro Ku, Tokyo 1538902, Japan
关键词
force field; protein folding; Monte Carlo simulation; glycine; alanine;
D O I
10.1002/jcc.10259
中图分类号
O6 [化学];
学科分类号
0703 ;
摘要
A simple strategy to compose a new all-atom protein force field (named as the SAAP force field), which utilizes the single amino acid potential (SAAP) functions obtained in various solvents by ab initio molecular orbital calculation applying the isodensity polarizable continuum model (IPCM), is presented. We considered that the total energy function of a protein force field (E-TOTAL) is divided into three components; a single amino acid potential term (E-SAAP), an interamino acid nonbonded interaction term (E-INTER), and a miscellaneous term (E-OTHERS), which is ignored (or considered to be constant) at the current version of the force field. The E INTER term consists of electrostatic interactions (E-ES') and van der Waals interactions (E-LJ'). Despite simplicity, the SAAP force field implicitly involves the correlation among individual terms of the Lifson's potential function within a single amino acid unit and can treat solvent effects unambiguously by choosing the SAAP function in an appropriate solvent and the dielectric constant (D) of medium. Application of the SAAP force field to the Monte Carlo simulation of For-Ala(2)-NH2 in vacuo reasonably reproduced the results of the extensive conformational search by ab initio molecular orbital calculation. In addition, the preliminary Monte Carlo simulations for For-Gly(10)-NH2 and For-Ala(10)-NH2 showed reversible transitions from the extended to, the pseudosecondary structures in water (D = 78.39) as well as in ether (D = 4.335). The result suggested that the new approach is efficient for fast modeling of protein structures in various environments. Decomposition analysis of the total energy function (E-TOTAL) by using the SAAP force field suggested that conformational propensities of single amino acids (i.e., the E-SAAP term) may play definitive roles on the topologies of protein secondary structures. (C) 2003 Wiley Periodicals, Inc.
引用
收藏
页码:1192 / 1200
页数:9
相关论文
共 61 条
[1]  
Akiyama S, 2000, NAT STRUCT BIOL, V7, P514
[2]   A combined discrete/continuum solvation model: Application to glycine [J].
Bandyopadhyay, P ;
Gordon, MS .
JOURNAL OF CHEMICAL PHYSICS, 2000, 113 (03) :1104-1109
[3]   Crystallographic evidence for C-alpha-H center dot center dot center dot O=C hydrogen bonds in a collagen triple helix [J].
Bella, J ;
Berman, HM .
JOURNAL OF MOLECULAR BIOLOGY, 1996, 264 (04) :734-742
[4]  
Berendsen H., 1981, INTERMOLECULAR FORCE, V331, P331, DOI [DOI 10.1007/978-94-015-7658-1_21, 10.1007/978-94-015-7658, DOI 10.1007/978-94-015-7658]
[5]  
Binder K., 1987, Applications of the Monte Carlo method in statistical physics
[6]   POTENTIAL FUNCTIONS AND CONFORMATIONS IN CYCLOALKANES [J].
BIXON, M ;
LIFSON, S .
TETRAHEDRON, 1967, 23 (02) :769-&
[7]   C-H•••π-interactions in proteins [J].
Brandl, M ;
Weiss, MS ;
Jabs, A ;
Sühnel, J ;
Hilgenfeld, R .
JOURNAL OF MOLECULAR BIOLOGY, 2001, 307 (01) :357-377
[8]   CHARMM - A PROGRAM FOR MACROMOLECULAR ENERGY, MINIMIZATION, AND DYNAMICS CALCULATIONS [J].
BROOKS, BR ;
BRUCCOLERI, RE ;
OLAFSON, BD ;
STATES, DJ ;
SWAMINATHAN, S ;
KARPLUS, M .
JOURNAL OF COMPUTATIONAL CHEMISTRY, 1983, 4 (02) :187-217
[9]   FUNNELS, PATHWAYS, AND THE ENERGY LANDSCAPE OF PROTEIN-FOLDING - A SYNTHESIS [J].
BRYNGELSON, JD ;
ONUCHIC, JN ;
SOCCI, ND ;
WOLYNES, PG .
PROTEINS-STRUCTURE FUNCTION AND BIOINFORMATICS, 1995, 21 (03) :167-195
[10]   Recent advances in molecular dynamics simulation towards the realistic representation of biomolecules in solution [J].
Cheatham, TE ;
Brooks, BR .
THEORETICAL CHEMISTRY ACCOUNTS, 1998, 99 (05) :279-288