Affiliation:
1. Department of Computer Science and Department of Physics, Virginia Tech, Blacksburg, VA 24060, USA
Abstract
We investigate the influence of three common definitions of the solute/solvent dielectric boundary (DB) on the accuracy of the electrostatic solvation energy ΔG el computed within the Poisson–Boltzmann (PB) and the generalized Born (GB) models of implicit solvation. The test structures include small molecules, peptides and small proteins; explicit solvent ΔG el are used as accuracy reference. For common atomic radii sets BONDI, PARSE (and ZAP9 for small molecules) the use of van der Waals (vdW) DB results, on average, in considerably larger errors in ΔG el than the molecular surface (MS) DB. The optimal probe radius ρw for which the MS DB yields the most accurate ΔG el varies considerably between structure types. The solvent accessible surface (SAS) DB becomes optimal at ρw ~ 0.2Å(exact value is sensitive to the structure and atomic radii), at which point the average accuracy of ΔG el is comparable to that of the MS-based boundary. The geometric equivalence of SAS to vdW surface based on the same atomic radii uniformly increased by ρw gives the corresponding optimal vdW DB. For small molecules, the optimal vdW DB based on BONDI +0.2 Å radii can yield ΔG el estimates at least as accurate as those based on the optimal MS DB. Also, in small molecules, pairwise charge–charge interactions computed with the optimal vdW DB are virtually equal to those computed with the MS DB, suggesting that in this case the two boundaries are practically equivalent by the electrostatic energy criteria. In structures other than small molecules, the optimal vdW and MS dielectric boundaries are not equivalent: the respective pairwise electrostatic interactions in the presence of solvent can differ by upto 5 kcal/mol for individual atomic pairs in small proteins, even when the total ΔG el are equal. For small proteins, the average decrease in pairwise electrostatic interactions resulting from the switch from optimal MS to optimal vdW DB definition can be mimicked within the MS DB definition by doubling of the solute dielectric constant. However, the use of the higher interior dielectric does not eliminate the large individual deviations between pairwise interactions computed within the two DB definitions. It is argued that while the MS-based definition of the DB is more physically correct in some types of practical calculations, the choice is not so clear in some other common scenarios.
Publisher
World Scientific Pub Co Pte Lt
Subject
Computational Theory and Mathematics,Physical and Theoretical Chemistry,Computer Science Applications
Cited by
22 articles.
订阅此论文施引文献
订阅此论文施引文献,注册后可以免费订阅5篇论文的施引文献,订阅后可以查看论文全部施引文献