 Research article
 Open access
 Published:
Examining the predictive accuracy of the novel 3D Nlinear algebraic molecular codifications on benchmark datasets
Journal of Cheminformatics volume 8, Article number: 10 (2016)
Abstract
Background
Recently, novel 3D alignmentfree molecular descriptors (also known as QuBiLSMIDAS) based on twolinear, threelinear and fourlinear algebraic forms have been introduced. These descriptors codify chemical information for relations between two, three and four atoms by using several (dis)similarity metrics and multimetrics. Several studies aimed at assessing the quality of these novel descriptors have been performed. However, a deeper analysis of their performance is necessary. Therefore, in the present manuscript an assessment and statistical validation of the performance of these novel descriptors in QSAR studies is performed.
Results
To this end, eight molecular datasets (angiotensin converting enzyme, acetylcholinesterase inhibitors, benzodiazepine receptor, cyclooxygenase2 inhibitors, dihydrofolate reductase inhibitors, glycogen phosphorylase b, thermolysin inhibitors, thrombin inhibitors) widely used as benchmarks in the evaluation of several procedures are utilized. Three to nine variable QSAR models based on Multiple Linear Regression are built for each chemical dataset according to the original division into training/test sets. Comparisons with respect to leaveoneout crossvalidation correlation coefficients \(\left( {Q_{loo}^{2} } \right)\) reveal that the models based on QuBiLSMIDAS indices possess superior predictive ability in 7 of the 8 datasets analyzed, outperforming methodologies based on similar or more complex techniques such as: Partial Least Square, Neural Networks, Support Vector Machine and others. On the other hand, superior external correlation coefficients \(\left( {Q_{ext}^{2} } \right)\) are attained in 6 of the 8 test sets considered, confirming the good predictive power of the obtained models. For the \(Q_{ext}^{2}\) values nonparametric statistic tests were performed, which demonstrated that the models based on QuBiLSMIDAS indices have the best global performance and yield significantly better predictions in 11 of the 12 QSAR procedures used in the comparison. Lastly, a study concerning to the performance of the indices according to several conformer generation methods was performed. This demonstrated that the quality of predictions of the QSAR models based on QuBiLSMIDAS indices depend on 3D structure generation method considered, although in this preliminary study the results achieved do not present significant statistical differences among them.
Conclusions
As conclusions it can be stated that the QuBiLSMIDAS indices are suitable for extracting structural information of the molecules and thus, constitute a promissory alternative to build models that contribute to the prediction of pharmacokinetic, pharmacodynamics and toxicological properties on novel compounds.
Background
Computational methods that employ statistical and/or artificial intelligence procedures are widely used in the drug discovery process, where the Quantitative Structure–Activity Relationship (QSAR) studies have an important role [1–4]. These studies are based on the principle that the biological activity (or property) of compounds depends on their structural and physicochemical features and thus, are primarily aimed at finding good correlations among molecular features and specific biological activities [5]. In this way, models with high external predictive ability in novel compounds could be built.
Right from the works developed by Hansch and Fujita in 1960s [6, 7], considered as the origins of the modern QSAR studies [8], several approaches have been reported in the literature with most of these being 2DQSAR methods, that is, they only consider the topological structural features of molecules often using matrix representations such as the connectivity and distance matrices [8]. However, with the introduction of the CoMFA [9] methodology in 1988, the 3DQSAR approaches become popular. These take into account the geometric (3D) features of molecules, which can be computed either from the information represented in a grid through an alignment process with respect to a reference compound or a pharmacophore [2, 10, 11], or using procedures based on Cartesian coordinates [8, 12, 13], molecular spectra [14, 15] and molecular transforms [16], or by the adaptation of 2D methods to take into account threedimensional (3D) aspects [17–21].
However, despite the number and variety of procedures defined up to date, there exists continued interest in creating or extending the current approaches to more generalized forms in order to codify more relevant chemical information with the aim of yielding QSAR models with better predictive ability. This assertion is in accordance with the Non Free Lunch Theorem [22], which could be interpreted as no single QSAR procedure yields superior predictions than all the others when its performance is averaged over all possible compound datasets. This can be confirmed in a report performed by Sutherland et al. [23], where it is observed how wellestablished procedures, assessed in eight diverse chemical datasets, present moderate predictions and without significant differences among them (see Additional file 1: Table S1 for a statistical analysis). The justification for this observation is that one family of molecular descriptors (MDs) may not suffice to codify all chemical information and/or molecular properties for different chemical datasets. In other words, the relevance of MDs depends on the nature of the compounds under study. It is therefore necessary to search for alternative methods/approaches to codify novel and orthogonal chemical information.
Inspired by the previous idea, recently the 3D Nlinear algebraic molecular descriptors have been introduced as a novel mathematical procedure for computing the structural features of chemical compounds [24–26]. These MDs employ the bilinear, quadratic and linear algebraic maps [27] to codify information between atompairs by using several (dis)similarity metrics [25]. Also, the Nlinear algebraic forms [28] were used as generalized expressions of the bilinear, quadratic and linear algebraic maps, when relations among three and four atoms are studied [26]. In this way, the geometric matrix [8] was extended to consider for the first time relations for more than two atoms.
Several studies aimed at assessing the quality of this novel descriptor family, also called QuBiLSMIDAS [acronym of Quadratic, Bilinear and NLinear Maps based on Ntuple Spatial Metric [(Dis)Similarity] Matrices and Atomic Weightings], were performed and these included an evaluation of the information content (variability) and linear independence using Shannon’s entropy based variability analysis [29] (using IMMAN software [30]) and the principal component analysis (PCA) technique [31], respectively. Also, comparisons with other MDs reported in the literature were performed [25, 26]. In general sense, the results demonstrated that the novel MDs have superior variability than 3D DRAGON indices and another approaches implemented in several software [32–35]. Furthermore, the results revealed that the novel 3D Nlinear indices not only do they codify all information contained in the 3D DRAGON MDs, but capture information orthogonal to the latter. Lastly, the QuBiLSMIDAS MDs were used for modeling the binding affinity to the corticosteroidbinding globulin (CBG), achieving superior results with respect to other QSAR methodologies (see Tables 8–9 in Ref. [25] and Tables 9–10 in Ref. [26]).
However, although the initial results with QuBiLSMIDAS MDs are promissory, it cannot be stated that these are most suitable for building QSAR models for all chemical datasets. It is thus necessary to evaluate the performance of the 3D Nlinear algebraic MDs in QSAR modeling with different molecular sets. Therefore, this paper is dedicated to the assessment of the utility of the QuBiLSMIDAS approach in the prediction of the biological activity in several compound datasets and the comparison of the obtained results with those of other QSAR procedures reported in the literature.
Mathematical overview of the 3D Nlinear algebraic molecular descriptors
In this report, the total and localfragment 3D Nlinear Algebraic indices [25, 26] (also known as QuBiLSMIDAS) are employed to assess the predictive accuracy of this approach in QSAR studies. These molecular descriptors (MDs) are calculated from the contribution of each atom in a molecule. That is, if a molecule is comprised of n atoms then the kth twolinear, threelinear and fourlinear algebraic indices for each atom “a” are computed as Nlinear (Multilinear) algebraic forms (maps) in \({\mathbb{R}}^{n}\), in a canonical basis set, when relations among two (N = 2), three (N = 3) and four (N = 4) atoms are considered, respectively. These descriptors are mathematically expressed as follows:
where, “a” is a specific atom (a = 1, 2,…,n), n is the number of atoms in a molecule, _{(F)} L _{ a } is the entry corresponding to the contribution of the atom “a” in the vector of atomlevel indices _{(F)} L, F is a localfragment (group or atomtype) that may or not be considered in the index computation, and x ^{1},…,x ^{n}, y ^{1},…,y ^{n}, z ^{1},…,z ^{n} and w ^{1},…,w ^{n} are the values (coordinates or components) of the molecular vectors \(\bar{x}\), \(\bar{y}\), \(\bar{z}\) and \(\bar{w}\), respectively. In addition, the coefficients \(g_{ij(F)}^{a,k}\), \(gt_{ijl(F)}^{a,k}\) and \(gq_{ijlh(F)}^{a,k}\) are the elements of the kth twotuple, threetuple and fourtuple atom  level total (or local  fragment) spatial(dis)similarity matrices [\({\mathbb{G}}_{\left( F \right)}^{a,k}\), \({{\mathbb{G}}{\mathbb{T}}}_{\left( F \right)}^{a,k}\) and \({{\mathbb{G}}{\mathbb{Q}}}_{\left( F \right)}^{a,k}\)], which are obtained from the corresponding kth twotuple, threetuple and fourtuple total (or local  fragment) spatial(dis)similarity matrices [\({\mathbb{G}}_{\left( F \right)}^{k}\), \({{\mathbb{G}}{\mathbb{T}}}_{\left( F \right)}^{k}\) and \({{\mathbb{G}}{\mathbb{Q}}}_{\left( F \right)}^{k}\)]. Lastly, k (±1, ±2,…,±12) is the power to which the matrix approaches are raised through the Hadamard product.
The molecular vectors (or property vectors) \(\bar{x}\), \(\bar{y}\), \(\bar{z}\) and \(\bar{w}\) are calculated by using the Chemistry Development Kit (CDK) library [36] considering the following fragment and atombased properties: atomic mass (m), the van der Waals volume (v), the atomic polarizability (p), atomic electronegativity in Pauling scale (e), atomic GhoseCrippen LogP (a), GasteigerMarsili atomic charge (c), atomic polar surface area (psa), atomic refractivity (r), atomic hardness (h) and atomic softness (s).
The total matrix approaches \({\mathbb{G}}_{{}}^{k}\), \({{\mathbb{G}}{\mathbb{T}}}_{{}}^{k}\) and \({{\mathbb{G}}{\mathbb{Q}}}_{{}}^{k}\) constitute the basis for the calculation of the twolinear, threelinear and fourlinear indices and these are employed to represent the chemical information codified on interactions among “N” atoms of a molecule. Specifically, for k = 1 (matrix of order 1) the coefficients \(g_{ij}^{1}\), \(gt_{ijl}^{1}\) and \(gq_{ijlh}^{1}\) corresponding to the matrices \({\mathbb{G}}_{{}}^{1}\), \({{\mathbb{G}}{\mathbb{T}}}_{{}}^{1}\) and \({{\mathbb{G}}{\mathbb{Q}}}_{{}}^{1}\) can be calculated by using several (dis)similarity metrics and multimetrics to capture the information on the relations between two, three and four atoms, respectively [25, 26]. To compute the atompair relations, metrics (see Table 1) derived from the general Minkowski definition (e.g. Manhattan, Euclidean) as well as others that have been successfully used in machine learning algorithms and similarity/dissimilarity studies (e.g. Canberra, Soergel, Clark) are employed. On the other hand, different multimetrics (see Table 2) to calculate the ternary (three) and quaternary (four) relations among atoms of a molecule can be utilized, such as: bond angle for relations among three atoms and dihedral angle for relations among four atoms. Table 3 shows examples of twotuple and threetuple total spatial(dis)similarity matrices calculated with some previously mentioned metrics and multimetrics.
From these total matrix approaches (\({\mathbb{G}}^{k}\), \({{\mathbb{G}}{\mathbb{T}}}^{k}\) and \({{\mathbb{G}}{\mathbb{Q}}}^{k}\)), localfragments matrices may be computed in order to consider atomtypes or chemical regions of interest and thus yielding the kth twotuple, threetuple and fourtuple local  fragment spatial(dis)similarity matrices, denoted by \({\mathbb{G}}_{F}^{k}\), \({{\mathbb{G}}{\mathbb{T}}}_{F}^{k}\) and \({{\mathbb{G}}{\mathbb{Q}}}_{F}^{k}\), respectively (see Eq. 13 in Ref. [25] and Eqs. 17–18 in Ref. [26]). Specifically, the localfragments (or atomtypes), F, that could be taken into account to compute these indices include: hydrogen bond acceptors (A), carbon atoms in aliphatic chains (C), hydrogen bond donors (D), halogens (G), terminal methyl groups (M), carbon atoms in aromatic portion (P) and heteroatoms (X) (see Table 4 for examples).
These total (or localfragment) matrix approaches (\({\mathbb{G}}_{(F)}^{k}\), \({{\mathbb{G}}{\mathbb{T}}}_{(F)}^{k}\) and \({{\mathbb{G}}{\mathbb{Q}}}_{(F)}^{k}\)) are also known as kth non  stochastic twotuple, threetuple and fourtuple total (or localfragment) spatial(dis)similarity matrices denoted by \(_{ns} {\mathbb{G}}_{(F)}^{k}\), \(_{ns} {{\mathbb{G}}{\mathbb{T}}}_{(F)}^{k}\) and \(_{ns} {{\mathbb{G}}{\mathbb{Q}}}_{(F)}^{k}\), respectively, because no normalizing procedure is used in their computation. Nonetheless, with the purpose of obtaining normalized matrix representations three probabilistic schemes may be employed to compute the QuBiLSMIDAS MDs. In this way, the following normalized matrix representations are obtained from the corresponding nonstochastic matrices: the kth simple  stochastic twotuple, threetuple and fourtuple total (or localfragment) spatial(dis)similarity matrices [\(_{ss} {\mathbb{G}}_{(F)}^{k}\), \(_{ss} {{\mathbb{G}}{\mathbb{T}}}_{(F)}^{k}\) and \(_{ss} {{\mathbb{G}}{\mathbb{Q}}}_{(F)}^{k}\)] (see Eq. 10 in Ref. [25] and Eqs. 13–14 in Ref. [26]), the kth double  stochastic twotuple total (or localfragment) spatial(dis)similarity matrix [\(_{ds} {\mathbb{G}}_{(F)}^{k}\)] (see Sinkhorn–Knopp algorithm in Ref. [37]) and the kth mutual probability twotuple, threetuple and fourtuple total (or localfragment) spatial(dis)similarity matrices [\(_{mp} {\mathbb{G}}_{(F)}^{k}\), \(_{mp} {{\mathbb{G}}{\mathbb{T}}}_{(F)}^{k}\) and \(_{mp} {{\mathbb{G}}{\mathbb{Q}}}_{(F)}^{k}\)] (see Eq. 12 in Ref. [25] and Eqs. 15–16 in Ref. [26]). Table 5 shows the results obtained with the three probabilistic transformations on a twotuple total spatial(dis)similarity matrix.
Finally, from the nonstochastic (simplestochastic, doublestochastic or mutualprobability) total (or localfragment) matrices [i.e. \({\mathbb{G}}_{(F)}^{k}\), \({{\mathbb{G}}{\mathbb{T}}}_{(F)}^{k}\) and \({{\mathbb{G}}{\mathbb{Q}}}_{(F)}^{k}\)], the corresponding atomlevel matrices [denoted as \({\mathbb{G}}_{\left( F \right)}^{a,k}\), \({{\mathbb{G}}{\mathbb{T}}}_{\left( F \right)}^{a,k}\) and \({{\mathbb{G}}{\mathbb{Q}}}_{\left( F \right)}^{a,k}\), respectively] are calculated and their coefficients are used in the descriptors calculation (see Eqs. 1–3). Each atomlevel matrix determines an atomlevel index for atom “a” of a molecule and this value constitutes a component (entry) of the vector _{(F)} L. Once the vector _{(F)} L is computed then the global definition of the kth twolinear, threelinear and fourlinear algebraic indices is obtained by applying over the entries of _{(F)} L one or several aggregation operators (see Additional file 1: Table S2 for mathematical definition) [25, 26], which have been successfully employed in other reports [38–40]. In the Scheme 1 a general flowchart regarding the calculation process of the QuBiLS MIDAS MDs detailed in this section may be observed, while Scheme 2 is a graphic representation of each step employed in the computation of a specific twolinear algebraic index.
In order to automatize the calculation of the 3D Nlinear algebraic indices used in the present manuscript the QuBiLSMIDAS software has been developed [41]. This software has as one of its main features the multicore processing of the MDs, as well as the option to carry out the distributed calculation of the indices by using the MultiServer Distributed Computing Platform known as Tarenal [42]. The latter is particularly useful for highthroughput calculation tasks. Both software are freely available via internet at: http://tomocomd.com/.
Methods
In order to assess the correlation ability of the QuBiLSMIDAS MDs for different biological activities eight wellknown chemical datasets were used. These were previously employed by Sutherland et al. in a comparative study of QSAR methods commonly used in chemoinformatics analysis [23] and since then, these have been utilized as “benchmarks” for comparing results obtained in other approaches [43–47]. These datasets are comprised by angiotensin converting enzyme (ACE) inhibitors, acetylcholinesterase (AchE) inhibitors, ligands for the benzodiazepine receptor (BZR), cyclooxygenase2 (COX2) inhibitors, dihydrofolate reductase inhibitors (DHFR), inhibitors of glycogen phosphorylase b (GPB), thermolysin inhibitors (THER) and thrombin inhibitors (THR). In this study the 3D coordinates were generated using CORINA software, and the same partitioning into training and test sets used in the initial study was considered in order to guarantee comparability of results.
For these datasets, several configurations based on 3D twolinear, threelinear and fourlinear algebraic indices were computed (see Additional file 1: Table S3) using the QuBiLS MIDAS software [41]. Due to the fact that numerous MDs are computed with this program yielding a highdimensional space, then strategies for data reduction are necessary. In this sense, the following workflow for each set of indices calculated for each chemical dataset was performed only considering the training set compounds:

The 1000 MDs with best variability behavior according to their Shannon’s Entropy values [29] were retained by using the IMMAN software [30].

The MDs with values represented as power of 10 (scientific notation) and whose exponents are greater or lesser than ±5 were removed.

Filters for removing the MDs with correlation equal or greater than 0.95 and standardized entropy lesser than 0.3 were applied.

The statistical method Multiple Linear Regression (MLR) implemented in the STATISTICA software was employed in order to select the MDs included in the model by using Forward Stepwise and Backward Stepwise selection procedures.

The MDs retained after applying the previous steps and computed for the same compounds were merged into a single dataset.
With the reduced data matrices for each chemical datasets, QSAR models were built with the MLR technique to determine the relationship between the response (activity) and predictor variables (MDs). The MLR technique is coupled with the Genetic Algorithm (GA) metaheuristic as the variable selection method [48]. This strategy (MLR + GA) is implemented in the MobyDigs software (version 1.0) which was utilized to carry out this study [49]. In this sense, to perform the search process several populations with 100 3D Nlinear MDs each were created, while the following configurations were used for the GA procedure: Number of iterations equal to 500,000, Population size equal to 100, Reproduction/mutation tradeoff equal to 0.5, Selection bias was initially set to 0 (indicative of random selection) until achieving the 80 % of the maximum number of iterations and was later set to 1 (indicates tournament selection) in order to increase the selection pressure. The values of the previous parameters were selected according to the study performed by Todeschini et al. in Ref. [49].
The search process was carried out by using the \(Q_{loo}^{2}\) (“leaveoneout” cross validation) statistical parameter as the fitness function. Once the exploration in each population was completed, then the MDs included in the built 9variable models were retained with the purpose of creating new populations until 100 MDs. This process is repeated until achieving an only one population with 100 MDs as maximum. Finally, from the final population and for each compound dataset, 3–9 variable regression models were built for the corresponding biological activity. However, as the MobyDigs software generates a set of MLR models then the choice of the most suitable model was performed according to the following steps:

