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.     

Conformational restraints and flexibility of 14-meric peptides in complex with HLA-B*3501

Human HLA-B*3501 binds an antigenic peptide of 14-aa length derived from an alternative reading frame of M-CSF with high affinity. Due to its extraordinary length, the exact HLA binding mode was unpredictable. The crystal structure of HLA-B*3501 at 1.5 A shows that the N and C termini of the peptide are embedded in the A and F pockets, respectively, similar to a peptide of normal length. The central part of the 14-meric peptide bulges flexibly out of the groove. Two variants of the alternative reading frame of M-CSF peptide substituted at P2 or P2 and P9 with Ala display weak or no T cell activation. Their structure differs mainly in flexibility and conformation from the agonistic peptide. Moreover, the variants induce subtle changes of MHC alpha-helical regions implicated as critical for TCR contact. The TCR specifically recognizing this peptide/MHC complex exhibits CDR3 length within the normal range, suggesting major conformational adaptations of this receptor upon peptide/MHC binding. Thus, the potential antigenic repertoire recognizable by CTLs is larger than currently thought.     

Cooperativity of hydrophobic anchor interactions: evidence for epitope selection by MHC class II as a folding process

Peptide binding to MHC class II (MHCII) molecules is stabilized by hydrophobic anchoring and hydrogen bond formation. We view peptide binding as a process in which the peptide folds into the binding groove and to some extent the groove folds around the peptide. Our previous observation of cooperativity when analyzing binding properties of peptides modified at side chains with medium to high solvent accessibility is compatible with such a view. However, a large component of peptide binding is mediated by residues with strong hydrophobic interactions that bind to their respective pockets. If these reflect initial nucleation events they may be upstream of the folding process and not show cooperativity. To test whether the folding hypothesis extends to these anchor interactions, we measured dissociation and affinity to HLA-DR1 of an influenza hemagglutinin-derived peptide with multiple substitutions at major anchor residues. Our results show both negative and positive cooperative effects between hydrophobic pocket interactions. Cooperativity was also observed between hydrophobic pockets and positions with intermediate solvent accessibility, indicating that hydrophobic interactions participate in the overall folding process. These findings point out that predicting the binding potential of epitopes cannot assume additive and independent contributions of the interactions between major MHCII pockets and corresponding peptide side chains.     

Predicting peptide binding to MHC pockets via molecular modeling, implicit solvation, and global optimization

Development of a computational prediction method based on molecular modeling, global optimization, and implicit solvation has produced accurate structure and relative binding affinity predictions for peptide amino acids binding to five pockets of the MHC molecule HLA-DRB1*0101. Because peptide binding to MHC molecules is essential to many immune responses, development of such a method for understanding and predicting the forces that drive binding is crucial for pharmaceutical design and disease treatment. Underlying the development of this prediction method are two hypotheses. The first is that pockets formed by the peptide binding groove of MHC molecules are independent, separating the prediction of peptide amino acids that bind within individual pockets from those that bind between pockets. The second hypothesis is that the native state of a system composed of an amino acid bound to a protein pocket corresponds to the system's lowest free energy. The prediction method developed from these hypotheses uses atomistic-level modeling, deterministic global optimization, and three methods of implicit solvation: solvent-accessible area, solvent-accessible volume, and Poisson-Boltzmann electrostatics. The method predicts relative binding affinities of peptide amino acids for pockets of HLA-DRB1*0101 by determining computationally an amino acid's global minimum energy conformation. Prediction results from the method are in agreement with X-ray crystallography data and experimental binding assays.     

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.     

Crystallographic structure of the human leukocyte antigen DRA, DRB3*0101: models of a directional alloimmune response and autoimmunity

We describe structural studies of the human leukocyte antigen DR52a, HLA-DRA/DRB3*0101, in complex with an N-terminal human platelet integrin alphaII(B)betaIII glycoprotein peptide which contains a Leu/Pro dimorphism. The 33:Leu dimorphism is the epitope for the T cell directed response in neonatal alloimmune thrombocytopenia and post-transfusion purpura in individuals with the alphaII(B)betaIII 33:Pro allele, and defines the unidirectional alloimmune response. This condition is always associated with DR52a. The crystallographic structure has been refined to 2.25 A. There are two alphabeta heterodimers to the asymmetric unit in space group P4(1)2(1)2. The molecule is characterized by two prominent hydrophobic pockets at either end of the peptide binding cleft and a deep, narrower and highly charged P4 opening underneath the beta 1 chain. Further, the peptide in the second molecule displays a sharp upward turn after pocket P9. The structure reveals the role of pockets and the distinctive basic P4 pocket, shared by DR52a and DR3, in selecting their respective binding peptide repertoire. We observe an interesting switch in a residue from the canonically assigned pocket 6 seen in prior class II structures to pocket 4. This occludes the P6 pocket helping to explain the distinctive "1-4-9" peptide binding motif. A beta57 Asp-->Val substitution abrogates the salt-bridge to alpha76 Arg and along with a hydrophobic beta37 is important in shaping the P9 pocket. DRB3*0101 and DRB1*0301 belong to an ancestral haplotype and are associated with many autoimmune diseases linked to antigen presentation, but whereas DR3 is susceptible to type 1 diabetes DR52a is not. This dichotomy is explored for clues to the disease.     

