Flavonoids as potent allosteric inhibitors of protein tyrosine phosphatase 1B: molecular dynamics simulation and free energy calculation

Protein tyrosine phosphatase 1B (PTP1B) is a member of the PTP superfamily which is considered to be a negative regulator of insulin receptor (IR) signaling pathway. PTP1B is a promising drug target for the treatment of type 2 diabetes, obesity, and cancer. The existence of allosteric site in PTP1B has turned the researcher’s attention to an alternate strategy for inhibition of this enzyme. Herein, the molecular interactions between the allosteric site of PTP1B with three non-competitive flavonoids, (MOR), (MOK), and (DPO) have been investigated. Three ligands were docked into allosteric site of the enzyme. The resulting protein–ligand complexes were used for molecular dynamics studies. Principal component and free-energy landscape (FEL) as well as cluster analyses were used to investigate the conformational and dynamical properties of the protein and identify representative enzyme substrates bounded to the inhibitors. Per residue energy decomposition analysis attributed dissimilar affinities of three inhibitors to the several hydrogen bonds and non-bonded interactions. In conclusion, our results exhibited an inhibitory pattern of the ligands against PTP1B.     

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.     

Identification of potential inhibitors for HCV NS5b of genotype 4a by combining dynamic simulation,protein–ligand interaction fingerprint, 3D pharmacophore,docking and 3D QSAR

Abstract

HCV NS5B polymerase has been one of the most attractive targets for developing new drugs for HCV infection and many drugs were successfully developed, but all of them were designed for targeting Hepatitis C Virus genotype 1 (HCV GT1). Hepatitis C virus genotype 4a (HCV GT4a) dominant in Egypt has paid less attention. Here, we describe our protocol of virtual screening in identification of novel potential potent inhibitors for HCV NS5B polymerase of GT4a using homology modeling, protein–ligand interaction fingerprint (PLIF), docking, pharmacophore, and 3D CoMFA quantitative structure activity relationship (QSAR). Firstly, a high-quality 3D model of HCV NS5B polymerase of GT4a was constructed using crystal structure of HCV NS5B polymerase of GT1 (PDB ID: 3hkw) as a template. Then, both the model and the template were simulated to compare conformational stability. PLIF was generated using five crystal structures of HCV NS5B (PDB ID: 4mia, 4mib, 4mk9, 4mka, and 4mkb), which revealed the most important residues and their interactions with the co-crystalized ligands. After that, a 3D pharmacophore model was developed from the generated PLIF data and then used as a screening filter for 17000328 drug-like zinc database compounds. 900 compounds passed the pharmacophore filter and entered the docking-based virtual screening stage. Finally, a 3D CoMFA QSAR was developed using 42 compounds as a training and 19 compounds as a test set. The 3D CoMFA QSAR was used to design and screen some potential inhibitors, these compounds were further evaluated by the docking stage. The highest ranked five hits from docking result (compounds (p1–p4) and compound q1) were selected for further analysis.

Communicated by Ramaswamy H. Sarma     

3D-QSAR pharmacophore-based virtual screening,molecular docking and molecular dynamics simulation toward identifying lead compounds for NS2B–NS3 protease inhibitors

NS2B–NS3 protease has been identified to serve as lead drug design target due to its significant role in West Nile viral (WNV) and dengue virus (DENV) reproduction and replication. There are currently no approved chemotherapeutic drugs and effective vaccines to inhibit DENV and WNV infections. In this work, 3D-QSAR pharmacophore model has been developed to discover potential inhibitory candidates. Validation through Fischer’s model and decoy test indicate that the developed 3D pharmacophore model is highly predictive for DENV inhibitors, which was then employed to screen ZINC chemical library to obtain reasonable hits. Following ADMET filtering, 15 hits were subjected to further filter through molecular docking and CoMFA modeling. Finally, top three hits were identified as lead compounds or potential inhibitory candidates with IC₅₀ values of ～0.4637?µM and fitness of ～57.73. It is implied from CoMFA modeling that substituents at the side site of benzotriazole such as a p-nitro group (e.g. biphenyl head) and a carbonyl (e.g. carboxylate function) at the side site of furan or amino group may improve bioactivity of ZINC85645245, respectively. Molecular dynamics simulations (MDS) were performed to discover new interactions and reinforce the binding modes from docking for the hits also. The QSAR and MDS results obtained from this work should be useful in determining structural requirements for inhibitor development as well as in designing more potential inhibitors for NS2B–NS3 protease.     

Stepwise development of structure–activity relationship of diverse PARP-1 inhibitors through comparative and validated in silico modeling techniques and molecular dynamics simulation

