Simulation of charge-mutant acetylcholinesterases

Computational Studies on Acetylcholinesterases

Functions of biomolecules, in particular enzymes, are usually modulated by structural fluctuations. This is especially the case in a gated diffusion-controlled reaction catalyzed by an enzyme such as acetylcholinesterase. The catalytic triad of acetylcholinesterase is located at the bottom of a long and narrow gorge, but it catalyzes the extremely rapid hydrolysis of the neurotransmitter, acetylcholine, with a reaction rate close to the diffusion-controlled limit. Computational modeling and simulation have produced considerable advances in exploring the dynamical and conformational properties of biomolecules, not only aiding in interpreting the experimental data, but also providing insights into the internal motions of the biomolecule at the atomic level. Given the remarkably high catalytic efficiency and the importance of acetylcholinesterase in drug development, great efforts have been made to understand the dynamics associated with its functions by use of various computational methods. Here, we present a comprehensive overview of recent computational studies on acetylcholinesterase, expanding our views of the enzyme from a microstate of a single structure to conformational ensembles, strengthening our understanding of the integration of structure, dynamics and function associated with the enzyme, and promoting the structure-based and/or mechanism-based design of new inhibitors for it.

Gorge Motions of Acetylcholinesterase Revealed by Microsecond Molecular Dynamics Simulations

Acetylcholinesterase, with a deep, narrow active-site gorge, attracts enormous interest due to its particularly high catalytic efficiency and its inhibitors used for treatment of Alzheimer’s disease. To facilitate the massive pass-through of the substrate and inhibitors, “breathing” motions to modulate the size of the gorge are an important prerequisite. However, the molecular mechanism that governs such motions is not well explored. Here, to systematically investigate intrinsic motions of the enzyme, we performed microsecond molecular dynamics simulations on the monomer and dimer of Torpedo californica acetylcholinesterase (TcAChE) as well as the complex of TcAChE bound with the drug E2020. It has been revealed that protein-ligand interactions and dimerization both keep the gorge in bulk, and opening events of the gorge increase dramatically compared to the monomer. Dynamics of three subdomains, S3, S4 and the Ω-loop, are tightly associated with variations of the gorge size while the dynamics can be changed by ligand binding or protein dimerization. Moreover, high correlations among these subdomains provide a basis for remote residues allosterically modulating the gorge motions. These observations are propitious to expand our understanding of protein structure and function as well as providing clues for performing structure-based drug design.

Internal Diffusion-Controlled Enzyme Reaction: The Acetylcholinesterase Kinetics

Acetylcholinesterase is an enzyme with a very high turnover rate; it quenches the neurotransmitter, acetylcholine, at the synapse. We have investigated the kinetics of the enzyme reaction by calculating the diffusion rate of the substrate molecule along an active site channel inside the enzyme from atomic-level molecular dynamics simulations. In contrast to the previous works, we have found that the internal substrate diffusion is the determinant of the acetylcholinesterase kinetics in the low substrate concentration limit. Our estimate of the overall bimolecular reaction rate constant for the enzyme is in good agreement with the experimental data. In addition, the present calculation provides a reasonable explanation for the effects of the ionic strength of solution and the mutation of surface residues of the enzyme. The study suggests that internal diffusion of the substrate could be a key factor in understanding the kinetics of enzymes of similar characteristics.

Redesigning protein pKa values

The ability to re-engineer enzymatic pH-activity profiles is of importance for industrial applications of enzymes. We theoretically explore the feasibility of re-engineering enzymatic pH-activity profiles by changing active site pK(a) values using point mutations. We calculate the maximum achievable DeltapK(a) values for 141 target titratable groups in seven enzymes by introducing conservative net-charge altering point mutations. We examine the importance of the number of mutations introduced, their distance from the target titratable group, and the characteristics of the target group itself. The results show that multiple mutations at 10A can change pK(a) values up to two units, but that the introduction of a requirement to keep other pK(a) values constant reduces the magnitude of the achievable DeltapK(a). The algorithm presented shows a good correlation with existing experimental data and is available for download and via a web server at http://enzyme.ucd.ie/pKD.

Protein complex formation by acetylcholinesterase and the neurotoxin fasciculin-2 appears to involve an induced-fit mechanism

