Effect of central-residue replacements on the helical stability of a monomeric peptide

Curcumin modulates the self-assembly of the islet amyloid polypeptide by disassembling α-helix

Understanding how small molecules affect amyloid formation is of major biomedical and pharmaceutical importance due to the association of amyloid with incurable diseases including Alzheimer's, Parkinson's, and type II diabetes. Using solution state (1)H NMR, we demonstrate that curcumin, a planar biphenolic compound found in the Indian spice turmeric, delays the self-assembly of islet amyloid polypeptide to NMR-invisible assemblies. Accompanying circular dichroism studies show that curcumin disassembles α-helix in maturing assemblies of IAPP. The amount of α-helix disassembled correlates with predicted and experimentally determined helical content of IAPP obtained by others. Taken together, these results indicate that curcumin modulates IAPP self-assembly by unfolding α-helix on pathway to amyloid. The implications of this work in the elucidation of the mechanism for amyloid formation by IAPP in the presence of curcumin are discussed.

Intrinsic amino acid side-chain hydrophilicity/hydrophobicity coefficients determined by reversed-phase high-performance liquid chromatography of model peptides: comparison with other hydrophilicity/hydrophobicity scales

An accurate determination of the intrinsic hydrophilicity/hydrophobicity of amino acid side-chains in peptides and proteins is fundamental in understanding many area of research, including protein folding and stability, peptide and protein function, protein-protein interactions and peptide/protein oligomerization, as well as the design of protocols for purification and characterization of peptides and proteins. Our definition of intrinsic hydrophilicity/hydrophobicity of side-chains is the maximum possible hydrophilicity/hydrophobicity of side-chains in the absence of any nearest-neighbor effects and/or any conformational effects of the polypeptide chain that prevent full expression of side-chain hydrophilicity/hydrophobicity. In this review, we have compared an experimentally derived intrinsic side-chain hydrophilicity/hydrophobicity scale generated from RP-HPLC retention behavior of de novo designed synthetic model peptides at pH 2 and pH 7 with other RP-HPLC-derived scales, as well as scales generated from classic experimental and calculation-based methods of octanol/water partitioning of Nalpha-acetyl-amino-acid amides or free energy of transfer of free amino acids. Generally poor correlation was found with previous RP-HPLC-derived scales, likely due to the random nature of the peptide mixtures in terms of varying peptide size, conformation and frequency of particular amino acids. In addition, generally poor correlation with the classical approaches served to underline the importance of the presence of a polypeptide backbone when generating intrinsic values. We have shown that the intrinsic scale determined here is in full agreement with the structural characteristics of amino acid side-chains.

Structure and energetics of the hydrogen-bonded backbone in protein folding

We seek to understand the link between protein thermodynamics and protein structure in molecular detail. A classical approach to this problem involves assessing changes in protein stability resulting from added cosolvents. Under any given conditions, protein molecules in aqueous buffer are in equilibrium between unfolded and folded states, U(nfolded) <==> N(ative). Addition of organic osmolytes, small uncharged compounds found throughout nature, shift this equilibrium. Urea, a denaturing osmolyte, shifts the equilibrium toward U; trimethylamine N-oxide (TMAO), a protecting osmolyte, shifts the equilibrium toward N. Using the Tanford Transfer Model, the thermodynamic response to many such osmolytes has been dissected into groupwise free energy contributions. It is found that the energetics involving backbone hydrogen bonding controls these shifts in protein stability almost entirely, with osmolyte cosolvents simply dialing between solvent-backbone versus backbone-backbone hydrogen bonds, as a function of solvent quality. This reciprocal relationship establishes the essential link between protein thermodynamics and the protein's hydrogen-bonded backbone structure.

Thermodynamics Of α‐Helix Formation

The alpha-helix was the first proposed and experimentally confirmed secondary structure. The elegant simplicity of the alpha-helical structure, stabilized by hydrogen bonding between the backbone carbonyl oxygen and the peptide amide four residues away, has captivated the scientific community. In proteins, alpha-helices are also stabilized by the so-called capping interactions that occur at both the C- and the N-termini of the helix. This chapter provides a brief historical overview of the thermodynamic studies of the energetics of helix formation, and reviews recent progress in our understanding of the thermodynamics of helix formation.

Assemblages of prion fragments: novel model systems for understanding amyloid toxicity

We report the conformational and toxic properties of two novel fibril-forming prion amyloid sequences, GAVVGGLG (PrP(119-126)) and VVGGLGG (PrP(121-127)). The conformational preferences of these fragments were studied in differing microenvironments of TFE/water mixtures and SDS solution. Interestingly, with an increase in TFE concentration, PrP(119-126) showed a helical conformational propensity, whereas PrP(121-127) adopted a more random coil structure. In 5% SDS, PrP(119-126) showed more alpha-helical content than in TFE solution, and PrP(121-127) exhibited a predominantly random coil conformation. However, both peptides took a random coil conformation in water, and over time the random coil transformed into a beta-sheet structure with a significant percentage of helical conformation and beta-turn structure in PrP(119-126) and PrP(121-127), respectively, as observed with CD spectroscopy. The aged fibrils of PrP(119-126) were insoluble in SDS, and PrP(121-127) was extractable with SDS solution. These fibrils were characterized by transmission electron microscopy. Both PrP(119-126) and PrP(121-127) formed stable monolayer's consisting of multimeric assemblages at the air-water interface. Monomeric PrP(119-126) was more toxic to astrocytes than the control Abeta peptide; however, the fibrillar form of PrP(119-126) was less toxic to astrocytes. PrP(121-127) elicited moderate toxicity in both soluble and fibrillar forms on astrocytes. Furthermore, quenching experiments using acroyl-labeled PrP(119-126) and PrP(121-127) with eosin-labeled synaptosomal membrane revealed that these prion fragments bind to anion-exchange protein. The binding of PrP(119-126) and PrP(121-127) with a membrane microdomain (lipid raft) was also analyzed using pyrenated derivatives. We conclude that the formation of PrP(119-126) and PrP(121-127) fibrils is a concentration-dependent process that involves coil to sheet conversion with aging. PrP(119-126), the sequence with intrinsic helical propensity, is more toxic in monomer form, and the fibril formation in this case seems to be protective to cells. For PrP(121-127), the SDS-soluble fibrils are more cytotoxic, indicating that a higher order assemblage structure is required for cytotoxic activity of this peptide.

