QSAR, Pharmacophore Mapping and Molecular Docking of 2,4-Diaminoquinazoline as Antitubercular Scaffold: A Computational Hybrid Approach

*Corresponding Author:

S. K. Kashaw
Department of Pharmaceutical Sciences, Dr. Harisingh Gour University (A Central University), Sagar-470 003, India
E-mail: sushilkashaw@gmail.com

This is an open access article distributed under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike 3.0 License, which allows others to remix, tweak, and build upon the work non-commercially, as long as the author is credited and the new creations are licensed under the identical terms

Abstract

The 2,4-diaminoquinazoline is a well-known scaffold with a great potential to generate a library. In the following research work different computational strategies were applied on 2,4-diaminoquinazoline moiety to estimate its efficiency as an antitubercular scaffold among its other versatile biological activities reported in several published studies. Each data has defined specific substitution mode on the scaffold to be active on different protein site or for disease condition. The following experiment regards of 2D-QSAR, 3D-QSAR, active site and blind docking, structural orientation, pharmacophore mapping and further designing of the data set with possible real active moiety. 2D-QSAR has shown good reliability with r² = 0.8190, q²_LOO = 0.7711 and external pred_r² = 0.5321, along with 3D-QSAR has good predictability, q² (r²cv) = 0.7601, pred_r² = 0.5567. Further ligand based pharmacophore mapping was carried out for estimating atomic contribution to chemical feature for the compound. The generated hypothesis established that hydrogen accepter, donor and aromatic ring with electronegative atom are the important features. Finally, outcomes of all results were recapitulated to design new compounds. Some compounds also were designed depending on the computational finding and concludes the suitability of the scaffold to be antiTB active.

Keywords

Computational strategy, QSAR, docking, pharmacophore mapping, molecular designing, antiTB activity, validation

With the increase in the statistics of the tuberculosis (TB) affected population and the multiplicity of the types and symptoms of TB each year, a scary picture has emerged in the worldwide scenario. In 2016, 4.1 % of new cases and 19 % cases of multidrug resistant or rifampicin-resistant TB have been listed in the report[1,2]. After the resistance to isoniazid, resistance to other first line drugs is also emerging with time reducing treatment and maintenance options for TB patients. In this context, a new versatile molecular structure is the requirement for generating new leads and finally for a drug molecule development. Among different well known scaffolds, 2,4-diaminoquinazoline has been an attractive system and reported to be effective against Mycobacterium tuberculosis in a promising way[3]. The diaminoquinazoline pharmacophore is a versatile moiety, as smallest changes in the functional subgroup surprisingly revealed new novel potential compounds on a number of protein family. The moiety with different group substitution have drawn much attention due to their wide range of biological activities for example, antimalarial[4,5], antileshmanial[6], antipsychotic[7], antitubercular, anticancer drugs[8,9]. The concept of being the scaffold to be antiTB template was first found in 2000 over a public data service system[10].

In 2008 in an US patent, it was first described that the modification over the quinazoline ring system specifically at 2,4 position on the quinazoline ring as well as the benzene ring can be exploited to develop potent antiTB drug candidate[11,12]. With that goal more molecules were reported in 2015 by Odingo et al.[3] and declared the scaffold to be a promising template which should be evaluated more. In this context, computational analysis is a logical approach for the development and discovery of drug molecules in more defined way. So, it has great potentials to generate drug like lead molecules, which are constructive and futuristic approach to draw a well evaluated results in terms of mathematical expression.

Materials and Methods

QSAR study using V-life MDS:

A set of 33, 2,4-diaminoquinazoline derivatives evaluated in solid medium assay for antiTB activity reported by Odingo et al.[3] was subjected to the QSAR analysis. All QSAR studies were performed in V-Life MDS software Version 4.4. The inhibitory activity data were reported as MIC₅₀ and converted to pMIC₅₀ by taking negative log of MIC₅₀ to reduce the skewness of data and activity range was taken in log 2.5 fold scatter. Table 1 shows the structure of the compounds along with their biological activity values.

