MHC class II binding prediction:a little help from a friend

Vaccines are the greatest single instrument of prophylaxis against infectious diseases, with immeasurable benefits to human wellbeing. The accurate and reliable prediction of peptide-MHC binding is fundamental to the robust identification of T-cell epitopes and thus the successful design of peptide- and protein-based vaccines. The prediction of MHC class II peptide binding has hitherto proved recalcitrant and refractory. Here we illustrate the utility of existing computational tools for in silico prediction of peptides binding to class II MHCs. Most of the methods, tested in the present study, detect more than the half of the true binders in the top 5% of all possible nonamers generated from one protein. This number increases in the top 10% and 15% and then does not change significantly. For the top 15% the identified binders approach 86%. In terms of lab work this means 85% less expenditure on materials, labour and time. We show that while existing caveats are well founded, nonetheless use of computational models of class II binding can still offer viable help to the work of the immunologist and vaccinologist.

HLA-DP2 binding prediction by molecular dynamics simulations

Major histocompatibility complex (MHC) II proteins bind peptide fragments derived from pathogen antigens and present them at the cell surface for recognition by T cells. MHC proteins are divided into Class I and Class II. Human MHC Class II alleles are grouped into three loci: HLA-DP, HLA-DQ, and HLA-DR. They are involved in many autoimmune diseases. In contrast to HLA-DR and HLA-DQ proteins, the X-ray structure of the HLA-DP2 protein has been solved quite recently. In this study, we have used structure-based molecular dynamics simulation to derive a tool for rapid and accurate virtual screening for the prediction of HLA-DP2-peptide binding. A combinatorial library of 247 peptides was built using the "single amino acid substitution" approach and docked into the HLA-DP2 binding site. The complexes were simulated for 1 ns and the short range interaction energies (Lennard-Jones and Coulumb) were used as binding scores after normalization. The normalized values were collected into quantitative matrices (QMs) and their predictive abilities were validated on a large external test set. The validation shows that the best performing QM consisted of Lennard-Jones energies normalized over all positions for anchor residues only plus cross terms between anchor-residues.

Peptide length significantly influences in vitro affinity for MHC class II molecules

Class II Major Histocompatibility Complex (MHC) molecules have an open-ended binding groove which can accommodate peptides of varying lengths. Several studies have demonstrated that peptide flanking residues (PFRs) which lie outside the core binding groove can influence peptide binding and T cell recognition. By using data from the AntiJen database we were able to characterise systematically the influence of PFRs on peptide affinity for MHC class II molecules.

Receptor activity-modifying protein-dependent effects of mutations in the calcitonin receptor-like receptor:implications for adrenomedullin and calcitonin gene-related peptide pharmacology

Background and Purpose Receptor activity-modifying proteins (RAMPs) define the pharmacology of the calcitonin receptor-like receptor (CLR). The interactions of the different RAMPs with this class B GPCR yield high-affinity calcitonin gene-related peptide (CGRP) or adrenomedullin (AM) receptors. However, the mechanism for this is unclear. Experimental Approach Guided by receptor models, we mutated residues in the N-terminal helix of CLR, RAMP2 and RAMP3 hypothesized to be involved in peptide interactions. These were assayed for cAMP production with AM, AM2 and CGRP together with their cell surface expression. Binding studies were also conducted for selected mutants. Key Results An important domain for peptide interactions on CLR from I32 to I52 was defined. Although I41 was universally important for binding and receptor function, the role of other residues depended on both ligand and RAMP. Peptide binding to CLR/RAMP3 involved a more restricted range of residues than that to CLR/RAMP1 or CLR/RAMP2. E101 of RAMP2 had a major role in AM interactions, and F111/W84 of RAMP2/3 was important with each peptide. Conclusions and Implications RAMP-dependent effects of CLR mutations suggest that the different RAMPs control accessibility of peptides to binding residues situated on the CLR N-terminus. RAMP3 appears to alter the role of specific residues at the CLR-RAMP interface compared with RAMP1 and RAMP2. © 2013 The Authors. British Journal of Pharmacology published by John Wiley &. Sons Ltd on behalf of The British Pharmacological Society.

