Abstract

Sulfa drugs are well-known antibacterial agents containing N-substituted sulfonamide group on para position of aniline ring (NH2RSO2NHR′). In this study 2,4-dichloro-1,3,5-triazine derivatives of sulfa drugs, sulfamerazine (1b), sulfaquinoxaline (2b), sulfadiazine (3b), sulfadimidine (4b), and sulfachloropyrazine (5b) (1a5a) were synthesized and characterized. Their carbonic anhydrase inhibition activity was evaluated against bovine cytosolic carbonic anhydrase isozyme II (bCA II). For the sake of comparison the CA inhibition activity of the parent sulfa drugs (1b5b) was also evaluated. A significant increase in CA inhibition activity of sulfa drugs was observed upon substitution with 2,4-dichloro-1,3,5-triazine moiety. Molecular docking studies were carried out to highlight binding site interactions. ADME properties were calculated to evaluate drug likeness of the compounds.

1. Introduction

Carbonic anhydrase (CA, EC 4.2.1.1) is a well-known ubiquitous zinc containing metalloenzyme that is found in animals including man, plants, bacteria, and archaea. The active site zinc ion (Zn2+) has been conserved in all classes of carbonic anhydrases [1]. This enzyme catalyzes an apparently simple yet physiologically important reaction of interconversion of water and carbon dioxide into bicarbonate and protons. Therefore CA has an important role to play in the transport of CO2 from metabolizing tissues to lungs. It is also responsible for maintaining acid/base and electrolyte balance in blood [2, 3]. Certain biosynthetic reactions are also assisted by CA such as lipogenesis, at the level of pyruvate carboxylation [4, 5], ureagenesis [6], and gluconeogenesis [6, 7]. In mammals carbonic anhydrase has sixteen different isozymes based on their distribution in tissues and subcellular localization. The cytosolic isozymes are CA I, CA II, CA III, CA VII, and CA XIII, whereas CA IV, CA IX, CA XII, CA XIV, and CA XV are membrane bound isozymes, and CA VA and CA VB are mitochondrial isozymes; CA VI is secreted isozyme mainly present in the saliva [8]. There are certain physiological disorders that are characterized by overexpression of CA [911], hence CAs have emerged as valuable drug targets. Many clinically established drugs are CA inhibitors and are used to treat disorders such as glaucoma, acidic ulcers, mountain/sea sickness, and epilepsy [12]. Carbonic anhydrase is also an important drug target for treating obesity and many sulfonamide inhibitors have proved to be efficient antiobesity agents [1316]. The transmembrane isozymes CA IX and CA XII have been found to be overexpressed in hypoxic tumors (having acidic environment) whereas their distribution in normal cells remains low [1722].

Sulfonamides and their derivatives are well-known inhibitors of carbonic anhydrase [18, 23]. Sulfa drugs are derived from sulfonamides; however all sulfonamides are not sulfa drugs, and the term sulfa drug is only used for clinically used antibacterial agents that are structurally derived from 4-aminobenzenesulfonamide, where the sulfonamide nitrogen is substituted (NH2RSO2NHR′) [24, 25]. The carbonic anhydrase inhibition activity of sulfa drugs has not been explored. Previously 1,3,5-triazine [26, 27] and 1,2,4-triazine [28] derivatives of different sulfonamides have been reported as efficient inhibitors of CA. Herein we report the synthesis of new 2,4,6-trichloro-1,3,5-triazine (TCT) derivatives of sulfa drugs (1a5a) and their carbonic anhydrase inhibition activity against bovine cytosolic carbonic anhydrase II (bCA II). For the purpose of comparison, the carbonic anhydrase inhibition activity of parent sulfa drugs (1b5b) is also reported.

2. Material and Methods

