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.

Mechanism of attenuation of angiotensin-II-induced protein degradation by insulin-like growth factor-I (IGF-I)

Insulin-like growth factor-I (IGF-I) has been shown to attenuate protein degradation in murine myotubes induced by angiotensin II through downregulation of the ubiquitin-proteasome pathway, although the mechanism is not known. Angiotensin II is known to upregulate this pathway through a cellular signalling mechanism involving release of arachidonic acid, activation of protein kinase Cα (PKCα), degradation of inhibitor-κB (I-κB) and nuclear migration of nuclear factor-κB (NF-κB), and all of these events were attenuated by IGF-I (13.2 nM). Induction of the ubiquitin-proteasome pathway has been linked to activation of the RNA-activated protein kinase (PKR), since an inhibitor of PKR attenuated proteasome expression and activity in response to angiotensin II and prevented the decrease in the myofibrillar protein myosin. Angiotensin II induced phosphorylation of PKR and of the eukaryotic initiation factor-2 (eIF2) on the α-subunit, and this was attenuated by IGF-I, by induction of the expression of protein phosphatase 1, which dephosphorylates PKR. Release of arachidonic acid and activation of PKCα by angiotensin II were attenuated by an inhibitor of PKR and IGF-I, and the effect was reversed by Salubrinal (15 μM), an inhibitor of eIF2α dephosphorylation, as was activation of PKCα. In addition myotubes transfected with a dominant-negative PKR (PKRΔ6) showed no release of arachidonate in response to Ang II, and no activation of PKCα. These results suggest that phosphorylation of PKR by angiotensin II was responsible for the activation of the PLA2/PKC pathway leading to activation of NF-κB and that IGF-I attenuates protein degradation due to an inhibitory effect on activation of PKR. © 2007 Elsevier Inc. All rights reserved.

Increased expression of the ubiquitin-proteasome pathway in murine myotubes by proteolysis-inducing factor (PIF) is associated with activation of the transcription factor NF-kappa B

Proteolysis-inducing factor (PIF), isolated from a cachexia-inducing murine tumour, has been shown to stimulate protein breakdown in C 2C12 myotubes. The effect was attenuated by the specific proteasome inhibitor lactacystin and there was an elevation of proteasome 'chymotrypsin-like' enzyme activity and expression of 205 proteasome α-subunits at concentrations of PIF between 2 and 16 nM. Higher concentrations of PIF had no effect. The action of PIF was attenuated by eicosapentaenoic acid (EPA) (50 μM). At a concentration of 4 nM, PIF induced a transient decrease in IκBα levels after 30 min incubation, while no effect was seen at 20 nM PIF. The level of IκBα, an NF-κB inhibitory protein, returned to normal after 60 min. Depletion of IκBα from the cytosol was not seen in myotubes pretreated with EPA, suggesting that the NF-κB/IκB complex was stabilised. At concentrations between 2 and 8 nM, PIF stimulated an increased nuclear migration of NF-κB, which was not seen in myotubes pretreated with EPA. The PIF-induced increase in chymotrypsin-like enzyme activity was also attenuated by the NF-κB inhibitor peptide SN50, suggesting that NF-κB may be involved in the PIF-induced increase in proteasome expression. The results further suggest that EPA may attenuate protein degradation induced by PIF, at least partly, by preventing NF-κB accumulation in the nucleus. © 2003 Cancer Research UK.

Induction of protein catabolism in myotubes by 15(S)-hydroxyeicosatetraenoic acid through increased expression of the ubiquitin-proteasome pathway

The potential role of 15(S)-hydroxyeicosatetraenoic acid (15(S)-HETE) as an intracellular signal for increased protein catabolism and induction of the expression of key components of the ubiquitin-proteasome proteolytic pathway induced by a tumour cachectic factor, proteolysis-inducing factor has been studied in murine C2C12 myotubes. 15(S)-HETE induced protein degradation in these cells with a maximal effect at concentrations between 78 and 312 nM. The effect was attenuated by the polyunsaturated fatty acid, eicosapentaenoic acid (EPA). There was an increase in 'chymotrypsin-like' enzyme activity, the predominant proteolytic activity of the proteasome, in the same concentration range as that inducing total protein degradation, and this effect was also attenuated by EPA. 15(S)-hydroxyeicosatetraenoic acid also increased maximal expression of mRNA for proteasome subunits C2 and C5, as well as the ubiquitin-conjugating enzyme, E214k, after 4 h incubation, as determined by quantitative competitive RT-PCR. The concentrations of 15-HETE affecting gene expression were the same as those inducing protein degradation. Western blotting of cellular supernatants of myotubes treated with 15(S)-HETE for 24 h showed increased expression of p42, an ATPase subunit of the regulatory complex at similar concentrations, as well as a decrease in expression of myosin in the same concentration range. 15(S)-hydroxyeicosatetraenoic acid activated binding of nuclear factor-κB (NF-κB) in the myotube nucleus and stimulated degradation of 1-κBα. The effect on the NF-κB/1-κBα system was attenuated by EPA. In addition, the NF-κB inhibitor peptide SN50 attenuated the increased chymotrypsin-like enzyme activity in the presence of 15(S)-HETE. These results suggest that 15(S)-HETE induces degradation of myofibrillar proteins in differentiated myotubes through an induction of an increased expression of the regulatory components of the ubiquitin-proteasome proteolytic pathway possibly through the intervention of the nuclear transcription factor NF-κB, and that this process is inhibited by EPA. © 2003 Cancer Research UK.

