Induction of changes in the secondary structure of globular proteins by a hydrophobic surface

Circular dichroism, ellipsometry and radiolabeling techniques were employed to study the induction of changes in the secondary structure of BSA, myoglobin and cytochrome C by a hydrophobic surface. The results showed that adsorbed protein molecules lose their ordered native structure in the initial stage of adsorption and the structure appears to be a random or disordered conformation. Protein molecules adsorbed in later stages adopt a more ordered secondary structure ( helix and structure). The changes of secondary structure of globular proteins induced by a hydrophobic surface can be explained by the steric interaction between adsorbed proteins as well as by hydrophobic interactions during the adsorption process. In addition, there is obviously an intermediate stage in which the protein molecules are mainly in the structure, indicating that for certain proteins, the structure may be a more stable secondary structure than helix on the hydrophobic surface. Correspondence to: S.-F. Sui     

Fuzzy cluster analysis of simple physicochemical properties of amino acids for recognizing secondary structure in proteins.

Fuzzy cluster analysis has been applied to the 20 amino acids by using 65 physicochemical properties as a basis for classification. The clustering products, the fuzzy sets (i.e., classical sets with associated membership functions), have provided a new measure of amino acid similarities for use in protein folding studies. This work demonstrates that fuzzy sets of simple molecular attributes, when assigned to amino acid residues in a protein''s sequence, can predict the secondary structure of the sequence with reasonable accuracy. An approach is presented for discriminating standard folding states, using near-optimum information splitting in half-overlapping segments of the sequence of assigned membership functions. The method is applied to a nonredundant set of 252 proteins and yields approximately 73% matching for correctly predicted and correctly rejected residues with approximately 60% overall success rate for the correctly recognized ones in three folding states: alpha-helix, beta-strand, and coil. The most useful attributes for discriminating these states appear to be related to size, polarity, and thermodynamic factors. Van der Waals volume, apparent average thickness of surrounding molecular free volume, and a measure of dimensionless surface electron density can explain approximately 95% of prediction results. hydrogen bonding and hydrophobicity induces do not yet enable clear clustering and prediction.     

Algorithmic computation of knot polynomials of secondary structure elements of proteins.

The classification of protein structures is an important and still outstanding problem. The purpose of this paper is threefold. First, we utilize a relation between the Tutte and homfly polynomial to show that the Alexander-Conway polynomial can be algorithmically computed for a given planar graph. Second, as special cases of planar graphs, we use polymer graphs of protein structures. More precisely, we use three building blocks of the three-dimensional protein structure--alpha-helix, antiparallel beta-sheet, and parallel beta-sheet--and calculate, for their corresponding polymer graphs, the Tutte polynomials analytically by providing recurrence equations for all three secondary structure elements. Third, we present numerical results comparing the results from our analytical calculations with the numerical results of our algorithm-not only to test consistency, but also to demonstrate that all assigned polynomials are unique labels of the secondary structure elements. This paves the way for an automatic classification of protein structures.     

Stabilization of secondary structure elements by specific combinations of hydrophilic and hydrophobic amino acid residues is more important for proteins encoded by GC-poor genes

Stabilization of secondary structure elements by specific combinations of hydrophobic and hydrophilic amino acids has been studied by the way of analysis of pentapeptide fragments from twelve partial bacterial proteomes. PDB files describing structures of proteins from species with extremely high and low genomic GC-content, as well as with average G + C were included in the study. Amino acid residues in 78,009 pentapeptides from alpha helices, beta strands and coil regions were classified into hydrophobic and hydrophilic ones. The common propensity scale for 32 possible combinations of hydrophobic and hydrophilic amino acid residues in pentapeptide has been created: specific pentapeptides for helix, sheet and coil were described. The usage of pentapeptides preferably forming alpha helices is decreasing in alpha helices of partial bacterial proteomes with the increase of the average genomic GC-content in first and second codon positions. The usage of pentapeptides preferably forming beta strands is increasing in coil regions and in helices of partial bacterial proteomes with the growth of the average genomic GC-content in first and second codon positions. Due to these circumstances the probability of coil-sheet and helix-sheet transitions should be increased in proteins encoded by GC-rich genes making them prone to form amyloid in certain conditions. Possible causes of the described fact that importance of alpha helix and coil stabilization by specific combinations of hydrophobic and hydrophilic amino acids is growing with the decrease of genomic GC-content have been discussed.     

Protein structure comparison by probability-based matching of secondary structure elements

