Pavel G. Polishchuk

Existing and Developing Approaches for QSAR Analysis of Mixtures

This review is devoted to the critical analysis of advantages and disadvantages of existing mixture descriptors and their usage in various QSAR/QSPR tasks. We describe good practices for the QSAR modeling of mixtures, data sources for mixtures, a discussion of various mixture descriptors and their application, recommendations about proper external validation specific for mixture QSAR modeling, and future perspectives of this field. The biggest problem in QSAR of mixtures is the lack of reliable data about the mixtures’ properties. Various mixture descriptors are used for the modeling of different endpoints. However, these descriptors have certain disadvantages, such as applicability only to 1 : 1 binary mixtures, and additive nature. The field of QSAR of mixtures is still under development, and existing efforts could be considered as a foundation for future approaches and studies. The usage of non‐additive mixture descriptors, which are sensitive to interaction effects, in combination with best practices of QSAR model development (e.g., thorough data collection and curation, rigorous external validation, etc.) will significantly improve the quality of QSAR studies of mixtures.

Interpretation of QSAR Models Based on Random Forest Methods.

A new algorithm for the interpretation of Random Forest models has been developed. It allows to calculate the contribution of each descriptor to the calculated property value. In case of the simplex representation of a molecular structure, contributions of individual atoms can be calculated, and thus it becomes possible to estimate the influence of separate molecular fragments on the investigated property. Such information can be used for the design of new compounds with a predefined property value. The proposed measure of descriptor contributions is not an alternative to the importance of Breiman’s variable, but it characterizes the contribution of a particular explanatory variable to the calculated response value.

QSPR Approach to Predict Nonadditive Properties of Mixtures. Application to Bubble Point Temperatures of Binary Mixtures of Liquids

This paper is devoted to the development of methodology for QSPR modeling of mixtures and its application to vapor/liquid equilibrium diagrams for bubble point temperatures of binary liquid mixtures. Two types of special mixture descriptors based on SiRMS and ISIDA approaches were developed. SiRMS‐based fragment descriptors involve atoms belonging to both components of the mixture, whereas the ISIDA fragments belong only to one of these components. The models were built on the data set containing the phase diagrams for 167 mixtures represented by different combinations of 67 pure liquids. Consensus models were developed using nonlinear Support Vector Machine (SVM), Associative Neural Networks (ASNN), and Random Forest (RF) approaches. For SVM and ASNN calculations, the ISIDA fragment descriptors were used, whereas Simplex descriptors were employed in RF models. The models have been validated using three different protocols: “Points out”, “Mixtures out” and “Compounds out”, based on the specific rules to form training/test sets in each fold of cross‐validation. A final validation of the models has been performed on an additional set of 94 mixtures represented by combinations of novel 34 compounds and modeling set chemicals with each other. The root mean squared error of predictions for new mixtures of already known liquids does not exceed 5.7 K, which outperforms COSMO‐RS models. Developed QSAR methodology can be applied to the modeling of any nonadditive property of binary mixtures (antiviral activities, drug formulation, etc.)

Universal Approach for Structural Interpretation of QSAR/QSPR Models

In this paper we offer a novel approach for the structural interpretation of QSAR models. The major advantage of our developed methodology is its universality, i.e., it can be applied to any QSAR/QSPR model irrespective of chemical descriptors and machine learning methods applied. This universality was achieved by using only the information obtained from substructures of the compounds of interest to interpret model outcomes. Reliability of the offered approach was confirmed by the results of three case studies, including end‐points of different types (continuous and binary classification) and nature (solubility, mutagenicity, and inhibition of Transglutaminase 2), various fragment and whole‐molecule descriptors (Simplex and Dragon), and multiple modeling techniques (partial least squares, random forest, and support vector machines). We compared the global contributions of molecular fragments obtained using our methodology with known SAR rules derived experimentally. In all cases high concordance between our interpretation and results published by others was observed. Although the proposed interpretation approach could be easily extended to any type of descriptors, we would recommend using Simplex descriptors to achieve a larger variety of investigated molecular fragments. The developed approach is a good tool for interpretation of such “black box” models like random forest, neural networks, etc. Analysis of fragment global contributions and their deviation across a dataset could be useful for the identification of key fragments and structural alerts. This information could be helpful to maximize the positive influence of structural surroundings on the given fragment and to decrease the negative effects.

Structure-Reactivity Relationships in Terms of the Condensed Graphs of Reactions

An approach for the prediction of rate constants of chemical reactions, based on the representation of a chemical reaction as a condensed graph, has been tested on more than 1000 bimolecular nucleophilic substitution reactions with neutral nucleophiles in 38 solvents. Molecular fragment descriptors, temperature, and solvent parameters characterizing solvation power have been used in the reaction modeling. The obtained models ensure a good correlation between the predicted and experimental values; the corresponding deviations are comparable with interlaboratory measurement errors.