Table 1: Biological Activity Data and Substituent Present In the Compounds at Different Positions

For 2D QSAR, all the compounds were subjected to energy minimization to get 3D structures, using Merck molecular force field (MMFF)[13], followed by considering distance-dependent dielectric constant of 1.0 and convergence criterion of 0.001 kcal/mol. The QSAR work sheet was generated using biological activity as dependent variable and various 2D descriptors as independent variables. The data set was divided into training and test sets were generated using 70 % search method followed by partial least squares regression coupled with stepwise forward-backward variable selection method (SW-PLS) was incorporate for estimation.

In 3D QSAR, conformers were generated by Monte Carlo conformational search method and the conformers of least energy were selected for the alignment[14]. All the compounds were aligned by template-based method, where a template was built by considering common substructures of the series. The lowest energy conformer of the most active compound was selected as a reference[15]. In the present study, all the compounds were aligned against minimum energy conformation and most active compound 10 (RMSE value 0.000110) using quinazoline nucleus as template as shown in figure 1a and the alignment of molecules is shown in figure 1b. The aligned molecules were used to generate QSAR work sheet using biological activity as dependent variable and 3D descriptors as independent variables. The data set was divided into training set and test set by considering the biological activity method along with 70 % data selection way of k-nearest neighbourmolecular field analysis coupled with stepwise forwardbackward variable selection method (known as SW kNN-MFA method) was incorporate for estimation.

The pharmacophore identification studies were carried out in Mol sign module of V-life MDS 4.4 to assign the geometrical representation of the features necessary for the molecule to be active[16]. A pharmacophore model is a set of three-dimensional features that are necessary for bioactive ligands. Thus, the series of antiTB molecular set was aligned as total set and pharmacophoric was developed which contains five pharmacophoric features keeping the tolerance distance at 10 Å with maximum limit of 10.

Docking study using GOLD:

The set of 33 molecules were docked with 4 different protein molecules reported to be important for antiTB activity and the crystalline protein structures were taken from the RSCB PDB site, in order to identify the enzyme for a potent inhibitor. The chemical nature of the macromolecule and ligands greatly influence the performance of docking routines. An evaluation criterion was based on the docking scores and interacting atoms. The docking study was carried out in GOLD software version 5.2[17,18] and the ligands were prepared using Hyper Chem 7 professional. The ligand set was first subjected to MM+, then PM3 and finally through abintio quantum mechanics and the energy diagram was prepared of the whole molecular set shown in figure 2. In GOLD ChemPLP scoring method was opt with early termination off, search efficiency to be 100 % and GA run 10 000.

Designing of a new molecule:

Depending on the results of the above mentioned studies as well as all considerations of the structural constituents of the molecular set, new molecular set was designed. The generated new molecular set first subjected to optimization and then conformers were generated by Monte Carlo conformational search method. Finally, conformers of least energy were selected for the alignment. All the compounds were aligned by template-based method with respect to the template moiety and evaluated through the 3D QSAR model prediction which estimated probable activity of the system.

Result and Discussion

In this statistical assessment, SW-PLS method was incorporated for searching the best 2D QSAR model with no randomness and biasness. The variation of the percentage selection is depicted in the Table 2. In order to construct and validate the QSAR models, both internally and externally, the data sets were divided into training (65, 70 and 75 %) and test sets (35, 30 and 25 %). Twenty nine trials were run leading to a training set with 75 % better r², q²_{_LOO}, pred_r² values and F test values than the others one. The selection criteria for training and test set in a way that biological activities of all compounds in the test set lie within the maximum and minimum value of biological activities of the training set compounds. 2D QSAR Eqns were selected by optimizing the statistical results developed with the variation of descriptors in these models. The frequency of use of a particular descriptor in the population of Eqns indicated relevant contributions of the descriptors. The n-value regards the sub domain of the descriptors from the whole set incorporated. Out of the applicable domains, only 3 descriptors proved to be valuable for the model defining the chemomathematical features of the molecular set. Among the several QSAR models constructed the best regression Eqn obtained is represented in Eqn. 1, pMIC₅₀=0.8084 (SssssCcount)–0.6067 (SulfursCount)+0.0418 (HydrogensCount)–2.4845.

