Ligand- and protein-based modeling studies of the inhibitors of human cytochrome P450 2D6 and a virtual screening for potential inhibitors from the Chinese herbal medicine, Scutellaria baicalensis (Huangqin,Baikal Skullcap)

We have previously examined the binding patterns of various substrates to human cytochrome P450 2D6 (CYP2D6) using a series of molecular modeling methods. In this study, we further explored the binding modes of various types of inhibitors to CYP2D6 using a combination of ligand- and protein-based modeling approaches. Firstly, we developed and validated a pharmacophore model for CYP2D6 inhibitors, which consisted of two hydrophobic features and one hydrogen bond acceptor feature. Secondly, we constructed and validated a quantitative structure-activity relationship (QSAR) model for CYP2D6 inhibitors which gave a poor to moderate prediction accuracy. Thirdly, a panel of CYP2D6 inhibitors were subject to molecular docking into the active site of wild-type and mutated CYP2D6 enzyme. We demonstrated that 8 residues in the active site (Leu213, Glu216, Ser217, Gln244, Asp301, Ser304, Ala305, and Phe483) played an important role in the binding to the inhibitors via hydrogen bond formation and/or π-π stacking interaction. Apparent changes in the binding modes of the inhibitors have been observed with Phe120Ile, Glu216Asp, Asp301Glu mutations in CYP2D6. Finally, we screened for potential binders/inhibitors from the Chinese herbal medicine Scutellaria baicalensis (Huangqin, Baikal Skullcap) using the established pharmacophore model for CYP2D6 inhibitors and molecular docking approach. Overall, 18 out of 40 compounds from S. baicalensis were mapped to the pharmacophore model of CYP2D6 inhibitors and most herbal compounds from S. baicalensis could be docked into the active site of CYP2D6. Our study has provided insights into the molecular mechanisms of interaction of synthetic and herbal compounds with human CYP2D6 and further benchmarking studies are needed to validate our modeling and virtual screening results.     

Pharmacophore based 3D-QSAR modeling,virtual screening and docking for identification of potential inhibitors of β-secretase

The enzyme β-secretase-1 is responsible for the cleavage of the amyloid precursor protein, a vital step in the process of the formation of amyloid-β peptides which are known to lead to neurodegeneration causing Alzheimer’s disease. Challenges associated with toxicity and blood brain permeation inability of potential inhibitors, continue to evade a successful therapy, thus demanding the search and development of highly active and effective inhibitors. Towards these efforts, we used a ligand based pharmacophore model generation from a dataset of known inhibitors whose activities against β-secretase hovered in the nano molar range. The identified 5 feature pharmacophore model, AHHPR, was validated via three dimensional quantitative structure activity relationship as indicated by r², q² and Pearson R values of 0.9013, 0.7726 and 0.9041 respectively. For a dataset of compounds with nano molar activity, the important pharmacophore features present in the current model appear to be similar with those observed in the models resulting from much wider activity range of inhibitors. Virtual screening of the ChemBridge CNS-Set™, a database having compounds with a better suitability for central nervous system based disorders followed by docking and analysis of the ligand protein interactions resulted in the identification of eight prospective compounds with considerable diversity. The current pharmacophore model can thus be useful for the identification, design and development of potent β-secretase inhibitors which by optimization can be potential therapeutics for Alzheimer’s disease.     

Development of predictive pharmacophore model for in silico screening,and 3D QSAR CoMFA and CoMSIA studies for lead optimization,for designing of potent tumor necrosis factor alpha converting enzyme inhibitors

A chemical feature-based pharmacophore model was developed for Tumor Necrosis Factor-α converting enzyme (TACE) inhibitors. A five point pharmacophore model having two hydrogen bond acceptors (A), one hydrogen bond donor (D) and two aromatic rings (R) with discrete geometries as pharmacophoric features was developed. The pharmacophore model so generated was then utilized for in silico screening of a database. The pharmacophore model so developed was validated by using four compounds having proven TACE inhibitory activity which were grafted into the database. These compounds mapped well onto the five listed pharmacophoric features. This validated pharmacophore model was also used for alignment of molecules in CoMFA and CoMSIA analysis. The contour maps of the CoMFA/CoMSIA models were utilized to provide structural insight for activity improvement of potential novel TACE inhibitors. The pharmacophore model so developed could be used for in silico screening of any commercial/in house database for identification of TACE inhibiting lead compounds, and the leads so identified could be optimized using the developed CoMSIA model. The present work highlights the tremendous potential of the two mutually complementary ligand-based drug designing techniques (i.e. pharmacophore mapping and 3D-QSAR analysis) using TACE inhibitors as prototype biologically active molecules.     

