Computer determination of peptide conformations in water: different roads to structure.

Fragments of proteins (short peptides) that "fold" suggest a mechanism of how complete conformational search in protein folding is avoided. We used a computational method to determine structures of two foldable peptides in explicit water: RVEW and CSVTC. The optimization starts from random structures and no experimental constraints are used. In agreement with NMR data, the simulations find a hydrophobic pair (Val/Trp) in REVW. The structure of CSVTC is induced by a surface water that bridges two amide hydrogens, a drive to structure hypothesized by Ben-Naim [Ben-Naim, A. (1990) J. Chem. Phys. 93, 8196-8210] that is largely ignored in studies of folding. Tendency to structure in short peptide chains suggests a mechanism for the formation of short-range nucleation sites in protein folding.

Dynamics and Conformational Studies of TOAC Spin Labeled Analogues of Ctx(Ile21)-Ha Peptide from Hypsiboas albopunctatus

Antimicrobial peptides (AMPs) isolated from several organisms have been receiving much attention due to some specific features that allow them to interact with, bind to, and disrupt cell membranes. The aim of this paper was to study the interactions between a membrane mimetic and the cationic AMP Ctx(Ile21)-Ha as well as analogues containing the paramagnetic amino acid 2,2,6,6-tetramethylpiperidine-1-oxyl-4-amino-4-carboxylic acid (TOAC) incorporated at residue positions n = 0, 2, and 13. Circular dichroism studies showed that the peptides, except for [TOAC13]Ctx(Ile21)-Ha, are unstructured in aqueous solution but acquire different amounts of α-helical secondary structure in the presence of trifluorethanol and lysophosphocholine micelles. Fluorescence experiments indicated that all peptides were able to interact with LPC micelles. In addition, Ctx(Ile21)-Ha and [TOAC13]Ctx(Ile21)-Ha peptides presented similar water accessibility for the Trp residue located near the N-terminal sequence. Electron spin resonance experiments showed two spectral components for [TOAC0]Ctx(Ile21)-Ha, which are most likely due to two membrane-bound peptide conformations. In contrast, TOAC2 and TOAC13 derivatives presented a single spectral component corresponding to a strong immobilization of the probe. Thus, our findings allowed the description of the peptide topology in the membrane mimetic, where the N-terminal region is in dynamic equilibrium between an ordered, membrane-bound conformation and a disordered, mobile conformation; position 2 is most likely situated in the lipid polar head group region, and residue 13 is fully inserted into the hydrophobic core of the membrane. © 2013 Vicente et al.

Alpha-helical, but not beta-sheet, propensity of proline is determined by peptide environment.

Proline is established as a potent breaker of both alpha-helical and beta-sheet structures in soluble (globular) proteins. Thus, the frequent occurrence of the Pro residue in the putative transmembrane helices of integral membrane proteins, particularly transport proteins, presents a structural dilemma. We propose that this phenomenon results from the fact that the structural propensity of a given amino acid may be altered to conform to changes imposed by molecular environment. To test this hypothesis on proline, we synthesized model peptides of generic sequence H2N-(Ser-LyS)2-Ala- Leu-Z-Ala-Leu-Z-Trp-Ala-Leu-Z-(Lys-Ser)3-OH (Z = Ala and/or Pro). Peptide conformations were analyzed by circular dichroism spectroscopy in aqueous buffer, SDS, lysophosphatidylglycerol micelles, and organic solvents (methanol, trifluoroethanol, and 2-propanol). The helical propensity of Pro was found to be greatly enhanced in the membrane-mimetic environments of both lipid micelles and organic solvents. Proline was found to stabilize the alpha-helical conformation relative to Ala at elevated temperatures in 2-propanol, an observation that argues against the doctrine that Pro is the most potent alpha-helix breaker as established in aqueous media. Parallel studies in deoxycholate micelles of the temperature-induced conformational transitions of the single-spanning membrane bacteriophage IKe major coat protein, in which the Pro-containing wild type was compared with Pro30 --> Ala mutant, Pro was found to protect the helix, but disrupt the beta-sheet structure as effectively as it does to model peptides in water. The intrinsic capacity of Pro to disrupt beta-sheets was further reflected in a survey of porins where Pro was found to be selectively excluded from the core of membrane-spanning beta-sheet barrels. The overall data provide a rationale for predicting and understanding the structural consequences when Pro occurs in the context of a membrane.

