e-Pharmacophore model-guided design of potential DprE1 inhibitors: synthesis, in vitro antitubercular assay and molecular modelling studies

Tuberculosis continues to wreak havoc worldwide and caused around 1.4 million deaths in 2019. Hence, in our pursuit of developing novel antitubercular compounds, we are reporting the e-Pharmacophore-based design of DprE1 (decaprenylphosphoryl-ribose 2′-oxidase) inhibitors. In the present work, we have developed a four-feature e-Pharmacophore model based on the receptor–ligand cavity of DprE1 protein (PDB ID 4P8C) and mapped our previous reported library of compounds against it. The compounds were ranked on phase screen score, and the insights obtained from their alignment were used to design some novel compounds. The designed compounds were docked with DprE1 protein in extra-precision mode using Glide module of Maestro, Schrodinger. Some derivatives like B1, B2, B4, B5 and B12 showed comparable docking score (docking score > − 6.0) with respect to the co-crystallized ligand. The designed compounds were synthesized and characterized. In vitro antitubercular activity was carried out on Mycobacterium tuberculosis H37Rv (ATCC27294) strain using the agar dilution method, and minimum inhibitory concentration (MIC) was determined. The compound B12 showed a MIC value of 1.56 μg/ml which was better than the standard drug ethambutol (3.125 μg/ml). Compounds B7 and B11 were found to be equipotent with ethambutol. Cytotoxicity studies against Vero cell lines proved that these compounds were non-cytotoxic. Molecular dynamic simulation study also suggests that compound B12 will form a stable complex with DprE1 protein and will show the crucial H-bond interaction with LYS418 residue. Further in vitro enzyme inhibition studies are required to validate these findings.


Introduction
More than 40 years have passed since the introduction of the quadruple drug therapy regimen for tuberculosis (TB), which has proven to be both effective and cost-effective. But despite this, TB continues to wreak havoc worldwide which is evident from the fact that worldwide around 1.4 million died of TB in 2019 (Raviglione et al. 2012; Global tuberculosis report 2020). This state of high mortality rate feeds on various factors like prolonged duration of antitubercular therapy, toxicity induced by the antitubercular drugs, lower patient compliance and the rise of MDR-TB (multidrugresistant). Bedaquiline, pretomanid and delamanid are the only drugs approved by the FDA (Food and Drug Administration) in the last four decades to treat MDR-TB which highlights the associated challenges in the field of antitubercular drug discovery (Cohen 2013). Hence, it is the urgent need Avinash Kumar and Revathi Rajappan have contributed equally to this work. of the hour to develop novel small molecules as antitubercular compounds acting on newer targets of the mycobacterium. Among many of the targets explored so far, DprE1 (Decaprenylphosphoryl-ribose 2′-oxidase) has emerged as a promising target for the development of antitubercular drugs. DprE1 is necessary for the synthesis of the mycobacterial cell wall as it facilitates the conversion of DPR (decaprenylphosphoryl-D-keto-erythro-pentofuranose) to DPX (decaprenylphosphoryl-2-ketoribose). DPX is further reduced to DPA (decaprenylphosporyl arabinose) by DprE2 as shown in Fig. 1 (Mikušová et al. 2005;Brecik et al. 2015).
The arabinosyltransferases use DPA as a substrate to synthesize arabinogalactan and lipoarabinomannan, which are major components of the mycobacterial cell wall. DprE1 is an oxidoreductase enzyme-containing FAD (flavin adenine dinucleotide) (Neres et al. 2012;Batt et al. 2012). DprE1 as a potential antitubercular drug target was validated by Makarov et al., who reported BTZ043 (benzothiazinone derivative) as the first covalent inhibitor of DprE1 (Makarov et al. 2009). After this, several authors have reported various heterocyclic derivatives like benzothiazole, azaindoles, triazoles, aminoquinolones, dinitrobenzamides, quinoxalines, pyrazolopyridines, etc. showing antitubercular activity through the inhibition of DprE1 (Landge et al. 2015;Mir et al. 2014;Magnet et al. 2010;Batt et al. 2016;Panda et al. 2014). Readers can refer to some excellent reviews by Chikhale et al. and Gawad et al. for more insights into the development of DprE1 inhibitors (Chikhale et al. 2018;Gawad and Bonde 2018).
This work is in continuation of our previous reported piperidone derivatives as potential anti-HIV/TB agents (Kumar et al. 2019). In the present work, we have designed non-covalent DprE1 inhibitors by developing an e-Pharmacophore model by using the receptor-ligand interactions of DprE1 with 2-carboxyquinoxalines as reported in its X-ray crystal structure (Neres et al. 2015). Several authors have reported e-Pharmacophore modelling combined with molecular docking and molecular dynamic simulations studies to design and identify potent inhibitors of several proteins (Arun et al. 2020;Kumar et al. 2020Kumar et al. , 2021Salam et al. 2009). We have used PDB ID 4P8C to develop a fourfeature e-Pharmacophore model and mapped it against our previous reported antitubercular compounds (Kumar et al. 2019). We found that only one aromatic feature matched, and our compounds had very low phase screen score. Based on the insights from the developed e-pharmacophore model, we added a benzimidazole group to the nitrogen of piperidone ring and the resultant derivatives had three out of four features matches with an improved phase screen score. Benzimidazole ring system was selected because historically its derivatives have been reported as antitubercular agents, and recently, Manjunath et al. have found a series Fig. 1 Epimerization of 2′-hydroxyl group of DPR by DprE1 and DprE2 and rationale behind selecting DprE1 as a potential target for antitubercular drug design of benzimidazole derivatives showing DprE1 inhibition (Manjunatha et al. 2019;Keri et al. 2016;Yeong et al. 2017;Surineni et al. 2019). Detailed design strategy is depicted in Fig. 2 and described under the materials and methods section. Herein, we have reported some potent antitubercular compounds with activity in the lower micromolar range and this design strategy can be used for targeting other proteins as well.

