Peptide length significantly influences in vitro affinity for MHC class II molecules

Background

Class II Major Histocompatibility Complex (MHC) molecules have an open-ended binding groove which can accommodate peptides of varying lengths. Several studies have demonstrated that peptide flanking residues (PFRs) which lie outside the core binding groove can influence peptide binding and T cell recognition. By using data from the AntiJen database we were able to characterise systematically the influence of PFRs on peptide affinity for MHC class II molecules.

Results

By analysing 1279 peptide elongation events covering 19 distinct HLA alleles it was observed that, in general, peptide elongation resulted in increased MHC class II molecule affinity. It was also possible to determine an optimal peptide length for MHC class II affinity of approximately 18–20 amino acids; elongation of peptides beyond this length resulted in a null or negative effect on affinity.

Conclusion

The observed relationship between peptide length and MHC class II affinity has significant implications for the design of vaccines and the study of the epitopic basis of immunological disease.     

Thermodynamic and kinetic analysis of a peptide-class I MHC interaction highlights the noncovalent nature and conformational dynamics of the class I heterotrimer

The class I major histocompatibility (MHC) molecule is a heterotrimer composed of a heavy chain, the small subunit beta(2)-microglobulin (beta(2)m), and a peptide. Fluorescence anisotropy has been used to assay the interaction of a labeled peptide with a recombinant, soluble form of the class I MHC HLA-A2. Consistent with earlier work, peptide binding is shown to be a two-step process limited by a conformational rearrangement in the heavy chain/beta(2)m heterodimer. However, we identify two pathways for peptide dissociation from the heterotrimer: (1) initial peptide dissociation leaving a heavy chain/beta(2)m heterodimer and (2) initial dissociation of beta(2)m, followed by peptide dissociation from the heavy chain. Eyring analyses of rate constants measured as a function of temperature permit for the first time a complete thermodynamic characterization of peptide binding. We find that in this case peptide binding is mostly entropically driven, likely reflecting the hydrophobic character of the peptide binding groove and the peptide anchor residues. Thermodynamic and kinetic analyses of peptide-MHC interactions as performed here may be of practical use in the engineering of peptides with desired binding properties and will aid in the interpretation of the effects of MHC and peptide substitutions on peptide binding and T cell reactivity. Finally, our data suggest a role for beta(2)m in dampening conformational dynamics in the heavy chain. Remaining conformational variability in the heavy chain once beta(2)m has bound may be a mechanism to promote promiscuity in peptide binding.     

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.     

Thermodynamic study of the binding of the tandem-SH2 domain of the Syk kinase to a dually phosphorylated ITAM peptide: evidence for two conformers

The cytosolic tyrosine kinase Syk is recruited to immune cell receptors via interactions of its tandem-SH2 domain with tyrosine-phosphorylated sequences called immune receptor tyrosine activation motifs (ITAMs). We have characterized the binding of the tandem-SH2 domain of Syk (Syk-tSH2) to a dually phosphorylated peptide derived from the ITAM of the T cell receptor CD3-epsilon subunit. The CD3-epsilon peptide binds with an affinity of 18-81 nM at 150 mM NaCl over the 4.5-30 degrees C temperature range that was studied. The enthalpy of binding, DeltaH degrees obs, shows an unusual nonlinear dependence on temperature, suggesting the possibility of a temperature-dependent conformational equilibrium coupled to binding. This hypothesis was tested and confirmed by examining the temperature dependence of (1) the on-rate constant for binding and (2) the fluorescence of Syk-tSH2 and its CD3-epsilon peptide complex. The DeltaH degrees obs, Kobs, fluorescence, and kinetic data are all well described by a model incorporating the hypothesized conformational equilibrium. Circular dichroism spectra at various temperatures indicate that the conformational change is not due to a partial unfolding of the protein. We suggest that the conformational equilibrium enables Syk-tSH2 to exhibit flexibility in its binding modality, which may be necessitated by Syk's involvement in a wide variety of signal tranduction pathways.     

Molecular modeling of the thyroid hormone interactions with alpha v beta 3 integrin