Modulation of the structural integrity of helix F in apomyoglobin by single amino acid replacements

The conformational features of native and mutant forms of sperm-whale apomyoglobin (apoMb) at neutral pH were probed by limited proteolysis experiments utilizing up to eight proteases of different substrate specificities. It was shown that all proteases selectively cleave apoMb at the level of chain segment 82-94 (HEAELKPLAQSHA), encompassing helix F in the X-ray structure of the holo form of the native protein; for example, thermolysin cleaves the Pro 88-Leu 89 peptide bond. These results indicate that helix F is highly flexible or largely disrupted in apoMb. Because helix F contains the helix-breaking Pro 88 residue, we propose that helix F is kept in place in the native holo protein by a variety of helix-heme stabilizing interactions. To modulate the stability of helix F, the Pro88Ala and Pro88Gly mutants were prepared by site-directed mutagenesis, and their conformational properties investigated by both far-UV circular dichroism spectroscopy and limited proteolysis. The helix content of the Pro88Ala mutant was somewhat enhanced with respect to that of both native and Pro88Gly mutant, as expected from the fact that Ala is the strongest helix inducer among the 20 amino acid residues. The rate of limited proteolysis of the three apoMb variants by thermolysin and proteinase K was in the order native > Pro88Gly >> Pro88Ala, in agreement with the scale of helix propensity of Ala, Gly, and Pro. The possible role of the flexible/unfolded chain segment 82-94 for the function and fate of apoMb at the cellular level is discussed.

Conformational polymorphism of the amyloidogenic peptide homologous to residues 113-127 of the prion protein

Conformational transitions are thought to be the prime mechanism of amyloid formation in prion diseases. The prion proteins are known to exhibit polymorphic behavior that explains their ability of "conformation switching" facilitated by structured "seeds" consisting of transformed proteins. Oligopeptides containing prion sequences showing the polymorphism are not known even though amyloid formation is observed in these fragments. In this work, we have observed polymorphism in a 15-residue peptide PrP (113-127) that is known to form amyloid fibrils on aging. To see the polymorphic behavior of this peptide in different solvent environments, circular dichroism (CD) spectroscopic studies on an aqueous solution of PrP (113-127) in different trifluoroethanol (TFE) concentrations were carried out. The results show that PrP (113-127) have sheet preference in lower TFE concentration whereas it has more helical conformation in higher TFE content (>40%). The structural transitions involved in TFE solvent were studied using interval-scan CD and FT-IR studies. It is interesting to note that the alpha-helical structure persists throughout the structural transition process involved in amyloid fibril formation implicating the involvement of both N- and C-terminal sequences. To unravel the role of the N-terminal region in the polymorphism of the PrP (113-127), CD studies on another synthetic peptide, PrP (113-120) were carried out. PrP(113-120) exhibits random coil conformation in 100% water and helical conformation in 100% TFE, indicating the importance of full-length sequence for beta-sheet formation. Besides, the influence of different chemico-physical conditions such as concentration, pH, ionic strength, and membrane like environment on the secondary structure of the peptide PrP (113-127) has been investigated. At higher concentration, PrP (113-127) shows features of sheet conformation even in 100% TFE suggesting aggregation. In the presence of 5% solution of sodium dodecyl sulfate, PrP (113-127) takes high alpha-helical propensity. The environment-dependent conformational polymorphism of PrP (113-127) and its marked tendency to form stable beta-sheet structure at acidic pH could account for its conformation switching behavior from alpha-helix to beta-sheet. This work emphasizes the coordinative involvement of N-terminal and C-terminal sequences in the self-assembly of PrP (113-127).

Lifson-Roig nucleation for alpha-helices in trifluoroethanol: context has a strong effect on the helical propensity of amino acids

We have investigated the effect of substituting each of 19 common amino acids (excluding P) at the X position in the peptide acetyl-Y-VAEAK-TSXSR-VAEAK-NH(2). This very different peptide is of interest because, in previous work, we showed that nucleation in the strong alpha-helix-forming pentamers VAEAK was unable to propagate the helix structure through the sequence TSDSR, which is neither a helix-forming sequence nor a breaker, but is indifferent to helix formation. Substitution in the center of the indifferent sequence reveals an interesting measure of the helical propensity for the 19 amino acids. CD spectra were measured in various mixtures of buffer and 2,2,2-trifluoroethanol (TFE), and then analyzed for helix propensity of the amino acids using the Lifson-Roig model. However, the nucleation parameter in the Lifson-Roig model has never been measured for TFE. We have empirically found that the nucleation parameter for a solvent can be determined from the data normally used to determine only the propagation parameters. The results of the analysis of the CD show that most amino acids are excellent or good helix formers in 90% TFE, while amino acids D, W, F and G are poor helix formers for the indifferent pentamer sequence. The helix propensity of the 19 amino acids is quite different from the helix propensity measured in other peptide sequences, demonstrating the context dependence of this property. The results as a function of alcohol concentration confirm that the relative order of helical propensity of amino acids changes with solvent environment. Clearly, the prediction of alpha-helical secondary structure from protein sequence requires more than a single helical propensity for each amino acid.

Helix-coil transitions re-visited

