Next Article in Journal
The Yeast Saccharomyces cerevisiae as a Model for Understanding RAS Proteins and their Role in Human Tumorigenesis
Previous Article in Journal
MiRNAs at the Crossroads between Innate Immunity and Cancer: Focus on Macrophages
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Linear Regression QSAR Models for Polo-Like Kinase-1 Inhibitors

by
Pablo R. Duchowicz
Instituto de Investigaciones Fisicoquímicas Teóricas y Aplicadas (INIFTA), CONICET, UNLP, Diag. 113 y 64, C.C. 16, Sucursal 4, La Plata 1900, Argentina
Cells 2018, 7(2), 13; https://doi.org/10.3390/cells7020013
Submission received: 27 November 2017 / Revised: 23 January 2018 / Accepted: 2 February 2018 / Published: 14 February 2018

Abstract

:
A structurally diverse dataset of 530 polo-like kinase-1 (PLK1) inhibitors is compiled from the ChEMBL database and studied by means of a conformation-independent quantitative structure-activity relationship (QSAR) approach. A large number (26,761) of molecular descriptors are explored with the main intention of capturing the most relevant structural characteristics affecting the bioactivity. The structural descriptors are derived with different freeware, such as PaDEL, Mold2, and QuBiLs-MAS; such descriptor software complements each other and improves the QSAR results. The best multivariable linear regression models are found with the replacement method variable subset selection technique. The balanced subsets method partitions the dataset into training, validation, and test sets. It is found that the proposed linear QSAR model improves previously reported models by leading to a simpler alternative structure-activity relationship.

Graphical Abstract

1. Introduction

Polo-like kinases (PLKs) are characterized by a multidomain structure consisting of a highly conserved N-terminal catalytic domain (KD) and a relatively divergent C-terminal polo-box domain (PBD), composed of either one or two polo boxes [1,2]. This serine/threonine kinase family is an important regulator of mitotic progression [3].
Among the different identified PLKs, PLK1 is the most investigated member of the family because it is highly expressed in proliferating cells and overexpressed in many cancers, thus resulting in an attractive target for anticancer therapeutic development [4].
Polo-like kinase 1 (PLK1) is involved in centrosome maturation, kinetochore function, spindle formation, chromosome segregation, and cytokinesis [5]. The PBD is critical to PLK1 localization and function and negatively regulates the kinase activity of the catalytic domain [1]. Inhibitors targeting KD, the so-called ATP-competitive PLK1 inhibitors, have attracted much attention over the last years [6].
Inhibiting PLK1 activity results in a potent antitumor effect both in vitro and in vivo [7]; therefore, much attention has been focused on characterizing PLK1 and synthesizing its inhibitors [5]. Figure 1 shows the molecular structures of some encouraging agents in current clinical trials.
Among the various methodologies available in the literature for predicting the biological activities of compounds, the quantitative structure-activity relationship (QSAR) theory [8,9,10,11,12] is considered as a useful and well-known strategy. The main hypothesis behind every QSAR study is that the chemical structure is responsible for the bioactivity of the compound.
Therefore, the structure is quantified through molecular descriptors; in other words, numerical quantities carrying specific and relevant information about the constitutional, topological, geometrical, hydrophobic, and/or electronic characteristics of the compounds are investigated [13,14,15,16]. Thousands of molecular descriptors are now available in the literature, and it has to be decided how to select those numerical variables that best characterize the experimental activity under consideration.
The most relevant descriptors are selected with an appropriate mathematical technique and statistically correlated to the experimental activity, resulting in a mathematical model that is used to find out useful structure-activity parallelisms. In this way, QSAR models constitute a fast and cost-effective alternative to experimental measurements.
The availability of newer and higher quality experimental measurements has encouraged and justified the development of newer and alternative QSAR models with improved statistical quality; therefore, this research field has continued to evolve over the years.
A large number of structure-activity relationship (SAR) studies has been reported in the past for analyzing the PLK1 inhibition although, however, very few account for quantitative approaches [17,18,19,20]. Therefore, in the present work, we develop a QSAR analysis for searching predictive models on a large and structurally diverse dataset of 530 PLK1 inhibitors. For this purpose, we resort to the conformation-independent QSAR approach and consider constitutional and topological representations of the inhibitors’ chemical structures for deriving the molecular descriptors.
The main advantage of not considering molecular conformations is that the only experimental data needed for establishing the QSAR models is the experimental inhibitory activity being analyzed [21,22,23,24]. No further experimental information is required, such as data on the experimental X-ray crystal structure of the PLK1 kinase domain with a given inhibitor in certain conformation [25]. In addition, it is known that the ligand–receptor complex has not been solved for all inhibitor types having different interaction modes—that is, in a heterogeneous dataset like the one considered in this work.
The inclusion of more specific experimental information is appropriate whenever a microscopic and more sophisticated modeling methodology is involved, such as CoMFA (comparative molecular field analysis) and CoMSIA (comparative molecular similarity indices analysis) [17,26]. However, as commented previously, such specific experimental details are usually unavailable for any chemical system of interest; therefore, the application of the conformation-independent QSAR approach can be considered as a useful and valid alternative.

