Allosteric crosstalk between peptide-binding, transport, and ATP hydrolysis of the ABC transporter TAP

The transporter associated with antigen processing (TAP) is essential for intracellular transport of protein fragments into the endoplasmic reticulum for loading of major histocompatibility complex (MHC) class I molecules. On the cell surface, these peptide–MHC complexes are monitored by cytotoxic T lymphocytes. To study the ATP hydrolysis of TAP, we developed an enrichment and reconstitution procedure, by which we fully restored TAP function in proteoliposomes. A TAP-specific ATPase activity was identified that could be stimulated by peptides and blocked by the herpes simplex virus protein ICP47. Strikingly, the peptide-binding motif of TAP directly correlates with the stimulation of the ATPase activity, demonstrating that the initial peptide-binding step is responsible for TAP selectivity. ATP hydrolysis follows Michaelis–Menten kinetics with a maximal velocity Vmax of 2 μmol/min per mg TAP, corresponding to a turnover number of approximately 5 ATP per second. This turnover rate is sufficient to account for the role of TAP in peptide loading of MHC molecules and the overall process of antigen presentation. Interestingly, sterically restricted peptides that bind but are not transported by TAP do not stimulate ATPase activity. These results point to coordinated dialogue between the peptide-binding site, the nucleotide-binding domain, and the translocation site via conformational changes within the TAP complex.

A Novel Quality Control Compartment Derived from the Endoplasmic Reticulum

Degradation of proteins that, because of improper or suboptimal processing, are retained in the endoplasmic reticulum (ER) involves retrotranslocation to reach the cytosolic ubiquitin-proteasome machinery. We found that substrates of this pathway, the precursor of human asialoglycoprotein receptor H2a and free heavy chains of murine class I major histocompatibility complex (MHC), accumulate in a novel preGolgi compartment that is adjacent to but not overlapping with the centrosome, the Golgi complex, and the ER-to-Golgi intermediate compartment (ERGIC). On its way to degradation, H2a associated increasingly after synthesis with the ER translocon Sec61. Nevertheless, it remained in the secretory pathway upon proteasomal inhibition, suggesting that its retrotranslocation must be tightly coupled to the degradation process. In the presence of proteasomal inhibitors, the ER chaperones calreticulin and calnexin, but not BiP, PDI, or glycoprotein glucosyltransferase, concentrate in the subcellular region of the novel compartment. The “quality control” compartment is possibly a subcompartment of the ER. It depends on microtubules but is insensitive to brefeldin A. We discuss the possibility that it is also the site for concentration and retrotranslocation of proteins that, like the mutant cystic fibrosis transmembrane conductance regulator, are transported to the cytosol, where they form large aggregates, the “aggresomes.”

Human cytomegalovirus inhibits antigen presentation by a sequential multistep process.

The human cytomegalovirus (HCMV) genomic unique short (US) region encodes a family of homologous genes essential for the inhibition of major histocompatibility complex (MHC) class I-mediated antigen presentation during viral infection. Here we show that US3, the only immediate early (IE) gene within the US region, encodes an endoplasmic reticulum-resident glycoprotein that prevents intracellular transport of MHC class I molecules. In contrast to the rapid degradation of newly synthesized MHC class I heavy chains mediated by the early gene product US11, we found that US3 retains stable MHC class I heterodimers in the endoplasmic reticulum that are loaded with peptides while retained in the ER. Consistent with the expression pattern of US3 and US11, MHC class I molecules are retained but not degraded during the IE period of infection. Our data identify the first nonregulatory role of an IE protein of HCMV and suggest that HCMV uses different T-cell escape strategies at different times during the infectious cycle.

Induction of cytotoxic T lymphocytes by intramuscular immunization with plasmid DNA is facilitated by bone marrow-derived cells.

Striated muscle is the predominant site of gene expression after i.m. immunization of plasmid DNA, but it is not clear if myocytes or professional antigen-presenting cells (APCs) of hematopoietic origin present the encoded antigens to class I major histocompatibility complex (MHC)-restricted cytotoxic T lymphocytes (CTL). To address this issue, CTL responses were assessed in mice engrafted with immune systems that were partially MHC matched with antigen-producing muscle cells. Spleen cells (sc) from immunocompetent F1 H-2bxd mice were infused into H-2b or H-2d mice carrying the severe combined immunodeficiency (scid) mutation, creating F1sc-->H-2b and F1sc-->H-2d chimeras, respectively. Immunization with DNA plasmids encoding the herpes simplex virus gB or the human immunodeficiency virus gp120 glycoproteins elicited antiviral CTL activity. F1sc-->H-2d chimeras responded to an H-2d-restricted gp120 epitope but not an H-2b restricted gB epitope, whereas F1sc-->H-2b chimeras responded to the H-2b but not the H-2d restricted epitope. This pattern of epitope recognition by the sc chimeras indicated that APCs of recipient (scid) origin were involved in initiation of CTL responses. Significantly, CTL responses against epitopes presented by the mismatched donor class I molecules were elicited if F1 bone marrow cells and sc were transferred into scid recipients before or several days to weeks after DNA immunization. Thus, bone marrow-derived APCs are sufficient for class I MHC presentation of viral antigens after i.m. immunization with plasmid DNA. Expression of plasmid DNA by these APCs is probably not a requirement for CTL priming. Instead, they appear to present proteins synthesized by other host cells.

