HLA-DP2 binding prediction by molecular dynamics simulations

Major histocompatibility complex (MHC) II proteins bind peptide fragments derived from pathogen antigens and present them at the cell surface for recognition by T cells. MHC proteins are divided into Class I and Class II. Human MHC Class II alleles are grouped into three loci: HLA-DP, HLA-DQ, and HLA-DR. They are involved in many autoimmune diseases. In contrast to HLA-DR and HLA-DQ proteins, the X-ray structure of the HLA-DP2 protein has been solved quite recently. In this study, we have used structure-based molecular dynamics simulation to derive a tool for rapid and accurate virtual screening for the prediction of HLA-DP2-peptide binding. A combinatorial library of 247 peptides was built using the "single amino acid substitution" approach and docked into the HLA-DP2 binding site. The complexes were simulated for 1 ns and the short range interaction energies (Lennard-Jones and Coulumb) were used as binding scores after normalization. The normalized values were collected into quantitative matrices (QMs) and their predictive abilities were validated on a large external test set. The validation shows that the best performing QM consisted of Lennard-Jones energies normalized over all positions for anchor residues only plus cross terms between anchor-residues.     

Membrane specializations and endosome maturation in dendritic cells and B cells

Interest in the cell biology of antigen presentation is centered on dendritic cells (DCs) as initiators of the immune response. The ability to examine primary antigen-presenting cells, as opposed to cell lines, has opened a new window for study of antigen processing and peptide acquisition by Class II major histocompatibility complex (MHC) products, especially where intracellular trafficking of peptide-Class-II complexes is concerned. Here, we review the dynamics of Class II MHC-positive intracellular structures in dendritic cells as well as B cells. We focus on the generation of multivesicular bodies, where Class II MHC products acquire antigenic peptide, on the endosomal transport of peptide-loaded Class II MHC to the cell surface and on the importance of Class II MHC localization in membrane microdomains.     

Empty class II major histocompatibility complex created by peptide photolysis establishes the role of DM in peptide association

DM catalyzes the exchange of peptides bound to Class II major histocompatibility complex (MHC) molecules. Because the dissociation and association components of the overall reaction are difficult to separate, a detailed mechanism of DM catalysis has long resisted elucidation. UV irradiation of DR molecules loaded with a photocleavable peptide (caged Class II MHC molecules) enabled synchronous and verifiable evacuation of the peptide-binding groove and tracking of early binding events in real time by fluorescence polarization. Empty DR molecules generated by photocleavage rapidly bound peptide but quickly resolved into species with substantially slower binding kinetics. DM formed a complex with empty DR molecules that bound peptide with even faster kinetics than empty DR molecules just having lost their peptide cargo. Mathematical models demonstrate that the peptide association rate of DR molecules is substantially higher in the presence of DM. We therefore unequivocally establish that DM contributes directly to peptide association through formation of a peptide-loading complex between DM and empty Class II MHC. This complex rapidly acquires a peptide analogous to the MHC class I peptide-loading complex.     

Peptide binding inhibits aggregation of soluble MHC class II in solution

Affinity-purified major histocompatability complex (MHC) class II molecules are known to bind antigenic peptide in vitro. This peptide-bound MHC class II is known to undergo a change in structure upon stable binding of antigenic peptide. Previous results from our, and other laboratories, have suggested a relationship between MHC class II structure and peptide association that enables class II to enter into a stable conformation upon peptide binding. In this report we describe that stable binding of high-affinity antigenic peptide to MHC class II molecule results in transition of aggregated purified MHC class II proteins to a stable heterodimeric state. Such transition was demonstrated by using purified human HLA-DR2 class II molecule and high-affinity myelin basic protein (MBP) 83-102)Y83 peptide. Highly aggregated purified DR2 (high molecular weight; HMW) was first separated from heterodimer (low molecular weight: LMW) in the presence of 50-fold molar excess of MBP(83-102)Y83 peptide. We then show that the aggregated HMW preparation can be successfully converted into a stable dimer by further incubation with MBP(83-102)Y83 and changing various binding parameters such as pH, temperature, reducing agent, and peptide concentrations. Under optimized conditions, the highly aggregated inactive DR2 molecules can be completely loaded with the antigenic peptide. The transformed heterodimers with bound peptide prepared by this method are biologically active, as shown by their ability to induce the production of gamma-interferon by SS8T-transformed human T cells. These results suggest that in solution, MHC class II molecules may be aggregated in the absence of bound peptide. Such aggregated MHC class II molecules can be converted to stable and biologically active heterodimers in the presence of high-affinity antigenic peptide.     

