The structure of the ends of α-helices in globular proteins: effect of additional hydrogen bonds and implications for helix formation

We prepared a set of about 2000 α-helices from a relational database of high-resolution three-dimensional structures of globular proteins, and identified additional main chain i ← i+3 hydrogen bonds at the ends of the helices (i.e., where the hydrogen bonding potential is not fulfilled by canonical i ← i+4 hydrogen bonds). About one-third of α-helices have such additional hydrogen bonds at the N-terminus, and more than half do so at the C-terminus. Although many of these additional hydrogen bonds at the C-terminus are associated with Schellman loops, the majority are not. We compared the dihedral angles at the termini of α-helices having or lacking the additional hydrogen bonds. Significant differences were found, especially at the C-terminus, where the dihedral angles at positions C2 and C1 in the absence of additional hydrogen bonds deviate substantially from those occurring within the α-helix. Using a novel approach we show how the structure of the C-terminus of the α-helix can emerge from that of constituent overlapping α-turns and β-turns, which individually show a variation in dihedral angles at different positions. We have also considered the direction of propagation of the α-helix using this approach. If one assumes that helices start as a single α-turn and grow by successive addition of further α-turns, the paths for growth in the N → C and C → N directions differ in a way that suggests that extension in the C → N direction is favored.     

The geometry of alpha-sheet: Implications for its possible function as amyloid precursor in proteins

alpha-sheet has been proposed as the main constituent of the prefibrillar intermediate during amyloid formation. Here the helical parameters of the alpha-sheet strand are calculated from average main-chain dihedral angles reported from molecular dynamics simulations. It is an almost linear polypeptide that forms a right-handed helix of about 100 A diameter, with 100 residues and a rise of 30 A per turn. The strands are curved but untwisted, which implies that neighboring strands need not coil to make interstrand hydrogen bonds. This suggests that compared to beta-sheets in native folded proteins, alpha-sheets can be larger and stack more easily to create extensive 3D blocks. It is shown that alpha-sheet is related to a category of structures termed "mirror" structures. Mirror structures have repetitive pairs of main-chain dihedral angles at residues i and i+1 that satisfy the condition phi(i) (+1) = -psi(i), psi(i) (+1) = -phi(i). They are uniquely identified by the two orientations of their peptide planes, specified by phi(i) and psi(i). Their side chains point alternately in opposite directions. Interestingly, their conformations are insensitive to phi(i) and psi(i) in that the pseudo dihedral angle formed by four consecutive C(alpha) atoms is always close to 180 degrees . There are two types: "beta-mirror" and "alpha-mirror" structure; beta-mirror structures relate to beta-sheet by small peptide plane rotations, of less than 90 degrees , while alpha-mirror structures are close to alpha-sheet and relate to beta-sheet by approximately 180 degrees peptide plane flips. Most mirror structures, and in particular the alpha-mirror, form wide helices with diameters 50-70 A. Their gentle curvature, and therefore that of the alpha-sheet, arises from the orientation of successive peptide units causing the difference in the bond angles at the C and N atoms of the peptide unit to gradually change the direction of the chain.     

Revisiting the Ramachandran plot: hard-sphere repulsion, electrostatics, and H-bonding in the alpha-helix

What determines the shape of the allowed regions in the Ramachandran plot? Although Ramachandran explained these regions in terms of 1–4 hard‐sphere repulsions, there are discrepancies with the data where, in particular, the α_R, α_L, and β‐strand regions are diagonal. The α_R‐region also varies along the α‐helix where it is constrained at the center and the amino terminus but diffuse at the carboxyl terminus. By analyzing a high‐resolution database of protein structures, we find that certain 1–4 hard‐sphere repulsions in the standard steric map of Ramachandran do not affect the statistical distributions. By ignoring these steric clashes (N···H_i+1 and O_i?1···C), we identify a revised set of steric clashes (C^β···O, O_i?1···N_i+1, C^β···N_i+1, O_i?1···C^β, and O_i?1···O) that produce a better match with the data. We also find that the strictly forbidden region in the Ramachandran plot is excluded by multiple steric clashes, whereas the outlier region is excluded by only one significant steric clash. However, steric clashes alone do not account for the diagonal regions. Using electrostatics to analyze the conformational dependence of specific interatomic interactions, we find that the diagonal shape of the α_R and α_L‐regions also depends on the optimization of the N···H_i+1 and O_i?1···C interactions, and the diagonal β‐strand region is due to the alignment of the CO and NH dipoles. Finally, we reproduce the variation of the Ramachandran plot along the α‐helix in a simple model that uses only H‐bonding constraints. This allows us to rationalize the difference between the amino terminus and the carboxyl terminus of the α‐helix in terms of backbone entropy.     

