Pathways of ligand clearance in acetylcholinesterase by multiple copy sampling

Organoselenium Compounds as Acetylcholinesterase Inhibitors: Evidence and Mechanism of Mixed Inhibition

Acetylcholinesterase (AChE) inhibitors are actively used for the effective treatment of Alzheimer's disease. In recent years, the neuroprotective effects of organoselenium compounds such as ebselen and diselenides on the AChE activity have been investigated as potential therapeutic agents. In this work, we have carried out systematic kinetic and intrinsic fluorescence assays in combination with docking and molecular dynamics (MD) simulations to elucidate the molecular mechanism of the mixed inhibition of AChE by ebselen and diphenyl diselenide (DPDSe) molecules. Our MD simulations demonstrate significant heterogeneity in the binding modes and allosteric hotspots for DPDSe on AChE due to non-specific interactions. We have further identified that both ebselen and DPDSe can strongly bind around the peripheral anionic site (PAS), leading to non-competitive inhibition similar to other PAS-binding inhibitors. We also illustrate the entry of the DPDSe molecule into the gorge through a "side door", which offers an alternate entry point for AChE inhibitors as compared to the usual substrate entry point of the gorge. Together with results from experiments, these simulations provide mechanistic insights into the mixed type of inhibition for AChE using DPDSe as a promising inhibitor for AChE.

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.

A wrench in the works of human acetylcholinesterase: soman induced conformational changes revealed by molecular dynamics simulations

Irreversible inactivation of human acetylcholinesterase (hAChE) by organophosphorous pesticides (OPs) and chemical weapon agents (CWA) has severe morbidity and mortality consequences. We present data from quantum mechanics/molecular mechanics (QM/MM) and 80 classical molecular dynamics (MD) simulations of the apo and soman-adducted forms of hAChE to investigate the effects on the dynamics and protein structure when the catalytic Serine 203 is phosphonylated. We find that the soman phosphonylation of the active site Ser203 follows a water assisted addition-elimination mechanism with the elimination of the fluoride ion being the highest energy barrier at 6.5 kcal/mole. We observe soman-dependent changes in backbone and sidechain motions compared to the apo form of the protein. These alterations restrict the soman-adducted hAChE to a structural state that is primed for the soman adduct to be cleaved and removed from the active site. The altered motions and resulting structures provide alternative pathways into and out of the hAChE active site. In the soman-adducted protein both side and back door pathways are viable for soman adduct access. Correlation analysis of the apo and soman adducted MD trajectories shows that the correlation of gorge entrance and back door motion is disrupted when hAChE is adducted. This supports the hypothesis that substrate and product can use two different pathways as entry and exit sites in the apo form of the protein. These alternative pathways have important implications for the rational design of medical countermeasures.

Evaluation of binding competition and group epitopes of acetylcholinesterase inhibitors by STD NMR, Tr-NOESY, DOSY and molecular docking: an old approach but new findings

We suggest that the small inhibitors of AChE like coumarin and gallic acid's structure-based scaffolds may help in controlling Alzheimer's disease.

Molecular dynamics simulation study and molecular docking descriptors in structure-based QSAR on acetylcholinesterase (AChE) inhibitors

In this study we present an approach for predicting the inhibitory activity of acetylcholinesterase (AChE) inhibitors by combining molecular dynamics (MD) simulation and docking studies in a structure-based quantitative structure-activity relationship (QSAR) model. The MD simulation was performed on AChE to obtain enzyme conformation in a water environment. The resulting conformation of the enzyme was used for docking with the most potent inhibitor (26a). Docking analysis revealed that hydrophobic interactions play important roles in the AChE-inhibitor complex. Then, all inhibitors that could bind simultaneously at the catalytic site and at the peripheral anionic site of AChE were docked into the enzyme and their interactions with AChE were used as new interpretable descriptors in a structure-based QSAR model. The least squares support vector regression was constructed using the four most relevant docking descriptors and one molecular structure descriptor. The Q(2) value of the model was found to be 0.790. Furthermore, to study the enzyme conformation stability, a second MD simulation was performed on AChE-inhibitor 26a complex. In MD simulation, the topological parameters of the inhibitor were derived from the PRODRG server, and partial atomic charges were modified using the B3LYP/6-31G level of theory. The radius of gyration for the complex showed that AChE conformation did not change in the presence of the inhibitors.

