Molecular dynamics simulations of interaction between protein-tyrosine phosphatase 1B and a bidentate inhibitor

An in silico analysis of the interaction of marine sponge-derived bioactive compounds with type 2 diabetes mellitus targets DPP-4 and PTP1B

Type 2 diabetes is a medical condition involving elevated blood glucose levels resulting from impaired or improper insulin utilization. As the number of type 2 diabetes cases increases each year, there is an urgent need to develop novel drugs having new targets and/or complementing existing therapeutic protocols. In this regard, marine sponge-derived compounds hold great potential due to their potent biological activity and structural diversity. In this study, a small library of 50 marine sponge-derived compounds were examined for their activity towards type 2 diabetes targets, namely dipeptidyl peptidase-4 (DPP-4) and protein tyrosine phosphatase 1B (PTP1B). The compounds were first subjected to molecular docking on protein models based on their respective co-crystal structures to assess binding free energies (BFE) and conformations. Clustering analysis yielded BFE that ranged from 24.54 kcal/mol to -9.97 kcal/mol for DPP-4, and from -4.98 kcal/mol to -8.67 kcal/mol for PTP1B. Interaction analysis on the top ten compounds with the most negative BFE towards each protein target showed similar intermolecular interactions and key interacting residues as in the previously solved co-crystal structure. These compounds were subjected to absorption, distribution, metabolism, excretion, and toxicity (ADMET) profiling to characterize drug-likeness and combining the results from these analyses, (S)-6'-debromohamacanthin B was identified as a potential multi-target inhibitor of DPP-4 and PTP1B, having favorable protein interaction, no Lipinski violations, good gastrointestinal (GI) tract absorption, blood-brain barrier (BBB) penetration, and no predicted toxicity. Finally, the interaction of (S)-6'-debromohamacanthin B with the two proteins was validated using molecular dynamics simulations over 100 ns through RMSD, radius of gyration, PCA, and molecular mechanics Poisson-Boltzmann surface area (MMPBSA) confirming favorable interactions with the respective proteins.Communicated by Ramaswamy H. Sarma.

A Novel PTP1B Inhibitor-Phosphate of Polymannuronic Acid Ameliorates Insulin Resistance by Regulating IRS-1/Akt Signaling

Protein tyrosine phosphatase 1B (PTP1B) is a critical negative modulator of insulin signaling and has attracted considerable attention in treating type 2 diabetes mellitus (T2DM). Low-molecular-weight polymannuronic acid phosphate (LPMP) was found to be a selective PTP1B inhibitor with an IC₅₀ of 1.02 ± 0.17 μM. Cellular glucose consumption was significantly elevated in insulin-resistant HepG2 cells after LPMP treatment. LPMP could alleviate oxidative stress and endoplasmic reticulum stress, which are associated with the development of insulin resistance. Western blot and polymerase chain reaction (PCR) analysis demonstrated that LPMP could enhance insulin sensitivity through the PTP1B/IRS/Akt transduction pathway. Furthermore, animal study confirmed that LPMP could decrease blood glucose, alleviate insulin resistance, and exert hepatoprotective effects in diabetic mice. Taken together, LPMP can effectively inhibit insulin resistance and has high potential as an anti-diabetic drug candidate to be further developed.

Application of machine learning methods for the development of antidiabetic drugs

Diabetes is a chronic non-communicable disease caused by several different routes, which has attracted increasing attention. In order to speed up the development of new selective drugs, machine learning (ML) technology has been applied in the process of diabetes drug development, which opens up a new blueprint for drug design. This review provides a comprehensive portrayal of the application of ML in antidiabetic drug use.

Conformation and dynamics of the C-terminal region in human phosphoglycerate mutase 1