Expanded turn conformations: characterization and sequence-structure correspondence in alpha-turns with implications in helix folding

Like the beta-turns, which are characterized by a limiting distance between residues two positions apart (i, i+3), a distance criterion (involving residues at positions i and i+4) is used here to identify alpha-turns from a database of known protein structures. At least 15 classes of alpha-turns have been enumerated based on the location in the phi,psi space of the three central residues (i+1 to i+3)-one of the major being the class AAA, where the residues occupy the conventional helical backbone torsion angles. However, moving towards the C-terminal end of the turn, there is a shift in the phi,psi angles towards more negative phi, such that the electrostatic repulsion between two consecutive carbonyl oxygen atoms is reduced. Except for the last position (i+4), there is not much similarity in residue composition at different positions of hydrogen and non-hydrogen bonded AAA turns. The presence or absence of Pro at i+1 position of alpha- and beta-turns has a bearing on whether the turn is hydrogen-bonded or without a hydrogen bond. In the tertiary structure, alpha-turns are more likely to be found in beta-hairpin loops. The residue composition at the beginning of the hydrogen bonded AAA alpha-turn has similarity with type I beta-turn and N-terminal positions of helices, but the last position matches with the C-terminal capping position of helices, suggesting that the existence of a "helix cap signal" at i+4 position prevents alpha-turns from growing into helices. Our results also provide new insights into alpha-helix nucleation and folding.     

A study of a Calpha,beta-didehydroalanine homo-oligopeptide series in the solid-state by 13C cross-polarization magic angle spinning NMR.

The fully extended peptide conformation (2.0(5)-helix) has been investigated for the first time in the solid-state by 13C cross-polarization magic angle spinning NMR. The compounds examined are members of a terminally protected, homo-oligopeptide series (from monomer through hexamer) based on Calpha,beta-didehydroalanine.     

Amplification,sequencing and characterization of pectin methyl esterase inhibitor 51 gene in Tectona grandis L.f.

Tectona grandis L.f. (Teak), a very important source of incomparable timber, withstands a wide range of tropical deciduous conditions. We achieved partial amplification of pectin methylesterase inhibitor 51 (PMEI) gene in teak by E. pilularis cinnamoyl Co-A reductase (CCR) gene specific primer. The amplified teak gene was of 750 bp, 79% identity and 97% query cover with PMEI of Sesamum indicum. The phylogenetic tree clustered the amplified gene with PMEI of database plant species, Erythranthe guttata and Sesamum indicum (87% bootstrap value). On conversion to amino acid sequence, the obtained protein comprised 237 amino acids. However, PMEI region spanned from 24 to 171 amino acids, 15.94 kDa molecular weight, 8.97 pI value and C₆₉₇H₁₁₁₇N₁₉₉O₂₁₁S₉ molecular formula with four conserved cysteine residues as disulfide bridges. 25.9 % protein residues were hydrophilic, 42.7% hydrophobic and 31.2% neutral. Teak 3D PMEI protein structure corresponded well with Arabidopsis thaliana and Actinidia deliciosa PMEIs. The gene maintains integrity of pectin component of middle lamella of primary cell wall and confers tolerance against various kinds of stresses. Teak conferred with overexpression of PMEI may secure a wide adaptability as well as luxuriant timber productivity and quality in adverse/ fluctuating/ scarce climatic and environmental conditions of tropical forests.     

Folding thermodynamics of three beta-sheet peptides: a model study

