# Steady-State Kinetics of Enzyme-Catalyzed Hydrolysis of Echothiophate, a P-S Bonded Organophosphorus as Monitored by Spectrofluorimetry

Irina Zueva, Sofya Lushchekina, David Daudé, Eric Chabriere, Patrick Masson

► To cite this version:

Irina Zueva, Sofya Lushchekina, David Daudé, Eric Chabriere, Patrick Masson. Steady-State Kinetics of Enzyme-Catalyzed Hydrolysis of Echothiophate, a P-S Bonded Organophosphorus as Monitored by Spectrofluorimetry. Molecules, 2020, 25 (6), pp.1371. 10.3390/molecules25061371. hal-02517839

HAL Id: hal-02517839

https://amu.hal.science/hal-02517839

Submitted on 27 Jun 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Article

# Steady-State Kinetics of Enzyme-Catalyzed Hydrolysis of Echothiophate, a P–S Bonded Organophosphorus as Monitored by Spectrofluorimetry

Irina V. Zueva<fnref n="1" />, Sofya V. Lushchekina<fnref n="2" /> ID, David Daudé<fnref n="3" />, Eric Chabrière<fnref n="4" /> and Patrick Masson<fnref n="5" />,*

<figure><img src="image_2.png" /></figure>

check for
updates

Received: 22 February 2020; Accepted: 16 March 2020; Published: 17 March 2020

**Abstract:** Enzyme-catalyzed hydrolysis of echothiophate, a P–S bonded organophosphorus (OP) model, was spectrofluorimetrically monitored, using Calbiochem Probe IV as the thiol reagent. OP hydrolases were: the G117H mutant of human butyrylcholinesterase capable of hydrolyzing OPs, and a multiple mutant of *Brevundimonas diminuta* phosphotriesterase, GG1, designed to hydrolyze a large spectrum of OPs at high rate, including V agents. Molecular modeling of interaction between Probe IV and OP hydrolases (G117H butyrylcholinesterase, GG1, wild types of *Brevundimonas diminuta* and *Sulfolobus solfataricus* phosphotriesterases, and human paraoxonase-1) was performed. The high sensitivity of the method allowed steady-state kinetic analysis of echothiophate hydrolysis by highly purified G117H butyrylcholinesterase concentration as low as 0.85 nM. Hydrolysis was Michaelian with K_m = 0.20 ± 0.03 mM and k_cat = 5.4 ± 1.6 min⁻<fnref n="1" />. The GG1 phosphotriesterase hydrolyzed echothiophate with a high efficiency (K_m = 2.6 ± 0.2 mM; k_cat = 53,400 min⁻<fnref n="1" />). With a k_cat/K_m = (2.6 ± 1.6) × 10⁷ M⁻¹min⁻<fnref n="1" />, GG1 fulfills the required condition of potential catalytic bioscavengers. Quantum mechanics/molecular mechanics (QM/MM) and molecular docking indicate that Probe IV does not interact significantly with the selected phosphotriesterases. Moreover, results on G117H mutant show that Probe IV does not inhibit butyrylcholinesterase. Therefore, Probe IV can be recommended for monitoring hydrolysis of P–S bonded OPs by thiol-free OP hydrolases.

**Keywords:** P–S bonded organophosphorus agents; echothiophate; Calbiochem Probe IV; organophosphate hydrolase; phosphotriesterase; cholinesterase; QM/MM

## 1. Introduction

Organophosphorus agents (OP) are among the most toxic chemicals. These compounds have been used for decades as pesticides or drugs, and the most dreadful are banned chemical warfare agents. They primarily inhibit acetylcholinesterase (AChE) and butyrylcholinesterase (BChE), causing a major cholinergic syndrome [1,2].

<footnote n="1">¹ Arbuzov Institute of Organic and Physical Chemistry, Federal Research Center “Kazan Scientific Center of the Russian Academy of Sciences”, Arbuzov str. 8, 420088 Kazan, Russia; zueva.irina.vladimirovna@gmail.com</footnote>

<footnote n="2">² Emanuel Institute of Biochemical Physics, Russian Academy of Sciences, Kosygin str 4, 119334 Moscow, Russia; sofya.lushchekina@gmail.com</footnote>

<footnote n="3">³ Gene&GreenTK, HU Méditerranée Infection, Jean Moulin Blvd 19-21, 13005 Marseille, France; david.daude@gene-greentk.com</footnote>

<footnote n="4">⁴ Aix-Marseille University, IRD, APHM, MEPHI, IHU-Méditerranée Infection, 15005 Marseille, France; eric.chabriere@univ-amu.fr</footnote>

<footnote n="5">⁵ Kazan Federal University, Neuropharmacology Laboratory, Kremlevskaya str 18, 480002 Kazan, Russia</footnote>

<footnote>* Correspondence: pym.masson@free.fr; Tel.: +7-96-5581-0473</footnote>
Several classes of OPs can be described according to atoms/chains bonded to the phosphorus atom. Among them, P–S bonded phosphoro- and phosphono-thioates are of interest (Scheme 1).

<figure><img src="image_3.png" /><figcaption>(a)</figcaption></figure>

<figure><img src="image_4.png" /><figcaption>(b)</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>(c)</figcaption></figure>

**Scheme 1.** Structure of three P–S bonded OPs: (a), echothiophate iodide; (b) malathion; (c) $S_p$(-)VX ([2-(diisopropylamino)ethyl]-O-ethyl methylphosphonothioate). Compound “a”, a phosphorothioate, was used as a P–S bonded model OP in the present study.

For example, echothiophate (phospholine iodide) is a drug used for treatment of glaucoma [3], malathion (carbophos) is an agriculture pesticide for eradication of mosquitos also used for treatment of pediculosis [4,5], and VX is a potent chemical warfare nerve agent [6]. Spontaneous hydrolysis or enzyme-catalyzed hydrolysis of phosphoro/phosphono-thioates and phosphorylation of ChEs by these compounds results in breakage of the P–S bond in the P–S-alkyl chain with release of thiol leaving group. There are two fundamentally different mechanisms for enzyme-catalyzed hydrolysis of OPs, for ChE mutants and phosphotriesterases (PTE) [7]. Because ChEs are serine hydrolases, OP hydrolysis catalyzed by ChE mutants involves a two chemical-step mechanisms, as for hydrolysis of carboxylesters: after formation of the productive enzyme-substrate complex the nucleophilic attack of the catalytic serine leads to formation of the phosphyl-intermediate accompanied by release of the oxo- or thio-leaving group, and then water-mediated hydrolysis of this intermediate takes place. On the other hand, PTEs are metalloenzymes and the hydrolysis mechanism involves only one chemical step. Although intimate mechanisms of PTE are still debated, there is a consensus, stating that hydrolysis results from a direct attack on P atom of a water-derived hydroxyl group bridging metal cations [8].