The “best” 50 models according to the \(Q_{loo}^{2}\) parameter were retained.

To each model retained the validation methods “bootstrapping” [50] \(\left( {Q_{boot}^{2} } \right)\) and “Yscrambling” [51] (a(Q ^{2})) were applied in order to assess the predictive power and the possible chance correlation with respect to the modeled biological activity, respectively. The former randomly creates training sets (with repeated objects) of the same size as the original and the objects left out constitute the test set, while the latter randomly changes the true response variables to determine the quality of the model. Both procedures were repeated 5000 and 300 times, respectively. These methods were applied due to the fact that \(Q_{loo}^{2}\) procedure does not suffice to validate the stability of a predictive model [52].

For each model the function \(f(x) = \left( {1  Q_{boot}^{2} } \right) + \left {a\left( {Q^{2} } \right)} \right\) was computed, which takes into account the results obtained with the two validation procedures employed and the model with the smallest f(x) value constitutes the “best” regression model.

The “best” regression model was assessed by using “external validation” \(\left( {Q_{ext}^{2} } \right)\) procedure in the corresponding test set in order to measure its generalization ability.
Results and discussion
Assessment of the QuBiLSMIDAS models versus other approaches
In this section the performance of the QuBiLSMIDAS models for the chemical datasets described in section “Methods” is compared with respect to 16 QSAR methodologies (or descriptor sets) reported in the literature. The Table 6 shows the statistical parameters and equations of the best regression model based on total and localfragment QuBiLSMIDAS indices corresponding to each chemical dataset used in this report. In general sense, it can be observed that the bootstrapping validation coefficient \(\left( {Q_{boot}^{2} } \right)\) calculated for each model presents a value greater than 0.6, indicative of the good predictive power of the built models. Also, the coefficients computed from scrambling tests (a(Q ^{2})) have in all cases values inferior to 0.4, indicating reduced propensity to chance correlation. Lastly, the values achieved in the external prediction \(\left( {Q_{ext}^{2} } \right)\) suggest that the models based on QuBiLSMIDAS MDs have appropriate generalization ability, given that all \(Q_{ext}^{2}\) parameters present values superior to 49 % of the total variance even when outlier compounds are retained in the validation set.
On the other hand, the Tables 7 and 8 show the comparisons with respect to other approaches reported in the literature, as well as the results obtained by the models based on total QuBiLSMIDAS MDs exclusively (see Additional file 1: Table S4 for information related with the best models from 3 to 9 variables). In this manner, the importance of considering localfragments (atomtypes or group) in the calculation of the QuBiLSMIDAS MDs and subsequently in the building of QSAR models can be analyzed. As can be observed in both tables, the performance of the QuBiLSMIDAS models is superior when localfragments are considered with respect to those QuBiLSMIDAS models that do not use them. Particularly, it can be noted that in 6 of the 8 datasets studied the \(Q_{loo}^{2}\) parameter is rather comparable, while better performances are attained according to \(Q_{ext}^{2}\). Both parameters for the COX2 dataset present the best improvements, achieving in the external prediction a value greater than 49 % of the total variance, while no other QSAR procedure outperforms this threshold. On the other hand, only in the DHFR and GPB datasets does the utilization of the localfragment QuBiLSMIDAS MDs not influence the performance of the developed QSAR models. It can thus be stated that considering a mixture of total and localfragment QuBiLSMIDAS MDs in building of QSAR models contributes to the improvement of the predictive ability.
Also, it can be observed from Table 7 that the crossvalidation performances achieved by the QuBiLSMIDAS models have comparabletosuperior behavior with respect to the approaches reported in the literature. Until now, the best \(Q_{loo}^{2}\) value for the datasets ACE, ACHE, BZR, COX2, GPB, THER and THR had been attained by the procedures HQSAR (and 2.5D) [\({\text{Q}}_{\text{loo}}^{2}\) = 0.72], SAMFARF (\({\text{Q}}_{\text{loo}}^{2}\) = 0.58), AllShortest Path [ASP] Fingerprint (\({\text{Q}}_{\text{loo}}^{2}\) = 0.56), COMSIA extra (\({\text{Q}}_{\text{loo}}^{2}\) = 0.57), HQSAR (and SAMFARF) [\({\text{Q}}_{\text{loo}}^{2}\) = 0.66], O3A/O3Q (\({\text{Q}}_{\text{loo}}^{2}\) = 0.67) and COMSIA extra (\({\text{Q}}_{\text{loo}}^{2}\) = 0.72), respectively, by using PLS, Random Forest (RF) or Support Vector Machine (SVM) techniques. However, all these previous results are clearly outperformed by the QuBiLSMIDAS models [(ACE, \({\text{Q}}_{\text{loo}}^{2}\) = 0.7756), (ACHE, \({\text{Q}}_{\text{loo}}^{2}\) = 0.6574), (BZR, \({\text{Q}}_{\text{loo}}^{2}\) = 0.6931), (COX2, \({\text{Q}}_{\text{loo}}^{2}\) = 0.6313), (GPB, \({\text{Q}}_{\text{loo}}^{2}\) = 0.8124), (THER, \({\text{Q}}_{\text{loo}}^{2}\) = 0.7530) and (THR, \({\text{Q}}_{\text{loo}}^{2}\) = 0.8149)], which were built with MLR that is a simpler method than those employed in the reported results. In the specific case of the DHFR dataset, although the attained value (\({\text{Q}}_{\text{loo}}^{2}\) = 0.7055) with the QuBiLSMIDAS approach is not better than the current best result (ASP fingerprint, \({\text{Q}}_{\text{loo}}^{2}\) = 0.76), the former is superior to the remaining QSAR procedures. However, it is important to remark that the best model (ASP fingerprint + SVM) for the DHFR dataset does not have the external prediction value (\(Q_{ext}^{2}\)) reported and thus the corresponding \(Q_{loo}^{2}\) could be overoptimistic.
According to the external predictions, it can be observed in the Table 8 that the models based on QuBiLSMIDAS indices yield comparabletosuperior performances with respect to the results reported in the literature. Specifically, the models for ACE (\({\text{Q}}_{\text{ext}}^{2}\) = 0.7422), BZR (\({\text{Q}}_{\text{ext}}^{2}\) = 0.5692), COX2 (\({\text{Q}}_{\text{ext}}^{2}\) = 0.4932), GPB (\({\text{Q}}_{\text{ext}}^{2}\) = 0.8283), THER (\({\text{Q}}_{\text{ext}}^{2}\) = 0.7248) and THR (\({\text{Q}}_{\text{ext}}^{2}\) = 0.7674) test sets outperform the best results reported up to date for each dataset previously mentioned, which correspond to COSMOsar3D (\({\text{Q}}_{\text{ext}}^{2}\) = 0.43) in COX2 and to the 2DFPT methodology in the other datasets [(ACE, \({\text{Q}}_{\text{ext}}^{2}\) = 0.713), (BZR, \({\text{Q}}_{\text{ext}}^{2}\) = 0.378), (GPB, \({\text{Q}}_{\text{ext}}^{2}\) = 0.667), (THER, \({\text{Q}}_{\text{ext}}^{2}\) = 0.649) and (THR, \({\text{Q}}_{\text{ext}}^{2}\) = 0.737)]. The 2DFPT models were developed by using SQS framework that determines linear and nonlinear models (see Table 8), while the model corresponding to COSMOsar3D is based on the PLS technique. Even so, the obtained MLR models have better predictive accuracy, even when these are compared with respect to more complex or similar procedures.
As for the ACHE and DHFR datasets, the predictive power obtained for models built with the QuBiLSMIDAS approach is inferior to the best results reported so far in the literature. In the former dataset, the methods 2DFPT (\({\text{Q}}_{\text{ext}}^{2}\) = 0.714), O3Q (\({\text{Q}}_{\text{ext}}^{2}\) = 0.67) and O3A/O3Q (\({\text{Q}}_{\text{ext}}^{2}\) = 0.65) offer better predictions than the proposed model (\({\text{Q}}_{\text{ext}}^{2}\) = 0.6309), albeit this can be considered as suitable (explains 63 % of total variance). Additionally, when the DHFR test set is taken into account the 2DFPT approach (\({\text{Q}}_{\text{ext}}^{2}\) = 0.683) has more predictive ability than the corresponding QuBiLSMIDAS model (\({\text{Q}}_{\text{ext}}^{2}\) = 0.6405), but the latter is superior to the remaining methodologies. Nonetheless, it is important to highlight that the procedures O3Q and O3A/O3Q are alignment dependent and thus their use is generally restricted to congeneric datasets [45]. In the specific case of the 2DFPT methodology for ACHE and DHFR datasets, the achieved results are based on nonlinear models while the proposed outcomes are determined with linear models.
The obtained results evidence that the QuBiLSMIDAS MDs properly codify structural information of the molecules considering interactions among N (N = 2, 3, 4) atoms and thus are suitable for developing QSAR models that contribute to the prediction of biological activity in novel structures. However, notwithstanding the comparabletosuperior predictions achieved by the proposed models, it is important to statistically validate these results.
Statistical analysis of the external predictive accuracy
To perform this analysis the values corresponding to the external predictions (\({\text{Q}}_{\text{ext}}^{2}\)) obtained by the QuBiLSMIDAS models were taken into consideration as well as the ones reported in the literature over the external compounds belonging to each dataset (see Table 8). Firstly, a descriptive analysis through boxplot graphics was performed (with SPSS software) and the obtained results are represented in Fig. 1. As can be observed, the QuBiLSMIDAS and 2DFTP models tend to have a similar behavior and superior to the remaining procedures. Also, it can be noted that the highest prediction among all procedures analyzed is achieved by the QuBiLSMIDAS models. In addition, taking into account the graphics corresponding to the QuBiLSMIDAS and 2DFPT approaches, it can be concluded that the predictions obtained by the former are less scattered than those attained by the latter and thus, the QuBiLSMIDAS models have a more suitable external predictive ability irrespective of the chemical dataset analyzed. However, these results are not enough to state that the models based on QuBiLSMIDAS MDs are statistically the best.
Therefore, an exploratory study was performed to analyze the normality of the data by using Kolmogorov–Smirnov (K–S) test corrected by Lilliefors [53] and the Shapiro–Wilk test [54]. This was done in order to guarantee that the variable \({\text{Q}}_{\text{ext}}^{2}\) is not normally distributed, at least for one model, and so to ensure that the nonparametric tests are the proper choice. As can be observed in Additional file 1: Table S5, the null hypotheses of normality can only be rejected with a high certainty for \({\text{Q}}_{\text{ext}}^{2}\) values in the 2DFTP and COSMOsar3D models, although with Shapiro–Wilk test the rejection of the null hypothesis is achieved for COMSIA basic as well. Therefore the nonparametric tests may be considered as suitable for this statistical analysis.
Subsequently, a Friedman test [55] for multiple comparisons was performed taking into consideration the results of all QSAR procedures. As can be seen in Additional file 1: Table S6A, there are global differences among the considered methods, with the QuBiLSMIDAS models being those with the best performance followed by the 2DFPT, O3Q and COSMOsar3D approaches, respectively, with a Kendall’s W [56] concordance level of 0.607 (see Additional file 1: Table S6B). In order to determine the specific statistical differences a Wilcoxon signedranks test [57] was carried out (with R software) by using Benjamini and Hochberg [58] (BH) as the adjustment method (onetailed p values calculation) for controlling the false discovery rate (FDR). The results of this analysis are shown in Table 9, where a significant p value (p value <0.05) means that the row approach is superior to the corresponding column. So, it can be noted that the QuBiLSMIDAS models yield statistically better predictions than the other methodologies considered, with the exception of the 2DFPT approach.
Analysis of the predictive ability according to conformer generation methods
The conformer generation constitutes an important step when chemoinformatics tasks are performed, particularly in the computeraided drug design, where the outcomes of a virtual screening process may depend on 3D structures employed to build the procedure to be used, e.g. a QSAR model [59]. Therefore, in this section an evaluation of the sensibility of the QuBiLSMIDAS MDs to the different conformer generation methods is performed in order to comprehend how these could affect in the performance of the indices. To this end, the software FROG2 [60], RDKit [61], BALLOON [62], OpenBabel [63] and Standardizer ChemAxon [64] were employed to generate the 3D structures, taking as starting point the SMILES representations corresponding to the eight compound datasets considered in this report.
Firstly, a study with the purpose of knowing if the models developed using the training structures generated with CORINA (see Table 6) are applicable to the test structures generated with the previously mentioned programs was performed. The external predictive abilities obtained after performing this study are graphically represented in Fig. 2. These results are significantly inferior to those achieved with the test sets based on CORINA (see Additional file 1: Table S8), with the exception of RDKIT. This demonstrates that QSAR models based on QuBiLSMIDAS MDs are not suitable to predict biological activity into compounds optimized with other procedure different from than used for the training structures. Thus, it can be stated that the performance of the QuBiLSMIDAS MDs depend on 3D conformations from which are computed.
It is important to highlight that the previous results do not mean that CORINA software is the most suitable to generate the 3D structures to be used in the development of the QSAR models based on QuBiLSMIDAS MDs. In this sense, in order to prove this assertion the following simple workflow was carried out considering the conformations generated by each previously mentioned program (including CORINA) for each chemical dataset:

8640 twolinear algebraic indices (Additional file 1: Table S9) were computed.

CfsSubsetEval feature selection procedure, implemented in WEKA software, was applied in order to retain those MDs with high correlation according to dependentvariable and with low intercorrelation among them.

The MLRGA procedure implemented in MobyDigs software was employed to build 9variable models performing 100,000 iterations and considering the tabu list options of removing MDs with correlation equal or greater than 0.95, fourth order moment greater than 8 and standardized entropy lesser than 0.3. The fitness function used was the statistical parameter \(Q_{loo}^{2}\).

The model with the highest \(Q_{loo}^{2}\) value was selected as the best model, to which the external predictive ability was determined.
Table 10 shows the external predictive power of the models developed from different 3D conformations, as well as the average of the rankings corresponding to the conformer generation methods considered in this study. As can be observed, the best predictions are achieved by the models built from 3D molecular structures generated by FROG2 procedure, followed by the results obtained from the methods CORINA, CHEMAXON, RDKIT, OPENBABEL and BALLOON, respectively. However, in Additional file 1: Table S10 is demonstrated through a Friedman test that there exists no global statistic differences among previous predictions, which proves, at least for this preliminary study that with QuBiLSMIDAS MDs can be developed QSAR models with good predictive accuracy irrespective of the procedure used to obtain optimized structures.
Note that for the forthcoming version of QuBiLSMIDAS software, RDKIT program will be incorporated in the QuBiLSMIDAS software as a builtin option for conformer generation. This is due to the fact that FROG2 procedure can only be accessed using a web browser, while CORINA and CHEMAXON software are not freely available for use. In addition, according to a study performed in Ref. [65] in order to assess the quality of the conformations generated by several free methods, RDKIT tends to generate the most similar conformations to the experimental structures, in addition to being the second fastest among all toolkits analyzed.
Conclusions
In this report the predictive accuracy of the novel alignmentfree geometric molecular descriptors based on Nlinear algebraic maps (so called QuBiLSMIDAS) has been examined. To this end, QSAR models for predicting the biological activity in eight molecular datasets were developed by using MLR as statistical technique. The results obtained with the QuBiLSMIDAS models were compared with respect to several QSAR procedures reported in the literature according to the correlation coefficients achieved with the leaveoneout crossvalidation \(\left( {Q_{loo}^{2} } \right)\) and external prediction \(\left( {Q_{ext}^{2} } \right)\) methods, and generally superior performances were observed with this QuBiLSMIDAS framework.
A few exceptions were observed: for the \(Q_{loo}^{2}\) parameter, the QuBiLSMIDAS approach is exclusively outperformed by the ASPbased (fingerprint) method in the DHFR dataset, while for the \(Q_{ext}^{2}\) parameter, the QuBiLSMIDAS method yields inferior results with respect to the 2DFPT methodology in the DHFR and ACHE test set, respectively. Also, inferior \(Q_{ext}^{2}\) values are yielded by the QuBiLSMIDAS approach with respect to the O3Q and O3A/O3Q procedures in the ACHE test set. However, these previous methodologies are based on techniques more complex than MLR and/or cannot be used in noncongeneric datasets because are alignmentdepend. Thus, considering the maximum parsimony principle (Ockham’s razor), the QuBiLSMIDAS approach seems to be more suitable than the other QSAR methods.
Additionally, several steps for statistically validating the obtained results are detailed. In this sense, the external predictive ability of the developed models was compared with respect to other methodologies by means of the multiple comparison tests. It was demonstrated that the QuBiLSMIDAS models yield the best predictions, and that these are significantly superior in 11 of the 12 methodologies compared. Therefore, it can be suggested that the 3D Algebraic Nlinear molecular descriptors (also known as QuBiLSMIDAS) are suitable for extracting structural information of the molecules and thus, constitute a promissory alternative to build models that contribute to the prediction of pharmacokinetic, pharmacodynamics and toxicological properties of novel compounds.
References
Norinder U (1996) Single and domain mode variable selection in 3D QSAR applications. J Chemom 10(2):95–105
SungSau S, Karplus M (1997) Threedimensional quantitative structure–activity relationships from molecular similarity matrices and genetic neural networks. 2. Applications. J Med Chem 40(26):4360–4371
AiresdeSousa J, Gasteiger J (2002) Prediction of enantiomeric selectivity in chromatography: application of conformationdependent and conformationindependent descriptors of molecular chirality. J Mol Graph Model 20(5):373–388
Chen H, Zhou J, Xie G (1998) PARM: a genetic evolved algorithm to predict bioactivity. J Chem Inf Comput Sci 38(2):243–250
Kubinyi H (1997) QSAR and 3D QSAR in drug design: 1. Methodology. Drug Discov Today 2(11):457–467
Fujita T, Iwasa J, Hansch C (1964) A new substituent constant, π, derived from partition coefficients. J Am Chem Soc 86(23):5175–5180
Hansch C et al (1962) Correlation of biological activity of phenoxyacetic acids with Hammett substituent constants and partition coefficients. Nature 194(4824):178–180
Todeschini R, Consonni V (2009) Molecular descriptors for chemoinformatics. In: Mannhold R, Kubinyi H, Folkers G (eds) Methods and principles in medicinal chemistry, 2nd edn. WileyVCH, Weinheim
Cramer RD, Patterson DE, Bunce JD (1988) Comparative molecular field analysis (CoMFA). 1. Effect of shape on binding of steroids to carrier proteins. J Am Chem Soc 110(18):5959–5967
Parretti MF et al (1997) Alignment of molecules by the Monte Carlo optimization of molecular similarity indices. J Comput Chem 18(11):1344–1353
Tominaga Y, Fujiwara I (1997) Novel 3D descriptors using excluded volume: application to 3D quantitative structure–activity relationships. J Chem Inf Comput Sci 37(6):1158–1161
Todeschini R, Lasagni M, Marengo E (1994) New molecular descriptors for 2D and 3D structures. Theory. J Chemom 8(4):263–272
Consonni V, Todeschini R, Pavan M (2002) Structure/response correlations and similarity/diversity analysis by GETAWAY descriptors. Part 1. Theory of the novel 3D molecular descriptors. J Chem Inf Comput Sci 42(3):682–692
Bursi R et al (1999) Comparative spectra analysis (CoSA): spectra as threedimensional molecular descriptors for the prediction of biological activities. J Chem Inf Comput Sci 39(5):861–867
Turner DB et al (1999) Evaluation of a novel molecular vibrationbased descriptor (EVA) for QSAR studies: 2. Model validation using a benchmark steroid dataset. J Comput Aided Mol Des 13(3):271–296
Gasteiger G et al (1996) Chemical information in 3D space. J Chem Inf Comput Sci 36(5):1030–1037
Balaban AT (1997) From chemical topology to threedimensional geometry. Springer, New York
Bogdanov B, Nikolic S, Trinajstic N (1990) On the threedimensional Wiener number: a comment. J Math Chem 5(3):305–306
Mekenyan O et al (1986) Modelling the interaction of small organic molecules with biomacromolecules. I. Interaction of substituted pyridines with anti3azopyridine antibody. Arzneim Forsch 36(2):176–183
Randić M (1995) Molecular profiles novel geometrydependent molecular descriptors. New J Chem 19:781–791
Pearlman RS, Smith KM (1998) Novel software tools for chemical diversity. In: Kubinyi H, Folkers G, Martin YC (eds) 3D QSAR in drug design. Kluwer/ESCOM, Dordrecht, pp 339–353
Wolpert DH, Macready WG (1997) No free lunch theorems for optimization. IEEE Trans Evolut Comput 1(1):67–82
Sutherland JJ, O’Brien LA, Weaver DF (2004) A comparison of methods for modeling quantitative structure–activity relationships. J Med Chem 47(22):5541–5554
Cubillán N et al (2015) Novel global and local 3D atombased linear descriptors of the Minkowski distance matrix: theory, diversity–variability analysis and QSPR applications. J Math Chem 53(9):2028–2064
MarreroPonce Y et al (2015) Optimum search strategies or novel 3D molecular descriptors: is there a stalemate? Curr Bioinf 10(5):533–564
GarcíaJacas CR et al (2014) Nlinear algebraic maps to codify chemical structures: is a suitable generalization to the atompairs approaches? Curr Drug Metab 15(4):441–469
Edwards CH, Penney DE (1988) Elementary linear algebra. Prentice Hall, Englewoods Cliffs
Johnson RW, Huang CH, Johnson JR (1991) Multilinear algebra and parallel programming. J Supercomput 5(2–3):189–217
Godden JW, Stahura FL, Bajorath J (2000) Variability of molecular descriptors in compound databases revealed by Shannon entropy calculations. J Chem Inf Comput Sci 40(3):796–800
Urias RWP et al (2015) IMMAN: free software for information theorybased chemometric analysis. Mol Divers 19(2):305–319
Somorjai RL (2010) Multivariate statistical methods. In: John L (ed) Encyclopedia of spectroscopy and spectrometry. Academic Press, Oxford, pp 1704–1709
Yap CW (2011) PaDELdescriptor: an open source software to calculate molecular descriptors and fingerprints. J Comput Chem 32(7):1466–1474
Georg H (2008) BlueDescmolecular descriptor calculator. University of Tübingen, Tübingen
Hong H et al (2008) Mold2, molecular descriptors from 2D structures for chemoinformatics and toxicoinformatics. J Chem Inf Comput Sci 48(7):1337–1344
Mauri A et al (2006) DRAGON software: an easy approach to molecular descriptor calculations. Match 56(2):237–248
Steinbeck C et al (2003) The Chemistry Development Kit (CDK): an opensource Java library for chemo and bioinformatics. J Chem Inf Comput Sci 43(2):493–500
Sinkhorn R, Knopp P (1967) Concerning nonnegative matrices and doubly stochastic matrices. Pacific J Math 21(2):343–348
Barigye SJ et al (2013) Shannon’s, mutual, conditional and joint entropybased information indices. Generalization of global indices defined from local vertex invariants. Curr Comput Aided Drug Des 9(2):164–183
Barigye SJ et al (2013) Relations frequency hypermatrices in mutual, conditional and joint entropybased information indices. J Comput Chem 34(4):259–274
MarreroPonce Y et al (2012) Derivatives in discrete mathematics: a novel graphtheoretical invariant for generating new 2/3D molecular descriptors. I. Theory and QSPR application. J Comput Aided Mol Des 26(11):1229–1246
GarcíaJacas CR et al (2014) QuBiLSMIDAS: a parallel freesoftware for molecular descriptors computation based on multilinear algebraic maps. J Comput Chem 35(18):1395–1409
GarcíaJacas CR et al (2015) Multiserver approach for highthroughput molecular descriptors calculation based on multilinear algebraic maps. Mol Inform 34(1):60–69
Manchester J, Czerminski R (2008) SAMFA: simplifying molecular description for 3DQSAR. J Chem Inf Model 48(6):1167–1173
Hinselmann G et al (2011) jCompoundMapper: An open source Java library and commandline tool for chemical fingerprints. J Cheminform 3(1):3
Tosco P, Balle T (2011) A 3DQSARdriven approach to binding mode and affinity prediction. J Chem Inf Model 52(2):302–307
Klamt A et al (2012) COSMOsar3D: molecular field analysis based on local COSMO σprofiles. J Chem Inf Model 52(8):2157–2164
Bonachéra F, Horvath D (2008) Fuzzy tricentric pharmacophore fingerprints. 2. Application of topological fuzzy pharmacophore triplets in quantitative structure–activity relationships. J Chem Inf Model 48(2):409–425
Leardi R, Boggia R, Terrile M (1992) Genetic algorithms as a strategy for feature selection. J Chemom 6(5):267–281
Todeschini R et al (2003) MobyDigs: software for regression and classification models by genetic algorithms. In: Leardi R (ed) Natureinspired methods in chemometrics: genetic algorithms and artificial neural networks. Elsevier, Amsterdam, pp 141–167
Wu CFJ (1986) Jackknife, bootstrap and other resampling methods in regression analysis. Ann Stat 14(4):1261–1295
Lindgren F et al (1996) Model validation by permutation tests: applications to variable selection. J Chemom 10(5–6):521–532
Elisseeff A, Pontil M (2003) Leaveoneout error and stability of learning algorithms with applications. NATO science series sub series III computer and systems sciences, vol 190, pp 111–130
Lilliefors HW (1967) On the Kolmogorov–Smirnov test for normality with mean and variance unknown. J Am Stat Assoc 62(318):399–402
Shapiro SS, Wilk MB (1965) An analysis of variance test for normality (complete samples). Biometrika 52(3/4):591–611
Friedman M (1940) A comparison of alternative tests of significance for the problem of m rankings. Ann Math Stat 11(1):86–92
Siegel S (1957) Nonparametric statistics. Am Stat 11(3):13–19
Wilcoxon F (1945) Individual comparisons by ranking methods. Biometrics 1(6):80–83
Benjamini Y, Hochberg Y (1995) Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc A 57(1):289–300
Hechinger M, Leonhard K, Marquardt W (2012) What is wrong with quantitative structureproperty relations models based on threedimensional descriptors? J Chem Inf Model 52(8):1984–1993
Miteva MA, Guyon F, Tufféry P (2010) Frog2: efficient 3D conformation ensemble generator for small compounds. Nucleic Acids Res 38(suppl 2):W622–W627
RDKit: cheminformatics and machine learning software. February 2, 2016; http://www.rdkit.org/
Vainio MJ, Johnson MS (2007) Generating conformer ensembles using a multiobjective genetic algorithm. J Chem Inf Model 47(6):2462–2474
O’Boyle N et al (2011) Open Babel: an open chemical toolbox. J Cheminform 3(1):33
Standardizer ChemAxon 5.9.0. February 2, 2016. https://www.chemaxon.com/products/standardizer/
Ebejer JP, Morris GM, Deane CM (2012) Freely available conformer generation methods: how good are they? J Chem Inf Model 52(5):1146–1158
Authors’ contributions
CRGJ proposed the theory of the QuBiLSMIDAS indices, supervised the QSAR modeling on the eight chemical datasets, performed the study about the performance of the indices according to several structure generation methods and prepared the manuscript. ECT worked in the QSAR modeling on the eight chemical datasets. YMP leaded the research related with this manuscript. MPM performed the statistical analysis. SJB worked in the definition of the QuBiLSMIDAS indices and prepared the manuscript. LCL worked in the QSAR modeling on the eight datasets. All authors read and approved the final manuscript.
Acknowledgements
GarcíaJacas, CR. thanks the program “International Professor” for a fellowship to work at “Pontificia Universidad Católica del Ecuador Sede Esmeraldas (PUCESE)” in 2015–2016. Barigye, S.J acknowledges support from Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq). MP. Y acknowledge also the partial financial support from Colegio de Medicina, USFQ. Finally, but not least, this work was supported in part by ISCUSFQ.
Competing interests
The authors declare that they have no competing interests.
Author information
Authors and Affiliations
Corresponding author
Additional file
13321_2016_122_MOESM1_ESM.rar
Additional file 1. The statistical parameters and equations of the 3–9 variable QSAR models developed and the corresponding outputs for the statistical analysis performed.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
GarcíaJacas, C.R., ContrerasTorres, E., MarreroPonce, Y. et al. Examining the predictive accuracy of the novel 3D Nlinear algebraic molecular codifications on benchmark datasets. J Cheminform 8, 10 (2016). https://doi.org/10.1186/s133210160122x
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s133210160122x