Comparison of protein solution structures refined by molecular dynamics simulation in vacuum, with a generalized Born model, and with explicit water

Bin Xia, Vickie Tsui, David A. Case, H. Jane Dyson, Peter E. Wright

Research output: Contribution to journalArticlepeer-review

112 Scopus citations

Abstract

The inclusion of explicit solvent water in molecular dynamics refinement of NMR structures ought to provide the most physically meaningful accounting for the effects of solvent on structure, but is computationally expensive. In order to evaluate the validity of commonly used vacuum refinements and of recently developed continuum solvent model methods, we have used three different methods to refine a set of NMR solution structures of a medium sized protein, Escherichia coli glutaredoxin 2, from starting structures calculated using the program DYANA. The three different refinement protocols used molecular dynamics simulated annealing with the program AMBER in vacuum (VAC), including a generalized Born (GB) solvent model, and a full calculation including explicit solvent water (WAT). The structures obtained using the three methods of refinements were very similar, a reflection of their generally well-determined nature. However, the structures refined with the generalized Born model were more similar to those from explicit water refinement than those refined in vacuum. Significant improvement was seen in the percentage of backbone dihedral angels in the most favored regions of φ, ψ space and in hydrogen bond pattern for structures refined with the GB and WAT models, compared with the structures refined in vacuum. The explicit water calculation took an average of 200 h of CPU time per structure on an SGI cluster, compared to 15-90 h for the GB calculation (depending on the parameters used) and 2 h for the vacuum calculation. The generalized Born solvent model proved to be an excellent compromise between the vacuum and explicit water refinements, giving results comparable to those of the explicit water calculation. Some improvement for φ and ψ angle distribution and hydrogen bond pattern can also be achieved by energy minimizing the vacuum structures with the GB model, which takes a much shorter time than MD simulations with the GB model.

Original languageEnglish (US)
Pages (from-to)317-331
Number of pages15
JournalJournal of biomolecular NMR
Volume22
Issue number4
DOIs
StatePublished - 2002
Externally publishedYes

All Science Journal Classification (ASJC) codes

  • Biochemistry
  • Spectroscopy

Keywords

  • Hydrogen bond
  • Solvent refinement
  • Structure determination method
  • Structure refinement

Fingerprint

Dive into the research topics of 'Comparison of protein solution structures refined by molecular dynamics simulation in vacuum, with a generalized Born model, and with explicit water'. Together they form a unique fingerprint.

Cite this