TY - JOUR
T1 - Gay-Berne and electrostatic multipole based coarse-grain potential in implicit solvent
AU - Wu, Johnny
AU - Zhen, Xia
AU - Shen, Hujun
AU - Li, Guohui
AU - Ren, Pengyu
PY - 2011/10/21
Y1 - 2011/10/21
N2 - A general, transferable coarse-grain (CG) framework based on the Gay-Berne potential and electrostatic point multipole expansion is presented for polypeptide simulations. The solvent effect is described by the Generalized Kirkwood theory. The CG model is calibrated using the results of all-atom simulations of model compounds in solution. Instead of matching the overall effective forces produced by atomic models, the fundamental intermolecular forces such as electrostatic, repulsion-dispersion, and solvation are represented explicitly at a CG level. We demonstrate that the CG alanine dipeptide model is able to reproduce quantitatively the conformational energy of all-atom force fields in both gas and solution phases, including the electrostatic and solvation components. Replica exchange molecular dynamics and microsecond dynamic simulations of polyalanine of 5 and 12 residues reveal that the CG polyalanines fold into alpha helix and beta sheet structures. The 5-residue polyalanine displays a substantial increase in the beta strand fraction relative to the 12-residue polyalanine. The detailed conformational distribution is compared with those reported from recent all-atom simulations and experiments. The results suggest that the new coarse-graining approach presented in this study has the potential to offer both accuracy and efficiency for biomolecular modeling.
AB - A general, transferable coarse-grain (CG) framework based on the Gay-Berne potential and electrostatic point multipole expansion is presented for polypeptide simulations. The solvent effect is described by the Generalized Kirkwood theory. The CG model is calibrated using the results of all-atom simulations of model compounds in solution. Instead of matching the overall effective forces produced by atomic models, the fundamental intermolecular forces such as electrostatic, repulsion-dispersion, and solvation are represented explicitly at a CG level. We demonstrate that the CG alanine dipeptide model is able to reproduce quantitatively the conformational energy of all-atom force fields in both gas and solution phases, including the electrostatic and solvation components. Replica exchange molecular dynamics and microsecond dynamic simulations of polyalanine of 5 and 12 residues reveal that the CG polyalanines fold into alpha helix and beta sheet structures. The 5-residue polyalanine displays a substantial increase in the beta strand fraction relative to the 12-residue polyalanine. The detailed conformational distribution is compared with those reported from recent all-atom simulations and experiments. The results suggest that the new coarse-graining approach presented in this study has the potential to offer both accuracy and efficiency for biomolecular modeling.
UR - http://www.scopus.com/inward/record.url?scp=80155157254&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=80155157254&partnerID=8YFLogxK
U2 - 10.1063/1.3651626
DO - 10.1063/1.3651626
M3 - Article
C2 - 22029338
AN - SCOPUS:80155157254
VL - 135
JO - Journal of Chemical Physics
JF - Journal of Chemical Physics
SN - 0021-9606
IS - 15
M1 - 155104
ER -