Determination of protein secondary structure using factor analysis of infrared spectra

A method is presented for determining the secondary structural composition of a protein in aqueous solution from its infrared spectrum. A factor analysis approach is used to analyze the infrared spectra of 18 proteins whose crystal structures are known from X-ray studies. Factor analysis followed by multiple linear regression identifies those eigenspectra that correlate with the variation in properties described by the calibration set. The properties of interest in this study are % alpha-helix, % beta-sheet, and % turns. In the analysis of an unknown, the factor loadings required to reproduce its spectrum are substituted in the regression equation for each property to predict its secondary structural composition. The accuracy of the method was determined by removing each standard, in turn, from the calibration set and using a calibration set generated from the remainder to predict its composition. By this method we obtain standard errors of prediction of 3.9% for alpha-helix, 8.3% for beta-sheet, and 6.6% for turns. The method may also be applied to the spectra of proteins in 2H2O. The method has important advantages over those currently in use for the quantitative analysis of the infrared spectra of proteins. Manipulation of the spectrum is kept to a minimum, no curve-fitting is necessary, and the several amide I band components need not be assigned.     

Sequence swapping does not result in conformation swapping for the beta4/beta5 and beta8/beta9 beta-hairpin turns in human acidic fibroblast growth factor

The beta-turn is the most common type of nonrepetitive structure in globular proteins, comprising ~25% of all residues; however, a detailed understanding of effects of specific residues upon beta-turn stability and conformation is lacking. Human acidic fibroblast growth factor (FGF-1) is a member of the beta-trefoil superfold and contains a total of five beta-hairpin structures (antiparallel beta-sheets connected by a reverse turn). beta-Turns related by the characteristic threefold structural symmetry of this superfold exhibit different primary structures, and in some cases, different secondary structures. As such, they represent a useful system with which to study the role that turn sequences play in determining structure, stability, and folding of the protein. Two turns related by the threefold structural symmetry, the beta4/beta5 and beta8/beta9 turns, were subjected to both sequence-swapping and poly-glycine substitution mutations, and the effects upon stability, folding, and structure were investigated. In the wild-type protein these turns are of identical length, but exhibit different conformations. These conformations were observed to be retained during sequence-swapping and glycine substitution mutagenesis. The results indicate that the beta-turn structure at these positions is not determined by the turn sequence. Structural analysis suggests that residues flanking the turn are a primary structural determinant of the conformation within the turn.     

A helix-turn-strand structural motif common in alpha-beta proteins.

By exhaustive structural comparisons, we have found that about one-third of the alpha-helix-turn-beta-strand polypeptides in alpha-beta barrel domains share a common structural motif. The chief characteristics of this motif are that first, the geometry of the turn between the alpha-helix and the beta-strand is somewhat constrained, and second, the beta-strand contains a hydrophobic patch that fits into a hydrophobic pocket on the alpha-helix. The geometry of the turn does not seem to be a major determinant of the alpha-beta unit, because the turns vary in length from four to six residues. However, the motif does not occur when there are few constraints on the geometry of the turn-for instance, when the turns between the alpha-helix and the beta-strands are very long. It also occurs much less frequently in flat-sheet alpha-beta proteins, where the topology is much less regular and the amount of twist on the sheet varies considerably more than in the barrel proteins. The motif may be one of the basic building blocks from which alpha-beta barrels are constructed.     

Structural and sequence features of two residue turns in beta‐hairpins

Beta‐turns in beta‐hairpins have been implicated as important sites in protein folding. In particular, two residue β‐turns, the most abundant connecting elements in beta‐hairpins, have been a major target for engineering protein stability and folding. In this study, we attempted to investigate and update the structural and sequence properties of two residue turns in beta‐hairpins with a large data set. For this, 3977 beta‐turns were extracted from 2394 nonhomologous protein chains and analyzed. First, the distribution, dihedral angles and twists of two residue turn types were determined, and compared with previous data. The trend of turn type occurrence and most structural features of the turn types were similar to previous results, but for the first time Type II turns in beta‐hairpins were identified. Second, sequence motifs for the turn types were devised based on amino acid positional potentials of two‐residue turns, and their distributions were examined. From this study, we could identify code‐like sequence motifs for the two residue beta‐turn types. Finally, structural and sequence properties of beta‐strands in the beta‐hairpins were analyzed, which revealed that the beta‐strands showed no specific sequence and structural patterns for turn types. The analytical results in this study are expected to be a reference in the engineering or design of beta‐hairpin turn structures and sequences. Proteins 2014; 82:1721–1733. © 2014 Wiley Periodicals, Inc.     

