Comparative QSAR studies using HQSAR,CoMFA, and CoMSIA methods on cyclic sulfone hydroxyethylamines as BACE1 inhibitors

The inhibition of β-secretase (BACE1) is currently the main pharmacological strategy available for Alzheimer’s disease (AD). 2D QSAR and 3D QSAR analysis on some cyclic sulfone hydroxyethylamines inhibitors against β-secretase (IC₅₀: 0.002–2.75 μM) were carried out using hologram QSAR (HQSAR), comparative molecular field analysis (CoMFA), and comparative molecular similarity indices analysis (CoMSIA) methods. The best model based on the training set was generated with a HQSAR q² value of 0.693 and r² value of 0.981; a CoMFA q² value of 0.534 and r² value of 0.913; and a CoMSIA q² value of 0.512 and r² value of 0.973. In order to gain further understand of the vital interactions between cyclic sulfone hydroxyethylamines and the protease, the analysis was performed by combining the CoMFA and CoMSIA field distributions with the active sites of the BACE1. The final QSAR models could be helpful in the design and development of novel active BACE1 inhibitors.     

QSAR analyses on avian influenza virus neuraminidase inhibitors using CoMFA, CoMSIA, and HQSAR

The recent wide spreading of the H5N1 avian influenza virus (AIV) in Asia, Europe and Africa and its ability to cause fatal infections in human has raised serious concerns about a pending global flu pandemic. Neuraminidase (NA) inhibitors are currently the only option for treatment or prophylaxis in humans infected with this strain. However, drugs currently on the market often meet with rapidly emerging resistant mutants and only have limited application as inadequate supply of synthetic material. To dig out helpful information for designing potent inhibitors with novel structures against the NA, we used automated docking, CoMFA, CoMSIA, and HQSAR methods to investigate the quantitative structure-activity relationship for 126 NA inhibitors (NIs) with great structural diversities and wide range of bioactivities against influenza A virus. Based on the binding conformations discovered via molecular docking into the crystal structure of NA, CoMFA and CoMSIA models were successfully built with the cross-validated q (2) of 0.813 and 0.771, respectively. HQSAR was also carried out as a complementary study in that HQSAR technique does not require 3D information of these compounds and could provide a detailed molecular fragment contribution to the inhibitory activity. These models also show clearly how steric, electrostatic, hydrophobicity, and individual fragments affect the potency of NA inhibitors. In addition, CoMFA and CoMSIA field distributions are found to be in well agreement with the structural characteristics of the corresponding binding sites. Therefore, the final 3D-QSAR models and the information of the inhibitor-enzyme interaction should be useful in developing novel potent NA inhibitors.     

QSAR and Pharmacophore Studies of Thiazolidine-4-carboxylic Acid Derivatives as Novel Influenza Neuraminidase Inhibitors Using HQSAR, Topomer CoMFA and CoMSIA

In order to understand the chemical-biological interactions governing their activities toward neuraminidase (NA), QSAR models of 28 thiazolidine-4-carboxylic acid derivatives with inhibitory influenza A virus were developed. The obtained HQSAR (hologram quantitative structure activity relationship), Topomer CoMFA and CoMSIA (comparative molecular similarity indices analysis) models were robust and had good exterior predictive capabilities. Moreover, QSAR modeling results elucidated that hydrogen bonds highly contributed to the inhibitory activity, then electrostatic and hydrophobic factors. Squared multiple correlation coefficients (R2) of HQSAR, Topomer CoMFA and CoMSIA models were 0.994, 0.978 and 0.996, respectively. Squared cross-validated correlation coefficients (Q2) of HQSAR, Topomer CoMFA and CoMSIA models were in turn 0.951, 919 and 0.820. Furthermore, squared multiple correlation coefficients for the test set (R2test) of HQSAR, CoMFA and CoMSIA models were 0.879, 0.912 and 0.953, respectively. Squared cross-validated correlation coefficients for the test set (Q2ext) of HQSAR, Topomer CoMFA and CoMSIA models were 0.867, 0.884 and 0.899, correspondingly.     

CoMFA,CoMSIA, Topomer CoMFA,HQSAR, molecular docking and molecular dynamics simulations study of triazine morpholino derivatives as mTOR inhibitors for the treatment of breast cancer

