Effect of the N2 residue on the stability of the alpha-helix for all 20 amino acids

Acyl Capping Group Identity Effects on α-Helicity: On the Importance of Amide·Water Hydrogen Bonds to α-Helix Stability

Acyl capping groups stabilize α-helices relative to free N-termini by providing one additional C═Oi···Hi+4-N hydrogen bond. The electronic properties of acyl capping groups might also directly modulate α-helix stability: electron-rich N-terminal acyl groups could stabilize the α-helix by strengthening both i/i + 4 hydrogen bonds and i/i + 1 n → π* interactions. This hypothesis was tested in peptides X-AKAAAAKAAAAKAAGY-NH2, where X = different acyl groups. Surprisingly, the most electron-rich acyl groups (pivaloyl and iso-butyryl) strongly destabilized the α-helix. Moreover, the formyl group induced nearly identical α-helicity to that of the acetyl group, despite being a weaker electron donor for hydrogen bonds and for n → π* interactions. Other acyl groups exhibited intermediate α-helicity. These results indicate that the electronic properties of the acyl carbonyl do not directly determine the α-helicity in peptides in water. In order to understand these effects, DFT calculations were conducted on α-helical peptides. Using implicit solvation, α-helix stability correlated with acyl group electronics, with the pivaloyl group exhibiting closer hydrogen bonds and n → π* interactions, in contrast to the experimental results. However, DFT and MD calculations with explicit water solvation revealed that hydrogen bonding to water was impacted by the sterics of the acyl capping group. Formyl capping groups exhibited the closest water-amide hydrogen bonds, while pivaloyl groups exhibited the longest. In α-helices in the PDB, the highest frequency of close amide-water hydrogen bonds is observed when the N-cap residue is Gly. The combination of experimental and computational results indicates that solvation (hydrogen bonding of water) to the N-terminal amide groups is a central determinant of α-helix stability.

Synthesis and conformational preferences of peptides and proteins with cysteine sulfonic acid

Cysteine sulfonic acid (Cys-SO₃H; cysteic acid) is an oxidative post-translational modification of cysteine, resulting from further oxidation from cysteine sulfinic acid (Cys-SO₂H). Cysteine sulfonic acid is considered an irreversible post-translational modification, which serves as a biomarker of oxidative stress that has resulted in oxidative damage to proteins. Cysteine sulfonic acid is anionic, as a sulfonate (Cys-SO₃^-; cysteate), in the ionization state that is almost exclusively present at physiological pH (pK_a ∼ -2). In order to understand protein structural changes that can occur upon oxidation to cysteine sulfonic acid, we analyzed its conformational preferences, using experimental methods, bioinformatics, and DFT-based computational analysis. Cysteine sulfonic acid was incorporated into model peptides for α-helix and polyproline II helix (PPII). Within peptides, oxidation of cysteine to the sulfonic acid proceeds rapidly and efficiently at room temperature in solution with methyltrioxorhenium (MeReO₃) and H₂O₂. Peptides containing cysteine sulfonic acid were also generated on solid phase using trityl-protected cysteine and oxidation with MeReO₃ and H₂O₂. Using methoxybenzyl (Mob)-protected cysteine, solid-phase oxidation with MeReO₃ and H₂O₂ generated the Mob sulfone precursor to Cys-SO₂^- within fully synthesized peptides. These two solid-phase methods allow the synthesis of peptides containing either Cys-SO₃^- or Cys-SO₂^- in a practical manner, with no solution-phase synthesis required. Cys-SO₃^- had low PPII propensity for PPII propagation, despite promoting a relatively compact conformation in ϕ. In contrast, in a PPII initiation model system, Cys-SO₃^- promoted PPII relative to neutral Cys, with PPII initiation similar to Cys thiolate but less than Cys-SO₂^- or Ala. In an α-helix model system, Cys-SO₃^- promoted α-helix near the N-terminus, due to favorable helix dipole interactions and favorable α-helix capping via a sulfonate-amide side chain-main chain hydrogen bond. Across all peptides, the sulfonate side chain was significantly less ordered than that of the sulfinate. Analysis of Cys-SO₃^- in the PDB revealed a very strong propensity for local (i/i or i/i + 1) side chain-main chain sulfonate-amide hydrogen bonds for Cys-SO₃^-, with >80% of Cys-SO₃^- residues exhibiting these interactions. DFT calculations conducted to explore these conformational preferences indicated that side chain-main chain hydrogen bonds of the sulfonate with the intraresidue amide and/or with the i + 1 amide were favorable. However, hydrogen bonds to water or to amides, as well as interactions with oxophilic metals, were weaker for the sulfonate than the sulfinate, due to lower charge density on the oxygens in the sulfonate.

Uncovering Streptomyces-Derived Compounds as Cosmeceuticals for the Development of Improved Skin Photoprotection Products: An In Silico Approach to Explore Multi-Targeted Agents

The search for novel photoprotective substances has become a challenge in cosmeceutical research. Streptomyces-derived compounds can serve as a promising source of photoprotective agents to formulate skin photoprotection products, such as sunscreens. This study aimed to identify specialized metabolites with the potential to modulate UV-induced cellular damage in the skin by identifying potential multi-target-directed ligands. Using a combination of ligand- and target-based virtual screening approaches, a public compound library comprising 6524 Streptomyces-derived specialized metabolites was studied for their photoprotective capability. The compounds were initially filtered by safety features and then examined for their ability to interact with key targets in the photodamage pathway by molecular docking. A set of 50 commercially available UV filters was used as the benchmark. The protein–ligand stability of selected Streptomyces-derived compounds was also studied by molecular dynamics (MD) simulations. From the compound library, 1981 compounds were found to meet the safety criteria for topically applied products, such as low skin permeability and low or non-toxicity-alerting substructures. A total of 34 compounds had promising binding scores against crucial targets involved in UV-induced photodamage, such as serotonin-receptor subtype 5-HT2A, platelet-activating factor receptor, IL-1 receptor type 1, epidermal growth factor receptor, and cyclooxygenase-2. Among these compounds, aspergilazine A and phaeochromycin F showed the highest ranked interactions with four of the five targets and triggered complex stabilization over time. Additionally, the predicted UV-absorbing profiles also suggest a UV-filtering effect. Streptomyces is an encouraging biological source of compounds for developing topical products. After in silico protein–ligand interactions, binding mode and stabilization of aspergilazine A and phaeochromycin F led to the discovery of potential candidates as photodamage multi-target inhibitors. Therefore, they can be further explored for the formulation of skin photoprotection products.

Structural preferences of cysteine sulfinic acid: The sulfinate engages in multiple local interactions with the peptide backbone

