Resurrection and Reactivation of Acetylcholinesterase and Butyrylcholinesterase.

Organophosphorus (OP) nerve agents and pesticides present significant threats to civilian and military populations. OP compounds include the nefarious G and V chemical nerve agents, but more commonly, civilians are exposed to less toxic OP pesticides, resulting in the same negative toxicological effects and thousands of deaths on an annual basis. After decades of research, no new therapeutics have been realized since the mid-1900s. Upon phosphylation of the catalytic serine residue, a process known as inhibition, there is an accumulation of acetylcholine (ACh) in the brain synapses and neuromuscular junctions, leading to a cholinergic crisis and eventually death. Oxime nucleophiles can reactivate select OP-inhibited acetylcholinesterase (AChE). Yet, the fields of reactivation of AChE and butyrylcholinesterase encounter additional challenges as broad-spectrum reactivation of either enzyme is difficult. Additional problems include the ability to cross the blood brain barrier (BBB) and to provide therapy in the central nervous system. Yet another complication arises in a competitive reaction, known as aging, whereby OP-inhibited AChE is converted to an inactive form, which until very recently, had been impossible to reverse to an active, functional form. Evaluations of uncharged oximes and other neutral nucleophiles have been made. Non-oxime reactivators, such as aromatic general bases and Mannich bases, have been developed. The issue of aging, which generates an anionic phosphylated serine residue, has been historically recalcitrant to recovery by any therapeutic approach-that is, until earlier this year. Mannich bases not only serve as reactivators of OP-inhibited AChE, but this class of compounds can also recover activity from the aged form of AChE, a process referred to as resurrection. This review covers the modern efforts to address all of these issues and notes the complexities of therapeutic development along these different lines of research.

Study of para-Quinone Methide Precursors toward the Realkylation of Aged Acetylcholinesterase.

Acetylcholinesterase (AChE) is an essential enzyme that can be targeted by organophosphorus (OP) compounds, including nerve agents. Following exposure to OPs, AChE becomes phosphylated (inhibited) and undergoes a subsequent aging process where the OP-AChE adduct is dealkylated. The aged AChE is unable to hydrolyze acetylcholine, resulting in accumulation of the neurotransmitter in the central nervous system (CNS) and elsewhere. Current therapeutics are only capable of reactivating inhibited AChE. There are no known therapeutic agents to reverse the aging process or treat aged AChE. Quinone methides (QMs) have been shown to alkylate phosphates under physiological conditions. In this study, a small library of novel quinone methide precursors (QMPs) has been synthesized and examined as potential alkylating agents against model nucleophiles, including a model phosphonate. Computational studies have been performed to evaluate the affinity of QMPs for the aged AChE active site, and preliminary testing with electric eel AChE has been performed.

Quantitative structure-activity relationship models of chemical transformations from matched pairs analyses.

The concepts of activity cliffs and matched molecular pairs (MMP) are recent paradigms for analysis of data sets to identify structural changes that may be used to modify the potency of lead molecules in drug discovery projects. Analysis of MMPs was recently demonstrated as a feasible technique for quantitative structure-activity relationship (QSAR) modeling of prospective compounds. Although within a small data set, the lack of matched pairs, and the lack of knowledge about specific chemical transformations limit prospective applications. Here we present an alternative technique that determines pairwise descriptors for each matched pair and then uses a QSAR model to estimate the activity change associated with a chemical transformation. The descriptors effectively group similar transformations and incorporate information about the transformation and its local environment. Use of a transformation QSAR model allows one to estimate the activity change for novel transformations and therefore returns predictions for a larger fraction of test set compounds. Application of the proposed methodology to four public data sets results in increased model performance over a benchmark random forest and direct application of chemical transformations using QSAR-by-matched molecular pairs analysis (QSAR-by-MMPA).

C2 hydroxyl group governs the difference in hydrolysis rates of methyl-α-D-glycero-D-guloseptanoside and methyl-ß-D-glycero-D-guloseptanoside.