mTOR has become a promising target for many types of cancer like breast, lung and renal cell carcinoma. CoMFA, CoMSIA, Topomer CoMFA and HQSAR were performed on the series of 39 triazine morpholino derivatives. CoMFA analysis showed q² value of 0.735, r²_cv value of 0.722 and r²_pred value of 0.769. CoMSIA analysis (SEHD) showed q² value of 0.761, r²_cv value of 0.775 and r²_pred value of 0.651. Topomer CoMFA analysis showed q² value of 0.693, r² (conventional correlation coefficient) value of 0.940 and r²_pred value of 0.720. HQSAR analysis showed q²,r²and r²_pred values of 0.694, 0.920 and 0.750, respectively. HQSAR analysis with the combination of atomic number (A), bond type (B) and atomic connections showed q² and r² values of 0.655 and 0.891, respectively. Contour maps from all studies provided significant insights. Molecular docking studies with molecular dynamics simulations were carried out on the highly potent compound 36. Furthermore, four acridine derivatives were designed and docking results of these designed compounds showed the same interactions as that of the standard PI-103 which proved the efficiency of 3D-QSAR and MD/MS study. In future, this study might be useful prior to synthesis for the designing of novel mTOR inhibitors.     

Structure-based CoMFA and CoMSIA study of indolinone inhibitors of PDK1

Phosphoinositide-dependent protein kinase-1 (PDK1) is a Ser/Thr kinase which phosphorylates and activates members of the AGC kinase group known to control processes such as tumor cell growth, protection from apoptosis, and tumor angiogenesis. In this paper, CoMFA and CoMSIA studies were carried out on a training set of 56 conformationally rigid indolinone inhibitors of PDK1. Predictive 3D QSAR models, established using atom fit alignment rule based on crystallographic-bound conformation, had cross-validated (r _cv²) values of 0.738 and 0.816 and non-cross-validated (r _ncv²) values of 0.912 and 0.949 for CoMFA and CoMSIA models, respectively. The predictive ability of the CoMFA and CoMSIA models was determined using a test set of 14 compounds, which gave predictive correlation coefficients (r _pred²) of 0.865 and 0.837, respectively. Structure-based interpretation of the CoMFA and CoMSIA field properties provided further insights for the rational design of new PDK1 inhibitors. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

CoMFA and CoMSIA studies on C-aryl glucoside SGLT2 inhibitors as potential anti-diabetic agents

SGLT2 has become a target of therapeutic interest in diabetes research. CoMFA and CoMSIA studies were performed on C-aryl glucoside SGLT2 inhibitors (180 analogues) as potential anti-diabetic agents. Three different alignment strategies were used for the compounds. The best CoMFA and CoMSIA models were obtained by means of Distill rigid body alignment of training and test sets, and found statistically significant with cross-validated coefficients (q ²) of 0.602 and 0.618, respectively, and conventional coefficients (r ²) of 0.905 and 0.902, respectively. Both models were validated by a test set of 36 compounds giving satisfactory predicted correlation coefficients (r ² _pred) of 0.622 and 0.584 for CoMFA and CoMSIA models, respectively. A comparison was made with earlier 3D QSAR study on SGLT2 inhibitors, which shows that our 3D QSAR models are better than earlier models to predict good inhibitory activity. CoMFA and CoMSIA models generated in this work can provide useful information to design new compounds and helped in prediction of activity prior to synthesis.     

Structural investigations of anthranilimide derivatives by CoMFA and CoMSIA 3D-QSAR studies reveal novel insight into their structures toward glycogen phosphorylase inhibition

In the present work, three-dimensional quantitative structure–activity relationship (3-D QSAR) studies on a set of 70 anthranilimide compounds has been performed using docking-based as well as substructure-based molecular alignments. This resulted in the selection of more statistically relevant substructure-based alignment for further studies. Further, molecular models with good predictive power were derived using CoMFA (r ^2?=?0.997; Q ^2?=?0.578) and CoMSIA (r ^2?=?0.976; Q ^2?=?0.506), for predicting the biological activity of new compounds. The so-developed contour plots identified several key features of the compounds explaining wide activity ranges. Based on the information derived from the CoMFA contour maps, novel leads were proposed which showed better predicted activity with respect to the already reported systems. Thus, the present study not only offers a highly significant predictive QSAR model for anthranilimide derivatives as glycogen phosphorylase (GP) inhibitors which can eventually assist and complement the rational drug-design attempts, but also proposes a highly predictive pharmacophore model as a guide for further development of selective and more potent GP inhibitors as anti-diabetic agents.     