Specific, rapid association of protein complexes is essential for all forms of cellular existence. The initial association of two molecules in diffusion-controlled reactions is often influenced by the electrostatic potential. Yet, the detailed binding mechanisms of proteins highly depend on the particular system. A complete protein complex formation pathway has been delineated by using structural information sampled over the course of the transformation reaction. The pathway begins at an encounter complex that is formed by one of the apo forms of neurotoxin fasciculin-2 (FAS2) and its high-affinity binding protein, acetylcholinesterase (AChE), followed by rapid conformational rearrangements into an intermediate complex that subsequently converts to the final complex as observed in crystal structures. Formation of the intermediate complex has also been independently captured in a separate 20-ns molecular dynamics simulation of the encounter complex. Conformational transitions between the apo and liganded states of FAS2 in the presence and absence of AChE are described in terms of their relative free energy profiles that link these two states. The transitions of FAS2 after binding to AChE are significantly faster than in the absence of AChE; the energy barrier between the two conformational states is reduced by half. Conformational rearrangements of FAS2 to the final liganded form not only bring the FAS2/AChE complex to lower energy states, but by controlling transient motions that lead to opening or closing one of the alternative passages to the active site of the enzyme also maximize the ligand's inhibition of the enzyme.

Modeling electrostatic effects in proteins

Electrostatic energies provide what is perhaps the most effective tool for structure-function correlation of biological molecules. This review considers the current state of simulations of electrostatic energies in macromolecules as well as the early developments of this field. We focus on the relationship between microscopic and macroscopic models, considering the convergence problems of the microscopic models and the fact that the dielectric 'constants' in semimacroscopic models depend on the definition and the specific treatment. The advances and the challenges in the field are illustrated considering a wide range of functional properties including pK(a)'s, redox potentials, ion and proton channels, enzyme catalysis, ligand binding and protein stability. We conclude by pointing out that, despite the current problems and the significant misunderstandings in the field, there is an overall progress that should lead eventually to quantitative descriptions of electrostatic effects in proteins and thus to quantitative descriptions of the function of proteins.

Alkaline phosphatase catalysis is ultrasensitive to charge sequestered between the active site zinc ions

Escherichia coli alkaline phosphatase (AP) is a prototypical bimetalloenzyme, facilitating catalysis of phosphate monoester hydrolysis with two Zn2+ metal ions that are only 4 A apart. In the reaction's transition state, one of the nonbridging oxygen atoms of the transferred group appears to interact directly with the Zn2+ ion metallocluster. To determine the importance and the energetic properties of this interaction, we systematically varied the charge on this oxygen atom, exploiting the ability of AP to catalyze reactions of different classes of substrates. We observed that the AP catalytic proficiency correlates very well (R2 = 0.98) with the charge on this oxygen atom, over 8 orders of magnitude of catalytic proficiency. The slope of this linear correlation (31 +/- 2 kcal/mol per unit charge) is extraordinarily steep, indicating that AP greatly discriminates between differentially charged substrates. We suggest that this discrimination arises via an electrostatic interaction with the bimetallocluster. The dependence of the AP catalytic proficiency on the nonbridging oxygen charge is much larger than charge perturbation effects observed previously for other proteins. We propose that AP uses folding energy to position the two Zn2+ metal ions in close proximity, thereby creating an active site with a high electrostatic potential that is extraordinarily sensitive to the charge that "solvates" the metallocluster. The sensitivity of enzyme energetics to systematic variation in electrostatic properties provides a powerful measure of the active site environment. Future work comparing the sensitivity of related enzymes that have been optimized to catalyze different reactions will help reveal how natural selection has tuned related active sites to favor different reactions.

Influence of the water structure on the acetylcholinesterase efficiency

We have studied the catalytic efficiency of acetylcholinesterase (AChE) in various solutions with ion-disturbed water structure to explore the role that the water structure plays in the substrate-enzyme encounter. The extent of water structuring in the different aqueous solutions was determined by near-infrared spectroscopy. The influence of water structure on the degree of solvation and on the intramolecular mobility of AChE was investigated for different aqueous ionic solutions by small-angle x-ray scattering technique and depolarization fluorescence spectroscopy. It was found that the encounter process between AChE and acetylthiocholine was promoted in solutions with less structured water. In these solutions it was also found that AChE is less solvated coinciding with higher intramolecular mobility. The found experimental results suggest that the water structure may influence the substrate-enzyme encounter process by diminishing the AChE solvation shell and may help diffusion of the substrate through the gorge by enhancing the intramolecular mobility of AChE.

Continuum diffusion reaction rate calculations of wild-type and mutant mouse acetylcholinesterase: adaptive finite element analysis