MOTIVATION: Protein structure comparison (PSC) has been used widely in studies of structural and functional genomics. However, PSC is computationally expensive and as a result almost all of the PSC methods currently in use look only for the optimal alignment and ignore many alternative alignments that are statistically significant and that may provide insight into protein evolution or folding. RESULTS: We have developed a new PSC method with efficiency to detect potentially viable alternative alignments in all-against-all database comparisons. The efficiency of the new PSC method derives from the ability to directly home in on a limited number of viable and ranked alignment solutions based on intuitively derived SSE (secondary structure element)-matching probabilities.     

An automatic method involving cluster analysis of secondary structures for the identification of domains in proteins.

With a growing number of structures available in the Brookhaven Protein Data Bank, automatic methods for domain identification are required for the construction of databases. Domains are considered to be clusters of secondary structure elements. Thus, helices and strands are first clustered using intersecondary structural distances between C alpha positions, and dendrograms based on this distance measure are used to identify domains. Individual domains are recognized by a disjoint factor, which enables the automatic identification and classification into disjoint, interacting, and conjoint domains. Application to a database of 83 protein families and 18 unique structures shows that the approach provides an effective delineation of boundaries and identifies those proteins that can be considered as a single domain. A quantitative estimate of the interaction between domains has been proposed. The database of protein domains is a useful tool for understanding protein folding, for recognizing protein folds, and for understanding structure-activity relationships.     

The hydrophobic cores of proteins predicted by wavelet analysis

