Exploring the mechanism of C473D mutation on CDC25B causing weak binding affinity with CDK2/CyclinA by molecular dynamics study.

CDC25B belongs to the CDC25 family, and it plays an important part in regulating the activity of CDK/CyclinA. Studies have shown that CDC25B is closely related to cancer development. When CYS473 on CDC25B is mutated into ASP, the affinity between CDC25B and CDK2/CyclinA weakens, and their dissociation speed is greatly improved. However, the mechanism by which the CDC25BC473D mutant weakens its binding to CDK2/CyclinA is unclear. In order to study the effect of CDC25BC473D mutants on CDK2/CyclinA substrates, we constructed and verified the rationality of the CDC25BWT:CDK2/CyclinA system and CDC25BC473D:CDK2/CyclinA system and conducted molecular dynamics (MD) simulation analysis. In the post-analysis, the fluctuations of residues ARG488-SER499, LYS541-TRP550 on CDC25B and residues ASP206-ASP210 on CDK2 were massive in the mutant CDC25BC473D:CDK2/CyclinA system. And the interactions between residue ARG492 and residue GLU208, residue ARG544 and residue GLU42, residue ARG544 and TRP550 were weakened in the mutant CDC25BC473D:CDK2/CyclinA system. The results showed that when CYS473 on CDC25B was mutated into ASP473, the mutant CDC25BC473D:CDK2/CyclinA system was less stable than the wild-type CDC25BWT:CDK2/CyclinA system. Finally, active site CYS473 of CDC25B was speculated to be the key residue, which had great effects on the binding between CDC25BCYS473 and CDK2 in the CDC25BC473D:CDK2/CyclinA system. Consequently, overall analyses appeared in this study ultimately provided a useful understanding of the weak interactions between CDC25BCYS473D and CDK2/CyclinA.Communicated by Ramaswamy H. Sarma.

Structure-based discovery of a specific SHP2 inhibitor with enhanced blood-brain barrier penetration from PubChem database.

The thiophene [2,3-d]pyrimidine structure-like small molecules were discovered from structure-based virtual screening of 1 billion compounds. Base on enzyme activity assay results, a SHP2-specific molecule inhibitor Comp#2 with IC50 of 1.174 µM, 85-fold more selective for SHP2 than the highly related SHP1 (IC50 > 100 µM). The compound can effectively inhibit SHP2-mediated cell signaling and cancer cell proliferation, including cervix cancer, human pancreatic cancer, large cell lung cancer, and mouse glioma cell. Moreover, the in vivo assay indicated that Comp#2 could inhibit cervix cancer tumors growth in BABL/c mice. This work has shown the specific SHP2 inhibitor can inhibit glioblastoma growth in vivo.

Design, synthesis, biological evaluation and molecular dynamics simulation studies of imidazolidine-2,4-dione derivatives as novel PTP1B inhibitors.

Protein tyrosine phosphatase 1B (PTP1B) is a member of the phosphotyrosine phosphatase family and plays an important role in the signal transduction of diabetes. Inhibition of PTP1B activity can increase insulin sensitivity and reduce blood sugar levels. Therefore, it is urgent to find compounds with novel structures that can inhibit PTP1B. This study designed imidazolidine-2,4-dione derivatives through the computer-aided drug design (CADD) strategy, and the Comp#10 showed outstanding inhibitory ability. (IC50 = 2.07 µM) and selectivity. The inhibitory mechanism at molecular level of Comp#10 on PTP1B was studied by molecular dynamics simulation. The results show that the catalytic region of PTP1B protein is more stable, which makes the catalytic sites unsuitable for exposure. Interestingly, the most obvious changes in the interaction between residues in the P-loop region (such as: His214, Cys215, and Ser216). In short, this study reported for the first time that imidazolidine-2,4-dione derivatives as novel PTP1B inhibitors had good inhibitory activity and selectivity, providing new ideas for the development of small molecule PTP1B inhibitors.

Molecular dynamics study of CDC25BR492L mutant causing the activity decrease of CDC25B.

