Highly flexible protein-peptide docking using CABS-dock

Protein-peptide molecular docking is a difficult modeling problem. It is even more challenging when significant conformational changes that may occur during the binding process need to be predicted. In this chapter, we demonstrate the capabilities and features of the CABS-dock server for flexible protein-peptide docking. CABS-dock allows highly efficient modeling of full peptide flexibility and significant flexibility of a protein receptor. During CABS-dock docking, the peptide folding and binding process is explicitly simulated and no information about the peptide binding site or its structure, is used. This chapter presents a successful CABS-dock use for docking a potentially therapeutic peptide to a protein target. Moreover, simulation contact maps, a new CABS-dock feature, are described and applied to the docking test case. Finally, a tutorial for running CABS-dock from the command line or command line scripts is provided. The CABS-dock web server is available from http://biocomp.chem.uw.edu.pl/CABSdock

Structure-based selection of human metabolite binding P4 pocket of DRB1*15:01 and DRB1*15:03, with implications for multiple sclerosis.

Binding of small molecules in the human leukocyte antigen (HLA) peptide-binding groove may result in conformational changes of bound peptide and an altered immune response, but previous studies have not considered a potential role for endogenous metabolites. We performed virtual screening of the complete Human Metabolite Database (HMDB) for docking to the multiple sclerosis (MS) susceptible DRB1*15:01 allele and compared the results to the closely related yet non-susceptible DRB1*15:03 allele; and assessed the potential impact on binding of human myelin basic peptide (MBP). We observed higher energy scores for metabolite binding to DRB1*15:01 than DRB1*15:03. Structural comparison of docked metabolites with DRB1*15:01 and DRB1*15:03 complexed with MBP revealed that PhenylalanineMBP92 allows binding of metabolites in the P4 pocket of DRB1*15:01 but ValineMBP89 abrogates metabolite binding in the P1 pocket. We observed differences in the energy scores for binding of metabolites in the P4 pockets of DRB1*15:01 vs. DRB1*15:03 suggesting stronger binding to DRB1*15:01. Our study confirmed that specific, disease-associated human metabolites bind effectively with the most polymorphic P4 pocket of DRB1*15:01, the primary MS susceptible allele in most populations. Our results suggest that endogenous human metabolites bound in specific pockets of HLA may be immunomodulatory and implicated in autoimmune disease

Whey-derived peptides interactions with ACE by molecular docking as a potential predictive tool of natural ACE inhibitors

Several milk/whey derived peptides possess high in vitro angiotensin I-converting enzyme (ACE) inhibitory activity. However, in some cases, poor correlation between the in vitro ACE inhibitory activity and the in vivo antihypertensive activity has been observed. The aim of this study is to gain insight into the structure-activity relationship of peptide sequences present in whey/milk protein hydrolysates with high ACE inhibitory activity, which could lead to a better understanding and prediction of their in vivo antihypertensive activity. The potential interactions between peptides produced from whey proteins, previously reported as high ACE inhibitors such as IPP, LIVTQ, IIAE, LVYPFP, and human ACE were assessed using a molecular docking approach. The results show that peptides IIAE, LIVTQ, and LVYPFP formed strong H bonds with the amino acids Gln 259, His 331, and Thr 358 in the active site of the human ACE. Interestingly, the same residues were found to form strong hydrogen bonds with the ACE inhibitory drug Sampatrilat. Furthermore, peptides IIAE and LVYPFP interacted with the amino acid residues Gln 259 and His 331, respectively, also in common with other ACE-inhibitory drugs such as Captopril, Lisinopril and Elanapril. Additionally, IIAE interacted with the amino acid residue Asp 140 in common with Lisinopril, and LIVTQ interacted with Ala 332 in common with both Lisinopril and Elanapril. The peptides produced naturally from whey by enzymatic hydrolysis interacted with residues of the human ACE in common with potent ACE-inhibitory drugs which suggests that these natural peptides may be potent ACE inhibitors

Molecular dynamics simulations and in silico peptide ligand screening of the Elk-1 ETS domain

Background: The Elk-1 transcription factor is a member of a group of proteins called ternary complex factors, which serve as a paradigm for gene regulation in response to extracellular signals. Its deregulation has been linked to multiple human diseases including the development of tumours. The work herein aims to inform the design of potential peptidomimetic compounds that can inhibit the formation of the Elk-1 dimer, which is key to Elk-1 stability. We have conducted molecular dynamics simulations of the Elk-1 ETS domain followed by virtual screening. Results: We show the ETS dimerisation site undergoes conformational reorganisation at the a1b1 loop. Through exhaustive screening of di- and tri-peptide libraries against a collection of ETS domain conformations representing the dynamics of the loop, we identified a series of potential binders for the Elk-1 dimer interface. The di-peptides showed no particular preference toward the binding site; however, the tri-peptides made specific interactions with residues: Glu17, Gln18 and Arg49 that are pivotal to the dimer interface. Conclusions: We have shown molecular dynamics simulations can be combined with virtual peptide screening to obtain an exhaustive docking protocol that incorporates dynamic fluctuations in a receptor. Based on our findings, we suggest experimental binding studies to be performed on the 12 SILE ranked tri-peptides as possible compounds for the design of inhibitors of Elk-1 dimerisation. It would also be reasonable to consider the score ranked tri-peptides as a comparative test to establish whether peptide size is a determinant factor of binding to the ETS domain

Structural investigation of nucleophosmin interaction with the tumor suppressor Fbw7γ

Nucleophosmin (NPM1) is a multifunctional nucleolar protein implicated in ribogenesis, centrosome duplication, cell cycle control, regulation of DNA repair and apoptotic response to stress stimuli. The majority of these functions are played through the interactions with a variety of protein partners. NPM1 is frequently overexpressed in solid tumors of different histological origin. Furthermore NPM1 is the most frequently mutated protein in acute myeloid leukemia (AML) patients. Mutations map to the C-terminal domain and lead to the aberrant and stable localization of the protein in the cytoplasm of leukemic blasts. Among NPM1 protein partners, a pivotal role is played by the tumor suppressor Fbw7γ, an E3-ubiquitin ligase that degrades oncoproteins like c-MYC, cyclin E, Notch and c-jun. In AML with NPM1 mutations, Fbw7γ is degraded following its abnormal cytosolic delocalization by mutated NPM1. This mechanism also applies to other tumor suppressors and it has been suggested that it may play a key role in leukemogenesis. Here we analyse the interaction between NPM1 and Fbw7γ, by identifying the protein surfaces implicated in recognition and key aminoacids involved. Based on the results of computational methods, we propose a structural model for the interaction, which is substantiated by experimental findings on several site-directed mutants. We also extend the analysis to two other NPM1 partners (HIV Tat and CENP-W) and conclude that NPM1 uses the same molecular surface as a platform for recognizing different protein partners. We suggest that this region of NPM1 may be targeted for cancer treatment

Import of honeybee prepromelittin into the endoplasmic reticulum

Honeybee prepromelittin is correctly processed and imported by dog pancreas microsomes. Insertion of prepromelittin into microsomal membranes, as assayed by signal sequence removal, does not depend on signal recognition particle (SRP) and docking protein. We addressed the question as to how prepromelittin bypasses the SRP/docking protein system. Hybrid proteins between prepromelittin, or carboxy-terminally truncated derivatives, and the cytoplasmic protein dihydrofolate reductase from mouse were constructed. These hybrid proteins were analysed for membrane insertion and sequestration into microsomes. The results suggest the following: (i) The signal sequence of prepromelittin is capable of interacting with the SRP/docking protein system, but this interaction is not mandatory for membrane insertion; this is related to the small size of prepromelittin. (ii) In prepromelittin a cluster of negatively charged amino acids must be balanced by a cluster of positively charged amino acids in order to allow membrane insertion. (iii) In general, a signal sequence can be sufficient to mediate membrane insertion independently of SRP and docking protein in the case of short precursor proteins; however, the presence and distribution of charged amino acids within the mature part of these precursors can play distinct roles

Substrate specificity provides insights into the sugar donor recognition mechanism of O-GlcNAc transferase (OGT).

O-Linked β-N-acetylglucosaminyl transferase (OGT) plays an important role in the glycosylation of proteins, which is involved in various cellular events. In human, three isoforms of OGT (short OGT [sOGT]; mitochondrial OGT [mOGT]; and nucleocytoplasmic OGT [ncOGT]) share the same catalytic domain, implying that they might adopt a similar catalytic mechanism, including sugar donor recognition. In this work, the sugar-nucleotide tolerance of sOGT was investigated. Among a series of uridine 5'-diphosphate-N-acetylglucosamine (UDP-GlcNAc) analogs tested using the casein kinase II (CKII) peptide as the sugar acceptor, four compounds could be used by sOGT, including UDP-6-deoxy-GlcNAc, UDP-GlcNPr, UDP-6-deoxy-GalNAc and UDP-4-deoxy-GlcNAc. Determined values of Km showed that the substitution of the N-acyl group, deoxy modification of C6/C4-OH or epimerization of C4-OH of the GlcNAc in UDP-GlcNAc decreased its affinity to sOGT. A molecular docking study combined with site-directed mutagenesis indicated that the backbone carbonyl oxygen of Leu653 and the hydroxyl group of Thr560 in sOGT contributed to the recognition of the sugar moiety via hydrogen bonds. The close vicinity between Met501 and the N-acyl group of GlcNPr, as well as the hydrophobic environment near Met501, were responsible for the selective binding of UDP-GlcNPr. These findings illustrate the interaction of OGT and sugar nucleotide donor, providing insights into the OGT catalytic mechanism

Structural Mechanism for the Specific Assembly and Activation of the Extracellular Signal Regulated Kinase 5 (ERK5) Module

Mitogen-activated protein kinase (MAPK) activation depends on a linear binding motif found in all MAPK kinases (MKK). In addition, the PB1 (Phox and Bem1) domain of MKK5 is required for extracellular signal regulated kinase 5 (ERK5) activation. We present the crystal structure of ERK5 in complex with an MKK5 construct comprised of the PB1 domain and the linear binding motif. We show that ERK5 has distinct protein-protein interaction surfaces compared with ERK2, which is the closest ERK5 paralog. The two MAPKs have characteristically different physiological functions and their distinct protein-protein interaction surface topography enables them to bind different sets of activators and substrates. Structural and biochemical characterization revealed that the MKK5 PB1 domain cooperates with the MAPK binding linear motif to achieve substrate specific binding, and it also enables co-recruitment of the upstream activating enzyme and the downstream substrate into one signaling competent complex. Studies on present day MAPKs and MKKs hint on the way protein kinase networks may evolve. In particular, they suggest how paralogous enzymes with similar catalytic properties could acquire novel signaling roles by merely changing the way they make physical links to other proteins