Steered Molecular Dynamics for Investigating the Interactions Between Insulin Receptor Tyrosine Kinase (IRK) and Variants of Protein Tyrosine Phosphatase 1B (PTP1B)

The aim of this study is to use steered molecular dynamics to investigate the dissociation process between IRK and PTP1Bs for wild type and five mutants (consisting of p.D181E, p.D181A, p.Q262A, p.D181A-Y46F, and p.D181A-Q262A). The gained results are observed not only the unbinding mechanism of IRK-PTP1B complexes came from pulling force profile, number of hydrogen bonds, and interaction energy between IRK and PTP1Bs but also described PTP1B’s point mutations could variably change its binding affinity towards IRK. Additionally, the binding free energy calculated by Molecular Mechanics/Poisson-Boltzmann Surface Area (MM-PBSA) is also revealed that electrostatic energy and polar solvation energy mainly made up the binding free energy of PTP1B-IRK complexes.     

Assessing the performance of the MM/PBSA and MM/GBSA methods. 1. The accuracy of binding free energy calculations based on molecular dynamics simulations

The Molecular Mechanics/Poisson-Boltzmann Surface Area (MM/PBSA) and the Molecular Mechanics/Generalized Born Surface Area (MM/GBSA) methods calculate binding free energies for macromolecules by combining molecular mechanics calculations and continuum solvation models. To systematically evaluate the performance of these methods, we report here an extensive study of 59 ligands interacting with six different proteins. First, we explored the effects of the length of the molecular dynamics (MD) simulation, ranging from 400 to 4800 ps, and the solute dielectric constant (1, 2, or 4) on the binding free energies predicted by MM/PBSA. The following three important conclusions could be observed: (1) MD simulation length has an obvious impact on the predictions, and longer MD simulation is not always necessary to achieve better predictions. (2) The predictions are quite sensitive to the solute dielectric constant, and this parameter should be carefully determined according to the characteristics of the protein/ligand binding interface. (3) Conformational entropy often show large fluctuations in MD trajectories, and a large number of snapshots are necessary to achieve stable predictions. Next, we evaluated the accuracy of the binding free energies calculated by three Generalized Born (GB) models. We found that the GB model developed by Onufriev and Case was the most successful model in ranking the binding affinities of the studied inhibitors. Finally, we evaluated the performance of MM/GBSA and MM/PBSA in predicting binding free energies. Our results showed that MM/PBSA performed better in calculating absolute, but not necessarily relative, binding free energies than MM/GBSA. Considering its computational efficiency, MM/GBSA can serve as a powerful tool in drug design, where correct ranking of inhibitors is often emphasized.     

Molecular modeling study of checkpoint kinase 1 inhibitors by multiple docking strategies and prime/MM-GBSA calculation

Developing chemicals that inhibit checkpoint kinase 1 (Chk1) is a promising adjuvant therapeutic to improve the efficacy and selectivity of DNA-targeting agents. Reliable prediction of binding-free energy and binding affinity of Chk1 inhibitors can provide a guide for rational drug design. In this study, multiple docking strategies and Prime/Molecular Mechanics Generalized Born Surface Area (Prime/MM-GBSA) calculation were applied to predict the binding mode and free energy for a series of benzoisoquinolinones as Chk1 inhibitors. Reliable docking results were obtained using induced-fit docking and quantum mechanics/molecular mechanics (QM/MM) docking, which showed superior performance on both ligand binding pose and docking score accuracy to the rigid-receptor docking. Then, the Prime/MM-GBSA method based on the docking complex was used to predict the binding-free energy. The combined use of QM/MM docking and Prime/MM-GBSA method could give a high correlation between the predicted binding-free energy and experimentally determined pIC(50) . The molecular docking combined with Prime/MM-GBSA simulation can not only be used to rapidly and accurately predict the binding-free energy of novel Chk1 inhibitors but also provide a novel strategy for lead discovery and optimization targeting Chk1.     

新型二氟甲基磷酸类酪氨酸蛋白磷酸酯酶1B抑制剂的分子动力学模拟和结合自由能计算