Slow dissociation of a charged ligand: analysis of the primary quinone Q(A) site of photosynthetic bacterial reaction centers

Reaction centers (RCs) are integral membrane proteins that undergo a series of electron transfer reactions during the process of photosynthesis. In the Q(A) site of RCs from Rhodobacter sphaeroides, ubiquinone-10 is reduced, by a single electron transfer, to its semiquinone. The neutral quinone and anionic semiquinone have similar affinities, which is required for correct in situ reaction thermodynamics. A previous study showed that despite similar affinities, anionic quinones associate and dissociate from the Q(A) site at rates ≈10(4) times slower than neutral quinones indicating that anionic quinones encounter larger binding barriers (Madeo, J.; Gunner, M. R. Modeling binding kinetics at the Q(A) site in bacterial reaction centers. Biochemistry 2005, 44, 10994-11004). The present study investigates these barriers computationally, using steered molecular dynamics (SMD) to model the unbinding of neutral ground state ubiquinone (UQ) and its reduced anionic semiquinone (SQ(-)) from the Q(A) site. In agreement with experiment, the SMD unbinding barrier for SQ(-) is larger than for UQ. Multi Conformational Continuum Electrostatics (MCCE), used here to calculate the binding energy, shows that SQ(-) and UQ have comparable affinities. In the Q(A) site, there are stronger binding interactions for SQ(-) compared to UQ, especially electrostatic attraction to a bound non-heme Fe(2+). These interactions compensate for the higher SQ(-) desolvation penalty, allowing both redox states to have similar affinities. These additional interactions also increase the dissociation barrier for SQ(-) relative to UQ. Thus, the slower SQ(-) dissociation rate is a direct physical consequence of the additional binding interactions required to achieve a Q(A) site affinity similar to that of UQ. By a similar mechanism, the slower association rate is caused by stronger interactions between SQ(-) and the polar solvent. Thus, stronger interactions for both the unbound and bound states of charged and highly polar ligands can slow their binding kinetics without a conformational gate. Implications of this for other systems are discussed.

Long route or shortcut? A molecular dynamics study of traffic of thiocholine within the active-site gorge of acetylcholinesterase

The principal role of acetylcholinesterase is termination of nerve impulse transmission at cholinergic synapses, by rapid hydrolysis of the neurotransmitter acetylcholine to acetate and choline. Its active site is buried at the bottom of a deep and narrow gorge, at the rim of which is found a second anionic site, the peripheral anionic site. The fact that the active site is so deeply buried has raised cogent questions as to how rapid traffic of substrate and products occurs in such a confined environment. Various theoretical and experimental approaches have been used to solve this problem. Here, multiple conventional molecular dynamics simulations have been performed to investigate the clearance of the product, thiocholine, from the active-site gorge of acetylcholinesterase. Our results indicate that thiocholine is released from the peripheral anionic site via random pathways, while three exit routes appear to be favored for its release from the active site, namely, along the axis of the active-site gorge, and through putative back- and side-doors. The back-door pathway is that via which thiocholine exits most frequently. Our results are in good agreement with kinetic and kinetic-crystallography studies. We propose the use of multiple molecular dynamics simulations as a fast yet accurate complementary tool in structural studies of enzymatic trafficking.

Influence of the acetylcholinesterase active site protonation on omega loop and active site dynamics

Existence of alternative entrances in acetylcholinesterase (AChE) could explain the contrast between the very high AChE catalytic efficiency and the narrow and long access path to the active site revealed by X-ray crystallography. Alternative entrances could facilitate diffusion of the reaction products or at least water and ions from the active site. Previous molecular dynamics simulations identified side door and back door as the most probable alternative entrances. The simulations of non-inhibited AChE suggested that the back door opening events occur only rarely (0.8% of the time in the 10ns trajectory). Here we present a molecular dynamics simulation of non-inhibited AChE, where the back door opening appears much more often (14% of the time in the 12ns trajectory) and where the side door opening was observed quite frequently (78% of trajectory time). We also present molecular dynamics, where the back door does not open at all, or where large conformational changes of the AChE omega loop occur together with alternative passage opening events. All these differences in AChE dynamical behavior are caused by different protonation states of two glutamate residues located on bottom of the active site gorge (Glu202 and G450 in Mus musculus AChE). Our results confirm the results of previous molecular dynamics simulations, expand the view and suggest the probable reasons for the overall conformational behavior of AChE omega loop.

