The organization of human leucocyte antigen class I epitopes in HIV genome products: implications for HIV evolution and vaccine design

Analysis of coevolution in nonstructural proteins of chikungunya virus

Background

RNA viruses are characterized by high rate of mutations mainly due to the lack of proofreading repair activities associated with its RNA-dependent RNA-polymerase (RdRp). In case of arboviruses, this phenomenon has lead to the existence of mixed population of genomic variants within the host called quasi-species. The stability of strains within the quasi-species lies on mutations that are positively selected which in turn depend on whether these mutations are beneficial in either or both hosts. Coevolution of amino acids (aa) is one phenomenon that leads to establishment of favorable traits in viruses and leading to their fitness.

Results

Fourteen CHIKV clinical samples collected over three years were subjected to RT-PCR, the four non-structural genes amplified and subjected to various genetic analyses. Coevolution analysis showed 30 aa pairs coevolving in nsP1, 23 aa pairs coevolving in nsP2, 239 in nsP3 and 46 aa coevolving pairs in nsP4 when each non-structural protein was considered independently. Further analysis showed that 705 amino acids pairs of the non-structural polyproteins coevolved together with a correlation coefficient of ≥0.5. Functional relevance of these coevolving amino acids in all the nonstructural proteins of CHIKV were predicted using Eukaryotic Linear Motifs (ELMs) of human.

Conclusions

The present study was undertaken to study co-evolving amino acids in the non-structural proteins of chikungunya virus (CHIKV), an important arbovirus. It was observed that several amino acids residues were coevolving and shared common functions.

Characterization of the Protective HIV-1 CTL Epitopes and the Corresponding HLA Class I Alleles: A Step towards Designing CTL Based HIV-1 Vaccine

Human immunodeficiency virus (HIV) possesses a major threat to the human life largely due to the unavailability of an efficacious vaccine and poor access to the antiretroviral drugs against this deadly virus. High mutation rate in the viral genome underlying the antigenic variability of the viral proteome is the major hindrance as far as the antibody based vaccine development is concerned. Although the exact mechanism by which CTL epitopes and the restricting HLA alleles mediate their action towards slow disease progression is still not clear, the important CTL restricted epitopes for controlling viral infections can be utilized in future vaccine design. This study was designed for the characterization the HIV-1 optimal CTL epitopes and their corresponding HLA alleles. CTL epitope cluster distribution analysis revealed only two HIV-1 proteins, namely, Nef and Gag, which have significant cluster forming capacity. We have found the role of specific HLA supertypes such as HLA B∗07, HLA B∗58, and HLA A∗03 in selecting the hydrophobic and conserved amino acid positions within Nef and Gag proteins, to be presented as epitopes. The analyses revealed that the clusters of optimal epitopes for Nef and p24 proteins of HIV-1 could potentially serve as a source of vaccine.

HLA supertypes contribute in HIV type 1 cytotoxic T lymphocyte epitope clustering in Nef and Gag proteins

Induction of HIV-1-specific cytotoxic T lymphocyte (CTL) responses largely depends upon the presentation of CTL epitopes to the CD8(+) T cells aided by a large number of different HLA class I alleles. Although several studies showed the clustering pattern of HIV-1 CTL epitopes, the underlying reason for this tendency remains unresolved. Moreover, the hypothesis that the CTL epitope clusters tend to coincide with the conserved and hydrophobic regions of HIV-1 proteins has been challenged in recent times. The present study aims to characterize and compare the HIV-1 CTL epitope clusters in terms of restricting HLA alleles, hydrophobicity, and sequence conservation in a proteome-wide manner by including a large number of experimentally validated CTL epitopes from the HIV Molecular Immunology Database. CTL epitope cluster distribution analysis in a proteome-wide manner revealed that only two HIV-1 proteins, namely Nef and Gag, have significant cluster-forming capacity where their epitope localization coincides with the hydrophobic and conserved regions. Furthermore, analyses of proteasomal cleavage sites and HLA anchoring motif frequencies in the epitope-dense regions highlighted the role of specific HLA supertypes such as HLA B*07, HLA B*58, HLA A*02, and HLA A*03 in selecting the hydrophobic and conserved amino acid positions within Nef and Gag proteins to be presented as epitopes. Based on our results, we hypothesize that the cluster-forming tendency of HIV-1 CTL epitopes is not a proteome-wide feature confined to Nef and Gag proteins. Their cluster-forming tendency largely depends on the host HLA alleles that contribute significantly in selecting functionally constrained hydrophobic regions within the HIV-1 proteome.