We study the folding thermodynamics of a beta-hairpin and two three-stranded beta-sheet peptides using a simplified sequence-based all-atom model, in which folding is driven mainly by backbone hydrogen bonding and effective hydrophobic attraction. The native populations obtained for these three sequences are in good agreement with experimental data. We also show that the apparent native population depends on which observable is studied; the hydrophobicity energy and the number of native hydrogen bonds give different results. The magnitude of this dependence matches well with the results obtained in two different experiments on the beta-hairpin.     

Influence of N-terminal residue stereochemistry on the prolyl amide geometry and the conformation of 5-tert-butylproline type VI beta-turn mimics.

The effects of N-terminal amino acid stereochemistry on prolyl amide geometry and peptide turn conformation were investigated by coupling both L- and D-amino acids to (2S, 5R)-5-tert-butylproline and L-proline to generate, respectively, N-(acetyl)dipeptide N'-methylamides 1 and 2. Prolyl amide cis- and trans-isomers were, respectively, favored for peptides 1 and 2 as observed by proton NMR spectroscopy in water, DMSO and chloroform. The influence of solvent composition on amide proton chemical shift indicated an intramolecular hydrogen bond between the N'-methylamide proton and the acetamide carbonyl for the major conformer of dipeptides (S)-1, that became less favorable in (R)-1 and 2. The coupling constant (3J(NH,alpha)) values for the cis-isomer of (R)-1 indicated a phi2 dihedral angle value characteristic of a type VIb beta-turn conformation in solution. X-ray crystallographic analysis of N-acetyl-D-leucyl-5-tert-butylproline N'-methylamide (R)-lb showed the prolyl residue in a type VIb beta-turn geometry possessing an amide cis-isomer and psi3-dihedral angle having a value of 157 degrees, which precluded an intramolecular hydrogen bond. Intermolecular hydrogen bonding between the leucyl residues of two turn structures within the unit cell positioned the N-terminal residue in a geometry where their phi2 and psi2 dihedral angle values were not characteristic of an ideal type VIb turn. The circular dichroism spectra of tert-butylprolyl peptides (S)- and (R)-1b were found not to be influenced by changes in solvent composition from water to acetonitrile. The type B spectrum exhibited by (S)-1b has been previously assigned to a type VIa beta-turn conformation [Halab L, Lubell WD. J. Org. Chem. 1999; 64: 3312-3321]. The type C spectrum exhibited by the (R)-lb has previously been associated with type II' beta-turn and alpha-helical conformations in solution and appears now to be also characteristic for a type VIb geometry.     

Transition state isosteres of the gamma-glutamyl peptide bond hydrolysis: synthesis and characterization of the psi (CH2NH) pseudopeptide analogue of glutathione.

The fully deprotected glutathione analogue containing the aminomethylene unit as transition state isostere of the gamma-Glu-Cys peptide bond was synthesized for the first time and characterized in both the reduced and oxidized forms.     

A conformational equilibrium in a protein fragment caused by two consecutive capping boxes: 1H-, 13C-NMR, and mutational analysis.

The conformational properties of an 18 residues peptide spanning the entire sequence, L1KTPA5QFDAD10ELRAA15MKG, of the first helix (A-helix) of domain 2 of annexin I, were thoroughly investigated. This fragment exhibits several singular features, and in particular, two successive potential capping boxes, T3xxQ6 and D8xxE11. The former corresponds to the native hydrogen bond network stabilizing the alpha helix N-terminus in the protein; the latter is a non-native capping box able to break the helix at residue D8, and is observed in the domain 2 partially folded state. Using 2D-NMR techniques, we showed that two main populations of conformers coexist in aqueous solution. The first corresponds to a single helix extending from T3 to K17. The second corresponds to a broken helix at residue Ds. Four mutants, T3A, F7A, D8A, and E11A, were designed to further analyze the role of key amino acids in the equilibrium between the two ensembles of conformers. The sensitivity of NMR parameters to account for the variations in the populations of conformers was evaluated for each peptide. Our data show the delta13Calpha chemical shift to be the most relevant parameter. We used it to estimate the population ratio in the various peptides between the two main ensembles of conformers, the full helix and the broken helix. For the WT, E11A, and F7A peptides, these ratios are respectively 35/65, 60/40, 60/40. Our results were compared to the data obtained from helix/coil transition algorithms.     

