PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of nihpaAbout Author manuscriptsSubmit a manuscriptNIH Public Access; Author Manuscript; Accepted for publication in peer reviewed journal;
 
Chem Biol Interact. Author manuscript; available in PMC Sep 25, 2009.
Published in final edited form as:
PMCID: PMC2572230
NIHMSID: NIHMS72523
Reactivation of model cholinesterases by oximes and intermediate phosphyloximes: A computational study
Shubham Vyas and Christopher M. Hadad*
Department of Chemistry, The Ohio State University, 100 W. 18th Ave. Columbus, OH 43210 USA
hadad.1/at/osu.edu
Phosphyloximes (POX) are generated upon the reactivation of organophosphorus (OP) inhibited cholinesterases (ChEs) by pyridinium oximes. These POXs are known to be potent inhibitors of the ChEs following reactivation. However, they can also decompose to give an OP derivative and a cyano derivative of the oxime when a base abstracts the benzylic proton. Using density functional theory, thermodynamic properties were calculated for the reactivation and decomposition pathways of three different oximes (2-PAM, 3-PAM and 4-PAM) with six different OPs (cyclosarin, paraoxon, sarin, tabun, VR and VX). For reactivation purposes, 2-PAM is predicted to be more efficient than 3- and 4-PAM. Based on atomic charges and relative energies, 2-POXs were found to be more inclined towards the decomposition process.
Keywords: Cholinesterase, phosphyloximes, inhibition, decomposition, density functional theory
Reactivation of organophosphorus (OP) inhibited cholinesterases (ChEs) can be accomplished by using pyridinium based oximes.1 Many oximes have been designed and tested against a variety of OPs in the past, but a general and efficacious remedy against all of the OPs has yet to be discovered. Phosphyloximes (POXs) are formed during the reactivation process by the reaction of the oxime with the OP, and these products are potent inhibitors of ChEs and can remain in the active site. In addition to inhibiting ChEs, POXs can also decompose to generate a cyano derivative as well as a phosphate analogue of the OP – a process initiated by base-induced abstraction of the benzylic hydrogen depicted in Figure 1.
Figure 1
Figure 1
Reactivation of the OP-inhibited ChEs by oxime and POX formation (top); POX decomposition pathway leading to cyano compound and OP analogue (bottom). (2-PAM and 2-POX are shown as the oxime in the top figure, 3-PAM and 3-POX, and 4-PAM and 4-POX are meta- (more ...)
To minimize the inhibition pathway for POXs, a thorough understanding of the thermodynamic and kinetic parameters governing the two competing processes is required. Doctor et al.2 have reported extensive data on the inhibition of ChEs by POXs formed upon reactivation of various OP-inhibited enzyme using 2-, 3- and 4- [hydroxyimino methyl]-1-methylpyridinium chloride (PAM) as shown in Figure 1. They concluded that POXs generated from 2-PAM (2-POXs) have a significantly higher rate of decomposition than the corresponding compounds formed from 3- and 4-PAM (3- and 4-POXs, respectively). However, the inhibition pathway was also highly facile for 2-POXs compared to the others. In addition, they performed some semi-empirical calculations, the results of which suggested that acidity of the benzylic hydrogen might play an important role in POX reactivity. However, at present there are no ab initio studies examining the relative energetics of the two pathways of POX activity which could be used to modulate the reactivation and inhibitory activity of novel oximes. In this study, we present a systematic computational investigation of the thermodynamic parameters for both pathways with six OPs (cyclosarin, paraoxon, sarin, tabun, VR and VX; Figure 2) and three oximes (2-, 3- and 4- PAM).
Figure 2
Figure 2
Organophosphorus (OP) compounds under study
All calculations were performed using the Gaussian 03 program suite3 at the Ohio Supercomputer Center. The geometries were optimized with density functional theory using Becke’s three parameter exchange functional and the correlation functional of Lee, Yang, and Parr (B3-LYP),4 with a 6-31+G(d) basis set.5 For computing the thermochemistries, the ChE was represented by an ethoxide group, simulating the catalytic serine. All structures were confirmed to be minima via vibrational frequency analyses; zero-point vibrational energy corrections were not scaled. Thermodynamic parameters discussed herein are ΔH values and were calculated at the B3LYP/6-311+G(d,p)//B3LYP/6-31+G(d) level of theory at 298K, unless noted otherwise. Electrostatic potential charges were calculated using the CHELPG6 method, and the inclusion of implicit solvation effects was performed using the polarizable continuum model (PCM) model for water.7
2-PAM is a widely used therapeutic for the reactivation of OP-inhibited ChEs. However, there are no comparative studies of the three PAM derivatives, which attempt to correlate computed properties with their reactivation efficacies. To this end, we have optimized the geometries of six energetically favorable conformers for 2-PAM and 3-PAM and three conformers of 4-PAM (Figure 3). It has been proposed that the ESA conformer of 2-PAM is the biologically active form.8 Recently, Harel et al.9 solved an x-ray crystal structure (2.71 Å resolution) with 2-PAM bound to AChE, and this study also suggested that 2-PAM is in the ESA or ESS conformation. However, the geometry of the oxime group about the exocyclic carbon atom is quite unusual with a C–C=N angle of ~172°. According to our calculations, the EAA and ESA conformers for 2- and 3-PAM were found to differ by only 0.5 kcal/mol (Table 1). The EAA (EA for 4-PAM) conformers were determined to be the most stable structures for all three PAMs. These results are consistent with the available X-ray structure10 and previous calculations performed for 2-PAM.11 Hence, the EAA conformers were used in calculating the overall reaction thermodynamics. All of the POXs prefer syn stereochemistry at the O–N=C(H)–C(Pyridine) bond, consistent with the conclusions from NMR experiments.12
Figure 3
Figure 3
Different conformers of 2-PAM considered for this study. (EAS conformer shows the naming convention of the conformers; EAA conformer shows the atomic numbering).
Table 1
Table 1
Relative energies (ΔE in kcal/mol) of different conformers of the PAMs at the B3LYP/6-311+G**//B3LYP/6-31+G* level. (Values in parentheses are ΔE including aqueous solvation at the PCM level.)
In the gas phase, all of the reactivation processes were found to be endothermic except for tabun, which is just slightly exothermic, as shown in Table 2 (the reverse process, inhibition, will thus be exothermic for all OPs except tabun). Although the enthalpy of reactivation for tabun is close to being thermoneutral, it is interesting to note that based on in vitro experiments performed by Szinicz et al.,13 these authors reported that reactivation was more favored for tabun relative to sarin. However, we should note that the present calculations are only thermodynamic predictions, and kinetic information will be needed to provide more insight to differentiate between this set of OPs. (Such studies are currently under investigation.) Further, these reactions became more endothermic upon solvation, due to an increase in charge delocalization for the POX compared to the parent oxime, which results in differential stabilization of the parent oximes over the POXs. 2-PAM was the best oxime for the reactivation process as suggested by the gas phase enthalpies. In contrast, enthalpic data for the reactivation in the aqueous phase suggested that 2-PAM was the least suitable oxime. This is likely the result of reduced solvent stabilization for 2-PAM and corresponding POXs, as they are less polar than their corresponding 3- and 4-substituted structural isomers. Magnitudes of the dipole moments of these species were found in the order 2-PAM < 3-PAM ≈ 4-PAM.
Table 2
Table 2
Calculated ΔH (298K) for the reactivation and decomposition processes in kcal/mol at the B3LYP/6-311+G**//B3LYP/6-31+G* level. (Values in parentheses are ΔH at 298K including aqueous solvation at the PCM level.)
Charge analyses are summarized in Table 3 for the oximes studied. Because the trend in the charge analysis is similar for both the gas phase and aqueous solution, we have shown only the aqueous phase results. A correlation was observed between increased charge separation at the oxime O–H bond and reaction exothermicity for the reactivation pathway. 2-PAM was found to have the highest charge separation at the oxime O–H bond, indicating that it is the most acidic among the three oximes studied. These results suggest that 2-PAM is also the most nucleophilic oxime in the set and therefore is a better candidate for increased reactivation than the others.
Table 3
Table 3
CHELPG charges (in units of electrons) at the B3LYP/6-311+G**//B3LYP/6-31+G* level in the aqueous phase for the oximes. (Atom labels are shown in Figure 1.)
The POXs can decompose in the presence of a base as shown in Figure 1. The reaction generates a negatively charged OP analogue and the conjugate acid of the base used, which is highly unfavorable from an energetic standpoint. Thermodynamic calculation were carried out using either water or ethoxide ion as a base, which yielded the reaction enthalpies for the decomposition pathway of about 160 kcal/mol and −55 kcal/mol respectively. These seemingly contradictory energetics for the decomposition process are consistent with the proton affinities of each nucleophile, water and ethoxide, respectively. Thus proton transfer to the OP analogue will take place simultaneously to the generation of the conjugate acid. Therefore, we hypothesize that the decomposition process occurs via a proton-shuttle mechanism as shown in Figure 4 and the corresponding reaction enthalpies are summarized in Table 2.
Figure 4
Figure 4
Hypothesized decomposition pathway for POXs via proton-shuttle mechanism.
All of the decomposition processes were found to be exothermic overall, confirming that decomposition is a thermodynamically viable pathway for POX species. Based on the decomposition enthalpies, 3-POXs were found to be more favorable toward the decomposition pathway. This result is supported by the fact that 2- and 4-POXs are resonance stabilized but 3-POXs are not. The exothermicities of the decomposition process were nearly doubled in the aqueous phase, because two molecules were solvated on the product side which yielded higher solvation energy compared to the parent POX at the reactant side. This was confirmed by comparing the gas phase and the solution phase energies.
By comparing the energetics among the sets of POXs, we found that 2-POXs were 3–5 kcal/mol less stable than 3- and 4-POXs. This could be explained by the steric bulk provided by the pyridinium methyl group around the oxime-OP linkage, which is not present in the case of 3-POXs and 4-POXs. This observation was also supported by previous work2 which showed that 3- and 4-POXs were 100–500 times more stable than the 2-POXs. Additionally, we performed CHELPG charge analysis for all of the POX species to gain insight into the parameters governing the decomposition pathway of the POXs (Table 4.) The trends in the charges were similar to the thermochemical results: charge separation at the O–N bond was about 3–5 fold higher in the 2-POXs than the others, indicating a higher propensity for decomposition, and thus lower stability. In addition, the benzylic proton had the highest positive charge concentration in the case of 2-PAM, which correlates with the increased propensity for decomposition in 2-POXs.
Table 4
Table 4
CHELPG charges (in units of electrons) at the B3LYP/6-311+G**//B3LYP/6-31+G* level in the aqueous phase for the POXs. (Atom labels are shown in Figure 1.)
In summary, the reactivation process was found to be endothermic for all systems and 2-PAM was concluded to be a more efficient reactivator than 3- and 4-PAM for the reactivation process, due to higher charge separation at the oxime bond. From an enthaplic point of view, the protein’s active site environment will provide an electrostatic environment which will be somewhat in between the gas phase and aqueous phase; quantum mechanics/molecular mechanics calculations are in process to obtain more refined results as to the electrostatic effect of the active site. The 2-POXs were found to be less stable than the others, which we attributed to the steric effects of the pyridinium methyl group and also reflected by the higher acidity of its benzylic proton which encouraged the decomposition. The calculations performed herein indicate that 2-substitution improves oxime performance with regard to reactivation and decomposition. Further calculations to characterize the transition states for the processes to complete a kinetic picture are currently in progress.
Acknowledgement
Support of this work by a generous grant from US Army Medical Research Institute of Chemical Defense (MRICD) and computational facilities from the Ohio Supercomputer Center is gratefully acknowledged.
Footnotes
This is a PDF file of an unedited manuscript that has been accepted for publication. As a service to our customers we are providing this early version of the manuscript. The manuscript will undergo copyediting, typesetting, and review of the resulting proof before it is published in its final citable form. Please note that during the production process errors may be discovered which could affect the content, and all legal disclaimers that apply to the journal pertain.
1. Worek F, Szinicz L, Eyer P, Thiermann H. Evaluation of oxime efficacy in nerve agent poisoning:Development of a kinetic-based dynamic model. Toxicol Appl Pharmacol. 2005;209:193–202. [PubMed]
2. Ashani Y, Bhattacharjee AK, Leader H, Saxena A, Doctor BP. Inhibition of cholinesterases with cationic phosphonyl oximes highlights distinctive properties of the charged pyridine groups of quaternary oxime reactivators. Biochem Pharmacol. 2003;66:191–202. [PubMed]
3. Gaussian 03, Revision C.02. Pittsburgh PA: Gaussian, Inc.; 2003. For full citation please refer to http://www.gaussian.com/citation_g03.htm.
4. (a) Becke AD. Density-functional exchange-energy approximation with correct asymptotic behavior. Phys Rev A. 1988;38:3098–3100. [PubMed] (b) Lee C, Yang W, Parr RG. Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density. Phys Rev B. 1988;37:785–789. [PubMed]
5. Hariharan PC, Pople JA. The effect of d-functions on molecular orbital energies for hydrocarbons. Chem Phys Lett. 1972;16:217–219.
6. Breneman CM, Wiberg KB. Determining atom-centered monopoles from molecular electrostatic potentials. The need for high sampling density in formamide conformational analysis. J Comp Chem. 1990;11:361–373.
7. Tomasi J, Mennucci B, Cammi R. Quantum mechanical continuum solvation models. Chem Rev. 2005;105:2999–3094. [PubMed]
8. Hagedorn I, Stark I, Lorenz HP. Reactivation of phosphorylated acetyl-cholinesterase. Dependence upon activator acidity. Angew Chem Int Ed. 1972;11:307–309. [PubMed]
9. Harel M, Silman I, Sussman JL. The structure of torpedo californica acetylcholinesterase complexed with 2-PAM (PDB ID: 2VB4) to be published.
10. Carlstrom D. Crystallographic study of N-methyl-pyridine-2-aldoxime (2-PAM) halides. Acta Chem Scand. 1966;20:1240–1246. [PubMed]
11. Castro AT, Figueroa-Villar JD. Molecular structure, conformational analysis and charge distribution of Pralidoxime: ab Initio and DFT studies. Int J Quantum Chem. 2002;89:135–146.
12. Crepaux D, Lehn JM. Nuclear spin-spin interactions. XII. Nitrogen-15-proton coupling constants in the determination of oxime configurations. Org Magn Reson. 1975;7:524–526.
13. Becker G, Kawan A, Szinicz L. Direct reaction of oximes with sarin, soman, or tabun in vitro. Arch Toxicol. 1997;71:714–718. [PubMed]