An implementation of hydrophobic force in implicit solvent molecular dynamics simulation for packed proteins

Li L. Duan, Tong Zhu, Ye Mei, Qing G. Zhang, Bo Tang, John Zhang

Research output: Contribution to journalArticle

Abstract

MD simulations of five proteins in which helical chains are held together by hydrophobic packing were carried out to investigate the effect of hydrophobic force on simulated structures of these protein complexes in implicit generalized Born (GB) model. The simulation study employed three different methods to treat hydrophobic effect: the standard GB method that does not include explicit hydrophobic force, the LCPO method that includes explicit hydrophobic force based directly on solvent accessible surface area (SASA), and a proposed packing enforced GB (PEGB) method that includes explicit hydrophobic force based on the radius of gyration of the protein complex. Our simulation study showed that all five protein complexes were unpacked in the standard GB simulation (without explicit hydrophobic force). In the LCPO method, three of the five protein systems remained well packed during the simulation, indicating the need for an explicit hydrophobic force in GB model for these packed protein systems. However, two of the five systems were still unpacked during LCPO simulation. For comparison, all five protein systems remain well packed in simulation using the new PEGB method. Analysis shows that the failure of the LCPO method in two cases is related to the way that SASA changes during the unpacking process for these two systems. These examples showed that standard GB method without explicit hydrophobic force is not suitable for MD simulation of protein systems involving hydrophobic packing. A similar problem remains but to a much lesser extent in the LCPO method for some packed protein systems. The proposed PEGB method seems quite promising for MD simulation of large, multi-domain packed proteins in implicit solvent model.

Original languageEnglish (US)
Pages (from-to)2605-2612
Number of pages8
JournalJournal of Molecular Modeling
Volume19
Issue number6
DOIs
StatePublished - Jun 2013

Fingerprint

Molecular dynamics
molecular dynamics
proteins
Proteins
Computer simulation
simulation
gyration
radii

Keywords

  • GB model
  • Hydrophobic force
  • Packed proteins
  • Radius of gyration
  • Solvent accessible surface area

ASJC Scopus subject areas

  • Physical and Theoretical Chemistry
  • Computer Science Applications
  • Computational Theory and Mathematics
  • Catalysis
  • Organic Chemistry
  • Inorganic Chemistry

Cite this

An implementation of hydrophobic force in implicit solvent molecular dynamics simulation for packed proteins. / Duan, Li L.; Zhu, Tong; Mei, Ye; Zhang, Qing G.; Tang, Bo; Zhang, John.

In: Journal of Molecular Modeling, Vol. 19, No. 6, 06.2013, p. 2605-2612.

Research output: Contribution to journalArticle

Duan, Li L. ; Zhu, Tong ; Mei, Ye ; Zhang, Qing G. ; Tang, Bo ; Zhang, John. / An implementation of hydrophobic force in implicit solvent molecular dynamics simulation for packed proteins. In: Journal of Molecular Modeling. 2013 ; Vol. 19, No. 6. pp. 2605-2612.
@article{0618207c2a3e4d8c86e6db9c0153df8b,
title = "An implementation of hydrophobic force in implicit solvent molecular dynamics simulation for packed proteins",
abstract = "MD simulations of five proteins in which helical chains are held together by hydrophobic packing were carried out to investigate the effect of hydrophobic force on simulated structures of these protein complexes in implicit generalized Born (GB) model. The simulation study employed three different methods to treat hydrophobic effect: the standard GB method that does not include explicit hydrophobic force, the LCPO method that includes explicit hydrophobic force based directly on solvent accessible surface area (SASA), and a proposed packing enforced GB (PEGB) method that includes explicit hydrophobic force based on the radius of gyration of the protein complex. Our simulation study showed that all five protein complexes were unpacked in the standard GB simulation (without explicit hydrophobic force). In the LCPO method, three of the five protein systems remained well packed during the simulation, indicating the need for an explicit hydrophobic force in GB model for these packed protein systems. However, two of the five systems were still unpacked during LCPO simulation. For comparison, all five protein systems remain well packed in simulation using the new PEGB method. Analysis shows that the failure of the LCPO method in two cases is related to the way that SASA changes during the unpacking process for these two systems. These examples showed that standard GB method without explicit hydrophobic force is not suitable for MD simulation of protein systems involving hydrophobic packing. A similar problem remains but to a much lesser extent in the LCPO method for some packed protein systems. The proposed PEGB method seems quite promising for MD simulation of large, multi-domain packed proteins in implicit solvent model.",
keywords = "GB model, Hydrophobic force, Packed proteins, Radius of gyration, Solvent accessible surface area",
author = "Duan, {Li L.} and Tong Zhu and Ye Mei and Zhang, {Qing G.} and Bo Tang and John Zhang",
year = "2013",
month = "6",
doi = "10.1007/s00894-013-1798-8",
language = "English (US)",
volume = "19",
pages = "2605--2612",
journal = "Journal of Molecular Modeling",
issn = "1610-2940",
publisher = "Springer Verlag",
number = "6",

}

