Conformational analysis of the C-terminal Gly-Leu-Met-NH2 tripeptide of substance P bound to the NK-1 receptor

We examined the effect of simultaneously incorporating proline or proline-amino acid chimeras in positions 9, 10, and/or 11 of substance P, on the affinity for the two NK-1 binding sites and on second-messenger activation. Because these 3-substituted prolines constrain not only the (phi,psi) values of the peptide backbone, but also the chi space of the amino acid side chain, we were able to gather data on the structural requirements for high-affinity binding to the NK-1 receptor. We were able to confirm that this C-terminal component is crucial and that it should adopt an extended conformation close to a polyproline II structure when bound to the receptor. The partial additivity of these constraints, more specifically, for the NK-1M site, suggests that the peptide backbone flexibility around the hinge-point residue Gly9 is essential to subtly position crucial side chains.     

Molecular and structural determinants of adamantyl susceptibility to HLA-DRs allelic variants: an in silico approach to understand the mechanism of MLEs

Class II major histocompatibility complex (MHC II) molecules as expressed by antigen-presenting cells are heterodimeric cell-surface glycoprotein receptors that are fundamental in initiating and propagating an immune response by presenting tumor-associated antigenic peptides to CD4⁺/T_H cells. The loading efficiency of such peptides can be improved by small organic compounds (MHC Loading Enhancers—MLEs), that convert the non-receptive peptide conformation of MHC II to a peptide-receptive conformation. In a reversible reaction, these compounds open up the binding site of MHC II molecules by specific interactions with a yet undefined pocket. Here, we performed molecular docking and molecular dynamics simulation studies of adamantyl compounds on the predicted cavity around the P1 pocket of 2 allelic variants of HLA-DRs. The purpose was to investigate the suitability of adamantyl compounds as MLEs at the dimorphic β86 position. Docking studies revealed that besides numerous molecular interactions formed by the adamantyl compounds, Asnβ82, Tyrβ83, and Thrβ90 are the crucial amino acid residues that are characterized as the “sensors” of peptide loading. Molecular dynamics simulation studies exposed the dynamical structural changes that HLA-DRs adopted as a response to binding of 3-(1-adamantyl)-5-hydrazidocarbonyl-1H-pyrazole (AdCaPy). The conformations of AdCaPy complexed with the Glyβ86 HLA-DR allelic variant are well correlated with the stabilized form of peptide-loaded HLA-DRs, further confirming the role of AdCaPy as a MLE. Hydrogen bonding interaction analysis clearly demonstrated that after making suitable contacts with AdCaPy, HLA-DR changes its local conformation. However, AdCaPy complexed with HLA-DR having Valβ86 at the dimorphic position did not accommodate AdCaPy as MLE due to steric hindrance caused by the valine.     

Structural characterization of a paramagnetic metal-ion-assembled three-stranded alpha-helical coiled coil

A helical peptide designed to present an all-leucine core upon folding has been shown to exhibit concentration-dependent helicity and to exist as an ill-defined equilibrium population of oligomers. In marked contrast, an identical peptide covalently modified with a 2,2'-bipyridyl group at the N terminus forms a stable three-stranded parallel coiled coil in the presence of transition metal ions. We have employed paramagnetic Ni(2+) and Co(2+) ions to stabilize the trimeric assembly and to exploit their shift and relaxation properties in NMR structural studies. We find that metal-ion binding and helix-bundle folding are tightly coupled. Surprisingly, the three-helix bundle exhibits a dynamic N-terminal region, and a well-structured C-terminal half. The spectra indicate the presence of a dual conformation for the bundle extending from the N terminus to residue 12. The structure of the two isomeric forms has been ascertained from interpretation of NOEs in the Ni(II) complex and (1)H pseudocontact shifts in the Co(II) complex. Two different facial isomers with distinct susceptibility tensors were identified. The bulky leucine side chain at position 3 in the peptide chain appears to play a role in the conformational variation at the N terminus.     

Conformational studies of immunodominant myelin basic protein 1�C11 analogues using NMR and molecular modeling