Structure-based virtual screening of dipeptidyl peptidase 4 inhibitors and their in vitro analysis

Type 2 diabetes mellitus (T2DM) is one of the most widely prevalent metabolic disorders with no cure to date thus remains the most challenging task in the current drug discovery. Therefore, the only strategy to control diabetes prevalence is to develop novel efficacious therapeutics. Dipeptidyl Peptidase 4 (DPP-4) inhibitors are currently used as anti-diabetic drugs for the inhibition of incretins. This study aims to construct the chemical feature based on pharmacophore models for dipeptidyl peptidase IV. The structure-based pharmacophore modeling has been employed to evaluate new inhibitors of DPP-4. A four-featured pharmacophore model was developed from crystal structure of DPP-4 enzyme with 4-(2-aminoethyl) benzenesulfonyl fluoride in its active site via pharmacophore constructing tool of Molecular Operating Environment (MOE) consisting F1 Hyd (hydrophobic region), F2 Hyd|Cat|Don (hydrophobic cationic and donor region), F3 Acc (acceptor region) and F4 Hyd (hydrophobic region). The generated pharmacophore model was used for virtual screening of in-house compound library (the available compounds which were used for initial screening to get the few compounds for the current studies). The resultant selected compounds, after virtual screening were further validated using in vitro assay. Furthermore, structure-activity relationship was carried out for the compounds possessing significant inhibition potential after docking studies. The binding free energy of analogs was evaluated via molecular mechanics generalized Born surface area (MM-GBSA) and Poisson-Boltzmann surface area (MM-PBSA) methods using AMBER 16 as a molecular dynamics (MD) simulation package. Based on potential findings, we report that selected candidates are more likely to be used as DPP-4 inhibitors or as starting leads for the development of novel and potent DPP-4 inhibitors.     

Virtual screening studies to design potent CDK2-cyclin A inhibitors

The cell division cycle is controlled by cyclin-dependent kinases (CDK), which consist of a catalytic subunit (CDK1-CDK8) and a regulatory subunit (cyclin A-H). Pharmacophore analysis indicates that the best inhibitor model consists of (1) two hydrogen bond acceptors, (2) one hydrogen bond donor, and (3) one hydrophobic feature. The HypoRefine pharmacophore model gave an enrichment factor of 1.31 and goodness of fit score of 0.76. Docking studies were carried out to explore the structural requirements for the CDK2-cyclin A inhibitors and to construct highly predictive models for the design of new inhibitors. Docking studies demonstrate the important role of hydrogen bond and hydrophobic interactions in determining the inhibitor-receptor binding affinity. The validated pharmacophore model is further used for retrieving the most active hits/lead from a virtual library of molecules. Subsequently, docking studies were performed on the hits, and novel series of potent leads were suggested based on the interaction energy between CDK2-cyclin A and the putative inhibitors.     

New leads for selective GSK-3 inhibition: pharmacophore mapping and virtual screening studies

Summary Glycogen Synthase Kinase-3 is a regulatory serine/threonine kinase, which is being targeted for the treatment of a number of human diseases including type-2 diabetes mellitus, neurodegenerative diseases, cancer and chronic inflammation. Selective GSK-3 inhibition is an important requirement owing to the possibility of side effects arising from other kinases. A pharmacophore mapping strategy is employed in this work to identify new leads for selective GSK-3 inhibition. Ligands known to show selective GSK-3 inhibition were employed in generating a pharmacophore map using distance comparison method (DISCO). The derived pharmacophore map was validated using (i) important interactions involved in selective GSK-3 inhibitions, and (ii) an in-house database containing different classes of GSK-3 selective, non-selective and inactive molecules. New Lead identification was carried out by performing virtual screening using validated pharmacophoric query and three chemical databases namely NCI, Maybridge and Leadquest. Further data reduction was carried out by employing virtual filters based on (i) Lipinski’s rule of 5 (ii) van der Waals bumps and (iii) restricting the number of rotatable bonds to seven. Final screening was carried out using FlexX based molecular docking study.     

Identification of potential HIV-1 integrase strand transfer inhibitors: In silico virtual screening and QM/MM docking studies

