Determination of the gas-phase basicities of proline and its di- and tripeptides with glycine: the enhanced basicity of prolylproline

Weak Acid-Base Interactions of Histidine and Cysteine Affect the Charge States, Tertiary Structure, and Zn(II)-Binding of Heptapeptides

Zinc fingers are proteins that are characterized by the coordination of zinc ions by an amino acid sequence that commonly contains two histidines and two cysteines (2His-2Cys motif). Investigations of oligopeptides that contain the 2His-2Cys motif, e.g., acetyl-His₁-Cys₂-Gly₃-Pro₄-Tyr₅-His₆-Cys₇, have discovered they exhibit pH-dependent Zn(II) chelation and have redox activities with Cu(I/II), forming a variety of metal complexes. To further understand how these 2His-2Cys oligopeptides bind these metal ions, we have undertaken a series of ion mobility-mass spectrometry and B3LYP/LanL2DZ computational studies of structurally related heptapeptides. Starting with the sequence above, we have modified the potential His, Cys, or C-terminus binding sites and report how these changes in primary structure affect the oligopeptides positive and negative charge states, conformational structure, collision-induced breakdown energies, and how effectively Zn(II) binds to these sequences. The results show evidence that the weak acid-base properties of Cys-His are intrinsically linked and can result in an intramolecular salt-bridged network that affects the oligopeptide properties.

Energetic switch of the proline effect in collision-induced dissociation of singly and doubly protonated peptide Ala-Ala-Arg-Pro-Ala-Ala

Suppression of the selective cleavage at N-terminal of proline is observed in the peptide cleavage by proteolytic enzyme trypsin and in the fragment ion mass spectra of peptides containing Arg-Pro sequence. An insight into the fragmentation mechanism of the influence of arginine residue on the proline effect can help in prediction of mass spectra and in protein structure analysis. In this work, collision-induced dissociation spectra of singly and doubly charged peptide AARPAA were studied by ESI MS/MS and theoretical calculation methods. The proline effect was evaluated by comparing the experimental ratio of fragments originated from cleavage of different amide bonds. The results revealed that the backbone amide bond cleavage was selected by the energy barrier height of the fragmentation pathway although the strong proton affinity of the Arg side chain affected the stereostructure of the peptide and the dissociation mechanism. The thermodynamic stability of the fragment ions played a secondary role in the abundance ratio of fragments generated via different pathways. Fragmentation studies of protonated peptide AACitPAA supported the energy-dependent hypothesis. The results provide an explanation to the long-term arguments between the steric conflict and the proton mobility mechanisms of proline effect.

Electron transfer dissociation mass spectrometry of acidic phosphorylated peptides cationized with trivalent praseodymium

The lanthanide ion praseodymium, Pr(III), was employed to study metallated ion formation and electron transfer dissociation (ETD) of 27 biological and model highly acidic phosphopeptides. All phosphopeptides investigated form metallated ions by electrospray ionization (ESI) that can be studied by ETD to yield abundant sequence information. The ions formed are [M + Pr - H]²⁺ , [M + Pr]³⁺ , and [M + Pr + H]⁴⁺ . All biological phosphopeptides with a chain length of seven or more residues generate [M + Pr]³⁺ . For biological phosphopeptides, [M + Pr]³⁺ undergoes more backbone cleavage by ETD than [M + Pr - H]²⁺ and, in some cases, full sequence coverage occurs. Acidic model phosphorylated hexa-peptides and octa-peptides, composed of alanine residues and one phosphorylated residue, form exclusively [M + Pr - H]²⁺ by ESI. Limited sequence information is obtained by ETD of [M + Pr - H]²⁺ with only metallated product ions being generated. For two biological phosphopeptides, [M + Pr + H]⁴⁺ is observed and may be due to the presence of at least one residue with a highly basic side chain that facilitates the addition of an extra proton. For the model phosphopeptides, more sequence coverage occurs when the phosphorylated residue is in the middle of the sequence than at either the N- or C-terminus. ETD of the metallated precursor ions formed by ESI generates exclusively metallated and nonmetallated c- and z-ions for the biological phosphopeptides, while metallated c-ions, z-ions, and a few y-ions form for the model phosphopeptides. Most of the product ions contain the phosphorylated residue indicating that the metal ion binds predominantly at the deprotonated phosphate group. The results of this study indicate that ETD is a promising tool for sequencing highly acidic phosphorylated peptides by metal adduction with Pr (III) and, by extension, all nonradioactive lanthanide metal ions.