The solution structure of a copper(II) compound of a new cyclic octapeptide by EPR spectroscopy and force field calculations

A new cyclic octapeptide, cyclo(Ile-Ser-(Gly)Thz-Ile-Thr-(Gly)Thz) (PatN), related to patellamide A, has been synthesized and reacted with copper(II) and base to form mono- and dinuclear complexes. The coordination environments around copper(TI) have been characterized by EPR spectroscopy. The solution structure of the thermodynamically most stable product, a purple dicopper(TI) compound, has been examined by simulating weakly dipole-dipole coupled EPR spectra based upon structural parameters obtained from force field (MM and MD) calculations. The MM-EPR method produces a saddle-shaped structure for [Cu-2(PatN)(OH2)(6)] that is similar to the known solution structure of patellamide A and the known solid-state structure of [Cu-2(AscidH(2))CO3(OH2)(2)]. Compared with the latter, [Cu-2(PatN)] has no carbonate bridge and a significantly flatter topology. The MM-EPR approach to solution-structure determination for paramagnetic metallopeptides may find wide applications to other metallopeptides and metalloproteins.

Conformational selection of inhibitors and substrates by proteolytic enzymes: Implications for drug design and polypeptide processing

Inhibitors of proteolytic enzymes (proteases) are emerging as prospective treatments for diseases such as AIDS and viral infections, cancers, inflammatory disorders, and Alzheimer's disease. Generic approaches to the design of protease inhibitors are limited by the unpredictability of interactions between, and structural changes to, inhibitor and protease during binding. A computer analysis of superimposed crystal structures for 266 small molecule inhibitors bound to 48 proteases (16 aspartic, 17 serine, 8 cysteine, and 7 metallo) provides the first conclusive proof that inhibitors, including substrate analogues, commonly bind in an extended beta-strand conformation at the active sites of all these proteases. Representative superimposed structures are shown for (a) multiple inhibitors bound to a protease of each class, (b) single inhibitors each bound to multiple proteases, and (c) conformationally constrained inhibitors bound to proteases. Thus inhibitor/substrate conformation, rather than sequence/composition alone, influences protease recognition, and this has profound implications for inhibitor design. This conclusion is supported by NMR, CD, and binding studies for HIV-1 protease inhibitors/ substrates which, when preorganized in an extended conformation, have significantly higher protease affinity. Recognition is dependent upon conformational equilibria since helical and turn peptide conformations are not processed by proteases. Conformational selection explains the resistance of folded/structured regions of proteins to proteolytic degradation, the susceptibility of denatured proteins to processing, and the higher affinity of conformationally constrained 'extended' inhibitors/substrates for proteases. Other approaches to extended inhibitor conformations should similarly lead to high-affinity binding to a protease.

Crystallographic analysis of endogenous peptides associated with HLA-DR1 suggests a common, polyproline II-like conformation for bound peptides.

The structure of the human major histocompatibility complex (MHC) class II molecule HLA-DR1 derived from the human lymphoblastoid cell line LG-2 has been determined in a complex with the Staphylococcus aureus enterotoxin B superantigen. The HLA-DR1 molecule contains a mixture of endogenous peptides derived from cellular or serum proteins bound in the antigen-binding site, which copurify with the class II molecule. Continuous electron density for 13 amino acid residues is observed in the MHC peptide-binding site, suggesting that this is the core length of peptide that forms common interactions with the MHC molecule. Electron density is also observed for side chains of the endogenous peptides. The electron density corresponding to peptide side chains that interact with the DR1-binding site is more clearly defined than the electron density that extends out of the binding site. The regions of the endogenous peptides that interact with DRI are therefore either more restricted in conformation or sequence than the peptide side chains or amino acids that project out of the peptide-binding site. The hydrogen-bond interactions and conformation of a peptide model built into the electron density are similar to other HLA-DR-peptide structures. The bound peptides assume a regular conformation that is similar to a polyproline type II helix. The side-chain pockets and conserved asparagine residues of the DR1 molecule are well-positioned to interact with peptides in the polyproline type II conformation and may restrict the range of acceptable peptide conformations.

Studies of the interaction of Potassium(I), Calcium(II), Magnesium(II), and Copper(II) with Cyclosporin A

