External and internal electrostatic potentials of cholinesterase models

Rational Design of Multifunctional Ferulic Acid Derivatives Aimed for Alzheimer's and Parkinson's Diseases

Ferulic acid has numerous beneficial effects on human health, which are frequently attributed to its antioxidant behavior. In this report, many of them are reviewed, and 185 new ferulic acid derivatives are computationally designed using the CADMA-Chem protocol. Consequently, their chemical space was sampled and evaluated. To that purpose, selection and elimination scores were used, which are built from a set of descriptors accounting for ADME properties, toxicity, and synthetic accessibility. After the first screening, 12 derivatives were selected and further investigated. Their potential role as antioxidants was predicted from reactivity indexes directly related to the formal hydrogen atom transfer and the single electron transfer mechanisms. The best performing molecules were identified by comparisons with the parent molecule and two references: Trolox and α-tocopherol. Their potential as polygenic neuroprotectors was investigated through the interactions with enzymes directly related to the etiologies of Parkinson's and Alzheimer's diseases. These enzymes are acetylcholinesterase, catechol-O-methyltransferase, and monoamine oxidase B. Based on the obtained results, the most promising candidates (FA-26, FA-118, and FA-138) are proposed as multifunctional antioxidants with potential neuroprotective effects. The findings derived from this investigation are encouraging and might promote further investigations on these molecules.

The Insight of In Silico and In Vitro evaluation of Beta vulgaris phytochemicals against Alzheimer's disease targeting acetylcholinesterase

B. vulgaris extracts possess antioxidant, anti-inflammatory along with its role in improving memory disorders. Subsequently, in vitro and in silico studies of its purified phytochemicals may expand complementary and alternative Alzheimer’s therapeutic option. Super activation of acetylcholinesterase enzyme is associated explicitly with Alzheimer’s disease (AD) ultimately resulting in senile dementia. Hence, acetylcholinesterase enzyme inhibition is employed as a promising approach for AD treatment. Many FDA approved drugs are unable to cure the disease progression completely. The Present study was devised to explore the potential bioactive phytochemicals of B. vulgaris as alternative therapeutic agents against AD by conducting in vitro and in silico studies. To achieve this, chemical structures of phytochemicals were recruited from PubChem. Further, these compounds were analyzed for their binding affinities towards acetylcholinesterase (AChE) enzyme. Pharmacophoric ligand-based models showed major characteristics like, HBA, HBD, hydrophobicity, aromaticity and positively ionizable surface morphology for receptor binding. Virtual screening identified three hit compounds including betanin, myricetin and folic acid with least binding score compared to the reference drug, donepezil (-17 kcal/mol). Further, in vitro studies for anti-acetylcholinesterase activity of betanin and glycine betaine were performed. Dose response analysis showed 1.271 μM and 1.203 μM 50% inhibitory concentration (IC₅₀) values for betanin and glycine betaine compounds respectively. Our findings indicate that phytoconstituents of B. vulgaris can be implicated as an alternative therapeutic drug candidate for cognitive disorders like Alzheimer’s disease.

A Comprehensive Review of Cholinesterase Modeling and Simulation

The present article reviews published efforts to study acetylcholinesterase and butyrylcholinesterase structure and function using computer-based modeling and simulation techniques. Structures and models of both enzymes from various organisms, including rays, mice, and humans, are discussed to highlight key structural similarities in the active site gorges of the two enzymes, such as flexibility, binding site location, and function, as well as differences, such as gorge volume and binding site residue composition. Catalytic studies are also described, with an emphasis on the mechanism of acetylcholine hydrolysis by each enzyme and novel mutants that increase catalytic efficiency. The inhibitory activities of myriad compounds have been computationally assessed, primarily through Monte Carlo-based docking calculations and molecular dynamics simulations. Pharmaceutical compounds examined herein include FDA-approved therapeutics and their derivatives, as well as several other prescription drug derivatives. Cholinesterase interactions with both narcotics and organophosphate compounds are discussed, with the latter focusing primarily on molecular recognition studies of potential therapeutic value and on improving our understanding of the reactivation of cholinesterases that are bound to toxins. This review also explores the inhibitory properties of several other organic and biological moieties, as well as advancements in virtual screening methodologies with respect to these enzymes.

Novel cholinesterase paralogs of Schistosoma mansoni have perceived roles in cholinergic signalling and drug detoxification and are essential for parasite survival

Cholinesterase (ChE) function in schistosomes is essential for orchestration of parasite neurotransmission but has been poorly defined with respect to the molecules responsible. Interrogation of the S. mansoni genome has revealed the presence of three ChE domain-containing genes (Smche)s, which we have shown to encode two functional acetylcholinesterases (AChE)s (Smache1 –smp_154600 and Smache2 –smp_136690) and a butyrylcholinesterase (BChE) (Smbche1 –smp_125350). Antibodies to recombinant forms of each SmChE localized the proteins to the tegument of adults and schistosomula and developmental expression profiling differed among the three molecules, suggestive of functions extending beyond traditional cholinergic signaling. For the first time in schistosomes, we identified ChE enzymatic activity in fluke excretory/secretory (ES) products and, using proteomic approaches, attributed this activity to the presence of SmAChE1 and SmBChE1. Parasite survival in vitro and in vivo was significantly impaired by silencing of each smche, either individually or in combination, attesting to the essential roles of these molecules. Lastly, in the first characterization study of a BChE from helminths, evidence is provided that SmBChE1 may act as a bio-scavenger of AChE inhibitors as the addition of recombinant SmBChE1 to parasite cultures mitigated the effect of the anti-schistosome AChE inhibitor 2,2- dichlorovinyl dimethyl phosphate—dichlorvos (DDVP), whereas smbche1-silenced parasites displayed increased sensitivity to DDVP.