Cysteine sulfinic acid (Cys-SO₂^-) is a non-enzymatic oxidative post-translational modification (PTM) that has been identified in hundreds of proteins. However, the effects of cysteine sulfination are in most cases poorly understood. Cys-SO₂^- is structurally distinctive, with long sulfur-carbon and sulfur-oxygen bonds, and with tetrahedral geometry around sulfur due to its lone pair. Cys-SO₂^- thus has a unique range of potential interactions with the protein backbone which could facilitate protein structural changes. Herein, the structural effects of cysteine oxidation to the sulfinic acid were investigated in model peptides and folded proteins using NMR spectroscopy, circular dichroism, bioinformatics, and computational studies. In the PDB, Cys-SO₂^- shows a greater preference for α-helix than Cys. In addition, Cys-SO₂^- is more commonly found in structures with φ > 0, including in multiple types of β-turn. Sulfinate oxygens engage in hydrogen bonds with adjacent (i or i + 1) amide hydrogens. Over half of sulfinates have at least one hydrogen bond with an adjacent amide, and several structures have hydrogen bonds with both adjacent amides. Alternately, sulfur or either oxygen can act as an electron donor for n→π* interactions with the backbone carbonyl of the same residue, as indicated by frequent S⋯CO or O⋯CO distances below the sums of their van der Waals radii in protein structures. In peptides, Cys-SO₂^- favored α-helical structure at the N-terminus, consistent with helix dipole effects and backbone hydrogen bonds with the sulfinate promoting α-helix. Cys-SO₂^- has only modestly greater polyproline II helix propensity than Cys-SH, likely due to competition from multiple side chain-backbone interactions. Cys-SO₂^- stabilizes the i+1 position of a β-turn relative to Cys-SH. Within proteins, the range of side chain-main chain interactions available to Cys-SO₂^- compared to Cys-SH provides a basis for potential changes in protein structure and function due to cysteine oxidation to the sulfinic acid.

Helix N-Cap Residues Drive the Acid Unfolding That Is Essential in the Action of the Toxin Colicin A

Numerous bacterial toxins and other virulence factors use low pH as a trigger to convert from water-soluble to membrane-inserted states. In the case of colicins, the pore-forming domain of colicin A (ColA-P) has been shown both to undergo a clear acidic unfolding transition and to require acidic lipids in the cytoplasmic membrane, whereas its close homologue colicin N shows neither behavior. Compared to that of ColN-P, the ColA-P primary structure reveals the replacement of several uncharged residues with aspartyl residues, which upon replacement with alanine induce an unfolded state at neutral pH. Here we investigate ColA-P's structural requirement for these critical aspartyl residues that are largely situated at the N-termini of α helices. As previously shown in model peptides, the charged carboxylate side chain can act as a stabilizing helix N-Cap group by interacting with free amide hydrogen bond donors. Because this could explain ColA-P destabilization when the aspartyl residues are protonated or replaced with alanyl residues, we test the hypothesis by inserting asparagine, glutamine, and glutamate residues at these sites. We combine urea (fluorescence and circular dichroism) and thermal (circular dichroism and differential scanning calorimetry) denaturation experiments with 1H-15N heteronuclear single-quantum coherence nuclear magnetic resonance spectroscopy of ColA-P at different pH values to provide a comprehensive description of the unfolding process and confirm the N-Cap hypothesis. Furthermore, we reveal that, in urea, the single domain ColA-P unfolds in two steps; low pH destabilizes the first step and stabilizes the second.

Computational Prediction and Design for Creating Iteratively Larger Heterospecific Coiled Coil Sets

A major biochemical goal is the ability to mimic nature in engineering highly specific protein-protein interactions (PPIs). We previously devised a computational interactome screen to identify eight peptides that form four heterospecific dimers despite 32 potential off-targets. To expand the speed and utility of our approach and the PPI toolkit, we have developed new software to derive much larger heterospecific sets (≥24 peptides) while directing against antiparallel off-targets. It works by predicting T_m values for every dimer on the basis of core, electrostatic, and helical propensity components. These guide interaction specificity, allowing heterospecific coiled coil (CC) sets to be incrementally assembled. Prediction accuracy is experimentally validated using circular dichroism and size exclusion chromatography. Thermal denaturation data from a 22-CC training set were used to improve software prediction accuracy and verified using a 136-CC test set consisting of eight predicted heterospecific dimers and 128 off-targets. The resulting software, qCIPA, individually now weighs core a-a' (II/NN/NI) and electrostatic g-e'⁺¹ (EE/EK/KK) components. The expanded data set has resulted in emerging sequence context rules for otherwise energetically equivalent CCs; for example, introducing intrahelical electrostatic charge blocks generated increased stability for designed CCs while concomitantly decreasing the stability of off-target CCs. Coupled with increased prediction accuracy and speed, the approach can be applied to a wide range of downstream chemical and synthetic biology applications, in addition more generally to impose specificity on structurally unrelated PPIs.

The severity of hereditary porphyria is modulated by the porphyrin exporter and Lan antigen ABCB6

Hereditary porphyrias are caused by mutations in genes that encode haem biosynthetic enzymes with resultant buildup of cytotoxic metabolic porphyrin intermediates. A long-standing open question is why the same causal porphyria mutations exhibit widely variable penetrance and expressivity in different individuals. Here we show that severely affected porphyria patients harbour variant alleles in the ABCB6 gene, also known as Lan, which encodes an ATP-binding cassette (ABC) transporter. Plasma membrane ABCB6 exports a variety of disease-related porphyrins. Functional studies show that most of these ABCB6 variants are expressed poorly and/or have impaired function. Accordingly, homozygous disruption of the Abcb6 gene in mice exacerbates porphyria phenotypes in the Fech^m1Pas mouse model, as evidenced by increased porphyrin accumulation, and marked liver injury. Collectively, these studies support ABCB6 role as a genetic modifier of porphyria and suggest that porphyrin-inducing drugs may produce excessive toxicities in individuals with the rare Lan(−) blood type.

Accumulation of intermediates of haem biosynthesis, porphyrins, is harmful and usually inherited, but it is unclear how the same mutation may make some individuals more ill than others. Here, the authors show that a porphyrin transporter ABCB6 is a modulator of porphyria, and that patients with functionally defective ABCB6 show more severe symptoms.

The importance of cooperative interactions and a solid-state paradigm to proteins: what Peptide chemists can learn from molecular crystals

Proteins and peptides in solution or in vivo share properties with both liquids and solids. More often than not, they are studied using the liquid paradigm rather than that of a solid. Studies of molecular crystals illustrate how the use of a solid paradigm may change the way that we consider these important molecules. Cooperative interactions, particularly those involving H-bonding, play much more important roles in the solid than in the liquid paradigms, as molecular crystals clearly illustrate. Using the solid rather than the liquid paradigm for proteins and peptides includes these cooperative interactions while application of the liquid paradigm tends to ignore or minimize them. Use of the solid paradigm has important implications for basic principles that are often implied about peptide and protein chemistry, such as the importance of entropy in protein folding and the nature of the hydrophobic effect. Understanding the folded states of peptides and proteins (especially alpha-helices) often requires the solid paradigm, whereas understanding unfolded states does not. Both theoretical and experimental studies of the energetics of protein and peptide folding require comparison to a suitable standard. Our perspective on these energetics depends on the reasonable choice of reference. The use of multiple reference states, particularly that of component amino acids in the gas phase, is proposed.

Deriving Heterospecific Self-Assembling Protein-Protein Interactions Using a Computational Interactome Screen

