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.     

Computational insights into allosteric interaction between benzoazepin-2-ones and lung cancer-associated PDK1: Implications for activator design

3-Phosphoinositide-dependent kinase-1 (PDK1) plays a key role in the regulation of physiological processes and its catalytic activity is tightly regulated by allosteric modulators which bind to the PDK1 Interacting Fragment (PIF) pocket. However, details on the allosteric modulators regulation of the PDK1 catalytic activity remain elusive. Here, molecular docking and molecular dynamics (MD) simulations were performed to investigate the allosteric regulation of PDK1 induced by one of the benzoazepin-2-ones, the most potent compound 17 (BAZ2O). Molecular docking and MD simulation revealed that BAZ2O was located in the PIF pocket formed by residues from β4 and β5 sheets and helices αB and αC. BAZ2O formed a hydrogen bond with Arg131 and participated in hydrophobic interactions with Ile119, Thr148, Gln150, Leu155 and Phe157. Further comparative analyses of PDK1 in its apo and BAZ2O-bound states unveiled that BAZ2O promoted the structural coupling between the important catalytic domains of PDK1, including the activation loop and the helices αB and αC, thereby stabilizing the PDK1 conformation for catalysis. Understanding the allosteric interaction of PDK1 with small molecules provides a potentially valuable possibility of designing more potent allosteric modulators with therapeutic implications for lung cancer.     

Identification of Two Distinct Sites for Antagonist and Biased Agonist Binding to the Human Chemokine Receptor CXCR3

The chemokine receptor CXCR3 is a G protein‐coupled receptor that conveys extracellular signals into cells by changing its conformation upon ligand binding. We previously hypothesized that small‐molecule allosteric CXCR3‐agonists do not bind to the same allosteric binding pocket as 8‐azaquinazolinone‐based negative allosteric modulators. We have now performed molecular‐dynamics (MD) simulations with metadynamics enhanced sampling on the CXCR3 system to refine structures and binding modes and to predict the CXCR3‐binding affinities of the biased allosteric agonist FAUC1036 and the negative allosteric modulator RAMX3. We have identified two distinct binding sites; a “shallow” and a second “deeper” pocket to which the biased allosteric agonist FAUC1036 and negative allosteric modulator RAMX3 bind, respectively.     

Apicidin选择性抑制ClassⅠ HDACs分子动力学研究

采用模拟方法研究组蛋白脱乙酰酶抑制剂(Apicidin)选择性抑制组蛋白去乙酰化酶(Histone deacetylases, HDACs)中的HDAC1和HDAC8. 通过HDAC8晶体结构同源模建HDAC1三维结构模型, 将Apicidin分别与HDAC1和HDAC8对接并进行分子动力学模拟, 结果表明, HDAC1活性口袋入口处的Arg270是Apicidin-HDAC1形成稳定结构的重要因素; HDAC1中Tyr303及His178与Apicidin形成2个持续存在的氢键, 而在HDAC8中未发现, 这是Apicidin选择性抑制HDAC1高于HDAC8的另一重要原因.     

Structure-based QSAR study on differential inhibition of human prostaglandin endoperoxide H synthase-2 (COX-2) by nonsteroidal anti-inflammatory drugs

The prostaglandin-endoperoxide H synthase-1 (PGHS-1) and prostaglandin-endoperoxide H synthase-2 (PGHS-2) are the targets of nonsteroidal anti-inflammatory drugs (NSAIDs).It appears that the high degree of selectivity for inhibition of PGHS-2 shown by certain compounds is the result of two mechanisms (time-dependent, time-independent inhibition), by which they interact with each isoform. Molecular models of the complexes formed by indomethacin, sulindac, fenamates, 2-phenylpropionic acids and selective cyclooxygenase-2 (COX-2) inhibitors with the cyclooxygenase active site of human PGHS-2 have been built, paying particular attention to water molecules that participate in the hydrogen-bonding network at the polar active site entrance. The stability of the complexes has been assessed by molecular dynamics simulations and interaction energy decomposition analysis, and their biological significance has been discussed in light of available X-ray crystallographic and kinetic results. The selective PGHS-2 inhibitors exploit the extra space of a side-pocket in the active site of PGHS-2 that is not found in PGHS-1. The results suggest that active site hydration together with residues Tyr³⁵⁵, Glu⁵²⁴, Arg¹²⁰ and Arg⁵¹³ are crucial to understand the time-dependent inhibition mechanism. A marked relationship between the isoform selectivity and tightly interactions with residues into the side pocket bordered by Val⁵²³ is also found.     

表皮生长因子受体和抑制剂之间分子对接的研究