Conservation of HIV-1 T cell epitopes across time and clades: validation of immunogenic HLA-A2 epitopes selected for the GAIA HIV vaccine

HIV genomic sequence variability has complicated efforts to generate an effective globally relevant vaccine. Regions of the viral genome conserved in sequence and across time may represent the "Achilles' heel" of HIV. In this study, highly conserved T-cell epitopes were selected using immunoinformatics tools combining HLA-A2 supertype binding predictions with relative global conservation. Analysis performed in 2002 on 10,803 HIV-1 sequences, and again in 2009, on 43,822 sequences, yielded 38 HLA-A2 epitopes. These epitopes were experimentally validated for HLA binding and immunogenicity with PBMCs from HIV-infected patients in Providence, Rhode Island, and/or Bamako, Mali. Thirty-five (92%) stimulated an IFNγ response in PBMCs from at least one subject. Eleven of fourteen peptides (79%) were confirmed as HLA-A2 epitopes in both locations. Validation of these HLA-A2 epitopes conserved across time, clades, and geography supports the hypothesis that such epitopes could provide effective coverage of virus diversity and would be appropriate for inclusion in a globally relevant HIV vaccine.

Further progress on defining highly conserved immunogenic epitopes for a global HIV vaccine: HLA-A3-restricted GAIA vaccine epitopes

Two major obstacles confronting HIV vaccine design have been the extensive viral diversity of HIV-1 globally and viral evolution driven by escape from CD8(+) cytotoxic T-cell lymphocyte (CTL)-mediated immune pressure. Regions of the viral genome that are not able to escape immune response and that are conserved in sequence and across time may represent the "Achilles' heel" of HIV and would be excellent candidates for vaccine development. In this study, T-cell epitopes were selected using immunoinformatics tools, combining HLA-A3 binding predictions with relative sequence conservation in the context of global HIV evolution. Twenty-seven HLA-A3 epitopes were chosen from an analysis performed in 2003 on 10,803 HIV-1 sequences, and additional sequences were selected in 2009 based on an expanded set of 43,822 sequences. These epitopes were tested in vitro for HLA binding and for immunogenicity with PBMCs of HIV-infected donors from Providence, Rhode Island. Validation of these HLA-A3 epitopes conserved across time, clades, and geography supports the hypothesis that epitopes such as these would be candidates for inclusion in our globally relevant GAIA HIV vaccine constructs.

Identification of immunogenic HLA-B7 "Achilles' heel" epitopes within highly conserved regions of HIV

Genetic polymorphisms in class I human leukocyte antigen molecules (HLA) have been shown to determine susceptibility to HIV infection as well as the rate of progression to AIDS. In particular, the HLA-B7 supertype has been shown to be associated with high viral loads and rapid progression to disease. Using a multiplatform in silico/in vitro approach, we have prospectively identified 45 highly conserved, putative HLA-B7 restricted HIV CTL epitopes and evaluated them in HLA binding and ELISpot assays. All 45 epitopes (100%) bound to HLA-B7 in cell-based HLA binding assays: 28 (62%) bound with high affinity, 6 (13%) peptides bound with medium affinity and 11 (24%) bound with low affinity. Forty of the 45 peptides (88%) stimulated a IFN-gamma response in PBMC from at least one subject. Eighteen of these 40 epitopes have not been previously described; an additional eight epitopes have not been previously described as restricted by B7. The HLA-B7 restricted epitopes discovered using this in silico screening approach are highly conserved across strains and clades of HIV as well as conserved in the HIV genome over the 20 years since HIV-1 isolates were first sequenced. This study demonstrates that it is possible to select a broad range of HLA-B7 restricted epitopes that comprise stable elements in the rapidly mutating HIV genome. The most immunogenic of these epitopes will be included in the GAIA multi-epitope vaccine.