2. Materials and Methods

2.1. Experimental Dataset

The structurally diverse PLK1 inhibitors were compiled from ChEMBL [27,28], an open data resource of binding, functional, and ADMET bioactivity data. The experimental inhibitory effectiveness is expressed as the half-maximal inhibitory concentration I C 50 (nM).
After removing duplicates, compounds with ambiguous data, compounds having molecular weights higher than 1000 g mol−1, and compounds without reported bioactivities, the dataset consisted of 530 compounds with I C 50 values ranging from 0.8 to 145,000 nM and molecular weights ranging from 164.2 to 949.97 g mol−1. The complete list of compounds studied here is provided in Table S1 as Supplementary Material.

2.2. Structural Representation and Molecular Descriptors Calculation

The 530 chemical structures studied here are provided as canonical SMILES notation in Table S1 of the Supplementary Material section. All file format conversions were performed with Open Babel for Windows [29]. The molecular structures were visualized with ACDLabs ChemSketch freeware [30].
The conformation-independent molecular descriptors were computed as follows. We used the Pharmaceutical Data Exploration Laboratory (PaDEL) freeware version 2.20 [31] because it has the advantage that it is a freely available and open-source program. PaDEL allowed us to calculate 1444 0D-2D descriptors and 12 fingerprint types (16,092 bits) [32]. The categorical (indicator) fingerprint descriptors involve the presence or count of specific chemical substructures: we treated the fingerprints like they were “constitutional descriptors” describing the molecular composition and, as such, they could be used for modeling any property of interest.
More molecular descriptors were calculated with the Molecular Descriptors from 2D structures (Mold2) freeware [33], which generated 777 1D-2D structural variables with molecules in MDL sdf format.
Finally, 2D molecular descriptors were calculated with the Quadratic, Bilinear and N-Linear MapS (QuBiLs) [34] suite by using the graph-theoretic electronic-density matrices and atomic weightings (MAS) module from the ToMoCoMD-CARDD free multi-platform freeware. The QuBiLs-MAS algebraic module calculated 8448 quadratic, bilinear, and linear maps based on pseudograph-theoretic electronic-density matrices and atomic weightings, when the program was used with the following options selected: ‘bilinear’, ‘linear’, and ‘quadratic’ algebraic forms; ‘atom-based’, ‘non-chiral’ and ‘duplex’ constraints; ‘non-stochastic’, ‘simple stochastic’, ‘double stochastic’ and ‘mutual probability’ matrix forms (maximum order 15); ‘keep all’ cut-off; ‘total’ groups; ‘Ghose-Crippen LogP’, ‘polarizability’, ‘charge’, ‘polar surface area’, ‘electronegativity’, ‘refractivity’, ‘mass’ and ‘van der Waals volume’ properties; ‘Euclidean distance’, ‘arithmetic mean’, and ‘standard deviation’ invariants (non-standardized option).
Through PaDEL, Mold2, and QuBiLs-MAS we derived 26,761 non-conformational molecular descriptors with the intention of exploring the most relevant structural characteristics affecting the studied PLK1 bioactivity.

2.3. Model Development

2.3.1. Molecular Descriptors Selection

First, the ‘collinear’ or linearly dependent descriptor pairs were identified, and only one variable from each pair was kept for further analysis. Non-informative descriptors not relevant to the QSAR analysis were excluded (i.e., descriptors with constant and near-constant values and descriptors with at least one missing value), leading to a pool of 11,565 linearly independent non-conformational descriptors.
We employed the replacement method (RM) technique [35] in order to generate multivariable linear regression (MLR) models on the training set (train) by searching in a pool having D = 11,565 descriptors for optimal subsets containing d descriptors (d is much lower than D), with the smallest values for the standard deviation ( S t r a i n ).
The main idea behind the RM is that one can approach the minimum of S t r a i n by judiciously taking into account the relative errors of the coefficients of the least-squares model given by a set of d descriptors. In other words, we should find the global minimum of S t r a i n ( d ) in a subspace of D ! / [ d ! ( D d ) ! ] points d, where D represents the total number of available descriptors.
The quality of the results achieved with this technique approaches that obtained by performing an exact (combinatorial) full search of molecular descriptors although, of course, requires much less computational work. The RM is computationally more expensive than the stepwise regression (SR) and genetics algorithm (GA) approaches, although produces similar or better results than GA and better results than SR [36].
Table S2 includes a list of mathematical equations involved in the present study. All the MatLab programmed algorithms used in our calculations are available upon request.

2.3.2. Model Validation

