Understand the interaction modes and reactivity of trimedoxime toward MmAChE inhibited by nerve agents: theoretical and experimental Aspects

doi:10.21203/rs.3.rs-71128/v1

The organophosphorus compounds (OP) are used as both chemical weapons and pesticides. However, these agents are very dangerous and toxic to humans, animals, and the environment. Thus, investigations with reactivators have been deeply developed in order to design new antidotes with better efficiency, as well as greater spectrum of action in the AChE reactivation process. With that in mind, in this work, we investigated the behavior of trimedoxime toward the Mus musculus Acetylcholinesterase (MmAChE) inhibited by a range of nerve agents, such as chemical weapons. From experimental essays, reactivation percentages were obtained for the reactivation of different AChE-OP complexes. On the other hand, theoretical calculations were performed to assess the differences of interaction modes and reactivity of trimedoxime within AChE active site. Comparing theoretical and experimental data, it is possible to notice that the oxime, in most case, showed better reactivation percentages at higher concentrations, with the best result for the reactivation of the AChE-VX adduct. From this work, it was revealed that the mechanistic process contributes most to the oxime efficiency than the interaction in the site. In this way, this study was important to better understand the reactivation process through trimedoxime, contributing to the proposal of novel antidotes.

Chemical Biology

Computational Biology

Biophysics

Medicinal Chemistry

Nerve agents

acetylcholinesterase

trimedoxime

reactivation

mechanistic studies

computational methods

Chemical weapons are defined as any chemical substance whose toxic properties are used for the purpose of killing, injuring or incapacitating an enemy in war or associated with military operations [1–3]. Even with the efforts from world entities to ban the use of chemical weapons, under the Chemical Weapons Convention, diverse countries still have an arsenal of these chemical substances [4–6]. Among these chemical weapons, the most toxic to humans are the well-known nerve agents, whose base structure consists of an organophosphorus compound (OP). In addition, a range of these toxic substances show potential for application in agricultural and industrial sectors [7].

From this class of OP substances, the pesticides are fundamental in agroindustrial applications [8]. Although the OP compounds are widely used to pest control, they are very dangerous and toxic to humans, animals, and the environment. These sorts of compounds act by inhibiting the acetylcholinesterase (AChE) enzyme causing a cholinergic neurotoxic effect. In the presence of OP, the Ser203 residue from the AChE active site covalently binds to the phosphorus atom forming a phosphorylated complex [2,9]. The AChE inhibition causes an accumulation of acetylcholine (ACh), once this enzyme is responsible for the hydrolysis of this neurotransmitter. This toxic framework results in the ACh accumulation, giving rise to the cholinergic syndrome, which is a set of symptoms associated with poisoning from certain toxic substances, such as OP nerve agents, caused by the overstimulation of muscarinic and nicotinic receptors [10]. Among the major symptoms of the intoxication, we can cite excessive salivation, lacrimation, urination, sweating, broncho-constriction and neuromuscular block, leading to death in severe cases of poisoning [10–12].

The current treatment protocol for OP poisoning consists mainly of the employment of a reactivating agent, commonly an oxime compound [13–16], which is capable of restoring the AChE catalytic activity through a nucleophilic attack, thus remediating the intoxication effects and reestablishing the ACh levels [17–19]. The general reactivation mechanism through oximes is represented in figure 1.

In view of the exposed so far, it is important to notice that there is no universal antidote to date, that is, a broad-spectrum oxime capable of reactivating all types of OP-inhibited AChE. In recent years, efforts have focused on the screening and identification of potent oximes, with sufficient permeability through the blood-brain barrier (BBB), maintaining a high reactivation rate [19–21]. In this study, we present a theoretical and experimental investigation to better understand the reactivation mechanism of the AChE inhibited by several kinds of OP agents. Based on interaction and mechanistic studies, we seek to explain the experimental data through molecular modeling in order to comprehend the reactivation process, by employing trimedoxime as the reactivating species (figure 2). We expect to understand the interaction modes and reactivity of trimedoxime in the reactivation process of the AChE-OP adduct.

2.1 In vitro test: experimental results

The results obtained through the experimental part of this work are summarized in

Table 1.

Table 1. Reactivation activity of trimedoxime (data obtained in triplicate experimental essays).

Trimedoxime
System	React. (%) Conc. 10-5 M	React. (%) Conc. 10-3 M
AChE-GA	0	30
AChE-GB	7	54
AChE-GF	0	0
AChE-GD	0	0
AChE-VX	9.8	85.3
AChE-POX	50	46
AChE-DDVP	17.3	31.5