Phosphoglycerate mutase 1 (PGAM1), an important enzyme in glycolysis, is overexpressed in a number of human cancers, thus has been proposed as a promising metabolic target for cancer treatments. The C-terminal portion of the available crystal structures of PGAM1 and its homologous proteins is partially disordered, as evidenced by weak electron density. In this study, we identified the conformational behavior of the C-terminal region of PGAM1 as well as its role during the catalytic cycle. Using the PONDR-FIT server, we demonstrated that the C-terminal region was intrinsically disordered. We applied the Monte Carlo (MC) method to explore the conformational space of the C-terminus and conducted a series of explicit-solvent molecular dynamics (MD) simulations, and revealed that the C-terminal region is inherently dynamic; large-scale conformational changes in the C-terminal segment led to the structural transition of PGAM1 from the closed state to the open state. Furthermore, the C-terminal segment influenced 2,3-bisphosphoglycerate (2,3-BPG) binding. The proposed swing model illustrated a critical role of the C-terminus in the catalytic cycle through the conformational changes. In conclusion, the C-terminal region induces large movements of PGAM1 from the closed state to the open state and influences cofactor binding during the catalytic cycle. This report describes the dynamic features of the C-terminal region in detail and should aid in design of novel and efficient inhibitors of PGAM1. A swing mechanism of the C-terminal region is proposed, to facilitate further studies of the catalytic mechanism and the physiological functions of its homologues.

BPN, a marine-derived PTP1B inhibitor, activates insulin signaling and improves insulin resistance in C2C12 myotubes

Insulin resistance is a key feature of type 2 diabetes mellitus (T2DM) and is characterized by defects in insulin signaling. Protein tyrosine phosphatase 1B (PTP1B) is a major negative regulator of insulin signaling cascade and has attracted intensive investigation in recent T2DM therapy study. BPN, a marine-derived bromophenol compound, was isolated from the red alga Rhodomela confervoides. This study investigated the effects of BPN on the insulin signaling pathway in insulin-resistant C2C12 myotubes by inhibiting PTP1B. Molecular docking study and analysis of small- molecule interaction with PTP1B all showed BPN inhibited PTP1B activity via binding to the catalytic site through hydrogen bonds. We then found that BPN permeated into C2C12 myotubes, on the one hand, activated insulin signaling in an insulin-independent manner in C2C12 cells; on the other hand, ameliorated palmitate-induced insulin resistance through augmenting insulin sensitivity. Moreover, our studies also showed that PTP1B inhibition by BPN increased glucose uptake in normal and insulin-resistant C2C12 myotubes through glucose transporter 4 (GLUT4) translocation. Taken together, BPN activates insulin signaling and alleviates insulin resistance and represents a potential candidate for further development as an antidiabetic agent.

Characterization of Protein Tyrosine Phosphatase 1B Inhibition by Chlorogenic Acid and Cichoric Acid

Protein tyrosine phosphatase 1B (PTP1B) is a known regulator of the insulin and leptin signaling pathways and is an active target for the design of inhibitors for the treatment of type II diabetes and obesity. Recently, cichoric acid (CHA) and chlorogenic acid (CGA) were predicted by docking methods to be allosteric inhibitors that bind distal to the active site. However, using a combination of steady-state inhibition kinetics, solution nuclear magnetic resonance experiments, and molecular dynamics simulations, we show that CHA is a competitive inhibitor that binds in the active site of PTP1B. CGA, while a noncompetitive inhibitor, binds in the second aryl phosphate binding site, rather than the predicted benzfuran binding pocket. The molecular dynamics simulations of the apo enzyme and cysteine-phosphoryl intermediate states with and without bound CGA suggest CGA binding inhibits PTP1B by altering hydrogen bonding patterns at the active site. This study provides a mechanistic understanding of the allosteric inhibition of PTP1B.

Insight into the interactive residues between two domains of human somatic Angiotensin-converting enzyme and Angiotensin II by MM-PBSA calculation and steered molecular dynamics simulation

Angiotensin-converting enzyme (ACE), a membrane-bound zinc metallopeptidase, catalyzes the formation of Angiotensin-II (AngII) and the deactivation of bradykinin in the renin-angiotensin-aldosterone and kallikrein-kinin systems. As a hydrolysis product of ACE, AngII is regarded as an inhibitor and displays stronger competitive inhibition in the C-domain than the N-domain of ACE. However, the AngII binding differences between the two domains and the mechanisms behind AngII dissociation from the C-domain are rarely explored. In this work, molecular docking, Molecular Mechanics/Poisson-Boltzmann Surface Area calculation, and steered molecular dynamics (SMD) are applied to explore the structures and interactions in the binding or unbinding of AngII with the two domains of human somatic ACE. Calculated free energy values suggest that the C-domain-AngII complex is more stable than the N-domain-AngII complex, consistent with available experimental data. SMD simulation results imply that electrostatic interaction is dominant in the dissociation of AngII from the C-domain. Moreover, Gln106, Asp121, Glu123, and Tyr213 may be the key residues in the unbinding pathway of AngII. The simulation results in our work provide insights into the interactions between the two domains of ACE and its natural peptide inhibitor AngII at a molecular level. Moreover, the results provide theoretical clues for the design of new inhibitors.