Interactions between naturally occurring proteins are highly specific, with protein-network imbalances associated with numerous diseases. For designed protein–protein interactions (PPIs), required specificity can be notoriously difficult to engineer. To accelerate this process, we have derived peptides that form heterospecific PPIs when combined. This is achieved using software that generates large virtual libraries of peptide sequences and searches within the resulting interactome for preferentially interacting peptides. To demonstrate feasibility, we have (i) generated 1536 peptide sequences based on the parallel dimeric coiled-coil motif and varied residues known to be important for stability and specificity, (ii) screened the 1,180,416 member interactome for predicted T_m values and (iii) used predicted T_m cutoff points to isolate eight peptides that form four heterospecific PPIs when combined. This required that all 32 hypothetical off-target interactions within the eight-peptide interactome be disfavoured and that the four desired interactions pair correctly. Lastly, we have verified the approach by characterising all 36 pairs within the interactome. In analysing the output, we hypothesised that several sequences are capable of adopting antiparallel orientations. We subsequently improved the software by removing sequences where doing so led to fully complementary electrostatic pairings. Our approach can be used to derive increasingly large and therefore complex sets of heterospecific PPIs with a wide range of potential downstream applications from disease modulation to the design of biomaterials and peptides in synthetic biology.

•

Naturally occurring protein–protein interactions (PPIs) are highly specific.

•

For designed PPIs, however, specificity can be notoriously difficult to engineer.

•

We have computationally screened a vast interactome to derive four heterospecific PPIs.

•

Eight peptides form four heterospecific coiled coils; all 32 off targets are disfavoured.

•

The method can derive larger and increasingly complex sets of heterospecific PPIs

Mapping side chain interactions at protein helix termini

Background

Interactions that involve one or more amino acid side chains near the ends of protein helices stabilize helix termini and shape the geometry of the adjacent loops, making a substantial contribution to overall protein structure. Previous work has identified key helix-terminal motifs, such as Asx/ST N-caps, the capping box, and hydrophobic and electrostatic interactions, but important questions remain, including: 1) What loop backbone geometries are favoured by each motif? 2) To what extent are multi-amino acid motifs likely to represent genuine cooperative interactions? 3) Can new motifs be identified in a large, recent dataset using the latest bioinformatics tools?

Results

Three analytical tools are applied here to answer these questions. First, helix-terminal structures are partitioned by loop backbone geometry using a new 3D clustering algorithm. Next, Cascade Detection, a motif detection algorithm recently published by the author, is applied to each cluster to determine which sequence motifs are overrepresented in each geometry. Finally, the results for each motif are presented in a CapMap, a 3D conformational heatmap that displays the distribution of the motif’s overrepresentation across loop geometries, enabling the rapid isolation and characterization of the associated side chain interaction. This work identifies a library of geometry-specific side chain interactions that provides a new, detailed picture of loop structure near the helix terminus. Highlights include determinations of the favoured loop geometries for the Asx/ST N-cap motifs, capping boxes, “big” boxes, and other hydrophobic, electrostatic, H-bond, and pi stacking interactions, many of which have not been described before.

Conclusions

This work demonstrates that the combination of structural clustering and motif detection in the sequence space can efficiently identify side chain motifs and map them to the loop geometries which they support. Protein designers should find this study useful, because it identifies side chain interactions which are good candidates for inclusion in synthetic helix-terminal loops with specific desired geometries, since they are used in nature to support these geometries. The techniques described here can also be applied to map side chain interactions associated with other structural components of proteins such as beta and gamma turns.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-015-0671-4) contains supplementary material, which is available to authorized users.

Local and macroscopic electrostatic interactions in single α-helices

The noncovalent forces that stabilize protein structures are not fully understood. One way to address this is to study equilibria between unfolded states and α-helices in peptides. Electrostatic forces-which include interactions between side chains, the backbone and side chains, and side chains and the helix macrodipole-are believed to contribute to these equilibria. Here we probe these interactions experimentally using designed peptides. We find that both terminal backbone-side chain and certain side chain-side chain interactions (which include both local effects between proximal charges and interatomic contacts) contribute much more to helix stability than side chain-helix macrodipole electrostatics, which are believed to operate at larger distances. This has implications for current descriptions of helix stability, the understanding of protein folding and the refinement of force fields for biomolecular modeling and simulations. In addition, this study sheds light on the stability of rod-like structures formed by single α-helices, which are common in natural proteins such as non-muscle myosins.

Effect of external pulling forces on the length distribution of peptides

BACKGROUND

The distribution of the length of a polypeptide, or that of the distance between any two of its atoms, is an important property as it can be analytically or numerically estimated for a number of polymer models. Importantly, it is directly measurable through a number of different experimental techniques. Length distributions can be straightforwardly assessed from molecular dynamics simulation; however, true convergence through full accurate coverage of the length range is difficult to achieve.

METHODS

The application of external constant force combined with the weighted-histogram analysis method (WHAM) is used to enhance sampling of unlikely 'long' or 'short' conformations and obtain the potential of mean force, while also collecting dynamic properties of the chain under variable tension.

RESULTS

We demonstrate the utility of constant force to enhance the sampling efficiency and obtain experimentally measurable quantities on a series of short peptides, including charge-rich sequences that are known to be highly helical but whose properties are distinct from those of helical peptides undergoing helix-coil transitions.

CONCLUSIONS

Force-enhanced sampling enhances the range and accuracy of the length-based potential of mean force of the peptide, in particular those sequences that contain increased numbers of charged residues.

GENERAL SIGNIFICANCE

This approach allows users to simultaneously probe the force-dependent behaviour of peptides directly, enhance the range and accuracy of the length-based PMF of the peptide and also test the convergence of simulations by comparing the overlap of PMF profiles from different constant forces. This article is part of a special issue entitled Recent developments of molecular dynamics.

Tyrosine phosphorylation within the intrinsically disordered cytosolic domains of the B-cell receptor: an NMR-based structural analysis

Intrinsically disordered proteins are found extensively in cell signaling pathways where they often are targets of posttranslational modifications e.g. phosphorylation. Such modifications can sometimes induce or disrupt secondary structure elements present in the modified protein. CD79a and CD79b are membrane-spanning, signal-transducing components of the B-cell receptor. The cytosolic domains of these proteins are intrinsically disordered and each has an immunoreceptor tyrosine-based activation motif (ITAM). When an antigen binds to the receptor, conserved tyrosines located in the ITAMs are phosphorylated which initiate further downstream signaling. Here we use NMR spectroscopy to examine the secondary structure propensity of the cytosolic domains of CD79a and CD79b in vitro before and after phosphorylation. The phosphorylation patterns are identified through analysis of changes of backbone chemical shifts found for the affected tyrosines and neighboring residues. The number of the phosphorylated sites is confirmed by mass spectrometry. The secondary structure propensities are calculated using the method of intrinsic referencing, where the reference random coil chemical shifts are measured for the same protein under denaturing conditions. Our analysis revealed that CD79a and CD79b both have an overall propensity for α-helical structure that is greatest in the C-terminal region of the ITAM. Phosphorylation of CD79a caused a decrease in helical propensity in the C-terminal ITAM region. For CD79b, the opposite was observed and phosphorylation resulted in an increase of helical propensity in the C-terminal part.

Effects of side chains in helix nucleation differ from helix propagation

