Anticooperativity in a Glu-Lys-Glu salt bridge triplet in an isolated alpha-helical peptide

Salt bridges between oppositely charged side chains are well-known to stabilize protein structure, though their contributions vary considerably. Here we study Glu-Lys and Lys-Glu salt bridges, formed when the residues are spaced i, i + 4 surface of an isolated alpha-helix in aqueous solution. Both are stabilizing by -0.60 and -1.02 kcal/mol, respectively, when the interacting residues are fully charged. When the side chains are spaced i, i + 4, i + 8, forming a Glu-Lys-Glu triplet, the second salt bridge provides no additional stabilization to the helix. We attribute this to the inability of the central Lys to form two salt bridges simultaneously. Analysis of these salt bridges in protein structures shows that the Lys-Glu interaction is dominant, with the side chains of the Glu-Lys pair far apart.     

Stabilization of Helical Peptides by Mixed Spaced Salt Bridges

Abstract

Whether or not surface salt bridges have a strong stabilizing effect on the native structure in proteins remains uncertain. Previous studies of model peptides have shown that salt bridges spaced at i,i+4 along the chain are more stabilizing than those spaced at i,i+3, with a preference for the order acid-base rather than base-acid from N to C terminus. An analysis of the effect of spacing the ion pairs in short helical peptides is presented, in which acidic and basic side chains spaced two or three residues apart alternate along the chain. The mixed spacing proves to be stabilizing relative to pure spacings. A control peptide in which salt bridges were spaced uniformly three residues apart proved to form a β-sheet structure rather than a-helix. This is due to formation of a silk-like apolar face consisting of alanine side chains; the mesoscopic structure formed by these sheets can be imaged by scanning microscopy.     

Folding simulations of alanine-based peptides with lysine residues.

The folding of short alanine-based peptides with different numbers of lysine residues is simulated at constant temperature (274 K) using the rigid-element Monte Carlo method. The solvent-referenced potential has prevented the multiple-minima problem in helix folding. From various initial structures, the peptides with three lysine residues fold into helix-dominated conformations with the calculated average helicity in the range of 60-80%. The peptide with six lysine residues shows only 8-14% helicity. These results agree well with experimental observations. The intramolecular electrostatic interaction of the charged lysine side chains and their electrostatic hydration destabilize the helical conformations of the peptide with six lysine residues, whereas these effects on the peptides with three lysine residues are small. The simulations provide insight into the helix-folding mechanism, including the beta-bend intermediate in helix initiation, the (i, i + 3) hydrogen bonds, the asymmetrical helix propagation, and the asymmetrical helicities in the N- and C-terminal regions. These findings are consistent with previous studies.     

Cooperative helix stabilization by complex Arg-Glu salt bridges

Among the interactions that stabilize the native state of proteins, the role of electrostatic interactions has been difficult to quantify precisely. Surface salt bridges or ion pairs between acidic and basic side chains have only a modest stabilizing effect on the stability of helical peptides or proteins: estimates are roughly 0.5 kcal/mol or less. On the other hand, theoretical arguments and the occurrence of salt bridge networks in thermophilic proteins suggest that multiple salt bridges may exert a stronger stabilizing effect. We show here that triads of charged side chains, Arg(+)-Glu(-)-Arg(+) spaced at i,i+4 or i,i+3 intervals in a helical peptide stabilize alpha helix by more than the additive contribution of two single salt bridges. The free energy of the triad is more than 1 kcal/mol in excess of the sum of the individual pairs, measured in low salt concentration (10 mM). The effect of spacing the three groups is severe; placing the charges at i,i+4 or i,i+3 sites has a strong effect on stability relative to single bridges; other combinations are weaker. A conservative calculation suggests that interactions of this kind between salt bridges can account for much of the stabilization of certain thermophilic proteins.     

Stabilizing nonpolar/polar side-chain interactions in the alpha-helix.

