Design, Synthesis and Biological Evaluation of 1,4-Benzenesulfonamide Derivatives as Glyoxalase I Inhibitors

Background Glyoxalase system is one of the defense cellular mechanisms that protect cells against endogenous harmful metabolites, mainly methylglyoxal (MG), through conversion of cytotoxic methylglyoxal into the non-toxic lactic acid. Glyoxalase system comprises of two enzymes glyoxalase I, glyoxalase II, and a catalytic amount of reduced glutathione. Cancerous cells overexpress glyoxalase I, making it a target for cancer therapy. Many studies have been conducted to identify potent Glx-I inhibitors. Methods Aiming to discover and develop novel Glx-I inhibitors, a series of 1,4-benzenesulfonamide derivatives were designed, synthesized, and biologically evaluated in vitro against human Glx-I enzyme. Seventeen compounds were designed based on the hit compound that was obtained from searching the National Cancer Institute (NCI) database. The synthesis of the target compounds (13–29) was accomplished utilizing an azo coupling reaction of aniline derivatives and activated substituted aromatic compounds. To understand the binding mode of the active compounds at the active site of Glx-I, docking studies were performed. Results Structure activity relationship (SAR) studies were accomplished which led to the identification of several compounds that showed potent inhibitory activity with IC50 values below 10 μM. Among the compounds tested, compounds (E)-2-hydroxy-5-((4-sulfamoylphenyl)diazenyl)benzoic acid (26) and (E)-4-((8-hydroxyquinolin-5-yl)diazenyl) benzenesulfonamide (28) displayed potent Glx-I inhibitory activity with IC50 values of 0.39 μM and 1.36 µM, respectively. Docking studies of compounds 26 and 28 were carried out to illustrate the binding mode of the molecules into the Glx-I active site. Conclusion Our results show that compounds 26 and 28 displayed potent Glx-I inhibitory activity and can bind the Glx-I well. These findings should lead us to discover new classes of compounds with better Glx-I inhibition.


Introduction
Advanced glycation end products (AGEs) are group of molecules produced from non-enzymatic glycation and oxidation of proteins and lipids. [1][2][3] It has been established that AGEs are involved in many diseases including aging, neurological disorders, diabetic complications and cancer. 4,5 α-Oxoaldehydes such as methylglyoxal (MG) are reactive dicarbonyl compounds contribute to the production of AGEs. 6,7 MG is produced endogenously through different enzymatic and nonenzymatic pathways including metabolic glycolysis processes and fragmentation of triosephosphates glyceraldehyde-3phosphate and dihydroxyacetone phosphate. [7][8][9][10][11][12] Glyoxalase system is one of the cellular mechanisms that protects cells against AGEs production through conversion of MG into the non-toxic lactic acid 13,14 and ultimately prevents the accumulation of reactive dicarbonyl compounds in cells (Figure 1).
Glyoxalase system comprises glyoxalase I (Glx-I), glyoxalase II (Glx-II), and a catalytic amount of reduced glutathione. 15 Under normal cellular conditions, thiohemiacetal 2 form from non-enzymatic condensation of glutathione (GSH) and MG 1. In the multistep pathway of glyoxalase system, Glx-I converts thiohemiacetal 2 into the S-Dlactoylglutathione 4 which is then hydrolyzed by Glx-II to produce D-lactate and glutathione. [16][17][18] In cancerous cells, where inducing apoptosis to the cells is the ultimate goal, Glx-I enzyme is abnormally overexpressed. [19][20][21][22][23][24] Inhibition of Glx-I enzyme should result in accumulation of MG 25 at cytotoxic levels eventually leading to apoptosis. 26 Therefore, inhibitors of Glx-I enzyme are considered effective anti-tumor agents.
Many studies have been conducted to identify potent Glx-I inhibitors. 27,28 Different drug discovery methodologies have been employed to identify such inhibitors including glutathione derivatives, [29][30][31][32] natural products and natural products derivatives, [33][34][35][36] high-throughput screening, 37,38 and in silico computer aided drug design. 35,[39][40][41][42][43][44][45] Recently, the active site of Glx-I was modeled using a structure-based pharmacophore approach that was utilized to construct a pharmacophore model, a tool that was used to search the NCI database to find new Glx-I inhibitors. Several promising hits were identified based on the search query, which were purchased and tested for their inhibitory potential, and showed good inhibitory activity. 46 This study selected compound 6 as a lead compound toward the development of more potent Glx-I inhibitors. The selected compound has a benzenesulfonamide moiety and exhibited an IC 50 value of 3.65 μM, as shown in Figure 2.
Recently, our research group has designed and synthesized 1,3-diazenylbenzenesulfonamide derivatives of the lead compound, that led to identification of promising Glx-I inhibitors. 47 The present study, as a continuation of the previous research, aims to prepared 1,4-benzenesulfonamide derivatives as structural analogues of compound 6. The compounds were synthesized and biologically evaluated in vitro against human Glx-I enzyme. Structure-activity relationship (SAR) studies were performed after structural modifications of the synthesized compounds.