Large-scale molecular dynamics simulations of HLA-A*0201 complexed with a tumor-specific antigenic peptide:can the α3 and β2m domains be neglected?

Large-scale massively parallel molecular dynamics (MD) simulations of the human class I major histo-compatibility complex (MHC) protein HLA-A*0201 bound to a decameric tumor-specific antigenic peptide GVY-DGREHTV were performed using a scalable MD code on high-performance computing platforms. Such computational capabilities put us in reach of simulations of various scales and complexities. The supercomputing resources available Large-scale massively parallel molecular dynamics (MD) simulations of the human class I major histocompatibility complex (MHC) protein HLA-A*0201 bound to a decameric tumor-specific antigenic peptide GVYDGREHTV were performed using a scalable MD code on high-performance computing platforms. Such computational capabilities put us in reach of simulations of various scales and complexities. The supercomputing resources available for this study allow us to compare directly differences in the behavior of very large molecular models; in this case, the entire extracellular portion of the peptide–MHC complex vs. the isolated peptide binding domain. Comparison of the results from the partial and the whole system simulations indicates that the peptide is less tightly bound in the partial system than in the whole system. From a detailed study of conformations, solvent-accessible surface area, the nature of the water network structure, and the binding energies, we conclude that, when considering the conformation of the α1–α2 domain, the α3 and β2m domains cannot be neglected. © 2004 Wiley Periodicals, Inc. J Comput Chem 25: 1803–1813, 2004

Differential impact of amino acid substitutions on critical residues of the human glucagon-like peptide-1 receptor involved in peptide activity and small-molecule allosterys

The glucagon-like peptide-1 receptor (GLP-1R) is a class B G protein-coupled receptor that has a critical role in the regulation of glucose homeostasis, principally through the regulation of insulin secretion. The receptor systemis highly complex, able to be activated by both endogenous [GLP-1(1-36)NH2, GLP-1(1-37), GLP-1(7-36)NH2, GLP-1(7-37), oxyntomodulin], and exogenous (exendin-4) peptides in addition to small-molecule allosteric agonists (compound 2 [6,7-dichloro-2-methylsulfonyl-3-tertbutylaminoquinoxaline], BETP [4-(3-benzyloxy)phenyl)-2-ethylsulfinyl-6-(trifluoromethyl)pyrimidine]). Furthermore, the GLP-1R is subject to single-nucleotide polymorphic variance, resulting in amino acid changes in the receptor protein. In this study, we investigated two polymorphic variants previously reported to impact peptidemediated receptor activity (M149) and small-molecule allostery (C333). These residues were mutated to a series of alternate amino acids, and their functionality was monitored across physiologically significant signaling pathways, including cAMP, extracellular signal-regulated kinase 1 and 2 phosphorylation, and intracellular Ca2+ mobilization, in addition to peptide binding and cell-surface expression. We observed that residue 149 is highly sensitive to mutation, with almost all peptide responses significantly attenuated at mutated receptors. However, most reductions in activity were able to be restored by the small-molecule allosteric agonist compound 2. Conversely, mutation of residue 333 had little impact on peptide-mediated receptor activation, but this activity could not be modulated by compound 2 to the same extent as that observed at the wild-type receptor. These results provide insight into the importance of residues 149 and 333 in peptide function and highlight the complexities of allosteric modulation within this receptor system.

MHC class II binding prediction by molecular docking

