A C-H triplebond O hydrogen bond stabilized polypeptide chain reversal motif at the C terminus of helices in proteins

The serendipitous observation of a C-H cdots, three dots, centered O hydrogen bond mediated polypeptide chain reversal in synthetic peptide helices has led to a search for the occurrence of a similar motif in protein structures. From a dataset of 634 proteins, 1304 helices terminating in a Schellman motif have been examined. The C-H triplebond O interaction between the T-4 C(alpha)H and T+1 Cz doublebond O group (C triplebond O< or =3.5A) becomes possible only when the T+1 residue adopts an extended beta conformation (T is defined as the helix terminating residue adopting an alpha(L) conformation). In all, 111 examples of this chain reversal motif have been identified and the compositional and conformational preferences at positions T-4, T, and T+1 determined. A marked preference for residues like Ser, Glu and Gln is observed at T-4 position with the motif being further stabilized by the formation of a side-chain-backbone O triplebond H-N hydrogen bond involving the side-chain of residue T-4 and the N-H group of residue T+3. In as many as 57 examples, the segment following the helix was extended with three to four successive residues in beta conformation. In a majority of these cases, the succeeding beta strand lies approximately antiparallel with the helix, suggesting that the backbone C-H triplebond O interactions may provide a means of registering helices and strands in an antiparallel orientation. Two examples were identified in which extended registry was detected with two sets of C-H cdots, three dots, centered O hydrogen bonds between (T-4) C(alpha)H triplebond O (T+1) and (T-8) C(alpha)H triplebondC doublebond O (T+3).     

Proline-induced constraints in alpha-helices

The disrupting effect of a prolyl residue on an α-helix has been analyzed by means of conformational energy computations. In the preferred, nearly α-helical conformations of Ac-Ala₄-Pro-NHMe and of Ac-Ala₇-Pro-Ala₇-NHMe, only the residue preceding Pro is not α-helical, while all other residues can occur in the α-helical A conformation; i.e., it is sufficient to introduce a conformational change of only one residue in order to accommodate proline in a distorted α-helix. Other low-energy conformations exist in which the conformational state of three residues preceding proline is altered considerably; on the other hand, another conformation in which these three residues retain the near-α-helical A-conformational state (with up to 26° changes of their dihedral angles ? and ψ, and a 48° change in one ω from those of the ideal α-helix) has a considerably higher energy. These conclusions are not altered by the substitution of other residues in the place of the Ala preceding Pro. The conformations of the peptide chain next to prolyl residues in or near an α-helix have been analyzed in 58 proteins of known structure, based on published atomic coordinates. Of 331 α-helices, 61 have a Pro at or next to their N-terminus, 21 have a Pro next to their C-terminus, and 30 contain a Pro inside the helix. Of the latter, 16 correspond to a break in the helix, 9 are located inside distorted first turns of the helix, and 5 are parts of irregular helices. Thus, the reported occurrence of prolyl residues next to or inside observed α-helices in proteins is consistent with the computed steric and energetic requirements of prolyl peptides.     

Variants of 3(10)-helices in proteins

An analysis of the shortest 3(10)-helices, containing three helical residues and two flanking capping residues that participate in two consecutive i + 3 --> i hydrogen bonds, shows that not all helices belong to the classic 3(10)-helix, where the three central residues adopt the right-handed helical conformation (alpha(R)). Three variants identified are: 3L10-helix with all residues in the left-handed helical region (alpha(L)), 3EL10-helix where the first residue is in the extended region followed by two residues in the alpha(L) conformation, and its mirror-image, the 3E'R10-helix. In the context of these helices, as well as the equivalent variants of alpha-helices, the length dependence of the handedness of secondary structures in protein structure is discussed. There are considerable differences in the amino acid preferences at different positions in the various types of 3(10)-helices. Each type of 3(10)-helix can be thought to be made up of an extension of a particular type of beta-turn (made up of residues i to i + 3) such that the (i + 3)th residue assumes the same conformation as the preceding residue. Distinct residue preferences at i and i + 3 positions seem to decide whether a particular stretch of four residues will be a beta-turn or a 3(10)-helix in the folded structure.     