Fluorogenic probes for monitoring peptide binding to class II MHC proteins in living cells

A crucial step in the immune response is the binding of antigenic peptides to major histocompatibility complex (MHC) proteins. Class II MHC proteins present their bound peptides to CD4(+) T cells, thereby helping to activate both the humoral and the cellular arms of the adaptive immune response. Peptide loading onto class II MHC proteins is regulated temporally, spatially and developmentally in antigen-presenting cells. To help visualize these processes, we have developed a series of novel fluorogenic probes that incorporate the environment-sensitive amino acid analogs 6-N,N-dimethylamino-2-3-naphthalimidoalanine and 4-N,N-dimethylaminophthalimidoalanine. Upon binding to class II MHC proteins these fluorophores show large changes in emission spectra, quantum yield and fluorescence lifetime. Peptides incorporating these fluorophores bind specifically to class II MHC proteins on antigen-presenting cells and can be used to follow peptide binding in vivo. Using these probes we have tracked a developmentally regulated cell-surface peptide-binding activity in primary human monocyte-derived dendritic cells.     

Kinetics of peptide binding to the class II MHC protein I-Ek

Class II MHC glycoproteins bind short (7-25 amino acid) peptides in an extended type II polyproline-like conformation and present them for immune recognition. Because empty MHC is unstable, measurement of the rate of the second-order reaction between peptide and MHC is challenging. In this report, we use dissociation of a pre-bound peptide to generate the active, peptide-receptive form of the empty class II MHC molecule I-Ek. This allows us to measure directly the rate of reaction between active, empty I-Ek and a set of peptides that vary in structure. We find that all peptides studied, despite having highly variable dissociation rates, bind with similar association rate constants. Thus, the rate-limiting step in peptide binding is minimally sensitive to peptide side-chain structure. An interesting complication to this simple model is that a single peptide can sometimes bind to I-Ek in two kinetically distinguishable conformations, with the stable peptide-MHC complex isomer forming much more slowly than the less-stable one. This demonstrates that an additional free-energy barrier limits the formation of certain specific MHC-peptide complex conformations.     

High variability in the MHC class II DA beta chain of the brushtail possum (<Emphasis Type="Italic">Trichosurus vulpecula</Emphasis>)

The diversity of class II major histocompatibility complex (MHC) loci was investigated in the brushtail possum, an important marsupial pest species in New Zealand. Immunocontraception, a form of fertility control that generates an autoimmune response, is being developed as a population control method for the possum. Because the immune response is partly under genetic control, an understanding of immunogenetics in possum will be crucial to the development of immunocontraceptive vaccines. MHC molecules are critical in the vertebrate immune response. Class II MHC molecules bind and present exogenously derived peptides to T lymphocytes and may be important in the presentation of immunocontraceptives. We used polymerase chain reaction primers designed to amplify the peptide binding region of possum class II MHC genes to isolate sequences from 49 animals. We have previously described 19 novel alleles from the DAB locus in the possum (Holland et al., Immunogenetics 60:449–460, 2008). Here, we report on another 11 novel alleles isolated from possum DAB, making this the most diverse marsupial locus described so far. This high level of diversity indicates that DAB is an important MHC locus in the possum and will need to be taken into account in the design of immunocontraceptive vaccines.     

Inhibition of endosomal proteolytic activity by leupeptin blocks surface expression of MHC class II molecules and their conversion to SDS resistance alpha beta heterodimers in endosomes.

The biosynthesis of MHC Class II molecules starts with the assembly of the alpha and beta subunits and the invariant chain. Intracellular transport of Class II molecules was followed in pulse-chase experiments of a human Epstein-Barr virus-transformed B lymphoblastoid cell line. Entry of Class II molecules into the endocytotic pathway and their cell surface appearance were monitored using neuraminidase as a fluid endocytotic marker and as a surface probe, respectively. In the course of intracellular transport, the Class II associated invariant chain is removed by proteases located in the endosomal pathway. Here, we show that leupeptin inhibits not only invariant chain breakdown, but also surface deposition of newly synthesized Class II molecules. Class II molecules display remarkable resistance to SDS at ambient temperature when occupied by peptide. We exploit this property to show that peptide binding precedes surface expression, and takes place in the course of intracellular transport through an endosomal compartment. Leupeptin blocks the conversion of Class II molecules to an SDS resistant complex.     

Engineering protein for X-ray crystallography: the murine Major Histocompatibility Complex class II molecule I-Ad.

