CORAL: probing the structural requirements for α-amylase inhibition activity of 5-(3-arylallylidene)-2-(arylimino)thiazolidin-4-one derivatives based on QSAR with correlation intensity index, molecular docking, molecular dynamics, and ADMET studies

The present study aims to examine the structural requirements governing α-amylase inhibitory activity of 5-(3-arylallylidene)-2-(arylimino)thiazolidin-4-one derivatives and their precursors by employing a multifaceted approach combining in vitro and in silico studies. The in vitro assay findings revealed strong inhibitory effect of this class of compounds against α-amylase and compound 20 exhibited maximum percentage inhibition of 88.54 ± 0.69, 84.98 ± 0.40, 77.26 ± 0.75, 67.80 ± 0.54, and 62.93 ± 1.17 at 200, 100, 50, 25, and 12.5 µg mL-1, respectively. Multiple CORAL QSAR models were developed from the randomly distributed eight splits by employing two target functions (TF1, TF2 with WCII = 0.0 and = 0.3, respectively), and the quality of predictions by the produced models was validated with the help of various statistical parameters. The model M-4 (R2Val = 0.8799) and model M-11 (R2Val = 0.9064) were the leading models developed by using TF1 and TF2. We designed five new congeneric inhibitors (D-1 to D-5) by incorporating SMILES features positively correlating with the activity. Molecular docking experiments were carried out to confirm the binding of these new inhibitors with the biological receptor α-amylase (PDB ID: 7TAA). Furthermore, molecular dynamic simulations provided a thorough knowledge of the binding process by shedding insight into the dynamic behavior and stability of the ligand-receptor complex over time. The results of this study highlight the key structural characteristics needed for improved α-amylase inhibitory efficacy and provide a rational basis for the development of more effective inhibitors.Communicated by Ramaswamy H. Sarma.

Identifications of good and bad structural fragments of hydrazone/2,5-disubstituted-1,3,4-oxadiazole hybrids with correlation intensity index and consensus modelling using Monte Carlo based QSAR studies, their molecular docking and ADME analysis

The application of QSAR along with other in silico tools like molecular docking, and molecular dynamics provide a lot of promise for finding new treatments for life-threatening diseases like Type 2 diabetes mellitus (T2DM). The present study is an attempt to develop Monte Carlo algorithm-based QSAR models using freely available CORAL software. The experimental data on the α-amylase inhibition by a series of benzothiazole-linked hydrazone/2,5-disubstituted-1,3,4-oxadiazole hybrids were selected as endpoint for the model generation. Initially, a total of eight QSAR models were built using correlation intensity index (CII) as a criterion of predictive potential. The model developed from split 6 using CII was the most reliable because of the highest numerical value of the determination coefficient of the validation set (r²_VAL = 0.8739). The important structural fragments responsible for altering the endpoint were also extracted from the best-built model. With the goal of improved prediction quality and lower prediction errors, the validated models were used to build consensus models. Molecular docking was used to know the binding mode and pose of the selected derivatives. Further, to get insight into their metabolism by living beings, ADME studies were investigated using internet freeware, SwissADME.

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.

The VEGAHUB Platform: The Philosophy and the Tools

VEGAHUB (www.vegahub.eu) is a repository of freely available, downloadable tools based on computational toxicology methodologies. The main software tool available in VEGAHUB is VEGA QSAR software encoding more than 90 quantitative structure-activity relationship (QSAR) models for tens of endpoints for human toxicology, ecotoxicology, environmental, physico-chemical and toxicokinetic properties. However, beyond VEGA QSAR, VEGAHUB offers several other tools. Here, we present these resources, the possibilities to fully exploit them and the ways in which to integrate results provided by different VEGAHUB tools. Read-across and weight-of-evidence represent a major advantage of VEGAHUB. Integration between hazard and exposure is provided within innovative tools, which are specific for well-defined scenarios, such as those for cosmetic products. Prioritisation can be achieved by integrating results from 48 models. Finally, we highlight how some tools may not only fit predefined endpoints but also could be applied to general problems and research applications in the QSAR field. A couple of examples are provided, in which a critical assessment of the predictions and the documentation associated with the prediction are considered, in order to properly assess the quality of the results. These results may be associated with different levels of uncertainty or even be conflicting.

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.

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.

Good and bad molecular fingerprints for human rhinovirus 3C protease inhibition: identification, validation, and application in designing of new inhibitors through Monte Carlo-based QSAR study

HRV 3 C protease (HRV 3Cpro) is an important target for common cold and upper respiratory tract infection. Keeping in view of the non-availability of drug for the treatment, newer computer-based modelling strategies should be applied to rationalize the process of antiviral drug discovery in order to decrease the valuable time and huge expenditure of the process. The present work demonstrates a structure wise optimization using Monte Carlo-based QSAR method that decomposes ligand compounds (in SMILES format) into several molecular fingerprints/descriptors. The current state-of-the-art in QSAR study involves the balance of correlation approach using four different sets: training, invisible training, calibration, and validation. The final models were also validated through mean absolute error, index of ideality of correlation, Y-randomization and applicability domain analysis. R2 and Q2 values for the best model were 0.8602, 0.8507 (training); 0.8435, 0.8331 (invisible training); 0.7424, 0.7020 (calibration); 0.5993, 0.5216 (validation), respectively. The process identified some molecular substructures as good and bad fingerprints depending on their effect to increase or decrease the HRV 3Cpro inhibition. Finally, new inhibitors were designed based on the fundamental concept to replace the bad fragments with the good fragments as well as including more good fragments into the structure. The study points out the importance of the fingerprint based drug design strategy through Monte Carlo optimization method in the modelling of HRV 3Cpro inhibitors.

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.

Design and development of novel focal adhesion kinase (FAK) inhibitors using Monte Carlo method with index of ideality of correlation to validate QSAR

Quantitative structure-activity relationship (QSAR) modelling of 55 focal adhesion kinase (FAK) (EC 2.7.10.2) inhibitors of triazinic nature was performed using the Monte Carlo method. The QSAR models were designed by CORAL software, and optimal descriptors were calculated with the simplified molecular input line entry system (SMILES). Four splits were made from the triazinic derivative data by random division into training, invisible training, calibration and validation sets. The QSAR results from these four random splits were robust, very simple, predictive and reliable. The best statistical parameters of the validation set (r² = 0.8398 and Q² = 0.7722) for the QSAR equation for split 3 with IIC = 0.9127 were obtained. The predictive potential of QSAR models of FAK inhibitors was explored by applying the index of ideality of correlation (IIC), which is a new criterion for the prediction of the potential for quantitative structure-property activity relationships (QSPRs/QSARs). The present method follows OECD principles.

Prediction of the adsorption coefficients of some aromatic compounds on multi-wall carbon nanotubes by the Monte Carlo method

In this investigation, quantitative structure-property relationship (QSPR) modelling of adsorption coefficients of 69 aromatic compounds on multi-wall carbon nanotubes (MWCNTs) was studied using the Monte Carlo method. QSPR models were calculated with CORAL software, and optimal descriptors were calculated with the simplified molecular input line entry system (SMILES) and hydrogen-suppressed molecular graphs (HSGs). The aromatic compound data set was randomly split into training, invisible training, calibration and validation sets. Analysis of three probes of the Monte Carlo optimization with three random splits was done. The results from three random splits displayed robust, very simple, predictable and reliable models for the training, invisible training, calibration and validation sets with a coefficient of determination (r²) equal to 0.9463-0.8528, 0.9020-0.8324, 0.9606-0.9178 and 0.9573-0.8228, respectively. As a result, the models obtained help to identify the hybrid descriptors for the increase and the decrease of the adsorption coefficient of aromatic compounds on MWCNTs. This simple QSPR model can be used for the prediction of the adsorption coefficient of numerous aromatic compounds on MWCNTs.

Towards the Development of Global Nano-Quantitative Structure-Property Relationship Models: Zeta Potentials of Metal Oxide Nanoparticles

Zeta potential indirectly reflects a charge of the surface of nanoparticles in solutions and could be used to represent the stability of the colloidal solution. As processes of synthesis, testing and evaluation of new nanomaterials are expensive and time-consuming, so it would be helpful to estimate an approximate range of properties for untested nanomaterials using computational modeling. We collected the largest dataset of zeta potential measurements of bare metal oxide nanoparticles in water (87 data points). The dataset was used to develop quantitative structure–property relationship (QSPR) models. Essential features of nanoparticles were represented using a modified simplified molecular input line entry system (SMILES). SMILES strings reflected the size-dependent behavior of zeta potentials, as the considered quasi-SMILES modification included information about both chemical composition and the size of the nanoparticles. Three mathematical models were generated using the Monte Carlo method, and their statistical quality was evaluated (R² for the training set varied from 0.71 to 0.87; for the validation set, from 0.67 to 0.82; root mean square errors for both training and validation sets ranged from 11.3 to 17.2 mV). The developed models were analyzed and linked to aggregation effects in aqueous solutions.

Comprehensive data on a 2D-QSAR model for Heme Oxygenase isoform 1 inhibitors

The data have been obtained from the Heme Oxygenase Database (HemeOxDB) and refined according to the 2D-QSAR requirements. These data provide information about a set of more than 380 Heme Oxygenase-1 (HO-1) inhibitors. The development of the 2D-QSAR model has been undertaken with the use of CORAL software using SMILES, molecular graphs and hybrid descriptors (SMILES and graph together). The 2D-QSAR model regressions for HO-1 half maximal inhibitory concentration (IC₅₀) expressed as pIC₅₀ (pIC₅₀=−LogIC₅₀) are here included. The 2D-QSAR model was also employed to predict the HO-1 pIC₅₀values of the FDA approved drugs that are herewith reported.

Nano-QSAR Model for Predicting Cell Viability of Human Embryonic Kidney Cells

Traditional Quantitative Structure-Activity Relationships (QSAR) models based on molecular descriptors as translators of chemical information show some drawbacks in predicting toxicity of nanomaterials due to their unique properties and to their nonhomogeneous structure.This chapter provides instructions on how to use CORAL, freely available software for building nano-QSAR models. CORAL makes use of descriptors based on "quasi-SMILES" representing physicochemical features and/or experimental conditions as an alternative to traditional SMILES encoding chemical structure to build up predictive nano-QSAR models for cytotoxicity.

New Quantitative Structure-Activity Relationship Models Improve Predictability of Ames Mutagenicity for Aromatic Azo Compounds

Existing Quantitative Structure-Activity Relationship (QSAR) models have limited predictive capabilities for aromatic azo compounds. In this study, 2 new models were built to predict Ames mutagenicity of this class of compounds. The first one made use of descriptors based on simplified molecular input-line entry system (SMILES), calculated with the CORAL software. The second model was based on the k-nearest neighbors algorithm. The statistical quality of the predictions from single models was satisfactory. The performance further improved when the predictions from these models were combined. The prediction results from other QSAR models for mutagenicity were also evaluated. Most of the existing models were found to be good at finding toxic compounds but resulted in many false positive predictions. The 2 new models specific for this class of compounds avoid this problem thanks to a larger set of related compounds as training set and improved algorithms.

QSPR modelling of dielectric constants of π-conjugated organic compounds by means of the CORAL software

Simplified molecular input line entry system (SMILES) notations of 116 π-conjugated organic compounds have been used in three random splits to develop single optimal descriptor based quantitative structure-property relationships (QSPR) models for the prediction of dielectric constants by CORAL (CORrelation And Logic). Four kinds of optimal descriptors were obtained based on SMILES, hydrogen suppressed graph (HSG), graph of atomic orbitals (GAO) and hybrid descriptors. The Monte Carlo optimization was carried out for each random split by three different methods: (i) classic scheme; (ii) balance of correlations; and (iii) balance of correlations with ideal slopes. The QSPR models gave reliable and accurate values of dielectric constant for all the π-conjugated organic compounds. SMILES and the hybrid-based QSPR model provided the best accuracy for the prediction of dielectric constants. Statistical characteristics of the QSPR model-1 based on classic scheme method are n = 110, r(2) = 0.860, Q(2) = 0.860, s = 1.84, MAE = 1.30 and F = 696 (training set), n = 6, r(2) = 0.947, Q(2) = 0.876, s = 0.955, MAE = 0.647 and F = 71 (test set). These QSPR models are further validated by an external validation set of 25 molecules and the robustness is checked by parameters like k, kk, rm(2), r(*)m(2), average rm(2), ∆rm(2)and randomization technique ([Formula: see text]).