Cooperativity during the formation of peptide/MHC class II complexes

To generate an effective immune response, class II major histocompatibility complex molecules (MHCII) must present a diverse array of peptide ligands for recognition by T lymphocytes. Peptide/MHCII complexes are stabilized by hydrophobic anchoring of peptide side chains to pockets in the MHCII protein and the formation of hydrogen bonds to the peptide backbone. Many current models of peptide/MHCII association assume an additive and independent contribution of the interactions between major MHCII pockets and corresponding side chains in the peptide. However, significant conformational rearrangements occur in both the peptide and MHCII during binding. Therefore, we hypothesize that peptide binding to MHCII could be viewed as a folding process in which both molecules cooperate to produce the final conformation. To directly test this hypothesis, we adapt a serial mutagenesis strategy to study cooperativity in the interaction of the human MHCII HLA-DR1 and a peptide derived from influenza hemagglutinin. Substitutions in either the peptide or HLA-DR1 that are predicted to interfere with hydrogen bond formation show cooperative effects on complex stability and affinity. Substitution of a peptide side chain that provides a hydrophobic contact also contributes to the cooperative effect, suggesting a role for all energetic sources in the folding process. We propose that cooperativity throughout the peptide-binding groove reflects the folding of segments of the MHCII molecule into helices around the peptide with a concomitant folding of the peptide into a polyproline helix. The implications of cooperativity for peptide/MHCII structure and epitope selection are discussed.     

Cutting edge: TCR contacts as anchors: effects on affinity and HLA-DM stability

Peptides presented via the class II MHC (MHCII) pathway are selected based on affinity for MHCII and stability in the presence of HLA-DM. Currently, epitope selection is thought to be controlled by the ability of peptide to sequester "anchor" residues into pockets in the MHCII. Residues exhibiting higher levels of solvent accessibility have been shown to contact TCR, but their roles in affinity and complex stability have not been directly studied. Using the HLA-DR1-binding influenza peptide, hemagglutinin (306-318), as a model, we show that side chain substitutions at these positions influence affinity and HLA-DM stability. Multiple substitutions reduce affinity to a greater extent than the loss of the major P1 anchor residue. We propose that these effects may be mediated through the H-bond network. These results demonstrate the importance of solvent-exposed residues in epitope selection and blur the distinctions between anchor and TCR contact residues.     

C-terminal anchoring of a peptide to class II MHC via the P10 residue is compatible with a peptide bulge.

The binding of antigenic peptide to class II MHC is mediated by hydrogen bonds between the MHC and the peptide, by salt bridges, and by hydrophobic interactions. The latter are confined to a number of deeper pockets within the peptide binding groove, and peptide side chains that interact with these pockets are referred to as anchor residues. T cell recognition involves solvent-accessible peptide residues along with minor changes in MHC helical pitch induced by the anchor residues. In class I MHC there is an added level of epitope complexity that results from binding of longer peptides that bulge out into the solvent-accessible, T cell contact area. Unlike class I MHC, class II MHC does not bind peptides of discrete length, and the possibility of peptide bulging has not been clearly addressed. A peptide derived from position 24-37 of integrin beta(3) can either bind or not bind to the class II MHC molecule HLA DRB3*0101 based on a polymorphism at the P9 anchor. We show that the loss of binding can be compensated by changes at the P10 position. We propose that this could be an example of a class II peptide bulge. Although not as efficient as P9 anchoring, the use of P10 as an anchor adds another possible mechanism by which T cell epitopes can be generated in the class II presentation system.     

Promiscuous antigen presentation by the nonclassical MHC Ib Qa-2 is enabled by a shallow, hydrophobic groove and self-stabilized peptide conformation.