Unexpected beta2-microglobulin sequence diversity in individual rainbow trout.

For mammals beta2-microglobulin (beta2m), the light chain of major histocompatibility complex (MHC) class I molecules, is invariant (or highly conserved) and is encoded by a single gene unlinked to the MHC. We find that beta2m of a salmonid fish, the rainbow trout (Oncorhynchus mykiss), does not conform to the mammalian paradigm. Ten of 12 randomly selected beta2m cDNA clones from an individual fish have different nucleotide sequences. A complex restriction fragment length polymorphism pattern is observed with rainbow trout, suggesting multiple beta2m genes in the genome, in excess of the two genes expected from the ancestral salmonid tetraploidy. Additional duplication and diversification of the beta2m genes might have occurred subsequently. Variation in the beta2m cDNA sequences is mainly at sites that do not perturb the structure of the mature beta2m protein, showing that the observed diversity of the trout beta2m genes is not primarily a result of pathogen selection.

Immunization with a circumsporozoite epitope fused to Bordetella pertussis adenylate cyclase in conjunction with cytotoxic T-lymphocyte-associated antigen 4 blockade confers protection against Plasmodium berghei liver-stage malaria

The adenylate cyclase toxoid (ACT) of Bordetella pertussis is capable of delivering its N-terminal catalytic domain into the cytosol of CD11b-expressing professional antigen-presenting cells such as myeloid dendritic cells. This allows delivery of CD8+ T-cell epitopes to the major histocompatibility complex (MHC) class I presentation pathway. Recombinant detoxified ACT containing an epitope of the Plasmodium berghei circumsporozoite protein (CSP), indeed, induced a specific CD8+ T-cell response in immunized mice after a single application, as detected by MHC multimer staining and gamma interferon (IFN-gamma) ELISPOT assay. This CSP-specific response could be significantly enhanced by prime-boost immunization with recombinant ACT in combination with anti-CTLA-4 during the boost immunization. This increased response was accompanied by complete protection in a number of mice after a challenge with P. berghei sporozoites. Transient blockade of CTLA-4 may overcome negative regulation and hence provide a strategy to enhance the efficacy of a vaccine by amplifying the number of responding T cells.

Understanding the phylogeographic patterns of European hedgehogs, Erinaceus concolor and E-europaeus using the MHC

The genome of the European hedgehog, Erinaceus concolor and E. europaeus, shows a strong signal of cycles of restriction to glacial refugia and postglacial expansion. Patterns of expansion, however, differ for mitochondrial DNA (mtDNA) and preliminary analysis of nuclear markers. In this study, we determine phylogeographic patterns in the hedgehog using two loci of the major histocompatibility complex (MHC), isolated for the first time in hedgehogs. These genes show long persistence times and high polymorphism in many species because of the actions of balancing selection. Among 84 individuals screened for variation, only two DQA alleles were identified in each species, but 10 DQB alleles were found in E. concolor and six in E. europaeus. A strong effect of demography on patterns of DQB variability is observed, with only weak evidence of balancing selection. While data from mtDNA clearly subdivide both species into monophyletic subgroups, the MHC data delineate only E. concolor into distinct subgroups, supporting the preliminary findings of other nuclear markers. Together with differences in variability, this suggests that the refugia history and/or expansion patterns of E. concolor and E. europaeus differ.

Novel visualization methods for protein data

Visualization of high-dimensional data has always been a challenging task. Here we discuss and propose variants of non-linear data projection methods (Generative Topographic Mapping (GTM) and GTM with simultaneous feature saliency (GTM-FS)) that are adapted to be effective on very high-dimensional data. The adaptations use log space values at certain steps of the Expectation Maximization (EM) algorithm and during the visualization process. We have tested the proposed algorithms by visualizing electrostatic potential data for Major Histocompatibility Complex (MHC) class-I proteins. The experiments show that the variation in the original version of GTM and GTM-FS worked successfully with data of more than 2000 dimensions and we compare the results with other linear/nonlinear projection methods: Principal Component Analysis (PCA), Neuroscale (NSC) and Gaussian Process Latent Variable Model (GPLVM).

Integrating in silico and in vitro analysis of peptide binding affinity to HLA-Cw*0102:a bioinformatic approach to the prediction of new epitopes