Inhibitors of poly (ADP-ribose) polymerase-1 (PARP-1) enzyme are useful for the treatment of various diseases including cancer. Comparative in silico studies were performed on different ligand-based (2D-QSAR, Kernel-based partial least square (KPLS) analysis, Pharmacophore Search Engine (PHASE) pharmacophore mapping), and structure-based (molecular docking, MM-GBSA analyses, Gaussian-based 3D-QSAR analyses on docked poses) modeling techniques to explore the structure–activity relationship of a diverse set of PARP-1 inhibitors. Two-dimensional (2D)-QSAR highlighted the importance of charge topological index (JGI7), fractional polar surface area (Jurs_FPSA3), and connectivity index (CIC2) along with different molecular fragments. Favorable and unfavorable fingerprints were demonstrated in KPLS analysis, whereas important pharmacophore features (one acceptor, one donor, and two ring aromatic) along with favorable and unfavorable field effects were demonstrated in PHASE-based pharmacophore model. MM-GBSA analyses revealed significance of different polar, non-polar, and solvation energies. Docking-based alignment of ligands was used to perform Gaussian-based 3D-QSAR study that further demonstrated importance of different field effects. Overall, it was found that polar interactions (hydrogen bonding, bridged hydrogen bonding, and pi–cation) play major roles for higher activity. Steric groups increase the total contact surface area but it should have higher fractional polar surface area to adjust solvation energy. Structure-based pharmacophore mapping spotted the positive ionizable feature of ligands as the most important feature for discriminating highly active compounds from inactives. Molecular dynamics simulation, conducted on highly active ligands, described the dynamic behaviors of the protein complexes and supported the interpretations obtained from other modeling analyses. The current study may be useful for designing PARP-1 inhibitors.     

Rational design,synthesis, and structure–activity relationships of 5-amino-1H-pyrazole-4-carboxylic acid derivatives as protein tyrosine phosphatase 1B inhibitors

A series of novel amino-carboxylic based pyrazole as protein tyrosine phosphatase 1B (PTP1B) inhibitors were designed on the basis of structure-based pharmacophore model and molecular docking. Compounds containing different hydrophobic tail (1,2-diphenyl ethanone, oxdiadizole and dibenzyl amines) were synthesized and evaluated in PTP1B enzymatic assay. Structure–activity relationship based optimization resulted in identification of several potent, metabolically stable and cell permeable PTP1B inhibitors.     

In silico approaches to identify novel myeloid cell leukemia-1 (Mcl-1) inhibitors for treatment of cancer

Myeloid cell leukemia-1 (Mcl-1) has been a validated and attractive target for cancer therapy. Over-expression of Mcl-1 in many cancers allows cancer cells to evade apoptosis and contributes to the resistance to current chemotherapeutics. Here, we identified new Mcl-1 inhibitors using a multi-step virtual screening approach. First, based on two different ligand-receptor complexes, 20 pharmacophore models were established by simultaneously using ‘Receptor-Ligand Pharmacophore Generation’ method and manual build feature method, and then carefully validated by a test database. Then, pharmacophore-based virtual screening (PB-VS) could be performed by using the 20 pharmacophore models. In addition, docking study was used to predict the possible binding poses of compounds, and the docking parameters were optimized before performing docking-based virtual screening (DB-VS). Moreover, a 3D QSAR model was established by applying the 55 aligned Mcl-1 inhibitors. The 55 inhibitors sharing the same scaffold were docked into the Mcl-1 active site before alignment, then the inhibitors with possible binding conformations were aligned. For the training set, the 3D QSAR model gave a correlation coefficient r² of 0.996; for the test set, the correlation coefficient r² was 0.812. Therefore, the developed 3D QSAR model was a good model, which could be applied for carrying out 3D QSAR-based virtual screening (QSARD-VS). After the above three virtual screening methods orderly filtering, 23 potential inhibitors with novel scaffolds were identified. Furthermore, we have discussed in detail the mapping results of two potent compounds onto pharmacophore models, 3D QSAR model, and the interactions between the compounds and active site residues.     

Human protein tyrosine phosphatase 1B inhibitors: QSAR by genetic function approximation