通过分子对接建立了一系列含二氟甲基磷酸基团(DFMP)或二氟甲基硫酸基团(DFMS)的抑制剂与酪氨酸蛋白磷酸酯酶1B(PTP1B)的相互作用模式, 并通过1 ns的分子动力学模拟和molecular mechanics/generalized Born surface area (MM/GBSA)方法计算了其结合自由能. 计算获得的结合自由能排序和抑制剂与靶酶间结合能力排序一致; 通过基于主方程的自由能计算方法, 获得了抑制剂与靶酶残基间相互作用的信息, 这些信息显示DFMP/DFMS基团的负电荷中心与PTP1B的221位精氨酸正电荷中心之间的静电相互作用强弱决定了此类抑制剂的活性, 进一步的分析还显示位于DFMP/DFMS基团中的氟原子或其他具有适当原子半径的氢键供体原子会增进此类抑制剂与PTP1B活性位点的结合能力.     

基于虚拟筛选和分子动力学模拟发现新型ULK1抑制剂

Unc-51样自噬激活激酶1（unc-51-like autophagy activating kinase 1,ULK1）作为自噬启动的重要调控因子,是肿瘤治疗的关键靶点之一。首先,以已知ULK1抑制剂为基础构建药效团模型,通过药效团模型筛选、分子对接以及分子力学广义波恩表面积（Molecular Mechanics/Generalized Born Surface Area,MM/GBSA）结合自由能计算等方法,对含有52万多个类药性小分子的数据库进行虚拟筛选,得到具有较高理论亲和力的化合物。随后,50ns的分子动力学模拟验证了蛋白质-配体复合物结合的稳定性,最后10ns的平均结合自由能的计算研究进一步验证了配体的结合能力。结果表明,6个化合物（F5258-0159、F3407-0428、F0529-1100、F0696-3531、F3222-5280、F6525-5596）具有骨架新颖、分子对接分数和结合自由能数值优异及与ULK1的结合状态稳定等特点,可以作为新型潜在的ULK1抑制剂用于肿瘤治疗的研究,也为新型ULK1抑制剂的设计和研发提供新的研究思路。     

Molecular dynamics simulations reveal the allosteric effect of F1174C resistance mutation to ceritinib in ALK-associated lung cancer

Anaplastic lymphoma kinase (ALK) has become as an important target for the treatment of various human cancers, especially non-small-cell lung cancer. A mutation, F1174C, suited in the C-terminal helix αC of ALK and distal from the small-molecule inhibitor ceritinib bound to the ATP-binding site, causes the emergence of drug resistance to ceritinib. However, the detailed mechanism for the allosteric effect of F1174C resistance mutation to ceritinib remains unclear. Here, molecular dynamics (MD) simulations and binding free energy calculations [Molecular Mechanics/Generalized Born Surface Area (MM/GBSA)] were carried out to explore the advent of drug resistance mutation in ALK. MD simulations observed that the exquisite aromatic-aromatic network formed by residues F1098, F1174, F1245, and F1271 in the wild-type ALK-ceritinib complex was disrupted by the F1174C mutation. The resulting mutation allosterically affected the conformational dynamic of P-loop and caused the upward movement of the P-loop from the ATP-binding site, thereby weakening the interaction between ceritinib and the P-loop. The subsequent MM/GBSA binding free energy calculations and decomposition analysis of binding free energy validated this prediction. This study provides mechanistic insight into the allosteric effect of F1174C resistance mutation to ceritinib in ALK and is expected to contribute to design the next-generation of ALK inhibitors.     

Computational and experimental studies of the interaction between phospho-peptides and the C-terminal domain of BRCA1

The C-terminal domain of BRCA1(BRCT) is involved in the DNA repair pathway by recognizing the pSXXF motif in interacting proteins. It has been reported that short peptides containing this motif bind to BRCA1(BRCT) in the micromolar range with high specificity. In this work, the binding of pSXXF peptides has been studied computationally and experimentally in order to characterize their interaction with BRCA1(BRCT). Elucidation of the contacts that drive the protein–ligand interaction is critical for the development of high affinity small-molecule BRCA1 inhibitors. Molecular dynamics (MD) simulations revealed the key role of threonine at the peptide P+2 position in providing structural rigidity to the ligand in the bound state. The mutation at P+1 had minor effects. Peptide extension at the N-terminal position with the naphthyl amino acid exhibited a modest increase in binding affinity, what could be explained by the dispersion interaction of the naphthyl side-chain with a hydrophobic patch. Three in silico end-point methods were considered for the calculation of binding free energy. The Molecular Mechanics Poisson–Boltzmann Surface Area and the Solvated Interaction Energy methods gave reasonable agreement with experimental data, exhibiting a Pearlman predictive index of 0.71 and 0.78, respectively. The MM-quantum mechanics-surface area method yielded improved results, which was characterized by a Pearlman index of 0.78. The correlation coefficients were 0.59, 0.61 and 0.69, respectively. The ability to apply a QM level of theory within an end-point binding free energy protocol may provide a way for a consistent improvement of accuracy in computer-aided drug design.     