Proteins of the Major Histocompatibility Complex (MHC) bind self and nonself peptide antigens or epitopes within the cell and present them at the cell surface for recognition by T cells. All T-cell epitopes are MHC binders but not all MCH binders are T-cell epitopes. The MHC class II proteins are extremely polymorphic. Polymorphic residues cluster in the peptide-binding region and largely determine the MHC's peptide selectivity. The peptide binding site on MHC class II proteins consist of five binding pockets. Using molecular docking, we have modelled the interactions between peptide and MHC class II proteins from locus DRB1. A combinatorial peptide library was generated by mutation of residues at peptide positions which correspond to binding pockets (so called anchor positions). The binding affinities were assessed using different scoring functions. The normalized scoring functions for each amino acid at each anchor position were used to construct quantitative matrices (QM) for MHC class II binding prediction. Models were validated by external test sets comprising 4540 known binders. Eighty percent of the known binders are identified in the best predicted 15% of all overlapping peptides, originating from one protein. © 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim.

Folding and Selective Exit of Reporter Proteins from the Yeast Endoplasmic Reticulum

The present study analyses the traffic of Hsp150 fusion proteins through the endoplasmic reticulum (ER) of yeast cells, from their post-translational translocation and folding to their exit from the ER via a selective COPI-independent pathway. The reporter proteins used in the present work are: Hsp150p, an O-glycosylated natural secretory protein of Saccharomyces cerevisiae, as well as fusion proteins consisting of a fragment of Hsp150 that facilitates in the yeast ER proper folding of heterologous proteins fused to it. It is thought that newly synthesized polypeptides are kept in an unfolded form by cytosolic chaperones to facilitate the post-translational translocation across the ER membrane. However, beta-lactamase, fused to the Hsp150 fragment, folds in the cytosol into bioactive conformation. Irreversible binding of benzylpenicillin locked beta-lactamase into a globular conformation, and prevented the translocation of the fusion protein. This indicates that under normal conditions the beta-lactamase portion unfolds for translocation. Cytosolic machinery must be responsible for the unfolding. The unfolding is a prerequisite for translocation through the Sec61 channel into the lumen of the ER, where the polypeptide is again folded into a bioactive and secretion-competent conformation. Lhs1p is a member of the Hsp70 family, which functions in the conformational repair of misfolded proteins in the yeast ER. It contains Hsp70 motifs, thus it has been thought to be an ATPase, like other Hsp70 members. In order to understand its activity, authentic Lhs1p and its recombinant forms expressed in E. coli, were purified. However, no ATPase activity of Lhs1p could be detected. Nor could physical interaction between Lhs1p and activators of the ER Hsp70 chaperone Kar2p, such as the J-domain proteins Sec63p, Scj1p, and Jem1p and the nucleotide exchange factor Sil1p, be demonstrated. The domain structure of Lhs1p was modelled, and found to consist of an ATPase-like domain, a domain resembling the peptide-binding domain (PBD) of Hsp70 proteins, and a C-terminal extension. Crosslinking experiments showed that Lhs1p and Kar2p interact. The interacting domains were the C-terminal extension of Lhs1p and the ATPase domain of Kar2p, and this interaction was independent of ATPase activity of Kar2p. A model is presented where the C-terminal part of Lhs1p forms a Bag-like 3 helices bundle that might serve in the nucleotide exchange function for Kar2p in translocation and folding of secretory proteins in the ER. Exit of secretory proteins in COPII-coated vesicles is believed to be dependent of retrograde transport from the Golgi to the ER in COPI-coated vesicles. It is thought that receptors escaping to the Golgi must be recycled back to the ER exit sites to recruit cargo proteins. We found that Hsp150 leaves the ER even in the absence of functional COPI-traffic from the Golgi to the ER. Thus, an alternative, COPI-independent ER exit pathway must exists, and Hsp150 is recruited to this route. The region containing the signature guiding Hsp150 to this alternative pathway was mapped.

Type III Secretion System of Phytopathogenic Bacterium Pseudomonas syringae:From Gene to Function

