This Article Full Text Full Text (PDF) All Versions of this Article: biophysj.106.102533v1 92/10/3397    most recent Alert me when this article is cited Alert me if a correction is posted Services Similar articles in this journal Similar articles in PubMed Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via HighWire Citing Articles via Google Scholar Google Scholar Articles by Cheng, Y. Articles by McCammon, J. A. Search for Related Content PubMed PubMed Citation Articles by Cheng, Y. Articles by McCammon, J. A.

Finite Element Analysis of the Time-Dependent Smoluchowski Equation for Acetylcholinesterase Reaction Rate Calculations

Yuhui Cheng ^* , Jason K. Suen , Deqiang Zhang , Stephen D. Bond , Yongjie Zhang ^¶, Yuhua Song ^||, Nathan A. Baker ^||, Chandrajit L. Bajaj ^¶ **, Michael J. Holst  and J. Andrew McCammon ^*  

^* Howard Hughes Medical Institute, Department of Chemistry and Biochemistry, and Center for Theoretical Biological Physics, University of California at San Diego, La Jolla, California; Accelrys, San Diego, California; Department of Computer Science, University of Illinois at Urbana-Champaign, Urbana, Illinois; ^¶ Institute of Computational Engineering and Sciences, Center for Computational Visualization, University of Texas at Austin, Texas; ^|| Department of Biochemistry and Molecular Biophysics, Center for Computational Biology, Washington University in St. Louis, Missouri; ^** Department of Computer Sciences, University of Texas at Austin, Texas; and Department of Mathematics, Department of Pharmacology, University of California at San Diego, La Jolla, California

Correspondence: Address reprint requests to Y. Cheng, Tel.: 858-822-2771; E-mail: ycheng{at}mccammon.ucsd.edu.

This article describes the numerical solution of the time-dependent Smoluchowski equation to study diffusion in biomolecular systems. Specifically, finite element methods have been developed to calculate ligand binding rate constants for large biomolecules. The resulting software has been validated and applied to the mouse acetylcholinesterase (mAChE) monomer and several tetramers. Rates for inhibitor binding to mAChE were calculated at various ionic strengths with several different time steps. Calculated rates show very good agreement with experimental and theoretical steady-state studies. Furthermore, these finite element methods require significantly fewer computational resources than existing particle-based Brownian dynamics methods and are robust for complicated geometries. The key finding of biological importance is that the rate accelerations of the monomeric and tetrameric mAChE that result from electrostatic steering are preserved under the non-steady-state conditions that are expected to occur in physiological circumstances.

This article has been cited by other articles:

				A. A. Gorfe, C.-e. A. Chang, I. Ivanov, and J. A. McCammon Dynamics of the Acetylcholinesterase Tetramer Biophys. J., February 15, 2008; 94(4): 1144 - 1154. [Abstract] [Full Text] [PDF]