Cell division cycle 25B (CDC25B) was responsible for regulating the various stages of cell division in the cell cycle. R492L was one of the common types of CDC25B mutants. Researches showed that compared to CDC25BWT, CDC25BR492L mutant had a ∼100-fold reduction in the rate constant for forming phosphatase intermediate (k2). However, the molecular basis of how the CDC25BR492L mutant influenced the process of binding between CDC25B and CDK2/CyclinA was not yet known. Therefore, the optimizations of three-dimensional structure of the CDC25BWT-CDK2/CyclinA system and the CDC25BR492L-CDK2/CyclinA system were constructed by ZDOCK and RDOCK, and five methods were employed to verify the reasonability of the docking structure. Then the molecular dynamics simulations on the two systems were performed to explore the reason why CDC25BR492L mutant caused the weak interactions between CDC25BR492L and CDK2/CyclinA, respectively. The remote docking site (Arg488-Tyr497) and the second active site (Lys538-Arg544) of CDC25B were observed to have high fluctuations in the CDC25BR492L-CDK2/CyclinA system with post-analysis, where the high fluctuation of these two regions resulted in weak interactions between CD25B and CDK2. In addition, Asp38-Glu42 and Asp206-Asp210 of CDK2 showed the slightly descending fluctuation, and CDK2 revealed an enhanced the self-interaction, which made CDK2 keep a relatively stable state in the CDC25BR492L-CDK2/CyclinA system. Finally, Leu492 of CDC25B was speculated to be the key residue, which had great effects on the binding between CDC25BR492L and CDK2 in the CDC25BR492L-CDK2/CyclinA system. Consequently, overall analyses appeared in this study ultimately offered a helpful understanding of the weak interactions between CDC25BR492L and CDK2.

Design, synthesis, biological evaluation and molecular dynamics of LAR inhibitors.

In this study, firstly, the pharmacophore model was established based on LAR inhibitors. ZINC database and drug-like database were screened by Hypo-1-LAR model, and the embryonic compound ZINC71414996 was obtained. Based on this compound, we designed 9 compounds. Secondly, the synthetic route of the compound was determined by consulting Reaxys and Scifinder databases, and 9 compounds (1a-1i) were synthesized by nucleophilic substitution, Suzuki reaction and so on. Meanwhile, their structures were confirmed by 1H NMR and 13C NMR. Thirdly, the Enzymatic assays was carried out, the biological evaluation of compounds 1a-1i led to the identification of a novel PTP-LAR inhibitor 1c, which displayed an IC50 value of 4.8 µM. At last, molecular dynamics simulation showed that compounds could interact strongly with the key amino acids LYS1350, LYS1352, ARG1354, TYR1355, LYS1433, ASP1435, TRP1488, ASP1490, VAL1493, SER1523, ARG1528, ARG1561, GLN1570, LYS1681, thereby inhibiting the protein activity. This study constructed the pharmacophore model of LAR protein, designed small-molecule inhibitors, conducted compound synthesis and enzyme activity screening, so as to provide a basis for searching for drug-capable lead compounds.

Exploring the dynamic mechanism of allosteric drug SHP099 inhibiting SHP2E69K.

The E69K mutation is one of the most frequent protein tyrosine phosphatase-2 (SHP2) mutations in leukemia, and it can cause the increase in the protein activity. Recent studies have shown that the E69K mutation was fairly sensitive to the allosteric inhibitor of SHP2 (SHP099). However, the molecular mechanism of the allosteric drug SHP099 inhibiting SHP2E69K remains unclear. Thus, the molecular dynamic simulations and the post-dynamics analyses (RMSF, PCA, DCCM, RIN and the binding free energies) for SHP2WT, SHP2WT-SHP099, SHP2E69K and SHP2E69K-SHP099 were carried out, respectively. Owing to the strong binding affinity of SHP099 to residues Thr219 and Arg220, the flexibility of linker region (residues Val209-Arg231) was reduced. Moreover, the presence of SHP099 kept the autoinhibition state of the SHP2 protein through enhancing the interactions between the linker region and Q loop in PTP domain, such as Thr219/Val490, Thr219/Asn491, Arg220/Ile488 and Leu254/Asn491. In addition, it was found that the residues (Thr219, Arg220, Leu254 and Asn491) might be the key residues responsible for the conformational changes of protein. Overall, this study may provide an important basis for understanding how the SHP099 effectively inhibited the SHP2E69K activity at the molecular level.

Exploring the mechanism of the potent allosteric inhibitor compound2 on SHP2 WT and SHP2F285S by molecular dynamics study.

