Abstract
Background
Identification of novel drug targets and their inhibitors is a major challenge in the field of drug designing and development. Diaminopimelic acid (DAP) pathway is a unique lysine biosynthetic pathway present in bacteria, however absent in mammals. This pathway is vital for bacteria due to its critical role in cell wall biosynthesis. One of the essential enzymes of this pathway is dihydrodipicolinate synthase (DHDPS), considered to be crucial for the bacterial survival. In view of its importance, the development and prediction of potent inhibitors against DHDPS may be valuable to design effective drugs against bacteria, in general.
Results
This paper describes a methodology for predicting novel/potent inhibitors against DHDPS. Here, quantitative structure activity relationship (QSAR) models were trained and tested on experimentally verified 23 enzyme's inhibitors having inhibitory value (K_{i}) in the range of 0.00522(mM). These inhibitors were docked at the active site of DHDPS (1YXD) using AutoDock software, which resulted in 11 energybased descriptors. For QSAR modeling, Multiple Linear Regression (MLR) model was engendered using best four energybased descriptors yielding correlation values R/q^{2 }of 0.82/0.67 and MAE of 2.43. Additionally, Support Vector Machine (SVM) based model was developed with three crucial descriptors selected using Fstepping removeone approach, which enhanced the performance by attaining R/q^{2 }values of 0.93/0.80 and MAE of 1.89. To validate the performance of QSAR models, external crossvalidation procedure was adopted which accomplished high training/testing correlation values (q^{2}/r^{2}) in the range of 0.780.83/0.930.95.
Conclusions
Our results suggests that ligandreceptor binding interactions for DHDPS employing QSAR modeling seems to be a promising approach for prediction of antibacterial agents. To serve the experimentalist to develop novel/potent inhibitors, a webserver "KiDoQ" has been developed http://crdd.osdd.net/raghava/kidoq webcite, which allows the prediction of K_{i }value of a new ligand molecule against DHDPS.
Background
An escalating magnitude of drug resistance among bacterial pathogens has been installing a serious threat on the public health and economy of the developed world. A survey report has suggested that the direct cost to US economy alone due to drug resistant bacterial infection is around $4$5 billion annually [13]. Even for pharmaceuticals companies, it turns out to be a heartdying situation that after investing ~$800 million and about 15 years of atrocious labor to introduce a drug in the market, the pathogens already attains resistance against the drug. Therefore, there is an urgent need to recognize new inhibitors against novel and/or known targets. Undoubtedly, wellestablished bacterial targets i.e. cell wall and membrane biosynthesis, protein biosynthesis, nucleic acid etc always the first choice for developing antibacterials. The recent trend in this direction indicates that researchers are looking for novel targets alongside to discover new classes of inhibitors/antibiotics.
The amino acids biosynthetic pathways specifically lysine pathway has gained special attention because of its potential role in bacterial cell wall and protein synthesis [4,5]. The D, Ldiaminopimelic acid (mesoDAP), an important intermediate in the biosynthetic pathway of lysine is crucial in crosslinking peptidoglycan chains to provide strength and rigidity to the bacterial cell wall (known as DAP pathway). The absence of this pathway in mammalian system suggests that specific inhibitors of this biosynthetic pathway may be a valuable for developing novel classes of antibacterial agents. In this study, we explored DHDPS enzyme of the pathway, which catalysis condensation of pyruvate and aspartate semialdehyde to form DHDP. Figure 1 shows the established DAP pathway for DAP and lysine biosynthesis. The enzyme is encoded by dapA gene, which has been cloned and expressed from several strains, including Thermatoga maritima, Corynebacterium glutamicum, Mycobacterium tuberculosis and Bacillus anthracis. The threedimensional structures of DHDPS enzyme from Escherichia coli, Staphylococcus aureus, M. tuberculosis and B. anthracis enzymes with substrate pyruvate and without have been reported [618].
Figure 1. Enzymatic action of DHDPS leads to the biosynthesis of bacterial cell wall and protein components. Figure 1 shows the action of DHDPS enzyme involved in protein and cell wall synthesis process.
The antibacterial identification using experimental techniques is invariably very expensive, requires extensive pains and labor. Therefore, in silico techniques, which have the power to cut down these unavoidable steps, would be valuable. In recent years, in silico techniques like quantitative structure activity relationship (QSAR) and molecular docking are gaining high popularity in the drug discovery [1921]. Both these methodologies allow the identification of probable lead candidates expeditiously prior to chemical synthesis and characterization, thereby, making the process more cost effective [22,23].
In the present study, we attempt to integrate power of two in silico potential techniques: QSAR and molecular docking by using docking generated energybased descriptors for building QSAR models. Using this strategy, the information regarding binding mode of ligands in the active site is accumulated which would in turn assist the accurate prediction of better inhibitor with improved K_{i }values. To facilitate this we also developed a webinterface to help experimentalist working in the field of designing novel inhibitors against DHDPS enzyme.
Results
For the docking of 23 inhibitors, E. coli DHDPS crystal structure stored in the PDB file 1YXD was retrieved. The crystal structure of DHDPS consisted of two similar chains (A and B) with inhibitor bound at allosteric site [13]. The water molecules and inhibitor were removed using PYMOL software and chain A was considered for the docking purpose. The python scripts were used for carrying out automated flexible docking of 23 inhibitors on the predefined and experimentally characterized binding pocket, where the residue LYS161 being particularly very important. Hence, it's important to consider the flexibility of LYS161 and the inhibitors, while performing docking. Figure 2 shows the docking of two inhibitors: Inh6 (having minimum K_{i }value) and Inh10 (with maximum K_{i }value) at the active site of DHDPS enzyme. In order to validate our docking methodology another crystal structure of E. coli DHDPS (3DU0) with substrate bound at the active site was obtained from PDB. The enzyme 1YXD could not be used as it enclosed bound conformation of an allosteric inhibitor (S)Lysine. Since crystals were remarkably similar (RMSD value of 0.15Å), therefore, the same procedure for the docking of pyruvate was adopted which resulted in very slight variation in the RMSD value of 0.31 Å. Hence, the docking protocol adopted in the present study was able to reproduce the conformation comparable to the crystal structure with substrate at active site. Additionally, analysis of 10 docked poses of substrate generated by AutoDock software was also carried out. In Additional File 1: table S1 we have shown the values of free binding energies and RMSDs in the inreasing order of ranking. It was observed that RMSD value for the fifth ranked pose was lesser in comparison with the pose with best and minimum free binding energy. We also calculated the pairwise corelation between free binding energy and RMSD, resulted in R value of 0.81, which reveals that there exists correlation between free binding energy and RMSD values, however not the ideal or perfect one. Therefore, it's not always true that the pose with the lowest binding energy is the one with the lowest RMSD to the crystal structure. Ofcourse, one can validate or check the RMSD values for a single ligand system with bound crystal structure known. However, during virtual screening procedure with large number of unknown structures to dock, it's practically impossible to obtain the RMSD values. Therefore, in such cases, it has been shown in the past that the compounds with the lowest binding energies are generally considered as potential hits.
Figure 2. View of docked Inh10 (A1 and A2) and Inh6 (B1 and B2) at the active binding site of DHDPS. Figure 2 shows the docked conformation of inh10 and inh6 in active site of DHDPS where protein is shown in secondary structure and inh10, inh6 is represented in ball and stick model.
Additional file 1. Textbox S1: Selection of descriptors based on the chemical structures and activities. Additional file shows the descriptor selection on the basis of similarity in chemical structure of inhibitors and their activity.
Format: DOC Size: 43KB Download file
This file can be viewed with: Microsoft Word Viewer
It's important to mention that in general, after docking, AutoDock computes 11 types of energy values i.e.  i) Estimated free energy of binding (E_{FreeBind}); ii) Final Intermolecular Energy (E_{InterMol}), which is the sum of 4 energies such as (iii) vdW + Hbond + desolv Energy (E_{VHD}), (iv) Electrostatic Energy (E_{Elec}), (v) Moving LigandFixed Receptor (E_{MLFR}), and (vi) Moving LigandMoving Receptor (E_{MLMR}); vii) Final Total Internal Energy (E_{FTot}), again the sum of 2 energy values such as (viii) Internal Energy Ligand (E_{IntL}), and (ix) Internal Energy Receptor (E_{IntR}); (x) Torsional Free Energy (E_{Tors}) and (xi) Unbound System's Energy (E_{Unb}). Finally, 11 types of energy values based descriptors were then used as independent variables for QSAR modeling.
To obtain significant and noncorrelated variables from the abovementioned 11 descriptors, a statistical package, STATISTICA, was used. Indeed all the descriptors were highly significant, showing the p < 0.05. To filter out correlated descriptors, pairwise correlation coefficient at the cutoff value of 0.9 was imposed. The two variables namely E_{Unb }and E_{IntL }yielded the pairwise correlation values > 0.9 with E_{FTot }and therefore filtered out from the further analysis (Table 1).
Table 1. Matrix showing the pairwise correlation values for docking generated 11 energybased descriptors
Using MLR, a QSAR based model was generated using 4 variables namely, E_{FreeBind}, E_{Elec}, E_{IntR}, and E_{Tors }which accomplished correlation (R/q^{2}) values of 0.81/0.65 with MAE of 2.61 (Table 2). Though model was able to obtain good correlation values however, q^{2 }value was observed to be very low (Figure 3). Next, SVM along with Fstepping variable selection approach was employed. During first cycle of Fstepping removeone, an elimination of fifth descriptor i.e. E_{MLFR }from the set of n = 9 and the development of SVM model using 8 remaining variables attained the best correlation R/q^{2 }values of 0.87/0.75 and MAE of 2.24 (C = 50, γ = 45) listed in Table 2. For the second cycle, the removal of 9^{th }descriptor i.e. E_{IntR }further improved the correlation value to 0.91/0.81 showing reduction in MAE to 2.01 (C = 50, γ = 50). The next cycle however did not enhance the correlation values significantly as exclusion of E_{MLMR }offered correlation values 0.90/0.81 and MAE of 2.16 (C = 75, γ = 50) therefore, making its absence or presence to elicit no influence on correlation values. The correlation between predicted and actual activity values is shown in Figure 4. In the subsequent cycles of variables selection, no improvement in correlation values was observed. Therefore, it can be deduced that 6 descriptors i.e. E_{FreeBind}, E_{InterMol}, E_{VHD}, E_{Elec}, E_{FTot }and E_{Tors }are important to predict the inhibitory activity values for the present dataset of 23 inhibitors against DHDPS. Further, an external crossvalidation was carried out by randomly dividing 23 inhibitors dataset into three different sizes of training and test sets such as 21 and 2; 19 and 4; 17 and 6 respectively. The highest correlation q^{2}/r^{2 }values of 0.81/0.97 (an average of 89 best models) was obtained for the largest training and smallest test sets of 21 and 2 inhibitors (Table 3). An increase in the size of test set with corresponding decrease in training set size reduced the r^{2 }values along with slight reduction in q^{2 }values. The notion behind this splitting was to appraise a high predictive correlation values on the test set even when the size of training set was very low.
Table 2. Correlation values for MLR and SVM based QSAR models developed using descriptors selected at pairwise correlation cutoff value 0.9
Table 3. Detailed results obtained during external crossvalidation procedure using six descriptors with pairwise correlation cutoff value below 0.9
Figure 3. Comparison between actual and predicted Ki values for MLR model generated using descriptors selected at pairwise correlation cutoff value 0.9. Figure 3 depict the experimental and predicted Ki value in X and Y direction respectively with q^{2 }value 0.6531 using MLR model.
Figure 4. The correlation between actual and predicted Ki values for SVM model generated using variables selected at pairwise correlation cutoff value 0.9. Figure 4 illustrate the experimental and predicted Ki value in X and Y direction respectively with q^{2 }value 0.811 using SVM model.
Besides, pairwise correlation coefficient values listed in Table 4 between K_{i }and energybased descriptors for 23 inhibitors were calculated. Surprisingly, three descriptors such as E_{FreeBind}, E_{InterMol}, and E_{VHD }from the finally selected 6 energybased descriptors (described earlier) showed high fluctuations with respect to K_{i }values (Figure 5), which in turn have higher pairwise correlation coefficient values (irrespective of signs). On the other hand, the variables i.e. E_{FTot }and E_{Tors }observed to be neutral, indeed provided low correlation values of 0.20 and 0.075 respectively. The E_{MLFR}, which showed high correlation value of 0.62 with K_{i}, was removed in the first cycle of variables selection. These deviations prompted us to carry out the clustering of the dataset of 23 inhibitors using JChem software http://www.chemaxon.com/ webcite. As shown in Figure 6, all 23 inhibitors clustered into two unique groups however, three compounds Inh8, Inh15, Inh17 were variable. Hence, to filter the noise that might be caused due to these 3 inhibitors, a QSAR model was again generated removing these 3 structures. Using Fstepping variable selection approach, 5 energybased input variables i.e. E_{FreeBind}, E_{InterMol}, E_{VHD}, E_{IntL }and E_{MLFR }generated a QSAR model attaining correlation R/q^{2 }values of 0.95/0.89 and MAE of 1.28 (Figure 7 and Table 2).
Table 4. Pairwise correlation values for 11 energybased descriptors with respect to Ki values
Figure 5. The variations in the values of 6 energybased input variables with respect to experimental Ki values. Figure 5 shows the variation in energy descriptors with respect to Ki values.
Figure 6. The clustering of 23 inhibitors dataset. Figure 6 illustrate the clustering graph of all 23 inhibitors where inh8, inh15 and inh17 exist as singleton.
Figure 7. Comparison between actual and predicted Ki values for SVM model developed using 20 inhibitors dataset. Figure 7 depict the experimental and predicted Ki value in X and Y direction respectively with q^{2 }value 0.89.
Additionally, QSAR modeling was also carried out using six noncorrelated descriptors i.e. E_{FreeBind}, E_{Elec}, E_{Ftot}, E_{IntR}, E_{MLMR }and E_{Tors }having pairwise correlation value less than 0.5. Using MLR, the QSAR model was developed with four types of input variables E_{FreeBind}, E_{MLMR}, E_{IntR}, and E_{Tors}, which accomplished correlation (R/q^{2}) values of 0.82/0.67 and MAE of 2.43 (Table 5). Hence, a small increase in correlation values in comparison to earlier MLR model (0.81/0.65 with MAE of 2.61), was observed. Further, SVM model was trained using six descriptors but the model attained poor correlation (R/q^{2}) values of 0.67/0.40 and MAE of 2.91 (C = 200, γ = 1) signifying the presence of some descriptors ideally may not required for robust model generation. Thus, employing the Fstepping variable selection approach, the removal E_{MLMR }from the set of n = 6 energybased descriptors and using the 5 remaining descriptors for the SVM model development achieved best correlation R/q^{2 }values of 0.83/0.67 and MAE of 2.63 (C = 20, γ = 55) (Table 5). Next, filtering of E_{Tors }enhanced the correlation (R/q^{2}) value of model to 0.84/0.69 with reduction in MAE to 2.51 (C = 75, γ = 125). Then, QSAR modeling was carried out on the remaining four descriptors and the removal of E_{FTot }augmented the correlation value to 0.93/0.80 with attenuation of MAE value to 1.89, a noteworthy enhancement. The correlation between predicted activities for the inhibitors and their actual experimental values is depicted in Figure 8A. There exists a good agreement between predicted and experimental activity values hence suggesting the robustness of QSAR model. Therefore, the three energybased descriptors such E_{FreeBind}, E_{Elec}, and E_{IntR }were imperative to predict the inhibitory activity values. Interestingly, the performance of three descriptor based QSAR model was found to be better in comparison to the six descriptors based SVM model (0.90/0.81 and MAE of 2.16) which was described earlier using cutoff value of 0.9. Additionally, the removal of three outliers like Inh7, Inh16 and Inh23, further enhanced the prediction efficiency of QSAR model by increasing the correlation (R/q^{2}) values to 0.94/0.87 and reduction in MAE to 1.45 (Figure 8B). Further, high external training/testing crossvalidated correlation values (q^{2}/r^{2}) in the range of 0.780.83/0.930.95 was attained by randomly splitting the dataset into several training sets for model building and independent testing on corresponding test sets (Table 6).
Table 5. Correlation values for MLR and SVM based QSAR models developed using noncorrelated variables selected at the cutoff value 0.5.
Table 6. Results obtained during external crossvalidation procedure using three noncorrelated descriptors
Figure 8. Scatter plot between experimental versus predicted Ki values provided by highly noncorrelated 3 energy values based SVM model for 23 inhibitors (A) and 20 inhibitors dataset (B). Figure 8 depict the experimental and predicted Ki value in X and Y direction respectively with q^{2 }value 0.80 and 0.87 for 23 and 20 inhibitors dataset respectively.
In order to assess robustness and validation of the finally developed three descriptors based QSAR model, a bootstrap analysis for 100 runs by statistical sampling of the original dataset was also performed which yielded a higher q^{2}_{bootstap }value of 0.88 ± 0.029. Thus, higher and lower values of q^{2}_{bootstap }and standard deviation parameters comprehensively support the statistical validity of the presently developed QSAR models. Further, Yrandomization test was also carried out using shuffled activity dataset which resulted in poor performance i.e. nearly all of the q^{2 }values were < zero (q^{2 }ranged from 0.15 to 0.41 and MAE from 4.15 to 5.13), thereby signifying the consistency of QSAR model.
2D descriptors based QSAR modeling
To compare the performance of three energybased QSAR model with simple 2D QSAR models, the 2D QSAR modeling with 14 noncorrelated descriptors was also performed. The study was commenced by MLR modeling, which tried to establish structureactivity relationship using five descriptors i.e. MSD, PJI2, Jhetm, ALOGP2, and Me by attaining correlation R/q^{2 }values of 0.78/0.61 and MAE of 3.08 (Table 7). The performance of 2D descriptors based model was found to be lower in comparison with four energybased MLR model described earlier (R/q^{2 }values of 0.82/0.67 and MAE of 2.43). Further, using all 14 noncorrelated descriptors for the training of SVM model, a very poor correlation (R) value of 0.23 was observed. The removal of 7 descriptors: nBm, BLI, Jhetm, GATS1v, nHAcc, ALOGP and MATS3 m (after employing 7 cycles of Fstepping removeone) and the training of SVM model with remaining 7 descriptors attained R/q^{2 }values of 0.77/0.57 and MAE of 3.26 (C = 25, γ = 25). During the next cycle, an exclusion of JGI2 (topological charge index) optimized the SVM model (C = 300, γ = 25) by achieving correlation R/q^{2 }values of 0.79/0.60 and MAE of 3.10. In the next cycle, the removal of nH (constitutional) descriptor improved the correlation value to 0.82/0.64 with reduction in MAE value to 2.68 (C = 200, γ = 25). Finally, elimination of Me constitutional descriptor and the development of QSAR model on the remaining 4 descriptors, which included two topological descriptorsMSD, PJI2, molecular propertyALOGP2 and Burden eigenvalues descriptorBEHm1, yielded correlation R/q^{2 }values of 0.84/0.67 and MAE of 2.61 (C = 300, γ = 25). Hence, the performance of SVM based 2D QSAR model was found to be very low in comparison with three energy values based QSAR model developed using SVM.
Table 7. Performance of 2D QSAR based MLR and SVM models
Implementation of webserver
We attempted to develop efficient QSAR model and based on these models, a web server "KiDoQ" (available at http://crdd.osdd.net/raghava/kidoq webcite) using CGIPERL and python scripts was developed. User can draw the structure of ligand molecule using JME editor incorporated on the server. The server also accepts input as mol/mol2 structure files pasted or uploaded on the server (Figure 9). The working flow of KiDoQ server is shown in Figure 10.
Figure 9. A snapshot of submission page of KiDoQ webserver. Figure 9 shows the screen shot of KiDoQ webserver.
Figure 10. Schematic flow for the working of KiDoQ webserver. Figure 10 illustrate the workflow diagram of KiDoQ webserver.
Discussion
The QSAR modeling has been accepted as a promising methodology for lead identification. Nevertheless, if highresolution target structure is available, then receptor structure based approach is often a first choice. The recent studies have shown the better performance of QSAR models even in the presence of target structure. However, simple QSAR approach can sometimes lead to false prediction if the collected data does not cover the complete property space or the selected 2D/3D descriptors are not reliable. Therefore, both techniques have their own advantages and limitations [1923]. Keeping in view, the importance of docking and better performance of QSAR, we integrated both approaches by using docking generated energybased scores as descriptors for QSAR modeling. The major benefit presumed by this integration would be an additional validation of the docking predicted inhibitors as bioactive or inactive by prediction of their bioactivity values using QSAR models henceforth, would facilitate in reduction of false positives.
In the present study, docking of the 23 experimentally known inhibitors of DHDPS at its active binding site resulted in 11 energybased descriptors. For valid statistical results, it was imperative to restrict the maximal number of descriptors or to remove highly correlated ones, as presence of redundancy reduced the discriminating power of input variables, thereby reducing their worth in model development. Ideally, a regression model with n training set compounds and k descriptors may be acceptable only if n > 4k and for any of the k descriptors i) the significance level p is < 0.05; ii) the pairwise correlation coefficient should be < 0.9 [24,25]. Therefore, we also looked for statistically significant and noncorrelated energybased descriptors. All 11 types of energybased descriptors were statistically significant but only 2 and 5 variables showed pairwise correlation value > 0.9 and 0.5 respectively, resulting their removal necessary for rigorous QSAR modeling. As QSAR modeling for DHDPS enzyme's inhibitors is being carried out for the first time, therefore both linear (MLR) and nonlinear (SVM) techniques were employed. In our study, the performance of SVM model was found to be much better in comparison with MLR model. The simple linear model was unable to handle the diversity of the present dataset; therefore, the cases where simple linear techniques fail, nonlinear techniques could provide a better option.
Further, it was also noticed that the structural diversity of 23 inhibitors and redundancy among finally selected six descriptors resulted in wrong selection of input variables. Therefore, we removed highly diverse three structures and used remaining 20 inhibitors for QSAR modeling. Interestingly, the performance of model was found to be enhanced in comparison to the model developed using 23 inhibitors. Further, QSAR modeling carried out with highly noncorrelated 3 descriptors i.e E_{FreeBind}, E_{Elec}, and E_{IntR }(selected at the cutoff pairwise correlation value of 0.5) provided better correlation values in comparison to the earlier six descriptors (selected at the cutoff pairwise correlation value of 0.9) based QSAR model. Therefore, removal of redundant descriptors reduced the noise and enabled the better training of QSAR models. The three noncorrelated descriptors appeared to be governing factors in establishing structure actvity relationship for DHDPS enzyme. One of the possible reasons for their selection was a higher pairwise correlation value with respect to K_{i }in comparison with other descriptors i.e. E_{Ftot}, E_{MLMR }and E_{Tors }removed during QSAR modeling.
Among three descriptors, the value of E_{FreeBind }based descriptor was found to be dependent on other two descriptors i.e. E_{Elec}, and E_{IntR }as well on other correlated energy based descriptors (such as E_{VHD}). Generally, in the absence of receptor's flexibility E_{IntR }remains constant and does not make any significant contribution to E_{FreeBind}, however, the flexibility incorporates transformations leading to internal energy changes. In the present study, the changes in the value of E_{IntR }were observed as the LYS161 was kept flexible. It was noticed that inhibitors with lower K_{i }value were characterized by high negative E_{IntR }values. These inhibitors included aliphatic compounds generally the pyruvate and aspartate semialdehyde analogues. On the other hand, inhibitors such as Inh17, Inh14, Inh23 and Inh10 with higher K_{i }values exhibited lower negative E_{IntR}. In view of this, we suggest E_{IntR }based descriptor is an important discriminating variable for developing robust QSAR models. Further, E_{Elec }was also found to be imperative as variations in the E_{Elec }values was highly dependent on the number and type of receptor residues involved in establishing charge interactions with inhibitors. We observed that inhibitors such as Inh9, Inh11, Inh18 and Inh20 with strong electrostatic interactions with receptor exhibited strong binding, resulting in higher negative E_{FreeBind }and E_{Elec }values that in turn provided lower K_{i }values. A few inhibitors such as Inh17 and Inh14 characterized by aromaticity, the strong electrostatic or πcationic interactions though provided higher negative E_{Elec }values however, at the cost of reasonable reduction in the E_{IntR }values, which in turn provided higher K_{i }values, in comparison to the aliphatic inhibitors and the ones with weak electrostatic interactions. In addition, it was also figured out that inhibitors i.e. Inh1, Inh3, Inh4, Inh5, Inh6, Inh12, Inh13, Inh15 and Inh19 exhibited strong affinity to receptor albeit no or very weak electrostatic interactions were observed. This suggest that binding of inhibitors to DHDPS is not specifically dependent on electrostatic interactions, however other bonded and nonbonded interactions appeared to be playing important role, which in turn provided higher negative E_{FreeBind }values and lower K_{i }values (For details see Additional File 1).
As we have employed a complex procedure of using docking generated energybased descriptors for QSAR modeling; therefore, it became imperative to compare the model performance with simple conventional 2D QSAR models. The SVM based 2D QSAR model achieved a poor correlation value of 0.84/0.67 in comparison with docking energybased SVM model (0.93/0.80) indicating inadequacy of 2D descriptors in providing acceptable and robust QSAR model for dataset of 23 inhibitors. This low performance of 2D QSAR models may be due to presence of high structural diversity among the inhibitors that was not easily captured using simple 2D descriptors.
To conclude, the present strategy of predicting K_{i }values using docking generated energybased descriptors for QSAR modeling is a promising approach to predict potent inhibitors against DHDPS enzyme.
Conclusions
In this study, we describe a new approach for prediction of antibacterial compounds that both take QSAR and docking strategy into its consideration. By using this approach, we get promising results instead of using these two strategies individually and develop a webserver called KiDoQ. This webserver will be helpful for better prediction of antibacterial compounds against dihydrodipicolinate synthase (DHDPS).
Methods
Inhibitors Dataset
The information regarding the experimentally known 23 inhibitors, classified as potent, moderate and slightly weak, was obtained from the literature [618]. The IUPAC names of these inhibitors along with K_{i }values are shown in Table 8. Chem3D Ultra (v11.0), windowsbased software was used for sketching the 2D structures for all inhibitors followed by cleaning and refinement in order to correct the accidentally distorted or unrealistic bond angles and lengths. The 2D structures were converted into 3D structures using CORINA software. Then each structure was energy minimized to give energetically preferred 3D structures.
Table 8. Dataset of 23 inhibitors along with their experimentally known K_{i }values
Docking energybased descriptors
were calculated using automated docking software AutoDock (v.4.0) (AD) [26]. It is a suite of three C programs: i) AutoTors, which facilitates input of ligand coordinates; ii) AutoGrid, which precalculates a threedimensional grid based on macromolecular coordinates; and iii) AutoDock, which performs a actual docking simulations. Before docking process, several separate predocking steps: ligand preparation, receptor preparation and grid map calculations were performed. The ligand and receptor preparation stage involved the addition of hydrogen atoms, computing charges, merging nonpolar hydrogen atoms and defining AD4 atom types to ensure that atoms conformed to the AutoDock atom types. The information about rotatable torsion bonds that defines the bond flexibility was acquired. The ligands and receptor molecule preparation was followed by grid construction using AutoGrid module. During grid construction, atom types of the ligand, which acted as probes in the calculation of grid maps, were identified. The grid with default volume of 40 × 40 × 40 Å with a spacing of 0.375Å centered on the receptor was prepared. For conformational searches, the docking calculations using the genetic algorithm (GA) procedure with default parameters was performed.
2D QSAR modeling
DRAGON software was used for the calculation of 2D descriptors. For our dataset, the software calculated ~848 types of 2D descriptors categorized into different descriptor blocks such as constitutional descriptors, topological descriptors, walk and path count descriptor, connectivity indices, information indices, 2D autocorrelation, edge adjacency, burden eigenvalues, topological charge indices, functional groups, molecular properties and eigenvalues based indices. Initially, the descriptors with zero or unassigned values were excluded and then pairwise correlation test to remove highly correlated descriptors at a cutoff value of 0.50 was executed. This procedure resulted in 14 descriptors for 2D QSAR modeling.
QSAR Model Construction
QSAR methodology quantitatively correlates the structural molecular properties (descriptors) with functions (biological activities) for a set of compounds by means of linear or nonlinear statistical methods. In the present study, we exploited both linear (MLR) and nonlinear (SVM) statistical methods for flourishing the robust QSAR models [27,28]. Retrospectively, for QSAR modeling, both linear and nonlinear models have been extensively used [2936]. MLR tries to model the relationship between two or more independent descriptors and dependent variable such as K_{i }by fitting a linear regression equation to the observed data with corresponding parameters (constants) and an error term. On the other hand, SVM based on statistical and optimization theory, handles complex structural features. In the present study, SVM_light http://www.cs.cornell.edu/People/tj/svm_light/ webcite, which is an implementation of SVM, was used for QSAR modeling.
Evaluation of QSAR models
To assess the predictive performance of QSAR models, different crossvalidation procedures were adopted. First, in leaveoneout strategy (LOOCV), one molecule was removed from the dataset as a test compound and the remaining 22 molecules were used to build the model. This process was repeated 23 times with each inhibitor as a test molecule. Once a regression model was constructed, goodness about the fit and statistical significance was assessed using the statistical parameters outlined below
where, x_{i }and y_{i }represents the actual and predicted K_{i }values for the i^{th }compound. N is the total number of compounds, represents the averaged value of the actual K_{i }for the entire dataset.
Here, it was equally important to use an independent test set to check the real predictive accuracy of trained QSAR models. However, 23 compounds were not expected to be sufficient for independent testing using existing QSAR models. Therefore, an alternative strategy, external crossvalidation, was adopted, where different number of inhibitors i.e. 2, 4, and 6 were randomly selected as independent test sets. The models were then trained on the remaining inhibitors i.e. 21, 19, and 17 using LOOCV procedure followed by independent testing on the corresponding test sets. This cycle of randomly separating test and training sets was repeated. Here, to determine the predictive accuracy of models on the test set, predictive r^{2 }value was used
where, SD is the sum of the squared deviations between the activities of the test set and mean activities of the training molecules.
Then, Yrandomization test was performed in order to appraise high training and testing correlation values observed during QSAR modeling, were not occurred incidentally. Here, the shuffled activity dataset was derived by randomly shuffling the dependent variables K_{i }and keeping the descriptors original, afterward using this randomly shuffled dataset to develop new QSAR models. The process of shuffling was carried out many times with subsequent generation of corresponding models nevertheless, with an assumption that the resulting models should give low performance, which would obviously imply the rigorous robustness of the original models.
Input variables selection
The selection of best descriptors that establish the relationship between chemical structure and an inhibitory property is crucial for the success of QSAR modeling. In the present study, we adopted Fstepping removeone approach for variable selection. Accordingly, each input variable was removed onebyone from the set of n variables followed by QSAR modeling using the remaining n1 variables. However, if the correlation value was increased, the particular variable was permanently removed from the analysis. These cycles were repeated until no further improvement in the correlation values was observed and stopped if n1 removal resulted in reduction of correlation values.
List of abbreviations used
The abbreviations used are: QSAR: Quantitative Structural Activity Relationship; DAP: Diaminopimelic Acid; CADD: Computed Aided Drug Designing; LYS161: Lysine161; SVM: Support Vector Machine; LOOCV: LeaveOneOut CrossValidation; MAE: Mean Absolute Error; MLR: Multiple Linear Regression.
Competing interests
The authors declare that they have no competing interests.
Authors' contributions
AG, GPSR and RT conceived and designed the experiments. AG performed the experiments, wrote perl scripts, developed server. GPSR and AG analyzed the data. AG wrote the manuscript. AG and GPSR carried out revision of the manuscript. This manuscript has been seen and approved by all authors.
Availability and requirements
Project name: A webserver for predicting DHDPS inhibitors
Project home page: http://crdd.osdd.net/raghava/kidoq webcite
Operating system(s): Platform independent;
Programming language: PERL, CGIPERL;
License: None;
Any restrictions to use by nonacademics: No restrictions
Acknowledgements
AG is thankful to CSIR for providing SRF. AG is also thankful to Manish Datt and Nitish Kumar for providing help in running AutoDock software. The authors are thankful to anonymous reviewers for their excellent suggestions and Dr. Kishore for critically editing our manuscript. The authors are also thankful to the Council of Scientific and Industrial Research (CSIR) and Department of Biotechnology, Government of India for financial assistance.
References