Structural prediction of peptides binding to MHC class I molecules

Peptide binding to class I major histocompatibility complex (MHCI) molecules is a key step in the immune response and the structural details of this interaction are of importance in the design of peptide vaccines. Algorithms based on primary sequence have had success in predicting potential antigenic peptides for MHCI, but such algorithms have limited accuracy and provide no structural information. Here, we present an algorithm, PePSSI (peptide-MHC prediction of structure through solvated interfaces), for the prediction of peptide structure when bound to the MHCI molecule, HLA-A2. The algorithm combines sampling of peptide backbone conformations and flexible movement of MHC side chains and is unique among other prediction algorithms in its incorporation of explicit water molecules at the peptide-MHC interface. In an initial test of the algorithm, PePSSI was used to predict the conformation of eight peptides bound to HLA-A2, for which X-ray data are available. Comparison of the predicted and X-ray conformations of these peptides gave RMSD values between 1.301 and 2.475 A. Binding conformations of 266 peptides with known binding affinities for HLA-A2 were then predicted using PePSSI. Structural analyses of these peptide-HLA-A2 conformations showed that peptide binding affinity is positively correlated with the number of peptide-MHC contacts and negatively correlated with the number of interfacial water molecules. These results are consistent with the relatively hydrophobic binding nature of the HLA-A2 peptide binding interface. In summary, PePSSI is capable of rapid and accurate prediction of peptide-MHC binding conformations, which may in turn allow estimation of MHCI-peptide binding affinity.

Evolution of nef variants in gut associated lymphoid tissue of rhesus macaques during primary simian immunodeficiency virus infection

We utilized the simian immunodeficiency virus model of AIDS to examine evolution of nef gene in gut-associated lymphoid tissue (GALT) during primary and early asymptomatic stages of infection. Macaques were infected with a cloned virus, SIVmac239/nef-stop harboring a premature stop codon in the nef gene. Restoration of the nef open reading frame occurred in GALT early at 3 days post-infection. Analysis of nef sequences by phylogenetic tools showed that evolution of nef was neutral thereafter, as evidenced by the ratio of synonymous to nonsynonymous substitutions, a star pattern in unrooted trees and distribution of amino acid replacements fitting a simple Poisson process. Two regions encoding for a nuclear localization signal and a CTL epitope were conserved. Thus, GALT was a site for strong positive selection of functional nef during initial stages of infection. However, evolution of the nef gene thereafter was neutral during early asymptomatic stage of infection.

Studies in macaques on cross-clade T cell responses elicited by a DNA/MVA AIDS vaccine, better conservation of CD8 than CD4 T cell responses

One of the unknowns faced by an HIV/AIDS vaccine is the ability of a single clade vaccine to protect against the multiple genetic subtypes and recombinant forms of HIV-1 present in the current pandemic. Here, we use a macaque model to investigate the ability of our clade B vaccine that consists of DNA priming and modified vaccinia Ankara (MVA) virus boosting to elicit T cell responses that recognize an A/G recombinant of HIV-1. To test for cross-reactive T cells, intracellular cytokine staining was conducted using five pools of Gag and six pools of Env peptides representing B or A/G sequences. Studies using the peptide pools revealed essentially complete conservation of the CD8 response but only approximately 50% conservation of the CD4 response. Thus, the ability of an HIV vaccine for one clade to protect against other clades may be more limited by the ability to provide CD4 T cell help than the ability to elicit CD8 effector functions.

Differential proteasomal processing of hydrophobic and hydrophilic protein regions: contribution to cytotoxic T lymphocyte epitope clustering in HIV-1-Nef