BACKGROUND: Qa-2 is a nonclassical MHC Ib antigen, which has been implicated in both innate and adaptive immune responses, as well as embryonic development. Qa-2 has an unusual peptide binding specificity in that it requires two dominant C-terminal anchor residues and is capable of associating with a substantially more diverse array of peptide sequences than other nonclassical MHC. RESULTS: We have determined the crystal structure, to 2.3 A, of the Q9 gene of murine Qa-2 complexed with a self-peptide derived from the L19 ribosomal protein, which is abundant in the pool of peptides eluted from the Q9 groove. The 9 amino acid peptide is bound high in a shallow, hydrophobic binding groove of Q9, which is missing a C pocket. The peptide makes few specific contacts and exhibits extremely poor shape complementarity to the MHC groove, which facilitates the presentation of a degenerate array of sequences. The L19 peptide is in a centrally bulged conformation that is stabilized by intramolecular interactions from the invariant P7 histidine anchor residue and by a hydrophobic core of preferred secondary anchor residues that have minimal interaction with the MHC. CONCLUSIONS: Unexpectedly, the preferred secondary peptide residues that exhibit tenuous contact with Q9 contribute significantly to the overall stability of the peptide-MHC complex. The structure of this complex implies a "conformational" selection by Q9 for peptide residues that optimally stabilize the large bulge rather than making intimate contact with the MHC pockets.     

Stabilizing roles of residual structure in the empty heme binding pockets and unfolded states of microsomal and mitochondrial apocytochrome b5

The microsomal (Mc) and mitochondrial (OM) isoforms of mammalian cytochrome b5 are the products of different genes, which likely arose via duplication of a primordial gene and subsequent functional divergence. Despite sharing essentially identical folds, heme-polypeptide interactions are stronger in OM b5s than in Mc b5s due to the presence of two conserved patches of hydrophobic amino acid side chains in the OM heme binding pockets. This is of fundamental interest in terms of understanding heme protein structure-function relationships, because stronger heme-polypeptide interactions in OM b5s in comparison to Mc b5s may represent a key source of their more negative reduction potentials. Herein we provide evidence that interactions amongst the amino acid side chains contributing to the hydrophobic patches in rat OM (rOM) b5 persist when heme is removed, rendering the empty heme binding pocket of rOM apo-b5 more compact and less conformationally dynamic than that in bovine Mc (bMc) apo-b5. This may contribute to the stronger heme binding by OM apo-b5 by reducing the entropic penalty associated with polypeptide folding. We also show that when bMc apo-b5 unfolds it adopts a structure that is more compact and contains greater nonrandom secondary structure content than unfolded rOM apo-b5. We propose that a more robust beta-sheet in Mc apo-b5s compensates for the absence of the hydrophobic packing interactions that stabilize the heme binding pocket in OM apo-b5s.     

Molecular Modeling Studies on Novel Open-chain and Cyclic Thia Compounds and their Ag(I) and Hg(II) Complexes

Molecular modeling techniques were used to generate structures of several HLA-DQ proteins associated with insulin-dependent diabetes mellitus (IDDM). A peptide fragment from glutamic acid decarboxylase (GAD), a known IDDM autoantigen, binds to certain HLA-DQ molecules positively associated with IDDM. Modeling studies were used to explore possible binding interactions between this GAD peptide and several HLA-DQ molecules. Based on the characterization of anchor pockets in the HLA-DQ binding groove and of peptide side chains, a novel binding mode was proposed. This binding mode predicts the GAD peptide is positioned in the binding groove in the direction opposite the orientation observed for class I proteins and the class II DR1, DR3, and I-E^k proteins. Peptide docking exercises were performed to construct models of the HLA-DQ/peptide complexes, and the resulting models have been used to design peptide binding experiments to test this reverse-orientation binding mode. A variety of experimental results are consistent with the proposed model and suggest that some peptide ligands of class II molecules may bind in a reversed orientation within the binding groove.Supplementary material to this paper is available in electronic form at http://dx.doi.org/10.1007/s0089460020205     

Flanking p10 contribution and sequence bias in matrix based epitope prediction: revisiting the assumption of independent binding pockets

