On the occurrence of linear groups in proteins

Linear groups—polypeptide conformations based on a single repeating φ,ψ-pair—are a foundational concept in protein structure, yet how they are presented in textbooks is based largely on theoretical studies from the early days of protein structure analysis. Now, ultra-high resolution protein structures provide a resource for an accurate empirical and systematic assessment of the linear groups that truly exist in proteins. Here, a purely conformation-based survey of linear groups shows that only three distinct φ,ψ-regions occur: a diverse set of extended conformations mostly present as β-strands, a broad population of polyproline-II-like spirals, and a tight cluster that includes the highly populated α-helix and the conformationally-similar but much less populated 3₁₀-helix. Rare, short left-handed α-/3₁₀-helical turns with repeating φ,ψ-angles occur, but none are longer than three residues. Misperceptions dispelled by this study are the existence of 2.2₇- and π-helices as linear groups, the existence of specific ideal φ,ψ-angles for each linear group, and the existence of a substantive difference in the φ,ψ-preferences for parallel versus antiparallel β-strands. This study provides a concrete basis for updating and enhancing how we think about and teach the basics of protein structure.     

Kinetic steps for α-helix formation

The kinetics of α-helix formation in polyalanine and polyglycine eicosamers (20-mers) were examined using torsional-coordinate molecular dynamics (MD). Of one hundred fifty-five MD experiments on extended (Ala)₂₀ carried out for 0.5 ns each, 129 (83%) formed a persistent α-helix. In contrast, the extended state of (Gly)₂₀ only formed a right-handed α-helix in two of the 20 MD experiments (10%), and these helices were not as long or as persistent as those of polyalanine. These simulations show helix formation to be a competition between the rates of (a) forming local hydrogen bonds (i.e. hydrogen bonds between any residue i and its i + 2, i + 3, i + 4, or i + 5th neighbor) and (b) forming nonlocal hydrogen bonds (HBs) between residues widely separated in sequence. Local HBs grow rapidly into an α-helix; but nonlocal HBs usually retard helix formation by “trapping” the polymer in irregular, “balled-up” structures. Most trajectories formed some nonlocal HBs, sometimes as many as eight. But, for (Ala)₂₀, most of these eventually rearranged to form local HBs that lead to α-helices. A simple kinetic model describes the rate of converting nonlocal HBs into α-helices. Torsional-coordinate MD speeds folding by eliminating bond and angle degrees of freedom and reducing dynamical friction. Thus, the observed 210 ps half-life for helix formation is likely to be a lower bound on the real rate. However, we believe the sequential steps observed here mirror those of real systems. Proteins 33:343–357, 1998. © 1998 Wiley-Liss, Inc.     

Circular dichroism of the parallel β helical proteins pectate lyase C and E

The pectate lyases, PelC and PelE, have an unusual folding motif, known as a parallel β-helix, in which the polypeptide chain is coiled into a larger helix composed of three parallel β-sheets connected by loops having variable lengths and conformations. Since the regular secondary structure consists almost entirely of parallel β-sheets these proteins provide a unique opportunity to study the effect of parallel β-helical structure on circular dichroism (CD). We report here the CD spectra of PelC and PelE in the presence and absence of Ca²⁺, derive the parallel β-helical components of the spectra, and compare these results with previous CD studies of parallel β-sheet structure. The shape and intensity of the parallel β-sheet spectrum is distinctive and may be useful in identifying other proteins that contain the parallel β-helical folding motif. © 1995 Wiley-Liss, Inc.     

Design and synthesis of novel nonpolar host peptides for the determination of the 310- and α-helix compatibilities of α-amino acid buildig blocks: An assessment of α,α-disubstituted glycines