Protein tyrosine phosphatase 1B (PTP 1B), a negative regulator of insulin receptor signaling system, has emerged as a highly validated, attractive target for the treatment of non-insulin dependent diabetes mellitus (NIDDM) and obesity. As a result there is a growing interest in the development of potent and specific inhibitors for this enzyme. This quantitative structure-activity relationship (QSAR) study for a series of formylchromone derivatives as PTP lB inhibitors was performed using genetic function approximation (GFA) technique. The QSAR models were developed using a training set of 29 compounds and the predictive ability of the QSAR model was evaluated against a test set of 7 compounds. The internal and external consistency of the final QSAR model was 0.766 and 0.785. The statistical quality of QSAR models was assessed by statistical parameters r2, r2 (crossvalidated r2), r2pred (predictive r2) and lack of fit (LOF) measure. The results indicate that PTP lB inhibitory activity of the formylchromone derivatives is strongly dependent on electronic, thermodynamic and shape related parameters.     

Identification of potential leukocyte antigen-related protein (PTP-LAR) inhibitors through 3D QSAR pharmacophore-based virtual screening and molecular dynamics simulation

Abstract

Owing to its negative regulatory role in insulin signaling, protein tyrosine phosphatase of leukocyte antigen-related protein (PTP-LAR) was widely thought as a potential drug target for diabetes. Now, it was urgent to search for potential LAR inhibitors targeting diabetes. Initially, the pharmacophore models of LAR inhibitors were established with the application of the HypoGen module. The cost analysis, test set validation, as well as Fischer’s test was used to verify the efficiency of pharmacophore model. Then, the best pharmacophore model (Hypo-1-LAR) was applied for the virtual screening of the ZINC database. And 30 compounds met the Lipinski’s rule of five. Among them, 10 compounds with better binding affinity than the known LAR inhibitor (BDBM50296375) were discovered by docking studies. Finally, molecular dynamics simulations and post-analysis experiments (RMSD, RMSF, PCA, DCCM and RIN) were conducted to explore the effect of ligands (ZINC97018474 and Compound 1) on LAR and preliminary understand why ZINC97018474 had better inhibitory activity than Compound 1 (BDBM50296375).

Communicated by Ramaswamy H. Sarma     

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     

Design and evaluation of non-carboxylate 5-arylidene-2-thioxo-4-imidazolidinones as novel non-competitive inhibitors of protein tyrosine phosphatase 1B

Protein tyrosine phosphatase 1B (PTP1B) acts as a negative regulator of insulin and leptin signalling and is crucially involved in the development of type 2 diabetes mellitus, obesity, cancer and neurodegenerative diseases. Pursuing our efforts to identify PTP1B inhibitors endowed with drug-like properties, we designed and evaluated 3-aryl-5-arylidene-2-thioxo-4-imidazolidinones (7) as a novel class of non-carboxylate PTP1B inhibitors. In agreement with our design, kinetic studies demonstrated that selected compounds 7 act as reversible, non-competitive inhibitors of the target enzyme at low micromolar concentrations. Accordingly, molecular docking experiments suggested that these inhibitors can fit an allosteric site of PTP1B that we previously individuated. Moreover, cellular assays demonstrated that compound 7e acts as a potent insulin-sensitizing agent in human liver HepG2 cells. Taken together, our results showed that these non-competitive PTP1B inhibitors can be considered promising lead compounds aimed to enhance druggability of the target enzyme and identify novel antidiabetic drugs.     

Polycyclic phloroglucinols as PTP1B inhibitors from Hypericum longistylum: Structures,PTP1B inhibitory activities,and interactions with PTP1B

Protein tyrosine phosphatase 1B (PTP1B) has been regarded as a target for the research and development of new drugs to treat type II diabetes and PTP1B inhibitors are potential lead compounds for this type of new drugs. A phytochemical investigation to obtain new PTP1B inhibitors resulted in the isolation of four new phloroglucinols, longistyliones A–D (1–4) from the aerial parts of Hypericum longistylum. The structures of 1–4 were elucidated on the basis of extensive 1D and 2D NMR spectroscopic data analysis, and the absolute configurations of these compounds were established by comparing their experimental electronic circular dichroism (ECD) spectra with those calculated by the time-dependent density functional theory method. Compounds 1–4 possess a rare polycyclic phloroglucinol skeleton. The following biological evaluation revealed that all of the compounds showed PTP1B inhibitory effects. The further molecular docking studies indicated the strong interactions between these bioactive compounds with the PTP1B protein, which revealed the possible mechanism of PTP1B inhibition of bioactive compounds. All of the results implied that these compounds are potentially useful for the treatment of type II diabetes.     

Molecular docking and virtual screening for novel protein tyrosine phosphatase 1B (PTP1B) inhibitors

