Amino acid distribution in protein secondary structures

Enhancement of sludge dewaterability using combined technology of bioleaching and Fenton: Microscopic structure and hydrophilic/hydrophobic properties of sludge particles

Bioleaching and Fenton technology are commonly used preconditioning techniques for sludge dewatering. This study compared the dewatering mechanisms of different conditioning technologies. The results showed that bound water, specific resistance to filtration (SRF), and capillary suction time decreased from 3.95 g/g, 6.16 × 1012 m/kg, and 130.6 s to 3.15 g/g, 2.81 × 1011 m/kg, and 33 s, respectively, under combined treatment condition. Moreover, the free radicals, including ·OH, O2-·and Fe (Ⅳ), further damaged the cell structure, thus increasing the concentration of DNA in the S-EPS layer. This intense degradation sludge particle size decreased by 15.6% and significantly increased zeta potential. Under the combined technology, the α-helix and β-sheet decreased by 42.2% and 56.5%, respectively, destabilizing the spatial structure of proteins and promoting the release of bound water. In addition, the combined technology decreased (Ala/Lys) ratio in the TB-EPS layer by 67.6%, indicating the weakening of protein water-holding capacity. Moreover, the conversion of oxygen-containing compounds to nonpolar hydrocarbons increased the hydrophobicity of the sludge under a combined treatment, thus enhancing dewatering performance.

The effect of the genomic GC content bias of prokaryotic organisms on the secondary structures of their proteins

One of the main characteristics of prokaryotic genomes is the ratio in which guanine-cytosine bases are used in their DNA sequences. This is known as the genomic GC content and varies widely, from values below 20% to values greater than 74%. It has been demonstrated that the genomic GC content varies in accordance with the phylogenetic distribution of organisms and influences the amino acid composition of their corresponding proteomes. This bias is particularly important for amino acids that are coded by GC content-rich codons such as alanine, glycine, and proline, as well as amino acids that are coded by AT-rich codons, such as lysine, asparagine, and isoleucine. In our study, we extend these results by considering the effect of the genomic GC content on the secondary structure of proteins. On a set of 192 representative prokaryotic genomes and proteome sequences, we identified through a bioinformatic study that the composition of the secondary structures of the proteomes varies in relation to the genomic GC content; random coils increase as the genomic GC content increases, while alpha-helices and beta-sheets present an inverse relationship. In addition, we found that the tendency of an amino acid to form part of a secondary structure of proteins is not ubiquitous, as previously expected, but varies according to the genomic GC content. Finally, we discovered that for some specific groups of orthologous proteins, the GC content of genes biases the composition of secondary structures of the proteins for which they code.

An electro-peroxone oxidation-Fe(III) coagulation sequential conditioning process for the enhanced waste activated sludge dewatering: Bound water release and organics multivariate change

As a by-product of wastewater treatment, waste activated sludge (WAS) has complex composition, strong hydrophilic extracellular polymeric substance (EPS), which make it difficult to dewater. In this study, an electro-peroxone oxidation-Fe(III) coagulation (E-peroxone-Fe(III)) sequential conditioning approach was developed to improve WAS dewaterability. At E-peroxone oxidation stage, hydrogen peroxide was generated through 2-electron path on a carbon polytetrafluoroethylene cathode, and reacted with the sparged O₃ to produce hydroxyl radicals. At the subsequent coagulation stage, Fe(III) was dosed to coagulate the small WAS fragments and release water from WAS. Along E-peroxone-Fe(III) subsequent conditioning process, the physicochemical properties of WAS, main components, functional groups and evolution of protein secondary structure, and typical amino acids in EPS, as well as the type and semi-quantitative of elements in WAS, were investigated. The results indicated that under the optimal conditions, the reductions of specific resistance to filterability (SRF) and capillary suction time (CST) for WAS equalled 78.18% and 71.06%, respectively, and its bound water content decreased from 8.87 g/g TSS to 7.67 g/g TSS. After E-peroxone oxidation, part of protein and polysaccharide migrated outside from TB-EPS to slime, the ratio of α-helix/(β-sheet + random coil) declined, even some of organic-N disintegrated to inorganic-N. At Fe(III) coagulation stage, re-coagulation of the dispersed WAS fragments and easy extraction from inner EPS for protein and polysaccharide occurred. Furthermore, the protein secondary structure of β-sheet increased by 13.48%, the contents of hydrophobic and hydrophilic amino acids also increased. In addition, a strong negative correlation between the hydrophobic amino acid content of Met in slime and CST or SRF (R²_CST = -0.999, p < 0.05 or R²_SRF = -0.948, p < 0.05) occurred, while a strong positive correlation between the hydrophilic amino acid content of Cys in TB-EPS and CST or SRF (R²_CST = 0.992, p < 0.05 or R²_SRF = 0.921, p < 0.05) occurred, which could be related to the WAS dewaterability.

