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

The coefficient of conformism of a correlative prediction (CCCP): Building up reliable nano-QSPRs/QSARs for endpoints of nanoparticles in different experimental conditions encoded via quasi-SMILES

Simulation of the physicochemical and biochemical behavior of nanomaterials has its own specifics. However, the main goal of modeling for both traditional substances and nanomaterials is the same. This is an ecologic risk assessment. The universal indicator of toxicity is the n-octanol/water partition coefficient. Mutagenicity indicates the possibility of future undesirable environmental effects, possibly greater than toxicity. Models have been proposed for the octanol/water distribution coefficient of gold nanoparticles and the mutagenicity of silver nanoparticles. Unlike the previous studies, here the models are built using an updated scheme, which includes two improvements. Firstly, the computing involves a new criterion for prediction potential, the so-called coefficient of conformism of a correlative prediction (CCCP); secondly, the Las Vegas algorithm is used to select the potentially most promising models from a group of models obtained by the Monte Carlo algorithm. Apparently, CCCP is a measure of the predictive potential (not only correlation). This can give an advantage in developing a model in comparison to using the classic determination coefficient. Likely, CCCP can be more informative than the classical determination coefficient. The Las Vegas algorithm is able to improve the model obtained by the Monte Carlo method.

Does the accounting of the local symmetry fragments in quasi-SMILES improve the predictive potential of the QSAR models of toxicity toward tadpoles?

Models of toxicity to tadpoles have been developed as single parameters based on special descriptors which are sums of correlation weights, molecular features, and experimental conditions. This information is presented by quasi-SMILES. Fragments of local symmetry (FLS) are involved in the development of the model and the use of FLS correlation weights improves their predictive potential. In addition, the index of ideality correlation (IIC) and correlation intensity index (CII) are compared. These two potential predictive criteria were tested in models built through Monte Carlo optimization. The CII was more effective than IIC for the models considered here.

In-silico activity prediction and docking studies of some flavonol derivatives as anti-prostate cancer agents based on Monte Carlo optimization