HIV-1 integrase (IN) is a retroviral enzyme that catalyses integration of the reverse-transcribed viral DNA into the host genome, which is necessary for efficient viral replication. In this study, we have performed an in silico virtual screening for the identification of potential HIV-1 IN strand transfer (ST) inhibitors. Pharmacophore modelling and atom-based 3D-QSAR studies were carried out for a series of compounds belonging to 3-Hydroxypyrimidine-2,4-diones. Based on the ligand-based pharmacophore model, we obtained a five-point pharmacophore with two hydrogen bond acceptors (A), one hydrogen bond donor (D), one hydrophobic group (H) and one aromatic ring (R) as pharmacophoric features. The pharmacophore hypothesis AADHR was used as a 3D query in a sequential virtual screening study to filter small molecule databases Maybridge, ChemBridge and Asinex. Hits matching with pharmacophore hypothesis AADHR were retrieved and passed progressively through Lipinski’s rule of five filtering, molecular docking and hierarchical clustering. The five compounds with best hits with novel and diverse chemotypes were subjected to QM/MM docking, which showed improved docking accuracy. We further performed molecular dynamics simulation and found three compounds that form stable interactions with key residues. These compounds could be used as a leads for further drug development and rational design of HIV-1 IN inhibitors.     

Hierarchical virtual screening: identification of potential high-affinity and selective β(3)-adrenergic receptor agonists

The hierarchical virtual screening (HVS) study, consisting of pharmacophore modelling, docking and VS of the generated focussed virtual library, has been carried out to identify novel high-affinity and selective β(3)-adrenergic receptor (β-AR) agonists. The best pharmacophore model, comprising one H-bond donor, two hydrophobes, one positive ionizable and one negative ionizable feature, was developed based on a training set of 51 β(3)-AR agonists using the pharmacophore generation protocol implemented in Discovery Studio. The model was further validated with the test set, external set and ability of the pharmacophoric features to complement the active site amino acids of the homology modelled β(3)-AR developed using MODELLER software. The focussed virtual library was generated using the structure-based insights gained from our earlier reported comprehensive study focussing on the structural basis of β-AR subtype selectivity of representative agonists and antagonists. The HVS with the sequential use of the best pharmacophore model and homology modelled β(3)-AR in the screening of the generated focussed library has led to the identification of potential virtual leads as novel high-affinity and selective β(3)-AR agonists.     

人类腺苷受体A3亚型拮抗剂的构效关系分析

构建人类腺苷受体A3亚型药效团模型和三维蛋白结构模型用于作用模式研究.以18个来源于文献具有腺苷受体A3亚型拮抗活性的化合物作为训练集,使用HypoGen方法构建药效团模型.通过同源模建和分子动力学模拟构建了人类腺苷受体A3亚型的三维蛋白模型,并利用PROCHECK方法评估该模型的合理性,对所得的结构使用分子对接程序进行作用模式分析,药效团模型和同源模建结果相互匹配较好.使用新药效团模型对MDL药物数据库(MDDR)中包含的约120000个化合物进行虚拟筛选,得到了8个候选化合物,用于进一步的生物学评价和活性测定.本工作对于人类腺苷受体A3亚型拮抗剂的设计和抗哮喘药物的研发具有一定的理论指导和应用价值.     

A virtual screening approach for thymidine monophosphate kinase inhibitors as antitubercular agents based on docking and pharmacophore models

Docking and pharmacophore screening tools were used to examine the binding of ligands in the active site of thymidine monophosphate kinase of Mycobacterium tuberculosis. Docking analysis of deoxythymidine monophosphate (dTMP) analogues suggests the role of hydrogen bonding and other weak interactions in enzyme selectivity. Water-mediated hydrogen-bond networks and a halogen-bond interaction seem to stabilize the molecular recognition. A pharmacophore model was developed using 20 dTMP analogues. The pharmacophoric features were complementary to the active site residues involved in the ligand recognition. On the basis of these studies, a composite screening model that combines the features from both the docking analysis and the pharmacophore model was developed. The composite model was validated by screening a database spiked with 47 known inhibitors. The model picked up 42 of these, giving an enrichment factor of 17. The validated model was used to successfully screen an in-house database of about 500,000 compounds. Subsequent screening with other filters gave 186 hit molecules.     

On the value of homology models for virtual screening: discovering hCXCR3 antagonists by pharmacophore-based and structure-based approaches

