Construction of pioneering quantitative structure activity relationship screening models for abuse potential of designer drugs using index of ideality of correlation in monte carlo optimization

Structural Features Influencing the Bioactive Conformation of Angiotensin II and Angiotensin A: Relationship between Receptor Desensitization, Addiction, and the Blood-Brain Barrier

The N-terminal portion of the octapeptide angiotensin II (DRVYIHPF; AngII), a vasopressor peptide that favorably binds to, and activates, AngII type 1 receptor (AT1R), has an important role in maintaining bioactive conformation. It involves all three charged groups, namely (i) the N-terminal amino group cation, (ii) the Asp sidechain anion and (iii) the Arg guanidino cation. Neutralization of any one of these three charged groups results in a substantial reduction (<5%) in bioactivity, implicating a specialized function for this cluster. In contrast, angiotensin A (ARVYIHPF; AngA) has reduced bioactivity at AT1R; however, replacement of Asp in AngII with sarcosine (N-methyl-glycine) not only restores bioactivity but increases the activity of agonist, antagonist, and inverse agonist analogues. A bend produced at the N-terminus by the introduction of the secondary amino acid sarcosine is thought to realign the functional groups that chaperone the C-terminal portion of AngII, allowing transfer of the negative charge originating at the C-terminus to be transferred to the Tyr hydroxyl-forming tyrosinate anion, which is required to activate the receptor and desensitizes the receptor (tachyphylaxis). Peptide (sarilesin) and nonpeptide (sartans) moieties, which are long-acting inverse agonists, appear to desensitize the receptor by a mechanism analogous to tachyphylaxis. Sartans/bisartans were found to bind to alpha adrenergic receptors resulting in structure-dependent desensitization or resensitization. These considerations have provided information on the mechanisms of receptor desensitization/tolerance and insights into possible avenues for treating addiction. In this regard sartans, which appear to cross the blood-brain barrier more readily than bisartans, are the preferred drug candidates.

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.

Quasi-SMILES as a basis to build up models of endpoints for nanomaterials

Simplified molecular input-line entry system (SMILES) is a format for representing of the molecular structure. Quasi-SMILES is an extended format for representing molecular structure data and some eclectic data, which in principle could be applied to improve a model's predictive potential. Nano-quantitative structure-property relationships (nano-QSPRs) for energy gap (Eg, eV) of the metals oxide nanoparticles based on the quasi-SMILES give a predictive model for Eg, characterized by the following statistical quality for external validation set n = 22, R2 = 0.83, RMSE = 0.267.

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.

Paradox of 'ideal correlations': improved model for air half-life of persistent organic pollutants

The persistence of organic pollutants is an important environmental property due to the extended possibility to have an impact of corresponding substances. In many cases, the experimental values of the thousands of contaminants are missing. The object of the study is novel computational modelling for air pollutions. Quantitative structure-property relationship (QSPR) for air half-life has been built using the Monte Carlo method with applying the index of ideality of correlation (IIC). The basis of the predictive model of air half-life is the representation of the molecular structure by simplifying molecular input-line entry system (SMILES) and numerical data on the above endpoint (expressed by hours) converted to a decimal logarithm. The statistical quality of the model has been checked up with different validation metrics and is quite good. Paradoxically, the improvement of the statistical quality via the IIC for the validation set is done in detriment to the training set. The new model has performed better than those obtained previously on the same set of compounds, for the prediction of new compounds in the validation set. Some semi-quantitative indicators for the mechanistic interpretation of the model are suggested.

Monte Carlo based QSGFEAR: prediction of Gibb's free energy of activation at different temperatures using SMILES based descriptors

Monte Carlo optimization based QSGFEAR model development using CII results in the formation of more reliable, robust and predictive models.

Unswerving modeling of hepatotoxicity of cadmium containing quantum dots using amalgamation of quasiSMILES, index of ideality of correlation, and consensus modeling