The discovery of potential acetylcholinesterase inhibitors: a combination of pharmacophore modeling, virtual screening, and molecular docking studies

BACKGROUND

Alzheimer's disease (AD) is the most common cause of dementia characterized by progressive cognitive impairment in the elderly people. The most dramatic abnormalities are those of the cholinergic system. Acetylcholinesterase (AChE) plays a key role in the regulation of the cholinergic system, and hence, inhibition of AChE has emerged as one of the most promising strategies for the treatment of AD.

METHODS

In this study, we suggest a workflow for the identification and prioritization of potential compounds targeted against AChE. In order to elucidate the essential structural features for AChE, three-dimensional pharmacophore models were constructed using Discovery Studio 2.5.5 (DS 2.5.5) program based on a set of known AChE inhibitors.

RESULTS

The best five-features pharmacophore model, which includes one hydrogen bond donor and four hydrophobic features, was generated from a training set of 62 compounds that yielded a correlation coefficient of R = 0.851 and a high prediction of fit values for a set of 26 test molecules with a correlation of R² = 0.830. Our pharmacophore model also has a high Güner-Henry score and enrichment factor. Virtual screening performed on the NCI database obtained new inhibitors which have the potential to inhibit AChE and to protect neurons from Aβ toxicity. The hit compounds were subsequently subjected to molecular docking and evaluated by consensus scoring function, which resulted in 9 compounds with high pharmacophore fit values and predicted biological activity scores. These compounds showed interactions with important residues at the active site.

CONCLUSIONS

The information gained from this study may assist in the discovery of potential AChE inhibitors that are highly selective for its dual binding sites.

The effect of Asp54 phosphorylation on the energetics and dynamics in the response regulator protein Spo0F studied by molecular dynamics

The response regulator protein Spo0F acts as an intermediate phospho-messenger in the signal transduction pathway that controls initiation of the differentiation process of sporculation in the bacterium Bacillus subtilis. The regulatory domain of Spo0F contains a triad of three conserved aspartate residues, whereof one aspartate (Asp54) is phosphorylated. Using molecular dynamics simulations, we have studied the changes in flexibility induced by phosphorylation and estimated the free energy cost of introducing a phosphate group at this position using alchemical free energy calculations. The deduced conformational flexibility compares well with experimental NMR results. We find that the apo-conformation of the protein explores a rough energy landscape resulting in a broad population of conformational substates. Phosphorylation of Spo0F reduces protein flexibility, and in particular, the so-called anchor and recognition regions exhibit lower mobility relative to the apo-conformation. Phosphorylation of Asp54 (P-Asp54), in which the apo-structure coordinates to the magnesium ion, results in extension of the sidechain, and depending on which carboxylate oxygen is phosphorylated, distinct interactions between P-Asp54 and magnesium ion as well as residues in its proximity are established. However, phosphorylation does not affect the coordination number of the magnesium ion yielding, within the statistical uncertainties, the same free energy change.

Pathways and mechanisms for product release in the engineered haloalkane dehalogenases explored using classical and random acceleration molecular dynamics simulations

Eight mutants of the DhaA haloalkane dehalogenase carrying mutations at the residues lining two tunnels, previously observed by protein X-ray crystallography, were constructed and biochemically characterized. The mutants showed distinct catalytic efficiencies with the halogenated substrate 1,2,3-trichloropropane. Release pathways for the two dehalogenation products, 2,3-dichloropropane-1-ol and the chloride ion, and exchange pathways for water molecules, were studied using classical and random acceleration molecular dynamics simulations. Five different pathways, denoted p1, p2a, p2b, p2c, and p3, were identified. The individual pathways showed differing selectivity for the products: the chloride ion releases solely through p1, whereas the alcohol releases through all five pathways. Water molecules play a crucial role for release of both products by breakage of their hydrogen-bonding interactions with the active-site residues and shielding the charged chloride ion during its passage through a hydrophobic tunnel. Exchange of the chloride ions, the alcohol product, and the waters between the buried active site and the bulk solvent can be realized by three different mechanisms: (i) passage through a permanent tunnel, (ii) passage through a transient tunnel, and (iii) migration through a protein matrix. We demonstrate that the accessibility of the pathways and the mechanisms of ligand exchange were modified by mutations. Insertion of bulky aromatic residues in the tunnel corresponding to pathway p1 leads to reduced accessibility to the ligands and a change in mechanism of opening from permanent to transient. We propose that engineering the accessibility of tunnels and the mechanisms of ligand exchange is a powerful strategy for modification of the functional properties of enzymes with buried active sites.