Materials and methods
The chemicals, reagents and solvents for the present work were sourced from Aldrich (Bengaluru, India), Spectrochem P. Ltd. (Mumbai, India), E. Merck (Mumbai, India), S. D. Fine Chem. P. Ltd. (Mumbai, India). TLC (thin-layer chromatography) was done on 0.25 mm silica gel plates purchased from E. Merck (silica gel 60 F254). The melting point was measured by using a laboratory melting point apparatus from Toshniwal P. Ltd., Mumbai, India. IR spectra were recorded on FT-IR Affinity-1 IR spectrometer from Shimadzu, Bengaluru, India. Mass spectra of the synthesized compounds were recorded on GC-MS-QP5050A from Shimadzu. NMR spectra were recorded on a 400 MHz spectrometer from Bruker, and dimethyl sulphoxide (DMSO)-d6 was used as the solvent. All in silico experiments were done on an HP computer with Linux (Ubuntu, version 18.04.1) as the operating system and Maestro 11.6 which is a small molecule drug discovery suite from Schrödinger LLC, 2018-3 (Schrödinger LLC, Bengaluru, India).

Design of non-covalent DprE1 inhibitors
For the present work, we have employed e-Pharmacophore modelling based on the receptor-ligand complex of the noncovalent inhibitor of DprE1 protein as reported by Neres et al.(Neres et al. 2015). The e-Pharmacophore method generates energetically optimized structure-based pharmacophores and has the advantages of both ligand-and structurebased pharmacophore modelling Loving et al. 2009). The e-Pharmacophore model was developed employing PHASE module of Maestro suite (Dixon et al. 2006). PHASE has a built-in set of six pharmacophore features like hydrogen bond donor (D), hydrogen bond acceptor (A), hydrophobic group (H), aromatic ring (R), positively ionizable (P) and negatively ionizable (N). The screened compounds were ranked based on the phase screen score. DprE1 protein (PDB ID 4P8C) was imported from the protein data bank followed by pre-processing, removal of water and energy minimization using the protein preparation wizard of Maestro (Madhavi Sastry et al. 2013). For all further studies, this energetically minimized protein structure of DprE1 was used. To further understand the binding of the designed compounds with DprE1 protein, molecular docking simulations were carried out.