Helix-coil transition theory connects observable properties of the α-helix to an ensemble of microstates and provides a foundation for analyzing secondary structure formation in proteins. Classical models account for cooperative helix formation in terms of an energetically demanding nucleation event (described by the σ constant) followed by a more facile propagation reaction, with corresponding s constants that are sequence dependent. Extensive studies of folding and unfolding in model peptides have led to the determination of the propagation constants for amino acids. However, the role of individual side chains in helix nucleation has not been separately accessible, so the σ constant is treated as independent of sequence. We describe here a synthetic model that allows the assessment of the role of individual amino acids in helix nucleation. Studies with this model lead to the surprising conclusion that widely accepted scales of helical propensity are not predictive of helix nucleation. Residues known to be helix stabilizers or breakers in propagation have only a tenuous relationship to residues that favor or disfavor helix nucleation.

OGlcNAcylation and phosphorylation have similar structural effects in α-helices: post-translational modifications as inducible start and stop signals in α-helices, with greater structural effects on threonine modification

OGlcNAcylation and phosphorylation are the major competing intracellular post-translational modifications of serine and threonine residues. The structural effects of both post-translational modifications on serine and threonine were examined within Baldwin model α-helical peptides (Ac-AKAAAAKAAAAKAAGY-NH₂ or Ac-YGAKAAAAKAAAAKAA-NH₂). At the N-terminus of an α-helix, both phosphorylation and OGlcNAcylation stabilized the α-helix relative to the free hydroxyls, with a larger induced structure for phosphorylation than for OGlcNAcylation, for the dianionic phosphate than for the monoanionic phosphate, and for modifications on threonine than for modifications on serine. Both phosphoserine and phosphothreonine resulted in peptides more α-helical than alanine at the N-terminus, with dianionic phosphothreonine the most α-helix-stabilizing residue here. In contrast, in the interior of the α-helix, both post-translational modifications were destabilizing with respect to the α-helix, with the greatest destabilization seen for threonine OGlcNAcylation at residue 5 and threonine phosphorylation at residue 10, with peptides containing either post-translational modification existing as random coils. At the C-terminus, both OGlcNAcylation and phosphorylation were destabilizing with respect to the α-helix, though the induced structural changes were less than in the interior of the α-helix. In general, the structural effects of modifications on threonine were greater than the effects on serine, because of both the lower α-helical propensity of Thr and the more defined induced structures upon modification of threonine than serine, suggesting threonine residues are particularly important loci for structural effects of post-translational modifications. The effects of serine and threonine post-translational modifications are analogous to the effects of proline on α-helices, with the effects of phosphothreonine being greater than those of proline throughout the α-helix. These results provide a basis for understanding the context-dependent structural effects of these competing protein post-translational modifications.

OGlcNAcylation and phosphorylation have opposing structural effects in tau: phosphothreonine induces particular conformational order

Phosphorylation and OGlcNAcylation are dynamic intracellular protein post-translational modifications that frequently are alternatively observed on the same serine and threonine residues. Phosphorylation and OGlcNAcylation commonly occur in natively disordered regions of proteins, and often have opposing functional effects. In the microtubule-associated protein tau, hyperphosphorylation is associated with protein misfolding and aggregation as the neurofibrillary tangles of Alzheimer's disease, whereas OGlcNAcylation stabilizes the soluble form of tau. A series of peptides derived from the proline-rich domain (residues 174-251) of tau was synthesized, with free Ser/Thr hydroxyls, phosphorylated Ser/Thr (pSer/pThr), OGlcNAcylated Ser/Thr, and diethylphosphorylated Ser/Thr. Phosphorylation and OGlcNAcylation were found by CD and NMR to have opposing structural effects on polyproline helix (PPII) formation, with phosphorylation favoring PPII, OGlcNAcylation opposing PPII, and the free hydroxyls intermediate in structure, and with phosphorylation structural effects greater than OGlcNAcylation. For tau196-209, phosphorylation and OGlcNAcylation had similar structural effects, opposing a nascent α-helix. Phosphomimic Glu exhibited PPII-favoring structural effects. Structural changes due to Thr phosphorylation were greater than those of Ser phosphorylation or Glu, with particular conformational restriction as the dianion, with mean (3)JαN = 3.5 Hz (pThr) versus 5.4 Hz (pSer), compared to 7.2, 6.8, and 6.2 Hz for Thr, Ser, and Glu, respectively, values that correlate with the backbone torsion angle ϕ. Dianionic phosphothreonine induced strong phosphothreonine amide protection and downfield amide chemical shifts (δmean = 9.63 ppm), consistent with formation of a stable phosphate-amide hydrogen bond. These data suggest potentially greater structural importance of threonine phosphorylation than serine phosphorylation due to larger induced structural effects.

The alphabet of intrinsic disorder: I. Act like a Pro: On the abundance and roles of proline residues in intrinsically disordered proteins

A significant fraction of every proteome is occupied by biologically active proteins that do not form unique three-dimensional structures. These intrinsically disordered proteins (IDPs) and IDP regions (IDPRs) have essential biological functions and are characterized by extensive structural plasticity. Such structural and functional behavior is encoded in the amino acid sequences of IDPs/IDPRs, which are enriched in disorder-promoting residues and depleted in order-promoting residues. In fact, amino acid residues can be arranged according to their disorder-promoting tendency to form an alphabet of intrinsic disorder that defines the structural complexity and diversity of IDPs/IDPRs. This review is the first in a series of publications dedicated to the roles that different amino acid residues play in defining the phenomenon of protein intrinsic disorder. We start with proline because data suggests that of the 20 common amino acid residues, this one is the most disorder-promoting.

Amino acid pair- and triplet-wise groupings in the interior of α-helical segments in proteins

A statistical approach has been applied to analyse primary structure patterns at inner positions of α-helices in proteins. A systematic survey was carried out in a recent sample of non-redundant proteins selected from the Protein Data Bank, which were used to analyse α-helix structures for amino acid pairing patterns. Only residues more than three positions apart from both termini of the α-helix were considered as inner. Amino acid pairings i, i+k (k=1, 2, 3, 4, 5), were analysed and the corresponding 20×20 matrices of relative global propensities were constructed. An analysis of (i, i+4, i+8) and (i, i+3, i+4) triplet patterns was also performed. These analysis yielded information on a series of amino acid patterns (pairings and triplets) showing either high or low preference for α-helical motifs and suggested a novel approach to protein alphabet reduction. In addition, it has been shown that the individual amino acid propensities are not enough to define the statistical distribution of these patterns. Global pair propensities also depend on the type of pattern, its composition and orientation in the protein sequence. The data presented should prove useful to obtain and refine useful predictive rules which can further the development and fine-tuning of protein structure prediction algorithms and tools.

Position-specific propensities of amino acids in the β-strand

Background

Despite the importance of β-strands as main building blocks in proteins, the propensity of amino acid in β-strands is not well-understood as it has been more difficult to determine experimentally compared to α-helices. Recent studies have shown that most of the amino acids have significantly high or low propensity towards both ends of β-strands. However, a comprehensive analysis of the sequence dependent amino acid propensities at positions between the ends of the β-strand has not been investigated.

Results