Cholinesterases—aceytlcholinesterases (AChE)s and butyrylcholinesterases (BChE)s—are multi-functional enzymes that play a pivotal role in the nervous system of parasites by regulating neurotransmission through acetylcholine hydrolysis. Herein, we provide a detailed characterization of schistosome cholinesterases using molecular, enzymatic and gene-silencing approaches and show evidence for these molecules having roles in addition to their neuronal function. Further, we demonstrate the importance of these proteins to parasite development and survival through gene knockdown experiments in laboratory animals, providing evidence for the use of these proteins in the development of novel intervention strategies against schistosomiasis.

Neurobiology and therapeutic utility of neurotoxins targeting postsynaptic mechanisms of neuromuscular transmission

The neuromuscular junction (NMJ) is the principal site for the translation of motor neurochemical signals to muscle activity. Therefore, the release and sensing machinery of acetylcholine (ACh) along with muscle contraction are two of the main targets of natural toxins and pathogens, causing paralysis. Given pharmacology and medical advances, the active ingredients of toxins that target postsynaptic mechanisms have become of major interest, showing promise as drug leads. Herein, we review key facets of prevalent toxins modulating the mechanisms of ACh sensing and generation of the postsynaptic response, with muscle contraction. We consider the correlation between their outstanding selectivity and potency plus effects on motor function, and discuss emerging data advocating their usage for the development of therapies alleviating neuromuscular dysfunction.

The impact of crystallization conditions on structure-based drug design: A case study on the methylene blue/acetylcholinesterase complex

Structure-based drug design utilizes apoprotein or complex structures retrieved from the PDB. >57% of crystallographic PDB entries were obtained with polyethylene glycols (PEGs) as precipitant and/or as cryoprotectant, but <6% of these report presence of individual ethyleneglycol oligomers. We report a case in which ethyleneglycol oligomers' presence in a crystal structure markedly affected the bound ligand's position. Specifically, we compared the positions of methylene blue and decamethonium in acetylcholinesterase complexes obtained using isomorphous crystals precipitated with PEG200 or ammonium sulfate. The ligands' positions within the active-site gorge in complexes obtained using PEG200 are influenced by presence of ethyleneglycol oligomers in both cases bound to W84 at the gorge's bottom, preventing interaction of the ligand's proximal quaternary group with its indole. Consequently, both ligands are ∼3.0Å further up the gorge than in complexes obtained using crystals precipitated with ammonium sulfate, in which the quaternary groups make direct π-cation interactions with the indole. These findings have implications for structure-based drug design, since data for ligand-protein complexes with polyethylene glycol as precipitant may not reflect the ligand's position in its absence, and could result in selecting incorrect drug discovery leads. Docking methylene blue into the structure obtained with PEG200, but omitting the ethyleneglycols, yields results agreeing poorly with the crystal structure; excellent agreement is obtained if they are included. Many proteins display features in which precipitants might lodge. It will be important to investigate presence of precipitants in published crystal structures, and whether it has resulted in misinterpreting electron density maps, adversely affecting drug design.

Molecular characterization of monoclonal antibodies that inhibit acetylcholinesterase by targeting the peripheral site and backdoor region

The inhibition properties and target sites of monoclonal antibodies (mAbs) Elec403, Elec408 and Elec410, generated against Electrophorus electricus acetylcholinesterase (AChE), have been defined previously using biochemical and mutagenesis approaches. Elec403 and Elec410, which bind competitively with each other and with the peptidic toxin inhibitor fasciculin, are directed toward distinctive albeit overlapping epitopes located at the AChE peripheral anionic site, which surrounds the entrance of the active site gorge. Elec408, which is not competitive with the other two mAbs nor fasciculin, targets a second epitope located in the backdoor region, distant from the gorge entrance. To characterize the molecular determinants dictating their binding site specificity, we cloned and sequenced the mAbs; generated antigen-binding fragments (Fab) retaining the parental inhibition properties; and explored their structure-function relationships using complementary x-ray crystallography, homology modeling and flexible docking approaches. Hypermutation of one Elec403 complementarity-determining region suggests occurrence of antigen-driven selection towards recognition of the AChE peripheral site. Comparative analysis of the 1.9Å-resolution structure of Fab408 and of theoretical models of its Fab403 and Fab410 congeners evidences distinctive surface topographies and anisotropic repartitions of charges, consistent with their respective target sites and inhibition properties. Finally, a validated, data-driven docking model of the Fab403-AChE complex suggests a mode of binding at the PAS that fully correlates with the functional data. This comprehensive study documents the molecular peculiarities of Fab403 and Fab410, as the largest peptidic inhibitors directed towards the peripheral site, and those of Fab408, as the first inhibitor directed toward the backdoor region of an AChE and a unique template for the design of new, specific modulators of AChE catalysis.

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.

Enzymatic toxins from snake venom: structural characterization and mechanism of catalysis

Snake venoms are cocktails of enzymes and non-enzymatic proteins used for both the immobilization and digestion of prey. The most common snake venom enzymes include acetylcholinesterases, l-amino acid oxidases, serine proteinases, metalloproteinases and phospholipases A(2) . Higher catalytic efficiency, thermal stability and resistance to proteolysis make these enzymes attractive models for biochemists, enzymologists and structural biologists. Here, we review the structures of these enzymes and describe their structure-based mechanisms of catalysis and inhibition. Some of the enzymes exist as protein complexes in the venom. Thus we also discuss the functional role of non-enzymatic subunits and the pharmacological effects of such protein complexes. The structures of inhibitor-enzyme complexes provide ideal platforms for the design of potent inhibitors which are useful in the development of prototypes and lead compounds with potential therapeutic applications.

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.

Acetylcholinesterase: from 3D structure to function

