Dynamic flexibility of a peptide-binding groove of human HLA-DR1 class II MHC molecules: normal mode analysis of the antigen peptide-class II MHC complex

Class II major histocompatibility complex (MHC) has tolerance for binding longer antigen peptides than those bound by class I MHC. In this paper, a normal mode analysis on HLA-DR1 class II MHC involving an antigen peptide indicated that the peptide-binding groove had some different dynamic characteristics from that of HLA-A2 class I MHC. The dynamic changes in the class I groove with removal of the bound peptide were limited primarily to the central region and the C-terminal side (corresponding to the C-terminal side of the bound peptide) of the groove, while the dynamic changes in the class II groove with removal of the bound peptide extended to the whole of the groove, and were especially remarkable around a strand located in the N-terminal side (corresponding to the N-terminal side of the bound peptide) of the groove. These results suggest that the N-terminal side of the class II groove is more flexible than the same side of the class I groove, and this flexibility may allow some N-terminal residues of the bound peptide to extend outside the class II groove. Definite anti-correlative motions with removal of the bound peptide appeared between two alpha-helical regions of class II MHC as in the case of class I MHC. These motions of the class II groove may play an important role in obtaining "a flexible dynamic fit" against diverse longer peptides both of whose terminals extend outside the groove.     

Synthesis and biological evaluation of two chemically modified peptide epitopes for the class I MHC protein HLA-B*2705

The T-cell receptor of a CD8(+) T-cell recognises peptide epitopes bound by class I major histocompatibility complex (MHC) glycoproteins presented in a groove on their upper surface. Within the groove of the MHC molecule are 6 pockets, two of which mostly display a high degree of specificity for binding amino acids capable of making conserved and energetically favourable contacts with the MHC. One type of MHC molecule, HLA-B*2705, preferentially binds peptides containing an arginine at position 2. In an effort to increase the affinity of peptides for HLA-B*2705, potentially leading to better immune responses to such a peptide, we synthesised two modified epitopes where the amino acid at position 2 involved in anchoring the peptide to the class I molecule was replaced with the alpha-methylated beta,gamma-unsaturated arginine analogue 2-(S)-amino-5-guanidino-2-methyl-pent-3-enoic acid. The latter was prepared via a multi-step synthetic sequence, starting from alpha-methyl serine, and incorporated into dipeptides which were fragment-coupled to resin-bound heptameric peptides yielding the target nonameric sequences. Biological characterisation indicated that the modified peptides were poorer than the native peptides at stabilising empty class I MHC complexes, and cells sensitised with these peptides were not recognised as well by cognate CD8(+) T-cells, where available, compared to those sensitised with the native peptide. We suggest that the modifications made to the peptide have decreased its ability to bind to the peptide binding groove of HLA-B*2705 molecules which may explain the decrease in recognition by cytotoxic T-cells when compared to the native peptide.     

Computational analysis on the binding of epitope peptide to human leukocyte antigen class I molecule A*2402 subtype

Immunological response induced by small amino peptide has attracted much recent attention in the field of immunotherapy. Wilms' tumor (WT1) protein is one of the potent tumor antigens inducing immunological response in mouse and human, because WT1 is over expressed in many types of leukemia and various kinds of solid tumors. A 9-mer peptide encoded in WT1 protein (CMTWNQMNL; amino acid 235-243) is known to serve as antigenic peptide for human leukocyte antigen (HLA)-A*2402 molecule. It was reported that the replacement of the second amino residue, which is deeply responsible for the peptide binding to HLA, induced strong immunological response compared to the natural peptide. In this study, 19 kinds of single amino substitutions were introduced at position 2 of this 9-mer WT1 peptide. We performed molecular dynamics simulation on the complex of each of WT1 epitope peptides and HLA-β2 micro globulin (β2m) molecule, and subsequently estimated the binding affinity using molecular mechanics/generalized-Born surface area method combined with normal mode analysis. Our computation indicated that the peptide containing M2Y or M2W mutation showed high binding affinity to the HLA-β2m molecule as well as the natural peptide. We have also examined the role of the residue at position 2 in peptide binding to HLA-β2m. The calculation showed that van der Waals interaction between the side chain of the residue at position 2 and hydrophobic residues inside B-pocket of HLA are important. These findings will be helpful to search other potent peptides that will enhance strong immunological response specific to HLA-A*2402 molecule.     