The thermally-induced helix-coil transition in polyamino acids is a good model for determining the helix-forming propensities of amino acids but not for the two-state folding/unfolding transition in globular proteins. The equilibrium and kinetic treatments of the helix-coil transition are summarized here together with a description of applications to various types of homopolymers and copolymers. Attention is then focused on the helix-coil transition in poly-L-alanine as an example of a non-polar polyamino acid. To render such a non-polar polymer water soluble, it is necessary to introduce polar amino acids such as lysines, but care must be taken as to the location of such polar residues. If they are attached as end groups, as in a triblock copolymer, they do not perturb the helix-forming tendency of the central poly-L-alanine block significantly, but if they are introduced within the sequence of alanine residues, then the hydration properties of the lysines dominate the behavior of the resulting copolymer, thereby leading to erroneous values of the parameters characterizing the helix-forming tendency of the alanines. Neutral but polar residues, such as glutamines, also exhibit hydration-dominating properties but less so than charged lysines. Some details of the calculations for an alanine/glutamine copolymer are presented here. It is concluded that random copolymers based on a neutral water-soluble host provide reliable information about the helix-forming tendencies of amino acid residues that are introduced as guests among such neutral host residues.

Photo-control of peptide helix content by an azobenzene cross-linker: steric interactions with underlying residues are not critical

Photo-control of protein conformation could prove useful for probing function in diverse biological systems. Recently, we reported photo-switching of helix content in a short peptide containing an azobenzene cross-linker between cysteine residues at positions i and i + 7 in the sequence. In the original sequence, underlying residues at positions i + 3 and i + 4 were made bulky as preliminary modelling suggested that this would enhance photo-control of helix content. To test this hypothesis, peptides with Val, Aib; Ile, Aib; and Ala, Ala at positions i + 3 and i + 4 were synthesized, cross-linked and characterized. Before cross-linking, the peptides show distinct conformational behaviours: two with differing helix/coil mixtures whereas the other has a circular dichroism (CD) spectrum characteristic of beta-sheet and a tendency to aggregate. However, upon cross-linking the peptides have very similar CD spectra: predominantly random coil in the dark but predominantly helical upon irradiation. These results refute the original hypothesis. Steric interactions between the linker and underlying residues do not appear to be critical for photo-switching behaviour. When the cross-linking bridge is lengthened by replacing the i, i + 7 cysteine residues with homocysteine, a lower degree of photo-control of helicity is observed. Furthermore, a non-cross-linking version of the azobenzene reagent is shown not to produce any photo-control of helicity. We conclude that the intramolecular cross-link is essential for photo-switching and that it should be applicable to a wide range of peptides and proteins.

Energetics of three-state unfolding of a protein: canine milk lysozyme

Thermodynamics of thermal transitions of a calcium-binding lysozyme, canine milk lysozyme (CML), was studied using differential scanning calorimetry and compared with those for homologous proteins, human alpha-lactalbumin (alpha-hLA) and equine milk lysozyme (EML). The results showed that CML and EML exhibit two clear heat absorption peaks in the absence of calcium ions (apo-form), although the cooperative thermal transition of alpha-hLA is apparently absent in this form. The first peak represents the unfolding transition from the native to an unfolding intermediate state (N-I transition) and the second peak represents that from the intermediate to the thermally unfolded state (I-U transition). We interpret that the cooperative thermal transition, which is observed between the intermediate and the thermally unfolded states of CML and EML, comes from the native-like packing interaction in their intermediate states. Furthermore, to examine the role of the stabilization mechanism of CML intermediate, we constructed four variant CMLs (H21G, I56L, A93S and V109K), in which the residues of CML are substituted for those of EML, and also investigated their thermal stability. Especially the His21 and Val109 of CML play a role in stabilization of the intermediate state and their contributions to the unfolding free energy are estimated to be 2.0 and 1.8 kJ/mol, respectively. From the results of the mutational analysis, a few differences in the local helical interactions within the alpha-domain are found to be predominant in stabilizing the intermediate state.

Design of a peptide inhibitor that blocks the cell fusion mediated by glycoprotein 41 of human immunodeficiency virus type 1

Fusion between the envelope of HIV and the plasma membrane of target cells is mediated by gp41, the envelope glycoprotein of HIV. Peptides derived from the membrane-proximal helical motif of the extracellular domain of gp41 effectively inhibit the infection of HIV, and their inhibitory activities are known to be correlated with the helical propensity of the peptides. We have designed small peptides that can form a stable alpha helix and thereby inhibit gp120/41-mediated cell fusion. A 19-mer peptide from the membrane-proximal helical motif of gp41 had no secondary structure in solution, and failed to block gp41-mediated cell fusion. When amino acids with low helical propensity were substituted, and helix-capping sequences were introduced at both ends of the peptides, the modified peptides formed a stable helical structure. They also bound to the coiled-coil motif of gp41 presented at the C terminus of thioredoxin and blocked the cell fusion mediated by gp120/41. These results implied that such modification was enough to change a short peptide derived from gp41 into a potent inhibitor against the infection of HIV.

Amino acid conformational preferences and solvation of polar backbone atoms in peptides and proteins

Amino acids in peptides and proteins display distinct preferences for alpha-helical, beta-strand, and other conformational states. Various physicochemical reasons for these preferences have been suggested: conformational entropy, steric factors, hydrophobic effect, and backbone electrostatics; however, the issue remains controversial. It has been proposed recently that the side-chain-dependent solvent screening of the local and non-local backbone electrostatic interactions primarily determines the preferences not only for the alpha-helical but also for all other main-chain conformational states. Side-chains modulate the electrostatic screening of backbone interactions by excluding the solvent from the vicinity of main-chain polar atoms. The deficiency of this electrostatic screening model of amino acid preferences is that the relationships between the main-chain electrostatics and the amino acid preferences have been demonstrated for a limited set of six non-polar amino acid types in proteins only. Here, these relationships are determined for all amino acid types in tripeptides, dekapeptides, and proteins. The solvation free energies of polar backbone atoms are approximated by the electrostatic contributions calculated by the finite difference Poisson-Boltzmann and the Langevin dipoles methods. The results show that the average solvation free energy of main-chain polar atoms depends strongly on backbone conformation, shape of side-chains, and exposure to solvent. The equilibrium between the low-energy beta-strand conformation of an amino acid (anti-parallel alignment of backbone dipole moments) and the high-energy alpha conformation (parallel alignment of backbone dipole moments) is strongly influenced by the solvation of backbone polar atoms. The free energy cost of reaching the alpha conformation is by approximately 1.5 kcal/mol smaller for residues with short side-chains than it is for the large beta-branched amino acid residues. This free energy difference is comparable to those obtained experimentally by mutation studies and is thus large enough to account for the distinct preferences of amino acid residues. The screening coefficients gamma(local)(r) and gamma(non-local)(r) correlate with the solvation effects for 19 amino acid types with the coefficients between 0.698 to 0.851, depending on the type of calculation and on the set of point atomic charges used. The screening coefficients gamma(local)(r) increase with the level of burial of amino acids in proteins, converging to 1.0 for the completely buried amino acid residues. The backbone solvation free energies of amino acid residues involved in strong hydrogen bonding (for example: in the middle of an alpha-helix) are small. The hydrogen bonded backbone is thus more hydrophobic than the peptide groups in random coil. The alpha-helix forming preference of alanine is attributed to the relatively small free energy cost of reaching the high-energy alpha-helix conformation. These results confirm that the side-chain-dependent solvent screening of the backbone electrostatic interactions is the dominant factor in determining amino acid conformational preferences.