As described previously, continuum models, such as the Smoluchowski equation, offer a scalable framework for studying diffusion in biomolecular systems. This work presents new developments in the efficient solution of the continuum diffusion equation. Specifically, we present methods for adaptively refining finite element solutions of the Smoluchowski equation based on a posteriori error estimates. We also describe new, molecular-surface-based models, for diffusional reaction boundary criteria and compare results obtained from these models with the traditional spherical criteria. The new methods are validated by comparison of the calculated reaction rates with experimental values for wild-type and mutant forms of mouse acetylcholinesterase. The results show good agreement with experiment and help to define optimal reactive boundary conditions.

Edge distribution method for solving elliptic boundary value problems with boundary singularities

Elliptic boundary value problems are difficult to treat in the vicinity of singularities, i.e., edges and corners, of the boundary. The concentration of electrical charge on the edges and corners of a conductor is perhaps the simplest example of such problems. Here we provide a rapid method for accurate treatment of these problems. It utilizes a Green's-function-based implementation of last-passage Monte Carlo diffusion methods. This is combined with a diffusion algorithm for the scaling of solutions to the Laplace equation near a corner singularity to yield the solution of a benchmark problem: the charge distribution near the edge and corner of a conducting cube.

The association between a negatively charged ligand and the electronegative binding pocket of its receptor

Many examples exist of charged amino acids that play a role in attracting or holding a charged ligand toward or inside an oppositely charged binding pocket of the protein. For example, the enzymes superoxide dismutase, triose-phosphate isomerase, and acetylcholinesterase can steer ligands toward their oppositely charged binding pockets or gorges. Interestingly, in our Brownian dynamics simulations of a phosphate-binding protein, we discovered that negatively charged phosphate (HPO(2-)(4)) could make its way into the negatively charged binding pocket. In fact, the phosphate-binding protein exhibits counterintuitive kinetics of association. That is, one would expect that the rate of association would increase on increases to the ionic strength since the interaction between the ligand, with a charge of -2, and the electronegative binding pocket would be repulsive and greater screening should reduce this repulsion and increase the rate of association. However, the opposite is seen-i.e., the rate of association decreases on increases in the ionic strength. We used Brownian dynamics techniques to compute the diffusion limited association rate constants between the negatively charged phosphate ligand and several open forms of PBP (wild-type and several mutants based on an x-ray structure of open-form PBP, mutant T141D). With the appropriate choices of reaction criteria and molecular parameters, the ligand was able to diffuse into the binding pocket. A number of residues influence binding of the ligand within the pocket via hydrogen bonds or salt bridges. Arg135 partially neutralizes the charges on the HPO(2-)(4) ligand in the binding pocket, allowing it to enter. It is also found that the positive electrostatic patches above and below the binding entrance of PBP contribute the major attractive forces that direct the ligand toward the surface of the protein near the binding site.

Kinetics of diffusion-assisted reactions in microheterogeneous systems

This review is focused on the basic theory of diffusion-assisted reactions in microheterogeneous systems, from porous solids to self-organized colloids and biomolecules. Rich kinetic behaviors observed experimentally are explained in a unified fashion using simple concepts of competing distance and time scales of the reaction and the embedding structure. We mainly consider pseudo-first-order reactions, such as luminescence quenching, described by the Smoluchowski type of equation for the reactant pair distribution function with a sink term defined by the reaction mechanism. Microheterogeneity can affect the microscopic rate constant. It also enters the evolution equation through various spatial constraints leading to complicated boundary conditions and, possibly, to the reduction of dimensionality of the diffusion space. The reaction coordinate and diffusive motion along this coordinate are understood in a general way, depending on the problem at hand. Thus, the evolution operator can describe translational and rotational diffusion of molecules in a usual sense, it can be a discrete random walk operator when dealing with hopping of adsorbates in solids, or it can correspond to conformational fluctuations in proteins. Mathematical formulation is universal but physical consequences can be different. Understanding the principal features of reaction kinetics in microheterogeneous systems enables one to extract important structural and dynamical information about the host environments by analyzing suitably designed experiments, it helps building effective strategies for computer simulations, and ultimately opens possibilities for designing systems with controllable reactivity properties.

Comparison of salt effects on the reactions of acetylcholinesterase with cationic and anionic inhibitors

The influence of inorganic salts on the inhibition of acetylcholinesterase by charged organophosphorous inhibitors has been studied. It has been shown that the salt effect on the reaction of acetylcholinesterase with anionic bis(p-nitrophenyl) phosphate is determined by the influence of added salts on the activity coefficient of the inhibitor. In contrast to the salt effects on the reaction of acetylcholinesterase with cationic compounds, it does not include contribution from the enzyme charges. The smaller salt effect in the case of anionic inhibitor can be explained assuming that the anionic inhibitor does not form a non-covalent complex with the enzyme before the phosphorylation step of the reaction. Comparison of salt effects on the substrate turnover showed that in the case of cholinesterases from natural sources they are larger than in the case of enzymes expressed in recombinant cell clones. The enhanced salt effects may result from post-translational modification of the enzyme.

