Conformations of Unsolvated Glycine-Based Peptides

Helical Peptoid Ions in the Gas Phase: Thwarting the Charge Solvation Effect by H-Bond Compensation

Folding and unfolding processes are key aspects that should be mastered for the design of foldamer molecules for targeted applications. In contrast to the solution phase, in vacuo conditions represent a well-defined environment to analyze the intramolecular interactions that largely control the folding/unfolding dynamics. Ion mobility mass spectrometry coupled to theoretical modeling represents an efficient method to decipher the spatial structures of gaseous ions, including foldamers. However, charge solvation typically compacts the ion structure in the absence of strong stabilizing secondary interactions. This is the case in peptoids that are synthetic peptide regioisomers whose side chains are connected to the nitrogen atoms of the backbone instead of α-carbon as in peptides, thus implying the absence of H-bonds among the core units of the backbone. A recent work indeed reported that helical peptoids based on Nspe units formed in solution do not retain their secondary structure when transferred to the gas phase upon electrospray ionization (ESI). In this context, we demonstrate here that the helical structure of peptoids bearing (S)-N-(1-carboxy-2-phenylethyl) bulky side chains (Nscp) is largely preserved in the gas phase by the creation of a hydrogen bond network, induced by the presence of carboxylic moieties, that compensates for the charge solvation process.

Gas-phase acidities of cysteine-polyglycine peptides: the effect of the cysteine position

The sequence and conformational effects on the gas-phase acidities of peptides have been studied by using two pairs of isomeric cysteine-polyglycine peptides, CysGly(3,4)NH(2) and Gly(3,4)CysNH(2). The extended Cooks kinetic method was employed to determine the gas-phase acidities using a triple quadrupole mass spectrometer with an electrospray ionization source. The ion activation was achieved via collision-induced dissociation experiments. The deprotonation enthalpies (Delta(acid)H) were determined to be 323.9 +/- 2.5 kcal/mol (CysGly(3)NH(2)), 319.2 +/- 2.3 kcal/mol (CysGly(4)NH(2)), 333.8 +/- 2.1 kcal/mol (Gly(3)CysNH(2)), and 321.9 +/- 2.8 kcal/mol (Gly(4)CysNH(2)), respectively. The corresponding deprotonation entropies (Delta(acid)S) of the peptides were estimated. The gas-phase acidities (Delta(acid)G) were derived to be 318.4 +/- 2.5 kcal/mol (CysGly(3)NH(2)), 314.9 +/- 2.3 kcal/mol (CysGly(4)NH(2)), 327.5 +/- 2.1 kcal/mol (Gly(3)CysNH(2)), and 317.4 +/- 2.8 kcal/mol (Gly(4)CysNH(2)), respectively. Conformations and energetic information of the neutral and anionic peptides were calculated through simulated annealing (Tripos), geometry optimization (AM1), and single point energy calculations (B3LYP/6-31+G(d)), respectively. Both neutral and deprotonated peptides adopt many possible conformations of similar energies. All neutral peptides are mainly random coils. The two C-cysteine anionic peptides, Gly(3,4)(Cys-H)(-)NH(2), are also random coils. The two N-cysteine anionic peptides, (Cys-H)(-)Gly(3,4)NH(2), may exist in both random coils and stretched helices. The two N-cysteine peptides, CysGly(3)NH(2) and CysGly(4)NH(2), are significantly more acidic than the corresponding C-terminal cysteine ones, Gly(3)CysNH(2) and Gly(4)CysNH(2). The stronger acidities of the former may come from the greater stability of the thiolate anion resulting from the interaction with the helix-macrodipole, in addition to the hydrogen bonding interactions.

Folding and unfolding of helix-turn-helix motifs in the gas phase

Ion mobility measurements and molecular dynamic simulations have been performed for a series of peptides designed to have helix-turn-helix motifs. For peptides with two helical sections linked by a short loop region: AcA(14)KG(3)A(14)K+2H(+), AcA(14)KG(5)A(14)K+2H(+), AcA(14)KG(7)A(14)K+2H(+), and AcA(14)KSar(3)A(14)K+2H(+) (Ac = acetyl, A = alanine, G = glycine, Sar = sarcosine and K = lysine); a coiled-coil geometry with two anti-parallel helices is the lowest energy conformation. The helices uncouple and the coiled-coil unfolds as the temperature is raised. Equilibrium constants determined as a function of temperature yield enthalpy and entropy changes for the unfolding of the coiled-coil. The enthalpy and entropy changes depend on the length and nature of the loop region. For a peptide with three helical sections: protonated AcA(14)KG(3)A(14)KG(3)A(14)K; a coiled-coil bundle with three helices side-by-side is substantially less stable than a geometry with two helices in an antiparallel coiled-coil and the third helix collinear with one of the other two.