Relative DNA binding affinity of helix 3 homeodomain analogues, major groove binders, can be rapidly screened by displacement of prebound ethidium bromide. A comparative study

The binding affinity for a 12-bp dsDNA of Antennapedia helix 3 analogues, major groove binders, has been measured by displacement of prebound ethidium bromide, a fluorescent displacement assay proposed for minor groove binders by Boger et al.(J. Am. Chem. Soc., 2000, 122, 6382-6394). Relative binding affinities determined by this method were compared to those obtained by gel mobility shift and footprinting assays for the 12-bp dsDNA and a 178-bp DNA fragment. The present work demonstrates that the fluorescence displacement assay is suitable for rapid screening of major groove binders, even though about 60 to 70% of the prebound ethidium bromide is displaced by these peptides. Total (100%) displacement of ethidium bromide was serendipitously achieved by addition in the peptide sequence, at the N-terminus, of a S-3-nitro-2-pyridinesulfenyl-N-acetyl-cysteine residue. S-3-nitro-2-pyridinesulfenylcysteine was shown to (i) bind to dsDNA with a micromolar affinity and (ii) direct within DNA grooves a peptide with no affinity for dsDNA.     

Predicting class I major histocompatibility complex (MHC) binders using multivariate statistics: comparison of discriminant analysis and multiple linear regression

The accurate in silico identification of T-cell epitopes is a critical step in the development of peptide-based vaccines, reagents, and diagnostics. It has a direct impact on the success of subsequent experimental work. Epitopes arise as a consequence of complex proteolytic processing within the cell. Prior to being recognized by T cells, an epitope is presented on the cell surface as a complex with a major histocompatibility complex (MHC) protein. A prerequisite therefore for T-cell recognition is that an epitope is also a good MHC binder. Thus, T-cell epitope prediction overlaps strongly with the prediction of MHC binding. In the present study, we compare discriminant analysis and multiple linear regression as algorithmic engines for the definition of quantitative matrices for binding affinity prediction. We apply these methods to peptides which bind the well-studied human MHC allele HLA-A*0201. A matrix which results from combining results of the two methods proved powerfully predictive under cross-validation. The new matrix was also tested on an external set of 160 binders to HLA-A*0201; it was able to recognize 135 (84%) of them.     

A polymorphic pocket at the P10 position contributes to peptide binding specificity in class II MHC proteins

Peptides bind to class II major histocompatibility complex (MHC) proteins in an extended conformation. Pockets in the peptide binding site spaced to accommodate peptide side chains at the P1, P4, P6, and P9 positions have been previously characterized and help to explain the obtained peptide binding specificity. However, two peptides differing only at P10 have significantly different binding affinities for HLA-DR1. The structure of HLA-DR1 in complex with the tighter binding peptide shows that the peptide binds in the usual polyproline type II conformation, but with the P10 residue accommodated in a shallow pocket at the end of the binding groove. HLA-DR1 variants with polymorphic residues at these positions were produced and found to exhibit different side chain specificity at the P10 position. These results define a new specificity position in HLA-DR proteins.     

Glycan side chains on naturally presented MHC class II ligands

The molecular characterization of unknown naturally presented major histocompatibility complex (MHC) class II glycopeptides carrying complex glycans has so far not been achieved, reflecting the different fragmentation characteristics of sugars and peptides in mass spectrometric analysis. Human leukocyte antigen (HLA)-DR-bound peptides were isolated by affinity purification, separated via high performance liquid chromatography and analyzed by matrix-assisted laser desorption/ionization and electrospray ionization mass spectrometry. We were able to identify two naturally processed MHC class II ligands, CD53(122-136) and CD53(121-136), carrying complex N-linked glycan side chains by a combination of in-source and collision-induced fragmentation on a quadrupole time-of-flight tandem mass spectrometer.     

In silico designing breast cancer peptide vaccine for binding to MHC class I and II: A molecular docking study

