Monte Carlo method based QSAR modeling of maleimide derivatives as glycogen synthase kinase-3β inhibitors

In silico aquatic toxicity prediction of chemicals toward Daphnia magna and fathead minnow using Monte Carlo approaches

The fast-increasing use of chemicals led to large numbers of chemical compounds entering the aquatic environment, raising concerns about their potential effects on ecosystems. Therefore, assessment of the ecotoxicological features of organic compounds on aquatic organisms is very important. Daphnia magna and Fathead minnow are two aquatic species that are commonly tested as standard test organisms for aquatic risk assessment and are typically chosen as the biological model for the ecotoxicology investigations of chemical pollutants. Herein, global quantitative structure-toxicity relationship (QSTR) models have been developed to predict the toxicity (pEC(LC)50) of a large dataset comprising 2106 chemicals toward Daphnia magna and Fathead minnow. The optimal descriptor of correlation weights (DCWs) is calculated using the notation of simplified molecular input line entry system (SMILES) and is used to construct QSTR models. Three target functions, TF1, TF2, and TF3 are utilized to generate 12 QSTR models from four splits, and their statistical characteristics are also compared. The designed QSTR models are validated using both internal and external validation criteria and are found to be reliable, robust, and excellently predictive. Among the models, those generated using the TF3 demonstrate the best statistical quality with R2 values ranging from 0.9467 to 0.9607, Q2 values ranging from 0.9462 to 0.9603 and RMSE values ranging from 0.3764 to 0.4413 for the validation set. The applicability domain and the mechanistic interpretations of generated models were also discussed.

Correlation intensity index-index of ideality of correlation: A hyphenated target function for furtherance of MAO-B inhibitory activity assessment

Monoamine oxidases are the enzymes involved in the management of brain homeostasis through oxidative deamination of monoamines such as neurotransmitters, tyramine etc. The excessive production of monoamine oxidase-B specifically results in numerous neurodegenerative disorders like Alzheimer's and Parkinson's diseases. Inhibitors of monoamine oxidase-B are applied in the management of these disorders. Here in this article we have developed robust hybrid descriptor based QSAR models related to 123 monoamine oxidase-B inhibitors through CORAL software by means of Monte Carlo optimization method. Three target functions were applied to prepare QSAR models and three splits were made for each target function. The most reliable, robust and better predictive QSAR models were developed with TF3 (correlation intensity index -index of ideality of correlation). Correlation intensity index showed positive effect on QSAR models. The structural features obtained from the QSAR modeling were incorporated in newly designed molecules and exhibited positive effect on their endpoint. Significant binding interactions were represented by these molecules in docking studies. Molecule B5 displayed prominent pIC50 (8.3) and binding affinity (-11.5 kcal mol-1) towards monoamine oxidase-B.

Molecular toxicity of nitrobenzene derivatives to tetrahymena pyriformis based on SMILES descriptors using Monte Carlo, docking, and MD simulations

It is challenging to model the toxicity of nitroaromatic compounds due to limited experimental data. Nitrobenzene derivatives are commonly used in industry and can lead to environmental contamination. Extensive research, including several QSPR studies, has been conducted to understand their toxicity. Predictive QSPR models can help improve chemical safety, but their limitations must be considered, and the molecular factors affecting toxicity should be carefully investigated. The latest QSPR methods, molecular modeling techniques, machine learning algorithms, and computational chemistry tools are essential for developing accurate and robust models. In this work, we used these methods to study a series of fifty compounds derived from nitrobenzene. The Monte Carlo approach was used for QSPR modeling by applying the SMILES molecular structure representation and optimal molecular descriptors. The correlation ideality index (CII) and correlation contradiction index (CCI) were further introduced as validation parameters to estimate the developed models' predictive ability. The statistical quality of the CII models was better than those without CII. The best QSPR model with the following statistical parameters (Split-3): (R2 = 0.968, CCC = 0.984, IIC = 0.861, CII = 0.979, Q2 = 0.954, QF12 = 0.946, QF22 = 0.938, QF32 = 0.947, Rm2 = 0.878, RMSE = 0.187, MAE = 0.151, FTraining = 390, FInvisible = 218, FCalibration = 240, RTest2 = 0.905) was selected to generate the studied promoters with increasing and decreasing activity.

In silico approaches to develop new phenyl-pyrimidines as glycogen synthase kinase 3 (GSK-3) inhibitors with halogen-bonding capabilities: 3D-QSAR CoMFA/CoMSIA, molecular docking and molecular dynamics studies

