Search tips
Search criteria 


Logo of springeropenLink to Publisher's site
Chromatographia. 2010 November; 72(9-10): 893–903.
Published online 2010 September 9. doi:  10.1365/s10337-010-1741-4
PMCID: PMC2965364

QSRR Models for Kováts’ Retention Indices of a Variety of Volatile Organic Compounds on Polar and Apolar GC Stationary Phases Using Molecular Connectivity Indexes


Quantitative structure-retention relationship (QSRR) approaches, based on molecular connectivity indices are useful to predict the gas chromatography of Kováts relative retention indices (GC-RRIs) of 132 volatile organic compounds (VOCs) on different 12 (4 apolar and 8 polar) stationary phases (C67, C103, C78, C, POH, TTF, MTF, PCL, PBR, TMO, PSH and PCN) at 130 °C. Full geometry optimization based on Austin model 1 semi-empirical molecular orbital method was carried out. The sets of 30 molecular descriptors were derived directly from the topological structures of the compounds from DRAGON program. By means of the final variable selection method, which is elimination selection stepwise regression algorithms, three optimal descriptors were selected to develop a QSRR model to predict the RRI of organic compounds on each stationary phase with a correlation coefficient between 0.9378 and 0.9673 and a leave-one-out cross-validation correlation coefficient between 0.9325 and 0.9653. The root mean squares errors over different 12 phases were within the range of 0.0333–0.0458. Furthermore, the accuracy of all developed models was confirmed using procedures of Y-randomization, external validation through an odd–even number and division of the entire dataset into training and test sets. A successful interpretation of the complex relationship between GC RRIs of VOCs and the chemical structures was achieved by QSRR. The three connectivity indexes in the models are also rationally interpreted, which indicated that all organic compounds’ RRI was precisely represented by molecular connectivity indexes.

Electronic supplementary material

The online version of this article (doi:10.1365/s10337-010-1741-4) contains supplementary material, which is available to authorized users.

Keywords: Gas chromatography, Quantitative structure–retention relationship, Kováts retention indices, Connectivity indices, Elimination selection stepwise regression algorithms


Retention is a phenomenon primarily dependent on the interactions between the solute and the stationary-phase molecules, which included directional force, induction force, dispersion force, hydrogen bond and so on. These forces can be related to the topological structures, geometric and electronic environments of the solute; therefore, it was possible to predict the solute retention from molecular parameters. Quantitative structure–retention relationships (QSRRs) have been demonstrated to be a powerful tool for the investigation of the chromatographic parameters. In QSRR approaches, the structural features of solutes encoded by non-empirical numerical descriptors and then the relationship between these descriptors and solute retention were studied [15]. The main steps in this method includes: data collection, molecular geometry optimization, molecular descriptors generation, descriptor selection, model development and finally model performance evaluation [5, 6].

The Kováts retention index is the most popular dependent variable in QSRR studies because of its reproducibility and accuracy [7]. The influence of the stationary phase polarity on the correlation and prediction of the retention of a variety of solutes by using different molecular descriptors has been the focus of several publications [814]. The development of QSRR studies is required for the suitable mathematical technique to make the models. In this regard, many available classical multivariate calibration techniques could be used to predict gas chromatographic retention index. These include multivariate linear regression (MLR), partial least-squares regression (PLS), and principal components regression (PCR) [13, 1520]. Junkes et al. [13] applied new topological index, called semi-empirical topological index (IET) by MLR to predict the chromatographic retention of aliphatic ketones and aldehydes on stationary phases of different polarities. Farkas and Héberger [15] compared ridge regression (RR), PLS, pairwise correlation (PCM), forward selection (FS), and best subset variable selection (BSS) methods for prediction of retention indices for 44 aliphatic alcohols. Initially Farkas et al. calculated 109 descriptors; they then reduced these to 17 by use of principal-component analysis. The 17 descriptors were ranked in five different ways by use of RR, PLS, PCM, FS, and BSS. Models for prediction of RIs were then built with multiple linear regression using the best three and four descriptors selected by RR, PLS, PCM, FS, and BSS. If the standard derivation of the models increase and the F values decrease, the model is unsatisfactory. Héberger et al. [1618] applied different linear multivariate techniques such as PCA, cluster, PLS and MLR in order to establish correlations between Kováts retention indices and different molecular descriptors for 35 aliphatic ketones and aldehydes, on different stationary phases at four temperatures. Tulasamma et al. [19] developed a unified QSRR with modified valence connectivity index, nχVm, based on the summation over inverse geometric mean terms by MLR to predict the retention properties of oxygen containing organic compounds on any new stationary phase. Good QSRR models were obtained to predict the retention index of 98 saturated esters on seven different polar stationary phases by MLR using the chemical descriptors proposed by Wang et al. [20].

Topological indices include valence and nonvalence molecular connectivity indices (MCIs); have been playing an important part in QSRR study for a long time. A large number of studies have demonstrated that many physicochemical and biological properties correlate with the connectivity index [2124]. The main advantage of the graph theoretical approach to the prediction of properties is that it permits the interpretation of results in terms of structurally related concepts. In spit of that, the most important criticism of the so-called topological indices is concerned with their physical meaning [25].