By rapid hydrolysis of the neurotransmitter, acetylcholine, acetylcholinesterase terminates neurotransmission at cholinergic synapses. Acetylcholinesterase is a very fast enzyme, functioning at a rate approaching that of a diffusion-controlled reaction. The powerful toxicity of organophosphate poisons is attributed primarily to their potent inhibition of acetylcholinesterase. Acetylcholinesterase inhibitors are utilized in the treatment of various neurological disorders, and are the principal drugs approved thus far by the FDA for management of Alzheimer's disease. Many organophosphates and carbamates serve as potent insecticides, by selectively inhibiting insect acetylcholinesterase. The determination of the crystal structure of Torpedo californica acetylcholinesterase permitted visualization, for the first time, at atomic resolution, of a binding pocket for acetylcholine. It also allowed identification of the active site of acetylcholinesterase, which, unexpectedly, is located at the bottom of a deep gorge lined largely by aromatic residues. The crystal structure of recombinant human acetylcholinesterase in its apo-state is similar in its overall features to that of the Torpedo enzyme; however, the unique crystal packing reveals a novel peptide sequence which blocks access to the active-site gorge.

The mannitol operon repressor MtlR belongs to a new class of transcription regulators in bacteria

Many bacteria express phosphoenolpyruvate-dependent phosphotransferase systems (PTS). The mannitol-specific PTS catalyze the uptake and phosphorylation of d-mannitol. The uptake system comprises several genes encoded in the single operon. The expression of the mannitol operon is regulated by a proposed transcriptional factor, mannitol operon repressor (MtlR) that was first studied in Escherichia coli. Here we report the first crystal structures of MtlR from Vibrio parahemeolyticus (Vp-MtlR) and its homolog YggD protein from Shigella flexneri (Sf-YggD). MtlR and YggD belong to the same protein family (Pfam05068). Although Vp-MtlR and Sf-YggD share low sequence identity (22%), their overall structures are very similar, representing a novel all alpha-helical fold, and indicate similar function. However, their lack of any known DNA-binding structural motifs and their unfavorable electrostatic properties imply that MtlR/YggD are unlikely to bind a specific DNA operator directly as proposed earlier. This structural observation is further corroborated by in vitro DNA-binding studies of E. coli MtlR (Ec-MtlR), which detected no interaction of Ec-MtlR with the well characterized mannitol operator/promoter region. Therefore, MtlR/YggD belongs to a new class of transcription factors in bacteria that may regulate gene expression indirectly as a part of a larger transcriptional complex.

Acetylcholinesterase: how is structure related to function?

In accordance with its biological role, termination of neurotransmission at cholinergic synapses by rapid hydrolysis of the neurotransmitter, acetylcholine, acetylcholinesterase is one of nature's most efficient enzymes. Solution of its three-dimensional structure revealed that its active site is located at the bottom of a deep and narrow gorge. Such an architecture was unanticipated in view of its high turnover number. The present review examines how the highly specialized structure of acetylcholinesterase, with its sequestered active site, contributes to its catalytic efficacy, and discusses how the traffic of substrate and products to and from the active site is controlled.

Attractive membrane domains control lateral diffusion

Lipid membranes play a fundamental role in vital cellular functions such as signal transduction. Many of these processes rely on lateral diffusion within the membrane, generally a complex fluid containing ordered microdomains. However, little attention has been paid to the alterations in transport dynamics of a diffusing species caused by long-range interactions with membrane domains. In this paper, we address the effect of such interactions on diffusive transport by studying lateral diffusion in a phase-separated Langmuir phospholipid monolayer via single-particle tracking. We find that attractive dipole-dipole interactions between condensed phase domains and diffusing probe beads lead to transient confinement at the phase boundaries, causing a transition from two- to one-dimensional diffusion. Using Brownian dynamics simulations, the long-term diffusion constant for such a system is found to have a sensitive, Boltzmann-like, dependence on the interaction strength. In addition, this interaction strength is shown to be a strong function of the ratio of domain to particle size. As similar interactions are expected in biological membranes, the modulation of diffusive transport dynamics by varying interaction strength and/or domain size may offer cells selective spatial and temporal control over signaling processes.

A server and database for dipole moments of proteins

An Internet server at http://bip.weizmann.ac.il/dipol calculates the net charge, dipole moment and mean radius of any 3D protein structure or its constituent peptide chains, and displays the dipole vector superimposed on a ribbon backbone of the protein. The server can also display the angle between the dipole and a selected list of amino acid residues in the protein. When the net charges and dipole moments of ∼12 000 non-homologous PDB biological units (PISCES set), and their unique chains of length 50 residues or longer, were examined, the great majority of both charges and dipoles fell into a very narrow range of values, with long extended tails containing a few extreme outliers. In general, there is no obvious relation between a protein's charge or dipole moment and its structure or function, so that its electrostatic properties are highly specific to the particular protein, except that the majority of chains with very large positive charges or dipoles bind to ribosomes or interact with nucleic acids.

Termination and beyond: acetylcholinesterase as a modulator of synaptic transmission

Termination of synaptic transmission by neurotransmitter hydrolysis is a substantial characteristic of cholinergic synapses. This unique termination mechanism makes acetylcholinesterase (AChE), the enzyme in charge of executing acetylcholine breakdown, a key component of cholinergic signaling. AChE is now known to exist not as a single entity, but rather as a combinatorial complex of protein products. The diverse AChE molecular forms are generated by a single gene that produces over ten different transcripts by alternative splicing and alternative promoter choices. These transcripts are translated into six different protein subunits. Mature AChE proteins are found as soluble monomers, amphipatic dimers, or tetramers of these subunits and become associated to the cellular membrane by specialized anchoring molecules or members of other heteromeric structural components. A substantial increasing body of research indicates that AChE functions in the central nervous system go far beyond the termination of synaptic transmission. The non-enzymatic neuromodulatory functions of AChE affect neurite outgrowth and synaptogenesis and play a major role in memory formation and stress responses. The structural homology between AChE and cell adhesion proteins, together with the recently discovered protein partners of AChE, predict the future unraveling of the molecular pathways underlying these multileveled functions.