Antigenic peptides or cancer peptide vaccines can be directly delivered to cancer patients to produce immunologic responses against cancer cells. Specifically, designed peptides can associate with Major Histocompatibility Complex (MHC) class I or II molecules on the cell surface of antigen presenting cells activating anti-tumor effector mechanisms by triggering helper T cell (Th) or cytotoxic T cells (CTL). In general, high binding to MHCs approximately correlates with in vivo immunogenicity. Consequently, a molecular docking technique was run on a library of novel discontinuous peptides predicted by PEPOP from Human epidermal growth factor receptor 2 (HER2 ECD) subdomain III. This technique is expected to improve the prediction accuracy in order to identify the best MHC class I and II binder peptides. Molecular docking analysis through GOLD identified the peptide 1412 as the best MHC binder peptide to both MHC class I and II molecules used in the study. The GOLD results predicted HLA-DR4, HLA-DP2 and TCR as the most often targeted receptors by the peptide 1412. These findings, based on bioinformatics analyses, can be exploited in further experimental analyses in vaccine design and cancer therapy to find possible proper approaches providing beneficial effects.     

抗原肽与MHC分子相互作用的QSAR模型研究

分别采用氨基酸序列及氨基酸结构描述符方法定量研究了短肽-MHC(major histocompatibility complex)分子结合亲合力的定量构效关系(QSAR)模型, 用这两个模型对短肽与HLA-A*0201分子结合的805个预测样本进行了预测, 预测准确度分别达到66.8%和65.5%. 采用了去除“劣点”方法检验模型的鲁棒性, 结果证明两个模型都具有良好的鲁棒性. 通过改变两个模型的参数, 对氨基酸残基个数为8或10的CTL(cytotoxicy T lymphocyte)表位进行预测.     

Efficient estimation of binding free energies between peptides and an MHC class II molecule using coarse‐grained molecular dynamics simulations with a weighted histogram analysis method

We estimate the binding free energy between peptides and an MHC class II molecule using molecular dynamics (MD) simulations with the weighted histogram analysis method (WHAM). We show that, owing to its more thorough sampling in the available computational time, the binding free energy obtained by pulling the whole peptide using a coarse‐grained (CG) force field (MARTINI) is less prone to significant error induced by inadequate‐sampling than using an atomistic force field (AMBER). We further demonstrate that using CG MD to pull 3–4 residue peptide segments while leaving the remaining peptide segments in the binding groove and adding up the binding free energies of all peptide segments gives robust binding free energy estimations, which are in good agreement with the experimentally measured binding affinities for the peptide sequences studied. Our approach thus provides a promising and computationally efficient way to rapidly and reliably estimate the binding free energy between an arbitrary peptide and an MHC class II molecule. © 2017 Wiley Periodicals, Inc.     

Design and synthesis of peptides that bind alpha-bungarotoxin with high affinity

BACKGROUND: Alpha-bungarotoxin (alpha-BTX) is a highly toxic snake venom alpha-neurotoxin that binds to acetylcholine receptor (AChR) at the neuromuscular junction, and is a potent inhibitor of this receptor. We describe the design and synthesis of peptides that bind alpha-BTX with high affinity, and inhibit its interaction with AChR with an IC(50) of 2 nM. The design of these peptides was based on a lead peptide with an IC(50) of 3x10(-7) M, previously identified by us [M. Balass et al., Proc. Natl. Acad. Sci. USA 94 (1997) 6054] using a phage-display peptide library. RESULTS: Employing nuclear magnetic resonance-derived structural information [T. Scherf et al., Proc. Natl. Acad. Sci. USA 94 (1997) 6059] of the complex of alpha-BTX with the lead peptide, as well as structure-function analysis of the ligand-binding site of AChR, a systematic residue replacement of the lead peptide, one position at a time, yielded 45 different 13-mer peptides. Of these, two peptides exhibited a one order of magnitude increase in inhibitory potency in comparison to the lead peptide. The design of additional peptides, with two or three replacements, resulted in peptides that exhibited a further increase in inhibitory potency (IC(50) values of 2 nM), that is more than two orders of magnitude better than that of the original lead peptide, and better than that of any known peptide derived from AChR sequence. The high affinity peptides had a protective effect on mice against alpha-BTX lethality. CONCLUSIONS: Synthetic peptides with high affinity to alpha-BTX may be used as potential lead compounds for developing effective antidotes against alpha-BTX poisoning. Moreover, the procedure employed in this study may serve as a general approach for the design and synthesis of peptides that interact with high affinity with any desired biological target.     