Homology modeling and molecular dynamic simulation of UDP-N-acetylmuramoyl-l-alanine-d-glutamate ligase (MurD) from Mycobacterium tuberculosis H37Rv using in silico approach

The present study aimed to identify the prospective inhibitors of MurD, a cytoplasmic enzyme that catalyzes the addition of d-glutamate to the UDP-N-acetylmuramoyl-l-alanine nucleotide precursor in Mycobacterium tuberculosis (MTB), using virtual screening, docking studies, pharmacokinetic analysis, Molecular Dynamic (MD) simulation, and Molecular Mechanics Generalized Born and Surface Area (MM-GBSA) analyses. The three dimensional (3D) structure was determined based on the homology technique using a template from Streptococcus agalactiae. The modeled structure had three binding sites, namely; substrate binding site (Val18, Thr19, Asp39, Asp40, Gly75, Asn147, Gln171 and His192), the ATP binding site (Gly123, Lys124, Thr125, Thr126, Glu166, Asp283, and Arg314) and the glutamic acid binding site (Arg382, Ser463, and Tyr470). These residues mentioned above play a critical role in the catalytic activity of the enzyme, and their inhibition could serve as a stumbling block to the normal function of the enzyme. A total of 10,344 obtained from virtual screened of Zinc and PubChem databases. These compounds further screened for Lipinski rule of five, docking studies and pharmacokinetic analysis. Four compounds with good binding energies (ZINC11881196 = −10.33 kcal/mol, ZINC12247644 = −8.90 kcal/mol, ZINC14995379 =−8.42 kcal/mol, and PubChem6185 = −8.20 kcal/mol), better than the binding energies of the ATP (−2.31 kcal/mol) and the ligand with known IC₅₀, Aminothiazole (−7.11 kcal/mol) were selected for the MD simulation and MM-GBSA analyses. The result of the analyses showed that all the four ligands formed a stable complex and had the binding free energies better than the binding energy of ATP. Therefore, these ligands considered as suitable prospective inhibitors of the MurD after experimental validation.     

Crystal molecular dynamics simulations to speed up MM/PB(GB)SA evaluation of binding free energies of di-mannose deoxy analogs with P51G-m4-Cyanovirin-N

Complexes of two Cyanovirin-N (CVN) mutants, m4-CVN and P51G-m4-CVN, with deoxy di-mannose analogs were employed as models to generate conformational ensembles using explicit water Molecular Dynamics (MD) simulations in solution and in crystal environment. The results were utilized for evaluation of binding free energies with the molecular mechanics Poisson-Boltzmann (or Generalized Born) surface area, MM/PB(GB)SA, methods. The calculations provided the ranking of deoxy di-mannose ligands affinity in agreement with available qualitative experimental evidences. This confirms the importance of the hydrogen-bond network between di-mannose 3'- and 4'-hydroxyl groups and the protein binding site B(M) as a basis of the CVN activity as an effective HIV fusion inhibitor. Comparison of binding free energies averaged over snapshots from the solution and crystal simulations showed high promises in the use of the crystal matrix for acceleration of the conformational ensemble generation, the most time consuming step in MM/PB(GB)SA approach. Correlation between energy values based on solution versus crystal ensembles is 0.95 for both MM/PBSA and MM/GBSA methods.     

Structure-based prediction of free energy changes of binding of PTP1B inhibitors

The goals were (1) to understand the driving forces in the binding of small molecule inhibitors to the active site of PTP1B and (2) to develop a molecular mechanics-based empirical free energy function for compound potency prediction. A set of compounds with known activities was docked onto the active site. The related energy components and molecular surface areas were calculated. The bridging water molecules were identified and their contributions were considered. Linear relationships were explored between the above terms and the binding free energies of compounds derived based on experimental inhibition constants. We found that minimally three terms are required to give rise to a good correlation (0.86) with predictive power in five-group cross-validation test (q2 = 0.70). The dominant terms are the electrostatic energy and non-electrostatic energy stemming from the intra- and intermolecular interactions of solutes and from those of bridging water molecules in complexes.     

Discovery of a potent,selective protein tyrosine phosphatase 1B inhibitor using a linked-fragment strategy

