Predicting class II MHC-Peptide binding:a kernel based approach using similarity scores

Background - Modelling the interaction between potentially antigenic peptides and Major Histocompatibility Complex (MHC) molecules is a key step in identifying potential T-cell epitopes. For Class II MHC alleles, the binding groove is open at both ends, causing ambiguity in the positional alignment between the groove and peptide, as well as creating uncertainty as to what parts of the peptide interact with the MHC. Moreover, the antigenic peptides have variable lengths, making naive modelling methods difficult to apply. This paper introduces a kernel method that can handle variable length peptides effectively by quantifying similarities between peptide sequences and integrating these into the kernel. Results - The kernel approach presented here shows increased prediction accuracy with a significantly higher number of true positives and negatives on multiple MHC class II alleles, when testing data sets from MHCPEP [1], MCHBN [2], and MHCBench [3]. Evaluation by cross validation, when segregating binders and non-binders, produced an average of 0.824 AROC for the MHCBench data sets (up from 0.756), and an average of 0.96 AROC for multiple alleles of the MHCPEP database. Conclusion - The method improves performance over existing state-of-the-art methods of MHC class II peptide binding predictions by using a custom, knowledge-based representation of peptides. Similarity scores, in contrast to a fixed-length, pocket-specific representation of amino acids, provide a flexible and powerful way of modelling MHC binding, and can easily be applied to other dynamic sequence problems.

Communication between the nucleotide binding domains of P-glycoprotein occurs via conformational changes that involve residue 508

Our aim is to provide molecular understanding of the mechanisms underlying the (i) interaction between the two nucleotide binding domains (NBDs) and (ii) coupling between NBDs and transmembrane domains within P-glycoprotein (Pgp) during a transport cycle. To facilitate this, we have introduced a number of unique cysteine residues at surface exposed positions (E393C, S452C, I500C, N508C, and K578C) in the N-terminal NBD of Pgp, which had previously been engineered to remove endogenous cysteines. Positions of the mutations were designed using a model based on crystallographic features of prokaryotic NBDs. The single cysteine mutants were expressed in insect cells using recombinant baculovirus and the proteins purified by metal affinity chromatography by virtue of a polyhistidine tag. None of the introduced cysteine residues perturbed the function of Pgp as judged by the characteristics of drug stimulated ATP hydrolysis. The role of residues at each of the introduced sites in the catalytic cycle of Pgp was investigated by the effect of covalent conjugation with N-ethyl-maleimide (NEM). All but one mutation (K578C) was accessible to labeling with [3H]-NEM. However, perturbation of ATPase activity was only observed for the derivitized N508C isoform. The principle functional manifestation was a marked inhibition of the "basal" rate of ATP hydrolysis. Neither the extent nor potency to which a range of drugs could affect the ATPase activity were altered in the NEM conjugated N508C isoform. The results imply that the accessibility of residue 508, located in the alpha-helical subdomain of NBD1 in Pgp, is altered by the conformational changes that occur during ATP hydrolysis.

Statistical deconvolution of enthalpic energetic contributions to MHC-peptide binding affinity

Background - MHC Class I molecules present antigenic peptides to cytotoxic T cells, which forms an integral part of the adaptive immune response. Peptides are bound within a groove formed by the MHC heavy chain. Previous approaches to MHC Class I-peptide binding prediction have largely concentrated on the peptide anchor residues located at the P2 and C-terminus positions. Results - A large dataset comprising MHC-peptide structural complexes was created by re-modelling pre-determined x-ray crystallographic structures. Static energetic analysis, following energy minimisation, was performed on the dataset in order to characterise interactions between bound peptides and the MHC Class I molecule, partitioning the interactions within the groove into van der Waals, electrostatic and total non-bonded energy contributions. Conclusion - The QSAR techniques of Genetic Function Approximation (GFA) and Genetic Partial Least Squares (G/PLS) algorithms were used to identify key interactions between the two molecules by comparing the calculated energy values with experimentally-determined BL50 data. Although the peptide termini binding interactions help ensure the stability of the MHC Class I-peptide complex, the central region of the peptide is also important in defining the specificity of the interaction. As thermodynamic studies indicate that peptide association and dissociation may be driven entropically, it may be necessary to incorporate entropic contributions into future calculations.

