New horizons in mouse immunoinformatics:reliable in silico prediction of mouse class I histocompatibility major complex peptide binding affinity

Quantitative structure–activity relationship (QSAR) analysis is a main cornerstone of modern informatic disciplines. Predictive computational models, based on QSAR technology, of peptide-major histocompatibility complex (MHC) binding affinity have now become a vital component of modern day computational immunovaccinology. Historically, such approaches have been built around semi-qualitative, classification methods, but these are now giving way to quantitative regression methods. The additive method, an established immunoinformatics technique for the quantitative prediction of peptide–protein affinity, was used here to identify the sequence dependence of peptide binding specificity for three mouse class I MHC alleles: H2–Db, H2–Kb and H2–Kk. As we show, in terms of reliability the resulting models represent a significant advance on existing methods. They can be used for the accurate prediction of T-cell epitopes and are freely available online (http://www.jenner.ac.uk/MHCPred).

SVRMHC prediction server for MHC-binding peptides

The binding between antigenic peptides (epitopes) and the MHC molecule is a key step in the cellular immune response. Accurate in silico prediction of epitope-MHC binding affinity can greatly expedite epitope screening by reducing costs and experimental effort. Recently, we demonstrated the appealing performance of SVRMHC, an SVR-based quantitative modeling method for peptide-MHC interactions, when applied to three mouse class I MHC molecules. Subsequently, we have greatly extended the construction of SVRMHC models and have established such models for more than 40 class I and class II MHC molecules. Here we present the SVRMHC web server for predicting peptide-MHC binding affinities using these models. Benchmarked percentile scores are provided for all predictions. The larger number of SVRMHC models available allowed for an updated evaluation of the performance of the SVRMHC method compared to other well- known linear modeling methods. SVRMHC is an accurate and easy-to-use prediction server for epitope-MHC binding with significant coverage of MHC molecules. We believe it will prove to be a valuable resource for T cell epitope researchers.

In silico prediction of peptide-MHC binding affinity using SVRMHC

The binding between peptide epitopes and major histocompatibility complex (MHC) proteins is a major event in the cellular immune response. Accurate prediction of the binding between short peptides and class I or class II MHC molecules is an important task in immunoinformatics. SVRMHC which is a novel method to model peptide-MHC binding affinities based on support rector machine regression (SVR) is described in this chapter. SVRMHC is among a small handful of quantitative modeling methods that make predictions about precise binding affinities between a peptide and an MHC molecule. As a kernel-based learning method, SVRMHC has rendered models with demonstrated appealing performance in the practice of modeling peptide-MHC binding.

Bayesian classification with Gaussian processes

We consider the problem of assigning an input vector class='mathrm'>bfx to one of class='mathrm'>m classes by predicting class='mathrm'>P(c|bfx) for class='mathrm'>c = 1, ldots, m. For a two-class problem, the probability of class 1 given class='mathrm'>bfx is estimated by class='mathrm'>s(y(bfx)), where class='mathrm'>s(y) = 1/(1 + e^-y). A Gaussian process prior is placed on class='mathrm'>y(bfx), and is combined with the training data to obtain predictions for new class='mathrm'>bfx points. We provide a Bayesian treatment, integrating over uncertainty in class='mathrm'>y and in the parameters that control the Gaussian process prior; the necessary integration over class='mathrm'>y is carried out using Laplace's approximation. The method is generalized to multi-class problems class='mathrm'>(m >2) using the softmax function. We demonstrate the effectiveness of the method on a number of datasets.

The clinical significance of medicines reconciliation in children admitted to hospital

Aims and Objectives: The NICE/NPSA guidance on Medicines Reconciliation in adults upon hospital admission excludes children under the age of 16.1 Hence the primary aim and objective of this study was to use medicines reconciliation to primarily identify if discrepancies occur upon hospital admission. Secondary objectives were to clinically assess for harm discrepancies that were identified in paediatric patients on long term medications at four hospitals across the UK. Method: Medicines reconciliation is a procedure where the current medication history of a patient prior to hospital admission would be taken and verifying the medication orders made at hospital admission against this history, addressing any discrepancies identified. Medicines reconciliation was carried out prospectively for 244 paediatric patients on chronic medication across four UK hospitals (Birmingham, London, Leeds and North Staffordshire) between January – May 2011. Medicines reconciliation was conducted by a clinical pharmacist using the following sources of information: 1) the patient's Pre-Admission Medication (PAM) from the patient's general practitioner 2) examination of the Patient's Own Medications brought into hospital, 3) a semi-structured interview with the parent-carers and 4) identification of admission medication orders written on the drug chart prior to clinical pharmacy input (Drug Chart). Discrepancies between the PAM and Drug Chart were documented and classified as intentional or unintentional. Intentional discrepancies were defined as changes that were made knowingly by the prescriber and confirmed. Unintentional discrepancies were assessed for clinical significance by an expert panel and assigned a significance score based on the likelihood of causing potential discomfort or clinical deterioration: class 1 unlikely, class 2 moderate and class 3 severe.2 Results: 1004 medication regimens were included from the 244 patients across the four sites. 588 of the 1004 (59%) medicines, had discrepancies between the PAM and Drug Chart; of these 36% (n = 209) were unintentional and included for clinically assessment. 189 drug discrepancies 30% were classified as class 1, 47% were class 2 and 23% were class 3 discrepancies. The remaining 20 discrepancies were cases where deviating from the PAM would have been the right thing to do, which might suggest that an intentional but undocumented discrepancy by the prescriber writing up the admission order may have occurred. Conclusion: The results suggest that medication discrepancies in paediatric patients do occur upon hospital admission, which do have a potential to cause harm and that medicines reconciliation is a potential solution to preventing such discrepancies. References: 1. National Institute for Health and Clinical Excellence. National Patient Safety Agency. PSG001. Technical patient safety solutions for medicines reconciliation on admission of adults to hospital. London: NICE; 2007. 2. Cornish, P. L., Knowles, S. R., Marchesano, et al. Unintended Medication Discrepancies at the Time of Hospital Admission. Archives of Internal Medicine 2005; 165:424–429

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.

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).

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.