Effects of acidic peptide size and sequence on trivalent praseodymium adduction and electron transfer dissociation mass spectrometry

Using the lanthanide ion praseodymium, Pr(III), metallated ion formation and electron transfer dissociation (ETD) were studied for 25 biological and model acidic peptides. For chain lengths of seven or more residues, even highly acidic peptides that can be difficult to protonate by electrospray ionization will metallate and undergo abundant ETD fragmentation. Peptides composed of predominantly acidic residues form only the deprotonated ion, [M + Pr - H]2+ ; this ion yields near complete ETD sequence coverage for larger peptides. Peptides with a mixture of acidic and neutral residues generate [M + Pr]3+ , which cleaves between every residue for many peptides. Acidic peptides that contain at least one residue with a basic side chain also produce the protonated ion, [M + Pr + H]4+ ; this ion undergoes the most extensive sequence coverage by ETD. Primarily metallated and non-metallated c- and z-ions form for all peptides investigated. Metal adducted product ions are only present when at least half of the peptide sequence can be incorporated into the ion; this suggests that the metal ion simultaneously attaches to more than one acidic site. The only site consistently lacking dissociation is at the N-terminal side of a proline residue. Increasing peptide chain length generates more backbone cleavage for metal-peptide complexes with the same charge state. For acidic peptides with the same length, increasing the precursor ion charge state from 2+ to 3+ also leads to more cleavage. The results of this study indicate that highly acidic peptides can be sequenced by ETD of complexes formed with Pr(III). Copyright © 2017 John Wiley & Sons, Ltd.

Gas-phase basicities of polyfunctional molecules. Part 4: Carbonyl groups as basic sites

This article constitutes the fourth part of a general review of the gas-phase protonation thermochemistry of polyfunctional molecules (Part 1: Theory and methods, Mass Spectrom Rev 2007, 26:775-835, Part 2: Saturated basic sites, Mass Spectrom Rev 2012, 31:353-390, Part 3: Amino acids, Mass Spectrom Rev 2012, 31:391-435). This fourth part is devoted to carbonyl containing polyfunctional molecules. After a short reminder of the methods of determination of gas-phase basicity and the underlying physicochemical concepts, specific examples are examined under two major chapters. In the first one, aliphatic and unsaturated (conjugated and cyclic) ketones, diketones, ketoalcohols, and ketoethers are considered. A second chapter describes the protonation energetic of gaseous acids and derivatives including diacids, diesters, diamides, anhydrides, imides, ureas, carbamates, amino acid derivatives, and peptides. Experimental data were re-evaluated according to the presently adopted basicity scale. Structural and energetic information given by G3 and G4 quantum chemistry computations on typical systems are presented.

Efficient and directed peptide bond formation in the gas phase via ion/ion reactions

Amide linkages are among the most important chemical bonds in living systems, constituting the connections between amino acids in peptides and proteins. We demonstrate the controlled formation of amide bonds between amino acids or peptides in the gas phase using ion/ion reactions in a mass spectrometer. Individual amino acids or peptides can be prepared as reagents by (i) incorporating gas phase-labile protecting groups to silence otherwise reactive functional groups, such as the N terminus; (ii) converting the carboxyl groups to the active ester of N-hydroxysuccinimide; and (iii) incorporating a charge site. Protonation renders basic sites (nucleophiles) unreactive toward the N-hydroxysuccinimide ester reagents, resulting in sites with the greatest gas phase basicities being, in large part, unreactive. The N-terminal amines of most naturally occurring amino acids have lower gas phase basicities than the side chains of the basic amino acids (i.e., those of histidine, lysine, or arginine). Therefore, reagents may be directed to the N terminus of an existing "anchor" peptide to form an amide bond by protonating the anchor peptide's basic residues, while leaving the N-terminal amine unprotonated and therefore reactive. Reaction efficiencies of greater than 30% have been observed. We propose this method as a step toward the controlled synthesis of peptides in the gas phase.