Peptoid residues and beta-turn formation.

A set of terminally protected tripeptoids containing a residue of either N-methylglycine or N-isobutylglycine in position i + 1/i + 2 were synthesized and tested for intramolecularly H-bonded beta-turn formation. By exploiting FT-IR absorption and 1H NMR techniques, their folding tendencies were compared with those of a variety of reference peptides. The amount of beta-turn induction and the relative extent of the various types of intramolecularly H-bonded beta-turn conformers were determined in chloroform solution.     

Beta-hairpin formation in aqueous solution and in the presence of trifluoroethanol: a (1)H and (13)C nuclear magnetic resonance conformational study of designed peptides

In order to check our current knowledge on the principles involved in beta-hairpin formation, we have modified the sequence of a 3:5 beta-hairpin forming peptide with two different purposes, first to increase the stability of the formed 3:5 beta-hairpin, and second to convert the 3:5 beta-hairpin into a 2:2 beta-hairpin. The conformational behavior of the designed peptides was investigated in aqueous solution and in 30% trifluoroethanol (TFE) by analysis of the following nuclear magnetic resonance (NMR) parameters: nuclear Overhauser effect (NOE) data, and C(alpha)H, (13)C(alpha), and (13)C(beta) conformational shifts. From the differences in the ability to adopt beta-hairpin structures in these peptides, we have arrived to the following conclusions: (i) beta-Hairpin population increases with the statistical propensity of residues to occupy each turn position. (ii) The loop length, and in turn, the beta-hairpin type, can be modified as a function of the type of turn favored by the loop sequence. These two conclusions reinforce previous results about the importance of beta-turn sequence in beta-hairpin folding. (iii) Side-chain packing on each face of the beta-sheet may play a major role in beta-hairpin stability; hence simplified analysis in terms of isolated pair interactions and intrinsic beta-sheet propensities is insufficient. (iv) Contributions to beta-hairpin stability of turn and strand sequences are not completely independent. (v) The burial of hydrophobic surface upon beta-hairpin formation that, in turn, depends on side-chain packing also contributes to beta-hairpin stability. (vi) As previously observed, TFE stabilizes beta-hairpin structures, but the extent of the contribution of different factors to beta-hairpin formation is sometimes different in aqueous solution and in 30% TFE.     

Conformational investigation of alpha,beta-dehydropeptides. XVI. Beta-turn tendency in Ac-Pro-DeltaXaa-NHMe: crystallographic and theoretical studies.

The crystal structures of two diastereomeric alpha,beta-dehydrobutyrine peptides Ac-Pro-(Z)-DeltaAbu-NHMe (I) and Ac-Pro-(E)-DeltaAbu-NHMe (II) have been determined. Both dehydropeptides adopt betaI-turn conformation characterized by the pairs of (phi(i+1), psi(i+1)) and (phi(i+2), psi(i+2)) angles as -66, -19, -97, 11 degrees for I and -59, -27, -119, 29 degrees for II. In each peptide, the betaI turn is stabilized by (i + 3) --> i intramolecular hydrogen bonds with N...O distance of 3.12 A for I and 2.93 A for II. These structures have been compared to the crystal structures of homologous peptides Ac-Pro-DeltaVal-NHMe and Ac-Pro-DeltaAla-NHMe. Theoretical analyses by DFT/B3LYP/6-31 + G** method of conformers formed by these four peptides and by the saturated peptide Ac-Pro-Ala-NHMe revealed that peptides with a (Z) substituent at the C(beta) (i+2) atom of dehydroamino acid, i.e. Ac-Pro-DeltaVal-NHMe and Ac-Pro-(Z)-DeltaAbu-NHMe, predominantly form beta turns, both in vacuo and in polar environment. The tendency to adopt beta-turn conformation is much weaker for the peptides lacking the (Z) substituent, Ac-Pro-(E)-DeltaAbu-NHMe and Ac-Pro-DeltaAla-NHMe. The latter adopts a semi-extended or an extended conformation in every polar environment, including a weakly polar solvent. The saturated peptide Ac-Pro-Ala-NHMe in vacuo prefers a beta-turn conformation, but in polar environment the differences between various conformers are small. The role of pi-electron correlation and intramolecular hydrogen bonds interaction in stabilizing the hairpin structures are discussed.     