Investigating the permanent electric dipole moment of β-lactoglobulin fibrils, using transient electric birefringence

Amyloid fibrils, which are polymeric assemblies of protein molecules, have been intensively studied on a structural level, yet due to complications such as the disorder within the molecules, several aspects of their structure remain mysterious. Similarly, the kinetics of assembly are not well understood. Here we investigate the electric dipole moment of beta-lactoglobulin fibrils, a model amyloid fibril system, by applying the technique of transient electric birefringence. This moment appears to be large, and comparable to the total moment of the constituent protein monomers if they were joined in a chain, head-to-tail, without changing conformation, suggesting an ordered joining of monomers in the fibril. Such an ordered assembly may have implications for the assembly mechanism of beta-lactoglobulin fibrils in particular, and amyloid fibrils in general.

Structural insights into substrate traffic and inhibition in acetylcholinesterase

Acetylcholinesterase (AChE) terminates nerve-impulse transmission at cholinergic synapses by rapid hydrolysis of the neurotransmitter, acetylcholine. Substrate traffic in AChE involves at least two binding sites, the catalytic and peripheral anionic sites, which have been suggested to be allosterically related and involved in substrate inhibition. Here, we present the crystal structures of Torpedo californica AChE complexed with the substrate acetylthiocholine, the product thiocholine and a nonhydrolysable substrate analogue. These structures provide a series of static snapshots of the substrate en route to the active site and identify, for the first time, binding of substrate and product at both the peripheral and active sites. Furthermore, they provide structural insight into substrate inhibition in AChE at two different substrate concentrations. Our structural data indicate that substrate inhibition at moderate substrate concentration is due to choline exit being hindered by a substrate molecule bound at the peripheral site. At the higher concentration, substrate inhibition arises from prevention of exit of acetate due to binding of two substrate molecules within the active-site gorge.

Inhibition of human acetyl- and butyrylcholinesterase by novel carbamates of (-)- and (+)-tetrahydrofurobenzofuran and methanobenzodioxepine

A new enantiomeric synthesis utilizing classical resolution provided two novel series of optically active inhibitors of cholinesterase: (-)- and (+)-O-carbamoyl phenols of tetrahydrofurobenzofuran and methanobenzodioxepine. An additional two series of (-)- and (+)-O-carbamoyl phenols of pyrroloindole and furoindole were obtained by known procedures, and their anticholinesterase actions were similarly quantified against freshly prepared human acetyl- (AChE) and butyrylcholinesterase (BChE). Both enantiomeric forms of each series demonstrated potent cholinesterase inhibitory activity (with IC(50) values as low as 10 nM for AChE and 3 nM for BChE), with the exception of the (+)-O-carbamoyl phenols of pyrroloindole, which lacked activity (IC(50) values >1 microM). Based on the biological data of these four series, a structure-activity relationship (SAR) analysis was provided by molecular volume calculations. In addition, a probable transition-state model was established according to the known X-ray structure of a transition-state complex of Torpedo californica AChE-m-(N,N,N-trimethylammonio)-2,2,2-trifluoroacetophenone (TcAChE-TMTFA). This model proved valuable in explaining the enantioselectivity and enzyme subtype selectivity of each series. These carbamates are more potent than, or similarly potent to, anticholinesterases in current clinical use, providing not only inhibitors of potential clinical relevance but also pharmacological tools to define drug-enzyme binding interactions within an enzyme crucial in the maintenance of cognition and numerous systemic physiological functions in health, aging, and disease.

Role of Ca2+ in the electrostatic stability and the functional activity of the globular domain of human C1q

C1q is the recognition subunit of the classical pathway of the complement system and a major connecting link between classical pathway-driven innate immunity and IgG- or IgM-mediated acquired immunity. The basic structural subunit of C1q is composed of an N-terminal triple-helical collagen-like region and a C-terminal heterotrimeric globular head domain (gC1q) that is made up of individual A, B, and C chains. Recent crystallographic studies have revealed that the gC1q domain, which is the main target-binding region of C1q, has a compact and spherical heterotrimeric assembly, held together by both electrostatic and nonpolar interactions, with quasi-3-fold symmetry. A characteristic feature of the gC1q domain is the presence of a exposed Ca(2+) located near the apex. We have investigated, using theoretical and experimental approaches, the role of Ca(2+) in the electrostatic stability and target-binding properties of the native C1q as well as recombinant monomeric forms of the C-terminal regions of the A, B, and C chains. Here, we report that Ca(2+) primarily influences the target recognition properties of C1q toward IgG, IgM, C-reactive protein, and pentraxin 3. At pH 7.4, the loss of Ca(2+) leads to changes in the direction of electric moment from coaxial (where the putative C-reactive protein-binding site is located) to perpendicular to the molecular axis (toward the most likely IgG-binding site), which appears important for target recognition by C1q and subsequent complement activation.

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].

Kinetic analysis of the inhibition of human butyrylcholinesterase with cymserine