Cohen ML: Changing patterns of infectious disease.
Nature 2000, 406:762767. PubMed Abstract  Publisher Full Text

Livermore DM: Bacterial resistance: origins, epidemiology and impact.
Clin Infect Dis 2003, 36:S1123. PubMed Abstract  Publisher Full Text

McGowan JE Jr: Economic impact of antimicrobial resistance.
Emerg Infect Dis 2001, 7:28692. PubMed Abstract  Publisher Full Text  PubMed Central Full Text

Hutton CA, Perugini MA, Gerrard JA: Inhibition of lysine biosynthesis: an evolving antibiotic strategy.
Mol BioSyst 2007, 3:458465. PubMed Abstract  Publisher Full Text

Cox RJ, Sutherland A, Vederas JC: Bacterial diaminopimelate metabolism as a target for antibiotic design.
Bioorg Med Chem 2000, 8:843871. PubMed Abstract  Publisher Full Text

Burgess BR, Dobson RC, Dogovski C, Jameson GB, Parker MW, Perugini MA: Purification, crystallization and preliminary Xray diffraction studies to nearatomic resolution of dihydrodipicolinate synthase from methicillinresistant Staphylococcus aureus.
Acta Crystallogr Sect F Struct Biol Cryst Commun 2008, 64:65961. PubMed Abstract  Publisher Full Text

