1
|
Azimi A, Ahmadi S, Javan MJ, Rouhani M, Mirjafary Z. QSAR models for the ozonation of diverse volatile organic compounds at different temperatures. RSC Adv 2024; 14:8041-8052. [PMID: 38454938 PMCID: PMC10918768 DOI: 10.1039/d3ra08805g] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/24/2023] [Accepted: 02/06/2024] [Indexed: 03/09/2024] Open
Abstract
In order to assess the fate and persistence of volatile organic compounds (VOCs) in the atmosphere, it is necessary to determine their oxidation rate constants for their reaction with ozone (kO3). However, given that experimental values of kO3 are only available for a few hundred compounds and their determination is expensive and time-consuming, developing predictive models for kO3 is of great importance. Thus, this study aimed to develop reliable quantitative structure-activity relationship (QSAR) models for 302 values of 149 VOCs across a broad temperature range (178-409 K). The model was constructed based on the combination of a simplified molecular-input line-entry system (SMILES) and temperature as an experimental condition, namely quasi-SMILES. In this study, temperature was incorporated in the models as an independent feature. The hybrid optimal descriptor generated from the combination of quasi-SMILES and HFG (hydrogen-filled graph) was used to develop reliable, accurate, and predictive QSAR models employing the CORAL software. The balance between the correlation method and four different target functions (target function without considering IIC or CII, target function using each IIC or CII, and target function based on the combination of IIC and CII) was used to improve the predictability of the QSAR models. The performance of the developed models based on different target functions was compared. The correlation intensity index (CII) significantly enhanced the predictability of the model. The best model was selected based on the numerical value of Rm2 of the calibration set (split #1, Rtrain2 = 0.9834, Rcalibration2 = 0.9276, Rvalidation2 = 0.9136, and calibration = 0.8770). The promoters of increase/decrease for log kO3 were also computed based on the best model. The presence of a double bond (BOND10000000 and $10 000 000 000), absence of halogen (HALO00000000), and the nearest neighbor codes for carbon equal to 321 (NNC-C⋯321) are some significant promoters of endpoint increase.
Collapse
Affiliation(s)
- Ali Azimi
- Department of Chemistry, Science and Research Branch, Islamic Azad University Tehran Iran
| | - Shahin Ahmadi
- Department of Pharmaceutical Chemistry, Faculty of Pharmaceutical Chemistry, Tehran Medical Sciences, Islamic Azad University Tehran Iran
| | - Marjan Jebeli Javan
- Department of Pharmaceutical Chemistry, Faculty of Pharmaceutical Chemistry, Tehran Medical Sciences, Islamic Azad University Tehran Iran
| | - Morteza Rouhani
- Department of Chemistry, Science and Research Branch, Islamic Azad University Tehran Iran
| | - Zohreh Mirjafary
- Department of Chemistry, Science and Research Branch, Islamic Azad University Tehran Iran
| |
Collapse
|
2
|
Ahmadi S, Lotfi S, Hamzehali H, Kumar P. A simple and reliable QSPR model for prediction of chromatography retention indices of volatile organic compounds in peppers. RSC Adv 2024; 14:3186-3201. [PMID: 38249679 PMCID: PMC10797599 DOI: 10.1039/d3ra07960k] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/21/2023] [Accepted: 01/03/2024] [Indexed: 01/23/2024] Open
Abstract
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.
Collapse
Affiliation(s)
- Shahin Ahmadi
- Department of Pharmaceutical Chemistry, Tehran Medical Sciences, Islamic Azad University Tehran Iran
| | - Shahram Lotfi
- Department of Chemistry, Payame Noor University (PNU) 19395-4697 Tehran Iran
| | - Hamideh Hamzehali
- Department of Chemistry, Islamic Azad University East Tehran Branch Tehran Iran
| | - Parvin Kumar
- Department of Chemistry, Kurukshetra University Kurukshetra Haryana 136119 India
| |
Collapse
|
3
|
Tajiani F, Ahmadi S, Lotfi S, Kumar P, Almasirad A. In-silico activity prediction and docking studies of some flavonol derivatives as anti-prostate cancer agents based on Monte Carlo optimization. BMC Chem 2023; 17:87. [PMID: 37496005 PMCID: PMC10373329 DOI: 10.1186/s13065-023-00999-y] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/02/2023] [Accepted: 06/30/2023] [Indexed: 07/28/2023] Open
Abstract
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 pIC50 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 pIC50 of eight natural flavonols with pIC50 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).
Collapse
Affiliation(s)
- Faezeh Tajiani
- Department of Medicinal Chemistry, Faculty of Pharmacy, Tehran Medical Sciences, Islamic Azad University, Tehran, Iran
| | - Shahin Ahmadi
- Department of Pharmaceutical Chemistry, Faculty of Pharmaceutical Chemistry, Tehran Medical Sciences, Islamic Azad University, Tehran, Iran.
| | - Shahram Lotfi
- Department of Chemistry, Payame Noor University (PNU), Tehran, 19395-4697, Iran
| | - Parvin Kumar
- Department of Chemistry, Kurukshetra University, Kurukshetra, Haryana, 136119, India
| | - Ali Almasirad
- Department of Medicinal Chemistry, Faculty of Pharmacy, Tehran Medical Sciences, Islamic Azad University, Tehran, Iran
| |
Collapse
|
4
|
Quantitative structure-activity relationship modeling for predication of inhibition potencies of imatinib derivatives using SMILES attributes. Sci Rep 2022; 12:21708. [PMID: 36522400 PMCID: PMC9755126 DOI: 10.1038/s41598-022-26279-8] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/11/2022] [Accepted: 12/13/2022] [Indexed: 12/23/2022] Open
Abstract
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 R2, Q2, 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 pIC50 are recognized and used for the mechanistic interpretation of structural attributes.
Collapse
|
5
|
|
6
|
QSAR modelling, molecular docking studies and ADMET predictions of polysubstituted pyridinylimidazoles as dual inhibitors of JNK3 and p38α MAPK. J Mol Struct 2022. [DOI: 10.1016/j.molstruc.2022.133504] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022]
|
7
|
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. J Mol Struct 2022. [DOI: 10.1016/j.molstruc.2022.133437] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023]
|
8
|
Lotfi S, Ahmadi S, Kumar P. Ecotoxicological prediction of organic chemicals toward Pseudokirchneriella subcapitata by Monte Carlo approach. RSC Adv 2022; 12:24988-24997. [PMID: 36199875 PMCID: PMC9434604 DOI: 10.1039/d2ra03936b] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/26/2022] [Accepted: 08/19/2022] [Indexed: 11/21/2022] Open
Abstract
In the ecotoxicological risk assessment, acute toxicity is one of the most significant criteria. Green alga Pseudokirchneriella subcapitata has been used for ecotoxicological studies to assess the toxicity of different toxic chemicals in freshwater. Quantitative Structure Activity Relationships (QSAR) are mathematical models to relate chemical structure and activity/physicochemical properties of chemicals quantitatively. Herein, Quantitative Structure Toxicity Relationship (QSTR) modeling is applied to assess the toxicity of a data set of 334 different chemicals on Pseudokirchneriella subcapitata, in terms of EC10 and EC50 values. The QSTR models are established using CORAL software by utilizing the target function (TF2) with the index of ideality of correlation (IIC). A hybrid optimal descriptor computed from SMILES and molecular hydrogen-suppressed graphs (HSG) is employed to construct QSTR models. The results of various statistical parameters of the QSTR model developed for pEC10 and pEC50 range from excellent to good and are in line with the standard parameters. The models prepared with IIC for Split 3 are chosen as the best model for both endpoints (pEC10 and pEC50). The numerical value of the determination coefficient of the validation set of split 3 for the endpoint pEC10 is 0.7849 and for the endpoint pEC50, it is 0.8150. The structural fractions accountable for the toxicity of chemicals are also extracted. The hydrophilic attributes like 1…n…(… and S…(…[double bond, length as m-dash]… exert positive contributions to controlling the aquatic toxicity and reducing algal toxicity, whereas attributes such as c…c…c…, C…C…C… enhance lipophilicity of the molecules and consequently enhance algal toxicity.
Collapse
Affiliation(s)
- Shahram Lotfi
- Department of Chemistry, Payame Noor University (PNU) 19395-4697 Tehran Iran
| | - Shahin Ahmadi
- Department of Pharmaceutical Chemistry, Faculty of Pharmaceutical Chemistry, Tehran Medical Sciences, Islamic Azad University Tehran Iran
| | - Parvin Kumar
- Department of Chemistry, Kurukshetra University Kurukshetra Haryana 136119 India
| |
Collapse
|
9
|
Kumar P, Kumar A, Singh D. 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. ENVIRONMENTAL TOXICOLOGY AND PHARMACOLOGY 2022; 93:103893. [PMID: 35654373 DOI: 10.1016/j.etap.2022.103893] [Citation(s) in RCA: 15] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/04/2022] [Revised: 05/21/2022] [Accepted: 05/26/2022] [Indexed: 06/15/2023]
Abstract
In the present study, ninety-five halogenated dioxins and related chemicals (dibenzo-p-dioxins, dibenzofurans, biphenyls, and naphthalene) with endpoint pEC50 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. TF1 (WIIC=WCII=0), TF2 (WIIC= 0.3 & WCII=0) and TF3 (WIIC= 0.0 &WCII=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 TF3 was found highest (RValid2=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 pEC50. Three random splits (Split 1, Split 2 and Split 4) were selected for the extraction of the promoter of increase/decrease for pEC50. 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.
Collapse
Affiliation(s)
- Parvin Kumar
- Department of Chemistry, Kurukshetra University, Kurukshetra, Haryana 136119, India.
| | - Ashwani Kumar
- Department of Pharmaceutical Sciences, Guru Jambheshwar University of Science and Technology, Hisar, Haryana 125001, India.
| | - Devender Singh
- Department of Chemistry, Maharshi Dayanand University, Rohtak, Haryana 124001, India
| |
Collapse
|
10
|
Singh R, Kumar P, Devi M, Lal S, Kumar A, Sindhu J, Toropova AP, Toropov AA, Singh D. Monte Carlo based QSGFEAR: prediction of Gibb's free energy of activation at different temperatures using SMILES based descriptors. NEW J CHEM 2022. [DOI: 10.1039/d2nj03515d] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
Monte Carlo optimization based QSGFEAR model development using CII results in the formation of more reliable, robust and predictive models.
Collapse
Affiliation(s)
- Rahul Singh
- Department of Chemistry, Kurukshetra University, Kurukshetra-136119, India
| | - Parvin Kumar
- Department of Chemistry, Kurukshetra University, Kurukshetra-136119, India
| | - Meena Devi
- Department of Chemistry, Kurukshetra University, Kurukshetra-136119, India
| | - Sohan Lal
- Department of Chemistry, Kurukshetra University, Kurukshetra-136119, India
| | - Ashwani Kumar
- Department of Pharmaceutical Sciences, GJUS&T, Hisar, 125001, India
| | - Jayant Sindhu
- Department of Chemistry, COBS&H, CCS Haryana Agricultural University, Hisar, 125004, India
| | - Alla P. Toropova
- Department of Environmental Health Science, Istituto di Ricerche Farmacologiche Mario Negri IRCCS, Via Mario Negri 2, 20156 Milano, Italy
| | - Andrey A. Toropov
- Department of Environmental Health Science, Istituto di Ricerche Farmacologiche Mario Negri IRCCS, Via Mario Negri 2, 20156 Milano, Italy
| | - Devender Singh
- Department of Chemistry, Maharshi Dayanand University, Rohtak, 124001, India
| |
Collapse
|
11
|
Ahmadi S, Lotfi S, Afshari S, Kumar P, Ghasemi E. CORAL: Monte Carlo based global QSAR modelling of Bruton tyrosine kinase inhibitors using hybrid descriptors. SAR AND QSAR IN ENVIRONMENTAL RESEARCH 2021; 32:1013-1031. [PMID: 34875951 DOI: 10.1080/1062936x.2021.2003429] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/03/2021] [Accepted: 11/03/2021] [Indexed: 06/13/2023]
Abstract
Global QSAR modelling was performed to predict the pIC50 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. TF1 (WIIC = 0) and TF2 (WIIC = 0.2) have been constructed. The statistical parameters of QSAR models computed by TF2 are more reliable and robust and were used to predict the pIC50 values. The model constructed for split 4 via TF2 is regarded as the best model and the numerical values of r2Train, r2Valid, Q2Train and Q2Valid 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 pIC50 of BTK inhibitors were also identified.
Collapse
Affiliation(s)
- S Ahmadi
- Department of Pharmaceutical Chemistry, Faculty of Pharmaceutical Chemistry, Tehran Medical Sciences, Islamic Azad University, Tehran, Iran
| | - S Lotfi
- Department of Chemistry, Payame Noor University (PNU), Tehran, Iran
| | - S Afshari
- Department of Medicinal Chemistry, Faculty of Pharmacy, Tehran Medical Sciences, Islamic Azad University, Tehran, Iran
| | - P Kumar
- Department of Chemistry, Kurukshetra University, Kurukshetra, Haryana, India
| | - E Ghasemi
- Department of Pharmaceutical Chemistry, Faculty of Pharmaceutical Chemistry, Tehran Medical Sciences, Islamic Azad University, Tehran, Iran
| |
Collapse
|