Liver toxicity of quantum dots varies with size, concentration, and other structural as well as experimental parameters. For modeling hepatotoxicity, the eclectic data associated with cadmium containing quantum dots have been used in the creation of quasiSMILES for their representation. The core diameter is normalized for wider applicability and the index of the ideality of correlation is applied to construct the better quantitative features toxicity relationship models. Total eight splits are created and the best model is obtained through split 1 with better prediction criteria of validation set objects. The values of all statistical criteria used in the quality determination of a QSAR model are within the specified range for all the eight toxicity models developed here. Factors like TGA ligand and 0.6-0.7 nm diameter are favorable for liver toxicity while L-cysteine ligand and 0.5-0.6 nm core diameter are helpful in the reduction of toxicity. Further, the intelligent consensus modeling process forms a total of 40 individual and 20 consensus models and the best individual and consensus models are 'Good' in MAE-based criteria. The consensus modeling enhances the prediction ability as well as the accuracy of the developed models and increases the applicability space of the built models for hepatotoxicity prediction of quantum dots.

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.

Exploring biological efficacy of novel benzothiazole linked 2,5-disubstituted-1,3,4-oxadiazole hybrids as efficient α-amylase inhibitors: Synthesis, characterization, inhibition, molecular docking, molecular dynamics and Monte Carlo based QSAR studies

In an effort to explore a class of novel antidiabetic agents, we have made an effort to synergize the α-amylase inhibitory potential of 1,3-benzothiazole and 1,3,4-oxadiazole scaffolds by combining the two into a single structure via an ether linkage. The structure of synthesized benzothiazole clubbed oxadiazole derivatives are established by different spectral techniques. The synthesized hybrids are evaluated for their in vitro inhibitory potential against α-amylase. Compound 8f is found to be the most potent with a significant inhibition (87.5 ± 0.74% at 50 μg/mL, 82.27 ± 1.85% at 25 μg/mL and 79.94 ± 1.88% at 12.5 μg/mL) when compared to positive control acarbose (77.96 ± 2.06%, 71.17 ± 0.60%, 67.24 ± 1.16% at 50 μg/mL, 25 μg/mL and 12.5 μg/mL concentration). Molecular docking of the most potent enzyme inhibitor, 8f, shows promising interaction with the binding site of biological macromolecule Aspergillus oryzae α-amylase (PDB ID: 7TAA) and human pancreatic α-amylase (PDB ID: 3BAJ). To a step further, in-depth QSAR studies show a significant correlation between the experimental and the predicted inhibitory activities with the best R_validation²= 0.8701. The developed QSAR model can provide ample information about the structural features responsible for the increase and decrease of inhibitory activity. The mechanistic interpretation of the structure-activity relationship (SAR) is done with the help of combined computational calculations i.e. molecular docking and QSAR. Finally, molecular dynamic simulations are performed to get an insight into the binding mode of the most potent derivative with α-amylase from A. oryzae (PDB ID: 7TAA) and human pancreas (PDB ID: 3BAJ).

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.

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.

Cytotoxicity of quantum dots: Use of quasiSMILES in development of reliable models with index of ideality of correlation and the consensus modelling

The assessment of cytotoxicity of quantum dots is very essential for environmental and health risk analysis. In the present work we have modelled HeLa cell cytotoxicity of sixty one CdSe quantum dots with ZnS shell as a function of its experimental conditions and molecular construction using quasiSMILES representations. The index of ideality of correlation helps in the building of ten statistically significant models having good fitting ability with value of R² ranging from 0.8414 to 0.9609 for the training set. The split 5 model is rated as the best model with values of R², Q²_F1, Q²_F2 and Q²_F3 as 0.8964, 0.8267, 0.8264 and 0.8777 respectively for the calibration set. The extraction of features causing increase and decrease of cytotoxicity of quantum dots indicates importance of neutral surface charge, surface modified with protein, 72 h exposure time, combination of MTT assay with surface protein in decreasing the cytotoxicity. Amphiphilic polymer, polyol ligand with neutral charge, 0.5 - 0.6 nm quantum dot diameter with lipid ligand and unmodified positively charged surface are grouped in toxicity enhancer features. Further, consensus modelling using split 5 and 8 patterns enhances the prediction quality by increasing the R²_val to 0.9361 and 0.9656 respectively.

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.