Sulfa Drugs as Inhibitors of Carbonic Anhydrase: New Targets for the Old Drugs

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) (1a–5a) 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 (1b–5b) 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.


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 (Zn 2+ ) 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 CO 2 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 [9][10][11], 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 [13][14][15][16]. 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 [17][18][19][20][21][22].
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 (NH 2 RSO 2 NHR ) [24,25]. The carbonic 2 BioMed Research International 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 (1a-5a) 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 (1b-5b) is also reported.

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 HF 254 , 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 1 H and 13 C-NMR analysis Bruker Avance DRX500 spectrometer was used with TMS as an internal standard and DMSO-d 6 as solvent.

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.       [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.  [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

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 dior tri-substituted products. The structures of sulfa drugs used 1b-5b are given in Figure 2.

bCA II Inhibition Studies and SAR.
The synthesized compounds (1a-5a) and their parent sulfa drugs (1b-5b) 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 (CO 2 ) 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 (1b-5b) 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 IC 50 value of 1.49 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 1b-5b, 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 (SO 2 NH 2 ), which is the typical zinc binding function that binds directly to the Zn 2+ 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  Figure 3: Compound 5a inside binding pocket of hCA II. 5a is shown in magenta color with sticks representation. The hydrophobic surfaces are drawn using Chimera and the color scheme is according to Kyte-Doolittle coloring scheme; blue color indicates most hydrophilic surface (hence polar residues), white for neutral, and orange for hydrophobic surface (nonpolar residues).
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 Zn 2+ 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 1a-5a, 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 1b-5b ( Figure 6). In most binding conformations for 1a-5b, 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 1a-5a are better able to inhibit CA II as compared to compounds 1b-5b, 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 ( + log and log ), the pH dependent octanol-water distribution coefficient ( + log ), 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 log value of less than 5 are considered to be orally bioavailable with favorable ADME profile [39,40]. The calculated ADME properties of compounds 1a-5a and 1b-5b are given in Table 3. All synthesized compounds showed favorable ADME properties.
No violation of Lipinski's rule of 5 was observed.

Conclusions
Sulfa drugs are well known antibacterial agents containing N-substituted sulfonamide group on para position of aniline ring (NH 2 RSO 2 NHR ). In this study 2,4dichloro-1,3,5-triazine derivatives of sulfa drugs (sulfamerazine, sulfaquinoxaline, sulfadiazine, sulfadimidine, and sulfachloropyrazine) (1a-5a) 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 (1b-5b) 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 (1a-5a) had significantly improved CAI activity (IC 50 = 1.49-24.9 M) and only compound 2a had IC 50 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 1b-5b, 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.