Table 2: Optimization for Choosing Training and Test Set By Manual Method Followed By PLS

The statistical data of the 2D QSAR model explains 81.9 % (r² = 0.8190) of the total variance in the training set as well as it has internal (q²_{_LOO}) and external (pred_r²) predicative ability of 77.11 and 53.21 %, respectively. The low standard error of r²_se=0.2633, q²_se=0.3321 and pred_r²se=0.2291 demonstrated accuracy of the model. The F-test=58.78 % showed the statistical significance of the model, which means that probability of failure of the model was 1 in 10 000. In addition randomization test showed confidence of 99.9 % that the generated model is not random and hence it is selected as the QSAR model.

Here, 3 descriptors, SssssCcount, SulfursCount, HydrogensCount have contributed for the estimation of 2D parameter of the QSAR model as shown in figure 3. Their effective contribution can be sustained with their correlation matrix (Table 3), which indicated that they were independent of each other, so reliability is good without the biased result. SssssCcount descriptor defines the total number of carbon connected with four single bonds and under the subclass of E-state number under the physicochemical class with positive coefficient value. It accounted for the bond formation of carbon atom within the system which can contribute to electropotential, steric hindrance as well hydrophobicity of the compounds. SulfursCount descriptor signifies number of sulphur atoms in a compound and defined subclass is element count with negative coefficient value. Sulfur atom positioned in space effects both hydrophobic and electronic orientation along with the steric effect in the molecular space by its atomic contribution to the molecular set.

Table 3: Correlation Matrix of Descriptors Used In 2d QSAR Model

HydrogensCount descriptor signifies number of hydrogen atoms in a compound and it is considered under the subclass of element count with positive coefficient value. Number of hydrogen greatly affects the hydrogen bond within the system both inter and intra hydrogen bond orientation will directly affects the hydrophobicity of the molecular set.

Several 3D QSAR models were generated using SW kNN-MFA method and among these models the best model was selected on the basis of values of statistical parameters. The best SW kNN-MFA 3D QSAR models with 23 training compounds have a q² (r²cv) value of 0.7601 and pred_r² value of 0.5567 was considered to be the best. The kNN-MFA QSAR method explores formally the active analogue approach which implies that compounds display similar profiles of pharmacological activities. In this method the activity of each compound was predicted as average activity of k most chemically similar compounds from that data set.

In 3D QSAR studies, 3D lattice generated around 2,4-diaminouinazoline pharmacophore were used to optimize the electrostatic and steric requirements of the nucleus for antiTB activity. The generated 3D grid box, assisted in designing new compounds, which are based on the variation of the field values at the chosen points using the most active molecule and its nearest neighbour set. The points generated in SW KNN MFA 3D QSAR model are H_1014 (0.550626, 0.550762), S_669 (-0.480982, -0.473374) (figure 4a). These are the interaction energy of the set with the virtual receptor. The groups on the molecules near the lattice points need to have the interaction energy with the lattice point within the given range to be active.

H_1014 (0.550626 0.550762), the positive range indicated that at lattice point 1014 positive hydrophobic potential is favourable for increase in the activity of the molecular set and hence a hydrophobic substituent would be favourable in that region to get more potent congeners. Alkyl group or benzene ring system as well as benzene ring substituted with appropriate functional group, which contributed to increase hydrophobic parameter will be suitable for getting more potent molecular set.