A simplistic, yet often used, view of protein stability is that amino acids attract other amino acids with similar polarity, whereas nonpolar and polar side chains repel. Here we show that nonpolar/polar interactions, namely Val or Ile bonding to Lys or Arg in alpha-helices, can in fact be stabilizing. Residues spaced i, i + 4 in alpha-helices are on the same face of the helix, with potential to favorably interact and stabilize the structure. We observe that the nonpolar/polar pairs Ile-Lys, Ile-Arg, and Val-Lys occur in protein helices more often than expected when spaced i, i + 4. Partially helical peptides containing pairs of nonpolar/polar residues were synthesized. Controls with i, i + 5 spacing have the residues on opposite faces of the helix and are less helical than the test peptides with the i, i + 4 interactions. Experimental circular dichroism results were analyzed with helix-coil theory to calculate the free energy for the interactions. All three stabilize the helix with DeltaG between -0.14 and -0.32 kcal x mol(-1). The interactions are hydrophobic with contacts between Val or Ile and the alkyl groups in Arg or Lys. Side chains such as Lys and Arg can thus interact favorably with both polar and nonpolar residues.     

Contribution of arginine-glutamate salt bridges to helix stability

Peptide side chain interactions were studied by molecular dynamics simulation using explicit solvent on a peptide with the sequence AAARAAAAEAAEAAAARA. Three different protonation states of the glutamic acid side chains were simulated for four 20 ns runs each, a total simulation time of 240 ns. Two different salt bridge geometries were observed and the preferred geometry was found to depend on Glu — Arg residue spacing. Stable charge clusters were also observed, particularly in the fully charged peptide. Salt bridges were selectively interrupted upon protonation, with concomitant changes in secondary structure. The fully charged peptide was highly helical between residues 9 and 13, although protonation increased helicity near the N-terminus. The contribution of salt bridges to helix stability therefore depends on both position and relative position of charged residues within a sequence.     

Helicity of short E-R/K peptides

Understanding the secondary structure of peptides is important in protein folding, enzyme function, and peptide‐based drug design. Previous studies of synthetic Ala‐based peptides (>12 a.a.) have demonstrated the role for charged side chain interactions involving Glu/Lys or Glu/Arg spaced three (i, i + 3) or four (i, i + 4) residues apart. The secondary structure of short peptides (<9 a.a.), however, has not been investigated. In this study, the effect of repetitive Glu/Lys or Glu/Arg side chain interactions, giving rise to E‐R/K helices, on the helicity of short peptides was examined using circular dichroism. Short E‐R/K–based peptides show significant helix content. Peptides containing one or more E‐R interactions display greater helicity than those with similar E‐K interactions. Significant helicity is achieved in Arg‐based E‐R/K peptides eight, six, and five amino acids long. In these short peptides, each additional i + 3 and i + 4 salt bridge has substantial contribution to fractional helix content. The E‐R/K peptides exhibit a strongly linear melt curve indicative of noncooperative folding. The significant helicity of these short peptides with predictable dependence on number, position, and type of side chain interactions makes them an important consideration in peptide design.     

Conformational study of two synthetic peptides with sequence analogies to the N-terminal fragment of RNase A

The conformational properties of two synthetic model peptides, AEAAHAAEAAHMG (PA) and AEAAHAFEAAHMG (PF), have been studied using CD and 1H-NMR methods. In both peptides, glutamate and histidine residues are situated in such a way that two salt bridges between Glu- (i) and His+ (i + 3) can be formed. A salt bridge of this type (Glu- 9-His+ 12) was postulated previously to stabilize, to a great extent, the alpha-helical conformation of isolated N-terminal fragments of RNase A: C-peptide and S-peptide (A. Bierzyński, P.S. Kim and R.L. Baldwin, Proc. Natl. Acad. Sci. U.S.A. 79 (1982) 2470). Although in both PA and PF salt bridges between glutamates and histidines are formed, as demonstrated by the pH-titration curves of the glutamate gamma-proton signals, no traces of helical conformation have been detected. Evidently, the Glu- (i)-His+ (i + 3) salt bridges do not stabilize the alpha-helical conformation. A comparative analysis of PA and PF NMR spectra provides strong evidence that the phenylalanyl ring in PF interacts not only with the hydrophobic methyl groups of almost all alanine residues but also with the histidine rings and the glutamate side chains in their protonated as well as deprotonated forms. Similar interactions, involving Phe 8, can be expected in the N-terminal fragments of RNase and should be taken into account as an important factor determining the conformational properties of C- and S-peptides.     