Background  

Eluted natural peptides from major histocompatibility molecules show patterns of conserved residues. Crystallographic structures show that the bound peptide in class II major histocompatibility complex adopts a near uniform polyproline II-like conformation. This way allele-specific favoured residues are able to anchor into pockets in the binding groove leaving other peptide side chains exposed for recognition by T cells. The anchor residues form a motif. This sequence pattern can be used to screen large sequences for potential epitopes. Quantitative matrices extend the motif idea to include the contribution of non-anchor peptide residues. This report examines two new matrices that extend the binding register to incorporate the polymorphic p10 pocket of human leukocyte antigen DR1. Their performance is quantified against experimental binding measurements and against the canonical nine-residue register matrix.     

Trifluoroethanol may form a solvent matrix for assisted hydrophobic interactions between peptide side chains

Several models for interactions between trifluoroethanol (TFE) and peptides and proteins have recently been proposed, but none have been able to rationalize the puzzling observations that on the one hand TFE can stabilize some hydrophobic interactions in secondary structures, but on the other can also melt the hydrophobic cores of globular proteins. The former is illustrated in this paper by the effect of TFE on a short elastin peptide, GVG(VPGVG)(3), which forms type II beta-turns stabilized by hydrophobic interactions between two intra-turn valine side chains. This folding, driven by increasing the entropy of bulk water, is stimulated in TFE-water mixtures and/or by raising the temperature. To explain these apparently contradictory observations, we propose a model in which TFE clusters locally assist the folding of secondary structures by first breaking down interfacial water molecules on the peptide and then providing a solvent matrix for further side chain--side chain interactions. This model also provides an explanation for TFE-induced transitions between secondary structures, in which the TFE clusters may redirect non-local to local interactions.     

The Carboxy Terminus of the Ligand Peptide Determines the Stability of the MHC Class I Molecule H-2Kb: A Combined Molecular Dynamics and Experimental Study

Major histocompatibility complex (MHC) class I molecules (proteins) bind peptides of eight to ten amino acids to present them at the cell surface to cytotoxic T cells. The class I binding groove binds the peptide via hydrogen bonds with the peptide termini and via diverse interactions with the anchor residue side chains of the peptide. To elucidate which of these interactions is most important for the thermodynamic and kinetic stability of the peptide-bound state, we have combined molecular dynamics simulations and experimental approaches in an investigation of the conformational dynamics and binding parameters of a murine class I molecule (H-2K^b) with optimal and truncated natural peptide epitopes. We show that the F pocket region dominates the conformational and thermodynamic properties of the binding groove, and that therefore the binding of the C terminus of the peptide to the F pocket region plays a crucial role in bringing about the peptide-bound state of MHC class I.     

Protein surface mapping by chemical oxidation: structural analysis by mass spectrometry

The solvent-accessible surface area of proteins is important in biological function for many reasons, including protein-protein interactions, protein folding, and catalytic sites. Here we present a chemical technique to oxidize amino acid side chains in a model protein, apomyoglobin, and subsequent elucidation of the effect of solvent accessibility on the sites of oxidation. Under conditions of low protein oxidation (zero to three oxygen atoms added per apomyoglobin molecule), we have positively identified five oxidation sites by liquid chromatography-tandem mass spectrometry and high-resolution Fourier transform mass spectrometry. Our results indicate that all oxidized amino acids, with the exception of methionine, have highly solvent-accessible side chains, but the rate of oxidation may not be dictated solely by solvent accessibility and amino acid identity.     

Energetics and cooperativity of the hydrogen bonding and anchor interactions that bind peptides to MHC class II protein