The complete molecular set of 530 inhibitors was split into three subsets: training (train), validation (val), and test sets. The training set was used to calibrate the model and to obtain its parameters through the RM technique, while the validation set helped to calibrate and partially validate the model by predicting the bioactivity of compounds not included in train. Finally, the test set contained compounds “never seen” during the calibration step with train and val, and demonstrates the real predictive capability of the QSAR.
The dataset partitioning has to achieve similar structure-activity relationships in the three subsets; in other words, the training set molecules should be representative of the validation and test set compounds. For this purpose, the split of the dataset was carried out by means of the balanced subsets method (BSM) [37,38], a procedure proposed by our group that ensures that balanced subsets are generated. The BSM is based on the k-means cluster analysis (k-MCA) method [39]: the essence of k-MCA is to create k-clusters or groups of compounds in such a way that compounds in the same cluster are very similar in terms of distance metrics (i.e., Euclidean distance), and compounds in different clusters are very distinct.
The linear regression models are also theoretically validated through the leave-one-out cross-validation (loo) procedure [40], and also through the more rigorous leave-30%-out cross-validation (l30‰) method.

2.3.3. Applicability Domain

A predictive QSAR model is only able to predict molecules falling within its applicability domain (AD), so that the predicted activity is not a result of substantial extrapolation (unreliable prediction) [41,42]. The AD definition is dependent on the model’s descriptors and the experimental activity.
In this work, we determined the AD through two alternative methodologies. The first one is based on the well-known leverage approach [43] where a test set compound i must have a calculated leverage h i smaller than the warning leverage h * . The second one is based on a simple standardization approach [42]: a given test set compound i having d standardized descriptor values s i k , k = 1 , , d must have a maximum value s i k max 3 . In the case that s i k max > 3 and its minimum value s i k min < 3 , then the s i n e w parameter has to be calculated and must fulfill the condition: s i n e w = s i + 1.28 . σ s i 3 , where s i is the mean of s i k values for i and σ s i is the standard deviation for such values.

3. Results and Discussion

After partitioning the dataset of 530 PLK1 inhibitors into train, val, and test sets using the BSM technique, we obtain balanced subsets with N t r a i n = 265 , N v a l = 133 and N t e s t = 132 compounds; in addition, Table S1 denotes the members of val (^) and test (*) sets. Therefore, the calibration compounds in train and val constitute 75% of the whole dataset.
The best MLR models, including the most representative 1–9 molecular descriptors, are presented in Table 1. A brief description of such descriptors is also supplied in Table S3. From the results of Table 1, it is clearly appreciated that the S t r a i n parameter continuously improves with the addition of molecular descriptors into the linear equation. However, according to the validation set results, the most predictive models (lowest S v a l ) have 8 and 9 descriptors. We kept the model’s dimension as small as possible and selected the following 8-descriptor model and associated statistical quality:
log 10 I C 50 = 0.46 m i n d s s C 0.85 m a x H C s a t s + 0.88 M 66 0.54 P C 494   2.76 P C 534 1.12 P C 686 + 2.68 K R 3577 1.44 K R 4268 + 4.37
N t r a i n = 265 ,   R t r a i n 2 = 0.69 ,   S t r a i n =   0.80
R i j 2 max = 0.14 ,   V I F max = 1.11 ,   o 3 = 1 ,   R r a n d 2 = 0.13 ,   S r a n d = 1.35
R l o o 2 = 0.67 ,   S l o o = 0.83 ,   R l 30 % o 2 = 0.58 ,   S l 30 % o = 0.95
N v a l = 133 ,   R v a l 2 = 0.75 ,   S v a l = 0.82
N t e s t = 132 ,   R t e s t 2 = 0.69 ,   S t e s t = 0.85
A plot for the log 10 I C 50 predictions given by Equation (1) as a function of the experimental values is provided in Figure 2. The dispersion plot of residuals in Figure 3 tends to obey a random pattern around the zero line, suggesting that Equation (1) predicts the whole dataset without systematic errors or residual bias.
The o 3 parameter indicates the number of outlier compounds in the training set having a residual (difference between experimental and predicted activity) greater than 3 times S t r a i n . The only outlier in the training set is 171, 1-{4-[(4-chlorophenyl)methoxy]-3-methoxyphenyl}-N-[(pyridin-4-yl)methyl]methanamine. After close inspection of this specific compound, it is easily concluded that the abnormal behavior can be completely attributed to the highly heterogeneous dataset being analyzed, involving molecular weights from 164.2 to 949.97 g mol−1 and bioactivities from 0.8 to 145,000 nM.
Our proposed 8-descriptor model approves the internal validation process of loo and l30‰ (500,000 cases) cross-validation procedures through the prediction of 1 or 80 molecules excluded at a time from the training set. According to the specialized literature [40], the cross-validation R l o o 2 and R l 30 2 explained variances should be greater than 0.5, although this is a necessary but not sufficient condition for the real predictive power.
As a way of demonstrating that the QSAR model is not a result of chance correlation, the experimental log 10 I C 50 activity values were scrambled with Y-randomization [44] (100,000 cases). When S r a n d (S for Y-randomization) is greater than S t r a i n , a valid and useful structure-activity relationship is found, as is the case for Equation (1).
The recommended external validation criteria to assure predictive capability [40] are also achieved: 1 R 0 2 / R t e s t 2 ( 1.18 × 10 3 ) < 0.1 or 1 R 0 2 / R t e s t 2 ( 0.12 ) < 0.1 ; 0.85 k ( 0.93 ) 1.15 or 0.85 k ( 1.00 ) 1.15 ; R m 2 ( 0.67 ) > 0.5 .
The R i j max parameter from Equation (1) is the maximum correlation coefficient between descriptor pairs: R i j 2 max = 0.14 indicates that there is no serious overlapping structural information. V I F max is the maximum variance inflation factor, a parameter that measures the maximum multicollinearity among descriptors. A VIF of 1 for a specific descriptor means that there is no correlation between this descriptor and all the remaining descriptors of the model, and a VIF exceeding 10 indicates that multicollinearity is a problem in the dataset [45]. For Equation (1), V I F max = 1.11 . The complete squared correlation matrix and VIF values are provided in Table S4.
It is known that a successful QSAR model is established only when it surpasses the validation process, in other words, by testing its ability to predict the experimental bioactivity of compounds that are not considered during the model calibration [46,47]. The QSAR of Equation (1) has an acceptable predictive capability for the external test set of 132 “never seen” experimental log 10 I C 50 values according to R t e s t 2 and S t e s t parameters and Figure 2 and Figure 3. This QSAR can thus be applied to predict new inhibitors with unknown experimental I C 50 .
The eight conformation-independent structural indices of Equation (1) have quite a straightforward structural interpretation:
  • Two electrotopological state atom-type descriptors: mindssC, the minimum atom-type E-state: =C<; and maxHCsats, the maximum atom-type H E-state: H bonded to B, Si, P, Ge, As, Se, Sn, or Pb.
  • A MACCS fingerprint descriptor: M66, the number of CC(C)(C)A fragments, where A is any valid periodic table element symbol.
  • Three PubChem fingerprint descriptors: PC494, the presence of O=C-C:N fragment, where ‘:’ denotes bond aromaticity; PC534, the presence of S-C:C-O fragment; and PC686, the presence of O=C-C-C-C-O fragment.
  • Two Klekota–Roth fingerprint descriptors: KR3577, the presence of SMARTS substructure Cc1cccc(C)c1NC=O; and KR4268, the presence of SMARTS substructure Nc1ccccc1O.