Profiling of amino acids and their interactions with proteinaceous compounds for sewage sludge dewatering by Fenton oxidation treatment

During advanced oxidation treatment for enhancing sludge dewaterability, the peptide chains of protein can be decomposed into amino acids. Protein exhibits a great impact on sewage sludge dewaterability. However, the role of amino acids in sludge dewatering remains unclear. In this study, among the 23 types of amino acids investigated, tryptophane (Trp) and lysine (Lys) were identified as the key amino acids affecting sludge dewaterability during Fenton oxidation treatment. The content of lysine showed positive correlations with capillary suction time (CST), specific resistance to filtration (SRF), and bound water content, and the concentrations of total protein, low molecular weight protein, amines and amides, and 3-turn helix of proteinaceous compounds in bound extracellular polymeric substances (EPS), while the content of tryptophane showed negative correlations with the above parameters. The amino acids may be sourced from damage of the membrane and ribosomal proteins by hydroxyl radicals, and the peptide bonds connected with tryptophane were more inclined to be decomposed than other amino acids. Particularly, more amino acids of tryptophane can result in more hydrophobic interaction, and less necessary energy barrier for aggregation of particles. As such, regulating protein degradation towards production of tryptophane may be related with enhanced sludge dewaterability by Fenton oxidation treatment.

Structure determination of UL49.5 transmembrane protein from bovine herpesvirus 1 by NMR spectroscopy and molecular dynamics

The transporter associated with antigen processing (TAP) directly participates in the immune response as a key component of the cytosolic peptide to major histocompatibility complex (MHC) class I protein loading machinery. This makes TAP an important target for viruses avoiding recognition by CD8+ T lymphocytes. Its activity can be suppressed by the UL49.5 protein produced by bovine herpesvirus 1, although the mechanism of this inhibition has not been understood so far.

Therefore, the main goal of our study was to investigate the 3D structure of bovine herpesvirus 1 - encoded UL49.5 protein. The final structure of the inhibitor was established using circular dichroism (CD), 2D nuclear magnetic resonance (NMR), and molecular dynamics (MD) in membrane mimetic environments. In NMR studies, UL49.5 was represented by two fragments: the extracellular region (residues 1–35) and the transmembrane-intracellular fragment (residues 36–75), displaying various functions during viral invasion. After the empirical structure determination, a molecular docking procedure was used to predict the complex of UL49.5 with the TAP heterodimer.

Our results revealed that UL49.5 adopted a highly flexible membrane-proximal helical structure in the extracellular part. In the transmembrane region, we observed two short α-helices. Furthermore, the cytoplasmic part had an unordered structure. Finally, we propose three different orientations of UL49.5 in the complex with TAP. Our studies provide, for the first time, the experimental structural information on UL49.5 and structure-based insight in its mechanism of action which might be helpful in designing new drugs against viral infections.

•

The UL49.5 viral protein forms a helical structure in the biological membrane

•

Our NMR-based 3D structure of UL49.5 differs from the theoretical predictions

•

Apart from the protruding N-terminal helix the structure is buried in the membrane

•

Attention should be paid to the turns in the external and transmembrane domains

•

Molecular docking proposes three possible structures of the UL49.5/TAP complexes

Helical ambivalency induced by point mutations

Background

Mutation of amino acid sequences in a protein may have diverse effects on its structure and function. Point mutations of even a single amino acid residue in the helices of the non-redundant database may lead to sequentially identical peptides which adopt different secondary structures in different proteins. However, various physico-chemical factors which govern the formation of these ambivalent helices generated by point mutations of a sequence are not clearly known.

Results

Sequences generated by point mutations of helices are mapped on to their non-helical counterparts in the SCOP database. The results show that short helices are prone to transform into non-helical conformations upon point mutations. Mutation of amino acid residues by helix breakers preferentially yield non-helical conformations, while mutation with residues of intermediate helix propensity display least preferences for non-helical conformations. Differences in the solvent accessibility of the mutating/mutated residues are found to be a major criteria for these sequences to conform to non-helical conformations. Even with minimal differences in the amino acid distributions of the sequences flanking the helical and non-helical conformations, helix-flanking sequences are found be more solvent accessible.

Conclusions

All types of mutations from helical to non-helical conformations are investigated. The primary factors attributing such changes in conformation can be: i) type/propensity of the mutating and mutant residues ii) solvent accessibility of the residue at the mutation site iii) context/environment dependence of the flanking sequences. The results from the present study may be used to design de novo proteins via point mutations.

(φ,ψ)₂ motifs: a purely conformation-based fine-grained enumeration of protein parts at the two-residue level

