Identification of novel, less toxic PTP-LAR inhibitors using in silico strategies: pharmacophore modeling, SADMET-based virtual screening and docking

Human leukocyte antigen-related (PTP-LAR) is a receptor-like transmembrane phosphatase and a potential target for diabetes, obesity and cancer. In the present study, a sequence of in silico strategies (pharmacophore mapping, a 3D database searching, SADMET screening, and docking and toxicity studies) was performed to identify eight novel nontoxic PTP-LAR inhibitors. Twenty different pharmacophore hypotheses were generated using two methods; the best (hypothesis 2) consisted of three hydrogen-bond acceptor (A), one ring aromatic (R), and one hydrophobic aliphatic (Z) features. This hypothesis was used to screen molecules from several databases, such as Specs, IBS, MiniMaybridge, NCI, and an in-house PTP inhibitor database. In order to overcome the general bioavailability problem associated with phosphatases, the hits obtained were filtered by Lipinski’s rule of five and SADMET properties and validated by molecular docking studies using the available crystal structure 1LAR. These docking studies suggested the ligand binding pattern and interactions required for LAR inhibition. The docking analysis also revealed that sulfonylurea derivatives with an isoquinoline or naphthalene scaffold represent potential LAR drugs. The screening protocol was further validated using ligand pharmacophore mapping studies, which showed that the abovementioned interactions are indeed crucial and that the screened molecules can be presumed to possess potent inhibitory activities.     

Pharmacophore modeling,molecular docking and molecular dynamics studies on natural products database to discover novel skeleton as non-purine xanthine oxidase inhibitors

Gout is a common inflammatory arthritis caused by the deposition of urate crystals within joints. It is increasingly in prevalence during the past few decades as shown by the epidemiological survey results. Xanthine oxidase (XO) is a key enzyme to transfer hypoxanthine and xanthine to uric acid, whose overproduction leads to gout. Therefore, inhibiting the activity of xanthine oxidase is an important way to reduce the production of urate. In the study, in order to identify the potential natural products targeting XO, pharmacophore modeling was employed to filter databases. Here, two methods, pharmacophore based on ligand and pharmacophore based on receptor-ligand, were constructed by Discovery Studio. Then GOLD was used to refine the potential compounds with higher fitness scores. Finally, molecular docking and dynamics simulations were employed to analyze the interactions between compounds and protein. The best hypothesis was set as a 3D query to screen database, returning 785 and 297 compounds respectively. A merged set of the above 1082 molecules was subjected to molecular docking, which returned 144 hits with high-fitness scores. These molecules were clustered in four main kinds depending on different backbones. What is more, molecular docking showed that the representative compounds established key interactions with the amino acid residues in the protein, and the RMSD and RMSF of molecular dynamics results showed that these compounds can stabilize the protein. The information represented in the study confirmed previous reports. And it may assist to discover and design new backbones as potential XO inhibitors based on natural products.     

Flexible ligand docking using conformational ensembles.

Molecular docking algorithms suggest possible structures for molecular complexes. They are used to model biological function and to discover potential ligands. A present challenge for docking algorithms is the treatment of molecular flexibility. Here, the rigid body program, DOCK, is modified to allow it to rapidly fit multiple conformations of ligands. Conformations of a given molecule are pre-calculated in the same frame of reference, so that each conformer shares a common rigid fragment with all other conformations. The ligand conformers are then docked together, as an ensemble, into a receptor binding site. This takes advantage of the redundancy present in differing conformers of the same molecule. The algorithm was tested using three organic ligand protein systems and two protein-protein systems. Both the bound and unbound conformations of the receptors were used. The ligand ensemble method found conformations that resembled those determined in X-ray crystal structures (RMS values typically less than 1.5 A). To test the method's usefulness for inhibitor discovery, multi-compound and multi-conformer databases were screened for compounds known to bind to dihydrofolate reductase and compounds known to bind to thymidylate synthase. In both cases, known inhibitors and substrates were identified in conformations resembling those observed experimentally. The ligand ensemble method was 100-fold faster than docking a single conformation at a time and was able to screen a database of over 34 million conformations from 117,000 molecules in one to four CPU days on a workstation.     

