Next Article in Journal
6-S-Trityl-mercapto-hexan-1-ol
Next Article in Special Issue
Headspace, Volatile and Semi-Volatile Organic Compounds Diversity and Radical Scavenging Activity of Ultrasonic Solvent Extracts from Amorpha fruticosa Honey Samples
Previous Article in Journal
Ethyl 2-(n-Butylamino)-7-hydroxy-5-methylpyrazolo[1,5-a]pyrimidine-3-carboxylate
 
 
Correction published on 26 January 2010, see Molecules 2010, 15(1), 604-605.
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On Two Novel Parameters for Validation of Predictive QSAR Models

Drug Theoretics and Cheminformatics Lab, Division of Medicinal and Pharmaceutical Chemistry, Department of Pharmaceutical Technology, Jadavpur University, Kolkata 700 032, India
*
Author to whom correspondence should be addressed.
Molecules 2009, 14(5), 1660-1701; https://doi.org/10.3390/molecules14051660
Submission received: 16 April 2009 / Accepted: 28 April 2009 / Published: 29 April 2009
(This article belongs to the Special Issue Molecular Diversity Feature Papers)

Abstract

:
Validation is a crucial aspect of quantitative structure–activity relationship (QSAR) modeling. The present paper shows that traditionally used validation parameters (leave-one-out Q2 for internal validation and predictive R2 for external validation) may be supplemented with two novel parameters rm2 and Rp2 for a stricter test of validation. The parameter rm2(overall) penalizes a model for large differences between observed and predicted values of the compounds of the whole set (considering both training and test sets) while the parameter Rp2 penalizes model R2 for large differences between determination coefficient of nonrandom model and square of mean correlation coefficient of random models in case of a randomization test. Two other variants of rm2 parameter, rm2(LOO) and rm2(test), penalize a model more strictly than Q2 and R2pred respectively. Three different data sets of moderate to large size have been used to develop multiple models in order to indicate the suitability of the novel parameters in QSAR studies. The results show that in many cases the developed models could satisfy the requirements of conventional parameters (Q2 and R2pred) but fail to achieve the required values for the novel parameters rm2 and Rp2. Moreover, these parameters also help in identifying the best models from among a set of comparable models. Thus, a test for these two parameters is suggested to be a more stringent requirement than the traditional validation parameters to decide acceptability of a predictive QSAR model, especially when a regulatory decision is involved.

1. Introduction

Quantitative structure-activity relationships (QSARs) are statistically derived models that can be used to predict the physicochemical and biological (including toxicological) properties of molecules from the knowledge of chemical structure. The structural features and properties are encoded within descriptors in numerical form. Descriptors support application of statistical tools generating relations which correlate activity data with descriptors (properties) in quantitative fashion. The description of QSAR models has been a topic for scientific research for more than 40 years and a topic within the regulatory framework for more than 20 years [1]. In the field of QSAR, the main objective is to investigate these relationships by building mathematical models that explain the relationship in a statistical way with ultimate goal of prediction and/or mechanistic interpretation. QSARs are being applied in many disciplines like drug discovery and lead optimization, risk assessment and toxicity prediction, regulatory decisions and agrochemicals [2,3,4]. One of the major applications of QSAR models is to predict the biological activity of untested compounds from their molecular structures [5]. The estimation of accuracy of predictions is a critical problem in QSAR modeling [6]. Only recently, validation of QSAR models has received considerable attention [7,8,9,10,11,12,13,14,15,16,17,18,19]. Four tools of assessing validity of QSAR models [20] are (i) randomization of the response data, (ii) cross-validation, (iii) bootstrapping, (iv) external validation by splitting of set of chemical compounds into a training and a test set and/or confirmation using an independent external validation set or external validation using a designed validation set. In order to be considered for regulatory use, especially in view of REACH (Registration, Evaluation, and Authorization of Chemicals) [1,21,22] legislation enforced in the European Union, it is widely agreed that QSARs need to be assessed in terms of their scientific validity, so that regulatory bodies have a sound scientific basis on which decisions regarding regulatory implementation can be taken. Several principles for assessing the validity of QSAR models were proposed at an International workshop held in Setubal (Portugal), which were subsequently modified in 2004 by the OECD Work Programme on QSARs [21,22]. Against this background, a review of the performance of the traditional validation parameters and the search for novel parameters which may be better metrics than the currently used ones appear to be of current need.
Recently the use of internal versus external validation has been a matter of great debate [23]. One group of QSAR workers supports internal validation, while the other group considers that internal validation is not a sufficient test for checking robustness of the models and external validation must be done. Hawkins et al., the major group of supporters of internal validation, are of the opinion that cross-validation is able to assess the model fit and to check whether the predictions will carry over to fresh data not used in the model fitting exercise. They have argued that when the sample size is small, holding a portion of it back for testing is wasteful and it is much better to use “computationally more burdensome” leave-one-out cross-validation [24,25].
An inconsistency between internal and external predictivity was reported in a few QSAR studies [26,27,28]. It was reported that, in general, there is no relationship between internal and external predictivity [29]: high internal predictivity may result in low external predictivity and vice versa.
Recently we have shown [15] that predictive R2 (R2pred) may not be a suitable measure to indicate external predictability, as it is highly dependent on training set mean. An alternative measure rm2 (based on observed and predicted data of the test set compounds) was suggested to be a better metric to indicate external predictability. But it can as well be applied for training set if one considers the correlation between observed and leave-one-out (LOO) predicted values of the training set compounds [30,31]. More interestingly, this can be used for the whole set considering LOO-predicted values for the training set and predicted values of the test set compounds. The advantages of such consideration are: (1) unlike external validation parameters (R2pred etc.), the rm2(overall) statistic is not based only on limited number of test set compounds. It includes prediction for both test set and training set (using LOO predictions) compounds. Thus, this statistic is based on prediction of comparably large number of compounds. In many cases, test set size is considerably small and regression based external validation parameter may be less reliable and highly dependent on individual test set observations. In such cases, the rm2(overall) statistic may be advantageous. (2) In many cases, comparable models are obtained where some models show comparatively better internal validation parameters and some other models show comparatively superior external validation parameters. This may create a problem in selecting the final model. The rm2(overall) statistic may be used for selection of the best predictive models from among comparable models.
Again, for an acceptable QSAR model, the average correlation coefficient (Rr) of randomized models should be less than the correlation coefficient (R) of the non-randomized model. No clear-cut recommendation was found in the literature for the difference between the average correlation coefficient (Rr) of randomized models and the correlation coefficient (R) of non-randomized model. We have used a parameter Rp2 [32] which penalizes the model R2 for the difference between squared mean correlation coefficient (Rr2) of randomized models and squared correlation coefficient (R2) of the non-randomized model.
In this paper, we demonstrate the usefulness of the parameters rm2 and Rp2 in deriving predictive QSAR models. For this task, we have chosen three different data sets of moderate to large size and developed multiple models to indicate the suitability of the parameters in QSAR studies. It may be noted here that the purpose of this paper is not to develop new QSAR models for the data sets but to explore suitability of the novel parameters rm2 and Rp2 in judging quality of predictive QSAR models.

2. Materials and Methods

2.1. The data sets and descriptors

In the present paper, three different data sets have been used for the QSAR model development: (1) CCR5 binding affinity data (IC50) of 119 piperidine derivatives [33,34,35,36]; (2) ovicidal activity data (LC50) of 90 2-(2′,6′-difluorophenyl)-4-phenyl-1,3-oxazoline derivatives [37] and (3) tetrahymena toxicity (IGC50) of 384 aromatic compounds [38]. For the three data sets (I, II and III), QSAR models were separately developed from genetic function approximation (GFA) technique [39] with 5,000 crossovers using Cerius2 version 4.10 software [40]. The descriptors used were from the classes of topological, structural, physicochemical and spatial types (vide infra).

2.1.1. Data set I

The CCR5 binding affinity data (IC50) of 119 piperidine derivatives [33,34,35,36] were converted to logarithmic scale [pIC50 = -logIC50 (mM)] and then used for the QSAR study. A total of 119 compounds were selected in our study, which are shown in Table 1. In cases of racemic compounds, only S configuration was considered for modeling because the R isomers are less potent [33,34]. For this data set, different classes of descriptors used were topological [Balaban index (Jx), kappa shape indices, Zagreb, Wiener, connectivity indices and E-state indices], structural [molecular weight (MW), numbers of rotatable bonds (Rotlbonds), number of hydrogen bond donors and acceptors and number of chiral centers], physicochemical [AlogP, AlogP98, LogP, MR and MolRef], spatial [RadOfGyration, Jurs, Shadow, Area, Density, Vm] and electronic [Apol, HOMO, LUMO and Sr] parameters. Definitions of all descriptors can be found at the Cerius2 tutorial available at the website http://www.accelrys.com.

2.1.2. Data set II

The ovicidal activity data (LC50) of 90 2-(2′,6′-difluorophenyl)-4-phenyl-1,3-oxazoline derivatives [37] were converted to reciprocal logarithmic values [pLC50 = -logLC50 (M)] which were used for the QSAR analysis. There is only one region of structural variations in the compounds, which is the R position of the phenyl ring. Thus the present QSAR study explores the impact of substitutional variations at the 4-phenyl ring of the 1,3-oxazoline nucleus on the ovicidal activity of the compounds. The structures of the compounds and associated ovicidal activities are listed in Table 2. The range of the ovicidal activity values is quite wide (6.1 log units). For this data set, only topological descriptors (Balaban J, kappa shape, flexibility, subgraph count, connectivity, Wiener, Zagreb and E-sate) along with structural parameters [molecular weight (MW), numbers of rotatable bonds (Rotlbonds), number of hydrogen bond donors and acceptors and number of chiral centers] and hydrophobic substituent constant π were used for the model development.

2.1.3. Data set III

Toxicity data (-log IGC50) (Table 3) determined against T. pyriformis [38] for 384 diverse compounds were used as the third data set. Different topological descriptors [ETA parameters [41,42] and non-ETA (Balaban J, kappa shape, flexibility, subgraph count, connectivity, Wiener, Zagreb, Hosoya and E-sate) parameters] were used to develop the models.

2.2. Model development

A model’s predictive accuracy and confidence for different unknown chemicals varies according to how well the training set represents the unknown chemicals and how robust the model is in extrapolating beyond the chemistry space defined by the training set. So, the selection of the training set is significantly important in QSAR analysis. Predictive potential of a model on the new data set is influenced by the similarity of chemical nature between training set and test set [43]. The test set molecules will be predicted well when these molecules are very similar to the training set compounds. The reason is that the model has represented all features common to the training set molecules. In this paper, for the development of models for a particular data set, standardized descriptor matrix was subjected to cluster analysis by K-nearest neighbour method [44]. After clustering, test set compounds were selected from each cluster so that both test set and training set could represent all clusters and characteristics of the whole dataset. This approach (clustering) ensures that the similarity principle can be employed for the activity prediction of the test set. Based on clustering, each data set was divided into 50 combinations of training and test sets. In each case, 75% of the total compounds were selected as training set and remaining 25% were selected as test set. Models were developed from a training set using genetic function approximation and the best model was selected from the population of models obtained based on lack-of-fit score. The selected model was then validated internally by leave-one-out method and then externally by predicting the activity values of the corresponding test set. Based on the results obtained from multiple models which are derived based on different combinations of training and test sets, we have tried to evaluate performance of different validation parameters.

2.3. Statistical methods

2.3.1. GFA