Molecular docking and free binding energy calculations
Molecular docking studies are useful in understanding the binding between a protein and ligand. For our study, we carried out extra-precision (XP) docking employing Glide module of Maestro 11.6 interface (Friesner et al. 2006). In the extra-precision mode, false positives are eliminated by extensive sampling and more advanced scoring functions. XP-docking was accomplished in three steps, viz. receptor grid generation, ligand minimization and docking. On the minimized protein structure of DprE1, a receptor grid was generated around the co-crystallized ligand as the site of binding of the bound ligand was the active site where our designed ligands would bind with DprE1. Then the designed compounds were energetically minimized using LigPrep tool of Maestro. In the final step, the structures generated from the LigPrep exercise were docked using the receptor grid file generated in the first step. The docked compounds were ranked based on the docking score. The docking protocol was validated by redocking the co-crystallized ligand and then computing the RMSD (root mean square deviation) value.
Free binding energy calculations were also carried out employing MM-GBSA (molecular mechanics, the generalized Born model and solvent accessibility) under the Prime module (Jacobson et al. 2004). Prime MM-GBSA generates many energy properties which report energies for the receptor, ligand and other complex structures. It can also compute the energy differences produced due to strain and binding. It employs VSGB 2.0 as an optimized implicit solvent model. Five fundamental energy are calculated in Prime, and from these energies, the strain and binding energies are computed (Li et al. 2011). For simplifying our work, we have considered MM-GBSA dG bind expressed in the unit of kcal/mol to predict the binding energy of our designed compounds with the DprE1 protein. For this study, the resultant docked complex files from Glide were directly put for MM-GBSA calculations. Based on the docking score and MM-GBSA dG bind score, twelve compounds were selected for synthesis. The docking score of the selected compounds ranged from high to low in comparison with the co-crystallized ligand.

Synthesis of 1-(1H-benzimidazol-2-ylmethyl) piperidin-4-one (3)
As shown in Scheme 1, 4-piperidone hydrochloride (0.7 gm, 0.005 mol) and 2-chloromethylbenzimidazole (0.81 gm, 0.005 mol) were taken in a round-bottomed flask containing 10 ml of DMF (dimethyl formamide). To this reaction mixture, 0.7 ml (0.005 mol) of TEA (triethylamine) was added and was refluxed for 7 h. The reaction progress was measured by taking TLC at periodic intervals. After the completion of the reaction, the contents were transferred into crushed ice. This mixture was further refrigerated for another 12 h. Finally, the precipitate was filtered under vacuum and recrystallized using methanol to obtain the pure product.

The general method of procedure for the synthesis of B2-12
The synthesis of compounds B2-B12 (Scheme 1) was carried out by using the method reported by Banerjee et al. (Banerjee et al. 2011). Different amines (0.01 mol) and sodium acetate (0.8 gm, 0.01 mol) were added to the ethanolic mixture of compound B1 (3.85 gm, 0.01 mol). The reaction mixture was kept for reflux with continuous stirring using a magnetic stirrer. The reaction was monitored using TLC, and after the completion (after 3.5 h), it was brought to room temperature and then transferred into icecold water. Finally, the precipitate was filtered under vacuum and recrystallized using methanol DMF (2:8) to obtain the pure product.

In vitro antitubercular evaluation
Antitubercular evaluation of the synthesized derivatives was carried out according the procedure reported by Addla et al. (Addla et al. 2014). In brief, serial dilution of the synthesized compounds and ethambutol (the standard drug used in this study) was carried out followed by inoculation into an agar medium (Middlebrook 7H11) supplemented with OADC (oleic acid, albumin, dextrose and catalase) procured from Difco (Kolkata, India). Inoculum; M. tuberculosis H37Rv ATCC27294 was made from freshly prepared Middlebrook 7H11 agar slants supplemented with OADC. An approximate concentration of 10 7 cfu/ml of mycobacterial suspension was obtained and 5 µl of this suspension was spotted into 7H11 agar tubes containing graded concentrations of the synthesized derivatives and ethambutol. It was incubated at 37 °C, and readings were taken after 28 days (Hall et al. 2012).

Cytotoxicity assay
Cell toxicity studies for these compounds were carried out on Vero cell line (African green monkey kidney cells) obtained from NCCS (National Centre for Cell Sciences), Pune, India. The assay was carried out using the well stablished protocol reported by Mosmann et al. and Denizot et al. (Mosmann 1983;Denizot and Lang 1986). The growth inhibition percentage was computed by the following formula: Further, SI (selectivity index) was computed using the formula ( SI = CC50∕MIC ) (Protopopova et al. 2005).