All chemicals used were purchased from either Sigma or Aldrich and used as such without further purification. Commercially available solvents were used. Ethanol was distilled and dried using standard methods and stored over molecular sieves. Reaction progress and product purity were checked via precoated TLC plates (silica gel, 0.2 mm, 60 HF254, Merck). TLC spots were visualized under short and long wavelength UV light. Bovine cytosolic carbonic anhydrase II (bCA II) was used. Melting points were taken on a Gallenkamp melting point apparatus and were uncorrected. FTIR spectra were taken on Perkin Elmer Spectrum BX-II. LECO CHNS 630 series elemental analyzer (model 630-200-200) was used for elemental analysis. For 1H and 13C-NMR analysis Bruker Avance DRX500 spectrometer was used with TMS as an internal standard and DMSO-d6 as solvent.

2.1. General Method of Synthesis

For the synthesis of TCT derived sulfa drugs, 2,4,6-trichloro-1,3,5-triazine (TCT, 0.01 mol), respective sulfa drug (0.01 mol), and sodium carbonate (0.1–0.2 g) were taken in a round bottom flask, 20 mL of ethanol and 5–7 mL acetone were added to it, and stirring was continued until a clear solution resulted. The reaction mixture was refluxed with constant stirring. After 2 hours solid precipitate began to appear in the reaction mixture; the reaction was allowed to continue for another hour after which the solid product was filtered, washed, and dried. Compounds were recrystallized with a mixture of acetone and acetonitrile.

2.1.1. Synthesis of 4-[(4,6-Dichloro-1,3,5-triazin-2-yl)amino]-N-(4-methylpyrimidin-2-yl)benzenesulfonamide (1a)

Yield 80%; m.p. 235–237°C; IR (ν, cm−1): 3438 (NH2), 1145 (), 1341 (), 779 (C–Cl); Anal. calcd for C14H11Cl2N7O2S: C(40.79%), H(2.69%), N(23.78%), S(7.78%); found C(39.89%), H(3.4%), N(23.18%), S(7.54%). 1H-NMR (500 MHz, DMSO-d6), δ (ppm): 2.29 (s, 3H, –CH3), 8.31 (1H, m, H4′′), 7.99 (2H, d,  Hz, H3′, H5′), 7.65 (2H, d,  Hz, H2′, H6′), 7.78 (1H, d,  Hz, H5′′), 11.13 (s, 1H, SO2NH), 8.34 (s, 1H, NH). 13C-NMR (125 MHz, DMSO-d6), δ (ppm): 157.70 (C1, C2), 156.99 (C3), 156.65 (C1′′), 150.02 (C1′), 130.08 (C3′′, C5′′), 129.11 (C3′, C5′), 120.74 (C4′), 119.59 (C2′, C6′), 114.85 (C4′′), 23.45 (CH3).

2.1.2. Synthesis of 4-[(4,6-Dichloro-1,3,5-triazin-2-yl)amino]-N-(quinoxalin-2-yl)benzenesulfonamide (2a)

Yield 87%; m.p. 268–270°C; IR (ν, cm−1): 3248 (NH2), 1151 (), 1314 (), 727 (C–Cl); Anal. calcd for C17H11Cl2N7O2S: C (45.55%) H(2.47%), N (21.87%), S(7.15%); found C(44.87%), H(2.71%), N(20.90%), S(7.54%). 1H-NMR (500 MHz, DMSO-d6), δ (ppm): 8.08 (2H, d,  Hz, H3′, H5′), 7.93 (2H, d,  Hz, H2′, H6′), 7.79 (1H, s, H9′′), 7.78 (2H, m, H4′′, H7′′), 7.72 (2H, m, H5′′, H6′′), 10.82 (s, 1H, SO2NH), 8.63 (s, 1H, NH). 13C-NMR (125 MHz, DMSO-d6), δ (ppm): 164.08 (C3), 154.14 (C1, C2), 150.02 (C1′′), 146.19 (C1′), 142.16 (C3′′), 131.09 (C8′′), 129.12 (C3′, C5′), 128.87 (C4′′, C7′′), 127.41 (C5′′, C6′′), 120.91 (C4′), 120.06 (C2′, C6′).

