Development of an effective polarizable bond method for biomolecular simulation

Xudong Xiao, Tong Zhu, Chang G. Ji, John Zhang

Research output: Contribution to journalArticle

Abstract

An effective polarizable bond (EPB) model has been developed for computer simulation of proteins. In this partial polarizable approach, all polar groups of amino acids are treated as polarizable, and the relevant polarizable parameters were determined by fitting to quantum calculated electrostatic properties of these polar groups. Extensive numerical tests on a diverse set of proteins (including 1IEP, 1MWE, 1NLJ, 4COX, 1PGB, 1K4C, 1MHN, 1UBQ, 1IGD) showed that this EPB model is robust in MD simulation and can correctly describe the structure and dynamics of proteins (both soluble and membrane proteins). Comparison of the computed hydrogen bond properties and dynamics of proteins with experimental data and with results obtained from the nonpolarizable force field clearly demonstrated that EPB can produce results in much better agreement with experiment. The averaged deviation of the simulated backbone N-H order parameter of the B3 immunoglobobulin-binding domain of streptococcal protein G from experimental observation is 0.0811 and 0.0332 for Amber99SB and EPB, respectively. This new model inherited the effective character of the classic force field and the fluctuating feature of previous polarizable models. Different from other polarizable models, the polarization cost energy is implicitly included in the present method. As a result, the present method avoids the problem of over polarization and is numerically stable and efficient for dynamics simulation. Finally, compared to the traditional fixed AMBER charge model, the present method only adds about 5% additional computational time and is therefore highly efficient for practical applications.

Original languageEnglish (US)
Pages (from-to)14885-14893
Number of pages9
JournalJournal of Physical Chemistry B
Volume117
Issue number48
DOIs
StatePublished - Dec 5 2013

Fingerprint

proteins
Proteins
simulation
field theory (physics)
Polarization
Computer simulation
polarization
amino acids
Amino acids
Electrostatics
Hydrogen bonds
Membrane Proteins
computerized simulation
electrostatics
hydrogen bonds
membranes
Membranes
costs
deviation
Amino Acids

ASJC Scopus subject areas

  • Physical and Theoretical Chemistry
  • Materials Chemistry
  • Surfaces, Coatings and Films

Cite this

Development of an effective polarizable bond method for biomolecular simulation. / Xiao, Xudong; Zhu, Tong; Ji, Chang G.; Zhang, John.

In: Journal of Physical Chemistry B, Vol. 117, No. 48, 05.12.2013, p. 14885-14893.

Research output: Contribution to journalArticle

Xiao, Xudong ; Zhu, Tong ; Ji, Chang G. ; Zhang, John. / Development of an effective polarizable bond method for biomolecular simulation. In: Journal of Physical Chemistry B. 2013 ; Vol. 117, No. 48. pp. 14885-14893.
@article{a869a1fd3d914bbd8114a39f1afdb590,
title = "Development of an effective polarizable bond method for biomolecular simulation",
abstract = "An effective polarizable bond (EPB) model has been developed for computer simulation of proteins. In this partial polarizable approach, all polar groups of amino acids are treated as polarizable, and the relevant polarizable parameters were determined by fitting to quantum calculated electrostatic properties of these polar groups. Extensive numerical tests on a diverse set of proteins (including 1IEP, 1MWE, 1NLJ, 4COX, 1PGB, 1K4C, 1MHN, 1UBQ, 1IGD) showed that this EPB model is robust in MD simulation and can correctly describe the structure and dynamics of proteins (both soluble and membrane proteins). Comparison of the computed hydrogen bond properties and dynamics of proteins with experimental data and with results obtained from the nonpolarizable force field clearly demonstrated that EPB can produce results in much better agreement with experiment. The averaged deviation of the simulated backbone N-H order parameter of the B3 immunoglobobulin-binding domain of streptococcal protein G from experimental observation is 0.0811 and 0.0332 for Amber99SB and EPB, respectively. This new model inherited the effective character of the classic force field and the fluctuating feature of previous polarizable models. Different from other polarizable models, the polarization cost energy is implicitly included in the present method. As a result, the present method avoids the problem of over polarization and is numerically stable and efficient for dynamics simulation. Finally, compared to the traditional fixed AMBER charge model, the present method only adds about 5{\%} additional computational time and is therefore highly efficient for practical applications.",
author = "Xudong Xiao and Tong Zhu and Ji, {Chang G.} and John Zhang",
year = "2013",
month = "12",
day = "5",
doi = "10.1021/jp4080866",
language = "English (US)",
volume = "117",
pages = "14885--14893",
journal = "Journal of Physical Chemistry B Materials",
issn = "1520-6106",
publisher = "American Chemical Society",
number = "48",

}