Protein tyrosine phosphatase 1B (PTP1B) functions as major negative regulator of insulin and leptin signaling pathways. In view of this, PTP1B is an significant target for drug development against cancer, diabetes and obesity. The aim of the current study is to identify PTP1B inhibitors by means of virtual screening with docking. 523,366 molecules from ZINC database have been screened and based on DOCK grid scores and hydrogen bonding interactions five new potential inhibitors were identified. ZINC12502589, ZINC13213457, ZINC25721858, ZINC31392733 and ZINC04096400 were identified as potential lead molecules for inhibition of PTP1B. The identified molecules were subjected to Lipinski''s rule of five parameters and found that they did not violate any rule. More specific analysis of pharmacological parameters may be scrutinized through a complete ADME/Tox evaluation. Pharma algorithm was used to Calculate ADME–Tox profiles for such molecules. In general, all the molecules presented advantages and as well as disadvantages when compared to each other. No marked difference in health effects and toxicity profiles were observed among these molecules.     

Molecular modeling of protein tyrosine phosphatase 1B (PTP 1B) inhibitors

Binding modes of a series of aryloxymethylphosphonates and monoanionic biosteres of phosphate group from a series of benzylic alpha,alpha-diflluoro phosphate and its biosteres as protein tyrosine phosphatase 1B (PTP 1B) inhibitors have been identified by molecular modeling techniques. We have performed docking and molecular dynamics simulations of these inhibitors with PTP 1B enzyme. The initial conformation of the inhibitors for docking was obtained from simulated annealing technique. Solvent accessible surface area calculations suggested that active site of PTP 1B is highly hydrophobic. The results indicate that for aryloxymethylphosphonates, in addition to hydrogen bonding interactions, Tyr46, Arg47, Asp48, Val49, Glu115, Lys116, Lys120 amino acid residues of PTP 1B are responsible for governing inhibitor potency of the compounds. The sulfonate and tetrazole functional groups have been identified as effective monoanionic biosteres of phosphate group and biphenyl ring system due to its favorable interactions with Glu115, Lys116, Lys120 residues of PTP 1B found to be more suitable aromatic functionality than naphthalene ring system for benzylic alpha,alpha-diflluoro phosphate and its biosteres. The information generated from the present study should be useful in the design of more potent PTP 1B inhibitors as anti diabetic agents.     

Discovery of novel PTP1B inhibitors via pharmacophore-oriented scaffold hopping from Ertiprotafib

An integrated molecular design strategy combining pharmacophore recognition and scaffold hopping was exploited to discover novel PTP1B inhibitors based on the known PTP1B inhibitor Ertiprotafib. A composite pharmacophore model was proposed from the interaction mode of Ertiprotafib, and 21 diverse molecules from five distinct structural classes were designed and synthesized accordingly. New compounds with considerable inhibition against PTP1B were identified from each series, and the most active compound 3a showed IC₅₀ value of 1.3 μmol L^?1 against human recombinant PTP1B. Docking study indicated that the new inhibitors assumed binding modes similar to that of Ertiprotafib.     

In silico study of M18 aspartyl amino peptidase (M18AAP) of Plasmodium vivax as an antimalarial drug target

Plasmodium vivax (Pv) is the second most malaria causing pathogen among Plasmodium species. M18 aspartic aminopeptidase (M18AAP) protein is a single gene copy present in Plasmodium. This protein is functional at the terminal stage of hemoglobin degradation of host and completes the hydrolysis process which makes it an important target for new chemotherapeutics. No experimental and structural study on M18AAP protein of P. vivax is reported till today. This paper advocates the application of multiple computational approaches like protein model prediction, ligand-based 3D QSAR study, pharmacophore, structure-based virtual screening and molecular docking simulation for identification of potent lead molecules against the enzyme. The 3D QSAR model was developed using known bioactive compounds against the PvM18AAP protein which statistically signify the k-NN model with q^2 = 0.7654. The study reports a lead molecule from ligand-centric approach with good binding affinity and possessing lowest docking score. The findings will be helpful for in-vivo and in-vitro validations and development of potent anti-malarial molecules against the drug resistant strains of malaria parasite.     

Discovery of potent PTP1B inhibitors via structure-based drug design,synthesis and in vitro bioassay of Norathyriol derivatives

