An integrated approach to epitope analysis I: Dimensional reduction, visualization and prediction of MHC binding using amino acid principal components and regression approaches

BACKGROUND: Operation of the immune system is multivariate. Reduction of the dimensionality is essential to facilitate understanding of this complex biological system. One multi-dimensional facet of the immune system is the binding of epitopes to the MHC-I and MHC-II molecules by diverse populations of individuals. Prediction of such epitope binding is critical and several immunoinformatic strategies utilizing amino acid substitution matrices have been designed to develop predictive algorithms. Contemporaneously, computational and statistical tools have evolved to handle multivariate and megavariate analysis, but these have not been systematically deployed in prediction of MHC binding. Partial least squares analysis, principal component analysis, and associated regression techniques have become the norm in handling complex datasets in many fields. Over two decades ago Wold and colleagues showed that principal components of amino acids could be used to predict peptide binding to cellular receptors. We have applied this observation to the analysis of MHC binding, and to derivation of predictive methods applicable on a whole proteome scale. RESULTS: We show that amino acid principal components and partial least squares approaches can be utilized to visualize the underlying physicochemical properties of the MHC binding domain by using commercially available software. We further show the application of amino acid principal components to develop both linear partial least squares and non-linear neural network regression prediction algorithms for MHC-I and MHC-II molecules. Several visualization options for the output aid in understanding the underlying physicochemical properties, enable confirmation of earlier work on the relative importance of certain peptide residues to MHC binding, and also provide new insights into differences among MHC molecules. We compared both the linear and non-linear MHC binding prediction tools to several predictive tools currently available on the Internet. CONCLUSIONS: As opposed to the highly constrained user-interaction paradigms of web-server approaches, local computational approaches enable interactive analysis and visualization of complex multidimensional data using robust mathematical tools. Our work shows that prediction tools such as these can be constructed on the widely available JMP(® )platform, can operate in a spreadsheet environment on a desktop computer, and are capable of handling proteome-scale analysis with high throughput

Patterns of Predicted T-Cell Epitopes Associated with Antigenic Drift in Influenza H3N2 Hemagglutinin

Antigenic drift allowing escape from neutralizing antibodies is an important feature of transmission and survival of influenza viruses in host populations. Antigenic drift has been studied in particular detail for influenza A H3N2 and well defined antigenic clusters of this virus documented. We examine how host immunogenetics contributes to determination of the antibody spectrum, and hence the immune pressure bringing about antigenic drift. Using uTOPE™ bioinformatics analysis of predicted MHC binding, based on amino acid physical property principal components, we examined the binding affinity of all 9-mer and 15-mer peptides within the hemagglutinin 1 (HA1) of 447 H3N2 virus isolates to 35 MHC-I and 14 MHC-II alleles. We provide a comprehensive map of predicted MHC-I and MHC-II binding affinity for a broad array of HLA alleles for the H3N2 influenza HA1 protein. Each HLA allele exhibited a characteristic predicted binding pattern. Cluster analysis for each HLA allele shows that patterns based on predicted MHC binding mirror those described based on antibody binding. A single amino acid mutation or position displacement can result in a marked difference in MHC binding and hence potential T-helper function. We assessed the impact of individual amino acid changes in HA1 sequences between 10 virus isolates from 1968–2002, representative of antigenic clusters, to understand the changes in MHC binding over time. Gain and loss of predicted high affinity MHC-II binding sites with cluster transitions were documented. Predicted high affinity MHC-II binding sites were adjacent to antibody binding sites. We conclude that host MHC diversity may have a major determinant role in the antigenic drift of influenza A H3N2

Effects of antiplatelet therapy on stroke risk by brain imaging features of intracerebral haemorrhage and cerebral small vessel diseases: subgroup analyses of the RESTART randomised, open-label trial

Background Findings from the RESTART trial suggest that starting antiplatelet therapy might reduce the risk of recurrent symptomatic intracerebral haemorrhage compared with avoiding antiplatelet therapy. Brain imaging features of intracerebral haemorrhage and cerebral small vessel diseases (such as cerebral microbleeds) are associated with greater risks of recurrent intracerebral haemorrhage. We did subgroup analyses of the RESTART trial to explore whether these brain imaging features modify the effects of antiplatelet therapy