Modulation of certain 5-HT-related behaviours by serotonergic and noradrenergic systems

The modulation of 5-hydroxytryptamine (5-HT)-related head-twitchbehaviour by antimigraine drugs and migraine triggers was examined inmice. The antimigraine drugs examined produced either inhibition or noeffect on 5-HT-related head-twitching. On the basis of these resultsit is suggested that 5-HT-related head-twitching is unlikely to beuseful in the preclinical screening and discovery of systemically-activeantimigraine agents. The migraine triggers examined, tyramineand beta-PEA initially produced a repeatable complex time-relatedeffect on 5-HT-related head-twitching, with both inhibition andpotentiation of this behaviour being observed, however, when furtherexamination of the effect of the migraine triggers on 5-HT-relatedhead-twitching was attempted some time later the effects seeninitially were no longer produced. The effect of (±)-1-<2, 5-dimethoxy-4-iodophenyl)-2-aminopropane,((±)DOl), on on-going behaviour of mice and rats was examined. Shakingbehaviour was observed in both species. In mice, excessive scratchingbehaviour was also present. (±)DOl-induced scratching and shakingbehaviour were found to be differentially modulated by noradrenergicand serotonergic agents, however, the fact that both behaviours wereblocked by ritanserin (5-HT2/5-HT1c receptor antagonist) and inhibitedby FLA-63 (a dopamine-beta-oxidase inhibitor which depletesnoradrenaline), suggests the pathways mediating these behaviours mustbe convergent in some manner, and that both behaviours require intact5-HT receptors, probably 5-HT2 receptors, for their production. Ingeneral, the behavioural profile of (±)DOI was as expected for anagent which exhibits high affinity binding to 5-HT2/5-HT1c receptors.Little sign of the 5-HTl-related '5-HT syndrome' was seen in eithermice or rats. The effect of a variety of noradrenergic agents on head-twitchinginduced by a variety of shake-inducing agents was examined. A patternof modulatory effect was seen whereby the modulatory effect of thenoradrenergic agents on 5-hydroxytryptophan <5-HTP) (and in some cases, 5-methoxy-N,N-dimethyltryptamine (5-MeODMT)) was found to be the opposite of that observed with quipazine and (±)DOI. The relationship between these effects, and their implications for understanding the pharmacology of centrally acting drugs is discussed.

Mechanism of attenuation of muscle protein degradation induced by tumor necrosis factor-alpha and angiotensin II by beta-hydroxy-beta-methylbutyrate

Both tumor necrosis factor-alpha (TNF-alpha)/interferon-gamma (IFN-gamma) and angiotensin II (ANG II) induced an increase in total protein degradation in murine myotubes, which was completely attenuated by treatment with beta-hydroxy-beta-methylbutyrate (HMB; 50 microM). There was an increase in formation of reactive oxygen species (ROS) within 30 min, as well as an increase in the activity of both caspase-3 and -8, and both effects were attenuated by HMB. Moreover, inhibitors of caspase-3 and -8 completely attenuated both ROS formation and total protein degradation induced by TNF-alpha/IFN-gamma and ANG II. There was an increased autophosphorylation of double-stranded RNA-dependent protein kinase (PKR), which was attenuated by the specific caspase-3 and -8 inhibitors. Neither ROS formation or protein degradation occurred in myotubes expressing a catalytically inactive PKR variant, PKRDelta6, in response to TNF-alpha/IFN-gamma, compared with myotubes expressing wild-type PKR, although there was still activation of caspase-3 and -8. HMB also attenuated activation of PKR, suggesting that it was important in protein degradation. Formation of ROS was attenuated by rotenone, an inhibitor of the mitochondrial electron transport chain, nitro-l-arginine methyl ester, an inhibitor of nitric oxide synthase, and SB 203580, a specific inhibitor of p38 mitogen-activated protein kinase (p38 MAPK), which also attenuated total protein degradation. Activation of p38 MAPK by PKR provides the link to ROS formation. These results suggest that TNF-alpha/IFN-gamma and ANG II induce muscle protein degradation by a common signaling pathway, which is attenuated by HMB, and that this involves the initial activation of caspase-3 and -8, followed by autophosphorylation and activation of PKR, which then leads to increased ROS formation via activation of p38 MAPK. Increased ROS formation is known to induce protein degradation through the ubiquitin-proteasome pathway.