Rapid and accurate structure-based therapeutic peptide design using GPU accelerated thermodynamic integration

Peptide-based therapeutics are an alternative to small molecule drugs as they offer superior specificity, lower toxicity, and easy synthesis. Here we present an approach that leverages the dramatic performance increase afforded by the recent arrival of GPU accelerated thermodynamic integration (TI). GPU TI facilitates very fast, highly accurate binding affinity optimization of peptides against therapeutic targets. We benchmarked TI predictions using published peptide binding optimization studies. Prediction of mutations involving charged side-chains was found to be less accurate than for non-charged, and use of a more complex 3-step TI protocol was found to boost accuracy in these cases. Using the 3-step protocol for non-charged side-chains either had no effect or was detrimental. We use the benchmarked pipeline to optimize a peptide binding to our recently discovered cancer target: EME1. TI calculations predict beneficial mutations using both canonical and non-canonical amino acids. We validate these predictions using fluorescence polarization and confirm that binding affinity is increased. We further demonstrate that this increase translates to a significant reduction in pancreatic cancer cell viability.     

Crystal structure of a myristoylated CAP-23/NAP-22 N-terminal domain complexed with Ca2+/calmodulin

A variety of viral and signal transduction proteins are known to be myristoylated. Although the role of myristoylation in protein-lipid interaction is well established, the involvement of myristoylation in protein-protein interactions is less well understood. CAP-23/NAP-22 is a brain-specific protein kinase C substrate protein that is involved in axon regeneration. Although the protein lacks any canonical calmodulin (CaM)-binding domain, it binds CaM with high affinity. The binding of CAP-23/NAP-22 to CaM is myristoylation dependent and the N-terminal myristoyl group is directly involved in the protein-protein interaction. Here we show the crystal structure of Ca2+-CaM bound to a myristoylated peptide corresponding to the N-terminal domain of CAP-23/NAP-22. The myristoyl moiety of the peptide goes through a hydrophobic tunnel created by the hydrophobic pockets in the N- and C-terminal domains of CaM. In addition to the myristoyl group, several amino-acid residues in the peptide are important for CaM binding. This is a novel mode of binding and is very different from the mechanism of binding in other CaM-target complexes.     

How a T cell receptor-like antibody recognizes major histocompatibility complex-bound peptide

We determined the crystal structures of the T cell receptor (TCR)-like antibody 25-D1.16 Fab fragment bound to a complex of SIINFEKL peptide from ovalbumin and the H-2K(b) molecule. Remarkably, this antibody directly "reads" the structure of the major histocompatibility complex (MHC)-bound peptide, employing the canonical diagonal binding mode utilized by most TCRs. This is in marked contrast with another TCR-like antibody, Hyb3, bound to melanoma peptide MAGE-A1 in association with HLA-A1 MHC class I. Hyb3 assumes a non-canonical orientation over its cognate peptide-MHC and appears to recognize a conformational epitope in which the MHC contribution is dominant. We conclude that TCR-like antibodies can recognize MHC-bound peptide via two different mechanisms: one is similar to that exploited by the preponderance of TCRs and the other requires a non-canonical antibody orientation over the peptide-MHC complex.     

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.     

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).     

Identification of high affinity polo-like kinase 1 (Plk1) polo-box domain binding peptides using oxime-based diversification

In an effort to develop improved binding antagonists of the polo-like kinase 1 (Plk1) polo-box domain (PBD), we optimized interactions of the known high affinity 5-mer peptide PLHSpT using oxime-based post solid-phase peptide diversification of the N-terminal Pro residue. This allowed us to achieve up to two orders of magnitude potency enhancement. An X-ray crystal structure of the highest affinity analogue in complex with Plk1 PBD revealed new binding interactions in a hydrophobic channel that had been occluded in X-ray structures of the unliganded protein. This study represents an important example where amino acid modification by post solid-phase oxime ligation can facilitate the development of protein-protein interaction inhibitors by identifying new binding pockets that would not otherwise be accessible to coded amino acid residues.     