The present work describes three novel nonpolar host peptide sequences that provide a ready assessment of the 3₁₀- and α-helix compatibilities of natural and unnatural amino acids at different positions of small- to medium-size peptides. The unpolar peptides containing Ala, Aib, and a C-terminal p-iodoanilide group were designed in such a way that the peptides could be rapidly assembled in a modular fashion, were highly soluble in solvent mixtures of triflouroethanol and H₂O for CD- and two-dimensional (2D) nmr spectroscopic analyses, and showed excellent crystallinity suited for x-ray structure analysis. To validate our approach we synthesized 9-mer peptides 79a–96 (Table IV), 12-mer peptides 99–110c (Table V), and 10-mer peptides 120a–125d and 129–133 (Table VI and Scheme 8) incorporating a series of optically pure cyclic and open-chain (R)- and (S)-α,α-disubstituted glycines 1–10 (Figure 2). These amino acids are known to significantly modulate the conformations of small peptides. Based on x-ray structures of 9-mers 79a, 80, and 87 (Figures 4–7), 10-mers 124c, 131, and 132 (Figures 9–12), and 12-mer peptide 102b (Figure 13), CD spectra of all peptides recorded in acidic, neutral, and basic media and detailed 2D-nmr analyses of 9-mer peptide 86 and 12-mer 102b, several interesting conformational observations were made. Especially interesting results were obtained using the convex constraint CD analysis proposed by Fasman on 9-mer peptides 79a–d, 80, 81, 86, and 87, which allowed us to determine the relative content of 3₁₀- and α-helical conformations. These results were fully supported by the corresponding x-ray and 2D-nmr analyses. As a striking example we found that the (S)- and (R)-β-tetralin derived amino acids (R)- and (S)-1 show excellent α-helix stabilisation, more pronounced than Aib and Ala. These novel reference peptide sequences should help establish a scale for natural and unnatural amino acids concerning their intrinsic 3₁₀- and α-helix compatibilities at different positions of medium-sized peptides and thus improve our understanding in the folding processes of peptides. © 1997 John Wiley & Sons, Inc. Biopoly 42: 575–626, 1997     

Folding type-specific secondary structure propensities of amino acids,derived from α-helical, β-sheet, α/β, and α+β proteins of known structures

Folding type-specific secondary structure propensities of 20 naturally occurring amino acids have been derived from α-helical, β-sheet, α/β, and α+β proteins of known structures. These data show that each residue type of amino acids has intrinsic propensities in different regions of secondary structures for different folding types of proteins. Each of the folding types shows markedly different rank ordering, indicating folding type-specific effects on the secondary structure propensities of amino acids. Rigorous statistical tests have been made to validate the folding type-specific effects. It should be noted that α and β proteins have relatively small α-helices and β-strands forming propensities respectively compared with those of α+β and α/β proteins. This may suggest that, with more complex architectures than α and β proteins, α+β and α/β proteins require larger propensities to distinguish from interacting α-helices and β-strands. Our finding of folding type-specific secondary structure propensities suggests that sequence space accessible to each folding type may have differing features. Differing sequence space features might be constrained by topological requirement for each of the folding types. Almost all strong β-sheet forming residues are hydrophobic in character regardless of folding types, thus suggesting the hydrophobicities of side chains as a key determinant of β-sheet structures. In contrast, conformational entropy of side chains is a major determinant of the helical propensities of amino acids, although other interactions such as hydrophobicities and charged interactions cannot be neglected. These results will be helpful to protein design, class-based secondary structure prediction, and protein folding. © 1998 John Wiley & Sons, Inc. Biopoly 45: 35–49, 1998     

In silico and in vitro studies to elucidate the role of Cu2+ and galanthamine as the limiting step in the amyloid beta (1–42) fibrillation process

The formation of fibrils and oligomers of amyloid beta (Aβ) with 42 amino acid residues (Aβ_1–42) is the most important pathophysiological event associated with Alzheimer''s disease (AD). The formation of Aβ fibrils and oligomers requires a conformational change from an α-helix to a β-sheet conformation, which is encouraged by the formation of a salt bridge between Asp 23 or Glu 22 and Lys 28. Recently, Cu²⁺ and various drugs used for AD treatment, such as galanthamine (Reminyl®), have been reported to inhibit the formation of Aβ fibrils. However, the mechanism of this inhibition remains unclear. Therefore, the aim of this work was to explore how Cu²⁺ and galanthamine prevent the formation of Aβ_1–42 fibrils using molecular dynamics (MD) simulations (20 ns) and in vitro studies using fluorescence and circular dichroism (CD) spectroscopies. The MD simulations revealed that Aβ_1–42 acquires a characteristic U-shape before the α-helix to β-sheet conformational change. The formation of a salt bridge between Asp 23 and Lys 28 was also observed beginning at 5 ns. However, the MD simulations of Aβ₁₋₄₂ in the presence of Cu²⁺ or galanthamine demonstrated that both ligands prevent the formation of the salt bridge by either binding to Glu 22 and Asp 23 (Cu²⁺) or to Lys 28 (galanthamine), which prevents Aβ₁₋₄₂ from adopting the U-characteristic conformation that allows the amino acids to transition to a β-sheet conformation. The docking results revealed that the conformation obtained by the MD simulation of a monomer from the 1Z0Q structure can form similar interactions to those obtained from the 2BGE structure in the oligomers. The in vitro studies demonstrated that Aβ remains in an unfolded conformation when Cu²⁺ and galanthamine are used. Then, ligands that bind Asp 23 or Glu 22 and Lys 28 could therefore be used to prevent β turn formation and, consequently, the formation of Aβ fibrils.     