Helices and Sheets in vacuo

The structures and properties of unsolvated peptides large enough to possess secondary structure have been examined by experiments and simulations. Some of the factors that stabilize unsolvated helices and sheets have been identified. The charge, in particular, plays a critical role in stabilizing alpha-helices and destabilizing beta-sheets. Some helices are much more stable in vacuum than in aqueous solution. Factors like helix propensity, context, and the incorporation of specific stabilizing interactions have been examined. The helix propensities in vacuum differ from those found in solution. Studies of the hydration of unsolvated peptides can be performed one water molecule at a time. The first few water molecules only bind weakly to unsolvated peptides, and they bind much more strongly to some conformations than to others. The most favorable binding locations are not the protonation sites, but clefts or pockets where a water molecule can establish a network of hydrogen bonds. Non-covalent interactions between secondary structure elements leads to the formation of tertiary structure. Helical peptides assemble into complexes with a variety of intriguing structures. The intramolecular coupling of helices to make antiparallel coiled-coil geometries has also been investigated with model peptides.

Left-Handed and Ambidextrous Helices in the Gas Phase

We have used ion mobility mass spectrometry to study the effect of d-residues on helix formation in unsolvated alanine-based peptides. The right-handed helix of AC-A15K + H+ is significantly disrupted when five or more of the natural L-residues are randomly replaced with D-residues. On the other hand, when a block of L-residues is replaced with D-residues, an unusual ambidextrous structure with helical segments of opposite chirality is formed. A peptide with all D-residues forms a left-handed helix.

Non-Covalent Interactions between Unsolvated Peptides: Helical Complexes Based on Acid−Base Interactions

We have used ion-mobility mass spectrometry to examine the conformations of the protonated complex formed between AcA(7)KA(6)KK and AcEA(7)EA(7), helical alanine-based peptides that incorporate glutamic acid (E) and lysine (K). Designed interactions between the acidic E and basic K residues help to stabilize the complex, which is generated by electrospray and studied in the gas phase. There are two main conformations: (1) a coaxial linear arrangement where the helices are tethered together by an EKK interaction between the pair of lysines at the C-terminus of the AcA(7)KA(6)KK peptide and a glutamic acid at the N-terminus of the AcEA(7)EA(7) peptide and (2) a coiled-coil arrangement with side-by-side antiparallel helices where there is an additional EK interaction between the E and K residues in the middle of the helices. The coiled-coil opens up to the coaxial linear structure as the temperature is raised. Entropy and enthalpy changes for the opening of the coiled-coil were derived from the measurements. The enthalpy change indicates that the interaction between the E and K residues in the middle of the helices is a weak neutral hydrogen bond. The EKK interaction is significantly stronger.

π-Helix Preference in Unsolvated Peptides

Ion mobility measurements have been used to examine helix formations in the gas phase for a series of alanine/glycine-based peptides that incorporate a glutamic acid (E) and lysine (K) at various positions along the backbone. Incorporation of an EK pair lowers the percent helix for all positions (presumably because hydrogen bonding between the backbone and the E and K side chains stabilize the nonhelical globular conformations). The largest percent helix is found when the EK pair is in an i,i+5 arrangement, which suggests that the preferred helical conformation for these peptides is a pi-helix. This conclusion is supported by comparison of cross sections deduced from the ion-mobility measurements to average cross sections calculated for conformations obtained from molecular dynamics simulations. The glutamic acid and lysine may form an ion pair that is stabilized by interactions with the helix macro-dipole.

Probing helix formation in unsolvated peptides

Ion mobility measurements have been used to examine helix formation in unsolvated glycine-based peptides containing three alanine residues. Nine sequence isomers of Ac-[12G3A]K+H(+) were studied (Ac = acetyl, G = glycine, A = alanine, and K = lysine). The amount of helix present for each peptide was examined using two metrics, and it is strongly dependent on the proximity and the location of the alanine residues. Peptides with three adjacent alanines have the highest helix abundances, and those with well-separated alanines have the lowest. The helix abundances for most of the peptides can be fit reasonably well using a modified Lifson-Roig theory. However, Lifson-Roig theory fails to account for several key features of the experimental results. The most likely explanation for the correlation between helix abundances and the number of adjacent alanines is that neighboring alanines promote helix nucleation.