Pathways of ligand clearance in acetylcholinesterase by multiple copy sampling

The clearance of seven different ligands from the deeply buried active-site of Torpedo californica acetylcholinesterase is investigated by combining multiple copy sampling molecular dynamics simulations, with the analysis of protein-ligand interactions, protein motion and the electrostatic potential sampled by the ligand copies along their journey outwards. The considered ligands are the cations ammonium, methylammonium, and tetramethylammonium, the hydrophobic methane and neopentane, and the anionic product acetate and its neutral form, acetic acid. We find that the pathways explored by the different ligands vary with ligand size and chemical properties. Very small ligands, such as ammonium and methane, exit through several routes. One involves the main exit through the mouth of the enzyme gorge, another is through the so-called back door near Trp84, and a third uses a side door at a direction of approximately 45 degrees to the main exit. The larger polar ligands, methylammonium and acetic acid, leave through the main exit, but the bulkiest, tetramethylammonium and neopentane, as well as the smaller acetate ion, remain trapped in the enzyme gorge during the time of the simulations. The pattern of protein-ligand contacts during the diffusion process is highly non-random and differs for different ligands. A majority is made with aromatic side-chains, but classical H-bonds are also formed. In the case of acetate, but not acetic acid, the anionic and neutral form, respectively, of one of the reaction products, specific electrostatic interactions with protein groups, seem to slow ligand motion and interfere with protein flexibility; protonation of the acetate ion is therefore suggested to facilitate clearance. The Poisson-Boltzmann formalism is used to compute the electrostatic potential of the thermally fluctuating acetylcholinesterase protein at positions actually visited by the diffusing ligand copies. Ligands of different charge and size are shown to sample somewhat different electrostatic potentials during their migration, because they explore different microscopic routes. The potential along the clearance route of a cation such as methylammonium displays two clear minima at the active and peripheral anionic site. We find moreover that the electrostatic energy barrier that the cation needs to overcome when moving between these two sites is small in both directions, being of the order of the ligand kinetic energy. The peripheral site thus appears to play a role in trapping inbound cationic ligands as well as in cation clearance, and hence in product release.

Electrostatic steering of substrate to acetylcholinesterase: analysis of field fluctuations

Based on previous molecular dynamics simulation results for acetylcholinesterase dimer, we calculate and analyse the electrostatic field fluctuations around the enzyme. The results show that dynamic features of the electrostatic field favor attraction of the positively-charged substrate. An Internet link to an animation of the results is also provided.

Territrem B, a tremorgenic mycotoxin that inhibits acetylcholinesterase with a noncovalent yet irreversible binding mechanism

Territrem B (TRB) is a fungal metabolite isolated from Aspergillus terreus shown previously to be a potent and irreversible inhibitor of acetylcholinesterase (AChE). In the present study, a number of binding and inhibition assays were carried out to further characterize the inhibitory effect of TRB. The results indicate that the binding of TRB (a) is much more selective than a well characterized selective inhibitor of AChE, BW284C51, (b) adopts a one-to-one stoichiometry with the enzyme, (c) cannot be undone by an AChE-regenerating oxime agent, which contrasts the ability of 8 M urea to release AChE-bound TRB, (d) is enhanced by high concentration NaCl but prevented, unless preincubated, by Triton X-100, and (e) exhibits quasi-first order kinetics with an overall inhibition constant of 0.01 nM(-1) min(-1). Together these results suggest a very different irreversible binding (a noncovalent type) from that of the covalent type, which involves typical irreversible AChE inhibitors such as diisopropylfluorophosphate and neostigmine. According to the prediction of a molecular modeling study, the distinct AChE inhibitory characteristics of TRB may arise from the inhibitor being noncovalently trapped within a unique active-site gorge structure of the enzyme. It was predicted that an optimal TRB. AChE binding would position a narrowing connection of the TRB structure at a constricted area near the entrance of the gorge, thereby providing a structural basis for the observed irreversible binding.

A modular treatment of molecular traffic through the active site of cholinesterase