Side-chain entropy effects on protein secondary structure formation

Loss of conformational entropy is one of the primary factors opposing protein folding. Both the backbone and side-chain of each residue in a protein will have their freedom of motion restricted in the final folded structure. The type of secondary structure of which a residue is part will have a significant impact on how much side-chain entropy is lost. Side-chain conformational entropies have previously been determined for folded proteins, simple models of unfolded proteins, alpha-helices, and a dipeptide model for beta-strands, but not for polyproline II (PII) helices. In this work, we present side-chain conformational estimates for the three regular secondary structure types: alpha-helices, beta-strands, and PII helices. Entropies are estimated from Monte Carlo computer simulations. Beta-strands are modeled as two structures, parallel and antiparallel beta-strands. Our data indicate that restraining a residue to the PII helix or antiparallel beta-strand conformations results in side-chain entropies equal to or higher than those obtained by restraining residues to the parallel beta-strand conformation. Side-chains in the alpha-helix conformation have the lowest side-chain entropies. The observation that extended structures retain the most side-chain entropy suggests that such structures would be entropically favored in unfolded proteins under folding conditions. Our data indicate that the PII helix conformation would be somewhat favored over beta-strand conformations, with antiparallel beta-strand favored over parallel. Notably, our data imply that, under some circumstances, residues may gain side-chain entropy upon folding. Implications of our findings for protein folding and unfolded states are discussed.     

Exploring the sequence patterns in the alpha-helices of proteins

This paper reports an extensive sequence analysis of the alpha-helices of proteins. alpha-Helices were extracted from the Protein Data Bank (PDB) and were divided into groups according to their sizes. It was found that some amino acids had differential propensity values for adopting helical conformation in short, medium and long alpha-helices. Pro and Trp had a significantly higher propensity for helical conformation in short helices than in medium and long helices. Trp was the strongest helix conformer in short helices. Sequence patterns favoring helical conformation were derived from a neighbor-dependent sequence analysis of proteins, which calculated the effect of neighboring amino acid type on the propensity of residues for adopting a particular secondary structure in proteins. This method produced an enhanced statistical significance scale that allowed us to explore the positional preference of amino acids for alpha-helical conformations. It was shown that the amino acid pair preference for alpha-helix had a unique pattern and this pattern was not always predictable by assuming proportional contributions from the individual propensity values of the amino acids. Our analysis also yielded a series of amino acid dyads that showed preference for alpha-helix conformation. The data presented in this study, along with our previous study on loop sequences of proteins, should prove useful for developing potential 'codes' for recognizing sequence patterns that are favorable for specific secondary structural elements in proteins.     

Helix folding simulations with various initial conformations.

Using a solvent-referenced energy calculation, a 16-residue peptide with alanine side chains folded into predominantly alpha-helical conformations during constant temperature (274 K) simulations. From different initial conformations, helical conformations were reached and the multiple energy minima did not become a serious problem. Under the same conditions, the simulation did not indiscriminately fold a sequence such as polyglycine into stable helices. Interesting observations from the simulations relate to the folding mechanism. The electrostatic interactions between the successive amides favored extended conformations (or beta strands) and caused energy barriers to helix folding. beta-bends were observed as intermediates during helix nucleation. The helix propagation toward the C-terminus seemed faster than that toward the N-terminus. In helical conformations, hydrogen bond oscillation between the i,i+ 4 and the i,i+3 patterns was observed. The i,i+3 hydrogen bonds occurred more frequently during helix propagation and deformation near both ends of the helical segment.     

Conformational effects of C(alpha,alpha)-dipropargylglycine as a constrained residue