According to the literature, the efficiency (reactivation percentage) of an oxime should be at least 10% to provide an appropriate remediation for the intoxicated patient [22,23]. In this regard, we observe different reactivation percentages through trimedoxime, taking into account its concentration as well as the type of OP-AChE complex. According to our experimental findings, note that trimedoxime demonstrated the best results at higher concentrations (10-3 M). At this concentration, the oxime showed a remarkable reactivation percentage of 85.3% for AChE-VX reactivation. At a concentration of 10-3 M, trimedoxime also provides a good performance in the reactivation of the AChE-GB (54%) and AChE-POX (46%) adducts. This result was more modest for AChE-GA (30%). An interesting outcome from this experimental investigation is the fact of trimedoxime does not reactivate the AChE-GF and AChE-GD adducts. These trends are more deeply approached in the next sections. From the experimental essays with trimedoxime at lower concentrations, we can observe that the experimental values indicate a significant reactivation percentage for the AChE-POX (50%), as well as a sufficient reactivation rate for AChE-DDVP (17.3%). Indeed, this oxime showed insufficient reactivating power for the AChE inhibited by the other OP agents investigated, such as GA, GF and GD, considering a concentration of 10-5 M.

2.2 Affinity and thermodynamics: docking results

According to the docking protocol, the calculations were performed in order to investigate the affinity between trimedoxime and inhibited AChE. For this, a cavity prediction algorithm based on a 3D box was used to find the binding sites in the inhibited enzyme active site. The active cavity presented a volume of 113.66 Å3, being appropriate to support the reactivator.

The molecular mechanics-based calculations generated diverse poses of trimedoxime within the cavity of the inhibited complexes, and the respective intermolecular interaction energy was computed to each system. As usual in these computations, the best oxime conformation was chosen for subsequent QM calculations, based on the lowest interaction energies as well as the most reactive conformations. Table 2 shows the values obtained from the docking calculations for the most appropriated poses of trimedoxime with different inhibited complexes.

Table 2. Docking results for trimedoxime inside different AChE-OP adducts.

Trimedoxime
System	*∆E (kcal mol**-1)	Residues
AChE-GA	-140.9	Ser298
AChE-GB	-154.7	Tyr124, Ser298, Arg296

AChE-GF	-161.3	Tyr124, Glu285
AChE-GD	-157.7	Tyr124, Glu285
AChE-VX	-115.0	Tyr124, Phe295, Arg296
AChE-POX	-144.1	Tyr124, Glu285
AChE-DDVP	-164.8	Arg296, Ser298, Trp286

*∆E = Intermolecular interaction energy.

According to the data reported in Table 2, note that trimedoxime showed stabilizing interactions within the inhibited enzyme complex site, for all OP agents investigated. From these results, the oxime demonstrated the lowest interaction energy in the AChE-DDVP (-164.8 kcal mol-1) adduct, followed by AChE-GF (-161.3 kcal mol-1) and AChE-GD (-157.7 kcal mol-1). In turn, the oxime showed a less stabilizing interaction energy within the AChE-VX cavity. As shown in the experimental section, at higher concentrations, the trimedoxime demonstrated to be more efficient in the reactivation of the AChE-VX adduct. This trend leads us to believe that the interaction energy is not the only factor responsible for the performance of this antidote in the reactivation, but others factors should be involved. In this regard, the results from the mechanistic study are presented in the next section.

From Table 2, the trimedoxime was stably docked in the inhibited AChE, with intermolecular interaction energy values in the range of -115.0 kcal mol-1 to -164.8 kcal mol-1. Diverse kinds of intermolecular interactions contribute to the stabilizing interaction in the site, such as hydrophobic interactions, electrostatic interactions and hydrogen bonds. It is important to mention that the AChE active site adopts distinct conformations according to the sort of OP agent. Thus, it is expected that trimedoxime performs different interactions with residues from the active site. These hydrogen bond-type interactions are generally the most important in studies of biological systems.

In most of the systems investigated, trimedoxime performed interaction with the Tyr124 amino acid residue, and according to the literature, this interaction is described as a possible π- π stacking, which takes place between Tyr124 residue and the pyridine ring of the oxime. This interaction is indicated to have an important role for helping in the transition state stabilization [2,24,25]. The hydrogen bonds revealed by the interaction of trimedoxime in each inhibited system are shown in figure 3.

From the discussed so far, it is important to notice that, together with the reactivation percentage, the interaction energy data do not explain thoroughly the experimental trends. The discussion of the mechanistic studies in the next section will rise new insights about the behavior of trimedoxime toward different AChE-OP systems.

2.3 Investigating kinetic parameters for biological activity: mechanistic studies