Guidelines for membrane protein engineering derived from de novo designed model peptides

Notwithstanding great advances in the engineering and structural analysis of globular proteins, relatively limited success has been achieved with membrane proteins--due largely to their intrinsic high insolubility and the concomitant difficulty in obtaining crystals. Progress with de novo synthesis of model membrane-interactive peptides presents an opportunity to construct simpler peptides with definable structures, and permits one to approach an understanding of the properties of the membrane proteins themselves. In the present article, we review how our laboratory and others have used peptide approaches to assess the detailed interactions of peptides with membranes, and primary folding at membrane surfaces and in membranes. Structural studies of model peptides identified the existence of a "threshold hydrophobicity," which controls spontaneous peptide insertion into membranes. Related studies of the relative helicity of peptides in organic media such as n-butanol indicate that the helical propensity of individual residues--not simply their hydrophobicity--may dictate the conformations of peptides in membranes. The overall experimental results provide fundamental guidelines for membrane protein engineering.

Structural cassette mutagenesis in a de novo designed protein: proof of a novel concept for examining protein folding and stability

The solution to the protein folding problem lies in defining the relative energetic contributions of short-range and long-range interactions. In other words, the tendency of a stretch of amino acids to adopt a final secondary structural fold is context dependent. Our approach to this problem is to address whether an amino acid sequence, a "cassette," with a defined secondary structure in the three-dimensional structure of a native protein, can adopt a different conformation when placed into a different protein environment. Thus, we designed de novo a disulfide-bridged two-stranded alpha-helical parallel coiled coil, where each polypeptide chain consisted of 39 residues, as a "cassette holder." The 11-residue cassette would be inserted into the center of each polypeptide chain between the two nucleating alpha-helices to replace the control sequence. This Structural Cassette Mutagenesis model permits the analysis of short-range interactions within the inserted cassette as well as long-range interactions between the nucleating helices and the cassette region. The cassette holder, with a control sequence as the cassette, had a GdnHCl transition midpoint during denaturation of 5.6M. To demonstrate the feasibility of our model, an 11-residue beta-strand cassette from an immunoglobulin fold was inserted. The cassette was fully induced into the alpha-helical conformation with a [GdnHCl]1/2 value of 3.2M. To demonstrate the importance of short-range interactions (beta-sheet/alpha-helical propensities of amino acid side chains) in modulating structure and stability, a series of 1-5 threonine residues (highest beta-sheet propensity) were substituted into the solvent-exposed portions of the cassette in the alpha-helical conformation. Each successive substitution systematically decreased the stability of the coiled coil with peptide T4b (4 Thr residues) having a [GdnHCl]1/2 value of 2.2M. The single substitution of Ile in the hydrophobic core of the cassette with Ala or Thr had the most dramatic effect on protein stability (peptide 120T, [GdnHCl]1/2 value of 1.4M). Though these substitutions were able to modulate stability, they were not able to disrupt the alpha-helical conformation of the cassette, showing the importance of the nucleating alpha-helices on either side of the cassette in controlling conformation of the cassette. We have demonstrated the feasibility of our model protein to accept a beta-strand cassette. The effect of cassettes containing other beta-strands, beta-turns, loops, regions of undefined structure, and helical segments on conformation and stability of our model protein will also be determined.

The relative order of helical propensity of amino acids changes with solvent environment

A model peptide of sequence Ac-Y-VAXAK-VAXAK-VAXAK-NH(2), where X is substituted with one of nineteen amino acids (P excluded), was synthesized and titrated with methanol to study helical propensity as a function of solvent environment. The CD spectra of these peptides are largely random coil in 2 mM sodium phosphate buffer (pH 5.5) and show a conformational change to alpha-helix with increasing methanol content. Singular value decomposition was used to correct the CD spectra for the absorbing side chains of W, Y, F, C, and M, and this correction can be substantial. With correction both W and F become good helix formers. The free energy for helix propagation was calculated using the Lifson-Roig statistical model for each of the nineteen amino acids at each point in their titration. The results show that the rank order of helical propensity for the nineteen amino acids changes with solvent environment. This result will be particularly important if proteins undergo hydrophobic collapse before secondary structures are formed, because amino acids can then see different solvent environments as the secondary structures are formed. Related amino acids are found to have interesting correlations in the shape of their titration curves. This finding provides one explanation for the limiting 70% accuracy in predicting secondary structure from sequence, since the helical propensities used are calculated for an average solvent environment. Proteins 2000;39:132-141.

Helix-enhancing propensity of fluoro and alkyl alcohols: influence of pH, temperature and cosolvent concentration on the helical conformation of peptides

