- Research article
- Open Access
- Published:
Examining the predictive accuracy of the novel 3D N-linear algebraic molecular codifications on benchmark datasets
Journal of Cheminformatics volume 8, Article number: 10 (2016)
Abstract
Background
Recently, novel 3D alignment-free molecular descriptors (also known as QuBiLS-MIDAS) based on two-linear, three-linear and four-linear 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 multi-metrics. 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, cyclooxygenase-2 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 leave-one-out cross-validation correlation coefficients \(\left( {Q_{loo}^{2} } \right)\) reveal that the models based on QuBiLS-MIDAS 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 non-parametric statistic tests were performed, which demonstrated that the models based on QuBiLS-MIDAS 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 QuBiLS-MIDAS 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 QuBiLS-MIDAS 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.

Comparative graphical representation of the performance of the novel QuBiLS-MIDAS 3D-MDs with respect to other methodologies in QSAR modeling of eight chemical datasets
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 2D-QSAR 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 3D-QSAR 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 three-dimensional (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 well-established 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 N-linear 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 atom-pairs by using several (dis-)similarity metrics [25]. Also, the N-linear 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 QuBiLS-MIDAS [acronym of Quadratic, Bilinear and N-Linear Maps based on N-tuple 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 N-linear indices not only do they codify all information contained in the 3D DRAGON MDs, but capture information orthogonal to the latter. Lastly, the QuBiLS-MIDAS MDs were used for modeling the binding affinity to the corticosteroid-binding 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 QuBiLS-MIDAS 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 N-linear algebraic MDs in QSAR modeling with different molecular sets. Therefore, this paper is dedicated to the assessment of the utility of the QuBiLS-MIDAS 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 N-linear algebraic molecular descriptors
In this report, the total and local-fragment 3D N-linear Algebraic indices [25, 26] (also known as QuBiLS-MIDAS) 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 two-linear, three-linear and four-linear algebraic indices for each atom “a” are computed as N-linear (Multi-linear) 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 atom-level indices (F) L, F is a local-fragment (group or atom-type) 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 two-tuple, three-tuple and four-tuple 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 two-tuple, three-tuple and four-tuple 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 atom-based properties: atomic mass (m), the van der Waals volume (v), the atomic polarizability (p), atomic electronegativity in Pauling scale (e), atomic Ghose-Crippen LogP (a), Gasteiger-Marsili 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 two-linear, three-linear and four-linear 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 multi-metrics to capture the information on the relations between two, three and four atoms, respectively [25, 26]. To compute the atom-pair 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 multi-metrics (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 two-tuple and three-tuple total spatial-(dis)similarity matrices calculated with some previously mentioned metrics and multi-metrics.
From these total matrix approaches (\({\mathbb{G}}^{k}\), \({{\mathbb{G}}{\mathbb{T}}}^{k}\) and \({{\mathbb{G}}{\mathbb{Q}}}^{k}\)), local-fragments matrices may be computed in order to consider atom-types or chemical regions of interest and thus yielding the kth two-tuple, three-tuple and four-tuple 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 local-fragments (or atom-types), 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 local-fragment) 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 two-tuple, three-tuple and four-tuple total (or local-fragment) 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 QuBiLS-MIDAS MDs. In this way, the following normalized matrix representations are obtained from the corresponding non-stochastic matrices: the kth simple - stochastic two-tuple, three-tuple and four-tuple total (or local-fragment) 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 two-tuple total (or local-fragment) spatial-(dis)similarity matrix [\(_{ds} {\mathbb{G}}_{(F)}^{k}\)] (see Sinkhorn–Knopp algorithm in Ref. [37]) and the kth mutual probability two-tuple, three-tuple and four-tuple total (or local-fragment) 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 two-tuple total spatial-(dis)similarity matrix.
Finally, from the non-stochastic (simple-stochastic, double-stochastic or mutual-probability) total (or local-fragment) matrices [i.e. \({\mathbb{G}}_{(F)}^{k}\), \({{\mathbb{G}}{\mathbb{T}}}_{(F)}^{k}\) and \({{\mathbb{G}}{\mathbb{Q}}}_{(F)}^{k}\)], the corresponding atom-level 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 atom-level matrix determines an atom-level 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 two-linear, three-linear and four-linear 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 two-linear algebraic index.
General workflow for calculating the QuBiLS MIDAS molecular descriptors. (1) Computation of the molecular vectors according to selected atomic properties; (2) Computation from 3D Cartesian coordinates of each atom of a molecule the non-stochastic two-tuple, three-tuple or four-tuple total spatial-(dis)similarity matrices for k = 1; (3) Consideration of atom-types or local-fragments (optional); (4) Computation of the simple-stochastic, double-stochastic and mutual probability matrices, as well as to determine the kth matrices through Hadamard product until the k value selected; (5) Splitting the calculated matrices into atom-level matrices; (6) Computation of the atom-level indices (descriptors) using the molecular vectors calculated in the step (1); and (7) Application of the selected aggregation operators over vector of atom-level descriptors
General workflow for the calculation of a two-linear descriptor based on the linear algebraic form, Euclidean metric, non-stochastic matrix approach, atomic mass as property and Manhattan aggregation operator. (1) Computation of the non-stochastic matrix for k = 1 \(\left( {{\mathbb{G}}^{1} } \right)\) from the 3D coordinates matrix and using the Euclidean metric; (2) Computation of the molecular vector based on the atomic mass property, \(\bar{\varvec{X}}_{\varvec{m}}\); (3) Splitting of the \({\mathbb{G}}^{1}\) matrix into “n” (number of atoms) atom-level matrices, \({\mathbb{G}}^{{\varvec{a},1}}\), where “a” is an atom of the molecule; (4) Computation of the atom-level descriptors and saving them into vector L; and (5) Application of the Manhattan aggregation operator over the entries of the vector L, being this value the molecular descriptor
In order to automatize the calculation of the 3D N-linear algebraic indices used in the present manuscript the QuBiLS-MIDAS software has been developed [41]. This software has as one of its main features the multi-core processing of the MDs, as well as the option to carry out the distributed calculation of the indices by using the Multi-Server Distributed Computing Platform known as T-arenal [42]. The latter is particularly useful for high-throughput calculation tasks. Both software are freely available via internet at: http://tomocomd.com/.
Methods
In order to assess the correlation ability of the QuBiLS-MIDAS MDs for different biological activities eight well-known chemical datasets were used. These were previously employed by Sutherland et al. in a comparative study of QSAR methods commonly used in chemo-informatics 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), cyclooxygenase-2 (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 two-linear, three-linear and four-linear 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 high-dimensional 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) meta-heuristic 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 N-linear 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 trade-off 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}\) (“leave-one-out” cross validation) statistical parameter as the fitness function. Once the exploration in each population was completed, then the MDs included in the built 9-variable 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 “Y-scrambling” [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 QuBiLS-MIDAS models versus other approaches
In this section the performance of the QuBiLS-MIDAS 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 local-fragment QuBiLS-MIDAS 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 QuBiLS-MIDAS 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 QuBiLS-MIDAS 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 local-fragments (atom-types or group) in the calculation of the QuBiLS-MIDAS MDs and subsequently in the building of QSAR models can be analyzed. As can be observed in both tables, the performance of the QuBiLS-MIDAS models is superior when local-fragments are considered with respect to those QuBiLS-MIDAS 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 local-fragment QuBiLS-MIDAS MDs not influence the performance of the developed QSAR models. It can thus be stated that considering a mixture of total and local-fragment QuBiLS-MIDAS MDs in building of QSAR models contributes to the improvement of the predictive ability.
Also, it can be observed from Table 7 that the cross-validation performances achieved by the QuBiLS-MIDAS models have comparable-to-superior 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], SAMFA-RF (\({\text{Q}}_{\text{loo}}^{2}\) = 0.58), All-Shortest Path [ASP] Fingerprint (\({\text{Q}}_{\text{loo}}^{2}\) = 0.56), COMSIA extra (\({\text{Q}}_{\text{loo}}^{2}\) = 0.57), HQSAR (and SAMFA-RF) [\({\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 QuBiLS-MIDAS 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 QuBiLS-MIDAS 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 QuBiLS-MIDAS indices yield comparable-to-superior 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 2D-FPT 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 2D-FPT models were developed by using SQS framework that determines linear and non-linear 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 QuBiLS-MIDAS approach is inferior to the best results reported so far in the literature. In the former dataset, the methods 2D-FPT (\({\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 2D-FPT approach (\({\text{Q}}_{\text{ext}}^{2}\) = 0.683) has more predictive ability than the corresponding QuBiLS-MIDAS 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 2D-FPT methodology for ACHE and DHFR datasets, the achieved results are based on non-linear models while the proposed outcomes are determined with linear models.
The obtained results evidence that the QuBiLS-MIDAS 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 comparable-to-superior 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 QuBiLS-MIDAS 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 QuBiLS-MIDAS and 2D-FTP 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 QuBiLS-MIDAS models. In addition, taking into account the graphics corresponding to the QuBiLS-MIDAS and 2D-FPT approaches, it can be concluded that the predictions obtained by the former are less scattered than those attained by the latter and thus, the QuBiLS-MIDAS 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 QuBiLS-MIDAS 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 non-parametric 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 2D-FTP and COSMOsar3D models, although with Shapiro–Wilk test the rejection of the null hypothesis is achieved for COMSIA basic as well. Therefore the non-parametric 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 QuBiLS-MIDAS models being those with the best performance followed by the 2D-FPT, 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 signed-ranks test [57] was carried out (with R software) by using Benjamini and Hochberg [58] (BH) as the adjustment method (one-tailed 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 QuBiLS-MIDAS models yield statistically better predictions than the other methodologies considered, with the exception of the 2D-FPT 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 computer-aided 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 QuBiLS-MIDAS 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 QuBiLS-MIDAS 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 QuBiLS-MIDAS MDs depend on 3D conformations from which are computed.
Boxplot graphic for the external predictive accuracy achieved by the QSAR models reported in this manuscript (see Table 6) and fitted using structures generated by CORINA software, over the corresponding test sets optimized by five different toolkits
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 QuBiLS-MIDAS 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 two-linear 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 dependent-variable and with low intercorrelation among them.
-
The MLR-GA procedure implemented in MobyDigs software was employed to build 9-variable 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 QuBiLS-MIDAS 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 QuBiLS-MIDAS software, RDKIT program will be incorporated in the QuBiLS-MIDAS software as a built-in 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 alignment-free geometric molecular descriptors based on N-linear algebraic maps (so called QuBiLS-MIDAS) 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 QuBiLS-MIDAS models were compared with respect to several QSAR procedures reported in the literature according to the correlation coefficients achieved with the leave-one-out cross-validation \(\left( {Q_{loo}^{2} } \right)\) and external prediction \(\left( {Q_{ext}^{2} } \right)\) methods, and generally superior performances were observed with this QuBiLS-MIDAS framework.
A few exceptions were observed: for the \(Q_{loo}^{2}\) parameter, the QuBiLS-MIDAS approach is exclusively outperformed by the ASP-based (fingerprint) method in the DHFR dataset, while for the \(Q_{ext}^{2}\) parameter, the QuBiLS-MIDAS method yields inferior results with respect to the 2D-FPT methodology in the DHFR and ACHE test set, respectively. Also, inferior \(Q_{ext}^{2}\) values are yielded by the QuBiLS-MIDAS 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 non-congeneric datasets because are alignment-depend. Thus, considering the maximum parsimony principle (Ockham’s razor), the QuBiLS-MIDAS 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 QuBiLS-MIDAS 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 N-linear molecular descriptors (also known as QuBiLS-MIDAS) 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
Sung-Sau S, Karplus M (1997) Three-dimensional quantitative structure–activity relationships from molecular similarity matrices and genetic neural networks. 2. Applications. J Med Chem 40(26):4360–4371
Aires-de-Sousa J, Gasteiger J (2002) Prediction of enantiomeric selectivity in chromatography: application of conformation-dependent and conformation-independent 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. Wiley-VCH, 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 three-dimensional 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 vibration-based 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 three-dimensional geometry. Springer, New York
Bogdanov B, Nikolic S, Trinajstic N (1990) On the three-dimensional 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 anti-3-azopyridine antibody. Arzneim Forsch 36(2):176–183
Randić M (1995) Molecular profiles novel geometry-dependent 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 atom-based linear descriptors of the Minkowski distance matrix: theory, diversity–variability analysis and QSPR applications. J Math Chem 53(9):2028–2064
Marrero-Ponce Y et al (2015) Optimum search strategies or novel 3D molecular descriptors: is there a stalemate? Curr Bioinf 10(5):533–564
García-Jacas CR et al (2014) N-linear algebraic maps to codify chemical structures: is a suitable generalization to the atom-pairs 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 theory-based 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) PaDEL-descriptor: an open source software to calculate molecular descriptors and fingerprints. J Comput Chem 32(7):1466–1474
Georg H (2008) BlueDesc-molecular 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 open-source 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 entropy-based 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 entropy-based information indices. J Comput Chem 34(4):259–274
Marrero-Ponce Y et al (2012) Derivatives in discrete mathematics: a novel graph-theoretical invariant for generating new 2/3D molecular descriptors. I. Theory and QSPR application. J Comput Aided Mol Des 26(11):1229–1246
García-Jacas CR et al (2014) QuBiLS-MIDAS: a parallel free-software for molecular descriptors computation based on multi-linear algebraic maps. J Comput Chem 35(18):1395–1409
García-Jacas CR et al (2015) Multi-server approach for high-throughput molecular descriptors calculation based on multi-linear algebraic maps. Mol Inform 34(1):60–69
Manchester J, Czerminski R (2008) SAMFA: simplifying molecular description for 3D-QSAR. J Chem Inf Model 48(6):1167–1173
Hinselmann G et al (2011) jCompoundMapper: An open source Java library and command-line tool for chemical fingerprints. J Cheminform 3(1):3
Tosco P, Balle T (2011) A 3D-QSAR-driven 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) Nature-inspired 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) Leave-one-out 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 structure-property relations models based on three-dimensional 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 J-P, 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 QuBiLS-MIDAS 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 QuBiLS-MIDAS 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ía-Jacas, 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). M-P. Y acknowledge also the partial financial support from Colegio de Medicina, USFQ. Finally, but not least, this work was supported in part by ISC-USFQ.
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ía-Jacas, C.R., Contreras-Torres, E., Marrero-Ponce, Y. et al. Examining the predictive accuracy of the novel 3D N-linear algebraic molecular codifications on benchmark datasets. J Cheminform 8, 10 (2016). https://doi.org/10.1186/s13321-016-0122-x
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s13321-016-0122-x
Keywords
- Multiple Linear Regression
- QuBiLS-MIDAS
- 3D-QSAR
- TOMOCOMD-CARDD