The QSAR models are employed to predict the anti-proliferative activity of 81 derivatives of flavonol against prostate cancer using the Monte Carlo algorithm based on the index of ideality of correlation (IIC) criterion. CORAL software is employed to design the QSAR models. The molecular structures of flavonols are demonstrated using the simplified molecular input line entry system (SMILES) notation. The models are developed with the hybrid optimal descriptors i.e. using both SMILES and hydrogen-suppressed molecular graph (HSG). The QSAR model developed for split 3 is selected as a prominent model ([Formula: see text]= 0.727, [Formula: see text]= 0.628, [Formula: see text]= 0.642, and [Formula: see text]=0.615). The model is interpreted mechanistically by identifying the characteristics responsible for the promoter of the increase or decrease. The structural attributes as promoters of increase of pIC₅₀ were aliphatic carbon atom connected to double-bound (C…=…, aliphatic oxygen atom connected to aliphatic carbon (O…C…), branching on aromatic ring (c…(…), and aliphatic nitrogen (N…). The pIC₅₀ of eight natural flavonols with pIC₅₀ more than 4.0, were predicted by the best model. The molecular docking is also performed for natural flavonols on the PC-3 cell line using the protein (PDB: 3RUK).

The System of Self-Consistent Models: QSAR Analysis of Drug-Induced Liver Toxicity

Removing a drug-like substance that can cause drug-induced liver injury from the drug discovery process is a significant task for medicinal chemistry. In silico models can facilitate this process. Semi-correlation is an approach to building in silico models representing the prediction in the active (1)-inactive (0) format. The so-called system of self-consistent models has been suggested as an approach for two tasks: (i) building up a model and (ii) estimating its predictive potential. However, this approach has been tested so far for regression models. Here, the approach is applied to building up and estimating a categorical hepatotoxicity model using the CORAL software. This new process yields good results: sensitivity = 0.77, specificity = 0.75, accuracy = 0.76, and Matthew correlation coefficient = 0.51 (all compounds) and sensitivity = 0.83, specificity = 0.81, accuracy = 0.83 and Matthew correlation coefficient = 0.63 (validation set).

QSAR and molecular docking studies of isatin and indole derivatives as SARS 3CLpro inhibitors

The 3C-like protease (3CL^pro), known as the main protease of SARS-COV, plays a vital role in the viral replication cycle and is a critical target for the development of SARS inhibitor. Comparative sequence analysis has shown that the 3CL^pro of two coronaviruses, SARS-CoV-2 and SARS-CoV, show high structural similarity, and several common features are shared among the substrates of 3CL^pro in different coronaviruses. The goal of this study is the development of validated QSAR models by CORAL software and Monte Carlo optimization to predict the inhibitory activity of 81 isatin and indole-based compounds against SARS CoV 3CL^pro. The models were built using a newer objective function optimization of this software, known as the index of ideality correlation (IIC), which provides favorable results. The entire set of molecules was randomly divided into four sets including: active training, passive training, calibration and validation sets. The optimal descriptors were selected from the hybrid model by combining SMILES and hydrogen suppressed graph (HSG) based on the objective function. According to the model interpretation results, eight synthesized compounds were extracted and introduced from the ChEMBL database as good SARS CoV 3CL^pro inhibitor. Also, the activity of the introduced molecules further was supported by docking studies using 3CL^pro of both SARS-COV-1 and SARS-COV-2. Based on the results of ADMET and OPE study, compounds CHEMBL4458417 and CHEMBL4565907 both containing an indole scaffold with the positive values of drug-likeness and the highest drug-score can be introduced as selected leads.

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.

In silico prediction of mosquito repellents for clothing application

Use of protective clothing is a simple and efficient way to reduce the contacts with mosquitoes and consequently the probability of transmission of diseases spread by them. This mechanical barrier can be enhanced by the application of repellents. Unfortunately the number of available repellents is limited. As a result, there is a crucial need to find new active and safer molecules repelling mosquitoes. In this context, a structure-activity relationship (SAR) model was proposed for the design of repellents active on clothing. It was computed from a dataset of 2027 chemicals for which repellent activity on clothing was measured against Aedes aegypti. Molecules were described by means of 20 molecular descriptors encoding physicochemical properties, topological information and structural features. A three-layer perceptron was used as statistical tool. An accuracy of 87% was obtained for both the training and test sets. Most of the wrong predictions can be explained. Avenues for increasing the performances of the model have been proposed.

The searching for agents for Alzheimer's disease treatment via the system of self-consistent models

Robust quantitative structure-activity relationships (QSARs) for hBACE-1 inhibitors (pIC50) for a large database (n = 1706) are established. New statistical criteria of the predictive potential of models are suggested and tested. These criteria are the index of ideality of correlation (IIC) and the correlation intensity index (CII). The system of self-consistent models is a new approach to validate the predictive potential of QSAR-models. The statistical quality of models obtained using the CORAL software (http://www.insilico.eu/coral) for the validation sets is characterized by the average determination coefficient R²_v= 0.923, and RMSE =0.345. Three new promising molecular structures which can become inhibitors hBACE-1 are suggested.

CORAL: Monte Carlo based global QSAR modelling of Bruton tyrosine kinase inhibitors using hybrid descriptors

Global QSAR modelling was performed to predict the pIC₅₀ values of 233 diverse heterocyclic compounds as BTK inhibitors with the Monte Carlo algorithm of CORAL software using the DCW hybrid descriptors extracted from SMILES notations of molecules. The dataset of 233 BTK inhibitors was randomly split into training, invisible training, calibration and validation sets. The index of ideality of correlation was also applied to build and judge the predictability of the QSAR models. Eight global QSAR models based on the hybrid optimal descriptor using two target functions, i.e. TF₁ (W_IIC = 0) and TF₂ (W_IIC = 0.2) have been constructed. The statistical parameters of QSAR models computed by TF₂ are more reliable and robust and were used to predict the pIC₅₀ values. The model constructed for split 4 via TF₂ is regarded as the best model and the numerical values of r²_Train, r²_Valid, Q²_Train and Q²_Valid are equal to 0.7981, 0.7429, 0.7898 and 0.6784, respectively. By internal and external validation techniques, the predictability and reliability of the designed models have been assessed. The structural attributes responsible for the increase and decrease of pIC₅₀ of BTK inhibitors were also identified.

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.

DFT based QSAR study on quinolone-triazole derivatives as antibacterial agents

QSAR modeling was performed on 39 quinolone-triazole derivatives against gram-positive Staphylococcus aureus and gram-negative Pseudomonas aeruginosa bacteria. The molecular structures were optimized using the DFT/B3LYP method and 6-31 G basis set. Molecular descriptors were extracted using quantum mechanical calculations. The hierarchical cluster analysis was performed for a rational subset division. The initial dataset was divided into calibration and validation sets, and modeling was done by stepwise MLR method for each of the two bacteria. Internal and external validation methods confirmed the robustness and predictability of the obtained models. According to the obtained model for S. aureus (R² = 0.889, R²_ext = 0.938, Q²_LOO = 0.853), the four descriptors- partial atomic charges for the N1 atom in triazole and C7 of the quinolone nucleus, 4-carbonyl bond length, and ¹³C-NMR chemical shift of 3-carboxylic acid- were found to be the descriptors controlling the activity. According to the obtained model for P. aeruginosa (R² = 0.957, R²_ext = 0.923, Q²_LOO = 0.909), the O atom's partial charge in carbonyl, LUMO-HOMO energy gap, and logP were found to be the descriptors having the highest correlation with the antibacterial activity. Finally, some new compounds with higher activities were designed and proposed.

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.

Monte-Carlo method-based QSAR model to discover phytochemical urease inhibitors using SMILES and GRAPH descriptors

Urease inhibitors are known to play a vital role in the field of medicine as well as agriculture. Special attention is attributed to the development of novel urease inhibitors with a view to treat the Helicobacter pylori infection. Amongst a number of urease inhibitors, a large number of molecules fail in vivo and in clinical trials due to their hydrolytic instability and toxicity profile. The search for potential inhibitors may require screening of large and diverse databases of small molecules and to design novel molecules. We developed a Monte-Carlo method-based QSAR model to predict urease inhibiting potency of molecules using SMILES and GRAPH descriptors on an existing diverse database of urease inhibitors. The QSAR model satisfies all the statistical parameters required for acceptance as a good model. The model is applied to identify urease inhibitors among the wide range of compounds in the phytochemical database, NPACT, as a test case. We combine the ligand-based and structure-based drug discovery methods to improve the accuracy of the prediction. The method predicts pIC₅₀ and estimates docking score of compounds in the database. The method may be applied to any other database or compounds designed in silico to discover novel drugs targeting urease.Communicated by Ramaswamy H. Sarma.

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.