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.     

Combined 2D and 3D-QSAR,molecular modelling and docking studies of pyrazolodiazepinones as novel phosphodiesterase 2 inhibitors

Selective inhibition of phosphodiesterase 2 (PDE2) in cells where it is located elevates cyclic guanosine monophosphate (cGMP) and acts as novel analgesic with antinociceptive activity. Three-dimensional quantitative structure–activity relationship (QSAR) studies for pyrazolodiazepinone inhibitors exhibiting PDE2 inhibition were performed using comparative molecular field analysis (CoMFA), comparative molecular similarity indices analysis (CoMSIA) and Topomer CoMFA, and two-dimensional QSAR study was performed using a Hologram QSAR (HQSAR) method. QSAR models were generated using training set of 23 compounds and were validated using test set of nine compounds. The optimum partial least squares (PLS) for CoMFA-Focusing, CoMSIA-SDH, Topomer CoMFA and HQSAR models exhibited good ‘leave-one-out’ cross validated correlation coefficient (q²) of 0.790, 0.769, 0.840 and 0.787, coefficient of determination (r²) of 0.999, 0.964, 0.979 and 0.980, and high predictive power (r²_pred) of 0.796, 0.833, 0.820 and 0.803 respectively. Docking studies revealed that those inhibitors able to bind to amino acid Gln859 by cGMP binding orientation called ‘glutamine-switch’, and also bind to the hydrophobic clamp of PDE2 isoform, could possess high selectivity for PDE2. From the results of all the studies, structure–activity relationships and structural requirements for binding to active site of PDE2 were established which provide useful guidance for the design and future synthesis of potent PDE2 inhibitors.     

Molecular docking and QSAR of aplyronine A and analogues: potent inhibitors of actin

Actin-binding natural products have been identified as a potential basis for the design of cancer therapeutic agents. We report flexible docking and QSAR studies on aplyronine A analogues. Our findings show the macrolide ‘tail’ to be fundamental for the depolymerisation effect of actin-binding macrolides and that it is the tail which forms the initial interaction with the actin rather than the macrocycle, as previously believed. Docking energy scores for the compounds were highly correlated with actin depolymerisation activity. The 3D-QSAR models were predictive, with the best model giving a q ² value of 0.85 and a r ² of 0.94. Results from the docking simulations and the interpretation from QSAR “coeff*stdev” contour maps provide insight into the binding mechanism of each analogue and highlight key features that influence depolymerisation activity. The results herein may aid the design of a putative set of analogues that can help produce efficacious and tolerable anti-tumour agents. Finally, using the best QSAR model, we have also made genuine predictions for an independent set of recently reported aplyronine analogues.     

Structure-based QSAR analysis of a set of 4-hydroxy-5,6-dihydropyrones as inhibitors of HIV-1 protease: an application of the receptor-dependent (RD) 4D-QSAR formalism

Receptor-dependent (RD) 4D-QSAR models were constructed for a set of 39 4-hydroxy-5,6-dihydropyrone analogue HIV-1 protease inhibitors. The receptor model used in this QSAR analysis was derived from the HIV-1 protease (PDB ID ) crystal structure. The bound ligand in the active site of the enzyme, also a 4-hydroxy-5,6-dihydropyrone analogue, was used as the reference ligand for docking the data set compounds. The optimized RD 4D-QSAR models are not only statistically significant (r(2) = 0.86, q(2) = 0.80 for four- and greater-term models) but also possess reasonable predictivity based on test set predictions. The proposed "active" conformations of the docked analogues in the active site of the enzyme are consistent in overall molecular shape with those suggested from crystallographic studies. Moreover, the RD 4D-QSAR models also "capture" the existence of specific induced-fit interactions between the enzyme active site and each specific inhibitor. Hydrophobic interactions, steric shape requirements, and hydrogen bonding of the 4-hydroxy-5,6-dihydropyrone analogues with the HIV-1 protease binding site model dominate the RD 4D-QSAR models in a manner again consistent with experimental conclusions. Some possible hypotheses for the development of new lead HIV-1 protease inhibitors can be inferred from the RD 4D-QSAR models.     

Identification of novel nt-MGAM inhibitors for potential treatment of type 2 diabetes: Virtual screening,atom based 3D-QSAR model,docking analysis and ADME study

In this study, a virtual screening procedure was applied to identify new potential nt-MGAM inhibitors as a possible medication for type 2 diabetes. To this aim, a series of salacinol analogues were first investigated by docking analysis for their binding to the X-ray structure of the biological target nt-MGAM. Key interactions for ligand binding into the receptor active site were identified which shared common features to those found for other known inhibitors, which strengthen the results of this study. 3D QSAR model was then built and showed to be statistically significant and with a good predictive power for the training (R² = 0.99, SD = 0.17, F = 555.3 and N = 27) and test set (Q² = 0.81, Pearson(r) = 0.92, RMSE = 0.52, N = 08). The model was then used to virtually screen the ZINC database with the aim of identifying novel chemical scaffolds as potential nt-MGAM inhibitors. Further, in silico predicted ADME properties were investigated for the most promising molecules. The outcome of this investigation sheds light on the molecular characteristics of the binding of salacinol analogues to nt-MGAM enzyme and identifies new possible inhibitors which have the potential to be developed into drugs, thus significantly contributing to the design and optimization of therapeutic strategies against type 2 diabetes.     