The numerical values for these descriptors are provided in Table S5: all of them have positive numerical values with the exception of mindssC, which has either positive or negative values. The sign of the regression coefficient in the linear model indicates when the descriptor contribution increases or decreases the predicted log 10 I C 50 values. Therefore, it is possible to propose the following useful QSAR guide for the chemical synthesis of new PLK1 inhibitors. Molecular structures of inhibitors simultaneously having higher positive values of m a x H C s a t s , P C 494 , P C 534 , P C 686 , and K R 4268 and lower values for m i n d s s C , M 66 , and K R 3577 would exhibit lower predicted log 10 I C 50 values, being predicted as more active PLK1 inhibitors.
In order to apply the proposed QSAR guide, the molecular structures to be predicted have to fall within the model’s applicability domain (AD). Within the leverage approach [43], a compound with high leverage ( h i ) would reinforce the model if the compound is in the train or val (good leverage) calibration sets; but such a compound in the test set could have unreliable predicted data, the result of substantial extrapolation of the model (bad leverage) [41]. Equation (1) reveals that most of the test set compounds have h i values falling under the h * limit (0.1019) with the exception of five test set compounds: 495, 508, 509, 511, and 516. The Williams plot (standardized residuals as a function of the h i values) is provided in Figure 4.
This result obtained with the leverage approach for the test set approximately coincides with the one obtained by using the standardization approach, as the two conditions s i k max 3 or s i n e w 3 are followed by all the test set compounds with the exception of seven compounds: the five previous test compounds and two more compounds lying near the h * limit: 453 and 457. Thus, the predicted log 10 I C 50 values for most of the test set compounds can be considered as reliable.
Finally, the obtained regression model in Equation (1) can be converted into a classification model by classifying compounds with experimental IC50 ≤ 1000 nM as highly active inhibitors and experimental IC50 > 1000 nM as poorly active inhibitors. Then, the Cooper statistics [48] related to accuracy (A%), sensitivity (SE), and specificity (SP) and the Matthews correlation coefficient (MCC) can be calculated. The classification results for Equation (1) in the test set are acceptable as A% = 83%, SE = 0.73, SP = 0.95, MCC = 0.69.
A previous study developed by Kong and Yan [20] employed the current ChEMBL database of PLK1 inhibitors for establishing various in silico classification models. The 16 single classifier models and one consensus Kohonen’s self-organizing map (SOM) model were applied to a dataset of 601 noncongeneric PLK1 inhibitors. For these 16 single classifier models, four machine learning methods were used: support vector machine (SVM), naive Bayes (NB), C4.5 decision tree (C4.5 DT), and random forest (RF), with MCC ranging from 0.609 to 0.864 and A% ranging from 78.7% to 93.1% for the test set. Then, a consensus SOM model was built based on four single classifier models to obtain a more reliable and robust model, outperforming all the single classifier models with MCC = 0.87 and A% = 93.6% on the test set.
The models reported in [20] achieved acceptable results. However, the linear QSAR model of Equation (1) represents an improved alternative model having the following characteristics:
i.
Our proposed model performs both regression and classification.
ii.
Dataset partitioning: three subsets are considered, such as train, val, and test instead of only two (train and test) in [20]. In this way, it is more convenient for analyzing the predictive performance of the model.
iii.
Model’s size: a fewer number of molecular descriptors are involved in the final selected model—i.e., 8 instead of 10–15. Therefore, the parsimony´s principle is accomplished (Ockham’s razor) [49] by following the common practice of keeping the model’s dimension as small as possible.
iv.
No energy or geometry optimization is performed on the inhibitor chemical structures. The conformation-independent QSAR approach considers only constitutional and topological representations for deriving the molecular descriptors.
v.
A simpler modeling methodology based on MLR analysis is applied in the present study.