Glycogen synthase kinase 3 (GSK-3) is involved in different diseases, such as manic-depressive illness, Alzheimer's disease and cancer. Studies have shown that insulin inhibits GSK-3 to keep glycogen synthase active. Inhibiting GSK-3 may have an indirect pro-insulin effect by favouring glycogen synthesis. Therefore, the development of GSK-3 inhibitors can be a useful alternative for the treatment of type II diabetes. Aminopyrimidine derivatives already proved to be interesting GSK-3 inhibitors. In the current study, comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) have been performed on a series of 122 aminopyrimidine derivatives in order to generate a robust model for the rational design of new compounds with promising antidiabetic activity. The q2 values obtained for the best CoMFA and CoMSIA models have been 0.563 and 0.598, respectively. In addition, the r2 values have been 0.823 and 0.925 for CoMFA and CoMSIA, respectively. The models were statistically validated, and from the contour maps analysis, a proposal of 10 new compounds has been generated, with predicted pIC50 higher than 9. The final contribution of our work is that: (a) we provide an extensive structure-activity relationship for GSK-3 inhibitory pyrimidines; and (b) these models may speed up the discovery of GSK-3 inhibitors based on the aminopyrimidine scaffold. Finally, we carried out docking and molecular dynamics studies of the two best candidates, which were shown to establish halogen-bond interactions with the enzyme.Communicated by Ramaswamy H. Sarma.

Comparison of various methods for validity evaluation of QSAR models

BACKGROUND

Quantitative structure-activity relationship (QSAR) modeling is one of the most important computational tools employed in drug discovery and development. The external validation of QSAR models is the main point to check the reliability of developed models for the prediction activity of not yet synthesized compounds. It was performed by different criteria in the literature.

METHODS

In this study, 44 reported QSAR models for biologically active compounds reported in scientific papers were collected. Various statistical parameters of external validation of a QSAR model were calculated, and the results were discussed.

RESULTS

The findings revealed that employing the coefficient of determination (r²) alone could not indicate the validity of a QSAR model. The established criteria for external validation have some advantages and disadvantages which should be considered in QSAR studies.

CONCLUSION

This study showed that these methods alone are not only enough to indicate the validity/invalidity of a QSAR model.

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.

Virtual Screening of Anti-Cancer Compounds: Application of Monte Carlo Technique

Possibility and necessity of standardization of predictive models for anti-cancer activity are discussed. The hypothesis about rationality of common quantitative analysis of anti-cancer activity and carcinogenicity is developed. Potential of optimal descriptors to be used as a tool to build up predictive models for anti-cancer activity is examined from practical point of view. Various perspectives of application of optimal descriptors are reviewed. Stochastic nature of phenomena which are related to carcinogenic potential of various substances can be successfully detected and interpreted by the Monte Carlo technique. Hypothesises related to practical strategy and tactics of the searching for new anticancer agents are suggested.

Use of the Monte Carlo Method for OECD Principles-Guided QSAR Modeling of SIRT1 Inhibitors

SIRT1 inhibitors offer therapeutic potential for the treatment of a number of diseases including cancer and human immunodeficiency virus infection. A diverse series of 45 compounds with reported SIRT1 inhibitory activity has been employed for the development of quantitative structure-activity relationship (QSAR) models using the Monte Carlo optimization method. This method makes use of simplified molecular input line entry system notation of the molecular structure. The QSAR models were built up according to OECD principles. Three subsets of three splits were examined and validated by respective external sets. All the three described models have good statistical quality. The best model has the following statistical characteristics: R²  = 0.8350, Q²_test  = 0.7491 for the test set and R²  = 0.9655, Q²_ext  = 0.9261 for the validation set. In the mechanistic interpretation, structural attributes responsible for the endpoint increase and decrease are defined. Further, the design of some prospective SIRT1 inhibitors is also presented on the basis of these structural attributes.

Monte Carlo-based quantitative structure-activity relationship models for toxicity of organic chemicals to Daphnia magna

Quantitative structure-activity relationships (QSARs) for toxicity of a large set of 758 organic compounds to Daphnia magna were built up. The simplified molecular input-line entry system (SMILES) was used to represent the molecular structure. The Correlation and Logic (CORAL) software was utilized as a tool to develop the QSAR models. These models are built up using the Monte Carlo method and according to the principle "QSAR is a random event" if one checks a group of random distributions in the visible training set and the invisible validation set. Three distributions of the data into the visible training, calibration, and invisible validation sets are examined. The predictive potentials (i.e., statistical characteristics for the invisible validation set of the best model) are as follows: n = 87, r²  = 0.8377, root mean square error = 0.564. The mechanistic interpretations and the domain of applicability of built models are suggested and discussed. Environ Toxicol Chem 2016;35:2691-2697. © 2016 SETAC.