Τwo dimensional nuclear magnetic resonance studies complimented by molecular dynamics simulations were conducted to investigate the conformation of the immunodominant epitope of acetylated myelin basic protein residues 1-11 (Ac-MBP(1-11)) and its altered peptide ligands, mutated at position 4 to an alanine (Ac-MBP(1-11)[4A]) or a tyrosine residue (Ac-MBP(1-11)[4Y]). Conformational analysis of the three analogues indicated that they adopt an extended conformation in DMSO solution as no long distance NOE connectivities were observed and seem to have a similar conformation when bound to the active site of the major histocompatibility complex (MHC II). The interaction of each peptide with MHC class II I-A(u) was further investigated in order to explore the molecular mechanism of experimental autoimmune encephalomyelitis induction/inhibition in mice. The present findings indicate that the Gln(3) residue, which serves as a T-cell receptor (TCR) contact site in the TCR/peptide/I-A(u) complex, has a different orientation in the mutated analogues especially in the Ac-MBP(1-11)[4A] peptide. In particular the side chain of Gln(3) is not solvent exposed as for the native Ac-MBP(1-11) and it is not available for interaction with the TCR.     

Revisiting the structural basis and energetic landscape of susceptibility difference between HLA isotypes to allergic rhinitis

The human leukocyte antigen class II (HLA II) molecules are implicated in the immunopathogenesis of allergic rhinitis (AR). The HLA II contains three allelic isotypes HLA-DR, -DQ, and -QP that exhibit considerably different susceptibility to AR. Here, we investigated the structural basis and energetic landscape of the susceptibility difference between the three HLA II isotypes to AR by combining computational analysis and experimental assay. Multiple sequence alignment revealed a low conservation among the three subtypes with sequence identity of ∼10% between them, suggesting that the peptide repertoires presented by HLA-DR, -DP and -DQ are not overlapped to each other, and they may be involved in different immune functions and dysfunctions. Structural analysis imparted that the antigenic peptides are rooted on the peptide-binding groove of HLA molecules and hold in a PPII-like helical conformation. Subsequently, the interaction behavior of 17 AR allergen-derived peptides with HLA-DR, −DP and −DQ was investigated using a statistics-based quantitative structure-activity relationship (QSAR) predictor. It was found a significant difference between the binding capabilities of these antigenic peptides to HLA-DR and to HLA-DP/-DQ; the former showed a generally higher affinity than the latter with p-value of 0.02 obtained from 2-tailed Student’s t-test. The computational findings were then confirmed by HLA II–peptide stability assay, which demonstrated that the AR allergen-derived peptides have a high in vitro selectivity for HLA-DR over HLA-DP/-DQ. Thus, the HLA-DR isotype, rather than HLA-DP and −DQ, is expected to associate with the pathological process of AR.     

Conformation study of HA(306-318) antigenic peptide of the haemagglutinin influenza virus protein

Several HLA-DR alleles present the immunodominant HA(306-318) peptide of haemagglutinin of the influenza virus to T cells. NMR data of the peptide in various water solutions exclude any alpha-helix or turn conformations. Circular dichroism and Fourier transform infrared spectroscopies indicate an estimated beta-extended structure in water of 31% and 28%, respectively, with spectra shape similar to the ones observed for beta-sheet containing proteins. The H/D amide exchange suggests a stable length-dependent interchain hydrogen-bonding. The partially beta-extended conformation of HA(306-318) in solution might be close to the one found in HA(306-318)-HLA-DR1 complex. These results suggest different interconverting extended conformations of HA(306-318), depending on the microenvironment of the solution medium. This flexibility emphasizes the ability of some peptides to fit more easily the binding site of several HLA-DR molecules. Similar results were obtained on the HIV P25(263-277) peptide which has been previously shown to be a good DR1 binder. From a vibrational point of view, infrared Amide I frequencies of secondary structures in peptides were ascertained. As previously demonstrated for proteins in solution, Fourier transform infrared and circular dichroism spectroscopies appear to be valuable tools for conformational properties of peptides. Their use may contribute to the detection of peptide conformation-binding relationship which has to be further tested by biochemical and biological studies.     

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.     