We present a model for the molecular traffic of ligands, substrates, and products through the active site of cholinesterases (ChEs). First, we describe a common treatment of the diffusion to a buried active site of cationic and neutral species. We then explain the specificity of ChEs for cationic ligands and substrates by introducing two additional components to this common treatment. The first module is a surface trap for cationic species at the entrance to the active-site gorge that operates through local, short-range electrostatic interactions and is independent of ionic strength. The second module is an ionic-strength-dependent steering mechanism generated by long-range electrostatic interactions arising from the overall distribution of charges in ChEs. Our calculations show that diffusion of charged ligands relative to neutral isosteric analogs is enhanced approximately 10-fold by the surface trap, while electrostatic steering contributes only a 1.5- to 2-fold rate enhancement at physiological salt concentration. We model clearance of cationic products from the active-site gorge as analogous to the escape of a particle from a one-dimensional well in the presence of a linear electrostatic potential. We evaluate the potential inside the gorge and provide evidence that while contributing to the steering of cationic species toward the active site, it does not appreciably retard their clearance. This optimal fine-tuning of global and local electrostatic interactions endows ChEs with maximum catalytic efficiency and specificity for a positively charged substrate, while at the same time not hindering clearance of the positively charged products.

Theoretical and experimental investigations of electrostatic effects on acetylcholinesterase catalysis and inhibition

The role of electrostatics in the function of acetylcholinesterase (AChE) has been investigated by both theoretical and experimental approaches. Second-order rate constants (kE = k(cat)/Km) for acetylthiocholine (ATCh) turnover have been measured as a function of ionic strength of the reaction medium for wild-type and mutant AChEs. Also, binding and dissociation rate constants have been measured as a function of ionic strength for the respective charged and neutral transition state analog inhibitors m-(N,N,N-trimethylammonio)trifluoroacetophenone (TMTFA) and m-(t-butyl)trifluoroacetophenone (TBTFA). Linear free-energy correlations between catalytic rate constants and inhibition constants indicate that kE for ATCh turnover is rate limited by terminal binding events. Comparison of binding rate constants for TMTFA and TBTFA attests to the sizable electrostatic discrimination of AChE. Free energy profiles for cationic ligand release from the active sites of wild-type and mutant AChEs have been calculated via a model that utilizes the structure of T. californica AChE, a spherical ligand, and energy terms that account for electrostatic and van der Waals interactions and chemical potential. These calculations indicate that EA and EI complexes are not bound with respect to electrostatic interactions, which obviates the need for a 'back door' for cationic ligand release. Moreover, the computed energy barriers for ligand release give linear free-energy correlations with log(kE) for substrate turnover, which supports the general correctness of the computational model.

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

It is known that anionic surface residues play a role in the long-range electrostatic attraction between acetylcholinesterase and cationic ligands. In our current investigation, we show that anionic residues also play an important role in the behavior of the ligand within the active site gorge of acetylcholinesterase. Negatively charged residues near the gorge opening not only attract positively charged ligands from solution to the enzyme, but can also restrict the motion of the ligand once it is inside of the gorge. We use Brownian dynamics techniques to calculate the rate constant kon, for wild type and mutant acetylcholinesterase with a positively charged ligand. These calculations are performed by allowing the ligand to diffuse within the active site gorge. This is an extension of previously reported work in which a ligand was allowed to diffuse only to the enzyme surface. By setting the reaction criteria for the ligand closer to the active site, better agreement with experimental data is obtained. Although a number of residues influence the movement of the ligand within the gorge, Asp74 is shown to play a particularly important role in this function. Asp74 traps the ligand within the gorge, and in this way helps to ensure a reaction.

Butyrylcholinesterase-catalysed hydrolysis of aspirin, a negatively charged ester, and aspirin-related neutral esters

