Acetylcholinesterase: role of the enzyme's charge distribution in steering charged ligands toward the active site.

The electrostatic steering of charged ligands toward the active site of Torpedo californica acetylcholinesterase is investigated by Brownian dynamics simulations of wild type enzyme and several mutated forms, in which some normally charged residues are neutralized. The simulations reveal that the total ligand influx through a surface of 42 A radius centered in the enzyme monomer and separated from the protein surface by 1-14 A is not significantly influenced by electrostatic interactions. Electrostatic effects are visible for encounters with a surface of 32 A radius, which is partially hidden inside the protein, but mostly within the solvent. A clear accumulation of encounter events for that sphere is observed in the area directly above the entrance to the active site gorge. In this area, the encounter events are increased by 40% compared to the case of a neutral ligand. However, the differences among the encounter rates for the various mutants considered here are not pronounced, all rate constants being within +/- 10% of the average value. The enzyme charge distribution becomes more important as the charged ligand moves toward the bottom of the gorge, where the active site is located. We show that neither the enzyme's total charge, nor its dipole moment, fully account for the electrostatic steering of ligand to the active site. Higher moments of the enzyme's charge distribution are also important. However, for a series of mutations for which the direction of the enzyme dipole moment is constant within a few degrees, one observes a gradual decrease in the diffusional encounter rate constant with the number of neutralized residues. On the other hand, for other mutants that change the direction of the dipole moment from that of the wild type, the calculated encounter rate constants can be very close to that of the wild type. The present work yields two new insights to the kinetics of acetylcholinesterase. First, evolution appears to have built a redundant electrostatic steering capability into this important enzyme through the overall distribution of its thousands of partially charged atoms. And second, roughly half of the rate enhancement due to electrostatics arises from steering of the substrate outside the enzyme; the other half of the rate enhancement arises from improved trapping of the substrate after it has entered the gorge. The computational results reproduce qualitatively, and help to rationalize, many surprising experimental results obtained recently for human acetylcholinesterase.

[1]  Kim A. Sharp,et al.  Electrostatic interactions in macromolecules , 1994 .

[2]  A. Goldman,et al.  Atomic structure of acetylcholinesterase from Torpedo californica: a prototypic acetylcholine-binding protein , 1991, Science.

[3]  M. Gilson,et al.  Acetylcholinesterase: Effects of Ionic Strength and Dimerization on the Rate Constants , 1994 .

[4]  K. Sharp,et al.  Electrostatic interactions in macromolecules: theory and applications. , 1990, Annual review of biophysics and biophysical chemistry.

[5]  J L Sussman,et al.  Acetylcholinesterase: electrostatic steering increases the rate of ligand binding. , 1993, Biochemistry.

[6]  J. Antosiewicz,et al.  Electrostatics of hemoglobins from measurements of the electric dichroism and computer simulations. , 1995, Biophysical journal.

[7]  T L Rosenberry,et al.  Effective charge on acetylcholinesterase active sites determined from the ionic strength dependence of association rate constants with cationic ligands. , 1980, Biochemistry.

[8]  Electrostatic attraction by surface charge does not contribute to the catalytic efficiency of acetylcholinesterase. , 1994 .

[9]  J. Mccammon,et al.  Stochastically gated diffusion‐influenced reactions , 1982 .

[10]  J. Andrew McCammon,et al.  Diffusional dynamics of ligand-receptor association , 1986 .

[11]  L. R. Scott,et al.  Electrostatics and diffusion of molecules in solution: simulations with the University of Houston Brownian dynamics program , 1995 .

[12]  M. Gilson Multiple‐site titration and molecular modeling: Two rapid methods for computing energies and forces for ionizable groups in proteins , 1993, Proteins.

[13]  J. Mccammon,et al.  Simulation of the diffusion‐controlled reaction between superoxide and superoxide dismutase. II. Detailed models , 1985, Biopolymers.

[14]  P Taylor,et al.  The cholinesterases: from genes to proteins. , 1994, Annual review of pharmacology and toxicology.

[15]  A. Szabó,et al.  Diffusion-controlled bimolecular reaction rates. The effect of rotational diffusion and orientation constraints. , 1981, Biophysical journal.

[16]  M. Gilson,et al.  Acetylcholinesterase: diffusional encounter rate constants for dumbbell models of ligand. , 1995, Biophysical journal.

[17]  John M. Deutch,et al.  Diffusion-Controlled Reactions , 1983 .

[18]  M. Gilson,et al.  Simulation of charge-mutant acetylcholinesterases. , 1995, Biochemistry.

[19]  Malcolm E. Davis,et al.  Electrostatics in biomolecular structure and dynamics , 1990 .