Predictive models of peptide-Major Histocompatibility Complex (MHC) binding affinity are important components of modern computational immunovaccinology. Here, we describe the development and deployment of a reliable peptide-binding prediction method for a previously poorly-characterized human MHC class I allele, HLA-Cw*0102.

Toward an atomistic understanding of the immune synapse:large-scale molecular dynamics simulation of a membrane-embedded TCR-pMHC-CD4 complex

T-cell activation requires interaction of T-cell receptors (TCR) with peptide epitopes bound by major histocompatibility complex (MHC) proteins. This interaction occurs at a special cell-cell junction known as the immune or immunological synapse. Fluorescence microscopy has shown that the interplay among one agonist peptide-MHC (pMHC), one TCR and one CD4 provides the minimum complexity needed to trigger transient calcium signalling. We describe a computational approach to the study of the immune synapse. Using molecular dynamics simulation, we report here on a study of the smallest viable model, a TCR-pMHC-CD4 complex in a membrane environment. The computed structural and thermodynamic properties are in fair agreement with experiment. A number of biomolecules participate in the formation of the immunological synapse. Multi-scale molecular dynamics simulations may be the best opportunity we have to reach a full understanding of this remarkable supra-macromolecular event at a cell-cell junction.

Toward the atomistic simulation of T cell epitopes automated construction of MHC:peptide structures for free energy calculations