Attenuation of depression of muscle protein synthesis induced by lipopolysaccharide, tumor necrosis factor, and angiotensin II by beta-hydroxy-beta-methylbutyrate

beta-Hydroxy-beta-methylbutyrate (HMB; 50 microM) has been shown to attenuate the depression in protein synthesis in murine myotubes in response to lipopolysaccharide (LPS), tumor necrosis factor-alpha (TNF-alpha) with or without interferon-gamma (IFN-gamma), and angiotensin II (ANG II). The mechanism for the depression of protein synthesis by all three agents was the same and was attributed to activation of double-stranded RNA-dependent protein kinase (PKR) with the subsequent phosphorylation of eukaryotic initiation factor 2 (eIF2) on the alpha-subunit as well as increased phosphorylation of the elongation factor (eEF2). Myotubes expressing a catalytically inactive PKR variant, PKRDelta6, showed no depression of protein synthesis in response to either LPS or TNF-alpha, confirming the importance of PKR in this process. There was no effect of any of the agents on phosphorylation of mammalian target of rapamycin (mTOR) or initiation factor 4E-binding protein (4E-BP1), and thus no change in the amount of eIF4E bound to 4E-BP1 or the concentration of the active eIF4E.eIF4G complex. HMB attenuated phosphorylation of eEF2, possibly by increasing phosphorylation of mTOR, and also attenuated phosphorylation of eIF2alpha by preventing activation of PKR. These results suggest that HMB may be effective in attenuating muscle atrophy in a range of catabolic conditions.

Peptide binding prediction for the human class II MHC allele HLA-DP2:a molecular docking approach

MHC class II proteins bind oligopeptide fragments derived from proteolysis of pathogen antigens, presenting them at the cell surface for recognition by CD4+ T cells. Human MHC class II alleles are grouped into three loci: HLA-DP, HLA-DQ and HLA-DR. In contrast to HLA-DR and HLA-DQ, HLA-DP proteins have not been studied extensively, as they have been viewed as less important in immune responses than DRs and DQs. However, it is now known that HLA-DP alleles are associated with many autoimmune diseases. Quite recently, the X-ray structure of the HLA-DP2 molecule (DPA*0103, DPB1*0201) in complex with a self-peptide derived from the HLA-DR a-chain has been determined. In the present study, we applied a validated molecular docking protocol to a library of 247 modelled peptide-DP2 complexes, seeking to assess the contribution made by each of the 20 naturally occurred amino acids at each of the nine binding core peptide positions and the four flanking residues (two on both sides).

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.

Static energy analysis of MHC class I and class II peptide-binding affinity

Antigenic peptide is presented to a T-cell receptor (TCR) through the formation of a stable complex with a major histocompatibility complex (MHC) molecule. Various predictive algorithms have been developed to estimate a peptide's capacity to form a stable complex with a given MHC class II allele, a technique integral to the strategy of vaccine design. These have previously incorporated such computational techniques as quantitative matrices and neural networks. A novel predictive technique is described, which uses molecular modeling of predetermined crystal structures to estimate the stability of an MHC class II-peptide complex. The structures are remodeled, energy minimized, and annealed before the energetic interaction is calculated.

Role of reactive oxygen species in protein degradation in murine myotubes induced by proteolysis-inducing factor and angiotensin II