Materials and Methods Chemistry
All organic reactions were performed with standard laboratory instruments with magnetic stirring. All chemicals or reagents and solvents used for syntheses were commercially available and were used as received without further purification. Analytical thin layer chromatography (TLC) was used to monitor all organic reactions in which the disappearance of the starting material and the formation of the product(s) indicated the completion of the reaction. TLC was performed on Merk silica gel 60 F254 aluminum backed plates. Merck silica gel 60 (0.2 mm) was used for column chromatography. Visualization of TLC was accomplished with UV light (254 nm). A Bruker AVANCE 400 MHz instrument was used to record 1 H/ 13 C NMR spectra at ambient temperature. The spectra were measured in ppm and calibrated with the residual proton peak of the solvent, DMSO-d 6 . The 1 H NMR data are presented as follows: chemical shift (δ ppm), multiplicity (s = singlet, d = doublet, t = triplet, q = quartet, m = multiplet), coupling constants, J (Hz), and integration. The 13 C NMR analyses were reported in terms of chemical shift. 1 H and 13 C NMR spectra for all compounds are provided in the Supplementary Material (S1-S34). High resolution mass spectrometry (HRMS) was performed on an Impact II by Brucker using UHR QTOF MS (Ultra-High-Resolution Q-Time-Of-Flight) with full scan MS and broadband CID (bbCID) MS/MS methods located at the NABA HIKMA industrial and testing services, Fuhais, Jordan.
Step 2: General procedure for preparation of compounds 13-29: Procedure 1: To a solution of a R 2 compound (1.0 equiv.) in a 2M sodium hydroxide (NaOH, 37.5 mL) at −5°C was added dropwise the solution prepared in step 1. Immediately after mixing, the solution changed color, and a precipitate was formed. Stirring continued at the same temperature until the reaction was completed. The pH was then adjusted to 7.0 after the mixture was warmed to room temperature. After filtration, the collected solid was washed with deionized water, air dried, then purified utilizing flash column chromatography (5% MeOH/DCM) or by washing with hot ethanol to yield the final products.
Procedure 2: Similar procedure as in procedure 1 except in procedure 2 a saturated aqueous solution of NaOAc (7.5 mL) was used as a base instead of 2M NaOH.

Biological Evaluation
The Glx-I inhibitory activity was conducted in vitro using established published procedure. 35,42,43,47 The IC 50 values of all compounds were calculated using GraphPad Prism 6 (2012). Myricetin was used as positive control with an IC 50 equals to 3.38 ± 0.41 (µM).

Molecular Modeling Study
Using discovery studio software from Biovia ® , the most active compounds were docked into the active site of Glx-I. 50 The enzyme was prepared according to the default parameters using "Prepare protein" protocol and minimized while the necessary modifications of the zinc atom were conducted to agree with the previous work set by our research group. The parameters for the minimize protein was left as default and the zinc atom at the active site was assigned to be +2 charge with octahedral geometry. The docking protocols were the LibDock and Flexible Docking, which are validated build-in software supported by Discovery studio ® . The LibDock parameters were set as default with the assignment of 10 A°as a docking sphere to enable the docking to reach all the active site three major areas (positively ionized entrance of the active site, zinc atom, deep hydrophobic pocket). Flexible docking was performed by allowing the amino acids within the active site to be flexible during docking. Parameters were left as default with 10 A°docking sphere.

