SMILES-based quantitative structure–property relationships for half-wave potential of N-benzylsalicylthioamides

A Graph Representation Approach Based on Light Gradient Boosting Machine for Predicting Drug-Disease Associations

Determining the association between drug and disease is important in drug development. However, existing approaches for drug-disease associations (DDAs) prediction are too homogeneous in terms of feature extraction. Here, a novel graph representation approach based on light gradient boosting machine (GRLGB) is proposed for prediction of DDAs. After the introduction of the protein into a heterogeneous network, nodes features were extracted from two perspectives: network topology and biological knowledge. Finally, the GRLGB classifier was applied to predict potential DDAs. GRLGB achieved satisfactory results on Bdataset and Fdataset through 10-fold cross-validation. To further prove the reliability of the GRLGB, case studies involving anxiety disorders and clozapine were conducted. The results suggest that GRLGB can identify novel DDAs.

Medicinal Chemistry Meets Electrochemistry: Redox Potential in the Role of Endpoint or Molecular Descriptor in QSAR/QSPR

Many biochemical reactions are based on redox reactions. Therefore, the redox potential of a chemical compound may be related to its therapeutic or physiological effects. The study of redox properties of compounds is a domain of electrochemistry. The subject of this review is the relationship between electrochemistry and medicinal chemistry, with a focus on quantifying these relationships. A summary of the relevant achievements in the correlation between redox potential and structure, therapeutic activity, resp., is presented. The first part of the review examines the applicability of QSPR for the prediction of redox properties of medically important compounds. The second part brings the exhaustive review of publications using redox potential as a molecular descriptor in QSAR of biological activity. Despite the complexity of medicinal chemistry and biological reactions, it is possible to employ redox potential in QSAR/QSPR. In many cases, this electrochemical parameter plays an essential but rarely absolute role.

The Monte Carlo Method as a Tool to Build up Predictive QSPR/QSAR

Background:

The Monte Carlo method has a wide application in various scientific researches. For the development of predictive models in a form of the quantitative structure-property / activity relationships (QSPRs/QSARs), the Monte Carlo approach also can be useful. The CORAL software provides the Monte Carlo calculations aimed to build up QSPR/QSAR models for different endpoints.

Methods:

Molecular descriptors are a mathematical function of so-called correlation weights of various molecular features. The numerical values of the correlation weights give the maximal value of a target function. The target function leads to a correlation between endpoint and optimal descriptor for the visible training set. The predictive potential of the model is estimated with the validation set, i.e. compounds that are not involved in the process of building up the model.

Results:

The approach gave quite good models for a large number of various physicochemical, biochemical, ecological, and medicinal endpoints. Bibliography and basic statistical characteristics of several CORAL models are collected in the present review. In addition, the extended version of the approach for more complex systems (nanomaterials and peptides), where behaviour of systems is defined by a group of conditions besides the molecular structure is demonstrated.

Conclusion:

The Monte Carlo technique available via the CORAL software can be a useful and convenient tool for the QSPR/QSAR analysis.

Development of prediction model for fructose- 1,6- bisphosphatase inhibitors using the Monte Carlo method

Fructose-1,6-bisphosphatase (FBPase) is an enzyme important for regulation of gluconeogenesis, which is a major process in the liver responsible for glucose production. Inhibition of FBPase enzyme causing blockage of the gluconeogenesis process represents a newer scheme in the progress of anti-diabetic drugs. The current research describes the development of hybrid optimal descriptors-based quantitative structure-activity relationship (QSAR) models intended for a set of 62 FBPase inhibitors with the Monte Carlo method. The molecular structures were expressed by the simplified molecular input line entry system (SMILES) notation. Three splits were prepared by random division of the molecules into training set, calibration set and validation set. Statistical parameters obtained from QSAR modelling were good for various designed splits. The best QSAR model showed the following parameters: the values of r² for calibration set and validation set of the best model were 0.6837 and 0.8623 and of Q² were 0.6114 and 0.8036, respectively. Based on the results obtained for correlation weights, different structural attributes were described as promoter of the endpoint. Further, these structural attributes were used in designing of new FBPase inhibitors and a molecular docking study was completed for the determination of interactions of the designed molecules with the enzyme.