Although aspirin (acetylsalicylic acid) is negatively charged, it is hydrolysed by butyrylcholinesterase (BuChE). Catalytic parameters were determined in 100 mM Tris buffer, pH 7.4, in the presence and absence of metal cations. The presence of Ca2+ or Mg2+ (<100 mM) in buffer did not change the Km, but accelerated the rate of hydrolysis of aspirin by wild-type or D70G mutant BuChE by 5-fold. Turnover numbers were of the order of 5000-12000 min-1 for the wild-type enzyme and the D70G and D70K enzymes in 100 mM Tris, pH 7.4, containing 50 mM CaCl2 at 25 degreesC; Km values were 6 mM for wild-type, 16 mM for D70G and 38 mM for D70K. People with 'atypical' BuChE have the D70G mutation. The apparent inhibition seen at high aspirin concentration was not due to inhibition by excess substrate but to spontaneous hydrolysis of aspirin, causing inhibition by salicylate. The wild-type and D70G enzymes were competitively inhibited by salicylic acid; the D70K enzyme showed a complex parabolic inhibition, suggesting multiple binding. The effect of salicylate was substrate-dependent, the D70K mutant being activated by salicylate with butyrylthiocholine as substrate. Km value for wild-type enzyme was lower than for D70 mutants, suggesting that residue 70 located at the rim of the active site gorge was not the major site for the initial encounter aspirin-BuChE complex. On the other hand, the virtual absence of affinity of the W82A mutant for aspirin indicated that W82 was the major residue involved in formation of the Michaelis complex. Molecular modelling of aspirin binding to BuChE indicated perpendicular interactions between the aromatic rings of W82 and aspirin. Kinetic study of BuChE-catalysed hydrolysis of different acetyl esters showed that the rate limiting step was acetylation. The bimolecular rate constants for hydrolysis of aspirin by wild-type, D70G and D70K enzymes were found to be close to 1x106 M-1 min-1. These results support the contention that the electrostatic steering due to the negative electrostatic field of the enzyme plays a role in substrate binding, but plays no role in the catalytic steps, i.e. in the enzyme acetylation.

Electrostatic steering and ionic tethering in enzyme-ligand binding: insights from simulations

To bind at an enzyme's active site, a ligand must diffuse or be transported to the enzyme's surface, and, if the binding site is buried, the ligand must diffuse through the protein to reach it. Although the driving force for ligand binding is often ascribed to the hydrophobic effect, electrostatic interactions also influence the binding process of both charged and nonpolar ligands. First, electrostatic steering of charged substrates into enzyme active sites is discussed. This is of particular relevance for diffusion-influenced enzymes. By comparing the results of Brownian dynamics simulations and electrostatic potential similarity analysis for triose-phosphate isomerases, superoxide dismutases, and beta-lactamases from different species, we identify the conserved features responsible for the electrostatic substrate-steering fields. The conserved potentials are localized at the active sites and are the primary determinants of the bimolecular association rates. Then we focus on a more subtle effect, which we will refer to as "ionic tethering." We explore, by means of molecular and Brownian dynamics simulations and electrostatic continuum calculations, how salt links can act as tethers between structural elements of an enzyme that undergo conformational change upon substrate binding, and thereby regulate or modulate substrate binding. This is illustrated for the lipase and cytochrome P450 enzymes. Ionic tethering can provide a control mechanism for substrate binding that is sensitive to the electrostatic properties of the enzyme's surroundings even when the substrate is nonpolar.

Calculated electrostatic gradients in recombinant human H-chain ferritin

Calculations to determine the electrostatic potential of the iron storage protein ferritin, using the human H-chain homopolymer (HuHF), reveal novel aspects of the protein. Some of the charge density correlates well with regions previously identified as active sites in the protein. The three-fold channels, the putative ferroxidase sites, and the nucleation sites all show expectedly negative values of the electrostatic potential. However, the outer entrance to the three-fold channels are surrounded by regions of positive potential, creating an electrostatic field directed toward the interior cavity. This electrostatic gradient provides a guidance mechanism for cations entering the protein cavity, indicating the three-fold channel as the major entrance to the protein. Pathways from the three-fold channels, indicated by electrostatic gradients on the inner surface, lead to the ferroxidase center, the nucleation center and to the interior entrance to the four-fold channel. Six glutamic acid residues at the nucleation site give rise to a region of very negative potential, surrounding a small positively charged center due to the presence of two conserved arginine residues, R63, in close proximity (4.9 A), suggesting that electrostatic fields could also play a role in the nucleation process. A large gradient in the electrostatic potential at the 4-fold channel gives rise to a field directed outward from the internal cavity, indicating the possibility that this channel functions to expel cations from inside the protein. The 4-fold channel could therefore provide an exit pathway for protons during mineralization, or iron leaving the protein cavity during de-mineralization.

Correlation between rate of enzyme-substrate diffusional encounter and average Boltzmann factor around active site

The utility of the average Boltzmann factor around the active site of an enzyme as the predictor of the electrostatic enhancement of the substrate binding rate is tested on a set of data on wild-type acetylcholinesterase and 18 charge mutants recently obtained by Brownian dynamics simulations. A good correlation between the average Boltzmann factors and the substrate binding rate constants is found. The effects of single charge mutations on both the Boltzmann factor and the substrate binding rate constant are modest, i.e., < 5 fold increase or decrease. This is consistent with the experimental results of Shafferman et al. but does not support their suggestion that the overall rate of the catalytic reaction is not limited by the diffusional encounter of acetylcholinesterase and its substrate.