Molecular dynamic simulations study
To better understand the interaction of the most potent antitubercular compound from this study with DprE1 protein, MD simulations studies were carried out with the help of Desmond module (Schrödinger Release 2018-3) (Bowers et al. 2006). MD simulations were carried out in three steps which included system builder, minimization and simulation. In the system builder step, the DprE1 protein-ligand complex was solvated using SPC as a solvent system in an orthorhombic box shape. This solvated complex was then minimized employing SD method with iterations fixed at 2000. The system was slowly equilibrated at 300 K (temperature) and 1.01325 bar (pressure) with the help of Nose-Hoover thermostat method and Martina-Tobias-Klein barostat method. In the last step, the solvated and minimized complex was simulated for 20,000 ps.

Design of non-covalent DprE1 inhibitors
Based on the non-covalent interactions between the amino acid residues in the active pocket of DprE1 protein (PDB ID 4P8C) and the co-crystallized ligand with the protein, i.e. 6-(trifluoromethyl)-3-{[4-(trifluoromethyl)benzyl]amino} quinoxaline-2-carboxylic acid, we developed a four-feature e-Pharmacophore model (ANRR) as shown in Fig. 2. Recently, Zhang et al. have also reported a three-feature pharmacophore (3-feature) model which had two hydrophobic and one hydrogen bond acceptor feature (Zhang et al. 2018). They had selected trifluoromethyl group substituted at the 6th position of quinoxaline ring as one of the hydrophobic features, while we have selected the fused phenyl ring of the quinoxaline ring system. The distance reported by Zhang et al., between the two hydrophobic groups, was different than ours. We found it to be 6.65 Å while Zhang et al. have reported it to be around 10.5 Å. The inter-site angle and inter-site distance are depicted in Fig. 2a and b. Since the developed model had features of a known non-covalent inhibitor of DprE1, it was assumed that the potential DprE1 inhibitors should have similar features. The developed model had two aromatic ring features (R7 and R9) at 6.65 Å and an acceptor feature (A2) at 6.41 Å from R7 and 2.77 Å from R9. A2 was at 2.97 Å from the negatively charged feature (N6) while N6 was at 5.64 Å from R9 and 6.49 Å from R7. The angle was computed between three features with the second feature at the centre. The angle between A2-R9-R7 was 72.9°, R9-R7-A2 was 24.4°, R9-R7-N6 was 50.8°, R7-N6-A2 was 75.1°, and R9-A2-N6 was 158.3°. Based on this insight, we propose that any ligand with all these four features at this specific distance and angle might be a potential DprE1 non-covalent inhibitor.
The developed e-Pharmacophore model was used to map ligands to assess whether they possessed pharmacophoric features essential for DprE1 inhibition. In our previous work, we had found some piperidone derivatives showing good antitubercular activity. We mapped all those 25 compounds on the developed e-Pharmacophore model and found none of the compounds had even two matching features. The best compound was R7 ((3E,5E)-3,5-Bis(furan-2-ylmethylidene)-4-(2-phenylhydrazinylidene)piperidine) with one match and a phase screen score of 0.245. The detailed design strategy is depicted in Fig. 3, where we found that if -NH group of piperidone ring could be substituted, it could increase the matches. As discussed in "Introduction", benzimidazole derivatives have been reported as DprE1 inhibitors with potent antitubercular activity; hence, we added benzimidazole ring system and this indeed led to increasing in the matches of features as well as phase screen score. All the 25 derivatives were substituted with 1H-benzo[d]imidazol-2-yl) methyl group. This exercise led to an increase in the number of matches to three out of four features. N6 pharmacophoric did not match, but there was a substantial increase in the phase screen score of the compounds. We decided to synthesize all those compounds which had at least three matches and a phase screen score more than 1.60 and twelve compounds passed this test. Derivative B1, i.e. (3E,5E)-1-((1Hbenzo[d]imidazol-2-yl)methyl)-3,5-bis(furan-2-ylmethylene)piperidin-4-one, showed the highest screen score of 1.84 while compound B4, i.e. 2-((3E,5E)-1-((1H-benzo[d] imidazol-2-yl)methyl)-3,5-bis(furan-2-ylmethylene)piperidin-4-ylidene)hydrazine-1-carbothioamide, showed the least screen score of 1.60. The rest of the compounds had a phase screen score between 1.84 and 1.60 as shown in Table 1.