A deep understanding of protein structure benefits from the use of a variety of classification strategies that enhance our ability to effectively describe local patterns of conformation. Here, we use a clustering algorithm to analyze 76,533 all-trans segments from protein structures solved at 1.2 Å resolution or better to create a purely φ,ψ-based comprehensive empirical categorization of common conformations adopted by two adjacent φ,ψ pairs (i.e., (φ,ψ)(2) motifs). The clustering algorithm works in an origin-shifted four-dimensional space based on the two φ,ψ pairs to yield a parameter-dependent list of (φ,ψ)(2) motifs, in order of their prominence. The results are remarkably distinct from and complementary to the standard hydrogen-bond-centered view of secondary structure. New insights include an unprecedented level of precision in describing the φ,ψ angles of both previously known and novel motifs, ordering of these motifs by their population density, a data-driven recommendation that the standard C(α(i))…C(α(i+3))<7 Å criteria for defining turns be changed to 6.5 Å, identification of β-strand and turn capping motifs, and identification of conformational capping by residues in polypeptide II conformation. We further document that the conformational preferences of a residue are substantially influenced by the conformation of its neighbors, and we suggest that accounting for these dependencies will improve protein modeling accuracy. Although the CUEVAS-4D(r(10)є(14)) 'parts list' presented here is only an initial exploration of the complex (φ,ψ)(2) landscape of proteins, it shows that there is value to be had from this approach, and it opens the door to more in-depth characterizations at the (φ,ψ)(2) level and at higher dimensions.

A polypeptide "building block" for the β-trefoil fold identified by "top-down symmetric deconstruction"

Fibroblast growth factor-1, a member of the 3-fold symmetric β-trefoil fold, was subjected to a series of symmetric constraint mutations in a process termed "top-down symmetric deconstruction." The mutations enforced a cumulative exact 3-fold symmetry upon symmetrically equivalent positions within the protein and were combined with a stability screen. This process culminated in a β-trefoil protein with exact 3-fold primary-structure symmetry that exhibited excellent folding and stability properties. Subsequent fragmentation of the repeating primary-structure motif yielded a 42-residue polypeptide capable of spontaneous assembly as a homotrimer, producing a thermostable β-trefoil architecture. The results show that despite pronounced reduction in sequence complexity, pure symmetry in the design of a foldable, thermostable β-trefoil fold is possible. The top-down symmetric deconstruction approach provides a novel alternative means to successfully identify a useful polypeptide "building block" for subsequent "bottom-up" de novo design of target protein architecture.

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.

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.

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.

Hidden Markov model-derived structural alphabet for proteins: the learning of protein local shapes captures sequence specificity

Understanding and predicting protein structures depend on the complexity and the accuracy of the models used to represent them. We have recently set up a Hidden Markov Model to optimally compress protein three-dimensional conformations into a one-dimensional series of letters of a structural alphabet. Such a model learns simultaneously the shape of representative structural letters describing the local conformation and the logic of their connections, i.e. the transition matrix between the letters. Here, we move one step further and report some evidence that such a model of protein local architecture also captures some accurate amino acid features. All the letters have specific and distinct amino acid distributions. Moreover, we show that words of amino acids can have significant propensities for some letters. Perspectives point towards the prediction of the series of letters describing the structure of a protein from its amino acid sequence.

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.

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.

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.

The interrelationships of side-chain and main-chain conformations in proteins

The accurate determination of a large number of protein structures by X-ray crystallography makes it possible to conduct a reliable statistical analysis of the distribution of the main-chain and side-chain conformational angles, how these are dependent on residue type, adjacent residue in the sequence, secondary structure, residue-residue interactions and location at the polypeptide chain termini. The interrelationship between the main-chain (phi, psi) and side-chain (chi 1) torsion angles leads to a classification of amino acid residues that simplify the folding alphabet considerably and can be a guide to the design of new proteins or mutational studies. Analyses of residues occurring with disallowed main-chain conformation or with multiple conformations shed some light on why some residues are less favoured in thermophiles.

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

N1 is the first residue in an alpha-helix. We have measured the contribution of all 20 amino acids to the stability of a small helical peptide CH(3)CO-XAAAAQAAAAQAAGY-NH(2) at the N1 position. By substituting every residue into the N1 position, we were able to investigate the stabilizing role of each amino acid in an isolated context. The helix content of each of the 20 peptides was measured by circular dichroism (CD) spectroscopy. The data were analyzed by our modified Lifson-Roig helix-coil theory, which includes the n1 parameter, to find free energies for placing a residue into the N1 position. The rank order for free energies is Asp(-), Ala > Glu(-) > Glu(0) > Trp, Leu, Ser > Asp(0), Thr, Gln, Met, Ile > Val, Pro > Lys(+), Arg, His(0) > Cys, Gly > Phe > Asn, Tyr, His(+). N1 preferences are clearly distinct from preferences for the preceding N-cap and alpha-helix interior. pK(a) values were measured for Asp, Glu, and His, and protonation-free energies were calculated for Asp and Glu. The dissociation of the Asp proton is less favorable than that of Glu, and this reflects its involvement in a stronger stabilizing interaction at the N terminus. Proline is not energetically favored at the alpha-helix N terminus despite having a high propensity for this position in crystal structures. The data presented are of value both in rationalizing mutations at N1 alpha-helix sites in proteins and in predicting the helix contents of peptides.

