Several methods for including the conformational flexibility of proteins in the calculation of titration curves are compared. The methods use the linearized Poisson-Boltzmann equation to calculate the electrostatic free energies of solvation and are applied to bovine pancreatic trypsin inhibitor (BPTI) and hen egg-white lysozyme (HEWL). An ensemble of conformations is generated by a molecular dynamics simulation of the proteins with explicit solvent. The average titration curve of the ensemble is calculated in three different ways: an average structure is used for the pK(a) calculation; the electrostatic interaction free energies are averaged and used for the pK(a) calculation; and the titration curve for each structure is calculated and the curves are averaged. The three averaging methods give very similar results and improve the pK(a) values to approximately the same degree. This suggests, in contrast to implications from other work, that the observed improvement of pK(a) values in the present studies is due not to averaging over an ensemble of structures, but rather to the generation of a single properly averaged structure for the pK(a) calculation. (C) 1998 Wiley-Liss, Inc.
%0 Journal Article
%1 hlwoodcock:T.1998a
%A van Vlijmen, H. W. T.
%A Schaefer, M.
%A Karplus, M.
%D 1998
%J Proteins-Structure Function and Genetics
%K poisson-boltzmann ph-dependence white ionizable egg- models electrostatistics solvent magnetic-resonance nuclear conformation dynamics dielectric trypsin-inhibitor bibtex-import average curves titration equation lysozyme groups pancreatic electrostatic protein molecular-dynamics model continuum molecular constant
%N 2
%P 145--158
%T Improving the accuracy of protein pk(a) calculations: conformational averaging versus the average structure
%V 33
%X Several methods for including the conformational flexibility of proteins in the calculation of titration curves are compared. The methods use the linearized Poisson-Boltzmann equation to calculate the electrostatic free energies of solvation and are applied to bovine pancreatic trypsin inhibitor (BPTI) and hen egg-white lysozyme (HEWL). An ensemble of conformations is generated by a molecular dynamics simulation of the proteins with explicit solvent. The average titration curve of the ensemble is calculated in three different ways: an average structure is used for the pK(a) calculation; the electrostatic interaction free energies are averaged and used for the pK(a) calculation; and the titration curve for each structure is calculated and the curves are averaged. The three averaging methods give very similar results and improve the pK(a) values to approximately the same degree. This suggests, in contrast to implications from other work, that the observed improvement of pK(a) values in the present studies is due not to averaging over an ensemble of structures, but rather to the generation of a single properly averaged structure for the pK(a) calculation. (C) 1998 Wiley-Liss, Inc.
@article{hlwoodcock:T.1998a,
abstract = {Several methods for including the conformational flexibility of proteins in the calculation of titration curves are compared. The methods use the linearized Poisson-Boltzmann equation to calculate the electrostatic free energies of solvation and are applied to bovine pancreatic trypsin inhibitor (BPTI) and hen egg-white lysozyme (HEWL). An ensemble of conformations is generated by a molecular dynamics simulation of the proteins with explicit solvent. The average titration curve of the ensemble is calculated in three different ways: an average structure is used for the pK(a) calculation; the electrostatic interaction free energies are averaged and used for the pK(a) calculation; and the titration curve for each structure is calculated and the curves are averaged. The three averaging methods give very similar results and improve the pK(a) values to approximately the same degree. This suggests, in contrast to implications from other work, that the observed improvement of pK(a) values in the present studies is due not to averaging over an ensemble of structures, but rather to the generation of a single properly averaged structure for the pK(a) calculation. (C) 1998 Wiley-Liss, Inc.},
added-at = {2006-06-16T05:03:46.000+0200},
author = {van Vlijmen, H. W. T. and Schaefer, M. and Karplus, M.},
biburl = {https://www.bibsonomy.org/bibtex/29804da57502df4d5d9596734b889b164/hlwoodcock},
citeulike-article-id = {569475},
comment = {125VE PROTEIN-STRUCT FUNCT GENET},
interhash = {8b1c09f4236ac6a23476965addcc17d3},
intrahash = {9804da57502df4d5d9596734b889b164},
journal = {Proteins-Structure Function and Genetics},
keywords = {poisson-boltzmann ph-dependence white ionizable egg- models electrostatistics solvent magnetic-resonance nuclear conformation dynamics dielectric trypsin-inhibitor bibtex-import average curves titration equation lysozyme groups pancreatic electrostatic protein molecular-dynamics model continuum molecular constant},
number = 2,
pages = {145--158},
priority = {2},
timestamp = {2006-06-16T05:03:46.000+0200},
title = {Improving the accuracy of protein pk(a) calculations: conformational averaging versus the average structure},
volume = 33,
year = 1998
}