Modeling the peptide-T cell receptor interaction by the comparative molecular similarity indices analysis-soft independent modeling of class analogy technique

A set of 38 epitopes and 183 non-epitopes, which bind to alleles of the HLA-A3 supertype, was subjected to a combination of comparative molecular similarity indices analysis (CoMSIA) and soft independent modeling of class analogy (SIMCA). During the process of T cell recognition, T cell receptors (TCR) interact with the central section of the bound nonamer peptide; thus only positions 4−8 were considered in the study. The derived model distinguished 82% of the epitopes and 73% of the non-epitopes after cross-validation in five groups. The overall preference from the model is for polar amino acids with high electron density and the ability to form hydrogen bonds. These so-called “aggressive” amino acids are flanked by small-sized residues, which enable such residues to protrude from the binding cleft and take an active role in TCR-mediated T cell recognition. Combinations of “aggressive” and “passive” amino acids in the middle part of epitopes constitute a putative TCR binding motif

Coupling in silico and in vitro analysis of peptide-MHC binding:a bioinformatic approach enabling prediction of superbinding peptides and anchorless epitopes

The ability to define and manipulate the interaction of peptides with MHC molecules has immense immunological utility, with applications in epitope identification, vaccine design, and immunomodulation. However, the methods currently available for prediction of peptide-MHC binding are far from ideal. We recently described the application of a bioinformatic prediction method based on quantitative structure-affinity relationship methods to peptide-MHC binding. In this study we demonstrate the predictivity and utility of this approach. We determined the binding affinities of a set of 90 nonamer peptides for the MHC class I allele HLA-A*0201 using an in-house, FACS-based, MHC stabilization assay, and from these data we derived an additive quantitative structure-affinity relationship model for peptide interaction with the HLA-A*0201 molecule. Using this model we then designed a series of high affinity HLA-A2-binding peptides. Experimental analysis revealed that all these peptides showed high binding affinities to the HLA-A*0201 molecule, significantly higher than the highest previously recorded. In addition, by the use of systematic substitution at principal anchor positions 2 and 9, we showed that high binding peptides are tolerant to a wide range of nonpreferred amino acids. Our results support a model in which the affinity of peptide binding to MHC is determined by the interactions of amino acids at multiple positions with the MHC molecule and may be enhanced by enthalpic cooperativity between these component interactions.

Septin6 and Septin7 GTP binding proteins regulate AP-3- and ESCRT-dependent multivesicular body biogenesis

Septins (SEPTs) form a family of GTP-binding proteins implicated in cytoskeleton and membrane organization, cell division and host/pathogen interactions. The precise function of many family members remains elusive. We show that SEPT6 and SEPT7 complexes bound to F-actin regulate protein sorting during multivesicular body (MVB) biogenesis. These complexes bind AP-3, an adapter complex sorting cargos destined to remain in outer membranes of maturing endosomes, modulate AP-3 membrane interactions and the motility of AP-3-positive endosomes. These SEPT-AP interactions also influence the membrane interaction of ESCRT (endosomal-sorting complex required for transport)-I, which selects ubiquitinated cargos for degradation inside MVBs. Whereas our findings demonstrate that SEPT6 and SEPT7 function in the spatial, temporal organization of AP-3- and ESCRT-coated membrane domains, they uncover an unsuspected coordination of these sorting machineries during MVB biogenesis. This requires the E3 ubiquitin ligase LRSAM1, an AP-3 interactor regulating ESCRT-I sorting activity and whose mutations are linked with Charcot-Marie-Tooth neuropathies.

APP controls the formation of PI(3,5)P2 vesicles through its binding of the PIKfyve complex