Pearce FG, Perugini MA, McKerchar HJ, Gerrard JA: Dihydrodipicolinate synthase from Thermotoga maritima.
Biochem J 2006, 400:35966. PubMed Abstract  Publisher Full Text  PubMed Central Full Text

Cremer J, Eggeling L, Sahm H: Cloning the dapA dapB cluster of the lysinesecreting bacterium Corynebacterium glutamicum.
Mol Gen Genet 1990, 229:478480. Publisher Full Text

Mirwaldt C, Korndorfer I, Huber R: The crystal structure of dihydrodipicolinate synthase from Escherichia coli at 2.5 A^{0 }resolution.
J Mol Biol 1995, 246:227239. PubMed Abstract  Publisher Full Text

Laber B, GomisRüth FX, Romão MJ, Huber R: Escherichia coli dihydrodipicolinate synthase. Identification of the active site and crystallization.
Biochem J 1992, 288:691695. PubMed Abstract  PubMed Central Full Text

Kefala G, Evans GL, Griffin MD, Devenish SR, Pearce FG, Perugini MA, Gerrard JA, Weiss MS, Dobson RC: Crystal structure and kinetic study of dihydrodipicolinate synthase from Mycobacterium tuberculosis.
Biochem J 2008, 411:35160. PubMed Abstract  Publisher Full Text