The propensities of the amino acids calculated from a large non-redundant database of proteins are found to be highly position-specific and vary continuously throughout the length of the β-strand. They follow an unexpected characteristic periodic pattern in inner positions with respect to the cap residues in both termini of β-strands; this periodic nature is markedly different from that of the α-helices with respect to the strength and pattern in periodicity. This periodicity is not only different for different amino acids but it also varies considerably for the amino acids belonging to the same physico-chemical group. Average hydrophobicity is also found to be periodic with respect to the positions from both termini of β-strands.

Conclusions

The results contradict the earlier perception of isotropic nature of amino acid propensities in the middle region of β-strands. These position-specific propensities should be of immense help in understanding the factors responsible for β-strand design and efficient prediction of β-strand structure in unknown proteins.

End-to-end and end-to-middle interhelical interactions: new classes of interacting helix pairs in protein structures

Helix-helix interactions are important for the structure, stability and function of alpha-helical proteins. Helices that either cross in the middle or show extensive contacts between each other, such as coiled coils, have been investigated in previous studies. Interactions between two helices can also occur only at the terminal regions or between the terminal region of one helix and the middle region of another helix. Examples of such helix pairs are found in aquaporin, H(+)/Cl(-) transporter and Bcl-2 proteins. The frequency of the occurrence of such ;end-to-end' (EE) and ;end-to-middle' (EM) helix pairs in protein structures is not known. Questions regarding the residue preferences in the interface and the mode of interhelical interactions in such helix pairs also remain unanswered. In this study, high-resolution structures of all-alpha proteins from the PDB have been systematically analyzed and the helix pairs that interact only in EE or EM fashion have been extracted. EE and EM helix pairs have been categorized into five classes (N-N, N-C, C-C, N-MID and C-MID) depending on the region of interaction. Nearly 13% of 5725 helix pairs belonged to one of the five classes. Analysis of single-residue propensities indicated that hydrophobic and polar residues prefer to occur in the C-terminal and N-terminal regions, respectively. Hydrophobic C-terminal interacting residues and polar N-terminal interacting residues are also highly conserved. A strong correlation exists between some of the residue properties (surface area/volume and length of side chains) and their preferences for occurring in the interface of EE and EM helix pairs. In contrast to interacting non-EE/EM helix pairs, helices in EE and EM pairs are farther apart. In these helix pairs, residues with large surface area/volume and longer side chains are preferred in the interfacial region.

Statistical mechanical model for helix-sheet-coil transitions in homopolypeptides

In this paper, we propose a simple statistical mechanical model to study the conformation transition between the alpha helix, beta sheet, and random coil in homopolypeptides. In our model, five parameters are introduced to obtain the partition function. There are two factors for helical propagation and initiation, which are the same as those used in the Zimm-Bragg model, and three newly introduced parameters for beta structures: the strand propagation factor for residues in beta strands and two correction factors for the initiation effect of the beta strand and beta sheet. Our model shows that the variation of these parameters may induce conformation transition from alpha helix or random coil to beta sheet. The sharpness of the transition depends on the initiation factors.

Lysine and arginine residues do not increase the helicity of alanine-rich peptide helices

The helix-disfavoring, versus alanine, propagation values of lysine (0.8) and arginine (1.0) residues placed centrally in an (Ala)(9) unit have been measured by (13)C NMR.

Amino acid pairing at the N- and C-termini of helical segments in proteins

A systematic survey was carried out in an unbiased sample of 815 protein chains with a maximum of 20% homology selected from the Protein Data Bank, whose structures were solved at a resolution higher than 1.6 A and with a R-factor lower than 25%. A set of 5556 subsequences with alpha-helix or 3(10)-helix motifs was extracted from the protein chains considered. Global and local propensities were then calculated for all possible amino acid pairs of the type (i, i + 1), (i, i + 2), (i, i + 3), and (i, i + 4), starting at the relevant helical positions N1, N2, N3, C3, C2, C1, and N-int (interior positions), and also at the first nonhelical positions in both termini of the helices, namely, N-cap and C-cap. The statistical analysis of the propensity values has shown that pairing is significantly dependent on the type of the amino acids and on the position of the pair. A few sequences of three and four amino acids were selected and their high prevalence in helices is outlined in this work. The Glu-Lys-Tyr-Pro sequence shows a peculiar distribution in proteins, which may suggest a relevant structural role in alpha-helices when Pro is located at the C-cap position. A bioinformatics tool was developed, which updates automatically and periodically the results and makes them available in a web site.

Structure-based interpretation of the mutagenesis database for the nucleotide binding domains of P-glycoprotein

P-glycoprotein (P-gp) is the most intensively studied eukaryotic ATP binding cassette (ABC) transporter, due to its involvement in the multidrug resistance phenotype of a number of cancers. In common with most ABC transporters, P-gp is comprised of two transmembrane domains (TMDs) and two nucleotide binding domains (NBD), the latter coupling ATP hydrolysis with substrate transport (efflux in the case of P-gp). Biochemical investigations over the past twenty years have attempted to unlock mechanistic aspects of P-glycoprotein through scanning and site-directed mutagenesis of both the TMDs and the NBDs. Contemporaneously, crystallographers have elucidated the atomic structure of numerous ABC transporter NBDs, as well as the intact structure (i.e. NBDs and TMDs) of a distantly related ABC-exporter Sav1866. Significantly, the structure of P-gp remains unknown, and only low resolution electron microscopy data exists. Within the current manuscript we employ crystallographic data for homologous proteins, and a molecular model for P-gp, to perform a structural interpretation of the existing "mutagenesis database" for P-gp NBDs. Consequently, this will enable testable predictions to be made that will result in further in-roads into our understanding of this clinically important drug pump.

Structural characterization of alpha-helices of implicitly solvated poly-alanine

The structural characteristics of alpha-helices in poly-alanine-based peptides have been investigated via molecular dynamics simulation with the goal of understanding the basic features of peptide simulations within the context of a model system, classical molecular dynamics with generalized Born (GB) solvation, and to shed insight into the formation and stabilization of alpha-helices in short peptides. The effects of peptide length, terminal charges, proline substitution, and temperature on the alpha-helical secondary structure have been studied. The simulations have shown that distinct secondary structure begins to develop in peptides with lengths approaching 10 residues while ambiguous structures occur in shorter peptides. The helical content of peptides with lengths > or =10 amino acids is observed to be nearly constant up to (Ala)(40). Interestingly, terminal charges and proline in the second position from the N-terminus alter the secondary structure locally with little effect on the overall alpha-helical content of the peptide. The free energy profile of helix formation was also investigated. A large increase in free energy accompanying the formation of helices with more than two consecutive hydrogen bonds in the (i, i + 4) pattern was observed while the free energy increases linearly with additional hydrogen bonds. Values for the change in enthalpy and entropy of helix nucleation and propagation are reported. Additionally the results obtained from the GB model are compared to explicit solvent simulations of two synthetic alanine-based peptides.

UV Raman spatially resolved melting dynamics of isotopically labeled polyalanyl peptide: slow alpha-helix melting follows 3(10)-helices and pi-bulges premelting