Human chemokine receptor CXCR3 (hCXCR3) antagonists have potential therapeutic applications as antivirus, antitumor, and anti-inflammatory agents. A novel virtual screening protocol, which combines pharmacophore-based and structure-based approaches, was proposed. A three-dimensional QSAR pharmacophore model and a structure-based docking model were built to virtually screen for hCXCR3 antagonists. The hCXCR3 antagonist binding site was constructed by homology modeling and molecular dynamics (MD) simulation. By combining the structure-based and ligand-based screenings results, 95% of the compounds satisfied either pharmacophore or docking score criteria and would be chosen as hits if the union of the two searches was taken. The false negative rates were 15% for the pharmacophore model, 14% for the homology model, and 5% for the combined model. Therefore, the consistency of the pharmacophore model and the structural binding model is 219/273 = 80%. The hit rate for the virtual screening protocol is 273/286 = 95%. This work demonstrated that the quality of both the pharmacophore model and homology model can be measured by the consistency of the two models, and the false negatives in virtual screening can be reduced by combining two virtual screening approaches.     

CoMFA 3D-QSAR analysis of HIV-1 RT nonnucleoside inhibitors, TIBO derivatives based on docking conformation and alignment

HIV-1 RT is one of the key enzymes in the duplication of HIV-1. Inhibitors of HIV-1 RT are classified as nonnucleoside RT inhibitors (NNRTIs) and nucleoside analogues. NNRTIs bind in a region not associated with the active site of the enzyme. Within the NNRTI category, there is a set of inhibitors commonly referred to as TIBO inhibitors. Fifty TIBO inhibitors were used in the work to build 3-D QSAR models. The two known crystal structures of complexes are used to investigate and validate the docking protocol. The results show that the docking simulations reproduce the crystal complexes very well with RMSDs of approximately 1 A and approximately 0.6 A for 1REV and 1COU, respectively. The alignment of molecules and "active" conformation selection are the key to a successful 3D-QSAR model by CoMFA. The flexible docking (Autodock3) was used on determination of "active" conformation and molecular alignment, and CoMFA and CoMSIA were used to develop 3D-QSAR models of 50 TIBOs in the work. The 3D-QSAR models demonstrate a good ability to predict the activity of studied compounds (r2 = 0.972, 0.944, q2 = 0.704, 0.776). It is shown that the steric and electrostatic properties predicted by CoMFA contours can be related to the binding structure of the complex. The results demonstrate that the combination of ligand-based and receptor-based modeling is a powerful approach to build 3D-QSAR models.     

Docking-based comparative intermolecular contacts analysis as new 3-D QSAR concept for validating docking studies and in silico screening: NMT and GP inhibitors as case studies

The significant role played by docking algorithms in drug discovery combined with their serious pitfalls prompted us to envisage a novel concept for validating docking solutions, namely, docking-based comparative intermolecular contacts analysis (dbCICA). This novel approach is based on the number and quality of contacts between docked ligands and amino acid residues within the binding pocket. It assesses a particular docking configuration on the basis of its ability to align a set of ligands within a corresponding binding pocket in such a way that potent ligands come into contact with binding site spots distinct from those approached by low-affinity ligands and vice versa. In other words, dbCICA evaluates the consistency of docking by assessing the correlation between ligands' affinities and their contacts with binding site spots. Optimal dbCICA models can be translated into valid pharmacophore models that can be used as 3-D search queries to mine structural databases for new bioactive compounds. dbCICA was implemented to search for new inhibitors of candida N-myristoyl transferase as potential antifungal agents and glycogen phosphorylase (GP) inhibitors as potential antidiabetic agents. The process culminated in five selective micromolar antifungal leads and nine GP inhibitory leads.     

Identification of Influenza PAN Endonuclease Inhibitors via 3D-QSAR Modeling and Docking-Based Virtual Screening

Structural and biochemical studies elucidate that PA_N may contribute to the host protein shutdown observed during influenza A infection. Thus, inhibition of the endonuclease activity of viral RdRP is an attractive approach for novel antiviral therapy. In order to envisage structurally diverse novel compounds with better efficacy as PA_N endonuclease inhibitors, a ligand-based-pharmacophore model was developed using 3D-QSAR pharmacophore generation (HypoGen algorithm) methodology in Discovery Studio. As the training set, 25 compounds were taken to generate a significant pharmacophore model. The selected pharmacophore Hypo1 was further validated by 12 compounds in the test set and was used as a query model for further screening of 1916 compounds containing 71 HIV-1 integrase inhibitors, 37 antibacterial inhibitors, 131 antiviral inhibitors and other 1677 approved drugs by the FDA. Then, six compounds (Hit01–Hit06) with estimated activity values less than 10 μM were subjected to ADMET study and toxicity assessment. Only one potential inhibitory ‘hit’ molecule (Hit01, raltegravir’s derivative) was further scrutinized by molecular docking analysis on the active site of PA_N endonuclease (PDB ID: 6E6W). Hit01 was utilized for designing novel potential PA_N endonuclease inhibitors through lead optimization, and then compounds were screened by pharmacophore Hypo1 and docking studies. Six raltegravir’s derivatives with significant estimated activity values and docking scores were obtained. Further, these results certainly do not confirm or indicate the seven compounds (Hit01, Hit07, Hit08, Hit09, Hit10, Hit11 and Hit12) have antiviral activity, and extensive wet-laboratory experimentation is needed to transmute these compounds into clinical drugs.     