Importance of context in protein folding: secondary structural propensities versus tertiary contact-assisted secondary structure formation

Molecular dynamics simulations can be used to reveal the detailed conformational behaviors of peptides and proteins. By comparing fragment and full-length protein simulations, we can investigate the role of each peptide segment in the folding process. Here, we take advantage of information regarding the helix formation process from our previous simulations of barnase and protein A as well as new simulations of four helical fragments from these proteins at three different temperatures, starting with both helical and extended structures. Segments with high helical propensity began the folding process by tethering the chain through side chain interactions involving either polar interactions, such as salt bridges, or hydrophobic staples. These tethers were frequently nonnative (i.e., not i --> i + 4 spacing) and provided a scaffold for other residues, thereby limiting the conformational search. The helical structure then propagated on both sides of the tether. Segments with low stability and propensity formed later in the folding process and utilized contacts with other portions of the protein when folding. These helices formed via a tertiary contact-assisted mechanism, primarily via hydrophobic contacts between residues distant in sequence. Thus, segments with different helical propensities appear to play different roles during protein folding. Furthermore, the active role of nonlocal side chains in helix formation highlights why we must move beyond simple hierarchical models of protein folding.     

A mechanism for the evolution of phosphorylation sites

Protein phosphorylation provides a mechanism for the rapid, reversible control of protein function. Phosphorylation adds negative charge to amino acid side chains, and negatively charged amino acids (Asp/Glu) can sometimes mimic the phosphorylated state of a protein. Using a comparative genomics approach, we show that nature also employs this trick in reverse by evolving serine, threonine, and tyrosine phosphorylation sites from Asp/Glu residues. Structures of three proteins where phosphosites evolved from acidic residues (DNA topoisomerase II, enolase, and C-Raf) show that the relevant acidic residues are present in salt bridges with conserved basic residues, and that phosphorylation has the potential to conditionally restore the salt bridges. The evolution of phosphorylation sites from glutamate and aspartate provides a rationale for why phosphorylation sometimes activates proteins, and helps explain the origins of this important and complex process.     

C-terminal anchoring of a peptide to class II MHC via the P10 residue is compatible with a peptide bulge.

The binding of antigenic peptide to class II MHC is mediated by hydrogen bonds between the MHC and the peptide, by salt bridges, and by hydrophobic interactions. The latter are confined to a number of deeper pockets within the peptide binding groove, and peptide side chains that interact with these pockets are referred to as anchor residues. T cell recognition involves solvent-accessible peptide residues along with minor changes in MHC helical pitch induced by the anchor residues. In class I MHC there is an added level of epitope complexity that results from binding of longer peptides that bulge out into the solvent-accessible, T cell contact area. Unlike class I MHC, class II MHC does not bind peptides of discrete length, and the possibility of peptide bulging has not been clearly addressed. A peptide derived from position 24-37 of integrin beta(3) can either bind or not bind to the class II MHC molecule HLA DRB3*0101 based on a polymorphism at the P9 anchor. We show that the loss of binding can be compensated by changes at the P10 position. We propose that this could be an example of a class II peptide bulge. Although not as efficient as P9 anchoring, the use of P10 as an anchor adds another possible mechanism by which T cell epitopes can be generated in the class II presentation system.     

A phosphoserine-lysine salt bridge within an alpha-helical peptide, the strongest alpha-helix side-chain interaction measured to date