Deamidation of α‐synuclein

The rates of deamidation of α-synuclein and single Asn residues in 13 Asn-sequence mutants have been measured for 5 × 10⁻⁵M protein in both the absence and presence of 10⁻²M sodium dodecyl sulfate (SDS). In the course of these experiments, 370 quantitative protein deamidation measurements were performed and 37 deamidation rates were determined by ion cyclotron resonance Fourier transform mass spectrometry, using an improved whole protein isotopic envelope method and a mass defect method with both enzymatic and collision-induced fragmentation. The measured deamidation index of α-synuclein was found to be 0.23 for an overall deamidation half-time of 23 days, without or with SDS micelles, owing primarily to the deamidation of Asn(103) and Asn(122). Deamidation rates of 15 Asn residues in the wild-type and mutant proteins were found to be primary sequence controlled without SDS. However, the presence of SDS micelles slowed the deamidation rates of nine N-terminal region Asn residues, caused by the known three-dimensional structures induced through protein binding to SDS micelles.     

TNF‐α‐mediated signal transduction pathway is a major determinant of apoptosis in dilated cardiomyopathy

Although J2N-k strain of cardiomyopathic hamsters is an excellent model of dilated cardiomyopathy, the presence and mechanisms of apoptosis in the hearts of these genetically modified animals have not been investigated. This study examined the hypothesis that cardiac dysfunction and apoptosis in the cardiomyopathic hamsters were associated with tumour necrosis factor-alpha (TNF-α)-mediated signalling pathway involving the activation of some pro-apoptotic proteins and/or deactivation of some antiapoptotic proteins. Echocardiographic assessment of 31-week-old hamsters indicated an increase in the internal dimension of the left ventricle as well as decreases in the ejection fraction, fractional shortening and cardiac output without any evidence of cardiac hypertrophy. Increased level of TNF-α and apoptosis in cardiomyopathic hearts were accompanied by increased protein content for protein kinase C (PKC) -α and -ɛ isozymes as well as caspases 3 and 9. Phosphorylated protein content for p38 MAPK and NFκB was increased whereas that for Erk1/2, BAD and Bcl-2 was decreased in cardiomyopathic hearts. These results support the view that TNF-α and PKC isozymes may promote apoptosis due to the activation of p38 MAPK and deactivation of Erk1/2 pathways, and these changes may contribute toward the development of cardiac dysfunction in dilated cardiomyopathy.     

Relationship of sidechain hydrophobicity and α-helical propensity on the stability of the single-stranded amphipathic α-helix

The aim of the present investigation is to determine the effect of α-helical propensity and sidechain hydrophobicity on the stability of amphipathic α-helices. Accordingly, a series of 18-residue amphipathic α-helical peptides has been synthesized as a model system where all 20 amino acid residues were substituted on the hydrophobic face of the amphipathic α-helix. In these experiments, all three parameters (sidechain hydrophobicity, α-helical propensity and helix stability) were measured on the same set of peptide analogues. For these peptide analogues that differ by only one amino acid residue, there was a 0.96 kcal/mole difference in α-helical propensity between the most (Ala) and the least (Gly) α-helical analogue, a 12.1-minute difference between the most (Phe) and the least (Asp) retentive analogue on the reversed-phase column, and a 32.3°C difference in melting temperatures between the most (Leu) and the least (Asp) stable analogue. The results show that the hydrophobicity and α-helical propensity of an amino acid sidechain are not correlated with each other, but each contributes to the stability of the amphipathic α-helix. More importantly, the combined effects of α-helical propensity and sidechain hydrophobicity at a ratio of about 2:1 had optimal correlation with α-helix stability. These results suggest that both α-helical propensity and sidechain hydrophobicity should be taken into consideration in the design of α-helical proteins with the desired stability.     

Energy minimization studies on α-turns