Accompanying the gradual rise in the average age of the population of most industrialized countries is a regrettable progressive rise in the number of individuals afflicted with age-related neurodegenerative disorders, epitomized by Alzheimer's disease (AD) but, additionally, including Parkinson's disease (PD) and stroke. The primary therapeutic strategy, to date, involves the use of cholinesterases inhibitors (ChEIs) to amplify residual cholinergic activity. The enzyme, acetylcholinesterase (AChE), along with other elements of the cholinergic system is depleted in the AD brain. In contrast, however, its sister enzyme, butyrylcholinesterase (BuChE), that likewise cleaves acetylcholine (ACh), is elevated and both AChE and BuChE co-localize in high amounts with the classical pathological hallmarks of AD. The mismatch between increased brain BuChE and depleted levels of both ACh and AChE, particularly late in the disease, has supported the design and development of new ChEIs with a preference for BuChE; exemplified by the novel agent, cymserine, whose binding kinetics are characterized for the first time. Specifically, as assessed by the Ellman method, cymserine demonstrated potent concentration-dependent binding with human BuChE. The IC50 was determined as 63 to 100 nM at the substrate concentration range of 25 to 800 microM BuSCh. In addition, the following new binding constants were investigated for human BuChE inhibition by cymserine: T(FPnubeta), K(nubeta), K(Bs), K(MIBA), M(IC50), D(Sc), R(f), (O)K(m), OIC100, K(sl), theta(max) and R(i). These new kinetic constants may open new avenues for the kinetic study of the inhibition of a broad array of other enzymes by a wide variety of inhibitors. In synopsis, cymserine proved to be a potent inhibitor of human BuChE in comparison to its structural analogue, phenserine.

Electric birefringence study of an amyloid fibril system: the short end of the length distribution

In this article, a system of amyloid fibrils, based on the protein beta-lactoglobulin, is studied by transient electric birefringence. Single pulses of an electric field were applied to the solution, and the initial rise and subsequent decay of birefringence analysed. The decay takes place on a range of relaxation times, and therefore contains information about the length distribution of fibrils in the system. The information can be extracted using theories of the electric polarisability of polyelectrolyte rods, since the fibrils are an example of these. Despite the long-standing complications of such theories, useful quantitative information about the system can still be obtained. Using the Fixman model of polyelectrolyte polarisability, we obtain a measurement of the short end of the length distribution which shows the fibril concentration as a function of length rising linearly from 0.02-2 microm. The short end of the length distribution was unobtainable in our previous study using rheo-optics (S.S. Rogers et al., Macromolecules 38, 2948 (2005)), but reasonable agreement between the two techniques shows they are complementary.

Protein electrostriction: a possibility of elastic deformation of the α-hemolysin channel by the applied field

While conformational flexibility of proteins is widely recognized as one of their functionally crucial features and enjoys proper attention for this reason, their elastic properties are rarely discussed. In ion channel studies, where the voltage-induced or ligand-induced conformational transitions, gating, are the leading topic of research, the elastic structural deformation by the applied electric field has never been addressed at all. Here we examine elasticity using a model channel of known crystal structure-Staphylococcus aureus alpha-hemolysin. Working with single channels reconstituted into planar lipid bilayers, we first show that their ionic conductance is asymmetric with voltage even at the highest salt concentration used where the static charges in the channel interior are maximally shielded. Second, choosing 18-crown-6 as a molecular probe whose size is close to the size of the narrowest part of the alpha-hemolysin pore, we analyze the blockage of the channel by the crown/K(+) complex. Analysis of the blockage within the framework of the Woodhull model in its generalized form demonstrates that the model is able to correctly describe the crown effect only if the parameters of the model are considered to be voltage-dependent. Specifically, one has to include either a voltage-dependent barrier for crown release to the cis side of the channel or voltage-dependent interactions between the binding site and the crown. We suggest that the voltage sensitivity of both the ionic conductance of the channel seen at the highest salt concentration and its blockage by the crown reflects a field-induced deformation of the pore.

Reversal of some viral IL-6 electrostatic properties compared to IL-6 contributes to a loss of alpha receptor component recruitment

Human interleukin-6 (hIL-6) is a pleiotropic mediator of activation and proliferation across a large number of different cell types. Human herpesvirus-8 (HHV-8) has been associated with classical and AIDS-related Kaposi's sarcoma (KS). HHV-8 encodes viral IL-6 (vIL-6), a functional homolog of human interleukin-6, that promotes the growth of KS and of some lymphoma cells. Signaling induced by human IL-6 requires recruitment of the glycoprotein gp130, which acts as the signal transducing chain, and of IL-6Ralpha, which is necessary for cognate recognition and high affinity receptor complex formation. In contrast, the formation of a functional complex between vIL-6 and gp130 does not require the presence of IL-6Ralpha. The physico-chemical properties of vIL-6 have been analyzed and compared to those of hIL-6 and of the receptor chains, gp130 and IL-6Ralpha. Interaction sites on vIL-6 involve more hydrophobic residues than those of hIL-6. The electrostatic fields induced by vIL-6 and IL-6Ralpha are repulsive and prevent interaction between vIL-6 and IL-6Ralpha, whereas the electrostatic field induced by hIL-6 steers the complex formation with IL-6Ralpha. Subsequently, electrostatic binding free energy in the vIL-6/IL-6Ralpha complex is destabilizing, whereas it is stabilizing in the complex comprising hIL-6. These properties result from charge reversals between viral and human IL-6, an unusual phenomenon of amino acid substitutions within a homologous protein family. This suggests a selection pressure for vIL-6 to by-pass the IL-6Ralpha control of host defense against virus infection. This selection pressure has yielded the reversal of electrostatic properties of vIL-6 when compared to hIL-6.

Structural and kinetic effects of mobile phone microwaves on acetylcholinesterase activity

The present study provides evidence that "in vitro" simple exposure of an aqueous solution of electric eel acetylcholinesterase (EeAChE; EC 3.1.1.7.) to cellular phone emission alters its enzymatic activity. This paper demonstrates, by combining different experimental techniques, that radio frequency (RF) radiations irreversibly affect the structural and biochemical characteristics of an important CNS enzyme. These results were obtained by using a commercial cellular phone to reproduce the reality of the human exposition. This experimental procedure provided surprising effects collected practically without experimental errors because they were obtained comparing native and irradiated sample of the same enzyme solution. Although these results cannot be used to conclude whether exposure to RF during the use of cellular phone can lead to any hazardous health effect, they may be a significant first step towards further verification of these effects on other "ex vivo" or "in vivo" biological systems.