S_669(-0.480982-0.473374), the negative range indicated that at lattice point 669 negative steric potential is favourable for increase in the activity and hence less bulkier substituent is preferred in this region. Hence, bulkier group like di- or tri- substituted hetero-phenyl ring system should be removed from that particular region or else benzene ring with single substitution, hetero-phenyl ring system or hetero atomic functional group should be more preferred for better molecular lead. 2D QSAR model validation was done by estimating various statistical parameter considerations. The q²_{_LOO} value showed the robustness of the model developed, which was within acceptable limits[19]. The Z score value of the model prepared were in acceptability range and also the residual error values were cross checked (Table 4). All the Z score values were noted when the model was prepared through the defined way and all the factors are within the acceptable limits, shown in the Table 5 which ensures the models reliability and acceptability that contributes directly to the validation of the prepared 2D QSAR model.

Table 4: Residual Error Values of 2d QSAR and 3d QSAR Models

Table 5: Validation Parameters of the 2d QSAR Model

For 3D QSAR, the residual error is given in Table 4. The actual vs. predicted and vice versa were plotted to estimate the reliability[20-21] of the 3D QSAR model’s predictability. The high value of cross validation parameter did not ensure reliable predictability of the model which was the key feature of a QSAR model[22]. Here, R² of the two regression lines were equal and the slope values have differed by less than 0.1 which ensured the reliability of the model prepared shown in the figure 4b. The residual error of each compound was less than 0.01 which the acceptable limit and confirmed that the model has a predictability of less error chance and failure.

Pharmacophore is a 3D description of the features needed for activity, which included hydrogen bond donors, hydrogen bond acceptors, aromatic ring, hydrophobic groups, electron donating group, electron accepting groups, which finally contributed to the chemical features of the moiety to be active. A set of pharmacophore hypothesis was generated using Mol Sign, each of the set contained the following features, hydrogen bond acceptors (blue colored cage), aromatic ring (yellow colored cage) and electronegative group (orange colored cage). The ideal bond distance of the crucial points in the structure was noted in Table 6 and figure 4c. The 1 and 3 positioned nitrogen atoms of the quinazoline ring system are important feature, which can generate electron donor and hydrogen bond acceptor regions. The benzene ring systems are useful as the electron cloud rich points for nucleophilic interaction and the electronegative fluorine substitution over the 4-aminophenolate group at position 4 of the quinazoline ring system is another crucial point for interaction as electronegative atomic region. The overlapped phamacophore map of the whole molecular set also concluded the same chemical features shown in figure 4d. The pharmacophoric features explored here are very much significant in contributing potent activity; these were considered for designing the new molecular set.

Table 6: Bond Distance Details of the Pharmacophoric Regional Points

The probable mechanism of action of the 2,4-diaminoquinazoline moiety has not been well defined or being reported, only has been reported to be antiTB activity in defined medium against virulent strain of M. tuberculosis. A set of the reported proteins for antiTB purpose was used for the docking to investigate the probable interacting points and the active site in macromolecule. The following proteins were under taken for docking study: enoyl ACPreductase (PDB ID 1ZID, 2IDZ)[23,24], caseinolytic proteases (PDB ID 2CE3)[25], malate synthase (PDB ID 3S9I)[26], MmpL membrane protein (PDB ID 4Y0L)[27]. The macromolecules are reported to have active sites for the antiTB molecules as well as can be serve as attractive drug target for the development of structure based drug design. Docking procedure includes both active site docking and blind docking. Only those water molecules were considered which are present within the active site. The procedure was carried out using Hermes 5.2 software. In case of blind docking no water molecules were considered as it was exceeding the limit of number of water molecule extracted by the software (i.e., 25 molecules).

Active site docking is performed with enoyl ACP reductase enzyme while caseinolytic proteases, MmpL membrane protein were subjected for the blind docking. But the molecular set of 2,4-diaminoquinazoline neither shown good docking score nor any hydrogen bond interaction was observed at the reported active site as well as in case of blind docking. Its interactive protein doesn’t inlcude the list used and demands for more research which open up for more possible ideas.