A useful synthon to approach artificial phenylalanyl peptides in a [2 + 2 + 2] cycloaddition reaction, C(alpha,alpha)-dipropargylglycine (Dprg) is examined for its conformational preferences as a constrained residue. Crystal structure analysis and preliminary NMR results establish possible preference of the residue for folded (alpha) rather than extended (beta) region of the straight phi,psi conformational space. Boc-Dprg-L-Leu-OMe (1) displays two molecular conformations within the same crystallographic asymmetric unit, with Dprg in the alpha(R) or alpha(L) conformation, participating in a type I beta-turn or an alpha(L)-alpha(R)-type fold, in which Leu(2) assumes the alpha(R) conformation stereochemically favored for an L-chiral residue. Boc-Dprg-D-Val-L-Leu-OMe (2) displays a type I' beta-turn conformation in crystal, with both Dprg(1) and D-Val(2) assuming the alpha(L) conformation stereochemically favored for a D-chiral residue, with 4 --> 1 type hydrogen bond linking L-Leu(3) NH with Boc CO. NMR analysis using temperature variation, solvent titration, and a spin probe study suggests a fully solvent-exposed nature of Dprg NH, ruling out a fully extended C(5)-type conformation for this residue, and solvent sequestered nature of L-Leu(3) NH, suggesting possibility of a beta-turn due to Dprg assuming a folded conformation.     

Pituitary adenylate cyclase activating polypeptide (PACAP) with 27 residues. Conformation determined by 1H NMR and CD spectroscopies and distance geometry in 25% methanol solution.

The conformation of pituitary adenylate cyclase activating polypeptide with 27 residues (PACAP27) has been determined by two-dimensional NMR and CD spectroscopies and distance geometry in 25% methanol. Residues 9-20 and 22-25 have well-defined conformations but other residues do not show ordered conformations. The conformation of residues 9-20 is composed of three distinct regions of beta turn-like conformation (residues 9-12), alpha helix (residues 12-14) and the looser helical conformation (residues 15-20), while residues 22-24 form alpha helix. PACAP27 has a 2 helices separated by a disordered region similar to a VIP analog reported by Fry et al. but is distinct from the VIP analog in the position of the first helix, which is shifted by 2 residues toward the C-terminus, and in the form of the second helix [Fry, D.C., Madison, V.S., Bolin, D.R., Greeley, D.N., Toome, V. and Wegrzynski, B.B. (1989) Biochemistry 28, 2399-2409].     

Spatial structure of gramicidin A in organic solvents. 1H-NMR analysis of conformation heterogeneity in ethanol

Structural features of double helices formed by polypeptides with alternating L- and D-amino acid residues were analysed. It was found that the map of short distances (less than 4 A) between protons of the two backbones is unique for each double helix type and even its fragment implies unambiguously parameters of the helix (i.e. parallel or antiparallel, handedness, pitch of helix, relative shift of polypeptide chains). By analysis of two-dimensional 1H-NMR spectra (COSY, RELSY, HOHAHA, NOESY), proton resonances of [Val1]gramicidin A (GA) in the ethanol solution were assigned. The results obtained show that the solution contains five stable conformations of GA in comparable concentrations. Monomer of GA is in a random coil conformation. Specific maps of short interproton distances for the other four species (1-4) were obtained by means of two dimensional nuclear Overhauser effect spectroscopy. The maps as well as spin-spin couplings of the H-NC alpha-H protons and solvent accessibilities of the individual amide groups correspond to four types of double helices pi pi LD 5,6 with 5.6 residues per turn. The double helices are related to the Veatch species 1-4 of GA. Species 1 and 2 are left-handed parallel double helices increase increase pi pi LD 5,6 with different relative shift of polypeptide chains. Species 3 is a left-handed antiparallel double helix increase decrease pi pi LD 5,6 and species 4 is a right-handed parallel double helix increase increase LD 5,6. In the dimers helices are fixed by the maximum number (28) of interbackbone hydrogen bonds NH...O = C possible for these structures. Species 1, 3 and 4 have C2 symmetry axes. Relationship between gramicidin A spatial structures induced by various media is discussed.     

Conformational features of a hexapeptide model Ac-TGAAKA-NH2 corresponding to a hydrated alpha helical segment from glyceraldehyde 3-phosphate dehydrogenase: implications for the role of turns in helix folding