The Complex of a Bivalent Derivative of Galanthamine withTorpedoAcetylcholinesterase Displays Drastic Deformation of the Active-Site Gorge: Implications for Structure-Based Drug Design

Bifunctional derivatives of the alkaloid galanthamine, designed to interact with both the active site of the enzyme acetylcholinesterase (AChE) and its peripheral cation binding site, have been assayed with Torpedo californica AChE (TcAChE), and the three-dimensional structures of their complexes with the enzyme have been solved by X-ray crystallography. Differences were noted between the IC(50) values obtained for TcAChE and those for Electrophorus electricus AChE. These differences are ascribed to sequence differences in one or two residues lining the active-site gorge of the enzyme. The binding of one of the inhibitors disrupts the native conformation of one wall of the gorge, formed by the loop Trp279-Phe290. It is proposed that flexibility of this loop may permit the binding of inhibitors such as galanthamine, which are too bulky to penetrate the narrow neck of the gorge formed by Tyr121 and Phe330 as seen in the crystal structure.

Identification of a structural site on acetylcholinesterase that promotes neurite outgrowth and binds laminin-1 and collagen IV

The cell adhesion and neurite outgrowth-promoting function of acetylcholinesterase has been localised to the area of the peripheral anionic site. In order to precisely determine the site involved, we used synthetic peptides representing sequences of the peripheral anionic site and its surrounds, and investigated their binding to a panel of monoclonal antibodies that inhibit cell adhesion/neurite outgrowth and/or to recognise the peripheral anionic site. Binding to laminin-1 and collagen IV was also investigated. A relationship between recognition of the sequence 37-50, representing a surface loop adjacent to the peripheral anionic site, and the degree of inhibition of cell adhesion was observed; both laminin-1 and collagen IV also bound this loop with high affinity. Neurite outgrowth on coverslips coated with this peptide was similar to those coated with acetylcholinesterase itself. Adhesion-inhibiting antibodies also recognised the omega loop 69-96, as did laminin-1 and collagen IV. Laminin also bound the sequences 55-66 and 340-353, recognised by the antibodies to varying degrees, but collagen did not. All these peptides were able to promote neurite outgrowth to some degree. No binding to the amyloid-binding omega loop 275-304 by the ligands was observed, nor did the antibodies recognise this consistently. No relationship was observed between the degree of inhibition of acetylcholinesterase and inhibition of neurite outgrowth by the antibodies from which we conclude that the neurite outgrowth function is non-cholinergic. In conclusion, we have identified a specific conformational structure on acetylcholinesterase, comprising adjacent surface loops between residues 37-50 and 69-96, with additional involvement of the sequences 55-66 and 340-353, that mediates cell adhesion and neurite outgrowth.

A paradigm for single nucleotide polymorphism analysis: The case of the acetylcholinesterase gene

Acetylcholinesterase (AChE) plays a crucial physiological role in termination of impulse transmission at cholinergic synapses through rapid hydrolysis of acetylcholine. It is a highly conserved molecule, and only a few naturally occurring genetic polymorphisms have been reported in the human gene. The goal of the present study was to make a systematic effort to identify natural single nucleotide polymorphisms (SNPs) in the human ACHE gene. To this end, the genomic coding sequences for acetylcholinesterase of 96 unrelated control individuals from three distinct ethnic groups were analyzed. A total of 13 ACHE SNPs were identified, 10 of which are newly described, and five that should produce amino acid substitutions [c.101G>A (p.Arg34Gln), c.169G>A (p.Gly57Arg), c.1031A>G (p.Glu344Gly), c.1057C>A (p.His353Asn), and c.1775C>G (p.Pro592Arg)]. Population frequencies of 11 of the 13 SNPs were established in four different populations: African Americans, Ashkenazi Jews, Sephardic Jews, and Israeli Arabs; 15 haplotypes and five ethnospecific alleles were identified. The low number of SNPs identified until now in the ACHE gene is ascribed to technical hurdles arising from the high GC content and the presence of numerous repeat sequences, and does not reflect its intrinsic heterozygosity. Among the SNPs resulting in an amino acid substitution, three are within the mature protein, mapping on its external surface: they are thus unlikely to affect its catalytic properties, yet could have antigenic consequences or affect putative protein-protein interactions. Furthermore, the newly identified SNPs open the door to a study of the possible association of AChE with deleterious phenotypes-such as adverse drug responses to AChE inhibitors employed in treatment of Alzheimer patients and hypersensitivity to pesticides.

A molecular mechanism of enantiorecognition of tertiary alcohols by carboxylesterases

Carboxylesterases containing the sequence motif GGGX catalyze the hydrolysis of esters of chiral tertiary alcohols, albeit with only low to moderate enantioselectivity, for three model substrates (linalyl acetate, methyl-1-pentin-1-yl acetate, 2-phenyl-3-butin-2-yl acetate). In order to understand the molecular mechanism of enantiorecognition and to improve enantioselectivity for this interesting substrate class, the interaction of both enantiomers with the substrate binding sites of acetylcholinesterases and p-nitrobenzyl esterase from Bacillus subtilis was modeled and correlated to experimental enantioselectivity. For all substrate-enzyme pairs, enantiopreference and ranking by enantioselectivity could be predicted by the model. In p-nitrobenzyl esterase, one of the key residues in determining enantioselectivity was G105: exchange of this amino acid for an alanine residue led to a sixfold increase of enantioselectivity (E = 19) towards 2-phenyl-3-butin-2-yl acetate. However, the effect of this mutation is specific: the same mutant had the opposite enantiopreference towards the substrate linalyl acetate. Thus, depending on the substrate structure, the same mutant has either increased enantioselectivity or opposite enantiopreference compared to the wild-type enzyme.