Structure-Activity Relationship (SAR) Studies
In this study compound 6 (IC 50 = 3.65 µM) was used as a lead compound to identify more potent Glx-I inhibitors. From the perspective of compound structure, compound 6 possesses 5-diazenyl-2-methylbenzenesulfonamide moiety attached to position 8 of the 5-amino-2-naphthalenol moiety, in which the diazenyl moiety is located at the meta position relative to the benzenesulfonamide moiety.
A series of structural analogues of compound 6 were designed to locate diazenyl moiety at the para position relative to the benzenesulfonamide moiety. The analogues were synthesized and biologically evaluated in vitro against human Glx-I enzyme to discover potential Glx-I inhibitors. Structure activity relationship (SAR) studies were accomplished which led to identification of a series of analogues (13)(14)(15)(16)(17)(18)(19)(20)(21)(22)(23)(24)(25)(26)(27)(28)(29) that showed a broad range of inhibition (potent to moderate) as shown in Table 1. SAR studies started with the synthesis of (E)-4-((4-amino-7-hydroxynaphthalen-1-yl)diazenyl) benzenesulfonamide (13) that possesses 4-diazenylbenzenesulfonamide moiety attached to the position 8 of the 5-amino-2-naphthalenol moiety. Removing the methyl group from compound 6 and shifting the diazenyl moiety to the para position relative to the benzenesulfonamide moiety resulted in reduced inhibitory activity with IC 50 39. 06 μM. Compounds 14 and 15 were synthesized to study the effect of naphthalene's substituents on the activity of compound 13. Moving the hydroxyl group in compound 13 from position 2 to position 3 led to the formation of compounds 14 with enhanced inhibitory activity IC 50 22.98 μM. Whereas removing the hydroxyl group from the naphthalene moiety (compound 15) significantly increased the inhibitory activity producing a potent compound with IC 50 1.85 μM. These results suggest that enhancement of the activity possibly due to the naphthalene moiety forming much stronger interactions with the surrounding hydrophobic amino acids in the active site of the enzyme. In addition, the presence of the hydroxyl group at the naphthalene ring may cause unfavorable interactions in the hydrophobic pocket at the active site of the enzyme. Also, the removal of the OH reduced the desolvation penalty contributed to presence of polar functional groups in the structure.
On the other hand, attaching the 4-diazenylbenzenesulfonamide moiety at position 1 of the 5-amino-2-naphthalenol moiety and moving the amino group from position 5 to position 8 produced compound 16 with good inhibitory activity (IC 50 8.61 μM). Removing the amino group from the naphthalene moiety of compound 16 led to the formation of compound 17 with increasing inhibitory activity (IC 50 5.30 μM). These results showed that the amino group at position 8 of the naphthalene ring reduced the inhibitory activity of compound 16, which may be due to the unfavorable interactions in the hydrophobic pocket at the active site of the enzyme.
Next, different analogues of compound 13 were synthesized in which the 5-amino-2-naphthalenol moiety was replaced with groups that possess the ability to interact with the zinc ion and the active site residues of the enzyme. First, replacement of the 5-amino-2-naphthalenol moiety with 7-hydroxyl-2H-chromen-2-one moiety led to moderately active compound (compound 18) with IC 50 33.11 μM. Second, the replacement of 5-amino-2-naphthalenol moiety with 2-hydroxy-1-naphthaldehyde moiety led to the formation of compound 19 with similar inhibitory activities (IC 50 39.00 μM), when compared to compound 13.
Furthermore, replacement of the formyl group in compound 19 with a carboxy group led to the formation of compound 20 with no noticeable improvement of the inhibitory activity. On the other hand, compounds 21-22 were synthesized to explore the SAR for compound 20. Isosteric replacement of the sulfonamide moiety with a carboxy group produced compound 21 (IC 50 32.09 μM) with slight increase in activity, whereas, removing the sulfonamide moiety from compound 20 resulted in the formation of compound 22 with good inhibitory activity (IC 50 12.46 μM). These results suggest that the presence of polar or ionizable groups at the para position of benzene ring reduces the inhibitory activity of the compounds, which may be attributed to an increase in the desolvation penalty. These groups have the potential to clasp water molecules tightly via a series of hydrogen bonds or via ion dipole interactions in the case of carboxylic acid. Most probably, the enthalpy gain from zinc chelation is not enough to compensate the desolvation penalty. Switching the position of the carboxyl and hydroxyl groups, led to compound 23 with significant enhancement in the inhibitory activity    Next, we synthesized compound 28 that contain 8-hydroxyquinoline moiety in the place of 5-amino-2-naphthalenol which provided a potent derivative with IC 50 1.36 μM. Compounds 27 and 29 were synthesized to explore the SAR for compound 28. Investigating the importance of the pyridine moiety in compound 28 was accomplished by synthesizing compound 27 that lack the pyridine moiety. It was found to have less inhibitory activity when compared to compound 28. This result indicates that the pyridine moiety has a role in binding interactions within the active site of the enzyme. Compound 29 was synthesized to explore the importance of the sulfonamide moiety to the activity of compound 28. The result showed that removing the sulfonamide moiety causes significant decrease in the inhibitory activity of the compound (IC 50 = 8.81 μM), possibly due to the loss of some binding interaction at the active site.