Recent analysis of alpha helices in protein crystal structures, available in literature, revealed hydrated alpha helical segments in which, water molecule breaks open helix 5-->1 hydrogen bond by inserting itself, hydrogen bonds to both C=O and NH groups of helix hydrogen bond without disrupting the helix hydrogen bond, and hydrogen bonds to either C=O or NH of helix hydrogen bond. These hydrated segments display a variety of turn conformations and are thought to be 'folding intermediates' trapped during folding-unfolding of alpha helices. A role for reverse turns is implicated in the folding of alpha helices. We considered a hexapeptide model Ac-1TGAAKA6-NH2 from glyceraldehyde 3-phosphate dehydrogenase, corresponding to a hydrated helical segment to assess its role in helix folding. The sequence is a site for two 'folding intermediates'. The conformational features of the model peptide have been investigated by 1H 2D NMR techniques and quantum mechanical perturbative configuration interaction over localized orbitals (PCILO) method. Theoretical modeling largely correlates with experimental observations. Based upon the amide proton temperature coefficients, the observed d alpha n(i, i + 1), d alpha n(i, i + 2), dnn(i, i + 1), d beta n(i, i + 1) NOEs and the results from theoretical modeling, we conclude that the residues of the peptide sample alpha helical and neck regions of the Ramachandran phi, psi map with reduced conformational entropy and there is a potential for turn conformations at N and C terminal ends of the peptide. The role of reduced conformational entropy and turn potential in helix formation have been discussed. We conclude that the peptide sequence can serve as a 'folding intermediate' in the helix folding of glyceraldehyde 3-phosphate dehydrogenase.     

Structures of the alpha L I domain and its complex with ICAM-1 reveal a shape-shifting pathway for integrin regulation

The structure of the I domain of integrin alpha L beta 2 bound to the Ig superfamily ligand ICAM-1 reveals the open ligand binding conformation and the first example of an integrin-IgSF interface. The I domain Mg2+ directly coordinates Glu-34 of ICAM-1, and a dramatic swing of I domain residue Glu-241 enables a critical salt bridge. Liganded and unliganded structures for both high- and intermediate-affinity mutant I domains reveal that ligand binding can induce conformational change in the alpha L I domain and that allosteric signals can convert the closed conformation to intermediate or open conformations without ligand binding. Pulling down on the C-terminal alpha 7 helix with introduced disulfide bonds ratchets the beta 6-alpha 7 loop into three different positions in the closed, intermediate, and open conformations, with a progressive increase in affinity.     

1H NMR studies of the solution conformations of an analogue of the C-peptide of ribonuclease A

Two-dimensional NMR experiments have been performed on a peptide, succinyl-AE-TAAAKFLRAHA-NH2, related to the amino-terminal sequence of ribonuclease A. This peptide contains 50-60% helix in 0.1 M NaCl solution, pH 5.2, 3 degrees C, as measured by circular dichroism. NOESY spectra of the peptide in aqueous solution at low temperatures show a number of NOE connectivities that are used to determine the highly populated conformations of the peptide in solution. Short-range dNN(i, i + 1) and d alpha N(i, i + 1) connectivities and medium-range d alpha beta(i, i + 3) and d alpha N(i, i + 3) connectivities are detected. The pattern of NOE connectivities unambiguously establishes the presence of helix in this peptide. The magnitudes of the 3JHN alpha coupling constants and the intensities of the dNN(i, i + 1) and d alpha N(i,i + 1) NOEs allow the evaluation of the position of the helix along the peptide backbone. These data indicate that the amino terminus of the peptide is less helical than the remainder of the peptide. The observation of several long-range NOEs that are atypical of helices indicates the presence of a high population of peptide molecules in which the first three residues are distorted out of the helical conformation. The absence of these NOEs in a related peptide, RN-31, in which Arg 10 has been changed to Ala, suggests that this distortion at the amino-terminal end of the peptide arises from the formation of a salt bridge between Glu 2 and Arg 10.(ABSTRACT TRUNCATED AT 250 WORDS)     

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.     

Analysis of the code relating sequence to conformation in globular proteins. An informational analysis of the role of the residue in determining the conformation of its neighbours in the primary sequence