4. Conclusions

Polo-like kinase-1 is an attractive target for anticancer therapeutic development so the prediction of its inhibitors has been of great interest during the last years. The linear regression QSAR model established in this work on a structurally diverse set of 530 PLK1 inhibitors has an acceptable predictive capability in the external test set and is based on eight non-conformational molecular descriptors.
For chemical structures falling within the applicability domain of this model, a QSAR guide for the chemical synthesis of new PLK1 inhibitors is provided as follows: molecular structures of inhibitors simultaneously having higher positive values of maxHCsats, PC494, PC534, PC686, and KR4268 and lower values for mindssC, M66, and KR3577 would exhibit lower predicted log10IC50 values, being predicted as more active PLK1 inhibitors.
The consideration of the constitutional and topological aspects of the molecular structures in the conformation-independent QSAR approach achieves acceptable results. New investigations on other physicochemical and biological properties of interest will be published soon elsewhere.

Supplementary Materials

The following are available online at https://www.mdpi.com/2073-4409/7/2/13/s1.

Acknowledgments

The author is a researcher from the National Scientific and Technical Research Council (CONICET) of Argentina and is grateful for the financial support of PIP11220130100311 project and for the electronic library facilities of the Ministry of Science, Technology and Productive Innovation.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Lowery, D.M.; Lim, D.; Yaffe, M.B. Structure and function of Polo-like kinases. Oncogene 2005, 24, 248–259. [Google Scholar] [CrossRef] [PubMed]
  2. McInnes, C.; Wyatt, M.D. PLK1 as an oncology target: Current status and future potential. Drug Discov. Today 2011, 16, 619–625. [Google Scholar] [CrossRef] [PubMed]
  3. Shakil, S.; Baig, M.H.; Tabrez, S.; Danish Rizvi, S.M.; Zaidi, S.K.; Ashraf, G.M.; Ansari, S.A.; Parwaz Khan, A.A.; Al-Qahtani, M.H.; Abuzenadah, A.M.; et al. Molecular and enzoinformatics perspectives of targeting Polo-like kinase 1 in cancer therapy. Semin. Cancer Biol. 2017. [Google Scholar] [CrossRef] [PubMed]
  4. Lee, K.S.; Burke, T.R., Jr.; Park, J.E.; Bang, J.K.; Lee, E. Recent Advances and New Strategies in Targeting PLK1 for Anticancer Therapy. Trends Pharmacol. Sci. 2015, 36, 858–877. [Google Scholar] [CrossRef] [PubMed]
  5. Strebhardt, K.; Becker, S.; Matthess, Y. Thoughts on the current assessment of Polo-like kinase inhibitor drug discovery. Expert Opin. Drug Discov. 2015, 10, 1–8. [Google Scholar] [CrossRef] [PubMed]
  6. Zhang, J.; Yang, P.Y.; Gray, N.S. Targeting cancer with small molecule kinase inhibitors. Nat. Rev. Cancer 2009, 9, 28–39. [Google Scholar] [CrossRef] [PubMed]
  7. Chopra, P.; Sethi, G.; Dastidar, S.G.; Ray, A. Polo-like kinase inhibitors: An emerging opportunity for cancer therapeutics. Expert Opin. Investig. Drugs 2010, 19, 27–43. [Google Scholar] [CrossRef] [PubMed]
  8. Hansch, C.; Leo, A. Exploring QSAR. Fundamentals and Applications in Chemistry and Biology; American Chemical Society: Washington, DC, USA, 1995. [Google Scholar]
  9. Esposito, E.; Hopfinger, A.J. Multi-Dimensional QSAR: Methods and Applications for Drug Discovery and Polymer Science; CRC Press: New York, NY, USA, 2012. [Google Scholar]
  10. Putz, M.V. QSAR and SPECTRAL-SAR in Computational Ecotoxicology; CRC Press: New York, NY, USA, 2012. [Google Scholar]
  11. Benfenati, E. Theory, Guidance and Applications on QSAR and REACH; Orchestra: Milan, Italy, 2012; Available online: http://ebook.insilico.eu/insilico-ebook-orchestra-benfenati-ed1_rev-June2013.pdf (accessed on 24 November 2017).
  12. Roy, K. Advances in QSAR Modeling. Applications in Pharmaceutical, Chemical, Food, Agricultural and Environmental Sciences; Springer International Publishing: Berlin, Germany, 2017. [Google Scholar]
  13. Katritzky, A.R.; Goordeva, E.V. Traditional Topological Indices vs. Electronic, Geometrical, and Combined Molecular Descriptors in QSAR/QSPR Research. J. Chem. Inf. Comput. Sci. 1993, 33, 835–857. [Google Scholar] [CrossRef] [PubMed]
  14. Diudea, M.V.E. QSPR/QSAR Studies by Molecular Descriptors; Nova Science Publishers: New York, NY, USA, 2001. [Google Scholar]
  15. Todeschini, R.; Consonni, V. Molecular Descriptors for Chemoinformatics (Methods and Principles in Medicinal Chemistry); Wiley-VCH: Weinheim, Germany, 2009. [Google Scholar]
  16. Masand, V.H.; Rastija, V. PyDescriptor: A new PyMOL plugin for calculating thousands of easily understandable molecular descriptors. Chemom. Intell. Lab. Syst. 2017, 169, 12–18. [Google Scholar] [CrossRef]
  17. Cao, S. QSAR, molecular docking studies of thiophene and imidazopyridine derivatives as polo-like kinase 1 inhibitors. J. Mol. Struct. 2012, 1020, 167–176. [Google Scholar] [CrossRef]
  18. Comelli, N.C.; Duchowicz, P.R.; Castro, E.A. QSAR models for thiophene and imidazopyridine derivatives inhibitors of the Polo-Like Kinase 1. Eur. J. Pharm. Sci. 2014, 62, 171–179. [Google Scholar] [CrossRef] [PubMed]
  19. Chekkara, R.; Kandakatla, N.; Gorla, V.R.; Tenkayala, S.R.; Susithra, E. Theoretical studies on benzimidazole and imidazo [1, 2-a] pyridine derivatives as Polo-like kinase 1 (Plk1) inhibitors: Pharmacophore modeling, atom-based 3D-QSAR and molecular docking approach. J. Saudi Chem. Soc. 2017, 21, S311–S321. [Google Scholar] [CrossRef]
  20. Kong, Y.; Yan, A. QSAR models for predicting the bioactivity of Polo-Like Kinase 1 inhibitors. Chemom. Intell. Lab. Syst. 2017, 167, 214–225. [Google Scholar] [CrossRef]
  21. Duchowicz, P.R.; Comelli, N.C.; Ortiz, E.V.; Castro, E.A. QSAR study for carcinogenicity in a large set of organic compounds. Curr. Drug Saf. 2012, 7, 282–288. [Google Scholar] [CrossRef] [PubMed]
  22. Duchowicz, P.R.; Bennardi, D.O.; Baselo, D.E.; Bonifazi, E.L.; Rios-Luci, C.; Padrón, J.M.; Burton, G.; Misico, R.I. QSAR on Antiproliferative Naphthoquinones Based on a Conformation-Independent Approach. Eur. J. Med. Chem. 2014, 77, 176–184. [Google Scholar] [CrossRef] [PubMed]
  23. Ortiz, E.V.; Bennardi, D.O.; Bacelo, D.E.; Fioressi, S.E.; Duchowicz, P.R. The conformation-independent QSPR approach for predicting the oxidation rate constant of water micropollutants. Environ. Sci. Pollut. Res. 2017. [Google Scholar] [CrossRef] [PubMed]
  24. Aranda, J.F.; Bacelo, D.E.; Leguizamón Aparicio, M.S.; Ocsachoque, M.A.; Castro, E.A.; Duchowicz, P.R. Predicting the Bioconcentration Factor through a Conformation-Independent QSPR Study. SAR&QSAR Environ. Res. 2017, 28, 749–763. [Google Scholar]
  25. Kothe, M.; Kohls, D.; Low, S.; Coli, R.; Cheng, A.C.; Jacques, S.L.; Johnson, T.L.; Lewis, C.; Loh, C.; Nonomiya, J.; et al. Structure of the catalytic domain of human polo-like kinase 1. Biochemistry 2007, 46, 5960–5971. [Google Scholar] [CrossRef] [PubMed]
  26. Kubinyi, H. (Ed.) Comparative Molecular Field Analysis (CoMFA); Wiley-VCH Verlag GmbH: Weinheim, Germany, 2003. [Google Scholar]
  27. Gaulton, A.; Bellis, L.J.; Bento, A.P.; Chambers, J.; Davies, M.; Hersey, A.; Light, Y.; McGlinchey, S.; Michalovich, D.; Al-Lazikani, B. ChEMBL: A large-scale bioactivity database for drug discovery. Nucleic Acids Res. 2012, 40, D1100–D1107. [Google Scholar] [CrossRef] [PubMed]
  28. ChEMBL. Available online: https://www.ebi.ac.uk/chembl/ (accessed on 24 November 2017).
  29. Open Babel for Windows. Available online: http://openbabel.org/wiki/Category:Installation (accessed on 24 November 2017).
  30. ACD/ChemSketch. Available online: www.acdlabs.com (accessed on 24 November 2017).
  31. Pharmaceutical Data Exploration Laboratory (PaDEL). Available online: http://www.yapcwsoft.com/ (accessed on 24 November 2017).
  32. Yap, C.W. PaDEL-Descriptor: An open source software to calculate molecular descriptors and fingerprints. J. Comput. Chem. 2011, 32, 1466–1474. [Google Scholar] [CrossRef] [PubMed]
  33. Hong, H.; Xie, Q.; Ge, W.; Qian, F.; Fang, H.; Shi, L.; Su, Z.; Perkins, R.; Tong, W. Mold2, molecular descriptors from 2D structures for chemoinformatics and toxicoinformatics. J. Chem. Inf. Model. 2008, 48, 1337–1344. [Google Scholar] [CrossRef] [PubMed]
  34. Valdes-Martini, J.R.; García Jacas, C.R.; Marrero-Ponce, Y.; Silveira Vaz‘d Almeida, Y.; Morrel, C. QuBiLS-MAS: Free Software for Molecular Descriptors Calculator from Quadratic, Bilinear and Linear Maps Based on Graph–Theoretic Electronic-Density Matrices and Atomic Weightings; Version 1.0; CAMD-BIR Unit, CENDA Number of Register: 2373-2012; Central University of Las Villas: Villa Clara, Cuba, 2012. [Google Scholar]
  35. Duchowicz, P.R.; Castro, E.A.; Fernández, F.M. Alternative Algorithm for the Search of an Optimal Set of Descriptors in QSAR-QSPR Studies. MATCH Commun. Math. Comput. Chem. 2006, 55, 179–192. [Google Scholar]
  36. Morales, A.H.; Duchowicz, P.R.; Cabrera Pérez, M.A.; Castro, E.A.; Cordeiro, M.N.D.S.; González, M.P. Application of the replacement method as a novel variable selection strategy in QSAR. 1. Carcinogenic potential. Chemom. Intell. Lab. Syst. 2006, 81, 180–187. [Google Scholar] [CrossRef]
  37. Rojas, C.; Duchowicz, P.R.; Tripaldi, P.; Pis Diez, R. Quantitative structure-property relationship analysis for the retention index of fragrance-like compounds on a polar stationary phase. J. Chromatogr. A 2015, 1422, 277–288. [Google Scholar] [CrossRef] [PubMed]
  38. Duchowicz, P.R.; Fioressi, S.E.; Castro, E.A.; Wróbel, K.; Ibezim, N.E.; Bacelo, D.E. Conformation-independent QSAR study on human epidermal growth factor receptor-2 (HER2) inhibitors. Chem. Sel. 2017, 2, 3725–3731. [Google Scholar] [CrossRef]
  39. Kaufman, L.; Rousseeuw, P.J. Finding Groups in Data: An Introduction to Cluster Analysis; Wiley: New York, NY, USA, 2005. [Google Scholar]
  40. Golbraikh, A.; Tropsha, A. Beware of q2! J. Mol. Graph. Model. 2002, 20, 269–276. [Google Scholar] [CrossRef]
  41. Gramatica, P. Principles of QSAR models validation: Internal and external. QSAR Comb. Sci. 2007, 26, 694–701. [Google Scholar] [CrossRef]
  42. Roy, K.; Kar, S.; Ambure, P. On a simple approach for determining applicability domain of QSAR models. Chemom. Intell. Lab. Syst. 2015, 145, 22–29. [Google Scholar] [CrossRef]
  43. Eriksson, L.; Jaworska, J.; Worth, A.P.; Cronin, M.T.; McDowell, R.M.; Gramatica, P. Methods for Reliability and Uncertainty Assessment and for Applicability Evaluations of Classification- and Regression-Based QSARs. Environ. Health Perspect. 2003, 111, 1361–1375. [Google Scholar] [CrossRef] [PubMed]
  44. Rücker, C.; Rücker, G.; Meringer, M. Y-Randomization and its variants in QSPR/QSAR. J. Chem. Inf. Model. 2007, 47, 2345–2357. [Google Scholar] [CrossRef] [PubMed]
  45. Roy, K.; Roy, P.P. Comparative chemometric modeling of cytochrome 3A4 inhibitory activity of structurally diverse compounds using stepwise MLR, FAMLR, PLS, GFA, G/PLS and ANN techniques. Eur. J. Med. Chem. 2009, 44, 2913–2922. [Google Scholar] [CrossRef] [PubMed]
  46. Gramatica, P.; Sangion, A. A Historical Excursus on the Statistical Validation Parameters for QSAR Models: A Clarification Concerning Metrics and Terminology. J. Chem. Inf. Model. 2016, 56, 1127–1131. [Google Scholar] [CrossRef] [PubMed]
  47. Roy, K.; Das, R.N.; Ambure, P.; Aher, R.B. Be aware of error measures. Further studies on validation of predictive QSAR models. Chemom. Intell. Lab. Syst. 2016, 152, 18–33. [Google Scholar] [CrossRef]
  48. Cooper, J.; Saracci, R.; Cole, P. Describing the validity of carcinogen screening tests. Br. J. Cancer 1979, 39, 87–89. [Google Scholar] [CrossRef] [PubMed]
  49. Hoffmann, R.; Minkin, V.I.; Carpenter, B.K. Ockham’s Razor and Chemistry. Bull. Soc. Chim. Fr. 1996, 133, 117–130. [Google Scholar]