Molecular connectivity is a method of molecular structure quantification based only on bonding and branching patterns rather than physical or chemical characteristics. Weighted counts of substructure fragments are incorporated into numerical indices and structural features (size, branching, unsaturation, heteroatoms content and cyclicity) are encoded. These indices are related to the number of atoms and how they are connected in a molecule. Only the carbon or heavy atoms are taken into consideration and the connectivity indices are derived from the hydrogen-suppressed graph of the molecule. Each atom is represented by a vertex in the graph, while the bonds become edges. Molecular valence connectivity index [26] uses the same invariant but modifies vertex degrees to account for heteroatoms by using the number of valence electrons in the corresponding atom. The details of their definition and the calculation method can be found elsewhere [2124]. The general expression for the mth-order molecular valence connectivity index is as follows:

equation M1

where m is the order of the molecular valence connectivity index, k denotes a contiguous path type of fragment, which is divided into paths (P), cluster (C), path-clusters (PC), and chains (cycles) (CH). nm is the number of the relevant paths, and δiv is the atomic valence connectivity index and is defined as:

equation M2

where Zi is the number of electrons of atom i, Ziv the number of valence electrons, and hi the number of hydrogen connected with atom i.

The study of the relationship between quantitative topological indices and chromatographic RI could be traced back to the 1970s. Many investigators have obtained good correlations between the experimental gas chromatographic RI and structural characteristics of molecules by using different topological indices as structural descriptors [2731]. Randic [27] first studied the relationship between molecular connection index and gas chromatographic RI of alcohols. To predict chromatographic RIs of alkyl benzenes, Sutter et al. [28] chose six from 182 variables, which were from topological, geometrical, and electronic descriptors. Liu et al. [29] modeled the relationship between novel topological indices, polarizability effect index (PEI), odd–even index (OEI) and steric effect index (SVij) of 90 saturated esters and their GC RIs on seven GC columns (SE-30, OV-7, DC-710, OV-25, XE-60, OV-225 and Silar-5CP) by the MLR method. The average prediction errors over seven phases are within the range of 0.5–0.7%. Katritzky et al. [30] chose four variables from 129 topological descriptors and built a good model of 178 methyalkanes with squared correlation coefficient of 0.9585 and standard deviation of 5.8, and they particularly interpreted the relationship between variables selected in the model and molecular structures. Recently, QSRR equations have been established to model gas chromatographic retention data of alkyl pyridines on apolar and polar stationary phases by Tulasamma and Reddy [31].

The purpose of the present study was to investigate the relationship between gas chromatographic Kováts RRIs of 132 VOCs having C, H, O, N, and halogen atoms and their valence and non-valence topological molecular descriptors on different 12 (4 apolar and 8 polar) stationary phases at 130 °C using elimination selection stepwise best multiple linear regression (BMLR) analysis. Moreover, molecular descriptors were discussed to explore the influence of structural features on the values of RI. This study provided a simple and straightforward way to predict the RIs of VOCs from their structures and gave some insight into structural features related to the retention of the molecules. It has been found that the QSRR models for each stationary phase have not only high estimation qualities and high stabilities but also good predictive potentials.

Experimental and Methodology

Kováts Retention Indices

The QSRR treatment started with the assembly of the dataset. The chromatographic data used were obtained from scientific resources omitted hexamethyldisiloxane [3236] and consisted of Kováts gas chromatographic RI of 132 solutes consisting of alkanes, alkens, ethers, amines, alcohols, alkylbenzenes, and alkylhalides on 12 diverse stationary phases with different polarities at 130 °C. The chemical structures of apolar and polar stationary phases are depicted in Fig. 1. The polar phases are all isochoric and isosteric with the C78 skeleton (Fig. 1 left side). All the polar phases have 78 heavy atoms (other than hydrogen), with the heavy atoms of the polar groups substituting for methylene or methyl groups in the hydrocarbon skeleton. The alkanes are more branched than the polar phases; hence they have a lower melting point. The alkane’s family represented in Fig. 1 (right part), only concerns compounds with an odd number of carbon atoms. The C78 and C alkanes have been derived from C67 and C103 (Eqs. 2, 3, Ref. [36]). Its detailed experimental conditions were listed in Table 1. Because of the large scale of RIs of dataset, organic compounds (max RI 1400 for tetradecane on all stationary phases), the RRIs of these compounds were recalculated by dividing the RI of a reference compound such as tetradecane as the internal standards. A complete list of the names and corresponding experimental Kováts RRI values of VOCs on each stationary phase has been categorized in Table S1 of the supplementary data.

