Physical reasons for the unusual α-helix stabilization afforded by charged or neutral polar residues in alanine-rich peptides

We have carried out conformational energy calculations on alanine-based copolymers with the sequence Ac-AAAAAXAAAA-NH2 in water, where X stands for lysine or glutamine, to identify the underlying source of stability of alanine-based polypeptides containing charged or highly soluble polar residues in the absence of charge–charge interactions. The results indicate that ionizable or neutral polar residues introduced into the sequence to make them soluble sequester the water away from the CO and NH groups of the backbone, thereby enabling them to form internal hydrogen bonds. This solvation effect dictates the conformational preference and, hence, modifies the conformational propensity of alanine residues. Even though we carried out simulations for specific amino acid sequences, our results provide an understanding of some of the basic principles that govern the process of folding of these short sequences independently of the kind of residues introduced to make them soluble. In addition, we have investigated through simulations the effect of the bulk dielectric constant on the conformational preferences of these peptides. Extensive conformational Monte Carlo searches on terminally blocked 10-mer and 16-mer homopolymers of alanine in the absence of salt were carried out assuming values for the dielectric constant of the solvent ɛ of 80, 40, and 2. Our simulations show a clear tendency of these oligopeptides to augment the α-helix content as the bulk dielectric constant of the solvent is lowered. This behavior is due mainly to a loss of exposure of the CO and NH groups to the aqueous solvent. Experimental evidence indicates that the helical propensity of the amino acids in water shows a dramatic increase on addition of certain alcohols, such us trifluoroethanol. Our results provide a possible explanation of the mechanism by which alcohol/water mixtures affect the free energy of helical alanine oligopeptides relative to nonhelical ones.

A host–guest system to study structure–function relationships of membrane fusion peptides

We designed a host–guest fusion peptide system, which is completely soluble in water and has a high affinity for biological and lipid model membranes. The guest sequences are those of the fusion peptides of influenza hemagglutinin, which are solubilized by a highly charged unstructured C-terminal host sequence. These peptides partition to the surface of negatively charged liposomes or erythrocytes and elicit membrane fusion or hemolysis. They undergo a conformational change from random coil to an obliquely inserted (≈33° from the surface) α-helix on binding to model membranes. Partition coefficients for membrane insertion were measured for influenza fusion peptides of increasing lengths (n = 8, 13, 16, and 20). The hydrophobic contribution to the free energy of binding of the 20-residue fusion peptide at pH 5.0 is −7.6 kcal/mol (1 cal = 4.18 J). This energy is sufficient to stabilize a “stalk” intermediate if a typical number of fusion peptides assemble at the site of membrane fusion. The fusion activity of the fusion peptides increases with each increment in length, and this increase strictly correlates with the hydrophobic binding energy and the angle of insertion.

A meanfield approach to the thermodynamics of a protein-solvent system with application to the oligomerization of the tumor suppressor p53

The thermodynamic stability and oligomerization status of the tumor suppressor p53 tetramerization domain have been studied experimentally and theoretically. A series of hydrophilic mutations at Met-340 and Leu-344 of human p53 were designed to disrupt the hydrophobic dimer–dimer interface of the tetrameric oligomerization domain of p53 (residues 325–355). Meanfield calculations of the free energy of the solvated mutants as a function of interdimer distance were compared with experimental data on the thermal stability and oligomeric state (tetramer, dimer, or equilibrium mixture of both) of each mutant. The calculations predicted a decreasing stability and oligomeric state for the following amino acids at residue 340: Met (tetramer) > Ser Asp, His, Gln, > Glu, Lys (dimer), whereas the experimental results showed the following order: Met (tetramer) > Ser > Gln > His, Lys > Asp, Glu (dimers). For residue 344, the calculated trend was Leu (tetramer) > Ala > Arg, Gln, Lys (dimer), and the experimental trend was Leu (tetramer) > Ala, Arg, Gln, Lys (dimer). The discrepancy for the lysine side chain at residue 340 is attributed to the dual nature of lysine, both hydrophobic and charged. The incorrect prediction of stability of the mutant with Asp at residue 340 is attributed to the fact that within the meanfield approach, we use the wild-type backbone configuration for all mutants, but low melting temperatures suggest a softening of the α-helices at the dimer–dimer interface. Overall, this initial application of meanfield theory toward a protein-solvent system is encouraging for the application of the theoretical model to more complex systems.

Fibroblast growth factor (FGF) homologous factors: new members of the FGF family implicated in nervous system development.

Four new members of the fibroblast growth factor (FGF) family, referred to as fibroblast growth factor homologous factors (FHFs), have been identified by a combination of random cDNA sequencing, data base searches, and degenerate PCR. Pairwise comparisons between the four FHFs show between 58% and 71% amino acid sequence identity, but each FHF shows less than 30% identity when compared with other FGFs. Like FGF-1 (acidic FGF) and FGF-2 (basic FGF), the FHFs lack a classical signal sequence and contain clusters of basic residues that can act as nuclear localization signals. In transiently transfected 293 cells FHF-1 accumulates in the nucleus and is not secreted. Each FHF is expressed in the developing and adult nervous systems, suggesting a role for this branch of the FGF family in nervous system development and function.

Direct measurement of salt-bridge solvation energies using a peptide model system: implications for protein stability.

The solvation energies of salt bridges formed between the terminal carboxyl of the host pentapeptide AcWL- X-LL and the side chains of Arg or Lys in the guest (X) position have been measured. The energies were derived from octanol-to-buffer transfer free energies determined between pH 1 and pH 9. 13C NMR measurements show that the salt bridges form in the octanol phase, but not in the buffer phase, when the side chains and the terminal carboxyl group are charged. The free energy of salt-bridge formation in octanol is approximately -4 kcal/mol (1 cal = 4.184 J), which is equal to or slightly larger than the sum of the solvation energies of noninteracting pairs of charged side chains. This is about one-half the free energy that would result from replacing a charge pair in octanol with a pair of hydrophobic residues of moderate size. Therefore, salt bridging in octanol can change the favorable aqueous solvation energy of a pair of oppositely charged residues to neutral or slightly unfavorable but cannot provide the same free energy decrease as hydrophobic residues. This is consistent with recent computational and experimental studies of protein stability.