Phosphoinositides are signalling lipids that are crucial for major signalling events as well as established regulators of membrane trafficking. Control of endosomal sorting and endosomal homeostasis requires phosphatidylinositol-3-phosphate (PI(3)P) and phosphatidylinositol-3,5-bisphosphate (PI(3,5)P2), the latter a lipid of low abundance but significant physiological relevance. PI(3,5)P2 is formed by phosphorylation of PI(3)P by the PIKfyve complex which is crucial for maintaining endosomal homeostasis. Interestingly, loss of PIKfyve function results in dramatic neurodegeneration. Despite the significance of PIKfyve, its regulation is still poorly understood. Here we show that the Amyloid Precursor Protein (APP), a central molecule in Alzheimer’s disease, associates with the PIKfyve complex (consisting of Vac14, PIKfyve and Fig4) and that the APP intracellular domain directly binds purified Vac14. We also show that the closely related APP paralogues, APLP1 and 2 associate with the PIKfyve complex. Whether APP family proteins can additionally form direct protein–protein interaction with PIKfyve or Fig4 remains to be explored. We show that APP binding to the PIKfyve complex drives formation of PI(3,5)P2 positive vesicles and that APP gene family members are required for supporting PIKfyve function. Interestingly, the PIKfyve complex is required for APP trafficking, suggesting a feedback loop in which APP, by binding to and stimulating PI(3,5)P2 vesicle formation may control its own trafficking. These data suggest that altered APP processing, as observed in Alzheimer’s disease, may disrupt PI(3,5)P2 metabolism, endosomal sorting and homeostasis with important implications for our understanding of the mechanism of neurodegeneration in Alzheimer’s disease.

Identifying key residues important for CGRP binding to its receptor

The calcitonin gene related peptide (CGRP) is a 37 amino acid neuropeptide. Its receptor is a heterodimeric complex of calcitonin receptor-like receptor (CLR) – a family B G-protein coupled receptor – and a single-pass transmembrane protein, receptoractivity modifying protein 1 (RAMP1). Here, we identify residues, within the N-terminal extracellular domain (ECD) of CLR, potentially involved in ligand binding.Certain residues presumed to be possible sites of contact for the CGRP were picked from the CLR/RAMP1 ECD crystal structure (PDB 3N7S). Residues were mutated to alanine (A) bysite-directed mutagenesis (QuikChangeTM, Stratagene). Mutants were analysed for their ability to stimulate cAMP and cell surface expression as previously described [1]. All mutants showed reduced potency, though to varying degrees as indicated by their pEC50 values. W69A and D70Ashowed significant reduction in cell surface expression.These findings suggest that these residues are important for the interaction of CGRP with its receptor. W69A and D70A, part of the WDG motif of family B GPCRs, are thought to rather play a role in receptor stability [2]. The data is consistent with CGRP binding in agroove between CLR and RAMP1. This project was supported byAston School of Life and Health Sciences.References1. Barwell J, Conner A & Poyner D (2011) Extracellular loops 1and 3 and their associated transmembrane regions of the calcitonin receptor-like receptor are needed for CGRP receptor function. Biochim Biophys Acta 1813, 1906–1916.2. Kumar S, Pioszak A, Zhang C et al. (2011) Crystal Structure of the PAC1R Extracellular Domain Unifies a Consensus Fold for Hormone Recognition by Class B G-Protein Cou-pled Receptors. PLoS One 6, e19682

Sequence-related structural DNA anomalies affect restriction enzyme activity and DNA polymerase I binding