We have analyzed the effects of trifluoroethanol (TFE) and three other alcohols(1-propanol, 2-propanol and hexafluoro-2-propanol) on S-peptide (residues 1-20) of ribonuclease A, an analog of S-peptide (QHM-->AAA, Sa-peptide) and TC-peptide (residues 295-316) of thermolysin to assess the helix-enhancing propensity of fluoro and alkyl alcohols under different environmental conditions of cosolvent concentration, pH and temperature by circular dichroism (CD). The dependence of cosolvent concentration on helix-induction showed a plateauing effect in all cases. 1-Propanol and 2-propanol were as effective as TFE in all the three peptides. Hexafluoro-2-propanol (HFIP) was a better helix enhancer in all cases however, the relative effectiveness varied with the peptide sequence. The alcohol transitions were analyzed assuming a two-state transition. The free energy decreased linearly in the cosolvent concentration range of 0-5 m for all the three peptides. The m-value (constant of proportionality) varied between peptides but was similar for any given peptide for TFE, 1-propanol or 2-propanol. The m-values of HFIP for all three peptides was much higher compared to other cosolvents. The isothermal cosolvent helix-induction curves for the three peptides exhibited similar features of shape and character for 1-propanol, 2-propanol and TFE. The additivity of cosolvent-induced helix formation was observed for different blends of alkyl and/or fluoro cosolvents. The pH-dependence of helix formation was observed in both TFE and 1-propanol solutions for S-peptide and TC-peptide, respectively, while in Sa-peptide, which was designed to perturb the pH-effect, helix formation was unaffected. The overall results provide some insight into the mechanism of cosolvent-mediated helix-enhancement in protein segments and are likely to facilitate optimization of conditions for cosolvent usage in chemistry and biology.

Role of hydrophobicity and solvent-mediated charge-charge interactions in stabilizing alpha-helices

A theoretical study to identify the conformational preferences of lysine-based oligopeptides has been carried out. The solvation free energy and free energy of ionization of the oligopeptides have been calculated by using a fast multigrid boundary element method that considers the coupling between the conformation of the molecule and the ionization equilibria explicitly, at a given pH value. It has been found experimentally that isolated alanine and lysine residues have somewhat small intrinsic helix-forming tendencies; however, results from these simulations indicate that conformations containing right-handed alpha-helical turns are energetically favorable at low values of pH for lysine-based oligopeptides. Also, unusual patterns of interactions among lysine side chains with large hydrophobic contacts and close proximity (5-6 A) between charged NH3+ groups are observed. Similar arrangements of charged groups have been seen for lysine and arginine residues in experimentally determined structures of proteins available from the Protein Data Bank. The lowest-free-energy conformation of the sequence Ac-(LYS)6-NMe from these simulations showed large pKalpha shifts for some of the NH3+ groups of the lysine residues. Such large effects are not observed in the lowest-energy conformations of oligopeptide sequences with two, three, or four lysine residues. Calculations on the sequence Ac-LYS-(ALA)4-LYS-NMe also reveal low-energy alpha-helical conformations with interactions of one of the LYS side chains with the helix backbone in an arrangement quite similar to the one described recently by (Proc. Natl. Acad. Sci. U.S.A. 93:4025-4029). The results of this study provide a sound basis with which to discuss the nature of the interactions, such as hydrophobicity, charge-charge interaction, and solvent polarization effects, that stabilize right-handed alpha-helical conformations.

Study of alpha-helix to beta-strand to beta-sheet transitions in amyloid: the role of segregated hydrophobic beta-strands

A conformational analysis including three polypeptide chains known to be amyloidogenic and neurotoxic has shown the occurrence of low probability hydrophobic beta-strands stabilized by intramolecular interactions. It is argued that by engaging in non-bonded and hydrophobic interactions these beta-strands seed the assembly of beta-sheets in amyloid fibrils following a non-cooperative mechanism dissimilar to beta-sheet folding in proteins. Molecular models of amyloid fibrils formed by such beta-strand templates were built. It is shown that the parallel alignment of beta-strands creates an extensive hydrophobic surface whereas the antiparallel alignment causes the formation of hydrophobic and hydrophilic aggregates. A comparison with experimental data and previous calculations is established.

Insight into signal transduction: structural alterations in transmembrane helices probed by multi-1 ns molecular dynamics simulations

The hypothesis of structural alteration in transmembrane helices for signal transduction process is viewed by molecular dynamics simulation techniques. For the c-erbB-2 transmembrane domain involved in oncogenicity, the occurrence of conformational changes has been previously described as transition from the alpha to pi helix. This dynamical feature is thoroughly analyzed for the wild phenotype and oncogenic sequences from a series of 18 simulations carried out on one nanosecond time scale. We show that these structural events do not depend upon the conditions of simulations like force field or starting helix coordinates. We demonstrate that the oncogenic mutations Val659 Glu, Gln and Asp do not prevent the transition. Furthermore, we show that beta branched residues, in conjunction with Gly residues in the c-erbB-2 sequence, act as destabilizers for the alpha helix structure, pi deformations are tightly related to other local structural motifs found in soluble and membrane proteins. These structural alterations are discussed in term of structure-activity relationships for the c-erbB-2 activating mechanism mediated by transmembrane domain dimerization.

The role of PII conformations in the calculation of peptide fractional helix content

Changes in the temperature, pH, ionic strength, or denaturant concentration of aqueous solutions of the monomeric non-alpha-helical peptide acetylYEAAAKEAPAKEAAAKAamide generate changes in its dichroic spectrum characteristic for a conformational transition. This transition has the characteristic features of a residue PII/unstructured conformational equilibrium in which PII denotes an extended left-handed helical conformation and unstructured denotes all the remaining conformations in a random coil ensemble. Replacement of the proline residue facilitates population of residues in an alpha-helical conformation. However, the ellipticity values for these non-proline peptides merge with the ellipticity of the proline peptide as the population of residues in the alpha-helix conformation is diminished. This convergence suggests that all residues in a host/guest peptide series of the same length share a common PII/unstructured conformational equilibrium in a given solvent. We propose that the fractional helix content of peptides within such a series may be estimated by using a two-state calculation in which the ellipticity for the non-alpha-helix conformations is provided by a peptide having a central proline guest residue.

Conformation of a water-soluble beta-sheet model peptide. A circular dichroism and Fourier-transform infrared spectroscopic study of double D-amino acid replacements