The complexity of the interaction between major histocompatibility complex class II (MHC II) proteins and peptide ligands has been revealed through structural studies and crystallographic characterization. Peptides bind through side-chain "anchor" interactions with MHC II pockets and an extensive array of genetically conserved hydrogen bonds to the peptide backbone. Here we quantitatively investigate the kinetic hierarchy of these interactions. We present results detailing the impact of single side-chain mutations of peptide anchor residues on dissociation rates, utilizing two I-A(d)-restricted peptides, one of which has a known crystal structure, and 24 natural and non-natural amino acid mutant variants of these peptides. We find that the N-terminal P1, P4 and P6 anchor-pocket interactions can make significant contributions to binding stability. We also investigate the interactions of these peptides with four I-A(d) MHC II proteins, each mutated to disrupt conserved hydrogen bonds to the peptide backbone. These complexes exhibit kinetic behavior suggesting that binding energy is disproportionately invested near the peptide N terminus for backbone hydrogen bonds. We then evaluate the effects of simultaneously modifying both anchor and hydrogen bonding interactions. A quantitative analysis of 71 double mutant cycles reveals that there is little apparent cooperativity between anchor residue interactions and hydrogen bonds, even when they are directly adjacent (<5A).     

Using metadynamics to understand the mechanism of calmodulin/target recognition at atomic detail

The ability of calcium-bound calmodulin (CaM) to recognize most of its target peptides is caused by its binding to two hydrophobic residues ('anchors'). In most of the CaM complexes, the anchors pack against the hydrophobic pockets of the CaM domains and are surrounded by fully conserved Met side chains. Here, by using metadynamics simulations, we investigate quantitatively the energetics of the final step of this process using the M13 peptide, which has a high affinity and spans the sequence of the skeletal myosin light chain kinase, an important natural CaM target. We established the accuracy of our calculations by a comparison between calculated and NMR-derived structural and dynamical properties. Our calculations provide novel insights into the mechanism of protein/peptide recognition: we show that the process is associated with a free energy gain similar to that experimentally measured for the CaM complex with the homologous smooth muscle MLCK peptide (Ehrhardt et al., 1995, Biochemistry 34, 2731). We suggest that binding is dominated by the entropic effect, in agreement with previous proposals. Furthermore, we explain the role of conserved methionines by showing that the large flexibility of these side chains is a key feature of the binding mechanism. Finally, we provide a rationale for the experimental observation that in all CaM complexes the C-terminal domain seems to be hierarchically more important in establishing the interaction.     

Determination of protein folding kinetic types using sequence and predicted secondary structure and solvent accessibility

Proteins fold through a two-state (TS), with no visible intermediates, or a multi-state (MS), via at least one intermediate, process. We analyze sequence-derived factors that determine folding types by introducing a novel sequence-based folding type predictor called FOKIT. This method implements a logistic regression model with six input features which hybridize information concerning amino acid composition and predicted secondary structure and solvent accessibility. FOKIT provides predictions with average Matthews correlation coefficient (MCC) between 0.58 and 0.91 measured using out-of-sample tests on four benchmark datasets. These results are shown to be competitive or better than results of four modern predictors. We also show that FOKIT outperforms these methods when predicting chains that share low similarity with the chains used to build the model, which is an important advantage given the limited number of annotated chains. We demonstrate that inclusion of solvent accessibility helps in discrimination of the folding kinetic types and that three of the features constitute statistically significant markers that differentiate TS and MS folders. We found that the increased content of exposed Trp and buried Leu are indicative of the MS folding, which implies that the exposure/burial of certain hydrophobic residues may play important role in the formation of the folding intermediates. Our conclusions are supported by two case studies.     

Smoothing of the GB1 Hairpin Folding Landscape by Interfacial Confinement

We study the effects of confinement between planar walls on the folding thermodynamics of a β-hairpin, using large-scale replica-exchange molecular-dynamics simulations with an all-atom model and explicit solvent. We find that the folding free-energy landscape of this peptide observed in bulk is significantly modified when the peptide is confined between the walls. Most notably, the propensity of the peptide to form a misfolded state observed in the bulk solution becomes negligible under confinement. The absence of the misfolded state under confinement can be explained by an increased tendency of hydrophobic aromatic side chains to stay near the walls, because the misfolded state is characterized by a nonnative arrangement of aromatic side chains. These results from a simple confinement model may provide clues about the role of chaperonin confinement in smoothing folding landscapes by avoiding trapped intermediates.     