Protein tyrosine phosphatase 1B (PTP1B) has recently been identified as a potential target of Norathyriol. Unfortunately, Norathyriol is not a potent PTP1B inhibitor, which somewhat hinders its further application. Based on the fact that no study on the relationship of chemical structure and PTP1B inhibitory activity of Norathyriol has been reported so far, we attempted to perform structural optimization so as to improve the potency for PTP1B. Via structure-based drug design (SBDD), a rational strategy based on the binding mode of Norathyriol to PTP1B, we designed 26 derivatives with substitutions at the four phenolic hydroxyl groups of Norathyriol. By chemical synthesis and in vitro bioassay, we identified seven PTP1B inhibitors that were more potent than Norathyriol, of which XWJ24 showed the highest potency (IC₅₀: 0.6 μM). We also found out that XWJ24 was a competitive inhibitor and showed the 4.5-fold selectivity over its close homolog, TC-PTP. Through molecular docking of XWJ24 against PTP1B, we highlighted the essential role of its hydrogen bond with Asp181 for PTP1B inhibition and identified a potential halogen bond with Asp48 that was not observed for Norathyriol. The current data indicate that our SBDD strategy is effective to discover potent PTP1B-targeted Norathyriol derivatives, and XWJ24 is a promising lead compound for further development.     

3D QSAR,pharmacophore and molecular docking studies of known inhibitors and designing of novel inhibitors for M18 aspartyl aminopeptidase of <Emphasis Type="Italic">Plasmodium falciparum</Emphasis>

Background

The Plasmodium falciparum M18 Aspartyl Aminopeptidase (PfM18AAP) is only aspartyl aminopeptidase which is found in the genome of P. falciparum and is essential for its survival. The PfM18AAP enzyme performs various functions in the parasite and the erythrocytic host such as hemoglobin digestion, erythrocyte invasion, parasite growth and parasite escape from the host cell. It is a valid target to develop antimalarial drugs. In the present work, we employed 3D QSAR modeling, pharmacophore modeling, and molecular docking to identify novel potent inhibitors that bind with M18AAP of P. falciparum.

Results

The PLSR QSAR model showed highest value for correlation coefficient r² (88 %) and predictive correlation coefficient (pred_r2) =0.6101 for external test set among all QSAR models. The pharmacophore modeling identified DHRR (one hydrogen donor, one hydrophobic group, and two aromatic rings) as an essential feature of PfM18AAP inhibitors. The combined approach of 3D QSAR, pharmacophore, and structure-based molecular docking yielded 10 novel PfM18AAP inhibitors from ChEMBL antimalarial library, 2 novel inhibitors from each derivative of quinine, chloroquine, 8-aminoquinoline and 10 novel inhibitors from WHO antimalarial drugs. Additionally, high throughput virtual screening identified top 10 compounds as antimalarial leads showing G-scores -12.50 to -10.45 (in kcal/mol), compared with control compounds(G-scores -7.80 to -4.70) which are known antimalarial M18AAP inhibitors (AID743024). This result indicates these novel compounds have the best binding affinity for PfM18AAP.

Conclusion

The 3D QSAR models of PfM18AAP inhibitors provided useful information about the structural characteristics of inhibitors which are contributors of the inhibitory potency. Interestingly, In this studies, we extrapolate that the derivatives of quinine, chloroquine, and 8-aminoquinoline, for which there is no specific target has been identified till date, might show the antimalarial effect by interacting with PfM18AAP.

     

Screening for the selective inhibitors of MMP-9 from natural products based on pharmacophore modeling and molecular docking in combination with bioassay experiment,hybrid QM/MM calculation,and MD simulation

Matrix metalloproteinase-9 (MMP-9) has been considered as an attractive target involving cancer therapy. In this study, the 3D QSAR pharmacophore model of MMP-9 inhibitors is built, and its reliability is subsequently validated based on different methods. The built pharmacophore model consists of the four chemical features, including two hydrogen bond acceptors (HBA), one hydrophobic (HY), and one ring aromatic (RA). Among them, both HY and RA are found to be especially important features because they involve the interactions of inhibitors with the S1′ pocket of MMP-9, which determines the selectivity of MMP-9 inhibitors. By combining pharmacophore model with molecular docking, the virtual screening is carried out to identify the selective MMP-9 inhibitors from natural products. The four potential selective MMP-9 inhibitors of natural products are found. One of them was used to carry out the bioassay experiment inhibiting MMP-9, and the estimated IC₅₀ value of only 26.94 µM clearly shows its strongly inhibitory activity; besides, both the hybrid quantum mechanics/molecular mechanics (QM/MM) calculation and the molecular dynamics simulation are performed to examine the reliability regarding the binding mode of this inhibitor with MMP-9 active sites predicted by molecular docking. All the screened four natural products are found to well bind with the MMP-9 active sites by different kinds of interactions. Finally, the ADMET properties of screened four natural products are assessed. These screened MMP-9 inhibitors of natural products could be used as the lead compounds to perform structural modifications and optimizations in the future work.

Communicated by Ramaswamy H. Sarma