In silico design of diacylglycerol acyltransferase-1 (DGAT1) inhibitors based on SMILES descriptors using Monte-Carlo method

QSAR study of tetrahydropteridin derivatives as polo-like kinase 1(PLK1) Inhibitors with molecular docking and dynamics study

PLK1 is the key target for dealing with different cancer because it plays an important role in cell proliferation. According to the regulation of OECD, a QSAR model was developed from a dataset of 68 tetrahydropteridin derivatives. Three descriptors (maxHaaCH, ATSC7i, AATS7m) were considered for the development of the QSAR model. The reliability and predictability of the developed QSAR model were evaluated by various statistical parameters (r² = 0.8213, r²_ext = 0.8771 and CCC_ext = 0.9364). The maxHaaCH descriptor is positively correlated to pIC₅₀ whereas, the ATSC7i and AATS7m are negatively correlated with pIC₅₀. The QSAR model explains all the structural features and shows a good correlation with the activity. Based on molecular modelling techniques, five compounds (D1-D5) were designed. Molecular docking and dynamics studies of the most active compound were performed with PDB ID: 2RKU. The results of the present investigation may be employed to identify and develop effective inhibitors for the treatment of PLK1-related pathophysiological disorders.

Development of Novel Therapeutics for Schizophrenia Treatment Based on a Selective Positive Allosteric Modulation of α1-Containing GABAARs-In Silico Approach

For the development of atypical antipsychotics, the selective positive allosteric modulation of the ionotropic GABAA receptor (GABAAR) has emerged as a promising approach. In the presented research, two unrelated methods were used for the development of QSAR models for selective positive allosteric modulation of 1-containing GABAARs with derivatives of imidazo [1,2-a]-pyridine. The development of conformation-independent QSAR models, based on descriptors derived from local molecular graph invariants and SMILES notation, was achieved with the Monte Carlo optimization method. From the vast pool of 0D, 1D, and 2D molecule descriptors, the GA-MLR method developed additional QSAR models. Various statistical methods were utilised for the determination of the developed models’ robustness, predictability, and overall quality, and according to the obtained results, all QSAR models are considered good. The molecular fragments that have a positive or negative impact on the studied activity were obtained from the studied molecules’ SMILES notations, and according to the obtained results, nine novel compounds were designed. The binding affinities to GABAAR of designed compounds were assessed with the application of molecular docking studies and the obtained results showed a high correlation with results obtained from QSAR modeling. To assess all designed molecules’ “drug-likeness”, their physicochemical descriptors were computed and utilised for the prediction of medicinal chemistry friendliness, pharmacokinetic properties, ADME parameters, and druglike nature.

CORAL: Development of a hybrid descriptor based QSTR model to predict the toxicity of dioxins and dioxin-like compounds with correlation intensity index and consensus modelling