Electrostatic effects in macromolecules: fundamental concepts and practical modeling

The past few years have seen an exponential growth in the calculations of electrostatic energies of macromolecules and an increased recognition of the crucial role of electrostatic effects. This review considers the current state of the field. Focus is placed on calculations of pKas, redox potentials and binding energies in macromolecules and clarification of the fact that the value of the dielectric 'constant' of a protein depends on its definition and that small dielectric constants should not be used in describing charge-charge interactions by current continuum models.

Design of fast enzymes by optimizing interaction potential in active site

The diffusional encounter between substrate and enzyme, and hence catalytic efficiency, can be enhanced by mutating charged residues on the surface of the enzyme. In this paper we present a simple method for screening such mutations. This is based on our earlier result that electrostatic enhancement of the enzyme-substrate binding rate constant can be accounted for just by the interaction potential within the active site. Assuming that catalytic and structural integrity is maintained, the catalytic efficiency can be optimized by surface charge mutations which lead to stronger interaction potential within the active site. Application of the screening method on superoxide dismutase shows that only charge mutations close to the active site will have practical effect on the catalytic efficiency. This rationalizes a large number of findings obtained in previous simulation and experimental studies.

External and internal electrostatic potentials of cholinesterase models

The electrostatic potentials for the three-dimensional structures of cholinesterases from various species were calculated, using the Delphi algorithm, on the basis of the Poisson-Boltzmann equation. We used structures for Torpedo californica and mouse acetylcholinesterase, and built homology models of the human, Bungarus fasciatus, and Drosophila melanogaster acetylcholinesterases and human butyrylcholinesterase. All these structures reveal a negative external surface potential, in the area around the entrance to the active-site gorge, that becomes more negative as the rim of the gorge is approached. Moreover, in all cases, the potential becomes increasingly more negative along the central axis running down the gorge, and is largest at the base of the gorge, near the active site. Ten key acidic residues conserved in the sequence alignments of AChE from various species, both in the surface area near the entrance of the active-site gorge and at its base, appear to be primarily responsible for these potentials. The potentials are highly correlated among the structures examined, down to sequence identities as low as 35%. This indicates that they are a conserved property of the cholinesterase family, could serve to attract the positively charged substrate into and down the gorge to the active site, and may play other roles important for cholinesterase function.

Electrostatic influence on the kinetics of ligand binding to acetylcholinesterase. Distinctions between active center ligands and fasciculin

To explore the role that surface and active center charges play in electrostatic attraction of ligands to the active center gorge of acetylcholinesterase (AChE), and the influence of charge on the reactive orientation of the ligand, we have studied the kinetics of association of cationic and neutral ligands with the active center and peripheral site of AChE. Electrostatic influences were reduced by sequential mutations of six surface anionic residues outside of the active center gorge (Glu-84, Glu-91, Asp-280, Asp-283, Glu-292, and Asp-372) and three residues within the active center gorge (Asp-74 at the rim and Glu-202 and Glu-450 at the base). The peripheral site ligand, fasciculin 2 (FAS2), a peptide of 6.5 kDa with a net charge of +4, shows a marked enhancement of rate of association with reduction in ionic strength, and this ionic strength dependence can be markedly reduced by progressive neutralization of surface and active center gorge anionic residues. By contrast, neutralization of surface residues only has a modest influence on the rate of cationic m-trimethylammoniotrifluoroacetophenone (TFK+) association with the active serine, whereas neutralization of residues in the active center gorge has a marked influence on the rate but with little change in the ionic strength dependence. Brownian dynamics calculations for approach of a small cationic ligand to the entrance of the gorge show the influence of individual charges to be in quantitative accord with that found for the surface residues. Anionic residues in the gorge may help to orient the ligand for reaction or to trap the ligand. Bound FAS2 on AChE not only reduces the rate of TFK+ reaction with the active center but inverts the ionic strength dependence for the cationic TFK+ association with AChE. Hence it appears that TFK+ must traverse an electrostatic barrier at the gorge entry imparted by the bound FAS2 with its net charge of +4.

Backbone makes a significant contribution to the electrostatics of alpha/beta-barrel proteins