The type III secretion system (T3SS) is an essential requirement for the virulence of many Gram-negative bacteria which infect plants, animals and men. Pathogens use the T3SS to deliver effector proteins from the bacterial cytoplasm to the eukaryotic host cells, where the effectors subvert host defenses. The best candidates for directing effector protein traffic are the bacterial type III-associated appendages, called needles or pili. In plant pathogenic bacteria, the best characterized example of a T3SS-associated appendage is the HrpA pilus of the plant pathogen Pseudomonas syringae pv. tomato DC3000. The components of the T3SS in plant pathogens are encoded by a cluster of hrp (hypersensitive reaction and pathogenicity) genes. Two major classes of T3SS-secreted proteins are: harpin proteins such as HrpZ which are exported into extracellular space, and avirulence (Avr) proteins such as AvrPto which are translocated directly to the plant cytoplasm. This study deals with the structural and functional characterization of the T3SS-associated HrpA pilus and the T3SS-secreted harpins. By insertional mutagenesis analysis of HrpA, we located the optimal epitope insertion site in the amino-terminus of HrpA, and revealed the potential application of the HrpA pilus as a carrier of antigenic determinants for vaccination. By pulse-expression of proteins combined with immuno-electron microscopy, we discovered the Hrp pilus assembly strategy as addition of HrpA subunits to the distal end of the growing pilus, and we showed for the first time that secretion of HrpZ occurs at the tip of the pilus. The pilus thus functions as a conduit delivering proteins to the extracellular milieu. By using phage-display and scanning-insertion mutagenesis methods we identified a conserved HrpZ-binding peptide and localized the peptide-binding site to the central domain of HrpZ. We also found that the HrpZ specifically interacts with a host bean protein. Taken together, the current results provide deeper insight into the molecular mechanism of T3SS-associated pilus assembly and effector protein translocation, which will be helpful for further studies on the pathogenic mechanisms of Gram-negative bacteria and for developing new strategies to prevent bacterial infection.

Structural plasticity and enzyme action: Crystal structures of Mycobacterium tuberculosis peptidyl-tRNA hydrolase

Peptidyl-tRNA hydrolase cleaves the ester bond between tRNA and the attached peptide in peptidyl-tRNA in order to avoid the toxicity resulting from its accumulation and to free the tRNA available for further rounds in protein synthesis. The structure of the enzyme from Mycobacteritan tuberculosis has been determined in three crystal forms. This structure and the structure of the enzyme frorn Escherichia coli in its crystal differ substantially on account of the binding of the C terminus of the E. coli enzyme to the peptide-binding site of a neighboring molecule in the crystal. A detailed examination of this difference led to an elucidation of the plasticity of the binding site of the enzyme. The peptide-binding site of the enzyme is a cleft between the body, of the molecule and a polypepticle Y stretch involving a loop and a helix. This stretch is in the open conformation when the enzyme is in the free state as in the crystals of M. tuberculosis peptidyl-tRNA hydrolase. Furthermore, there is no physical continuity between the tRNA and the peptide-binding sites. The molecule in the E. coli crystal mimics the peptide-bound enzyme molecule. The peptide stretch referred to earlier now closes on the bound peptide. Concurrently, a channel connecting the tRNA and the peptide-binding site opens primarily through the concerted movement of two residues. Thus, the crystal structure of M. tuberculosis peptidyl-tRNA hydrolase when compared with the crystal structure of the E. coli enzyme, leads to a model of structural changes associated with enzyme action on the basis of the plasticity of the molecule. (c) 2007 Elsevier Ltd. All rights reserved.

Structural plasticity and enzyme action: Crystal structures of Mycobacterium tuberculosis peptidyl-tRNA hydrolase