A simple electrostatic criterion for predicting the thermal stability of proteins

The enhancement of protein thermostability is an important issue for both basic science and biotechnology purposes. We have developed a thermostability criterion for a protein in terms of a quasi-electric dipole moment (contributed by its charged residues) defined for an electric charge distribution whose total charge is not zero. It was found that minimization of the modulus of this dipole moment increased its thermal stability, as demonstrated by surveying these values in pairs of mesostable-thermostable homologous proteins and in mutations described in the literature. The analysis of these observations provides criteria for thermostabilization of a protein, by computing its dipole profile. This profile is obtained by direct substitution of each amino acid of the sequence by either a positive, negative or neutral amino acid, followed by a recalculation of the dipole moment. As an experimental example, these criteria were applied to a beta-glucanase to enhance its thermal stability.

Acetylcholinesterase in motion: visualizing conformational changes in crystal structures by a morphing procedure

In order to visualize and appreciate conformational changes between homologous three-dimensional (3D) protein structures or protein/inhibitor complexes, we have developed a user-friendly morphing procedure. It enabled us to detect coordinated conformational changes not easily discernible by analytic methods or by comparison of static images. This procedure was applied to comparison of native Torpedo californica acetylcholinesterase and of complexes with reversible inhibitors and conjugates with covalent inhibitors. It was likewise shown to be valuable for the visualization of conformational differences between acetylcholinesterases from different species. The procedure involves generation, in Cartesian space, of 25 interpolated intermediate structures between the initial and final 3D structures, which then serve as the individual frames in a QuickTime movie.

Acetylcholinesterase complexed with bivalent ligands related to huperzine a: experimental evidence for species-dependent protein-ligand complementarity

Acetylcholinesterase (AChE) inhibitors improve the cognitive abilities of Alzheimer patients. (-)-Huperzine A [(-)-HupA], an alkaloid isolated from the club moss, Huperzia serrata, is one such inhibitor, but the search for more potent and selective drugs continues. Recently, alkylene-linked dimers of 5-amino-5,6,7,8-tetrahydroquinolinone (hupyridone, 1a), a fragment of HupA, were shown to serve as more potent inhibitors of AChE than (-)-HupA and monomeric 1a. We soaked two such dimers, (S,S)-(-)-bis(10)-hupyridone [(S,S)-(-)-2a] and (S,S)-(-)-bis(12)-hupyridone [(S,S)-(-)-2b] containing, respectively, 10 and 12 methylenes in the spacer, into trigonal TcAChE crystals, and solved the X-ray structures of the resulting complexes using the difference Fourier technique, both to 2.15 A resolution. The structures revealed one HupA-like 1a unit bound to the "anionic" subsite of the active-site, near the bottom of the active-site gorge, adjacent to Trp84, as seen for the TcAChE/(-)-HupA complex, and the second 1a unit near Trp279 in the "peripheral" anionic site at the top of the gorge, both bivalent molecules thus spanning the active-site gorge. The results confirm that the increased affinity of the dimeric HupA analogues for AChE is conferred by binding to the two "anionic" sites of the enzyme. Inhibition data show that (-)-2a binds to TcAChE approximately 6-7- and > 170-fold more tightly than (-)-2b and (-)-HupA, respectively. In contrast, previous data for rat AChE show that (-)-2b binds approximately 3- and approximately 2-fold more tightly than (-)-2a and (-)-HupA, respectively. Structural comparison of TcAChE with rat AChE, as represented by the closely related mouse AChE structure (1maa.pdb), reveals a narrower gorge for rat AChE, a perpendicular alignment of the Tyr337 ring to the gorge axis, and its conformational rigidity, as a result of hydrogen bonding between its hydroxyl group and that of Tyr341, relative to TcAChE Phe330. These structural differences in the active-site gorge explain the switch in inhibitory potency of (-)-2a and 2b and the larger dimer/(-)-HupA potency ratios observed for TcAChE relative to rat AChE. The results offer new insights into factors affecting protein-ligand complementarity within the gorge and should assist the further development of improved AChE inhibitors.

A structure-based design approach to the development of novel, reversible AChE inhibitors

Chimeras of tacrine and m-(N,N,N-Trimethylammonio)trifluoroacetophenone (1) were designed as novel, reversible inhibitors of acetylcholinesterase. On the basis of the X-ray structure of the apoenzyme, a molecular modeling study determined the favored attachment positions on the 4-aminoquinoline ring (position 3 and the 4-amino nitrogen) and the favored lengths of a polymethylene link between the two moieties (respectively 5-6 and 4-5 sp(3) atoms). Seven compounds matching these criteria were synthesized, and their inhibitory potencies were determined to be in the low nanomolar range. Activity data for close analogues lacking some of the postulated key features showed that our predictions were correct. In addition, a subsequent crystal structure of acetylcholinesterase complexed with the most active compound 27 was in good agreement with our model. The design strategy is therefore validated and can now be developed further.

Predicted Michaelis-Menten complexes of cocaine-butyrylcholinesterase. Engineering effective butyrylcholinesterase mutants for cocaine detoxication

Butyrylcholinesterase (BChE) is important in cocaine metabolism, but it hydrolyzes (-)-cocaine only one-two thousandth as fast as the unnatural (+)-stereoisomer. A starting point in engineering BChE mutants that rapidly clear cocaine from the bloodstream, for overdose treatment, is to elucidate structural factors underlying the stereochemical difference in catalysis. Here, we report two three-dimensional Michaelis-Menten complexes of BChE liganded with natural and unnatural cocaine molecules, respectively, that were derived from molecular modeling and supported by experimental studies. Such complexes revealed that the benzoic ester group of both cocaine stereoisomers must rotate toward the catalytic Ser(198) for hydrolysis. Rotation of (-)-cocaine appears to be hindered by interactions of its phenyl ring with Phe(329) and Trp(430). These interactions do not occur with (+)-cocaine. Because the rate of (-)-cocaine hydrolysis is predicted to be determined mainly by the re-orientation step, it should not be greatly influenced by pH. In fact, measured rates of this reaction were nearly constant over the pH range from 5.5 to 8.5, despite large rate changes in hydrolysis of (+)-cocaine. Our models can explain why BChE hydrolyzes (+)-cocaine faster than (-)-cocaine, and they suggest that mutations of certain residues in the catalytic site could greatly improve catalytic efficiency and the potential for detoxication.