Determination of alpha-helix N1 energies after addition of N1, N2, and N3 preferences to helix/coil theory

Surveys of protein crystal structures have revealed that amino acids show unique structural preferences for the N1, N2, and N3 positions in the first turn of the alpha-helix. We have therefore extended helix-coil theory to include statistical weights for these locations. The helix content of a peptide in this model is a function of N-cap, C-cap, N1, N2, N3, C1, and helix interior (N4 to C2) preferences. The partition function for the system is calculated using a matrix incorporating the weights of the fourth residue in a hexamer of amino acids and is implemented using a FORTRAN program. We have applied the model to calculate the N1 preferences of Gln, Val, Ile, Ala, Met, Pro, Leu, Thr, Gly, Ser, and Asn, using our previous data on helix contents of peptides Ac-XAKAAAAKAAGY-CONH2. We find that Ala has the highest preference for the N1 position. Asn is the most unfavorable, destabilizing a helix at N1 by at least 1.4 kcal mol(-1) compared to Ala. The remaining amino acids all have similar preferences, 0.5 kcal mol(-1) less than Ala. Gln, Asn, and Ser, therefore, do not stabilize the helix when at N1.

Hidden Markov model approach for identifying the modular framework of the protein backbone

The hidden Markov model (HMM) was used to identify recurrent short 3D structural building blocks (SBBs) describing protein backbones, independently of any a priori knowledge. Polypeptide chains are decomposed into a series of short segments defined by their inter-alpha-carbon distances. Basically, the model takes into account the sequentiality of the observed segments and assumes that each one corresponds to one of several possible SBBs. Fitting the model to a database of non-redundant proteins allowed us to decode proteins in terms of 12 distinct SBBs with different roles in protein structure. Some SBBs correspond to classical regular secondary structures. Others correspond to a significant subdivision of their bounding regions previously considered to be a single pattern. The major contribution of the HMM is that this model implicitly takes into account the sequential connections between SBBs and thus describes the most probable pathways by which the blocks are connected to form the framework of the protein structures. Validation of the SBBs code was performed by extracting SBB series repeated in recoding proteins and examining their structural similarities. Preliminary results on the sequence specificity of SBBs suggest promising perspectives for the prediction of SBBs or series of SBBs from the protein sequences.

Periodicity in alpha-helix lengths and C-capping preferences

We surveyed 299 high resolution, non-homologous protein crystal structures for alpha-helix lengths and capping preferences. We find that helices show a preference to have close to an integral number of turns. Helices can be usefully subdivided into either "favoured length" with 6, 7, 10, 11, 13, 14, 17, 18, 21, 22, 24, 25, 28, 29 or 31 residues, or "disfavoured length" with 8, 9, 12, 15, 16, 19, 20, 23, 26, 27 or 30 residues. Favoured length helices have their N and C-caps on the same side of the helix so they can lie on the protein surface. There is no significant difference in amino acid preferences at the N terminus between favoured and disfavoured length helices. At the C terminus, favoured length helices prefer non-polar side-chains at C4 and polar amino acid residues at C2, while disfavoured length helices prefer non-polar amino acid residues at C2. There are strong periodic trends in the likelihood of terminating a helix with a Schellman or alphaL C-capping motif. These can be rationalised by the preference for a non-polar side-chain at C3 with these motifs, favouring placing C3 on the buried side of the helix. We suggest that algorithms aiming to predict helices or C-capping in proteins should include a weight for the helix length.

Side-chain structures in the first turn of the alpha-helix

The first three residues at the N terminus of the alpha-helix are called N1, N2 and N3. We surveyed 2102 alpha-helix N termini in 298 high-resolution, non-homologous protein crystal structures for N1, N2 and N3 amino acid and side-chain rotamer propensities and hydrogen-bonding patterns. We find strong structural preferences that are unique to these sites. The rotamer distributions as a function of amino acid identity and position in the helix are often explained in terms of hydrogen-bonding interactions to the free N1, N2 and N3 backbone NH groups. Notably, the "good N2" amino acid residues Gln, Glu, Asp, Asn, Ser, Thr and His preferentially form i, i or i,i+1 hydrogen bonds to the backbone, though this is hindered by good N-caps (Asp, Asn, Ser, Thr and Cys) that compete for these hydrogen bond donors. We find a number of specific side-chain to side-chain interactions between N1 and N2 or between the N-cap and N2 or N3, such as Arg(N-cap) to Asp(N2). The strong energetic and structural preferences found for N1, N2 and N3, which differ greatly from positions within helix interiors, suggest that these sites should be treated explicitly in any consideration of helical structure in peptides or proteins.