Protein tyrosine phosphatase 1B (PTP1B) is an enzyme that downregulates the insulin receptor. Inhibition of PTP1B is expected to improve insulin action, and the design of small molecule PTP1B inhibitors to treat type II diabetes has received considerable attention. In this work, NMR-based screening identified a nonselective competitive inhibitor of PTP1B. A second site ligand was also identified by NMR-based screening and then linked to the catalytic site ligand by rational design. X-ray data confirmed that the inhibitor bound with the catalytic site in the native, "open" conformation. The final compound displayed excellent potency and good selectivity over many other phosphatases. The modular approach to drug design described in this work should be applicable for the design of potent and selective inhibitors of other therapeutically relevant protein tyrosine phosphatases.     

Molecular Dynamics Approach to Probe the Allosteric Inhibition of PTP1B by Chlorogenic and Cichoric Acid

Protein tyrosine phosphatase 1B (PTP1B), a major negative regulator of the insulin and leptin signaling pathway, is a potential target for therapeutic intervention against diabetes and obesity. The recent discovery of an allosteric site in PTP1B has created an alternate strategy in the development of PTP1B targeted therapy. The current study investigates the molecular interactions between the allosteric site of PTP1B with two caffeoyl derivatives, chlorogenic acid (CGA) and cichoric acid (CHA), using computational strategies. Molecular docking analysis with CGA and CHA at the allosteric site of PTP1B were performed and the resulting protein-ligand complexes used for molecular dynamics simulation studies for a time scale of 10 ns. Results show stable binding of CGA and CHA at the allosteric site of PTP1B. The flexibility of the WPD loop was observed to be constrained by CGA and CHA in the open (inactive), providing molecular mechanism of allosteric inhibition. The allosteric inhibition of CGA and CHA of PTP1B was shown to be favorable due to no restriction by the α-7 helix in the binding of CGA and CHA at the allosteric binding site. In conclusion, our results exhibit an inhibitory pattern of CGA and CHA against PTP1B through potent binding at the allosteric site.     

Structural Bases for Hesperetin Derivatives: Inhibition of Protein Tyrosine Phosphatase 1B,Kinetics Mechanism and Molecular Docking Study

In the present study, we investigated the structure-activity relationship of naturally occurring hesperetin derivatives, as well as the effects of their glycosylation on the inhibition of diabetes-related enzyme systems, protein tyrosine phosphatase 1B (PTP1B) and α-glycosidase. Among the tested hesperetin derivatives, hesperetin 5-O-glucoside, a single-glucose-containing flavanone glycoside, significantly inhibited PTP1B with an IC₅₀ value of 37.14 ± 0.07 µM. Hesperetin, which lacks a sugar molecule, was the weakest inhibitor compared to the reference compound, ursolic acid (IC₅₀ = 9.65 ± 0.01 µM). The most active flavanone hesperetin 5-O-glucoside suggested that the position of a sugar moiety at the C-5-position influences the PTP1B inhibition. It was observed that the ability to inhibit PTP1B is dependent on the nature, position, and number of sugar moieties in the flavonoid structure, as well as conjugation. In the kinetic study of PTP1B enzyme inhibition, hesperetin 5-O-glucoside led to mixed-type inhibition. Molecular docking studies revealed that hesperetin 5-O-glucoside had a higher binding affinity with key amino residues, suggesting that this molecule best fits the PTP1B allosteric site cavity. The data reported here support hesperetin 5-O-glucoside as a hit for the design of more potent and selective inhibitors against PTP1B in the search for a new anti-diabetic treatment.     

Structural and energetic insight into the interactions between the benzolactam inhibitors and tumor marker HSP90α

The heat shock protein 90α (HSP90α) provides a promising molecular target for cancer therapy. A series of novel benzolactam inhibitors exhibited distinct inhibitory activity for HSP90α. However, the structural basis for the impact of distinct R₁ substituent groups of nine benzolactam inhibitors on HSP90α binding affinities remains unknown. In this study, we carried out molecular docking, molecular dynamics (MD) simulations, and molecular mechanics and generalized Born/surface area (MM–GBSA) binding free energy calculations to address the differences. Molecular docking studies indicated that all nine compounds presented one conformation in the ATP-binding site of HSP90α N-terminal domain. MD simulations and subsequent MM–GBSA calculations revealed that the hydrophobic interactions between all compounds and HSP90α contributed the most to the binding affinity and a good linear correlation was obtained between the calculated and the experimental binding free energies (R = 0.88). The per residue decomposition revealed that the most remarkable differences of residue contributions were found in the residues Ala55, Ile96, and Leu107 defining a hydrophobic pocket for the R₁ group, consistent with the analysis of binding modes. This study may be helpful for the future design of novel HSP90α inhibitors.     