CoMFA and CoMSIA 3D-quantitative structure-activity relationship model on benzodiazepine derivatives, inhibitors of phosphodiesterase IV

Recently, we reported structurally novel PDE4 inhibitors based on 1,4-benzodiazepine derivatives. The main interest in developing bezodiazepine-based PDE4 inhibitors is in their lack of adverse effects of emesis with respect to rolipram-like compounds. A large effort has thus been made toward the structural optimization of this series. In the absence of structural information on the inhibitor binding mode into the PDE4 active site, 2D-QSAR (H-QSAR) and two 3D-QSAR (CoMFA and CoMSIA) methods were applied to improve our understanding of the molecular mechanism controlling the PDE4 affinity of the benzodiazepine derivatives. As expected, the CoMSIA 3D contour maps have provided more information on the benzodiazepine interaction mode with the PDE4 active site whereas CoMFA has built the best tool for activity prediction. The 2D pharmacophoric model derived from CoMSIA fields is consistent with the crystal structure of the PDE4 active site reported recently. The combination of the 2D and 3D-QSAR models was used not only to predict new compounds from the structural optimization process, but also to screen a large library of bezodiazepine derivatives.     

QSAR studies on hepatitis C virus NS5A protein tetracyclic inhibitors in wild type and mutants by CoMFA and CoMSIA

ABSTRACT

Several 3D-QSAR models were built based on 196 hepatitis C virus (HCV) NS5A protein inhibitors. The bioactivity values EC₉₀ for three types of inhibitors, the wild type (GT1a) and two mutants (GT1a Y93H and GT1a L31V), were collected to build three datasets. The programs OMEGA and ROCS were used for generating conformations and aligning molecules of the dataset, respectively. Each dataset was randomly divided into a training set and a test set three times to reduce the contingency of only one random selection. QSAR models were computed by comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA). For the datasets GT1a, GT1a Y93H, and GT1a L31V, the best models CoMFA-INDX, CoMSIA-SEHA, and CoMSIA-SEHA showed an r² value of 0.682 ± 0.033, 0.779 ± 0.036, and 0.782 ± 0.022 on the test sets, respectively. From the contour maps of the three best models, we summarized the favourable and unfavourable substituents on the tetracyclic core, the Z group, the proline group, and the valine group of inhibitors. We guessed the mutants could change the electrostatic surfaces of the wild type active pocket. In addition, we used ECFP analyses to find important substructures and could intuitively understand the results from QSAR models.     

3D-QSAR and molecular docking study of LRRK2 kinase inhibitors by CoMFA and CoMSIA methods

Three-dimensional quantitative structure–activity relationship (3D-QSAR) modelling was conducted on a series of leucine-rich repeat kinase 2 (LRRK2) antagonists using CoMFA and CoMSIA methods. The data set, which consisted of 37 molecules, was divided into training and test subsets by using a hierarchical clustering method. Both CoMFA and CoMSIA models were derived using a training set on the basis of the common substructure-based alignment. The optimum PLS model built by CoMFA and CoMSIA provided satisfactory statistical results (q² = 0.589 and r² = 0.927 and q² = 0.473 and r² = 0.802, respectively). The external predictive ability of the models was evaluated by using seven compounds. Moreover, an external evaluation set with known experimental data was used to evaluate the external predictive ability of the porposed models. The statistical parameters indicated that CoMFA (after region focusing) has high predictive ability in comparison with standard CoMFA and CoMSIA models. Molecular docking was also performed on the most active compound to investigate the existence of interactions between the most active inhibitor and the LRRK2 receptor. Based on the obtained results and CoMFA contour maps, some features were introduced to provide useful insights for designing novel and potent LRRK2 inhibitors.     

QSAR models of cytochrome P450 enzyme 1A2 inhibitors using CoMFA,CoMSIA and HQSAR