Validation of the method was done using rescoring method for the blind docking and for the active site docking or direct docking superimpose fitting method using Pymol 2.0.7. The scoring vs rescoring of the total molecular set and vice versa were plotted to estimate the reliability of the docking method used. The R² of the two regression lines are equal and the slope values have differed by less than 0.1, which ensured the reliability of the method. On the other side, the fitting was done over 42 atoms of the 81 atoms constituted co-crystal system extracted from the reported pdb file over the docked conformer co-crystal system, which contributes over 50 % of the atom fitted within 2Ǻ of acceptability range (figure 5a).

The results of 2DQSAR such as SssssCcount, SulfursCount and hydrogen count, pointed towards the hydrophobicity and steric potential of the molecules, which strongly supported the results of 3D QSAR.

So, an optimum balance of hyrdophobicity and steric potential distribution is required in the molecular system for developing better moieties. For the development of better molecular set, the following major points of the QSAR results was considered; instead of huge bulky group, substitution was done with those functional groups which contributed hydrophobicity to the whole molecular system with less steric hindrance, such as, benzene ring, heterocyclic ring, hetero-alicyclic ring with very less further substitution on them. Among the molecular set, compound 10 was reported to be the most active molecule (reported MIC₅₀= 469, -log MIC₅₀= -2.6707), which was considered as pharmacophoric reference moiety. Taking all considerations, a set of 1000 molecular system was generated. The activity of designed compounds was predicted using the developed 3D QSAR model. Only the good and efficient molecule having activity greater than the reported value is reported in Table 4.

The most active compound predicted has all the quality ensured in the result of the pharmacophore and the QSAR results. Its basic template is 2,4-diaminoquinazoline ring with substitution at amino groups with piperidine and 4-aminophenolate group over the other amino position. This type of compound has nowhere being reported nor being evaluated. In our study it has estimated to be a potent moiety and structurally is possible to be a real molecule.

The molecular structure was subjected to the MM+ method then, PM3 and finally using DFT b3lyp/6-311g method. So the optimized structure (figure 5b) of the designed molecule is a stable and real possible molecule. The IR spectrum given in figure 6, is the ideal one, which provides about the functional groups present in the system. This concluded the molecule is not an arbitrary outcome and imaginary, but could exist in reality. So, further development of the molecular set could reveal whether this moiety could be a drug like molecule or further molecular set development will be required for more suitable pharmacologically active molecule.

This scaffold was indicated in a number of studies for various activities and established to be a potent lead molecule. The notable reports were cited here to summaries the importance of the scaffold and also the small variations that developed moiety to be an attractive one. A set of substituted 2,4 diaminoquinazoline reported to be antimalarial in which derivatives were prepared by substitution on both the benzenoid ring system and 2 and 4 positioned amino groups with piperizine or pyridine or closed heterocyclic ring systems simultaneously. It revealed that these small variation in the substituents and different combination of substitutions over the positions on the basic structure would lead to biologically active moieties[4].

Hao et al.[8] focused on anticancer activity of the scaffold on PAK4 protein, which is a cancer pathway protein family in which N2 and N4 substituted with either 1H-indazol-5-yl or, (1H-indol-3-yl)ethyl simultaneously with the substitution at the benzenoid ring with methoxy or halogen substitution at different position turning the basic scaffold to be a PAK4 inhibitor molecule. Horn et al.[28] reported that substitution at the N2 and N4 either with benzyl or phenyl simultaneously produced molecular set with different permutations and combinations with antibacterial activity. Besides, branched chain substitution on same position has played some role along with some more moieties were reported which were substituted at the quinazoline benzenoid ring of the system; leading to make the scaffold active against multi-drug resistant Stephylococcus aureus[28]. 4-substituted-piperazine-1-carbodithioate derivatives of 2,4-diaminoquinazoline have been reported by Cao et al., who screened those against a panel of five human tumour cell lines, A549, MCF-7, HeLa, HT29 and HCT-116 for antitumor activity and the results clearly depicted a molecular set from the scaffold to be antitumor[29].