1. The effect exerted by a residue on the conformation of neighbouring residues was analysed by using data from nine globular proteins of known sequence and conformation. 2. An information measure was used which estimated the role of a residue in influencing neighbouring conformations and also its tendency to influence the lengths of runs of residues in that conformation. This measure was estimated for each residue in all conformations defined by domains on the varphi, psi diagram. 3. Plots of the information measure yielded an intercept, which was a measure of intra-residue information for a residue. The slope was a measure of the statistical co-operativity or tendency of the residue to influence the occurrence of its neighbours in runs of a particular conformation. Both parameters are a function of the residue type. Statistical co-operativity is found in the alpha(1)-helical (H(1)) and beta-pleated-sheet (P(2)) conformations and, to a lesser extent, in their distorted variants H(2) and P(1). 4. The directional nature of these influences for H(1) and P(2) conformations is illustrated by plots of the information measure against the distance m from the residue, for m=-10 to +10. 5. The results for statistical co-operativity are discussed in relation to theories of helix-coil and pleated-sheet-coil transitions. The value of the information-theory-derived parameters in obtaining s parameters for the Zimm & Bragg (1959) equations is illustrated. 6. Directional effects are discussed with particular relation to mechanisms of the termination of helices and the involvement of the alpha(II) conformation and also to discontinuities in pleated-sheet conformations.     

Role of azaamino acid residue in beta-turn formation and stability in designed peptide.

The structural perturbation induced by C(alpha)-->N(alpha) exchange in azaamino acid-containing peptides was predicted by ab initio calculation of the 6-31G* and 3-21G* levels. The global energy-minimum conformations for model compounds, For-azaXaa-NH2 (Xaa=Gly, Ala, Leu) appeared to be the beta-turn motif with a dihedral angle of phi= +/- 90 degrees, psi=0 degrees. This suggests that incorporation of the azaXaa residue into the i+2 position of designed peptides could stabilize the beta-turn structure. The model azaLeu-containing peptide, Boc-Phe-azaLeu-Ala-OMe, which is predicted to adopt a beta-turn conformation was designed and synthesized in order to experimentally elucidate the role of the azaamino acid residue. Its structural preference in organic solvents was investigated using 1H NMR, molecular modelling and IR spectroscopy. The temperature coefficients of amide protons, the characteristic NOE patterns, the restrained molecular dynamics simulation and IR spectroscopy defined the dihedral angles [ (phi i+1, psi i+1) (phi i+2, psi i+2)] of the Phe-azaLeu fragment in the model peptide, Boc-Phe-azaLeu-Ala-OMe, as [(-59 degrees, 127 degrees) (107 degrees, -4 degrees)]. This solution conformation supports a betaII-turn structural preference in azaLeu-containing peptides as predicted by the quantum chemical calculation. Therefore, intercalation of the azaamino acid residue into the i+2 position in synthetic peptides is expected to provide a stable beta-turn formation, and this could be utilized in the design of new peptidomimetics adopting a beta-turn scaffold.     

Analysis of protein main-chain solvation as a function of secondary structure

We have analysed the hydration of main-chain carbonyl and amide groups in 24 high-resolution well-refined protein structures as a function of the secondary structure in which these polar groups occur. We find that main-chain atoms in beta-sheets are as hydrated as those in alpha-helices, with most interactions involving "free" amide and carbonyl groups that do not participate in secondary structure hydrogen bonds. The distributions of water molecules around these non-bonded carbonyl groups reflect specific steric interactions due to the local secondary structure. Approximately 20% and 4%, respectively of bonded carbonyl and amide groups interact with solvent. These include interactions with carbonyl groups on the exposed faces of alpha-helices that have been correlated previously with bending of the helix. Water molecules interacting with alpha-helices occur mainly at the amino and carbonyl termini of the helices, in which case the solvent sites maintain the hydrogen bonding by bridging between residues i and i-3 or i-4 at the amino terminus and between i and i+3 or i+4 at the carbonyl terminus. We also see a number of solvent-mediated Ncap and Ccap interactions. The water molecules interacting with beta-sheets occur mainly at the edges, in which case they extend the sheet structure, or at the ends of strands, in which case they extend the beta-ladder. In summary, the solvent networks appear to extend the hydrogen-bonding structure of the secondary structures. In beta-turns, which usually occur at the surface of a protein, exposed amide and carbonyl groups are often hydrated, especially close to glycine residues. Occasionally water molecules form a bridge between residues i and i+3 in the turn and this may provide extra stabilization.     