Using a grid search technique, the entire conformational space of a system of four linked peptide units (tetrapeptide) was scanned to pick out geometrically possible 5→1 type hydrogen-bonded conformations defined as an α-turn. The energy minimization of these conformations led to 23 distinct minimum energy conformations (MECs) falling in 13 different classes. The presence of β and γ turn type hydrogen bonds along with 5→1 type hydrogen bond gave conformational variability in a given class. The occurrence of bifurcated hydrogen bonding network was a characteristic feature of most of the MECs. In many prototype MECs non-glycyl residues such as Ala and Pro could be accommodated. Comparison of MECs with the α-turn examples that are observed in proteins showed that the conformationally worked out MECs occurred in isolation in proteins, with the α-helical α-turn being distinctly the most predominant. © 1998 European Peptide Society and John Wiley & Sons, Ltd.     

Characterization of proline-containing α-helix (helix F model of bacteriorhodopsin) by molecular dynamics studies

Many of the bilayer spanning segments of membrane transport proteins contain proline residues, and most of them are believed to occur in α-helical form. A proline residue in the middle of an α-helix is known to produce a bend in the helix, and recent studies have focused on characterizing such a bend at atomic level. In the present case, molecular dynamics (MD) studies are carried out on helix F model of bacteriorhodopsin (BR) Ace-(Ala)₇-Trp-(Ala)₂-Tyr-Pro-(Ala)₂-Trp-(Ala)₈-NHMe and compared with Ace-(Ala)₇-Trp-(Ala)₂-Tyr-(Ala)₃-Trp-(Ala)₈-NHMe in which the proline is replaced by alanine. The bend in the helix is characterized by structural parameters such as kink angle (α), wobble angle (θ), virtual torsion angle (ρ), and the hydrogen bond distance d (O_p?3 … N_p+1). The average values and the flexibility involved in these parameters are evaluated. The correlation among the bend related parameters are estimated. The equilibrium side chain orientations of tryptophan and tyrosine residues are discussed and compared with those found in the recently proposed model of bacteriorhodopsin. Finally, a detailed characterization of the bend in terms of secondary structures such as α_I, α_II and goniometric helices are discussed, which can be useful in the interpretation of the experimental results on the secondary structures of membrane proteins involving the proline residue. © 1993 Wiley-Liss, Inc.     

The power of hard-sphere models: explaining side-chain dihedral angle distributions of Thr and Val

The energy functions used to predict protein structures typically include both molecular-mechanics and knowledge-based terms. In contrast, our approach is to develop robust physics- and geometry-based methods. Here, we investigate to what extent simple hard-sphere models can be used to predict side-chain conformations. The distributions of the side-chain dihedral angle χ₁ of Val and Thr in proteins of known structure show distinctive features: Val side chains predominantly adopt χ₁ = 180°, whereas Thr side chains typically adopt χ₁ = 60° and 300° (i.e., χ₁ = ±60° or g− and g⁺ configurations). Several hypotheses have been proposed to explain these differences, including interresidue steric clashes and hydrogen-bonding interactions. In contrast, we show that the observed side-chain dihedral angle distributions for both Val and Thr can be explained using only local steric interactions in a dipeptide mimetic. Our results emphasize the power of simple physical approaches and their importance for future advances in protein engineering and design.     

Age‐related sensitivity to endotoxin‐induced liver inflammation: Implication of inflammasome/IL‐1β for steatohepatitis

Aging is associated with increased vulnerability to inflammatory challenge. However, the effects of altered inflammatory response on the metabolic status of tissues or organs are not well documented. In this study, we present evidence demonstrating that lipopolysaccharide (LPS)-induced upregulation of the inflammasome/IL-1β pathway is accompanied with an increased inflammatory response and abnormal lipid accumulation in livers of aged rats. To monitor the effects of aging on LPS-induced inflammation, we administered LPS (2 mg kg⁻¹) to young (6-month old) and aged (24-month old) rats and found abnormal lipid metabolism in only aged rats with increased lipid accumulation in the liver. This lipid accumulation in the liver was due to the dysregulation of PPARα and SREBP1c. We also observed severe liver inflammation in aged rats as indicated by increased ALT levels in serum and increased Kupffer cells in the liver. Importantly, among many inflammation-associated factors, the aged rat liver showed chronically increased IL-1β production. Increased levels of IL-1β were caused by the upregulation of caspase-1 activity and inflammasome activation. In vitro studies with HepG2 cells demonstrated that treatment with IL-1β significantly induced lipid accumulation in hepatocytes through the regulation of PPARα and SREBP1c. In summary, we demonstrated that LPS-induced liver inflammation and lipid accumulation were associated with a chronically overactive inflammasome/IL-1β pathway in aged rat livers. Based on the present findings, we propose a mechanism of aging-associated progression of steatohepatitis induced by endotoxin, delineating a pathogenic role of the inflammasome/IL-1β pathway involved in lipid accumulation in the liver.     