DNA serves as a target molecule for several types of enzymes and may assume a wide variety of structural motifs depending upon the local sequence. The BssHII restriction site (GC)3 resides in a 9bp region of alternating pyrimidine and purine residues within the &phis;X174 genome. Such sequences are known to demonstrate non-canonical helical behavior under the appropriate conditions. The kinetics of BssHII cleavage was investigated in supercoiled and linear plasmid DNA, and in a 323bp DNA fragment obtained via amplification of &phis;X174. The rate of enzyme cleavage was enhanced in the supercoiled form and in the presence of 50μM cobalt hexamine. Similarly, cobalt hexamine was also found to enhance TaqI activity directly adjacent to the (GC)3 region. ^ Initial DNA polymerase I binding studies (including a gel mobility shift assay and a protection assay) indicated a notable interaction between DNA polymerase I and the BssHII site. An in-depth study revealed that equilibrium binding of DNA polymerase I to the T7 RNA polymerase promoter was comparable to that of the (GC)3 site, however the strongest interaction was observed with a cruciform containing region. Increasing the ionic strength of the solution environment, including the addition of DNA polymerase I reaction buffer significantly decreased the equilibrium dissociation constant values. ^ It is suggested that the region within or around the BssHII site experiences a conformational change generating a novel structure under the influence of supercoiled tension or 50μM cobalt hexamine. It is proposed that this transition may enhance enzyme activity and binding by providing an initial enzyme-docking site—the rate-limiting step in restriction enzyme kinetics. The high binding potential of DNA polymerase I for each of the motifs described, is hypothesized to be due to recognition of the structural DNA anomalies by the 3′–5′ exonuclease domain. ^

beta-Arrestin1 modulates lymphoid enhancer factor transcriptional activity through interaction with phosphorylated dishevelled proteins.

One aspect of the function of the beta-arrestins is to serve as scaffold or adapter molecules coupling G-protein coupled receptors (GPCRs) to signal transduction pathways distinct from traditional second messenger pathways. Here we report the identification of Dishevelled 1 and Dishevelled 2 (Dvl1 and Dvl2) as beta-arrestin1 (betaarr1) interacting proteins. Dvl proteins participate as key intermediates in signal transmission from the seven membrane-spanning Frizzled receptors leading to inhibition of glycogen synthase kinase-3beta (GSK-3beta), stabilization of beta-catenin, and activation of the lymphoid enhancer factor (LEF) transcription factor. We find that phosphorylation of Dvl strongly enhances its interaction with betaarr1, suggesting that regulation of Dvl phosphorylation and subsequent interaction with betaarr1 may play a key role in the activation of the LEF transcription pathway. Because coexpression of the Dvl kinases, CK1epsilon and PAR-1, with Dvl synergistically activates LEF reporter gene activity, we reasoned that coexpression of betaarr1 with Dvl might also affect LEF-dependent gene activation. Interestingly, whereas betaarr1 or Dvl alone leads to low-level stimulation of LEF (2- to 5-fold), coexpression of betaarr1 with either Dvl1 or Dvl2 leads to a synergistic activation of LEF (up to 16-fold). Additional experiments with LiCl as an inhibitor of GSK-3beta kinase activity indicate that the step affected by betaarr1 is upstream of GSK-3beta and most likely at the level of Dvl. These results identify betaarr1 as a regulator of Dvl-dependent LEF transcription and suggest that betaarr1 might serve as an adapter molecule that can couple Frizzled receptors and perhaps other GPCRs to these important transcription pathways.

Characterisation of an ubiquitin binding CUE domain in presenilin-1