Quantitative structure–activity relationship (QSAR) studies were conducted on an in-house database of cytochrome P450 enzyme 1A2 inhibitors using the comparative molecular field analysis (CoMFA), comparative molecular similarity analysis (CoMSIA) and hologram QSAR (HQSAR) approaches. The database consisted of 36 active molecules featuring varied core structures. The model based on the naphthalene substructure alignment incorporating 19 molecules yielded the best model with a CoMFA cross validation value q² of 0.667 and a Pearson correlation coefficient r² of 0.976; a CoMSIA q² value of 0.616 and r² value of 0.985; and a HQSAR q² value of 0.652 and r² value of 0.917. A second model incorporating 34 molecules aligned using the benzene substructure yielded an acceptable CoMFA model with q² value of 0.5 and r² value of 0.991. Depending on the core structure of the molecule under consideration, new CYP1A2 inhibitors will be designed based on the results from these models.     

CoMFA and CoMSIA studies on a new series of xanthone derivatives against the oral human epidermoid carcinoma (KB) cancer cell line

Abstract  A new series of xanthone derivatives against the oral human epidermoid carcinoma (KB) cancer cell line is examined to determine the relationship between the structural properties and the biological activity of these compounds—the 3-D quantitative structure–activity relationship (3D-QSAR)—using comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA). The best CoMFA and CoMSIA models were obtained using the atom-based alignment of 33 compounds, 22 training compounds and 11 tested compounds, and these give desirable statistics; those for the CoMFA standard model were: r _cv² = 0.691, r ² = 0.998, S _press = 0.178, s = 0.014 and F = 1080.765, while CoMSIA combined steric, electrostatic, hydrophobic and hydrogen-bond acceptor fields: r _cv² = 0.600, r ² = 0.988, S _press = 0.206, s = 0.034 and F = 284.433. The 3D-QSAR models calculated satisfactory test set activities. The 3D-QSAR contour plots correlated strongly with the experimental data for the binding topology. For this reason, these results would be beneficial for predicting affinities with the compounds of interest, and they are advantageous for guiding the design and synthesis of new and more effective anticancer agents. Graphical abstract   A new and more effective anticancer agent of xanthone derivatives against the oral human epidermoid carcinoma (KB) cell line, as investigated by CoMFA and CoMSIA analysis     

Docking-based CoMFA and CoMSIA study of azaindole carboxylic acid derivatives as promising HIV-1 integrase inhibitors

Three-dimensional quantitative structure–activity relationship (3D-QSAR) studies were performed based on a series of azaindole carboxylic acid derivatives that had previously been reported as promising HIV-1 integrase inhibitors. Docking studies to explore the binding mode were performed based on the highly active molecule 36. The best docked conformation of molecule 36 was used as template for alignment. The comparative molecular field analysis (CoMFA) model (including steric and electrostatic fields) yielded the cross validation q ² = 0.655, non-cross validation r ² = 0.989 and predictive r ² _pred = 0.979. The best comparative molecular similarity indices analysis (CoMSIA) model (including steric, electrostatic, hydrophobic and hydrogen-bond acceptor fields) yielded the cross validation q ² = 0.719, non-cross validation r ² = 0.992 and predictive r ² _pred = 0.953. A series of new azaindole carboxylic acid derivatives were designed and the HIV-1 integrase inhibitory activities of these designed compounds were predicted based on the CoMFA and CoMSIA models.     

3D-QSAR (CoMFA,CoMSIA, HQSAR and topomer CoMFA), MD simulations and molecular docking studies on purinylpyridine derivatives as B-Raf inhibitors for the treatment of melanoma cancer

