×

Poisson-Boltzmann model for protein-surface electrostatic interactions and grid-convergence study using the PyGBe code. (English) Zbl 1348.92008

Summary: Interactions between surfaces and proteins occur in many vital processes and are crucial in biotechnology: the ability to control specific interactions is essential in fields like biomaterials, biomedical implants and biosensors. In the latter case, biosensor sensitivity hinges on ligand proteins adsorbing on bioactive surfaces with a favorable orientation, exposing reaction sites to target molecules. Protein adsorption, being a free-energy-driven process, is difficult to study experimentally. This paper develops and evaluates a computational model to study electrostatic interactions of proteins and charged nanosurfaces, via the Poisson-Boltzmann equation. We extended the implicit-solvent model used in the open-source code PyGBe to include surfaces of imposed charge or potential. This code solves the boundary integral formulation of the Poisson-Boltzmann equation, discretized with surface elements. PyGBe has at its core a treecode-accelerated Krylov iterative solver, resulting in \(O(N \log N)\) scaling, with further acceleration on hardware via multi-threaded execution on gpus. It computes solvation and surface free energies, providing a framework for studying the effect of electrostatics on adsorption. We derived an analytical solution for a spherical charged surface interacting with a spherical dielectric cavity, and used it in a grid-convergence study to build evidence on the correctness of our approach. The study showed the error decaying with the average area of the boundary elements, i.e., the method is \(O(1 / N)\), which is consistent with our previous verification studies using PyGBe. We also studied grid-convergence using a real molecular geometry (protein GB1D4’), in this case using Richardson extrapolation (in the absence of an analytical solution) and confirmed the \(O(1 / N)\) scaling. With this work, we can now access a completely new family of problems, which no other major bioelectrostatics solver, e.g. APBS, is capable of dealing with. PyGBe is open-source under an MIT license and is hosted under version control at https://github.com/barbagroup/pygbe. To supplement this paper, we prepared “reproducibility packages” consisting of running and post-processing scripts in Python for replicating the grid-convergence studies, all the way to generating the final plots, with a single command.

MSC:

92-08 Computational methods for problems pertaining to biology
92C05 Biophysics
78A70 Biological applications of optics and electromagnetic theory
78A30 Electro- and magnetostatics
78M15 Boundary element methods applied to problems in optics and electromagnetic theory
PDFBibTeX XMLCite
Full Text: DOI arXiv

References:

[1] Gray, J. J., The interaction of proteins with solid surfaces, Curr. Opin. Struct. Biol., 14, 110-115 (2004)
[2] Rabe, M.; Verdes, D.; Seeger, S., Understanding protein adsorption phenomena at solid surfaces, Adv. Colloid Interface Sci., 162, 87-106 (2011)
[3] Mijajlovic, M.; Penna, M. J.; Biggs, M. J., Free energy of adsorption for a peptide at a liquid/solid interface via nonequilibrium molecular dynamics, Langmuir, 29, 9, 2919-2926 (2013), pMID: 23394469
[4] Roux, B.; Simonson, T., Implicit solvent models, Biophys. Chem., 78, 1-20 (1999)
[5] Bardhan, J. P., Biomolecular electrostatics—I want your solvation (model), Comput. Sci. Disc., 5, 013001 (2012)
[6] Yoon, B. J.; Lenhoff, A., Computation of the electrostatic interaction energy between a protein and a charged surface, J. Phys. Chem., 96, 3130-3134 (1992)
[7] Roth, C. M.; Lenhoff, A. M., Electrostatic and van der Waals contributions to protein adsorption: computation of equilibrium constants, Langmuir, 9, 962-972 (1993)
[8] Asthagiri, D.; Lenhoff, A. M., Influence of structural details in modeling electrostatically driven protein adsorption, Langmuir, 13, 6761-6768 (1997)
[9] Yao, Y.; Lenhoff, A. M., Electrostatic contributions to protein retention in ion-exchange chromatography. 1. Cytochrome c variants, Anal. Chem., 76, 6743-6752 (2004)
[10] Yao, Y.; Lenhoff, A. M., Electrostatic contributions to protein retention in ion-exchange chromatography. 2. Proteins with various degrees of structural differences, Anal. Chem., 77, 2157-2165 (2005)
[11] Roth, C. M.; Neal, B. L.; Lenhoff, A. M., Van der Waals interactions involving proteins, Biophys. J., 70, 977-987 (1996)
[13] Cooper, C. D.; Bardhan, J. P.; Barba, L. A., A biomolecular electrostatics solver using Python, GPUs and boundary elements that can handle solvent-filled cavities and Stern layers, Comput. Phys. Comm., 185, 3, 720-729 (2014), Preprint on arXiv:/1309.4018 · Zbl 1360.78047
[14] Tajima, N.; Takai, M.; Ishihara, K., Significance of antibody orientation unraveled: Well-oriented antibodies recorded high binding affinity, Anal. Chem., 83, 1969-1976 (2011)
[15] Trilling, A. K.; Beekwilder, J.; Zuilhof, H., Antibody orientation on biosensor surfaces: a minireview, Analyst, 138, 1619-1627 (2013)
[16] Cooper, C. D.; Clementi, N. C.; Barba, L. A., Probing protein orientation near charged nanosurfaces for simulation-assisted biosensor design, J. Chem. Phys., 143, 12, Article 124709 pp. (2015), Preprint on arxiv:1503.08150
[17] Yoon, B. J.; Lenhoff, A. M., A boundary element method for molecular electrostatics with electrolyte effects, J. Comput. Chem., 11, 9, 1080-1086 (1990)
[18] Juffer, A. H.; Botta, E. F.F.; van Keulen, B. A.M.; van der Ploeg, A.; Berendsen, H. J.C., The electric potential of a macromolecule in a solvent: A fundamental approach, J. Comput. Phys., 97, 144-171 (1991) · Zbl 0743.65094
[19] Lu, B.; Cheng, X.; Huang, J.; McCammon, J. A., Order \(N\) algorithm for computation of electrostatic interactions in biomolecular systems, Proc. Natl. Acad. Sci. USA, 103, 51, 19314-19319 (2006)
[20] Bajaj, C.; Chen, S. C.; Rand, A., An efficient higher-order fast multipole boundary element solution for Poisson-Boltzmann-based molecular electrostatics, SIAM J. Sci. Comput., 33, 2, 826-848 (2011) · Zbl 1227.92005
[21] Altman, M. D.; Bardhan, J. P.; White, J. K.; Tidor, B., Accurate solution of multi-region continuum electrostatic problems using the linearized Poisson-Boltzmann equation and curved boundary elements, J. Comput. Chem., 30, 132-153 (2009)
[22] Geng, W. H.; Krasny, R., A treecode-accelerated boundary integral Poisson-Boltzmann solver for solvated biomolecules, J. Comput. Phys., 247, 62-78 (2013) · Zbl 1349.78084
[23] Hess, J. L.; Smith, A., Calculation of potential flow about arbitrary bodies, Prog. Aerosp. Sci., 8, 1-138 (1967) · Zbl 0204.25602
[25] Greengard, L.; Rokhlin, V., A fast algorithm for particle simulations, J. Comput. Phys., 73, 2, 325-348 (1987) · Zbl 0629.65005
[26] Barnes, J.; Hut, P., A hierarchical \(O(N \log N)\) force-calculation algorithm, Nature, 324, 446-449 (1986)
[27] Phillips, J. R.; White, J. K., A precorrected-FFT method for electrostatic analysis of complicated 3-D structures, IEEE Trans. Comput.-Aided Des. Integr. Circuits Syst., 16, 10, 1059-1072 (1997)
[28] Li, P.; Johnston, H.; Krasny, R., A Cartesian treecode for screened Coulomb interactions, J. Comput. Phys., 228, 3858-3868 (2009) · Zbl 1165.78304
[29] Chan, D. Y.; Mitchell, D. J., The free energy of an electrical double layer, J. Colloid Interface Sci., 95, 1, 193-197 (1983)
[30] Carnie, S. L.; Chan, D. Y., Interaction free energy between identical spherical colloidal particles: The linearized Poisson-Boltzmann theory, J. Colloid Interface Sci., 155, 297-312 (1993)
[31] Carnie, S. L.; Chan, D. Y.; Gunning, J. S., Electrical double layer interaction between dissimilar spherical colloidal particles and betwen a sphere and a plate: The linearized Poisson-boltzmann theory, Langmuir, 10, 2993-3009 (1994)
[32] Lotan, I.; Head-Gordon, T., An analytical electrostatic model for salt screened interactions between multiple proteins, J. Chem. Theory Comput., 2, 3, 541-555 (2006)
[33] Marcelja, S.; Mitchell, D. J.; Ninham, B. W.; Sculley, M. J., Role of solvent structure in solution theory, J.  Chem.  Soc. Faraday Trans. II, 73, 630-648 (1977)
[34] Sanner, M. F.; Olson, A. J.; Spehner, J.-C., Fast and robust computation of molecular surfaces, (Proceedings of the Eleventh Annual Symposium on Computational Geometry (1995), ACM), 406-407
[35] Dolinsky, T. J.; Nielsen, J. E.; McCammon, J. A.; Baker, N. A., PDB2PQR: an automated pipeline for the setup of Poisson-Boltzmann electrostatics calculations, Nucleic Acids Res., 32, W665-W667 (2004)
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.