Blickling S, Renner C, Laber B, Pohlenz H, Holak T, Huber R: Reaction mechanism of Escherichia coli dihydrodipicolinate synthase investigated by Xray crystallography and NMR spectroscopy.
Biochem J 1997, 36:2433. Publisher Full Text

Dobson RCJ, Griffin MDW, Jameson GB, Gerrard JA: The crystal structures of native and (S)lysinebound dihydrodipicolinate synthase from Escherichia coli with improved resolution show new features of biological significance.

Blagova E, Levdikov V, Milioti N, Fogg MJ, Kalliomaa AK, Brannigan JA, Wilson KS, Wilkinson AJ: Crystal structure of dihydrodipicolinate synthase (BA3935) from Bacillus anthracis at 1.94 A^{0 }resolution.
Proteins 2006, 62:297301. PubMed Abstract  Publisher Full Text

Karsten WE: Dihydrodipicolinate Synthase from Escherichia coli: pH dependent changes in the kinetic mechanism and kinetic mechanism of allosteric inhibition by LLysine.
Biochem J 1997, 36:17301739. Publisher Full Text

Turner JJ, Gerrard JA, Hutton CA: Heterocyclic inhibitors of dihydrodipicolinate synthase are not competitive.
Bioorg Med Chem 2005, 13:213340. PubMed Abstract  Publisher Full Text