The presenilins are the catalytic component of the gamma-secretase protease complex, involved in the regulated intramembrane proteolysis of numerous type-1 transmembrane proteins, including Amyloid precursor protein (APP) and Notch. In addition to their role in the γ-secretase complex the presenilins are involved in a number of γ-secretase independent functions such as calcium homeostasis, apoptosis, inflammation and protein trafficking. Presenilin function is known to be regulated through posttranslational modifications like endoproteolysis, phosphorylation and ubiquitination. Using a bioinformatics and protein sequence analysis approach this lab has identified a putative ubiquitin binding CUE domain in the presenilins. The aim of this project was to characterise the function of the presenilin CUE domains. Firstly, the presenilins are shown to contain a functional ubiquitin-binding CUE domain that preferentially binds to K63-linked polyubiquitin chains. The PS1 CUE domain is shown to be dispensable for PS1 endoproteolysis and γ-secretase mediated cleavage of APP, Notch and IL-1R1. This suggests the PS1 CUE domain is involved in a γ-secretase independent PS1 function. Our hypothesis is that the PS1 CUE domain is involved in regulating PS1’s intermolecular protein-protein interactions or intramolecular PS1:PS1 interactions. Here the PS1 CUE domain is shown to be dispensable for the interaction of PS1 and the K63-linked polyubiquitinated PS1 interacting proteins P75NTR, IL-1R1, TRAF6, TRAF2 and RIP1. To further investigate PS1 CUE domain function a mass spectrometry proteomics based approach is used to identify PS1 CUE domain interacting proteins. This proteomics approach demonstrated that the PS1 CUE domain is not required for PS1 dimerization. Instead a number of proteins thatinteract with the PS1 CUE domain are identified as well as proteins whose interaction with PS1 is downregulated by the presence of the PS1 CUE domain. Bioinformatic analysis of these proteins suggests possible roles for the PS1 CUE domain in regulating cell signalling, ubiquitination or cellular trafficking.

Interaction between Candida albicans and Pseudomonas aeruginosa

Fungal pathogen Candida albicans causes serious nosocomial infections in patients, in part, due to formation of drug-resistant biofilms. Protein kinases (PK) and transcription factors (TF) mediate signal transduction and transcription of proteins involved in biofilm development. To discover biofilm-related PKs, a collection of 63 C. albicans PK mutants was screened twice independently with microtiter plate-based biofilm assay (XTT). Thirty-eight (60%) mutants showed different degrees of biofilm impairment with the poor biofilm formers additionally possessing filamentation defects. Most of these genes were already known to encode proteins associated with Candida morphology and biofilms but VPS15, PKH3, PGA43, IME2 and CEX1, were firstly associated with both processes in this study. Previous studies of Holcombe et al. (2010) had shown that bacterial pathogen, Pseudomonas aeruginosa can impair C. albicans filamentation and biofilm development. To investigate their interaction, the good biofilm former PK mutants of C. albicans were assessed for their response to P. aeruginosa supernatants derived from two strains, wildtype PAO1 and homoserine lactone (HSL)-free mutant ΔQS, without finding any nonresponsive mutants. This suggested that none of the PKs in this study was implicated in Candida-Pseudomonas signaling. To screen promoter sequences for overrepresented TFs across C. albicans gene sets significantly up/downregulated in presence of bacterial supernatants from Holcombe et al. (2010) study, TFbsST database was created online. The TFbsST database integrates experimentally verified TFs of Candida to analyse promoter sequences for TF binding sites. In silico studies predicted that Efg1p was overrepresented in C. albicans and C. parapsilosis RBT family genes.

Exploring several nonconventional interactions of ice-binding proteins

Traditionally, ice-binding proteins (IBPs), also known as antifreeze proteins (AFPs), have been defined by two universal activities: ice recrystallization inhibition and thermal hysteresis. However, there remains the possibility IBPs have other complementary functions given the diversity found within this protein group. This thesis explores some of these in both natural and applied settings, in the hopes of furthering our understanding of this remarkable group of proteins. Plant IBPs could function as part of a defensive strategy against ice nucleators produced by certain pathogens. To assess this hypothesis, recombinant IBPs from perennial ryegrass and purple false brome were combined with the ice nucleation protein (INP) from the plant pathogen, Pseudomonas syringae. Strikingly, the plant proteins depressed the freezing point of the bacterial INP, while a fish AFP could not, nor did the INPs have any effect on IBP activity. Thus, the interaction between these two different proteins suggests a role in plant defensive strategies against pathogenic bacteria as another IBP function. In addition, the potential use of hyperactive insect IBPs in organ preservation was investigated. Current kidney preservation techniques involve storing the organ at 4 °C for a maximum of 24 h prior to transplantation. Extending this “safe” time would have profound effects on renal transplants, however, ischemic injury is prevalent when storage periods are prolonged. Experiments described here allowed subzero preservation for 72 h with the addition of a beetle IBP to CryoStasis® solution. Kidneys stored using the traditional technique for 24 h and the method developed here for 72 h showed similar levels of biomarker enzymes, underscoring the potential utility of insect IBPs for future transplant purposes. Finally, IBP function in the freeze-tolerant gall fly, Eurosta solidaginis, was examined. Larvae representing the mid-autumn stage displayed ice-binding activity, suggesting an IBP is being expressed, possibly as a protective measure against freezing damage when fall temperatures can unpredictably drop. IBP activity was also observed in the larvae’s host plant, Solidago spp. Mass spectrometry analysis of ice-affinity purified plant extracts provided three candidate pathogenesis-related proteins that could be responsible for the detected activity, further demonstrating additional functions of IBPs.

