Design of 1-(furan-2-yl)-N-(5-substituted phenyl-1, 3, 4-thiadiazol -2-yl) methanimine derivatives as Enoyl-ACP reductase inhibitors: Synthesis, molecular docking studies and anti-tubercular activity

combination of above these pharmacophore in the designed scaffold would exhibit better antitubercular nature. Abstract A series of 5-phenylsubstiuted 1, 3, 4 thiadiazoles clubbed with furan moiety ( Fa-Fe ) by means of azomethine linkage have been synthesized. All the newly synthesized compounds were characterized by IR, 1 HNMR and Mass analyses. All the synthesized molecules have been predicted as anti-tubercular in nature by PASS in silico approach. In vitro anti-tubercular screening was performed by alamar blue assay method on Mycobacterium tuberculosis H37Rv strain. Among the synthesized derivatives Fb and Fe were active at 3.125 µg/mL against Mycobacterium tuberculosis H37Rv strain. The mechanism of action of the active compounds was carried out by docking of receptor enoyl-ACP reductase. It has been concluded that both Fb and Fe posses a significant interaction of hydrogen bonding and electrostatic attraction with Tyr 158 and Met103 in the active site of enzyme.


Introduction
Tuberculosis (TB), a lung infection is one of the contagious and deadly diseases which have added to the woes of the mankind.It is the most dreadful infectious disease caused by Mycobacterium tuberculosis.The widespread of this disease is primarily due to the population growth, emergence of multi-drug resistant TB strains, financial burden in the developing countries and unsuccessful attempt to synthesize a new drug with novel mechanism of action (Thomas et al., 2011).Mycolic acid biosynthesis is essential for the building up of cell wall in mycobacterial and related species (Kolattukudy et al., 1997).InhA, the enoyl acyl carrier protein reductase from Mycobacterium tuberculosis is one of the key enzymes involved in the mycobacterial fatty acid elongation cycle and has been validated as an effective antimicrobial target (Banerjee et al., 1994).A series of pyrrolidine carboxamides have been recog-nized as potent direct class of InhA inhibitors (He et al., 2006).Recently Oxadiazolo pyrrolidine carboxamides were designed as enol-ACP reductase inhibitors (Sonia et al., 2012).Considering the structural features of the inhibitor of enol-ACP reductase enzyme, we had undertaken the design a novel class of 1,3,4-thiadiazole clubbed with the furan moiety by means an azomethine linkage (Figure 1).

Materials and Methods
Melting points were determined by using open capillary tube method and the values were uncorrected.IR spectra were recorded on JASCO FT/IR-140 spectrophotometer by using KBr pellets technique.PMR spectra were recorded using BRUCKER FT-NMR-500 MHz spectrophotometer by using DMSO as solvent and TMS as internal standard.The chemical shift was expressed in  ppm.Mass spectra were recorded on a JEOL GCmate mass spectrometer.

Synthesis of thiosemicarbazones
Aromatic aldehyde (0.2 M) in warm alcohol (300 mL) and thiosemicarbazide (0.2 M) in warm water (300 mL) were mixed slowly with continuous stirring.The product separated immediately on cooling which was filtered with suction, dried and recrystallized in ethanol: water (7:3).

Synthesis of 5-substituted phenyl-1, 3, 4-thiadiazol-2amines
Thiosemicarbazone (0.05 M) was suspended in 300 mL warm water, FeCl3 (0.15 M) in 300 mL water was added quantitatively, slowly with constant stirring.The contents were heated at 80-90°C for 45 min.Solution was filtered hot and then citric acid (0.11 M) and sodium citrate (0.05 M) were added.The resulting mixture was divided into 4 parts and each part was neutralized separately with ammonia (10%).The required amine separated out, filtered with suction, dried and recrystallized in ethanol:water (5:5).

PASS prediction
Prediction of this spectrum by PASS is based on SAR analysis of the training set containing more than 35,000 compounds which have more than 500 kinds of biological activity.If Pa>0.7 the chance to find the activity in experiment is high, but in many cases the compound may occur to be the close analogue of known pharmaceutical agents.If 0.5<Pa<0.7 the chance to find the activity in experiment is less, but the compound is not so similar to known pharmaceutical agents.If Pa<0.5 the chance to find the activity in experiment is even more less, but if it will be confirmed the compound might occur to be a New Chemical Entity.It was shown that the numbers of experiments necessary to find the known effects for the heterogenous set of substances by using PASS prediction and without it are 1650 and 316 respectively.(Poroikov et al., 2002;Mathew et al., 2013).