In this work, all models were developed using genetic function approximation (GFA) technique. Genetic algorithms are derived from an analogy with the evolution of DNA [39]. The genetic function approximation algorithm was initially anticipated by: 1) Holland’s genetic algorithm and 2) Friedman’s multivariate adaptive regression splines (MARS) algorithm. In this algorithm an individual or model is represented as one-dimensional string of bits. A distinctive feature of GFA is that it produces a population of models (e.g. 100), instead of generating a single model, as do most other statistical methods. Genetic algorithm makes superior models to those developed using stepwise regression techniques because it selects the basis functions genetically. Descriptors, which were selected by this algorithm, were subjected to multiple linear regression for generation of models. A “fitness function” or lack of fit (LOF) was used to estimate the quality of a model, so that best model receives the best fitness score. The error measurement term LOF is determined by the following equation:
L O F = L S E ( 1 c + d p M ) 2
In Eq. (1), ‘c’ is the number of basis functions (other than constant term); ‘d’ is smoothing parameter (adjustable by the user); ‘M’ is number of samples in the training set; LSE is least squares error and ‘p’ is total numbers of features contained in all basis functions.
Once models in the population have been rated using the LOF score, the genetic cross over operation is repeatedly performed. Initially two good models are probabilistically selected as parents and each parent is randomly cut into two pieces and a new model (child) is generated using a piece from each parents. After many mating steps, i.e., genetic crossover type operation, average fitness of models in the population increases as good combinations of genes are discovered and spread through the population. It can build not only linear models but also higher-order polynomials, splines and Gaussians. In our present work, only linear terms have been used. For the development of genetic function approximation (GFA) model, Cerius2 version 4.10 [38] has been used. The mutation probabilities were kept at 5,000 iterations. Smoothness (d) was kept at 1.00. Initial equation length value was selected as 4 and the length of the final equation was not fixed.

2.3.2. Validation parameters

2.3.2.1. Q2

In case of leave-one-out (LOO) cross-validation, each member of the sample in turn is removed, the full modeling method is applied to the remaining n-1 members, and the fitted model is applied to the holdback member. The LOO approach perturbs the data structure by removing 1/Nth compound in each crossvalidation round, thus, accomplishing an increasingly smaller perturbation with increasing N. Hence, the Q2 value of LOO approaches to that of R2, which is highly unsatisfactory [20].
Cross-validated squared correlation coefficient R2 (LOO-Q2) is calculated according to the formula:
Q 2 = 1 ( Y p r e d Y ) 2 ( Y Y ¯ ) 2
In Eq. (2), Ypred and Y indicate predicted and observed activity values respectively and Y ¯ indicate mean activity value. A model is considered acceptable when the value of Q2 exceeds 05.

2.3.2.2. R2pred

Cross validation provides a reasonable approximation of ability with which the QSAR predicts the activity values of new compounds. However, external validation gives the ultimate proof of the true predictability of a model. In many cases, truly external data points being unavailable for prediction purpose, original data set compounds are divided into training and test sets [45], thus enabling external validation. This subdivision of the data set can be accomplished in many ways, but approximately similar ranges of the biological responses and structural properties and all available structural and/or physicochemical features should be represented in both training and test sets.
Equations are generated based on training set compounds and predictive capacity of the models is judged based on the predictive R2 (R2pred) values calculated according to the following equation:
R 2 p r e d = 1 ( Y p r e d ( t e s t ) Y ( t e s t ) ) 2 ( Y ( t e s t ) Y ¯ t r a i n i n g ) 2
In Eq. (3), Ypred(test) and Y(test) indicate predicted and observed activity values respectively of the test set compounds and Y ¯ training indicates mean activity value of the training set. For a predictive QSAR model, the value of R2pred should be more than 0.5.

2.3.2.3. rm2

It has been previously shown [15] that R2pred may not be sufficient to indicate external predictivity of a model. The value of R2pred is mainly controlled by ( Y o b s ( t e s t ) Y ¯ t r a i n i n g ) 2 , i.e., sum of squared differences between observed values of test set compounds and mean observed activity values of training data set. Thus, it may not truly reflect the predictive capability of the model on a new dataset. Besides this, a good value of squared correlation coefficient (r2) between observed and predicted values of the test set compounds does not necessarily mean that the predicted values are very near to corresponding observed activity (there may be considerable numerical difference between the values though maintaining an overall good intercorrelation). So, for better external predictive potential of the model, a modified r2 [rm2(test)] was introduced by the following equation [15]:
r m ( t e s t ) 2 = r 2 ( 1 r 2 r 0 2 )
In Eq. (4), r02 is squared correlation coefficient between the observed and predicted values of the test set compounds with intercept set to zero. The value of r2m(test) should be greater than 0.5 for an acceptable model.
Initially, the concept of rm2 was applied only to the test set prediction [15], but it can as well be applied for training set if one considers the correlation between observed and leave-one-out (LOO) predicted values of the training set compounds [39,40]. More interestingly, this can be used for the whole set considering LOO-predicted values for the training set and predicted values of the test set compounds. The rm2(overall) statistic may be used for selection of the best predictive models from among comparable models.

2.3.2.4. Rp2

Further statistical significance of the relationship between activity and the descriptors can be checked by randomization test (Y-randomization) of the models. This method is of two types: process randomization and model randomization. In case of process randomization, the values of the dependent variable are randomly scrambled and variable selection is done freshly from the whole descriptor matrix. In case of model randomization, the Y column entries are scrambled and new QSAR models are developed using same set of variables as present in the unrandomized model. For an acceptable QSAR model, the average correlation coefficient (Rr) of randomized models should be less than the correlation coefficient (R) of non-randomized model. We have used a parameter Rp2 [32] in the present paper, which penalizes the model R2 for the difference between squared mean correlation coefficient (Rr2) of randomized models and squared correlation coefficient (R2) of the non-randomized model. The above mentioned novel parameter can be calculated by the following equation:
R p 2 = R 2 R 2 R r 2
This novel parameter Rp2 ensures that the models thus developed are not obtained by chance. We have assumed that the value of Rp2 should be greater than 0.5 for an acceptable model.

3. Results and Discussion

3.1. Data set I

The dataset (n = 119) was divided into training set of 89 compounds and test set of 30 compounds in 50 different combinations. Each of the 50 different training sets was then used for developing QSAR models using the genetic function approximation (GFA) technique. Each of the selected QSAR models was validated internally using the leave-one-out technique and externally using the corresponding test set compounds. All the models were also validated by the process randomization technique. From the internal validation technique, the value of Q2 was determined and from the external validation technique the value of R2pred was calculated which were then used as the parameters for determining the model predictivity. Using the process randomization technique, the average of the correlation coefficients of the randomized models (Rr) was compared with the correlation coefficient (R) of the non-randomized model. To penalize a model for the difference between the squared correlation coefficients of the randomized and the non-randomized models, the value Rp2 was also calculated.
An illustration of the results obtained for each combination studied is given in Table 4. The Q2 values obtained for all the models are well above the stipulated value of 0.5 with model no. 39 showing the highest Q2 value of 0.701. However, external validation of the models showed a wide range of variation in the values of R2pred. A very low value of R2pred is obtained for models showing high values of Q2 while models with moderate values of Q2 showed a similarly moderate values of R2pred. The value of R2pred for model no. 39 is only 0.240 which is far below the stipulated acceptable value of 0.5 although the model gives the maximum value of Q2. Similarly model no. 12 gives the lowest value of R2pred (0.117) in spite of having a quite acceptable value of Q2 (0.632). On the contrary, only model nos. 3, 6, 10, 11, 15, 18, 29, 37, 41 and 42 having Q2 values just exceeding 0.5 give values of R2pred above 0.5. Again for model nos. 41 and 42, the value of R2pred is greater the value of Q2. Thus it may be inferred that very a high value of Q2 does not indicate the model to be highly predictive while determining the activity of external dataset and also a model with high external predictivity may be poorly predictive internally. Thus the parameter, rm2(overall), was used which penalizes a model for large differences in observed and predicted activity values of the congeners. A model may be considered satisfactory when rm2(overall) is greater than 0.5.
As we know, high or acceptable values of the two parameters, Q2 and R2pred, may be obtained as long as a moderate overall correlation is maintained between the observed and predicted activity values even if there is a considerable difference between them. The parameter rm2(overall) determines whether the predicted activities are really close to the observed values or not since high values of Q2 and R2pred does not necessarily mean that the predicted values are very close to the observed ones. The value of rm2(overall) is a good compromise between a high value of Q2 and a low value of R2pred and vice versa. For models showing high acceptable values of Q2 but very low values of R2pred (below 0.5) and vice versa, it becomes difficult to conclude whether the model is well predictive or not. Similarly, the results obtained here show that some of the models give high Q2 values while others give high R2pred values. So, the selection of the best model becomes difficult. The value of rm2(overall) takes into consideration predictions for both training and test set compounds and maintains a balance between the values of Q2 and R2pred. This fact can be well established from the Figure 1 showing a comparative plot of the values of Q2, R2pred and rm2(overall) for the 50 different models (trial nos. in x axis). The line showing the values of rm2(overall) indicates that it can penalize a model with high Q2 but low R2pred. Furthermore, models with rm2(overall) values greater than 0.5 may be considered acceptable. Thus, in this dataset, although some of the models are acceptable considering the values of the conventional parameters (Q2 and R2pred), none of the models satisfy the value of rm2(overall). So none of the models obtained using the present descriptor matrix appears to be truly predictive.
In all the models developed for this dataset, there is a difference of at least 0.15 or more between the values of Q2 and rm2(LOO), the latter parameter showing lower values. Model no. 8 having an acceptable value of Q2 (0.620) may appear to be quite good at a first glance, but this model bears the maximum difference between the values of Q2 and rm2(LOO) (0.204). The rm2(LOO) parameter for a given model indicates the extent of deviation of the LOO predicted activity values from the observed ones for the training set compounds. This implies that model 8, despite having an acceptable Q2, is not capable of accurately predicting the activities of some training set molecules (7 out of 89 training set compounds have LOO predicted residuals of more than 1 log unit) and this is reflected in the value of rm2(LOO). Similar results are also obtained for model nos. 2, 9, 16, 28 and 39. Interestingly, model 39 has the maximum Q2 value (0.701) while the rm2(LOO) value of this model is only 0.517. Figure 2 shows a comparative plot of the values of Q2 and rm2(LOO) for the 50 different models.
The rm2(test) parameter determines the extent of deviation of the predicted activity from the observed activity values of test set compounds where the predicted activity is calculated on the basis of the model developed using the corresponding training set. Model nos. 3, 6, 10, 11, 15, 18 and 41 show acceptable values of R2pred and rm2(test).
Moreover, for these models the difference between the value of R2pred and rm2(test) is very low (less than 0.1) indicating that the predicted activity values of the test set compounds obtained from the corresponding models are very close to the corresponding observed activities of the compounds. Figure 3 shows a comparative plot of the values of R2pred and rm2(test) for the 50 different models.
The developed models were further validated by the process randomization technique. The values of Rr2 and R2 were determined which were then used for calculating the value of Rp2. Models with Rp2 values greater than 0.5 are considered to be statistically robust. If the value of Rp2 is less than 0.5, then it may be concluded that the outcome of the models is merely by chance and they are not at all well predictive for truly external datasets. Figure 4 shows a comparative plot of the values of R2, Rr2 and Rp2 for the 50 different models. In this work although some of the models satisfy the requirement for Rp2, they do not achieve the stipulated value of rm2(overall). Model nos. 9, 13, 24, 33, 39, 46 show acceptable values of Rp2 (above 0.5) but at the same time none of them achieve the required value (0.5) of rm2(overall). Thus it may be concluded that the different models obtained for this dataset using the given descriptor matrix do not appear to be truly predictive as none of them fulfills the requirements of both the parameters, rm2(overall) and Rp2, though many of them satisfy the conventional parameters, Q2 and R2pred.

3.2. Data set II