In the last part of this investigation, theoretical calculations were carried out to determine the relative activation energy (∆∆E#) through the hybrid QM/MM for the reactivation of each inhibited AChE system. The ∆E# values were computed based on the energy difference between the transition states and the initial system configurations from the reactants.

For the reaction mechanism simulation, steric and electronic effects of the chemical reactions are important aspects over the reaction pathway. In addition, the strain and interaction energies are significant contributing factors that dictate the reaction course. The interaction energy is responsible for stabilizing the reaction. On the other hand, the strain energy is responsible for distorting the reactants to adopt a pentacoordinate transition state. The relation between interaction and strain energies results in the height of the reaction barrier (∆E#), the so-called activation energy. This parameter was elucidated for some of these reactions in order to better comprehend the trimedoxime’s behavior in the reactivation process. For this, a combined procedure of docking and DFT calculations at the QM/MM interface for the mechanism was carried out. The transition states were characterized through potential energy curves. Table 3 shows the kinetic parameters ∆∆E#, as well as the experimental values of reactivation at the concentration of 10-3 M.

Table 3. Experimental reactivation percentage and relative activation energy for trimedoxime in the reactivation process.

Trimedoxime
System	∆∆E# * (kcal mol-1)	React. (%) Conc. 10-3 M
AChE-GA	46.83	30
AChE-GB	33.43	54
AChE-GF	-	0
AChE-GD	-	0
AChE-VX	0	85.3
AChE-POX	41.59	46
AChE-DDVP	47.75	31.5

*∆∆E# = Relative activation energy.

According to Table 3, these quantum theoretical results corroborate our experimental findings. Trimedoxime has shown itself to be very efficient in reactivating the inhibited AChE- VX at a concentration of 10^-3 M, which is according to the reactional barrier observed in the reactivation of this inhibited complex. From Table 3, the reactivation of the AChE-VX adduct revealed the lowest barrier. This fact helps explain the higher experimental reactivation percentage of the AChE inhibited by VX, which was 85.3%. This fact suggests that its transition state is better stabilized, allowing for the oxime to interact stronger with the nerve agent.

As we can see from Table 3, the reactivation of the AChE-GB complex showed the second more stabilizing barrier (33.43 kcal mol-1), which corroborates the second best reactivation percentage found in our experimental essays (54%). In addition, it was computed for the AChE- POX a barrier of 41.59 kcal mol-1, and for the AChE-GA and AChE-DDVP, our computations indicate very close barriers, 46.83 kcal mol^-1 and 47.75 kcal mol^-1 respectively, which correlate very well with the close experimental reactivation percentages found for these respective systems. It is worth mentioning that our simulations for the reactivation of the AChE-GF and AChE-GD were not succeeded in our study, that is, the AChE inhibited by these OP agents did not provide a feasible conformation for the nucleophilic attack by trimedoxime in the active site.

This fact could be explained from the interaction modes of these toxic agents in the site, mostly due to steric hindrance effects, as well as intermolecular interactions.

From the ΔΔE^# values in Table 3, we performed a multiple linear regression (MLR) between this parameter and the reactivation percentage, as well as the interaction energy. Our results revealed that the combination of interaction energy (ΔE) and activation energy (ΔΔE#) is able to efficiently explain the experimental outcomes. By increasing the number of system descriptors, a better correlation between theory and experiment is expected. Based on this, the MLR between the experimental and theoretical parameters resulted in the equation below. The regression was obtained with an excellent correlation value of 0.97.

By analyzing equation 1, we can observe some important trends about the studied systems. Starting with the correlation value from the MLR, it shows that the docking conjugated to QM/MM calculations result in a better representation of the systems investigated. According to the coefficients of the equation, the importance of each stage for the AChE reactivation process by trimedoxime can be evaluated. Note that the highest modulus of the coefficient of the term ΔΔE^# (relative activation energy) indicates that the reaction step presents a greater contribution for the AChE reactivation than the interaction energy [23,25]. This means that trimedoxime can more easily fit to the transition state structure in the reactivation process. In addition, the binding mode of trimedoxime in the site is not a critical step for activity. With the exposed in this investigation, we observe that trimedoxime stands for a significant advance in the development of more efficient reactivators for the remediation of the intoxication caused by neurotoxic nerve agents.

Previous studies have shown that there is not a direct correlation in oxime-mediated reactivation between species, and comparative studies in one species may not truly reflect the reactivation effects in humans. Due to structural differences, the active site of both enzymes from rat and human may adopt distinct conformations in the presence of the neurotoxic agent, and the antidote might be led to specific reactional behaviors. In this context, in silico and in vitro investigations with the human AChE are equally important. These aspects will be considered in future investigations [26].

3.1 Experimental details

Trimedoxime was prepared at the department of Toxicology in School of Military Health Sciences (Czech Republic), according to the synthesis route described earlier [27]. The purity of the reactivator was detected through the TLC and HPLC technique and NMR [28]. All compounds were obtained from the Brno Military Facility (95% purity and higher).

The animals employed in this experiment were handled under the supervision of the Ethics Committee of the School of Military Health Sciences in Hradec Kralove, Czech Republic. As a source of cholinesterases, a 10% rat brain homogenate (w/v) was used. The homogenate was prepared as described: Ether-narcotized rats (n=6) were killed by bleeding from a carotid artery. The brain was removed, washed with saline and homogenized using an Ultra-Turrax homogenizer, in distilled water.

For the in vitro test, 0.5 mL of brain homogenate was mixed with 20 μL of isopropanol solution of the selected nerve agent and distilled water (0.5 mL). The mixture was incubated for 30 minutes at 25oC to achieve 95% inhibition of AChE. 2.5 mL of sodium chloride (3 M) and distilled water were added to a volume of 23 mL. Finally, 2 mL of the substrate – ACh iodide (0.02 M) were added. The enzyme activity (analyzed by potentiometric titration of decomposed ACh iodide) was measured at pH 7.6 and 25°C on an autotitrator RTS 822 (Radiometer, Denmark). The same procedure was undertaken with the inhibited enzyme and further treatment with 10 min incubation with an aqueous solution of the reactivator (0.2 mL of 10-3 M), which replaced 0.2 mL of water. Activities of intact AChE (a₀), inhibited AChE (a_i) and reactivated AChE (a_r) were deduced from the consumption of NaOH solution (0.01 M) over time; NaOH reacted with acetate released from the decomposed ACh iodide. The reactivation percentage (%) was calculated from the measured data according to the formula (Equation 2):

The entire method is described in details in the work from Kuca and Cabal [29]. This same methodology was successfully employed in the work from Polisel et al (2019) [23].

3.2 Docking procedure

In the docking studies, the affinity of trimedoxime with the AChE inhibited by diverse OP agents was investigated. The oxime chemical structure was constructed and optimized at the DFT level, with B3LYP density functional and 6-31g(d,p) basis set, as implemented in the Gaussian 09 package [30]. The oxime was then docked inside the crystallographic structure of Mus musculus AChE (PDB code 3ZLU; resolution = 2.60 Å) [31] inhibited by GA (Tabun), GB (Sarin), GF (Cyclosarin), GD (Soman), VX, POX (Paraoxon) and DDVP (Dichlorvos), using the Molegro Virtual Docker program (MVD®) [32], according to similar procedures employed previously [33–35]. From our calculation protocol, it was considered a radius of about 20 Å, where the residues of the catalytic triad were kept as flexible. Due to the nature of the docking methods, the calculations were carried out, generating approximately 50 poses (hence such as conformation and orientation) for each ligand studied.

In the MVD program, the MolDock score algorithm method used as a scoring function is based on the piecewise linear potential, which fundamentally is a simplified potential whose parameters are in turn fitted to protein-ligand structures, binding data scoring functions and further extended in Generic Evolutionary Method for molecular docking, including a new hydrogen bonding term as well as new charge schemes [32]. Along this line, the docking scoring function values, E_score, are usually defined by Eq. 3:

Note that the E_PLP stands for ‘‘piecewise linear potential’’, which consists of the use of two different parameter sets, as described forward: one for the approximation of the steric term (i.e., Van der Waals) among atoms, as well as the other potential for the hydrogen bonding. As can be seen, the second term is, of course, related to the electrostatic interactions among overloaded atoms. Typically, it is a Coulomb potential with a dielectric constant dependent on the distance (which can be approximately described as D(r) = 4r). Hence, for this, the numerical value of 332.0 is responsible for the electrostatic energy unit to be given in kilocalories per molecule, as well [32].

E_intra is defined as the internal energy of each ligand. That is:

Note that the first part of the equation (double summation) is among all pairs of atoms in the ligand, taking off those connected by two bonds. Thus, in this equation, the second term denotes the torsional energy, where θ is the torsional angle of the bond. Hence, if several torsions could be determined, then, each torsional energy is considered as an average among them. Being that the last term, E_clash, assigns a penalty of about 1.000 if the distance between two heavy atoms (e.g., more than two bonds apart) is smaller than 2.0 Å, but not taking into account infeasible ligand conformations [32]. Thus, the docking search algorithm that is applied in the MVD program considers an evolutionary algorithm, that is, based on the interactive optimization techniques (inspired by Darwinian evolution theory), which implies a new hybrid search algorithm conveniently called guided differential evolution. As such, this hybrid combines the differential evolution optimization technique with a cavity prediction algorithm during the search process, allowing that way a fast and accurate identification of potential binding modes (poses) [32,36,37].

3.3 QM/MM procedure

In line with the large number of atoms present in the investigated systems, a quantum mechanics (QM)-based treatment becomes infeasible due to the high computational demand. However, the covalent bond re-arrangements in the reactional process cannot be ignored and treated exclusively through molecular mechanics (MM). In this context, the hybrid quantum mechanics-molecular mechanics (QM/MM) was employed in this investigation in order to study the reaction pathway involved in the reactivation process [38]. From this protocol, the AChE active site was treated through QM methods, DFT in this case, and the rest of the system was treated with MM-based methods [39]. From these calculations, the energetic barrier of the reactivation process of each enzyme-OP complex with trimedoxime was determined. This theoretical strategy has been previously employed in other works [26,39–44]. The QM part of the calculations was performed through the Gaussian 09 package, at the DFT level and 6- 31g(d,p) basis set [45,46]. The delimited QM region includes: Ser203 residue bound to the respective OP, the residues Tyr124, Phe295, Arg296, Glu285, Ser298 and Trp286, in addition to trimedoxime. In this simulation, all precursors, transition states and intermediates were calculated and characterized identifying imaginary frequencies [25,47,48]. Each system was fully optimized at the DFT level with conjugate gradient and quasi-Newton-Raphson algorithms. The final geometries were obtained with the density functional Becke’s three- parameter exchange functional and the gradient-corrected functional of Lee, Yang and Paar (B3LYP) [36,49], by using 6-31g(d, p) basis set.

In this work, we tested the in vitro efficiency of trimedoxime and applied computational techniques to evaluate the interaction modes and reactivity of this antidote, in the reactivation process of the AChE inhibited by a range of OP nerve agents. Thus, kinetic factors and interactions that govern the AChE enzyme reactivation process were investigated. With that in mind, our theoretical outcomes show that the active site of the inhibited AChE adopts different conformations according to the kind of neurotoxic agent. Therefore, these conformational changes in the site result in the different interactions and reactivity of trimedoxime in the active cavity.

Our findings indicate that the performance of trimedoxime enhances by increasing its concentration, being the best result found for the reactivation of the AChE-VX adduct. On the other hand, our experimental results show that trimedoxime was inefficient in the reactivation of the AChE-GF and AChE-GD complexes. Interestingly, appropriate conformations were not found for simulating the reactivation mechanisms with these complexes, which can be explained, for instance, by the steric hindrance observed in the site, thus causing a significant conformational change in the cavity.

Through MLR analysis, we can observe that the combination of interaction energy and reaction energy is sufficient to well explain the experimental data, with a high correlation. However, the mechanistic part has a greater weight and contributes most to the reactivation process through trimedoxime. Therefore, this work will bring about important contributions to the field of drug design and therapies, assisting in the development of broad spectrum and more efficient reactivators.

AChE	Acetylcholinesterase
MmAChE	Mus musculus Acetylcholinesterase
OP	Organophosphorus compounds
ACh BBB GA GB GF GD POX DDVP SER Tyr ARG GLU PHE	Acetylcholine Blood Brain Barrier Tabun Sarin Cyclosarin Soman Paraoxon Dichlorvos Serine Tyrosine Arginine Glutamic acid Phenylalanine

TRP

QM/MM MLR TLC HPLC NMR

MVD

Tryptophan

Quantum Mechanics/Molecular Mechanics Multiple Linear Regression

Thin Layer Chromatography

High-performance liquid chromatography Nuclear Magnetic Resonance

Molegro Virtual Docker

Author Contributions: Alexandre A. de Castro and Bruna T. L. Pereira performed the theoretical calculations, data analysis, elaboration of initial versions of this manuscript and figures preparation; Teodorico C. Ramalho, Elaine F. F. da Cunha and Kamil Kuca contributed in the technical-scientific evaluation of the final version and adjustments of language requirements.

Acknowledgments: The authors wish to thank the Brazilian financial agencies Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq), Fundação de Amparo ao Ensino e Pesquisa de Minas Gerais (FAPEMIG) and Coordenação de Aperfeiçoamento de Pessoal de Nível Superior/Ministério da Defesa (CAPES/MD) for financial support, and the Federal University of Lavras (UFLA) for providing the physical infrastructure and working space. Also supported by excellence project FIM UHK.

Conflicts of Interest: The authors declare no conflict of interest.

Sadik, O.A.; Land, W.H.; Wang, J. Targeting Chemical and Biological Warfare Agents at the Molecular Level. Electroanalysis 2003, 15, 1149–1159.
Ramalho, T.C.; de Castro, A.A.; Silva, D.R.; Silva, M.C.; Franca, T.C.C.; Bennion, B.J.; Kuca, K. Computational Enzymology and Organophosphorus Degrading Enzymes: Promising Approaches Toward Remediation Technologies of Warfare Agents and Pesticides. Med. Chem. 2016, 23, 1041-1061.
Spiers, J. New perspectives on vulnerability using emic and etic approaches. Adv. Nurs. 2000, 31, 715–721.
França, T.C.C.; Silva, G.R.; Castro, A.T. De Chemical Defense: A new subject in the Chemical Teaching. Virtual Química 2010, 2, 84–104.
Chauhan, S.; Chauhan, S.; D’Cruz, R.; Faruqi, S.; Singh, K.K.; Varma, S.; Singh, M.; Karthik, V. Chemical warfare agents. Toxicol. Pharmacol. 2008, 26, 113–122.
Gravett, M.R.; Hopkins, F.B.; Self, A.J.; Webb, A.J.; Timperley, C.M.; Baker, M.J. Evidence of VX nerve agent use from contaminated white mustard plants. R. Soc. A Math. Phys. Eng. Sci. 2014, 470, 20140076–20140076.
Ganesan, K.; Raza, S.; Vijayaraghavan, R. Chemical warfare agents. Pharm. Bioallied Sci. 2010, 2, 166–178.

El-Ebiary, A.A.; Elsharkawy, R.E.; Soliman, N.A.; Soliman, M.A.; Hashem, A.A. N - acetylcysteine in Acute Organophosphorus Pesticide Poisoning: A Randomized, Clinical Trial. Basic Clin. Pharmacol. Toxicol. 2016, 119, 222–227.
Dos Santos, V.M.R.; Donnici, C.L.; DaCosta, J.B.N.; Caixeiro, J.M.R. Organophosphorus pentavalent compounds: history, synthetic methods of preparation and application as insecticides and antitumor agents. Nova 2007, 30, 159–170.
Black, R.M.; Read, R.W. Biological markers of exposure to organophosphorus nerve agents. Toxicol. 2013, 87, 421–437.
Dong, H.; Weng, Y.B.; Zhen, G.S.; Li, F.J.; Jin, A.C.; Liu, ; Pany, S. Clinical emergency treatment of 68 critical patients with severe organophosphorus poisoning and prognosis analysis after rescue. Med. (United States) 2017, 96, 9–12.
Kassa, J.; Korabecny, J.; Nepovimova, E.; Jun, D. The influence of modulators of acetylcholinesterase on the resistance of mice against soman and on the effectiveness of antidotal treatment of soman poisoning in mice. Appl. Biomed. 2016, 1–5.
Worek, F.; Thiermann, H.; Wille, T. Organophosphorus compounds and oximes: a critical review. Toxicol. 2020, 94, 2275–2292.
de Castro, A.A.; Assis, L.C.; Soares, F. V; Kuca, K.; Polisel, D.A.; da Cunha, E.F.F.; Ramalho, T.C. Trends in the Recent Patent Literature on Cholinesterase Reactivators (2016-2019). Biomolecules 2020, 10,
Chambers, J.E.; Meek, E.C. Novel centrally active oxime reactivators of acetylcholinesterase inhibited by surrogates of sarin and VX. Dis. 2020, 133, 104487.
Lorke, D.E.; Petroianu, G.A. The Experimental Oxime K027-A Promising Protector From Organophosphate Pesticide Poisoning. A Review Comparing K027, K048, Pralidoxime, and Obidoxime. Front. Neurosci. 2019, 13, 427.

Herbert, J.; Thiermann, H.; Worek, F.; Wille, T. Precision cut lung slices as test system for candidate therapeutics in organophosphate poisoning. Toxicology 2017, 389, 94–100.
Alencar Filho, E.B.; Santos, A.A.; Oliveira, B.G. A quantum chemical study of molecular properties and QSPR modeling of oximes, amidoximes and hydroxamic acids with nucleophilic activity against toxic organophosphorus agents. Mol. Struct. 2017, 1133, 338–347.
Malfatti, M.A.; Enright, H.A.; Be, N.A.; Kuhn, A.; Hok, S.; McNerney, M.W.; Lao, V.; Nguyen, T.H.; Lightstone, F.C.; Carpenter, T.S.; et al. The biodistribution and pharmacokinetics of the oxime acetylcholinesterase reactivator RS194B in guinea pigs. Chem. Biol. Interact. 2017, 277, 159–167.
Kuča, K.; Kassa, J. A comparison of the ability of a new bispyridinium oxime - 1-(4- hydroxyiminomethylpyridinium)-4-(4-carbamoylpyridinium)butane dibromide and currently used oximes to reactivate nerve agent-inhibited rat brain acetylcholinesterase by in vitro methods. Enzyme Inhib. Med. Chem. 2003, 18, 529–535.
Delfino, R.T.; Ribeiro, T.S.; Figueroa-Villar, J.D. Organophosphorus compounds as chemical warfare agents: a review. Title Toxicol. 2009, 20, 407–428.
Bajgar, J.; Fusek, J.; Kuca, K.; Bartosova, L.; Jun, Treatment of Organophosphate Intoxication Using Cholinesterase Reactivators:Facts and Fiction. Mini-Reviews Med. Chem. 2007, 7, 461–466.
Polisel, D.A.; de Castro, A.A.; Mancini, D.T.; da Cunha, E.F.F.; França, T.C.C.; Ramalho, T.C.; Kuca, K. Slight difference in the isomeric oximes K206 and K203 makes huge difference for the reactivation of organophosphorus-inhibited AChE: Theoretical and experimental aspects. Biol. Interact. 2019, 309, 108671.
Matos, K.S.; Mancini, D.T.; Cunha, E.F.F. da; Kuča, K.; França, T.C.C.; Ramalho, T.C. Molecular aspects of the reactivation process of acetylcholinesterase inhibited by J. Braz. Chem. Soc. 2011, 6, 286–289.
Giacoppo, J.O.S.; França, T.C.C.; Kuča, K.; Da Cunha, E.F.F.; Abagyan, R.; Mancini, D.T.; Ramalho, T.C. Molecular modeling and in vitro reactivation study between the oxime BI-6 and acetylcholinesterase inhibited by different nerve agents. Biomol. Struct. Dyn. 2015, 33, 2048–2058.
Kuca, K.; Musilek, K.; Jun, D.; Zdarova-Karasova, J.; Nepovimova, E.; Soukup, O.; Hrabinova, M.; Mikler, J.; Franca, T.C.C.; Da Cunha, E.F.F.; et al. A newly developed oxime K203 is the most effective reactivator of tabun-inhibited acetylcholinesterase. BMC Pharmacol. Toxicol. 2018, 19, 1–10.
Musilek, K.; Holas, O.; Kuca, K.; Jun, D.; Dohnal, V.; Opletalova, V.; Dolezal, M. Synthesis of monooxime-monocarbamoyl bispyridinium compounds bearing (E)-but-2- ene linker and evaluation of their reactivation activity against tabun- and paraoxon- inhibited acetylcholinesterase. Enzyme Inhib. Med. Chem. 2008, 23, 70–76.
Jun, D.; Stodulka, P.; Kuca, K.; Dolezal, B. High-performance liquid chromatography analysis of by-products and intermediates arising during the synthesis of the acetylcholinesterase reactivator HI-6. Chromatogr. Sci. 2010, 48, 694–696.
Kuča, K.; Cabal, J. Evaluation of Newly Synthesized Reactivators of the Brain Cholinesterase Inhibited by Sarin Nerve Agent. Mech. Methods 2005, 15, 247–252.
J. Frisch, G. W. Trucks, H. B. Schlegel, G. E. Scuseria, M. A. Robb, J. R. Cheeseman, G. Scalmani, V. Barone, G. A. Petersson, H. Nakatsuji, X. Li, M. Caricato, A. Marenich, J. Bloino, B. G. Janesko, R. Gomperts, B. Mennucci, H. P. Hratchian, J. V. Ort, W.C. Gaussian, Inc., Wallingford CT 2009.
Artursson, E.; Andersson, P.O.; Akfur, C.; Linusson, A.; Börjegren, S.; Ekström, F. Catalytic-site conformational equilibrium in nerve-agent adducts of acetylcholinesterase: Possible implications for the HI-6 antidote substrate specificity. Pharmacol. 2013, 85, 1389–1397.

Thomsen, R.; Christensen, M.H. MolDock: A New Technique for High-Accuracy Molecular Docking. Med. Chem. 2006, 49, 3315–3321.
Silva, C.; Pires, M. dos S.; de Castro, A.A.; da Cunha, E.F.F.; Caetano, M.S.; Ramalho,T.C. Molecular insight into the inhibition mechanism of plant and rat 4- hydroxyphenylpyruvate dioxygenase by molecular docking and DFT calculations. Med. Chem. Res. 2015, 24, 3958–3971.

Guimaraes, A.P.; Oliveira, A.A.; da Cunha, E.F.F.; Ramalho, T.C.; Franco, T.C.C. Analysis of Bacillus anthracis nucleoside hydrolase via in silico docking with inhibitors and molecular dynamics simulation. Mol. Model. 2011, 17, 2939–2951.
Matos, K.S.; Mancini, D.T.; da Cunha, E.F.F.; Kuca, K.; Franca, T.C.C.; Ramalho, T.C. Molecular Aspects of the Reactivation Process of Acetylcholinesterase Inhibited by J. Braz. Chem. Soc. 2011, 22, 1999–2004.
da Cunha, E.F.F.; Barbosa, E.F.; Oliveira, A.A.; Ramalho, T.C. Molecular Modeling of Mycobacterium Tuberculosis DNA Gyrase and its Molecular Docking Study with Gatifloxacin Inhibitors. Biomol. Struct. Dyn. 2010, 27, 619–625.
Souza, T.C.S.; Josa, D.; Ramalho, T.C.; Caetano, M.S.; da Cunha, E.F.F. Molecular modelling of Mycobacterium tuberculosis acetolactate synthase catalytic subunit and its molecular docking study with inhibitors. Simul. 2008, 34, 707–713.
Nemukhin, A. V.; Grigorenko, B.L.; Morozov, D.I.; Kochetov, M.S.; Lushchekina, S. V.; Varfolomeev, S.D. On quantum mechanical-molecular mechanical (QM/MM) approaches to model hydrolysis of acetylcholine by acetylcholinesterase. Biol. Interact. 2013, 203, 51–56.
Da Silva Gonçalves, A.; França, T.C.C.; Caetano, M.S.; Ramalho, T.C. Reactivation steps by 2-PAM of tabun-inhibited human acetylcholinesterase: Reducing the computational cost in hybrid QM/MM methods. Biomol. Struct. Dyn. 2014, 32, 301–307.
Matos, K.; Cunha, E.; Abagyan, R.; Ramalho, T. Computational Evidence for the Reactivation Process of Human Acetylcholinesterase Inhibited by Carbamates. Chem. High Throughput Screen. 2013, 17, 554–564.
Heyden, A.; Lin, H.; Truhlar, D.G. Adaptive partitioning in combined quantum mechanical and molecular mechanical calculations of potential energy functions for multiscale simulations. Phys. Chem. B 2007, 111, 2231–41.
Ramalho, T.C.; Caetano, M.S.; da Cunha, E.F.F.; Souza, T.C.S.; Rocha, M.V.J. Construction and Assessment of Reaction Models of Class I EPSP Synthase: Molecular Docking and Density Functional Theoretical Calculations. Biomol. Struct. Dyn. 2009, 27, 195–207.
Kuca, K.; Musilek, K.; Jun, ; Nepovimova, E.; Soukup, O.; Korabecny, J.; França, T.C.C.; de Castro, A.A.; Krejcar, O.; da Cunha, E.F.F.; et al. Oxime K074 – in vitro and in silico reactivation of acetylcholinesterase inhibited by nerve agents and pesticides. Toxin Rev. 2020, 39, 157–166.
da Silva, J.A.; Pereira, A.F.; LaPlante, S.R.; Kuca, K.; Ramalho, T.C.; França, T.C. Reactivation of VX-Inhibited Human Acetylcholinesterase by Deprotonated Pralidoxime. A Complementary Quantum Mechanical Study. 2020, 10, 192.
Besler, B.H.; Merz, K.M.; Kollman, P.A. Atomic charges derived from semiempirical methods. Comput. Chem. 1990, 11, 431–439.
Singh, U.C.; Kollman, P.A. An approach to computing electrostatic charges for molecules. Comput. Chem. 1984, 5, 129–145.
Li, R.; Liu, ; Zhang, J.; Chen, K.; Li, S.; Jiang, J. An isofenphos-methyl hydrolase (Imh) capable of hydrolyzing the P-O-Z moiety of organophosphorus pesticides containing an aryl or heterocyclic group. Appl. Microbiol. Biotechnol. 2012, 94, 1553–1564.
Gorecki, L.; Korabecny, J.; Musilek, K.; Malinak, D.; Nepovimova, E.; Dolezal, R.; Jun, ;Soukup, O.; Kuca, K. SAR study to find optimal cholinesterase reactivator against organophosphorous nerve agents and pesticides. Arch. Toxicol. 2016, 90, 2831–2859.

da Cunha, E.F.F.; Ramalho, T.C.; Reynolds, R.C. Binding mode analysis of 2, 4-diamino- 5-methyl-5-deaza-6-substituted pteridines with mycobacterium tuberculosis and human dihydrofolate reductases. Biomol. Struct. Dyn. 2008, 25, 377–385.
da Silva, J.A.V.; Pereira, A.F.; LaPlante, S.R.; Kuca, K.; Ramalho, T.C.; França, T. C. C. Biomolecules. 2020, 10,

Understand the interaction modes and reactivity of trimedoxime toward MmAChE inhibited by nerve agents: theoretical and experimental Aspects

Status:

Version 1

Abstract

Figures

Introduction

Results And Discussion

2.2 Affinity and thermodynamics: docking results

2.3 Investigating kinetic parameters for biological activity: mechanistic studies

Materials And Methods

3.2 Docking procedure

3.3 QM/MM procedure

Conclusions

Abbreviations

Declarations

References

Status:

Version 1