A cell surface receptor for thyroid hormone has recently been identified on the extracellular domain of integrin alphavbeta3. In a variety of human and animal cell lines this hormone receptor mediates activation by thyroid hormone of the cellular mitogen-activated protein kinase (MAPK) signal transduction cascade. An arginine-glycine-aspartate (RGD) recognition site on the heterodimeric integrin is essential to the binding of a variety of extracellular matrix proteins. Recent competition data reveal that RGD peptides block hormone-binding by the integrin and consequent MAPK activation, suggesting that the hormone interaction site is located at or near the RGD recognition site on integrin alphavbeta3. A deaminated thyroid hormone (l-thyroxine, T4) analogue, tetraiodothyroacetic acid (tetrac, T4ac), inhibits binding of T4 and 3,5,3'-triiodo-l-thyronine (T3) to alphavbeta3, but does not activate MAPK. Structural data show that the RGD cyclic peptide binds at the interface of the propeller of the alphav and the B domains on the integrin head [Xiong JP, Stehle T, Zhang R, Joachimiack A, Frech M, Goodman SL, et al. Crystal structure of the extracellular segment of integrin alphavbeta3 in complexing with an Arg-Gly-Asp ligand. Science 2002;296:151-5]. To model potential interactions of thyroid hormone analogues with integrin, we mapped T4 and T4ac to the binding site of the RGD peptide. Modeling studies indicate that there is sufficient space in the cavity for the thyroid hormone to bind. Since the hormone is smaller in overall length than the RGD peptide, the hormone does not interact with the Arg recognition site in the propeller domain from alphav. In this model, most of the hormone interactions are with betaA domain of the integrin. Mutagenic studies can be carried out to validate the role of these residues in directing hormone interactions.     

Model of calcium channel inactivation: a qualitative analysis

A simple model of calcium channel inactivation has been developed, based on the accumulation of calcium ions at the inner mouth of the channel and on their binding to a receptor which inactivates the channel. A qualitative analysis has shown that upon an appropriate choice of parameters corresponding to the cell structure and to kinetic properties of its components, the calcium dependent inactivation and that assumed to be voltage dependent can both be emulated. The model suggests that the supposed variety of calcium channels might be explained by quantitative differences in nonlinear interactions of the channels with other cell components.     

Learning from Decoys to Improve the Sensitivity and Specificity of Proteomics Database Search Results

The statistical validation of database search results is a complex issue in bottom-up proteomics. The correct and incorrect peptide spectrum match (PSM) scores overlap significantly, making an accurate assessment of true peptide matches challenging. Since the complete separation between the true and false hits is practically never achieved, there is need for better methods and rescoring algorithms to improve upon the primary database search results. Here we describe the calibration and False Discovery Rate (FDR) estimation of database search scores through a dynamic FDR calculation method, FlexiFDR, which increases both the sensitivity and specificity of search results. Modelling a simple linear regression on the decoy hits for different charge states, the method maximized the number of true positives and reduced the number of false negatives in several standard datasets of varying complexity (18-mix, 49-mix, 200-mix) and few complex datasets (E. coli and Yeast) obtained from a wide variety of MS platforms. The net positive gain for correct spectral and peptide identifications was up to 14.81% and 6.2% respectively. The approach is applicable to different search methodologies- separate as well as concatenated database search, high mass accuracy, and semi-tryptic and modification searches. FlexiFDR was also applied to Mascot results and showed better performance than before. We have shown that appropriate threshold learnt from decoys, can be very effective in improving the database search results. FlexiFDR adapts itself to different instruments, data types and MS platforms. It learns from the decoy hits and sets a flexible threshold that automatically aligns itself to the underlying variables of data quality and size.     

Thermodynamic database for protein-nucleic acid interactions (ProNIT).