Epitopes mediated by T cells lie at the heart of the adaptive immune response and form the essential nucleus of anti-tumour peptide or epitope-based vaccines. Antigenic T cell epitopes are mediated by major histocompatibility complex (MHC) molecules, which present them to T cell receptors. Calculating the affinity between a given MHC molecule and an antigenic peptide using experimental approaches is both difficult and time consuming, thus various computational methods have been developed for this purpose. A server has been developed to allow a structural approach to the problem by generating specific MHC:peptide complex structures and providing configuration files to run molecular modelling simulations upon them. A system has been produced which allows the automated construction of MHC:peptide structure files and the corresponding configuration files required to execute a molecular dynamics simulation using NAMD. The system has been made available through a web-based front end and stand-alone scripts. Previous attempts at structural prediction of MHC:peptide affinity have been limited due to the paucity of structures and the computational expense in running large scale molecular dynamics simulations. The MHCsim server (http://igrid-ext.cryst.bbk.ac.uk/MHCsim) allows the user to rapidly generate any desired MHC:peptide complex and will facilitate molecular modelling simulation of MHC complexes on an unprecedented scale.

Multi-level visualisation using Gaussian process latent variable models

Projection of a high-dimensional dataset onto a two-dimensional space is a useful tool to visualise structures and relationships in the dataset. However, a single two-dimensional visualisation may not display all the intrinsic structure. Therefore, hierarchical/multi-level visualisation methods have been used to extract more detailed understanding of the data. Here we propose a multi-level Gaussian process latent variable model (MLGPLVM). MLGPLVM works by segmenting data (with e.g. K-means, Gaussian mixture model or interactive clustering) in the visualisation space and then fitting a visualisation model to each subset. To measure the quality of multi-level visualisation (with respect to parent and child models), metrics such as trustworthiness, continuity, mean relative rank errors, visualisation distance distortion and the negative log-likelihood per point are used. We evaluate the MLGPLVM approach on the ‘Oil Flow’ dataset and a dataset of protein electrostatic potentials for the ‘Major Histocompatibility Complex (MHC) class I’ of humans. In both cases, visual observation and the quantitative quality measures have shown better visualisation at lower levels.

Large-scale molecular dynamics simulations of HLA-A*0201 complexed with a tumor-specific antigenic peptide:can the α3 and β2m domains be neglected?

Large-scale massively parallel molecular dynamics (MD) simulations of the human class I major histo-compatibility complex (MHC) protein HLA-A*0201 bound to a decameric tumor-specific antigenic peptide GVY-DGREHTV were performed using a scalable MD code on high-performance computing platforms. Such computational capabilities put us in reach of simulations of various scales and complexities. The supercomputing resources available Large-scale massively parallel molecular dynamics (MD) simulations of the human class I major histocompatibility complex (MHC) protein HLA-A*0201 bound to a decameric tumor-specific antigenic peptide GVYDGREHTV were performed using a scalable MD code on high-performance computing platforms. Such computational capabilities put us in reach of simulations of various scales and complexities. The supercomputing resources available for this study allow us to compare directly differences in the behavior of very large molecular models; in this case, the entire extracellular portion of the peptide–MHC complex vs. the isolated peptide binding domain. Comparison of the results from the partial and the whole system simulations indicates that the peptide is less tightly bound in the partial system than in the whole system. From a detailed study of conformations, solvent-accessible surface area, the nature of the water network structure, and the binding energies, we conclude that, when considering the conformation of the α1–α2 domain, the α3 and β2m domains cannot be neglected. © 2004 Wiley Periodicals, Inc. J Comput Chem 25: 1803–1813, 2004

Visualisation of bioinformatics datasets

Analysing the molecular polymorphism and interactions of DNA, RNA and proteins is of fundamental importance in biology. Predicting functions of polymorphic molecules is important in order to design more effective medicines. Analysing major histocompatibility complex (MHC) polymorphism is important for mate choice, epitope-based vaccine design and transplantation rejection etc. Most of the existing exploratory approaches cannot analyse these datasets because of the large number of molecules with a high number of descriptors per molecule. This thesis develops novel methods for data projection in order to explore high dimensional biological dataset by visualising them in a low-dimensional space. With increasing dimensionality, some existing data visualisation methods such as generative topographic mapping (GTM) become computationally intractable. We propose variants of these methods, where we use log-transformations at certain steps of expectation maximisation (EM) based parameter learning process, to make them tractable for high-dimensional datasets. We demonstrate these proposed variants both for synthetic and electrostatic potential dataset of MHC class-I. We also propose to extend a latent trait model (LTM), suitable for visualising high dimensional discrete data, to simultaneously estimate feature saliency as an integrated part of the parameter learning process of a visualisation model. This LTM variant not only gives better visualisation by modifying the project map based on feature relevance, but also helps users to assess the significance of each feature. Another problem which is not addressed much in the literature is the visualisation of mixed-type data. We propose to combine GTM and LTM in a principled way where appropriate noise models are used for each type of data in order to visualise mixed-type data in a single plot. We call this model a generalised GTM (GGTM). We also propose to extend GGTM model to estimate feature saliencies while training a visualisation model and this is called GGTM with feature saliency (GGTM-FS). We demonstrate effectiveness of these proposed models both for synthetic and real datasets. We evaluate visualisation quality using quality metrics such as distance distortion measure and rank based measures: trustworthiness, continuity, mean relative rank errors with respect to data space and latent space. In cases where the labels are known we also use quality metrics of KL divergence and nearest neighbour classifications error in order to determine the separation between classes. We demonstrate the efficacy of these proposed models both for synthetic and real biological datasets with a main focus on the MHC class-I dataset.

THE IMPORTANCE OF MHC-INDEPENDENT NATURAL KILLER CELL RECEPTORS IN THE IMMUNOREGULATION OF MURINE PREGNANCY

Numerous leukocyte populations are essential for pregnancy success. Uterine natural killer (uNK) cells are chief amongst these leukocytes and represent a unique lineage with limited cytotoxicity but abundant angiokine production. They possess a distinct phenotype of activating and inhibitory receptors that recognize major histocompatibility complex (MHC) molecules, such as the killer immunoglobulin like receptors (KIRs; mouse Ly49), and MHC-independent activating receptors, including the aryl hydrocarbon receptor (AHR) and natural cytotoxicity receptor 1 (NCR1). While the roles of MHC-dependent receptors are widely addressed in pregnancy, MHC-independent receptors are relatively unstudied. This thesis investigated the roles of MHC-independent receptors in promotion of mouse pregnancy and characterized early leukocyte interactions in the presence and absence of NCR1. It was hypothesized that loss of MHC-independent receptors impairs uNK cell development resulting in aberrations in leukocyte function and decidual vasculature. Implantation sites from Ahr-/- and Ncr1Gfp/Gfp mice were assessed using whole mount in situ immunohistochemistry (WM-IHC) and histochemical techniques. Leukocyte interactions identified during preliminary WM-IHC studies were confirmed as immune synapses. The novel identification of immune synapses in early mouse pregnancy compelled further examination of leukocyte conjugates in wildtype C57BL/6 and Ncr1Gfp/Gfp mice. In Ahr-/- and Ncr1Gfp/Gfp mice, receptor loss resulted in reduced uNK cell diameters, impaired decidual vasculature, and failures in spiral artery remodeling. Ahr-/- mice had severe fertility deficits whereas Ncr1Gfp/Gfp mice had increased fetal resorption indicating differing receptor requirements in pregnancy success. NCR1 loss primarily affected uNK cell maturation and function as identified by alterations in granule ultrastructure, lytic protein expression, and angiokine production. Leukocyte conjugates were frequent in early C57BL/6 decidua basalis and included uNK cells conjugating first with antigen presenting cells and then with T cells. Overall conjugate formation was reduced in the absence of NCR1, but specific uNK cell conjugations were unaffected by receptor loss. While KIR-MHC interactions are associated with numerous pregnancy complications in humans, the role of other uNK cell receptors are not well characterized. These results illustrate the importance of MHC-independent receptors in uNK cell activation during early pregnancy in mice and encourage further studies of pregnancy complications that may occur independently of maternal KIR-MHC contributions.