Flexibility of aromatic residues in the active-site gorge of acetylcholinesterase: X-ray versus molecular dynamics

The high aromatic content of the deep and narrow active-site gorge of acetylcholinesterase (AChE) is a remarkable feature of this enzyme. Here, we analyze conformational flexibility of the side chains of the 14 conserved aromatic residues in the active-site gorge of Torpedo californica AChE based on the 47 three-dimensional crystal structures available for the native enzyme, and for its complexes and conjugates, and on a 20-ns molecular dynamics (MD) trajectory of the native enzyme. The degree of flexibility of these 14 aromatic side chains is diverse. Although the side-chain conformations of F330 and W279 are both very flexible, the side-chain conformations of F120, W233, W432, Y70, Y121, F288, F290 and F331 appear to be fixed. Residues located on, or adjacent to, the Omega-loop (C67-C94), namely W84, Y130, Y442, and Y334, display different flexibilities in the MD simulations and in the crystal structures. An important outcome of our study is that the majority of the side-chain conformations observed in the 47 Torpedo californica AChE crystal structures are faithfully reproduced by the MD simulation on the native enzyme. Thus, the protein can assume these conformations even in the absence of the ligand that permitted their experimental detection. These observations are pertinent to structure-based drug design.

Peripheral site acetylcholinesterase blockade induces RACK1-associated neuronal remodeling

BACKGROUND

Peripheral anionic site (PAS) blockade of acetylcholinesterase (AChE) notably affects neuronal activity and cyto-architecture, however, the mechanism(s) involved are incompletely understood.

OBJECTIVE

We wished to specify the PAS extracellular effects on specific AChE mRNA splice variants, delineate the consequent cellular remodeling events, and explore the inhibitory effects on interchanging RACK1 interactions.

METHODS

We exposed rat hippocampal cultured neurons to BW284C51, the peripheral anionic site inhibitor of AChE, and to the non-selective AChE active site inhibitor, physostigmine for studying the neuronal remodeling of AChE mRNA expression and trafficking.

RESULTS

BW284C51 induced overexpression of both AChE splice variants, yet promoted neuritic translocation of the normally rare AChE-R, and retraction of AChE-S mRNA in an antisense-suppressible manner. BW284C51 further caused modest decreases in the expression of the scaffold protein RACK1 (receptor for activated protein kinase betaII), followed by drastic neurite retraction of both RACK1 and the AChE homologue neuroligin1, but not the tubulin-associated MAP2 protein. Accompanying BW284C51 effects involved decreases in the Fyn kinase and membrane insertion of the glutamate receptor NR2B variant and impaired glutamatergic activities of treated cells. Intriguingly, molecular modeling suggested that direct, non-catalytic competition with Fyn binding by the RACK1-interacting AChE-R variant may be involved.

CONCLUSIONS

Our findings highlight complex neuronal AChE-R/RACK1 interactions and are compatible with the hypothesis that peripheral site AChE inhibitors induce RACK1-mediated neuronal remodeling, promoting suppressed glutamatergic neurotransmission.

Molecular docking study on the “back door” hypothesis for product clearance in acetylcholinesterase

Acetylcholinesterase (AChE) is one of the fastest enzymes known, even though the active site is buried inside the protein at the end of a 20-A deep narrow gorge. Among the great variety of crystal structures of this enzyme, both in the absence and presence of various ligands and proteins, the structure of a complex of AChE with the pseudo-irreversible inhibitor Mf268 is of particular interest, as it assists in the proposal of a back door for product clearance from the active site. Binding of Mf268 to AChE results in the carbamoylation of Ser200 and liberation of an eseroline-fragment as the leaving group. The crystal structure of the AChE-Mf268 complex, however, proves that eseroline has escaped from the enzyme, despite the fact that the Ser-bound inhibitor fragment blocks the gorge entrance. The existence of alternative routes other than through the gorge for product clearance has been postulated but is still controversially discussed in the literature, as an experimental proof for such a back door is still missing. We have used Monte Carlo-based molecular docking methods in order to examine possible alternative pathways that could allow eseroline to be released from the protein after being cleaved from the substrate by Ser200. Based on our results, a short channel at the bottom of the gorge seems to be the most probable back-door site, which begins at amino acid Trp84 and ends at the enzyme surface in a cavity close to amino acid Glu445. [Figure: see text].