Accurate computer-based design of a new backbone conformation in the second turn of protein L.

The rational design of loops and turns is a key step towards creating proteins with new functions. We used a computational design procedure to create new backbone conformations in the second turn of protein L. The Protein Data Bank was searched for alternative turn conformations, and sequences optimal for these turns in the context of protein L were identified using a Monte Carlo search procedure and an energy function that favors close packing. Two variants containing 12 and 14 mutations were found to be as stable as wild-type protein L. The crystal structure of one of the variants has been solved at a resolution of 1.9 A, and the backbone conformation in the second turn is remarkably close to that of the in silico model (1.1 A RMSD) while it differs significantly from that of wild-type protein L (the turn residues are displaced by an average of 7.2 A). The folding rates of the redesigned proteins are greater than that of the wild-type protein and in contrast to wild-type protein L the second beta-turn appears to be formed at the rate limiting step in folding.     

Analysis of gammabeta, betagamma, gammagamma, betabeta continuous turns in proteins.

We report the observation of continuous turns in proteins which comprise individual gamma-turns or beta-turns or both that are situated immediately one after the other along the polypeptide chain. The continuous turns were identified from a representative data set of three-dimensional protein crystal structures. The gammabeta/betagamma, gammagamma and betabeta continuous turns represent peptides of varying amino acid residue lengths and conformations. The continuous turns frequently observed in proteins were: gammabeta, between a coil and a strand; betagamma, between a helix and a strand; gammagamma, between coils; and betabeta, either between a strand and a coil or between strands or coils. We determined the statistically significant amino acid residue preferences at individual positions in the turn, calculated amino acid positional potentials and analyzed main chain hydrogen bonds and side-chain interactions likely to stabilize the continuous turns. The data on continuous turns have been integrated in the database of structural motifs in proteins (DSMP) on our web server at (http://www.cdfd.org.in/dsmp.html). This is useful to make queries on sequences compatible with different continuous turns.     

Sequence and structural parameters enhancing adaptation of proteins to low temperatures

A systematic analysis compared sequence and structural parameters distributions between 13 pairs of psychrophilic and mesophilic proteins for elucidating the cold adaptation parameters. The results of statistical test (t-test) revealed that helical content, tight turn content, disulfide bonds and hydrogen bonds do not show significant difference between psychrophilic and mesophilic proteins. However, it was demonstrated in this study that a larger proportion of open beta-turn in psychrophilic proteins is an effective parameter in specific activity at low temperature. In addition, substitution of amino acids of charged and aliphatic groups with amino acids of tiny and small groups in protein chains, tight turns and alpha-helices in the direction from mesophilic to psychrophilic proteins is one of the mechanisms of low temperature adaptation. Such sequence and structural parameter differences would help to develop a strategy for designing cold-adapted proteins.     

Prediction of tight turns and their types in proteins

A tight turn in protein structure is defined as a site where (i) a polypeptide chain reverses its overall direction, i.e., leads the chain to fold back on itself by nearly 180 degrees, and (ii) the amino acid residues directly involved in forming the turn are no more than six. Tight turns are generally categorized as delta-turn, gamma-turn, beta-turn, alpha-turn, and pi-turn, which are formed by two-, three-, four-, five-, and six-amino-acid residues, respectively. According to the folding mode, each of such tight turns can be further classified into several different types. Tight turns play an important role in globular proteins from both the structural and functional points of view. In view of this, various efforts have been made to predict tight turns and their types. This Review summarizes the development in this area, with an emphasis focused on the most recent work concerned that is featured by the sequence-coupled model. Meanwhile, the future challenge in this area has also been briefly addressed.     

Prediction and classification of alpha-turn types

Tight turns play an important role in globular proteins from both the structural and functional points of view. Of tight turns, beta-turns and gamma-turns have been extensively studied, but alpha-turns were little investigated. Recently, a systematic search for alpha-turns was conducted by V. Pavone et al. [(1996) Biopolymers, Vol. 38, pp. 705-721] from 190 proteins (221 protein chains). They found 356 alpha-turns that were classified into nine different types according to their backbone trajectory features. In view of this new discovery, a sequence-coupled model based on Markov chain theory is proposed for predicting the alpha-turn types in proteins. The high rates of correct prediction by resubstitution test and jackknife test imply that that the formation of different alpha-turn types is evidently correlated with the sequence of a pentapeptide, and hence can be approximately predicted based on the sequence information of the pentapeptide alone, although the role of its interaction with the other part of a protein cannot be completely ignored. The algorithm presented here can also be used to conduct the prediction in which a distinction between alpha-turns and non-alpha-turns is also required.     

Continuous and discontinuous domains: an algorithm for the automatic generation of reliable protein domain definitions.

An algorithm is presented for the fast and accurate definition of protein structural domains from coordinate data without prior knowledge of the number or type of domains. The algorithm explicitly locates domains that comprise one or two continuous segments of protein chain. Domains that include more than two segments are also located. The algorithm was applied to a nonredundant database of 230 protein structures and the results compared to domain definitions obtained from the literature, or by inspection of the coordinates on molecular graphics. For 70% of the proteins, the derived domains agree with the reference definitions, 18% show minor differences and only 12% (28 proteins) show very different definitions. Three screens were applied to identify the derived domains least likely to agree with the subjective definition set. These screens revealed a set of 173 proteins, 97% of which agree well with the subjective definitions. The algorithm represents a practical domain identification tool that can be run routinely on the entire structural database. Adjustment of parameters also allows smaller compact units to be identified in proteins.     

Formation of cytoplasmic turns between two closely spaced transmembrane helices during membrane protein integration into the ER membrane

The helical hairpin, two closely spaced transmembrane helices separated by a short turn, is a recurring structural element in integral membrane proteins, and may serve as a compact unit that inserts into the membrane en bloc. Previously, we have determined the propensities of the 20 natural amino acids, when present in the middle of a long hydrophobic stretch, to induce the formation of a helical hairpin with a lumenally exposed turn during membrane protein assembly into the endoplasmic reticulum membrane. Here, we present results from a similar set of measurements, but with the turn placed on the cytoplasmic side of the membrane. We find that a significantly higher number of turn-promoting residues need to be present to induce a cytoplasmic turn compared to a lumenal turn, and that, in contrast to the lumenal turn, the positively charged residues Arg and Lys are the strongest turn-promoters in cytoplasmic turns. These results suggest that the process of turn formation between transmembrane helices is different for lumenal and cytoplasmic turns.     

Polypeptide secondary structure determination by nuclear magnetic resonance observation of short proton-proton distances

The use of proton-proton nuclear Overhauser enhancement (NOE) distance information for identification of polypeptide secondary structures in non-crystalline proteins was investigated by stereochemical studies of standard secondary structures and by statistical analyses of the secondary structures in the crystal conformations of a group of globular proteins. Both regular helix and beta-sheet secondary structures were found to contain a dense network of short 1H-1H distances. The results obtained imply that the combined information on all these distances obtained from visual inspection of the two-dimensional NOE (NOESY) spectra is sufficient for determination of the helical and beta-sheet secondary structures in small globular proteins. Furthermore, cis peptide bonds can be identified from unique, short sequential proton-proton distances. Limitations of this empirical approach are that the exact start or end of a helix may be difficult to define when the adjoining residues form a tight turn, and that unambiguous identification of tight turns can usually be obtained only in the hairpins of antiparallel beta-structures. The short distances between protons in pentapeptide segments of the different secondary structures have been tabulated to provide a generally applicable guide for the analysis of NOESY spectra of proteins.     

FAMCS： Finding All Maximal Common Substructures in Proteins

Finding the common substructures shared by two proteins is considered as one of the central issues in computational biology because of its usefulness in understanding the structure-function relationship and application in drug and vaccine design. In this paper, we propose a novel algorithm called FAMCS （Finding All Maximal Common Substructures） for the common substructure identification problem. Our method works initially at the protein secondary structural element （SSE） level and starts with the identification of all structurally similar SSE pairs. These SSE pairs are then merged into sets using a modified Apriori algorithm, which will test the similarity of various sets of SSE pairs incrementally until all the maximal sets of SSE pairs that deemed to be similar are found. The maximal common substructures of the two proteins will be formed from these maximal sets. A refinement algorithm is also proposed to fine tune the alignment from the SSE level to the residue level. Comparison of FAMCS with other methods on various proteins shows that FAMCS can address all four requirements and infer interesting biological discoveries.     

Analysis of gammabeta, betagamma, gammagamma, betabeta multiple turns in proteins.

The number of gamma-turns in a representative protein dataset selected from the current Protein Data Bank has increased almost seven times during the past decade. Eighty percent classic gamma-turns and 57% inverse gamma-turns are associated as multiple turns with either another y-turn or a beta-turn. We refer to these as multiple turns of the (gammabeta)1,2,3 or (betagamma)1,2,3 type, depending upon whether the gamma-turn is before or after the beta-turn along the protein chain, respectively. However, for multiple turns involving only gamma-turns, we follow the nomenclature analogous to that proposed earlier for the multiple (or double) beta-turns. Fifty-eight per cent beta-turns are associated as multiple turns with another beta-turn. We extracted multiple turns from the protein dataset and classified them on the basis of individual gamma- or beta-turn types and the number of overlapping residues. Furthermore, we evaluated the amino acid positional potentials and determined the statistically significant amino acid preferences, hydrogen bond/side-chain interaction preferences in the multiple turns and secondary structure preferences for residues immediately flanking these turns. The results of our analysis would be useful in the modeling, prediction or design of multiple turns in proteins. The amino acid sequence corresponding to the multiple turn, position in the protein chain, PDB Code/chain in which multiple turn is present and the individual turn types constituting the multiple turns are available from our website and this information would also be integrated in our Database of Structural Motifs in Proteins (http://www.cdfd.org.in/dsmp.html).     

Stabilization of globular proteins via introduction of temperature-activated elastin-based switches

To investigate whether swapping native turns of a globular protein with an elastin-based turn sequence (VPGVG) can increase its thermostability, we have performed molecular dynamics simulations of wild-type chymotrypsin inhibitor 2 (CI2) and variants containing elastin-based turns at 10 degrees C and 40 degrees C. Wild-type CI2 is more stable at 10 degrees C, while both of the variant forms are more stable at 40 degrees C. Detailed analyses indicate that the elastin-based turns do indeed contribute to the inverse temperature behavior of the modified proteins. Therefore, swapping a wild-type turn sequence with an elastin-based turn provides a novel way to both improve stability of target proteins at body temperature and to possibly introduce a temperature-sensitive switch.     

Role of amino acid residues at turns in the conformational stability and folding of human lysozyme

To clarify the role of amino acid residues at turns in the conformational stability and folding of a globular protein, six mutant human lysozymes deleted or substituted at turn structures were investigated by calorimetry, GuHCl denaturation experiments, and X-ray crystal analysis. The thermodynamic properties of the mutant and wild-type human lysozymes were compared and discussed on the basis of their three-dimensional structures. For the deletion mutants, Delta47-48 and Delta101, the deleted residues are in turns on the surface and are absent in human alpha-lactalbumin, which is homologous to human lysozyme in amino acid sequence and tertiary structure. The stability of both mutants would be expected to increase due to a decrease in conformational entropy in the denatured state; however, both proteins were destabilized. The destabilizations were mainly caused by the disappearance of intramolecular hydrogen bonds. Each part deleted was recovered by the turn region like the alpha-lactalbumin structure, but there were differences in the main-chain conformation of the turn between each deletion mutant and alpha-lactalbumin even if the loop length was the same. For the point mutants, R50G, Q58G, H78G, and G37Q, the main-chain conformations of these substitution residues located in turns adopt a left-handed helical region in the wild-type structure. It is thought that the left-handed non-Gly residue has unfavorable conformational energy compared to the left-handed Gly residue. Q58G was stabilized, but the others had little effect on the stability. The structural analysis revealed that the turns could rearrange the main-chain conformation to accommodate the left-handed non-Gly residues. The present results indicate that turn structures are able to change their main-chain conformations, depending upon the side-chain features of amino acid residues on the turns. Furthermore, stopped-flow GuHCl denaturation experiments on the six mutants were performed. The effects of mutations on unfolding-refolding kinetics were significantly different among the mutant proteins. The deletion/substitutions in turns located in the alpha-domain of human lysozyme affected the refolding rate, indicating the contribution of turn structures to the folding of a globular protein.     

Solution structure of the loops of bacteriorhodopsin closely resembles the crystal structure

Bacteriorhodopsin is one of very few transmembrane proteins for which high resolution structures have been solved. The structure shows a bundle of seven helices connected by six turns. Some turns in proteins are stabilized by short range interactions and can behave as small domains. These observations suggest that peptides containing the sequence of the turns in a membrane protein such as bacteriorhodopsin may form stable turn structures in solution. To test this hypothesis, we determined the solution structure of three peptides each containing the sequence of one of the turns in bacteriorhodopsin. The solution structures of the peptides closely resemble the structures of the corresponding turns in the high resolution structures of the intact protein.     

Analysis and design of turns in alpha-helical hairpins

Although the analysis and design of turns that connect the strands in antiparallel beta-hairpins has reached an advanced state, much less is known concerning turns between antiparallel helices in helical hairpins. We have conducted an analysis of the structures and sequence preferences of two types of interhelical turns, each of which connects the two helices by a two-residue linker in an alphaL-beta conformation. Based on this analysis, it became apparent that the turn introduced into a designed four-helix bundle protein, DF1, did not occur within an optimal structural context. DF1 is a dimeric model for the diiron class of proteins. A longer loop with a beta-alphaR-beta conformation was inserted between two helices in the protein, and a sequence was chosen to stabilize its conformation. X-ray crystallography and NMR analysis of the protein showed the structure to be in excellent agreement with design.     

Roles of beta-turns in protein folding: from peptide models to protein engineering

Reverse turns are a major class of protein secondary structure; they represent sites of chain reversal and thus sites where the globular character of a protein is created. It has been speculated for many years that turns may nucleate the formation of structure in protein folding, as their propensity to occur will favor the approximation of their flanking regions and their general tendency to be hydrophilic will favor their disposition at the solvent-accessible surface. Reverse turns are local features, and it is therefore not surprising that their structural properties have been extensively studied using peptide models. In this article, we review research on peptide models of turns to test the hypothesis that the propensities of turns to form in short peptides will relate to the roles of corresponding sequences in protein folding. Turns with significant stability as isolated entities should actively promote the folding of a protein, and by contrast, turn sequences that merely allow the chain to adopt conformations required for chain reversal are predicted to be passive in the folding mechanism. We discuss results of protein engineering studies of the roles of turn residues in folding mechanisms. Factors that correlate with the importance of turns in folding indeed include their intrinsic stability, as well as their topological context and their participation in hydrophobic networks within the protein's structure.     

In silico platform for predicting and initiating β‐turns in a protein at desired locations

Numerous studies have been performed for analysis and prediction of β‐turns in a protein. This study focuses on analyzing, predicting, and designing of β‐turns to understand the preference of amino acids in β‐turn formation. We analyzed around 20,000 PDB chains to understand the preference of residues or pair of residues at different positions in β‐turns. Based on the results, a propensity‐based method has been developed for predicting β‐turns with an accuracy of 82%. We introduced a new approach entitled “Turn level prediction method,” which predicts the complete β‐turn rather than focusing on the residues in a β‐turn. Finally, we developed BetaTPred3, a Random forest based method for predicting β‐turns by utilizing various features of four residues present in β‐turns. The BetaTPred3 achieved an accuracy of 79% with 0.51 MCC that is comparable or better than existing methods on BT426 dataset. Additionally, models were developed to predict β‐turn types with better performance than other methods available in the literature. In order to improve the quality of prediction of turns, we developed prediction models on a large and latest dataset of 6376 nonredundant protein chains. Based on this study, a web server has been developed for prediction of β‐turns and their types in proteins. This web server also predicts minimum number of mutations required to initiate or break a β‐turn in a protein at specified location of a protein. Proteins 2015; 83:910–921. © 2015 Wiley Periodicals, Inc.