We used UV resonance Raman (UVRR) to examine the spatial dependence of the T-jump secondary structure relaxation of an isotopically labeled 21-residue mainly Ala peptide, AdP. The AdP penultimate Ala residues were perdeuterated, leaving the central residues hydrogenated, to allow separate monitoring of melting of the middle versus the end peptide bonds. For 5 to 30 degrees C T-jumps, the central peptide bonds show a approximately 2-fold slower relaxation time (189 +/- 31 ns) than do the exterior peptide bonds (97 +/- 15 ns). In contrast, for a 20 to 40 degrees C T-jump, the central peptide bond relaxation appears to be faster (56 +/- 6 ns) than that of the penultimate peptide bonds (131 +/- 46 ns). We show that, if the data are modeled as a two-state transition, we find that only exterior peptide bonds show anti-Arrhenius folding behavior; the middle peptide bonds show both normal Arrhenius-like folding and unfolding. This anti-Arrhenius behavior results from the involvement of pi-bulges/helices and 3(10)-helix states in the melting. The unusual temperature dependence of the (un)folding rates of the interior and exterior peptide bonds is due to the different relative (un)folding rates of 3(10)-helices, alpha-helices, and pi-bulges/helices. Pure alpha-helix unfolding rates are approximately 12-fold slower (approximately 1 micros) than that of pi-bulges and 3(10)-helices. In addition, we also find that the alpha-helix is most stable at the AdP N-terminus where eight consecutive Ala occur, whereas the three hydrophilic Arg located in the middle and at the C-terminus destabilize the alpha-helix in these regions and induce defects such as pi-bulges and 3(10)-helices.

A dynamic N-capping motif in cytochrome b5: evidence for a pH-controlled conformational switch

Apocytochrome b5 is a marginally stable protein exhibiting under native conditions a slow conformational exchange in its C-terminal region. The affected elements of secondary structure include a 3(10)-helix containing at its N-terminus a histidine Ncap and a subsequent proline. Participation of the neutral histidine side-chain in backbone amide capping lowers the imidazole pKa. To explore the nature of the conformational exchange in the protein and determine whether it is related to cis-trans isomerization of the His-Pro bond, three octapeptides encompassing the helix were synthesized and studied by NMR spectroscopy. One corresponded to the wild-type sequence, the second was the D-histidine epimer, and the third contained an alanine in place of the proline. It was found that the rates of cis-trans interconversion in the proline-containing peptides were slower than the rates of the conformational exchange in the protein. In addition, the wild-type peptide hinted at a predisposition for Ncap formation when in the trans configuration. Analysis of the pH response of the peptides and protein suggested that at pH near neutral, the conformational exchange detected in the protein involved only species with a trans His-Pro bond and could be approximated with a three-state model by which the terminal helix sampled a locally unfolded state. This state, which contained an uncapped histidine with a normal pKa, partitioned into neutral and protonated populations according to pH. The intrinsic conformational bias of the wild-type peptide and the pH-driven equilibria illustrated how a 3(10)-element could serve as a nucleation site for structural rearrangement.

The CXXC motif at the N terminus of an alpha-helical peptide

An active site containing a CXXC motif is always found in the thiol-disulphide oxidoreductase superfamily. A survey of crystal structures revealed that the CXXC motif had a very high local propensity (26.3 +/- 6.2) for the N termini of alpha-helices. A helical peptide with the sequence CAAC at the N terminus was synthesized to examine the helix-stabilizing capacity of the CXXC motif. Circular dichroism was used to confirm the helical nature of the peptide and study behavior under titration with various species. With DTT, a redox potential of E(o) = -230 mV was measured, indicating that the isolated peptide is reducing in nature and similar to native human thioredoxin. The pK(a) values of the individual Cys residues could not be separated in the titration of the reduced state, giving a single transition with an apparent pK(a) of 6.74 (+/-0.06). In the oxidized state, the N-terminal pK(a) is 5.96 (+/-0.05). Analysis of results with the modified helix-coil theory indicated that the disulfide bond stabilized the alpha-helical structure by 0.5 kcal/mol. Reducing the disulfide destabilizes the helix by 0.9 kcal/mol.

Protein secondary structure prediction for a single-sequence using hidden semi-Markov models

Background

The accuracy of protein secondary structure prediction has been improving steadily towards the 88% estimated theoretical limit. There are two types of prediction algorithms: Single-sequence prediction algorithms imply that information about other (homologous) proteins is not available, while algorithms of the second type imply that information about homologous proteins is available, and use it intensively. The single-sequence algorithms could make an important contribution to studies of proteins with no detected homologs, however the accuracy of protein secondary structure prediction from a single-sequence is not as high as when the additional evolutionary information is present.

Results

In this paper, we further refine and extend the hidden semi-Markov model (HSMM) initially considered in the BSPSS algorithm. We introduce an improved residue dependency model by considering the patterns of statistically significant amino acid correlation at structural segment borders. We also derive models that specialize on different sections of the dependency structure and incorporate them into HSMM. In addition, we implement an iterative training method to refine estimates of HSMM parameters. The three-state-per-residue accuracy and other accuracy measures of the new method, IPSSP, are shown to be comparable or better than ones for BSPSS as well as for PSIPRED, tested under the single-sequence condition.

Conclusions

We have shown that new dependency models and training methods bring further improvements to single-sequence protein secondary structure prediction. The results are obtained under cross-validation conditions using a dataset with no pair of sequences having significant sequence similarity. As new sequences are added to the database it is possible to augment the dependency structure and obtain even higher accuracy. Current and future advances should contribute to the improvement of function prediction for orphan proteins inscrutable to current similarity search methods.

Functional role of a protein foldon-An Ω-loop foldon controls the alkaline transition in ferricytochrome c

Hydrogen exchange results for cytochrome c and several other proteins show that they are composed of a number of foldon units which continually unfold and refold and account for some functional properties. Previous work showed that one Omega-loop foldon controls the rate of the structural switching and ligand exchange behavior of cytochrome c known as the alkaline transition. The present work tests the role of foldons in the alkaline transition equilibrium. We measured the effects of denaturant and 14 destabilizing mutations. The results show that the ligand exchange equilibrium is controlled by the stability of the same foldon unit implicated before. In addition, the results obtained confirm the epsilon-amino group of Lys79 and Lys73 as the alkaline replacement ligands and bear on the search for a triggering group.

Disruption of an intermonomer salt bridge in the p53 tetramerization domain results in an increased propensity to form amyloid fibrils

We describe in molecular detail how disruption of an intermonomer salt bridge (Arg337-Asp352) leads to partial destabilization of the p53 tetramerization domain and a dramatically increased propensity to form amyloid fibrils. At pH 4.0 and 37 degrees C, a p53 tetramerization domain mutant (p53tet-R337H), associated with adrenocortical carcinoma in children, readily formed amyloid fibrils, while the wild-type (p53tet-wt) did not. We characterized these proteins by equilibrium denaturation, 13C(alpha) secondary chemical shifts, (1H)-15N heteronuclear NOEs, and H/D exchange. Although p53tet-R337H was thermodynamically less stable, NMR data indicated that the two proteins had similar secondary structure and molecular dynamics. NMR derived pK(a) values indicated that at low pH the R337H mutation partially disrupted an intermonomer salt bridge. Backbone H/D exchange results showed that for at least a small population of p53tet-R337H molecules disruption of this salt bridge resulted in partial destabilization of the protein. It is proposed that this decrease in p53tet-R337H stability resulted in an increased propensity to form amyloid fibrils.