Class II Major Histocompatibility (MHC) molecules are cell surface heterodimeric glycoproteins that play a central role in the immune response by presenting peptide antigens for surveillance by T cells. Due to the inherent instability of the class II MHC heterodimer, and its dependence on bound peptide for proper assembly, the production of electrophoretically pure samples of class II MHC proteins in complex with specific peptides has been problematic. A soluble form of the murine class II MHC molecule, I-Ad, with a leucine zipper tail added to each chain to enhance dimer assembly and secretion, has been produced in Drosophila melanogaster SC2 cells. To facilitate peptide loading, a high affinity ovalbumin peptide was covalently engineered to be attached by a six-residue linker to the amino terminus of the I-Adbeta chain. This modified I-Ad molecule was purified using preparative IEF and one fraction, after removal of the leucine zipper tails, produced crystals suitable for X-ray crystallographic analysis. The protein engineering and purification methods described here should be of general value for the expression of I-A and other class II MHC-peptide complexes.     

Prediction of MHC class II binding affinity using SMM-align,a novel stabilization matrix alignment method

Background  

Antigen presenting cells (APCs) sample the extra cellular space and present peptides from here to T helper cells, which can be activated if the peptides are of foreign origin. The peptides are presented on the surface of the cells in complex with major histocompatibility class II (MHC II) molecules. Identification of peptides that bind MHC II molecules is thus a key step in rational vaccine design and developing methods for accurate prediction of the peptide:MHC interactions play a central role in epitope discovery. The MHC class II binding groove is open at both ends making the correct alignment of a peptide in the binding groove a crucial part of identifying the core of an MHC class II binding motif. Here, we present a novel stabilization matrix alignment method, SMM-align, that allows for direct prediction of peptide:MHC binding affinities. The predictive performance of the method is validated on a large MHC class II benchmark data set covering 14 HLA-DR (human MHC) and three mouse H2-IA alleles.     

Stable peptide binding to MHC class II molecule is rapid and is determined by a receptive conformation shaped by prior association with low affinity peptides

Formation of stable class II MHC/peptide complex involves conformational changes and proceeds via an intermediate. Although this intermediate complex forms and dissociates in minutes, its conversion to a stable complex is a very slow process, taking up to a few days to reach completion. Here, we investigate the different steps of this binding and demonstrate that the conformational changes necessary to generate a receptive molecule is the rate-determining slow step in the process, while formation of the stable MHC/peptide complex is very rapid. With HLA-DR1 as our model class II molecule, we first used low affinity variants of hemagglutinin peptide (HA306-318), which lack the principal anchor, to shape the conformation of the MHC and then studied the kinetics of stable binding of HA306-318 to such an induced conformation. We found that the apparent association rate of HA306-318 is equivalent to the dissociation rate of the low affinity peptide. A 4- to 18-fold enhancement in the binding rates of HA306-318 was observed depending on the dissociation rates of the low affinity peptides. These results establish that 1) formation of stable MHC/peptide complexes is very rapid and 2) prior binding of low affinity peptide induces a receptive conformation in MHC for efficient stable peptide binding. Furthermore, in the absence of any free peptide, this receptive molecule rapidly reverts to slow binding behavior toward the subsequently offered peptide. These results have important implications for the roles of low affinity MHC/peptide complexes in Ag presentation.     

Endosomal sorting of MHC class II determines antigen presentation by dendritic cells

Dendritic cells (DCs) initiate primary immune responses by presenting pathogen-derived antigens in association with major histocompatibility Class II molecules (MHC II) to T cells. In DCs, MHC II is constitutively synthesized and loaded at endosomes with peptides from hydrolyzed endogenous proteins or exogenously acquired antigens. Whether peptide loaded MHC II (MHC II-p) is subsequently recruited to and stably expressed at the plasma membrane or degraded in lysosomes is determined by the status of the DC. In immature DCs, MHC II-p is ubiquitinated after peptide loading, driving its sorting to the luminal vesicles of multivesicular bodies. These luminal vesicles, and the MHC II-p they carry, are delivered to lysosomes for degradation. MHC II-p is inefficiently ubiquitinated in DCs that are activated by pathogens or inflammatory stimuli, thus allowing its transfer to and stable expression at the plasma membrane.     

Prediction of MHC binding peptide using Gibbs motif sampler,weight matrix and artificial neural network