2.1.3. Synthesis of 4-[(4,6-Dichloro-1,3,5-triazin-2-yl)amino]-N-(pyrimidin-2-yl)benzenesulfonamide (3a)

Yield 88%; m.p. 221–223°C; IR (ν, cm−1): 3286 (NH2), 1143 (), 1335 (), 688 (C–Cl); Anal. calcd for C13H9Cl2N7O2S: C(39.21%), H(2.28%), N(24.62%), S(8.05%); found C(39.19%), H(2.46%), N(24.17%), S(8.51%). 1H-NMR (500 MHz, DMSO-d6), δ (ppm): 8.50 (2H, d,  Hz, H3′, H5′), 7.96 (2H, d,  Hz, H2′, H6′), 7.77 (2H, d,  Hz, H3′′, H5′′), 7.04 (1H, t,  Hz, H4′′), 8.0 (1H, s, NH), 10.82 (1H, s, SO2NH). 13C-NMR (125 MHz, DMSO-d6), δ (ppm): 158.48 (C1, C2), 157.03 (C3), 154.21 (C1′′), 141.81 (C1′), 134.89 (C3′′, C5′′), 128.92 (C3′, C5′), 120.81 (C4′), 120.07 (C2′, C6′), 115.91 (C4′′).

2.1.4. Synthesis of 4-[(4,6-Dichloro-1,3,5-triazin-2-yl)amino]-N-(4,6-dimethylpyrimidin-2-yl)benzenesulfonamide (4a)

Yield 84%; m.p. 327–329°C; IR (ν, cm−1): 3211 (NH2), 1161 (), 1334 (), 697 (C–Cl); Anal. calcd for C15H13Cl2N7O2S: C(42.26%), H(3.07%), N(23.00%), S(7.52%); found C(42.87%), H(2.91%), N(23.20%), S(7.56). 1H-NMR (500 MHz, DMSO-d6), δ (ppm): 2.24 (6H, s, CH3), 6.74 (1H, s, H4′′), 7.95 (2H, d,  Hz, H3′, H5′), 7.85 (2H, d,  Hz, H2′, H6′), 7.83 (1H, s, NH), 10.61 (1H, s, SO2NH). 13C-NMR (125 MHz, DMSO-d6), δ (ppm): 23.02 (CH3), 164.0 (C1, C2), 156.33 (C3), 142.20 (C1′), 129.21 (C3′′, C5′′), 129.09 (C3′, C5′), 120.02 (C4′), 119.71 (C4′′).

2.1.5. Synthesis of N-(6-Chloropyrazin-2-yl)-4-[(4,6-dichloro-1,3,5-triazin-2-yl)amino]benzenesulfonamide (5a)

Yield 85%; m.p. 236–238°C; IR (ν, cm−1): 3255 (NH2), 1158 (), 1336 (), 692 (C–Cl); Anal. calcd for C13H8Cl3N7O2S: C(36.09%), H(1.86%), N(22.66%), S(7.41%); found C(36.15%), H(1.94%), N(22.03%), S(7.44%). 1H-NMR (500 MHz, DMSO-d6), δ (ppm): 8.01 (2H, d,  Hz, H3′, H5′), 7.62 (2H, d,  Hz, H2′, H6′), 7.58 (1H, s, H4′′), 7.32 (1H, s, H6′′), 10.71 (s, 1H, SO2NH), 8.54 (s, 1H, NH). 13C-NMR (125 MHz, DMSO-d6), δ (ppm): 164.08 (C3), 154.14 (C1, C2), 150.02 (C1′′), 146.19 (C1′), 152.10 (C3′′), 129.12 (C3′, C5′), 132.87 (C4′′), 128.41 (C6′′), 121.21 (C4′), 120.13 (C2′, C6′).