Frequency patterns of T-cell exposed motifs in immunoglobulin heavy chain peptides presented by MHCs

Immunoglobulins are highly diverse protein sequences that are processed and presented to T-cells by B-cells and other antigen presenting cells. We examined a large dataset of immunoglobulin heavy chain variable regions (IGHV) to assess the diversity of T-cell exposed motifs (TCEM). TCEM comprise those amino acids in a MHC-bound peptide which face outwards, surrounded by the MHC histotope, and which engage the T-cell receptor. Within IGHV there is a distinct pattern of predicted MHC class II binding and a very high frequency of re-use of the TCEMs. The re-use frequency indicates that only a limited number of different cognate T-cells are required to engage many different clonal B-cells. The amino acids in each outward-facing TCEM are intercalated with the amino acids of inward-facing MHC groove-exposed motifs (GEM). Different GEM may have differing, allele-specific, MHC binding affinities. The intercalation of TCEM and GEM in a peptide allows for a vast combinatorial repertoire of epitopes, each eliciting a different response. Outcome of T-cell receptor binding is determined by overall signal strength, which is a function of the number of responding T-cells and the duration of engagement. Hence, the frequency of T-cell exposed motif re-use appears to be an important determinant of whether a T-cell response is stimulatory or suppressive. The frequency distribution of TCEMs implies that somatic hypermutation is followed by clonal expansion that develop along repeated pathways. The observations of TCEM and GEM derived from immunoglobulins suggest a relatively simple, yet powerful, mechanism to correlate T-cell polyspecificity, through re-use of TCEMs, with a very high degree of specificity achieved by combination with a diversity of GEMs. The frequency profile of TCEMs also points to an economical mechanism for maintaining T-cell memory, recall, and self-discrimination based on an endogenously generated profile of motifs

A Role for Epitope Networking in Immunomodulation by Helminths

Helminth infections, by nematodes, trematodes, or cestodes, can lead to the modulation of host immune responses. This allows long-duration parasite infections and also impacts responses to co-infections. Surface, secreted, excreted, and shed proteins are thought to play a major role in modulation. A commonly reported feature of such immune modulation is the role of T regulatory (Treg) cells and IL-10. Efforts to identify helminth proteins, which cause immunomodulation, have identified candidates but not provided clarity as to a uniform mechanism driving modulation. In this study, we applied a bioinformatics systems approach, allowing us to analyze predicted T-cell epitopes of 17 helminth species and the responses to their surface proteins. In addition to major histocompatibility complex (MHC) binding, we analyzed amino acid motifs that would be recognized by T-cell receptors [T-cell-exposed motifs (TCEMs)]. All the helminth species examined have, within their surface proteins, peptides, which combine very common TCEMs with predicted high affinity binding to many human MHC alleles. This combination of features would result in large cognate T cell and a high probability of eliciting Treg responses. The TCEMs, which determine recognition by responding T-cell clones, are shared to a high degree between helminth species and with Plasmodium falciparum and Mycobacterium tuberculosis, both common co-infecting organisms. The implication of our observations is not only that Treg cells play a significant role in helminth-induced immune modulation but also that the epitope specificities of Treg responses are shared across species and genera of helminth. Hence, the immune response to a given helminth cannot be considered in isolation but rather forms part of an epitope ecosystem, or microenvironment, in which potentially immunosuppressive peptides in the helminth network via their common T-cell receptor recognition signals with T-cell epitopes in self proteins, microbiome, other helminths, and taxonomically unrelated pathogens. Such a systems approach provides a high-level view of the antigen-immune system signaling dynamics that may bias a host’s immune response to helminth infections toward immune modulation. It may indicate how helminths have evolved to select for peptides that favor long-term parasite host coexistence

Recognition of Higher Order Patterns in Proteins: Immunologic Kernels