The total data set (n=90) was divided into training set (n=68) and test (external evaluation) set (n=22) (75% and 25% respectively of the total number of compounds) in 50 different combinations, based on clusters obtained from K-means clustering applied on standardized topological, structural and physicochemical descriptor matrix. Models were generated with topological, structural and physicochemical descriptors of each of the training sets using GFA. The predictive potentials of those models were determined on the corresponding test sets. Each of the models were validated both internally (using Q2) and externally (using R2pred). The models were further validated using process randomization technique. A comparison of statistical quality parameters and validation parameters of the models are listed in Table 5. The Q2 values of model nos. 8, 37 and 42 did not cross the stipulated value, i.e., 0.5. But, the rest 47 models successfully crossed that threshold value. A very low value of R2pred was obtained for models showing a high value of Q2 and vice versa, while models with a moderate value of Q2 showed a similarly moderate value of R2pred. As for example, model number 44 has the maximum leave-one-out (LOO) predicted variance (Q2 = 0.723), but the external predictive power of that model is very poor (R2pred = 0.136), which is far less than the threshold value, i.e., 0.5. Similarly, model number 35 has also high internal predictive variance (Q2 = 0.704), but the external predictive potential of that model is very poor (R2pred = -0.002). However, in case of model number 8, internal predictive variance (Q2 = 0.468) is quite less than the stipulated value, but the external predictive potential of that model (R2pred = 0.714) is very good. However, the models with acceptable moderate values (greater than 0.5) of LOO predicted variance (Q2) like the model nos. 4, 6, 9, 13, 15, 17, 20, 22, 25, 28, 29, 34, 36, 46, 47, 50 showed satisfactory moderate values (higher than 0.5) of external predictive variance (R2pred). This dataset also implies that very high value of Q2 does not indicate the model to be highly predictive while determining the activity of external dataset and also a model with high external predictivity may be poorly predictive internally. Thus the values of rm2(overall) were also calculated to penalize the models for large differences between observed and predictive values of the congeners.
Due to the wide distribution of the ovicidal activity among the congeners (range: 6.1 log units) acceptable values of the two parameters, Q2 and R2pred, were obtained in spite of bearing a considerable difference in numerical values of the observed and predicted activities. To penalize a model for large predicted residuals, rm2(overall) was calculated. The results obtained here show that some of the models give high Q2 values while others give high R2pred values, so for selecting the best model the values of rm2(overall) were compared. The fact that the value of r2m(overall) takes into consideration predictions for the whole dataset and maintains a compromise between the values of Q2 and R2pred is established from the Figure 5 showing a comparative plot of the values of Q2, R2pred and rm2(overall) for the 50 different models. The line showing the values of rm2(overall) indicates that it penalizes a model for large difference between Q2 and R2pred values. Models with rm2(overall) values greater than (or, at least near to) 0.5 may be considered acceptable. Thus, in this dataset, although some of the models are acceptable considering the values of the conventional parameters (Q2 and R2pred), yet none of the models satisfy the value of r2m(overall). But, the value of rm2(overall) of the model no. 22 (0.488) is very close to the predetermined criterion.
The rm2(LOO) parameter for a given model is a measure of the extent of deviation of the LOO predicted activity values from the observed ones for the training set compounds. In all the models developed for this dataset, there is a difference of at least 0.111 or more between the values of Q2 and rm2(LOO) and value of the latter parameter is always lower than the former. A very high value of Q2 may indicate the model to be well predictive internally but at the same time low value of rm2(LOO) (below 0.5) for that model indicates that there exists a considerable difference between the observed and LOO predicted activity values. Hence, it may be considered that a model predictivity improves as the difference between these two parameters [Q2 and rm2(LOO)] reduces. Model number 44 has a considerably high value of Q2 (0.723) and thus the predictive potential of the model may appear to be a highly acceptable but the LOO predicted residuals of 13 compounds (out of 68) in the training set are more than 1 log unit. This has not been reflected in the Q2 value while rm2(LOO) value of the model is comparatively much lower (0.551). Thus the parameter rm2(LOO) has been able to capture the information on deviation of LOO predicted values from the observed ones for the training set compounds more efficiently and it may serve as a more strict parameter than Q2 for internal validation. Figure 6 shows a comparative plot of the values of Q2 and rm2(LOO) for the 50 different models. Similarly, rm2(test) parameter determines the extent of deviation of the predicted activity from the observed activity values for the test set compounds. Model number 25 has an acceptable value of R2pred (0.525) but the predicted residuals of 6 compounds (out of 22 compounds) in the test set are more than 1 log unit. Though the model bears an acceptable value of R2pred (0.525), the model can not be concluded to be truly predictive externally and it has not been reflected in the value of R2pred. However, the value of rm2(test) (0.484) has not crossed the threshold value of 0.5. Thus rm2(test) appears to be a more stringent parameter than R2pred for external validation. Figure 7 shows a comparative plot of the values of R2pred and rm2(test) for the 50 different models.
Robustness of the models relating the ovicidal activity with selected descriptors was judged by randomization (Y-randomization) of the model development process. To penalize the model R2 for the difference between Rr2 and R2, Rp2 was also determined. Figure 8 shows a comparative plot of the values of R2 and Rp2 for the 50 different models. In this data set, the values of Rp2 of 23 models out of 50 models crossed the threshold value of 0.5 and thus those models may be considered to be statistically robust. But, at the same time if the value of rm2(overall) is considered then those models are not acceptable since none of them achieve the required value (0.5) of rm2(overall). But, we mentioned previously that the value of rm2(overall) of the model number 22 (0.488) is very close to the required value (0.5) and that model has also acceptable value of Rp2 (0.522). These results thus suggest that this combination of training and test sets is the best one out of the 50 combinations.

3.3. Data set III

Based on cluster analysis applied on standardized descriptor matrix, the dataset (n=384) was divided into training set of 288 compounds and test set of 96 compounds in 50 different combinations. Each of the 50 different training sets was then used for developing QSAR models using the genetic function approximation (GFA) technique. Each of the best QSAR models obtained from training set was validated internally using the leave-one-out technique and externally using the corresponding test set compounds to determine the values of Q2 and R2pred respectively which were used for determining model predictivity. The models were also validated by the process randomization technique and the values of Rr and R were calculated to obtain the value of Rp2 which penalizes the models for differences in the values of Rr2 and R2.
The results of the above-mentioned 50 different trials are shown in Table 6. For this dataset all the 50 models passed the critical value (0.5) for Q2 (Q2 ranging from 0.660 to 0.774) while only two models (37, 23) failed to cross the 0.5 limit for R2pred (R2pred ranging from 0.384 to 0.834). For all the models the difference between R2 and Q2 values is not very high (less than 0.3). As illustrated in Table 6 that models with maximum internal predictive variance do not correspond to model with maximum external prediction power and vice versa. Trial 50 has the highest Q2 value (0.774) but the corresponding predictive R2 value is 0.596. On the other hand trial 45 shows the maximum value of R2pred (0.834) and the corresponding Q2 value is 0.677. Models with small differences in the above two parameters values are observed in the trials (6, 10, 13, 18, 27, 33, 35, 37 and 40). Large differences in the values of the parameters are observed in trials 1, 9, 15, 20, 25, 42 and 50. Except models 37 and 23 all the other models are statistically acceptable (Q2> 0.5 and R2pred> 0.5). Thus for selecting the best model, values of rm2(overall) for all the models was determined. As shown above, this parameter penalizes a model for large differences in observed and predicted activity values of the congeners.
Similar to the results obtained for the two datasets mentioned above, Table 6 also corresponds to the fact that the parameter, rm2(overall) penalizes a model for wide difference in the values of Q2 and R2pred. This fact can be further established from the Figure 9 showing a comparative plot of the values of Q2, R2pred and rm2(overall) for the 50 different models. For this data set all the models have the rm2(overall) value above 0.5 (0.631-0.699). The best model according to r2m(overall) is obtained from trial 30 and the corresponding Q2 and R2pred values are 0.771 and 0.692 respectively. It is obvious none of the parameter (Q2 and R2pred ) has its maximum value for this trial, however the overall parameter, rm2(overall), shows a maximum.
Besides rm2(overall), we have calculated rm2(test) and rm2(LOO) values for all the 50 trials. These two parameters signify the differences between the observed and predicted activities of the test and training set compounds in that order. For an ideal predictive model, the difference between R2pred and rm2(test) and difference between Q2 and rm2(LOO) should be low. Large difference between the values of R2pred and rm2(test) and that between Q2 and rm2(LOO) will ultimately lead to poor values of rm2(overall) parameter. Figure 10 shows a comparative plot of the values of Q2 and rm2(LOO) for the 50 different models while Figure 11 shows a comparative plot of the values of R2pred and rm2(test) for the 50 different models. For this data set, the difference between Q2 and rm2(LOO) is quite less (-0.008 to 0.057) and that between R2pred and rm2(test) is also very less (-0.019 to 0.099). Thus indicates that the models obtained for this data set using the topological descriptors are quite robust and predictive.
Further validation of the developed models by the randomization technique and the subsequent calculation of the value of Rp2 yielded results showing that none of the models developed were by chance only and the models were statistically robust. Figure 12 shows a comparative plot of the values of R2 and Rp2 for the 50 different models. In this dataset, values of Rp2 for all the models are well above the stipulated value of 0.5 (Rp2: 0.574-0.695) as shown in Table 6. Moreover since all the models showed acceptable values of r2m(overall), it can be concluded that besides being robust all the models developed are well predictive.

3.4. Overview

The QSAR models obtained for all the datasets considered in this work and their subsequent validation show that the parameters which are traditionally calculated during internal and external validation of models (Q2 and R2pred) are not enough for determining whether the model obtained is acceptable or not from the view point of predictability. Thus, additional parameters are needed for selecting the best model and confirming that the model obtained is robust and not by mere chance. These criteria are fulfilled by the parameters r2m(overall) and Rp2. The value of r2m(overall) determines whether the range of predicted activity values for the whole dataset of molecules are really close to the observed activity or not. Since the value of r2m(overall) takes into consideration the whole dataset, it penalizes models for differences between the values of Q2 and R2pred enabling one to select the best predictive model. The value of Rp2, on the contrary, determines whether the model obtained is really robust or obtained as a result of chance only. Hence it can be inferred that if the values of rm2(overall) and Rp2 are equal to or above 0.5 (or at least near 0.5), a QSAR model can be considered acceptable. Finally it can be inferred that selection of QSAR models on the basis of Q2 and R2pred may mislead the search for the ideally predictive model. The selection of robust and well predictive QSAR models may be done merely on the basis of the two parameters, rm2(overall) and Rp2, in addition to the conventional parameters. Consideration of these parameters helps one to develop more stringent models which can be successfully applied to predict the activities of molecules in a truly external dataset.
The results obtained from the present study on the three data sets show that only the third data set gives Q2 values very close to corresponding rm2(LOO) values (Figure 10) while other two data sets show large fluctuations of Q2 values from the corresponding rm2(LOO) values, the latter being always less than the former (Figure 2 and Figure 6). The reason may be the quality of the biological activity data, apart from the performance of the selected descriptors to explain a particular biological activity in relation to the structural features. In case of data sets I and III, the biological activity data are satisfactorily distributed (Figure 13), while in case of data set II the distribution is not satisfactory. Thus, for data set I, the differences between Q2 and corresponding rm2(LOO) values may be attributed to the inability of the selected descriptors to satisfactorily explain the change of biological activity values with changes in structural features while in case of the second data set, it may be due to unsatisfactory distribution of the biological activity values.
It may be noted here that rm2 values do not take into account the number of predictor variables included in a model. When different models, having different number of predictor variables are compared then it may be very difficult to determine which one is the best model as rm2 does not consider the number of predictor variables used. To solve this problem, another parameter [rm2(overall)(adjusted)] may be calculated in a manner similar to the adjusted R2 (R2a):
r m ( o v e r a l l ) 2 ( a d j u s t e d ) = ( n 1 ) r m ( o v e r a l l ) 2 p n p 1
In Eq. (6), n is the total number of compounds and p is the number of predictor variables. The values of the parameter rm2(overall)(adjusted) for all the models of data sets I, II and III have been shown in Table 4, Table 5 and Table 6 respectively.

4. Conclusions

QSAR models have been traditionally tested for their predictive potential using internal (Q2) and external validation (R2pred) parameters. The present study shows that even in presence of considerable differences between observed and LOO predicted values of the training set compounds, Q2 value may be considerably high thus not reflecting bad predictions for some compounds. The parameter rm2(LOO) is a stricter metric for internal validation than Q2. Similarly rm2(test) appears to be a better metric to denote external predictivity than the traditional parameter R2pred. The parameter rm2(overall) is unique in that it considers predictions for both training and test set compounds and its value is not obtained from prediction of limited number of test set compounds as is the case for R2pred. In addition to this, rm2(overall) helps to identify the best model from among comparable models, especially when different models show different patterns in internal and external predictivity. The parameter Rp2 penalizes model R2 for large differences between determination coefficient of nonrandom model and square of mean correlation coefficient of random models in case of a randomization test and thus confirms whether a model has been obtained by chance or not. A model can be considered robust, truly predictive and not obtained by chance when the parameters rm2 (all three variants) and Rp2 cross the minimum limit of 0.5 (or at least near 0.5). Thus, in addition to the traditional validation parameters, tests for rm2 and Rp2 should be carried out for a more stringent test of validation of predictive QSAR models, especially when a regulatory decision is involved.