研究了表皮生长因子受体(EGFR)和4-苯胺喹唑啉类抑制剂之间的相互作用模式，表皮生长因子受体的三维结构通过同源蛋白模建的方法得到，而抑制剂和靶酶结合复合物结构则通过分子力学和分子动力学结合的方法计算得到。从模拟结果得到的抑制剂和靶酶之间的相互作用模式表明范德华相互作用、疏水相互作用以及氢键相互作用对抑制剂的活性都有重要的影响，抑制剂的苯胺部分位于活性口袋的底部，能够与受体残基的非极性侧链产生很强的范德华和疏水相互作用，抑制剂双环上的取代基团也能和活性口袋外部的部分残基形成一定的范德华和疏水性相互作用，而抑制剂喹唑啉环上的氮原子能和周围的残基形成较强的氢键相互作用，对抑制剂的活性有较大的影响，计算得到抑制剂和靶酶之间的非键相互作用能以及抑制剂和靶酶之间的相互作用信息能够很好地解释抑制剂活性和结构的关系，为全新抑制剂的设计提供了重要的结构信息。     

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

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

A Comprehensive Computational Screening of Phytochemicals Derived from Saudi Medicinal Plants against Human CC Chemokine Receptor 7 to Identify Potential Anti-Cancer Therapeutics

Homeostatic trafficking of immune cells by CC chemokine receptor 7 (CCR7) keeps immune responses and tolerance in a balance. The involvement of this protein in lymph node metastasis in cancer marks CCR7 as a penitential drug target. Using the crystal structure of CCR7, herein, a comprehensive virtual screening study is presented to filter novel strong CCR7 binding phytochemicals from Saudi medicinal plants that have a higher binding affinity for the intracellular allosteric binding pocket. By doing so, three small natural molecules named as Hit-1 (1,8,10-trihydroxy-3-methoxy-6-methylanthracen-9(4H)-one), Hit-2 (4-(3,4-dimethoxybenzyl)-3-(4-hydroxy-3-methoxybenzyl)dihydrofuran-2(3H)-one), and Hit-3 (10-methyl-12,13-dihydro-[1,2]dioxolo[3,4,5-de]furo[3,2-g]isochromeno[4,3-b]chromen-8-ol) are predicted showing strong binding potential for the CC chemokine receptor 7 allosteric pocket. During molecular dynamics simulations, the compounds were observed in the formation of several chemical bonding of short bond distances. Additionally, the molecules remained in strong contact with the active pocket residues and experienced small conformation changes that seemed to be mediated by the CCR7 loops to properly engage the ligands. Two types of binding energy methods (MM/GBPBSA and WaterSwap) were additionally applied to further validate docking and simulation findings. Both analyses complement the good affinity of compounds for CCR7, the electrostatic and van der Waals energies being the most dominant in intermolecular interactions. The active pocket residue’s role in compounds binding was further evaluated via alanine scanning, which highlighted their importance in natural compounds binding. Additionally, the compounds fulfilled all drug-like rules: Lipinski, Ghose, Veber, Egan, and Muegge passed many safety parameters, making them excellent anti-cancer candidates for experimental testing.     

Computational Characterization of Binding of Small Molecule Inhibitors to HIV‐1 gp41

Developing orally available small molecule inhibitors of HIV‐1 fusion has attracted significant interest over many years. Frey had recently reported several synthetic compounds which are experimentally shown to inhibit cell‐cell fusion in the low micromolar range. We carried out computational study to help identify possible binding modes by docking these compounds onto the hydrophobic pocket on gp41 and to characterize structures of binding complexes. The detailed gp41‐molecule binding interactions and free energies of binding are obtained through molecular dynamics simulation and MM‐PBSA calculation. Specific molecular interactions in the gp41‐inhibitor complexes are identified. The present computational study complements the corresponding experimental investigation and helps establish a good starting point for further refinement of small molecular gp41 inhibitors.     

A Molecular Tool Targeting the Base-Flipping Activity of Human UHRF1

During DNA replication, ubiquitin-like, containing PHD and RING fingers domains 1 (UHRF1) plays key roles in the inheritance of methylation patterns to daughter strands by recognizing through its SET and RING-associated domain (SRA) the methylated CpGs and recruiting DNA methyltransferase 1 (DNMT1). Herein, our goal is to identify UHRF1 inhibitors targeting the 5′-methylcytosine (5mC) binding pocket of the SRA domain to prevent the recognition and flipping of 5mC and determine the molecular and cellular consequences of this inhibition. For this, we used a multidisciplinary strategy combining virtual screening and molecular modeling with biophysical assays in solution and cells. We identified an anthraquinone compound able to bind to the 5mC binding pocket and inhibit the base-flipping process in the low micromolar range. We also showed in cells that this hit impaired the UHRF1/DNMT1 interaction and decreased the overall methylation of DNA, highlighting the critical role of base flipping for DNMT1 recruitment and providing the first proof of concept of the druggability of the 5mC binding pocket. The selected anthraquinone appears thus as a key tool to investigate the role of UHRF1 in the inheritance of methylation patterns, as well as a starting point for hit-to-lead optimizations.     