Molecular docking and free binding energy calculations
Molecular docking study (XP mode) was done to understand the binding of the designed compounds with DprE1 protein.
The docking protocol was validated by redocking the cocrystallized ligand at its reported binding site and then by calculating the RMS value by superimposing both the poses. As shown in Fig. 4, both poses of the co-crystallized ligand were similar to an RMS value of 0.5149. The compounds were ranked by their docking score, and non-bonding amino acid interactions were compared with the co-crystallized ligand. As shown in Fig. 5f, the co-crystallized ligand showed H-bond interaction with LYS418 amino acid residue. It also showed hydrophobic interactions with TYR60, TRP230, PRO316, LEU317, LUE363, VAL365, PHE369 and CYS387 residues. It showed a docking score of − 6.728 and free binding energy of − 47.92 kcal/mol as computed by prime MM-GBSA and shown by MM-GBSA dG bind. Compound B1 which had the highest phase screen score also showed the highest docking score of − 7.644, even better than the co-crystallized ligand. The dG bind value of B1 was also more than − 49.64 kcal/mol. It also showed possible H-bond interactions with LYS418 residue and its other interactions were also like the co-crystallized ligand. The docking score for other five compounds like B2, B3, B4, B5 and B12 was more than − 6.0 and was comparable to the co-crystallized ligand. Compound B3 showed potential H-bond interaction with GLN336 residue and not with LYS418 residue. The rest of the compounds had docking scores below − 6.0 and more than − 2.50. Docking study and MM-GBSA free binding studies suggest that few of the designed compounds like B1, B2, B3, B4, B5 and B12 might form a stable complex with DprE1 protein. Ligand interaction diagram (2D) is depicted for the top five ranked compounds in Fig. 5a-e.

Chemistry
The reaction between piperidin-4-one (1) and 2-chloromethylbenzimidazole (2) resulted in the synthesis of 1-(1H-benzimidazol-2-ylmethyl)piperidin-4-one (3). In this reaction, chloro group of the 2-chloromethylbenzimidazole acts as a good leaving group. This leads to a nucleophilic attack on the methylene carbon of compound 2 by the secondary amine of 4-piperidone. Hydrochloric acid formed as a by-product of this reaction is neutralized by TEA thus driving the reaction forward. In the next step, intermediate 3 undergoes aldol condensation when treated with 2-furfural and NaOH to form an aldol product which upon loss of water molecule forms the first derivative B1. In the last step, B1 was treated with different substituted amines to yield compounds B2-B12 through sodium acetate catalysed Schiff's reaction.

In vitro antitubercular evaluation
Agar dilution method was employed for the in vitro antitubercular assessment of the synthesized derivatives against M. tuberculosis H37Rv (ATCC27294) strain, and minimum inhibitory concentration (MIC) was determined in triplicates. The MIC values (μg/ml) are given in Table 1.
The compound B12 showed the most potent antitubercular activity with a MIC of 1.56 μg/ml which was better than the standard drug ethambutol (3.125 μg/ml). Compounds B7 and B11 were found to be equipotent with ethambutol and had a MIC value of 3.125 μg/ml. Compound B1 had a MIC value of 50 μg/ml which was the least potent among all the tested compounds. Compounds B2, B4, B5, B6 and B10 showed a MIC value of 6.25 μg/ml which was comparable with ethambutol activity. It was observed that nonsubstitution at the 4th position of piperidone ring like in B1 resulted in the least potent compound of the series. Hence, substitution at 4th position was favourable for the activity. The most active compound B12 had o,p-difluoro substitution at the phenyl ring attached to the 4th position of the piperidine ring. Similarly, B11 with o,p-dichloro substitution was equipotent with ethambutol. From this, we can conclude that electron-withdrawing groups at the phenyl ring on o,p-position led to an increase in the activity. Single fluoro substitution at para position (B6) was less potent than o,p-difluoro substitution (B12). Hence, di-substitution at o,p-position was more favourable. Electron-donating groups like methyl (B9) and methoxy (B8) at the phenyl ring led to a decrease in the activity as both compounds showed a MIC value of 25 μg/ ml. From this study, we can conclude that the piperidine ring substituted with furylidene rings and methyl benzimidazole ring showed potent antitubercular activity. But if we compare the antitubercular activity of parent compounds without substitution at the 4th position of piperidine ring like R2, R3, R8, R9 and R11 for which the corresponding compounds are B2, B3, B8, B9 and B11, we found that the activity was same in both the cases (Kumar et al. 2019). The compounds like B4, B5, B6 and B12 were twice less potent than their unsubstituted derivatives. Hence, we can conclude that the substitution at the 4th position of piperidine ring was well tolerated.