Abnormal activation of Ras/MAPK signaling pathway could trigger excessive cell division. Src-homology 2 (SH2) domain-containing protein tyrosine phosphatase (SHP2) could promote Ras/MAPK activation by integrating growth factor signals. Thus, SHP2 inhibitors had become a hot topic in the treatment of cancer. SHP2F285S, mutation in SHP2, was detected in leukemia variants. The compound 2 (3-benzyl-8-chloro-2-hydroxy-4H-benzo[4,5]thiazolo[3,2-a]pyrimidin-4-one) had been reported that it was a potent allosteric inhibitor of both SHP2 wild type (SHP2WT) and the F285S mutant (SHP2F285S). However, the mechanism of inhibition remained to be further discovered. Herein, molecular docking and molecular dynamic (MD) simulation were performed to explain the inhibition mechanism of compound 2 on SHP2WT and SHP2F285S. Overall, the molecular docking analysis revealed that compound 2 maintained the "close" structure of SHP2 protein probably by locking the C-SH2 and PTP domain. Next, post-analysis demonstrated that compound 2 might make TYR66-GLU76 of D'E-loop in N-SH2 and GLU258-LYS266 of B'C-loop, HIS458-ARG465 of P-loop, VAL497-THR507 of Q-loop in PTP domain regions tightly connect and much easier maintain "self-inhibited" conformation of SHP2F285S-compound2 than that of SHP2WT-compound2. Importantly, GLU76 of D'E-loop could play a vital role in inhibition of SHP2WT-compound2 and SHP2F285S-compound2. This work provided a reliable clue to develop novel inhibitors for leukemia related to SHP2F285S.

Synthesis, evaluation, molecular dynamics simulation and targets identification of novel pyrazole-containing imide derivatives.