The 2,4-diaminoquinazoline is a β-catenin inhibitor, which could disrupt the transcriptional activity of β-catenin and ultimately interfere many of the downstream genes expression. TGFβ-induced PPARγ expression was significantly inhibited by 2,4-diaminoquinazoline, which clearly demonstrated that β-catenin was an upstream regulator of PPARγ and critical for its expression[30]. The 7-(benzimidazol-1-yl) substituted derivatives of 2,4-diaminoquinazoline are DHFR inhibitors, among them compounds containing thiazol-2-yl group in the 2-position of the benzimidazole is highly potent in case of S. aureus DHFR as well as highly selective[31]. Chao et al. reported the scaffold to be active against dengue virus. The new derivatives were prepared by the substitution on the benzenoid ring at position 5,6,8 mainly with methoxy or OPh group substitution, i.e., viral target has a spatially restricted conformation that accommodates a narrow range of C-5 substituent[32].

A molecular set of substituting the amino group and benzenoid ring has been reported as antileishmanial scaffold. In some moieties the di-amino same group substitution with the benzenoid substitution over C5 and C6 positions have produced active moiety and some are outcome of the multiple substitution on the benzenoid ring system[33]. A report of 2017 revealed the potent antibacterial activities of N2- benzyl-N4-methylquinazoline-2,4-diamines against Gram-negative A. baumannii and multidrug resistant infectious species. The 6 or 7 substituted molecules displayed promising activity, with MICs ranging from 0.5 to 30 μM. The most potent in vitro activities were obtained with quinazoline-2,4-diamines bearing a N2- benzyl moiety and a N4-methyl group[34].

A well-defined protocol of computational experiments were followed in the present study, which managed to reach a conclusion that 2,4-diaminoquinazoline scaffold has high potential to be developed as a lead molecular set which could be subjected to lead optimization. In accordance to QSAR analysis an optimum balance of hydrophobicity and steric potential distribution is required in the molecular system for better moieties. Similar findings were observed in pharmacophore mapping, where 1 and 3 positioned nitrogen atoms of the quinazoline ring are important along with the benzene ring and the electronegative fluorine substitution over the 4-aminophenolate group at position 4 of the quinazoline ring system. Designing and evaluation of the activity through 3D QSAR model revealed that the substitution of piperidine ring at 2-amino position and 4-aminophenolate group over the 4-amino position is prime requirement as observed from QSAR and pharmacophore mapping results. The final molecule reported included all the features estimated and believed to be an active moiety. Possibly larger library generation might be required for identifying a potential candidate drug molecule. As the biological cell interaction functionality has not been checked in a conventional way, but has been checked through the protein interaction docking method, a huge probability existed that in the long run a compound could be identified, which could serve as a potential antiTB drug.

Acknowledgement:

One of the authors PB wishes to thanks AICTE for the grant of fellowship.

Conflict of interest:

The authors confirm that this article content has no conflict of interest.

References