Solution and structural characterization of iron(II) complexes with ortho-halogenated phenolates: insights into potential substrate binding modes in hydroquinone dioxygenases

The new ligand cis,cis-1,3,5-tris-(E)-(tolylideneimino)cyclohexane (TACH-o-tolyl) forms a 1:1 complex with iron(II). Addition of substituted phenolates forms 1:1:1 ligand:iron:phenolate complexes, which have been characterized both in the solid state and in solution. There is complete binding of the phenolate to the complex only when there are ortho-halogens on the phenolate. The tertiary complexes with ortho-halo-substituted phenolates exhibit short Fe-halogen distances, and the complex containing a non-coordinating but similarly sized ortho-methyl phenolate has a significantly different conformation and coordination geometry. Therefore, it is likely that the metal-halogen interaction stabilizes the complexes. The iron(II)-halogen interaction in these complexes may explain the substrate specificity of PcpA and LinE, enzymes that preferentially bind phenols and hydroquinones containing halogen substituents in ortho positions.     

Designed high affinity Cu2+-binding alpha-helical foldamer

The design, synthesis, and characterization of a folded high-affinity metal-binding peptide is described. Based on the previously described folded peptide NTH-18, in which an alpha-helix was constrained through two disulfide bonds to a C-terminal extension of noncanonical secondary structure, a peptide (1) was designed to contain two histidine residues in positions 3 and 7. Air oxidation of 1 led to the formation of peptide 2, which contained two intramolecular disulfide bonds. The presence of the two histidines significantly destabilized the alpha-helical structure of 2 when compared to NTH-18. However, CD spectroscopy revealed that the addition of certain transition metal ions allowed the reformation of a stable alpha-helix. CD, NMR, and EPR spectroscopy as well as MALDI-TOF mass spectrometry indicated that 2 bound to Cu2+ to form a 1:1 complex via the imidazoles of the two histidine side chains. A glycine displacement assay revealed a dissociation constant for this complex of 5 nM at pH 8, which is the lowest reported value for a designed Cu2+-binding peptide. This peptide displayed more than 100-fold selectivity for Cu2+ over Zn2+, Ni2+, and Co2+. The 1.05 A crystal structure of the Cu(II)-complex of 2 revealed a square-pyramidal coordination geometry and confirmed that 2 bound to copper in an alpha-helical conformation via its two histidine side chains. The high affinity metal binding of peptide 2 demonstrates that metals can be used for the selective nucleation of alpha-helices.     

Non-radioactive detection of MHC class II-peptide antigen complexes in the sub-picomole range by high-performance size-exclusion chromatography with fluorescence detection

In order to avoid chemical or structural modification of T-cell epitopes by labelling, a high-performance size-exclusion chromatographic fluorescence binding assay was developed, based on the intrinsic Trp fluorescence of major histocompatibility complex (MHC) proteins. The increase in Trp fluorescence intensity of the isolated human MHC product HLA-DR 1 on complex formation with unlabelled influenza matrix peptide[18-29] (IM[18-29]) was examined. Binding of IM[18-29] to the heterodimeric form of HLA-DR 1 (Kd = 4.8 mM) and to the disassembled alpha-and beta-subunits (Kd = 9.2 mM) could be demonstrated. In addition, the assay showed the peptide-induced formation of a dimeric conformer of HLA-DR 1, the nature of which is still undefined. Detection of HLA-DR 1 subunit-peptide complexes was possible in amounts of 25 ng in 10 microliter (80 fmol/microliter). The technique proved to be reproducible and less time consuming than common methods that need fluorescence or radioactive labelling.     

Site-selective metal binding by designed alpha-helical peptides