How Does Huperzine A Enter and Leave the Binding Gorge of Acetylcholinesterase? Steered Molecular Dynamics Simulations

The entering and leaving processes of Huperzine A (HupA) binding with the long active-site gorge of Torpedo californica acetylcholinesterase (TcAChE) have been investigated by using steered molecular dynamics simulations. The analysis of the force required along the pathway shows that it is easier for HupA to bind to the active site of AChE than to disassociate from it, which for the first time interprets at the atomic level the previous experimental result that unbinding process of HupA is much slower than its binding process to AChE. The direct hydrogen bonds, water bridges, and hydrophobic interactions were analyzed during two steered molecular dynamics (SMD) simulations. Break of the direct hydrogen bond needs a great pulling force. The steric hindrance of bottleneck might be the most important factor to produce the maximal rupture force for HupA to leave the binding site but it has a little effect on the binding process of HupA with AChE. Residue Asp72 forms a lot of water bridges with HupA leaving and entering the AChE binding gorge, acting as a clamp to take out HupA from or put HupA into the active site. The flip of the peptide bond between Gly117 and Gly118 has been detected during both the conventional MD and SMD simulations. The simulation results indicate that this flip phenomenon could be an intrinsic property of AChE and the Gly117-Gly118 peptide bond in both HupA bound and unbound AChE structures tends to adopt the native enzyme structure. At last, in a vacuum the rupture force is increased up to 1500 pN while in water solution the greatest rupture force is about 800 pN, which means water molecules in the binding gorge act as lubricant to facilitate HupA entering or leaving the binding gorge.

Essential motions in a fungal lipase with bound substrate, covalently attached inhibitor and product

As an aid to understanding the influence of dynamic fluctuations during esterolytic catalysis, we follow protein flexibility at three different steps along the catalytic pathway from substrate binding to product clearance via a covalently attached inhibitor, which represents a transition-state mimic. We have applied a classical approach, using molecular dynamics simulations to monitor protein dynamics in the nanosecond regime. We filter out small amplitude fluctuations and focus on the anharmonic contributions to the overall dynamics. This 'essential dynamics' analysis reveals different modes of response along the pathway suggesting that binding, catalysis and product clearance occur along different energy surfaces. Motions in the enzyme with a covalently attached ligand are more complex and occur along several eigenvectors. The magnitudes of the fluctuations in these individual subspaces are significantly smaller than those observed for the substrate and product molecules, indicating that the energy surface is shallow and that a relatively large number of conformational substates are accessible. On the other hand, substrate binding and product release occur at distinct modes of the protein flexibility suggesting that these processes occur along rough energy surfaces with only a few minima. Detailed energetic analyses along the trajectories indicated that in all cases binding is dominated by van der Waals interactions. The carboxylate form of the product is stabilized by a tight hydrogen bond network involving in particular Ser82, which may be a potential cause of product inhibition. Considerations such as these should aid the understanding of mechanisms of substrate, inhibitor or product recognition and could become of importance in the design of new substrates or inhibitors for enzymes.

Normal mode analysis of macromolecular motions in a database framework: developing mode concentration as a useful classifying statistic

We investigated protein motions using normal modes within a database framework, determining on a large sample the degree to which normal modes anticipate the direction of the observed motion and were useful for motions classification. As a starting point for our analysis, we identified a large number of examples of protein flexibility from a comprehensive set of structural alignments of the proteins in the PDB. Each example consisted of a pair of proteins that were considerably different in structure given their sequence similarity. On each pair, we performed geometric comparisons and adiabatic-mapping interpolations in a high-throughput pipeline, arriving at a final list of 3,814 putative motions and standardized statistics for each. We then computed the normal modes of each motion in this list, determining the linear combination of modes that best approximated the direction of the observed motion. We integrated our new motions and normal mode calculations in the Macromolecular Motions Database, through a new ranking interface at http://molmovdb.org. Based on the normal mode calculations and the interpolations, we identified a new statistic, mode concentration, related to the mathematical concept of information content, which describes the degree to which the direction of the observed motion can be summarized by a few modes. Using this statistic, we were able to determine the fraction of the 3,814 motions where one could anticipate the direction of the actual motion from only a few modes. We also investigated mode concentration in comparison to related statistics on combinations of normal modes and correlated it with quantities characterizing protein flexibility (e.g., maximum backbone displacement or number of mobile atoms). Finally, we evaluated the ability of mode concentration to automatically classify motions into a variety of simple categories (e.g., whether or not they are "fragment-like"), in comparison to motion statistics. This involved the application of decision trees and feature selection (particular machine-learning techniques) to training and testing sets derived from merging the "list" of motions with manually classified ones.