Evolutional and functional analysis of a serine protease in Spodoptera litura

Spodoptera litura is a threatening agricultural insect in tropical and subtropical areas and accounts for tremendous annual crop losses. As seen in virtually all insect species, serine proteases (SPs) are crucial to S. litura. The expression pattern of SPs from the midgut of S. litura was studied through expressed sequence tags (ESTs) analysis. One of SP (SlSP1) was chosen for detailed study, because the expression of the gene was midgut and larvae specific. SlSP1 was conducted as a model of its evolution, structure, and potential binding activity with corresponding substrates. SlSP1 is composed of 255 amino acids including a signal peptide at N-terminal followed by a putative activation peptide and the mature protein along with five putative phosphorylation sites, three disulphide bridges, and two N-glycosylation positions. At least nine conserved motifs were obtained in multiple sequence alignments. Some conserved residues, such as the catalytic triad His84, Asp127, and Ser229 as well as six cysteines at position 66, 82, 194, 211, 223, and 247, were examined. After homology modeling and molecular dynamics simulation, the resultant three-dimensional (3D) structure of SlSP1 was docked with the substrates 2PTC-Arg and 2PTC-Lys, respectively. Molecular Mechanic/Poisson-Boltzmann surface area analysis was applied to anticipate optimal binding mode and crucial active sites of this enzyme. The residues Trp28, Gly187, Aso188, Arg249, Ile250, Lys246, and Lys278 are crucial for the substrate binding and molecule process. This information can be used in logical design of SPs inhibitors. New inhibitors may be a basis for development of a new pest control technology.

Protein tyrosine phosphatases in cancer: friends and foes!

Tyrosine phosphorylation of proteins serves as an exquisite switch in controlling several key oncogenic signaling pathways involved in cell proliferation, apoptosis, migration, and invasion. Since protein tyrosine phosphatases (PTPs) counteract protein kinases by removing phosphate moieties on target proteins, one may intuitively think that PTPs would act as tumor suppressors. Indeed, one of the most described PTPs, namely, the phosphatase and tensin homolog (PTEN), is a tumor suppressor. However, a growing body of evidence suggests that PTPs can also function as potent oncoproteins. In this chapter, we provide a broad historical overview of the PTPs, their mechanism of action, and posttranslational modifications. Then, we focus on the dual properties of classical PTPs (receptor and nonreceptor) and dual-specificity phosphatases in cancer and summarize the current knowledge of the signaling pathways regulated by key PTPs in human cancer. In conclusion, we present our perspective on the potential of these PTPs to serve as therapeutic targets in cancer.

Molecular dynamics simulation of the interaction between protein tyrosine phosphatase 1B and aryl diketoacid derivatives

The protein tyrosine phosphatase 1B (PTP-1B) is acknowledged as an outstanding therapeutic target for the treatment of diabetes, obesity and cancer. In this work, six aryl diketoacid compounds have been studied on the basis of molecular dynamics simulations. Hydrogen bonds, binding energies and conformation changes of the WPD loop have been analyzed. The results indicated that their activation model falls into two parts: the target region of the monomeric aryl diketoacid compounds is the active site, whereas the target region of the dimeric aryl diketoacid compounds is the WPD loop or the R loop. The van der Waals interactions exhibit stronger effects than the short-range electrostatic interactions. The van der Waals interaction energy and the IC50 values exhibit an approximately exponential relationship. Furthermore, the van der Waals interactions cooperate with the hydrogen bond interactions. This study provides a more thorough understanding of the PTP-1B inhibitor binding processes.