As per the World Health Organization (WHO), cancer is the second most leading cause of death after cardiovascular diseases in worldwide with around 9.88 million total new cases and 1.08 million were observed due to skin cancer in 2018. Amongst two types of skin cancer, progression of melanoma cancer is increasing day by day due to the environmental changes than non-melanoma cancer. Most of B-Raf mutation, specifically B-RafV600E, is responsible for the progression of the melanoma cancer. Here, various 3D-QSAR techniques like comparative molecular field analysis (CoMFA), comparative molecular similarity indices analysis (CoMSIA), molecular hologram QSAR (HQSAR) and topomer CoMFA were used to design novel B-Raf inhibitors by using 28 synthetic B-Raf inhibitors. Except for topomer CoMFA model, remaining models were generated by three different alignment methods in which distil-based alignment method was found best and gave prominent statistical values. After performing N-fold statistical validation, in CoMFA, q², r² and r²_pred values were found to be 0.638, 0.969 and 0.848, respectively. Similarly, q², r² and r²_pred values were found to be 0.796, 0.978 and 0.891 in CoMSIA (SHD) and 0.761, 0.973 and 0.852 in CoMSIA (SH) by N-fold statistical validation. In HQSAR analysis, statistical values were found for q² as 0.984, r² as 0.999 and r²_pred as 0.634 with 97 as best hologram length (BHL). The results of topomer CoMFA showed the q² value of 0.663 and the r² value of 0.967. Important features of purinylpyridine were identified by contour map analysis of all 3D-QSAR techniques, which could be useful to design the novel molecules as B-Raf inhibitors for the treatment of melanoma cancer.

     

A combined CoMFA and CoMSIA 3D-QSAR study of benzamide type antibacterial inhibitors of the FtsZ protein in drug-resistant Staphylococcus aureus

A major problem today is bacterial resistance to antibiotics and the small number of new therapeutic agents approved in recent years. The development of new antibiotics capable of acting on new targets is urgently required. The filamenting temperature-sensitive Z (FtsZ) bacterial protein is a key biomolecule for bacterial division and survival. This makes FtsZ an attractive new pharmacological target for the development of antibacterial agents. There have been several attempts to develop ligands able to inhibit FtsZ. Despite the large number of synthesized compounds that inhibit the FtsZ protein, there are no quantitative structure–activity relationships (QSAR) that allow for the rational design and synthesis of promising new molecules. We present the first 3D-QSAR study of a large and diverse set of molecules that are able to inhibit the FtsZ bacterial protein. We summarize a set of chemical changes that can be made in the steric, electrostatic, hydrophobic and donor/acceptor hydrogen-bonding properties of the pharmacophore, to generate new bioactive molecules against FtsZ. These results provide a rational guide for the design and synthesis of promising new antibacterial agents, supported by the strong statistical parameters obtained from CoMFA (r²_pred = 0.974) and CoMSIA (r²_pred = 0.980) analyses.     

Comparative molecular field analysis and comparative molecular similarity index analysis studies on 1H NMR chemical shift of NH group of diaryl triazene derivatives

Comparative molecular field analysis (CoMFA), comparative molecular field analysis region focusing (CoMFA‐RF) for optimizing the region for the final partial least square analysis, and comparative molecular similarity indices analysis (CoMSIA) methods were employed to develop three‐dimensional quantitative structure–activity relationship (3D‐QSAR) models of ¹H NMR chemical shift of NH proton of diaryl triazene derivatives. The best orientation was searched by all‐orientation search (AOS) strategy to minimize the effect of the initial orientation of the structures. The predictive abilities of CoMFA‐RF and CoMSIA models were determined using a test set of ten compounds affording predictive correlation coefficients of 0.721 and 0.754, respectively, indicating good predictive power. For further model validation, cross validation (leave one out), progressive scrambling, and bootstrapping were also applied. The accuracy and speed of obtained 3D‐QSAR models for the prediction of ¹H NMR chemical shifts of NH group of diaryl triazene derivatives were greater compared to some computational well‐known procedures. Copyright © 2013 John Wiley & Sons, Ltd.     

3D－QSAR Study on Diindolylmethane and Its Analogues with Comparative Molecular Field Analysis(CoMFA)

Comparative molecular field analysis (CoMFA),a three dimensional quantitative structure-activity relationship (3D-QSAR) method was applied to a series of diindolylmethane(DIM) analogs to study the relationship between their structure and their induction of CYP 1A1-associated ethoxyresorufin-O-deethylase(EROD) activity.A DISCO model of pharmacophore was derved to guide the superposition of the compounds.The coefficient of cross-validation (q^2) and non cross-validation(r^2) for the model established by the study are 0.827 and 0.988 respectively,the value of variance ratio (F) is 103.53 and standard error estimate (SEE)is 0.044.These values indicate that the CoMFA model derived is significant and might have a good prediction for the catalytic activity of DIM compounds.As a consequence,the predicted activity values of new designed compounds were all higher than that of the reported value.