Molecular modeling of class I and II alleles of the major histocompatibility complex in <Emphasis Type="Italic">Salmo salar</Emphasis>

Knowledge of the 3D structure of the binding groove of major histocompatibility (MHC) molecules, which play a central role in the immune response, is crucial to shed light into the details of peptide recognition and polymorphism. This work reports molecular modeling studies aimed at providing 3D models for two class I and two class II MHC alleles from Salmo salar (Sasa), as the lack of experimental structures of fish MHC molecules represents a serious limitation to understand the specific preferences for peptide binding. The reliability of the structural models built up using bioinformatic tools was explored by means of molecular dynamics simulations of their complexes with representative peptides, and the energetics of the MHC-peptide interaction was determined by combining molecular mechanics interaction energies and implicit continuum solvation calculations. The structural models revealed the occurrence of notable differences in the nature of residues at specific positions in the binding groove not only between human and Sasa MHC proteins, but also between different Sasa alleles. Those differences lead to distinct trends in the structural features that mediate the binding of peptides to both class I and II MHC molecules, which are qualitatively reflected in the relative binding affinities. Overall, the structural models presented here are a valuable starting point to explore the interactions between MHC receptors and pathogen-specific interactions and to design vaccines against viral pathogens.     

Fine specificity of antigen binding to two class I major histocompatibility proteins (B*2705 and B*2703) differing in a single amino acid residue

Starting from the X-ray structure of a class I majorhistocompatibility complex (MHC)-encoded protein (HLA-B*2705), a naturallypresented self-nonapeptide and two synthetic analogues were simulated in thebinding groove of two human leukocyte antigen (HLA) alleles (B*2703 andB*2705) differing in a single amino acid residue. After 200 ps moleculardynamics simulations of the solvated HLA–peptide pairs, some molecularproperties of the complexes (distances between ligand and protein center ofmasses, atomic fluctuations, buried versus accessible surface areas,hydrogen-bond frequencies) allow a clear discrimination of potent from weakMHC binders. The binding specificity of the three nonapeptides for the twoHLA alleles could be explained by the disruption of one hydrogen-bondingnetwork in the binding pocket of the HLA-B*2705 protein where the singlemutation occurs. Rearrangements of interactions in the B pocket, which bindsthe side chain of peptidic residue 2, and a weakening of interactionsinvolving the C-terminal end of the peptide also took place. In addition,extension of the peptide backbone using a -Ala analogue did notabolish binding to any of the two HLA-B27 subtypes, but increased theselectivity for B*2703, as expected from the larger peptide binding groovein this subtype. A better understanding of the atomic details involved inpeptide selection by closely related HLA alleles is of crucial importancefor unraveling the molecular features linking particular HLA alleles toautoimmune diseases, and for the identification of antigenic peptidestriggering such pathologies.     

A critical cross-validation of high throughput structural binding prediction methods for pMHC

T-cells recognize antigens via their T-cell receptors. The major histocompatibility complex (MHC) binds antigens in a specific way, transports them to the surface and presents the peptides to the TCR. Many in silico approaches have been developed to predict the binding characteristics of potential T-cell epitopes (peptides), with most of them being based solely on the amino acid sequence. We present a structural approach which provides insights into the spatial binding geometry. We combine different tools for side chain substitution (threading), energy minimization, as well as scoring methods for protein/peptide interfaces. The focus of this study is on high data throughput in combination with accurate results. These methods are not meant to predict the accurate binding free energy but to give a certain direction for the classification of peptides into peptides that are potential binders and peptides that definitely do not bind to a given MHC structure. In total we performed approximately 83,000 binding affinity prediction runs to evaluate interactions between peptides and MHCs, using different combinations of tools. Depending on the tools used, the prediction quality ranged from almost random to around 75% of accuracy for correctly predicting a peptide to be either a binder or a non-binder. The prediction quality strongly depends on all three evaluation steps, namely, the threading of the peptide, energy minimization and scoring.     