Gas phase basicities of polyfunctional molecules. Part 3: Amino acids

The present article is the third part of a general overview of the gas-phase protonation thermochemistry of polyfunctional molecules (first part: Mass Spectrom. Rev., 2007, 26:775-835, second part: Mass Spectrom. Rev., 2011, in press). This review is devoted to the 20 proteinogenic amino acids and is divided in two parts. In the first one, the experimental data obtained during the last 30 years using the equilibrium, thermokinetic and kinetic methods are presented. A general re-assignment of the values originating from these various experiments has been done on the basis of the commonly accepted Hunter & Lias 1998 gas-phase basicity scale in order to provide an homogeneous set of data. In the second part, theoretical investigations on gaseous neutral and protonated amino acids are reviewed. Conformational landscapes of both types of species were examined in order to provide theoretical protonation thermochemistry based on the truly identified most stable conformers. Proton affinities computed at the presently highest levels of theory (i.e. composite methods such as Gn procedures) are presented. Estimates of thermochemical parameters calculated using a Boltzmann distribution of conformers at 298K are also included. Finally, comparison between experiment and theory is discussed and a set of evaluated proton affinities, gas-phase basicities and protonation entropies is proposed.

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.

Collision-induced dissociative chemical cross-linking reagent for protein structure characterization: applied Edman chemistry in the gas phase

Chemical cross-linking combined with a subsequent enzymatic digestion and mass spectrometric analysis of the created cross-linked products presents an alternative approach to assess low-resolution protein structures and to gain insight into protein interfaces. In this contribution, we report the design of an innovative cross-linker based on Edman degradation chemistry, which leads to the formation of indicative mass shifted fragment ions and constant neutral losses (CNLs) in electrospray ionization (ESI)-tandem-mass spectrometry (MS/MS) product ion mass spectra, allowing an unambiguous identification of cross-linked peptides. Moreover, the characteristic neutral loss reactions facilitate automated analysis by multiple reaction monitoring suited for high throughput studies with good sensitivity and selectivity. The functioning of the novel cross-linker relies on the presence of a highly nucleophilic sulfur in a thiourea moiety, safeguarding for effective intramolecular attack leading to predictive and preferred cleavage of a glycyl-prolyl amide bond. Our innovative analytical concept and the versatile applicability of the collision-induced dissociative chemical cross-linking reagent are exemplified for substance P, luteinizing hormone releasing hormone LHRH and lysozyme. The novel cross-linker is expected to have a broad range of applications for probing protein tertiary structures and for investigating protein-protein interactions.

MS/MS of protonated polyproline peptides: the influence of N-terminal protonation on dissociation

A unique collision-induced dissociation pattern was observed for protonated polyproline peptides of length n in which y(n-2) and/or y(n-4) ions were formed in much higher abundance than any other product ions. Cleavage occurs only at every other amide bond, such that product ions are formed only from the losses of even numbers of proline residues. Exclusive losses of even numbers of proline residues were not observed from sodiated peptides. Further study of the tandem mass spectrometry (MS/MS) patterns of protonated proline-rich peptides showed that the substitution of alanine in the second position of polyproline peptides did not prevent the dominant formation of y(n-2) and y(n-4) ions. The loss of ProAla to form the y(8) ion from (ProAlaPro(8)NH(2)+H)(+) was as abundant as the loss of ProPro from (Pro(10)NH(2)+H)(+). However, modification of the peptides that presumably affected the location of the proton on the peptide did alter the MS/MS spectra. Pro(10) and Pro(5) with blocked N-termini or with arginine substituted for the first proline residue did not form abundant y(n-2) or y(n-4) ions. MS(3) and double resonance experiments showed that dissociation of intermediate y(n) product ions can produce y(n-2) ions, but are not necessary dissociation pathway intermediates. This analysis suggests that the ionizing proton must be located at the N-terminus for the peptide ion to dissociate in this manner.

Prediction of Low-Energy Collision-Induced Dissociation Spectra of Peptides

