[Chimera-users] APBS electrostatic potential
Elaine Meng
meng at cgl.ucsf.edu
Fri Jun 5 11:14:14 PDT 2015
Hi Andrés,
I would not recommend making any quantitative comparisons between Coulombic and Poisson-Boltzmann potentials: not only are the equations different, but you can easily change the magnitudes by changing the dielectric (e.g. increasing the dielectric constant from 1 to 10 will decrease Coulombic potential magnitudes by a factor of 10).
Nevertheless, for visual rendering they should (and in my experience, they do) produce similar qualitative conclusions of which surface regions are more negative or positive.
We have made the Coulombic calculations very easy to set up, whereas for Poisson-Boltzmann calculations, more preparations of the structure are required. For example, yes, you do need to add the hydrogens and charges (including on the hydrogens) before running APBS or else you will have a ridiculously negative net charge. This is because the force field choices in PDB2PQR for adding charges are derived for such all-atom structures,.
As per the APBS manpage,
<http://www.rbvi.ucsf.edu/chimera/docs/ContributedSoftware/apbs/apbs.html>
"A structure should be prepared for APBS calculations by reconstructing missing heavy atoms, adding hydrogens, and assigning atomic charges and radii. These tasks can be done with PDB2PQR alone or in combination with parts of Dock Prep. Atomic charges can be assigned with Add Charge or PDB2PQR, although the latter may be preferred because it includes force fields developed specifically for Poisson-Boltzmann calculations.”
For additional details on the Chimera PDB2PQR and Coulombic tools:
<http://www.rbvi.ucsf.edu/chimera/docs/ContributedSoftware/apbs/pdb2pqr.html>
<http://www.rbvi.ucsf.edu/chimera/docs/ContributedSoftware/coulombic/coulombic.html>
I hope this helps,
Elaine
-----
Elaine C. Meng, Ph.D.
UCSF Computer Graphics Lab (Chimera team) and Babbitt Lab
Department of Pharmaceutical Chemistry
University of California, San Francisco
On Jun 3, 2015, at 5:45 PM, Felipe Vasquez <anfelvas at gmail.com> wrote:
> Hi,
> I am interested on calculating the electrostatic potential of a molecular surface on a protein (specifically visualizing that potential on a previously cavity determined by CASTp). However, althought I did get reasonable results using Coulombic potential, I could not obtain similar results using APBS potential (I obtained a surface too electrostatic, charged negatively and the values have a huge range, not only from -10 to 10 as in the case of Coulombic potential). Should I modify some default values out of the calculation parameters in order to get a similar result to that obtain by Coulombic potential? And it is absolutely indispensable add hydrogens in an APBS estimation?
> Thanks in advance for your help.
> Best regards,
> Andrés Felipe Vásquez J., BSc, MSc.
More information about the Chimera-users
mailing list