Stabilization of the fibrous structure of an α-helix-forming peptide by sequence reversal

The alpha3-peptide, which comprises three repeats of the sequence Leu-Glu-Thr-Leu-Ala-Lys-Ala and forms an amphipathic alpha-helix, is unique among various alpha-helix-forming peptides in that it assembles into fibrous structures that can be observed by transmission electron microscopy. As part of our investigation of the structure-stability relationships of the alpha3-peptide, we synthesized the r3-peptide, whose amino acid sequence is the reverse of that of the alpha3-peptide, and we investigated the effects of sequence reversal on alpha-helix stability and the formation of fibrous structures. Unexpectedly, the r3-peptide formed a more-stable alpha-helix and longer fibers than did the alpha3-peptide. The stability of the r3-peptide helix decreased when the ionic strength of the buffer was increased and when the pH of the buffer was adjusted to 2 or 12. These results suggest that the r3-peptide underwent a "magnet-like" oligomerization and that an increase in the charge-distribution inequality may be the driving force for the formation of fibrous structures.

Improved prediction for N-termini of alpha-helices using empirical information

The prediction of the secondary structure of proteins from their amino acid sequences remains a key component of many approaches to the protein folding problem. The most abundant form of regular secondary structure in proteins is the alpha-helix, in which specific residue preferences exist at the N-terminal locations. Propensities derived from these observed amino acid frequencies in the Protein Data Bank (PDB) database correlate well with experimental free energies measured for residues at different N-terminal positions in alanine-based peptides. We report a novel method to exploit this data to improve protein secondary structure prediction through identification of the correct N-terminal sequences in alpha-helices, based on existing popular methods for secondary structure prediction. With this algorithm, the number of correctly predicted alpha-helix start positions was improved from 30% to 38%, while the overall prediction accuracy (Q3) remained the same, using cross-validated testing. Although the algorithm was developed and tested on multiple sequence alignment-based secondary structure predictions, it was also able to improve the predictions of start locations by methods that use single sequences to make their predictions. Furthermore, the residue frequencies at N-terminal positions of the improved predictions better reflect those seen at the N-terminal positions of alpha-helices in proteins. This has implications for areas such as comparative modeling, where a more accurate prediction of the N-terminal regions of alpha-helices should benefit attempts to model adjacent loop regions. The algorithm is available as a Web tool, located at http://rocky.bms.umist.ac.uk/elephant.

The helical alanine controversy: an (Ala)6 insertion dramatically increases helicity

Employing chemical shift melts and hydrogen/deuterium exchange NMR techniques, we have determined the stabilization of the Trp-cage miniprotein due to multiple alanine insertions within the N-terminal alpha-helix. Alanine is shown to be uniquely helix-stabilizing and this stabilization is reflected in the global fold stability of the Trp-cage. The associated free energy change per alanine can be utilized to calculate the alanine propagation value. From the Lifson-Roig formulation, the calculated value (wAla = 1.6) is comparable to those obtained for short, solubilized, alanine-rich helices and is much larger than the values obtained by prior host-guest techniques or in N-terminally templated helices and peptides bearing long contiguous strings of alanines with no capping or solubilizing units present.

Effect of the N3 residue on the stability of the alpha-helix

N3 is the third position from the N terminus in the alpha-helix with helical backbone dihedral angles. All 20 amino acids have been placed in the N3 position of a synthetic helical peptide (CH(3)CO-[AAX AAAAKAAAAKAGY]-NH(2)) and the helix content measured by circular dichroism spectroscopy at 273 K. The dependence of peptide helicity on N3 residue identity has been used to determine a free energy scale by analysis with a modified Lifson-Roig helix coil theory that includes a parameter for the N3 energy (n3). The most stabilizing residues at N3 in rank order are Ala, Glu, Met/Ile, Leu, Lys, Ser, Gln, Thr, Tyr, Phe, Asp, His, and Trp. Free energies for the most destabilizing residues (Cys, Gly, Asn, Arg, and Pro) could not be fitted. The results correlate with N1, N2, and helix interior energies and not at all with N-cap preferences. This completes our work on studying the structural and energetic preferences of the amino acids for the N-terminal positions of the alpha-helix. These results can be used to rationally modify protein stability, help design helices, and improve prediction of helix location and stability.

Amino acid propensities are position-dependent throughout the length of alpha-helices

The 20 commonly occurring amino acids have been shown to have distinct position-dependent, helix-forming propensities near the ends of alpha-helices. Here, we show that the amino acids also have very strong position-dependent propensities throughout the length of a helix. Most helices are amphiphilic, and they have a strong tendency to both begin and end on the solvent-inaccessible face of the helix. These position-specific propensities should provide valuable parameters to guide de novo protein design, and should allow more precise prediction of helical topology in natural proteins.

De Novo Design of a Molecular Switch: Phosphorylation-dependent Association of Designed Peptides

The de novo design of peptides that switch their oligomerization state in response to a chemical stimulus is of interest, both as a tool for understanding the basis of molecular switching as well as development of reagents for the study of signal transduction in cells. The target of the current study is the design of a series of peptides that undergo a transition from an unstructured monomer to a four-helical bundle upon phosphorylation by the enzyme cyclic AMP-dependent protein kinase (PKA). The designed peptides are based on the 20-residue Lac repressor tetramerization domain. Beginning with this structure, we introduced a phosphorylation site near the N terminus. Phosphorylation leads to a 2-4.6 kcal/mol increase in the stability of the tetramer, depending on the design. The most successful switches were designed such that phosphorylation would increase the stability of the individual helices and also relieve an unfavorable electrostatic interaction in the tetramer.

Exploration of structural features of monomeric helical peptides designed with a genetic algorithm

A genetic algorithm (GA)-based strategy to dissect the determinants of peptide folding into alpha-helix was developed. The structural information of helical peptides was obtained with respect to patterns of sequence variability. In many previously reported studies the intrinsic alpha-helical propensities of amino acids although sequence-dependent are apparently independent of the amino acid position. In this research, monomeric helical peptides selected from possible sequences produced by a GA-chemical synthesis were analyzed to identify possible influential structural features. These hexadeca-peptides were obtained after four successive generations. A total of 128 synthetic peptides were evaluated via circular dichroism (CD) measurements in aqueous solution, while the mean ellipticity at 222 nm confirmed the monomeric state of the peptides. The results presented here show that our GA-based strategy may be useful in the design of proteins with increased alpha-helix content.

Importance of alpha-helix N-capping motif in stabilization of betabetaalpha fold