Preparation of proteins and ligands
The crystal structure of Enoyl Acyl Carrier Protein Reductase (InhA) from Mycobacterium tuberculosis (PDB The geometry optimization was done by using UFF molecular mechanics method.Clean hybridization option was done and make sure with exact hybridization pattern of the molecule.Final geometry optimization was done by using semi empirical Quantum mechanics PM3 method.The energy minimized structures were saved in pdb format for further docking studies.(Mathew et al., 2012).

Determination of the active site of the enzyme
The active sites of the (PDB entry: 2HYM) were identified by using Q-Site Finder: an energy-based method for the prediction of protein-ligand binding sites (Laurie et al., 2005) VAL 189,ALA 191,GLY 192,PRO 193,ILE 194,THR 196,MET 199,ILE 202,LEU 207 were the interacting residues.These predicted amino acid residues were selected and saved as the binding site for the docking study of final designed derivatives.

Docking calculation
Argus lab 4.0.1.versionwas employed in the docking between protein and ligand using argus dock with a fast and simplified potential of mean force (PMF).Dock a ligand in to binding site was done by 'Argus Dock' as the docking engine.'Regular precision' was selected in docking precision menu, 'Dock' was selected as calculation type, 'Flexible' for the ligand and 'Ascore' was used for the scoring function.The binding site box was set in to (24 x 22 x 22) angstroms for enclosing the entire active binding site of the protein with a grid resolution of 0.4Å.The process of docking is repeated until a constant value of docking score is obtained.The resulted final docked structures were saved in pdb format and the docking snap shots were imported in to Molegro virtual viewer version 2012 2.5.0.This program can generate the amino acid residue around each segment of the final synthesized derivatives and interpreting the different drug-receptor interaction such as hydrogen bonding, hydrophobic interaction and electrostatic force of attraction.

Validation of potential means force method
It is necessary to perform the docking of 1-cyclohexyl-N -(3,5-dichlorophenyl)-5-oxopyrrolidine-3-carboxamide into the active site of enoyl-ACP reductase (2H7M), before docking the test compounds for validate the molecular docking program prescribed in the current study.This inhibitor of binds in the active site region of enoyl-ACP reductase with a binding score of -11.7 with a root mean square of 0.23 was observed.105number of final unique configurations appeared during the docking calculation.The major interaction of the inhibitor is shown in (Figure 2).

Anti-tubercular activity
The anti-mycobacterial activity of synthesized compounds (Fa-Fe) were assessed against M. tuberculosis using microplate Alamar Blue assay (MABA).This methodology is non-toxic, uses a thermally stable reagent and shows good correlation with proportional and BACTEC radiometric method.Briefly, 200 µL l of sterile deionzed water was added to all outer perimeter wells of sterile 96 wells plate to minimized evaporation of medium in the test wells during incubation.The 96 wells plate received 100 µL of the Middlebrook 7H9 broth and serial dilution of compounds were made directly on plate.The final drug concentrations tested were 100 to 0.2 µg/mL.Plates were covered and sealed with parafilm and incubated at 37ºC for five days.After this time, 25 µL of freshly prepared 1:1 mixture of Almar Blue reagent and 10% tween 80 was added to the plate and incubated for 24 hours.A blue colour in the well was interpreted as no bacterial growth, and pink colour was scored as growth.The MIC was defined as lowest drug concentration which prevented the colour change from blue to pink (Collins et al., 1997).
Thiosemicarbazone was obtained by the reaction between aromatic aldehyde and thiosemicarbazide.This thiosemicarbazones undergo oxidative cyclization with FeCl3 in presence of citric acid medium to form 5phenyl 2-amino 1,3,4-thiadiazole.The amino group of the thiadiazole ring undergoes nucleophilic addition followed by dehydration with furfural in presence of conc.H2SO4 with DMF medium (Mathew et al., 2011).The presence of a band at 1618 cm-1 in Fc strongly recommend the presence of azomethine linkage.A band of 721 and 1350 cm-1 indicated the presence of C-S -C and C-O groups.In 1 HNMR a singlet peak at 8.55 showed the presence of N=CH unit in Fc.The deshielding of this CH signal is due to the strong inductive effect of neighbouring electronegative nitrogen atom.The appearance of a singlet peak in the region of 3.55 showed the presence of CH3 group in the methoxyl group of Fc.The multiplet peaks ranges from 7.0-7.90revealed the presence of aryl protons in the final synthesised derivative.The mass fragment of Fe showed a molecular ion peak M + at m/z 300 correspond to molecular formula C13H8N4O3S.A peak of m/z 301 of M + +1 was also observed.The splitting of a oxygen atom group from the Fe derivative showed a peak of 285 which correspond to the molecular formula of C13H8N4O2S gave a full agreement of the incorporation of the imines linkage between furfural and 2-amino 5(3-nitro phenyl) 1,3,4-thiadiazole.254, 233, 206, 106 and 94 are the other prominent peaks obtained.The mass fragmentation pattern of Fe is shown in (Figure 4).
The predicted activity spectra of the 1-(furan-2-yl)-N-(5substiuted-phenyl-1,3,4-thiadiazol-2-yl) methanimines (Fa-Fe) in the current study was done to identify its anti -tubercular nature with a Pa scores more than 0.68.The activity score of the all designed molecules were shown in Table I.
The prediction of the pharmacological profile was established by the in vitro studies.The rank order of PASS showed as Fe>Fa>Fd >Fb>Fc.The PASS program significantly predicted the Fe compound has the highest probability and was most active in the in vitro anti-tubercular studies, while Fc with lowest Pa value in the series was practically less active.It has been noted that the PASS program predicts Fb with a low Pa value of 0.691.This low score cannot favours its activity ratio because it showed significant activity with MIC of 3.3 µg/mL.
The active pocket was considered to be the site where 1cyclohexyl-N-(3, 5-dichlorophenyl)-5-oxopyrrolidine-3carboxamide (pyrrolidine carboxamides) complexed with enoyl-ACP reductase in 2H7M.Hydrogen bonding network of ligand and Tyr158 seems to be a prominent feature among all the InhA-inhibitor complexes identified so far (Kuo et al., 2003).All the newly designed molecules showed good binding interaction towards the active site of enoyl-ACP reductase.Molecular docking study suggested that compound Fb and Fe showed the best docking score of-9.crucial π-π stacking interaction between the azomethine unit and phenyl system of Tyr158 (Figure 5).
The in silico study of the Fb and Fe revealed that its binding affinity towards the enzyme is due to the hydrogen bonding interaction with Tyr 158.The ligand map of Fb (Figure 6) suggested that a steric interaction of furan part in the scaffold to the Met 103 also contribute an interesting result in the activity ratio.Tyr 158 and Met 103 are the major interacting residues of the known inhibitor of 1-cyclohexyl-N-(3,5dichlorophenyl)-5-oxopyrrolidine-3-carboxamide towards enoyl-ACP reductase.
Another interesting factor raised from our finding is that the poor affinity and MIC of Fc is due to its lack of hydrogen bonding between Tyr158.Estimated binding energy (Kcal/mol) of all the designed molecules and the amino acid residues enveloped by the designed molecules was represented in Table II.

Figure 1 :
Figure 1: Design of enoyl-ACP reductase inhibitor based up on known pyrrolidine carboxamides

Figure 3 :
Figure 3: Synthetic route of the titled compounds

FigureFigure 5 :
Figure 4: Mass fragmentation pattern of Fe Journal homepage: www.banglajol.infoAbstracted/indexed in Academic Search Complete, Agroforestry Abstracts, Asia Journals Online, Bangladesh Journals Online, Biological Abstracts, BIOSIS Previews, CAB Abstracts, Current Abstracts, Directory of Open Access Journals, EMBASE/Excerpta Medica, Google Scholar, HINARI (WHO), International Pharmaceutical Abstracts, Open J-gate, Science Citation Index Expanded, SCOPUS and Social Sciences Citation Index ISSN: 1991-0088