Structure and therapeutic uses of butyrylcholinesterase: Application in detoxification, Alzheimer's disease, and fat metabolism

Structural information of butyrylcholinesterase (BChE) and its variants associated with several diseases are discussed here. Pure human BChE has been proved safe and effective in treating organophosphorus (OPs) poisoning and has completed Phase 1 and 2 pharmacokinetic (PK) and safety studies. The introduction of specific mutations into native BChE to endow it a self-reactivating property has gained much progress in producing effective OPs hydrolases. The hydrolysis ability of native BChE on cocaine has been confirmed but was blocked to clinical application due to poor PK properties. Several BChE mutants with elevated cocaine hydrolysis activity were published, some of which have shown safety and efficiency in treating cocaine addiction of human. The increased level of BChE in progressed Alzheimer's disease patients made it a promising target to elevate acetylcholine level and attenuate cognitive status. A variety of selective BChE inhibitors with high inhibitory activity published in recent years are reviewed here. BChE could influence the weight and insulin secretion and resistance of BChE knockout (KO) mice through hydrolyzing ghrelin. The BChE-ghrelin pathway could also regulate aggressive behaviors of BChE-KO mice.

Enantioseparation, in vitro testing, and structural characterization of triple-binding reactivators of organophosphate-inhibited cholinesterases

The enantiomers of racemic 2-hydroxyimino-N-(azidophenylpropyl)acetamide-derived triple-binding oxime reactivators were separated, and tested for inhibition and reactivation of acetylcholinesterase (AChE) and butyrylcholinesterase (BChE) inhibited with tabun (GA), cyclosarin (GF), sarin (GB), and VX. Both enzymes showed the greatest affinity toward the methylimidazole derivative (III) of 2-hydroxyimino-N-(azidophenylpropyl)acetamide (I). The crystal structure was determined for the complex of oxime III within human BChE, confirming that all three binding groups interacted with active site residues. In the case of BChE inhibited by GF, oximes I (kr = 207 M-1 min-1) and III (kr = 213 M-1 min-1) showed better reactivation efficiency than the reference oxime 2-PAM. Finally, the key mechanistic steps in the reactivation of GF-inhibited BChE with oxime III were modeled using the PM7R6 method, stressing the importance of proton transfer from Nε of His438 to Oγ of Ser203 for achieving successful reactivation.

Detection of Chemical Weapon Agents Using Spectroscopic Probes: A Computational Study

Organophosphorus compounds are organic compounds widely employed in agriculture as well as in chemical weapons. The use in agriculture is due to their insecticidal properties. However, in chemical warfare, the use of organophosphorus is associated with acetylcholinesterase inhibition, which promotes the cholinergic syndromes. In this line, the fast detection of this class of compound is crucial for the determination of environmental exposure. This improved detection will naturally allow for more prompt courses of treatment depending on the contaminant findings. In this perspective, the dipyrrinone oxime (1) was employed for the detection of organophosphorus compounds that are employed as nerve agents, such as cyclosarin, sarin, soman, diethyl chlorophosphate, diisopropylfluorophosphate, 2-(dimethylamino)ethyl N,N-dimethylphosphoramidofluoridate, O-ethyl-S-[2-(diethylamino)ethyl]methylphosphonothioate, O-ethyl-S-[2(diisopropylamino)ethyl] methylphosphonothioate, and O,O-diethyl-S-[2-(diethylamino)ethyl] phosphorothioate, through fluorescent emission. The thermodynamics and kinetic parameters as well as spectroscopic properties of the complexes formed for 1 and all organophosphorus compounds previously cited were investigated by means of theoretical calculations. From our findings, only the diethyl chlorophosphate, 2-(dimethylamino)ethyl N,N-dimethylphosphoramidofluoridate, and O,O-diethyl-S-[2-(diethylamino)ethyl] phosphorothioate emitted fluorescence in the hexane, toluene, chloroform, dichloromethane, methanol, acetonitrile, water, and dimethyl sulfoxide solvents. The study of the absorption wavelength with the most polar solvent showed higher values compared to apolar solvents. In the same solvent, for instance, soman in hexane showed the lowest absorption wavelength value, 324.5 nm, and DCP the highest value, 330.8 nm. This behavior was observed in other tested solvents. The thermodynamic parameters indicate negative Gibbs free energy (ΔG) values for the O-ethyl-S-[2(diisopropylamino)ethyl] methylphosphonothioate with 1 reaction. On the other hand, the sarin and cyclosarin revealed the lowest Gibbs free energy (ΔG‡) values, being kinetically favorable and presenting more reactivity.

Unexpected protonation state of Glu197 discovered from simulations of tacrine in butyrylcholinesterase

Butyrylcholinesterase (BChE) has been actively involved in drug discoveries from many fields for decades. In the crystal structure of the BChE-tacrine complex, there is an unanticipated formyl-proline molecule resolved very close to tacrine, raising an essential question on how reliable it is to apply the binding pose in a crystal structure to analyze related experimental observations, in which no formyl-proline is actually involved. In this study, by performing a series of 100 ns molecular dynamics simulations, we demonstrate that it is safe to employ the structural information from this crystal structure to analyze related experimental observations. Surprisingly, Glu197 needs to be protonated to have the structures simulated appropriately. It should be noted that Glu197 has been commonly considered as deprotonated in diverse analyses due to its low pKa in aqueous solution, for which some interpretations are inconsistent or unclear. Our further investigation shows that the protonated Glu197 plays a very important role in preserving His438 within the catalytic triad through stabilizing a highly conserved water molecule. Interestingly, the catalytic triad and Glu197 have been long recognized for possibly deviating largely from the crystal structure, which might be catalytically deficient and is generally considered to result from the difference between the crystal and aqueous environment. Herein, our results suggest that the large deviations of the catalytic triad and Glu197 from the crystal structure are caused by the inappropriate protonation state of Glu197. This finding shall provide an important clue that has been long missing for a better understanding of BChE-related puzzles or even reconsideration of some BChE-catalyzed reaction mechanisms.

