Numerical calculation of protein-ligand binding rates through solution of the Smoluchowski equation using smoothed particle hydrodynamics
Abstract
Background: The calculation of diffusion-controlled ligand binding rates is important for understanding enzyme mechanisms as well as designing enzyme inhibitors. Methods: We demonstrate the accuracy and effectiveness of a Lagrangian particle-based method, smoothed particle hydrodynamics (SPH), to study diffusion in biomolecular systems by numerically solving the time-dependent Smoluchowski equation for continuum diffusion. Unlike previous studies, a reactive Robin boundary condition (BC), rather than the absolute absorbing (Dirichlet) BC, is considered on the reactive boundaries. This new BC treatment allows for the analysis of enzymes with “imperfect” reaction rates. Results: The numerical method is first verified in simple systems and then applied to the calculation of ligand binding to a mouse acetylcholinesterase (mAChE) monomer. Rates for inhibitor binding to mAChE are calculated at various ionic strengths and compared with experiment and other numerical methods. We find that imposition of the Robin BC improves agreement between calculated and experimental reaction rates. Conclusions: Although this initial application focuses on a single monomer system, our new method provides a framework to explore broader applications of SPH in larger-scale biomolecular complexes by taking advantage of its Lagrangian particle-based nature.
- Authors:
-
- Pacific Northwest National Lab. (PNNL), Richland, WA (United States)
- Publication Date:
- Research Org.:
- Pacific Northwest National Laboratory (PNNL), Richland, WA (United States)
- Sponsoring Org.:
- USDOE Office of Science (SC), Advanced Scientific Computing Research (ASCR)
- OSTI Identifier:
- 1203555
- Grant/Contract Number:
- AC06-76RL01830
- Resource Type:
- Accepted Manuscript
- Journal Name:
- BMC Biophysics
- Additional Journal Information:
- Journal Volume: 8; Journal Issue: 1; Journal ID: ISSN 2046-1682
- Publisher:
- BioMed Central
- Country of Publication:
- United States
- Language:
- English
- Subject:
- 59 BASIC BIOLOGICAL SCIENCES; Diffusion, Smoluchowski equation, Smoothed particle hydrodynamics, Protein-ligand interactions, Binding rates, Acetylcholinesterase
Citation Formats
Pan, Wenxiao, Daily, Michael, and Baker, Nathan A. Numerical calculation of protein-ligand binding rates through solution of the Smoluchowski equation using smoothed particle hydrodynamics. United States: N. p., 2015.
Web. doi:10.1186/s13628-015-0021-y.
Pan, Wenxiao, Daily, Michael, & Baker, Nathan A. Numerical calculation of protein-ligand binding rates through solution of the Smoluchowski equation using smoothed particle hydrodynamics. United States. https://doi.org/10.1186/s13628-015-0021-y
Pan, Wenxiao, Daily, Michael, and Baker, Nathan A. Thu .
"Numerical calculation of protein-ligand binding rates through solution of the Smoluchowski equation using smoothed particle hydrodynamics". United States. https://doi.org/10.1186/s13628-015-0021-y. https://www.osti.gov/servlets/purl/1203555.
@article{osti_1203555,
title = {Numerical calculation of protein-ligand binding rates through solution of the Smoluchowski equation using smoothed particle hydrodynamics},
author = {Pan, Wenxiao and Daily, Michael and Baker, Nathan A.},
abstractNote = {Background: The calculation of diffusion-controlled ligand binding rates is important for understanding enzyme mechanisms as well as designing enzyme inhibitors. Methods: We demonstrate the accuracy and effectiveness of a Lagrangian particle-based method, smoothed particle hydrodynamics (SPH), to study diffusion in biomolecular systems by numerically solving the time-dependent Smoluchowski equation for continuum diffusion. Unlike previous studies, a reactive Robin boundary condition (BC), rather than the absolute absorbing (Dirichlet) BC, is considered on the reactive boundaries. This new BC treatment allows for the analysis of enzymes with “imperfect” reaction rates. Results: The numerical method is first verified in simple systems and then applied to the calculation of ligand binding to a mouse acetylcholinesterase (mAChE) monomer. Rates for inhibitor binding to mAChE are calculated at various ionic strengths and compared with experiment and other numerical methods. We find that imposition of the Robin BC improves agreement between calculated and experimental reaction rates. Conclusions: Although this initial application focuses on a single monomer system, our new method provides a framework to explore broader applications of SPH in larger-scale biomolecular complexes by taking advantage of its Lagrangian particle-based nature.},
doi = {10.1186/s13628-015-0021-y},
journal = {BMC Biophysics},
number = 1,
volume = 8,
place = {United States},
year = {Thu May 07 00:00:00 EDT 2015},
month = {Thu May 07 00:00:00 EDT 2015}
}
Web of Science
Works referenced in this record:
Multiple time step diffusive Langevin dynamics for proteins
journal, February 1998
- Eastman, Peter; Doniach, Sebastian
- Proteins: Structure, Function, and Genetics, Vol. 30, Issue 3
Continuum Diffusion Reaction Rate Calculations of Wild-Type and Mutant Mouse Acetylcholinesterase: Adaptive Finite Element Analysis
journal, September 2004
- Song, Yuhua; Zhang, Yongjie; Bajaj, Chandrajit L.
- Biophysical Journal, Vol. 87, Issue 3
Non-Boltzmann Rate Distributions in Stochastically Gated Reactions
journal, January 1999
- Baker, Nathan A.; McCammon, J. Andrew
- The Journal of Physical Chemistry B, Vol. 103, Issue 4
Electrostatics and diffusion of molecules in solution: simulations with the University of Houston Brownian dynamics program
journal, March 1991
- Davis, Malcolm E.; Madura, Jeffry D.; Luty, Brock A.
- Computer Physics Communications, Vol. 62, Issue 2-3
Smoothed particle hydrodynamics continuous boundary force method for Navier–Stokes equations subject to a Robin boundary condition
journal, February 2014
- Pan, Wenxiao; Bao, Jie; Tartakovsky, Alexandre M.
- Journal of Computational Physics, Vol. 259
Finite concentration effects on diffusion-controlled reactions
journal, January 2004
- Senapati, Sanjib; Wong, Chung F.; McCammon, J. Andrew
- The Journal of Chemical Physics, Vol. 121, Issue 16
Smoothed Particle Hydrodynamics and Its Diverse Applications
journal, January 2012
- Monaghan, J. J.
- Annual Review of Fluid Mechanics, Vol. 44, Issue 1
Self-consistent generalized Langevin equation for colloid dynamics
journal, November 2001
- Yeomans-Reyna, Laura; Medina-Noyola, Magdaleno
- Physical Review E, Vol. 64, Issue 6
A new smoothed particle hydrodynamics non-Newtonian model for friction stir welding: Process modeling and simulation of microstructure evolution in a magnesium alloy
journal, September 2013
- Pan, Wenxiao; Li, Dongsheng; Tartakovsky, Alexandre M.
- International Journal of Plasticity, Vol. 48
Finite Element Solution of the Steady-State Smoluchowski Equation for Rate Constant Calculations
journal, April 2004
- Song, Yuhua; Zhang, Yongjie; Shen, Tongye
- Biophysical Journal, Vol. 86, Issue 4
Very fast empirical prediction and rationalization of protein pKa values
journal, October 2005
- Li, Hui; Robertson, Andrew D.; Jensen, Jan H.
- Proteins: Structure, Function, and Bioinformatics, Vol. 61, Issue 4
Poisson–Nernst–Planck equations for simulating biomolecular diffusion–reaction processes I: Finite element solutions
journal, September 2010
- Lu, Benzhuo; Holst, Michael J.; Andrew McCammon, J.
- Journal of Computational Physics, Vol. 229, Issue 19
Diffusion Control in Biochemical Reactions
book, January 1974
- Eigen, Manfred
- Quantum Statistical Mechanics in the Natural Sciences
Finite Element Analysis of the Time-Dependent Smoluchowski Equation for Acetylcholinesterase Reaction Rate Calculations
journal, May 2007
- Cheng, Yuhui; Suen, Jason K.; Zhang, Deqiang
- Biophysical Journal, Vol. 92, Issue 10
A commentary on gating of the active site of triose phosphate isomerase: Brownian dynamics simulations of flexible peptide loops in the enzyme, by R. C. Wade, M. E. Davis, B. A. Luty, J. D. Madura, and J. A. McCammon
journal, January 1993
- Allison, S. A.
- Biophysical Journal, Vol. 64, Issue 1
Kinetics of diffusion-controlled enzymatic reactions with charged substrates
journal, January 2010
- Lu, Benzhuo; McCammon, J. Andrew
- PMC Biophysics, Vol. 3, Issue 1
Gating of the active site of triose phosphate isomerase: Brownian dynamics simulations of flexible peptide loops in the enzyme
journal, January 1993
- Wade, R. C.; Davis, M. E.; Luty, B. A.
- Biophysical Journal, Vol. 64, Issue 1
[22] Diffusion-controlled enzymatic reactions
book, January 1991
- Davis, Malcolm E.; Madura, Jeffry D.; Sines, Jacqueline
- Methods in Enzymology
PDB2PQR: an automated pipeline for the setup of Poisson-Boltzmann electrostatics calculations
journal, July 2004
- Dolinsky, T. J.; Nielsen, J. E.; McCammon, J. A.
- Nucleic Acids Research, Vol. 32, Issue Web Server, p. W665-W667
Fractional diffusion-limited component of reactions catalyzed by acetylcholinesterase
journal, January 1986
- Bazelyansky, Michael; Robey, Ellen; Kirsch, Jack F.
- Biochemistry, Vol. 25, Issue 1
Effective charge on acetylcholinesterase active sites determined from the ionic strength dependence of association rate constants with cationic ligands
journal, August 1980
- Nolte, Hans-Juergen; Rosenberry, Terrone L.; Neumann, Eberhard
- Biochemistry, Vol. 19, Issue 16
A Monte Carlo simulation study of the influence of internal motions on the molecular conformation deduced from two-dimensional nmr experiments
journal, November 1989
- Genest, D.
- Biopolymers, Vol. 28, Issue 11
Improving implicit solvent simulations: a Poisson-centric view
journal, April 2005
- Baker, Nathan A.
- Current Opinion in Structural Biology, Vol. 15, Issue 2
A novel method for modeling Neumann and Robin boundary conditions in smoothed particle hydrodynamics
journal, December 2010
- Ryan, Emily M.; Tartakovsky, Alexandre M.; Amon, Cristina
- Computer Physics Communications, Vol. 181, Issue 12
PDB2PQR: expanding and upgrading automated preparation of biomolecular structures for molecular simulations
journal, May 2007
- Dolinsky, T. J.; Czodrowski, P.; Li, H.
- Nucleic Acids Research, Vol. 35, Issue Web Server
Hydrolysis of low concentrations of the acetylthiocholine analogs acetyl(homo)thiocholine and acetyl(nor)thiocholine by acetylcholinesterase may be limited by selective gating at the enzyme peripheral site
journal, March 2013
- Beri, Veena; Auletta, Jeffrey T.; Maharvi, Ghulam M.
- Chemico-Biological Interactions, Vol. 203, Issue 1
Fast Monte Carlo Simulation Methods for Biological Reaction-Diffusion Systems in Solution and on Surfaces
journal, January 2008
- Kerr, Rex A.; Bartol, Thomas M.; Kaminsky, Boris
- SIAM Journal on Scientific Computing, Vol. 30, Issue 6
Lateral diffusion and aggregation. A Monte Carlo study
journal, January 1992
- Saxton, M. J.
- Biophysical Journal, Vol. 61, Issue 1
Stochastic Simulation of Chemical Kinetics
journal, May 2007
- Gillespie, Daniel T.
- Annual Review of Physical Chemistry, Vol. 58, Issue 1
Diffusion Control in Biochemical Reactions
book, January 1974
- Eigen, Manfred
- Quantum Statistical Mechanics in the Natural Sciences
Fractional diffusion-limited component of reactions catalyzed by acetylcholinesterase
journal, January 1986
- Bazelyansky, Michael; Robey, Ellen; Kirsch, Jack F.
- Biochemistry, Vol. 25, Issue 1
Effective charge on acetylcholinesterase active sites determined from the ionic strength dependence of association rate constants with cationic ligands
journal, August 1980
- Nolte, Hans-Juergen; Rosenberry, Terrone L.; Neumann, Eberhard
- Biochemistry, Vol. 19, Issue 16
Electrostatic Influence on the Kinetics of Ligand Binding to Acetylcholinesterase
journal, September 1997
- Radić, Zoran; Kirchhoff, Paul D.; Quinn, Daniel M.
- Journal of Biological Chemistry, Vol. 272, Issue 37
Acetylcholinesterase: electrostatic steering increases the rate of ligand binding
journal, January 1993
- Tan, Raymond C.; Truong, Thanh N.; McCammon, J. Andrew
- Biochemistry, Vol. 32, Issue 2
Continuum Diffusion Reaction Rate Calculations of Wild-Type and Mutant Mouse Acetylcholinesterase: Adaptive Finite Element Analysis
journal, September 2004
- Song, Yuhua; Zhang, Yongjie; Bajaj, Chandrajit L.
- Biophysical Journal, Vol. 87, Issue 3
Rapid binding of a cationic active site inhibitor to wild type and mutant mouse acetylcholinesterase: Brownian dynamics simulation including diffusion in the active site gorge
journal, December 1998
- Tara, Sylvia; Elcock, Adrian H.; Kirchhoff, Paul D.
- Biopolymers, Vol. 46, Issue 7
Electrostatics and diffusion of molecules in solution: simulations with the University of Houston Brownian dynamics program
journal, March 1991
- Davis, Malcolm E.; Madura, Jeffry D.; Luty, Brock A.
- Computer Physics Communications, Vol. 62, Issue 2-3
[22] Diffusion-controlled enzymatic reactions
book, January 1991
- Davis, Malcolm E.; Madura, Jeffry D.; Sines, Jacqueline
- Methods in Enzymology
A Monte Carlo simulation study of the influence of internal motions on the molecular conformation deduced from two-dimensional nmr experiments
journal, November 1989
- Genest, D.
- Biopolymers, Vol. 28, Issue 11
Fast Monte Carlo Simulation Methods for Biological Reaction-Diffusion Systems in Solution and on Surfaces
journal, January 2008
- Kerr, Rex A.; Bartol, Thomas M.; Kaminsky, Boris
- SIAM Journal on Scientific Computing, Vol. 30, Issue 6
Brownian dynamics of cytochrome c and cytochrome c peroxidase association
journal, July 1988
- Northrup, S.; Boles, J.; Reynolds, J.
- Science, Vol. 241, Issue 4861
Improving implicit solvent simulations: a Poisson-centric view
journal, April 2005
- Baker, Nathan A.
- Current Opinion in Structural Biology, Vol. 15, Issue 2
Finite concentration effects on diffusion-controlled reactions
journal, January 2004
- Senapati, Sanjib; Wong, Chung F.; McCammon, J. Andrew
- The Journal of Chemical Physics, Vol. 121, Issue 16
A Lattice Relaxation Algorithm for Three-Dimensional Poisson-Nernst-Planck Theory with Application to Ion Transport through the Gramicidin A Channel
journal, February 1999
- Kurnikova, Maria G.; Coalson, Rob D.; Graf, Peter
- Biophysical Journal, Vol. 76, Issue 2
Analysis of Synaptic Transmission in the Neuromuscular Junction Using a Continuum Finite Element Model
journal, October 1998
- Smart, Jason L.; McCammon, J. Andrew
- Biophysical Journal, Vol. 75, Issue 4
Finite Element Analysis of the Time-Dependent Smoluchowski Equation for Acetylcholinesterase Reaction Rate Calculations
journal, May 2007
- Cheng, Yuhui; Suen, Jason K.; Zhang, Deqiang
- Biophysical Journal, Vol. 92, Issue 10
Poisson–Nernst–Planck equations for simulating biomolecular diffusion–reaction processes I: Finite element solutions
journal, September 2010
- Lu, Benzhuo; Holst, Michael J.; Andrew McCammon, J.
- Journal of Computational Physics, Vol. 229, Issue 19
Kinetics of diffusion-controlled enzymatic reactions with charged substrates
journal, January 2010
- Lu, Benzhuo; McCammon, J. Andrew
- PMC Biophysics, Vol. 3, Issue 1
Smoothed Particle Hydrodynamics and Its Diverse Applications
journal, January 2012
- Monaghan, J. J.
- Annual Review of Fluid Mechanics, Vol. 44, Issue 1
Smoothed particle hydrodynamics non-Newtonian model for ice-sheet and ice-shelf dynamics
journal, June 2013
- Pan, W.; Tartakovsky, A. M.; Monaghan, J. J.
- Journal of Computational Physics, Vol. 242
A new smoothed particle hydrodynamics non-Newtonian model for friction stir welding: Process modeling and simulation of microstructure evolution in a magnesium alloy
journal, September 2013
- Pan, Wenxiao; Li, Dongsheng; Tartakovsky, Alexandre M.
- International Journal of Plasticity, Vol. 48
A novel method for modeling Neumann and Robin boundary conditions in smoothed particle hydrodynamics
journal, December 2010
- Ryan, Emily M.; Tartakovsky, Alexandre M.; Amon, Cristina
- Computer Physics Communications, Vol. 181, Issue 12
Smoothed particle hydrodynamics continuous boundary force method for Navier–Stokes equations subject to a Robin boundary condition
journal, February 2014
- Pan, Wenxiao; Bao, Jie; Tartakovsky, Alexandre M.
- Journal of Computational Physics, Vol. 259
Stochastic Simulation of Chemical Kinetics
journal, May 2007
- Gillespie, Daniel T.
- Annual Review of Physical Chemistry, Vol. 58, Issue 1
Hydrolysis of low concentrations of the acetylthiocholine analogs acetyl(homo)thiocholine and acetyl(nor)thiocholine by acetylcholinesterase may be limited by selective gating at the enzyme peripheral site
journal, March 2013
- Beri, Veena; Auletta, Jeffrey T.; Maharvi, Ghulam M.
- Chemico-Biological Interactions, Vol. 203, Issue 1
Conformation gating as a mechanism for enzyme specificity
journal, August 1998
- Zhou, Huan-Xiang; Wlodek, Stanislaw T.; McCammon, J. Andrew
- Proceedings of the National Academy of Sciences, Vol. 95, Issue 16
Non-Boltzmann Rate Distributions in Stochastically Gated Reactions
journal, January 1999
- Baker, Nathan A.; McCammon, J. Andrew
- The Journal of Physical Chemistry B, Vol. 103, Issue 4
On the calculation of diffusive reaction rates using Brownian dynamics simulations
journal, March 1990
- Zhou, Huan‐Xiang
- The Journal of Chemical Physics, Vol. 92, Issue 5
Divergent Structure–Activity Relationships of Structurally Similar Acetylcholinesterase Inhibitors
journal, September 2013
- Andersson, C. David; Forsgren, Nina; Akfur, Christine
- Journal of Medicinal Chemistry, Vol. 56, Issue 19
Very fast empirical prediction and rationalization of protein pKa values
journal, October 2005
- Li, Hui; Robertson, Andrew D.; Jensen, Jan H.
- Proteins: Structure, Function, and Bioinformatics, Vol. 61, Issue 4
PDB2PQR: expanding and upgrading automated preparation of biomolecular structures for molecular simulations
journal, May 2007
- Dolinsky, T. J.; Czodrowski, P.; Li, H.
- Nucleic Acids Research, Vol. 35, Issue Web Server
PDB2PQR: an automated pipeline for the setup of Poisson-Boltzmann electrostatics calculations
journal, July 2004
- Dolinsky, T. J.; Nielsen, J. E.; McCammon, J. A.
- Nucleic Acids Research, Vol. 32, Issue Web Server, p. W665-W667
Electrostatics of nanosystems: Application to microtubules and the ribosome
journal, August 2001
- Baker, N. A.; Sept, D.; Joseph, S.
- Proceedings of the National Academy of Sciences, Vol. 98, Issue 18, p. 10037-10041
Potential theory, path integrals and the Laplacian of the indicator
journal, November 2012
- Lange, Rutger-Jan
- Journal of High Energy Physics, Vol. 2012, Issue 11