Turner JJ, Healy JP, Dobson RC, Gerrard JA, Hutton CA: Two new irreversible inhibitors of dihydrodipicolinate synthase: diethyl (E, E)4oxo2,5heptadienedioate and diethyl (E)4oxo2heptenedioate.
Bioorg Med Chem Lett 2005, 15:9958. PubMed Abstract  Publisher Full Text

Turner JJ, Gerrard JA, Hutton CA: Conformationally constrained diketopimelic acid analogues as inhibitors of dihydrodipicolinate synthase.
Bioorg Med Chem 2005, 13:213340. PubMed Abstract  Publisher Full Text

Schneider G, Fechner U: Computerbased de novo design of druglike molecules.
Nat Rev Drug Discov 2005, 8:649663. Publisher Full Text

Buolamwini JK, Assefa H: CoMFA and CoMSIA 3D QSAR and docking studies on conformationallyrestrained cinnamoyl HIV1 integrase inhibitors: Exploration of a binding bode at the active Site.
J Med Chem 2002, 45:841852. PubMed Abstract  Publisher Full Text

Lyne PD: Structurebased virtual screening: an overview.
Drug Discov Today 2002, 7:10471055. PubMed Abstract  Publisher Full Text

Muegge I, Oloff S: Advances in virtual screening.
Drug Discov Today 2006, 3:405411. Publisher Full Text