Discover potential inhibitors for PFKFB3 using 3D-QSAR,virtual screening,molecular docking and molecular dynamics simulation

Abstract

The 6-phosphofructo-2-kinase/fructose-2,6-bisphosphatase-3 (PFKFB3) is a master regulator of glycolysis in cancer cells by synthesizing fructose-2,6-bisphosphate (F-2,6-BP), a potent allosteric activator of phosphofructokinase-1 (PFK-1), which is a rate-limiting enzyme of glycolysis. PFKFB3 is an attractive target for cancer treatment. It is valuable to discover promising inhibitors by using 3D-QSAR pharmacophore modeling, virtual screening, molecular docking and molecular dynamics simulation. Twenty molecules with known activity were used to build 3D-QSAR pharmacophore models. The best pharmacophore model was ADHR called Hypo1, which had the highest correlation value of 0.98 and the lowest RMSD of 0.82. Then, the Hypo1 was validated by cost value method, test set method and decoy set validation method. Next, the Hypo1 combined with Lipinski's rule of five and ADMET properties were employed to screen databases including Asinex and Specs, total of 1,048,159 molecules. The hits retrieved from screening were docked into protein by different procedures including HTVS, SP and XP. Finally, nine molecules were picked out as potential PFKFB3 inhibitors. The stability of PFKFB3-lead complexes was verified by 40?ns molecular dynamics simulation. The binding free energy and the energy contribution of per residue to the binding energy were calculated by MM-PBSA based on molecular dynamics simulation.     

Novel virtual lead identification in the discovery of hematopoietic cell kinase (HCK) inhibitors: application of 3D QSAR and molecular dynamics simulation

High level of hematopoietic cell kinase (Hck) is associated with drug resistance in chronic myeloid leukemia. Additionally, Hck activity has also been connected with the pathogenesis of HIV-1 and chronic obstructive pulmonary disease. In this study, three-dimensional (3D) QSAR pharmacophore models were generated for Hck based on experimentally known inhibitors. A best pharmacophore model, Hypo1, was developed with high correlation coefficient (0.975), Low RMS deviation (0.60) and large cost difference (49.31), containing three ring aromatic and one hydrophobic aliphatic feature. It was further validated by the test set (r?=?0.96) and Fisher’s randomization method (95%). Hypo 1 was used as a 3D query for screening the chemical databases, and the hits were further screened by applying Lipinski’s rule of five and ADMET properties. Selected hit compounds were subjected to molecular docking to identify binding conformations in the active site. Finally, the appropriate binding modes of final hit compounds were revealed by molecular dynamics (MD) simulations and free energy calculation studies. Hence, we propose the final three hit compounds as virtual candidates for Hck inhibitors.     

Pharmacoinformatics-based identification of chemically active molecules against Ebola virus

Abstract

Ebola is a dangerous virus transmitted by animals and humans and to date there is no curable agent for such a deadly infectious disease. In this study, pharmacoinformatics-based methods were adopted to find effective novel chemical entities against Ebola virus. A well predictive and statistical robust pharmacophore model was developed from known Ebola virus inhibitors collected from the literature. The model explained the significance of each of hydrogen bond acceptor and donor, and two hydrophobic regions for activity. The National Cancer Institute and Asinex (Antiviral library) databases were screened using the final validated pharmacophore model. Initial hits were further screened with a set of criteria and finally eight molecules from both databases were proposed as promising anti Ebola agents. Further molecular docking and molecular dynamics studies were carried out and it was found that the proposed molecules possessed capability to interact with amino residues of Ebola protein as well as retaining equilibrium of protein-ligand systems. Finally, the binding energies were calculated using molecular mechanics Poisson–Boltzmann surface area approach and all proposed molecules showed strong binding affinity towards the Ebola protein receptor.

