Accuracy of continuum electrostatic calculations based on three common dielectric boundary definitions

Author:

Onufriev Alexey V.1,Aguilar Boris1

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篇论文的施引文献,订阅后可以查看论文全部施引文献

同舟云学术

1.学者识别学者识别

2.学术分析学术分析

3.人才评估人才评估

"同舟云学术"是以全球学者为主线,采集、加工和组织学术论文而形成的新型学术文献查询和分析系统,可以对全球学者进行文献检索和人才价值评估。用户可以通过关注某些学科领域的顶尖人物而持续追踪该领域的学科进展和研究前沿。经过近期的数据扩容,当前同舟云学术共收录了国内外主流学术期刊6万余种,收集的期刊论文及会议论文总量共计约1.5亿篇,并以每天添加12000余篇中外论文的速度递增。我们也可以为用户提供个性化、定制化的学者数据。欢迎来电咨询!咨询电话:010-8811{复制后删除}0370

www.globalauthorid.com

TOP

Copyright © 2019-2024 北京同舟云网络信息技术有限公司
京公网安备11010802033243号  京ICP备18003416号-3