Fig. 1
Schematic representation of the stationary phases synthesized and used by Kováts and co-workers [3235], and Laffort and co-workers [36]. The polar phases, with X and Y representing different functional groups (see also Table  ...
Table 1
Name and structure details of the stationary phases synthesized and used by Kováts and co-workersa and Laffort and co-workers (experimentallyb and interpolated or extrapolatedc)

Computer Hardware and Software

All calculations were run on a Pentium IV personal computer (CPU at 2.6 MB) under Windows XP operating system. The ISIS/Draw version 2.3 software was used for drawing the molecular structures [37]. Molecular modeling and geometry optimization were employed by HyperChem (version 7.1, HyperCube) [38]. Dragon software [39] was employed for calculation of theoretical topological connectivity indices. SPSS software (version 13.0, SPSS) was used for stepwise MLR analysis and other calculations were performed in the MATLAB (version 7.0, Math Works) environment.

Descriptor Generation and Models Developing

To obtain QSRR models, solutes must be represented using molecular descriptors. Descriptors are generated solely from the molecular structures and aimed to numerically encode meaningful features of each molecule. The calculation process of the molecular descriptors is described as below: all the two-dimensional structures of the molecules were drawn using ISIS/Draw 2.3 program [37]. Then the 3D geometry structures of the molecules were pre-optimized using MM+ molecular mechanics force filed and precisely optimized with semi-empirical Austin Model 1 (AM1) method implemented in HyperChem software package (HyperCube, version 7.1) [38]. All calculations were carried out at restricted Hartree–Fock level with no configuration interaction. The molecular structures were optimized using the Polak-Ribiere algorithm until the root mean square gradient was 0.01 kcal Å mol−1 (200 K, gas phase). In order to prevent the structures locating at local minima, geometry optimization was run many times with different stating points for each molecule. The output files exported from HyperChem software were transferred into software Dragon, developed by Todeschini et al. [39], to calculate Randic topological indices as the mostly used topological indices in the literature Randic indices [23, 40]. We computed 30 different MCIs by the Dragon software including a set of 12 connectivity and average connectivity index (6 connectivity indices 0X − 5X, 6 average connectivity indices 0XA − 5XA), 12 valence and average valence connectivity index (6 valence connectivity indices 0Xv − 5Xv, 6 average valence connectivity indices 1XAv − 5XAv), and 6 solvation connectivity index (0Xsol − 5Xsol).

The calculated molecular descriptors were collected in a data matrix (D) whose number of rows and columns were the number of molecules and descriptors, respectively. At the beginning, in order to minimize the information overlap in descriptors and to reduce the number of descriptors required in regression equation, the concept of non-redundant descriptors (NRD) [41] was used in our study. That is, when two descriptors are correlated by a linear correlation coefficient value >0.85, both descriptors are correlated with the dependent variables, the better correlation is used for the actual analysis, leaving out the descriptors showing a lower correlation. This objective-based feature selection left reduced and predictive descriptors for the studied compounds. By using these criteria, 23 out of 30 original descriptors were eliminated. These descriptors are not correlating with each other as revealed from the correlation matrix presented in Table 2. In addition, these descriptors can give some information on the affecting degree for RI of different descriptors and well understand the correlation between the experimental and calculated values. Therefore, 1Xv, 3Xv, 4Xv, 1Xsol, 2Xsol, 3Xsol and 4Xsol set of descriptors has been used in the QSRR model for all stationary phases.

Table 2
Correlation matrix for the inter-correlation of various connectivity indices of 132 solutes

In order to select the subset of descriptors that best explain Kováts RI for each stationary phase, we have used elimination selection stepwise regression (ES-SWR) algorithm to select the most appropriate [4244]. This method can be regarded as a combination of the forward and backward approaches. Stepwise model-building techniques for regression designs with a single dependent variable involve identifying an initial model, repeatedly altering the model from the previous step by adding (forward stepwise) or removing (back stepwise) a predictor variable and terminating the search when stepping does not further improve the model. The forward stepwise method employs a combination of the forward entry of independent variables and backward removal of insignificant variables. The best single predictor, which is the most significant variable, was used for the initial linear regression step. Next, descriptors were added one at a time, always adding the one that most improved the fit, until the fit was not significantly improved. Once all the significant variables were determined, the regression equation was constructed. The number of variables retained in the model is based on the levels of significance assumed for inclusion and exclusion of variables from the model for each stationary phase column.

An MLR model assumes that there is a linear relationship between the molecular descriptors of a compound, which is usually expressed as a feature vector x (with each descriptor as a component of this vector), and its target property, y. An MLR model can be described using the following equation:

equation M3

where {X1,…, Xk} are molecular descriptors, β0 is the regression model constant, β1 to βk are the coefficients corresponding to the descriptors X1 to Xk and y is dependent variable. The values for β0 to βk are chosen by minimizing the sum of squares of the vertical distances of the points from the hyperplane so as to give the best prediction of y from X. Regression coefficients represent the independent contributions of each calculated molecular descriptor. In matrix notation, we will write the MLR model is defined in Eq. (4) as:

equation M4

where X(n×k) is of full column rank, including a column of 1 s for the intercept if the intercept is included in the mean function y. We will further assume that we have selected a parameterization for the y so that X has full column rank, meaning that the inverse (XTX)−1 exists; this is not an important limitation on regression models because we can always delete terms from the y, or equivalently delete columns from X, until we have full rank. The k × 1 vector b is the unknown parameter vector. The vector e consists of unobservable errors that we assume are equally variable and uncorrelated, unless stated otherwise. In appropriate Eq. (4), the least squares solution estimate b by equation M5 and the fitted values y corresponding to the experimental RIs are then given by:

equation M6

where H is the n × n matrix defined by H = X(XTX)1XT, which is called the hat matrix because it transforms the vector of experimental responses y into the vector of fitted responses ŷ. The vector of residuals ê is defined by equation M7

The advantages of MLR are that it is simple to use and the derived models are easy to interpret. The sign of the coefficients β0 to βk shows whether the molecular descriptors contribute positively or negatively to the target property and their magnitudes indicates the relative importance of the descriptors to the target property. However, the molecular descriptors should be mathematically independent (orthogonal) of one another and the number of compounds in the training set should exceed the number of molecular descriptors by at least a factor of 5 [45]. Studies have shown that collinear descriptors may result in the coefficients β0 to βk being larger than expected or have the wrong sign [46]. To examine multicollinearity, the variance inflation factor (VIF) was calculated for each variable in the regressions, which is defined as:

equation M8

where Rj2 is the squared correlation coefficient between the jth coefficient regressed against all the other descriptors in the model [47]. The VIF is uniform and equal to 1.00 if there is no linear correlation between a given variable and rest of the variables in the regression equations. Higher values of VIF indicate a more serious multicollinearity problem (popular cut off value is 10). In addition, for inspection of the relative importance and contribution of each descriptors in the QSRR models, the value of mean effect (MF) was calculated for each descriptors by the following equation and it is shown in Table 3:

equation M9

where MFj is the mean effect for considered descriptor j, βj is the coefficient of descriptor j, dij denotes the value of descriptor j of molecule i, m is the number of descriptors in the model and n is the number of molecules in the data sets. The value and sign of the mean effect shows the relative contribution and direction of influence of each descriptor on the RI.

Table 3
The value of mean effect for each descriptor of QSRR models

Results and Discussion

Selecting seven molecular theoretical descriptors including 3 valence connectivity index and 4 solvation connectivity index and the GC-RRI values of 132 solutes on 12 column stationary phases as a dependent one, the best combination of variables is selected by the ES-SWR algorithm to build the relationship between the molecular structure and RRIs using BMLR analysis. The influence of the best number of variables for the subset of all solutes was selected by ES-SWR on the calibration correlation coefficient (R) and square standard error (SE2) for each stationary phase were included in Fig. 2. As can be seen in Fig. 2 for the former 12 GC stationary phases, the R values increases gradually with increasing the number of variables until reaches a plateau while SE2 declines until it drops to a lower limit value. We used the best correlation equation with three optimal variables including solvation connectivity index of order 1 (1Xsol), valence connectivity index of order one (1Xv) and a solvation connectivity index of order 4 (4Xsol) for the analysis of all QSRR models. The descriptors that appear in the BMLR equations for the 12 different stationary phases are identical. A complete list of solutes and the calculated values of the molecular connectivity indexes appearing in the QSRR models are summarized in Table S2 of the supplementary data. These descriptors are related to shape, and the degree of branching of the molecules. This indicates that dispersion interactions and the extent of branching of the molecules affected the retention behavior of organic compounds on the polar and apolar stationary phase columns. The solvation connectivity index (1Xsol) shows an average mean effect of 0.409 for all columns, which is the largest among the descriptors appearing in the QSRR models. This descriptor can be considered as entropy of solvation and somehow indicates the dispersion interactions occurring in the solution. 1Xsol also is a measure of branching of the molecules. The large contributions of this descriptor in the retention behavior of organic molecules is in agreement with the contribution that one would expect for the interaction of nonpolar stationary phases such as C67, C103, C, and C78 with the nonpolar organic molecules. The coefficient of correlation between the GC-RRIs of 132 solutes and 1Xsol index was 0.8339 > 0.8411 > 0.8442 > 0.8463 > 0.8716 > 0.8733 > 0.8736 > 0.8743 > 0.8766 > 0.8828 > 0.8863 > 0.8886, respectively, for the TMO, PCN, TTF, POH, C, PBR, PSH, PCL, MTF, C103, C78 and C67 stationary phases. The presence of 1Xv as a connectivity index with the average mean effect −0.084 in all models indicates that the retention indices depend on the presence and the position of the heteroatoms in the organic molecules. The descriptor shows the negative effect on the retention indices, consequently, the RI decreases with increasing of 1Xv. The solvation connectivity index 4Xsol with the average mean effect 0.039 in all models is a measure of branching of the molecules. Thus, the emergence of the 4Xsol in all QSRR models reflects the influence of the degree of branching on the values of RIs. The positive sign of the corresponding coefficient indicates that the higher ramifications in the solutes, the smaller is 4Xsol and therefore, the bigger is the molecule. This descriptor could also be considered as entropy of solvation and somehow indicates the influence of the dispersion interaction occurring in the stationary phases on the values of RIs. The regression coefficients and the statistical results of the resulted QSRR models for each stationary phase are given in Table 4. The value after the symbol “±” in the parenthesis is the standard deviation related to the regression coefficient. The results in Table 4 indicate that GC-RRIs of 132 solutes on all column stationary phases are strongly dependent on the first order solvation connectivity index (1Xsol) with the maximum mean effect (Table 3) for each stationary phase. The qualities of the models derived from various subsets are evaluated using some statistics, such as the correlation coefficient (R), root mean-square error (RMS), relative error of prediction (REP) and the Fisher’s criterion at the 95% level probability (F) are included in Table 4, also the best fitted equations and results of the root mean square error (RMSE) along with correlation coefficient (R) are plotted in Fig. 3. From Table 4 and Fig. 3, the predicted correlation coefficients over 0.9378, the overall F values higher than 310, and the RMS and REP below 0.0458 and 7.8109, respectively, indicated that the BMLR models have good statistical qualities with low prediction error and demonstrated an excellent predictive power of the obtained QSRR models for all stationary phases. From Table 4, the VIF values are lower than 10 (4.32, 2.54, and 2.34) for three variables 1Xsol, 1Xv, and 4Xsol, respectively, indicating that the QSRR models could have some multicollinearity but it was not serious. According to the statement of Mihalic and Trinajstic [48], the models we have constructed represent good QSRR models judging from the statistics.

Fig. 2
The influences of the number of descriptors on the correlation coefficient (R) and the standard deviation (SE2)
Table 4
QSRR models and statistical parameters of GLC-RRI values versus tetradecane for the total sets (n = 132) of solutes on 12 column stationary phases
Fig. 3
Plot of the root mean square error (♦ RMS and ○ RMScv) and correlation coefficient (♦ R and ○ Rcv) for the resulted 12 QSRR models versus GC stationary phases

Model Prediction-Validation

Model validation is a critical component of QSRR development. A number of procedures have been established to determine the quality of QSRR models. Therefore, a leave-one-out cross-validation (LOO-CV), Y-randomization, and external validation (EV) procedures through an odd–even number and division of the entire dataset into training and test sets are used to validate the predictive ability and check the statistical significance of the developed 12 QSRR models.


The most popular validation method is cross-validation (CV), known as jack-knifing or leave-one-out (LOO). This method systematically removes one data point at a time from the training set, and constructs a model with the reduced dataset. Subsequently, the model is used to predict the data point that has been left out. By repeating the procedure for the entire dataset, a complete set of predicted properties and cross-validated statistics can be obtained. It has been argued that the LOO procedure often overestimates the predictivity of the model and that, subsequently, the QSRR models are overoptimistic [49]. For cross-validated statistics, it has been suggested that prediction residual error sum of squares (PRESS), cross-validated square correlation coefficient (Rcv2) and root mean square error in cross-validation (RMScv) are good estimates of the real prediction error of a model:

equation M10
equation M11
equation M12

where N is the number of training patterns, yobs,i and ypred,i are the experimental, and predicted RRIs of the left-out compound i, respectively and equation M13 is the average experimental RRI of left-in compounds different from i. Values of Rcv2 can range from 1 to <0. A value of one indicates a perfect prediction, and a value of 0 means that the QSRR derived has no modeling power. Negative values arise from a situation where the derived QSRR is a poorer description of data than no model at all. The Rcv2 values can be considered as a measure of the predictive power of a model: whereas R2 can always be increased artificially by adding more parameters, Rcv2 decreases if a model is over parameterized [50], and is therefore a more meaningful summary statistic for predictive models. The correlation coefficients (Rcv) and RMScv for each subset are presented in Table 4 and the resulted values are plotted in Fig. 3. The cross-validation results shows that the Rcv are higher than 0.9325 and RMScv lower than 0.0477 for all GC stationary phases, respectively. Furthermore, in all cases, the cross-validated Rcv values are very close to the corresponding R values and the cross-validated RMScv values are only slightly larger than the corresponding RMS values. Clearly, the cross-validation demonstrates the final models to be statistically significant.

This method is not a very rigorous model predictivity test and suffers from two other major deficiencies: the time to carry out the cross-validation increases as the square of the size of training set; the method produces n final models (each corresponding to one of the training set molecules being left out) and it is not clear which is the ‘best’ model. To further check the prediction ability of the resulting QSRR models two better methods are applied here, one by Hawkins [51] namely as odd–even external validation and the other better method is to remove a percentage of the training set into a prediction set [50, 52].

Odd–Even External Validation

To validate and develop a credible QSRR model, it is not enough to build a model for the whole dataset. So, the 132 dataset solutes for all stationary phases were sorted in the ascending order of GC Kováts RRI values and then divided into two sets namely ‘‘odd set’’ and ‘‘even set’’ RRIs [50, 52]. This way of splitting ensures that the distribution of RRI values of the two subsets were very similar. The QSRR models were fitted to the odd set and even set samples separately and the resulted fitness were assessed by applying QSRR models to both samples. To compare the estimation abilities of the models, two statistical parameters namely root mean squares error (RMSE) and R, were calculated. The same dataset (i.e., ‘calibration set’) that was already used to fit the models was employed to determine resubstitution parameters, i.e. RMSERS and RRS, also to determine holdout parameters, i.e. RMSEHO and RHO for the other dataset, which was not involved in the fitting. The resubstitution statistical parameters of the samples base their predictions on the regression fitted to those samples and this is while the holdout statistical parameters base their predictions on the regression fitted to the other samples. The plots of RRIs estimated by odd- and even-set QSRR models (holdout prediction) versus the RRIs observed experimentally are given in Fig. 4, also Table 5 summarizes these statistical parameters achieved by this approach. As can be seen, in the odd and even-set samples, the resubstitution and holdout RMSE are very similar, indicating that the same sample and other sample predictions are equally precise for all stationary phases.

Fig. 4
Plots of the RRIs estimated for the odd set (♦ A) and even set (○ B) samples by holdout model versus that observed RRIs experimentally for all stationary phases of 132 solutes
Table 5
Statistical parameters of the over-fitting and predictive ability of the models

Y-Randomization Test

Another procedure that is easy to perform is a randomization test called Y-randomization (randomization of response, i.e., in our case RRI). In this method for each column stationary phase, the output RRIs values of the compounds are shuffled randomly, and the resulting dataset is examined by the QSRR method against real (unscrambled) input descriptors to determine the correlation and predictivity of the resulting “model” [5356]. The whole procedure is repeated on many different scrambled datasets. The rationale behind this test is that the significance of the real QSRR model would be suspected if there is a strong correlation between the selected descriptors and the randomized response variables. The randomization was repeated ten times. If the statistical qualities of these models are much lower than the original model, it can be considered that the model is reasonable and had not been obtained by chance. The results are shown in last column of Table 4. Very low level of Rmax (in the interval of 0.0282 for C103 and 0.0551 for MTF) indicates good results in our original models and is not due to a chance correlation or structural dependency of the training set for each stationary phase of the GC column.

Calibration and Prediction Sets

In this investigation, for further testing the predictive ability of the models for the external compounds without the models, part of the congeners are picked up from 132 solutes to construct a training set which is used to develop a prediction model and then predict the values of Kováts RRIs in the remaining congeners. How to pick up the compounds in the training set is very important for developing of the predictive QSRR models. In this case, before each training run, all datasets were split randomly into two separate sub-matrices: the training set matrix and external testing set matrix. Out of 132 organic compounds, 82 solutes (62%) were used for the training set and 50 congeners (38%) were used as external validation. The solutes constituting the training and testing sets are clearly presented in Table S1 of the supplementary data. Moreover, the same divisions were repeated with corresponding RRIs values. The test examples are marked as bold font and training set was also used to obtain the best fit equation of MLR with three molecular descriptors. Furthermore, the testing set was used to monitor overfitting the MLR models. The resulted MLR models for training set congeners were the same as those obtained for the entire set of all solute in each subset subject to use descriptors of all solute models supporting sufficient ability for the prediction set of 50 solutes. The resulting regressions equations of the training set for individual GC column stationary phases with the optimal three molecular descriptors are indexed in Table 6, and results obtained are plotted in Fig. 5. Statistical parameters for the best-fitted models are also presented in Table 7. The correlation coefficients (R) of the obtained models are >0.94 for all the stationary phases and the highest one is 0.9687 for stationary phase C67. The RMS and relative error prediction (REP) of estimation ranged from 0.0331, 5.8983 of C67 stationary phases to 0.0452, 7.6901 of TMO stationary phase, respectively, also the F statistic values are >203.5. The LOO-CV method was used to examine the stability of QSRR models, and the values of Rcv and RMScv for the models were above 0.9359 and in the range of 0.0347 for C67 stationary phase and 0.0496 for TMO stationary phase, respectively. The predicted Kováts RRIs versus the observed Kováts RRIs of the 82 solute training sets are plotted in Fig. 5. As shown in Table 7 and Fig. 5, the QSRR statistical results exhibit good estimation capacity and stability for internal training set solute samples to individual stationary phases. High predictive ability of QSRR models for external examples is another criterion of a good QSRR model. The predicted RRIs of 50 solutes in the external testing set by the models in Table 6 are also demonstrated in Fig. 5 versus the observed RRIs of 12 GC stationary phases. For all 12 GC stationary phases, the regression of the observed and predicted RRIs had a high agreement with the diagonal of each chart. The predicted correlation coefficients (R) over 0.9310 and the RMS REP below 0.0346 and 8.2626, respectively, demonstrated an excellent predictive power of the obtained QSRR models.

Table 6
QSRR models of GC-RRI values versus tetradecane for the training sets (n = 82) of solutes on 12 GC stationary phases
Fig. 5
Plots of the RRIs estimated by the QSPR models in Table 7 versus that observed for 82 training set solutes (♦) and 50 testing set ones (○) for all stationary phases
Table 7
Statistical parameters of the QSRR models obtained using different molecular descriptors


In this study, a novel QSRR tool of BMLR is performed to describe the GC Kováts RRI values of 132 solutes on 12 polar and apolar stationary phase columns based on connectivity molecular descriptors. MLR analysis produced more predictive, informative and significantly improved QSRR models. The use of connectivity indices molecular descriptors revealed to be a completely successful strategy. The effectiveness of the stepwise forward and backward elimination algorithm is demonstrated by the selection of the best set of molecular descriptors. All QSRR models provide a reasonably good calibrated correlation coefficient. The validation and predictive ability of the models were examined by the leave one-out cross-validation, Y-randomization, and external validation. The three methods indicated that the resulting multiparametric QSRR models possess high prediction ability and low overfitting.

Electronic supplementary material

Below is the link to the electronic supplementary material.

Open Access

This article is distributed under the terms of the Creative Commons Attribution Noncommercial License which permits any noncommercial use, distribution, and reproduction in any medium, provided the original author(s) and source are credited.


1. Kaliszan R (2000) In: Valko K (ed) Recent advances in quantitative structure-retention relationships; separation methods in drug synthesis and purification. Elsevier, Amsterdam
2. Kaliszan R (1997) Structure and retention in chromatography. A chemometric approach. Harwood, Amsterdam
3. Karelson M, Lovanov M, Katritzky AR. Chem Rev. 1996;96:1027–1044. doi: 10.1021/cr950202r. [PubMed] [Cross Ref]
4. Kaliszan R. Chem Rev. 2007;107:3212–3246. doi: 10.1021/cr068412z. [PubMed] [Cross Ref]
5. Liu F, Liang Y, Cao C, Zhou N. Anal Chim Acta. 2007;594:279–289. doi: 10.1016/j.aca.2007.05.023. [PubMed] [Cross Ref]
6. Hu RJ, Liu HX, Zhang RS, Xue CX, Yao XJ, Liu MC, Hu ZD, Fan BT. Talanta. 2005;68:31–39. doi: 10.1016/j.talanta.2005.04.034. [PubMed] [Cross Ref]
7. Héberger K. J Chromatogr A. 2007;1158:273–305. doi: 10.1016/j.chroma.2007.03.108. [PubMed] [Cross Ref]
8. Junkes BS, Amboni RDMC, Yunes RA, Heinzen VEF. Anal Chim Acta. 2003;477:29–39. doi: 10.1016/S0003-2670(02)01413-7. [Cross Ref]
9. Ren B. Chemom Intell Lab Syst. 2003;66:29–39. doi: 10.1016/S0169-7439(03)00004-2. [Cross Ref]
10. Estrada E, Gutierrez Y. J Chromatogr A. 1999;858:187–199. doi: 10.1016/S0021-9673(99)00808-0. [PubMed] [Cross Ref]
11. Liu S, Yin C, Cai S, Li Z. Chemom Intell Lab Syst. 2002;61:3–15. doi: 10.1016/S0169-7439(01)00146-0. [Cross Ref]
12. Tello AM, Lebrón-Aguilar R, Quintanilla-López JE, Santiuste JM. J Chromatogr A. 2009;1216:1630–1639. doi: 10.1016/j.chroma.2008.10.025. [PubMed] [Cross Ref]
13. Junkes BS, Amboni RDMC, Yunes RA, Heinzen VEF (2004) J Braz Chem Soc 15:183–189. doi:10.1590/S0103-50532004000200005
14. Garkani-Nejad Z, Karlovits M, Demuth M, Stimpfl T, Vycudilik W, Jalali-Heravi M, Varmuza K. J Chromatogr A. 2004;1028:287–295. doi: 10.1016/j.chroma.2003.12.003. [PubMed] [Cross Ref]
15. Farkas O, Héberger K. J Chem Inf Model. 2005;45:339–346. doi: 10.1021/ci049827t. [PubMed] [Cross Ref]
16. Héberger K, Görgényi M. J Chromatogr A. 1999;845:21–31. doi: 10.1016/S0021-9673(99)00323-4. [Cross Ref]
17. Héberger K, Görgényi M, Sjöström M. Chromatographia. 2000;51:595–600. doi: 10.1007/BF02490818. [Cross Ref]
18. Körtvélyesi T, Görgényi M, Héberger K. Anal Chim Acta. 2001;428:73–82. doi: 10.1016/S0003-2670(00)01220-4. [Cross Ref]
19. Tulasamma P, Reddy KS. Internet Electron J Mol Des. 2007;6:345–362.
20. Wang Y, Yao X, Zhang X, Zhang R, Liu M, Hu Z, Fan B. Talanta. 2002;57:641–652. doi: 10.1016/S0039-9140(02)00078-4. [PubMed] [Cross Ref]
21. Kier LB, Hall LH. Molecular connectivity in structure-activity analysis. New York: Wiley; 1986.
22. Hall LH, Kier LB. J Mol Graph Model. 2001;20:4–18. doi: 10.1016/S1093-3263(01)00097-3. [PubMed] [Cross Ref]
23. Randić M. J Mol Graph Model. 2001;20:19–35. doi: 10.1016/S1093-3263(01)00098-5. [PubMed] [Cross Ref]
24. Todeschini R, Consonni V. Handbook of molecular descriptors. Weinheim, Germany: Wiley-VCH; 2000.
25. Mihalic Z, Nikolic S, Trinajstic N. J Chem Inf Comput Sci. 1992;32:28–37.
26. Estrada E. Chem Phys Lett. 2001;336:248–252. doi: 10.1016/S0009-2614(01)00127-0. [Cross Ref]
27. Randic M. J Chromatogr. 1978;161:1–14. doi: 10.1016/S0021-9673(01)85209-2. [Cross Ref]
28. Sutter JM, Peterson TA, Jurs PC. Anal Chim Acta. 1997;342:113–122. doi: 10.1016/S0003-2670(96)00578-8. [Cross Ref]
29. Liu F, Liang Y, Cao C, Zhou N. Talanta. 2007;72:1307–1315. doi: 10.1016/j.talanta.2007.01.038. [PubMed] [Cross Ref]
30. Katritzky AR, Cheng K, Maran Y, Carlson DA. Anal Chem. 2000;72:101–109. doi: 10.1021/ac990800w. [PubMed] [Cross Ref]
31. Tulasamma P, Reddy KS. J Mol Graph Model. 2006;25:507–513. doi: 10.1016/j.jmgm.2006.04.003. [PubMed] [Cross Ref]
32. Reddy KS, Dutoit JC, Kováts Esz. J Chromatogr. 1992;609:229–259. doi: 10.1016/0021-9673(92)80167-S. [Cross Ref]
33. Reddy KS, Cloux R, Kováts Esz. J Chromatogr A. 1994;673:181–209. doi: 10.1016/0021-9673(94)85037-2. [Cross Ref]
34. Defayes G, Reddy KS, Dallos A, Kováts Esz. J Chromatogr A. 1995;699:131–154. doi: 10.1016/0021-9673(95)00023-G. [Cross Ref]
35. Reddy KS, Cloux R, Kováts Esz. J Chromatogr A. 1995;704:387–436. doi: 10.1016/0021-9673(95)93206-B. [Cross Ref]
36. Laffort P, Chauvin F, Dallos A, Callegari P, Valentin D. J Chromatogr A. 2005;1100:90–107. doi: 10.1016/j.chroma.2005.09.022. [PubMed] [Cross Ref]
37. ISIS Draw 2.3 (1990–2000) MDL Information Systems, Inc
38. HyperChem Release 7.1 for windows molecular modeling system program package, HyperCube, 2002
39. Todeschini R, Consonni V (2002) Dragon software version 2.1.
40. Randic M, Zupan J. J Chem Inf Comput Sci. 2001;41:550–560. [PubMed]
41. Olivero J, Garcia T, Payares P, Vivas R, Diaz D, Daza E, Geerliger P. J Pharm Sci. 1997;86:625–630. doi: 10.1021/js960196u. [PubMed] [Cross Ref]
42. Brown S, Tauler R, Walczak B. Comprehensives chemometrics. Amsterdam: Elsevier; 2009.
43. Afantitis A, Melagraki G, Sarimveis H, Koutentis PA, Markopoulos J, Markopoulou OI. QSAR Comb Sci. 2006;25:928–935. doi: 10.1002/qsar.200530208. [Cross Ref]
44. Rencher AC. Methods of multivariate analysis. 2. Canada: Wiley; 2002.
45. Topliss JG, Edwards RP. J Med Chem. 1979;22:1238–1244. doi: 10.1021/jm00196a017. [PubMed] [Cross Ref]
46. Eriksson L, Jaworska J, Cronin M, Worth A, Gramatica P, McDowell R (2003) Environ Health Perspect 111:1361–75 (PMID: 12896860) [PMC free article] [PubMed]
47. Chatterjee S, Hadi AS, Price B (2006) Analysis of collinear data. In: Regression analysis by example, 4th edn. Wiley, New York, pp 221–258 (Chap. 9)
48. Mihalic M, Trinajstic N. J Chem Educ. 1992;69:701–712. doi: 10.1021/ed069p701. [Cross Ref]
49. Shao J. J Am Stat Assoc. 1993;88:486–494. doi: 10.2307/2290328. [Cross Ref]
50. Hawkins DM, Basak SC, Mills D. J Chem Inf Comput Sci. 2003;43:579–586. [PubMed]
51. Hawkins DM. J Chem Inf Comput Sci. 2004;44:1–12. [PubMed]
52. Tetko IV, Livingstone DJ, Luik AI. J Chem Inf Comput Sci. 1995;35:826–833.
53. Golbraikh A, Tropsha A. Mol Divers. 2000;5:231–243. doi: 10.1023/A:1021372108686. [PubMed] [Cross Ref]
54. Melagraki G, Afantitis A, Sarimveis H, Koutentis PA, Markopoulos J, Igglessi-Markopoulou O. Bioorg Med Chem. 2007;15:7237–7247. doi: 10.1016/j.bmc.2007.08.036. [PubMed] [Cross Ref]
55. Ghavami R, Sadehgi F. Chromatographia. 2009;70:851–868. doi: 10.1365/s10337-009-1233-6. [Cross Ref]
56. Hemmateenejad B, Sanchooli M. J Chemometr. 2007;21:96–107. doi: 10.1002/cem.1039. [Cross Ref]

Articles from Springer Open Choice are provided here courtesy of Springer