Geometrical and sequence characteristics of alpha-helices in globular proteins

Understanding the sequence-structure relationships in globular proteins is important for reliable protein structure prediction and de novo design. Using a database of 1131 alpha-helices with nonidentical sequences from 205 nonhomologous globular protein chains, we have analyzed structural and sequence characteristics of alpha-helices. We find that geometries of more than 99% of all the alpha-helices can be simply characterised as being linear, curved, or kinked. Only a small number of alpha-helices ( approximately 4%) show sharp localized bends in their middle regions, and thus are classified as kinked. Approximately three-fourths (approximately 73%) of the alpha-helices in globular proteins show varying degrees of smooth curvature, with a mean radius of curvature of 65 +/- 33 A; longer helices are less curved. Computation of helix accessibility to the solvent indicates that nearly two-thirds of the helices ( approximately 66%) are largely buried in the protein core, and the length and geometry of the helices are not correlated with their location in the protein globule. However, the amino acid compositions and propensities of individual amino acids to occur in alpha-helices vary with their location in the protein globule, their geometries, and their lengths. In particular, Gln, Glu, Lys, and Arg are found more often in helices near the surface of globular proteins. Interestingly, kinks often seem to occur in regions where amino acids with low helix propensities (e.g., beta-branched and aromatic residues) cluster together, in addition to those associated with the occurrence of proline residues. Hence the propensities of individual amino acids to occur in a given secondary structure depend not only on conformation but also on its length, geometry, and location in the protein globule.

Dissecting alpha-helices: position-specific analysis of alpha-helices in globular proteins

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

Deciphering the structural code for proteins: helical propensities in domain classes and statistical multiresidue information in alpha-helices

We made several statistical analyses in a large sample of nearly 4,000 helices (from 546 redundancy-controlled PDB protein subunits), which give new insights into the helical properties of globular proteins. In a first experiment, the amino acid composition of the whole sample was compared with the composition of two helical sample subgroups (the "mainly-alpha" and the "(alpha/beta)8 barrel" domain classes); we reached the conclusion that composition-based helical propensities for secondary structure prediction do not depend on the structural class. Running a five-residue window through the whole sample, the positional composition revealed that positive and negative residues are located throughout the helices and tend to neutralize the macrodipole effect. On this basis, we analyzed charged triplets using a running five-residue window. The conclusion was that only mixed charged residues [positive (+) and negative (-)] located at positions 1-2-5 and 1-4-5 are clearly favored. In these locations the most abundant are (- -..+) and (-..+ +), and this shows the existence of side chain microdipoles, which neutralize the large macrodipole of the helix. We made a systematic statistical analysis of charged, dipolar, and hydrophobic + aromatic residues, which enabled us to work out rules that should be useful for modeling and design purposes. Finally, we analyzed the relative abundance of all the different amphipathic double-arcs that are present in helices formed by octapeptides (8) and nonapeptides (18). All of the double-arcs that make up Schiffer and Edmundson's classical helical wheel are found in abundance in the sample.

The refined crystal structure of Bacillus cereus oligo-1,6-glucosidase at 2.0 A resolution: structural characterization of proline-substitution sites for protein thermostabilization