Communicated by Ramaswamy H. Sarma     

3D-QSAR pharmacophore modelling,virtual screening and docking studies for lead discovery of a novel scaffold for VEGFR 2 inhibitors: Design,synthesis and biological evaluation

A series of novel 6,7-dihydro-5H-cyclopenta[d]pyrimidine derivatives was successfully designed, synthesized and evaluated as a new chemical scaffold with vascular endothelial growth factor receptor (VEGFR 2) inhibitory activity. Compounds 6c and 6b showed enzyme inhibition of 97% and 87% at 10 µM, respectively, and exhibited potent dose-related VEGFR 2 inhibition with IC₅₀ values of 0.85 µM and 2.26 µM, respectively. The design of the 6,7-dihydro-5H-cyclopenta[d]pyrimidine scaffold was implemented via consecutive molecular modelling protocols prior to the synthesis and biological evaluation of the derivatives. First, sorafenib was docked in the binding site of VEGFR 2 to study its binding orientation and affinity, followed by the generation of a valid 3D QSAR pharmacophore model for use in the virtual screening of different 3D databases. Structures with promising pharmacophore-based virtual screening results were refined using molecular docking studies in the binding site of VEGFR 2. A novel scaffold was designed by incorporating the results of the pharmacophore model generation and molecular docking studies. The new scaffold showed hydrophobic interactions with the kinase front pocket that may be attributed to increasing residence time in VEGFR 2, which is a key success factor for ligand optimization in drug discovery. Different derivatives of the novel scaffold were validated using docking studies and pharmacophore mapping, where they exhibited promising results as VEGFR 2 inhibitors to be synthesized and biologically evaluated. 6,7-dihydro-5H-cyclopenta[d]pyrimidine is a new scaffold that can be further optimized for the synthesis of promising VEGFR 2 inhibitors.     

Discovery of promising FtsZ inhibitors by E-pharmacophore, 3D-QSAR,molecular docking study,and molecular dynamics simulation

Filamentous temperature-sensitive protein Z (FtsZ), playing a key role in bacterial cell division, is regarded as a promising target for the design of antimicrobial agent. This study is looking for potential high-efficiency FtsZ inhibitors. Ligand-based pharmacophore and E-pharmacophore, virtual screening and molecular docking were used to detect promising FtsZ inhibitors, and molecular dynamics simulation was used to study the stability of protein-ligand complexes in this paper. Sixty-three inhibitors from published literatures with pIC₅₀ ranging from 2.483 to 5.678 were collected to develop ligand-based pharmacophore model. 4DXD bound with 9PC was selected to develop the E-pharmacophore model. The pharmacophore models validated by test set method and decoy set were employed for virtual screening to exclude inactive compounds against ZINC database. After molecular docking, ADME analysis, IFD docking and MM-GBSA, 8 hits were identified as potent FtsZ inhibitors. A 50?ns molecular dynamics simulation was implemented on the compounds to assess the stability between potent inhibitors and FtsZ. The results indicated that the candidate compounds had a high docking score and were strongly combined with FtsZ by forming hydrogen bonding interactions with key amino acid residues, and van der Waals forces and hydrophobic interactions had significant contribution to the stability of the binding. Molecular dynamics simulation results showed that the protein-ligand compounds performed well in both the stability and flexibility of the simulation process.     

Identification of novel inhibitor of protein tyrosine phosphatases delta: structure-based pharmacophore modeling,virtual screening,flexible docking,molecular dynamics simulation,and post-molecular dynamics analysis

Abstract