The surface of β-sheet proteins contains amphiphilic regions which may provide clues about protein folding

A major bottleneck in the field of biochemistry is our limited understanding of the processes by which a protein folds into its native conformation. Much of the work on this issue has focused on the conserved core of the folded protein. However, one might imagine that a ubiquitous motif for unaided folding or for the recognition of chaperones may involve regions on the surface of the native structure. We explore this possibility by an analysis of the spatial distribution of regions with amphiphilic α-helical potential on the surface of β-sheet proteins. All proteins, Including β-sheet proteins, contain regions with amphiphilic α-helical potential. That is, any α-helix formed by that region would be amphiphilic, having both hydrophobic and hydrophilic surfaces. In the three-dimensional structure of all β-sheet proteins analyzed, we have found a distinct pattern in the spatial distribution of sequences with amphiphilic α-helical potential. The amphiphilic regions occur in ring shaped clusters approximately 20 to 30 Å in diameter on the surface of the protein. In addition, these regions have a strong preference for positively charged amino acids and a lower preference for residues not favorable to α-helix formation. Although the purpose of these amphiphilic regions which are not associated with naturally occurring α-helix is unknown, they may play a critical role in highly conserved processes such as protein folding. © 1996 Wiley-Liss, Inc.     

Refolding simulations of an isolated fragment of barnase into a native-like β hairpin: Evidence for compactness and hydrogen bonding as concurrent stabilizing factors

Experimental evidence and theoretical models both suggest that protein folding is initiated within specific fragments intermittently adopting conformations close to that found in the protein native structure. These folding initiation sites encompassing short portions of the protein are ideally suited for study in isolation by computational methods aimed at peering into the very early events of folding. We have used Molecular Dynamics (MD) technique to investigate the behavior of an isolated protein fragment formed by residues 85 to 102 of barnase that folds into a β hairpin in the protein native structure. Three independent MD simulations of 1.3 to 1.8 ns starting from unfolded conformations of the peptide portrayed with an all-atom model in water were carried out at gradually decreasing temperature. A detailed analysis of the conformational preferences adopted by this peptide in the course of the simulations is presented. Two of the unfolded peptide conformations fold into a hairpin characterized by native and a larger bulk of nonnative interactions. Both refolding simulations substantiate the close relationship between interstrand compactness and hydrogen bonding network involving backbone atoms. Persistent compactness witnessed by side-chain interactions always occurs concomitantly with the formation of backbone hydrogen bonds. No highly populated conformations generated in a third simulation starting from the remotest unfolded conformer relative to the native structure are observed. However, nonnative long-range and medium-range contacts with the aromatic moiety of Trp94 are spotted, which are in fair agreement with a former nuclear magnetic resonance study of a denaturing solution of an isolated barnase fragment encompassing the β hairpin. All this lends reason to believe that the 85–102 barnase fragment is a strong initiation site for folding. Proteins 29:212–227, 1997. © 1997 Wiley-Liss, Inc.     

Dynamic properties of pH-dependent structural organization of the amyloidogenic ��-protein (1�C40)

The structural organization of the amyloidogenic β-protein containing 40 amino acid residues (Aβ40) was studied by the high temperature molecular dynamics simulations in the acidic (pH ∼ 3) and basic (pH ∼ 8) pH regions. The obtained data suggest that the central Ala21-Gly29 segment of Aβ40 can adopt folded and partially unfolded structures. At the basic pH, this segment forms folded structures stabilized by electrostatic interactions and hydrogen bonds. At the acidic pH, it forms partially unfolded structures. Two other segments flanking to the central segment exhibit the propensity to adopt unstable interconverting α-helical, 3₁₀-helical and turn-like structures. One of these segments is comprised of the Ala30-Val36 residues at both of the considered pHs. The second segment is comprised of the Glu11-Phe20 at the basic pH and of the Glu11-Val24 residues at the acidic pHs. The revealed pH-dependent structuration of the Aβ40 allowed us to suggest a possible scenario for initial Aβ aggregation. According to this scenario, the occurrence of the partially unfolded states of the Ala21-Gly29 segment plays main role in the Aβ oligomerization process.Key words: amyloid-β protein, Alzheimer disease, oligomerization, fibril, electrostatic interactions, molecular dynamics simulations