Phosphorylation is ubiquitous in control of protein activity, yet its effects on protein structure are poorly understood. Here we investigate the effect of serine phosphorylation in the interior of an alpha-helix when a salt bridge is present between the phosphate group and a positively charged side chain (in this case lysine) at i,i + 4 spacing. The stabilization of the helix is considerable and can overcome the intrinsically low preference of phosphoserine for the interior of the helix. The effect is pH dependent, as both the lysine and phosphate groups are titratable, and so calculations are given for several charge combinations. These results, with our previous work, highlight the different, context-dependent effects of phosphorylation in the alpha-helix. The interaction between the phosphate(2)(-) group and the lysine side chain is the strongest yet recorded in helix-coil studies. The results are of interest both in de novo design of peptides and in understanding the structural modes of control by phosphorylation.     

Surface electrostatic interactions contribute little of stability of barnase

Electrostatic interactions are believed to play an important role in stabilizing the native structure of proteins. We have quantified the contribution to stability of an interaction between two oppositely charged side-chains on the surface of barnase. Using site-directed mutagenesis, glutamate 28 and lysine 32 were introduced onto the solvent-accessible side of the second alpha-helix in barnase. These two residues are separated by one turn of the helix, and so are ideally situated for their opposite charges to interact. Double mutant cycle analysis reveals that the interaction between Glu28 and Lys32 contributes only approximately 0.2 kcal/mol to stability of the protein. All other interactions between exposed charged side-chains in barnase examined so far also contribute little to stability. We explain this low value by their location on the surface, rather than in the interior, of the protein.     

Structural and functional implications of a proline residue in the antimicrobial peptide gaegurin.

Although it is commonly known as a helix breaker, proline residues have been found in the alpha-helical regions of many peptides and proteins. The antimicrobial peptide gaegurin displays alpha-helical structure and has a central proline residue (P14). The structure and activity of gaegurin and its alanine derivative (P14A) were determined by various spectroscopic methods, restrained molecular dynamics, and biological assays. Both P14 and P14A exhibited cooperative helix formation in solution, but the helical stability of P14 was reduced substantially when compared to that of P14A. Chemical-shift analysis indicated that both of the peptides formed curved helices and that P14 showed diminished stability in the region around the central proline. However, hydrogen-exchange data revealed remarkable differences in the location of stable amide protons. P14 showed a stable region in the concave side of the curved helix, while P14A exhibited a stable region in the central turn of the helix. The model structure of P14 exhibited a pronounced kink, in contrast to the uniform helix of P14A. Both peptides showed comparable binding affinities for negatively charged lipids, while P14 had a considerably reduced affinity for a neutral lipid. With its destabilized alpha-helix, P14 exhibited greater antibacterial activity than did P14A. Hence, electrostatic interaction between helical peptides and lipid membranes is believed to be the dominant factor for antibacterial activity. Moreover, helical stability can modulate peptide binding to membranes that is driven by electrostatic interactions. The observation that P14 is a more potent antibacterial agent than P14A implies that the helical kink of P14 plays an important role in the disruption of bacterial membranes.     

Addition of side chain interactions to modified Lifson-Roig helix-coil theory: application to energetics of phenylalanine-methionine interactions.

We introduce here i, i + 3 and i, i + 4 side chain interactions into the modified Lifson-Roig helix-coil theory of Doig et al. (1994, Biochemistry 33:3396-3403). The helix/coil equilibrium is a function of initiation, propagation, capping, and side chain interaction parameters. If each of these parameters is known, the helix content of any isolated peptide can be predicted. The model considers every possible conformation of a peptide, is not limited to peptides with only a single helical segment, and has physically meaningful parameters. We apply the theory to measure the i, i + 4 interaction energies between Phe and Met side chains. Peptides with these residues spaced i, i + 4 are significantly more helical than controls where they are spaced i, i + 5. Application of the model yields delta G for the Phe-Met orientation to be -0.75 kcal.mol-1, whereas that for the Met-Phe orientation is -0.54 kcal.mol-1. These orientational preferences can be explained, in part, by rotamer preferences for the interacting side chains. We place Phe-Met i, i + 4 at the N-terminus, the C-terminus, and in the center of the host peptide. The model quantitatively predicts the observed helix contents using a single parameter for the side chain-side chain interaction energy. This result indicates that the model works well even when the interaction is at different locations in the helix.     

