Encyclopedia
The cholinesterase enzyme family has but two members: acetylcholinesterase (AChE) and butyrylcholinesterase (BChE). Computational efforts to better understand these enzymes, with a focus on their structural and catalytic properties and dynamics, includes docking and other Monte Carlo based calculations, as well as dynamic simulations of varying rigor and resource-demand including molecular dynamics, Langevin dynamics, quantum mechanical/molecular mechanics, and related algorithms.
- acetylcholinesterase
- butyrylcholinesterase
- docking
- molecular dynamics
- hydrolysis
- molecular recognition
1. Introduction
The cholinesterase enzyme family has but two members: acetylcholinesterase (AChE) and butyrylcholinesterase (BChE). The former’s primary biological purpose is regulating acetylcholine, a neurotransmitter, via hydrolysis at neuromuscular junctions, thus proving itself to be an essential component in the maintenance and performance of nervous systems. AChE, which has also been referred to as “true cholinesterase”, is created in muscle, nerve, and hematopoietic cells and is considered to be one of the most efficient enzymes due to its rapid rate of catalysis [1]. AChE’s plasma analog, BChE, previously referred to as “pseudocholinesterase”, is produced in the liver and, unlike AChE, has been viewed as having a much more ambiguous biological purpose, as it was long believed to be vestigial [2]. BChE has a higher concentration in plasma than AChE, is present in many vertebrates, and tolerates several mutations, which has prompted the theory that BChE evolved from AChE to be a general detoxifier while still retaining some function in the process of neurotransmission [3]. This theory seems apt given the structural similarity of their binding sites, their shared affinities for certain substrates and ligands, and their sequence homology of approximately 65% [4]. Figure 1 depicts both AChE (bottom) and BChE (top) from a gorge-centric view and after a 90° rotation.
Figure 1. Visualizations of (a) butyrylcholinesterase (BChE) (PDBID 1P0I) and (b) acetylcholinesterase (AChE) (PDBID 1B41) in white ribbon mode with residues in notable binding sites shown as semi-transparent van der Waals surfaces colored according to the key. Top panels present views directly into the active site gorges with structures in the center panel rotated 90° about the vertical axis. Bottom panels present magnified views of the binding site regions of (c) BChE and (d) AChE with key residues labeled.
One difference between the two enzymes is the relative size of the binding pocket, with BChE and AChE having approximate gorge volumes of 1500 Å3 and 1300 Å3, respectively [5]. X-ray crystallography of Torpedo californica (Tc, pacific ray) AChE revealed the enzyme to have a deep hydrophobic gorge with residues that stabilize substrates in the pocket [6], as well as a bottleneck region in the active site [7] that narrows to approximately 4 Å in width [8]. Common models of AChE, including human, mouse, and Torpedo californica (TcAChE), demonstrate conserved active sites, save for a few residues that participate in ligand binding [9], and both have negative surface potentials that become more negative deeper within the gorge. This negative potential, which is high near the catalytic site at the “bottom” of the gorge, seems to have evolved to facilitate electrostatic attraction of positively charged choline substrates by both enzymes [8].
Although BChE is structurally similar to AChE, 6 of the 14 aromatic amino acids that line the active site gorge in AChE are substituted with aliphatic residues in BChE [10]. In particular, the substitution of Phe288 and Phe290 in TcAChE with the smaller Leu286 and Val288 in human BChE lead to conformational changes that result in a deeper gorge in BChE, thereby allowing BChE to interact with, and potentially hydrolize, a much wider range of substrates and inhibitors than AChE [11]. BChE is thus characterized as the promiscuous, or non-specific, bigger sibling to the smaller and much more specific AChE.
The tremendous growth and improvement in computational resources and modeling techniques over the past few decades have led to an exponential increase in computational studies of biomolecular systems. In 2003, the Protein Data Bank (PDB) launched its online presence, thereby making myriad biomolecules and complexes available for structural and computational research, and by 2016 there were 178 AChE structures available to the public, many of those including bound substrates or inhibitors of medical and pharmacological significance [9]. In fact, some of the earliest and most significant contributions of cholinesterase models were reported by Sussman et al., who performed an X-ray analysis of AChE at 2.8 resolution in 1991 [12], and Nicolet et al., who examined several crystal structures of BChE in 2003 [13]. Modeling and simulation-based studies of the cholinesterase enzymes have advanced in tandem with this growth.
The application of molecular dynamics (MD) simulations, Monte Carlo (MC) based docking calculations, and more sophisticated quantum mechanical/molecular mechanics (QM/MM) simulations have proven highly insightful in probing cholinesterase structure and activity. Notable areas of interest include the mechanisms of cholinesterase catalysis, reversible and irreversible inhibition of both enzymes to manage Alzheimer’s Disease (AD) and other ailments, BChE-specific inhibition, and the reactivation of phosphorylated cholinesterases, a process that normally follows nerve agent attacks or pesticide poisoning [14]. Figure 2a,b depict the common paths of substrate or ligand binding and phosphorylation, respectively.
Figure 2. (a) Depiction of the enzyme catalytic mechanism (black) and reversible inhibition mechanisms including competitive inhibition (blue), uncompetitive inhibition (red), and noncompetitive inhibition (both blue and red). (b) Representation of irreversible (covalent) inhibition by organophosphorus ligands and the reactivation mechanisms to revert the aged enzyme (blue) back to the inhibited complex (green). L and R represent the first and second leaving groups of the organophosphate, respectively, following the paths of inhibition and reactivation.
Most therapeutic treatments of AD, and other illnesses to which the cholinesterases have been linked, are reversible inhibitors that form non-covalent molecular recognition (MR) complexes with cholinesterases and can leave the active site, thus existing in an equilibrium between bound and unbound states characterized by KI and/or IC50 values. The three mechanisms that reversible inhibition can follow are competitive, noncompetitive, and uncompetitive, as depicted in Figure 2a. During competitive inhibition, the substrate and inhibitor are competing for the same binding site, and the binding of a competitive inhibitor within the active site blocks entrance of the substrate, thereby hindering formation of the enzyme substrate complex. By contrast, uncompetitive inhibitors function by binding to the enzyme-substrate complex to prevent product formation, and noncompetitive inhibitors can bind to either the enzyme or the enzyme-substrate complex to regulate catalytic activity.
Many organic compounds, including organophosphates (OPs), are categorized as irreversible inhibitors, which covalently bond to residues in the gorge and thus cannot leave the active site. In the cholinesterases, OPs covalently bind to the serine residue of the catalytic triad. Irreversibly inhibited cholinesterases, however, can be reactivated via various pathways with compounds such as oximes, as illustrated in Figure 2b. Oximes reactivate an inhibited cholinesterase to its native structure by nucleophilic substitution of the phosphorylated serine. If untreated, the inhibited enzyme can undergo the dealkylation process, also known as aging, where the loss of the second leaving group produces an oxyanion on the phosphoryl group.
Aged cholinesterase is highly stable due to the strong electrostatic interactions between the oxyanion and the positively charged catalytic histidine. Known reactivators of OP-inhibited AChE, such as oximes, are ineffectual against aged AChE [15]. Though a successful attempt to reactivate aged AChE by a class of compounds called “quinone methide precursors” (QMPs) was reported by Hadad and coworkers [16], the reaction between QMPs and aged AChE was rendered too slow to be useful [17].
2. Structure and Dynamics
While understanding the structure, function, and behavior of cholinesterase binding sites seems a practical starting point for modeling studies, early researchers initially wanted to understand how substrates and inhibitors entered the active sites of these enzymes. Molecular dynamics (MD) simulations of human BChE lasting for 5 and 10 ns indicated that inhibitors can access the binding subsites in the catalytic cavity due to the highly flexible entrance (or “mouth”) of the gorge, a portion of which is formed by the flexible omega loop (Ω-loop) region, as well as the peripheral aromatic site (PAS, formerly mistakenly called the peripheral anionic site). In these simulations, the Asp70 residue in PAS showed significant deviation from the crystal structure with root-mean-square deviation (RMSD) values of 2 to 6 Å [18]. AChE has also been shown to experience such fluctuations at the gorge entrance, with a similarly flexible Ω-loop region that is thought to increase enzyme specificity by making it more difficult for large molecules to enter the AChE gorge without hindering the productivity of the enzyme [19].
Apart from the fluctuations at the mouth of the gorge, AChE also experiences what is known as bottleneck fluctuations, with the bottleneck being the narrowest part of the gorge. These fluctuations, which have come to be known as the “breathing” of the enzyme, can help inhibitors or substrates move from the surface to deeper regions of the gorge [20]. Cheng et al. recently defined this breathing by monitoring the varying distance between the Cε2 atom of Phe330 (CBS) and the OH atom of Tyr121 (PAS) in TcAChE, which suggested that a number of subdomains within the enzyme, particularly the Ω-loop, contribute to modulating the size of the gorge bottleneck [21]. A comprehensive comparison between 47 crystal structures of AChE (in its native form and in complex with small molecules), as well as a 20 ns simulation of TcAChE, provided by Sussman and coworkers, suggested that the 14 aromatic residues lining the AChE gorge and creating over half of the gorge surface area contribute greatly to the overall flexibility of the enzyme, the observed bottleneck breathing motions, and the resulting ability to perform its catalytic function [22].
Some of these aromatic residues play significant roles in primary subsites within the cholinesterase gorges, such as the catalytic active site (CAS) and the peripheral aromatic site (PAS) [22]. While these sites, which are addressed in more detail below, are integral to cholinergic activity, cholinesterases do not only perform cholinergic functions. As discussed by Chinnadurai et al., the aryl acylamidase activity (AAA) of AChE, which also involves hydrolysis, only requires the CAS and does not interact with the PAS at all [23]. This has prompted the theory that AAA substrates enter from a side-door into the enzyme, rather than via the mouth of the gorge and, indeed, researchers have suggested that there are a number of doors through which substrates can enter the gorge including a back door [24], an Ω-loop door, and the suggested side door [23]. It was suggested that the side door may open more frequently than the other doors to mediate AAA activity, and MD simulations of side door probing emphasize the importance of hydrophobic interactions, hydrogen bonding, and water mediated interactions (“water bridges”) in moving the substrate towards the CAS [23]. To be sure, simulations of AChE in explicit solvent sans substrate [25], as well as analyses of TcAChE crystal structures in its native and several inhibited forms [26], have underscored the importance of the presence of molecular water in enzyme structure.
At higher concentrations, dimerization and the further dimerization of dimers to form tetramers is known to affect the structure and function of cholinesterase enzymes [27]. While MD simulations suggested that two of the four binding sites in tetramerized cholinesterases are sterically blocked, thereby becoming less active [28], as reflected by a 15% decrease in catalytic activity [29], a recently elucidated CryoEM structure of hBChE shows distinct structural variance from the simulated tetramer, with the active site gorges being fully solvent accessible [27]. Tetramerization, however, increases the half-life of the enzymes, which is a desirable result when cholinesterases, and BChE in particular, are being used to counteract drug overdoses. For example, the addition of proline-rich attachment domains (PRADs) to BChE increases tetramer stability, leading to an extended circulation time [30]. Interestingly, glycosylated models of BChE increase the enzyme’s flexibility and half-life without hindering its ability to bind to glycans, which cannot be said of all therapeutic protein targets [31]. On the other hand, cholinesterase phosphonylation, or the irreversible binding of organophosphate to the active site, which will be discussed in more detail in the organophosphate inhibition section below, severely restricts the flexibility of both AChE and BChE, as reported by Bennion et al. [32]. Experimental findings of AChE covalently bound to the nerve agent soman agree that OP-poisoned AChE is significantly stiffer [33].
2.1. Important Binding Sites
There has been substantial past effort to study the sites responsible for molecular recognition (MR) and binding affinity within the active site gorge of both enzymes, and it is clear that those binding sites, specific chemical subsites within each active site gorge, mirror each other and perform similar functions respective to each enzyme. While it is important to note that these binding sites have been examined experimentally with X-ray and kinetic studies, including a recent study by Rosenberry et al. [34], the present review focuses on the unique perspective provided by computational investigations. As expected, one of the most important sites for both cholinesterases is the catalytic active site (CAS), and the peripheral aromatic site (PAS) also plays an indispensable role in cholinesterase or ligand binding, while the Ω-loop (OML), acyl binding site (ABS), and oxyanion hole (OAH) sites are more essential in contributing to binding affinity and complex stability. Alvarado et al. have provided a method of succinct graphical tabulation of BChE-ligand contacts and interactions, referred to as contact tables, that include these five sites and additional protein residues of interest [35], as discussed below. A detailed analysis of these binding sites is provided here in the same order that they are encountered by substrates and inhibitors upon entering and moving into the gorge.
2.1.1. Peripheral Aromatic Site
The peripheral aromatic site (PAS, red in Figure 1) is located near the mouth of gorge [36] and plays a prominent role in substrate and ligand binding [37]. For decades, peer-reviewed studies have used PAS to denote the “peripheral anionic site”. In recent years, however, the aromatic properties of this binding site that are vital to cholinesterase function have driven the community to instead refer to this region as the “peripheral aromatic site”. Important amino acids in the PAS of AChE include serine, tyrosine, aspartic acid, and tryptophan [36,38], while notable PAS residues in BChE include asparagine, aspartic acid, glutamine, serine, and tyrosine [35], highlighting the polar, negatively charged, and electron-rich nature of residues in this site. As previously mentioned, one distinction between the cholinesterases is the aromatic nature of the residues surrounding the PAS in AChE, which is more aliphatic in BChE [10]. The PAS makes contact with many loops and secondary structural elements at the surface of the protein, including the Ω-loop, which contributes to the needed flexibility discussed above. Although steric and electrostatic interactions may slow the catalytic efficiency of AChE, the PAS is valuable for trafficking ligands into the gorge [39], particularly positively charged species such as cholines. MD simulations have emphasized the importance of cation-π interactions, which stabilize the ligand at the rim of the gorge entrance prior to entering the gorge [40], and it has been proposed that non-cholinergic activity of the PAS could include the deposition of amyloids, adhesion to cells, and outgrowth of neurites [41].
2.1.2. Acyl and Choline Binding Sites
Once a ligand has entered the gorge, the acyl and choline binding sites (ABS and CBS, shown as blue and green in Figure 1, respectively), which are located near the catalytic triad, assist in positioning the ligand for catalysis. The ABS and CBS are hydrophobic regions composed primarily of tryptophan, tyrosine, and phenylalanine in human AChE. Tyr337 in the choline binding site of AChE is replaced by Ala328 in that of BChE; Phe295 and Phe297 in the acyl binding site of AChE are replaced by Leu286 and Val288, respectively, in BChE. The replacement of aromatic residues in the ABS and CBS of BChE enable it to bind larger substrates and inhibitors than AChE [42]. In addition, the ABS and CBS are largely responsible for the specificity of these enzymes and are thus primary targets studied when synthesizing inhibitors such as imidazole or pyridine derivatives [43].
2.1.3. Catalytic Active Site
The catalytic active site (CAS, yellow in Figure 1) is located approximately 20 Å deep at the bottom of both the AChE and BChE gorges [12,24] and is made up of serine, glutamic acid, and histidine residues, prompting the name “catalytic triad” [24,35]. The CAS is surrounded by numerous aromatic and acidic residues [44] and is observed to engage in shorter hydrogen bonds in crystal and NMR structures than observed in simulation [45,46]. More importantly, the CAS is highly conserved [44], emphasizing historically vital biological roles of these enzymes and their cholinergic activities. QM/MM simulations at the MP2(6-31 + G*) level reveal a potential energy barrier of 10.5 kcal/mol, which agrees with experimental data [47], and MD simulations of AChE bound to acetylcholine (ACh) show that ACh stabilizes the CAS and improves the binding ability of the peripheral aromatic site [48].
It has been postulated that a back-door exists in AChE, just behind the CAS and controlled by a tryptophan residue, which was theorized after a single water molecule exited the active site gorge from a direction opposite that of the gorge entrance in an MD simulation [49]. This “back door” was later thought to open three to four Å wide such that catalysis products could exit the gorge of the enzyme without blocking the gorge entrance, and thus contributing to a high catalytic rate [50]. Aromatic residues surrounding the CAS histidine also largely influence the productivity and efficiency of the enzyme [51], which decreased approximately 600-fold when disrupted or replaced by aliphatic side chains [52]. More recently, Xu et al. used MD to study TcAChE and, from 27 of their 40 trajectories, observed thiocholine to frequently exit the active site via a back-door created by cooperative motions of CAS residue Trp86 along with Val132 and Gly448 [53], for which previous experimental support was noted [54,55].
Indeed, other mutations in or near the CAS are known to have effects on the structure, and subsequent function, of the enzyme. This research has naturally focused on, and is more applicable to, BChE due to its much greater natural affinity for mutations than AChE [3]. For example, prolonged use of muscle relaxers led to the discovery of the “silent phenotype” in which an alanine is mutated to a valine near the CAS of BChE. This mutation was studied in silico and observed to severely disrupt interactions between the histidine and serine in the catalytic triad [56], leading to a dysfunctional CAS, regardless of the inhibitor, and increases in the volume of the enzyme, indicating that this may be a pre-denaturation state [57].
The mutation of a nearby alanine in silico to cystine in wild-type (WT) BChE causes the histidine in the catalytic triad to flip, an event that is largely guided by local water molecules [58]. This man-made mutation, while possibly slowing the speed of binding, ultimately still allows for substrate binding to the active site; the naturally occurring mutation of that alanine to aspartic acid, however, is claimed to be catalytically inactive due to strong disruptive interactions between aspartic acid and the CAS histidine residue [59]. Both mutations showcase the possible hysteretic behavior of BChE, or its reliance on past-states, which can likely be attributed to its toxicological and pharmaceutical functions [58,59]. For instance, a man-made mutant of BChE was recently modeled and examined by Masson and coworkers, using QM/MM and Markov state analysis, and was determined to be a template for future investigations into organophosphate hydrolase functions [60].
2.1.4. Oxyanion Hole
The oxyanion hole (OAH, orange in Figure 1), made of two glycines and one alanine [35,47], is generally a two-pronged site in many proteases and hydrolases; in the case of AChE, and subsequently BChE, the OAH is a three-pronged hole [61]. Early MD simulations of AChE phosphonylation suggested that the OAH exerts a pulling force on leaving groups during the alkylation step [45], and the OAH is known to lower the energy barrier for ACh hydrolysis in both cholinesterases [62]. QM/MM simulations exhibited consistent, tight hydrogen bonding between the OAH and the carbonyl carbon of the substrate, suggesting that the OAH facilitates stabilizing interactions in intermediate and transition states [61].
2.1.5. Ω-loop (Omega Loop)
The Ω-loop (OML, purple in Figure 1), consisting of a series of nearly 30 residues, is located along one side of the active site gorge wall. In the presence of a substrate or inhibitor, the Ω-loop plays an important role in modulating enzyme “breathing”, regulating the size of the gorge, and thereby enabling the passage of the ligand to the active site [21]. The OML undergoes conformational changes, such as gorge enlargement, through torsional motion and segmental fluctuations [63,64]. Unregulated motions and decreased electrostatic interactions, however, can significantly decrease the binding affinity of these enzymes, such as the case of the atypical mutation from Asp70 to Gly70 in the OML of BChE, which Masson et al. reported could increase the Km values 10- to 100-fold [65]. Moreover, the Ω-loop is speculated to facilitate an alternative entrance, the proposed side-door model noted above, to the active site of the enzyme; MD simulations performed by Wiesner et al. suggested that protonation of the AChE active site leads to conformational changes within the Ω-loop at Asn87 and Glu84 that result in the opening of this side door [66]. A similar observation was recently reported by the Rydzewski laboratory for TcAChE, where opening and closing of this side door due to the displacement of the OML favored alternative dissociate routes of the substrate and inhibitor [67]. To put these computational results into perspective, a number of experimental investigations into the omega loop and backdoor of AChE from various species have suggested that the back door opens in some cases [50,68,69], but also that this opening is likely not relevant functionally [70].
This entry is adapted from the peer-reviewed paper 10.3390/biom11040580