Willet P: Similaritybased virtual screening using 2D fingerprints.
Drug Discov Today 2006, 11:10461053. PubMed Abstract  Publisher Full Text

Golbraikh A, Tropsha A: Beware of q2!
J Mol Graph Model 2002, 20:269276. PubMed Abstract  Publisher Full Text

Tropsha A, Gramatica P, Gombar V: The importance of being earnest: validation is the absolute essential for successful application and interpretation of QSPR models.

Morris GM, Goodsell DS, Halliday RS, Huey R, Hart WE, Belew RK, Olson AJ: Automated docking using a Lamarckian Genetic algorithm and empirical binding free energy function.
J Computational Chemistry 1998, 19:16391662. Publisher Full Text

Jurs PC: Computer Software Applications in Chemistry. 2nd edition. John Wiley and Sons, New York; 1996.

Joachims T: Making largescale SVM learning particle. MIT Press, Cambridge, MA and London; 1999.

Warmuth MK, Liao J, Ratsch G, Mathieson M, Putta S, Lemmen C: Active Learning with Support Vector Machines in the drug discovery process.

Oloff S, Mailman RB, Tropsha A: Application of validated QSAR models of D1 dopaminergic antagonists for database mining.
Med Chem 2005, 48:73227332. Publisher Full Text

