Surendar Chitti1, Kevin Van Calster2, Davie Cappoen2, Adinarayana Nandikolla1, Yogesh Mahadu Khetmalis1, Paul Cos2, Banoth Karan Kumar3, Sankaranarayanan Murugesan3, Kondapalli Venkata Gowri Chandra Sekhar1. 1. Department of Chemistry, Birla Institute of Technology and Science Pilani, Hyderabad Campus, Jawahar Nagar Hyderabad-500 078 Telangana India kvgc@hyderabad.bits-pilani.ac.in kvgcs.bits@gmail.com +91 40 66303527. 2. Laboratory of Microbiology, Parasitology and Hygiene, Department of Pharmaceutical Sciences, University of Antwerp Universiteitsplein 1, 2610 Wilrijk Belgium davie.cappoen@uantwerpen.be. 3. Medicinal Chemistry Research Laboratory, Department of Pharmacy, Birla Institute of Technology and Science Pilani Pilani Campus Pilani-333031 Rajasthan India.
According to the report submitted by the world health organization, globally, 1.4 million people died from tuberculosis (TB) in 2019; almost one-fourth of the world's population is infected with its etiological agent, the weakly Gram-positive Mycobacterium tuberculosis (Mtb). TB is a pulmonary infectious disease spread by aerosols and ranks among the top 10 causes of mortality worldwide.[1] TB is curable and preventable, if the patient takes appropriate treatment and precautions. Patients with smear-positive pulmonary TB (without HIV co-infection) have a 50–60% risk to die within 5 years without treatment while 20–25% achieve spontaneous resolution.[2] The remainder of these patients will develop a Latent TB Infection (LTBI), a state of persistent Mtb infection without evidence of clinically manifested active TB. However, people with LTBI risk developing active TB when they become immunocompromised later in life. Hence, development of a super-selective drug candidate to cure or prevent TB is imperative.[3] Nitrogen, oxygen, and sulphur nucleus-based heterocycles are important compounds in medicinal chemistry.[4-9] Many biologically active pharmaceutical ingredients are constructed with fused five-membered bicyclic rings. Hence fused heterocyclic moieties became valuable in novel drug discovery and medicinal chemistry. Imidazo [2,1-b]-thiazoles are reported to exhibit antitubercular,[10] antibacterial,[11] antifungal,[12,13] anticancer,[14] antiviral,[15] anti-inflammatory,[16] anthelmintic,[17] antihypertensive,[18] cardiotonic,[19,20] and insecticidal,[21] properties. Benzo-[d]-imidazo-[2,1-b] scaffold is often observed in several active pharmaceutical ingredients and bioactive molecules. Compounds containing benzo-[d]-imidazo-[2,1-b] skeleton have antimicrobial,[22] anticancer,[23] antibacterial,[24] and antiallergic properties.[25] In addition, benzo-[d]-imidazo-[2,1-b]-thiazole products exhibit a potential role as various kinase inhibitors.[26-29] Imidazo-[2,1-b]-thiazole carboxamide derivatives (A) also have notable anti-TB activity.[29,30] Imidazo-[2,1-b]-thiazole-bearing aza spiro carboxamide derivative B of imidazo-[2,1-b]-thiazole showed strong antitubercular activity.[31] Samala et al. reported anti-TB activity of a few imidazo-[2,1-b]-thiazole and benzo-[d]-imidazo-[2,1-b]-thiazole carboxamide derivatives (C) with MIC 3.53 μM and IC50 0.53 ± 0.13 μM.[32] Moraski and group recently reported that novel carboxamide derivatives of imidazo-[2,1-b]-thiazole have promising antitubercular activity, among them compound D displayed remarkable activity with MIC 0.061 μM and toxicity >100 μM (ref. 33) (Fig. 1).
Fig. 1
Structure of imidazo-[2,1-b]-thiazole and benzo-[d]-imidazo-[2,1-b]-thiazole based active antitubercular agents.
In medicinal, agrochemical research, and material sciences the triazole analogues play a key role.[34-36] Triazole is an unsaturated five membered heterocyclic compound containing three nitrogen's and two carbons. In recent years, I-A09 small molecule was identified as a potent mycobacterium protein tyrosine phosphate B inhibitor. It carries both benzofuran salicylic acid and 1,2,3-triazole analogue, which have been in clinical use since the last decade and have provided an innovative drug candidate for the treatment of Mtb.[37-39] Recently reported notable active molecules (E, F, G) exhibiting antitubercular activity based on the 1,2,3-triazole skeleton are depicted in Fig. 2.[40-42]
Fig. 2
Structure of 1,2,3-triazole based antitubercular molecules.
Pantothenate synthetase (PS) of Mtb is a homodimer protein which catalyzes the synthesis of pantothenate from d-pantoate and β-alanine. It is the last enzyme present in the pantothenate biosynthesis pathway and has been considered as crucial for the pathogenicity, persistence and survival of Mtb in the host. Generally, pantothenate can be synthesized by micro-organisms. In contrast, mammals cannot produce the molecule and need to obtain it from dietary sources. For these reasons, in recent years the enzyme made for an interesting target in developing potential agents against various forms of Mtb.[43-46] Recently, many Mtb PS inhibitors are reported by the researchers, including 2-methylimidazo-[1,2-a]-pyridine-3-carboxamide, pyrazole-[4,3-c]-pyridine carboxamide, 1-benzoyl-N-(4-nitrophenyl)-3-phenyl-6,7-dihydro-1H-pyrazolo-[4,3-c]-pyridine-5-(4H)-carboxamide, imidazo-[2,1-b]-thiazole, benzo-[d]-imidazo-[2,1-b]-thiazole derivatives and one of the widely used first-line antitubercular drug, pyrazinamide.[43,44] Moreover, our group recently offered a critical review on inhibitors of Mtb PS from a medicinal chemist's perspective.[47] Bearing in mind these findings as well as our continuous interest in the field of anti-TB drug discovery, we identified our titled compounds as potential inhibitors of Mtb PS during the molecular docking study.Piperazine moiety is also well explored for anti-tubercular activity.[32,48-50] Aiming to discover novel antitubercular compounds[35,51] and owing to the importance of these heterocyclics, in this manuscript, we focused on the synthesis of piperazine bearing benzo-[d]-imidazo-[2,1-b]-thiazole and imidazo-[2,1-b]-thiazole scaffolds with various 1,2,3-triazole analogues. Our design strategy is depicted in Fig. 3.
Fig. 3
Scaffold design strategy for the proposed novel antitubercular agents.
Results and discussion
In silico prediction of drug-likeness properties
In silico physicochemical and pharmacokinetic parameters of imidazo-[2,1-b]-thiazole (IT01–05, IT19–IT27) and benzo-[d]-imidazo[2,1-b]-thiazole (IT06–IT18) analogues were predicted using the knowledge-based FAF Drugs4 and Swiss ADME software tools to enable the selection of designed analogues that are suitable for further biological evaluation.[52,53] Narrowing the collection to a set of compounds including a high drug-likeness saves wasting money on compounds that could have negative side effects later in the drug research and development process.[54]Table 1 represents the anticipated drug-likeness parameters of each designed hybrid analogues; the prescribed cut-off ranges of each and every parameter that are being followed by 95% of known marketed drugs were added in the footnotes of the table,[49] Other than molecular weight, most of the designed imidazo-[2,1-b]-thiazole and benzo-[d]-imidazo-[2,1-b]-thiazole derivatives met the predicted n-octanol values, and the water partition coefficient, hydrogen bond acceptors and hydrogen bond donors' parameters. The values were all within the established limits and validating Lipinski's rule of five. Additionally, the number of atoms was within the acceptable range, enhancing drug-likeness. All the titled imidazo-[2,1-b]-thiazole and benzo-[d]-imidazo-[2,1-b]-thiazole's predicted aqueous solubility was good. Based on Veber's rule, the imidazo-[2,1-b]-thiazole and benzo-[d]-imidazo-[2,1-b]-thiazole analogues were also considered to have good oral bioavailability.[55] Overall outcomes disclose that the designed imidazo-[2,1-b]-thiazole and benzo-[d]-imidazo-[2,1-b]-thiazole derivatives generally possess drug-like behaviour. Hence, all derivatives were selected for further research.
In silico predicted physicochemical and pharmacokinetic parameters of the imidazo-[2,1-b]-thiazole and benzo-[d]-imidazo-[2,1-b]-thiazole derivatives
Cmpd.
MWa
log Pb
HBDc
HBAd
Heavy atomse
Rotatable bondsf
log Swg
Solubility forecast index
Oral bioavailability (Veber's rule)
IT01
401.53
3.38
8
5
28
8
−3.92
Good
Good
IT02
477.62
4.25
7
5
34
7
−5.51
Good
Good
IT03
471.66
4.88
13
5
33
13
−5.69
Good
Good
IT04
493.61
3.38
7
6
34
7
−5.03
Good
Good
IT05
472.95
3.60
7
7
32
7
−4.53
Good
Good
IT06
526.44
4.14
5
5
35
5
−6.60
Good
Good
IT07
487.58
4.21
6
6
35
6
−5.46
Good
Good
IT08
502.55
3.48
6
7
36
6
−5.48
Good
Good
IT09
485.60
4.31
5
5
35
5
−5.99
Good
Good
IT10
502.55
3.40
6
7
36
6
−5.48
Good
Good
IT11
492.00
4.17
5
5
34
5
−5.99
Good
Good
IT12
451.59
4.13
8
5
32
8
−5.16
Good
Good
IT13
527.68
4.55
7
5
38
7
−6.74
Good
Good
IT14
521.72
5.10
13
5
37
13
−6.93
Good
Good
IT15
471.58
3.87
6
5
34
6
−5.45
Good
Good
IT16
543.67
3.89
7
6
38
7
−6.25
Good
Good
IT17
523.01
4.00
7
7
36
7
−5.76
Good
Good
IT18
435.55
3.57
6
5
31
6
−4.48
Good
Good
IT19
385.49
3.07
6
5
27
6
−3.24
Good
Good
IT20
421.52
3.29
6
5
30
6
−4.22
Good
Good
IT21
486.39
3.82
5
5
30
5
−5.07
Good
Good
IT22
459.93
3.66
5
6
31
5
−4.93
Good
Good
IT23
459.93
3.66
5
6
31
5
−4.24
Good
Good
IT24
435.55
3.92
5
5
31
5
−4.75
Good
Good
IT25
452.49
2.98
6
7
32
6
−4.24
Good
Good
IT26
441.94
3.43
5
5
30
5
−4.76
Good
Good
IT27
437.52
3.61
6
6
31
6
−4.23
Good
Good
Molecular Weight (MW; 130 to 725).
Log partition coefficient between n-octanol and water (−2 to 6.5).
# hydrogen bond donors (0–6).
# hydrogen acceptors (2–20).
# Heavy atoms (non-hydrogen) (20–70).
Number of rotatable bonds (0–15).
Log aqueous solubility (−6.5 to 0.5).
Molecular Weight (MW; 130 to 725).Log partition coefficient between n-octanol and water (−2 to 6.5).# hydrogen bond donors (0–6).# hydrogen acceptors (2–20).# Heavy atoms (non-hydrogen) (20–70).Number of rotatable bonds (0–15).Log aqueous solubility (−6.5 to 0.5).
Chemistry
To achieve the title derivatives (IT01–IT27), we followed a synthetic strategy represented in Scheme 1. The cyclization of compound 1 with ethyl 2-chloro-3-oxobutanoate (2) in 1,4-dioxane under reflux condition produced the ethyl 2-methylbenzo-[d]-imidazo-[2,1-b]-thiazole-3-carboxylate (3). The obtained compound 3 was established by 1H NMR. Earlier literature reports also confirm the formation of compound 3.[29-32,51] A specific singlet peak for the methyl group protons at δ 2.28 ppm was observed. In mass spectroscopy we observed [M + 1] peak at m/z = 261.2, confirming the presence of compound 3. Later, compound 3 was further subjected to ester hydrolysis in presence of LiOH to yield corresponding acid 4; structure of 4 was confirmed by loss of the ethoxy proton signals and appearance of a broad exchangeable peak at δ 12.08 ppm in 1H NMR. The acid 4 was subjected to amide coupling with n-Boc piperazine using HATU as coupling agent and Hunig's base to get the carboxamide 5, which was confirmed by MS analysis. The carboxamide (5) was subjected to Boc deprotection in acidic medium to obtain key intermediate 6; the absence of Boc protons at 1.49 ppm in 1H NMR confirmed the formation of 6. N-alkylation of 6 with 2-bromoethanol and propargyl bromide in the presence of K2CO3 yielded the corresponding N-alkylated products (8 and 11). Compound 8 was treated with methyl silyl chloride followed by sodium azide to get the compound 9. The alkylated products (9 and 11) finally were reacted with various acetylenes and azides using well-known click chemistry to yield titled 1,2,3-triazole derivatives. In the ESI† we discussed the detailed procedure for each step, while the corresponding analytical data of the final compounds is in Section 1.
Scheme 1
Synthetic strategy for the designed molecules. Reagents and conditions: (a) 1,4-dioxane, reflux, 24 h (b) LiOH–H2O, THF–MeOH, 0 °C to rt, 16 h (c) HATU, DIPEA, DMF, rt, 12 h (d) 4 M HCl in 1,4-dioxane, CH2Cl2, 0 °C to rt, 4 h (e) K2CO3, ACN, rt, 6 h (f) i. MsCl, Et3N, CH2Cl2, 0 °C to rt, 2 h; ii. NaN3, DMSO, 70 °C, 5 h (g) CuSO4·5H2O, sodium ascorbate, DMF-t-BuOH-H2O (5 : 3 : 2), 0 °C to rt, 12–14 h.
TLC (thin-layer chromatography) and ESI mass spectrometry were used to monitor all reactions. The obtained crude product from every step was purified via column chromatography using silica gel (mesh size is 100–200) and (10 to 100%) EtOAc : petroleum ether as eluent. All title compounds were subjected to ESI-MS, 1H NMR and 13C NMR for further confirmation of structures.
Biological evaluation
The titled synthesized compounds, imidazo-[2,1-b]-thiazole (IT01–IT05, IT19–IT27) and benzo-[d]-imidazo-[2,1-b]-thiazole (IT06–IT18) were examined for their antimycobacterial activity against Mtb. A prototype based on a luminescent Mtb H37Ra lab strain (H37Ralux) was used. As previously reported in the literature, this technique provides a sensitive and repeatable tool, able to substitute enumeration of bacteria by careful plating on agar.[57,58] Antitubercular properties of the synthesized compounds were elucidated by comparing the reduction of luminescence emitted by a culture exposed to the compound to a negative control culture. After seven days of exposure of Mtb H37Ralux to successive dilutions of the compounds, the effectiveness was calculated as the minimal inhibitory concentration (MIC) at which the mycobacterial growth is reduced by 90% as well as the IC50, in which 50% of mycobacterial growth was reduced, as reported in Table 2. Within the series, hits were defined during profiling as test compounds that reached an MIC against Mtb H37Ralux below 15 μM. Albeit arbitrary by nature, the selected stringency was chosen based on the overall potency of the series. Out of the compound library, nine hits were selected, meeting the outlined criteria. The nine compounds include three derivatives (IT02–IT04) from the imidazo-[2,1-b]-thiazole cluster with MIC between 10.22 μM and 14.59 μM and six derivatives from the benzo-[d]-imidazo-[2,1-b]-thiazole cluster (IT08–IT11, IT13 and IT16) with MIC between 7.05 μM and 14.49 μM. These nine compounds were selected for a study of potential acute cellular toxicity as reported in Table 2. As a cellular model for acute toxicity, the human MRC-5 fibroblast cell line was selected.[59] The cell line was exposed to the synthesized derivatives for 24 h after which viability was measured. Viability was detected by the addition of Resazurin, a non-toxic, cell-permeable purple dye which is reduced in living cells to a pink, highly fluorescent resorufin.[60] The acute cytotoxic concentration of the test compounds was defined as the concentration at which the cells failed to reduce 50% of Resazurin compared to the non-exposed negative control cells. For the tested derivatives, acute toxicity could only be detected at concentration tenfold higher than the IC50 towards Mtb H37Ralux. For three benzo-[d]-imidazo-[2,1-b]-thiazole derivatives (IT08–IT10) and one imidaz-[2,1-b]-thiazole derivative IT04, a reduction in viability following exposure was not observed altogether at 128 μM. By dividing the CC50 with the IC50, the Therapeutic Index (T.I.) can be calculated to allow the selection of the derivative with the most favourable activity. From this analysis, the 4-nitro phenyl substituted benzo-[d]-imidazo-[2,1-b]-thiazole derivative IT10 showed best activity with an observed IC90 of 7.05 μM and an acute cytotoxicity CC50 >128 μM.
Antimycobacterial results of the synthesized imidazo-[2,1-b]-thiazole and benzo-[d]-imidazo-[2,1-b]-thiazole derivatives against H37Ra and MRC-5 strains
Cmpd.
R
R1
Mtb H37Raluxa
MRC-5b
T.I.c
IC90 (μM)
IC50 (μM)
CC50 (μM)
CC50/IC50
IT01
—
52.03
6.12
N.D.(d)
N.A.e
IT02
—
12.83
5.20
46.01
8.8
IT03
—
10.22
2.88
45.49
15.8
IT04
—
14.59
3.94
>128
32.5
IT05
—
>128
49.6
N.D.
N.A.
IT06
—
15.22
2.03
N.D.
N.A.
IT07
—
15.36
7.27
N.D.
N.A.
IT08
—
11.81
2.96
>128
43.2
IT09
—
13.49
3.46
>128
37.0
IT10
—
7.05
2.32
>128
55.2
IT11
—
14.46
3.86
29.42
7.6
IT12
—
62.76
43.67
N.D.
N.A.
IT13
—
7.24
3.19
43.67
13.7
IT14
—
27.96
13.07
N.D.
N.A.
IT15
—
27.21
12.22
N.D.
N.A.
IT16
—
14.49
7.09
58.09
8.2
IT17
—
115.57
19.74
N.D.
N.A.
IT18
—
63.77
29.13
N.D.
N.A.
IT19
—
>128
>128
N.D.
N.A.
IT20
—
54.94
20.59
N.D.
N.A.
IT21
—
53.02
20.67
N.D.
N.A.
IT22
—
57.82
5.28
N.D.
N.A.
IT23
—
84.83
43.82
N.D.
N.A.
IT24
—
48.78
19.83
N.D.
N.A.
IT25
—
51.36
23.05
N.D.
N.A.
IT26
—
62.14
30.16
N.D.
N.A.
IT27
—
98.08
46.03
N.D.
N.A.
INHd
0.09
0.04
N.D.
N.A.
Tamoxifene
N.D.
N.D.
12.40
MIC and IC50 concentrations calculated from the growth inhibition of 90% and 50%, respectively, based on luminometry of cultures exposed to a 10-point serial dilution in triplicate. The experiment was followed by two independent repeats. Percentage inhibitions >10% were repeated.
Viability of MRC-5 fibroblast cell line cultures, measured as a 50% reduction in fluorescence in a resazurin assay.
T.I., the therapeutic index calculated by dividing the CC50 by the IC50. INH, isoniazid a first-line TB drug used as reference drug. Tamoxifen, a reference compound for cellular toxicity.
N.D = Not determined.
N.A = Not applicable.
MIC and IC50 concentrations calculated from the growth inhibition of 90% and 50%, respectively, based on luminometry of cultures exposed to a 10-point serial dilution in triplicate. The experiment was followed by two independent repeats. Percentage inhibitions >10% were repeated.Viability of MRC-5 fibroblast cell line cultures, measured as a 50% reduction in fluorescence in a resazurin assay.T.I., the therapeutic index calculated by dividing the CC50 by the IC50. INH, isoniazid a first-line TB drug used as reference drug. Tamoxifen, a reference compound for cellular toxicity.N.D = Not determined.N.A = Not applicable.To discriminate between the potency of the compounds against slow versus fast growing mycobacteria, the synthesized benzo-[d]-imidazo-[2,1-b]-thiazole derivatives with an observed MIC >15 μM against Mtb H37Ra were tested against a panel of non-tuberculous mycobacteria (NTM). The panel of NTM included two members of the Mycobacterium avium complex (MAC), M. avium subsp. avium and M. intracellulare.[61] Both MAC species cause TB-like pulmonary infections in patients that are immuno-compromised or suffer from severe lung disease. Both species are considered slow growing mycobacteria with a generation time of (10–12 h) similar to the generation time of Mtb (20–24 h). As a model for fast growing mycobacteria, Mycobacterium smegmatis and Mycobacterium abscessus were selected.[62,63]M. smegmatis is considered a non-pathogenic microorganism that infrequently causes opportunistic infections while M. abscessus can cause chronic lung disease in vulnerable hosts as well as skin infections in immune deficient patients. M. abscessus infections appear to be on the rise, and outbreaks have been observed in hospitals and clinical settings around the world. Both species are considered fast-growing mycobacteria with a shorter generation time of 3 to 4 h. In general, we observed a significantly reduced potency of the (benzo-[d])-imidazo-[2,1-b]-thiazole derivatives against all NTMs as shown in Table 3. Apart from the 4-t-butyl phenyl substituted benzo-imidazo-[2,1-b]-thiazole derivative IT13, no IC50 was observed for any tested compound below 100 μM against the slow growing MAC species. Against the fast-growing NTMs, most derivatives showed activity against the fast-growing M. smegmatis but not against M. abscessus. Surprisingly, the most active compound out of the Mtb screen (benzo-[d])-imidazo-[2,1-b]-thiazole derivative IT10 showed no activity >128 μM against any NTM tested within the panel.
Antimycobacterial results of the synthesized imidazo-[2,1-b]-thiazole and benzo-[d]-imidazo-[2,1-b]-thiazole derivatives against a panel of NTMsa
Cmpd.
IC50 (μM)
M. avium subsp. avium
M. intracellulare
M. smegmatis
M. abscessus
IT02
125.8
119.3
>128.0
>128.0
IT03
>128.0
63.4
28.2
>128.0
IT04
>128.0
>128.0
16.4
>128.0
IT08
>128.0
>128.0
44.1
>128.0
IT09
>128.0
>128.0
5.7
>128.0
IT10
>128.0
>128.0
>128.0
>128.0
IT11
100.2
98.8
52.7
>128.0
IT13
17.1
>128.0
34.6
>128.0
IT16
123.5
>128.0
103.2
>128.0
IC50 concentrations calculated from the growth inhibition 50% based on luminometry of cultures exposed to a 10-point serial dilution in triplicate. The experiment was followed by two independent repeats. Percentage inhibitions >10% were repeated.
IC50 concentrations calculated from the growth inhibition 50% based on luminometry of cultures exposed to a 10-point serial dilution in triplicate. The experiment was followed by two independent repeats. Percentage inhibitions >10% were repeated.
Molecular docking studies
The crystal structure of Mtb pantothenate synthetase (PDB-3IUB) was retrieved from the RCSB Protein Data Bank[56] (https://www.rcsb.org/structure/3iub) with a resolution of 1.5 Å. After importation, the protein was prepared by Schrödinger Protein Preparation Wizard and a grid was generated. The co-crystal ligand FG-2 was extracted/removed from the PDB structure, re-docked in the generated grid, and checked for superimposition. The Root-mean-square deviation (RMSD) was found to be 0.22 Å, indicating that the docking protocol was validated and suitable for the docking study (Fig. 4) of the test molecules.
Fig. 4
Superimposed view of the native orientation of ligand FG-2 (X-ray crystallized pose) and docked orientation of the same ligand in the active site of the protein (3IUB) (Root mean square deviation 0.22 Å) (color interpretation—pink color—binding pose after docking, white color—X-ray native pose of ligand).
Molecular docking is one of the widely used structure based drug design (SBDD) techniques in computer-aided drug design approach. It is useful in drug discovery to predict the predominant binding mode(s) of a test ligand with a target protein of a known three-dimensional structure. In the current study also, in order to identify the variety of orientations, conformations, and interactions exhibited by the significantly active molecule, a docking study was initially carried out on an exemplar molecules, analog IT06 and IT10, followed by the docked poses were critically analyzed using an XP visualizer. The docking scores and the nature of interactions exhibited by the co-crystal ligand and test molecules are depicted in Table 4. By scrutinizing the co-crystal ligand's (FG-2) 3D and 2D (Fig. 5) poses, amino acid residues MET-40, HIS-47 and VAL-187, GLN-164 contributed their part in hydrogen bond formation with the co-crystal ligand and showed the docking score of −4.96 kcal mol−1. Apart from the hydrogen bond interactions, the co-crystal ligand also exhibited three aromatic bonds with HIS-47, MET-195, and GLN-164 amino-acid residues. HIS-44 amino-acid residue revealed pi–pi stacking interaction with the co-crystal ligand. Interestingly, one water molecule was also involved in the hydrogen bond formation with the ligand in the active site of the target.
Docking results and interactions of the test compounds (FG-2, IT06, IT10)
S. no.
Code
Hydrogen-bond
Aromatic bond
Pi–Pi interactions
Glide score (kcal mol−1)
1
Co-crystal ligand (FG-2)
MET-40
HIS-47
HIS-44 (Pi–Pi stacking)
−4.96
GLN-164
GLN-164
HIS-47
MET-195
VAL-187
2
IT06
HIS-44
GLN-164
VAL-187 (halogen bond)
−1.95
MET-195
HIS-47
3
IT10
VAL-187
MET-195
ARG 198 (2)
−1.29
SER 197
HIS-44
Fig. 5
(A) Docked orientations (3D and 2D) of the co-crystal ligand (FG-2) in the active site of the protein (3IUB). (B) Docked orientations (3D and 2D) of the significantly active compound (IT06) in the active site of the protein (3IUB) C. Docked orientations (3D and 2D) of the significantly active compound (IT10) in the active site of the protein (3IUB) (color interpretation of amino acid residues—yellow hydrogen bond. Blue aromatic bond, pink Pi interactions, yellow halogen bond interaction (2D)).
When the selected model compound IT06 was subjected to docking at the active site of the target protein, many interactions were observed. The amino acid HIS-44 was involved in the hydrogen bond interactions with the test compound IT06. In terms of aromatic interactions, GLN-164 and MET-195 were found common as that of co-crystal ligand. Additionally, one halogen bond was exhibited with a chlorine atom of phenyl triazole nucleus. Another compound IT10 also revealed two hydrogen bond interactions with VAL187, SER-197 amino acid residues. Additionally, MET-195 exposed an aromatic bond interaction with the target. Two Pi–Pi interactions were observed with compound IT10 with the amino acid residue ARG-198, and another with HIS-44. All these interactions made the complex best fit in the active site of the target PS (PDB-3IUB) of Mtb.
Molecular dynamic studies
The best conformations of each docked pose of co-crystal ligand and significantly active compounds IT06, IT10 were subjected to a molecular dynamics (MD) simulation study to examine the stability of the protein and ligand complex (PLC).Here, we analyzed the binding energy between target enzyme and ligand, the structural changes during complex docking, and interactions of salt bridges and hydrogen bonds. An excellent root mean square deviation (RMSD) was found for the co-crystal ligand (Fig. 6). Critical analysis was carried out to check the stability of the enzyme and co-crystal ligand; at the initial stage of the molecular dynamics, around 15 ns, minor instability was revealed. After the simulation started around 20 ns, the enzyme and co-crystal ligand were in close contact. Molecular docking and molecular dynamics studies both revealed the same amino acids involved in various types of interactions. The same amino acids, for instance, VAL-187, were wholly involved in the hydrogen bond interaction throughout the simulation. Apart from this, HIS-44, HIS-47, LYS-160, VAL-187, MET-195, SER-196 and SER-197 residues actively participated in hydrogen bond interactions with the co-crystal ligand at the active site of the target PS (PDB-3IUB) of Mtb (Fig. 7). In the same way, 49% to 99% of the MD simulated contacts were found with the co-crystal ligand FG-2. Furthermore, the SER-197 residue contributed a continuous contact of 99% throughout with an oxygen atom of the SO2 group of the ligand, and the indole nitrogen atom of the ligand collaborated 98% with the MET-195 residue of the protein. During amino-acid residue interactions analysis, water-mediated interactions were also identified. All these interactions made the complex stable during the entire duration of the study.
Fig. 6
Root mean square deviation (RMSD) of the co-crystal ligand FG-2, compound IT06 and compound IT10 complexes.
Fig. 7
Protein ligand histogram of the co-crystal ligand FG-2, compound IT06 and compound IT10 complexes.
The significantly active compound IT06 was also analyzed with the molecular dynamics simulation to confirm its conformational changes occurring during the simulation. The protein and the compound IT06 ligand complex were analyzed by checking RMSD. After critical analysis of the compound IT06, results were acceptable as a 2.8 Å distance was observed between protein and IT06 ligand complex. The amino acids HIS-44 and GLN-164 highly contributed to the hydrogen bond interaction with the protein and IT06 ligand. The same amino acids were involved in the interactions during molecular docking studies also. With these inferences, the protein and compound IT06 ligand fits perfectly in the active site of the target PS (PDB-3IUB) of Mtb. GLN-164 amino-acid residue exposed 94% of the interaction with the carbonyl oxygen of piperidinoyl group. The phenyl triazole moiety divulged two hydrogen bond interactions with HIS-44 (46%) and SER-197 (29%) amino acid residues (Fig. 8). All these simulations gave insights for future studies.
Fig. 8
Percentage of amino-acid residue interactions with the co-crystal ligand FG-2, compound IT06 and compound IT10 complexes.
The investigation was continued to another protein-ligand complex IT10. The conformational changes of the ligand was observed up to 20 ns. After 20 ns, the IT10 complex was observed stable with slight fluctuations that were within the range. The hydrogen bond interaction observed in the molecular docking studies with the amino acid residue VAL 187 revealed the same in the histogram of the molecular dynamics studies as well. The nitro group of the IT10 was contributed to the pi interaction with the amino acid residue HIS44 (46%). The keto group also involved in the active hydrogen bond formation with MET40 (53%). Another amino acid residue of the target LYS 160 (49%) was actively collaborated with the triazine moiety of the ligand IT10 with pi interaction. Overall, the molecular dynamics studies of the selected ligands were perfectly fitted in the active site of the target of PS from Mycobacterium tuberculosis.
Conclusion
A panel of around twenty-seven imidazo-[2,1-b]-thiazole (IT01-IT05 and IT19-IT27) and benzo-[d]-imidazo-[2,1-b]-thiazole (IT06-IT18) derivatives were designed, in silico ADMET predicted, synthesized in a straightforward approach by the combination of piperazine and various 1,2,3-triazole analogues. Admittedly, the potency of the synthesized panel towards Mtb was fairly low which led to arbitrary hit selection of >15 μM with the best compound selected, derivative IT10. Substitution at position 4 with a nitro phenyl moiety, as accomplished in compound IT10, lowered the observed IC90 of the benzo-[d]-imidazo-[2,1-b]-thiazole scaffold to 7.05 μM and IC50 of 2.32 μM and removed observable acute cellular toxicity >128 μM against MRC-5 lung fibroblast cell line. Another benzo-[d]-imidazo-[2,1-b]-thiazole compound (IT06) holding 2,4-dichloro phenyl moiety also exhibited significant activity (IC50 2.03 μM) against the tested strain of Mycobacterium tuberculosis. The observed activity confirms the antitubercular activity and the versatility of the imidazo-[2,1-b]-thiazoles as a potential pharmacophore. The flexibility of the scaffold should enable a straightforward library enrichment that could allow an in-depth structure activity relation (SAR) study and lead to more active library members based on the scaffolding design strategy. The low or absent toxicity towards the MRC-5 fibroblast cell line and the significantly reduced activity towards various NTM species potentially points towards a highly selective inhibitory activity exerted by the tested panel. In order to find the putative binding mode of the significantly active molecule and stability of protein-ligand complex, molecular docking and molecular dynamics studies were carried out at the active site of Pantothenate synthetase (PDB-3IUB).
Experimental section
Chemistry: general procedure for preparation of intermediates and final compounds
The general procedure of preparation of intermediates and final compounds is detailed in the ESI† section.
Authors: Maíra Bidart de Macedo; Roman Kimmel; Damijana Urankar; Martin Gazvoda; Antonio Peixoto; Freya Cools; Eveline Torfs; Luc Verschaeve; Emerson Silva Lima; Antonín Lyčka; David Milićević; Antonín Klásek; Paul Cos; Stanislav Kafka; Janez Košmrlj; Davie Cappoen Journal: Eur J Med Chem Date: 2017-06-29 Impact factor: 6.514