A cellular automata model of ligand passage over a protein hydrodynamic landscape

The subject of ligand passage to an active site on a protein is addressed. Current views on the mechanism and the possible role of surface water are discussed. A theory is presented in which the pattern of hydropathic states of protein surface amino acid side chains is invoked as the influence on the relative hydrophobic effects of nearby water. The theory describes a ligand passage through the hydrodynamic near-surface water which exhibits temporary organized cavities resembling the chreodes introduced by Waddington. The passage of the ligand to the active site is facilitated by this dynamic mechanism. Cellular automata models of preferential directional diffusion through these chreodes support the theory. The theory may be invoked to explain a number of ligand-active site observations and serves as an idea for further studies.

Concentration-dependent reversible activation-inhibition of human butyrylcholinesterase by tetraethylammonium ion

Tetraalkylammonium (TAA) salts are well known reversible inhibitors of cholinesterases. However, at concentrations around 10 mm, they have been found to activate the hydrolysis of positively charged substrates, catalyzed by wild-type human butyrylcholinesterase (EC 3.1.1.8) [Erdoes, E.G., Foldes, F.F., Zsigmond, E.K., Baart, N. & Zwartz, J.A. (1958) Science 128, 92]. The present study was undertaken to determine whether the peripheral anionic site (PAS) of human BuChE (Y332, D70) and/or the catalytic substrate binding site (CS) (W82, A328) are involved in this phenomenon. For this purpose, the kinetics of butyrylthiocholine (BTC) hydrolysis by wild-type human BuChE, by selected mutants and by horse BuChE was carried out at 25 degreeC and pH 7.0 in the presence of tetraethylammonium (TEA). It appears that human enzymes with more intact structure of the PAS show more prominent activation phenomenon. The following explanation has been put forward: TEA competes with the substrate at the peripheral site thus inhibiting the substrate hydrolysis at the CS. As the inhibition by TEA is less effective than the substrate inhibition itself, it mimics activation. At the concentrations around 40 mm, well within the range of TEA competition at both substrate binding sites, it lowers the activity of all tested enzymes.

Molecular dynamics study of active-site interactions with tetracoordinate transients in acetylcholinesterase and its mutants

The role of active-site residues in the dealkylation reaction in the P(S)C(S) diastereomer of 2-(3,3-dimethylbutyl)methylphosphonofluoridate (soman)-inhibited Torpedo californica acetylcholinesterase (AChE) was investigated by full-scale molecular dynamics simulations using CHARMM: >400 ps equilibration was followed by 150-200 ps production runs with the fully solvated tetracoordinate phosphonate adduct of the wild-type, Trp84Ala and Gly199Gln mutants of AChE. Parallel simulations were carried out with the tetrahedral intermediate formed between serine-200 Ogamma of AChE and acetylcholine. We found that the NepsilonH in histidine H(+)-440 is positioned to protonate the oxygen in choline and thus promote its departure. In contrast, NepsilonH in histidine H(+)-440 is not aligned for a favourable proton transfer to the pinacolyl O to promote dealkylation, but electrostatic stabilization by histidine H(+)-440 of the developing anion on the phosphonate monoester occurs. Destabilizing interactions between residues and the alkyl fragment of the inhibitor enforce methyl migration from Cbeta to Calpha concerted with C-O bond breaking in soman-inhibited AChE. Tryptophan-84, phenyalanine-331 and glutamic acid-199 are within 3.7-3.9 A (1 A=10(-10) m) from a methyl group in Cbeta, 4.5-5.1 A from Cbeta and 4.8-5.8 A from Calpha, and can better stabilize the developing carbenium ion on Cbeta than on Calpha. The Trp84Ala mutation eliminates interactions between the incipient carbenium ion and the indole ring, but also reduces its interactions with phenylalanine-331 and aspartic acid-72. Tyrosine-130 promotes dealkylation by interacting with the indole ring of tryptophan-84. Glutamic acid-443 can influence the orientation of active-site residues through tyrosine-421, tyrosine-442 and histidine-440 in soman-inhibited AChE, and thus facilitate dealkylation.