Niu B, Lu WC, Yang SS, Cai YD, Li GZ: Support vector machine for SAR/QSAR of phenethylamines.
Acta Pharmacol Sin 2007, 28:107586. PubMed Abstract  Publisher Full Text

Han LY, Ma XH, Lin HH, Jia J, Zhu F, Xue Y, Li ZR, Cao ZW, Ji ZL, Chen YZ: A support vector machines approach for virtual screening of active compounds of single and multiple mechanisms from large libraries at an improved hitrate and enrichment factor.
J Mol Graph Model 2008, 26:127686. PubMed Abstract  Publisher Full Text

Ashek A, Cho SJ: A combined approach of docking and 3D QSAR study of betaketoacylacyl carrier protein synthase III (FabH) inhibitors.
Bioorg Med Chem 2006, 14:147482. PubMed Abstract  Publisher Full Text

Kim HJ, Chae CH, Yi KY, Park KL, Yoo SE: Computational studies of COX2 inhibitors: 3DQSAR and docking.
Bioorg Med Chem 2004, 12:162941. PubMed Abstract  Publisher Full Text

Doble M, Karthikeyan S, Padmaswar PA, Akamanchi KG: QSAR studies of paeonol analogues for inhibition of platelet aggregation.
Bioorg Med Chem 2005, 13:59966001. PubMed Abstract  Publisher Full Text

Daszykowski M, Stanimirova I, Walczak B, Daeyaert F, de Jonge MR, Heeres J, Koymans LM, Lewi PJ, Vinkers HM, Janssen PA, Massart DL: Improving QSAR models for the biological activity of HIV Reverse Transcriptase inhibitors: Aspects of outlier detection and uninformative variable elemination.
Talanta 2005, 68:5460. PubMed Abstract  Publisher Full Text