A kinetic model, based on the "mobile proton" model of peptide fragmentation, was developed to quantitatively simulate the low-energy collision-induced dissociation (CID) spectra of peptides dissociated in a quadrupole ion trap mass spectrometer. The model includes most fragmentation pathways described in the literature, plus some additional pathways based on the author's observations. The model was trained by optimizing parameters within the model for predictions of CID spectra of known peptides. A best set of parameters was optimized to obtain best match between the simulated spectra and the experimental spectra in a training data set. The performance of the mathematical model and the associated optimized parameter set used in the CID spectra simulation was evaluated by generating predictions for a large number of known peptides, which were not included in the training data set. It was shown that the model is able to predict peptide CID spectra with reasonable accuracy in fragment ion intensities for both singly and doubly charged peptide parent ions up to 2000 u in mass. The optimized parameter set was evaluated to gain insight into the collision-induced peptide fragmentation process.

Differentiation of a pair of diastereomeric tertiarybutoxycarbonylprolylproline ethyl esters by collision-induced dissociation of sodium adduct ions in electrospray ionization mass spectrometry and evidence for chiral recognition by ab initio molecular orbital calculations

The fragmentation of the sodium adduct ions for tert-butoxycarbonyl-L-prolyl-L-proline ethyl ester (Boc-L-Pro-L-Pro-OEt) was compared with that for Boc-D-Pro-L-Pro-OEt in positive-ion electrospray ionization (ESI) mass spectrometry. In the collision-induced dissociation (CID) mass spectra of the [M + Na](+) ions, the abundance of the [M + Na - C(CH(3))(3) + H](+) ion, which is due to the loss of a tert-butyl group from the [M + Na](+) ion for Boc-D-Pro-L-Pro-OEt, was about eight times higher than that for Boc-L-Pro-L-Pro-OEt. In addition, in the CID spectra of the sodium adduct fragment ion ([M + Na - Boc + H](+)), the abundance of the [M + Na - Boc - prolylresidue + H](+) ion, which is due to the loss of prolyl residue from the [M + Na - Boc + H](+) ion for Boc-L-Pro-L-Pro-OEt, was about five times higher than that for Boc-D-Pro-L-Pro-OEt. These results indicate that Boc-L-Pro-L-Pro-OEt was distinguished from Boc-D-Pro-L-Pro-OEt by the CID mass spectra of the sodium adduct ions in ESI mass spectrometry. The optimized geometries of the [M + Na](+) and the [M + Na - Boc + H](+) ions calculated by ab initio molecular orbital calculations suggest that the chiral recognition of these diastereomers was due to the difference of the orientation of a sodium ion to the oxygen and nitrogen atoms in dipeptide derivatives, and to the difference of the total energies between them.

Gas-phase basicity of glycine, alanine, proline, serine, lysine, histidine and some of their peptides by the thermokinetic method

The thermokinetic method is applied to a set of six amino acids (glycine, alanine, proline, serine, lysine, histidine) and 30 of their di- and tri-peptides for which experimental proton transfer rate constants were available. The comparison between the presently determined gas-phase basicities, GBs, of the amino acids with values obtained from equilibrium constant determination is generally good (a mean deviation of appoximately 3 kJ mol(-1) is observed). Derived proton affinities values are discussed. The gas-phase basicities of peptides provided by the present study correct several previously estimated values thus offering a more firm basis for structural discussion. Composite reaction efficiency curves indicate the existence, for several peptides, of at least two non-interconverting populations of protonated forms.

Cis-trans signatures of proline-containing tryptic peptides in the gas phase

High-resolution ion mobility/time-of-flight techniques were used to measure collision cross sections for 968 tryptic digest peptide ions obtained from digestion of common proteins. Here, we report a mobility signature that aids in identifying proline-containing peptides containing 4-10 residues. Of 129 peptides (< or = 10 residues in length) in the database that contain proline residues, 57% show multiple resolved features in the ion mobility distribution for at least one of the [M + H]+ or [M + 2H]2+ ions. These multiple features are attributed to different conformations that arise from populations of cis and trans forms of proline. The number of resolved peaks in the ion mobility distribution appears to be correlated with the peptide ion charge state and the number of proline residues in the peptide.