Predicting class I major histocompatibility complex (MHC) binders using multivariate statistics:comparison of discriminant analysis and multiple linear regression

The accurate in silico identification of T-cell epitopes is a critical step in the development of peptide-based vaccines, reagents, and diagnostics. It has a direct impact on the success of subsequent experimental work. Epitopes arise as a consequence of complex proteolytic processing within the cell. Prior to being recognized by T cells, an epitope is presented on the cell surface as a complex with a major histocompatibility complex (MHC) protein. A prerequisite therefore for T-cell recognition is that an epitope is also a good MHC binder. Thus, T-cell epitope prediction overlaps strongly with the prediction of MHC binding. In the present study, we compare discriminant analysis and multiple linear regression as algorithmic engines for the definition of quantitative matrices for binding affinity prediction. We apply these methods to peptides which bind the well-studied human MHC allele HLA-A*0201. A matrix which results from combining results of the two methods proved powerfully predictive under cross-validation. The new matrix was also tested on an external set of 160 binders to HLA-A*0201; it was able to recognize 135 (84%) of them.

Quantitative prediction of mouse class I MHC peptide binding affinity using support vector machine regression (SVR) models

Background - The binding between peptide epitopes and major histocompatibility complex proteins (MHCs) is an important event in the cellular immune response. Accurate prediction of the binding between short peptides and the MHC molecules has long been a principal challenge for immunoinformatics. Recently, the modeling of MHC-peptide binding has come to emphasize quantitative predictions: instead of categorizing peptides as "binders" or "non-binders" or as "strong binders" and "weak binders", recent methods seek to make predictions about precise binding affinities. Results - We developed a quantitative support vector machine regression (SVR) approach, called SVRMHC, to model peptide-MHC binding affinities. As a non-linear method, SVRMHC was able to generate models that out-performed existing linear models, such as the "additive method". By adopting a new "11-factor encoding" scheme, SVRMHC takes into account similarities in the physicochemical properties of the amino acids constituting the input peptides. When applied to MHC-peptide binding data for three mouse class I MHC alleles, the SVRMHC models produced more accurate predictions than those produced previously. Furthermore, comparisons based on Receiver Operating Characteristic (ROC) analysis indicated that SVRMHC was able to out-perform several prominent methods in identifying strongly binding peptides. Conclusion - As a method with demonstrated performance in the quantitative modeling of MHC-peptide binding and in identifying strong binders, SVRMHC is a promising immunoinformatics tool with not inconsiderable future potential.

MHC class I bound to an immunodominant Theileria parva epitope demonstrates unconventional presentation to T cell receptors