The antioxidants butylated hydroxytoluene (BHT, 1 mM) and d-α-tocopherol (10 μM) completely attenuated protein degradation in murine myotubes in response to both proteolysis-inducing factor (PIF) and angiotensin II (Ang II), suggesting that the formation of reactive oxygen species (ROS) plays an important role in this process. Both PIF and Ang II induced a rapid and transient increase in ROS formation in myotubes, which followed a parabolic dose-response curve, similar to that for total protein degradation. Antioxidant treatment attenuated the increase in expression and activity of the ubiquitin-proteasome proteolytic pathway by PIF and Ang II, by preventing the activation of the transcription factor nuclear factor-κB (NF-κB), through inhibition of phosphorylation of the NF-κB inhibitor protein (I-κB) and its subsequent degradation. ROS formation by both PIF and Ang II was attenuated by diphenyleneiodonium (10 μM), suggesting that it was mediated through the NADPH oxidase system. ROS formation was also attenuated by trifluoroacetyl arachidonic acid (10 μM), a specific inhibitor of cytosolic phospholipase A2, U-73122 (5 μM) and D609 (200 μM), inhibitors of phospholipase C and calphostin C (300 nM), a highly specific inhibitor of protein kinase C (PKC), all known activators of NADPH oxidase. Myotubes containing a dominant-negative mutant of PKC did not show an increase in ROS formation in response to either PIF or Ang II. The two Rac1 inhibitors W56 (200 μM) and NSC23766 (10 μM) also attenuated both ROS formation and protein degradation induced by both PIF and Ang II. Rac1 is known to mediate signalling between the phosphatidylinositol-3 kinase (PI-3K) product and NADPH oxidase, and treatment with LY24002 (10 μM), a highly selective inhibitor of PI-3K, completely attenuated ROS production in response to both PIF and Ang II, and inhibited total protein degradation, while the inactive analogue LY303511 (100 μM) had no effect. ROS formation appears to be important in muscle atrophy in cancer cachexia, since treatment of weight losing mice bearing the MAC16 tumour with d-α-tocopherol (1 mg kg- 1) attenuated protein degradation and increased protein synthesis in skeletal muscle. © 2007 Elsevier Inc. All rights reserved.

Mechanism of induction of muscle protein degradation by angiotensin II

Angiotensin I and II have been shown to directly induce protein degradation in skeletal muscle through an increased activity and expression of the ubiquitin-proteasome proteolytic pathway. This investigation determines the role of the nuclear transcription factor nuclear factor-κB (NF-κB) in this process. Using murine myotubes as a surrogate model system both angiotensin I and II were found to induce activation of protein kinase C (PKC), with a parabolic dose-response curve similar to the induction of total protein degradation. Activation of PKC was required for the induction of proteasome expression, since calphostin C, a highly specific inhibitor of PKC, attenuated both the increase in total protein degradation and in proteasome expression and functional activity increased by angiotensin II. PKC is known to activate I-κB kinase (IKK), which is responsible for the phosphorylation and subsequent degradation of I-κB. Both angiotensin I and II induced an early decrease in cytoplasmic I-κB levels followed by nuclear accumulation of NF-κB. Using an NF-κB luciferase construct this was shown to increase transcriptional activation of NF-κB regulated genes. Maximal luciferase expression was seen at the same concentrations of angiotensin I/II as those inducing protein degradation. Total protein degradation induced by both angiotensin I and II was attenuated by resveratrol, which prevented nuclear accumulation of NF-κB, confirming that activation of NF-κB was responsible for the increased protein degradation. These results suggest that induction of proteasome expression by angiotensin I/II involves a signalling pathway involving PKC and NF-κB. © 2005 Elsevier Inc. All rights reserved.

Toward the prediction of class I and II mouse major histocompatibility complex-peptide-binding affinity:in silico bioinformatic step-by-step guide using quantitative structure-activity relationships