Supplementary Files

Supplementary File 1

Acknowledgements

The authors thank Gopinath Ghosh and Asim Sattwa Mandal for their help in computation of the descriptors. Financial support under a Major Research Grant of University Grant Commission (UGC), New Delhi is thankfully acknowledged. One of the authors (P. P. Roy) thanks the UGC, New Delhi for a fellowship.

References and Notes

  1. Zvinavashe, E.; Murk, A.J.; Rietjens, I.M.C.M. Promises and pitfalls of quantitative structure-activity relationship approaches for predicting metabolism and toxicity. Chem. Res. Toxicol. 2008, 21, 2229–2236. [Google Scholar] [CrossRef] [PubMed]
  2. Perkins, R.; Fang, H.; Tong, W.; Welsh, W.J. Quantitative structure-activity relationship methods: perspectives on drug discovery and toxicology. Environ. Toxicol. Chem. 2003, 22, 1666–1679. [Google Scholar] [CrossRef] [PubMed]
  3. Yang, G.F.; Huang, F. Development of Quantitative Structure-Activity Relationships and Its Application in Rational Drug Design. Curr. Pharm. Des. 2006, 12, 4601–4611. [Google Scholar] [CrossRef] [PubMed]
  4. Mazzatorta, P.; Benfenati, E.; Lorenzini, P.; Vighi, M. QSAR in ecotoxicity: an overview of modern classification techniques. J. Chem. Inf. Comput. Sci. 2004, 44, 105–112. [Google Scholar] [CrossRef] [PubMed]
  5. Konovalov, D.A.; Llewellyn, L.E.; Heyden, Y.V.; Coomans, D.J. Robust cross-validation of linear regression QSAR models. Chem. Inf. Model. 2008, 48, 2081–2094. [Google Scholar] [CrossRef] [PubMed]
  6. Tetko, I.V.; Sushko, I.; Pandey, A.K.; Zhu, H.; Tropsha, A.; Papa, E.; Oberg, T.; Todeschini, R.; Fourches, D.; Varnek, A. Critical assessment of QSAR models of environmental toxicity against Tetrahymena pyriformis: focusing on applicability domain and overfitting by variable selection. J. Chem. Inf. Model. 2008, 48, 1733–1746. [Google Scholar] [CrossRef] [PubMed]
  7. Golbraikh, A.; Tropsha, A. Beware of q2! J. Mol. Graphics Mod. 2002, 20, 269–276. [Google Scholar] [CrossRef]
  8. Tropsha, A.; Gramatica, P.; Gombar, V.K. The Importance of Being Earnest: Validation is the Absolute Essential for Successful Application and Interpretation of QSPR Models. QSAR Comb. Sci. 2003, 22, 69–77. [Google Scholar] [CrossRef]
  9. Tong, W.; Xie, Q.; Hong, H.; Shi, L.; Fang, H.; Perkins, R. Assessment of prediction confidence and domain extrapolation of two structure-activity relationship models for predicting estrogen receptor binding activity. Environ. Health Perspect. 2004, 112, 1249–1254. [Google Scholar] [PubMed]
  10. Aptula, A.O.; Jeliazkova, N.G.; Schultz, T.W.; Cronin, M.T.D. The better predictive model: High q2 for the training set or low root mean square error of prediction for the test set? QSAR Comb. Sci. 2005, 24, 385–396. [Google Scholar] [CrossRef]
  11. He, L.; Jurs, P.C. Assessing the reliability of a QSAR model’s predictions. J. Mol. Graphics Mod. 2005, 23, 503–523. [Google Scholar] [CrossRef] [PubMed]
  12. Ghafourian, T.; Cronin, M.T.D. The impact of variable selection on the modelling of oestrogenicity. SAR QSAR Environ. Res. 2005, 16, 171–190. [Google Scholar] [CrossRef] [PubMed]
  13. Roy, K.; Leonard, J.T. On selection of training and test sets for the development of predictive QSAR models. QSAR Comb. Sci. 2006, 25, 235–251. [Google Scholar]
  14. Kolossov, E.; Stanforth, R. The quality of QSAR models: problems and solutions. SAR and QSAR Environ. Res. 2007, 18, 89–100. [Google Scholar] [CrossRef] [PubMed]
  15. Roy, P.P.; Roy, K. On some aspects of variable selection for partial least squares regression models. QSAR Comb. Sci. 2008, 27, 302–313. [Google Scholar] [CrossRef]
  16. Roy, P.P.; Leonard, J.T.; Roy, K. Exploring the impact of the size of training sets for the development of predictive QSAR models. Chemom. Intell. Lab. Sys. 2008, 90, 31–42. [Google Scholar] [CrossRef]
  17. Schuurmann, G.; Ebert, R.U.; Chen, J.; Wang, B.; Kuhne, R. External validation and prediction employing the predictive squared correlation coefficient test set activity mean vs training set activity mean. J. Chem. Inf. Model. 2008, 48, 2140–2145. [Google Scholar] [CrossRef] [PubMed]
  18. Hawkins, D.M.; Kraker, J.J.; Basak, S.C.; Mills, D. QSPR checking and validation: a case study with hydroxy radical reaction rate constant. SAR and QSAR Environ. Res. 2008, 19, 525–539. [Google Scholar] [CrossRef] [PubMed]
  19. Benigni, R.; Bossa, C. Predictivity of QSAR. J. Chem. Inf. Model. 2008, 48, 971–980. [Google Scholar] [CrossRef] [PubMed]
  20. Wold, S.; Eriksson, L. Chemometrics Methods in Molecular Design; van de Waterbeemd, H., Ed.; VCH: Weinheim, Germany, 1995; pp. 309–318. [Google Scholar]
  21. Gramatica, P. Principles of QSAR models validation: internal and external. QSAR Comb. Sci. 2007, 26, 694–701. [Google Scholar] [CrossRef]
  22. Available online: http://ec.europa.eu/environment/chemicals/reach/reach_intro.htm accessed on 28 April 2009.
  23. Roy, K. On some aspects of validation of predictive QSAR models. Expert Opin. Drug Discov. 2007, 2, 1567–1577. [Google Scholar] [CrossRef] [PubMed]
  24. Hawkins, D.M.; Basak, S.C.; Mills, D. Assessing model fit by crossvalidation. J. Chem. Inf. Comput. Sci. 2003, 43, 579–586. [Google Scholar] [CrossRef] [PubMed]
  25. Hawkins, D.M. The problem of overfitting. J. Chem. Inf. Comput. Sci. 2003, 44, 1–12. [Google Scholar] [CrossRef] [PubMed]
  26. Novellino, E.; Fattorusso, C.; Greco, G. Use of comparative molecular field analysis and cluster analysis in series design. Pharm. Acta Helv. 1995, 70, 149–154. [Google Scholar] [CrossRef]
  27. Norinder, U. Single and domain variable selection in 3D QSAR applications. J. Chemom. 1996, 10, 95–105. [Google Scholar] [CrossRef]
  28. Kubinyi, H. A general view on similarity and QSAR studies. In Computer-Assisted Lead Finding and Optimization; van de Waterbeemd, H., Testa, B., Folkers, G., Eds.; VHChA and VCH: Basel, Weinheim, 1997; pp. 9–28. [Google Scholar]
  29. Kubinyi, H.; Hamprecht, F.A.; Mietzner, T. Three-dimensional quantitative similarity-activity relationships (3D QSiAR) from SEAL similarity matrices. J. Med. Chem. 1998, 41, 2553–2564. [Google Scholar] [CrossRef] [PubMed]
  30. Roy, K.; Roy, P.P. Comparative QSAR studies of CYP1A2 inhibitor flavonoids using 2D and 3D descriptors. Chem. Biol. Drug Des. 2008, 5, 370–382. [Google Scholar] [CrossRef] [PubMed]
  31. Roy, K.; Ghosh, G. QSTR with Extended Topochemical Atom (ETA) Indices. 10. Modeling of Toxicity of Organic Chemicals to Humans Using Different Chemometric Tools. Chem. Biol Drug Des. 2008, 5, 383–394. [Google Scholar] [CrossRef] [PubMed]
  32. Roy, K.; Paul, S. Exploring 2D and 3D QSARs of 2,4-diphenyl-1,3-oxazolines for ovicidal activity against Tetranychus urticae. QSAR Comb. Sci. 2008, 28, 406–425. [Google Scholar] [CrossRef]
  33. Dorn, C.P.; Finke, P.E.; Oates, B.; Budhu, R.J.; Mills, S.G.; MacCoss, M.; Malkowitz, L.; Springer, M.S.; Daugherty, B.L.; Gould, S.L.; DeMartino, J.A.; Siciliano, S.J.; Carella, A.; Carver, G.; Holmes, K.; Danzeisen, R.; Hazuda, D.; Kessler, J.; Lineberger, J.; Miller, M.; Schleif, W.A.; Emini, E.A. Antagonists of the human CCR5 receptor as anti-HIV-1 agents. Part 1: discovery and initial structure-activity relationships for 1-amino-2-phenyl-4-(piperidin-1-yl) butanes. Bioorg. Med. Chem. Lett. 2001, 11, 259–264. [Google Scholar] [CrossRef]
  34. Finke, P.E.; Meurer, L.C.; Oates, B.; Mills, S.G.; MacCoss, M.; Malkowitz, L.; Springer, M.S.; Daugherty, B.L.; Gould, S.L.; DeMartino, J.A.; Sicilino, S.J.; Carella, A.; Carver, G.; Holmes, K.; Danzeisen, R.; Hazuda, D.; Kessler, J.; Lineberger, J.; Miller, M.; Schleif, W.A.; Emini, E.A. Antagonists of the human CCR5 receptor as anti-HIV-1 agents. Part 2: structure-activity relationships for substituted 2-aryl-1-[N-(methyl)-N-(phenylsulfonyl) amino]-4-(piperidin-1-yl) butanes. Bioorg. Med. Chem. Lett. 2001, 11, 265–270. [Google Scholar] [CrossRef]
  35. Finke, P.E.; Meurer, L.C.; Oates, B.; Shah, S.K.; Loebach, J.L.; Mills, S.G.; MacCoss, M.; Castonguay, L.; Malkowitz, L.; Springer, M.S.; Gould, S.L.; DeMartino, J.L. Antagonists of the human CCR5 receptor as anti-HIV-1 agents. Part 3: a proposed pharmacophore model for 1-[N-(methyl)-N-(phenylsulfonyl) amino]-2-(phenyl)-4-[4-(substituted)piperidin-1-yl] butanes. Bioorg. Med. Chem. Lett. 2001, 11, 2469–2473. [Google Scholar] [CrossRef]
  36. Finke, P.E.; Oates, B.; Mills, S.G.; MacCoss, M.; Malkowitz, L.; Springer, M.S.; Gould, S.L.; DeMartino, J.A.; Carella, A.; Carver, G.; Holmes, K.; Danzeisen, R.; Hazuda, D.; Kessle, J.; Lineberger, J.; Miller, M.; Schleif, W.A.; Emini, E.A. Antagonists of the human CCR5 receptor as anti-HIV-1 agents. Part 4: synthesis and structure-activity relationships for 1-[N-(methyl)-N-(phenylsulfonyl)amino]-2-(phenyl)-4-(4-(N-(alkyl)-N-(benzyloxycarbonyl)amino)piperidin-1-yl)-butanes. Bioorg. Med. Chem. Lett. 2001, 11, 2475–2479. [Google Scholar] [CrossRef]
  37. Suzuki, J.; Tanji, I.; Ota, Y.; Toda, K.; Nakagawa, Y. QSAR of 2,4-diphenyl-1,3-oxazolines for ovicidal activity against the two-spotted spider mite Tetranychus urticae. J. Pestic. Sci. 2006, 31, 409–416. [Google Scholar] [CrossRef]
  38. Schultz, T.W.; Netzeva, T.I.; Cronin, M.T.D. Selection of data sets for QSARs: Analyses of Tetrahymena toxicity from aromatic compounds. SAR and QSAR Environ. Res. 2003, 14, 59–81. [Google Scholar] [CrossRef] [PubMed]
  39. Rogers, D.; Hopfinger, A. J. Application of genetic function approximation to quantitative structure-activity relationships and quantitative structure-property relationships. J. Chem. Inf. Comput. Sci. 1994, 34, 854–866. [Google Scholar] [CrossRef]
  40. Cerius2 Version 4.10. Accelrys Inc.: San Diego, CA, USA.
  41. Roy, K.; Ghosh, G. QSTR with extended topochemical atom (ETA) indices. 9. Comparative QSAR for the toxicity of diverse functional organic compounds to Chlorella vulgaris using chemometric tools. Chemosphere 2007, 70, 1–12. [Google Scholar] [CrossRef] [PubMed]
  42. Roy, K.; Ghosh, G. QSTR with extended topochemical atom (ETA) indices. 8. QSAR for the inhibition of substituted phenols on germination rate of Cucumis sativus using chemometric tools. QSAR Comb. Sci. 2006, 25, 846–859. [Google Scholar] [CrossRef]
  43. Eriksson, L.; Jaworska, J.; Worth, A.P.; Cronin, M.T.D.; McDowell, R.M.; Gramatica, R.P. Methods for reliability and uncertainty assessment and for applicability evaluations of classification and regression-based QSARs. Environ. Health Perspect. 2003, 111, 1361–1375. [Google Scholar] [CrossRef] [PubMed]
  44. Dougherty, E.R.; Barrera, J.; Brun, M.; Kim, S.; Cesar, R.M.; Chen, Y.; Bittner, M.; Trent, J.M. Inference from clustering with application to gene-expression microarrays. J. Comput. Biol. 2002, 9, 105–126. [Google Scholar] [CrossRef] [PubMed]
  45. Wu, W.; Walczak, B.; Massart, D.L.; Heuerding, S.; Erni, F.; Last, I.R.; Prebble, K.A. Artificial neural networks in classification of NIR spectral data: Design of the training set. Chemom. Intell. Lab. Syst. 1996, 33, 35–46. [Google Scholar] [CrossRef]
