Grid convergence of PyGBe with protein G B1 D4'
Reproducibility package containing data, running script, plotting script and final plot of a grid-convergence study on protein G B1 D4' near a charged surface.
The running script invokes the (open-source) bioelectrostatics solver PyGBe with the given input data, computes the extrapolated value of solvation plus surface energy (with Richardson extrapolation), and plots the final result of convergence.
This result is part of the following publication:
—"Poisson–Boltzmann model for protein-surface electrostatic interactions and grid-convergence study using the PyGBe code," Christopher D. Cooper and Lorena A. Barba, Comput. Phys. Comm. (2016) doi:10.1016/j.cpc.2015.12.019
PyGBe solves biomolecular electrostatics problems using an implicit-solvent model (Poisson-Boltzmann) and it uses GPU hardware for fast execution. It is written in Python, PyCUDA and CUDA.
More information about the PyGBe code in:
—"Validation of the PyGBe code for Poisson-Boltzmann equation with boundary element methods," Christopher Cooper, Lorena A. Barba. figshare.
http://dx.doi.org/10.6084/m9.figshare.154331
—"A biomolecular electrostatics solver using Python, GPUs and boundary elements that can handle solvent-filled cavities and Stern layers," Christopher D. Cooper, Jaydeep P. Bardhan, L. A. Barba. Comput. Phys. Comm., 185(3):720–729 (March 2014). 10.1016/j.cpc.2013.10.028 // Preprint arXiv:1309.4018
Acknowledgement:
This research is made possible by support from the Office of Naval Research, Applied Computational Analysis Program, N00014-11-1-0356. LAB also acknowledges support from NSF CAREER award OCI-1149784.