MOTIVATION: Protein-nucleic acid interactions are fundamental to the regulation of gene expression. In order to elucidate the molecular mechanism of protein-nucleic acid recognition and analyze the gene regulation network, not only structural data but also quantitative binding data are necessary. Although there are structural databases for proteins and nucleic acids, there exists no database for their experimental binding data. Thus, we have developed a Thermodynamic Database for Protein-Nucleic Acid Interactions (ProNIT). RESULTS: We have collected experimentally observed binding data from the literature. ProNIT contains several important thermodynamic data for protein-nucleic acid binding, such as dissociation constant (K(d)), association constant (K(a)), Gibbs free energy change (DeltaG), enthalpy change (DeltaH), heat capacity change (DeltaC(p)), experimental conditions, structural information of proteins, nucleic acids and the complex, and literature information. These data are integrated into a relational database system together with structural and functional information to provide flexible searching facilities by using combinations of various terms and parameters. A www interface allows users to search for data based on various conditions, with different display and sorting options, and to visualize molecular structures and their interactions. AVAILABILITY: ProNIT is freely accessible at the URL http://www.rtc.riken.go.jp/jouhou/pronit/pronit.html.     

KDBI: Kinetic Data of Bio-molecular Interactions database

Understanding of cellular processes and underlying molecular events requires knowledge about different aspects of molecular interactions, networks of molecules and pathways in addition to the sequence, structure and function of individual molecules involved. Databases of interacting molecules, pathways and related chemical reaction equations have been developed. The kinetic data for these interactions, which is important for mechanistic investigation, quantitative study and simulation of cellular processes and events, is not provided in the existing databases. We introduce a new database of Kinetic Data of Bio-molecular Interactions (KDBI) aimed at providing experimentally determined kinetic data of protein-protein, protein-RNA, protein-DNA, protein-ligand, RNA-ligand, DNA-ligand binding or reaction events described in the literature. KDBI contains information about binding or reaction event, participating molecules (name, synonyms, molecular formula, classification, SWISS-PROT AC or CAS number), binding or reaction equation, kinetic data and related references. The kinetic data is in terms of one or a combination of the following quantities as given in the literature of a particular event: association/dissociation or on/off rate constant, first/second/third/. order rate constant, equilibrium rate constant, catalytic rate constant, equilibrium association/dissociation constant, inhibition constant and binding affinity constant. Each entry can be retrieved through protein or nucleic acid or ligand name, SWISS-PROT AC number, ligand CAS number and full-text search of a binding or reaction event. KDBI currently contains 8273 entries of biomolecular binding or reaction events involving 1380 proteins, 143 nucleic acids and 1395 small molecules. Hyperlinks are provided for accessing references in Medline and available 3D structures in PDB and NDB. This database can be accessed at http://xin.cz3.nus.edu.sg/group/kdbi/kdbi.asp.     

Molecular interaction analysis in ligand design using mass transport,kinetic and thermodynamic methods

Ligand design in biotechnology is underpinned by the control of molecular affinity. Hence, measuring binding interactions is a key component in designing ligands for such uses as therapeutics, diagnostics, biomaterials and separation science. Mass transport, kinetic and thermodynamic methods have been used for macromolecular interaction analysis but also have potential applicability as direct methods for measuring small molecular interactions. They can enhance the ligand design process by providing the ability to choose ligands based on both their kinetic and thermodynamic binding properties.     

Kinetic and thermodynamic analyses of adhesion of a peptide,Trp-Lys-Tyr-Met-Val-D-Met (WKYMVm), and human formyl peptide receptor (hFPR)

The kinetic and thermodynamic properties of a peptide–receptor interaction was investigated by measuring the adhesion force in the reaction via atomic force microscopy (AFM). Trp-Lys-Tyr-Met-Val-D-Met (WKYMVm), considered as a model system in the present study, is a potent neutrophil chemo-attractant. Since being identified as an agonist for formyl peptide receptor (FPR), WKYMVm’s high affinity to FPR has been verified through investigation of its kinetic and physiological behaviors via conventional methods. However, there have been no reports on the adhesion force of WKYMVm-FPR. In this research, we measured the adhesion force of WKYMVm-FPR using AFM. Kinetic parameters obtained from the relationship between the adhesion force and loading rate were used to characterize the thermodynamic properties of WKYMVm-hFPR binding.     

Structural basis of neurophysin hormone specificity: Geometry, polarity, and polarizability in aromatic ring interactions