Helical screw sense of peptide molecules: The pentapeptide system (Aib)4/L-Val[L-(αMe)Val] in the crystal state

The crystal-state preferred conformations of six N^α-blocked pentapeptide esters, each containing four helicogenic, achiral α-aminoisobutyric acid (Aib) residues followed by one chiral L -valine (L -Val) or C^α-methyl-L -valine [(αMe)Val] residue at the C-terminus, have been assessed by x-ray diffraction analysis. In all of the compounds the  (Aib)₄ sequence is folded in a regular 3₁₀-helical conformation. In the four pentapeptides characterized by the L -(αMe)Val residue two conformationally distinct molecules occur in the asymmetric unit. Conversely, only one molecule is observed in the asymmetric unit of two pentapeptides with the C-terminal L -Val residue. In the L -Val based peptides the helical screw sense of the  (Aib)₄ sequence is right-handed, whereas in the L  (αMe)Val analogues both right- and left-handed helical screw senses concomitantly occur in the two crystallographically independent molecules. © 1998 John Wiley & Sons, Inc. Biopoly 46: 433–443, 1998     

Dissecting α-helices: Position-specific analysis of α-helices in globular proteins

An analysis of the amino acid distributions at 15 positions, viz., N“, N′, Ncap, N1, N2, N3, N4, Mid, C4, C3, C2, C1, Ccap, C′, and C” in 1,131 α-helices reveals that each position has its own unique characteristics. In general, natural helix sequences optimize by identifying the residues to be avoided at a given position and minimizing the occurrence of these avoided residues rather than by maximizing the preferred residues at various positions. Ncap is most selective in its choice of residues, with six amino acids (S, D, T, N, G, and P) being preferred at this position and another 11 (V, I, F, A, K, L, Y, R, E, M, and Q) being strongly avoided. Ser, Asp, and Thr are all more preferred at Ncap position than Asn, whose role at helix N-terminus has been highlighted by earlier analyses. Furthermore, Asn is also found to be almost equally preferred at helix C-terminus and a novel structural motif is identified, involving a hydrogen bond formed by N^δ2 of Asn at Ccap or C1 position, with the backbone carbonyl oxygen four residues inside the helix. His also forms a similar motif at the C-terminus. Pro is the most avoided residue in the main body (N4 to C4 positions) and at C-ter-minus, including Ccap of an α-helix. In 1,131 α-helices, no helix contains Pro at C3 or C2 positions. However, Pro is highly favoured at N1 and C′. The doublet X-Pro, with Pro at C′ position and extended backbone conformation for the X residue at Ccap, appears to be a common structural motif for termination of α-helices, in addition to the Schellman motif. Main body of the helix shows a high preference for aliphatic residues Ala, Leu, Val, and Ile, while these are avoided at helix termini. A propensity scale for amino acids to occur in the middle of helices has been obtained. Comparison of this scale with several previously reported scales shows that this scale correlates best with the experimentally determined values. Proteins 31:460–476, 1998. © 1998 Wiley-Liss, Inc.     

An electronic effect on protein structure

The well-known preference of the peptide bond for the trans conformation has been attributed to steric effects. Here, we show that a proline residue with an N-formyl group (H(i-1)-C'(i-1)=O(i-1)), in which H(i-1) presents less steric hindrance than does O(i-1), likewise prefers a trans conformation. Thus, the preference of the peptide bond for the trans conformation cannot be explained by steric effects alone. Rather, an n --> pi* interaction between the oxygen of the peptide bond (O(i-1)), and the subsequent carbonyl carbon in the polypeptide chain (C'(i)) also contributes to this preference. The O(i-1) and C'(i) distance and O(i-1).C'(i)=O(i) angle are especially favorable for such an n --> pi* interaction in a polyproline II helix. We propose that this electronic effect provides substantial stabilization to this and other elements of protein structure.