Fragment-similarity-based QSAR (FS-QSAR) algorithm for ligand biological activity predictions

Quantitative structure–activity relationship (QSAR) studies are useful computational tools often used in drug discovery research and in many scientific disciplines. In this study, a robust fragment-similarity-based QSAR (FS-QSAR) algorithm was developed to correlate structures with biological activities by integrating fragment-based drug design concept and a multiple linear regression method. Similarity between any pair of training and testing fragments was determined by calculating the difference of lowest or highest eigenvalues of the chemistry space BCUT matrices of corresponding fragments. In addition to the BCUT-similarity function, molecular fingerprint Tanimoto coefficient (Tc) similarity function was also used as an alternative for comparison. For validation studies, the FS-QSAR algorithm was applied to several case studies, including a dataset of COX2 inhibitors and a dataset of cannabinoid CB2 triaryl bis-sulfone antagonist analogues, to build predictive models achieving average coefficient of determination (r ²) of 0.62 and 0.68, respectively. The developed FS-QSAR method is proved to give more accurate predictions than the traditional and one-nearest-neighbour QSAR methods and can be a useful tool in the fragment-based drug discovery for ligand activity prediction.     

Monte Carlo method based QSAR modelling of natural lipase inhibitors using hybrid optimal descriptors

Obesity is one of the most provoking health burdens in the developed countries. One of the strategies to prevent obesity is the inhibition of pancreatic lipase enzyme. The aim of this study was to build QSAR models for natural lipase inhibitors by using the Monte Carlo method. The molecular structures were represented by the simplified molecular input line entry system (SMILES) notation and molecular graphs. Three sets – training, calibration and test set of three splits – were examined and validated. Statistical quality of all the described models was very good. The best QSAR model showed the following statistical parameters: r² = 0.864 and Q² = 0.836 for the test set and r² = 0.824 and Q² = 0.819 for the validation set. Structural attributes for increasing and decreasing the activity (expressed as pIC₅₀) were also defined. Using defined structural attributes, the design of new potential lipase inhibitors is also presented. Additionally, a molecular docking study was performed for the determination of binding modes of designed molecules.     

Development of docking-based 3D QSAR models for the design of substituted quinoline derivatives as human dihydroorotate dehydrogenase (hDHODH) inhibitors

Abstract

This study has investigated docking-based 3D quantitative structure–activity relationships (QSARs) for a range of quinoline carboxylic acid derivatives by comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA). A docking study has shown that most of the compounds formed H-bonds with Arg136 and Gln47, which have already been shown to be essential for the binding of ligands at the active site of the hydroorotate dehydrogenase adenovirus (hDHODH). Bioactive conformations of all the molecules obtained from the docking study were used for the 3D QSAR study. The best CoMFA and CoMSIA models were obtained for the training set and were found to be statistically significant, with cross-validated coefficients (q² ) of 0.672 and 0.613, r² _cv of 0.635 and 0.598 and coefficients of determination (r² ) of 0.963 and 0.896, respectively. Both models were validated by a test set of 15 compounds, giving satisfactory predicted correlation coefficients (r² _pred) of 0.824 and 0.793 for the CoMFA and CoMSIA models, respectively. From the docking-based 3D QSAR study we designed 34 novel quinoline-based compounds and performed structure-based virtual screening. Finally, in silico pharmacokinetics and toxicities were predicted for 24 of the best docked molecules. The study provides valuable information for the understanding of interactions between hDHODH and the novel compounds.     

Design,synthesis and computational approach to study novel pyrrole scaffolds as active inhibitors of enoyl ACP reductase (InhA) and Mycobacterium tuberculosis antagonists

Novel 54 pyrrolyl acetohydrazide analogues were designed, synthesized and screened for antitubercular activity against InhA. Enoyl-ACP reductase/InhA one of the significant enzymes implicated in type II FAS (fatty acid synthase) biosynthetic pathway of bacterial outer cell membrane, in addition Mycobacterium tuberculosis H37Rv inhibition potency has proven to be one of the most promising drug target used for designing and testing against TB. In silico molecular modeling was achieved using Surflex-docking method to recognize important binding sites of the enoyl-ACP reductase. 3D-QSAR studies like CoMFA and CoMSIA approaches were studied to create 3D-QSAR depictions for InhA inhibitors. Based on docking results the synthesized molecules are oriented towards the core of the active site. The pyrolyl acetohydrazides exhibited one or two H-bonding connections with InhA enzyme. Molecule 8l (MIC 0.4 μg/mL, InhA- 70% at 50 μM) showed H-bonding connections with Tyr158 and NAD⁺ in a similar mode to that of ligand pyrrolidine carboxamide. The tested molecules also showed good antibacterial (Escherichia coli, Staphylococcus aureus) activity (MIC 0.4–25 μg/mL), while the study against A549 lung adenocarcinoma cell line confirms nontoxic nature of the reported molecules. The QSAR model from CoMFA and CoMSIA through the database configuration showed the most excellent data. The prognostic ability of CoMFA and CoMSIA representations was developed by a test set of 15 molecules that produced cross validated correlation coefficients (q²) of 0.642 and 0.701, respectively. This study gives insight into structural requirement needed for the development of more active InhA inhibitors through in silico approach.     

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.     

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.