Metal ion binding properties of the immunosuppressant drug cyclosporin A have been investigated. Complexation studies in acetonitrile solution using H-1 NMR and CD spectroscopy yielded 1:1 metal-peptide binding constants (log(10)K) for potassium(l), < 1, magnesium(II), 4.8 +/- 0.2. and calcium(II), 5.0 +/- 1.0. The interaction of copper(II) with cyclosporin A in methanol was investigated with UV/visible and electron paramagnetic resonance (EPR) spectroscopy. No complexation of copper(II) was observed in neutral solution. In the presence of base, monomeric copper(II) complexes were detected. These results support the possibility that cyclosporin A has ionophoric properties for biologically important essential metal ions. (C) 2003 Elsevier Inc. All rights reserved.

Conformations of platypus venom C-type natriuretic peptide in aqueous solution and sodium dodecyl sulfate micelles

Nuclear magnetic resonance spectroscopy was used to investigate the conformations of the platypus venom C-type natriuretic peptide A (OvCNPa) in aqueous solutions and in solutions containing sodium dodecyl sulfate (SDS) micelles. The chemically synthesized OvCNPa showed a substantial decrease in flexibility in aqueous solution at 10 degreesC, allowing the observation of medium- and long-range nuclear Overhauser enhancement (NOE) connectivities. Three-dimensional structures calculated using these data showed flexible and reasonably well-defined regions, the locations of which were similar in the two solvents. In aqueous solution, the linear part that spans residues 3-14 was basically an extended conformation while the cyclic portion, defined by residues 23-39, contained a series of beta-turns. The overall shape of the cyclic portion was similar to that observed for an atrial natriuretic peptide (ANP) variant in aqueous solution. OvCNPa adopted a different conformation in SDS micelles wherein the N-terminal region, defined by residues 2-10, was more compact, characterised by turns and a helix, while the cyclic region had turns and an overall shape that was fundamentally different from those structures observed in aqueous solution. The hydrophobic cluster, situated at the centre of the ring of the structure in aqueous solution, was absent in the structure in the presence of SDS micelles. Thus, OvCNPa interacts with SDS micelles and can possibly form ion-channels in cell membranes. (C) 2002 Elsevier Science Ltd. All rights reserved.

Peptide transporter DtpA has two alternate conformations, one of which is promoted by inhibitor binding

Peptide transporters (PTRs) of the large PTR family facilitate the uptake of di- and tripeptides to provide cells with amino acids for protein synthesis and for metabolic intermediates. Although several PTRs have been structurally and functionally characterized, how drugs modulate peptide transport remains unclear. To obtain insight into this mechanism, we characterize inhibitor binding to the Escherichia coli PTR dipeptide and tripeptide permease A (DtpA), which shows substrate specificities similar to its human homolog hPEPT1. After demonstrating that Lys[Z-NO2]-Val, the strongest inhibitor of hPEPT1, also acts as a high-affinity inhibitor for DtpA, we used single-molecule force spectroscopy to localize the structural segments stabilizing the peptide transporter and investigated which of these structural segments change stability upon inhibitor binding. This characterization was done with DtpA embedded in the lipid membrane and exposed to physiologically relevant conditions. In the unbound state, DtpA adopts two main alternate conformations in which transmembrane α-helix (TMH) 2 is either stabilized (in ∼43% of DtpA molecules) or not (in ∼57% of DtpA molecules). The two conformations are understood to represent the inward- and outward-facing conformational states of the transporter. With increasing inhibitor concentration, the conformation characterized by a stabilized TMH 2 becomes increasingly prevalent, reaching ∼92% at saturation. Our measurements further suggest that Lys[Z-NO2]-Val interacts with discrete residues in TMH 2 that are important for ligand binding and substrate affinity. These interactions in turn stabilize TMH 2, thereby promoting the inhibited conformation of DtpA.

Metal Clips Induce Folding of a Short Unstructured Peptide into an a-Helix via Turn Conformations in Water. Kinetic versus Thermodynamic Products