Structural resiliency of an EGF-like subdomain bound to its target protein, thrombin.

The thrombin-bound structures of native peptide fragments from the fifth EGF-like domain of thrombomodulin were determined by use of NMR and transferred NOE spectroscopy. The bound peptides assume an EGF-like structure of an antiparallel beta-sheet, a novel structural motif observed for a bound peptide in protein-peptide complexes. There is a remarkable structural resiliency of this structure motif manifested in its ability to accommodate a different number of residues within the disulfide loop. Docking experiments revealed that the key contacts with thrombin are hydrophobic interactions between the side chains of residues Ile 414 and Ile 424 of thrombomodulin and a hydrophobic pocket on the thrombin surface. Residues Leu 415, Phe 419, and Ile 420, which would have been buried in intact EGF-like domains, are unfavorably exposed in the complex of thrombin with the EGF-like thrombomodulin fragment, thus providing a rationale for the enhancement of binding affinity upon the deletion of Ile 420. The unique beta-sheet structures of the bound peptides are specified by the presence of disulfide bridges in the peptides because a corresponding linear thrombomodulin fragment folds into a sheet structure with a different backbone topology. The different bound conformations for the linear and the cyclized peptides indicate that side-chain interactions within a specific environment may dictate the folding of bound peptides in protein-peptide complexes.     

A novel target recognition revealed by calmodulin in complex with Ca2+-calmodulin-dependent kinase kinase.

The structure of calcium-bound calmodulin (Ca2+/CaM) complexed with a 26-residue peptide, corresponding to the CaM-binding domain of rat Ca2+/CaM-dependent protein kinase kinase (CaMKK), has been determined by NMR spectroscopy. In this complex, the CaMKK peptide forms a fold comprising an alpha-helix and a hairpin-like loop whose C-terminus folds back on itself. The binding orientation of this CaMKK peptide by the two CaM domains is opposite to that observed in all other CaM-target complexes determined so far. The N- and C-terminal hydrophobic pockets of Ca2+/CaM anchor Trp 444 and Phe 459 of the CaMKK peptide, respectively. This 14-residue separation between two key hydrophobic groups is also unique among previously determined CaM complexes. The present structure represents a new and distinct class of Ca2+/CaM target recognition that may be shared by other Ca2+/CaM-stimulated proteins.     

Changes at the floor of the peptide-binding groove induce a strong preference for proline at position 3 of the bound peptide: molecular dynamics simulations of HLA-A*0217

We report on molecular dynamics simulations of major histocompatibility complex (MHC)-peptide complexes. Class I MHC molecules play an important role in cellular immunity by presenting antigenic peptides to cytotoxic T cells. Pockets in the peptide-binding groove of MHC molecules accommodate anchor side chains of the bound peptide. Amino acid substitutions in MHC affect differences in the peptide-anchor motifs. HLA-A*0217, human MHC class I molecule, differs from HLA-A*0201 only by three amino acid residues substitutions (positions 95, 97, and 99) at the floor of the peptide-binding groove. A*0217 showed a strong preference for Pro at position 3 (p3) and accepted Phe at p9 of its peptide ligands, but these preferences have not been found in other HLA-A2 ligands. To reveal the structural mechanism of these observations, the A*0217-peptide complexes were simulated by 1000 ps molecular dynamics at 300 K with explicit solvent molecules and compared with those of the A*0201-peptide complexes. We examined the distances between the anchor side chain of the bound peptide and the pocket, and the rms fluctuations of the bound peptides and the HLA molecules. On the basis of the results from our simulations, we propose that Pro at p3 serves as an optimum residue to lock the dominant anchor residue (p9) tightly into pocket F and to hold the peptide in the binding groove, rather than a secondary anchor residue fitting optimally the complementary pocket. We also found that Phe at p9 is used to occupy the space created by replacements of three amino acid residues at the floor within the groove. These findings would provide a novel understanding in the peptide-binding motifs of class I MHC molecules.     