2.2. In Vitro Carbonic Anhydrase Inhibition Assay

Carbonic anhydrase inhibition was measured by the reported method [29] after standardization of reaction conditions such as concentration of enzyme and substrate, buffer pH, and duration of reaction. The method is based on spectrophotometric determination p-nitrophenol, the CA catalyzed hydrolysis product of p-nitrophenyl acetate. The reaction mixture consisted of 60 μL of 50 mM Tris-sulfate buffer (pH 7.6, containing 0.1 mM ZnCl2), 10 μL (0.5 mM) test compound in 1% DMSO, and 10 μL (50 U) bovine enzyme per well. Contents were mixed together and preincubated at 25°C for 10 minutes. Plates were preread at 348 nm using a 96-well plate reader. Substrate, p-nitrophenyl acetate (6 mM stock using <5% acetonitrile in buffer), 20 μL was freshly prepared and added per well to achieve 0.6 mM concentration per well. Total reaction volume was made to 100 μL. After incubating for 30 minute at 25°C, the contents were mixed and read at 348 nm. Suitable controls with DMSO and standard inhibitor acetazolamide (AZM) were included in the assay. Results reported are mean of three independent experiments (±SEM) and expressed as percent inhibitions calculated by the following formula: Inhibition (%) = [100 − (abs of test comp/abs of control) × 100)]. IC50 values of compounds (0.5 mM concentration) exhibiting >50% inhibition activity were calculated after suitable dilutions.

2.3. Molecular Docking Studies
2.3.1. Receptor and Ligand Preparation

High resolution crystal structures of hCA II (0.90 Å) [30] and bCA II (1.95 Å) [31] were downloaded from the Protein Data Bank (PDB-IDs: 3K34 and 1V9E, resp.). Structural comparison and percent similarity (87.3%) and percent identity (80%) between bovine and human CA II have already been reported by our research group [32]. The RMSD with respect to the active site residues of the two proteins is 0.22 Å [32]. Method validation was done by redocking the bound ligand extracted from the active site of hCA-II (3K34). The docking method was able to reproduce the experimentally bound conformation of ligand in the active site with an RMSD of <2 Å. Method validation via self-docking could not be carried out for bCA II, since this enzyme did not contain any cocrystallized ligand. On this basis hCA-II was selected to carry out further molecular docking studies. For docking the receptor (hCA II) was prepared using DockPrep utility of Chimera software [33], which includes standard preparation steps such as adding hydrogen atoms and adding gasteiger charges using ANTECHAMBER [34] and repairing incomplete side chains (if any) using Dunbrack rotamer library. A charge of +2 was added on the zinc atom. Before docking the structures of all molecules were drawn using ACD/ChemSketch [35]. Gasteiger charges were added on each ligand using ANTECHAMBER [34] and the energy of each molecule was minimized through 100 steepest descents and 100 conjugate gradient steps using a step size of 0.02 each using Chimera [33]. Before docking, compounds were protonated (as in aqueous environment) using LeadIT 2.1.7 [36].

2.3.2. Docking Methodology

For docking studies the binding site was defined, using LeadIT [36] with the aid of reference ligand which had cocrystallized with the enzyme. The binding site included amino acid residues within 8.0 Å of the active site. Docking studies were carried out using LeadIT [36] program that incorporates FlexX for docking calculations. The default docking and scoring parameters were kept for the docking calculations; first top 10 docked conformations of lowest energy were retained for detailed analysis of binding site interactions and binding modes. In order to validate the docking protocol, method validation was carried out by redocking the ligand extracted from 3K34; the docking was successfully able to reproduce the correct binding mode for the native ligand with an RMSD of <2.

3. Results and Discussion

3.1. Chemistry