Short peptides corresponding to two to four a-helical turns of proteins are not thermodynamically stable helices in water. Unstructured octapeptide Ac-His1*-Ala2-Ala3-His4*-His5*-Glu6-Leu7-His8*-NH2 (1) reacts with two [Pd ((NH2)-N-15(CH2)(2) (NH2)-N-15)(NO3)(2)] in water to form a kinetically stable intermediate, [{Pden}(2)-{(1,4)(5,8)-peptide}](2), in which two 19-membered metallocyclic rings stabilize two peptide turns. Slow subsequent folding to a thermodynamically more stable two-turn a-helix drives the equilibrium to [{Pden}(2)-{(1,5)(4,8)-peptide}] (3), featuring two 22-membered rings. This transformation from unstructured peptide via turns to an a-helix suggests that metal clips might be useful probes for investigating peptide folding.

Conformations of cyclic octapeptides and the influence of heterocyclic ring constraints upon calcium binding

A comparison is made between the structures and calcium binding properties of four cyclic octapeptides that differ in the number of heterocyclic thiazole and oxazoline ring constraints. The conformations of the naturally occurring cyclic octapeptides ascidiacyclamide 1 and patellamide D 2, which each contain two oxazoline and two thiazole rings, are compared by H-1 NMR spectroscopy with the analogues cyclo(Thr-D-Val(Thz)-Ile)(2) 3 with just two thiazoles, and cyclo(Thr-D-Val-alpha Abu-Ile)(2) 4, with no 5-membered rings. The conformations observed in the solid state for ascidiacyclamide (saddle) and patellamide D (twisted figure of eight) were retained in solution, whilst peptide 3 was found to have a chair shape and peptide 4 displayed a range of conformations. The solid state structure of 4 revealed that the peptide takes a relatively planar conformation with a number of transannular hydrogen bonds, which are apparently retained in solution. Complexation studies utilising H-1 NMR and CD spectroscopy yielded 1∶1 calcium-peptide binding constants (log K) for the four peptides (2.9 (1), 2.8 (2), 4.0 (3) and 5.5 (4)) as well as a 1 : 2 metal-peptide binding constant for 3 (log K = 4.5). The affinity for Ca2+ thus decreases with increasing number of 5-membered ring constraints in the macrocycle (4 > 3 > 2 approximate to 1).

Identification of initiation sites for T4 lysozyme folding using CD and NMR spectroscopy of peptide fragments

Using CD and 2D H-1 NMR spectroscopy, we have identified potential initiation sites for the folding of T4 lysozyme by examining the conformational preferences of peptide fragments corresponding to regions of secondary structure. CD spectropolarimetry showed most peptides were unstructured in water, but adopted partial helical conformations in TFE and SDS solution. This was also consistent with the H-1 NMR data which showed that the peptides were predominantly disordered in water, although in some cases, nascent or small populations of partially folded conformations could be detected. NOE patterns, coupling constants, and deviations from random coil Her chemical shift values complemented the CD data and confirmed that many of the peptides were helical in TFE and SDS micelles. In particular, the peptide corresponding to helix E in the native enzyme formed a well-defined helix in both TFE and SDS, indicating that helix E potentially forms an initiation site for T4 lysozyme folding. The data for the other peptides indicated that helices D, F, G, and H are dependent on tertiary interactions for their folding and/or stability. Overall, the results from this study, and those of our earlier studies, are in agreement with modeling and IID-deuterium exchange experiments, and support an hierarchical model of folding for T4 lysozyme.

Synthesis of new peptide derivatives that self-assemble into nanostructured hydrogels for biomedical applications

Tese de Doutoramento em Ciências (área de especialização em Química)

Modeling of the TCR-MHC-peptide complex.

The methodology for generating a homology model of the T1 TCR-PbCS-K(d) class I major histocompatibility complex (MHC) class I complex is presented. The resulting model provides a qualitative explanation of the effect of over 50 different mutations in the region of the complementarity determining region (CDR) loops of the T cell receptor (TCR), the peptide and the MHC's alpha(1)/alpha(2) helices. The peptide is modified by an azido benzoic acid photoreactive group, which is part of the epitope recognized by the TCR. The construction of the model makes use of closely related homologs (the A6 TCR-Tax-HLA A2 complex, the 2C TCR, the 14.3.d TCR Vbeta chain, the 1934.4 TCR Valpha chain, and the H-2 K(b)-ovalbumine peptide), ab initio sampling of CDR loops conformations and experimental data to select from the set of possibilities. The model shows a complex arrangement of the CDR3alpha, CDR1beta, CDR2beta and CDR3beta loops that leads to the highly specific recognition of the photoreactive group. The protocol can be applied systematically to a series of related sequences, permitting the analysis at the structural level of the large TCR repertoire specific for a given peptide-MHC complex.