The SAAP force field. A simple approach to a new all-atom protein force field by using single amino acid potential (SAAP) functions in various solvents

A simple strategy to compose a new all-atom protein force field (named as the SAAP force field), which utilizes the single amino acid potential (SAAP) functions obtained in various solvents by ab initio molecular orbital calculation applying the isodensity polarizable continuum model (IPCM), is presented. We considered that the total energy function of a protein force field (E(TOTAL)) is divided into three components; a single amino acid potential term (E(SAAP)), an interamino acid nonbonded interaction term (E(INTER)), and a miscellaneous term (E(OTHERS)), which is ignored (or considered to be constant) at the current version of the force field. The E(INTER) term consists of electrostatic interactions (E(ES')) and van der Waals interactions (E(LJ')). Despite simplicity, the SAAP force field implicitly involves the correlation among individual terms of the Lifson's potential function within a single amino acid unit and can treat solvent effects unambiguously by choosing the SAAP function in an appropriate solvent and the dielectric constant (D) of medium. Application of the SAAP force field to the Monte Carlo simulation of For-Ala(2)-NH(2) in vacuo reasonably reproduced the results of the extensive conformational search by ab initio molecular orbital calculation. In addition, the preliminary Monte Carlo simulations for For-Gly(10)-NH(2) and For-Ala(10)-NH(2) showed reversible transitions from the extended to the pseudosecondary structures in water (D = 78.39) as well as in ether (D = 4.335). The result suggested that the new approach is efficient for fast modeling of protein structures in various environments. Decomposition analysis of the total energy function (E(TOTAL)) by using the SAAP force field suggested that conformational propensities of single amino acids (i.e., the E(SAAP) term) may play definitive roles on the topologies of protein secondary structures.

Helix-turn-helix motifs in unsolvated peptides

The conformations of unsolvated Ac-A14KG3A14K + 2H+ (Ac = acetyl, A = alanine, K = lysine, G = glycine) have been examined by ion mobility measurements and molecular dynamics simulations. This peptide was designed as a model helix-turn-helix motif. It was found to adopt three distinct geometries which were assigned to an extended helical conformation which is only stable at low temperatures (<230 K), a relatively high energy but metastable structure with exchanged lysines, and a coiled-coil. The coiled coil (which consists of an antiparallel arrangement of two helical alanine sections linked by a flexible glycine loop) is the dominant conformation. For temperatures >350 K, the experimental results indicate the helices uncouple and the loop randomizes. From equilibrium constants determined for this helix coupling right arrow over left arrow uncoupling transition, we found DeltaH degrees = -45 kJ mol-1 and DeltaS degrees = 114 J K-1 mol-1. -DeltaH degrees is essentially the enthalpy change for docking the two helices together while DeltaS degrees is essentially the entropy change for freeing up the glycine loop.

The energy landscape of unsolvated peptides: the role of context in the stability of alanine/glycine helices

Ion mobility measurements have been used to examine the conformations present for unsolvated Ac-(AG)(7)A+H(+) and (AG)(7)A+H(+) peptides (Ac = acetyl, A = alanine, and G = glycine) over a broad temperature range (100-410 K). The results are compared to those recently reported for Ac-A(4)G(7)A(4)+H(+) and A(4)G(7)A(4)+H(+), which have the same compositions but different sequences. Ac-(AG)(7)A+H(+) shows less conformational diversity than Ac-A(4)G(7)A(4)+H(+); it is much less helical than Ac-A(4)G(7)A(4)+H(+) at the upper end of the temperature range studied, and at low temperatures, one of the two Ac-A(4)G(7)A(4)+H(+) features assigned to helical conformations is missing for Ac-(AG)(7)A+H(+). Molecular dynamics simulations suggest that the different conformational preferences are not due to differences in the stabilities of the helical states, but differences in the nonhelical states: it appears that Ac-(AG)(7)A+H(+) is more flexible and able to adopt lower energy globular conformations (compact random looking three-dimensional structures) than Ac-A(4)G(7)A(4)+H(+). The helix to globule transition that occurs for Ac-(AG)(7)A+H(+) at around 250-350 K is not a direct (two-state) process, but a creeping transition that takes place through at least one and probably several intermediates.