Structural implications of placing cationic residues at either the NH2- or COOH-terminus in a pore-forming synthetic peptide

Restoration of chloride conductance via introduction of an anion-selective pore, formed by a channel-forming peptide, has been hypothesized as a novel treatment modality for patients with cystic fibrosis. Delivery of these peptides from an aqueous environment in the absence of organic solvents is paramount. M2GlyR peptides, designed based on the glycine receptor, insert into lipid bilayers and polarized epithelial cells and assemble spontaneously into chloride-conducting pores. Addition of 4 lysine residues to either terminus increases the solubility of M2GlyR peptides. Both orientations of the helix within the membrane form an anion-selective pore, however, differences in solubility, associations and channel-forming activity are observed. To determine how the positioning of the lysine residues affects these properties, structural characteristics of the lysyl-modified peptides were explored utilizing chemical cross-linking, NMR and molecular modeling. Initial model structures of the a-helical peptides predict that lysine residues at the COOH-terminus form a capping structure by folding back to form hydrogen bonds with backbone carbonyl groups and hydroxyl side chains of residues in the helical segment of the peptide. In contrast, lysine residues at the NH2-terminus form fewer H-bonds and extend away from the helical backbone. Results from NMR and chemical cross-linking support the model structures. The C-cap formed by H-bonding of lysine residues is likely to account for the different biophysical properties observed between NH2- and COOH-terminal-modified M2GlyR peptides.     

Structural requirements for N-trimethylation of lysine 115 of calmodulin

Calmodulin is trimethylated at lysine 115 by a highly specific methyltransferase that utilizes S-adenosylmethionine as a co-substrate. Lysine 115 is found within a highly conserved six-amino acid loop (LGEKLT) that forms a 90 degrees turn between EF-hand III and EF-hand IV in the carboxyl-terminal lobe. In the present work a mutagenesis approach was used to investigate the structural features of the carboxyl-terminal lobe that lead to the specificity of calmodulin methylation. Three structural regions within the carboxyl-terminal lobe appear to be involved in methyltransferase recognition: the highly conserved six-amino acid loop-turn region that contains lysine 115 as well as the adjacent alpha-helices (helix 6 and helix 7) from EF-hands III and IV. Site-directed mutagenesis of residues in the loop show that three residues, glycine 113, glutamate 114, and leucine 116 are essential for methylation. In addition, subdomain (individual helix or Ca(2+) binding loop) exchange mutants show that the substitutions of either helix 6 (EF-hand III) with helix 2 (EF-hand I) or helix 7 (EF-hand IV) with helix 3 (EF-hand II) compromises methylation. Charge-to-alanine mutations in helix 7 show that substitution of conserved charged residues at positions 118, 120, 122, 126, and 127 reduced lysine 115 methylation rates, suggesting possible electrostatic interactions between this helix and the methyltransferase. Single substitutions in helix 6 did not affect calmodulin methylation, suggesting this region may play a more indirect role in stabilizing the conformation of the methyltransferase recognition sequence.     

Effects of charged amino acids at b and c heptad positions on specificity and stability of four-chain coiled coils

An understanding of the balance of chemical forces responsible for protein stability and specificity of structure is essential for the success of efforts in protein design. Specifically, electrostatic interactions between charged amino acids have been explored extensively to understand the contribution of this force to protein stability. Much research on the importance of electrostatic interactions as specificity and stability determinants in two-stranded coiled coils has been done, but there remains significant controversy about the magnitude of the attractive forces using such systems. We have developed a four-stranded coiled-coil system with charged residues incorporated at b and c heptad positions to explore the role of charge interactions. Here, we test quantitatively the effects of varying sidechain length on the magnitude of such electrostatic interactions. We synthesized peptides containing either aspartate or ornithine at both b and c heptad positions and tested their ability to self-associate and to hetero-associate with one another and with peptides containing glutamate or lysine at the same positions. We find that interactions between glutamate and either lysine or ornithine are more favorable than the corresponding interactions involving aspartate. In each case, charged interactions provide additional stability to coiled coils, although helix propensity effects may play a significant role in determining the overall stability of these structures.     