The identification of MHC restricted epitopes is an important goal in peptide based vaccine and diagnostic development. As wet lab experiments for identification of MHC binding peptide are expensive and time consuming, in silico tools have been developed as fast alternatives, however with low performance. In the present study, we used IEDB training and blind validation datasets for the prediction of peptide binding to fourteen human MHC class I and II molecules using Gibbs motif sampler, weight matrix and artificial neural network methods. As compare to MHC class I predictor based on sequence weighting (Aroc=0.95 and CC=0.56) and artificial neural network (Aroc=0.73 and CC=0.25), MHC class II predictor based on Gibbs sampler did not perform well (Aroc=0.62 and CC=0.19). The predictive accuracy of Gibbs motif sampler in identifying the 9-mer cores of a binding peptide to DRB1 alleles are also limited (40¢), however above the random prediction (14¢). Therefore, the size of dataset (training and validation) and the correct identification of the binding core are the two main factors limiting the performance of MHC class-II binding peptide prediction. Overall, these data suggest that there is substantial room to improve the quality of the core predictions using novel approaches that capture distinct features of MHC-peptide interactions than the current approaches.     

Retroviral antibody binding of the MHC class II molecule: a biochemical influence on CD4 T cell differentiation in HIV infection?

Retroviral antibody capable of binding to the major histocompatibility complex (MHC) Class II molecule has been documented in human immunodeficiency virus-1 (HIV-1)-infected patients. Interactions between the MHC Class II receptor and the T-cell receptor (TCR) are central to the immune response. Importantly, retroviral antibody possesses a much higher binding affinity for the MHC Class II receptor, when compared to the TCR. Experiments have manipulated a number of factors related to antigen-presenting cell (APC) interaction with differentiating T-cells. These studies have observed the effects of lowering antigen dose and reducing ligand density on precursor Th (T helper) cell differentiation. Studies have also examined the effect of downregulated MHC Class II receptors and co-stimulatory molecules on APC-Th cell interaction. In addition, the sequestration of antigens away from the Class II processing pathway has been studied. These investigations reveal a general trend that can determine whether a naive CD4 T-cell becomes a Th1 or Th2-like cell. If the experimental manipulation weakens the APC-Th cell interaction, a weak ligating TCR signal results. Consequently, a weak ligating TCR signal can influence precursor Th cells to become Th2-like cells. Retroviral antibody binding of MHC Class II receptors may mimic a number of experimental conditions responsible for creating a weak ligating TCR signal.     

Improved prediction of MHC class I and class II epitopes using a novel Gibbs sampling approach

MOTIVATION: Prediction of which peptides will bind a specific major histocompatibility complex (MHC) constitutes an important step in identifying potential T-cell epitopes suitable as vaccine candidates. MHC class II binding peptides have a broad length distribution complicating such predictions. Thus, identifying the correct alignment is a crucial part of identifying the core of an MHC class II binding motif. In this context, we wish to describe a novel Gibbs motif sampler method ideally suited for recognizing such weak sequence motifs. The method is based on the Gibbs sampling method, and it incorporates novel features optimized for the task of recognizing the binding motif of MHC classes I and II. The method locates the binding motif in a set of sequences and characterizes the motif in terms of a weight-matrix. Subsequently, the weight-matrix can be applied to identifying effectively potential MHC binding peptides and to guiding the process of rational vaccine design. RESULTS: We apply the motif sampler method to the complex problem of MHC class II binding. The input to the method is amino acid peptide sequences extracted from the public databases of SYFPEITHI and MHCPEP and known to bind to the MHC class II complex HLA-DR4(B1*0401). Prior identification of information-rich (anchor) positions in the binding motif is shown to improve the predictive performance of the Gibbs sampler. Similarly, a consensus solution obtained from an ensemble average over suboptimal solutions is shown to outperform the use of a single optimal solution. In a large-scale benchmark calculation, the performance is quantified using relative operating characteristics curve (ROC) plots and we make a detailed comparison of the performance with that of both the TEPITOPE method and a weight-matrix derived using the conventional alignment algorithm of ClustalW. The calculation demonstrates that the predictive performance of the Gibbs sampler is higher than that of ClustalW and in most cases also higher than that of the TEPITOPE method.     

Kinetics of registry selection of chimeric peptides binding to MHC II