2,4,6-Trichloro-1,3,5-triazine is an important and versatile molecule that provides the possibility of nucleophilic replacement at all three chlorine atoms. The synthetic utility and biological activities of this molecule are well known [37]. The ease of nucleophilic chloride replacement decreases as the number of chlorine atoms to be replaced decreases; thus nucleophilic replacement that results in synthesis of (a) proceeds far more swiftly than subsequent reactions that result in synthesis of molecules of types (b) and (c) (Figure 1).

Thus by reacting equimolar amounts of 2,4,6-trichloro-1,3,5-triazine and sulfa drugs, sulfamerazine (1b), sulfaquinoxaline (2b), sulfadiazine (3b), sulfadimidine (4b), and sulfachloropyrazine (5b), products of type (a) (see Scheme 1) were easily obtained without interference from di- or tri-substituted products. The structures of sulfa drugs used 1b5b are given in Figure 2.

3.2. bCA II Inhibition Studies and SAR

The synthesized compounds (1a5a) and their parent sulfa drugs (1b5b) were screened against bovine cytosolic carbonic anhydrase II (bCA II). The CA active site residues involved in direct zinc binding are His94, His96, and His119. In addition the substrate (CO2) association pocket, Thr199 loop, and the histidine proton shuttle mechanism (for regeneration of active site Zn bound hydroxide ion) are essential for CA catalytic activity and are highly conserved in all alpha class CAs. For this reason initial screening against bCA-II was carried out. For comparison the bCA II inhibition activity of clinically used standard CA inhibitor acetazolamide (AZM) was also carried out (Table 1). The CA inhibition activity for all compounds is given in Table 1. All chloro triazine derived compounds were active inhibitors of bCA II. A significant increase in CA inhibition activity was observed for sulfa drugs containing chloro triazine moieties. Sulfa drugs (1b5b) also exhibited the potential to inhibit bCA II; however their inhibitor potency is not significant. The most active sulfa drug was 2b showing CA inhibition up to 51.8% followed by 5b (42.4%). The corresponding chloro triazine derivative of 5b, that is, compound 5a, was the most active inhibitor of bCA II having IC50 value of 1.49 μM. For compounds 1a5a, the bCA II inhibition activity decreases in the order 5a > 1a > 4a > 3a > 2a, whereas for the parent sulfa drugs 1b5b, the inhibition activity decreases as 2b > 5b > 4b > 3b > 1b.

Molecular docking studies were carried out to investigate the mode of binding as well as important binding site interactions that could possibly explain the increased CA inhibition observed for chloro triazine derived compounds. Molecular docking studies were performed using LeadIT software. Before starting the molecular docking studies, the docking methodology was validated by redocking the ligand (that had cocrystallized with the hCA II enzyme 3K34). The docking method successfully reproduced the binding mode of native ligand with an RMSD of <2 Å. Although none of the compounds tested here contained a free sulfonamide group (SO2NH2), which is the typical zinc binding function that binds directly to the Zn2+ metal center of CA, still all compounds were able to inhibit this enzyme. Molecular docking studies reveal that, in addition to forming stable hydrogen bond interactions with surrounding amino acid residues as well as certain favorable hydrophobic interactions, these molecules were able to inhibit the enzyme by fitting inside the entrance of the active site such that the aryl group substituted on the sulfonamide nitrogen fits inside the hydrophobic pocket (the typical arylsulfonamide binding pocket) and the dichloro-triazine moiety fills the adjacent hydrophilic pocket (Figures 3 and 4). In all compounds studied, the sulfonamide group did not form a direct bond with the Zn2+ metallic centre, as expected, since the sulfonamide nitrogen atom is not free and is substituted. The binding modes of compounds were similar to one another. In general, for compounds 1a5a, the sulfonamide nitrogen makes hydrogen bond interaction with Gln92 and sulfonamide oxygen makes hydrogen bond interactions with Asn62 and Asn67 (Figure 5). Similar interactions were observed for sulfonamide group of sulfa drugs 1b5b (Figure 6). In most binding conformations for 1a5b, the nitrogen atom of the triazine ring was found to be involved in hydrogen bond contact with Trp5. Common residues involved in hydrophobic and van der Waals interactions are Leu198, Val121, His94, Thr200, and Phe131. As indicated by molecular docking studies compounds 1a5a are better able to inhibit CA II as compared to compounds 1b5b, due to the presence of additional 2,4-dichloro-1,3,5-triazine moiety that fits in the adjacent hydrophilic pocket as seen in Figures 3 and 4.