<div><p>By applying analysis of the principal components of amino acid physical properties we predicted cathepsin cleavage sites, MHC binding affinity, and probability of B-cell epitope binding of peptides in tetanus toxin and in ten diverse additional proteins. Cross-correlation of these metrics, for peptides of all possible amino acid index positions, each evaluated in the context of a ±25 amino acid flanking region, indicated that there is a strongly repetitive pattern of short peptides of approximately thirty amino acids each bounded by cathepsin cleavage sites and each comprising B-cell linear epitopes, MHC–I and MHC-II binding peptides. Such “immunologic kernel” peptides comprise all signals necessary for adaptive immunologic cognition, response and recall. The patterns described indicate a higher order spatial integration that forms a symbolic logic coordinating the adaptive immune system.</p></div

Data_Sheet_1_A Role for Epitope Networking in Immunomodulation by Helminths.PDF

<p>Helminth infections, by nematodes, trematodes, or cestodes, can lead to the modulation of host immune responses. This allows long-duration parasite infections and also impacts responses to co-infections. Surface, secreted, excreted, and shed proteins are thought to play a major role in modulation. A commonly reported feature of such immune modulation is the role of T regulatory (Treg) cells and IL-10. Efforts to identify helminth proteins, which cause immunomodulation, have identified candidates but not provided clarity as to a uniform mechanism driving modulation. In this study, we applied a bioinformatics systems approach, allowing us to analyze predicted T-cell epitopes of 17 helminth species and the responses to their surface proteins. In addition to major histocompatibility complex (MHC) binding, we analyzed amino acid motifs that would be recognized by T-cell receptors [T-cell-exposed motifs (TCEMs)]. All the helminth species examined have, within their surface proteins, peptides, which combine very common TCEMs with predicted high affinity binding to many human MHC alleles. This combination of features would result in large cognate T cell and a high probability of eliciting Treg responses. The TCEMs, which determine recognition by responding T-cell clones, are shared to a high degree between helminth species and with Plasmodium falciparum and Mycobacterium tuberculosis, both common co-infecting organisms. The implication of our observations is not only that Treg cells play a significant role in helminth-induced immune modulation but also that the epitope specificities of Treg responses are shared across species and genera of helminth. Hence, the immune response to a given helminth cannot be considered in isolation but rather forms part of an epitope ecosystem, or microenvironment, in which potentially immunosuppressive peptides in the helminth network via their common T-cell receptor recognition signals with T-cell epitopes in self proteins, microbiome, other helminths, and taxonomically unrelated pathogens. Such a systems approach provides a high-level view of the antigen-immune system signaling dynamics that may bias a host’s immune response to helminth infections toward immune modulation. It may indicate how helminths have evolved to select for peptides that favor long-term parasite host coexistence.</p

Conceptual model of an immunologic kernel.

<p>Relationships of the components are shown based on the cross correlations conducted. Two-headed arrows indicate there will be minor positional differences based on the host MHC alleles. Cathepsin cleavage is a requirement at the C-terminal of the MHC peptides; a high frequency of cathepsin cleavage occurs on the proximal side of the B-cell epitope, but no functional requirement for such cleavage has been demonstrated. Each cleavage site comprises an octomer, with the central P1-P1’ scissile bond indicated by the vertical arrow and the octomer amino acids shown as beads. We have characterized a kernel to comprise both B-cell epitope and T-cell epitope components; as shown T-independent and B-independent epitopes comprise subunits of the whole.</p

Predicted cleavage of tetanus toxin by human cathepsin L and S.

<p>A: Shows the distribution of the distance between successive cleavage probabilities of ≥0.5 for the two cathepsins. λ = expected value (mean) and σ = over dispersion (variance) of the fitted gamma Poisson distribution. B: Cross correlation of cleavage by cathepsin L and cathepsin S cleavage probabilities. A high correlation centered at zero indicates that the two cathepsins have a tendency to cut at the same site within the protein. This is flanked by high correlation of low probability of cleavage at ±5 amino acids of the initial cleavage. The red dashed lines at ±0.04 indicate the 95^th percentile confidence limits for non-significant correlation; values outside this band are significant p<0.05.</p