Owing to their unique functions in regulating the synapse activity of protein tyrosine phosphatases delta (PTPδ) that has drawn special attention for developing drugs to autism spectrum disorders (ASDs). In this study, the PTPδ pharmacophore was first established by the structure-based pharmacophore method. Subsequently, 10 compounds contented Lipinski’s rule of five was acquired by the virtual screening of the PTPδ pharmacophore against ZINC and PubChem databases. Then, the 10 identified molecules were discovered that had better binding affinity than a known PTPδ inhibitors compound SCHEMBL16375396. Two compounds SCHEMBL16375408 and ZINC19796658 with high binding score, low toxicity were gained. They were observed by docking analysis and molecular dynamics simulations that the novel potential inhibitors not only possessed the same function as SCHEMBL16375396 did in inhibiting PTPδ, but also had more favorable conformation to bind with the catalytic active regions. This study provides a new method for identify PTPδ inhibitor for the treatment of ASDs disease.

Communicated by Ramaswamy H. Sarma     

Modeling correlated main-chain motions in proteins for flexible molecular recognition

We describe a new method for modeling protein and ligand main-chain flexibility, and show its ability to model flexible molecular recognition. The goal is to sample the full conformational space, including large-scale motions that typically cannot be reached in molecular dynamics simulations due to the computational intensity, as well as conformations that have not been observed yet by crystallography or NMR. A secondary goal is to assess the degree of flexibility consistent with protein-ligand recognition. Flexibility analysis of the target protein is performed using the graph-theoretic algorithm FIRST, which also identifies coupled networks of covalent and noncovalent bonds within the protein. The available conformations of the flexible regions are then explored with ROCK by random-walk sampling of the rotatable bonds. ROCK explores correlated motions by only sampling dihedral angles that preserve the coupled bond networks in the protein and generates conformers with good stereochemistry, without using a computationally expensive potential function. A representative set of the conformational ensemble generated this way can be used as targets for docking with SLIDE, which handles the flexibility of protein and ligand side-chains. The realism of this protein main-chain conformational sampling is assessed by comparison with time-resolved NMR studies of cyclophilin A motions. ROCK is also effective for modeling the flexibility of large cyclic and polycyclic ligands, as demonstrated for cyclosporin and zearalenol. The use of this combined approach to perform docking with main-chain flexibility is illustrated for the cyclophilin A-cyclosporin complex and the estrogen receptor in complex with zearalenol, while addressing the question of how much flexibility is allowed without hindering molecular recognition.     

Pharmacophore models generation by catalyst and phase consensus-based virtual screening protocol against PI3Kα inhibitors

Phosphoinositide 3-kinases (PI3Ks) family has emerged as promising targets for novel therapeutic agents against neoplastic diseases. Pharmacophore and 3D-quantitative structure–activity relationship modelling were applied to study the structure–activity relationship of PI3K inhibitors. The best HypoGen pharmacophore hypothesis Hypo1 with a correlation coefficient of 0.961 consists of one hydrogen-bond acceptor, one hydrogen-bond donor and two hydrophobic features, whereas the best phase hypothesis AADRRR.378 with favourable statistics (q² = 0.7368, r² = 0.9863) has two hydrogen-bond acceptors, one hydrogen-bond donor and three ring aromatic features. Multiple methods, such as Fischer validation, molecular docking and mapping of test set molecules, were carried out to validate these pharmacophore models. Furthermore, a comparative molecular similarity indices analysis candidate hypothesis model was generated as a supplement of pharmacophore hypothesis. Detailed protein–ligand binding information obtained by Glide was utilised in compound optimisation and virtual screening. A molecular database of 133 known inhibitors and 6179 decoys was built for a screening test to quantitatively analyse various hypotheses and scoring parameters. Finally, we designed a workflow integrating HypoGen pharmacophore searching, phase pharmacophore searching and molecular docking for screening the database. With an improved criterion of enrichment factor (EF = 17.43) and ROC curve (AUC = 0.946), this workflow would provide us an original method for novel PI3K inhibitors.     

A hierarchical method for molecular docking using cloud computing

