Recently, the use of the Bayesian network as an alternative to existing tools for similarity-based virtual screening has received noticeable attention from researchers in the chemoinformatics field. The main aim of the Bayesian network model is to improve the retrieval effectiveness of similarity-based virtual screening. To this end, different models of the Bayesian network have been developed. In our previous works, the retrieval performance of the Bayesian network was observed to improve significantly when multiple reference structures or fragment weightings were used. In this article, the authors enhance the Bayesian inference network (BIN) using the relevance feedback information. In this approach, a few high-ranking structures of unknown activity were filtered from the outputs of BIN, based on a single active reference structure, to form a set of active reference structures. This set of active reference structures was used in two distinct techniques for carrying out such BIN searching: reweighting the fragments in the reference structures and group fusion techniques. Simulated virtual screening experiments with three MDL Drug Data Report data sets showed that the proposed techniques provide simple ways of enhancing the cost-effectiveness of ligand-based virtual screening searches, especially for higher diversity data sets.
The mammalian hyaluronidase degrades hyaluronic acid by the cleavage of the β-1,4-glycosidic bond furnishing a tetrasaccharide molecule as the main product which is a highly angiogenic and potent inducer of inflammatory cytokines. Ursolic acid 1, isolated from Prismatomeris tetrandra, was identified as having the potential to develop inhibitors of hyaluronidase. A series of ursolic acid analogues were either synthesized via structure modification of ursolic acid 1 or commercially obtained. The evaluation of the inhibitory activity of these compounds on the hyaluronidase enzyme was conducted. Several structural, topological and quantum chemical descriptors for these compounds were calculated using semi empirical quantum chemical methods. A quantitative structure activity relationship study (QSAR) was performed to correlate these descriptors with the hyaluronidase inhibitory activity. The statistical characteristics provided by the best multi linear model (BML) (R² = 0.9717, R²cv = 0.9506) indicated satisfactory stability and predictive ability of the developed model. The in silico molecular docking study which was used to determine the binding interactions revealed that the ursolic acid analog 22 had a strong affinity towards human hyaluronidase.
A quantitative structure-human intestinal absorption relationship was developed using artificial neural network (ANN) modeling. A set of 86 drug compounds and their experimentally-derived intestinal absorption values used in this study was gathered from the literature and a total of 57 global molecular descriptors, including constitutional, topological, chemical, geometrical and quantum chemical descriptors, calculated for each compound. A supervised network with radial basis transfer function was used to correlate calculated molecular descriptors with experimentally-derived measures of human intestinal absorption. A genetic algorithm was then used to select important molecular descriptors. Intestinal absorption values (IA%) were used as the ANN's output and calculated molecular descriptors as the inputs. The best genetic neural network (GNN) model with 15 input descriptors was chosen, and the significance of the selected descriptors for intestinal absorption examined. Results obtained with the model that was developed indicate that lipophilicity, conformational stability and inter-molecular interactions (polarity, and hydrogen bonding) have the largest impact on intestinal absorption.
Aromatherapy with essential oils (EOs) has been linked to improvement of cognitive function in patients with dementia. In order to act systemically, active EO components must be absorbed through the skin, enter the systemic circulation, and cross the blood brain barrier (BBB). Thus, the aim of this work was to develop quantitative structure activity relationships (QSARs), to predict skin and blood barrier penetrative abilities of 119 terpenoids from EOs used in aromatherapy. The first model was based on experimentally measured skin permeability for 162 molecules, and the second model on BBB permeability for 138 molecules. Each molecule was encoded with 63 calculated molecular descriptors and an artificial neural network was used to correlate molecular descriptors to permeabilities. Developed QSAR models confirm that EOs components penetrate through the skin and across the BBB. Some well-known descriptors, such as log P (lipophilicity), molecular size and shape, dominated the QSAR model for BBB permeability. Compounds with the highest predicted BBB penetration were hydrocarbon terpenes with the smallest molecular size and highest lipophilicity. Thus, molecular size is a limiting factor for penetration. Compounds with the highest skin permeability have slightly higher molecular size, high lipophilicity and low polarity. Our work shows that a major disadvantage of novel multitarget compounds developed for the treatment of Alzheimer's disease is the size of molecules, which cause problems in their delivery to the brain. Therefore, there is a need for smaller compounds, which possess more desirable physicochemical properties and pharmacokinetics, in addition to targeted biological effects.Communicated by Ramaswamy H. Sarma.
Many of the similarity-based virtual screening approaches assume that molecular fragments that are not related to the biological activity carry the same weight as the important ones. This was the reason that led to the use of Bayesian networks as an alternative to existing tools for similarity-based virtual screening. In our recent work, the retrieval performance of the Bayesian inference network (BIN) was observed to improve significantly when molecular fragments were reweighted using the relevance feedback information. In this paper, a set of active reference structures were used to reweight the fragments in the reference structure. In this approach, higher weights were assigned to those fragments that occur more frequently in the set of active reference structures while others were penalized. Simulated virtual screening experiments with MDL Drug Data Report datasets showed that the proposed approach significantly improved the retrieval effectiveness of ligand-based virtual screening, especially when the active molecules being sought had a high degree of structural heterogeneity.
The Plasmodium falciparum Lactate Dehydrogenase enzyme (PfLDH) catalyzes inter-conversion of pyruvate to lactate during glycolysis producing the energy required for parasitic growth. The PfLDH has been studied as a potential molecular target for development of anti-malarial agents. In an attempt to find the potent inhibitor of PfLDH, we have used Discovery studio to perform molecular docking in the active binding pocket of PfLDH by CDOCKER, followed by three-dimensional quantitative structure-activity relationship (3D-QSAR) studies of tricyclic guanidine batzelladine compounds, which were previously synthesized in our laboratory. Docking studies showed that there is a very strong correlation between in silico and in vitro results. Based on docking results, a highly predictive 3D-QSAR model was developed with q(2) of 0.516. The model has predicted r(2) of 0.91 showing that predicted IC50 values are in good agreement with experimental IC50 values. The results obtained from this study revealed the developed model can be used to design new anti-malarial compounds based on tricyclic guanidine derivatives and to predict activities of new inhibitors.
Time-varying binary gravitational search algorithm (TVBGSA) is proposed for predicting antidiabetic activity of 134 dipeptidyl peptidase-IV (DPP-IV) inhibitors. To improve the performance of the binary gravitational search algorithm (BGSA) method, we propose a dynamic time-varying transfer function. A new control parameter, μ , is added in the original transfer function as a time-varying variable. The TVBGSA-based model was internally and externally validated based on
Q
int
2
,
Q
L G O
2
,
Q
B o o t
2
,
M S
E
t r a i n
,
Q
e x t
2
,
M S
E
t e s t
, Y-randomization test, and applicability domain evaluation. The validation results indicate that the proposed TVBGSA model is robust and not due to chance correlation. The descriptor selection and prediction performance of TVBGSA outperform BGSA method. TVBGSA shows higher
Q
int
2
of 0.957,
Q
L G O
2
of 0.951,
Q
B o o t
2
of 0.954,
Q
e x t
2
of 0.938, and lower
M S
E
t r a i n
and
M S
E
t e s t
compared to obtained results by BGSA, indicating the best prediction performance of the proposed TVBGSA model. The results clearly reveal that the proposed TVBGSA method is useful for constructing reliable and robust QSARs for predicting antidiabetic activity of DPP-IV inhibitors prior to designing and experimental synthesizing of new DPP-IV inhibitors.
An improved binary differential search (improved BDS) algorithm is proposed for QSAR classification of diverse series of antimicrobial compounds against Candida albicans inhibitors. The transfer functions is the most important component of the BDS algorithm, and converts continuous values of the donor into discrete values. In this paper, the eight types of transfer functions are investigated to verify their efficiency in improving BDS algorithm performance in QSAR classification. The performance was evaluated using three metrics: classification accuracy (CA), geometric mean of sensitivity and specificity (G-mean), and area under the curve. The Kruskal-Wallis test was also applied to show the statistical differences between the functions. Two functions, S1 and V4, show the best classification achievement, with a slightly better performance of V4 than S1. The V4 function takes the lowest iterations and selects the fewest descriptors. In addition, the V4 function yields the best CA and G-mean of 98.07% and 0.977%, respectively. The results prove that the V4 transfer function significantly improves the performance of the original BDS.
A robust screening approach and a sparse quantitative structure-retention relationship (QSRR) model for predicting retention indices (RIs) of 169 constituents of essential oils is proposed. The proposed approach is represented in two steps. First, dimension reduction was performed using the proposed modified robust sure independence screening (MR-SIS) method. Second, prediction of RIs was made using the proposed robust sparse QSRR with smoothly clipped absolute deviation (SCAD) penalty (RSQSRR). The RSQSRR model was internally and externally validated based on [Formula: see text], [Formula: see text], [Formula: see text], [Formula: see text], Y-randomization test, [Formula: see text], [Formula: see text], and the applicability domain. The validation results indicate that the model is robust and not due to chance correlation. The descriptor selection and prediction performance of the RSQSRR for training dataset outperform the other two used modelling methods. The RSQSRR shows the highest [Formula: see text], [Formula: see text], and [Formula: see text], and the lowest [Formula: see text]. For the test dataset, the RSQSRR shows a high external validation value ([Formula: see text]), and a low value of [Formula: see text] compared with the other methods, indicating its higher predictive ability. In conclusion, the results reveal that the proposed RSQSRR is an efficient approach for modelling high dimensional QSRRs and the method is useful for the estimation of RIs of essential oils that have not been experimentally tested.
A penalized quantitative structure-property relationship (QSPR) model with adaptive bridge penalty for predicting the melting points of 92 energetic carbocyclic nitroaromatic compounds is proposed. To ensure the consistency of the descriptor selection of the proposed penalized adaptive bridge (PBridge), we proposed a ridge estimator ([Formula: see text]) as an initial weight in the adaptive bridge penalty. The Bayesian information criterion was applied to ensure the accurate selection of the tuning parameter ([Formula: see text]). The PBridge based model was internally and externally validated based on [Formula: see text], [Formula: see text], [Formula: see text], [Formula: see text], [Formula: see text], [Formula: see text], the Y-randomization test, [Formula: see text], [Formula: see text], [Formula: see text], [Formula: see text] and the applicability domain. The validation results indicate that the model is robust and not due to chance correlation. The descriptor selection and prediction performance of PBridge for the training dataset outperforms the other methods used. PBridge shows the highest [Formula: see text] of 0.959, [Formula: see text] of 0.953, [Formula: see text] of 0.949 and [Formula: see text] of 0.959, and the lowest [Formula: see text] and [Formula: see text]. For the test dataset, PBridge shows a higher [Formula: see text] of 0.945 and [Formula: see text] of 0.948, and a lower [Formula: see text] and [Formula: see text], indicating its better prediction performance. The results clearly reveal that the proposed PBridge is useful for constructing reliable and robust QSPRs for predicting melting points prior to synthesizing new organic compounds.
One of the recently developed metaheuristic algorithms, the coyote optimization algorithm (COA), has shown to perform better in a number of difficult optimization tasks. The binary form, BCOA, is used in this study as a solution to the descriptor selection issue in classifying diverse antifungal series. Z-shape transfer functions (ZTF) are evaluated to verify their efficiency in improving BCOA performance in QSAR classification based on classification accuracy (CA), the geometric mean of sensitivity and specificity (G-mean), and the area under the curve (AUC). The Kruskal-Wallis test is also applied to show the statistical differences between the functions. The efficacy of the best suggested transfer function, ZTF4, is further assessed by comparing it to the most recent binary algorithms. The results prove that ZTF, especially ZTF4, significantly improves the performance of the original BCOA. The ZTF4 function yields the best CA and G-mean of 99.03% and 0.992%, respectively. It shows the fastest convergence behaviour compared to other binary algorithms. It takes the fewest iterations to reach high classification performance and selects the fewest descriptors. In conclusion, the obtained results indicate the ability of the ZTF4-based BCOA to find the smallest subset of descriptors while maintaining the best classification accuracy performance.
Benzyl-o-vanillin and benzimidazole nucleus serve as important pharmacophore in drug discovery. The benzyl vanillin (2-(benzyloxy)-3-methoxybenzaldehyde) compound shows anti-proliferative activity in HL60 leukemia cancer cells and can effect cell cycle progression at G2/M phase. Its apoptosis activity was due to disruption of mitochondrial functioning. In this study, we have studied a series of compounds consisting of benzyl vanillin and benzimidazole structures. We hypothesize that by fusing these two structures we can produce compounds that have better anticancer activity with improved specificity particularly towards the leukemia cell line. Here we explored the anticancer activity of three compounds namely 2-(2-benzyloxy-3-methoxyphenyl)-1H-benzimidazole, 2MP, N-1-(2-benzyloxy-3-methoxybenzyl)-2-(2-benzyloxy-3-methoxyphenyl)-1H-benzimidazole, 2XP, and (R) and (S)-1-(2-benzyloxy-3-methoxyphenyl)-2, 2, 2-trichloroethyl benzenesulfonate, 3BS and compared their activity to 2-benzyloxy-3-methoxybenzaldehyde, (Bn1), the parent compound. 2XP and 3BS induces cell death of U937 leukemic cell line through DNA fragmentation that lead to the intrinsic caspase 9 activation. DNA binding study primarily by the equilibrium binding titration assay followed by the Viscosity study reveal the DNA binding through groove region with intrinsic binding constant 7.39 µM/bp and 6.86 µM/bp for 3BS and 2XP respectively. 2XP and 3BS showed strong DNA binding activity by the UV titration method with the computational drug modeling showed that both 2XP and 3BS failed to form any electrostatic linkages except via hydrophobic interaction through the minor groove region of the nucleic acid. The benzylvanillin alone (Bn1) has weak anticancer activity even after it was combined with the benzimidazole (2MP), but after addition of another benzylvanillin structure (2XP), stronger activity was observed. Also, the combination of benzylvanillin with benzenesulfonate (3BS) significantly improved the anticancer activity of Bn1. The present study provides a new insight of benzyl vanillin derivatives as potential anti-leukemic agent.
Peroxisome Proliferator-Activated Receptor γ (PPARγ) activators have drawn great recent attention in the clinical management of type 2 diabetes mellitus, prompting several attempts to discover and optimize new PPARγ activators. With this in mind, we explored the pharmacophoric space of PPARγ using seven diverse sets of activators. Subsequently, genetic algorithm and multiple linear regression analysis were employed to select an optimal combination of pharmacophoric models and 2D physicochemical descriptors capable of accessing self-consistent and predictive quantitative structure-activity relationship (QSAR) (r2(71)=0.80, F=270.3, r2LOO=0.73, r2PRESS against 17 external test inhibitors=0.67). Three orthogonal pharmacophores emerged in the QSAR equation and were validated by receiver operating characteristic (ROC) curves analysis. The models were then used to screen the national cancer institute (NCI) list of compounds. The highest-ranking hits were tested in vitro. The most potent hits illustrated EC50 values of 15 and 224 nM.
In high-dimensional quantitative structure-activity relationship (QSAR) modelling, penalization methods have been a popular choice to simultaneously address molecular descriptor selection and QSAR model estimation. In this study, a penalized linear regression model with L1/2-norm is proposed. Furthermore, the local linear approximation algorithm is utilized to avoid the non-convexity of the proposed method. The potential applicability of the proposed method is tested on several benchmark data sets. Compared with other commonly used penalized methods, the proposed method can not only obtain the best predictive ability, but also provide an easily interpretable QSAR model. In addition, it is noteworthy that the results obtained in terms of applicability domain and Y-randomization test provide an efficient and a robust QSAR model. It is evident from the results that the proposed method may possibly be a promising penalized method in the field of computational chemistry research, especially when the number of molecular descriptors exceeds the number of compounds.
A high-dimensional quantitative structure-activity relationship (QSAR) classification model typically contains a large number of irrelevant and redundant descriptors. In this paper, a new design of descriptor selection for the QSAR classification model estimation method is proposed by adding a new weight inside L1-norm. The experimental results of classifying the anti-hepatitis C virus activity of thiourea derivatives demonstrate that the proposed descriptor selection method in the QSAR classification model performs effectively and competitively compared with other existing penalized methods in terms of classification performance on both the training and the testing datasets. Moreover, it is noteworthy that the results obtained in terms of stability test and applicability domain provide a robust QSAR classification model. It is evident from the results that the developed QSAR classification model could conceivably be employed for further high-dimensional QSAR classification studies.
High-dimensionality is one of the major problems which affect the quality of the quantitative structure-activity relationship (QSAR) modelling. Obtaining a reliable QSAR model with few descriptors is an essential procedure in chemometrics. The binary grasshopper optimization algorithm (BGOA) is a new meta-heuristic optimization algorithm, which has been used successfully to perform feature selection. In this paper, four new transfer functions were adapted to improve the exploration and exploitation capability of the BGOA in QSAR modelling of influenza A viruses (H1N1). The QSAR model with these new quadratic transfer functions was internally and externally validated based on MSEtrain, Y-randomization test, MSEtest, and the applicability domain (AD). The validation results indicate that the model is robust and not due to chance correlation. In addition, the results indicate that the descriptor selection and prediction performance of the QSAR model for training dataset outperform the other S-shaped and V-shaped transfer functions. QSAR model using quadratic transfer function shows the lowest MSEtrain. For the test dataset, proposed QSAR model shows lower value of MSEtest compared with the other methods, indicating its higher predictive ability. In conclusion, the results reveal that the proposed QSAR model is an efficient approach for modelling high-dimensional QSAR models and it is useful for the estimation of IC50 values of neuraminidase inhibitors that have not been experimentally tested.
One of the most challenging issues when facing a Quantitative structure-activity relationship (QSAR) classification model is to deal with the descriptor selection. Penalized methods have been adapted and have gained popularity as a key for simultaneously performing descriptor selection and QSAR classification model estimation. However, penalized methods have drawbacks such as having biases and inconsistencies that make they lack the oracle properties. This paper proposes an adaptive penalized logistic regression (APLR) to overcome these drawbacks. This is done by employing a ratio (BWR) of the descriptors between-groups sum of squares (BSS) to the within-groups sum of squares (WSS) for each descriptor as a weight inside the L1-norm. The proposed method was applied to one dataset that consists of a diverse series of antimicrobial agents with their respective bioactivities against Candida albicans. By experimental study, it has been shown that the proposed method (APLR) was more efficient in the selection of descriptors and classification accuracy than the other competitive methods that could be used in developing QSAR classification models. Another dataset was also successfully experienced. Therefore, it can be concluded that the APLR method had significant impact on QSAR analysis and studies.
The horse herd optimization algorithm (HOA), one of the more contemporary metaheuristic algorithms, has demonstrated superior performance in a number of challenging optimization tasks. In the present work, the descriptor selection issue is resolved by classifying different essential oil retention indices using the binary form, BHOA. Based on internal and external prediction criteria, Z-shape transfer functions (ZTF) were tested to verify their efficiency in improving BHOA performance in QSPR modelling for predicting retention indices of essential oils. The evaluation criteria involved the mean-squared error of the training and testing datasets (MSE), and leave-one-out internal and external validation (Q2). The degree of convergence of the proposed Z-shaped transfer functions was compared. In addition, K-fold cross validation with k = 5 was applied. The results show that ZTF, especially ZTF1, greatly improves the performance of the original BHOA. Comparatively speaking, ZTF, especially ZTF1, exhibits the fastest convergence behaviour of the binary algorithms. It chooses the fewest descriptors and requires the fewest iterations to achieve excellent prediction performance.
A series of 3-amino-2-hydroxybenzofused 2-phosphalactones (4a-l) has been synthesized from the Kabachnik-Fields reaction via a facile route from a one-pot three-component reaction of diphenylphosphite with various 2-hydroxybenzaldehyes and heterocyclic amines in a new way of expansion. The in vitro anti-cell proliferation studies by MTT assay have revealed them as potential Panc-1, Miapaca-2, and BxPC-3 pancreatic cell growth inhibitors, and the same is supported by molecular docking, QSAR, and ADMET studies. The MTT assay of their SAHA derivatives against the same cell lines evidenced them as potential HDAC inhibitors and identified 4a, 4b, and 4k substituted with 1,3-thiazol, 1,3,4-thiadiazol, and 5-sulfanyl-1,3,4-thiadiazol moieties on phenyl and diethylamino phenyl rings as potential ones. Additionally, the flow cytometric analyses of 4a, 4b, and 4k against BxPC-3 cells revealed compound 4k as a lead compound that arrests the S phase cell cycle growth at low micromolar concentrations. The ADMET properties have ascertained their inherent pharmacokinetic potentiality, and the wholesome results prompted us to report it as the first study on anti-pancreatic cancer activity of cyclic α-aminophosphonates. Ultimately, this study serves as a good contribution to update the existing knowledge on the anticancer organophosphorus heterocyclic compounds and elevates the scope for generation of new anticancer drugs. Further, the studies like QSAR, drug properties, toxicity risks, and bioactivity scores predicted for them have ascertained the synthesized compounds as newer and potential drug candidates. Hence, this study had augmented the array of α-aminophosphonates by adding a new collection of 3-amino-2-hydroxybenzofused 2-phosphalactones, a class of cyclic α-aminophosphonates, to it, which proved them as potential anti-pancreatic cancer agents.
In advanced oxidation processes (AOPs), the aqueous hydroxyl radical (HO) acts as a strong oxidant to react with organic contaminants. The hydroxyl radical rate constant (kHO) is important for evaluating and modelling of the AOPs. In this study, quantitative structure-property relationship (QSPR) method is applied to model the hydroxyl radical rate constant for a diverse dataset of 457 water contaminants from 27 various chemical classes. The constricted binary particle swarm optimization and multiple-linear regression (BPSO-MLR) are used to obtain the best model with eight theoretical descriptors. An optimized feed forward neural network (FFNN) is developed to investigate the complex performance of the selected molecular parameters with kHO. Although the FFNN prediction results are more accurate than those obtained using BPSO-MLR, the application of the latter is much more convenient. Various internal and external validation techniques indicate that the obtained models could predict the logarithmic hydroxyl radical rate constants of a large number of water contaminants with less than 4% absolute relative error. Finally, the above-mentioned proposed models are compared to those reported earlier and the structural factors contributing to the AOP degradation efficiency are discussed.