Heme Oxygenase Database (HemeOxDB) and QSAR Analysis of Isoform 1 Inhibitors

Due to increasing interest in the field of heme oxygenases (HOs), we built a ligand database called HemeOxDB that includes the entire set of known HO-1 and HO-2 inhibitors, resulting in more than 400 compounds. The HemeOxDB is available online at http://www.researchdsf.unict.it/hemeoxdb/, and having a robust search engine allows end users to build complex queries, sort tabulated results, and generate color-coded two- and three-dimensional graphs. This database will grow to be a tool for the design of potent and selective HO-1 or HO-2 inhibitors. We were also interested in virtually searching for alternative inhibitors, and, for the first time in the field of HOs, a quantitative structure-activity relationship (QSAR) model was built using half-maximal inhibitory concentration (IC₅₀ ) values of the whole set of known HO-1 inhibitors, taken from the HemeOxDB and employing the Monte Carlo technique. The statistical quality suggested that the model is robust and possesses desirable predictive potential. The screening of US Food and Drug Administration (FDA)-approved drugs, external to our dataset, suggested new predicted inhibitors, opening the way for replacing imidazole groups. The HemeOxDB and the QSAR model reported herein may help in prospectively identifying or repurposing new drugs with optimal structural attributes for HO enzyme inhibition.

QSPR/QSAR Analyses by Means of the CORAL Software

In this chapter, the methodology of building up quantitative structure—property/activity relationships (QSPRs/QSARs)—by means of the CORAL software is described. The Monte Carlo method is the basis of this approach. Simplified Molecular Input-Line Entry System (SMILES) is used as the representation of the molecular structure. The conversion of SMILES into the molecular graph is available for QSPR/QSAR analysis using the CORAL software. The model for an endpoint is a mathematical function of the correlation weights for various features of the molecular structure. Hybrid models that are based on features extracted from both SMILES and a graph also can be built up by the CORAL software. The conceptually new ideas collected and revealed through the CORAL software are: (1) any QSPR/QSAR model is a random event; and (2) optimal descriptor can be a translator of eclectic information into an endpoint prediction.

Searching for anthranilic acid-based thumb pocket 2 HCV NS5B polymerase inhibitors through a combination of molecular docking, 3D-QSAR and virtual screening

A combination of the following computational methods: (i) molecular docking, (ii) 3-D Quantitative Structure Activity Relationship Comparative Molecular Field Analysis (3D-QSAR CoMFA), (iii) similarity search and (iv) virtual screening using PubChem database was applied to identify new anthranilic acid-based inhibitors of hepatitis C virus (HCV) replication. A number of known inhibitors were initially docked into the "Thumb Pocket 2" allosteric site of the crystal structure of the enzyme HCV RNA-dependent RNA polymerase (NS5B GT1b). Then, the CoMFA fields were generated through a receptor-based alignment of docking poses to build a validated and stable 3D-QSAR CoMFA model. The proposed model can be first utilized to get insight into the molecular features that promote bioactivity, and then within a virtual screening procedure, it can be used to estimate the activity of novel potential bioactive compounds prior to their synthesis and biological tests.

QSAR as a random event: a case of NOAEL

Quantitative structure-activity relationships (QSAR) for no observed adverse effect levels (NOAEL, mmol/kg/day, in logarithmic units) are suggested. Simplified molecular input line entry systems (SMILES) were used for molecular structure representation. Monte Carlo method was used for one-variable models building up for three different splits into the "visible" training set and "invisible" validation. The statistical quality of the models for three random splits are the following: split 1 n = 180, r (2) = 0.718, q (2) = 0.712, s = 0.403, F = 454 (training set); n = 17, r (2) = 0.544, s = 0.367 (calibration set); n = 21, r (2) = 0.61, s = 0.44, r m (2) = 0.61 (validation set); split 2 n = 169, r (2) = 0.711, q (2) = 0.705, s = 0.409, F = 411 (training set); n = 27, r (2) = 0.512, s = 0.461 (calibration set); n = 22, r (2) = 0.669, s = 0.360, r m (2) = 0.63 (validation set); split 3 n = 172, r (2) = 0.679, q (2) = 0.672, s = 0.420, F = 360 (training set); n = 19, r (2) = 0.617, s = 0.582 (calibration set); n = 21, r (2) = 0.627, s = 0.367, r m (2) = 0.54 (validation set). All models are built according to OCED principles.