CD studies on films of amyloid proteins and polypeptides: quantitative g-factor analysis indicates a common folding motif

Irrespective of the constituent protein, all amyloid fibrils show similar morphology in the electron microscope and x-ray diffraction patterns characteristic of a "cross-beta" structure, with extended beta-strands perpendicular to the fibril's long axis. Little is known about the amount or type of this structure. I have measured CD spectra of films formed from a number of amyloid proteins and polypeptides, and estimated their contents of extended secondary structure, by analysis of their g-factor spectra, the ratio of the CD and absorbance signals (P. McPhie, Analytical Biochemistry, 2001, Vol. 293, pp. 109-119). Amyloid films of Abeta-(1-40) peptide, beta-2-microglobulin, insulin, and three homopolypeptides show very intense CD spectra, compatible with the presence of a beta-helix-like structure, arranged in a common framework in the fibrils. The extent of this structure was estimated as 45-80% in the protein fibrils and 30-80% in the polypeptide fibrils.     

Conversion of type I 4:6 to 3:5 beta-turn types in human acidic fibroblast growth factor: effects upon structure, stability, folding, and mitogenic function

Human acidic fibroblast growth factor (FGF-1) is a member of the beta-trefoil superfold, a protein architecture that exhibits a characteristic threefold axis of structural symmetry. FGF-1 contains 11 beta-turns, the majority being type I 3:5; however, a type I 4:6 turn is also found at three symmetry-related locations. The relative uniqueness of the type I 4:6 turn in the FGF-1 structure suggests it may play a key role in the stability, folding, or function of the protein. To test this hypothesis a series of deletion mutations were constructed, the aim of which was to convert existing type I 4:6 turns at two locations into type I 3:5 turns. The results show it is possible to successfully substitute the type I 4:6 turn by a type I 3:5 turn with minimal impact upon protein stability or folding. Thus, these different turn structures, even though they differ in length, exhibit similar energetic properties. Additional sequence swapping mutations within the introduced type I 3:5 turns suggests that the turn sequence primarily affects stability but not turn structure (which appears dictated primarily by the local environment). Although the results suggest that a stable, foldable beta-trefoil protein may be designed utilizing a single turn type (type I 3:5), a type I 4:6 turn at turn 1 of FGF-1 appears essential for efficient mitogenic function.     

Folding of beta/alpha-unit scrambled forms of S. cerevisiae triosephosphate isomerase: Evidence for autonomy of substructure formation and plasticity of hydrophobic and hydrogen bonding interactions in core of (beta/alpha)8-barrel

We show that residues at the interfaces of protein-protein complexes have higher side-chain energy than other surface residues. Eight different sets of protein complexes were analyzed. For each protein pair, the complex structure was used to identify the interface residues in the unbound monomer structures. Side-chain energy was calculated for each surface residue in the unbound monomer using our previously developed scoring function.1 The mean energy was calculated for the interface residues and the other surface residues. In 15 of the 16 monomers, the mean energy of the interface residues was higher than that of other surface residues. By decomposing the scoring function, we found that the energy term of the buried surface area of non-hydrogen-bonded hydrophilic atoms is the most important factor contributing to the high energy of the interface regions. In spite of lacking hydrophilic residues, the interface regions were found to be rich in buried non-hydrogen-bonded hydrophilic atoms. Although the calculation results could be affected by the inaccuracy of the scoring function, patch analysis of side-chain energy on the surface of an isolated protein may be helpful in identifying the possible protein-protein interface. A patch was defined as 20 residues surrounding the central residue on the protein surface, and patch energy was calculated as the mean value of the side-chain energy of all residues in the patch. In 12 of the studied monomers, the patch with the highest energy overlaps with the observed interface. The results are more remarkable when only three residues with the highest energy in a patch are averaged to derive the patch energy. All three highest-energy residues of the top energy patch belong to interfacial residues in four of the eight small protomers. We also found that the residue with the highest energy score on the surface of a small protomer is very possibly the key interaction residue.     