QM/MM Description of Newly Selected Catalytic Bioscavengers Against Organophosphorus Compounds Revealed Reactivation Stimulus Mediated by Histidine Residue in the Acyl-Binding Loop

Butyrylcholinesterase (BChE) is considered as an efficient stoichiometric antidote against organophosphorus (OP) poisons. Recently we utilized combination of calculations and ultrahigh-throughput screening (uHTS) to select BChE variants capable of catalytic destruction of OP pesticide paraoxon. The purpose of this study was to elucidate the molecular mechanism underlying enzymatic hydrolysis of paraoxon by BChE variants using hybrid quantum mechanical/molecular mechanical (QM/MM) calculations. Detailed analysis of accomplished QM/MM runs revealed that histidine residues introduced into the acyl-binding loop are always located in close proximity with aspartate residue at position 70. Histidine residue acts as general base thus leading to attacking water molecule activation and subsequent SN2 inline hydrolysis resulting in BChE reactivation. This combination resembles canonical catalytic triad found in active centers of various proteases. Carboxyl group activates histidine residue by altering its pK_a, which in turn promotes the activation of water molecule in terms of its nucleophilicity. Observed re-protonation of catalytic serine residue at position 198 from histidine residue at position 438 recovers initial configuration of the enzyme’s active center, facilitating next catalytic cycle. We therefore suggest that utilization of uHTS platform in combination with deciphering of molecular mechanisms by QM/MM calculations may significantly improve our knowledge of enzyme function, propose new strategies for enzyme design and open new horizons in generation of catalytic bioscavengers against OP poisons.

Optimization of Cholinesterase-Based Catalytic Bioscavengers Against Organophosphorus Agents

Organophosphorus agents (OPs) are irreversible inhibitors of acetylcholinesterase (AChE). OP poisoning causes major cholinergic syndrome. Current medical counter-measures mitigate the acute effects but have limited action against OP-induced brain damage. Bioscavengers are appealing alternative therapeutic approach because they neutralize OPs in bloodstream before they reach physiological targets. First generation bioscavengers are stoichiometric bioscavengers. However, stoichiometric neutralization requires administration of huge doses of enzyme. Second generation bioscavengers are catalytic bioscavengers capable of detoxifying OPs with a turnover. High bimolecular rate constants (k_cat/K_m > 10⁶ M⁻¹min⁻¹) are required, so that low enzyme doses can be administered. Cholinesterases (ChE) are attractive candidates because OPs are hemi-substrates. Moderate OP hydrolase (OPase) activity has been observed for certain natural ChEs and for G117H-based human BChE mutants made by site-directed mutagenesis. However, before mutated ChEs can become operational catalytic bioscavengers their dephosphylation rate constant must be increased by several orders of magnitude. New strategies for converting ChEs into fast OPase are based either on combinational approaches or on computer redesign of enzyme. The keystone for rational conversion of ChEs into OPases is to understand the reaction mechanisms with OPs. In the present work we propose that efficient OP hydrolysis can be achieved by re-designing the configuration of enzyme active center residues and by creating specific routes for attack of water molecules and proton transfer. Four directions for nucleophilic attack of water on phosphorus atom were defined. Changes must lead to a novel enzyme, wherein OP hydrolysis wins over competing aging reactions. Kinetic, crystallographic, and computational data have been accumulated that describe mechanisms of reactions involving ChEs. From these studies, it appears that introducing new groups that create a stable H-bonded network susceptible to activate and orient water molecule, stabilize transition states (TS), and intermediates may determine whether dephosphylation is favored over aging. Mutations on key residues (L286, F329, F398) were considered. QM/MM calculations suggest that mutation L286H combined to other mutations favors water attack from apical position. However, the aging reaction is competing. Axial direction of water attack is not favorable to aging. QM/MM calculation shows that F329H+F398H-based multiple mutants display favorable energy barrier for fast reactivation without aging.

Understanding the non-catalytic behavior of human butyrylcholinesterase silent variants: Comparison of wild-type enzyme, catalytically active Ala328Cys mutant, and silent Ala328Asp variant

Conformational dynamics of wild-type human butyrylcholinesterase (BChE), two mutants of residue Ala328, the catalytically active Ala328Cys, and the catalytically inactive (silent) Ala328Asp, and their interactions with butyrylcholine were studied. The aim was to understand the molecular mechanisms by which point mutations may lead to silent BChE variant or alter catalytic activity. Importance of BChE natural variants is due to medical consequences, i.e. prolonged apnea, following administration of the myorelaxant esters, succinylcholine and mivacurium. Comparison of molecular dynamics (MD) simulations for the three model systems showed that: 1) the active mutant Ala328Cys mutant has some changes in configuration of catalytic residues, which do not prevent binding of butyrylcholine to the active site; 2) in the naturally-occurring silent variant Ala328Asp, the Asp328 carboxylate may either form a salt bridge with Lys339 or a H-bond with His438. In the first case, the Ω-loop swings off the gorge, disrupting the π-cation binding site and the catalytic triad. In the second case, binding of cationic substrates in the catalytic center is also impaired. MD simulations carried out in 0.15 M NaCl, close to physiological ionic strength conditions, favored the second situation. It was seen that Asp328 forms a H-bond with the catalytic triad His438, which in turn disrupts the catalytic machinery. Therefore, we concluded that the Ala328Asp variant is not catalytically active because of that dramatic event. Computational results, consistent with in vitro biochemical data and clinical observations, validate our MD approach.