T cell receptor (TCR) recognition of peptide-MHC class I (pMHC) complexes is a crucial event in the adaptive immune response to pathogens. Peptide epitopes often display a strong dominance hierarchy, resulting in focusing of the response on a limited number of the most dominant epitopes. Such T cell responses may be additionally restricted by particular MHC alleles in preference to others. We have studied this poorly understood phenomenon using Theileria parva, a protozoan parasite that causes an often fatal lymphoproliferative disease in cattle. Despite its antigenic complexity, CD8+ T cell responses induced by infection with the parasite show profound immunodominance, as exemplified by the Tp1(214-224) epitope presented by the common and functionally important MHC class I allele N*01301. We present a high-resolution crystal structure of this pMHC complex, demonstrating that the peptide is presented in a distinctive raised conformation. Functional studies using CD8+ T cell clones show that this impacts significantly on TCR recognition. The unconventional structure is generated by a hydrophobic ridge within the MHC peptide binding groove, found in a set of cattle MHC alleles. Extremely rare in all other species, this feature is seen in a small group of mouse MHC class I molecules. The data generated in this analysis contribute to our understanding of the structural basis for T cell-dependent immune responses, providing insight into what determines a highly immunogenic p-MHC complex, and hence can be of value in prediction of antigenic epitopes and vaccine design.

A novel requirement for C. elegans Alix/ALX-1 in RME-1-mediated membrane transport

BACKGROUND: Alix/Bro1p family proteins have recently been identified as important components of multivesicular endosomes (MVEs) and are involved in the sorting of endocytosed integral membrane proteins, interacting with components of the ESCRT complex, the unconventional phospholipid LBPA, and other known endocytosis regulators. During infection, Alix can be co-opted by enveloped retroviruses, including HIV, providing an important function during virus budding from the plasma membrane. In addition, Alix is associated with the actin cytoskeleton and might regulate cytoskeletal dynamics. RESULTS: Here we demonstrate a novel physical interaction between the only apparent Alix/Bro1p family protein in C. elegans, ALX-1, and a key regulator of receptor recycling from endosomes to the plasma membrane, called RME-1. The analysis of alx-1 mutants indicates that ALX-1 is required for the endocytic recycling of specific basolateral cargo in the C. elegans intestine, a pathway previously defined by the analysis of rme-1 mutants. The expression of truncated human Alix in HeLa cells disrupts the recycling of major histocompatibility complex class I, a known Ehd1/RME-1-dependent transport step, suggesting the phylogenetic conservation of this function. We show that the interaction of ALX-1 with RME-1 in C. elegans, mediated by RME-1/YPSL and ALX-1/NPF motifs, is required for this recycling process. In the C. elegans intestine, ALX-1 localizes to both recycling endosomes and MVEs, but the ALX-1/RME-1 interaction appears to be dispensable for ALX-1 function in MVEs and/or late endosomes. CONCLUSIONS: This work provides the first demonstration of a requirement for an Alix/Bro1p family member in the endocytic recycling pathway in association with the recycling regulator RME-1.

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.

A hydrogen-bonded polar network in the core of the glucagon-like peptide-1 receptor is a fulcrum for biased agonism:lessons from class B crystal structuress

The glucagon-like peptide 1 (GLP-1) receptor is a class B G protein-coupled receptor (GPCR) that is a key target for treatments for type II diabetes and obesity. This receptor, like other class B GPCRs, displays biased agonism, though the physiologic significance of this is yet to be elucidated. Previous work has implicated R2.60190 , N3.43240 , Q7.49394 , and H6.52363 as key residues involved in peptide-mediated biased agonism, with R2.60190 , N3.43240 , and Q7.49394 predicted to form a polar interaction network. In this study, we used novel insight gained from recent crystal structures of the transmembrane domains of the glucagon and corticotropin releasing factor 1 (CRF1) receptors to develop improved models of the GLP-1 receptor that predict additional key molecular interactions with these amino acids. We have introduced E6.53364 A, N3.43240 Q, Q7.49493N, and N3.43240 Q/Q7.49 Q/Q7.49493N mutations to probe the role of predicted H-bonding and charge-charge interactions in driving cAMP, calcium, or extracellular signal-regulated kinase (ERK) signaling. A polar interaction between E6.53364 and R2.60190 was predicted to be important for GLP-1- and exendin-4-, but not oxyntomodulin-mediated cAMP formation and also ERK1/2 phosphorylation. In contrast, Q7.49394 , but not R2.60190 /E6.53364 was critical for calcium mobilization for all three peptides. Mutation of N3.43240 and Q7.49394 had differential effects on individual peptides, providing evidence for molecular differences in activation transition. Collectively, this work expands our understanding of peptide-mediated signaling from the GLP-1 receptor and the key role that the central polar network plays in these events.