The crystal structure of oligo-1,6-glucosidase (dextrin 6-alpha-glucanohydrolase, EC 3.2.1.10) from Bacillus cereus ATCC7064 has been refined to 2.0 A resolution with an R-factor of 19.6% for 43,328 reflections. The final model contains 4646 protein atoms and 221 ordered water molecules with respective root-mean-square deviations of 0.015 A for bond lengths and of 3.166 degrees for bond angles from the ideal values. The structure consists of three domains: the N-terminal domain (residues 1 to 104 and 175 to 480), the subdomain (residues 105 to 174) and the C-terminal domain (residues 481 to 558). The N-terminal domain takes a (beta/alpha)8-barrel structure with additions of an alpha-helix (N alpha6') between the sixth strand Nbeta6 and the sixth helix N alpha6, an alpha-helix (N alpha7') between the seventh strand Nbeta7 and the seventh helix N alpha7 and three alpha-helices (N alpha8', N alpha8" and N alpha8'" between the eighth strand Nbeta8 and the eighth helix N alpha8. The subdomain is composed of an alpha-helix, a three-stranded antiparallel beta-sheet, and long intervening loops. The C-terminal domain has a beta-barrel structure of eight antiparallel beta-strands folded in double Greek key motifs, which is distorted in the sixth strand Cbeta6. Three catalytic residues, Asp199, Glu255 and Asp329, are located at the bottom of a deep cleft formed by the subdomain and a cluster of the two additional alpha-helices N alpha8' and N alpha8" in the (beta/alpha)8-barrel. The refined structure reveals the locations of 21 proline-substitution sites that are expected to be critical to protein thermostabilization from a sequence comparison among three Bacillus oligo-1,6-glucosidases with different thermostability. These sites lie in loops, beta-turns and alpha-helices, in order of frequency, except for Cys515 in the fourth beta-strand Cbeta4 of the C-terminal domain. The residues in beta-turns (Lys121, Glu208, Pro257, Glu290, Pro443, Lys457 and Glu487) are all found at their second positions, and those in alpha-helices (Asn109, Glu175, Thr261 and Ile403) are present at their N1 positions of the first helical turns. Those residues in both secondary structures adopt phi and phi values favorable for proline substitution. Residues preceding the 21 sites are mostly conserved upon proline occurrence at these 21 sites in more thermostable Bacillus oligo-1,6-glucosidases. These structural features with respect to the 21 sites indicate that the sites in beta-turns and alpha-helices have more essential prerequisites for proline substitution to thermostabilize the protein than those in loops. This well supports the previous finding that the replacement at the appropriate positions in beta-turns or alpha-helices is the most effective for protein thermostabilization by proline substitution.

Structures of N-termini of helices in proteins

We have surveyed 393 N-termini of alpha-helices and 156 N-termini of 3(10)-helices in 85 high resolution, non-homologous protein crystal structures for N-cap side-chain rotamer preferences, hydrogen bonding patterns, and solvent accessibilities. We find very strong rotamer preferences that are unique to N-cap sites. The following rules are generally observed for N-capping in alpha-helices: Thr and Ser N-cap side chains adopt the gauche - rotamer, hydrogen bond to the N3 NH and have psi restricted to 164 +/- 8 degrees. Asp and Asn N-cap side chains either adopt the gauche - rotamer and hydrogen bond to the N3 NH with psi = 172 +/- 10 degrees, or adopt the trans rotamer and hydrogen bond to both the N2 and N3 NH groups with psi = 1-7 +/- 19 degrees. With all other N-caps, the side chain is found in the gauche + rotamer so that the side chain does not interact unfavorably with the N-terminus by blocking solvation and psi is unrestricted. An i, i + 3 hydrogen bond from N3 NH to the N-cap backbone C = O in more likely to form at the N-terminus when an unfavorable N-cap is present. In the 3(10)-helix Asn and Asp remain favorable N-caps as they can hydrogen bond to the N2 NH while in the trans rotamer; in contrast, Ser and Thr are disfavored as their preferred hydrogen bonding partner (N3 NH) is inaccessible. This suggests that Ser is the optimum choice of N-cap when alpha-helix formation is to be encouraged while 3(10)-helix formation discouraged. The strong energetic and structural preferences found for N-caps, which differ greatly from positions within helix interiors, suggest that N-caps should be treated explicitly in any consideration of helical structure in peptides or proteins.

Non-random usage of 'degenerate' codons is related to protein three-dimensional structure

We report an analysis of a novel sequence-structure database of mammalian proteins incorporating nucleotide sequences of the exon regions of their genes together with protein sequence and structural information. We find that synonymous codon families (i.e. coding the same residue) have non-random codon distribution frequencies between protein secondary structure types. Their structural preferences are related to the third, 'silent' nucleotide position in a codon. We also find that some synonymous codons show very different or even opposite structural preferences at the N- or C-termini of structure fragments, relative to those observed for their amino acid residues.

Patterns in ionizable side chain interactions in protein structures

In a selected set of 44 high-resolution, non-homologous protein structures, the intramolecular hydrogen bonds or salt bridges formed by ionizable amino acid side chains were identified and analyzed. The analysis was based on the investigation of several properties of the involved residues such as their solvent exposure, their belonging to a certain secondary structural element, and their position relative to the N- and C-termini of their respective structural element. It was observed that two-thirds of the interactions made by basic or acidic side chains are hydrogen bonds to polar uncharged groups. In particular, the majority (78%) of the hydrogen bonds between ionizable side chains and main chain polar groups (sch:mch bonds) involved at least one buried atom, and in 42% of the cases both interacting atoms were buried. In alpha-helices, the sch:mch bonds observed in the proximity of the C- and N-termini show a clear preference for acidic and basic side chains, respectively. This appears to be due to the partial charges of peptide group atoms at the termini of alpha-helices, which establish energetically favorable electrostatic interactions with side chain carrying opposite charge, at distances even greater than 4.5 angstrom. The sch:mch interactions involving ionizable side chains that belong either to beta-strands or to the central part of alpha-helices are based almost exclusively on basic residues. This results from the presence of main chain carbonyl oxygen atoms in the protein core which have unsatisfied hydrogen bonding capabilities.

The alpha-helix as seen from the protein tertiary structure: a 3-D structural classification

Helices are selected from globular protein structures defined at high resolution by X-ray analysis. We cluster alpha-helices in two ways: according to their position in the tertiary structure by considering patterns of solvent inaccessible residues and according to the arc of the solvent inaccessible face. For each class of helices we have defined propensities for amino acids at each position; these can be used to calculate templates for recognition of a member of that class. The analysis provides a basis for the prediction of alpha-helices and estimation of their approximate position in a protein tertiary structure. It also provides an approach to estimating the probability of finding amino acid sequences as helices in solution and in a folded protein, thus indicating those helices that might be involved in nucleation of protein folding.

Multiple proline substitutions cumulatively thermostabilize Bacillus cereus ATCC7064 oligo-1,6-glucosidase. Irrefragable proof supporting the proline rule

Nine residues of Bacillus cereus ATCC7064 oligo-1,6-glucosidase were replaced stepwise with proline residues. Of the nine residues, Lys121, Glu208 and Glu290 were at second sites of beta turns; Asn109, Glu175 and Thr261 were at N-caps of alpha helices; Glu216, Glu270 and Glu378 were in coils within loops. The replacements were carried out in the order, Lys121-->Pro, Glu175-->Pro, Glu290-->Pro, Glu208-->Pro, Glu270-->Pro, Glu378-->Pro, Thr261-->Pro, Glu216-->Pro and Asn109-->Pro. The resultant nine active mutant enzymes contained 1-9 more proline residues than B. cereus oligo-1,6-glucosidase. The thermostability of these mutants was additively enhanced with the increase in the number of proline residues introduced. The increase in the thermostability was most remarkable when proline residues were introduced at second sites of beta turns or at N-caps of alpha helices. The above results afforded irrefragable proof for the proline rule as an effective principle for increasing protein thermostability [Suzuki, Y., Oishi, K., Nakano, H. & Nagayama, T. (1987) Appl. Microbiol. Biotechnol. 26, 546-551].

Helix stop signals in proteins and peptides: the capping box

The alpha-helix [Pauling, L., Corey, R. B., & Branson, H. R. (1951) Proc. Natl. Acad. Sci. U.S.A. 37, 205-211] is a common motif in both proteins and peptides. Despite intense investigation, predictive understanding of helices is still lacking. A recent hypothesis [Presta, L. G., & Rose, G. D. (1988) Science 240, 1632-1641] proposed that the structural specificity of helices resides, in part, in those residues that flank helix termini. If so, then signals that arrest helix propagation--i.e., helix stop signals--should be found among these flanking residues. Evidence is presented for the existence of one such signal, a reciprocal backbone-side-chain hydrogen-bonding interaction, dubbed the capping box. In proteins, the capping box is found uniquely at helix N-termini. In peptides, the capping box can function as a helix stop signal, as shown in the work of Kallenbach and co-workers.

Studies of the zein-like alpha-prolamins based on an analysis of amino acid sequences: implications for their evolution and three-dimensional structure

alpha-Prolamins are the major seed storage proteins of species of the grass tribe Andropogonea. They are unusually rich in glutamine, proline, alanine, and leucine residues and their sequences show a series of tandem repeats presumed to be the result of multiple intragenic duplication. Two new sequences of alpha-prolamin clones from Coix (pBCX25.12 and pBCX25.10) are compared with similar clones from maize and Sorghum in order to investigate evolutionary relationships between the repeat motifs and to propose a schematic model for their three-dimensional structure based on hydrophobic membrane-helix propensities and helical "wheels." A scheme is proposed for the most recent events in the evolution of the central part of the molecule (repeats 3 to 8) which involves two partial intragenic duplications and in which contemporary odd-numbered and even-numbered repeats arise from common ancestors, respectively. Each pair of repeats is proposed to form an antiparallel alpha-helical hairpin and that the helices of the molecule as a whole are arranged on a hexagonal net. The majority of helices show six faces of alternating hydrophobic and polar residues, which give rise to intersticial holes around each helix which alternate in chemical character. The model is consistent with proteins which contain different numbers of repeats, with oligomerization and with the dense packaging of alpha-prolamins within the protein body of the seed endosperm.

Beta-breakers: an aperiodic secondary structure

We have studied the architecture of parallel beta-sheets in proteins and focused on the residues that initiate and terminate the beta-strands. These beta-breaker residues are at the origin of the kink between the beta-strand and the turn that precedes or follows it. beta-Breakers can be located automatically using a consensus approach based on algorithmic secondary structure assignment, solvent accessibility and backbone dihedral angles. These beta-breakers are conformationally homogeneous with respect to side-chain solvent accessibility and backbone dihedral angle profile. A sequence-structure correlation is noted: a restricted subset of amino acids is observed at these positions. Analysis of homologous protein sequences shows that these residues are more highly conserved than other residues in the loop. We conclude that beta-breakers are the structural analogs of the N and C-terminal caps of alpha-helices. The identification of this aperiodic substructure suggests a strategy for improving secondary structure prediction and may guide site-directed mutagenesis experiments.

Influence of proline residues on protein conformation

To study the influence of proline residues on three-dimensional structure, an analysis has been made of all proline residues and their local conformations extracted from the Brookhaven Protein Data bank. We have considered the conformation of the proline itself, the relative occurrence of cis and trans peptides preceding proline residues, the influence of proline on the conformation of the preceding residue and the conformations of various proline patterns (Pro-Pro, Pro-X-Pro, etc.). The results highlight the unique role of proline in determining local conformation.

Nonenzymatic deamidation of asparaginyl and glutaminyl residues in proteins

Some asparagine and glutamine residues in proteins undergo deamidation to aspartate and glutamate with rates that depend upon the sequence and higher-order structure of the protein. Functional groups within the protein can catalyze this reaction, acting as general acids, bases, or stabilizers of the transition state. Information from specific proteins that deamidate and analysis of protein sequence and structure data bases suggest that asparagine and glutamine lability has been a selective pressure in the evolution of protein sequence and folding. Asparagine and glutamine deamidation can affect protein structure and function in natural and engineered mutant sequences, and may play a role in the regulation of protein folding, protein breakdown, and aging.

Tertiary structural constraints on protein evolutionary diversity: templates, key residues and structure prediction

The pattern of residue substitution in divergently evolving families of globular proteins is highly variable. At each position in a fold there are constraints on the identities of amino acids from both the three-dimensional structure and the function of the protein. To characterize and quantify the structural constraints, we have made a comparative analysis of families of homologous globular proteins. Residues are classified according to amino acid type, secondary structure, accessibility of the sidechain, and existence of hydrogen bonds from sidechain to other sidechains or peptide carbonyl or amide functions. There are distinct patterns of substitution especially where residues are both solvent inaccessible and hydrogen bonded through their sidechains. The patterns of residue substitution can be used to construct templates or to identify 'key' residues if one or more structures are known. Conversely, analysis of conversation and substitution across a large family of aligned sequences in terms of substitution profiles can allow prediction of tertiary environment or indicate a functional role. Similar analyses can be used to test the validity of putative structures if several homologous sequences are available.

On the appearance and role of a spacer group in the protein amino acids

Of the 20 protein amino acids, 16 have a methylene group at the beta position, and a further three bear a methine group. No aromatic, carboxamido, carboxylic carbon, or hetero atoms are attached directly to the alpha carbon, but they are separated by this methylene or occasionally by a longer n-alkylene spacer group. Therefore, the structure of the protein amino acids should rather be formulated as H2N-CH((CH2)n-R')-COOH instead of the generally accepted H2N-CH(R)-COOH. The appearance of and the role played by the spacer group are discussed in an evolutionary context. It is suggested that the spacer group appeared as a result of prebiotic selection, based on the relative abundance, racemization rate, and suitability for thermal polymerization of the protein amino acids and their homologs with various spacer group lengths. At the biotic level of evolution the requirements for ribosomal polymerization, as well as the abilities of polypeptides to maintain a stable and flexible three-dimensional structure and to bind ligands are considered and are proposed to have been responsible for the possible exclusion of longer spacer groups. It is concluded that the general role of the spacer group is to ensure the uniformity of the constant regions H2N-CH(-)-COOH and the individuality of the R' contact groups by spatially separating them.

Analysis of sequence-similar pentapeptides in unrelated protein tertiary structures. Strategies for protein folding and a guide for site-directed mutagenesis

From the most recent Brookhaven Protein Co-ordinate Databank, 229 sequence-identical pentapeptide pairs, each found in two unrelated protein structures, were collected; 9115 such pairs differing in only one residue were also gathered. For both samples the main-chain fold was conserved about 20% of the time, despite the different atomic environments presented by the unrelated protein architectures. An analysis of the substituted residues as well as the composition of the sequence-similar pentapeptides allowed several suggestions regarding protein folding mechanisms. An examination of the most frequently observed residue substitutions and their correlation with structural changes in the oligopeptide pairs yields a possible guide for site-directed mutagenesis experiments, especially when no tertiary structural information is at hand.