Sample availability: Not available.
Figure 1. Comparative plots of Q2, R2pred and rm2(overall) values of 50 models (data set I).
Figure 1. Comparative plots of Q2, R2pred and rm2(overall) values of 50 models (data set I).
Molecules 14 01660 g001
Figure 2. Comparative plots of Q2 and rm2(LOO) values of 50 models (data set I).
Figure 2. Comparative plots of Q2 and rm2(LOO) values of 50 models (data set I).
Molecules 14 01660 g002
Figure 3. Comparative plots of R2pred and rm2(test) values of 50 models (data set I).
Figure 3. Comparative plots of R2pred and rm2(test) values of 50 models (data set I).
Molecules 14 01660 g003
Figure 4. Comparative plots of R2, Rr2 and Rp2 values of 50 models (data set I).
Figure 4. Comparative plots of R2, Rr2 and Rp2 values of 50 models (data set I).
Molecules 14 01660 g004
Figure 5. Comparative plots of Q2, R2pred and rm2(overall) values of 50 models (data set II).
Figure 5. Comparative plots of Q2, R2pred and rm2(overall) values of 50 models (data set II).
Molecules 14 01660 g005
Figure 6. Comparative plots of Q2 and rm2(LOO) values of 50 models (data set II).
Figure 6. Comparative plots of Q2 and rm2(LOO) values of 50 models (data set II).
Molecules 14 01660 g006
Figure 7. Comparative plots of R2pred and rm2(test) values of 50 models (data set II).
Figure 7. Comparative plots of R2pred and rm2(test) values of 50 models (data set II).
Molecules 14 01660 g007
Figure 8. Comparative plots of R2, Rr2 and Rp2 values of 50 models (data set II).
Figure 8. Comparative plots of R2, Rr2 and Rp2 values of 50 models (data set II).
Molecules 14 01660 g008
Figure 9. Comparative plots of Q2, R2pred and rm2(overall) values of 50 models (data set III).
Figure 9. Comparative plots of Q2, R2pred and rm2(overall) values of 50 models (data set III).
Molecules 14 01660 g009
Figure 10. Comparative plots of Q2 and rm2(LOO) values of 50 models (data set III).
Figure 10. Comparative plots of Q2 and rm2(LOO) values of 50 models (data set III).
Molecules 14 01660 g010
Figure 11. Comparative plots of R2pred and rm2(test) values of 50 models (data set III).
Figure 11. Comparative plots of R2pred and rm2(test) values of 50 models (data set III).
Molecules 14 01660 g011
Figure 12. Comparative plots of R2, Rr2 and Rp2 values of 50 models (data set III).
Figure 12. Comparative plots of R2, Rr2 and Rp2 values of 50 models (data set III).
Molecules 14 01660 g012
Figure 13. Frequency distribution of compounds for different relative ranges of biological activity data (from low to high in log units): (a) data set I, (b) data set II, (c) data set III.
Figure 13. Frequency distribution of compounds for different relative ranges of biological activity data (from low to high in log units): (a) data set I, (b) data set II, (c) data set III.
Molecules 14 01660 g013
Table 1. Structural features and CCR5 binding affinities of piperidine containing compounds.
Table 1. Structural features and CCR5 binding affinities of piperidine containing compounds.
Molecules 14 01660 i001
Sl. No. Structural FeaturesCCR5 binding affinity (-logIC50(mM))
Number of oxygen atoms (n)R1R2YXY-ZObserved [33,34,35,36]
10(S)-3,4-Cl2-phenylPhenyl---3.000
21(S)-3,4-Cl2-phenylPhenyl---4.456
32(S)-3,4-Cl2-phenylPhenyl---4.000
41(S)-3,4-Cl2-phenyl2-Thienyl---4.222
52(S)-3,4-Cl2-phenyl 2-Thienyl---3.921
61(S)-3,4-Cl2-phenylDimethylamino---3.469
71(S)-3,4-Cl2-phenylBenzyl---3.229
81(S)-3,4-Cl2-phenylMethyl---3.071
91(S)-3,4-Cl2-phenyln-Octyl---2.854
101(S)-3,4-Cl2-phenylCyclopentyl---4.000
111(S)-3,4-Cl2-phenylCyclohexyl---4.000
121(S)-3,4-Cl2-phenyl2-Cl-phenyl---4.097
131(S)-3,4-Cl2-phenyl3-Cl-phenyl---4.155
141(S)-3,4-Cl2-phenyl 4-Cl-phenyl---4.398
152(S)-3,4-Cl2-phenyl3-NO2-phenyl---3.824
162(S)-3,4-Cl2-phenyl4-NO2-phenyl---4.222
171(S)-3,4-Cl2-phenyl4-MeO-phenyl---4.398
181(S)-3,4-Cl2-phenyl4-Phenyl-phenyl---4.398
191(S)-3,4-Cl2-phenylNaphth-1-yl---3.444
201(S)-3,4-Cl2-phenylNaphth-2-yl---4.222
211(S)-3,4-Cl2-phenylIndan-5-yl---4.155
221(S)-3,4-Cl2-phenylPyridin-3-yl---4.000
231(S)-3,4-Cl2-phenylQuinolin-8-yl---4.046
241(S)-3,4-Cl2-phenylQuinolin-3-yl---3.921
251(S)-3,4-Cl2-phenyl1-Me-imidazol-4-yl---3.469
260(R/S)-phenylPhenyl---3.347
271(R/S)-phenylPhenyl---4.456
282(R/S)-phenylPhenyl---4.523
291(R/S)-2-Cl-phenylPhenyl---2.699
302(R/S)-2-Cl-phenylPhenyl---2.886
310(S)-3-Cl-phenylPhenyl---3.569
321(S)-3-Cl-phenylPhenyl---5.000
332(S)-3-Cl-phenylPhenyl---4.824
341(S)-4-Cl-phenylPhenyl---3.569
351(S)-4-F-phenylPhenyl---3.244
361(R/S)-3,5- Cl2-phenylPhenyl---4.046
372(R/S)-3,5- Cl2-phenylPhenyl---3.959
38-Phenyl(R/S)-Phenyl-CH---3.921
39-Phenyl(R/S)-2-Cl-phenyl-CH---2.523
40-Phenyl(S)-3-Cl-phenyl-CH---4.523
41-Phenyl(S)-4-F-phenyl-CH---3.000
42-Phenyl(R/S)-3,5-Cl2-phenyl-CH---3.523
43-Phenyl(R/S)-3-F-phenyl-CH---4.000
44-Phenyl(R/S)-3-Me-phenyl-CH---4.097
45-Phenyl(R/S)-3-Et-phenyl-CH---3.959
46-Phenyl(R/S)-3-CF3-phenyl-CH---3.301
47-Phenyl(R/S)-4-Me-phenyl-CH---3.699
48-Phenyl(R/S)-3,5-Me2-phenyl-CH---3.796
49-Phenyl(R/S)-3,4-F2-phenyl-CH---3.244
50-Phenyl(R/S)-3,4-Me2-phenyl-CH---4.222
51-Phenyl(R/S)-3-Me-4-F-phenyl-CH---3.745
52-Phenyl(R/S)-3-F-4-Me-phenyl-CH---3.959
53-Phenyl3-Cl-phenyl-N---3.155
54-2-Methyl-phenyl3-Cl-phenyl-N---2.620
55-2-Methyl-phenyl3-Cl-phenyl-CH---3.398
56-2-MeO-phenyl3-Cl-phenyl-CH---4.155
57-3-CF3-phenyl3-Cl-phenyl-CH---3.921
58-4-Cl-phenyl3-Cl-phenyl-CH---3.699
59-4-F-phenyl3-Cl-phenyl-CH---4.602
60-Benzyl3-Cl-phenyl-CH---3.602
61-C6H5CH2CH23-Cl-phenyl-CH---4.187
62-C6H5CH2CH2CH23-Cl-phenyl-CH---5.301
63-----a-CH2CH2-3.745
64-----a-NHCH2-4.301
65-----a-C(O)CH2-5.301
66-----a-C(O)NH-4.347
67-----a-C(O)N(Me)4.000
68-----a-C(O)NHCH2-4.456
69-----a-NHC(O)CH2-4.456
70-----a-CH(OH)CH2-4.000
71-----CH2--O-3.585
72-MeHHO-3.000
73-t-BuHHO-3.000
74-t-BuEtHO-4.523
75-MeMeHO-3.824
76-MeEtHO-4.398
77-Men-PrHO-4.699
78-Men-BuHO-4.824
79-Men-C6H13HO-5.000
80-Mec-C6H11-CH2HO-5.222
81-MeBnHO-4.000
82-Etc-C6H11-CH2HO-4.456
83-Bnc-C6H11-CH2HO-3.097
84-EtEtHO-4.398
85-t-BuEtHO-4.602
86-c-C6H11-CH2EtHO-4.824
87-PhEtHO-5.000
88-BnEtHO-5.699
89-BnEtClO-5.699
90-BnMeHO-5.301
91-Bnn-PrHO-5.699
92-Bnn-PrClO-5.398
93-Bnn-BuHO-5.301
94-BnAllylHO-5.824
95-2-Me-C6H4-CH2n-PrHO-5.398
96-3-Me-C6H4-CH2n-PrHO-5.523
97-4-Me-C6H4-CH2n-PrHO-5.523
98-4-CF3-C6H4-CH2n-PrHO-5.222
99-4-NO2-C6H4-CH2n-PrHO-5.824
100-4-NO2-C6H4-CH2AllylHO-5.699
101-4-NO2-C6H4-CH2AllylClO-5.699
102-3-NH2COC6H4-CH2n-PrHO-6.097
103-4-NH2COC6H4-CH2n-PrHO-5.699
104-4-NH2COC6H4-CH2n-PrClO-5.523
105-Bnn-PrHO-5.699
106-MeHHNH-3.000
107-MeEtHNH-3.921
108-BnHHNH-4.000
109-Bnn-PrHNH-5.602
110-Phn-PrHNH-5.398
111-Bnn-PrHN-Me-4.699
112-(S)-α-Me-Bnn-PrHNH-4.125
113-4-NO2-BnAllylHNH-6.125
114-MeEtH--3.921
115-Phn-PrH--4.000
116-Bnn-PrH--5.523
117-PhOCH2n-PrH--5.398
118-PhCH2CH2n-PrH--4.699
119-4-NO2-BnAllylH--5.699
aThe X feature in these structures is a single bond.
Table 2. Structural features and ovicidal activity of 2-(2′,6′-difluorophenyl)-4-phenyl-1,3-oxazoline derivatives.
Table 2. Structural features and ovicidal activity of 2-(2′,6′-difluorophenyl)-4-phenyl-1,3-oxazoline derivatives.
Molecules 14 01660 i002
Sl. No.Substitution (R)Ovicidal activity
Observed [37]
1H4.71
22-CH33.74
32-Et4.76
42-OCH33.76
52-OEt3.78
62-F4.74
72-Cl5.77
83-CH33.74
93-Et3.76
103-OCH34.76
113-OEt4.78
123-F4.74
133-Cl4.77
144-CH35.74
154-Et7.76
164-i-Pr7.78
174-n-Bu8.8
184-i-Bu8.8
194-t-Bu8.8
204-n-C6H138.84
214-n-C8H178.87
224-n-C10H218.9
234-n-C12H258.93
244-n-C15H317.97
254-OH3.74
264-OCH34.76
274-OEt7.78
284-O-iPr7.8
294-n-Bu8.82
304-O-n-C8H178.89
314-O-n-C10H218.92
324-O-n-C13H277.96
334-O-n-C14H296.97
344-OCF37.84
354-OCH2CF38.85
364-SCH35.79
374-S-i-Pr5.82
384-S-NC9H196.92
394-S(=O)CH33.81
404-SO2CH32.83
414-F5.74
424-Cl7.77
434-Br7.83
444-CF36.82
454-N(CH3)23.78
464-Si(CH3)38.82
472-CH3, 4-CH33.76
482-CH3, 4-n-C8H178.89
492-CH3, 4-Cl5.79
502-OCH3, 4-t-Bu7.84
512-OCH3, 4-n-C8H176.9
522-OCH3, 4-n-C9H197.92
532-OCH3, 4-n-C10H216.93
542-OCH3, 4-F5.79
552-OCH3, 4-Cl5.81
562-OEt, 4-i-Pr6.84
572-OEt, 4-t-Bu7.86
582-OEt, 4-n-C5H118.87
592-OEt, 4-F7.81
602-OEt, 4-Cl5.83
612-OEt, 4-Br5.88
622-O-n-Pr, 4-i-Pr8.86
632-O-n-Pr, 4-t-Bu7.87
642-O-n-Pr, 4-n-C5H117.89
652-O-n-Bu, 4-t-Bu6.89
662-O-n-Bu, 4-F8.84
672-O-n-Hex, 4-t-Bu5.92
682-F, 4-Et5.79
692-F, 4-n-C6H138.86
702-F, 4-n-C7H158.88
712-F, 4-n-C8H178.89
722-F, 4-n-C10H217.92
732-F, 4-n-C12H256.95
742-F, 4-F6.77
752-F, 4-Cl8.79
762-Cl, 4-Et7.81
772-Cl, 4-i-Bu8.84
782-Cl, 4-n-C6H138.88
792-Cl, 4-n-C8H178.91
802-Cl, 4-n-C10H215.94
812-Cl, 4-n-C12H255.97
822-Cl, 4-F5.79
832-Cl, 4-Cl6.82
843-CH3, 4-CH34.76
853-F, 4-n-C6H135.86
863-F, 4-F5.77
873-F, 4-Cl6.79
883-Cl, 4-n-C6H135.88
893-Cl, 4-F5.79
903-Cl, 4-Cl5.82
Table 3. Toxicity (-log IGC50) of diverse compounds against T. Pyriformis.
Table 3. Toxicity (-log IGC50) of diverse compounds against T. Pyriformis.
Sl. NoNameToxicity [38]
13-Aminobenzyl alcohol-1.13
22-Aminobenzyl alcohol-1.07
3Benzyl alcohol-0.83
44-Hydroxyphenethyl alcohol-0.83
54-Aminobenzyl cyanide-0.76
62-Nitrobenzamide-0.72
74-Hydroxy-3-methoxybenzyl alcohol-0.7
82-Methoxyaniline-0.69
9(sec)-Phenethyl alcohol-0.66
101,3-Dihydroxybenzene-0.65
111-Phenyl-2-propanol-0.62
12Phenethyl alcohol-0.59
132-Phenyl-2-propanol-0.57
143-Amono-2-cresol-0.55
152,4,6-tris-(Dimethylaminomethyl)phenol-0.52
164-Methylbenzyl alcohol-0.49
17Phenylacetic acid hydrazide-0.48
183-Cyanoaniline-0.47
19Acetophenone-0.46
202-Methylbenzyl alcohol-0.43
21(±)1-Phenyl-1-propanol-0.43
222,3-Dimethylaniline-0.43
232,6-Dimethylaniline-0.43
242-Methyl-1-phenyl-2-propanol-0.41
25N-Methylphenethylamine-0.41
262-Phenyl-1-propanol-0.4
273-Fluorobenzyl alcohol-0.39
284-Hydroxybenzyl cyanide-0.38
294-Cyanobenzamide-0.38
302-Fluoroaniline-0.37
313,5-Dimethylaniline-0.36
32Benzyl cyanide-0.36
33Phenol-0.35
343-Methoxyphenol-0.33
352,5-Dimethylaniline-0.33
362-Methylphenol-0.29
372,4-Dimethylaniline-0.29
383-Methylaniline-0.28
39β- Methylphenethylamine-0.28
404-Methylphenethyl alcohol-0.26
41Benzylamine-0.24
422-Tolunitrile-0.24
433-Methylbenzyl alcohol-0.24
44Aniline-0.23
452-Ethylaniline-0.22
463-Nitrobenzyl alcohol-0.22
473-Phenyl-1-propanol-0.21
48Benzaldehyde-0.2
492-Phenyl-3-butyn-2-ol-0.18
501-Phenylethylamine-0.18
512-Chloroaniline-0.17
521-Phenyl-2-butanol-0.16
533,4-Dimethylaniline-0.16
542-Methylaniline-0.16
554-Methylphenol-0.16
563-Phenylpropionitrile-0.16
573-Acetamidophenol-0.16
584-Methoxyphenol-0.14
59Phenetole-0.14
603-Hydroxy-4-methoxybenzaldehyde-0.14
61Chlorobenzene-0.13
62Benzene-0.12
632-Phenyl-1-butanol-0.11
64Benzaldoxime-0.11
65Anisole-0.1
663-Fluoroaniline-0.1
672,4,5-Trimethoxybenzaldehyde-0.1
68(S±)-1-Phenyl-1-butanol-0.09
693,5-Dimethoxyphenol-0.09
703-Methylphenol-0.08
713-Phenyl-2-propen-1-ol-0.08
72α,α-Dimethylbenzenepropanol-0.07
73Propiophenone-0.07
742-Nitroanisole-0.07
754-Methylaniline-0.05
762,4,6-Trimethylaniline-0.05
772-(4-Tolyl)-ethylamine-0.04
783-Ethylaniline-0.03
793-Methoxy-4-hydroxybenzaldehyde-0.03
804-Hydroxy-3-methoxybenzonitrile-0.03
81Ethyl phenylcyanoacetate-0.02
82(R±)-1-Phenyl-1-butanol-0.01
834-Methylbenzylamine-0.01
84Thioacetanilide-0.01
853-Phenyl-1-butanol0.01
86α-Methylbenzyl cyanide0.01
874-Ethoxyphenol0.01
883-Ethoxy-4-hydroxybenzaldehyde0.02
894-Fluorophenol0.02
904-Ethylaniline0.03
913-Nitroaniline0.03
924-Chloroaniline0.05
93(±)-2-Phenyl-2-butanol0.06
94Benzyl chloride0.06
95N-Methylaniline0.06
964-Ethylbenzyl alcohol0.07
97N-Ethylaniline0.07
98Bromobenzene0.08
992-Nitroaniline0.08
1002-Propylaniline0.08
1013-Hydroxybenzaldehyde0.08
102Thiobenzamide0.09
1031-Fluoro-4-nitrobenzene0.1
1042-Bromobenzyl alcohol0.1
1054-Methoxybenzonitrile0.1
1063,5-Dimethylphenol0.11
1073-Nitrobenzaldehyde0.11
1084-Phenyl-1-butanol0.12
1094/-Hydroxypropiophenone0.12
1102-iso-Propylaniline0.12
1113,4-Dimethylphenol0.12
1122,3-Dimethylphenol0.12
1134-Chlororesorcinol0.13
1142,4-Dimethylphenol0.14
1152-(4-Chlorophenyl)-ethylamine0.14
116Nitrobenzene0.14
1172,5-Dimethylphenol0.14
1184-Phenylbutyronitrile0.15
1193-Chlorobenzyl alcohol0.15
1202-Anisaldehyde0.15
1212-Ethylphenol0.16
1224-Chlorobenzylamine0.16
123(±)-1-Phenyl-2-pentanol0.16
124Cinnamonitrile0.16
1252-Nitrobenzaldehyde0.17
126Thioanisole0.18
1272-Chloro-4-methylaniline0.18
1284-iso-Propylbenzyl alcohol0.18
129Phenyl-1,3-dialdehyde0.18
1302-Fluorophenol0.19
1314-Nitrobenzaldehyde0.2
1324-Ethylphenol0.21
133Butyrophenone0.21
1344-iso-propylaniline0.22
1353-Chloroaniline0.22
1364-(Dimethylamino)-benzaldehyde0.23
1373-Anisaldehyde0.23
1381-Fluoro-2-nitrobenzene0.23
1394-Xylene0.25
140Toluene0.25
1414-Methylanisole0.25
1424-Chlorobenzyl alcohol0.25
1432,4-Dihydroxyacetophenone0.25
1442-Nitrotoluene0.26
145Pentafluoroaniline0.26
1462-Phenylpyridine0.27
1473-Hydroxy-4-nitrobenzaldehyde0.27
1482,3,6-Trimethylphenol0.28
1493-Ethylphenol0.29
1502,6-Diethylaniline0.31
151Methyl-4-methylaminobenzoate0.31
152Benzoyl cyanide0.31
1534-Chlorophenethyl alcohol0.32
1543/-Nitroacetophenone0.32
1552-Allylphenol0.33
1565-Hydroxy-2-nitrobenzaldehyde0.33
1572-Bromophenol0.33
1582,5-Difluoronitrobenzene0.33
1594-Chloro-2-methylaniline0.35
1602-Iodoaniline0.35
1612,3,5-trimethylphenol0.36
162Iodobenzene0.36
1634-(tert)-Butylaniline0.36
1644-methyl-2-nitroaniline0.37
1652-Amino-4-(tert)-butylphenol0.37
1662-Benzylpyridine0.38
1673-Chloro-2-methylaniline0.38
1683-Chloro-4-methylaniline0.39
169Methyl-4-nitrobenzoate0.39
1704-Chlorobenzaldehyde0.4
1715-Phenyl-1-pentanol0.42
172(2-Bromoethyl)-benzene0.42
1732,4,6-Trimethylphenol0.42
1743-Nitrotoluene0.42
1752-Hydroxybenzaldehyde0.42
1761-Chloro-4-nitrobenzene0.43
177Dimethylnitroterephthalate0.43
1782-Amino-5-chlorobenzonitrile0.44
1793-Nitrobenzonitrile0.45
1804-Bromotoluene0.47
1813-Phenylpyridine0.47
1824-iso-Propylphenol0.47
1834-(tert)-Butylbenzyl alcohol0.48
184Benzhydrol0.5
1855-Chloro-2-methylaniline0.5
1863-Nitrophenol0.51
1871,2-Dichlorobenzene0.53
1882-Chloro-5-nitrobenzaldehyde0.53
1894-Chlorophenol0.54
190Phenyl propargyl sulfide0.54
1912-Chloro-5-methylphenol0.54
1922-Hydroxy-4-methoxyacetophenone0.55
1932,4-Dichloroaniline0.56
1941,2-Dimethyl-3-nitrobenzene0.56
195Valerophenone0.56
1964-Methyl-2-nitrophenol0.57
1972,5-Dichloroaniline0.58
198trans-Methyl cinnamate0.58
1991,2-Dimethyl-4-nitrobenzene0.59
2005-Chloro-2-hydroxybenzamide0.59
2015-Methyl-2-nitrophenol0.59
2024-Chloroanisole0.6
2032-Bromo-4-methylphenol0.6
2044-Bromophenyl acetonitrile0.6
2054-Butoxyaniline0.61
2064-sec-Butylaniline0.61
2073-iso-Propylphenol0.61
2082-iso-Propylphenol0.61
2093-Methyl-2-nitrophenol0.61
2104-Hydroxy-3-nitrobenzaldehyde0.61
2115-Bromovanillin0.62
212α,α,α-Trifluoro-4-cresol0.62
2134-Benzylpyridine0.63
2144-Propylphenol0.64
215Benzylidine malononitrile0.64
2164-Nitrotoluene0.65
2173-Iodoaniline0.65
218Benzyl methacrylate0.65
2194-Chlorobenzylcyanide0.66
2202-Methyl-5-nitrophenol0.66
2212-Nitroresorcinol0.66
2221-Bromo-4-ethylbenzene0.67
2234-iso-Propylbenzaldehyde0.67
2242-Nitrophenol0.67
2251,4-Dibromobenzene0.68
2262-Chloro-6-nitrotoluene0.68
2271-Chloro-2-nitrobenzene0.68
2284-Bromophenol0.68
2294-Benzoylaniline0.68
230iso-Propylbenzene0.69
2312-Chloro-4,5-dimethylphenol0.69
2324-Butoxyphenol0.7
2334-Chloro-2-methylphenol0.7
2343,5-Dichloroaniline0.71
2352-Hydroxy-4,5-dimethylacetophenone0.71
236Ethyl-4-nitrobenzoate0.71
2373-Nitroanisole0.72
2382,4-Dinitroaniline0.72
2391-Chloro-3-nitrobenzene0.73
2402,6-Dichlorophenol0.73
2413-tert-Butylphenol0.74
2421,1-Diphenyl-2-propanol0.75
2432-Chloro-4-nitroaniline0.75
2441-Bromo-2-nitrobenzene0.75
2452-Methoxy-4-propenylphenol0.75
2462-Chloromethyl-4-nitrophenol0.75
2474,5-Difluoro-2-nitroaniline0.75
2482,6-Diisopropylaniline0.76
2493-Chloro-5-methoxyphenol0.76
2504-Ethoxy-2-nitroaniline0.76
2511,3-Dinitrobenzene0.76
252α,α,α-4-Tetrafluoro-3-touidine0.77
253Ethyl-4-methoxybenzoate0.77
254(±)-1,2-Diphenyl-2-propanol0.8
2554-Chloro-3-methylphenol0.8
2563-Chloro-4-fluoronitrobenzene0.8
257Methyl-2,5-dichlorobenzoate0.81
2584-Chloro-2-nitrotoluene0.82
259Pentafluorobenzaldehyde0.82
2604-Bromophenyl-3-pyridyl ketone0.82
261Methyl-4-chloro-2-nitrobenzoate0.82
2624-Nitrophenetole0.83
2632,6-Dinitrophenol0.83
2642,6-Dinitroaniline0.84
2654-Iodophenol0.85
2661,3,5-Trimethyl-2-nitrobenzene0.86
2676-Phenyl-1-hexanol0.87
2683-Chlorophenol0.87
269Benzophenone0.87
2701,3,5-Trichlorobenzene0.87
2712,4-Dinitrotoluene0.87
2724-(tert)-Butylphenol0.91
2734-Biphenylmethanol0.92
2743,4,5-Trimethylphenol0.93
2752,2/,4,4/-Tetrahydroxybenzophenone0.96
2764-Pentyloxyaniline0.97
2772,4-Dichloronitrobenzene0.99
278(trans)-Ethyl cinnamate0.99
2794-Benzoylphenol1.02
2801-Bromo-3-nitrobenzene1.03
2812,4-Dichlorophenol1.04
2822,5-Dinitrophenol1.04
2832,4-Dichlorobenzaldehyde1.04
284Biphenyl1.05
2852,4-Dinitrophenol1.06
2864-Butylaniline1.07
2873,4-Dichlorotoluene1.07
2882,3-Dichloronitrobenzene1.07
289Benzyl-4-hydroxylphenyl ketone1.07
2901,2,4-Trichlorobenzene1.08
2914-Chloro-3-ethylphenol1.08
2921-Fluoro-3-iodo-5-nitrobenzene1.09
293Resorcinol monobenzoate1.11
2946-Chloro-2,4-dinitroaniline1.12
2954-Biphenylcarboxaldehyde1.12
2963,5-Dichloronitrobenzene1.13
2972,5-Dichloronitrobenzene1.13
2982-Bromo-5-nitrotoluene1.16
2993,4-Dichloronitrobenzene1.16
3006-tert-butyl-2,4-dimethylphenol1.16
3014-Bromo-2,6-dimethylphenol1.16
3022,2/-Dihydroxybenzophenone1.16
3033,5-Dibromo-4-hydroxybenzonitrile1.16
3044-(Pentyloxy)-benzaldehyde1.18
3054-Nitrobenzyl chloride1.18
306Hexanophenone1.19
3074-Chloro-3,5-dimethylphenol1.2
3084-tert-Pentylphenol1.23
309n-Propyl cinnamate1.23
3102-Bromo-4,6-dinitroaniline1.24
311n-Butylbenzene1.25
3121,2-Dinitrobenzene1.25
3134-Bromobenzophenone1.26
3142,4-Dichloro-6-nitroaniline1.26
3154-Phenoxybenzaldehyde1.26
3164-Chloro-3-nitrophenol1.27
3174-Bromo-6-chloro-2-cresol1.28
3182,4,5-Trichloroaniline1.3
3191,4-Dinitrobenzene1.3
3202-Nitrobiphenyl1.3
3215-Pentylresorcinol1.31
322Ethyl-4-bromobenzoate1.33
3232/,3/,4/-Trichloroacetophenone1.34
324Phenyl benzoate1.35
325Phenyl-4-hydroxybenzoate1.37
3262,5-Dibromonitrobenzene1.37
3274-Hexyloxyaniline1.38
3282,4-Dibromophenol1.4
3292,4,6-Trichlorophenol1.41
330Phenyl isothiocyanate1.41
3312-Hydroxy-4-methoxybenzophenone1.42
3321,3,5-Trichloro-2-nitrobenzene1.43
333Benzyl benzoate1.45
334iso-Amyl-4-hydroxybenzoate1.48
3352,5-Diphenyl-1,4-benzoquinone1.48
3364-Chlorobenzophenone1.5
3371,2,3-Trichloro-4-nitrobenzene1.51
3381,2,4-Trichloro-5-nitrobenzene1.53
339n-Butyl cinnamate1.53
3403-Chlorobenzophenone1.55
3413,5-Dichlorosalicylaldehyde1.55
342Heptanophenone1.56
3433,5-Dichlorophenol1.56
3444-Nitrophenyl phenyl ether1.58
3452,4-Dibromo-6-nitroaniline1.62
3464-Chloro-6-nitro-3-cresol1.63
347Pentafluorophenol1.63
3483,5-Di-tert-butylphenol1.64
3493,5-Dibromosalicylaldehyde1.65
3503-Trifluoromethyl-4-nitrophenol1.65
3514,5-Dichloro-2-nitroaniline1.66
3522,4-Dinitro-1-fluorobenzene1.71
3532-(Benzylthio)-3-nitropyridine1.72
3544,6-Dinitro-2-methylphenol1.73
3552,4-Dichloro-6-nitrophenol1.75
3562,3,5,6-Tetrachloroaniline1.76
3574-Bromo-2,6-dichlorophenol1.78
3582,3,4,5-Tetrachloronitrobenzene1.78
359n -Amylbenzene1.79
3604-Hexylresorcinol1.8
3614-(tert)-Butyl-2,6-dinitrophenol1.8
3622,6-Diiodo-4-nitrophenol1.81
3632,3,5,6- Tetrachloronitrobenzene1.82
3642,3,4,6- Tetrachloronitrobenzene1.87
365Octanophenone1.89
3661,2,3-Trifluoro-4-nitrobenzene1.89
3672,4,6-Tribromophenol1.91
3682,3,4,5-Tetrachloroaniline1.96
3694-Ethylbiphenyl1.97
3701,2,4,5-Tetrachlorobenzene2
371Pentachlorophenol2.07
3722,4,5-Trichlorophenol2.1
3732,4-Dinitro-1-iodobenzene2.12
3741-Chloro-2,4-dinitrobenzene2.16
3752,3,4,6-Tetrachlorophenol2.18
3761,3,5-Trichloro-2,4-dinitrobenzene hemihydrate2.19
3771,2-Dichloro-4,5-dinitrobenzene2.21
3781,5-Dichloro-2,3-dinitrobenzene2.42
379Nonylphenol2.47
3803,4,5,6-Tetrabromo-2-cresol2.57
3811,3-Dinitro-2,4,5-trichlorobenzene2.60
382Pentabromophenol2.66
3832,3,4,5-Tetrachlorophenol2.72
3841,4-Dinitrotetrachlorobenzene2.82
Table 4. Comparison of statistical qualities and validation parameters of different models (Data set I).
Table 4. Comparison of statistical qualities and validation parameters of different models (Data set I).
Trial No.No. of predictor variablesLOFR2Q2R2predrm2(LOO)rm2(test)rm2 (overall)rm2(overall)(adjusted)Rr2Rp2
0140.2760.6770.6420.3290.4660.3130.4440.4180.1430.495
0220.3360.5960.5690.3580.3780.3480.3830.3690.0980.421
0330.3800.5520.5110.595*0.3670.5660.4000.3790.1180.364
0440.3490.6210.5690.4380.4090.3910.4070.3790.1080.445
0540.3230.6340.5670.3670.4070.3390.4020.3740.0780.473
0620.3570.5420.5110.5420.3660.5080.3990.3850.1160.354
0730.3510.5970.5600.4360.4020.4660.4160.3950.1000.421
0840.3040.6600.6200.3460.4140.3030.4000.3710.1170.486
0930.2560.6970.6750.0800.4940.1420.4170.3960.1020.538
1040.3470.5960.5490.5300.3940.5090.4310.4040.0870.425
1130.3590.5670.5190.5560.3720.5060.4070.3860.1020.387
1230.2940.6630.6320.1170.4580.1330.4050.3840.1380.480
1340.2730.6780.6400.3260.4630.3240.4410.4140.0890.520
1430.3450.6040.5680.3900.4080.3640.4100.3890.1760.395
1530.3690.5580.5020.5230.3600.5000.3860.3640.1300.365
1630.3180.6270.5840.2820.4010.3100.3730.3510.1260.444
1740.3300.6220.5620.4620.4050.4450.4170.3890.1000.449
1840.3700.5810.5310.5420.3810.5290.4150.3870.0910.407
1940.3460.6150.5640.4470.4060.4270.4110.3830.1110.437
2030.2890.6570.6250.3010.4520.2680.4200.4000.1080.487
2130.2990.6480.6140.2540.4430.2410.4120.3910.1240.469
2230.3240.6100.5730.4260.4100.4260.4180.3970.1430.417
2330.3470.5810.5190.4710.3730.4400.3980.3770.1160.396
2440.2900.6730.6360.2380.4610.2540.4250.3980.0920.513
2530.3130.6220.5910.3430.4250.3240.4110.3900.1220.440
2640.2570.6860.6450.2330.4670.1790.4050.3770.1080.521
2740.2990.6590.6150.2120.4450.2190.4040.3760.0950.495
2840.3420.6030.5580.4970.3690.4680.3960.3670.1220.418
2940.3850.5930.5360.5440.3860.4740.3990.3700.1180.409
3040.3240.6270.5800.3940.4180.3610.4140.3860.0950.457
3140.3530.5920.5440.2860.3890.2600.3680.3380.1110.411
3230.3140.6360.6020.2640.4340.2720.4110.3900.1130.460
3350.2950.6850.6440.1790.4680.2010.4130.3780.1060.521
3420.2710.6520.6290.2630.4540.2440.4150.4010.1320.470
3540.3400.6150.5660.3030.4060.2730.3870.3580.0880.447
3640.3350.6410.6040.2860.4360.3210.4250.3980.0960.473
3730.3410.5850.5470.5170.3920.5030.4130.3920.0950.409
3830.2790.6590.6280.2530.4540.2660.4190.3980.1660.463
3940.2100.7310.701*0.2400.5170.2850.452*0.4260.1350.564
4030.3020.6300.5970.3590.4290.3670.4220.4020.1580.433
4130.3800.5580.5100.5650.3670.5420.3990.3780.1070.375
4230.4040.5570.5170.595*0.3740.5780.4030.3820.1060.374
4330.2850.6320.5850.2840.4200.2820.3960.3750.1340.446
4440.3370.6110.5650.4240.4050.4530.4210.3930.1370.421
4530.3600.6020.5670.3200.4080.2990.3980.3770.0940.429
4660.3020.6970.6460.2390.4710.2620.4310.3890.0760.549
4730.3120.6150.5730.3690.4110.3650.4110.3900.1340.427
4830.2980.6530.6170.1670.4460.1790.4040.3830.1340.470
4930.2900.6230.5890.4000.4210.4120.4240.4040.0970.452
5020.3110.5900.5610.4200.4280.4010.4150.4010.1060.410
*Models with maximum Q2, R2pred and rm2(overall) values are shown in bold.
Table 5. Comparison of statistical qualities and validation parameters of different models (Data set II).
Table 5. Comparison of statistical qualities and validation parameters of different models (Data set II).
Trial No.No. of predictor variablesLOFR2Q2R2predrm2(LOO)rm2(test)rm2(overall)rm2(overall)(adjusted)Rr2Rp2
0141.3060.6730.6170.3250.4620.2800.4260.3900.0760.520
0241.6960.5770.5100.4790.3840.4330.3930.3540.0780.408
0341.5290.6120.5590.3470.4180.3260.4080.3700.0780.447
0461.6200.6070.5170.5400.3850.4730.4150.3570.0790.441
0541.3470.6460.6060.4410.4490.4300.4440.4090.0710.490
0641.5340.6060.5480.6000.4080.5850.4370.4010.0590.448
0741.4960.6420.5850.0240.4400.1490.3720.3320.1070.470
0841.6440.5530.4680.714*0.3570.6840.4080.3700.0500.392
0941.5930.5880.5210.6330.3910.5350.4230.3860.0660.425
1021.5140.5470.5130.3250.3810.2910.3670.3480.1040.364
1151.4570.6580.5890.4480.4390.4720.4480.4030.0510.513
1241.4360.6420.5960.4700.4430.4350.4390.4030.0750.483
1341.5170.5900.5290.6130.3940.5770.4330.3970.0740.424
1441.3180.6540.6090.4430.4520.4330.4490.4140.0760.497
1541.5230.5860.5230.6520.3900.5730.4340.3980.1030.407
1641.4660.6220.5670.2030.4220.2430.3970.3590.0940.452
1761.4090.6810.6130.5970.4570.5970.4710.4190.0720.531
1851.2530.7050.6560.3510.4930.3280.4480.4030.0720.561
1951.1730.7110.6650.3310.4990.3120.4550.4110.1000.556
2051.5460.6300.5580.5070.4160.4680.4250.3790.0600.476
2141.2880.6810.636-0.0280.4770.1290.3820.3430.0560.538
2261.3490.6750.6120.6080.4570.5380.488*0.4380.0770.522
2351.3920.6600.6000.4880.4490.4670.4470.4020.0460.517
2451.3210.6800.6370.4090.4750.3740.4510.4070.0860.524
2561.3600.7010.6350.5250.4760.4840.4750.4230.0750.555
2661.2310.7220.6660.4030.5040.3630.4640.4110.0680.584
2741.1160.7080.6720.2820.5030.2540.4510.4160.0630.569
2851.3630.6480.5820.5880.4320.5520.4550.4110.0970.481
2951.4140.6270.5640.6140.4180.5720.4470.4020.1100.451
3041.2670.6730.6300.2130.4700.2600.4360.4000.0580.528
3141.4540.6260.5770.3300.4300.3020.4110.3740.0840.461
3251.5950.6130.5400.4330.4070.3490.3910.3420.0810.447
3341.4080.6330.5770.2490.4290.2480.3920.3530.0680.476
3441.5220.5860.5170.6560.3870.6350.4340.3980.0700.421
3561.0750.7580.704-0.0020.5360.1080.4220.3650.0830.623
3641.4460.5980.5350.6160.3980.5450.4450.4100.0740.433
3741.6950.5520.4860.6140.3680.5590.4090.3710.0980.372
3841.3050.6500.5960.3680.4420.4500.4430.4080.0800.491
3951.2980.6870.6160.3610.4630.3220.4370.3920.0900.531
4041.3300.6630.6170.1250.4600.1490.3970.3590.0780.507
4151.3190.6820.6200.0770.4650.1400.3930.3440.0930.523
4241.6010.5560.4850.6560.3650.6340.4130.3760.0470.396
4341.2180.6510.5880.4960.4360.4820.4440.4090.0600.500
4460.9930.7700.723*0.1360.5510.1690.4620.4090.0750.642
4541.0970.7050.6630.2000.4960.1730.4270.3910.0780.558
4651.4940.6330.5580.6360.4180.5500.4390.3940.1030.461
4751.3920.6490.5750.5450.4270.5360.4390.3940.0590.498
4851.2540.6820.6230.0770.4660.1340.3880.3390.0700.533
4941.2520.6840.6360.1510.4760.1730.4110.3740.0730.535
5051.2700.6570.5830.5560.4330.5480.4470.4020.0570.509
*Models with maximum Q2, R2pred and rm2(overall) values are shown in bold.
Table 6. Comparison of statistical qualities and validation parameters of different models (Data set III).
Table 6. Comparison of statistical qualities and validation parameters of different models (Data set III).
Trial No.No. of predictor variablesLOFR2Q2R2predrm2(LOO)rm2(test)rm2(overall)rm2(overall)(adjusted)Rr2Rp2
01080.1320.7740.7580.5510.7110.5590.6750.6660.0420.662
02080.1470.7530.7210.6410.6940.6470.6930.6840.0370.637
03080.1670.7210.6600.7500.6680.7210.6570.6470.0250.601
04070.1390.7640.7440.6850.7230.5860.6670.6590.0450.648
05060.1350.7600.6710.6810.6590.6530.6310.6230.0520.640
06070.1480.7470.7270.7030.7040.6610.6800.6720.0370.629
07060.1590.7310.7080.6120.6940.6200.6690.6620.0350.610
08070.1440.7580.7030.6410.6810.6280.6500.6410.0310.646
09070.1230.7720.7590.5720.7120.5770.6800.6720.0360.662
10090.1370.7650.7340.7420.7010.7520.6770.6670.0420.651
11090.1450.7480.7130.5830.6930.5900.6570.6460.0360.631
12080.1500.7380.6720.7340.6690.7120.6690.6600.0370.618
13120.1290.7800.7380.7160.6980.6690.6910.6780.0320.675
14090.1430.7590.7030.6220.6790.5950.6390.6270.0380.645
15090.1220.7890.7690.5450.7240.5180.6580.6470.0290.688
16070.1490.7340.6920.7530.6760.7280.6880.6800.0320.615
17070.1230.7700.7550.5950.7060.5940.6720.6640.0370.659
18090.1380.7560.7310.7410.6990.6710.6880.6780.0250.646
19070.1620.7260.6760.6780.6730.6740.6430.6340.0270.607
20070.1380.7690.7520.5770.7200.5360.6590.6500.0280.662
21080.1470.7330.6900.7310.6690.6430.6700.6610.0470.607
22080.1600.7300.6930.7310.6790.6880.6660.6560.0440.605
23060.1310.7690.7550.4970.7100.4780.6540.6470.0350.659
24090.1540.7510.7210.6350.6970.6100.6760.6660.0380.634
25060.1080.7840.7720.5750.7150.5940.6740.6670.0230.684
26080.1530.7230.6970.7810.6830.7520.6880.6790.0320.601
27080.1580.7320.7060.7440.6920.7420.6870.6780.0250.615
28080.1640.7260.6960.7360.6960.6860.6640.6540.0520.596
29070.1650.7200.6900.7460.6810.7270.6830.6750.0380.594
30090.1230.7920.7710.6920.7200.6870.699*0.6890.0520.682
31080.1180.7830.7660.5800.7160.5590.6650.6550.0320.678
32070.1620.7090.6850.7120.6790.6780.6810.6730.0400.580
33090.1440.7590.7300.7300.7050.6990.6830.6730.0340.646
34130.1540.7580.7180.6780.6990.6380.6740.6590.0250.649
35130.1300.7950.7570.7040.7150.7010.6810.6660.0330.694
36080.1460.7540.7280.5790.7030.5100.6410.6310.0350.639
37050.1350.7690.7570.3820.7200.3850.6460.6400.0320.660
38100.1510.7480.7190.6010.6930.5680.6590.6470.0330.632
39060.1640.7090.6870.7390.6810.7140.6730.6660.0340.583
40080.1530.7390.7100.7580.6920.7220.6910.6820.0370.619
41080.1640.7270.6920.6800.6840.6640.6590.6490.0320.606
42090.1390.7660.7340.5220.6970.4730.6340.6220.0360.655
43070.1470.7480.7270.6430.6990.6380.6610.6530.0390.630
44080.1670.7260.6990.6560.6840.6000.6550.6450.0310.605
45070.1680.7000.6770.834*0.6760.7530.6850.6770.0270.574
46080.1620.7080.6760.7530.6790.7250.6760.6670.0390.579
47070.1510.7360.7120.6590.6890.6690.6740.6660.0420.613
48070.1590.7230.6950.7370.6850.7140.6850.6770.0210.606
49080.1300.7810.7640.5960.7190.6100.6930.6840.0350.675
50090.1230.7920.774*0.5960.7260.5870.6780.6680.0230.695
*Models with maximum Q2, R2pred and rm2(overall) values are shown in bold.

Share and Cite

MDPI and ACS Style

Pratim Roy, P.; Paul, S.; Mitra, I.; Roy, K. On Two Novel Parameters for Validation of Predictive QSAR Models. Molecules 2009, 14, 1660-1701. https://doi.org/10.3390/molecules14051660

AMA Style

Pratim Roy P, Paul S, Mitra I, Roy K. On Two Novel Parameters for Validation of Predictive QSAR Models. Molecules. 2009; 14(5):1660-1701. https://doi.org/10.3390/molecules14051660

Chicago/Turabian Style

Pratim Roy, Partha, Somnath Paul, Indrani Mitra, and Kunal Roy. 2009. "On Two Novel Parameters for Validation of Predictive QSAR Models" Molecules 14, no. 5: 1660-1701. https://doi.org/10.3390/molecules14051660

APA Style

Pratim Roy, P., Paul, S., Mitra, I., & Roy, K. (2009). On Two Novel Parameters for Validation of Predictive QSAR Models. Molecules, 14(5), 1660-1701. https://doi.org/10.3390/molecules14051660

Article Metrics

Back to TopTop