The proton affinity of proline analogs using the kinetic method with full entropy analysis

The proton affinity of proline analogs, L-azetidine-2-carboxylic acid (Aze), L-proline (Pro), and L-pipecolic acid (Pip), have been measured using the Armentrout modification of the extended kinetic method in a quadrupole ion trap instrument. Experimental values of 223.0 +/- 1.5, 224.9 +/- 1.6, and 225.6 +/- 1.6 kcal/mol have been determined for the 298K proton affinities of Aze, Pro, and Pip respectively. High level theoretical calculations using both MP2 and B3LYP methods at a variety of basis sets were carried out in order to give theoretical predictions for the 298 K proton affinity and gas phase basicity of all three analogs. Recommended values for the gas phase basicity and proton affinity for proline based on our work and other recent determinations are 216 +/- 2 and 224 +/- 2 kcal/mol.

Gas-phase basicities for ions from bradykinin and its des-arginine analogues

Apparent gas-phase basicities (GB(app)s) for [M + H]+ of bradykinin, des-Arg1-bradykinin and des-Arg9-bradykinin have been assigned by deprotonation reactions of [M + 2H]2+ in a Fourier transform ion cyclotron resonance mass spectrometer. With a GB(app) of 225.8 +/- 4.2 kcal x mol(-1), bradykinin [M + H]+ is the most basic of the ions studied. Ions from des-Arg1-bradykinin and des-Arg9-bradykinin have GB(app) values of 222.8 +/- 4.3 kcal x mol(-1) and 214.9 +/- 2.3 kcal x mol(-1), respectively. One purpose of this work was to determine a suitable reaction efficiency 'break point' for assigning GB(app) values to peptide ions using the bracketing method. An efficiency value of 0.1 (i.e. approximately 10% of all collisions resulting in a deprotonation reaction) was used to assign GB(app)s. Support for this criterion is provided by the fact that our GB(app) values for des-Arg1-bradykinin and des-Arg9-bradykinin are identical, within experimental error, to literature values obtained using a modified kinetic method. However, the GB(app)s for bradykinin ions from the two studies differ by 10.3 kcal x mol(-1). The reason for this is not clear, but may involve conformation differences produced by experimental conditions. The results may be influenced by salt-bridge conformers and/or by conformational changes caused by the use of a proton-bound heterodimer in the kinetic method. Factors affecting the basicities of these peptide ions are also discussed, and molecular modeling is used to provide information on protonation sites and conformations. The presence of two highly basic arginine residues on bradykinin results in its high GB(app), while the basicity of des-Arg1-bradykinin ions is increased by the presence of two proline residues at the N-terminus. The proline residue in the second position folds the peptide chain in a manner that increases intramolecular hydrogen bonding to the protonated N-terminal amino group of the proline at the first position.

Density functional computations of proton affinity and gas-phase basicity of proline

The proton affinity and gas-phase basicity of proline were evaluated by using density functional theory coupling the B3-LYP hybrid functional with the extended 6--311++G** basis set. Cis and trans conformations of the carboxyl moiety for both exo and endo ring structures were considered for the neutral proline. The results show that the most stable structure of proline has the endo ring conformation with the carboxyl group in the cis position. The structure at the global minimum is stabilized by an intramolecular hydrogen bond. The nitrogen of the ring in the exo form is the preferred protonation site. The calculated proton affinity (924.3 kJ mol(-1)) and gas-phase basicity (894.4 kJ mol(-1)) are in very good agreement with the experimental counterparts.

Estimation of proton affinity of proline and tryptophan under electrospray ionization conditions using the extended kinetic method

The relative order of the proton affinity (PA) of 20 naturally occurring amino acids has been determined under electrospray ionization conditions and compared with earlier studies of different research workers. The order we obtained is similar to that reported by other groups except in three cases viz., valine--aspartic acid, asparagine--glutamic acid and tryptophan--proline. The PA values of proline and tryptophan are determined by the extended kinetic method using amino acids themselves as reference bases. The PA values we thus obtained for proline and tryptophan are 219.9 and 221.6 kcal/mol, respectively.