The electrostatic properties of seven alpha/beta-barrel enzymes selected from different evolutionary families were studied: triose phosphate isomerase, fructose-1,6-bisphosphate aldolase, pyruvate kinase, mandelate racemase, trimethylamine dehydrogenase, glycolate oxidase, and narbonin, a protein without any known enzymatic activity. The backbone of the alpha/beta-barrel has a distinct electrostatic field pattern, which is dipolar along the barrel axis. When the side chains are included in the calculations the general effect is to modulate the electrostatic pattern so that the electrostatic field is generally enhanced and is focused into a specific area near the active site. We use the electrostatic flux through a square surface near the active site to gauge the functionally relevant magnitude of the electrostatic field. The calculations reveal that in six out of the seven cases the backbone itself contributes greater than 45% of the total flux. The substantial electrostatic contribution of the backbone correlates with the known preference of alpha/beta-barrel enzymes for negatively charged substrates.

Incorporating protein conformational flexibility into the calculation of pH-dependent protein properties

A method for combining calculations of residue pKa's with changes in the position of polar hydrogens has been developed. The Boltzmann distributions of proton positions in hydroxyls and neutral titratable residues are found in the same Monte Carlo sampling procedure that determines the amino acid ionization states at each pH. Electrostatic, Lennard-Jones potentials, and torsion angle energies are considered at each proton position. Many acidic and basic residues are found to have significant electrostatic interactions with either a water- or hydroxyl-containing side chain. Protonation state changes are coupled to reorientation of the neighboring hydroxyl dipoles, resulting in smaller free energy differences between neutral and ionized residues than when the protein is held rigid. Multiconformation pH titration gives better agreement with the experimental pKa's for triclinic hen egg lysozyme than conventional rigid protein calculations. The hydroxyl motion significantly increases the protein dielectric response, making it sensitive to the composition of the local protein structure. More than one conformer per residue is often found at a given pH, providing information about the distribution of low-energy lysozyme structures.

Structure-based drug design: computational advances

Structure-based computational methods continue to enhance progress in the discovery and refinement of therapeutic agents. Several such methods and their applications are described. These include molecular visualization and molecular modeling, docking, fragment methods, 3-D database techniques, and free-energy perturbation. Related issues that are discussed include the use of simplified potential energy functions and the determination of the positions of tightly bound waters. Strengths and weaknesses of the various methods are described.

Role of ionic interactions in cholinesterase catalysis

Acetylcholinesterase (AChE, EC 3.1.1.7) is an enzyme terminating the transmission of nerve impulse in synapses by rapid and selective hydrolysis of the neurotransmitter acetylcholine. Recent years have added a considerable amount of structural knowledge about this protein as well as opened new perspectives to the study of the molecular mechanism of cholinesterase catalysis. In this paper the current state of understanding the molecular recognition by cholinesterases is critically surveyed with particular emphasis on the role of electrostatic interactions.

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.

Crystal structure of an acetylcholinesterase-fasciculin complex: interaction of a three-fingered toxin from snake venom with its target

BACKGROUND

Fasciculin (FAS), a 61-residue polypeptide purified from mamba venom, is a three-fingered toxin which is a powerful reversible inhibitor of acetylcholinesterase (AChE). Solution of the three-dimensional structure of the AChE/FAS complex would provide the first structure of a three-fingered toxin complexed with its target.

RESULTS

The structure of a complex between Torpedo californica AChE and fasciculin-II (FAS-II), from the venom of the green mamba (Dendroaspis angusticeps) was solved by molecular replacement techniques, and refined at 3.0 A resolution to an R-factor of 0.231. The structure reveals a stoichiometric complex with one FAS molecule bound to each AChE subunit. The AChE and FAS conformations in the complex are very similar to those in their isolated structures. FAS is bound at the 'peripheral' anionic site of AChE, sealing the narrow gorge leading to the active site, with the dipole moments of the two molecules roughly aligned. The high affinity of FAS for AChE is due to a remarkable surface complementarity, involving a large contact area (approximately 2000 A2) and many residues either unique to FAS or rare in other three-fingered toxins. The first loop, or finger, of FAS reaches down the outer surface of the thin aspect of the gorge. The second loop inserts into the gorge, with an unusual stacking interaction between Met33 in FAS and Trp279 in AChE. The third loop points away from the gorge, but the C-terminal residue makes contact with the enzyme.

CONCLUSIONS

Two conserved aromatic residues in the AChE peripheral anionic site make important contacts with FAS. The absence of these residues from chicken and insect AChEs and from butyrylcholinesterase explains the very large reduction in the affinity of these enzymes for FAS. Several basic residues in FAS make important contacts with AChE. The complementarity between FAS and AChE is unusual, inasmuch as it involves a number of charged residues, but lacks any intermolecular salt linkages.