A computational investigation into the hydrolysis of two methyl septanosides, methyl-α-D-glycero-D-guloseptanoside and methyl-ß-D-glycero-D-guloseptanoside was undertaken. These septanosides were chosen as model compounds for comparison to methyl pyranosides and allowed direct comparison of α versus ß hydrolysis rates for a specific septanoside isomer. Results suggest that hydrolysis takes place without proceeding through a transition state, an observation that was suggested in previous computational studies on exocyclic bond cleavage of carbohydrates. A conformational analysis of α- and ß-anomers 1 and 2 and their corresponding oxocarbenium 3, coupled with relaxed potential energy surface (PES) scans (M06-2X/6-311+G**, implicit methanol), indicated that hydrolysis of the α-anomer is favored by 1-2 kcal/mol over the ß-anomer, consistent with experiment. Model systems revealed that the lowest energy conformations of the septanoside ring system destabilize the ß-anomer by 2-3 kcal/mol relative to the α-anomer, and the addition of a single hydroxyl group at the C2-position on a minimal oxepane acetal can reproduce the PES for the septanoside 1. These results suggest that the C2 hydroxyl plays a unique role in the hydrolysis mechanism, destabilizing the septanoside via its proximity to the anomeric carbon and also through its interaction with the departing methanol from the α-anomer via hydrogen-bonding interactions.

Butyrylcholinesterase and G116H, G116S, G117H, G117N, E197Q and G117H/E197Q mutants: a molecular dynamics study.

Butyrylcholinesterase (BuChE) is a stoichiometric bioscavenger against organophosphorus (OP) nerve agent poisoning, and efforts to make BuChE variants that are catalytically active against a wide spectrum of nerve agents have been ongoing for the last decade. In order to understand the structural consequences for BuChE, we carried out extensive molecular dynamics (MD) simulations on wild-type BuChE (PDB ID: 1P0I) and several known and new variants of this enzyme, but without the presence of any ligand in the active site. The MD simulations on WT-BuChE identified two labile orientations for the catalytic serine, and also showed the likelihood of a backdoor. Upon changes at the G116 position, severe alterations around the active site region were identified. Simulations on both G117H and G117N variants showed the existence of a bound water molecule that is in close proximity to S198. Modeling of the E197Q mutant suggested that Q197 can be in two distinct orientations, one similar to the E202Q-AChE crystal structure and another in proximity to G439 and E441. The double mutant, G117H/E197Q, was found to have structural characteristics of both G117H and E197Q. In light of the computational results, previous experimental observations are discussed.

Reaction profiles of the interaction between sarin and acetylcholinesterase and the S203C mutant: model nucleophiles and QM/MM potential energy surfaces.

The phosphonylation mechanism of AChE and the S203C mutation by sarin (GB) is evaluated using two reaction schemes: a small model nucleophile (ethoxide, CH(3)CH(2)O(-)) and quantum mechanical/molecular mechanical (QM/MM) simulations. Calculations utilizing small model nucleophiles indicate that the reaction barrier for addition to GB is the rate-limiting step for both ethoxide and ethyl thiolate (CH(3)CH(2)S(-)); moreover, the activation barrier for addition to the phosphorus center of GB by ethyl thiolate is significantly larger (13.2 kcal/mol) than for ethoxide (8.3 kcal/mol). The decomposition transition state for both nucleophiles was determined to be approximately 1 kcal/mol. QM/MM simulations for AChE suggest a similar reaction mechanism for phosphonylation of the catalytic S203; however, the relative energetics are altered significantly compared to the isolated system. QM/MM results indicate that formation of the penta-coordinate intermediate is the rate-limiting step in the enzymatic system, with an activation barrier of 3.6 kcal/mol. Hydrogen-bonding interactions between the fluoride leaving group of GB with Y124 in AChE are observed throughout the reaction profile. The S203C mutation alters the relative energetics of the reaction, increasing the energy barrier for formation of the penta-coordinate intermediate to a value of 4.5 kcal/mol; moreover, the penta-coordinate intermediate (as product) is stabilized by an additional 6 kcal/mol when compared to wild-type AChE.