Binding Mode of 26 and 28 into the Active Site of Glx-I Enzyme
Molecular docking was performed using flexible docking and Libdock protocols. The poses of the most active compounds (26 and 28) were analyzed and studied to investigate their binding pattern within the active site of the enzyme and to explore the most important ligand-receptor interactions such as the role of zinc atom. Binding mode of the most active compound 26 into the catalytic site of Glx-I is shown in (Figure 3), which explains the good activity of the compound. The binding model showed that the ionizable carboxyl group chelates zinc atom cofactor in a bidentate pattern at the active site of the enzyme. Additionally, the sulfonamide group provided significant hydrogen bond interactions with Arg37, Asn103, and Arg122 that enhanced compound's potency. Furthermore, the potency of compound 26 is also due to hydrophobic interactions between the aromatic benzene moieties of compound 26 and the hydrophobic amino acid residues (Phe162, Phe67, Phe160, and Met157) at the active site of the enzyme. On the other hand, binding model of compound 28 was generated by its docking into the active site of Glx-I. The probable binding mode is illustrated in (Figure 4). The binding model showed that the hydroxyl group chelates the zinc atom at the active site of the enzyme. Moreover, the sulfonamide group forms hydrogen bond interactions with Arg122, Arg37 and Asn103. Additionally, hydrophobic interactions occurred between the quinoline moiety and the hydrophobic amino acid residues (Phe162, Met157, Leu69, and Leu160) at the active site, whereas the benzene moiety forms hydrophobic interactions with the hydrophobic amino acid residue Phe67at the active site of the enzyme.

Conclusion
To sum up, the previously reported Glx-I inhibitor (Compound 6, IC 50 3.65 μM) was used as a lead compound to identify more potent Glx-I inhibitors. To achieve our goals, a series of structural analogues were synthesized and biologically evaluated in vitro against human Glx-I enzyme. Several analogues were found to exhibit potent activity with IC 50 below 10 μM. Structural optimization of the lead compound (compound 6) produced potent compounds, compound 26 with IC 50 0.39 μM and compound 28 with IC 50 1.36 μM. SAR studies revealed that the carboxy moieties in compound 26 plays pivotal role in the inhibitory activity of the compound. Additionally, SAR evaluation of compound 28 signifies the importance of the sulfonamide and pyridine moieties for the inhibitory activity of the compound. Furthermore, molecular docking was performed to give the probable binding modes of compounds 26 and 28 into the catalytic site of Glx-I enzyme. These findings will open the doors for better discoveries of potent Glx I inhibitors.