Major histocompatability complex type II proteins (MHC II) are alphabeta-heterodimeric glycoproteins that present peptides to the T cell receptor (TCR) of CD4(+) T-cells. This presentation may result in activation of these T-cells, depending on the nature of the peptide. Peptides interact specifically with MHC II with nine peptide amino acid positions, and the corresponding MHC II pocket positions are usually labeled P1-P9. However, the length of peptides binding to MHC II may be greater than nine amino acids, and therefore these peptides may potentially bind to the MHC II in more than one registry. To investigate the mechanism by which a long peptide binds to I-E(k), a murine MHC II, a chimeric peptide with two nonoverlapping registries, f-IAYLKQATKQLRMATPLLMR was designed. The IAYLKQATK peptide segment is based on moth cytochrome c 95-103 (MCC 95-103), and the QLRMATPLLMR segment is based on murine Ii CLIP 89-99 M90L (Ii CLIP 89-99 M90L). This chimeric peptide forms two isomeric complexes. The MCC and Ii CLIP registries dissociate from I-E(k) with t(1/2) values of >800 and 4.94 h, respectively. The registry composition of this MHC II/chimeric peptide complex was found to change as a function of time in approaching thermodynamic equilibrium: the results are consistent with a kinetic model that involves no intramolecular isomer interconversion. The model depicts uncorrelated binding to the MHC II determined by relative association rates to the two registries. This is followed by dissociation and subsequent rebinding, leading ultimately to a preponderance of the most stable complex. Similar results were obtained at pH 5.3. The behavior of this chimeric peptide approximates the binding of a 1:1 solution mixture of two peptides to MHC II, where the more stable complex is selected over time. We have also found that a chimeric peptide and a human MHC II, HLA-DR40401, form isomers with relative association rates to DR0401 at pH 5.3 of 15% for one isomer and 85% for the second isomer.     

Flexibility of the MHC class II peptide binding cleft in the bound, partially filled, and empty states: a molecular dynamics simulation study

Major histocompatibility (MHC) Class II cell surface proteins present antigenic peptides to the immune system. Class II structures in complex with peptides but not in the absence of peptide are known. Comparative molecular dynamics (MD) simulations of a Class II protein (HLA-DR3) with and without CLIP (invariant chain-associated protein) peptide were performed starting from the CLIP-bound crystal structure. Depending on the protonation of acidic residues in the P6 peptide-binding pocket the simulations stayed overall close to the start structure. The simulations without CLIP showed larger conformational fluctuations especially of alpha-helices flanking the binding cleft. Largest fluctuations without CLIP were observed in a helical segment near the peptide C-terminus binding region matching a segment recognized by antibodies specific for empty Class II proteins. Simulations on a Val86Tyr mutation that fills the peptide N-terminus binding P1 pocket or of a complex with a CLIP fragment (dipeptide) bound to P1 showed an unexpected long range effect. In both simulations the mobility not only of P1 but also of the entire binding cleft was reduced compared to simulations without CLIP. It correlates with the experimental finding that the CLIP fragment binding to P1 is sufficient to prevent antibody recognition specific for the empty form at a site distant from P1. The results suggest a mechanism how a local binding event of small peptides or of an exchange factor near P1 may promote peptide binding and exchange through a long range stabilization of the whole binding cleft in a receptive (near bound) conformation.     

Reconstructing an ancestral mammalian immune supercomplex from a marsupial major histocompatibility complex

The first sequenced marsupial genome promises to reveal unparalleled insights into mammalian evolution. We have used theMonodelphis domestica (gray short-tailed opossum) sequence to construct the first map of a marsupial major histocompatibility complex (MHC). The MHC is the most gene-dense region of the mammalian genome and is critical to immunity and reproductive success. The marsupial MHC bridges the phylogenetic gap between the complex MHC of eutherian mammals and the minimal essential MHC of birds. Here we show that the opossum MHC is gene dense and complex, as in humans, but shares more organizational features with non-mammals. The Class I genes have amplified within the Class II region, resulting in a unique Class I/II region. We present a model of the organization of the MHC in ancestral mammals and its elaboration during mammalian evolution. The opossum genome, together with other extant genomes, reveals the existence of an ancestral “immune supercomplex” that contained genes of both types of natural killer receptors together with antigen processing genes and MHC genes.     

Monoclonal antibodies specific for the empty conformation of HLA-DR1 reveal aspects of the conformational change associated with peptide binding

Class II major histocompatibility complex (MHC) proteins bind peptides and present them at the cell surface for interaction with CD4+ T cells as part of the system by which the immune system surveys the body for signs of infection. Peptide binding is known to induce conformational changes in class II MHC proteins on the basis of a variety of hydrodynamic and spectroscopic approaches, but the changes have not been clearly localized within the overall class II MHC structure. To map the peptide-induced conformational change for HLA-DR1, a common human class II MHC variant, we generated a series of monoclonal antibodies recognizing the beta subunit that are specific for the empty conformation. Each antibody reacted with the empty but not the peptide-loaded form, for both soluble recombinant protein and native protein expressed at the cell surface. Antibody binding epitopes were characterized using overlapping peptides and alanine scanning substitutions and were localized to two distinct regions of the protein. The pattern of key residues within the epitopes suggested that the two epitope regions undergo substantial conformational alteration during peptide binding. These results illuminate aspects of the structure of the empty forms and the nature of the peptide-induced conformational change.