Among peptide secondary structures beta-sheet domains have been much less intensively studied than alpha-helical conformations, mainly because of the lack of well characterized model peptides. In the present paper the secondary structure of a water-soluble de novo peptide consisting of 26 amino acids (DPKGDPKGVTVTVTVTVTGKGDPKPD-NH2) and the corresponding double D-amino acid replacement set have been studied by circular dichroism and Fourier-transform infrared spectroscopy. The model peptide was found to be unstructured in aqueous solution at peptide concentrations < 10(-3) mol/L but to adopt a predominantly beta-sheet structure in the presence of 15 mM sodium dodecyl sulfate or at apolar/water interfaces. Although the peptide is composed of amino acids with low helical propensity, it formed a single-stranded helical structure in aqueous trifluoroethanol. The D-amino acid replacement set was synthesized in order to study the conformational stability of the model peptide selectively in distinct regions. The data show that both the alpha-helix present in 50% trifluoroethanol as well as the beta-sheet domain formed in the presence of sodium dodecyl sulfate or at apolar/water interfaces, are located in the region between Val9 and Thr18. Pairwise substitution of adjacent amino acids by their corresponding D-amino acids provides a pronounced beta-sheet disturbance. These findings demonstrate that double D-amino acid replacements may be used to locate beta-sheet domains in peptides.

Hydrophobic solvation in aqueous trifluoroethanol solution

The titration of an aqueous solution of a de novo designed peptide with trifluoroethanol (TFE) shows complete helix formation with the addition of only 30% TFE. A molecular simulation of the peptide, in which a single shell of TFE molecules initially surrounds the peptide, reveals preferred sites of solvent interaction. The TFE molecules show greater preference for the hydrophobic compared with hydrophilic side chains. The helix-enhancing ability of TFE in aqueous solution may be rationalized in terms of stabilizing the hydrophobic collapse of apolar side chains of the formed helix.

Relationship of sidechain hydrophobicity and alpha-helical propensity on the stability of the single-stranded amphipathic alpha-helix

The aim of the present investigation is to determine the effect of alpha-helical propensity and sidechain hydrophobicity on the stability of amphipathic alpha-helices. Accordingly, a series of 18-residue amphipathic alpha-helical peptides has been synthesized as a model system where all 20 amino acid residues were substituted on the hydrophobic face of the amphipathic alpha-helix. In these experiments, all three parameters (sidechain hydrophobicity, alpha-helical propensity and helix stability) were measured on the same set of peptide analogues. For these peptide analogues that differ by only one amino acid residue, there was a 0.96 kcal/mole difference in alpha-helical propensity between the most (Ala) and the least (Gly) alpha-helical analogue, a 12.1-minute difference between the most (Phe) and the least (Asp) retentive analogue on the reversed-phase column, and a 32.3 degrees C difference in melting temperatures between the most (Leu) and the least (Asp) stable analogue. The results show that the hydrophobicity and alpha-helical propensity of an amino acid sidechain are not correlated with each other, but each contributes to the stability of the amphipathic alpha-helix. More importantly, the combined effects of alpha-helical propensity and sidechain hydrophobicity at a ratio of about 2:1 had optimal correlation with alpha-helix stability. These results suggest that both alpha-helical propensity and sidechain hydrophobicity should be taken into consideration in the design of alpha-helical proteins with the desired stability.

Helix design, prediction and stability

Recent work revealing that our knowledge is now sufficient to build a reasonable quantitative model for the helix/coil transition in heteropolypeptides represents a watershed in research into alpha-helix stability, prediction and design. The opportunity is presented to design specific alpha-helix propensity patterns that may be used both to modify thermodynamic properties of target proteins and peptides, and for de novo protein design. Despite these advances, the picture is not yet complete and further studies of still poorly characterized factors are required to obtain a more precise understanding of alpha-helix stability.

Structural analysis of the N- and C-termini in a peptide with consensus sequence

We present a structural analysis of a peptide, the sequence of which includes amino acids that show preferences for specific positions near the N- and C-termini in protein helices. This peptide has the sequence ac-YMSEDELKAAEAAFKRHGVP-amide, which includes a strong version of an N-terminal Harper-Rose capping box structure as well as a Gly located close to the C-terminus designed to elucidate its role in C-terminal capping. The sequence of five residues at the middle is inserted to separate effects at the two ends via a helix-stabilizing linker. Application of a simulated annealing procedure using interproton distance constraints derived from 1H NOESY experiments in water reveals the presence of a C-terminal structure in this model. The C-terminus forms a folded back structure in a significant fraction of structures generated by the annealing, in most of which Gly assumes an alpha L conformation. This structure occurs within a highly flexible region of the molecule and hence is occupied only a fraction of the time.

Periodicity of polar and nonpolar amino acids is the major determinant of secondary structure in self-assembling oligomeric peptides

The tendency of a polypeptide chain to form alpha-helical or beta-strand secondary structure depends upon local and nonlocal effects. Local effects reflect the intrinsic propensities of the amino acid residues for particular secondary structures, while nonlocal effects reflect the positioning of the individual residues in the context of the entire amino acid sequence. In particular, the periodicity of polar and nonpolar residues specifies whether a given sequence is consistent with amphiphilic alpha-helices or beta-strands. The importance of intrinsic propensities was compared to that of polar/nonpolar periodicity by a direct competition. Synthetic peptides were designed using residues with intrinsic propensities that favored one or the other type of secondary structure. The polar/nonpolar periodicities of the peptides were designed either to be consistent with the secondary structure favored by the intrinsic propensities of the component residues or in other cases to oppose these intrinsic propensities. Characterization of the synthetic peptides demonstrated that in all cases the observed secondary structure correlates with the periodicity of the peptide sequence--even when this secondary structure differs from that predicted from the intrinsic propensities of the component amino acids. The observed secondary structures are concentration dependent, indicating that oligomerization of the amphiphilic peptides is responsible for the observed secondary structures. Thus, for self-assembling oligomeric peptides, the polar/nonpolar periodicity can overwhelm the intrinsic propensities of the amino acid residues and serves as the major determinant of peptide secondary structure.