HIV proteins contain a multitude of naturally processed cytotoxic T lymphocyte (CTL) epitopes that concentrate in clusters. The molecular basis of epitope clustering is of interest for understanding HIV immunogenicity and for vaccine design. We show that the CTL epitope clusters of HIV proteins predominantly coincide with hydrophobic regions, whereas the noncluster regions are predominantly hydrophilic. Analysis of the proteasomal degradation products of full-length HIV-Nef revealed a differential sensitivity of cluster and noncluster regions to proteasomal processing. Compared with the epitope-scarce noncluster regions, cluster regions are digested by proteasomes more intensively and with greater preference for hydrophobic P1 residues, resulting in substantially greater numbers of fragments with the sizes and COOH termini typical of epitopes and their precursors. Indeed, many of these fragments correspond to endogenously processed Nef epitopes and/or their potential precursors. The results suggest that differential proteasomal processing contributes importantly to the clustering of CTL epitopes in hydrophobic regions.

Magnitude and frequency of cytotoxic T-lymphocyte responses: identification of immunodominant regions of human immunodeficiency virus type 1 subtype C

A systematic analysis of immune responses on a population level is critical for a human immunodeficiency virus type 1 (HIV-1) vaccine design. Our studies in Botswana on (i) molecular analysis of the HIV-1 subtype C (HIV-1C) epidemic, (ii) frequencies of major histocompatibility complex class I HLA types, and (iii) cytotoxic T-lymphocyte (CTL) responses in the course of natural infection allowed us to address HIV-1C-specific immune responses on a population level. We analyzed the magnitude and frequency of the gamma interferon ELISPOT-based CTL responses and translated them into normalized cumulative CTL responses. The introduction of population-based cumulative CTL responses reflected both (i) essentials of the predominant virus circulating locally in Botswana and (ii) specificities of the genetic background of the Botswana population, and it allowed the identification of immunodominant regions across the entire HIV-1C. The most robust and vigorous immune responses were found within the HIV-1C proteins Gag p24, Vpr, Tat, and Nef. In addition, moderately strong responses were scattered across Gag p24, Pol reverse transcriptase and integrase, Vif, Tat, Env gp120 and gp41, and Nef. Assuming that at least some of the immune responses are protective, these identified immunodominant regions could be utilized in designing an HIV vaccine candidate for the population of southern Africa. Targeting multiple immunodominant regions should improve the overall vaccine immunogenicity in the local population and minimize viral escape from immune recognition. Furthermore, the analysis of HIV-1C-specific immune responses on a population level represents a comprehensive systematic approach in HIV vaccine design and should be considered for other HIV-1 subtypes and/or different geographic areas.

Clustering patterns of cytotoxic T-lymphocyte epitopes in human immunodeficiency virus type 1 (HIV-1) proteins reveal imprints of immune evasion on HIV-1 global variation

The human cytotoxic T-lymphocyte (CTL) response to human immunodeficiency virus type 1 (HIV-1) has been intensely studied, and hundreds of CTL epitopes have been experimentally defined, published, and compiled in the HIV Molecular Immunology Database. Maps of CTL epitopes on HIV-1 protein sequences reveal that defined epitopes tend to cluster. Here we integrate the global sequence and immunology databases to systematically explore the relationship between HIV-1 amino acid sequences and CTL epitope distributions. CTL responses to five HIV-1 proteins, Gag p17, Gag p24, reverse transcriptase (RT), Env, and Nef, have been particularly well characterized in the literature to date. Through comparing CTL epitope distributions in these five proteins to global protein sequence alignments, we identified distinct characteristics of HIV amino acid sequences that correlate with CTL epitope localization. First, experimentally defined HIV CTL epitopes are concentrated in relatively conserved regions. Second, the highly variable regions that lack epitopes bear cumulative evidence of past immune escape that may make them relatively refractive to CTLs: a paucity of predicted proteasome processing sites and an enrichment for amino acids that do not serve as C-terminal anchor residues. Finally, CTL epitopes are more highly concentrated in alpha-helical regions of proteins. Based on amino acid sequence characteristics, in a blinded fashion, we predicted regions in HIV regulatory and accessory proteins that would be likely to contain CTL epitopes; these predictions were then validated by comparison to new sets of experimentally defined epitopes in HIV-1 Rev, Tat, Vif, and Vpr.