TY - JOUR

T1 - Development of an effective polarizable bond method for biomolecular simulation

AU - Xiao, Xudong

AU - Zhu, Tong

AU - Ji, Chang G.

AU - Zhang, John

PY - 2013/12/5

Y1 - 2013/12/5

N2 - An effective polarizable bond (EPB) model has been developed for computer simulation of proteins. In this partial polarizable approach, all polar groups of amino acids are treated as polarizable, and the relevant polarizable parameters were determined by fitting to quantum calculated electrostatic properties of these polar groups. Extensive numerical tests on a diverse set of proteins (including 1IEP, 1MWE, 1NLJ, 4COX, 1PGB, 1K4C, 1MHN, 1UBQ, 1IGD) showed that this EPB model is robust in MD simulation and can correctly describe the structure and dynamics of proteins (both soluble and membrane proteins). Comparison of the computed hydrogen bond properties and dynamics of proteins with experimental data and with results obtained from the nonpolarizable force field clearly demonstrated that EPB can produce results in much better agreement with experiment. The averaged deviation of the simulated backbone N-H order parameter of the B3 immunoglobobulin-binding domain of streptococcal protein G from experimental observation is 0.0811 and 0.0332 for Amber99SB and EPB, respectively. This new model inherited the effective character of the classic force field and the fluctuating feature of previous polarizable models. Different from other polarizable models, the polarization cost energy is implicitly included in the present method. As a result, the present method avoids the problem of over polarization and is numerically stable and efficient for dynamics simulation. Finally, compared to the traditional fixed AMBER charge model, the present method only adds about 5% additional computational time and is therefore highly efficient for practical applications.

AB - An effective polarizable bond (EPB) model has been developed for computer simulation of proteins. In this partial polarizable approach, all polar groups of amino acids are treated as polarizable, and the relevant polarizable parameters were determined by fitting to quantum calculated electrostatic properties of these polar groups. Extensive numerical tests on a diverse set of proteins (including 1IEP, 1MWE, 1NLJ, 4COX, 1PGB, 1K4C, 1MHN, 1UBQ, 1IGD) showed that this EPB model is robust in MD simulation and can correctly describe the structure and dynamics of proteins (both soluble and membrane proteins). Comparison of the computed hydrogen bond properties and dynamics of proteins with experimental data and with results obtained from the nonpolarizable force field clearly demonstrated that EPB can produce results in much better agreement with experiment. The averaged deviation of the simulated backbone N-H order parameter of the B3 immunoglobobulin-binding domain of streptococcal protein G from experimental observation is 0.0811 and 0.0332 for Amber99SB and EPB, respectively. This new model inherited the effective character of the classic force field and the fluctuating feature of previous polarizable models. Different from other polarizable models, the polarization cost energy is implicitly included in the present method. As a result, the present method avoids the problem of over polarization and is numerically stable and efficient for dynamics simulation. Finally, compared to the traditional fixed AMBER charge model, the present method only adds about 5% additional computational time and is therefore highly efficient for practical applications.

UR - http://www.scopus.com/inward/record.url?scp=84890053414&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=84890053414&partnerID=8YFLogxK

U2 - 10.1021/jp4080866

DO - 10.1021/jp4080866

M3 - Article

C2 - 24251550

AN - SCOPUS:84890053414

VL - 117

SP - 14885

EP - 14893

JO - Journal of Physical Chemistry B Materials

JF - Journal of Physical Chemistry B Materials

SN - 1520-6106

IS - 48

ER -