It is known that the designed alpha-helical peptide family TRI [(Ac-G(LKALEEK)4G-CONH2)], containing single site substitution of a cysteine for a leucine, is capable of binding Cd(II) within a three-stranded coiled coil. The binding affinity of cadmium is dependent upon the site of substitution, with cysteine incorporated at the a site leading to cadmium complexes of higher affinity than when a d site was modified. In this work we have examined whether this differential binding affinity can be expressed in a di-cysteine-substituted peptide in order to develop site specificity within a designed system. The peptide TRI L9CL19C was used to determine whether significant differences in binding affinities at nearly proximal sites could be achieved in a short designed peptide. On the basis of 113Cd, 1H NMR, and circular dichroic spectroscopies, we have shown that 1 equiv of Cd(II) binds exclusively at the a site. Only after that position is filled does the second site become populated. Thus, the TRI system represents the first example where stoichiometrically equivalent peptides with different sequences form the framework for designing molecular assemblies that show site-specific ion recognition. We propose that the distinct metal affinities are due to the cysteine conformers at different substitution points along the peptide. Furthermore, we have shown that site selectivity in biomolecules can be encoded into relatively short peptides with helical sequences and, therefore, do not necessarily require the extensive protein scaffolds found in natural systems.     

Tuning the Binding Affinity and Selectivity of Perfluoroaryl-Stapled Peptides by Cysteine-Editing

A growing number of approaches to “staple” α-helical peptides into a bioactive conformation using cysteine cross-linking are emerging. Here, the replacement of l -cysteine with “cysteine analogues” in combinations of different stereochemistry, side chain length and beta-carbon substitution, is explored to examine the influence that the thiol-containing residue(s) has on target protein binding affinity in a well-explored model system, p53–MDM2/MDMX, which is constituted by the interaction of the tumour suppressor protein p53 and proteins MDM2 and MDMX, which regulate p53 activity. In some cases, replacement of one or more l -cysteine residues afforded significant changes in the measured binding affinity and target selectivity of the peptide. Computationally constructed homology models indicate that some modifications, such as incorporating two d -cysteine residues, favourably alter the positions of key functional amino acid side chains, which is likely to cause changes in binding affinity, in agreement with measured surface plasmon resonance data.     

Side-chain conformational space analysis (SCSA): A multi conformation-based QSAR approach for modeling and prediction of protein–peptide binding affinities

In this article, the concept of multi conformation-based quantitative structure-activity relationship (MCB-QSAR) is proposed, and based upon that, we describe a new approach called the side-chain conformational space analysis (SCSA) to model and predict protein-peptide binding affinities. In SCSA, multi-conformations (rather than traditional single-conformation) have received much attention, and the statistical average information on multi-conformations of side chains is determined using self-consistent mean field theory based upon side chain rotamer library. Thereby, enthalpy contributions (including electrostatic, steric, hydrophobic interaction and hydrogen bond) and conformational entropy effects to the binding are investigated in terms of occurrence probability of residue rotamers. Then, SCSA was applied into the dataset of 419 HLA-A 0201 binding peptides, and nonbonding contributions of each position in peptide ligands are well determined. For the peptides, the hydrogen bond and electrostatic interactions of the two ends are essential to the binding specificity, van der Waals and hydrophobic interactions of all the positions ensure strong binding affinity, and the loss of conformational entropy at anchor positions partially counteracts other favorable nonbonding effects.     

Simultaneous prediction of binding free energy and specificity for PDZ domain–peptide interactions

Interactions between protein domains and linear peptides underlie many biological processes. Among these interactions, the recognition of C-terminal peptides by PDZ domains is one of the most ubiquitous. In this work, we present a mathematical model for PDZ domain–peptide interactions capable of predicting both affinity and specificity of binding based on X-ray crystal structures and comparative modeling with Rosetta. We developed our mathematical model using a large phage display dataset describing binding specificity for a wild type PDZ domain and 91 single mutants, as well as binding affinity data for a wild type PDZ domain binding to 28 different peptides. Structural refinement was carried out through several Rosetta protocols, the most accurate of which included flexible peptide docking and several iterations of side chain repacking and backbone minimization. Our findings emphasize the importance of backbone flexibility and the energetic contributions of side chain-side chain hydrogen bonds in accurately predicting interactions. We also determined that predicting PDZ domain–peptide interactions became increasingly challenging as the length of the peptide increased in the N-terminal direction. In the training dataset, predicted binding energies correlated with those derived through calorimetry and specificity switches introduced through single mutations at interface positions were recapitulated. In independent tests, our best performing protocol was capable of predicting dissociation constants well within one order of magnitude of the experimental values and specificity profiles at the level of accuracy of previous studies. To our knowledge, this approach represents the first integrated protocol for predicting both affinity and specificity for PDZ domain–peptide interactions.     