Effect of a novel proteoglycan PTP1B inhibitor from Ganoderma lucidum on the amelioration of hyperglycaemia and dyslipidaemia in db/db mice

Protein tyrosine phosphatase 1B (PTP1B) is implicated in the negative regulation of the insulin signalling pathway by dephosphorylating the insulin receptor (IR) and IR substrates. Ganoderma lucidum has traditionally been used for the treatment of diabetes in Chinese medicine; however, its anti-diabetic potency and mechanism in vivo is still unclear. Our previously published study reported a novel proteoglycan PTP1B inhibitor, named Fudan-Yueyang-Ganoderma lucidum (FYGL) from G. lucidum, with a half-maximal inhibitory concentration (IC₅₀) value of 5·12 (sem 0·05) μg/ml, a protein:polyglycan ratio of 17:77 and 78 % glucose in polysaccharide, and dominant amino acid residues of aspartic acid, glycine, glutamic acid, alanine, serine and threonine in protein. FYGL is capable of decreasing plasma glucose in streptozotocin-induced diabetic mice with a high safety of median lethal dose (LD₅₀) of 6 g/kg. In the present study, C57BL/6 db/db diabetic mice were trialed further using FYGL as well as metformin for comparison. Oral treatment with FYGL in db/db diabetic mice for 4 weeks significantly (P < 0·01 or 0·05) decreased the fasting plasma glucose level, serum insulin concentration and the homeostasis model assessment of insulin resistance. FYGL also controlled the biochemistry indices relative to type 2 diabetes-accompanied lipidaemic disorders. Pharmacology research suggests that FYGL decreases the plasma glucose level by the mechanism of inhibiting PTP1B expression and activity, consequently, regulating the tyrosine phosphorylation level of the IR β-subunit and the level of hepatic glycogen, thus resulting in the improvement of insulin sensitivity. Therefore, FYGL is promising as an insulin sensitiser for the therapy of type 2 diabetes and accompanied dyslipidaemia.

Molecular dynamics studies on the interactions of PTP1B with inhibitors: from the first phosphate-binding site to the second one

Protein tyrosine phosphatases 1B (PTP1B) is a major negative regulator of both insulin and leptin signaling pathways. In view of this, it becomes an important target for drug development against cancers, diabetes and obesity. The aim of the current study is to use the long time-scale molecular dynamics (MD) simulations to investigate the structural and dynamic factors that cause its inhibition by INTA and INTB, the two most potent and highly selective PTP1B inhibitors known so far. In order to investigate the mode of collective motions that is vitally important to the biological function, the covariance matrix of C(alpha) atoms was introduced for performing the dynamic analysis of the inhibition systems. It has been observed that the conformational and dynamic features of WPD-Loop, R-Loop and S-Loop play a key role in providing a smooth entrance for the inhibitors moving into the binding pocket as well as a favorable microenvironment to stabilize them. Furthermore, the hydrogen bonding networks formed around the active site with INTA and INTB may be the main reason of why the inhibition of PTP1B by the two ligands is so potent and selective. All these findings might provide useful insights for developing novel and effective drugs to treat cancer, diabetes and obesity.

Molecular modeling and dynamics studies of Shikimate Kinase from Bacillus anthracis

Bacillus anthracis has been used as weapon in bioterrorist activities, with high mortality, despite anti-microbial treatment, which strongly indicates a need of new drugs to treat anthrax. Shikimate Pathway is a seven-step biosynthetic route which generates chorismic acid. The shikimate pathway is essential for many pathological organisms, whereas it is absent in mammals. Therefore, these enzymes are potential targets for the development of non-toxic anti-microbial agents and herbicides and have been submitted to intensive structural studies. Shikimate Kinase is the fifth enzyme of shikimate pathway and catalyzes the specific phosphorylation of the 3-hydroxyl group of shikimate using ATP as a co-substrate, resulting in shikimate-3-phosphate and ADP. The present work describes for the first time a structural model for the Shikimate Kinase from B. anthracis using molecular modeling approach and molecular dynamics simulations. This study was able to identify the main residues of the ATP-binding and the shikimate pockets responsible for ligand affinities. Analysis of the molecular dynamics simulations indicates the structural features responsible for the stability of the structure. This study may help in the identification of new inhibitors for this enzyme.