1. Who TB sheet for drug resistant TB [cited 2018 March 24]. Available from: http://www.who.int/tb/areas-of-work/drug-resistant-tb/en.
2. WHO Global Tuberculosis report 2017 [cited 2018 February 21]. Available from: http://www.who.int/tb/publications/globalreports /en.
3. Odingo J, Malley T, Kesicki EA, Alling T, Bailey MA, Early J, et al. Synthesis and evaluation of the 2,4-diaminoquinazoline series as anti-tubercular agents. Bioorg Med Chem 2014;22:6965-79.
4. Sundriyal S, Malmquist NA, Carona J, Blundelld S, Liue F, Chene X, et al. Development of Diaminoquinazoline Histone Lysine Methyltransferase Inhibitors as Potent Blood-stage Anti-Malarial Compounds. Chem Med 2014;9(10):2360-73.
5. Mishra M, Mishra VK, Senger P, Pathak AK, Kashaw SK. Exploring QSAR studies on 4-substituted quinazoline derivatives as antimalarial compounds for the development of predictive models. Med Chem Res 2014;23:1397-405.
6. Zhu X, Horn KSV, Barber M, Yang S, Wang MZ, Manetsch R, et al. SAR refinement of antileishmanial N 2,N 4-disubstituted quinazoline-2,4-diamines. Bioorg Med Chem 2015;23(16):5182-9.
7. Mohamed T, Shakeri A, Tin G, Rao PPN. Structure−Activity Relationship Studies of Isomeric 2,4-Diaminoquinazolines on β-Amyloid Aggregation Kinetics. ACS Med Chem Lett 2016;7:502-7.
8. Hao C, Huang W, Li X, Guo J, Chen M, Yan Z, et al. Development of 2, 4-diaminoquinazoline derivatives as potent PAK4 inhibitors by the core refinement strategy. Eur J Med Chem 2017;131:1-13.
9. Lin LC, Hsu SL, Wu CL, Hsueh CM. TGFβ can stimulate the p38/b-catenin/PPARc signaling pathway to promote the EMT, invasion and migration of non-small cell lung cancer (H460 cells). Clin Exp Metastasis 2014;31:881-95.
10. Wishart DS, Knox C, Guo AC, Cheng D, Shrivastava S, Tzur D, et al. Drug Bank: a knowledgebase for drugs, drug actions and drug targets. Nucleic Acids Res 2008;36:D901-6.
11. Wynne G, DeMoor O, Johnson P, Vickers R, inventor; Infectious Disease Research Institute, assignee. Use of Compounds for Preparing Anti-tuberculosis Agents. WPO 20100317607. 2008 Dec 16.
12. Ananthan S, Faaleolea ER, Goldman RC, Hobrath JV, Kwong CD, Laughon BE, et al. High-throughput screening for inhibitors of Mycobacterium tuberculosis H37Rv. Tuberculosis 2009;89(5):334-53.
13. Halgren TA. Merck molecular force field. III. Molecularn geometries and vibrational frequencies for MMFF94. J Comput Chem 1996;17:553-86
14. Ajmani S, Jadhav K, Kulkarni SA, Three-Dimensional QSAR Using the k-Nearest Neighbor Method and Its Interpretation. J Chem Inf Model 2006;46:24-31.
15. Xu M, Zhang A, Han S, Wang L. Studies of 3d-quantitative structure–activity relationships on a set of nitroaromatic compounds: CoMFA, advanced CoMFA and CoMSIA, Chemosphere 2002;48:707-15.
16. Vitthal UB, Pratapro GS. 3D QSAR and pharmacophore modeling on substituted cyanopyrrolidines as type-II anti-diabetic agents potential dipeptidyl peptidase-IV Inhibitors. Pharmacophore 2016;7(5):342-8.
17. Jones G, Willett P, Glen RC, Leach AR, Taylor R. Development and validation of a genetic algorithm for flexible docking. J Mol Biol 1997;267:727-48.
18. Verdonk ML, Cole JC, Hartshorn HJ, Murray CW, Taylor RD. Improved Protein-Ligand Docking using GOLD. Proteins 2003;52:609-23.
19. Gramatica P. Principles of QSAR models validation: internal and external. QSAR Comb Sci 2007;26(5):694-701.
20. Zheng W, Tropsha A. Novel Variable Selection Quantitative Structure-Property Relationship Approach Based on the k-Nearest-Neighbor Principle. J Chem Inf Comput Sci 2000;40:185-94
21. Shen M, Xiao Y, Golbraikh A, Gombar VK, Tropsha A. Development and Validation of k-Nearest-Neighbor QSPR Models of Metabolic Stability of Drug Candidates. J Med Chem 2003;46:3013-3020.
22. Golbraikh A, Tropsha A. Beware of q2. J Mol Graphics Model 2002;20:269-76.
23. Rozwarski DA, Grant GA, Barton DH, Jacobs WR, Sacchettini JC. Modification of the NADH of the isoniazid target (InhA) from Mycobacterium tuberculosis. Science 1998;279(5347):98-102.
24. Dias MV, Vasconcelos IB, Prado AM, Fadel V, Basso LA, Azevedo WF, et al. Crystallographic studies on the binding of isonicotinyl-NAD adduct to wild-type and isoniazid resistant 2-trans-enoyl-ACP (CoA) reductase from Mycobacterium tuberculosis. J Struct Biol 2007;159(3):369-80.
25. Ingvarsson H, Maté MJ, Högbom M, Portnoï D, Benaroudj N, Alzari PM, et al. Insights into the inter-ring plasticity of caseinolytic proteases from the X-ray structure of Mycobacterium tuberculosis ClpP1. Acta Crystallogr D Biol Crystallogr 2007;63(2):249-59.
26. Krieger IV, Freundlich JS, Gawandi VB, Roberts JP, Gawandi VB, Sun Q, et al. Structure-guided discovery of phenyl-diketo acids as potent inhibitors of M. tuberculosis malate synthase. Chem Biol 2012;19(12):1556-67.
27. Chim N, Torres R, Liu Y, Capri J, Batot G, Whitelegge JP, Goulding CW. The Structure and Interactions of Periplasmic Domains of Crucial MmpL Membrane Proteins from Mycobacterium tuberculosis. Chem Biol 2015;22(8):1098-107.
28. Horn KSV, Burda WN, Fleeman R, Shaw LN, Manetsch R. Antibacterial Activity of a Series of N2,N4-Disubstituted Quinazoline-2,4-diamines. J Med Chem 2014;57(7):3075-93.
29. Cao SL, Han Y, Yuan CZ, Wang Y, Xiahou ZK, Liao J, et al. Synthesis and antiproliferative activity of 4-substituted-piperazine- 1-carbodithioate derivatives of 2,4-diaminoquinazoline. Eur J Med Chem 2013;64:401-9.
30. Chen Z, Venkatesan AM, Dehnhardt CM, Dos Santos O, Delos Santos E, Ayral-Kaloustian S, et al. 2,4-Diamino-quinazolines as inhibitors of β-catenin/Tcf-4 pathway: potential treatment for colorectal cancer. Bioorg Med Chem Lett 2009;19:4980-3.
31. Lam T, Hilgers MT, Cunningham ML, Kwan BP, Nelson KJ, Driver VB, et al. Structure-Based Design of New Dihydrofolate Reductase Antibacterial Agents: 7-(Benzimidazol-1-yl)-2,4-diaminoquinazolines. J Med Chem 2014;57:651-68.
32. Chao B, Tong XK, Tang W, Li DW, He PL, Garcia JM, et al. Discovery and Optimization of 2,4-Diaminoquinazoline Derivatives as a New Class of Potent Dengue Virus Inhibitors. J Med Chem 2012;55:3135-43.
33. Horn KSV, Zhu X, Pandharkar T, Yang S, Vesely B, Vanaerschot M, et al. Antileishmanial Activity of a Series of N2,N4-Disubstituted Quinazoline-2,4-diamines. J Med Chem 2014;57(12):5141-56.
34. Fleeman R, Horn KSV, Barber MM, Burda WN, Flanigan DL, Manetsch R, et al. Characterizing the Antimicrobial Activity of N2,N4-Disubstituted Quinazoline-2,4-Diamines Towards Multidrug Resistant Acinetobacter baumannii. Antimicrob Agents Chemother 2017;61(6):59-17.