Therefore, hydrolytic degradation of phosphoro/phosphono-thioates can be monitored by measuring the rate of formation of the thiol leaving group. There are numerous chromogenic thiol probes, one of the most widely used is the Ellman's reagent, DTNB, (5,5'-dithiobis-2-nitro) benzoic acid, that releases a yellow anion, 5-thio-2-nitrobenzoate, upon reaction with thiol compounds. This reaction is implemented in the most popular method for assay of ChEs with thiocholine esters as substrates [9]. Because hydrolysis of echothiophate releases also thiocholine (cf. Scheme 1, a), this reaction has been used to monitor the hydrolysis of this OP by the G117H mutant of human BChE and by the mutant HQF of *Bungarus fasciatus* AChE, mutants capable of hydrolyzing OPs [10,11].

However, rather high concentrations of enzymes are needed for kinetic analysis of P–S bonded hydrolysis of OPs. A fluorogenic thiol probe would be more sensitive than DTNB. The use of the fluorogenic Calbiochem Probe IV ((3-(7-Hydroxy-2-oxo-2H-chromen-3-ylcarbamoyl) acrylic acid
methylester) (Scheme 2) was proven to be 2 orders of magnitude more sensitive than the Ellman assay for kinetic analysis of ChEs-catalyzed hydrolysis of thiocholine esters [12].

<figure><img src="image_6.png" /><figcaption>**Scheme 2.** Reaction of Probe IV with a thiol compound (R-SH) giving a highly fluorescent conjugate.</figcaption></figure>

In the present work, we report steady-state kinetic analysis of degradation a P-S-bonded OP by two OP hydrolases, using Probe IV as thiol probe. OP hydrolases are of great biotechnological and medical interest. Phosphotriesterases (PTE) of various origins, in particular, can be used for detection of OPs [13], decontamination and remediation [8,14–18] and for therapy of OP poisoning as catalytic bioscavengers [7,19–22]. Because OPs are hemi-substrates of ChEs, research of novel human ChE mutants capable of hydrolyzing OPs is of great interest for improving catalytic bioscavenger-based medical countermeasures of OP poisoning [23–25] and implementation of pseudocatalytic bioscavenger systems composed of ChE-reactivator [26–28]. Enzymatic hydrolysis of P-S bonded OPs, such as V-agents, is currently monitored by the reaction of the OP thiol leaving group with DTNB [16,29]. However, DTNB is two orders of magnitude less sensitive than Probe IV.

The use of Probe IV makes it possible to push the limit of detection of P-S bonded OPs. The P-S bonded model OP in the present study was echothiophate iodide, and the enzymes were the reference G117H mutant of human BChE that hydrolyze OPs, and a new multiple mutant of *Brevundimonas diminuta* PTE, GG1, specially designed for rapid hydrolytic detoxification of V-agents [30].

## 2. Results and Discussion

### 2.1. *Interaction of Probe IV With Enzymes*

#### Molecular Docking Studies

Molecular docking was used to provide information about possible reversible binding of Probe IV to ChEs in our previous work [12]. However, molecular docking techniques have intrinsic limitations [31]. In particular, in the case of ChEs, due to the complex architecture of the deep active site gorge, direct transformation of docking binding affinities into inhibition constants should be used with caution. In the case of metallo-PTEs, a similar warning has to be pointed out. Nevertheless, molecular docking often provides useful clues about molecular interactions and it is helpful for comparative studies, as in the present report. Furthermore, there is a body of docking studies about numerous ligands into ChEs, coupled with experimental measurements, serving refinement of this approach [32,33]. Also, substantial efforts have been made to parametrize metal-containing systems [34,35] to provide adequate docking score.

Molecular docking showed that Probe IV can bind to active sites of all considered enzymes. For human BChE, mutation G117H slightly alters the position of Probe IV in the BchE active site (Figure 1a compared to the wild-type enzyme [12]). This causes a moderate increase of estimated weak binding affinity (−5.76 kcal/mol vs. −6.61 kcal/mol). Thus, we may consider that 1 µM Probe IV does not inhibit G117H.

Very close docked binding affinities were obtained for Probe IV with PON-1 (−6.75 kcal/mol) (Figure 1b) and wild-type *Brevundimonas diminuta* PTE (−6.41 kcal/mol). However, for the PTE mutant GG1 binding was much weaker (−4.16 kcal/mol). Comparison of binding position of Probe IV in active sites of wild-type and GG1 *B. diminuta* PTE (Figure 1c,d) shows that in the latter case Probe IV forms
fewer hydrogen bonds. This explains drop in estimated binding free energies. Thus, inhibition of these PTEs would certainly be significant above 10–20 µM Probe IV.

Reaction between Probe IV and Cys106 in GG1 was found to be unlikely, though possible: the minimal distance between sulfur atom and reacting atom of Probe IV among all docked poses was 3.8 Å. This is too long for covalent reaction, but could be got over due to thermal motions.

Binding position of Probe IV in the active site of *Sulfolobus solfataricus* PTE with QM/MM-derived point charges (Figure 1e,f) was closer to position in *B. diminuta* PTE, but binding affinity was weaker (−5.34 kcal/mol). Docking with formal charges +2 on metals shows the importance of accurate parameterization of bi-cations in active sites (without this step, binding affinities were strongly overestimated: −13.4 kcal/mol for *B. diminuta* PTE and −17.10 kcal/mol for *S. solfataricus* PTEs).

<figure><img src="image_7.png" /><figcaption>Figure 1. *Cont.*</figcaption></figure>
<figure><img src="image_8.png" /></figure>

**Figure 1.** Docked positions of Probe IV inside active sites of (a) G117H butyrylcholinesterase (BChE), (b) PON-1, (c) wt *B. diminuta* phosphotriesterases (PTE), and (d) GG1 mutant of *B. diminuta* PTE. In panel (d) reacting atom of Probe IV is highlighted yellow and (e) *S. solfataricus* PTE.

## 2.2. Steady-State Kinetics of Enzyme-Catalyzed Hydrolysis of Echothiophate

### 2.2.1. Hydrolysis of Echothiophate by the G117H Mutant of Human Butyrylcholinesterase

Steady-state kinetics analysis of G117H-catalyzed hydrolysis of echothiophate up to 0.6 mM (Figure 2), under pseudo-first order conditions ([E]<<[OP]) was interpreted in terms of Michaelis–Menten model. However, at higher concentration, binding of a second echothiophate molecule on the enzyme peripheral anionic site would be expected, thus, altering the reactivity of the catalytic center as do positively charged substrates and OPs at high concentration [36,37].

<figure><img src="image_9.png" /></figure>

**Figure 2.** Steady state kinetics of echothiophate hydrolysis by the G117H mutant of human BChE ([E] = 8.5 nM) in 0.1M phosphate buffer pH 7.0 at 25 °C. Fluorescence change/min was converted in micromoles of released thiol/min, using calibration plot established in proper buffer.

The calculated catalytic parameters are: $K_m$ = 0.20 ± 0.03mM, $k_{cat}$ = 5.4 ± 1.6 min⁻¹, and $k_{cat}/K_m$ = 2.7 ± 1.6 × 10⁴ M⁻¹min⁻¹. Though slightly higher, these values are in accordance with previously reported parameters for this mutant with echothiophate under the same conditions: $K_m$ = 0.074 mM, $k_{cat}$ = 0.75 min⁻¹, and $k_{cat}/K_m$ = 1.01 × 10⁴ M⁻¹min⁻¹ [10]. It should be noted that these old results were obtained, using a partially purified enzyme (<90% pure) at enzyme concentration
of 3 × 10⁻⁶ M per assay. The high sensitivity of Probe IV allowed us to determine the catalytic parameters of G117H with an enzyme more than three–four orders of magnitude less concentrated (0.85–28.5 × 10⁻⁹ M).

However, our results showed that 1% dimethylsulfoxide (DMSO) inhibits the enzyme. This is in accordance with previously reported data on AChE and BChE [12,38]. Because solvent inhibiting effects are independent on the nature of substrate, the effect of DMSO on G117H was investigated, using the classical Ellman's method with butyrylthiocholine (BTC) as the substrate. Thus, it was possible to compare the present results on G117H mutant with results previously obtained on wild type and other mutants of human BChE [12,38]. Kinetic analysis of DMSO-induced inhibition of G117H activity toward BTC showed that DMSO is a weak reversible competitive inhibitor with Kᵢ = 0.6%, i.e., 76 mM (Figure 3). Thus, 1% DMSO (= 128 mM) leads to a K_m increment 1 + [I] / Kᵢ = 2.68. This value is consistent with the observed 2.7-times increase of apparent K_m of G117H for echothiophate, and three-times increase of apparent K_m of wild-type BChE and E197G and E197D mutants for BTC [12]. Also, the fact that weak reversible inhibition of the enzyme by DMSO fully explains the slight increase in K_m of our results compared to previously reported results [10], indicates that ProbeIV does not compete with echothiophate in G117H enzyme active site.

<figure><img src="image_10.png" /><figcaption>**Figure 3.** Reversible inhibition of G117H by DMSO, (a), Dixon plot; (b), Cornish-Bowden plot.</figcaption></figure>

### 2.2.2. Hydrolysis of Echothiophate by the GG1 Mutant of Brevundimonas Diminuta PTE

The GG1 mutant of *Br. d.* PTE was found to be highly effective to hydrolyze echothiophate with K_m = 2.6 ± 0.2 mM; *k*_{cat} = 55400 min⁻¹ and *k*_{cat}/K_m = (2.1 ± 1.6) × 10⁷ M⁻¹min⁻¹ (Figure 4). These values are in agreement with reported values for related multiple mutants of *Br. d.* PTE with other P–S bonded OPs.

Moreover, assay of the enzyme with 500 µM echothiophate in the presence of various concentrations of DMSO (1%, 1.5%, 2%, and 2.5%) showed no activity change. This indicates that the catalytic activity of GG1 is not altered by the presence of DMSO up to 2.5% (not shown).
<figure><img src="image_11.png" /><figcaption>Figure 4. Steady-state GG1-catalyzed hydrolysis of echothiophate. [E] = 0.024 µM. Fluorescence change/min was converted in micromoles of released thiol/min, using calibration plot established in proper buffer.</figcaption></figure>

Although no data are available with echothiophate for other Br.d. PTE mutants that cleave P–S bond, it was reported by Raushel’s group that different multiple mutants of this enzyme capable of hydrolyzing P–S bonded pesticides and V-agents display similar properties, showing the highest specificity constants for malathion: 1 × 10⁷ M⁻¹min⁻¹ [16], Sₚ(VX): 5.1 × 10⁷ M⁻¹min⁻¹, and Sₚ(VR): 1.2 × 10⁷ M⁻¹min⁻¹ [8]. These reported kinetic values were determined by using diluted enzymes whose concentrations typically ranged between 8 × 10⁻⁸ M (for the most active mutants) and 3 × 10⁻⁵ M (for wt) [39]. Thus, with a high specificity constant (1.8 × 10⁷ M⁻¹min⁻¹) for echothiophate as the substrate, the multiple mutant GG1 almost fulfill the minimal required condition (kₐt/Kₘ > 5 × 10⁷ M⁻¹min⁻¹) to be considered as a potential catalytic bioscavenger [7,19].

This statement has important practical therapeutic applications, regardless galenic formulation and safety issues. Indeed, even in the most severe cases of poisoning the OP concentration in blood is extremely low (e.g., ≈ 3 × 10⁻⁶ M), in vivo enzymatic detoxification of OP by a catalytic bioscavenger operates under first order conditions, according to Equation (3), with the first order rate constant equal to (kₐt/Kₘ) [E], and:

$$[OP]_t = [OP]_0 \exp\left(-\frac{k_{cat}[E]_0}{K_m} \cdot t\right). \qquad (1)$$

Thus, to drop a concentration in echothiophate of 3 × 10⁻⁶ M by 100-fold (ln[OP]₀/[OP]ₜ = 4.6) in human blood in 1min (median time for a full turn of blood circulation), an amount of 9 mg of GG1 would be needed to inject to a human of 70 kg. Such a dose is about 20-fold lower than the equivalent dose that would be needed with a stoichiometric bioscavenger like human BChE [40,41]. Decontamination of biological samples or fragile materials by GG1 would also be very effective and fast. Although PTE-based prophylactic and post-exposure treatments would require further developments, including encapsulation in nanoparticles and pharmacokinetics studies, such an enzymatic variant is of prime interest for external decontamination (e.g., skin, materials) for limiting poisoning effects of OPs in humans and environment. On the other hand, the specificity constant (kₐt/Kₘ) of ChE-based mutants capable of self-reactivation after phosphorylation would have to be increased by four orders of magnitude to reach the >5 × 10⁷ M⁻¹min⁻¹ condition. This challenge has already been pointed out [24]. However, progress in directed evolution and computer re-design of enzymes make possible to take up the challenge.
# 3. Materials and Methods

## 3.1. *Chemicals*

Echothiophate iodide was from Biobasal AG (Basel, Switzerland). Calbiochem Probe IV (3-(7-Hydroxy-2-oxo-2H-chromen-3-ylcarbamoyl) acrylic acid methylester) was from Merck Millipore (Darmstadt, Germany). Probe IV is a light sensitive chemical. Glutathione (GSH) used for calibration curve was purchased from Sigma-Aldrich (Saint Louis, MO, USA). All other chemicals were of biochemical grade.

## 3.2. *Enzymes*

G117 H mutant of human BChE: Low glycosylated and truncated monomeric recombinant G117H mutant of human BChE (MW = 70,000) was produced by O. Lockridge in Chinese hamster ovary (CHO) cells and highly purified for X ray structure determination [42]. The enzyme appeared 100% pure on SDS-4-30% polyacrylamide gradient gel; it was concentrated to 6 mg/mL ([E] = 0.85 × 10⁻⁴ M) in 10 mM MES pH 6.5, containing 0.02% sodium azide. The enzyme was stored for years at -80 °C without loss of ChE activity as assayed with 1 mM butyrylthiocholine iodide as the substrate (200 unit/mg) [30].

GG1 mutant of *Brevundimonas diminuta* PTE: Multiple mutated PTE (MW = 36,000) was designed to hydrolyze V agents in addition to other OPs. The variant GG1 was previously optimized through ancestral engineering and shown to degrade V-agent surrogates [30]. The gene coding for the variant GG1 from *Brevundimonas diminuta* containing N-terminal Strep-TEV tag was cloned in a pET22b+ vector using NdeI/NotI restriction sites. Heterologous production was performed in *Escherichia coli* BL21 (DE3)-pGro7/GroEL (TaKaRa) chaperone expressing strain as previously described. Briefly, production was performed in ZYP medium (10 g/L Tryptone, 5 g/L Yeast Extract, 25 mM (NH4)2SO4, 100 mM KH2PO4, 50 mM Na2HPO4, 0.5% (w/v) glycerol, 0.05% (w/v) glucose 25 g, 0.2% (w/v) α-lactose, 100 µg/mL ampicillin and 34 µg/mL chloramphenicol) at 37 °C and 140 rpm shaking. When OD600nm reached 0.8–1, CoCl2 (final concentration 0.2 mM) and L-arabinose (final concentration 2 g·L⁻¹) were added, and the temperature was decreased to 16 °C for 16–20 h. Cells were harvested by centrifugation (4400 g, 10 °C, 30 min). Pellets were re-suspended in buffer (50 mM Tris, 300 mM NaCl, pH 8.0). Then, cell lysis was performed by adding DNAseI, lysozyme and PMSF to reach final concentrations of 10 µg·mL⁻¹, 0.25 mg·mL⁻¹ and 0.1 mM, respectively. Subsequently, three steps of 30 s of sonifications (Qsonica, Q700; Amplitude 45) were applied in ice. Cell debris were removed by centrifugation (12,000 g, 10 °C, 20 min) and supernatant was filtered (0.8 µM) prior to purification. Purification was performed by affinity chromatography (StrepTrapTM HP GE Healthcare; ÄKTA Avant) and eluted in 50 mM Tris, 300 mM NaCl, 2.5 mM desthiobiotin, pH 8.0. Protein concentration was measured by using a NanoDrop 2000 (Thermo Scientific) spectrophotometer. Then, the protein solution was adjusted to 1 mg/mL with 20 mM Tris buffer, pH 8.0 supplemented with 0.2 mM CoCl2. Enzyme purity (85%, corresponding to ([E] = 2.43 × 10⁻⁵ M) was checked by SDS-PAGE (T = 12.5%) followed by Coomasie Brilliant Blue staining (Figure 5).
<figure><img src="image_12.png" /><figcaption>**Figure 5.** SDS-PAGE of purified GG1 mutant of *Brevundimonas diminuta* PTE. The band of 28 kDa was excised from the gel manually and identified as a degraded form of GG1 by peptide mass fingerprinting using a MALDIT-OF MS Bruker Ultraflex I spectrometer (Bruker Daltonics) after trypsin digestion as previously described [43].</figcaption></figure>

### 3.3. Reaction of Probe IV With Thiols and Calibration

“Calbiochem probe IV”, is a coumarinyl derivative that reacts with thiol-chemicals to form a highly fluorescent conjugate. It was found to be the fastest and the most sensitive thiol reagent with a large signal-to-noise ratio [44]. It was, in particular, applied for screening human BChE clones with butyrylthiocholine as a substrate in microfluidic droplets [45]. More recently, Probe IV was used for quantitative measurements of low activity of BChE and AChE and for kinetic analysis of ChE-catalyzed hydrolysis of thioesters [12].

Stock solution of 0.1 M echothiophate was in water and stored at -20 °C. Stock solution of Probe IV was 1mM in DMSO and stored at -20 °C. Because fluorescence of Probe IV-thio-conjugate was found dependent on buffer, 60 µM stock solution of glutathione (GSH), the calibration thiol compound, was in 20 mM Tris/HCl pH 8.0 for calibration of velocity data of GG1 and in 0.1 M sodium phosphate buffer pH 7.0 for calibration of velocity data of G117H mutant.

Kinetic measurements were performed in a Peltier thermostated spectrofluorimeter F-7100 (Hitachi Ltd., Japan), using standard spectrofluorimetric cuvettes of 1 cm-path length in a total volume of 1.5 mL. The fluorescence emission of Probe IV-thiol conjugate has $\lambda_{ex}$ = 400 nm and $\lambda_{em}$ = 465 nm. $\lambda_{ex}$ slit was 5.0 nm and $\lambda_{em}$ slit 10.0 nm. Calibration plot was made with glutathione up to 5 µM final in cuvette containing either 0.1 M phosphate buffer pH 7.0 (calibration for G117H activity) or 20 mM Tris/HCl pH 8.0 (calibration for GG1 activity). Fluorescence intensity plot, FI vs [GSH] was built after subtraction of the intrinsic fluorescence of Probe IV.

The buffer for the G117H mutant of human ChE was 0.1 M sodium phosphate buffer pH 7.0, and 20 mM Tris/HCl pH 8.0 for GG1. The final concentration of Probe IV in cuvette was 1 µM, and the final DMSO concentration was 1% *v/v*. The final echothiophate concentration in cuvette ranged from 0.05 to 0.8 mM. The final concentration in active sites per assay ranged from 8.51 × 10⁻¹⁰ M to 2.85 × 10⁻⁸ M for G117H mutant and from 1.6 × 10⁻⁸ to 9.4 × 10⁻⁸ M for GG1. The temperature for kinetic measurements was set up at 25 °C for the G117H mutant of BChE and at 30 °C for the GG1 mutant of *Brevundimonas diminuta*. PTE.

The rate of GG1-catalyzed hydrolysis (ΔIF/dt) of echothiophate was monitored for 8 min, while for G117H mutant of BChE it was 2 min. The fluorescence background of Probe IV was subtracted, as well as the fluorescence background due to spontaneous hydrolysis of echothiophate.
### 3.4. Steady-State Kinetic Study of Enzyme-Catalyzed Hydrolysis of Echthiophate

Steady-state kinetic measurements of initial velocity for enzyme-catalyzed hydrolysis of echthiophate were repeated in triplicate for each enzyme concentration used. Values of ΔF/min were converted in µmol of thiocholine released/min using the GSH calibration plot.

Previously reported results for G117H showed that at least up to 1 mM echthiophate, the enzyme obeys the Michaelis–Menten kinetic model [10,11]. Also, all reported studies performed on wild type and mutants of *Brevundimonas diminuta* showed that these enzymes obey the Michaelis–Menten model regardless the substrate [39]. Thus, catalytic parameters $V_{max}$ and $K_m$ were determined by non-linear fitting of initial velocity to the Michaelis–Menten rate equation (Equation (2)), using Origin (OriginLab Co, Northampton, MA, USA).

$$v = \frac{k_{cat}[E][S]}{K_m + [S]} \qquad (2)$$

At very low echthiophate concentration, [OP] << $K_m$, the bimolecular rate constants $k_{cat}/K_m$ were also determined by linear regression of Equation (3).

$$v = \frac{k_{cat}[E][S]}{K_m} \qquad (3)$$

### 3.5. Possibility of Unwanted Interactions of DMSO and Probe IV With Enzymes

As previously suggested for wild-type BChE [12], the possibility that 1% DMSO, the solvent of Probe IV, could reversibly inhibit the G117H mutant was checked. The enzyme was assayed by the method of Ellman [9] with BTC (0.1, 0.2, and 0.5 mM) as the substrate and DTNB as the chromogene, in the presence of different concentrations of DMSO (0.5%, 1%, 1.5%, and 2% final). The effect of DMSO on the enzyme was analyzed by Dixon and Cornish-Bowden plots [46].

The possibility that Probe IV could be a reversible inhibitor of G117H mutant was checked. Indeed, Probe IV bears structural components in common with certain fluorescent probes that have been recently recognized has reversible inhibitors of BChE [47].

The possibility that Probe IV could alkylate free thiols was investigated. G117H BChE has no solvent exposed cysteine [42]. On the other hand, GG1 has two naturally occurring cysteines and a residue I106 mutated in cysteine (I106C). This later is located in the active site pocket. For this purpose, 0.047 µM GG1 were pre-incubated from 5 to 20 min in dark with 1 µM in 1% DMSO. After the different pre-incubation times, 0.5 nM echthiophate was added and the enzyme activity measured.

### 3.6. Molecular Modeling

In a previous work, we showed that Probe IV does not inhibit significantly human BChE [12]. Thus, to estimate the possibility of reversible inhibition and/or irreversible inhibition of PTE by Probe IV molecular docking of Probe IV was performed on *Brevundimonas diminuta* PTE. It should be noted that wt PTE possesses 2 free Cysteines in its structure, and the GG1 mutant has an additional Cys in its active center (I106C). A reported result on closely related mutants showed that the two natural Cys do not react with DTNB as the thiol probe, and that C106 is also inaccessible to DTNB [16].

Although the different PTEs has common features in their active site, they have different 3D structures (TIM-barrel fold such as *Brevenundimonas diminuta* PTE, pita bread-fold, 6-propeller-fold: [48], our work was expanded to docking of Probe IV to paraoxonase, a 6-propeller-fold enzyme and *Sulfolobus solfataricus*, another PLL- PTE of TIM-barrel type. Docking into active sites of PTEs is challenging due to the presence of bi-cations and their coordination by active site residues, leading to charge redistribution [49]. To overcome this problem, one solution is to use specially developed charge-independent AutoDock4Zn force field [34] for Zn-containing *B. diminuta* PTE. In case of Co and Fe-containing *S. solfataricus* PTE solution is to use QM-derived partial charges on active site atoms. Due the protein size, structure optimization requires combined QM/MM method [50].
The structure of Probe IV was quantum-mechanically optimized with GameSS-US [51] software (B3LYP/6-31G*). According to previous findings [33], partial atomic charges assigned according to the Löwdin scheme [52] from quantum mechanics data was used for molecular docking.

X-ray structures of G117H BChE (PDB ID 2XMD [42]); *Brevundimonas diminuta* PTE (PDB ID 1I0D [53] with both Zn²⁺ cations) and GG1 mutant model based on it; PON-1 (PDB ID 1V04 [54]); *Sulfolobus solfataricus* PTE (PDB ID 3UF9 [55]) were used. All water and co-crystallized molecules were removed from all structures; active site bi-cations were kept in PTEs. Hydrogen atoms were added using Reduce software [56].

For QM/MM optimization of *Sulfolobus solfataricus* PTE, we proceeded as follows: protein was surrounded with 1006 water molecules, forming solvation shell, total size of the system was 8059 atoms (Figure 6). Active site residues and 7 water molecules were included into quantum subsystem 137 atoms total). NwChem 6.5 [57] software was used (PBE0-D3/cc-pvdz/Amber). The effective core potentials (ECPs) with LanL2DZ basis set were used to describe Co²⁺ and Fe²⁺ cations. For optimized active site geometry, partial atomic charges were derived in the same fashion as for the ligand.

<figure><img src="image_13.png" /><figcaption>**Figure 6.** Full *Sulfolobus solfataricus* PTE structure in solvation shell used for QM/MM modeling, QM subsystem with active site is shown in the circle.</figcaption></figure>

For molecular docking AutoDock 4.2.6 software [58] was used with AutoDock4Zn force field for *B. diminuta* PTE [34]. For all proteins grid boxes were 22.5 Å × 22.5 Å × 22.5 Å size (grid spacing 0.375 Å) and embraced active sites of the enzymes. The main Lamarckian Genetic Algorithm (LGA) [59] parameters were set as follows: 256 runs, 25 × 10⁶ evaluations, 27 × 10⁴ generations, and a population size of 3000.

**4. Conclusions**

The simplicity, accuracy, and sensitivity of Probe IV-based spectrofluorimetric assay for investigating steady-state kinetics of enzyme catalyzed hydrolysis of echothiophate as a P-S bonded organophosphorus model make this new analytical method of great interest for research of novel OP hydrolases capable of hydrolyzing V-agents and easy determination of their catalytic properties.
In particular, it could advantageously be used for high throughput direct microfluidic screening of recombinant OP hydrolase clones [45]. The high sensitivity of Probe IV is an asset to screen low activities at the beginning of directed evolution cycles. Thus, Probe IV-based assays would be particularly suitable for rapid screening of evolved enzyme libraries, enzyme mining in various biotopes, and for enzyme-based biosensors with improved detection limit of compounds of interest. However, possible covalent binding and non-covalent interactions of Probe IV with enzymes may affect the catalytic parameters of OP hydrolases. Thus, quantitative studies must be supported by thorough molecular modeling and non-spectro(fluorimetric) determinations of binding affinity, e.g., isothermal calorimetry, as well as kinetic study of the reactivity of free cysteines present in enzyme structures. In particular, it should be demonstrated that concentrations of the thiol probe used in activity assays does not affect the measured values of P-S bonded OP hydrolyzing enzyme activity.

**Author Contributions:** I.V.Z. performed kinetic studies; S.V.L. performed molecular modeling studies; D.D. and E.C. produced the GG1 mutant of *Brevundimonas diminuta*; P.M. wrote the manuscript. All authors have read and agreed to the published version of the manuscript.

**Funding:** Biochemical experiments performed by I.Z. were supported by the Russian government assignment for FRC Kazan Scientific Center of RAS. Molecular modeling part of this study was funded by RFBR according to the research project #19-03-00043 to S.L.

**Acknowledgments:** The authors are indebted to O. L. (UNMC, Omaha, NE, USA) for the G117H mutant of human BChE and for valuable discussions. The authors are grateful to Pauline Jacquet and Lucile Pinault for their help with GG1 production and purification. Computer modeling was carried out using equipment from the shared research facilities of the HPC computing resources at Lomonosov Moscow State University [60].

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

## References

1. Eddleston, M.; Buckley, N.A.; Eyer, P.; Dawson, A.H. Management of acute organophosphorus pesticide poisoning. *Lancet* **2008**, *371*, 597-607. [CrossRef]

2. Chowdhary, S.; Bhattacharyya, R.; Banerjee, D. Acute organophosphorus poisoning. *Clin. Chim. Acta* **2014**, *431*, 66-76. [CrossRef]

3. Schmidt, K.G.; Horowitz, Y.; Buckman, G.; Segev, E.; Levinger, E.; Geyer, O. Lowering of IOP by echothiophate iodide in pseudophakic eyes with glaucoma. *Curr. Eye. Res.* **2010**, *35*, 698-702. [CrossRef] [PubMed]

4. Bonner, M.R.; Coble, J.; Blair, A.; Beane Freeman, L.E.; Hoppin, J.A.; Sandler, D.P.; Alavanja, M.C. Malathion exposure and the incidence of cancer in the agricultural health study. *Am. J. Epidemiol.* **2007**, *166*, 1023-1034. [CrossRef] [PubMed]

5. Idriss, S.; Levitt, J. Malathion for head lice and scabies: treatment and safety considerations. *J. Drugs Dermatol.* **2009**, *8*, 715-720. [PubMed]

6. Reiter, G.; Muller, S.; Hill, I.; Weatherby, K.; Thiermann, H.; Worek, F.; Mikler, J. In vitro and in vivo toxicological studies of V nerve agents: Molecular and stereoselective aspects. *Toxicol. Lett.* **2015**, *232*, 438-448. [CrossRef] [PubMed]

7. Masson, P.; Luschekina, S. Catalytic Bioscavengers: The second Generation of Bioscavenger-Based Medical Countermeasures. In *Handbook of Toxicology of Chemical Warfare Agents*, 3rd ed.; Gupta, R.D., Ed.; Academic Press: London, UK, 2020; in press.

8. Bigley, A.N.; Raushel, F.M. The evolution of phosphotriesterase for decontamination and detoxification of organophosphorus chemical warfare agents. *Chem. Biol. Interact.* **2019**, *308*, 80-88. [CrossRef]

9. Ellman, G.L.; Courtney, K.D.; Andres, V., Jr.; Feather-Stone, R.M. A new and rapid colorimetric determination of acetylcholinesterase activity. *Biochem. Pharmacol.* **1961**, *7*, 88-95. [CrossRef]

10. Lockridge, O.; Blong, R.M.; Masson, P.; Froment, M.T.; Millard, C.B.; Broomfield, C.A. A single amino acid substitution, Gly117His, confers phosphotriesterase (organophosphorus acid anhydride hydrolase) activity on human butyrylcholinesterase. *Biochemistry* **1997**, *36*, 786-795. [CrossRef]

11. Poyot, T.; Nachon, F.; Froment, M.T.; Loiodice, M.; Wieseler, S.; Schopfer, L.M.; Lockridge, O.; Masson, P. Mutant of *Bungarus fasciatus* acetylcholinesterase with low affinity and low hydrolase activity toward organophosphorus esters. *Biochim. Biophys. Acta* **2006**, *1764*, 1470-1478. [CrossRef]
12. Mukhametgalieva, A.R.; Zueva, I.V.; Aglyamova, A.R.; Lushchekina, S.V.; Masson, P. A new sensitive spectrofluorimetric method for measurement of activity and kinetic study of cholinesterases. *Biochim Biophys. Acta Proteins Proteom.* **2020**, *1868*, 140270. [CrossRef] [PubMed]

13. Latip, W.; Knight, V.F.; Abdul Halim, N.; Ong, K.K.; Mohd Kassim, N.A.; Wan Yunus, W.M.Z.; Mohd Noor, S.A.; Mohamad Ali, M.S. Microbial Phosphotriesterase: Structure, Function, and Biotechnological Applications. *Catalysts* **2019**, *9*, 671. [CrossRef]

14. Jacquet, P.; Daude, D.; Bzdrenga, J.; Masson, P.; Elias, M.; Chabriere, E. Current and emerging strategies for organophosphate decontamination: special focus on hyperstable enzymes. *Environ. Sci. Pollut. Res. Int.* **2016**, *23*, 8200–8218. [CrossRef] [PubMed]

15. Ramalho, T.C.; de Castro, A.A.; Silva, D.R.; Cristina Silva, M.; 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. *Curr. Med. Chem.* **2016**, *23*, 1041–1061. [CrossRef]

16. Bigley, A.N.; Desormeaux, E.; Xiang, D.F.; Bae, S.Y.; Harvey, S.P.; Raushel, F.M. Overcoming the Challenges of Enzyme Evolution To Adapt Phosphotriesterase for V-Agent Decontamination. *Biochemistry* **2019**, *58*, 2039–2053. [CrossRef] [PubMed]

17. Despotovic, D.; Aharon, E.; Dubovetskyi, A.; Leader, H.; Ashani, Y.; Tawfik, D.S. A mixture of three engineered phosphotriesterases enables rapid detoxification of the entire spectrum of known threat nerve agents. *Protein Eng. Des. Sel.* **2019**. [CrossRef]

18. Cardozo, M.; de Almeida, J.; Cavalcante, S.F.A.; Salgado, J.R.S.; Goncalves, A.S.; Franca, T.C.C.; Kuca, K.; Bizzo, H.R. Biodegradation of Organophosphorus Compounds Predicted by Enzymatic Process Using Molecular Modelling and Observed in Soil Samples Through Analytical Techniques and Microbiological Analysis: A Comparison. *Molecules* **2019**, *25*, 58. [CrossRef]

19. Worek, F.; Thiermann, H.; Wille, T. Catalytic bioscavengers in nerve agent poisoning: A promising approach? *Toxicol. Lett.* **2016**, *244*, 143–148. [CrossRef]

20. Goldsmith, M.; Ashani, Y. Catalytic bioscavengers as countermeasures against organophosphate nerve agents. *Chem. Biol. Interact.* **2018**, *292*, 50–64. [CrossRef]

21. Zhang, P.; Liu, E.J.; Tsao, C.; Kasten, S.A.; Boeri, M.V.; Dao, T.L.; DeBus, S.J.; Cadieux, C.L.; Baker, C.A.; Otto, T.C.; et al. Nanoscavenger provides long-term prophylactic protection against nerve agents in rodents. *Sci. Transl. Med.* **2019**, *11*. [CrossRef]

22. Job, L.; Kohler, A.; Escher, B.; Worek, F.; Skerra, A. A catalytic bioscavenger with improved stability and reduced susceptibility to oxidation for treatment of acute poisoning with neurotoxic organophosphorus compounds. *Toxicology letters* **2020**, *321*, 138–145. [CrossRef] [PubMed]

23. Masson, P.; Nachon, F.; Broomfield, C.A.; Lenz, D.E.; Verdier, L.; Schopfer, L.M.; Lockridge, O. A collaborative endeavor to design cholinesterase-based catalytic scavengers against toxic organophosphorus esters. *Chem. Biol. Interact.* **2008**, *175*, 273–280. [CrossRef] [PubMed]

24. Lushchekina, S.V.; Schopfer, L.M.; Grigorenko, B.L.; Nemukhin, A.V.; Varfolomeev, S.D.; Lockridge, O.; Masson, P. Optimization of Cholinesterase-Based Catalytic Bioscavengers Against Organophosphorus Agents. *Front. Pharmacol.* **2018**, *9*, 211. [CrossRef]

25. Grigorenko, B.L.; Novichkova, D.A.; Lushchekina, S.V.; Zueva, I.V.; Schopfer, L.M.; Nemukhin, A.V.; Varfolomeev, S.D.; Lockridge, O.; Masson, P. Computer-designed active human butyrylcholinesterase double mutant with a new catalytic triad. *Chem. Biol. Interact.* **2019**, *306*, 138–146. [CrossRef] [PubMed]

26. Katalinic, M.; Macek Hrvat, N.; Baumann, K.; Morasi Pipercic, S.; Makaric, S.; Tomic, S.; Jovic, O.; Hrenar, T.; Milicevic, A.; Jelic, D.; et al. A comprehensive evaluation of novel oximes in creation of butyrylcholinesterase-based nerve agent bioscavengers. *Toxicol. Appl. Pharmacol.* **2016**, *310*, 195–204. [CrossRef] [PubMed]

27. Zhang, L.; Baker, S.L.; Murata, H.; Harris, N.; Ji, W.; Amitai, G.; Matyjaszewski, K.; Russell, A.J. Tuning Butyrylcholinesterase Inactivation and Reactivation by Polymer-Based Protein Engineering. *Adv. Sci. (Weinh)* **2020**, *7*, 1901904. [CrossRef]

28. Zorbaz, T.; Malinak, D.; Kuca, K.; Musilek, K.; Kovarik, Z. Butyrylcholinesterase inhibited by nerve agents is efficiently reactivated with chlorinated pyridinium oximes. *Chem. Biol. Interact.* **2019**, *307*, 16–20. [CrossRef]

29. Goldsmith, M.; Aggarwal, N.; Ashani, Y.; Jubran, H.; Greisen, P.J.; Ovchinnikov, S.; Leader, H.; Baker, D.; Sussman, J.L.; Goldenzweig, A.; et al. Overcoming an optimization plateau in the directed evolution of highly efficient nerve agent bioscavengers. *Protein Eng. Des. Sel.* **2017**, *30*, 333–345. [CrossRef]
30. Chabrière, E.; Daudé, D.; Elias, M. Nouvelles enzymes pte mutées. WO2019016468A1, 17 July 2019.

31. Chen, Y.C. Beware of docking! *Trends Pharmacol. Sci.* **2015**, *36*, 78–95. [CrossRef]

32. Xu, Y.; Cheng, S.; Sussman, J.L.; Silman, I.; Jiang, H. Computational Studies on Acetylcholinesterases. *Molecules* **2017**, *22*, 1324. [CrossRef] [PubMed]

33. Lushchekina, S.V.; Makhaeva, G.F.; Novichkova, D.A.; Zueva, I.V.; Kovaleva, N.V.; Richardson, R.J. Supercomputer Modeling of Dual-Site Acetylcholinesterase (AChE) Inhibition. *Supercomput. Front. Innov.* **2018**, *5*, [CrossRef]

34. Santos-Martins, D.; Forli, S.; Ramos, M.J.; Olson, A.J. AutoDock4(Zn): an improved AutoDock force field for small-molecule docking to zinc metalloproteins. *J. Chem. Inf. Model.* **2014**, *54*, 2371–2379. [CrossRef] [PubMed]

35. Hu, X.; Balaz, S.; Shelver, W.H. A practical approach to docking of zinc metalloproteinase inhibitors. *J. Mol. Graph. Model.* **2004**, *22*, 293–307. [CrossRef] [PubMed]

36. Friboulet, A.; Rieger, F.; Goudou, D.; Amitai, G.; Taylor, P. Interaction of an organophosphate with a peripheral site on acetylcholinesterase. *Biochemistry* **1990**, *29*, 914–920. [CrossRef]

37. Masson, P.; Legrand, P.; Bartels, C.F.; Froment, M.-T.; Schopfer, L.M.; Lockridge, O. Role of Aspartate 70 and Tryptophan 82 in binding of succinyldithiocholine to human butyrylcholinesterase. *Biochemistry* **1997**, *36*, 2266–2277. [CrossRef]

38. Kumar, A.; Darreh-Shori, T. DMSO: A Mixed-Competitive Inhibitor of Human Acetylcholinesterase. *ACS Chem. Neurosci.* **2017**, *8*, 2618–2625. [CrossRef]

39. Bigley, A.N.; Xu, C.; Henderson, T.J.; Harvey, S.P.; Raushel, F.M. Enzymatic neutralization of the chemical warfare agent VX: evolution of phosphotriesterase for phosphorothiolate hydrolysis. *J. Am. Chem. Soc.* **2013**, *135*, 10426–10432. [CrossRef]

40. Reed, B.A.; Sabourin, C.L.; Lenz, D.E. Human butyrylcholinesterase efficacy against nerve agent exposure. *J. Biochem. Mol. Toxicol.* **2017**, *31*. [CrossRef]

41. Cerasoli, D.M.; Armstrong, S.J.; Reeves, T.E.; Hodgins, S.M.; Kasten, S.A.; Lee-Stubbs, R.B.; Cadieux, C.L.; Otto, T.C.; Capacio, B.R.; Lenz, D.E. Butyrylcholinesterase, a stereospecific in vivo bioscavenger against nerve agent intoxication. *Biochem. Pharmacol.* **2020**, *171*, 113670. [CrossRef]

42. Nachon, F.; Carletti, E.; Wandhammer, M.; Nicolet, Y.; Schopfer, L.M.; Masson, P.; Lockridge, O. X-ray crystallographic snapshots of reaction intermediates in the G117H mutant of human butyrylcholinesterase, a nerve agent target engineered into a catalytic bioscavenger. *Biochem. J.* **2011**, *434*, 73–82. [CrossRef] [PubMed]

43. Kowalczewska, M.; Fenollar, F.; Villard, C.; Azza, S.; Roux, M.; Raoult, D. An immunoproteomic approach for identification of clinical biomarkers of Whipple's disease. *Proteomics Clin. Appl.* **2008**, *2*, 504–516. [CrossRef] [PubMed]

44. Yi, L.; Li, H.; Sun, L.; Liu, L.; Zhang, C.; Xi, Z. A highly sensitive fluorescence probe for fast thiol-quantification assay of glutathione reductase. *Angew. Chem. Int. Ed. Engl.* **2009**, *48*, 4034–4037. [CrossRef] [PubMed]

45. Terekhov, S.S.; Smirnov, I.V.; Stepanova, A.V.; Bobik, T.V.; Mokrushina, Y.A.; Ponomarenko, N.A.; Belogurov, A.A., Jr.; Rubtsova, M.P.; Kartseva, O.V.; Gomzikova, M.O.; et al. Microfluidic droplet platform for ultrahigh-throughput single-cell screening of biodiversity. *Proc. Natl. Acad. Sci. USA* **2017**, *114*, 2550–2555. [CrossRef] [PubMed]

46. Cornish-Bowden, A. A simple graphical method for determining the inhibition constants of mixed, uncompetitive and non-competitive inhibitors. *Biochem. J.* **1974**, *137*, 143–144. [CrossRef] [PubMed]

47. Pajk, S.; Knez, D.; Kosak, U.; Zorovic, M.; Brazzolotto, X.; Coquelle, N.; Nachon, F.; Colletier, J.P.; Zivin, M.; Stojan, J.; et al. Development of potent reversible selective inhibitors of butyrylcholinesterase as fluorescent probes. *J. Enzyme. Inhib. Med. Chem.* **2020**, *35*, 498–505. [CrossRef]

48. Bigley, A.N.; Raushel, F.M. Catalytic mechanisms for phosphotriesterases. *Biochim. Biophys. Acta* **2013**, *1834*, 443–453. [CrossRef]

49. Chen, D.; Menche, G.; Power, T.D.; Sower, L.; Peterson, J.W.; Schein, C.H. Accounting for ligand-bound metal ions in docking small molecules on adenylyl cyclase toxins. *Proteins* **2007**, *67*, 593–605. [CrossRef]

50. Nemukhin, A.V.; Grigorenko, B.L.; Lushchekina, S.V.; Varfolomeev, S.D. Quantum chemical modelling in the research of molecular mechanisms of enzymatic catalysis. *Russ. Chem. Rev.* **2012**, *81*, 1011–1025. [CrossRef]

51. Schmidt, M.W.; Baldridge, K.K.; Boatz, J.A.; Elbert, S.T.; Gordon, M.S.; Jensen, J.H.; Koseki, S.; Matsunaga, N.; Nguyen, K.A.; Su, S.J.; et al. General Atomic and Molecular Electronic-Structure System. *J. Comput. Chem.* **1993**, *14*, 1347–1363. [CrossRef]
52. Löwdin, P.-O. On the Nonorthogonality Problem. In *Adv. in Quantum Chem.*; Per-Olov, L., Ed.; Academic Press: Cambridge, MA, USA, 1970; Volume 5, pp. 185–199.

53. Benning, M.M.; Shim, H.; Raushel, F.M.; Holden, H.M. High resolution X-ray structures of different metal-substituted forms of phosphotriesterase from Pseudomonas diminuta. *Biochem. 2001*, *40*, 2712–2722. [CrossRef] [PubMed]

54. Harel, M.; Aharoni, A.; Gaidukov, L.; Brumshtein, B.; Khersonsky, O.; Meged, R.; Dvir, H.; Ravelli, R.B.; McCarthy, A.; Toker, L.; et al. Structure and evolution of the serum paraoxonase family of detoxifying and anti-atherosclerotic enzymes. *Nat. Struct. Mol. Biol.* **2004**, *11*, 412–419. [CrossRef] [PubMed]

55. Hiblot, J.; Gotthard, G.; Chabriere, E.; Elias, M. Characterisation of the organophosphate hydrolase catalytic activity of SsoPox. *Sci. Rep.* **2012**, *2*, 779. [CrossRef] [PubMed]

56. Word, J.M.; Lovell, S.C.; Richardson, J.S.; Richardson, D.C. Asparagine and glutamine: using hydrogen atom contacts in the choice of side-chain amide orientation. *J. Mol. Biol.* **1999**, *285*, 1735–1747. [CrossRef]

57. Valiev, M.; Bylaska, E.J.; Govind, N.; Kowalski, K.; Straatsma, T.P.; Van Dam, H.J.J.; Wang, D.; Nieplocha, J.; Apra, E.; Windus, T.L.; et al. NWChem: A comprehensive and scalable open-source solution for large scale molecular simulations. *Comput. Phys. Commun.* **2010**, *181*, 1477–1489. [CrossRef]

58. Morris, G.M.; Huey, R.; Lindstrom, W.; Sanner, M.F.; Belew, R.K.; Goodsell, D.S.; Olson, A.J. AutoDock4 and AutoDockTools4: Automated docking with selective receptor flexibility. *J. Comput. Chem.* **2009**, *30*, 2785–2791. [CrossRef]

59. Morris, G.M.; Goodsell, D.S.; Halliday, R.S.; Huey, R.; Hart, W.E.; Belew, R.K.; Olson, A.J. Automated docking using a Lamarckian genetic algorithm and an empirical binding free energy function. *J. Comput. Chem.* **1998**, *19*, 1639–1662. [CrossRef]

60. Voevodin, V.; Antonov, A.; Nikitenko, D.; Shvets, P.; Sobolev, S.; Sidorov, I.; Stefanov, K.; Voevodin, V.; Zhumatiy, S. Supercomputer Lomonosov-2: Large Scale, Deep Monitoring and Fine Analytics for the User Community. *Supercomput. Front. Innov.* **2019**, *6*, 4–11. [CrossRef]

<figure><img src="image_14.png" /></figure>

© 2020 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).