FSD-1 (full sequence design 1) is a protein folded in a betabetaalpha motif, designed on the basis of the second zinc finger domain of Zif268 by a substitution of its metal coordination site with a hydrophobic core. In this work, we analyzed the possibility of introducing the DNA recognition motif of the template zinc finger (S(13)RSDH(17)) into FSD-1 sequence in order to obtain a small DNA-binding module devoid of cross-link(s) or metal cofactors. The hybrid protein was unfolded, as judged by CD and NMR criteria. To reveal the role of each of the five amino acids, which form the N-capping motif of the alpha-helix, we analyzed conformational and stability properties of eight FSD-1 mutants. We used a shielded methyl group of Leu 18 and a CD signal at 215 nm as a convenient measure of the folded state. Glu 17-->His substitution at the N(3) in S(13)NEKE(17) variant decreased the folded structure content from 90% to 25% (equivalent to 1.8 kcal * mole(-1) destabilization) by disruption of N-capping interactions, and had the most significant effect among single mutants studied here. The N(cap) Asn 14 substitution with Arg considerably decreased stability, reducing structure content from 90% to 40% (1.4 kcal * mole(-1) destabilization) by disruption of a helix-capping hydrogen bond and destabilization of a helix macrodipole. The N(1) Glu 15-->Ser mutation also produced a considerable effect (1.0 kcal * mole(-1) destabilization), again emphasizing the significance of electrostatic interactions in alpha-helix stabilization.

Exceptional pairs of amino acid neighbors in alpha-helices

Amino acids seem to have specific preferences for various locations in alpha-helices. These specific preferences, called singlet local propensity (SLP), have been determined by calculating the preference of occurrence of each amino acid in different positions of the alpha-helix. We have studied the occurrence of amino acids, single or pairs, in different positions, singlet or doublet, of alpha-helices in a database of 343 non-homologous proteins representing a unique superfamily from the SCOP database with a resolution better than 2.5 A from the Protein Data Bank. The preference of single amino acids for various locations of the helix was shown by the relative entropy of each amino acid with respect to the background. Based on the total relative entropy of all amino acids occurring in a single position, the N(cap) position was found to be the most selective position in the alpha-helix. A rigorous statistical analysis of amino acid pair occurrences showed that there are exceptional pairs for which, the observed frequency of occurrence in various doublet positions of the alpha-helix is significantly different from the expected frequency of occurrence in that position. The doublet local propensity (DLP) was defined as the preference of occurrences of amino acid pairs in different doublet positions of the alpha-helix. For most amino acid pairs, the observed DLP (DLP(O)) was nearly equal to the expected DLP (DLP(E)), which is the product of the related SLPs. However, for exceptional pairs of amino acids identified above, the DLP(O) and DLP(E) values were significantly different. Based on the relative values of DLP(O) and DLP(E), exceptional amino acid pairs were divided into two categories. Those, for which the DLP(O) values are higher than DLP(E), should have a strong tendency to pair together in the specified position. For those pairs which the DLP(O) values are less than DLP(E), there exists a hindrance in neighboring of the two amino acids in that specific position of the alpha-helix. These cases have been identified and listed in various tables in this paper. The amount of mutual information carried by the exceptional pairs of amino acids was significantly higher than the average mutual information carried by other amino acid pairs. The average mutual information conveyed by amino acid pairs in each doublet position was found to be very small but non-zero.

Recent advances in helix-coil theory

Peptide helices in solution form a complex mixture of all helix, all coil or, most frequently, central helices with frayed coil ends. In order to interpret experiments on helical peptides and make theoretical predictions on helices, it is therefore essential to use a helix-coil theory that takes account of this equilibrium. The original Zimm-Bragg and Lifson-Roig helix-coil theories have been greatly extended in the last 10 years to include additional interactions. These include preferences for the N-cap, N1, N2, N3 and C-cap positions, capping motifs, helix dipoles, side chain interactions and 3(10)-helix formation. These have been applied to determine energies for these preferences from experimental data and to predict the helix contents of peptides. This review discusses these newly recognised structural features of helices and how they have been included in helix-coil models.

Helix-coil transitions re-visited

The thermally-induced helix-coil transition in polyamino acids is a good model for determining the helix-forming propensities of amino acids but not for the two-state folding/unfolding transition in globular proteins. The equilibrium and kinetic treatments of the helix-coil transition are summarized here together with a description of applications to various types of homopolymers and copolymers. Attention is then focused on the helix-coil transition in poly-L-alanine as an example of a non-polar polyamino acid. To render such a non-polar polymer water soluble, it is necessary to introduce polar amino acids such as lysines, but care must be taken as to the location of such polar residues. If they are attached as end groups, as in a triblock copolymer, they do not perturb the helix-forming tendency of the central poly-L-alanine block significantly, but if they are introduced within the sequence of alanine residues, then the hydration properties of the lysines dominate the behavior of the resulting copolymer, thereby leading to erroneous values of the parameters characterizing the helix-forming tendency of the alanines. Neutral but polar residues, such as glutamines, also exhibit hydration-dominating properties but less so than charged lysines. Some details of the calculations for an alanine/glutamine copolymer are presented here. It is concluded that random copolymers based on a neutral water-soluble host provide reliable information about the helix-forming tendencies of amino acid residues that are introduced as guests among such neutral host residues.

A critical assessment of the secondary structure alpha-helices and their termini in proteins

Secondary structure prediction from amino acid sequence is a key component of protein structure prediction, with current accuracy at approximately 75%. We analysed two state-of-the-art secondary structure prediction methods, PHD and JPRED, comparing predictions with secondary structure assigned by the algorithms DSSP and STRIDE. The specific focus of our study was alpha-helix N-termini, as empirical free energy scales are available for residue preferences at N-terminal positions. Although these prediction methods perform well in general at predicting the alpha-helical locations and length distributions in proteins, they perform less well at predicting the correct helical termini. For example, although most predicted alpha-helices overlap a real alpha-helix (with relatively few completely missed or extra predicted helices), only one-third of JPRED and PHD predictions correctly identify the N-terminus. Analysis of neighbouring N-terminal sequences to predicted helical N-termini shows that the correct N-terminus is often within one or two residues. More importantly, the true N-terminal motif is, on average, more favourable as judged by our experimentally measured free energies. This suggests a simple, but powerful, strategy to improve secondary structure prediction using empirically derived energies to adjust the predicted output to a more favourable N-terminal sequence.

Amino acid intrinsic alpha-helical propensities III: positional dependence at several positions of C terminus

In this study, we have analyzed experimentally the helical intrinsic propensities of non-charged and non-aromatic residues at different C-terminal positions (C1, C2, C3) of an Ala-based peptide. The effect was found to be complex, resulting in extra stabilization or destabilization, depending on guest amino acid and position under consideration. Polar (Ser, Thr, Cys, Asn, and Gln) amino acids and Gly were found to have significantly larger helical propensities at several C-terminal positions compared with the alpha-helix center (-1.0 kcal/mole in some cases). Some of the nonpolar residues, especially beta-branched ones (Val and Ile) are significantly more favorable at position C3 (-0.3 to -0.4 kcal/mole), although having minor differences at other C-terminal positions compared with the alpha-helix center. Leu has moderate (-0.1 to -0.2 kcal/mole) stabilization effects at position C2 and C3, whereas being relatively neutral at C1. Finally, Met was found to be unfavorable at C1 and C2 ( +0.2 kcal/mole) and favorable at C3 (-0.2 kcal/mole). Thus, significant differences found between the intrinsic helical propensities at the C-terminal positions and those in the alpha-helix center must be accounted for in helix/coil transition theories and in protein design.