Molecular modeling of peptide-MHC class I complexes : applications to peptide based immunotherapy

Summary The specific CD8+ T cell immune response against tumors relies on the recognition by the T cell receptor (TCR) on cytotoxic T lymphocytes (CTL) of antigenic peptides bound to the class I major histocompatibility complex (MHC) molecule. Such tumor associated antigenic peptides are the focus of tumor immunotherapy with peptide vaccines. The strategy for obtaining an improved immune response often involves the design of modified tumor associated antigenic peptides. Such modifications aim at creating higher affinity and/or degradation resistant peptides and require precise structures of the peptide-MHC class I complex. In addition, the modified peptide must be cross-recognized by CTLs specific for the parental peptide, i.e. preserve the structure of the epitope. Detailed structural information on the modified peptide in complex with MHC is necessary for such predictions. In this thesis, the main focus is the development of theoretical in silico methods for prediction of both structure and cross-reactivity of peptide-MHC class I complexes. Applications of these methods in the context of immunotherapy are also presented. First, a theoretical method for structure prediction of peptide-MHC class I complexes is developed and validated. The approach is based on a molecular dynamics protocol to sample the conformational space of the peptide in its MHC environment. The sampled conformers are evaluated using conformational free energy calculations. The method, which is evaluated for its ability to reproduce 41 X-ray crystallographic structures of different peptide-MHC class I complexes, shows an overall prediction success of 83%. Importantly, in the clinically highly relevant subset of peptide-HLAA*0201 complexes, the prediction success is 100%. Based on these structure predictions, a theoretical approach for prediction of cross-reactivity is developed and validated. This method involves the generation of quantitative structure-activity relationships using three-dimensional molecular descriptors and a genetic neural network. The generated relationships are highly predictive as proved by high cross-validated correlation coefficients (0.78-0.79). Together, the here developed theoretical methods open the door for efficient rational design of improved peptides to be used in immunotherapy. Résumé La réponse immunitaire spécifique contre des tumeurs dépend de la reconnaissance par les récepteurs des cellules T CD8+ de peptides antigéniques présentés par les complexes majeurs d'histocompatibilité (CMH) de classe I. Ces peptides sont utilisés comme cible dans l'immunothérapie par vaccins peptidiques. Afin d'augmenter la réponse immunitaire, les peptides sont modifiés de façon à améliorer l'affinité et/ou la résistance à la dégradation. Ceci nécessite de connaître la structure tridimensionnelle des complexes peptide-CMH. De plus, les peptides modifiés doivent être reconnus par des cellules T spécifiques du peptide natif. La structure de l'épitope doit donc être préservée et des structures détaillées des complexes peptide-CMH sont nécessaires. Dans cette thèse, le thème central est le développement des méthodes computationnelles de prédiction des structures des complexes peptide-CMH classe I et de la reconnaissance croisée. Des applications de ces méthodes de prédiction à l'immunothérapie sont également présentées. Premièrement, une méthode théorique de prédiction des structures des complexes peptide-CMH classe I est développée et validée. Cette méthode est basée sur un échantillonnage de l'espace conformationnel du peptide dans le contexte du récepteur CMH classe I par dynamique moléculaire. Les conformations sont évaluées par leurs énergies libres conformationnelles. La méthode est validée par sa capacité à reproduire 41 structures des complexes peptide-CMH classe I obtenues par cristallographie aux rayons X. Le succès prédictif général est de 83%. Pour le sous-groupe HLA-A*0201 de complexes de grande importance pour l'immunothérapie, ce succès est de 100%. Deuxièmement, à partir de ces structures prédites in silico, une méthode théorique de prédiction de la reconnaissance croisée est développée et validée. Celle-ci consiste à générer des relations structure-activité quantitatives en utilisant des descripteurs moléculaires tridimensionnels et un réseau de neurones couplé à un algorithme génétique. Les relations générées montrent une capacité de prédiction remarquable avec des valeurs de coefficients de corrélation de validation croisée élevées (0.78-0.79). Les méthodes théoriques développées dans le cadre de cette thèse ouvrent la voie du design de vaccins peptidiques améliorés.