Computational Modeling of Human Paraoxonase 1: Preparation of Protein Models, Binding Studies, and Mechanistic Insights.

The enzyme human paraoxonase 1 (huPON1) has demonstrated significant potential for use as a bioscavenger for treatment of exposure to organophosphorus (OP) nerve agents. Herein we report the development of protein models for the human isoform derived from a crystal structure of a chimeric version of the protein (pdb ID: 1V04) and a homology model derived from the related enzyme diisopropylfluorophosphatase (pdb ID: 1XHR). From these structural models, binding modes for OP substrates are predicted, and these poses are found to orient substrates in proximity to residues known to modulate specificity of the enzyme. Predictions are made with regard to the role that residues play in altering substrate binding and turnover, in particular with regard to the stereoselectivity of the enzyme, and the known differences in activity related to a natural polymorphism in the enzyme. Potential mechanisms of action of the protein for catalytic hydrolysis of OP substrates are also evaluated in light of the proposed binding modes.

Chemical synthesis of two series of nerve agent model compounds and their stereoselective interaction with human acetylcholinesterase and human butyrylcholinesterase.

Both G and V type nerve agents possess a center of chirality about phosphorus. The S(p) enantiomers are generally more potent inhibitors than their R(p) counterparts toward acetylcholinesterase (AChE) and butyrylcholinesterase (BChE). To develop model compounds with defined centers of chirality that mimic the target nerve agent structures, we synthesized both the S(p) and the R(p) stereoisomers of two series of G type nerve agent model compounds in enantiomerically enriched form. The two series of model compounds contained identical substituents on the phosphorus as the G type agents, except that thiomethyl (CH(3)-S-) and thiocholine [(CH(3))(3)NCH(2)CH(2)-S-] groups were used to replace the traditional nerve agent leaving groups (i.e., fluoro for GB, GF, and GD and cyano for GA). Inhibition kinetic studies of the thiomethyl- and thiocholine-substituted series of nerve agent model compounds revealed that the S(p) enantiomers of both series of compounds showed greater inhibition potency toward AChE and BChE. The level of stereoselectivity, as indicated by the ratio of the bimolecular inhibition rate constants between S(p) and R(p) enantiomers, was greatest for the GF model compounds in both series. The thiocholine analogues were much more potent than the corresponding thiomethyl analogues. With the exception of the GA model compounds, both series showed greater potency against AChE than BChE. The stereoselectivity (i.e., S(p) > R(p)), enzyme selectivity, and dynamic range of inhibition potency contributed from these two series of compounds suggest that the combined application of these model compounds will provide useful research tools for understanding interactions of nerve agents with cholinesterase and other enzymes involved in nerve agent and organophosphate pharmacology. The potential of and limitations for using these model compounds in the development of biological therapeutics against nerve agent toxicity are also discussed.

Hydrolysis of nerve agents by model nucleophiles: a computational study.

Density functional theory calculations were employed to study the reaction of five nerve agents with model nucleophiles, including EtX(-) and EtXH (X=O, S, Se) for serine, cysteine and selenocysteine, respectively. Calculations at the B3LYP/6-311++G(2d,p) level of theory predict an exothermic reaction between ethoxide and all of the nerve agents studied. As compared to EtO(-) as a nucleophile, these reactions become approximately 30 kcal/mol more endothermic for EtS(-), and by approximately 40 kcal/mol for EtSe(-). The equivalent reactions with the neutral nucleophiles (EtXH) were more endothermic. The effect of solvation on the reaction thermochemistry was determined using a polarizable continuum model simulating the dielectric constant of chloroform. While there was a large exothermic shift for reactions involving charged nucleophiles with solvation modeling, the corresponding shift was minimal for the reaction with neutral nucleophiles.