The structural origins of the specificity of the neurophysin hormone-binding site for an aromatic residue in peptide position 2 were explored by analyzing the binding of a series of peptides in the context of the crystal structure of liganded neurophysin. A new modeling method for describing the van der Waals surface of binding sites assisted in the analysis. Particular attention was paid to the unusually large (5 kcal/mol) difference in binding free energy between Phe and Leu in position 2, a value representing more than three times the maximum expected based on hydrophobicity alone, and additionally remarkable since modeling indicated that the Leu side chain was readily accommodated by the binding pocket. Although evidence was obtained of a weak thermodynamic linkage between the binding interactions of the residue 2 side chain and of the peptide alpha-amino group, two factors are considered central. (1) The bound Leu side chain can establish only one-third of the van der Waals contacts available to a Phe side chain. (2) The bound Phe side chain appears to be additionally stabilized relative to Leu by more favorable dipole and induced dipole interactions with nonaromatic polar and sulfur ligands in the binding pocket, as evidenced by examination of its interactions in the pocket, analysis of the detailed energetics of transfer of Phe and Leu side chains from water to other phases, and comparison with thermodynamic and structural data for the binding of residue 1 side chains in this system. While such polar interactions of aromatic rings have been previously observed, the present results suggest their potential for significant thermodynamic contributions to protein structure and ligand recognition.     

In-depth biophysical analysis of interactions between therapeutic antibodies and the extracellular domain of the epidermal growth factor receptor

Targeting of the epidermal growth factor receptor (EGFR) with monoclonal antibodies has become an established antitumor strategy in clinical use or in late stages of drug development. The mAbs effector mechanisms have been widely analyzed based on in vivo or cell studies. Hereby we intend to complement these functional studies by investigating the mAb-EGFR interactions on a molecular level. Surface plasmon resonance, isothermal titration calorimetry, and static light scattering were employed to characterize the interactions of matuzumab, cetuximab, and panitumumab with the extracellular soluble form ecEGFR. The kinetic and thermodynamic determinants dissected the differences in mAbs binding mechanism toward ecEGFR. The quantitative stoichiometric data clearly demonstrated the bivalent binding of the mAbs to two ecEGFR molecules. Our results complement earlier studies on simultaneous binding of cetuximab and matuzumab. The antibodies retain their bivalent binding mode achieving a 1:2:1 complex formation. Interestingly the binding parameters remain nearly constant for the individual antibodies in this ternary assembly. In contrast the binding of panitumumab is almost exclusive either by directly blocking the accessibility for the second antibody or by negative allosteric modulation. Overall we provide a comprehensive biophysical dataset on binding parameters, the complex assembly, and relative epitope accessibility for therapeutic anti-EGFR antibodies.     

Identifying Protein-protein Interaction Sites Using Peptide Arrays

Protein-protein interactions mediate most of the processes in the living cell and control homeostasis of the organism. Impaired protein interactions may result in disease, making protein interactions important drug targets. It is thus highly important to understand these interactions at the molecular level. Protein interactions are studied using a variety of techniques ranging from cellular and biochemical assays to quantitative biophysical assays, and these may be performed either with full-length proteins, with protein domains or with peptides. Peptides serve as excellent tools to study protein interactions since peptides can be easily synthesized and allow the focusing on specific interaction sites. Peptide arrays enable the identification of the interaction sites between two proteins as well as screening for peptides that bind the target protein for therapeutic purposes. They also allow high throughput SAR studies. For identification of binding sites, a typical peptide array usually contains partly overlapping 10-20 residues peptides derived from the full sequences of one or more partner proteins of the desired target protein. Screening the array for binding the target protein reveals the binding peptides, corresponding to the binding sites in the partner proteins, in an easy and fast method using only small amount of protein.In this article we describe a protocol for screening peptide arrays for mapping the interaction sites between a target protein and its partners. The peptide array is designed based on the sequences of the partner proteins taking into account their secondary structures. The arrays used in this protocol were Celluspots arrays prepared by INTAVIS Bioanalytical Instruments. The array is blocked to prevent unspecific binding and then incubated with the studied protein. Detection using an antibody reveals the binding peptides corresponding to the specific interaction sites between the proteins.