Peptide tic-tac-toe: heterotrimeric coiled-coil specificity from steric matching of multiple hydrophobic side chains

Specific coiled-coil heterotrimers result from steric matching of hydrophobic core side chains. A 2:1 heterotrimer is formed by peptides containing alanine or cyclohexylalanine, respectively, at a central core residue. Detailed thermodynamic analysis reveals that the designed complex is considerably more stable than the corresponding alanine homotrimer (deltaT(m) = 25 degrees C, deltadeltaG(unf) = 4.5 kcal/mol), while control complexes with naphthylalanine or cyclopropylalanine peptides are much less stable. However, the cyclohexylalanine homotrimer is of comparable stability to the 2:1 complex, prompting an investigation of multiply substituted peptides. A specific 1:1:1 heterotrimer is formed from three independent peptide strands, each bearing one large (cyclohexylalanine) and two small (alanine) side chains at the same three core positions but in different order. The combined impact of three substitutions improves specificity to the point where each pure peptide and all pairwise equimolar mixtures form significantly less stable complexes (deltaTm = 22-24 degrees C). The capacity for specific complex formation governed by multiple unnatural core side chains should facilitate design of numerous new peptide assemblies.     

Controlling the DNA binding specificity of bHLH proteins through intramolecular interactions

Reversible control of the conformation of proteins was employed to probe the relationship between flexibility and specificity of the basic helix-loop-helix protein MyoD. A fusion protein (apaMyoD) was designed where the basic DNA binding helix of MyoD was stablized by an amino-terminal extension with a sequence derived from the bee venom peptide apamin. The disulfide-stabilized helix from apamin served as a nucleus for a helix that extended for a further ten residues, thereby holding apaMyoD's DNA recognition helix in a predominantly alpha-helical conformation. The thermal stability of the DNA complexes of apaMyoD was increased by 13 degrees C relative to MyoD-bHLH. Measurements of the fluorescence anisotropy change on DNA binding indicated that apaMyoD bound to E-box-containing DNA sequences with enhanced affinity relative to MyoD-bHLH. Consequently, the DNA binding specificity of apaMyoD was increased 10-fold.     

Design of Peptides with α,β-Dehydro Residues: Synthesis,Crystal Structure and Molecular Conformation of a Peptide N-Boc-Phe-ΔPhe-Ile-OCH<Subscript>3</Subscript>

To develop a complete set of design rules with α,β-dehydro residues, a tripeptide N-Boc-Phe-ΔPhe-Ile-OCH₃ was synthesized. The synthesis was carried out in solution phase using azlactone procedure. The three-dimensional structure of the peptide was determined by X-ray diffraction method and refined to an R-factor of 0.085. The structure contains three peptide molecules in the asymmetric unit. In all the three crystallographically independent molecules ΔPhe residue adopts one of the three conformations that have been reported for a ΔPhe residue. The overall conformations of three peptide molecules in the asymmetric unit are not similar. Two out of three crystallographically independent molecules adopt type II β-turn conformations whereas the third molecule is found having the characteristic S-shaped conformation in which the values of dihedral angles φ, ψ have opposite signs alternately. One of these two types of conformations has been observed when a ΔPhe is introduced at (i+2) position of a tetrapeptide. The β-turn conformation is stabilized by a 4→1 hydrogen bond where the hydrophobic side chains of residues at (i+1) and (i+3) positions stabilized the unfolded conformation with van der Waals interactions. The three independent molecules are locked together by three hydrogen bonds between molecules A and B and two hydrogen bonds between molecules B and C.     

The Impact of Amino Acid Side Chain Mutations in Conformational Design of Peptides and Proteins