The energy landscape of unsolvated peptides: helix formation and cold denaturation in Ac-A4G7A4 + H+

Ion mobility measurements and molecular dynamics simulations were performed for unsolvated A4G7A4 + H+ and Ac-A4G7A4 + H+ (Ac = acetyl, A = alanine, G = glycine) peptides. As expected, A4G7A4 + H+ adopts a globular conformation (a compact, random-looking, three-dimensional structure) over the entire temperature range examined (100-410 K). Ac-A4G7A4 + H+ on the other hand is designed to have a flat energy landscape with a marginally stable helical state. This peptide shows at least four different conformations at low temperatures (<230 K). The two conformations with the largest cross sections are attributed to - and partial -helices, while the one with the smallest cross section is globular. The other main conformation may be partially helical. Ac-A4G7A4 + H+ becomes predominantly globular at intermediate temperatures and then becomes more helical as the temperature is raised further. This unexpected behavior may be due to the helix having a higher vibrational entropy than the globular state, as predicted by some recent calculations (Ma, B.; Tsai, C.-J.; Nussinov, R. Biophys. J. 2000, 79, 2739-2753).

Helix unfolding in unsolvated peptides

The conformations of unsolvated Ac-K(AGG)(5)+H(+) and Ac-(AGG)(5)K+H(+) peptides (Ac = acetyl, A = alanine, G = glycine, and K = lysine) have been examined by ion mobility measurements over a wide temperature range (150-410 K). The Ac-K(AGG)(5)+H(+) peptide remains a globule (a compact, roughly spherical structure) over the entire temperature range, while both an alpha-helix and a globule are found for Ac-(AGG)(5)K+H(+) at low temperature. As the temperature is raised the alpha-helix unfolds. Rate constants for loss of the helix (on a millisecond time scale) have been determined as a function of temperature and yield an Arrhenius activation energy and preexponential factor of 38.2 +/- 1.0 kJ mol(-1) and 6.5 +/- 3.7 x 10(9) s(-1), respectively. The alpha-helix apparently does not unfold directly into the globule, but first converts into a long-lived intermediate which survives to a significantly higher temperature before converting. According to molecular dynamics simulations, there is a partially untwisted helical conformation that has both a low energy and a well-defined geometry. This special structure lies between the helix and globule and may be the long-lived intermediate.

A theoretical study on the origin of cooperativity in the formation of 3(10)- and alpha-helices

By using a simple repeating unit method, we have conducted a theoretical study which delineates the preferences for beta-strand, 2(7)-ribbon, 3(10)-helix, and alpha-helix formation for a series of polyglycine models up to 14 amino acid residues (Ac-(Gly)(n), n = 0, 1, 2,., 14). Interactions among residues, which result in cooperativity, are clearly indicated by variations in calculated energies of the residues. Whereas no cooperativity is found in the formation of beta-strands and 2(7)-ribbons, there is a significant cooperativity in the formation of 3(10)- and alpha-helices, especially for the latter. In the case of alpha-helices, the 14th residue is more stable than the 3rd by about 3 kcal/mol. A good correlation between calculated residue energy and residue dipole moment was uncovered, indicating the importance of long-range electrostatic interactions to the cooperativity. The results of our calculations are compared with those of the AMBER and PM3 methods, and indicate that both methods, AMBER and PM3, need further development in the cooperative view of electrostatic interactions. The result should be of importance in providing insight into protein folding and formation of helical structures in a variety of polymeric compounds. This also suggests a strategy for the development of more consistent molecular mechanics force fields.

Peptides and proteins in the vapor phase

This article provides a review of recent studies of the properties of unsolvated (and partially solvated) peptides and proteins. The methods used to produce vapor-phase peptide and protein ions are described along with some of the techniques used to study them, such as H/D exchange, blackbody infrared radiative dissociation, and ion mobility measurements. Studies of unsolvated peptides and proteins provide information about their intrinsic intramolecular interactions. The topics covered include the role of zwitterions and salt bridges in the vapor phase, Coulomb interactions in multiply charged ions, the unfolding and refolding of vapor-phase proteins, and the stability of unsolvated helices and sheets. Finally, dehydration and rehydration studies of proteins in the vapor phase are described. These can provide exquisitely detailed information about hydration interactions, such as the enthalpy and entropy changes associated with adsorbing individual water molecules.