Cytotoxicity studies
As antitubercular drugs should be non-toxic to the normal cells; hence, all the synthesized derivatives were tested for their cytotoxicity potential using 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyl tetrazolium bromide (MTT) assay against Vero cell lines. MTT is a colorimetric assay where yellowcoloured MTT dye is reduced to insoluble purple-coloured formazan in living cells. We selected Vero cell lines as several authors have reported such cell safety study using Vero cell lines (Sukheja et al. 2017;Logu et al. 2002

MD simulations studies
Molecular docking studies and MM-GBSA analysis also showed that compound B12 might form a stable complex with DprE1 enzyme, but crucial H-bond interaction with LYS418 residue was missing. Hence, to better understand the interaction of compound B12 with DprE1 protein, MD simulations were carried out. The trajectory file obtained after 20,000 ps simulation was analysed, the simulation interaction diagram report was generated, and the detailed result is shown in Fig. 6. One frame was captured every 20 ps, and thus, a total of 1000 frames were generated for the study which was superimposed on the first frame, and RMSD value was calculated as shown in Fig. 6a. DprE1-B12 RMSD plot suggests that the complex was stable throughout the simulation period as the deviation was less than 2 Å at any given point of time. However, a slight drift was observed for the first half of simulation, i.e. first 10 ns; after that, the complex stabilized and drifted together. The most interesting finding was observed in the 2D protein-ligand contact plot as shown in Fig. 6b. We observed that B12 showed crucial H-bond interaction with LYS418 residue for 75% of the simulations time. This interaction was not observed in the docking studies. Also, it was evident from the histogram plot ( Fig. 6c) that B12 showed similar non-bonding interactions with DprE1 enzyme as shown by its co-crystallized ligand. It showed hydrophobic interactions with LYS60, TRP230, LUE283, TYR314, PRO316, LEU317, PHE362, LEU363, VAL385 and CYS387 residues. A water bridge interaction was also observed with TYR314 residue. Thus, MD simulation analysis supports our hypothesis that B12 will inhibit DprE1 protein.

Conclusion
In the present work, we have developed a four-feature e-Pharmacophore model based on the receptor-ligand cavity of DprE1 protein and mapped our previous reported library of compounds against it. The compounds were ranked on phase screen score, and the insights obtained from their alignment was used to design some novel benzimidazole-substituted compounds which might act as DprE1 inhibitors. These designed compounds had three matches out of four pharmacophoric features and favourable phase screen score up to 1.84. The designed compounds were also put for docking studies in XP mode, and some derivatives like B1, B2, B4, B5 and B12 showed comparable docking score with respect to the co-crystallized ligand. The designed compounds were synthesized in good yield, purified and characterized using various techniques like IR, NMR, LC-MS and elemental analysis. In vitro antitubercular activity was carried out on M. tuberculosis H37Rv (ATCC27294) strain using agar dilution method, and minimum inhibitory concentration (MIC) was determined in triplicates. All the compounds showed potent antitubercular activity in comparison with ethambutol. The compound B12 showed the most potent antitubercular activity with a MIC of 1.56 μg/ml which was better than the standard drug ethambutol (3.125 μg/ml). Compounds B7 and B11 were found to be equipotent with ethambutol and had a MIC value of 3.125 μg/ml. From this study, we can conclude that the piperidine ring substituted with furylidene rings and methyl benzimidazole ring showed potent antitubercular activity. Further substitution at the 4th position of piperidine ring was well tolerated with electron-withdrawing groups at the phenyl ring on o,p-position that led to an increase in the activity. Cytotoxicity studies against Vero cell lines proved that these compounds can be considered to be safe which is reflected in their high SI values. MD simulation study was done for the most potent compound B12 in complex with DprE1 protein and it also suggests that B12 will form a stable complex with DprE1 protein and showed the crucial H-bond interaction with LYS418 residue. To the best of our knowledge, this is the first work reported where e-Pharmacophore modelling has resulted in