Local energetic effects of amino acid replacements are often considered to have only a moderate influence on the backbone conformation of proteins or peptides. As these effects are difficult to determine experimentally, no comparison has yet been performed. However, knowledge of the influence of side chain mutations is essential in protein homology modeling and in optimizing biologically active peptide ligands in medicinal chemistry. Furthermore, the tool of N‐methylation of peptides is of increasing importance for the design of peptidic drugs to gain oral availability or receptor selectivity. However, N‐methylation is often accompanied by considerable population of cis‐peptide bond structures, resulting in completely different conformations compared with the parent peptide. To retain a favored structure, it might be important to understand the effect of different side chains on the backbone conformation and to enable the introduction of an N‐methylation at the right position without disturbing a biologically active conformation. In order to detect even small energetic effects due to side chain mutations, we employed a trick to investigate the structural equilibrium of a selected cyclic pentapeptide in which two conformations are equally populated. Very small energetic differences between both conformations could easily be determined experimentally by identifying shifts in the population of both isomers.     

Oxime Side‐Chain Cross‐Links in an α‐Helical Coiled‐Coil Protein: Structure,Thermodynamics, and Folding‐Templated Synthesis of Bicyclic Species

Covalent side‐chain cross‐links are a versatile method to control peptide folding, particularly when α‐helical secondary structure is the target. Here, we examine the application of oxime bridges, formed by the chemoselective reaction between aminooxy and aldehyde side chains, for the stabilization of a helical peptide involved in a protein–protein complex. A series of sequence variants of the dimeric coiled coil GCN4‐p1 bearing oxime bridges at solvent‐exposed positions were prepared and biophysically characterized. Triggered unmasking of a side‐chain aldehyde in situ and subsequent cyclization proceed rapidly and cleanly at pH 7 in the folded protein complex. Comparison of folding thermodynamics among a series of different oxime bridges show that the cross links are consistently stabilizing to the coiled coil, with the extent of stabilization sensitive to the exact size and structure of the macrocycle. X‐ray crystallographic analysis of a coiled coil with the best cross link in place and a second structure of its linear precursor show how the bridge is accommodated into an α‐helix. Preparation of a bicyclic oligomer by simultaneous formation of two linkages in situ demonstrates the potential use of triggered oxime formation to both trap and stabilize a particular peptide folded conformation in the bound state.     

Substrate specificity and novel selective inhibitors of TNF-alpha converting enzyme (TACE) from two-dimensional substrate mapping

We report a systematic analysis of the P1' and P2' substrate specificity of TNF-alpha converting enzyme (TACE) using a peptide library and a novel analytical method, and we use the substrate specificity information to design novel reverse hydroxamate inhibitors. Initial truncation studies, using the amino acid sequence around the cleavage site in precursor-TNF-alpha, showed that good turnover was obtained with the peptide DNP-LAQAVRSS-NH2. Based on this result, 1000 different peptide substrates of the form Biotin-LAQA-P1'-P2'-SSK(DNP)-NH2 were prepared, with 50 different natural and unnatural amino acids at P1' in combination with 20 different amino acids at P2'. The peptides were pooled, treated with purified microsomal TACE, and the reaction mixtures were passed over a streptavidin affinity column to remove unreacted substrate and the N-terminal biotinylated product. C-terminal cleavage products not binding to streptavidin were subjected to liquid chromatography/mass spectrometry analysis where individual products were identified and semiquantitated. 25 of the substrates were resynthesized as discrete peptides and assayed with recombinant TACE. The experiments show that recombinant TACE prefers lipophilic amino acids at the P1' position, such as phenylglycine, homophenylalanine, leucine and valine. At the P2' position, TACE can accommodate basic amino acids, such as arginine and lysine, as well as certain non-basic amino acids such as citrulline, methionine sulfoxide and threonine. These substrate preferences were used in the design of novel reverse hydroxamate TACE inhibitors with phenethyl and 5-methyl-thiophene-methyl side-chains at P1', and threonine and nitro-arginine at P2'.