A new series of novel pyrazole-containing imide derivatives were synthesized and evaluated for their anticancer activities against A-549, Bel7402, and HCT-8 cell lines. Among these compounds A2, A4, A11 and A14 possessed high inhibition activity against A-549 cell lines with IC50 values at 4.91, 3.22, 27.43 and 18.14 µM, respectively, better than that of 5-fluorouracil (IC50=59.27 µM). A2, A4, and A11 also exhibited significant inhibitory activity towards HCT-8 and Bel7402 cell lines. Interestingly, the Heat Shock Protein 90α (Hsp90α, PDB ID: 1UYK) was found to be the potential drug target of these synthesized compounds with the aid of PharmMapper server (http://lilab.ecust.edu.cn/pharmmapper/) and docking module of Schrödinger (Maestro 10.2). Additionally, molecular dynamics simulation was performed out to explore the most likely binding mode of compound A2 with Hsp90α.Communicated by Ramaswamy H. Sarma.

Design, synthesis, biological evaluation, common feature pharmacophore model and molecular dynamics simulation studies of ethyl 4-(phenoxymethyl)-2-phenylthiazole-5-carboxylate as Src homology-2 domain containing protein tyrosine phosphatase-2 (SHP2) inhibitors.

SHP2 is a non-receptor protein tyrosine phosphatase (PTP) encoded by the PTPN11 gene involved in cell death pathway (PD-1/PD-L1) and cell growth and differentiation pathway (MAPK). Moreover, mutations in SHP2 have been implicated in Leopard syndrome (LS), Noonan syndrome (NS), juvenile myelomonocytic leukemia (JMML) and several types of cancer and solid tumors. Thus, SHP2 inhibitors are much needed reagents for evaluation of SHP2 as a therapeutic target. A series of novel ethyl 4-(phenoxymethyl)-2-phenylthiazole-5-carboxylate derivatives were designed and synthesized, and their SHP2 inhibitory activities (IC50) were determined. Among the desired compounds, 1d shares the highest inhibitory activity (IC50 = 0.99 µM) against SHP2. Additionally, a common feature pharmacophore model was established to explain the structure activity relationship of the desired compounds. Finally, molecular dynamics simulation was carried out to explore the most likely binding mode of compound 1d with SHP2. In brief, the findings reported here may at least provide a new strategy or useful insights in discovering novel effective SHP2 inhibitors.Communicated by Ramaswamy H. Sarma.

Identification of the potential dual inhibitor of protein tyrosine phosphatase sigma and leukocyte common antigen-related phosphatase by virtual screen, molecular dynamic simulations and post-analysis.

Owing to their inhibitory role in regulating oligodendrocyte differentiation and apoptosis, protein tyrosine phosphatase sigma (PTPσ) and leukocyte common antigen-related phosphatase (LAR) play a crucial potential role in treating spinal cord injury (SCI) disease. In this research, the computer aided drug design (CADD) methods were applied to discover the potential dual-target drug involving virtual screen, molecular docking and molecular dynamic simulation. Initially, the top 20 compounds with higher docking score than the positive controls (ZINC13749892, ZINC14516161) were virtually screened out from NCI and ZINC databases, and then were submitted in ADMET to predict their drug properties. Among these potential compounds, ZINC72417086 showed a higher docking score and satisfied Lipinski's rule of five. In addition, the post-analysis demonstrated that when ZINC72417086 bound to PTPσ and LAR, it could stable proteins conformations and destroy the residues interactions between P-loop and other loop regions in active pocket. Meanwhile, residue ARG1595 and ARG1528 could play a crucial role in in the inhibition of PTPσ and LAR, respectively. This research offered a novel approach for rapid discovery of dual-target leads compounds to treat SCI.Communicated by Ramaswamy H. Sarma.

Allosteric Inhibitors of SHP2: An Updated Patent Review (2015-2020).

Srchomology-2-domain-containing PTP 2 (SHP2) is a nonreceptor phosphatase encoded by the PTPN11 gene. Over expression of SHP2 is associated with various human diseases, such as Noonan syndrome, LEOPARD syndrome, and cancers. To overcome the shortcomings of existing orthosteric inhibitors, novel inhibitors targeting the allosteric site of SHP2 with high selectivity and low toxicity are under development. This paper reviews allosteric inhibitors of SHP2 published in patents from 2015 to 2020. The molecules are classified according to the chemical structure of the central core. SHP2 has long been considered as an 'undruggable' protein. Fortunately, a critical breakthrough was made by researchers from Novartis AG Ltd., who identified SHP099 as a highly potent, selective, soluble, and orally bioavailable SHP2 allosteric inhibitor. Currently, there are several allosteric inhibitors of SHP2 in clinical development. However, drug resistance is still a major challenge. The combination of SHP2 allosteric inhibitors and immunotherapy drugs or molecular targeted drugs is emerging as a promising therapeutic strategy against drug resistance.

Scaffold-based selective SHP2 inhibitors design using core hopping, molecular docking, biological evaluation and molecular simulation.

PTPN11 (coding the gene of SHP2), a classic non-receptor protein tyrosine phosphatase, is implicated in multiple cell signaling pathway. Abnormal activation of SHP2 has been shown to contribute to a variety of human diseases, including Juvenile myelomonocytic leukemia (JMML), Noonan syndrome and tumors. Thus, the SHP2 inhibitors have important therapeutic value. Here, based on the compound PubChem CID 8,478,960 (IC50 = 45.01 µM), a series of thiophene [2,3-d] pyrimidine derivatives (IC50 = 0.4-37.87 µM) were discovered as novel and efficient inhibitors of SHP2 through powerful "core hopping" and CDOCKER technology. Furthermore, the SHP2-PTP phosphatase activity assay indicated that Comp#5 (IC50 = 0.4 µM) was the most active SHP2 inhibitor. Subsequently, the effects of Comp#5 on the structure and function of SHP2 were investigated through molecular dynamics (MD) simulation and post-kinetic analysis. The result indicated that Comp#5 enhanced the interaction of residues THR357, ARG362, LYS366, PRO424, CYS459, SER460, ALA461, ILE463, ARG465, THR507 and GLN510 with the surrounding residues, improving the stability of the catalytic active region and the entrance of catalytic active region. In particular, the Comp#5 conjugated with residue ARG362, elevating the efficient and selectivity of SHP2 protein. The study here may pave the way for discovering the novel SHP2 inhibitors for suffering cancer patients.

Design, synthesis and biological evaluation of pyridine derivatives as selective SHP2 inhibitors.

SHP2 is a non-receptor protein tyrosine phosphatase encoded by the PTPN11 gene, which affects the transduction of multiple signaling pathways, including RAS-ERK, PI3K-AKT and JAK-STAT. SHP2 also plays an important role in the programmed cell death pathway (PD-1/PD-L1). Studies have shown that SHP2 is associated with a variety of cancers, including breast, liver and gastric cancers. Therefore, the development of SHP2 inhibitors has attracted extensive attention. In this study, based on the known inhibitor 1 (SHP099), novel SHP2 inhibitors were designed by means of scaffold hopping, and 35 pyridine derivatives as SHP2 inhibitors were found. The in vitro enzyme activity assay was performed on these compounds, and multiple selective SHP2 inhibitors with activity potency similar to that of SHP099 were obtained. Among them, compound (2-(4-(aminomethyl)piperidin-1-yl)-5-(2,3-dichlorophenyl)pyridin-3-yl)methanol (11a) was the most potent and highly selective SHP2 inhibitor with an in vitro enzyme activity IC50 value of 1.36 µM. Fluorescence titration assay verified that 11a bound directly to SHP2 protein. Subsequently, cell assay of representative compounds showed that these compounds could effectively inhibit the proliferation of Ba/F3 cells. In addition, the pharmacokinetic characteristics of the designed compounds were analyzed by the in silico ADMET prediction. Molecular docking study provided more detailed information on the binding mode of compounds and SHP2 protein. In brief, this study reported for the first time that pyridine derivatives as novel SHP2 inhibitors had good inhibitory activity and selectivity, providing new clues for the development of small molecule SHP2 inhibitors.

Identification of protein tyrosine phosphatase 1B (PTP1B) inhibitors through De Novo Evoluton, synthesis, biological evaluation and molecular dynamics simulation.

Protein tyrosine phosphatase 1B (PTP1B) is a widely expressed 50 kDa enzyme and the first intracellular PTP to be purified from human placental tissue. It has been proved that protein tyrosine phosphatase 1B played a significant role in the negative regulation of insulin signaling pathway and overexpression of PTP1B could lead to the decrease of insulin resistance. Therefore PTP1B has emerged as a novel promising therapeutic target for the treatment of type-2 diabetes mellitus. Computer aided drug design (CADD), chemical synthesis and biological activity assay resulted in the identification of a novel potent PTP1B inhibitor, compound 1a, which shared an IC50 value of 4.46 µM. Finally, the analysis of molecular dynamics simulation provided the theoretical basis for favorable activity of compound 1a.

Investigating the reason for loss-of-function of Src homology 2 domain-containing protein tyrosine phosphatase 2 (SHP2) caused by Y279C mutation through molecular dynamics simulation.

Noonan syndrome with multiple lentigines (NSML), formerly known as LEOPARD syndrome (LS), is an autosomal dominant inherited multisystemic disorder. Most patients involve mutation in SHP2 encoded by tyrosine-protein phosphatase non-receptor type 11 (PTPN11) gene. Studies have shown that NSML-associated Y279C mutation exhibited the reduced phosphatase activity, leading to loss-of-function (LOF) of SHP2. However, the effect of the Y279C mutation on the SHP2 at the molecular level is unclear. In this study, molecular dynamics simulations of SHP2 wild-type (SHP2WT) and Y279C mutant (SHP2Y279C) were performed to investigate the structural differences in proteins after Y279C mutation and to find out the reason for loss-of-function of SHP2. Through a series of post-dynamic analyses, it was found that the protein occupied a smaller phase space after Y279C mutation, showing reduced flexibility. Specifically, due to the mutation of Y279C, the secondary structures of these two regions (residues Lys70-Ala72 and Gly462-Arg465) were significantly transformed from Turn to α-helix and ß-strand. Furthermore, by calculating the residue interaction network, hydrogen bond occupancy and binding free energy, it was further revealed that the conformational differences between SHP2WT and SHP2Y279C systems were mainly caused by the differences in the interaction between Arg465-Phe469, Ile463-Gly467, Cys279-Lys70, Cys459-Ala72, Gly464-Phe71, Phe71-Ile463, Ile463-Ala505 and Arg465-Glu361. Consequently, this finding is expected to provide a new insight into the reason for loss-of-function of SHP2 caused by Y279C mutation.Communicated by Ramaswamy H. Sarma.

Probing the acting mode and advantages of RMC-4550 as an Src-homology 2 domain-containing protein tyrosine phosphatase (SHP2) inhibitor at molecular level through molecular docking and molecular dynamics.

The over-activation of Ras/mitogen-activated protein kinase (MAPK) signaling pathway associated with a variety of cancers is usually related with abnormal activation of Src-homology 2 domain-containing protein tyrosine phosphatase (SHP2). For this purpose, SHP2 has attracted extensive interest as a potential target for cancer treatment. RMC-4550, as a newly developed selective inhibitor of SHP2, possesses an overwhelming advantage over the previous generation inhibitor SHP099 in terms of in vitro activity. However, the binding mode of SHP2 with RMC-4550 and the reason for the high efficiency of RMC-4550 as SHP2 inhibitor at molecular level are still unclear. Therefore, in this study, the binding mode of RMC-4550 with SHP2 and the superiorities of RMC-4550 as inhibitor at binding affinity and dynamic interactive behavior with SHP2 were probed by molecular docking and molecular dynamics (MD) simulations. By comparing the results of molecular docking, it was found that SHP2 formed more tight interaction with RMC-4550 than that with SHP099. Subsequently, a series of post-dynamic analyses on three simulation trajectories (SHP2WT, SHP2SHP099 and SHP2RMC-4550) were performed and found that the SHP2 protein bound with RMC-4550 maintained a firmer interaction between N-Src-homology 2 (N-SH2) and PTP domain throughout the MD simulation, leading to a more stable protein conformation. The finding here provides new clues for the design of SHP2 inhibitor against the over-activation of Ras/MAPK pathway.Communicated by Ramaswamy H. Sarma.

Design, synthesis, biological evaluation and molecular dynamics simulation studies of (R)-5-methylthiazolidin-4-One derivatives as megakaryocyte protein tyrosine phosphatase 2 (PTP-MEG2) inhibitors for the treatment of type 2 diabetes.

PTP-MEG2 plays a significant role in insulin production and is able to enhance insulin signaling and improve insulin sensitivity. So, PTP-MEG2 inhibitors are closely associated with type 2 diabetes therapy. A series of novel (R)-5-methylthiazolidin-4-one derivatives were designed and synthesized, and their PTP-MEG2 inhibitory activities (IC50) were determined. Among the desired compounds, 1h shares the highest inhibitory activity (IC50 = 1.34 µM) against PTP-MEG2. Additionally, various post-dynamic analyses confirmed that when compound 1h bound to the PTP-MEG2, the protein conformations became unstable and the function of the pTyr recognition loop (Asn331-Cys338) would be disturbed. And thus, the ideal conformations needed for the catalytic activity was difficult to be maintained. In brief, these might be how the compound 1h worked. Furthermore, we also found that the key residues Arg332 would play a critical role in disturbing the residue interactions. AbbreviationsDCCMdynamic cross-correlation mappingDMFN,N-dimethylformamideDSSPdefinition of secondary structure of proteinsFOXOforkhead transcription factorsMDmolecular dynamicsPCAprincipal component analysisPDBprotein data bankPTKsprotein tyrosine kinasesPTPsprotein tyrosine phosphatasesPTP-MEG2megakaryocyte protein tyrosine phosphatase 2RINresidue interaction networkRINGResidue Interaction Network GeneratorRMSDroot means square deviationRMSFroot mean square fluctuationCommunicated by Ramaswamy H. Sarma.

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

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.

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, synthesis, biological evaluation and molecular dynamics studies of 4-thiazolinone derivatives as protein tyrosine phosphatase 1B (PTP1B) inhibitors.

Protein tyrosine phosphatase 1B (PTP1B) is a key negative regulator of insulin signaling pathway, and more and more studies have shown that it is a potential target for the treatment of type 2 diabetes mellitus (T2DM). In this study, 17 new 4-thiazolinone derivatives were designed and synthesized as novel PTP1B inhibitors, and ADMET prediction confirmed that these compounds were to be drug-like. In vitro enzyme activity experiments were performed on these compounds, and it was found that a plurality of compounds had good inhibitory activity and high selectivity against PTP1B protein. Among them, compound 7p exhibited the best inhibitory activity with an IC50 of 0.92 µM. The binding mode of compound 7p and PTP1B protein was explored, revealing the reason for its high efficiency. In addition, molecular dynamics simulations for the PTP1BWT and PTP1Bcomp#7p systems revealed the effects of compound 7p on PTP1B protein at the molecular level. In summary, the study reported for the first time that 4-thiazolinone derivatives as a novel PTP1B inhibitor had good inhibitory activity and selectivity for the treatment of T2DM, providing more options for the development of PTP1B inhibitors. AbbreviationsBBBblood-brain barrierCDC25Bcell division cycle 25 homolog BCYP2D6Cytochrome P450 2D6 bindingDCCMdynamic cross-correlation mapDSDiscovery StudioH bondhydrogen bondHIAhuman intestinal absorptionLARleukocyte antigen-related phosphataseMDmolecular dynamicsMEG-2maternal-effect germ-cell defective 2MM-PBSAmolecular mechanics Poisson Boltzmann surface area)PCAprincipal component analysisPDBProtein Data BankpNPPp-nitrophenyl phosphatePPBplasma protein bindingPTP1Bprotein tyrosine phosphotase 1BRMSDroot mean square deviationRMSFroot mean square fluctuationSHP-1src homologous phosphatase-1SHP-2src homologous phosphatase-2SPCsingle-point chargeTCPTPT cell protein tyrosine phosphataseT2DMType 2 diabetes mellitusVDWvan der WaalsCommunicated by Ramaswamy H. Sarma.