Identification of a new subfamily of HNH nucleases and experimental characterization of a representative member, HphI restriction endonuclease

The restriction endonuclease (REase) R. HphI is a Type IIS enzyme that recognizes the asymmetric target DNA sequence 5'-GGTGA-3' and in the presence of Mg(2+) hydrolyzes phosphodiester bonds in both strands of the DNA at a distance of 8 nucleotides towards the 3' side of the target, producing a 1 nucleotide 3'-staggered cut in an unspecified sequence at this position. REases are typically ORFans that exhibit little similarity to each other and to any proteins in the database. However, bioinformatics analyses revealed that R.HphI is a member of a relatively big sequence family with a conserved C-terminal domain and a variable N-terminal domain. We predict that the C-terminal domains of proteins from this family correspond to the nuclease domain of the HNH superfamily rather than to the most common PD-(D/E)XK superfamily of nucleases. We constructed a three-dimensional model of the R.HphI catalytic domain and validated our predictions by site-directed mutagenesis and studies of DNA-binding and catalytic activities of the mutant proteins. We also analyzed the genomic neighborhood of R.HphI homologs and found that putative nucleases accompanied by a DNA methyltransferase (i.e. predicted REases) do not form a single group on a phylogenetic tree, but are dispersed among free-standing putative nucleases. This suggests that nucleases from the HNH superfamily were independently recruited to become REases in the context of RM systems multiple times in the evolution and that members of the HNH superfamily may be much more frequent among the so far unassigned REase sequences than previously thought.     

C(alpha)-hydroxymethyl methionine: synthesis, optical resolution and crystal structure of its (+)-N(alpha)-benzoyl derivative.

(R, S)-Methionine was transformed into C(alpha)-hydroxymethyl methionine by a route involving C(alpha)-hydroxymethylation of 2-phenyl-4-methylthioethyl-5-oxo-4,5-dihydro-1,3-oxazole. The absolute configuration of (-)-C(alpha)-hydroxymethyl methionine was elucidated to be (S) by chemical correlation with (S) (-)-C(alpha)-ethyl serine. Absolute structure determination (by single crystal X-ray diffraction) on N(alpha)-benzoyl-C(alpha)-hydroxymethyl methionine confirmed the (R)-configuration for the (+)-enantiomer. In addition, the X-ray diffraction analysis showed that the C(alpha,alpha)-disubstituted glycyl residue adopts the fully extended (C5) conformation.     

Mechanism-based protein design: attempted "nucleation-condensation" approach to a possible minimal helix-bundle protein

In an intended mechanism-based de novo approach, a 22-mer peptide was so designed as to make it both a stereochemically nucleatable and hydrophobically condensable minimal globular protein. Framework-like nucleation of a triple-helix bundle was targeted by employing as folding nucleators composite beta-turns that could both nucleate helices and place them in close juxtaposition for possible interhelical interaction. To promote the targeted triple-helix bundle to condense as a globular protein, an amphipathic sequence pattern was adopted for possible hydrophobic interhelical interaction. A predominantly helicogenic 22-mer amphipathic peptide was thus designed, punctuating it with composite type II'-III and type II-Asx type beta-turns as the helix nucleators cum chain reversal elements. The peptide made by solid-phase synthesis was shown by NMR and CD to be a nascent and distorted triple-helix bundle in a trifluoroethanol (TFE)-water mixture, but more or less a random coil in water. A fold nucleation effect is evident in the TFE-water mixture, but apparently the hydrophobic effect cannot sustain the peptide conformational order in water. A lack of synergy between folding nucleation and hydrophobic condensation of the peptide is possible. Indeed, a mismatch between the sequential H,P pattern of the peptide and its nascent-type globular fold in a TFE-water mixture is evident based on a simulated annealing study guided by NMR.