Molecular dynamics simulations of pro-apoptotic BH3 peptide helices in aqueous medium: relationship between helix stability and their binding affinities to the anti-apoptotic protein Bcl-X<Subscript>L</Subscript>

The B-cell lymphoma 2 (Bcl-2) family of proteins regulates the intrinsic pathway of apoptosis. Interactions between specific anti- and pro-apoptotic Bcl-2 proteins determine the fate of a cell. Anti-apoptotic Bcl-2 proteins have been shown to be over-expressed in certain cancers and they are attractive targets for developing anti-cancer drugs. Peptides from the BH3 region of pro-apoptotic proteins have been shown to interact with anti-apoptotic Bcl-2 proteins and induce biological activity similar to that observed in parent proteins. However, the specificity of BH3 peptides derived from different pro-apoptotic proteins differ for different anti-apoptotic Bcl-2 proteins. In this study, we have investigated the relationship between the stable helical nature of BH3 peptides and their affinities to Bcl-X_L, an anti-apoptotic Bcl-2 protein. We have carried out molecular dynamics simulations of six BH3 peptides derived from Bak, Bad and Bim pro-apoptotic proteins for a period of 50 ns each in aqueous medium. Due to the amphipathic nature of BH3 peptides, the hydrophobic residues on the hydrophobic face tend to cluster together in all BH3 peptides. While this process resulted in a complete loss of helical structure in 16-mer Bak and 16-mer Bad wild type peptides, stabilizing interactions in the hydrophilic face of the BH3 peptides and capping interactions helped to maintain partial helical character in 16-mer Bad mutant and 16-mer Bim peptides. The latter two 16-mer peptides exhibit higher affinity for Bcl-X_L. Similarly the longer BH3 peptides, 25-mer Bad and 33-mer Bim, also resulted in smaller and stable helical fragments and their helical conformation is stabilized by interactions between residues in the solvent-exposed hydrophilic half of the peptide. The stable nature of helical segment in a BH3 peptide can be directly correlated to its binding affinity and the helical region encompassed the highly conserved Leu residue. We propose that upon approaching the hydrophobic groove of anti-apoptotic proteins, a longer helix will be induced in high affinity BH3 peptides by extending the smaller stable helical segments around the conserved Leu residue in both N- and C-terminal regions. The results reported in this study will have implications in developing peptide-based inhibitors for anti-apoptotic Bcl-2 proteins.     

Does adsorption at hydroxyapatite surfaces induce peptide folding? Insights from large-scale B3LYP calculations

Large-scale periodic quantum mechanical calculations (509 atoms, 7852 atomic orbitals) based on the hybrid B3LYP functional focused on the peptide folding induced by the adsorption on the (001) and (010) hydroxyapatite (HA) surfaces give interesting insights on the role of specific interactions between surface sites and the peptide, which stabilize the helix conformation over the "native" random coil ones for in silico designed model peptides. The two peptides were derived from the 12-Gly oligomer, with one (P1, C-tGGKGGGGGGEGGN-t) and two (P2, C-tGGKGGKEGGEGGN-t) glutamic acid (E) and lysine (K) residue mutations. The most stable gas-phase "native" conformation for both peptides resulted in a random coil (RC) structure, with the helix (H) conformation being ≈100 kJ mol(-1) higher in free energy. The two peptide conformations interact with the HA (001) and (010) surfaces by C═O groups via Ca(2+) ions, by hydrogen bond between NH(2) groups and the basic PO(4)(3-) groups and by a relevant fraction due to dispersion forces. Peptide adsorption was studied on the dry (001) surface, the wet one envisaging 2 H(2)O per surface Ca(2+) and, on the latter, also considering the adsorption of microsolvated peptides with 4 H(2)O molecules located at sites responsible of the interaction with the surface. The P1 mutant does prefer to be adsorbed as a random coil by ≈160 kJ/mol, whereas the reverse is computed for P2, preferring the helix conformation by ≈50 kJ/mol. Adsorption as helix of both P1 and P2 mutants brings about proton transfer toward the HA surfaces with a large charge transfer component to the interaction energy.     