Figure 1. Some polo-like kinase-1 (PLK1) inhibitors involved in current clinical trials [5].
Figure 1. Some polo-like kinase-1 (PLK1) inhibitors involved in current clinical trials [5].
Cells 07 00013 g001
Figure 2. Predicted and experimental log 10 I C 50 values according to the quantitative structure-activity relationship (QSAR) of Equation (1).
Figure 2. Predicted and experimental log 10 I C 50 values according to the quantitative structure-activity relationship (QSAR) of Equation (1).
Cells 07 00013 g002
Figure 3. Dispersion plot of residuals for Equation (1).
Figure 3. Dispersion plot of residuals for Equation (1).
Cells 07 00013 g003
Figure 4. Williams plot for Equation (1).
Figure 4. Williams plot for Equation (1).
Cells 07 00013 g004
Table 1. Molecular descriptors involved in the best linear regression quantitative structure-activity relationship (QSAR) models for polo-like kinase-1 (PLK1) inhibitors. The selected model appears in bold.
Table 1. Molecular descriptors involved in the best linear regression quantitative structure-activity relationship (QSAR) models for polo-like kinase-1 (PLK1) inhibitors. The selected model appears in bold.
dDescriptors R t r a i n 2 S t r a i n R v a l 2 S v a l R t e s t 2 S t e s t
1Sub990.311.180.391.250.281.31
2PC534; AP1700.491.020.561.080.521.06
3PC534; KR4261; AP1700.520.990.680.950.570.98
4nHBAcc3; PC534; KR4261; AP1700.570.940.710.900.620.93
5PC534; KR3577; KR4268; AP170; KRC38970.610.900.710.890.710.83
6maxHCsats; M66; PC534; KR3577; KR4268; KRC38970.640.870.740.850.690.84
7maxHCsats; M66; PC534; PC686; KR3577; KR4268; AP1590.660.840.740.840.660.89
8mindssC; maxHCsats; M66; PC494; PC534; PC686; KR3577; KR42680.690.800.750.820.690.85
9mindssC; maxHCsats; M66; PC494; PC534; PC686; KR3577; KR4268; APC5100.700.790.750.820.700.85

Share and Cite

MDPI and ACS Style

Duchowicz, P.R. Linear Regression QSAR Models for Polo-Like Kinase-1 Inhibitors. Cells 2018, 7, 13. https://doi.org/10.3390/cells7020013

AMA Style

Duchowicz PR. Linear Regression QSAR Models for Polo-Like Kinase-1 Inhibitors. Cells. 2018; 7(2):13. https://doi.org/10.3390/cells7020013

Chicago/Turabian Style

Duchowicz, Pablo R. 2018. "Linear Regression QSAR Models for Polo-Like Kinase-1 Inhibitors" Cells 7, no. 2: 13. https://doi.org/10.3390/cells7020013

APA Style

Duchowicz, P. R. (2018). Linear Regression QSAR Models for Polo-Like Kinase-1 Inhibitors. Cells, 7(2), 13. https://doi.org/10.3390/cells7020013

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop