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

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

On the hydrophobicity of peptides:comparing empirical predictions of peptide log P values

Peptides are of great therapeutic potential as vaccines and drugs. Knowledge of physicochemical descriptors, including the partition coefficient logP, is useful for the development of predictive Quantitative Structure-Activity Relationships (QSARs). We have investigated the accuracy of available programs for the prediction of logP values for peptides with known experimental values obtained from the literature. Eight prediction programs were tested, of which seven programs were fragment-based methods: XLogP, LogKow, PLogP, ACDLogP, AlogP, Interactive Analysis's LogP and MlogP; and one program used a whole molecule approach: QikProp. The predictive accuracy of the programs was assessed using r(2) values, with ALogP being the most effective (r( 2) = 0.822) and MLogP the least (r(2) = 0.090). We also examined three distinct types of peptide structure: blocked, unblocked, and cyclic. For each study (all peptides, blocked, unblocked and cyclic peptides) the performance of programs rated from best to worse is as follows: all peptides - ALogP, QikProp, PLogP, XLogP, IALogP, LogKow, ACDLogP, and MlogP; blocked peptides - PLogP, XLogP, ACDLogP, IALogP, LogKow, QikProp, ALogP, and MLogP; unblocked peptides - QikProp, IALogP, ALogP, ACDLogP, MLogP, XLogP, LogKow and PLogP; cyclic peptides - LogKow, ALogP, XLogP, MLogP, QikProp, ACDLogP, IALogP. In summary, all programs gave better predictions for blocked peptides, while, in general, logP values for cyclic peptides were under-predicted and those of unblocked peptides were over-predicted.

Structure-function analysis of amino acid 74 of human RAMP1 and RAMP3 and its role in peptide interactions with adrenomedullin and calcitonin gene-related peptide receptors

The receptors for calcitonin gene-related peptide (CGRP) and adrenomedullin (AM) are complexes of the calcitonin receptor-like receptor (CLR) and receptor activity-modifying proteins (RAMP). The CGRP receptor is a CLR/RAMP1 pairing whereas CLR/RAMP2 and CLR/RAMP3 constitute two subtypes of AM receptor: AM(1) and AM(2), respectively. Previous studies identified Glu74 in RAMP3 to be important for AM binding and potency. To further understand the importance of this residue and its equivalent in RAMP1 (Trp74) we substituted the native amino acids with several others. In RAMP3, these were Trp, Phe, Tyr, Ala, Ser, Thr, Arg and Asn; in RAMP1, Glu, Phe, Tyr, Ala and Asn substitutions were made. The mutant RAMPs were co-expressed with CLR in Cos7 cells; receptor function in response to AM, AM(2)/intermedin and CGRP was measured in a cAMP assay and cell surface expression was determined by ELISA. Phe reduced AM potency in RAMP3 but had no effect in RAMP1. In contrast, Tyr had no effect in RAMP3 but enhanced AM potency in RAMP1. Most other substitutions had a small effect on AM potency in both receptors whereas there was little impact on CGRP or AM(2) potency. Overall, these data suggest that the geometry and charge of the residue at position 74 contribute to how AM interacts with the AM(2) and CGRP receptors and confirms the role of this position in dictating differential AM pharmacology at the AM(2) and CGRP receptors.

Structure-activity relationships for α-calcitonin gene-related peptide

Calcitonin gene-related peptide (CGRP) is a member of the calcitonin (CT) family of peptides. It is a widely distributed neuropeptide implicated in conditions such as neurogenic inflammation. With other members of the CT family, it shares an N-terminal disulphide-bonded ring which is essential for biological activity, an area of potential α-helix, and a C-terminal amide. CGRP binds to the calcitonin receptor-like receptor (CLR) in complex with receptor activity-modifying protein 1 (RAMP1), a member of the family B (or secretin-like) GPCRs. It can also activate other CLR or calcitonin-receptor/RAMP complexes. This 37 amino acid peptide comprises the N-terminal ring that is required for receptor activation (residues 1-7); an α-helix (residues 8-18), a region incorporating a β-bend (residues 19-26) and the C-terminal portion (residues 27-37), that is characterized by bends between residues 28-30 and 33-34. A few residues have been identified that seem to make major contributions to receptor binding and activation, with a larger number contributing either to minor interactions (which collectively may be significant), or to maintaining the conformation of the bound peptide. It is not clear if CGRP follows the pattern of other family B GPCRs in binding largely as an α-helix. Linked Articles This article is part of a themed section on Neuropeptides. To view the other articles in this section visit http://dx.doi.org/10.1111/bph.2013.170.issue-7 © 2012 The Authors. British Journal of Pharmacology © 2012 The British Pharmacological Society.

Structure-activity relationships of the N-terminus of calcitonin gene-related peptide:key roles of alanine-5 and threonine-6 in receptor activation

Background and purpose - The N-terminus of calcitonin gene-related peptide (CGRP) is important for receptor activation, especially the disulphide-bonded ring (residues 1-7). However, the roles of individual amino acids within this region have not been examined and so the molecular determinants of agonism are unknown. This study has examined the role of residues 1, 3-6 and 8-9, excluding Cys-2 and Cys-7. Experimental approach - CGRP derivatives were substituted with either cysteine or alanine; further residues were introduced at position 6. Their affinity was measured by radioligand binding and their efficacy by measuring cAMP production in SK-N-MC cells and ß-arrestin 2 translocation in CHO-K1 cells at the CGRP receptor. Key results - Substitution of Ala-5 by cysteine reduced affinity 270-fold and reduced efficacy for production of cAMP in SK-N-MCs. Potency at ß-arrestin translocation was reduced by 9-fold. Substitution of Thr-6 by cysteine destroyed all measurable efficacy of both cAMP and ß-arrestin responses; substitution with either alanine or serine impaired potency. Substitutions at positions 1, 4, 8 and 9 resulted in approximately 10-fold reductions in potency at both responses. Similar observations were made at a second CGRP-activated receptor, the AMY1(a) receptor. Conclusions and implications - Ala-5 and Thr-6 are key determinants of agonist activity for CGRP. Ala-5 is also very important for receptor binding. Residues outside of the 1-7 ring also contribute to agonist activity.

Structure-activity relationship of calcitonin gene-related peptide and its receptor

DUE TO COPYRIGHT RESTRICTIONS ONLY AVAILABLE FOR CONSULTATION AT ASTON UNIVERSITY LIBRARY AND INFORMATION SERVICES WITH PRIOR ARRANGEMENT

Toward the prediction of class I and II mouse major histocompatibility complex-peptide-binding affinity:in silico bioinformatic step-by-step guide using quantitative structure-activity relationships

Quantitative structure-activity relationship (QSAR) analysis is a cornerstone of modern informatics. Predictive computational models of peptide-major histocompatibility complex (MHC)-binding affinity based on QSAR technology have now become important components of modern computational immunovaccinology. Historically, such approaches have been built around semiqualitative, classification methods, but these are now giving way to quantitative regression methods. We review three methods--a 2D-QSAR additive-partial least squares (PLS) and a 3D-QSAR comparative molecular similarity index analysis (CoMSIA) method--which can identify the sequence dependence of peptide-binding specificity for various class I MHC alleles from the reported binding affinities (IC50) of peptide sets. The third method is an iterative self-consistent (ISC) PLS-based additive method, which is a recently developed extension to the additive method for the affinity prediction of class II peptides. The QSAR methods presented here have established themselves as immunoinformatic techniques complementary to existing methodology, useful in the quantitative prediction of binding affinity: current methods for the in silico identification of T-cell epitopes (which form the basis of many vaccines, diagnostics, and reagents) rely on the accurate computational prediction of peptide-MHC affinity. We have reviewed various human and mouse class I and class II allele models. Studied alleles comprise HLA-A*0101, HLA-A*0201, HLA-A*0202, HLA-A*0203, HLA-A*0206, HLA-A*0301, HLA-A*1101, HLA-A*3101, HLA-A*6801, HLA-A*6802, HLA-B*3501, H2-K(k), H2-K(b), H2-D(b) HLA-DRB1*0101, HLA-DRB1*0401, HLA-DRB1*0701, I-A(b), I-A(d), I-A(k), I-A(S), I-E(d), and I-E(k). In this chapter we show a step-by-step guide into predicting the reliability and the resulting models to represent an advance on existing methods. The peptides used in this study are available from the AntiJen database (http://www.jenner.ac.uk/AntiJen). The PLS method is available commercially in the SYBYL molecular modeling software package. The resulting models, which can be used for accurate T-cell epitope prediction, will be made are freely available online at the URL http://www.jenner.ac.uk/MHCPred.

Structure-function analysis of RAMP1 by Alanine Mutagenesis

Receptor activity modifying protein 1 (RAMP1) is an integral component of several receptors including the calcitonin gene-related peptide (CGRP) receptor. It forms a complex with the calcitonin receptor-like receptor (CLR) and is required for receptor trafficking and ligand binding. The N-terminus of RAMP1 comprises three helices. The current study investigated regions of RAMP1 important for CGRP or CLR interactions by alanine mutagenesis. Modeling suggested the second and third helices were important in protein-protein interactions. Most of the conserved residues in the N-terminus (M48, W56, Y66, P85, N66, H97, F101, D113, P114, P115), together with a further 13 residues spread throughout three helices of RAMP1, were mutated to alanine and coexpressed with CLR in Cos 7 cells. None of the mutations significantly reduced RAMP expression. Of the nine mutants from helix 1, only M48A had any effect, producing a modest reduction in trafficking of CLR to the cell surface. In helix 2 Y66A almost completely abolished CLR trafficking; L69A and T73A reduced the potency of CGRP to produce cAMP. In helix 3, H97A abolished CLR trafficking; P85A, N86A, and F101A had caused modest reductions in CLR trafficking and also reduced the potency of CGRP on cAMP production. F93A caused a modest reduction in CLR trafficking alone and L94A increased cAMP production. The data are consistent with a CLR recognition site particularly involving Y66 and H97, with lesser roles for adjacent residues in helix 3. L69 and T73 may contribute to a CGRP recognition site in helix 2 also involving nearby residues.

Structure-function analysis of RAMP1-RAMP3 chimeras

The role of receptor activity modifying protein 1 (RAMP1) in forming receptors with the calcitonin receptor-like receptor (CLR) and the calcitonin receptor (CTR) was examined by producing chimeras between RAMP1 and RAMP3. RAMPs have three extracellular helices. Exchange of helix 1 of the RAMPs or residues 62-69 in helix 2 greatly reduced CLR trafficking (a marker for CLR association). Modeling suggests that these exchanges alter the CLR recognition site on RAMP1, which is more exposed than on RAMP3. Exchange of residues 86-89 of RAMP1 had no effect on the trafficking of CLR but reduced the potency of human (h) alphaCGRP and adrenomedullin. However, these alterations to RAMP1 had no effect on the potency of hbetaCGRP. These residues of RAMP1 lie at the junction of helix 3 and its connecting loop with helix 2. Modeling suggests that the loop is more exposed in RAMP1 than RAMP3; it may play an important role in peptide binding, either directly or indirectly. Exchange of residues 90-94 of RAMP1 caused a modest reduction in CLR expression and a 15-fold decrease in CGRP potency. It is unlikely that the decrease in expression is enough to explain the reduction in potency, and so these may have dual roles in recognizing CLR and CGRP. For CTR, only 6 out of 26 chimeras covering the extracellular part of RAMP1 did not reduce agonist potency. Thus the association of CTR with RAMP1 seems more sensitive to changes in RAMP1 structure induced by the chimeras than is CLR.

Tethering a laminin peptide to a crosslinked collagen scaffold for biofunctionality

Cell adhesion peptide regulates various cellular functions like proliferation, attachment, and spreading. The cellular response to laminin peptide (PPFLMLLKGSTR), a motif of laminin-5 alpha3 chain, tethered to type I collagen, crosslinked using microbial transglutaminase (mTGase) was investigated. mTGase is an enzyme that initiates crosslinking by reacting with the glutamine and lysine residues on the collagen fibers stabilizing the molecular structure. In this study that tethering of the laminin peptide in a mTGase crosslinked collagen scaffold enhanced cell proliferation and attachment. Laminin peptide tethered crosslinked scaffold showed unaltered cell morphology of 3T3 fibroblasts when compared with collagen and crosslinked scaffold. The triple helical structure of collagen remained unaltered by the addition of laminin peptide. In addition a dose-dependent affinity of the laminin peptide towards collagen was seen. The degree of crosslinking was measured by amino acid analysis, differential scanning calorimeter and fourier transform infrared spectroscopy. Increased crosslinking was observed in mTGase crosslinked group. mTGase crosslinking showed higher shrinkage temperature. There was alteration in the fibrillar architecture due to the crosslinking activity of mTGase. Hence, the use of enzyme-mediated linking shows promise in tethering cell adhesive peptides through biodegradable scaffolds.

Peptide binding prediction for the human class II MHC allele HLA-DP2:a molecular docking approach

MHC class II proteins bind oligopeptide fragments derived from proteolysis of pathogen antigens, presenting them at the cell surface for recognition by CD4+ T cells. Human MHC class II alleles are grouped into three loci: HLA-DP, HLA-DQ and HLA-DR. In contrast to HLA-DR and HLA-DQ, HLA-DP proteins have not been studied extensively, as they have been viewed as less important in immune responses than DRs and DQs. However, it is now known that HLA-DP alleles are associated with many autoimmune diseases. Quite recently, the X-ray structure of the HLA-DP2 molecule (DPA*0103, DPB1*0201) in complex with a self-peptide derived from the HLA-DR a-chain has been determined. In the present study, we applied a validated molecular docking protocol to a library of 247 modelled peptide-DP2 complexes, seeking to assess the contribution made by each of the 20 naturally occurred amino acids at each of the nine binding core peptide positions and the four flanking residues (two on both sides).

Characterisation of calcitonin generelated peptide (CGRP) and amylin receptors

The aims of this study were to examine the binding characteristics of the rat CGRP receptor and to further the classification of CGRP and amylin receptors in guinea-pig tissue preparations. Binding characteristics of CGRP were investigated on rat splenic, cerebellar and liver membrane preparations. Human-α-CGRP, rat-α-CGRP and the CGRP receptor analogues Tyrº -CGRPC28-37) and [Cys (ACM)2,7 ]-human CGRP and the CGRP receptor antagonist CGRPC8-37) were utilised in competitive radioligand binding experiments to identify possible CGRP receptor subtypes in these tissues. There appeared to be no significant differences between the rat CGRP receptors examined. A panel of monoclonal antibodies (Mabs) raised against CGRP were employed to investigate the structure-activity relationships of CGRP and its receptor. No differences between the tissue receptors were observed using this panel of Mabs. The effects of human-α, human-β, rat-α-CGRP, human and rat amylin and adrenomedullin(13-52) were examined on the spontaneously beating right atria and on electrically evoked twitch contractions of isolated guinea-pig ileum, vas deferens and left atria. All of the peptides caused concentration-dependent inhibition of twitch amplitude in the ileum and vas deferens. CGRP produced positive inotropic effects in the right and left atria and positive chronotropic effects in the right atria. A variety of CGRP receptor antagonists and putative amylin receptor antagonists were used to antagonise these effects. CGRP(8-37) is currently used as a basis for CGRP receptor classification (Dennis, et al., 1989). Based upon results obtained using CGRP(8-37) it has been shown that the guinea-pig ileum contains mainly CGRP 1 receptors and the vas deferens contain CGRP2 receptors. Amylin was shown to act at receptors distinct from those for CGRP and it is postulated that amylin has its own receptors in these preparations. Experiments using CGRP (19-37) and Tyrº -CGRP(28-37) indicate that human and rat CGRP act at distinct receptors in guinea-pig ileum and vas deferens. The amylin receptor antagonist amylin(8-37) and the putative antagonist AC187 provide evidence to suggest human and rat amylin also act at receptors able to distinguish between the two types of amylin.

An insight into the activation mechanism of the calcatonin-gene-related peptide receptor

The calcitonin-gene- related peptide (CGRP) receptor is unique among G-protein coupled receptors (GPCRs) as it consists of at least three proteins: calcitonin receptor like receptor (CLR), receptor activity modifying protein (RAMP)1 and receptor component protein (RCP). An endogenous agonist for this curious receptor is aCGRP, which is a sensory nerve-derived peptide made up of 37 amino acids. aCGRP acts as a potent vasodilator having pronounced effects on arterioles and capillaries. Understanding the pharmacodynamics of the CGRP receptor may have pharmaceutical benefit as the receptor has been associated with the onset of migraines and implicated in Raynauds syndrome. The primary aim of this thesis was to identify functionally important residues in the extracellular face of the CGRP receptor. Three areas of interest were selected including the extreme N-terminus of the CLR, extracellular loop 1 (ECL1) of the CLR and its associated transmembrane (TM) regions, and finally extracellular loop 3 (ECL3) of the CLR and its juxtamembrane regions. A site-directed mutagenesis (SDM) strategy was used to investigate these regions, primarily substituting the innate residues of CLR with alanine and assessing the mutation on multiple criteria including a functional cAMP assay, cell-surface expression, total expression, agonist-mediated internalisation and aCGRP binding. The results are interpreted and discussed taking into consideration contemporary concepts surrounding Secretin-like GPCRs. Moreover, the thesis also contains details of RAMP purification. Overall the thesis provides novel data that furthers insight into the complex phenomenon of CGRP receptor activation. Site-directed mutants have been identified that affect aCGRP binding, receptor signal transduction, the CLR/RAMP1 interface and the integrity of the protein complex structure.

