CORAL: Quantitative Structure Retention Relationship (QSRR) of flavors and fragrances compounds studied on the stationary phase methyl silicone OV-101 column in gas chromatography using correlation intensity index and consensus modelling

Synthesis of thiazolidine-2,4-dione tethered 1,2,3-triazoles as α-amylase inhibitors: In vitro approach coupled with QSAR, molecular docking, molecular dynamics and ADMET studies

A new series of thiazolidine-2,4-dione tethered 1,2,3-triazole derivatives were designed, synthesized and screened for their α-amylase inhibitory potential employing in vitro and in silico approaches. The target compounds were synthesized with the help of Cu (I) catalyzed [3 + 2] cycloaddition of terminal alkyne with numerous azides, followed by unambiguously characterizing the structure by employing various spectroscopic approaches. The synthesized derivatives were assessed for their in vitro α-amylase inhibition and it was found that thiazolidine-2,4-dione derivatives 6e, 6j, 6o, 6u and 6x exhibited comparable inhibition with the standard drug acarbose. The compound 6e with a 7-chloroquinolinyl substituent on the triazole ring exhibited significant inhibition potential with IC50 value of 0.040 μmol mL-1 whereas compound 6c (IC50 = 0.099 μmol mL-1) and 6h (IC50 = 0.098 μmol mL-1) were poor inhibitors. QSAR studies revealed the positively correlating descriptors that aid in the design of novel compounds. Molecular docking was performed to investigate the binding interactions with the active site of the biological receptor and the stability of the complex over a period of 100 ns was examined using molecular dynamics studies. The physiochemical properties and drug-likeliness behavior of the potent derivatives were investigated by carrying out the ADMET studies.

A general procedure for finding potentially erroneous entries in the database of retention indices

BACKGROUND

The NIST retention index database is one the most widely used sources of retention indices. In both untargeted analysis and machine learning studies filtering for potential errors is rather lacking or nonexistent. According to our estimates about 80% of the compounds from both NIST 17 and NIST 20 retention index databases have only one RI value per stationary phase, which makes searching for erroneous values with statistical methods impossible. Manual inspection is also impractical because the database contains more than 300 000 entries.

RESULTS

We suggest a two-step procedure to find potentially erroneous retention indices based on machine learning. The first step is to use five predictive models to obtain predicted retention index values for the whole database. The second one is to compare these predicted values against the experimental ones. We consider a retention index erroneous if its accuracy (the difference between predicted and experimental value) is in the bottom 5% for each of the five models simultaneously. Using this method, we were able to detect 2093 outlier entries for standard and semi-standard non-polar stationary phases in the NIST 17 retention index database, 566 of those were corrected or removed by the developers in the NIST 20.

SIGNIFICANCE

This is a novel approach to find potentially erroneous entries in a large-scale database with mostly unique entries, which can be applied not only to retention indices. The procedure can help filter and report mishandled data to improve the quality of the dataset for machine learning applications and experimental use.

A simple and reliable QSPR model for prediction of chromatography retention indices of volatile organic compounds in peppers

Worldwide, various types of pepper are used in food as an additive due to their unique pungency, aroma, taste, and color. This spice is valued for its pungency contributed by the alkaloid piperine and aroma attributed to volatile essential oils. The essential oils are composed of volatile organic compounds (VOCs) in different concentrations and ratios. In chromatography, the identification of compounds is done by comparing obtained peaks with a reference standard. However, there are cases where reference standards are either unavailable or the chemical information of VOCs is not documented in reference libraries. To overcome these limitations, theoretical methodologies are applied to estimate the retention indices (RIs) of new VOCs. The aim of the present work is to develop a reliable QSPR model for the RIs of 273 identified VOCs of different types of pepper. Experimental retention indices were measured using comprehensive two-dimensional gas chromatography coupled to quadrupole mass spectrometry (GC × GC/qMS) using a coupled BPX5 and BP20 column system. The inbuilt Monte Carlo algorithm of CORAL software is used to generate QSPR models using the hybrid optimal descriptor extracted from a combination of SMILES and HFG (hydrogen-filled graph). The whole dataset of 273 VOCs is used to make ten splits, each of which is further divided into four sets: active training, passive training, calibration, and validation. The balance of correlation method with four target functions i.e. TF0 (WIIC = WCII = 0), TF1 (WIIC = 0.5 & WCII = 0), TF2 (WIIC = 0 & WCII = 0.3) and TF3 (WIIC = 0.5 & WCII = 0.3) is used. The results of the statistical parameters of each target function are compared with each other. The simultaneous application of the index of ideality of correlation (IIC) and correlation intensity index (CII) improves the predictive potential of the model. The best model is judged on the basis of the numerical value of R2 of the validation set. The statistical result of the best model for the validation set of split 6 computed with TF3 (WIIC = 0.5 & WCII = 0.3) is R2 = 0.9308, CCC = 0.9588, IIC = 0.7704, CII = 0.9549, Q2 = 0.9281 and RMSE = 0.544. The promoters of increase/decrease for RI are also extracted using the best model (split 6). Moreover, the proposed model was used for an external validation set.

The System of Self-Consistent Models: The Case of Henry's Law Constants

Data on Henry's law constants make it possible to systematize geochemical conditions affecting atmosphere status and consequently triggering climate changes. The constants of Henry's law are desired for assessing the processes related to atmospheric contaminations caused by pollutants. The most important are those that are capable of long-term movements over long distances. This ability is closely related to the values of Henry's law constants. Chemical changes in gaseous mixtures affect the fate of atmospheric pollutants and ecology, climate, and human health. Since the number of organic compounds present in the atmosphere is extremely large, it is desirable to develop models suitable for predictions for the large pool of organic molecules that may be present in the atmosphere. Here, we report the development of such a model for Henry's law constants predictions of 29,439 compounds using the CORAL software (2023). The statistical quality of the model is characterized by the value of the coefficient of determination for the training and validation sets of about 0.81 (on average).

The enhancement scheme for the predictive ability of QSAR: A case of mutagenicity

Mutagenicity is one of the most dangerous properties from the point of view of medicine and ecology. Experimental determination of mutagenicity remains a costly process, which makes it attractive to identify new hazardous compounds based on available experimental data through in silico methods or quantitative structure-activity relationships (QSAR). A system for constructing groups of random models is proposed for comparing various molecular features extracted from SMILES and graphs. For mutagenicity (mutagenicity values were expressed by the logarithm of the number of revertants per nanomole assayed by Salmonella typhimurium TA98-S9 microsomal preparation) models, the Morgan connectivity values are more informative than the comparison of quality for different rings in molecules. The resulting models were tested with the previously proposed model self-consistency system. The average value of the determination coefficient for the validation set is 0.8737 ± 0.0312.

Parsing structural fragments of thiazolidin-4-one based α-amylase inhibitors: A combined approach employing in vitro colorimetric screening and GA-MLR based QSAR modelling supported by molecular docking, molecular dynamics simulation and ADMET studies

α-Amylase (EC.3.2.1.1) is a ubiquitous digestive endoamylase. The abrupt rise in blood glucose levels due to the hydrolysis of carbohydrates by α-amylase at a faster rate is one of the main reasons for type 2 diabetes. The inhibitors prevent the action of digestive enzymes, slowing the digestion of carbs and eventually assisting in the management of postprandial hyperglycemia. In the course of developing α-amylase inhibitors, we have screened 2-aryliminothiazolidin-4-one based analogs for their in vitro α-amylase inhibitory potential and employed various in silico approaches for the detailed exploration of the bioactivity. The DNSA bioassay revealed that compounds 5c, 5e, 5h, 5j, 5m, 5o and 5t were more potent than the reference drug (IC₆₀ value = 22.94 ± 0.24 μg mL^-1). The derivative 5o with -NO₂ group at both the rings was the most potent analog with an IC₆₀ value of 19.67 ± 0.20 μg mL^-1 whereas derivative 5a with unsubstituted aromatic rings showed poor inhibitory potential with an IC₆₀ value of 33.40 ± 0.15 μg mL^-1. The reliable QSAR models were developed using the QSARINS software. The high value of R²_ext = 0.9632 for model IM-9 showed that the built model can be applied to predict the α-amylase inhibitory activity of the untested molecules. A consensus modelling approach was also employed to test the reliability and robustness of the developed QSAR models. Molecular docking and molecular dynamics were employed to validate the bioassay results by studying the conformational changes and interaction mechanisms. A step further, these compounds also exhibited good ADMET characteristics and bioavailability when tested for in silico pharmacokinetics prediction parameters.

Quantitative structure-activity relationship modeling for predication of inhibition potencies of imatinib derivatives using SMILES attributes

Chronic myelogenous leukemia (CML) which is resulted from the BCR-ABL tyrosine kinase (TK) chimeric oncoprotein, is a malignant clonal disorder of hematopoietic stem cells. Imatinib is used as an inhibitor of BCR-ABL TK in the treatment of CML patients. The main object of the present manuscript is focused on constructing quantitative activity relationships (QSARs) models for the prediction of inhibition potencies of a large series of imatinib derivatives against BCR-ABL TK. Herren, the inbuilt Monte Carlo algorithm of CORAL software is employed to develop QSAR models. The SMILES notations of chemical structures are used to compute the descriptor of correlation weights (CWs). QSAR models are established using the balance of correlation method with the index of ideality of correlation (IIC). The data set of 306 molecules is randomly divided into three splits. In QSAR modeling, the numerical value of R², Q², and IIC for the validation set of splits 1 to 3 are in the range of 0.7180-0.7755, 0.6891-0.7561, and 0.4431-0.8611 respectively. The numerical result of [Formula: see text] > 0.5 for all three constructed models in the Y-randomization test validate the reliability of established models. The promoters of increase/decrease for pIC₅₀ are recognized and used for the mechanistic interpretation of structural attributes.

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.

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.