Direct evidence for an adhesive function in the noncholinergic role of acetylcholinesterase in neurite outgrowth

Acetylcholinesterase (AChE) can promote neurite outgrowth through a mechanism that is independent of its role in hydrolyzing the neurotransmitter acetylcholine. It has been proposed that this neuritogenic capacity of AChE may result from its intrinsic capacity to function in adhesion. In this study we directly tested this hypothesis using neuroblastoma cell lines that have been engineered for altered cell-surface expression of AChE. Using a microtiter-plate adhesion assay and the electrical cell-substrate impedance-sensing (ECIS) method, we demonstrate that the level of cell-substratum adhesion of these cells directly correlates with their level of AChE expression. Furthermore, this adhesion is blocked by either an anti-AChE antibody or a highly specific AChE inhibitor (BW284c51), both of which have also been shown to block neurite outgrowth. In addition, cells that overexpress AChE showed enhanced neurite initiation. By employing cell lines with different levels of AChE expression in two types of cell-substratum adhesion assays, our current studies provide evidence for an adhesive function for AChE. These results, together with the fact that AChE shares sequence homology and structural similarities with several known cell adhesion molecules, support the hypothesis that AChE promotes neurite outgrowth, at least in part, through an adhesive function.

Three-dimensional structures of Drosophila melanogaster acetylcholinesterase and of its complexes with two potent inhibitors

We have crystallized Drosophila melanogaster acetylcholinesterase and solved the structure of the native enzyme and of its complexes with two potent reversible inhibitors, 1,2,3,4-tetrahydro-N-(phenylmethyl)-9-acridinamine and 1,2,3,4-tetrahydro-N-(3-iodophenyl-methyl)-9-acridinamine--all three at 2.7 A resolution. The refined structure of D. melanogaster acetylcholinesterase is similar to that of vertebrate acetylcholinesterases, for example, human, mouse, and fish, in its overall fold, charge distribution, and deep active-site gorge, but some of the surface loops deviate by up to 8 A from their position in the vertebrate structures, and the C-terminal helix is shifted substantially. The active-site gorge of the insect enzyme is significantly narrower than that of Torpedo californica AChE, and its trajectory is shifted several angstroms. The volume of the lower part of the gorge of the insect enzyme is approximately 50% of that of the vertebrate enzyme. Upon binding of either of the two inhibitors, nine aromatic side chains within the active-site gorge change their conformation so as to interact with the inhibitors. Some differences in activity and specificity between the insect and vertebrate enzymes can be explained by comparison of their three-dimensional structures.

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.

Evidence for the direct role of acetylcholinesterase in neurite outgrowth in primary dorsal root ganglion neurons

Dorsal root ganglion (DRG) neurons show a transient peak expression of acetylcholinesterase (AChE) during periods of axonal outgrowth prior to synaptogenesis, suggesting that AChE has a non-enzymatic role during development. We have previously shown that perturbation of cell surface AChE in cultured embryonic rat DRG neurons results in decreased neurite outgrowth and neurite detachment. In this report, we demonstrate a direct correlation between endogenous AChE content and neurite outgrowth in primary DRG neurons. Adenoviral vectors were constructed using full-length rat AChE(T) cDNA in either the sense or antisense orientations to overexpress or knock down AChE expression, respectively. Treatment with the sense-expressing vector produced a 2.5-fold increase in AChE expression and a 2-fold increase in neurite length compared with either untreated or null virus-treated control cells. Conversely, treatment with the antisense-expressing vector reduced AChE expression by 40% and resulted in a reduction in neurite length of similar magnitude. We also observed that overexpression of AChE resulted in greater branching at the distal tips of each primary neurite as well as an increase in cell body size. These findings further indicate that AChE expressed on the axonal surface of developing DRG neurons may modulate their adhesive properties and thereby support axonal development.

Active-site gorge and buried water molecules in crystal structures of acetylcholinesterase from Torpedo californica

Buried water molecules and the water molecules in the active-site gorge are analyzed for five crystal structures of acetylcholinesterase from Torpedo californica in the resolution range 2.2-2.5 A (native enzyme, and four inhibitor complexes). A total of 45 buried hydration sites are identified, which are populated with between 36 and 41 water molecules. About half of the buried water is located in a distinct region neighboring the active-site gorge. Most of the buried water molecules are very well conserved among the five structures, and have low displacement parameters, B, of magnitudes similar to those of the main-chain atoms of the central beta-sheet structure. The active-site gorge of the native enzyme is filled with over 20 water molecules, which have poor hydrogen-bond coordination with an average of 2.9 polar contacts per water molecule. Upon ligand binding, distinct groups of these water molecules are displaced, whereas the others remain in positions similar to those that they occupy in the native enzyme. Possible roles of the buried water molecules are discussed, including their possible action as a lubricant to allow large-amplitude fluctuations of the loop structures forming the gorge wall. Such fluctuations are required to facilitate traffic of substrate, products and water molecules to and from the active-site. Because of their poor coordination, the gorge water molecules can be considered as "activated" as compared to bulk water. This should allow their easy displacement by incoming substrate. The relatively loose packing of the gorge water molecules leaves numerous small voids, and more efficient space-filling by substrates and inhibitors may be a major driving force of ligand binding.