酪氨酸蛋白磷酸酯酶1B抑制剂的分子对接和三维定量构效关系研究

用一种柔性分子对接方法(FlexX)将12个2-草酰胺苯甲酸类抑制剂和酪氨酸蛋白磷酸酯酶(PTP1B)活性口袋进行分子对接,对接程序预测的抑制剂和酶之间的相互作用能与抑制活性之间有很好的相关性(非线性相关系数R²达0.859),这说明对接结果可以比较准确地预测抑制剂和PTP1B之间的结合模式.然后,将33个同类抑制剂的骨架叠合在分子对接预测的结合构象上,用比较分子力场分析方法(CoMFA)对其进行三维定量活性构效关系研究,得到的CoMFA模型具有很好的统计相关性(交互验证回归系数q²为0.650),并可以准确地预测测试集6个化合物的活性(平均标准偏差为0.177).同时,由CoMFA模型得出的抑制剂改造信息与用FlexX预测的结合模式是一致的,进一步证明我们预测的结合模式是正确的.为研究这类抑制剂和PTP1B的结合模式及对抑制剂进行结构改造提供了信息.     

Binding affinities of factor Xa inhibitors estimated by thermodynamic integration and MM/GBSA

We present free energy estimates of nine 3-amidinobenzyl-1H-indole-2-carboxamide inhibitors of factor Xa. Using alchemical thermodynamic integration (TI) calculations, we estimate the difference in binding free energies with high accuracy and precision, except for mutations involving one of the amidinobenzyl rings. Crystal studies show that the inhibitors may bind in two distinct conformations, and using TI, we show that the two conformations give a similar binding affinity. Furthermore, we show that we can reduce the computational demand, while still retaining a high accuracy and precision, by using fewer integration points and shorter protein-ligand simulations. Finally, we have compared the TI results to those obtained with the simpler MM/GBSA method (molecular-mechanics with generalized Born surface-area solvation). MM/GBSA gives better results for the mutations that involve a change of net charge, but if a precision similar to that of the TI method is required, the MM/GBSA method is actually slightly more expensive. Thus, we have shown that TI could be a valuable tool in drug design.     

Selectivity of neutral/weakly basic P1 group inhibitors of thrombin and trypsin by a molecular dynamics study

Molecular dynamics (MD) simulations followed by molecular mechanics generalized Born surface area (MM-GBSA) analyses have been carried out to study the selectivity of two neutral and weakly basic P1 group inhibitors (177 and CDA) to thrombin and trypsin. Detailed binding free energies between these inhibitors and individual protein residues are calculated by using a per-residue basis decomposition method. The analysis of the detailed interaction energies provides insight on the protein-inhibitor-binding mechanism and helps to elucidate the basis for achieving selectivity through interpretation of the structural and energetic results from the simulations. The study shows that the dominant factor of selectivity for both inhibitors is van der Waals energy, which suggests better shape complementarity and packing with thrombin. Nonpolar solvation free energy and total entropy contribution are also in favor of selectivity, but the contributions are much smaller. Binding mode and structural analysis show that 177 binds to thrombin and trypsin in a similar binding mode. In contrast, the CDA binds to thrombin and trypsin in very different modes.     

Synthesis and crystal structure of a chiral lactam and three amino alcohols as potential protein tyrosine phosphates 1B inhibitors

Chiral lactam 2 and three chiral β-amino alcohols 3–5 have been synthesized and characterized by spectroscopic techniques. Regioselective ring opening reaction of chiral styrene oxide by an amine nucleophile was confirmed by X-ray diffraction data. Ligand 2–4 crystallizes in the tetragonal, orthorhombic and tetragonal crystal lattice system respectively. Ligands 2–6 have been used as potential inhibitors for protein tyrosine phosphatase 1B enzyme (PTP1B). The potential inhibitor effect of these molecules to the target protein was investigated by Dock and molecular dynamics calculations. Dock score analysis and Lipinski parameters suggested that ligands 1, 2, 4–6 are potential inhibitors towards PTP1B, thus indicating that the residues Arg24, Arg254 and Met258, Asp29 in the second active site of PTP1B are essential for the high selectivity of inhibitors. The results indicate that the polar hydrogen bonding interacts with Asp29, Gln102, and the amino acid residues of PTP1B are responsible for governing inhibitory potency of ligands 1–6.