The quest for an AIDS vaccine: is the CD8+ T-cell approach feasible?

The rationale for developing anti-HIV vaccines that stimulate cytotoxic T-lymphocyte responses is given. We argue that such vaccines will work, provided that attention is paid to the development of memory T-cell responses that are strong and preferably activated. Furthermore, the vaccine should match the prevailing virus clade as closely as possible. Vaccines will have to stimulate a wide range of responses, but it is not clear how this can be achieved.

Definition of the Mamu A*01 peptide binding specificity: application to the identification of wild-type and optimized ligands from simian immunodeficiency virus regulatory proteins

Single amino acid substitution analogs of the known Mamu A*01 binding peptide gag 181-190 and libraries of naturally occurring sequences of viral or bacterial origin were used to rigorously define the peptide binding motif associated with Mamu A*01 molecules. The presence of S or T in position 2, P in position 3, and hydrophobic or aromatic residues at the C terminus is associated with optimal binding capacity. At each of these positions, additional residues are also tolerated but associated with significant decreases in binding capacity. The presence of at least two preferred and one tolerated residues at the three anchor positions is necessary for good Mamu A*01 binding; optimal ligand size is 8-9 residues. This detailed motif has been used to map potential epitopes from SIVmac239 regulatory proteins and to engineer peptides with increased binding capacity. A total of 13 wild type and 17 analog candidate epitopes were identified. Furthermore, our analysis reveals a significantly lower than expected frequency of epitopes in early regulatory proteins, suggesting a possible evolutionary- and/or immunoselection directed against variants of viral products that contain CTL epitopes.

The effect of human beta2-microglobulin on major histocompatibility complex I peptide loading and the engineering of a high affinity variant. Implications for peptide-based vaccines

The ability to directly load cell surface major histocompatibility complex (MHC) class I molecules with peptides provides a potentially powerful approach toward the development of vaccines to generate cell-mediated immunity. We demonstrate that exogenous beta2-microglobulin (beta2m) stabilizes human cell surface MHC I molecules and facilitates their loading with exogenous peptides. Additionally, using three-dimensional crystal structures and known interaction sites between MHC I heavy chains and beta2m, we engineered variants of human beta2m (hbeta2m) with a single serine substitution at residue 55. This alteration was predicted to promote hydrophobic interactions at the MHC I heavy chain/beta2m interface and displace an ordered water molecule. Compared with hbeta2m, the serine to valine substitution at residue 55 had improved ability to bind to cell surface HLA-A1, HLA-A2, and HLA-A3 molecules, facilitate exogenous peptide loading, and promote recognition by peptide-specific T cells. The inclusion of hbeta2m or higher affinity variants when pulsing cells with MHC-restricted peptides increases the efficiency of peptide loading 50-80-fold. Therefore, the inclusion of hbeta2m in peptide-based vaccines may increase cell surface antigen densities above thresholds that allow recognition of peptide antigens by the immune system, particularly for cryptic, subdominant, or marginally antigenic peptides.

Structural principles that govern the peptide-binding motifs of class I MHC molecules

The peptides that bind class I MHC molecules are restricted in length and often contain key amino acids, anchor residues, at particular positions. The side-chains of peptide anchor residues interact with the polymorphic complementary pockets in MHC peptide-binding grooves and provide the molecular basis for allele-specific recognition of antigenic peptides. We establish correlations between class I MHC specificities for anchor residues and class I MHC sequence markers that occur at the polymorphic positions lining the structural pockets. By analyzing the pocket structures of nine crystallized class I MHC molecules and the modeled structures of another 39 class I MHC molecules, we show that class I pockets can be classified into families that are distinguishable by their common physico-chemical properties and peptide side-chain selectivities. The identification of recurrent structural principles among class I pockets makes it possible to greatly expand the repertoire of known peptide-binding motifs of class I MHC molecules. The evolutionary strategies underlying the emergence of pocket families is briefly discussed.