Discovering small molecules that interact with protein targets will be a key part of future drug discovery efforts. Molecular docking of drug-like molecules is likely to be valuable in this field; however, the great number of such molecules makes the potential size of this task enormous. In this paper, a method to screen small molecular databases using cloud computing is proposed. This method is called the hierarchical method for molecular docking and can be completed in a relatively short period of time. In this method, the optimization of molecular docking is divided into two subproblems based on the different effects on the protein–ligand interaction energy. An adaptive genetic algorithm is developed to solve the optimization problem and a new docking program (FlexGAsDock) based on the hierarchical docking method has been developed. The implementation of docking on a cloud computing platform is then discussed. The docking results show that this method can be conveniently used for the efficient molecular design of drugs.     

Pharmacophore modeling study based on known Spleen tyrosine kinase inhibitors together with virtual screening for identifying novel inhibitors

In this investigation, chemical features based 3D pharmacophore models were developed based on the known inhibitors of Spleen tyrosine kinase (Syk) with the aid of hiphop and hyporefine modules within catalyst. The best quantitative pharmacophore model, Hypo1, was used as a 3D structural query for retrieving potential inhibitors from chemical databases including Specs, NCI, MayBridge, and Chinese Nature Product Database (CNPD). The hit compounds were subsequently subjected to filtering by Lipinski’s rule of five and docking studies to refine the retrieved hits. Finally 30 compounds were selected from the top ranked hit compounds and conducted an in vitro kinase inhibitory assay. Six compounds showed a good inhibitory potency against Syk, which have been selected for further investigation.     

3D QSAR pharmacophore-based virtual screening for the identification of potential inhibitors of tyrosinase

Abstract

Tyrosinase plays an important role in melanin biosynthesis and protects skin against ultraviolet radiations. Functional deficiency of tyrosinase results in serious dermatological diseases. Tyrosinase also participates in neuromelanin formation in the human brain, which leads to neurodegeneration resulting in Parkinson’s disease. In fruits and vegetables, tyrosinase plays a critical role in senescence, causing undesired browning that results in faster deterioration and shorter shelf lines. The only commercially available tyrosinase is mushroom tyrosinase and it shows the highest homology to the mammalian tyrosinase. Although kojic acid is currently used as a tyrosinase inhibitor, they have serious side effects such as dermatitis, carcinogenesis and hepatotoxicity. Therefore, in order to develop a more active and safer tyrosinase inhibitor, 3D QSAR pharmacophore models were generated based on experimentally known inhibitors. The pharmacophore model, Hypo1, was developed with a large cost difference, high correlation coefficient and low RMS deviation. Hypo1 showed a good spatial arrangement; consisting of five-point features including two hydrogen bond acceptor, one hydrogen bond donor and two hydrophobic features. Hypo1 was further validated by cost analysis, test set and Fisher’s randomisation method. Hypo1 was used as a 3D query for screening the in-house drug-like databases, and the hits were further selected by applying ADMET, Lipinski’s rule of five and fit value criteria. To identify binding conformations, the obtained hits were subjected to molecular docking. Finally, molecular dynamics simulations revealed the appropriate binding modes of hit compounds. To conclude, we propose the final three hit compounds with new structural scaffolds as a virtual candidate as tyrosinase inhibitors.

Communicated by Ramaswamy H. Sarma     

E-pharmacophore-based virtual screening to identify GSK-3β inhibitors