Docking scores and their corresponding ranks are given in Table 2. The ranks of docked ligand conformations with most favorable binding mode were decided after running Hyde affinity assessment of first 10 lowest energy docked conformations. Hyde Affinity Assessment is a utility included in the LeadIT software; it allows selection of best docked conformation on the basis of calculated binding free energy. The docked conformation which has the least binding free energy (ΔG, KJ/mol) is most favorable. Thus, out of 10 lowest energy conformations, the conformation that had the lowest binding free energy ΔG, after Hyde assessment, was chosen and its corresponding rank and docking score are given in Table 2. The 2D interaction diagrams of binding site interactions are given in Figures 5 and 6.

In drug discovery process, the evaluation of ADME properties (Absorption, Distribution, Metabolism, and Excretion) is of considerable importance, since the likelihood of success of a molecule as a drug greatly depends on its favorable ADME. Accordingly the parameters selected for this study include octanol-water distribution coefficients ( and ), the pH dependent octanol-water distribution coefficient (), number of hydrogen bond donors (HBDH), hydrogen bond acceptor (sum of nitrogen and oxygen atoms MNO), and topological polar surface area (TPSA). The evaluation of TPSA is used as a model to assess the ability of compounds to cross the blood brain barrier (BBB). Molecules with TPSA < 60 Å2 are usually expected to be completely absorbed whereas molecules with TPSA > 140 Å2 are generally unwanted and cannot be expected to exhibit sufficient bioavailability [38]. Typically, compounds having molecular weight less than 500, number of hydrogen bond acceptor (HBDH) less than 10, number of hydrogen bond donor (MNO) less than 5, and a value of less than 5 are considered to be orally bioavailable with favorable ADME profile [39, 40]. The calculated ADME properties of compounds 1a5a and 1b5b are given in Table 3. All synthesized compounds showed favorable ADME properties. No violation of Lipinski’s rule of 5 was observed.

4. Conclusions

Sulfa drugs are well known antibacterial agents containing N-substituted sulfonamide group on para position of aniline ring (NH2RSO2NHR′). In this study 2,4-dichloro-1,3,5-triazine derivatives of sulfa drugs (sulfamerazine, sulfaquinoxaline, sulfadiazine, sulfadimidine, and sulfachloropyrazine) (1a5a) were synthesized and characterized. Their carbonic anhydrase inhibition activity was evaluated against bCA II isozyme. For the sake of comparison the CA inhibition activity of parent sulfa drugs (1b5b) was also evaluated. It is interesting to note here that although all compounds including sulfa drugs contained an N-substituted sulfonamide group, they were still able to inhibit CA activity. Sulfa drugs exhibited weak CA inhibition (19.6–51.8%). However 2,4-dichloro-1,3,5-triazine derivatives of sulfa drugs (1a5a) had significantly improved CAI activity (IC50 = 1.49–24.9 μM) and only compound 2a had IC50 value of 105.3 μM. For compounds 1a-5a, the bCA II inhibition activity decreases in the order 5a > 1a > 4a > 3a > 2a, whereas for the parent sulfa drugs 1b5b, the inhibition activity decreases as 2b > 5b > 4b > 3b > 1b. Molecular docking studies were carried out to give insights into the binding site interactions. ADME properties were calculated to evaluate the drug likeness of synthesized compounds. All compounds indicated a favorable ADME profile.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This work was financially supported by the German-Pakistani Research Collaboration Programme. The authors are grateful to Dr. Mohammad Shahid for help with docking studies.