Recognition of protein allosteric states and residues: Machine learning approaches

Allostery is a process by which proteins transmit the effect of perturbation at one site to a distal functional site upon certain perturbation. As an intrinsically global effect of protein dynamics, it is difficult to associate protein allostery with individual residues, hindering effective selection of key residues for mutagenesis studies. The machine learning models including decision tree (DT) and artificial neural network (ANN) models were applied to develop classification model for a cell signaling allosteric protein with two states showing extremely similar tertiary structures in both crystallographic structures and molecular dynamics simulations. Both DT and ANN models were developed with 75% and 80% of predicting accuracy, respectively. Good agreement between machine learning models and previous experimental as well as computational studies of the same protein validates this approach as an alternative way to analyze protein dynamics simulations and allostery. In addition, the difference of distributions of key features in two allosteric states also underlies the population shift hypothesis of dynamics‐driven allostery model. © 2018 Wiley Periodicals, Inc.     

激酶ABL与肉豆蔻酰位点小分子作用机理的分子动力学模拟及自由能计算

采用分子动力学方法研究激酶ABL 与ATP 位点小分子imatinib、P16 及变构位点小分子STJ、MS7、MS9、3YY、MYR等的结合, 并用GBSA (generalized Born surface area)方法将结合自由能分解到各残基. 自由能计算结果表明, 小分子STJ、MS7、MS9 有利于imatinib 与ABL 结合; 小分子STJ、MS7、MS9 与激酶ABL的结合自由能接近, 绝对值均大于ABL 与3YY、MYR 的结合自由能. 能量分解表明, ABL 残基ILE502、VAL506、LEU510与STJ和MYR的相互作用是αI 螺旋处于弯曲状态的重要原因. 模拟过程中ABL肉豆蔻酰口袋残基均方根偏差(RMSD)变化值表明, STJ等小分子抑制剂与ABL结合后降低了肉豆蔻酰口袋残基的柔性.     

Synthesis,Spectral, Thermal and Biological Studies of Mn(II), Co(II), Ni(II) and Cu(II) Complexes with 1‐(((5‐Mercapto‐1H‐1,2,4‐triazol‐3‐yl)imino)‐methyl)naphthalene‐2‐ol

Mn(II), Co(II), Ni(II) and Cu(II) complexes of 5‐mercapto‐1,2,4‐triazol‐3‐imine‐2′‐hydroxynaphthaline have been synthesized and characterized by elemental analysis, IR, ¹H NMR, EI‐mass, UV‐Vis, and ESR (electron spin resonance) spectra, molar conductance, magnetic moment measurements, DC conductivity and thermogravimetric analysis. IR spectra confirm that the ligand molecule existed in both thione and thiole forms. The molar conductance values indicate the complexes are nonelectrolyte. The magnetic moment values of the complexes display paramagnetic behavior. All studies confirm the formation of an octahedral geometry for complex 1 and the other complexes have tetrahedral geometrical structures. The structures of the complexes have also been theoretically studied by using the molecular mechanic calculations by the hyperchem. 8.03 molecular modeling program which confirm the proposed structures. The Schiff‐base ligand and its metal complexes have also been screened for their antimicrobial activities.     

Binding Modes of Small-Molecule Inhibitors to the EED Pocket of PRC2

Polycomb Polycomb repressive complex 2 (PRC2) plays a key role in silencing epigenetic gene through trimethylation of lysine 27 on histone 3 (H3K27). Dysregulations of PRC2 caused by overexpression and mutations of the core subunits of PRC2 have been implicated in many cancers. The core subunits EZH1/2 are histone-lysine N-methyltransferases that function as the enzymatic component of PRC2. While the core subunit EED is a scaffolding protein to support EZH1/2 and binds JARID2K116me3/H3K27me3 to enhance the enzymatic activity of PRC2 through allosteric activation. Recently, several small molecules that compete with JARI2K116me3 and H3K27me3 have been reported. These molecules selectively bind to the JARID2K116me3/H3K27me3-binding pocket of EED, thereby preventing the allosteric regulation of PRC2. These first-in-class PRC2 inhibitors show robust suppression in DLBCL cell lines, demonstrating anticancer drugs that target the EED subunit of PRC2 are viable. In this study, we used the recently developed MM/GBSA_IE and the alanine scanning method to analyze the hot spots in EED/inhibitor interactions. The analysis of these hot and warm spots helps us to understand the fundamental differences between inhibitors. Our results give a quantitative explanation on why the binding affinities of EED/A-395 interactions are stronger than that of EED/EED226 while their binding modes are similar and provide valuable insights for rational design of novel EED inhibitors.     