The stability of thermophilic proteins: a study based on comprehensive genome comparison

We address the question of the thermal stability of proteins in thermophiles through comprehensive genome comparison, focussing on the occurrence of salt bridges. We compared a set of 12 genomes (from four thermophilic archaeons, one eukaryote, six mesophilic eubacteria, and one thermophilic eubacteria). Our results showed that thermophiles have a greater content of charged residues than mesophiles, both at the overall genomic level and in alpha helices. Furthermore, we found that in thermophiles the charged residues in helices tend to be preferentially arranged with a 1–4 helical spacing and oriented so that intra-helical charge pairs agree with the helix dipole. Collectively, these results imply that intra-helical salt bridges are more prevalent in thermophiles than mesophiles and thus suggest that they are an important factor stabilizing thermophilic proteins. We also found that the proteins in thermophiles appear to be somewhat shorter than those in mesophiles. However, this later observation may have more to do with evolutionary relationships than with physically stabilizing factors. In all our statistics we were careful to controls for various biases. These could have, for instance, arisen due to repetitive or duplicated sequences. In particular, we repeated our calculation using a variety of random and directed sampling schemes. One of these involved making a "stratified sample," a representative cross-section of the genomes derived from a set of 52 orthologous proteins present roughly once in each genome. For another sample, we focused on the subset of the 52 orthologs that had a known 3D structure. This allowed us to determine the frequency of tertiary as well as main-chain salt bridges. Our statistical controls supported our overall conclusion about the prevalence of salt bridges in thermophiles in comparison to mesophiles. Electronic Publication     

Solution structures of cyclic and dicyclic analogues of growth hormone releasing factor as determined by two-dimensional NMR and CD spectroscopies and constrained molecular dynamics.

Solution structures were determined for a linear analogue of growth hormone releasing factor (GRF), and cyclic and dicyclic analogues in which the side chains of aspartyl and lysyl residues spaced at positions i-(i + 4) were joined to form a lactam. The four analogues were [Ala15]-GRF-(1-29)-NH2 and its cyclo8-12, cyclo21-25, and dicyclo8-12;21-25 derivatives. The peptides were studied in two solvent systems: 75% methanol/25% water at pH 6.0; and 100% water at pH 3.0. CD spectroscopy was used to assess the overall alpha-helical content. Nuclear magnetic resonance spectroscopy was used to determine the structures in more detail. Nearly complete proton resonance assignments were made for each of the peptides, in both solvents. Nuclear Overhauser effects were converted into distance constraints and applied in the molecular dynamics program CHARMM to evaluate the range of low-energy structures that satisfied the nmr data. In 75% methanol, all of the peptides are comprised of a single alpha-helical segment with fraying of one to three residues at each end. The linear analogue has a tendency to kink. In water, the analogues have two helical segments with flexible regions between them and at the termini of the peptides. The linear analogue is helical at residues 7-14 and 21-28. In the cyclo8-12 analogue, the N-terminal helical region extends to include residues 7-19, while the other helical region is slightly shortened. In the cyclo21-25 analogue, the C-terminal helical region is extended to include residues 19-28, while the N-terminal helical region is destabilized. The dicyclic analogue has the largest N-terminal helix, spanning residues 7-20, but its helical segment at residues 21-28 is not well ordered. All of the analogues exhibit substantial biological activity. The cyclic and dicyclic analogues show dramatically increased resistance to degradation during incubation with human plasma. The i-(i + 4) lactam, therefore, appears to be a synthetic means of stabilizing a local alpha-helical conformation, which may be of general use in the design of active, stable peptides.