Peptidyl-tRNA hydrolase cleaves the ester bond between tRNA and the attached peptide in peptidyl-tRNA in order to avoid the toxicity resulting from its accumulation and to free the tRNA available for further rounds in protein synthesis. The structure of the enzyme from Mycobacteritan tuberculosis has been determined in three crystal forms. This structure and the structure of the enzyme frorn Escherichia coli in its crystal differ substantially on account of the binding of the C terminus of the E. coli enzyme to the peptide-binding site of a neighboring molecule in the crystal. A detailed examination of this difference led to an elucidation of the plasticity of the binding site of the enzyme. The peptide-binding site of the enzyme is a cleft between the body, of the molecule and a polypepticle Y stretch involving a loop and a helix. This stretch is in the open conformation when the enzyme is in the free state as in the crystals of M. tuberculosis peptidyl-tRNA hydrolase. Furthermore, there is no physical continuity between the tRNA and the peptide-binding sites. The molecule in the E. coli crystal mimics the peptide-bound enzyme molecule. The peptide stretch referred to earlier now closes on the bound peptide. Concurrently, a channel connecting the tRNA and the peptide-binding site opens primarily through the concerted movement of two residues. Thus, the crystal structure of M. tuberculosis peptidyl-tRNA hydrolase when compared with the crystal structure of the E. coli enzyme, leads to a model of structural changes associated with enzyme action on the basis of the plasticity of the molecule. (c) 2007 Elsevier Ltd. All rights reserved.

Functional role of the Mso1p-Sec1p complex in membrane fusion regulation

Sec1/Munc18 (SM) protein family members are evolutionary conserved proteins. They perform an essential, albeit poorly understood function in SNARE complex formation in membrane fusion. In addition to the SNARE complex components, only a few SM protein binding proteins are known. Typically, their binding modes to SM proteins and their contribution to the membrane fusion regulation is poorly characterised. We identified Mso1p as a novel Sec1p interacting partner. It was shown that Mso1p and Sec1p interact at sites of polarised secretion and that this localisation is dependent on the Rab GTPase Sec4p and its GEF Sec2p. Using targeted mutagenesis and N- and C-terminal deletants, it was discovered that the interaction between an N-terminal peptide of Mso1p and the putative Syntaxin N-peptide binding area in Sec1p domain 1 is important for membrane fusion regulation. The yeast Syntaxin homologues Sso1p and Sso2p lack the N-terminal peptide. Our results show that in addition to binding to the putative N-peptide binding area in Sec1p, Mso1p can interact with Sso1p and Sso2p. This result suggests that Mso1p can mimic the N-peptide binding to facilitate membrane fusion. In addition to Mso1p, a novel role in membrane fusion regulation was revealed for the Sec1p C-terminal tail, which is missing in its mammalian homologues. Deletion of the Sec1p-tail results in temperature sensitive growth and reduced sporulation. Using in vivo and in vitro experiments, it was shown that the Sec1p-tail mediates SNARE complex binding and assembly. These results propose a regulatory role for the Sec1p-tail in SNARE complex formation. Furthermore, two novel interaction partners for Mso1p, the Rab GTPase Sec4p and plasma membrane phospholipids, were identified. The Sec4p link was identified using Bimolecular Fluorescence Complementation assays with Mso1p and the non-SNARE binding Sec1p(1-657). The assay revealed that Mso1p can target Sec1p(1-657) to sites of secretion. This effect is mediated via the Mso1p C-terminus, which previously has been genetically linked to Sec4p. These results and in vitro binding experiments suggest that Mso1p acts in cooperation with the GTP-bound form of Sec4p on vesicle-like structures prior to membrane fusion. Mso1p shares homology with the PIP2 binding domain of the mammalian Munc18 binding Mint proteins. It was shown both in vivo and in vitro that Mso1p is a phospholipid inserting protein and that this insertion is mediated by the conserved Mso1p amino terminus. In vivo, the Mso1p phospholipid binding is needed for sporulation and Mso1p-Sec1p localisation at the sites of secretion at the plasma membrane. The results reveal a novel layer of membrane fusion regulation in exocytosis and propose a coordinating role for Mso1p in connection with membrane lipids, Sec1p, Sec4p and SNARE complexes in this process.

Structures of new crystal forms of Mycobacterium tuberculosis peptidyl-tRNA hydrolase and functionally important plasticity of the molecule