Simplified molecular input line entry system-based: QSAR modelling for MAP kinase-interacting protein kinase (MNK1)

Quantitative structure-activity relationship (QSAR) models were built for the prediction of inhibition (pIC50, i.e. negative logarithm of the 50% effective concentration) of MAP kinase-interacting protein kinase (MNK1) by 43 potent inhibitors. The pIC50 values were modelled with five random splits, with the representations of the molecular structures by simplified molecular input line entry system (SMILES). QSAR model building was performed by the Monte Carlo optimisation using three methods: classic scheme; balance of correlations; and balance correlation with ideal slopes. The robustness of these models were checked by parameters as rm(2), r(*)m(2), [Formula: see text] and randomisation technique. The best QSAR model based on single optimal descriptors was applied to study in vitro structure-activity relationships of 6-(4-(2-(piperidin-1-yl) ethoxy) phenyl)-3-(pyridin-4-yl) pyrazolo [1,5-a] pyrimidine derivatives as a screening tool for the development of novel potent MNK1 inhibitors. The effects of alkyl group, -OH, -NO2, F, Cl, Br, I, etc. on the IC50 values towards the inhibition of MNK1 were also reported.

Quasi-QSAR for mutagenic potential of multi-walled carbon-nanotubes

Available on the Internet, the CORAL software (http://www.insilico.eu/coral) has been used to build up quasi-quantitative structure-activity relationships (quasi-QSAR) for prediction of mutagenic potential of multi-walled carbon-nanotubes (MWCNTs). In contrast with the previous models built up by CORAL which were based on representation of the molecular structure by simplified molecular input-line entry system (SMILES) the quasi-QSARs based on the representation of conditions (not on the molecular structure) such as concentration, presence (absence) S9 mix, the using (or without the using) of preincubation were encoded by so-called quasi-SMILES. The statistical characteristics of these models (quasi-QSARs) for three random splits into the visible training set and test set and invisible validation set are the following: (i) split 1: n=13, r(2)=0.8037, q(2)=0.7260, s=0.033, F=45 (training set); n=5, r(2)=0.9102, s=0.071 (test set); n=6, r(2)=0.7627, s=0.044 (validation set); (ii) split 2: n=13, r(2)=0.6446, q(2)=0.4733, s=0.045, F=20 (training set); n=5, r(2)=0.6785, s=0.054 (test set); n=6, r(2)=0.9593, s=0.032 (validation set); and (iii) n=14, r(2)=0.8087, q(2)=0.6975, s=0.026, F=51 (training set); n=5, r(2)=0.9453, s=0.074 (test set); n=5, r(2)=0.8951, s=0.052 (validation set).

QSPR/QSAR Analyses by Means of the CORAL Software

In this chapter, the methodology of building up quantitative structure—property/activity relationships (QSPRs/QSARs)—by means of the CORAL software is described. The Monte Carlo method is the basis of this approach. Simplified Molecular Input-Line Entry System (SMILES) is used as the representation of the molecular structure. The conversion of SMILES into the molecular graph is available for QSPR/QSAR analysis using the CORAL software. The model for an endpoint is a mathematical function of the correlation weights for various features of the molecular structure. Hybrid models that are based on features extracted from both SMILES and a graph also can be built up by the CORAL software. The conceptually new ideas collected and revealed through the CORAL software are: (1) any QSPR/QSAR model is a random event; and (2) optimal descriptor can be a translator of eclectic information into an endpoint prediction.