Application of Random Forest and Multiple Linear Regression Techniques to QSPR Prediction of an Aqueous Solubility for Military Compounds

The relationship between the aqueous solubility of more than two thousand eight hundred organic compounds and their structures was investigated using a QSPR approach based on Simplex Representation of Molecular Structure (SiRMS). The dataset consists of 2537 diverse organic compounds. Multiple Linear Regression (MLR) and Random Forest (RF) methods were used for statistical modeling at the 2D level of representation of molecular structure. Statistical characteristics of the best models are quite good (MLR method: R2=0.85, Q2=0.83; RF method: R2=0.99, R2oob=0.88). The external validation set of 301 compounds (including 47 nitro‐, nitroso‐ and nitrogen‐rich compounds of military interest) which were not included in the training set and modeling process, was used for evaluation of the models predictivity. Thus, well‐fitted and robust (R2test(MLR)=0.76 and R2test(RF)=0.82) models were obtained for both statistical techniques using descriptors based on the topological structural information only. The predicted solubility values for military compounds are in good agreement with experimental ones. Developed QSPR models represent powerful and easy‐to‐use virtual screening tool that can be recommended for prediction of aqueous solubility.

Design, Virtual Screening, and Synthesis of Antagonists of αIIbβ3 as Antiplatelet Agents

This article describes design, virtual screening, synthesis, and biological tests of novel αIIbβ3 antagonists, which inhibit platelet aggregation. Two types of αIIbβ3 antagonists were developed: those binding either closed or open form of the protein. At the first step, available experimental data were used to build QSAR models and ligand- and structure-based pharmacophore models and to select the most appropriate tool for ligand-to-protein docking. Virtual screening of publicly available databases (BioinfoDB, ZINC, Enamine data sets) with developed models resulted in no hits. Therefore, small focused libraries for two types of ligands were prepared on the basis of pharmacophore models. Their screening resulted in four potential ligands for open form of αIIbβ3 and four ligands for its closed form followed by their synthesis and in vitro tests. Experimental measurements of affinity for αIIbβ3 and ability to inhibit ADP-induced platelet aggregation (IC50) showed that two designed ligands for the open form 4c and 4d (IC50 = 6.2 nM and 25 nM, respectively) and one for the closed form 12b (IC50 = 11 nM) were more potent than commercial antithrombotic Tirofiban (IC50 = 32 nM).

Structural and Physico-Chemical Interpretation (SPCI) of QSAR Models and Its Comparison with Matched Molecular Pair Analysis

This paper describes the Structural and Physico-Chemical Interpretation (SPCI) approach, which is an extension of a recently reported method for interpretation of quantitative structure-activity relationship (QSAR) models. This approach can efficiently be used to reveal structural motifs and the major physicochemical factors affecting the investigated properties. Its efficacy was demonstrated both on the classical Free-Wilson data set and on several data sets with different end points (permeability of the blood-brain barrier, fibrinogen receptor antagonists, acute oral toxicity). Structure-activity patterns extracted from QSAR models with SPCI were in good correspondence with experimentally observed relationships and molecular docking, regardless of the machine learning method used. Comparison of SPCI with the matched molecular pair (MMP) method clearly shows an advantage of our approach over MMP, especially for small or structurally diverse data sets. The developed approach has been implemented in the SPCI software tool with a graphical user interface, which is publicly available at http://qsar4u.com/pages/sirms_qsar.php .

Synthesis, biological evaluation, X-ray molecular structure and molecular docking studies of RGD mimetics containing 6-amino-2,3-dihydroisoindolin-1-one fragment as ligands of integrin αIIbβ3

A series of novel RGD mimetics containing phthalimidine fragment was designed and synthesized. Their antiaggregative activity determined by Borns method was shown to be due to inhibition of fibrinogen binding to αIIbβ₃. Molecular docking of RGD mimetics to αIIbβ₃ receptor showed the key interactions in this complex, and also some correlations have been observed between values of biological activity and docking scores. The single crystal X-ray data were obtained for five mimetics.

Structure–reactivity modeling using mixture-based representation of chemical reactions

We describe a novel approach of reaction representation as a combination of two mixtures: a mixture of reactants and a mixture of products. In turn, each mixture can be encoded using an earlier reported approach involving simplex descriptors (SiRMS). The feature vector representing these two mixtures results from either concatenated product and reactant descriptors or the difference between descriptors of products and reactants. This reaction representation doesn’t need an explicit labeling of a reaction center. The rigorous “product-out” cross-validation (CV) strategy has been suggested. Unlike the naïve “reaction-out” CV approach based on a random selection of items, the proposed one provides with more realistic estimation of prediction accuracy for reactions resulting in novel products. The new methodology has been applied to model rate constants of E2 reactions. It has been demonstrated that the use of the fragment control domain applicability approach significantly increases prediction accuracy of the models. The models obtained with new “mixture” approach performed better than those required either explicit (Condensed Graph of Reaction) or implicit (reaction fingerprints) reaction center labeling.