The X-ray structures of new crystal forms of peptidyl-tRNA hydrolase from M.similar to tuberculosis reported here and the results of previous X-ray studies of the enzyme from different sources provide a picture of the functionally relevant plasticity of the protein molecule. The new X-ray results confirm the connection deduced previously between the closure of the lid at the peptide-binding site and the opening of the gate that separates the peptide-binding and tRNA-binding sites. The plasticity of the molecule indicated by X-ray structures is in general agreement with that deduced from the available solution NMR results. The correlation between the lid and the gate movements is not, however, observed in the NMR structure.

Cloning and expression analysis of a HSP70 gene from Pacific abalone (Haliotis discus hannai)

Heat shock protein 70 (HSP70), the primary member of HSPs that are responsive of thermal stress, is found in all multicellular organisms and functions mostly as molecular chaperon. The inducible HSP70 cDNA cloned from Pacific abalone (Haliotis discus hannai) using rapid amplification of cDNA ends (RACE), was highly homologous to other HSP70 genes. The full-length cDNA of the Pacific abalone HSP70 was 2631 bp, consisting of a 5'-terminal untranslated region (UTR) of 90 bp, a 3'-terminal UTR of 573 by with a canonical polyadenylation signal sequence AATAAA and a poly (A) tail, and an open reading frame of 1968 bp. The HSP70 cDNA encoded a polypeptide of 655 amino acids with an ATPase domain of 382 amino acids, the substrate peptide binding domain of 161 amino acids and a C-terminus domain of 112 amino acids. The temporal expression of HSP70 was measured by semi-quantitative RT-PCR after heat shock and bacterial challenge. Challenge of Pacific abalone with heat shock or the pathogenic bacteria Vibrio anguillarum resulted in a dramatic increase in the expression of HSP70 mRNA level in muscle, followed by a recovery to normal level after 96 h. Unlike the muscle, the levels of HSP70 expression in gills reached the top at 12 h and maintained a relatively high level compared with the control after thermal and bacterial challenge. The upregulated mRNA expression of HSP70 in the abalone following heat shock and infection response indicates that the HSP70 gene is inducible and involved in immune response. (c) 2006 Elsevier Ltd. All rights reserved.

Comparative immunogenicity of HIV-1 clade C envelope proteins for prime/boost studies.

BACKGROUND: Previous clinical efficacy trials failed to support the continued development of recombinant gp120 (rgp120) as a candidate HIV vaccine. However, the recent RV144 HIV vaccine trial in Thailand showed that a prime/boost immunization strategy involving priming with canarypox vCP1521 followed by boosting with rgp120 could provide significant, although modest, protection from HIV infection. Based on these results, there is renewed interest in the development of rgp120 based antigens for follow up vaccine trials, where this immunization approach can be applied to other cohorts at high risk for HIV infection. Of particular interest are cohorts in Africa, India, and China that are infected with clade C viruses. METHODOLOGY/PRINCIPAL FINDINGS: A panel of 10 clade C rgp120 envelope proteins was expressed in 293 cells, purified by immunoaffinity chromatography, and used to immunize guinea pigs. The resulting sera were collected and analyzed in checkerboard experiments for rgp120 binding, V3 peptide binding, and CD4 blocking activity. Virus neutralization studies were carried out with two different assays and two different panels of clade C viruses. A high degree of cross reactivity against clade C and clade B viruses and viral proteins was observed. Most, but not all of the immunogens tested elicited antibodies that neutralized tier 1 clade B viruses, and some sera neutralized multiple clade C viruses. Immunization with rgp120 from the CN97001 strain of HIV appeared to elicit higher cross neutralizing antibody titers than the other antigens tested. CONCLUSIONS/SIGNIFICANCE: While all of the clade C antigens tested were immunogenic, some were more effective than others in eliciting virus neutralizing antibodies. Neutralization titers did not correlate with rgp120 binding, V3 peptide binding, or CD4 blocking activity. CN97001 rgp120 elicited the highest level of neutralizing antibodies, and should be considered for further HIV vaccine development studies.