Novel peptide antagonists of adrenomedullin and calcitonin gene-related peptide receptors:identification, pharmacological characterization, and interactions with position 74 in receptor activity-modifying protein 1/3

Human adrenomedullin (AM) is a 52-amino acid peptide belonging to the calcitonin peptide family, which also includes calcitonin gene-related peptide (CGRP) and AM2. The two AM receptors, AM(1) and AM(2), are calcitonin receptor-like receptor (CL)/receptor activity-modifying protein (RAMP) (RAMP2 and RAMP3, respectively) heterodimers. CGRP receptors comprise CL/RAMP1. The only human AM receptor antagonist (AM(22-52)) is a truncated form of AM; it has low affinity and is only weakly selective for AM(1) over AM(2) receptors. To develop novel AM receptor antagonists, we explored the importance of different regions of AM in interactions with AM(1), AM(2), and CGRP receptors. AM(22-52) was the framework for generating further AM fragments (AM(26-52) and AM(30-52)), novel AM/alphaCGRP chimeras (C1-C5 and C9), and AM/AM(2) chimeras (C6-C8). cAMP assays were used to screen the antagonists at all receptors to determine their affinity and selectivity. Circular dichroism spectroscopy was used to investigate the secondary structures of AM and its related peptides. The data indicate that the structures of AM, AM2, and alphaCGRP differ from one another. Our chimeric approach enabled the identification of two nonselective high-affinity antagonists of AM(1), AM(2), and CGRP receptors (C2 and C6), one high-affinity antagonist of AM(2) receptors (C7), and a weak antagonist selective for the CGRP receptor (C5). By use of receptor mutagenesis, we also determined that the C-terminal nine amino acids of AM seem to be responsible for its interaction with Glu74 of RAMP3. We provide new information on the structure-activity relationship of AM, alphaCGRP, and AM2 and how AM interacts with CGRP and AM(2) receptors.

Similarity between class A and class B G-protein-coupled receptors exemplified through calcitonin gene-related peptide receptor modelling and mutagenesis studies

Modelling class B G-protein-coupled receptors (GPCRs) using class A GPCR structural templates is difficult due to lack of homology. The plant GPCR, GCR1, has homology to both class A and class B GPCRs. We have used this to generate a class A-class B alignment, and by incorporating maximum lagged correlation of entropy and hydrophobicity into a consensus score, we have been able to align receptor transmembrane regions. We have applied this analysis to generate active and inactive homology models of the class B calcitonin gene-related peptide (CGRP) receptor, and have supported it with site-directed mutagenesis data using 122 CGRP receptor residues and 144 published mutagenesis results on other class B GPCRs. The variation of sequence variability with structure, the analysis of polarity violations, the alignment of group-conserved residues and the mutagenesis results at 27 key positions were particularly informative in distinguishing between the proposed and plausible alternative alignments. Furthermore, we have been able to associate the key molecular features of the class B GPCR signalling machinery with their class A counterparts for the first time. These include the [K/R]KLH motif in intracellular loop 1, [I/L]xxxL and KxxK at the intracellular end of TM5 and TM6, the NPXXY/VAVLY motif on TM7 and small group-conserved residues in TM1, TM2, TM3 and TM7. The equivalent of the class A DRY motif is proposed to involve Arg(2.39), His(2.43) and Glu(3.46), which makes a polar lock with T(6.37). These alignments and models provide useful tools for understanding class B GPCR function.