Unraveling the Compositional and Molecular Features Involved in Lysozyme-Benzothiazole Derivative Interactions

In this work we present a computational analysis together with experimental studies, focusing on the interaction between a benzothiazole (BTS) and lysozyme. Results obtained from isothermal titration calorimetry, UV-vis, and fluorescence were contrasted and complemented with molecular docking and machine learning techniques. The free energy values obtained both experimentally and theoretically showed excellent similarity. Calorimetry, UV-vis, and 3D/2D-lig-plot analysis revealed that the most relevant interactions between BTS and lysozyme are based on a predominance of aromatic, hydrophobic Van der Waals interactions, mainly aromatic edge-to-face (T-shaped) π-π stacking interactions between the benzene ring belonging to the 2-(methylthio)-benzothiazole moiety of BTS and the aromatic amino acid residue TRP108 of the lysozyme receptor. Next, conventional hydrogen bonding interactions contribute to the stability of the BTS-lysozyme coupling complex. In addition, mechanistic approaches performed using elastic network models revealed that the BTS ligand theoretically induces propagation of allosteric signals, suggesting non-physiological conformational flexing in large blocks of lysozyme affecting α-helices. Likewise, the BTS ligand interacts directly with allosteric residues, inducing perturbations in the conformational dynamics expressed as a moderate conformational softening in the α-helices H1, H2, and their corresponding β-loop in the lysozyme receptor, in contrast to the unbound state of lysozyme.     

Synthetic hosts for molecular recognition of ureas

Four hosts (7-10) containing 2,6-bisamidopyridine- and 2,5-bisamidopyrrole-bearing pyridyl or 1,8-naphthyridyl groups have been prepared and their structures studied by a combination of multinuclear NMR spectroscopy and X-ray crystallography. Their behavior in molecular recognition of urea derivatives, including (+)-biotin methyl ester, has been approached by molecular modeling (Monte Carlo conformational search, AMBER force field). The minimum energy values for the complexes are correlated with the experimental binding energies determined by means of (1)H NMR titrations.     

Relationship between chiroptical properties, structural changes and interactions in enzymes: a computational study on beta-lactamases from class A

Enzyme functions such as catalysis, binding and regulation are directly related to a variety of conformational changes. A sensitive and useful method for their investigation is circular dichroism (CD) and a rotational strength (R) is its fundamental characteristic. In this study, we show how the sensitivity of the mechanisms of rotational strengths to important conformational changes depends on the chromophore environment in two beta-lactamases from class A (from Escherichia coli and B. licheniformis). Rotational strengths have been calculated using the matrix method and including the effects of local environment (LE). X-ray structures (of protein components) of several enzyme-ligand complexes from the catalytic cycle of the TEM-1 enzyme and for both crystallographic monomers of the enzyme from B. licheniformis were used. An analysis of the relative degree of perturbation of the rotational strengths upon local interactions is performed.     

Two-track virtual screening approach to identify both competitive and allosteric inhibitors of human small C-terminal domain phosphatase 1

Despite a wealth of persuasive evidence for the involvement of human small C-terminal domain phosphatase 1 (Scp1) in the impairment of neuronal differentiation and in Huntington’s disease, small-molecule inhibitors of Scp1 have been rarely reported so far. This study aims to the discovery of both competitive and allosteric Scp1 inhibitors through the two-track virtual screening procedure. By virtue of the improvement of the scoring function by implementing a new molecular solvation energy term and by reoptimizing the atomic charges for the active-site Mg²⁺ ion cluster, we have been able to identify three allosteric and five competitive Scp1 inhibitors with low-micromolar inhibitory activity. Consistent with the results of kinetic studies on the inhibitory mechanisms, the allosteric inhibitors appear to be accommodated in the peripheral binding pocket through the hydrophobic interactions with the nonpolar residues whereas the competitive ones bind tightly in the active site with a direct coordination to the central Mg²⁺ ion. Some structural modifications to improve the biochemical potency of the newly identified inhibitors are proposed based on the binding modes estimated with docking simulations.     

Structure-function relationship of amino acid-[2]rotaxanes

Synthetic methodology was developed to construct amino acid-[2]rotaxanes that have phenylalanine and 3,5-di-tert-butylbenzene as blocking groups and dibenzo-24-crown-8, derivatized with either N-acetylargininyl or a carboxylic group, as the ring. A relative measure of the intramolecular interaction energies between the functional groups in DMSO/water mixtures is obtained by comparing their pK(a) values. Rotaxane structures were investigated through 2D NMR analysis and molecular dynamics simulations. Association constants for complexes of amino acids and rotaxanes in various protonation states were determined in a variety of solvent systems by (1)H NMR analysis. The unique intracomponent interactions that exist in the rotaxanes and their ability to act as artificial receptors are discussed.