TY - JOUR

T1 - An implementation of hydrophobic force in implicit solvent molecular dynamics simulation for packed proteins

AU - Duan, Li L.

AU - Zhu, Tong

AU - Mei, Ye

AU - Zhang, Qing G.

AU - Tang, Bo

AU - Zhang, John

PY - 2013/6

Y1 - 2013/6

N2 - MD simulations of five proteins in which helical chains are held together by hydrophobic packing were carried out to investigate the effect of hydrophobic force on simulated structures of these protein complexes in implicit generalized Born (GB) model. The simulation study employed three different methods to treat hydrophobic effect: the standard GB method that does not include explicit hydrophobic force, the LCPO method that includes explicit hydrophobic force based directly on solvent accessible surface area (SASA), and a proposed packing enforced GB (PEGB) method that includes explicit hydrophobic force based on the radius of gyration of the protein complex. Our simulation study showed that all five protein complexes were unpacked in the standard GB simulation (without explicit hydrophobic force). In the LCPO method, three of the five protein systems remained well packed during the simulation, indicating the need for an explicit hydrophobic force in GB model for these packed protein systems. However, two of the five systems were still unpacked during LCPO simulation. For comparison, all five protein systems remain well packed in simulation using the new PEGB method. Analysis shows that the failure of the LCPO method in two cases is related to the way that SASA changes during the unpacking process for these two systems. These examples showed that standard GB method without explicit hydrophobic force is not suitable for MD simulation of protein systems involving hydrophobic packing. A similar problem remains but to a much lesser extent in the LCPO method for some packed protein systems. The proposed PEGB method seems quite promising for MD simulation of large, multi-domain packed proteins in implicit solvent model.

AB - MD simulations of five proteins in which helical chains are held together by hydrophobic packing were carried out to investigate the effect of hydrophobic force on simulated structures of these protein complexes in implicit generalized Born (GB) model. The simulation study employed three different methods to treat hydrophobic effect: the standard GB method that does not include explicit hydrophobic force, the LCPO method that includes explicit hydrophobic force based directly on solvent accessible surface area (SASA), and a proposed packing enforced GB (PEGB) method that includes explicit hydrophobic force based on the radius of gyration of the protein complex. Our simulation study showed that all five protein complexes were unpacked in the standard GB simulation (without explicit hydrophobic force). In the LCPO method, three of the five protein systems remained well packed during the simulation, indicating the need for an explicit hydrophobic force in GB model for these packed protein systems. However, two of the five systems were still unpacked during LCPO simulation. For comparison, all five protein systems remain well packed in simulation using the new PEGB method. Analysis shows that the failure of the LCPO method in two cases is related to the way that SASA changes during the unpacking process for these two systems. These examples showed that standard GB method without explicit hydrophobic force is not suitable for MD simulation of protein systems involving hydrophobic packing. A similar problem remains but to a much lesser extent in the LCPO method for some packed protein systems. The proposed PEGB method seems quite promising for MD simulation of large, multi-domain packed proteins in implicit solvent model.

KW - GB model

KW - Hydrophobic force

KW - Packed proteins

KW - Radius of gyration

KW - Solvent accessible surface area

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

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

U2 - 10.1007/s00894-013-1798-8

DO - 10.1007/s00894-013-1798-8

M3 - Article

VL - 19

SP - 2605

EP - 2612

JO - Journal of Molecular Modeling

JF - Journal of Molecular Modeling

SN - 1610-2940

IS - 6

ER -