Design, synthesis and evaluation of β-lactam antigenic peptide hybrids; unusual opening of the β-lactam ring in acidic media

β-Lactam peptides were envisioned as conformational constraints in antigenic peptides (APs). Three different β-lactam tripeptides of varying flexibility were prepared in solution and incorporated in place of the central part of the altered melanoma associated antigenic peptide Leu(27)-Melan-A(26-35) using solid phase synthesis techniques. Upon TFA cleavage from the solid support, an unexpected opening of the β-lactam ring occurred with conservation of the amide bond. After adaptation of the solid phase synthesis strategy, β-lactam peptides were successfully obtained and both opened and closed forms were evaluated for their capacity to bind to the antigen-presenting class-I MHC HLA-A2 protein system. None of the closed β-lactam peptides bound to HLA-A2, but their opened variants were shown to be moderate to good HLA-A2 ligands, one of them being even capable of stimulating a Melan-A-specific T cell line.     

Evaluation of the membrane-binding properties of the proximal region of the angiotensin II receptor (AT1A) carboxyl terminus by surface plasmon resonance.

The proximal region of the angiotensin II receptor (AT1A) carboxyl-terminus (known as helix VIII) is important for receptor function. In this study, we used surface plasmon resonance (SPR) to examine the interaction of helix VIII-derived peptides with three model lipid membranes. The membrane-binding properties of these synthetic peptides, as well as a series of peptide analogues with modified amino acid sequences, could be explained by both amino acid sequence and kinetic binding data by SPR. The helix VIII peptides showed a higher affinity for lipid membranes that contained negatively charged phospholipid, rather than zwitterionic phospholipid. The findings of an SPR study may be useful for estimating the cooperative binding of intracellular receptor domains with G proteins and the components of the lipid bilayer.     

Probing the role of the C-H...O hydrogen bond stabilized polypeptide chain reversal at the C-terminus of designed peptide helices. Structural characterization of three decapeptides

The structural characterization in crystals of three designed decapeptides containing a double d-segment at the C-terminus is described. The crystal structures of the peptides Boc-Leu-Aib-Val-Xxx-Leu-Aib-Val-(D)Ala-(D)Leu-Aib-OMe, (Xxx = Gly 2, (D)Ala 3, Aib 4) have been determined and compared with those reported earlier for peptide 1 (Xxx = Ala) and the all l analogue Boc-Leu-Aib-Val-Ala-Leu-Aib-Val-Ala-Leu-Aib-OMe, which yielded a perfect right-handed alpha-helical structure. Peptides 1 and 2 reveal a right-handed helical segment spanning residues 1 to 7, ending in a Schellman motif with (D)Ala(8) functioning as the terminating residue. Polypeptide chain reversal occurs at residue 9, a novel feature that appears to be the consequence of a C-H.O hydrogen bond between residue 4 C(alpha)H and residue 9 CO groups. The structures of peptides 3 and 4, which lack the pro R hydrogen at the C(alpha) atom of residue 4, are dramatically different. Peptide 3 adopts a right-handed helical conformation over the 1 to 7 segment. Residues 8 and 9 adopt alpha(L) conformations forming a C-terminus type I' beta-turn, corresponding to an incipient left-handed twist of the polypeptide chain. In peptide 4, helix termination occurs at Aib(6), with residues 6 to 9 forming a left-handed helix, resulting in a structure that accommodates direct fusion of two helical segments of opposite twist. Peptides 3 and 4 provide examples of chiral residues occurring in the less favored sense of helical twist; (D)Ala(4) in peptide 3 adopts an alpha(R) conformation, while (L)Val(7) in 4 adopts an alpha(L) conformation. The structural comparison of the decapeptides reported here provides evidence for the role of specific C-H.O hydrogen bonds in stabilizing chain reversals at helix termini, which may be relevant in aligning contiguous helical and strand segments in polypeptide structures.