Crystal structure of a non-canonical low-affinity peptide complexed with MHC class I: a new approach for vaccine design

Peptides bind with high affinity to MHC class I molecules by anchoring certain side-chains (anchors) into specificity pockets in the MHC peptide-binding groove. Peptides that do not contain these canonical anchor residues normally have low affinity, resulting in impaired pMHC stability and loss of immunogenicity. Here, we report the crystal structure at 1.6 A resolution of an immunogenic, low-affinity peptide from the tumor-associated antigen MUC1, bound to H-2Kb. Stable binding is still achieved despite small, non-canonical residues in the C and F anchor pockets. This structure reveals how low-affinity peptides can be utilized in the design of novel peptide-based tumor vaccines. The molecular interactions elucidated in this non-canonical low-affinity peptide MHC complex should help uncover additional immunogenic peptides from primary protein sequences and aid in the design of alternative approaches for T-cell vaccines.     

Molecular dynamics study of the human insulin B peptide SHLVEALYLVCGERGG complexed with HLA-DQ8 reveals important hydrogen bond interactions

The basis of proper recognition of pathogens and tumours is provided by adaptive immunity. This immunological reaction of the recognition function of T-cell receptors on T lymphocytes detects antigenic peptides bound to major histocompatibility complex (MHC) molecules. Structural insight into this process has few grown considerably in the last years. In some of the cases, antigens are self-protein fragments causing autoimmunity diseases. Type 1 diabetes is such a disease connected with the human leukocyte antigen-DQ8 molecule, a class II MHC glycoprotein. Its crystal structure, complexed with LVEALYLVCGERGG peptide (insulin B peptide), has been solved, and important information about the significance of P1, P4 and P9 binding pockets has been discovered. The complex structure also revealed an unusual large number of intermolecular hydrogen bonds between insulin B peptide and MHC molecule. To further investigate the dynamics of peptide/MHC interactions, we perform molecular dynamic simulations in explicit water. Analysis of the results provided useful information of the binding of the peptide antigen to MHC molecule, which is supported by numerous hydrogen bonds besides the electrostatic (P1 and P9 pockets) or hydrophobic interactions (P4). Results also allowed some implications to be drawn for the role of residues located outside of the binding groove.     

Adaptive recognition in RNA complexes with peptides and protein modules

Recently, progress has been made towards the structural characterization of the novel folds of RNA-bound arginine-rich peptides and the architecture of their peptide-binding RNA pockets in viral and phage systems. These studies are based on an approach whereby the peptide and RNA components are minimalist modular domains that undergo adaptive structural transitions upon complex formation. Such complexes are characterized by recognition alignments in which the tertiary fold of the RNA generates binding pockets with the potential to envelop minimal elements of protein secondary structure. Strikingly, the peptides fold as isolated alpha-helical or beta-hairpin folds within their RNA major-groove targets, without the necessity of additional appendages for anchorage within the binding pocket. The RNA peptide-binding pocket architectures are sculptured through precisely positioned mismatches, triples and looped-out bases, which accommodate amino acid sidechains through hydrophobic, hydrogen bonding and ionic intermolecular contacts. By contrast, protein modules associated with the HIV-1 nucleocapsid and MS2 phage coat target their RNA binding sites through the insertion of specificity-determining RNA base residues within conserved hydrophobic pockets and crevices on the protein surface, with the bases anchored through hydrogen bonding interactions. These alternative strategies of RNA recognition at the peptide and protein module level provide novel insights into the principles, patterns and diversity of the adaptive transitions associated with the recognition process.