State Key Laboratory of Scientific/Engineering Computing, Institute of Computational Mathematics and Scientific/Engineering Computing, Academy of Mathematics and Systems Science, Chinese Academy of Sciences, Beijing 100190, China

Department of Chemistry and Biochemistry, Center for Theoretical Biological Physics, Department of Pharmacology, Howard Hughes Medical Institute, University of California at San Diego, La Jolla, CA, 92093-0365, USA

Abstract

The Debye-Hückel limiting law (DHL) has often been used to estimate rate constants of diffusion-controlled reactions under different ionic strengths. Two main approximations are adopted in DHL: one is that the solution of the linearized Poisson-Boltzmann equation for a spherical cavity is used to estimate the excess electrostatic free energy of a solution; the other is that details of electrostatic interactions of the solutes are neglected. This makes DHL applicable only at low ionic strengths and dilute solutions (very low substrate/solute concentrations). We show in this work that through numerical solution of the Poisson-Nernst-Planck equations, diffusion-reaction processes can be studied at a variety of conditions including realistically concentrated solutions, high ionic strength, and certainly with non-equilibrium charge distributions. Reaction rate coefficients for the acetylcholine-acetylcholinesterase system are predicted to strongly depend on both ionic strength and substrate concentration. In particular, they increase considerably with increase of substrate concentrations at a fixed ionic strength, which is open to experimental testing. This phenomenon is also verified on a simple model, and is expected to be general for electrostatically attracting enzyme-substrate systems.

**PACS Codes**: 82.45.Tv, 87.15.Vv

**MSC Codes**: 92C30

1 Background

Electrostatically steered diffusion-reaction processes exist widely in chemistry and biochemistry

where _{on}, _{E }and _{I }are the charges of the enzyme and substrate involved in the interaction. The DHL says that the rate constant (in an electrostatically-steered process) decays exponentially with the increase of the square root of ionic strength, as is observed under some conditions

The finite concentration effect was recently studied using Brownian dynamics simulation

2 Methods

We use a continuum model to simulate the electrodiffusion processes. The theoretical background is introduced in ^{+ }and Cl^{-}), and the atomic charges of the enzyme. Usually, we use three NP equations to describe the diffusion of three mobile species (counterion, coion, and substrate) respectively in the PNP model:

where ^{i}(^{i }and charge ^{i}, ρ^{f }is the fixed atomic charge distribution on the enzyme, β is the inverse Boltzmann energy, ε is the dielectric coefficient, ϕ is the electrostatic potential determined by the Poisson equation. The flux is

The reaction between the enzyme and the diffusing substrate is modeled by an absorbing boundary condition on a reactive site represented by a molecular surface patch. It is worth noting that this treatment is due to the fact that acetylcholinesterase is considered a fast enzyme. But in the context of high concentration of acetylcholine, ca. 500 mM, the simultaneously absorbing boundary condition may not be proper due to limited diffusion speed. In such case, a more complicated boundary condition such as Robin boundary condition, or inclusion of coupled ordinary partial differential equations can serve as a better description of the reaction event. The implementation of these considerations would be a future direction. For consistency and convenience in the setup of the computational model, a same absorbing boundary condition is used in this work. The reaction rate

Calculations of steady-state rate coefficients are performed for the enzyme catalyzed degradation of acetylcholine (ACh), which is an electrostatically steered diffusion-controlled reaction _{+ }and _{- }are the total bulk concentrations of cation and anion respectively, and that _{subs }is the bulk concentration of substrate ACh. These bulk concentrations are used as the outer boundary conditions of the diffusion domain in solving the PNP equations _{+ }+ _{subs }- _{- }= 0. As a comparison, the condition _{+ }- _{- }= 0 will lead to quite different results, which will be addressed later.

3 Results

The reaction rate coefficient is shown as a function of ionic strength (= "spectator" + bulk substrate) for different prescribed substrate concentrations in Figure _{+ }+ _{subs }= 300 mM), the rate coefficient is 1.36 × 10^{11 }M^{-1}min^{-1 }for _{subs }= 1 mM and is increased to 3.28 × 10^{11 }M^{-1}min^{-1 }for _{subs }= 300 mM. The physical origins of the observed behavior can be explained as follows. If substrate concentration is not considered, as in most previous work based on the DHL, the concentration of the counter ion of the enzyme, i.e., _{+ }here, is equal to the concentration of the coion, i.e., _{+ }= _{-}. The counter ions are attracted and concentrated around the negatively charged active site, which serves to screen the Coulomb interaction between ACh molecules and AChE, hence slowing the association. When _{subs }is considered in the PNP model, to maintain the same ionic strength, _{+ }needs to be reduced by _{subs }compared with that in the familiar Debye-Hückel theory. This leads to a thinner counter-ion atmosphere around the active site, and it can not be compensated by the additional substrate (ACh) density that is relatively low due to reactant depletion that results from the absorbing boundary condition. In other words, in the resulting non-equilibrium state, the sum of counter-ion density and ACh density near the active site is lower than that obtained with the Boltzmann distribution for a +1e particle. The consequences are a reduced overall screening effect and thereby an enhanced reaction rate.

Reaction rate affected by both ionic strength and substrate concentration

**Reaction rate affected by both ionic strength and substrate concentration**. Reaction rate coefficients for (a-b) ACh-AChE reaction system, (c) an absorbing unit sphere model system. _{0 }is bulk substrate concentration;

The ionic atmosphere always screens the electrostatic interactions, and hence reduces the rate coefficients. At very high ionic strength, due to strong ionic screening effects, the electrostatic interactions become very weak. This is close to the pure diffusion case, and all the rate constants for different substrate concentrations are close to the pure diffusion-reaction rate constant.

4 Discussion

The phenomena observed above are expected to be general for attractive substrate-enzyme systems, which can be illustrated with an idealized sphere model. Figure

As a comparison, if we use an aforementioned neutrality condition _{+ }- _{- }= 0 and the substrate concentration is not counted into ionic strength as was done in our former work

However, it is worth pointing out an issue in this model that the reaction products, choline and acetic acid that will ionize and generate acetate and a proton, are also charged species. The distribution and diffuse of these added ionic species will definitely affect the local ionic strength, hence the reaction rate coefficient. Therefore, a more complete treatment of some enzymes that catalyze reactions of charged substrates should include additional species such as charged products of the reaction in the model. But this brings the methodology a new issue that is how to elaborate the current PNP model to include the product diffusion originated from the reactive site, which will lead to some additional lines of work in the future. Therefore, the present work is limited to the case in which product concentrations are small, so that experimental tests of the current model would need to be in the early steady-state regime.

5 Conclusions

To summarize, the DHL only applies to very dilute situations. Our numerical results show that for electrostatically steered diffusion-controlled reaction processes, the rate coefficients strongly depend on both ionic strength and substrate concentration. At the same ionic strength, the current model predicts that increasing substrate concentration results in significant increase in rate coefficients for the attractive substrate-enzyme systems in case the product concentration can be ignored (like in the early steady-state regime). We are extending the theory and simulation methods to account for finite product concentrations, which will allow for easier comparison with experiments.

Acknowledgements

This work is supported in part by the NIH, NSF, the Howard Hughes Medical Institute, National Biomedical Computing Resource, the NSF Center for Theoretical Biological Physics, the W. M. Keck Foundation, and Accelrys, Inc. BZ is partially funded by the Academy of Mathematics and Systems Science of Chinese Academy of Sciences, the State Key Laboratory of Scientific/Engineering Computing, and NSFC (NSFC10971218).