Interactions between hydrophobic side chains within alpha-helices

The thermodynamic basis of helix stability in peptides and proteins is a topic of considerable interest. Accordingly, we have computed the interactions between side chains of all hydrophobic residue pairs and selected triples in a model helix, using Boltzmann-weighted exhaustive modeling. Specifically, all possible pairs from the set Ala, Cys, His, Ile, Leu, Met, Phe, Trp, Tyr, and Val were modeled at spacings of (i, i + 2), (i, i + 3), and (i, i + 4) in the central turn of a model poly-alanyl alpha-helix. Significant interactions--both stabilizing and destabilizing-- were found to occur at spacings of (i, i + 3) and (i, i + 4), particularly in side chains with rings (i.e., Phe, Tyr, Trp, and His). In addition, modeling of leucine triples in a helix showed that the free energy can exceed the sum of pairwise interactions in certain cases. Our calculated interaction values both rationalize recent experimental data and provide previously unavailable estimates of the constituent energies and entropies of interaction.

The effect of environment on the stability of an integral membrane helix: molecular dynamics simulations of surfactant protein C in chloroform, methanol and water

A series of three molecular dynamics simulations at 300 K in explicit solvent environments of chloroform, methanol and water has been performed on the pulmonary surfactant lipoprotein, SP-C, comprising several consecutive valine residues in order to investigate the stability of the alpha-helical conformation. Two additional simulations were performed on truncated SP-C with a five-residue N-terminal deletion at 300 K and 500 K in water, the high temperature run in order to increase the rate of peptide denaturation. Indications of destabilization appear in chloroform during 1 ns while the SP-C alpha-helix is remarkably stable during 1 ns in methanol and water. In particular the polyvalyl part comprising residues Val15 to Val21 remains intact even at elevated temperature, and the valines do not disrupt the alpha-helical conformation. The valyl-rotamer sampling is partly restricted. Unfolding takes place successively along the primary sequence starting from the C-terminal end. Factors affecting polypeptide stability in molecular dynamics simulations are addressed. The intrinsic helix-forming tendency of valine residues and its dependence on the sequence context, and the role of the solvent environment in stabilizing or destabilizing an alpha-helical fold, are discussed.

Competing interactions contributing to alpha-helical stability in aqueous solution

The stability of a 15-residue peptide has been investigated using CD spectroscopy and molecular simulation techniques. The sequence of the peptide was designed to include key features that are known to stabilize alpha-helices, including ion pairs, helix dipole capping, peptide bond capping, and aromatic interactions. The degree of helicity has been determined experimentally by CD in three solvents (aqueous buffer, methanol, and trifluoroethanol) and at two temperatures. Simulations of the peptide in the aqueous system have been performed over 500 ps at the same two temperatures using a fully explicit solvent model. Consistent with the CD data, the degree of helicity is decreased at the higher temperature. Our analysis of the simulation results has focused on competition between different side-chain/side-chain and side-chain/main-chain interactions, which can, in principle, stabilize the helix. The unfolding in aqueous solution occurs at the amino terminus because the side-chain interactions are insufficient to stabilize both the helix dipole and the peptide hydrogen bonds. Loss of capping of the peptide backbone leads to water insertion within the first peptide hydrogen bond and hence unfolding. In contrast, the carboxy terminus of the alpha-helix is stable in both simulations because the C-terminal lysine residue stabilizes the helix dipole, but at the expense of an ion pair.

The role of amphipathicity in the folding, self-association and biological activity of multiple subunit small proteins

The effect that altering amphipathicity has on the folding process and self association of melittin, a small model protein, has been investigated using single amino acid substitutions of lysine 7, a residue distant from the contact residues involved in the hydrophobic core of tetrameric melittin. While substitutions of such a residue were not expected to interfere with the packing process, the largest alterations in the potential overall amphipathicity of melittin were found to prevent the folding into an alpha-helical conformation to occur and, in turn, to prevent the self association. Amphipathic alpha-helices were found to be a key determining feature in the early folding process of the self association of peptides and protein segments. Those substitutions, which prevented the inducible amphipathic folding ability, were also found to result in a loss in hemolytic and antimicrobial activity. These results, combined with studies of the binding to artificial liposomes and to polysialic acids, indicate that the losses in activity were due to an initial inability to be induced into an amphipathic alpha-helix and to self associate. Ultimately, melittin's self association is proposed to be required to penetrate the carbohydrate barrier present in biological membranes.

Dissecting titin into its structural motifs: identification of an alpha-helix motif near the titin N-terminus

Titin, also known as connectin, is a giant modular protein specifically found in vertebrate striated muscle. Since the huge size of titin does not allow a direct structure determination, we have started a long-term project to characterize the protein by cutting it into smaller domains or structural units. The major part of the titin sequence is assembled by modules approximately 100 amino acids long that belong to two major protein superfamilies. Most of these modules are linked together by stretches of variable length with unique sequence. No direct structural characterization has been achieved so far for any of these linkers. We present here a study of a stretch located in the titin N-terminus and part of a linker between two modules. Our attention was drawn toward this region because it shows 100% probability to form a coiled coil when analyzed by a prediction program. A synthetic 38 amino acid peptide spanning such a sequence was studied in aqueous solution by circular dichroism, nuclear magnetic resonance, and analytical ultracentrifugation at various pH, salt, and peptide concentrations. Under all conditions, it shows a strong tendency to form alpha-helical structures. In the presence of salt, this conformation is associated with the formation of helical bundles below pH 5. Above pH 5, any aggregate breaks, and the titin peptide is a monomeric helix in equilibrium with its random coil conformation. We discuss the factors which stabilize the helical conformation and the possible role of this stretch in vivo.

Designed polyanionic coiled-coil proteins: acceleration of heparin cofactor II inhibition of thrombin