THE HERPESVIRUS NUCLEAR EGRESS COMPLEX COMPONENT, UL31, IS RECRUITED TO SITES OF DNA DAMAGE THROUGH POLY(ADP-RIBOSE) BINDING

The herpes simplex virus (HSV) UL31 gene encodes a conserved member of the herpesvirus nuclear egress complex that not only functions in the egress of DNA-containing capsids from the nucleus, but is also required for optimal viral genome expression, replication and packaging into capsids. Here, we report that the UL31 protein from HSV-2 and the orthologous protein, ORF69, from Kaposi's sarcoma-associated herpesvirus (KSHV) are recruited to sites of DNA damage. Recruitment of UL31 to sites of DNA damage occurred in HSV-2 infected cells, but did not require other viral proteins. The N-terminus of UL31 contains sequences resembling a poly(ADP-ribose) (PAR) binding motif. As protein poly-ADP ribosylation (PARylation) is a hallmark of the DNA damage response we examined the relationship between PARylation and UL31 recruitment to DNA damage. While the PAR polymerase (PARP)1/2 inhibitor, olaparib, prevented UL31 recruitment to damaged DNA, KU55933 inhibition of signaling through the ataxia telangiectasia mutated (ATM) DNA damage response pathway had no effect. These findings were further supported by experiments demonstrating direct and specific interaction between HSV-2 UL31 and PAR using purified components. Co-transfection with the viral kinase Us3, known to phosphorylate UL31, inhibited UL31 recruitment to DNA damage but also prevented the recruitment of other proteins recruited to DNA damage sites. The viral E3 ubiquitin ligase ICP0 was observed to co-localize with UL31 in transfected cells in a manner that is independent of the PAR-binding ability of UL31. However, inhibition of PARP1/2/3 did not reduce the ability of HSV-2 to replicate and we observed reduced PAR levels in the nuclei of infected cells. This study reveals a previously unrecognized function for UL31 orthologs and may suggest that the recognition of PAR by UL31 is coupled to the nuclear egress of herpesvirus capsids, influences viral DNA replication and packaging, or possibly modulates the DNA damage response mounted by virally infected cells.

von Hippel Lindau binding protein 1-mediated degradation of integrase affects HIV-1 gene expression at a postintegration step

HIV-1 integrase, the viral enzyme responsible for provirus integration into the host genome, can be actively degraded by the ubiquitin-proteasome pathway. Here, we identify von Hippel-Lindau binding protein 1(VBP1), a subunit of the prefoldin chaperone, as an integrase cellular binding protein that bridges interaction between integrase and the cullin2 (Cul2)-based von Hippel-Lindau (VHL) ubiquitin ligase. We demonstrate that VBP1 and Cul2/VHL are required for proper HIV-1 expression at a step between integrase-dependent proviral integration into the host genome and transcription of viral genes. Using both an siRNA approach and Cul2/VHL mutant cells, we show that VBP1 and the Cul2/VHL ligase cooperate in the efficient polyubiquitylation of integrase and its subsequent proteasome-mediated degradation. Results presented here support a role for integrase degradation by the prefoldin-VHL-proteasome pathway in the integration-transcription transition of the viral replication cycle.