In the present study, ninety-five halogenated dioxins and related chemicals (dibenzo-p-dioxins, dibenzofurans, biphenyls, and naphthalene) with endpoint pEC₅₀ were used to develop twelve quantitative structure toxicity relationship (QSTR) models using inbuilt Monte Carlo algorithm of CORAL software. The hybrid optimal descriptor of correlation weights (DCW) using a combination of SMILES and HSG (hydrogen suppressed graph) was employed to generate QSTR models. Three target functions i.e. TF₁ (W_IIC=W_CII=0), TF₂ (W_IIC= 0.3 & W_CII=0) and TF₃ (W_IIC= 0.0 &W_CII=0.3) were employed to develop robust QSTR models and the statistical outcomes of each target function were compared with each other. The correlation intensity index (CII) was found a reliable benchmark of the predictive potential for QSTR models. The numerical value of the determination coefficient of the validation set of split 1 computed by TF₃ was found highest (R_Valid²=0.8438). The fragments responsible for the toxicity of dioxins and related chemicals were also identified in terms of the promoter of increase/decrease for pEC₅₀. Three random splits (Split 1, Split 2 and Split 4) were selected for the extraction of the promoter of increase/decrease for pEC₅₀. In the last, consensus modelling was performed using the intelligent consensus tool of DTC lab (https://dtclab.webs.com/software-tools). The original consensus model, which was created by combining four distinct models employing the split 4 arrangement, was more predictive for the validation set and the numerical value of the determination coefficient of the test set (validation set) was increased from 0.8133 to 0.9725. For the validation set of split 4, the mean absolute error (MAE 100%) was also lowered from 0.513 to 0.2739.

Carcinogenicity prediction using the index of ideality of correlation

Carcinogenicity testing is necessary to protect human health and comply with regulations, but testing it with the traditionally used two-year rodent studies is time-consuming and expensive. In certain cases, such as for impurities, alternative methods may be convenient. Thus there is an urgent need for alternative approaches for reliable and robust assessments of carcinogenicity. The Monte Carlo technique with CORAL software is a tool to tackle this task for unknown compounds using available experimental data for a representative set of compounds. The models can be constructed with the simplified molecular input line entry system without additional physicochemical descriptors. We describe here a model based on a data set of 1167 substances. Matthew's correlation coefficient values for calibration and validation sets are 0.747 and 0.577, respectively. Double bonds between carbon atoms and double bonds of oxygen atoms are the molecular features that indicate the carcinogenic potential of a compound.

Development of novel antipsychotic agents by inhibiting dopamine transporter – in silico approach

Various in silico methods were employed for the development of antipsychotic agents by dopamine transporter inhibition.

Quantitative structure-toxicity relationship models for predication of toxicity of ionic liquids towards Leukemia rat cell line IPC-81 based on index of ideality of correlation

The application of ion liquids (ILs) as green solvents has attracted the attention of the scientific community. However, ILs may play the role of toxins. Even though ionic liquids may assist to minimise air pollution, but their discharge into aquatic ecosystems might result in significant water pollution due to their potential toxicity and inaccessibility to biodegradation. Recently, more attention has been paid to the toxicity of ILs on plants, bacteria, and humans. Here, a quantitative structure-toxicity relationship study (QSTR) based on the Monte Carlo method of CORAL software has been applied to estimate the logarithm of the half-maximal effective concentration of toxicity of ILs against leukemia rat cell line IPC-81 (logEC₅₀). A hybrid optimal descriptor is used to build QSTR models for a large set of 304 diverse ILs including ammonium, imidazolium, morpholinium, phosphonium, piperidinium, pyridinium, pyrrolidinium, quinolinium, sulfonium, and protic ILs. The SMILES notations of Ils are utilized to compute the descriptor correlation weight (DCW). Four splits are made from the whole dataset and each split is randomly divided into four sets (training subsets and validation set). The index of ideality of correlation (IIC) is applied to evaluate the authenticity and robustness of the QSTR models. A QSTR model with statistical parameters R²=0.85, CCC =0.92, Q²=0.84, and MAE =0.25 for the validation set of the best split is considered as a prime model. The outliers and promoters of increase/decrease of logEC₅₀ are extracted and the mechanistic interpretation of effective descriptors for the model is also offered.HighlightsGlobal SMILES-based QSAR model was developed to predict the toxicity of ILs.The CORAL software is used to model the ILs toxicity on IPC-81 leukemia rat cell lineIIC is tested as a criterion of predictive potential.The toxicological effects of ILs are discussed based on the proposed model.

The Monte Carlo approach to model and predict the melting point of imidazolium ionic liquids using hybrid optimal descriptors

Ionic liquids (ILs) have captured intensive attention owing to their unique properties such as high thermal stability, negligible vapour pressure, high dissolution capacity and high ionic conductivity as well as their wide applications in various scientific fields including organic synthesis, catalysis, and industrial extraction processes. Many applications of ionic liquids (ILs) rely on the melting point (T_m). Therefore, in the present manuscript, the melting points of imidazolium ILs are studied employing a quantitative structure–property relationship (QSPR) approach to develop a model for predicting the melting points of a data set of imidazolium ILs. The Monte Carlo algorithm of CORAL software is applied to build up a robust QSPR model to calculate the values T_m of 353 imidazolium ILs. Using a combination of SMILES and hydrogen-suppressed molecular graphs (HSGs), the hybrid optimal descriptor is computed and used to generate the QSPR models. Internal and external validation parameters are also employed to evaluate the predictability and reliability of the QSPR model. Four splits are prepared from the dataset and each split is randomly distributed into four sets i.e. training set (≈33%), invisible training set (≈31%), calibration set (≈16%) and validation set (≈20%). In QSPR modelling, the numerical values of various statistical features of the validation sets such as R_Validation², Q_Validation², and IIC_Validation are found to be in the range of 0.7846–0.8535, 0.7687–0.8423 and 0.7424–0.8982, respectively. For mechanistic interpretation, the structural attributes which are responsible for the increase/decrease of T_m are also extracted.

The melting points of imidazolium ILs are studied employing a quantitative structure–property relationship (QSPR) approach to develop a model for predicting the melting points of a data set of imidazolium ILs.

Prediction of power conversion efficiency of phenothiazine-based dye-sensitized solar cells using Monte Carlo method with index of ideality of correlation

Simplified molecular-input line-entry system (SMILES) notation and inbuilt Monte Carlo algorithm of CORAL software were employed to construct generative and prediction QSPR models for the analysis of the power conversion efficiency (PCE) of 215 phenothiazine derivatives. The dataset was divided into four splits and each split was further divided into four sets. A hybrid descriptor, a combination of SMILES and hydrogen suppressed graph (HSG), was employed to build reliable and robust QSPR models. The role of the index of ideality of correlation (IIC) was also studied in depth. We performed a comparative study to predict PCE using two target functions (TF₁ without IIC and TF₂ with IIC). Eight QSPR models were developed and the models developed with TF₂ was shown robust and reliable. The QSPR model generated from split 4 was considered a leading model. The different statistical benchmarks were computed for the lead model and these were rtraining set2=0.7784; rinvisible training set2=0.7955; rcalibration set2=0.7738; rvalidation set2=0.7506; Qtraining set2=0.7691; Qinvisible training set2=0.7850; Qcalibration set2=0.7501;  Qvalidation set2=0.7085; IIC_{training set} = 0.8590; IIC_{invisible training set} = 0.8297; IIC_{calibration set} = 0.8796; IIC_{validation set} = 0.8293, etc. The promoters of increase and decrease of endpoint PCE were also extracted.

Semi-correlations as a tool to model for skin sensitization

Semi-correlation specifically assesses the correlation between a binary variable and a continuous variable. Semi-correlations were applied to develop binary models for various endpoints. We applied the semi-correlation to develop models of two kinds of skin sensitization one related to animals (local lymph node assay LLNA) and one to human beings (direct peptide reactivity assay DPRA and/or human cell line activation test h-CLAT). The models refer to binary classification for a two-level strategy: the first level (analysis of all compounds) is used in the format "sensitizer or non-sensitizer", and the second level (only sensitizers) is a further classification in the format "strong or weak sensitizer". The ranges of statistical characteristics of the models depend on the endpoint, LLNA or DPRA/h-CLAT: for the first level, sensitivity: 0.69-0.88, specificity: 0.75-0.89, accuracy: 0.77-0.87, Matthew's correlation coefficient (MCC): 0.54-0.57 and for the second level, sensitivity: 0.70-1.0, specificity: 0.78-0.83, accuracy: 0.77-0.87, MCC: 0.54-0.76. Thus, the described approach can be applied to building up models of the skin sensitization potency.

SMILES-based QSAR and molecular docking study of xanthone derivatives as α-glucosidase inhibitors

Increasing diabetic population is one of the major health concerns all over the world. Inhibition of α-glucosidase is a clinically proved and attractive strategy to manage diabetes. In this study, robust and reliable QSAR models to predict α-glucosidase inhibitory potential of xanthone derivatives are developed by the Monte Carlo technique. The chemical structures are represented by SMILES notation without any 3D-optimization. The significance of the index of ideality correlation (IIC) with applicability domain (AD) is also studied in depth. The models developed using CORAL software by considering IIC criteria are found to be statistically more significant and robust than simple balance of correlation. The QSAR models are validated by both internal and external validation methods. The promoters of increase and decrease of activity are also extracted and interpreted in detail. The interpretation of developed models explains the role of different structural attributes in predicting the pIC₅₀ of xanthone derivatives as α-glucosidase inhibitors. Based on the results of model interpretation, modifications are done on some xanthone derivatives and 15 new molecules are designed. The α-glucosidase inhibitory activity of novel molecules is further supported by docking studies.

The index of ideality of correlation: QSAR studies of hepatitis C virus NS3/4A protease inhibitors using SMILES descriptors

Robust and reliable QSAR models were developed to predict half-maximal inhibitory concentration (IC₅₀) values of hepatitis C virus NS3/4A protease inhibitors from the Monte Carlo technique. 524 HCV NS3/4A protease inhibitors were extracted from the scientific literature to create a reasonably large set. The models were developed using CORAL software by using two target functions namely target function 1 (TF1) without applying the index of ideality of correlation (IIC) and target function 2 (TF2) that uses IIC. The constructed models based on TF2 were statistically more significant and robust than the models based on TF1. The determination coefficients (r²) of training and test sets were 0.86 and 0.88 for the best split based on TF2. The promoters of the increase/decrease of activity were also extracted and interpreted in detail. The model interpretation results explain the role of different structural attributes in predicting the pIC₅₀ values of hepatitis C virus NS3/4A protease inhibitors. Based on the mechanistic model interpretation results, eight new compounds were designed and their pIC₅₀ values were predicted based on the average prediction of ten models.

The Monte Carlo method to build up models of the hydrolysis half-lives of organic compounds

The hydrolysis of organic chemicals such as pesticides, pollutants, or drugs can affect the fate and behaviour of environmental contaminants, so it is of interest to evaluate the stability of substances in water for various purposes. For the registration of organic compounds in Europe, information on hydrolysis must be presented. However, the experimental measurements of all chemicals would require enormous resources, and computational models may become attractive. Applying the CORAL software (http://www.insilico.eu/coral) quantitative structure-property relationships (QSPRs) were built up to model hydrolysis. The 2D-optimal descriptor is calculated with so-called correlation weights for attributes of simplified molecular input-line entry systems (SMILES). The correlation weights are obtained as results of the special Monte Carlo optimization. The nature of (five- and six-member) rings is an important component of this approach. Another important component is the atom pair proportions for nitrogen, oxygen, and sulphur. The statistical quality of the best model is: n = 44, r² = 0.74 (training set); n = 14, r² = 0.75 (calibration set); and n = 12, r² = 0.80 (validation set).

Quantitative structure activity relationship studies of novel hydrazone derivatives as α-amylase inhibitors with index of ideality of correlation

The present manuscript describes the synthesis, α-amylase inhibition, in silico studies and in-depth quantitative structure-activity relationship (QSAR) of a library of aroyl hydrazones based on benzothiazole skeleton. All the compounds of the developed library are characterized by various spectral techniques. α-Amylase inhibitory potential of all compounds has been explored, where compound 7n exhibits remarkable α-amylase inhibition of 87.5% at 50 µg/mL. Robust QSAR models are made by using the balance of correlation method in CORAL software. The chemical structures at different concentration with optimal descriptors are represented by SMILES. A data set of 66 SMILES of 22 hydrazones at three distinct concentrations are prepared. The significance of the index of ideality of correlation (IIC) with applicability domain (AD) is also studied at depth. A QSAR model with best Rvalidation2 = 0.8587 for split 1 is considered as a leading model. The outliers and promoters of increase and decrease of endpoint are also extracted. The binding modes of the most active compound, that is, 7n in the active site of Aspergillus oryzae α-amylase (PDB ID: 7TAA) are also explored by in silico molecular docking studies. Compound 7n displays high resemblance in binding mode and pose with the standard drug acarbose. Molecular dynamics simulations performed on protein-ligand complex for 100 ns, the protein gets stabilised after 20 ns and remained below 2 Å for the remaining simulation. Moreover, the deviation observed in RMSF during simulation for each amino acid residue with respect to Cα carbon atom is insignificant.

Prediction of No Observed Adverse Effect Concentration for inhalation toxicity using Monte Carlo approach

Ideal correlation is one variable model based on so-called optimal descriptors calculated with simplified molecular input-line entry systems (SMILES). The optimal descriptor is calculated according to the index of ideality of correlation, a new criterion of predictive potential of quantitative structure-property/activity relationships (QSPRs/QSARs). The aim of the present study was the building and estimation of models for inhalation toxicity as No Observed Adverse Effect Concentration (NOAEC) based on the OECD guidelines 413. Three random distributions into the training set and validation set were examined. In practice, a structured training set that contains active training set, passive training set and calibration set is used as the training set. The statistical characteristics of the best model for negative logarithm of NOAEC (pNOAEC) are for training set n = 108, average r ² = 0.52 + 0.62 + 0.76/3 = 0.63 and for validation set n = 35, r ² = 0.73.

A Monte Carlo method based QSPR model for prediction of reaction rate constants of hydrated electrons with organic contaminants

The Monte Carlo algorithm was applied to formulate a robust quantitative structure-property relationship (QSPR) model to compute the reactions rate constants of hydrated electron values for a data set of 309 water contaminants containing 125 aliphatic and 184 phenyl-based chemicals. The QSPR models were computed with the hybrid optimal descriptors which were procured by combining the SMILES and hydrogen-suppressed molecular graph for both classes of compounds. Approximately 75% of the total experimental data set was randomly divided into training and invisible training sets, while approximately 25% was divided into calibration and validation sets. The authenticity and robustness of the developed QSPR models were also judged by the Index of Ideality of Correlation. In QSPR modelling of aliphatic compounds, the numerical values of r T r a i n i n g 2 , r V a l i d a t i o n 2 , Q T r a i n i n g 2 and Q V a l i d a t i o n 2 were in the range of 0.852-0.905, 0.815-0.894, 0.839-0.897 and 0.737-0.867, respectively. Whereas, in the QSPR modelling of phenyl-based compounds, the numerical values of r T r a i n i n g 2 , r V a l i d a t i o n 2 , Q T r a i n i n g 2 and Q V a l i d a t i o n 2 were in the range of 0.867-0.896, 0.852-0.865, 0.816-0.850 and 0.760-0.762, respectively. The structural attributes, which are promoters of l o g K e a q - increase/decrease are also extracted from the SMILES notation for mechanistic interpretation. These QSPR models can also be applied to compute the reaction rate constants of organic contaminants.

Correlation intensity index: Building up models for mutagenicity of silver nanoparticles

Nanomaterials become significant component of economics. Consequently, nanomaterials become object of environmental sciences. There is a traditional list of endpoints which are indicators of the ecological risk. Mutagenicity is one of important component in this list. The quasi-SMILES approach, that in contrast to majority of work dedicated to modelling behaviour of nanomaterials gives possibility to consider experimental conditions as well as other circumstances which can impact the behaviour of nanomaterials is suggested. This is carried out via so-called quasi-SMILES. The quasi-SMILES is a line on of codes that contains all the above available eclectic data. Modelling process aimed to build up a model involves Correlation Intensity Index (CII) that is a new criterion of predictive potential of models. The scheme of calculation of CII is described in this work in the first time. The applying of CII together with Index of Ideality Correlation (IIC) in modelling of mutagenicity of silver nanoparticles by the Monte Carlo method using the CORAL software (http://www.insilico.eu/coral) indicates that application of the CII improves the predictive potential of these models for three random splits into the training set (75%) and validation set (25%).

QSAR modelling of larvicidal phytocompounds against Aedes aegypti using index of ideality of correlation

Aedes aegypti is the primary vector of several infectious viruses that cause yellow, dengue, chikungunya, and Zika fevers. Recently, plant-derived products have been tested as safe and eco-friendly larvicides against Ae. aegypti. The present study aimed to improve QSAR models for 62 larvicidal phytocompounds against Ae. aegypti via the Monte Carlo method based on the index of the ideality of correlation (IIC) criterion. The representation of structures was done with SMILES. Three splits were prepared randomly and three QSAR models were constructed using IIC target function. The molecular descriptors were selected from SMILES descriptors and the hydrogen-filled molecular graphs. The predictability of three models was evaluated on the validation sets, the r ² of which was 0.9770, 0.8660, and 0.8565 for models 1 to 3, respectively. The statistical results of three randomized splits indicated that robust, simple, predictive, and reliable models were obtained for different sets. From the modelling results, important descriptors were identified to enhance and reduce the larvicidal activity of compounds. Based on the identified important descriptors, some new structures of larvicidal compounds were proposed. The larvicidal activity of novel molecules designed further was supported by docking studies. Using the simple QSAR model, one can predict pLC₅₀ of new similarity larvicidal phytocompounds.

In silico enhancement of azo dye adsorption affinity for cellulose fibre through mechanistic interpretation under guidance of QSPR models using Monte Carlo method with index of ideality correlation

Azo dyes are a group of chemical moieties joined by azo (-N=N-) group with potential usefulness in different industrial applications. But these dyes are not devoid of hazardous consequence because of poor affinity for the fibre and discharge into the water stream. The chemical aspects of 72 azo dyes towards cellulose fibre in terms of their affinity by QSPR have been explored in the present work. We have employed two approaches, namely balance of correlation without IIC (TF₁) and balance of correlation with IIC (TF₂), to generate 16 QSAR models from 8 splits. The determination coefficient of calibration and validation set was found higher when the QSPR models were developed using the index of ideality correlation (IIC) parameter (TF₂). The model developed with TF₂ for split 3 was considered as a prominent model because the determination coefficient of the validation set was maximum (r ² = 0.9468). The applicability domain (AD) was also analysed based on 'statistical defect', d(A) for a SMILES attribute. The mechanistic interpretation was done by identifying the SMILES attributes responsible for the promoter of endpoint increase and promoter of endpoint decrease. These SMILES attributes were applied to design 15 new dyes with higher affinity for cellulose fibre.

QSPR/QSAR: State-of-Art, Weirdness, the Future

Ability of quantitative structure-property/activity relationships (QSPRs/QSARs) to serve for epistemological processes in natural sciences is discussed. Some weirdness of QSPR/QSAR state-of-art is listed. There are some contradictions in the research results in this area. Sometimes, these should be classified as paradoxes or weirdness. These points are often ignored. Here, these are listed and briefly commented. In addition, hypotheses on the future evolution of the QSPR/QSAR theory and practice are suggested. In particular, the possibility of extending of the QSPR/QSAR problematic by searching for the "statistical similarity" of different endpoints is suggested and illustrated by an example for relatively "distanced each from other" endpoints, namely (i) mutagenicity, (ii) anticancer activity, and (iii) blood-brain barrier.