Novel polyanionic proteins were designed to increase the rate of heparin cofactor II (HC) inhibition of alpha-thrombin, an essential protease in the coagulation cascade. Two alpha-helical coiled-coil proteins, a 62-residue dimer containing 8 Glu residues (E8C) and a 104-residue dimer containing 14 Glu residues (E14C), plus two 31-residue control peptides containing 8 Glu residues each (E8A and E8B), were chemically synthesized, structurally characterized and enzymatically assayed. Circular dichroic spectrophotometry indicated that both E8C and E14C formed stable two-chain alpha-helical coiled coils at pH 7 and 25 degrees C. The control peptides were only partially alpha-helical. E14C remained folded at 90 degrees C but E8C was half unfolded at 49 degrees C. Coiled-coil proteins E8C and E14C maximally accelerated by 35- and 33-fold, respectively, the rate of HC inhibition of alpha-thrombin. None of these compounds accelerated antithrombin inhibition of alpha-thrombin, and neither control peptide accelerated HC inhibition of alpha-thrombin. Acceleration of the HC inhibition of alpha-thrombin showed bimodal dependence on the concentration of the polyanionic protein, which is consistent with formation of a HC-coiled-coil-thrombin ternary complex. The results suggest that antithrombotic polyanionic alpha-helical coiled-coil proteins can be designed and synthesized and that the occurrence of secondary structure can be correlated with biological activity.

Structural stability of short subsequences of the tropomyosin chain

The native tropomyosin molecule is a parallel, registered, alpha-helical coiled coil made from two 284-residue chains. Long excised subsequences (> or = 95 residues) form the same structure with comparable thermal stability. Here, we investigate local stability using shorter subsequences (20-50 residues) that are chemically synthesized or excised from various regions along the protein chain. Thermal unfolding studies of such shorter peptides by CD in the same solvent medium used in extant studies of the parent protein indicate very low helix content, almost no coiled-coil formation, and high thermal lability of such secondary structure as does form. This behavior is in stark contrast to extant data on leucine-zipper peptides and short "designed" synthetic peptides, many of which have high alpha-helix content and form highly stable coiled coils. The existence of short coiled coils calls into question the older idea that short subsequences of a protein have little structure. The present study supports the older view, at least in its application to tropomyosin. The intrinsic local alpha-helical propensity and helix-helix interaction in this prototypical alpha-helical protein is sufficiently weak as to require not only dimerization, but macro-molecular amplification in order to attain its native conformation in common benign media near neutral pH.

Exploring the peptide 3(10)-helix reversible alpha-helix equilibrium with double label electron spin resonance

Over the last several years we have used spin labeling as a means for exploring the structure of helical peptides. Two nitroxide labels are engineered into a peptide sequence and distances are ranked with electron spin resonance (ESR). We have found that there is a significant amount of 3(10)-helix in 16-residue model peptides containing only L-amino acids. This review covers several facets of the methodology including spin labeling strategy, interpretation of ESR spectra and the influence of molecular dynamics on the spectral line shapes. Also covered are recent findings of a length-dependent 3(10)-helix-->alpha-helix transition and the role of Arg+ in the stabilization of specific helix structures.

Peptides in membranes: helicity and hydrophobicity

Synthetic model membrane-interactive peptides--both of natural and designed sequence--have become convenient and systematic tools for determination of how the membrane-spanning segments within integral membrane proteins confer protein structure and biology. Conformational studies on these peptides demonstrate that the alpha-helix is the natural choice of conformation for a peptide segment in a membrane, and that a helical conformation will arise "automatically" in a peptide above a threshold hydrophobicity that allows it to associate stably with the membrane. Environmental and sequential contexts thus impart conformational versatility to many of the amino acids, thereby providing a mechanism for producing the diverse structural and functional properties of proteins.

Stabilizing and destabilizing effects of placing beta-branched amino acids in protein alpha-helices

In order to gain greater insight into the effects of beta-branched amino acids on protein alpha-helices, hydrophobic amino acids with varying degrees of beta-branching, including the fully beta-substituted L-2-amino-3,3-dimethylbutanoic acid (ADBA), were incorporated into the protein T4 lysozyme. The unnatural and natural amino acids were substituted at two solvent-exposed alpha-helical sites, Ser 44 and Asn 68, in the protein using the technique of unnatural amino acid mutagenesis. The stabilities of the mutant proteins were determined by using a heat of inactivation assay and from their circular dichroism thermal denaturation curves. Surprisingly, while substitution of the amino acid with the greatest degree of beta-branching, ADBA, destabilizes the protein by 2.5 +/- 0.1 degrees C (0.69 +/- 0.03 kcal/mol) relative to Ala at site 44, the same substitution stabilizes the protein by 1.0 +/- 0.1 degree C (0.27 +/- 0.03 kcal/mol) at site 68. The difference observed at these two positions illustrates the extent to which the local context can mediate the impact of a particular mutation. Molecular dynamics simulations were carried out in parallel to model the structures of the mutant proteins and to examine the energetic consequences of incorporating ADBA. Together, these results suggest that the conformationally restricted beta-branched amino acids are destabilizing, in part, because the beta-branched methyl groups can cause distortions in the local helix backbone. In addition, it is proposed that in some contexts the conformational rigidity of beta-branched amino acids may be stabilizing because it lowers the entropic cost of forming favorable side-chain van der Waals interactions.

Sequence determinants of the capping box, a stabilizing motif at the N-termini of alpha-helices

The capping box, a recurrent hydrogen bonded motif at the N-termini of alpha-helices, caps 2 of the initial 4 backbone amide hydrogen donors of the helix (Harper ET, Rose GD, 1993, Biochemistry 32:7605-7609). In detail, the side chain of the first helical residue forms a hydrogen bond with the backbone of the fourth helical residue and, reciprocally, the side chain of the fourth residue forms a hydrogen bond with the backbone of the first residue. We now enlarge the earlier definition of this motif to include an accompanying hydrophobic interaction between residues that bracket the capping box sequence on either side. The expanded box motif--in which 2 hydrogen bonds and a hydrophobic interaction are localized within 6 consecutive residues--resembles a glycine-based capping motif found at helix C-termini (Aurora R, Srinivasan R, Rose GD, 1994, Science 264:1126-1130).