Computer-aided molecular design of (E)-N-Aryl-2-ethene-sulfonamide analogues as microtubule targeted agents in prostate cancer

Microtubules are tube-shaped, filamentous and cytoskeletal proteins that are essential in all eukaryotic cells. Microtubule is an attractive and promising target for anticancer agents. In this study, three-dimensional quantitative structure activity relationships (3D-QSAR) including comparative molecular field analysis, CoMFA, and comparative molecular similarity indices analysis, CoMSIA, were performed on a set of 45 (E)-N-Aryl-2-ethene-sulfonamide analogues as microtubule-targeted anti-prostate cancer agents. Automated grid potential analysis, AutoGPA module in Molecular Operating Environment 2009.10 (MOE) as a new 3D-QSAR approach with the pharmacophore-based alignment was carried out on the same dataset. AutoGPA-based 3D-QSAR model yielded better prediction parameters than CoMFA and CoMSIA. Based on the contour maps generated from the models, some key features were identified in (E)-N-Aryl-2-arylethene-sulfonamide analogues that were responsible for the anti-cancer activity. Virtual screening was performed based on pharmacophore modeling and molecular docking to identify the new inhibitors from ZINC database. Seven top ranked compounds were found based on Gold score fitness function. In silico ADMET studies were performed on compounds retrieved from virtual screening in compliance with the standard ranges.     

Sphingosine kinase 1 (SK1) allosteric inhibitors that target the dimerization site

The sphingosine kinase 1 (SK1)/sphingosine-1-phosphate (S1P) signaling pathway is a crucial target for numerous human diseases from cancer to cardiovascular diseases. However, available SK1 inhibitors that target the active site suffer from poor potency, selectivity and pharmacokinetic properties. The selectivity issue of the kinases, which share a highly-conserved ATP-pocket, can be overcome by targeting the less-conserved allosteric sites. SK1 is known to function minimally as a dimer; however, the crystal structure of the SK1 dimer has not been determined. In this study, a template-based algorithm implemented in PRISM was used to predict the SK1 dimer structure and then the possible allosteric sites at the dimer interface were determined via SiteMap. These sites were used in a virtual screening campaign that includes an integrated workflow of structure-based pharmacophore modeling, virtual screening, molecular docking, re-screening of common scaffolds to propose a series of compounds with different scaffolds as potential allosteric SK1 inhibitors. Finally, the stability of the SK1-ligand complexes was analyzed by molecular dynamics simulations. As a final outcome, ligand 7 having a 4,9-dihydro-1H-purine scaffold and ligand 12 having a 2,3,4,9-tetrahydro-1H-β-carboline scaffold were found to be potential selective inhibitors for SK1.     

Site-Identification by Ligand Competitive Saturation (SILCS) assisted pharmacophore modeling

Database screening using receptor-based pharmacophores is a computer-aided drug design technique that uses the structure of the target molecule (i.e. protein) to identify novel ligands that may bind to the target. Typically receptor-based pharmacophore modeling methods only consider a single or limited number of receptor conformations and map out the favorable binding patterns in vacuum or with a limited representation of the aqueous solvent environment, such that they may suffer from neglect of protein flexibility and desolvation effects. Site-Identification by Ligand Competitive Saturation (SILCS) is an approach that takes into account these, as well as other, properties to determine 3-dimensional maps of the functional group-binding patterns on a target receptor (i.e. FragMaps). In this study, a method to use the FragMaps to automatically generate receptor-based pharmacophore models is presented. It converts the FragMaps into SILCS pharmacophore features including aromatic, aliphatic, hydrogen-bond donor and acceptor chemical functionalities. The method generates multiple pharmacophore hypotheses that are then quantitatively ranked using SILCS grid free energies. The pharmacophore model generation protocol is validated using three different protein targets, including using the resulting models in virtual screening. Improved performance and efficiency of the SILCS derived pharmacophore models as compared to published docking studies, as well as a recently developed receptor-based pharmacophore modeling method is shown, indicating the potential utility of the approach in rational drug design.