Glycogen synthase kinase-3β (GSK-3β) is a serine/threonine kinase which has attracted significant attention during recent years in drug design studies. The deregulation of GSK-3β increased the loss of hippocampal neurons by triggering apoptosis-mediating production of neurofibrillary tangles and alleviates memory deficits in Alzheimer’s disease (AD). Given its role in the formation of neurofibrillary tangles leading to AD, it has been a major therapeutic target for intervention in AD, hence was targeted in the present study. Twenty crystal structures were refined to generate pharmacophore models based on energy involvement in binding co-crystal ligands. Four common e-pharmacophore models were optimized from the 20 pharmacophore models. Shape-based screening of four e-pharmacophore models against nine established small molecule databases using Phase v3.9 had resulted in 1800 compounds having similar pharmacophore features. Rigid receptor docking (RRD) was performed for 1800 compounds and 20 co-crystal ligands with GSK-3β to generate dock complexes. Interactions of the best scoring lead obtained through RRD were further studied with quantum polarized ligand docking (QPLD), induced fit docking (IFD) and molecular mechanics/generalized Born surface area. Comparing the obtained leads to 20 co-crystal ligands resulted in 18 leads among them, lead1 had the lowest docking score, lower binding free energy and better binding orientation toward GSK-3β. The 50?ns MD simulations run confirmed the stable nature of GSK-3β-lead1 docking complex. The results from RRD, QPLD, IFD and MD simulations confirmed that lead1 might be used as a potent antagonist for GSK-3β.     

Homology modeling and virtual screening studies of FGF-7 protein—a structure-based approach to design new molecules against tumor angiogenesis

Keratinocyte growth factor (KGF) protein is a member of the fibroblast growth factor (FGF) family, which is also known as FGF-7. The FGF-7 plays an important role in tumor angiogenesis. In the present work, FGF-7 is treated as a potential therapeutic target to prevent angiogenesis in cancerous tissue. Computational techniques are applied to evaluate and validate the 3D structure of FGF-7 protein. The active site region of the FGF-7 protein is identified based on hydrophobicity calculations using CASTp and Q-site Finder active site prediction tools. The protein–protein docking study of FGF-7 with its natural receptor FGFR2b is carried out to confirm the active site region in FGF-7. The amino acid residues Asp34, Arg67, Glu116, and Thr194 in FGF-7 interact with the receptor protein (FGFR2b). A grid is generated at the active site region of FGF-7 using Glide module of Schrödinger suite. Subsequently, a virtual screening study is carried out at the active site using small molecular structural databases to identify the ligand molecules. The binding interactions of the ligand molecules, with piperazine moiety as a pharmacophore, are observed at Arg67 and Glu149 residues of the FGF-7 protein. The identified ligand molecules against the FGF-7 protein show permissible pharmacokinetic properties (ADME). The ligand molecules with good docking scores and satisfactory pharmacokinetic properties are prioritized and identified as novel ligands for the FGF-7 protein in cancer therapy.     

Similarity-driven flexible ligand docking

A similarity-driven approach to flexible ligand docking is presented. Given a reference ligand or a pharmacophore positioned in the protein active site, the method allows inclusion of a similarity term during docking. Two different algorithms have been implemented, namely, a similarity-penalized docking (SP-DOCK) and a similarity-guided docking (SG-DOCK). The basic idea is to maximally exploit the structural information about the ligand binding mode present in cases where ligand-bound protein structures are available, information that is usually ignored in standard docking procedures. SP-DOCK and SG-DOCK have been derived as modified versions of the program DOCK 4.0, where the similarity program MIMIC acts as a module for the calculation of similarity indices that correct docking energy scores at certain steps of the calculation. SP-DOCK applies similarity corrections to the set of ligand orientations at the end of the ligand incremental construction process, penalizing the docking energy and, thus, having only an effect on the relative ordering of the final solutions. SG-DOCK applies similarity corrections throughout the entire ligand incremental construction process, thus affecting not only the relative ordering of solutions but also actively guiding the ligand docking. The performance of SP-DOCK and SG-DOCK for binding mode assessment and molecular database screening is discussed. When applied to a set of 32 thrombin ligands for which crystal structures are available, SG-DOCK improves the average RMSD by ca. 1 A when compared with DOCK. When those 32 thrombin ligands are included into a set of 1,000 diverse molecules from the ACD, DIV, and WDI databases, SP-DOCK significantly improves the retrieval of thrombin ligands within the first 10% of each of the three databases with respect to DOCK, with minimal additional computational cost. In all cases, comparison of SP-DOCK and SG-DOCK results with those obtained by DOCK and MIMIC is performed.