Quantitative structure-activity relationship (QSAR) analysis is a cornerstone of modern informatics. Predictive computational models of peptide-major histocompatibility complex (MHC)-binding affinity based on QSAR technology have now become important components of modern computational immunovaccinology. Historically, such approaches have been built around semiqualitative, classification methods, but these are now giving way to quantitative regression methods. We review three methods--a 2D-QSAR additive-partial least squares (PLS) and a 3D-QSAR comparative molecular similarity index analysis (CoMSIA) method--which can identify the sequence dependence of peptide-binding specificity for various class I MHC alleles from the reported binding affinities (IC50) of peptide sets. The third method is an iterative self-consistent (ISC) PLS-based additive method, which is a recently developed extension to the additive method for the affinity prediction of class II peptides. The QSAR methods presented here have established themselves as immunoinformatic techniques complementary to existing methodology, useful in the quantitative prediction of binding affinity: current methods for the in silico identification of T-cell epitopes (which form the basis of many vaccines, diagnostics, and reagents) rely on the accurate computational prediction of peptide-MHC affinity. We have reviewed various human and mouse class I and class II allele models. Studied alleles comprise HLA-A*0101, HLA-A*0201, HLA-A*0202, HLA-A*0203, HLA-A*0206, HLA-A*0301, HLA-A*1101, HLA-A*3101, HLA-A*6801, HLA-A*6802, HLA-B*3501, H2-K(k), H2-K(b), H2-D(b) HLA-DRB1*0101, HLA-DRB1*0401, HLA-DRB1*0701, I-A(b), I-A(d), I-A(k), I-A(S), I-E(d), and I-E(k). In this chapter we show a step-by-step guide into predicting the reliability and the resulting models to represent an advance on existing methods. The peptides used in this study are available from the AntiJen database (http://www.jenner.ac.uk/AntiJen). The PLS method is available commercially in the SYBYL molecular modeling software package. The resulting models, which can be used for accurate T-cell epitope prediction, will be made are freely available online at the URL http://www.jenner.ac.uk/MHCPred.

Oxidised LDL-lipids increase beta amyloid production by SH-SY5Y cells through glutathione depletion and lipid raft formation

Elevated total cholesterol in midlife has been associated with increased risk of dementia in later life. We have previously shown that low-density lipoprotein (LDL) is more oxidized in the plasma of dementia patients, although total cholesterol levels are not different from those of age-matched controls. β-Amyloid (Aβ) peptide, which accumulates in Alzheimer disease (AD), arises from the initial cleavage of amyloid precursor protein by β-secretase-1 (BACE1). BACE1 activity is regulated by membrane lipids and raft formation. Given the evidence for altered lipid metabolism in AD, we have investigated a mechanism for enhanced Aβ production by SH-SY5Y neuronal-like cells exposed to oxidized LDL (oxLDL). The viability of SH-SY5Y cells exposed to 4 μg oxLDL and 25 μM 27-hydroxycholesterol (27OH-C) was decreased significantly. Lipids, but not proteins, extracted from oxLDL were more cytotoxic than oxLDL. In parallel, the ratio of reduced glutathione (GSH) to oxidized glutathione was decreased at sublethal concentrations of lipids extracted from native and oxLDL. GSH loss was associated with an increase in acid sphingomyelinase (ASMase) activity and lipid raft formation, which could be inhibited by the ASMase inhibitor desipramine. 27OH-C and total lipids from LDL and oxLDL independently increased Aβ production by SH-SY5Y cells, and Aβ accumulation could be inhibited by desipramine and by N-acetylcysteine. These data suggest a mechanism whereby oxLDL lipids and 27OH-C can drive Aβ production by GSH depletion, ASMase-driven membrane remodeling, and BACE1 activation in neuronal cells. © 2014 The Authors.

Toward prediction of class II mouse major histocompatibility complex peptide binding affinity:in silico bioinformatic evaluation using partial least squares, a robust multivariate statistical technique

The accurate identification of T-cell epitopes remains a principal goal of bioinformatics within immunology. As the immunogenicity of peptide epitopes is dependent on their binding to major histocompatibility complex (MHC) molecules, the prediction of binding affinity is a prerequisite to the reliable prediction of epitopes. The iterative self-consistent (ISC) partial-least-squares (PLS)-based additive method is a recently developed bioinformatic approach for predicting class II peptide−MHC binding affinity. The ISC−PLS method overcomes many of the conceptual difficulties inherent in the prediction of class II peptide−MHC affinity, such as the binding of a mixed population of peptide lengths due to the open-ended class II binding site. The method has applications in both the accurate prediction of class II epitopes and the manipulation of affinity for heteroclitic and competitor peptides. The method is applied here to six class II mouse alleles (I-Ab, I-Ad, I-Ak, I-As, I-Ed, and I-Ek) and included peptides up to 25 amino acids in length. A series of regression equations highlighting the quantitative contributions of individual amino acids at each peptide position was established. The initial model for each allele exhibited only moderate predictivity. Once the set of selected peptide subsequences had converged, the final models exhibited a satisfactory predictive power. Convergence was reached between the 4th and 17th iterations, and the leave-one-out cross-validation statistical terms - q2, SEP, and NC - ranged between 0.732 and 0.925, 0.418 and 0.816, and 1 and 6, respectively. The non-cross-validated statistical terms r2 and SEE ranged between 0.98 and 0.995 and 0.089 and 0.180, respectively. The peptides used in this study are available from the AntiJen database (http://www.jenner.ac.uk/AntiJen). The PLS method is available commercially in the SYBYL molecular modeling software package. The resulting models, which can be used for accurate T-cell epitope prediction, will be made freely available online (http://www.jenner.ac.uk/MHCPred).