Proton affinities of the N- and C-terminal segments arising upon the dissociation of the amide bond in protonated peptides

Dissociation of the amide bonds in a protonated peptide leads to N-terminal sequence fragments with cyclic structures and C-terminal sequence fragments with linear structures. The ionic fragments containing the N-terminus (bn) have been shown to be protonated oxazolones, whereas those containing the C-terminus (Yn) are protonated linear peptides. The coproduced neutral fragments are cyclic peptides from the N-terminus and linear peptides from the C-terminus. A likely determinant of these structural choices is the proton affinity (PA) of the described peptide segments. This study determines the PA values of such segments (Pep), i.e., cyclic and linear dipeptides and a relevant oxazolone, based on the dissociations of proton-bound dimers [Pep + Bi]H+ in which Bi is a reference base of known PA value (Cooks kinetic method). The dissociations are assessed at different internal energies to thereby obtain both proton affinities as well as entropies of protonation. For species with comparable amino acid composition, the proton affinity (and gas phase basicity) follows the order cyclic peptide << oxazolone approximately linear peptide. This ranking is consistent with dissociation of the protonated peptide via interconverting proton-bound complexes involving N-terminal oxazolone (O) or cyclopeptide (C) segments and C-terminal linear peptide segments (L), viz. O...H+...L reversible C...H+...L. N-terminal sequence ions (bn) are formed with oxazolone structures which can efficiently compete for the proton with the linear segments. On the other hand, N-terminal neutral fragments detach as cyclic peptides, with H+ now being retained by the more basic linear segment from the C-terminus to yield Yn.

Conformational interconversions in peptide beta-turns: analysis of turns in proteins and computational estimates of barriers

The two most important beta-turn features in peptides and proteins are the type I and type II turns, which differ mainly in the orientation of the central peptide unit. Facile conformational interconversion is possible, in principle, by a flip of the central peptide unit. Homologous crystal structures afford an opportunity to structurally characterize both possible conformational states, thus allowing identification of sites that are potentially stereochemically mobile. A representative data set of 250 high-resolution (</=2.0 A), non-homologous protein crystal structures and corresponding variant and homologous entries, obtained from the Brookhaven Protein Data Bank, was examined to identify turns that are assigned different conformational types (type I/type II) in related structures. A total of 55 examples of beta-turns were identified as possible candidates for a stereochemically mobile site. Of the 55 examples, 45 could be classified as a potential site for interconversion between type I and type II beta-turns, while ten correspond to flips from type I' to type II' structures. As a further check, the temperature factors of the central peptide unit carbonyl oxygen atom of the 55 examples were examined. The analysis reveals that the turn assignments are indeed reliable. Examination of the secondary structures at the flanking positions of the flippable beta-turns reveals that seven examples occur in the loop region of beta-hairpins, indicating that the formation of ordered secondary structures on either side of the beta-turn does not preclude local conformational variations. In these beta-turns, Pro (11 examples), Lys (nine examples) and Ser (seven examples) were most often found at the i+1 position. Glycine was found to occur overwhelmingly at position i+2 (28 examples), while Ser (seven examples) and Asn (six examples) were amongst the most frequent residues. Activation energy barriers for the interconversion between type I and type II beta-turns were computed using the peptide models Ac-Pro-Aib-NHMe and Ac-Pro-Gly-NHMe within the framework of the AM1 semi-empirical molecular orbital procedure. In order to have a uniform basis for comparison and to eliminate the distracting influence of the deviation of backbone dihedral angles from that expected for ideal beta-turns, the dihedral angles phii+1 and psii+2 were fixed at the ideal values (phii+1=-60 degrees and psii+2=0 degrees). The other two angles (psii+1 and phii+2) were varied systematically to go from type II to type I beta-turn structures. The computational results suggest that there exists one stereospecific, concerted flip of the central peptide unit involving correlated single bond rotation that can occur with an activation barrier of the order of 3 kcal/mol. The results presented here suggest that conformational variations in beta-turns are observed in protein crystal structures and such changes may be an important dynamic feature in solution.