High resolution structure of streptavidin in complex with a novel high affinity peptide tag mimicking the biotin binding motif

A novel peptide was designed which possesses nanomolar affinity of less than 20 nM for streptavidin. Therefore it was termed Nano-tag and has been used as an affinity tag for recombinant proteins. The minimized version of the wild type Nano-tag is a seven-amino acid peptide with the sequence fMDVEAWL. The three-dimensional structure of wild type streptavidin in complex with the minimized Nano-tag was analyzed at atomic resolution of 1.15 A and the details of the binding motif were investigated. The peptide recognizes the same pocket of streptavidin where the natural ligand biotin is bound, but the peptide requires significantly more space than biotin. Therefore the binding loop adopts an "open" conformation in order to release additional space for the peptide. The conformation of the bound Nano-tag corresponds to a 3(10) helix. However, the analysis of the intermolecular interactions of the Nano-tag with residues of the binding pocket of streptavidin reveals astonishing similarities to the biotin binding motif. In principle the three-dimensional conformation of the Nano-tag mimics the binding mode of biotin. Our results explain why the use of the Nano-tag in fusion with recombinant proteins is restricted to their N-terminus and we describe the special significance of the fMet residue for the high affinity binding mode.     

Constrained peptides mimic a viral suppressor of RNA silencing

The design of high-affinity, RNA-binding ligands has proven very challenging. This is due to the unique structural properties of RNA, often characterized by polar surfaces and high flexibility. In addition, the frequent lack of well-defined binding pockets complicates the development of small molecule binders. This has triggered the search for alternative scaffolds of intermediate size. Among these, peptide-derived molecules represent appealing entities as they can mimic structural features also present in RNA-binding proteins. However, the application of peptidic RNA-targeting ligands is hampered by a lack of design principles and their inherently low bio-stability. Here, the structure-based design of constrained α-helical peptides derived from the viral suppressor of RNA silencing, TAV2b, is described. We observe that the introduction of two inter-side chain crosslinks provides peptides with increased α-helicity and protease stability. One of these modified peptides (B3) shows high affinity for double-stranded RNA structures including a palindromic siRNA as well as microRNA-21 and its precursor pre-miR-21. Notably, B3 binding to pre-miR-21 inhibits Dicer processing in a biochemical assay. As a further characteristic this peptide also exhibits cellular entry. Our findings show that constrained peptides can efficiently mimic RNA-binding proteins rendering them potentially useful for the design of bioactive RNA-targeting ligands.     

Identification of model peptides as affinity ligands for the purification of humanized monoclonal antibodies by means of phage display.

A proof-of-principle study was initiated to determine whether phage-display technology could be used to identify peptides as leads in the customization of ligands for affinity chromatography and to identify a peptide or peptidomimetic for use as a Protein A alternative in the affinity purification of monoclonal antibodies. The constant region of humanized anti-Tac (HAT), prepared by pepsin digestion and receptor-affinity chromatography, was used as the target for phage display in this study. As such, 20 phage-derived peptide sequences were identified from four rounds of biopanning with two linear phage-display libraries (7-mer, containing 100 copies of 2 x 10(9) sequences and 12-mer, containing 70 copies of 1.4 x 10(9) sequences). Five peptides were synthesized for use as affinity ligands, based on sequence homology to Protein A, sequence redundancy, and amino acid motifs. The best HAT binding immobilized peptide was EPIHRSTLTALL. The best-fit analysis of this peptide sequence with Protein A yielded an alignment well within the Fc binding domain of Protein A. These results suggest that phage display can serve as a tool in the identification of peptides as model ligands for affinity chromatography.     

Evidence for a new subclass of calcitonin/ calcitonin gene-related peptide binding site in rat brain