MOTIVATION: In the process of protein construction, buried hydrophobic residues tend to assemble in a core of a protein. Methods used to predict these cores involve use or no use of sequential alignment. In the case of a close homology, prediction was more accurate if sequential alignment was used. If the homology was weak, predictions would be unreliable. A hydrophobicity plot involving the hydropathy index is useful for purposes of prediction, and smoothing is essential. However, the proposed methods are insufficient. We attempted to predict hydrophobic cores with a low frequency extracted from the hydrophobicity plot, using wavelet analysis. RESULTS: The cores were predicted at a rate of 68.7%, by cross-validation. Using wavelet analysis, the cores of non-homologous proteins can be predicted with close to 70% accuracy, without sequential alignment. AVAILABILITY: The program used in this study is available from Intergalactic Reality (http://www.intergalact.com). CONTACT: hirakawa@grt.kyushu-u.ac.jp, kuhara@grt.kyushu-u.ac.jp     

Conformational analysis of peptides corresponding to all the secondary structure elements of protein L B1 domain: secondary structure propensities are not conserved in proteins with the same fold.

The solution conformation of three peptides corresponding to the two beta-hairpins and the alpha-helix of the protein L B1 domain have been analyzed by circular dichroism (CD) and nuclear magnetic resonance spectroscopy (NMR). In aqueous solution, the three peptides show low populations of native and non-native locally folded structures, but no well-defined hairpin or helix structures are formed. In 30% aqueous trifluoroethanol (TFE), the peptide corresponding to the alpha-helix adopts a high populated helical conformation three residues longer than in the protein. The hairpin peptides aggregate in TFE, and no significant conformational change occurs in the NMR observable fraction of molecules. These results indicate that the helical peptide has a significant intrinsic tendency to adopt its native structure and that the hairpin sequences seem to be selected as non-helical. This suggests that these sequences favor the structure finally attained in the protein, but the contribution of the local interactions alone is not enough to drive the formation of a detectable population of native secondary structures. This pattern of secondary structure tendencies is different to those observed in two structurally related proteins: ubiquitin and the protein G B1 domain. The only common feature is a certain propensity of the helical segments to form the native structure. These results indicate that for a protein to fold, there is no need for large native-like secondary structure propensities, although a minimum tendency to avoid non-native structures and to favor native ones could be required.     

The graduation of secondary structure elements

A method of graduating (i.e., least-squares fitting) a smooth polynomial curve through long elements of protein secondary structure is described. It uses the Chebyshev polynomials of a discrete (integer) variable with several restraints to prevent artifactual curvatures. A new recursion formula is given which allows the evaluation of the polynomials on rational-number points as well as on the integer points. High-order splines suitable for interpolation between integer points are also discussed. The new method finds applications in graphics and in structural analysis.     

Refinement by shifting secondary structure elements improves sequence alignments

Constructing a model of a query protein based on its alignment to a homolog with experimentally determined spatial structure (the template) is still the most reliable approach to structure prediction. Alignment errors are the main bottleneck for homology modeling when the query is distantly related to the template. Alignment methods often misalign secondary structural elements by a few residues. Therefore, better alignment solutions can be found within a limited set of local shifts of secondary structures. We present a refinement method to improve pairwise sequence alignments by evaluating alignment variants generated by local shifts of template‐defined secondary structures. Our method SFESA is based on a novel scoring function that combines the profile‐based sequence score and the structure score derived from residue contacts in a template. Such a combined score frequently selects a better alignment variant among a set of candidate alignments generated by local shifts and leads to overall increase in alignment accuracy. Evaluation of several benchmarks shows that our refinement method significantly improves alignments made by automatic methods such as PROMALS, HHpred and CNFpred. The web server is available at http://prodata.swmed.edu/sfesa . Proteins 2015; 83:411–427. © 2014 Wiley Periodicals, Inc.     

Secondary structure in ribonuclease. II. Relations between folding kinetics and secondary structure elements

The kinetics of regain of 2′-CMP binding are monitored during renaturation of RNAase S. Experiments were performed by mixing equimolar amounts of S-peptide with S-protein. The S-protein fragment was incubated initially (i.e. before mixing with S-peptide) at pH 6.2 or 1.7 and various guanidine hydrochloride (GuHCl) concentrations. Three well-resolved phases are observed. The fastest phase is second-order. The reciprocal half-time increases linearly with fragment concentration and is independent of initial conditions for the S-protein fragment. An apparent on rate of k_on = 2 × 10⁵m^?1s^?1 is measured in 0.5 m-GuHCl (pH 6.2) and 20 ° C. Identical association kinetics are observed by changes in tyrosine absorbance. The fraction of native RNAase S formed in this second-order reaction strictly equals the fraction of S-protein molecules with intact β-sheet in initial conditions. The relation holds for different pH values, GuHCl concentrations and temperatures. The fraction of apparent helical content of S-protein in initial conditions may also vary but this is not reflected by the association reaction. We interpret this to mean that the β-sheet but not the α-helices must be preformed in initial conditions in order to generate the high-affinity peptide binding site of S-protein. Furthermore, it is concluded that the S-protein moiety β-sheet forms or unfolds in a single one-step reaction. 2′-CMP binding reports, additionally, two slower phases of renaturation. These are produced by S-protein molecules that have their β-sheet unfolded in initial conditions. It is observed that a unique dependence of these two folding rates exists for RNAase A, RNAase S and S-protein as function of t_m, the temperature of half-completion of thermal denaturation as measured by unfolding of the β-sheet in the respective compound in final conditions. The t_m value varies with changing pH, with GuHCl concentration and (for RNAase S) with changing fragment concentration. The findings are interpreted to argue in favor of a sequential mechanism of folding, where the stability of a structural precursor determines the rate of folding.     

Role of loops connecting secondary structure elements in the stabilization of proteins isolated from thermophilic organisms

It has been recently discovered that the connection of secondary structure elements (ββ‐unit, βα‐ and αβ‐units) in proteins follows quite stringent principles regarding the chirality and the orientation of the structural units (Koga et al., Nature 2012;491:222–227). By exploiting these rules, a number of protein scaffolds endowed with a remarkable thermal stability have been designed (Koga et al., Nature 2012;491:222–227). By using structural databases of proteins isolated from either mesophilic or thermophilic organisms, we here investigate the influence of supersecondary associations on the thermal stability of natural proteins. Our results suggest that β‐hairpins of proteins from thermophilic organisms are very frequently characterized by shortenings of the loops. Interestingly, this shortening leads to states that display a very strong preference for the most common connectivity of the strands observed in native protein hairpins. The abundance of selective states in these proteins suggests that they may achieve a high stability by adopting a strategy aimed to reduce the possible conformations of the unfolded ensemble. In this scenario, our data indicate that the shortening is effective if it increases the adherence to these rules. We also show that this mechanism may operate in the stabilization of well‐known protein folds (thioredoxin and RNase A). These findings suggest that future investigations aimed at defining mechanism of protein stabilization should also consider these effects.     

Structural homology among the peroxidase enzyme family revealed by hydrophobic cluster analysis

A number of peroxidase amino acid sequences show limited homology to short regions comprising the known active site cleft of yeast cytochrome c peroxidase. Otherwise no clear homology is visible in linear alignments between this enzyme and other peroxidases. We have subjected eight peroxidase sequences to hydrophobic cluster analysis. Our results suggest that these peroxidases are evolutionary related and that they share many folding characteristics.     

Phylogenetic analysis of tmRNA secondary structure.

The bacterial tmRNA acts with dual tRNA-like and mRNA-like character to tag incomplete translation products for degradation. Comparative analysis of 17 tmRNA genes (including eight new sequences) has allowed us to deduce conserved features of the tmRNA secondary structure. Except in a segment that includes the first codon of the tag reading frame, tmRNA is highly structured, with four pseudoknots and a total of 11 conserved base pairing regions. The previously identified tRNA minihelix structure is connected by a long base paired region to a large structured domain composed of a pseudoknot, followed by the tag reading frame and a string of three rather similar pseudoknots. The conservation of numerous structural elements among diverse eubacterial species indicates that these elements have important function beyond simply forming an endonuclease-resistant link between the reading frame and the tRNA-like domain.     

Contribution of proline to the pre-structuring tendency of transient helical secondary structure elements in intrinsically disordered proteins

Background

IDPs function without relying on three-dimensional structures. No clear rationale for such a behavior is available yet. PreSMos are transient secondary structures observed in the target-free IDPs and serve as the target-binding “active” motifs in IDPs. Prolines are frequently found in the flanking regions of PreSMos. Contribution of prolines to the conformational stability of the helical PreSMos in IDPs is investigated.

Methods

MD simulations are performed for several IDP segments containing a helical PreSMo and the flanking prolines. To measure the influence of flanking-prolines on the structural content of a helical PreSMo calculations were done for wild type as well as for mutant segments with Pro → Asp, His, Lys, or Ala. The change in the helicity due to removal of a proline was measured both for the PreSMo region and for the flanking regions.

Results

The α-helical content in ~ 70% of the helical PreSMos at the early stage of simulation decreases due to replacement of an N-terminal flanking proline by other residues whereas the helix content in nearly all PreSMos increases when the same replacements occur at the C-terminal flanking region. The helix destabilizing/terminating role of the C-terminal flanking prolines is more pronounced than the helix promoting effect of the N-terminal flanking prolines.

General significance

This work represents a novel example demonstrating that a proline is encoded in an IDP with a defined purpose. The helical PreSMos presage their target-bound conformations. As they most likely mediate IDP-target binding via conformational selection their helical content can be an important feature for IDP function.     

Identification of secondary structure elements in intermediate-resolution density maps

An increasing number of structural studies of large macromolecular complexes, both in X-ray crystallography and cryo-electron microscopy, have resulted in intermediate-resolution (5-10 A) density maps. Despite being limited in resolution, significant structural and functional information may be extractable from these maps. To aid in the analysis and annotation of these complexes, we have developed SSEhunter, a tool for the quantitative detection of alpha helices and beta sheets. Based on density skeletonization, local geometry calculations, and a template-based search, SSEhunter has been tested and validated on a variety of simulated and authentic subnanometer-resolution density maps. The result is a robust, user-friendly approach that allows users to quickly visualize, assess, and annotate intermediate-resolution density maps. Beyond secondary structure element identification, the skeletonization algorithm in SSEhunter provides secondary structure topology, which is potentially useful in leading to structural models of individual molecular components directly from the density.     

Automatic identification of secondary structure in globular proteins

A computer program is used to analyse automatically and objectively the atomic co-ordinates of a large number of globular proteins in order to identify the regions of α-helix, β-sheet and reverse-turn secondary structure. Several different criteria for the assignment of secondary structure are tested for accuracy, reproducibility and efficiency. The most successful criterion, which is based on patterns of peptide hydrogen bonds, inter-C^α distances and inter-C^α torsion angles, is used to find the secondary structure of all the proteins studied. The accuracy of the derived assignments is assessed by comparing them with the secondary structure reported in the literature for each protein. The reliability of the methods is assessed by comparing the secondary structures derived from the independently determined sets of co-ordinates available for some proteins.We provide the first objective and consistent compilation of α-helix, β-sheet and reverse-turn secondary structure in almost all globular proteins of known tertiary structure. These data will be invaluable for analysing the relative tendencies of different amino acids to occur in different types of secondary structure, for analysing the regularity of the secondary structure itself, and for analysing how the pieces of secondary structure fit together to form the globular tertiary structure of each protein.     

Metals in proteins: cluster analysis studies

We have conducted a prospective analysis of the Protein Data Bank in order to study certain constituents of proteins: elements that are neither halogens nor phosphorus nor part of the biological amino acid set. A sample of 5749 structures was analyzed and classified according to the 56 elements encountered. Fifteen metals (Na, Mg, K, Ca, Mn, Fe, Co, Ni, Cu, Zn, As, Mo, Cd, W, Hg) are involved in almost half of the structures, with each metal figuring in more than 100 structures. We analyzed this subsample in more detail by computing the amino acid residues occurring within a coordination sphere of 5 Å centered on the element, and using methods of cluster analysis to group the elements. The analyses undertaken here are able to distinguish between real components of proteins and elements inserted by artefacts of the crystallization process or experimental techniques.