Binding sites for calcitonin and calcitonin gene-related peptide are widely distributed in the central nervous system. In this study, binding of [¹²⁵I]-alpha-rat calcitonin gene-related peptide and [¹²⁵I]-salmon calcitonin in adjacent sections of rat brain revealed clearly distinct patterns of binding in most regions although in some restricted areas such as parts of the ventral striatum, including the nucleus accumbens, there was some overlap in the patterns of binding. In the primary olfactory cortex, which bound only calcitonin gene-related peptide, salmon calcitonin was very weak in inhibiting the binding of calcitonin gene-related peptide. In the nucleus accumbens, high affinity binding of calcitonin and calcitonin gene-related peptide at their homologous receptors was observed, with affinity constants for calcitonin and calcitonin gene-related peptide of 1.4 × 10⁹ M⁻¹ and 1.2 × 10⁹ M⁻¹ respectively. Cross competition studies in this nucleus demonstrated that salmon calcitonin was able to compete for [¹²⁵I]-rat calcitonin gene-related peptide labelled sites with high affinity, with an affinity constant of 0.8 × 10⁹ M⁻¹. However, rat calcitonin gene-related peptide was less potent in inhibiting the binding of [¹²⁵I]-salmon calcitonin labelled sites with only 28% inhibition at 10⁻⁶M. Further characterization of the calcitonin sensitive calcitonin gene-related peptide labelled sites demonstrated that a range of calcitonin analogs inhibited the binding of [¹²⁵I]-rat calcitonin gene-related peptide with the same order of potency as the analogs competed for [¹²⁵I]-salmon calcitonin labelled sites. Digital substraction mapping revealed calcitonin-sensitive calcitonin gene-related peptide binding sites over parts of the ventral striatum, including mid-caudal nucleus accumbens and fundus striati; over the lateral border of the lateral bed nucleus of the stria terminalis; part of the central amygdaloid nucleus; the organum vasculosum of the lamina terminalis and area postrema and over the wings of the dorsal raphe.These results demonstrate the existence of a new subtype of calcitonin/calcitonin gene-related peptide binding site, which has high affinity for the two otherwise biochemically distinct peptides.     

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     

Importance of N- and C-terminal regions of gastrin-Gly for preferential binding to high and low affinity gastrin-Gly receptors

G17-Gly has been shown to stimulate the growth of DLD-1 human colon cancer cells in a biphasic manner via high and low affinity receptors. In the current study, the existence of heterogeneous receptor populations for G17-Gly on the HT-29 human colon cancer cell line was investigated. The effect of either N- or C-terminal peptide truncation on receptor binding and cell growth stimulation was also explored. [Leu15]G17-Gly bound to both high (nM) and low (microM) affinity sites on HT-29 cells. The peptide stimulated cell growth in a dose-dependent and biphasic manner with maximal stimulation at 10(-9) M peptide concentration, suggesting that, as in the case of DLD-1 cells, it is the high affinity receptor which is responsible for the growth-promoting effects. In contrast, G17(1-12) stimulated the growth of HT-29 cells in a sigmoidal fashion with an EC50 of 4.6x10(-9) M. Sequential N-terminal truncation of [Leu15]G17-Gly results in decreased binding to the high affinity G17-Gly receptor on DLD-1 cells. [Leu15]G17(11-17)Gly bound to the low affinity G17-Gly receptor with an affinity similar to that of the full sequence peptide but was unable to displace the radioligand from high affinity sites. G17(1-6)-NH2 was unable to displace [3H]G17-Gly from either site. These results suggest that the important residues for binding to the low affinity receptor are in the C-terminal region of the peptide while those required for interaction with the high affinity receptor lie further towards the N-terminus.     

Isolation and biochemical characterization of fusicoccin-insensitive variant lines of Arabidopsis thaliana (L.) Heynh

Arabidopsis thaliana lines have been isolated that are insensitive to the fungal toxin fusicoccin (FC). Initial screening was done by selecting for plants that either grew well on high concentrations of FC or did not respond to FC by increases in H⁺-extrusion. All selected plants were tested, in several additional rounds of screening, for binding to microsomal proteins of a ³H-labeled radioligand of fusicoccin. A novel assay allowing for the direct selection of individual plants exhibiting reduced binding of FC was developed and used as screening procedure. Independent variant lines (43) with stably expressed, reduced binding of FC were isolated and subjected to a detailed characterization of their binding sites. The lines could be subdivided into several distinct classes with respect to these characteristics. In class-I lines, the data indicate a partial conversion of high-affinity binding sites to a low-affinity state. In class-II lines, the affinity of the binding site to FC is strongly reduced while the number of sites, as well as several other biochemical parameters, is completely unchanged, suggesting a specific alteration in the properties of the fusicoccin-binding protein. In class-III lines, the ligand-binding protein complex, while retaining its high affinity, is destabilized at supraoptimal concentrations of FC (such as those used for screening). In wild-type plants, only the high-affinity binding site was detected. Combined, these data prove that the high-affinity sites represent the plant's FC receptor.Abbreviations A_o binding site concentration - FC fusicoccin - FCBP fusicoccin-binding protein - FCol 9-nor-8-hydroxyfusicoccin - K_D dissociation constant of the FCBP-radioligand complex We are grateful to Iris Sandorf and Gudrun Henrichs for excellent technical assistance. This work was supported by the Deutsche Forschungsgemeinschaft, Bonn, Germany and by Fonds der Chemischen Industrie (literature provision).