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

A hydrophobic spine stabilizes a surface-exposed α-helix according to analysis of the solvent-accessible surface area

Background

Most of hydrophilic and hydrophobic residues are thought to be exposed and buried in proteins, respectively. In contrast to the majority of the existing studies on protein folding characteristics using protein structures, in this study, our aim was to design predictors for estimating relative solvent accessibility (RSA) of amino acid residues to discover protein folding characteristics from sequences.

Methods

The proposed 20 real-value RSA predictors were designed on the basis of the support vector regression method with a set of informative physicochemical properties (PCPs) obtained by means of an optimal feature selection algorithm. Then, molecular dynamics simulations were performed for validating the knowledge discovered by analysis of the selected PCPs.

Results

The RSA predictors had the mean absolute error of 14.11% and a correlation coefficient of 0.69, better than the existing predictors. The hydrophilic-residue predictors preferred PCPs of buried amino acid residues to PCPs of exposed ones as prediction features. A hydrophobic spine composed of exposed hydrophobic residues of an α-helix was discovered by analyzing the PCPs of RSA predictors corresponding to hydrophobic residues. For example, the results of a molecular dynamics simulation of wild-type sequences and their mutants showed that proteins 1MOF and 2WRP_H16I (Protein Data Bank IDs), which have a perfectly hydrophobic spine, have more stable structures than 1MOF_I54D and 2WRP do (which do not have a perfectly hydrophobic spine).

Conclusions

We identified informative PCPs to design high-performance RSA predictors and to analyze these PCPs for identification of novel protein folding characteristics. A hydrophobic spine in a protein can help to stabilize exposed α-helices.

Electronic supplementary material

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

Construction and characterization of protein libraries composed of secondary structure modules

Only a minute fraction of all possible protein sequences can exist in the genomes of all life forms. To explore whether physicochemical constraints or a lack of need causes the paucity of different protein folds, we set out to construct protein libraries without any restriction of topology. We generated different libraries (all alpha-helix, all beta-strand, and alpha-helix plus beta-strand) with an average length of 100 amino acid residues, composed of designed secondary structure modules (alpha-helix, beta-strand, and beta-turn) in various proportions, based primarily on the patterning of polar and nonpolar residues. We wished to explore that part of sequence space that is rich in secondary structure. The analysis of randomly chosen clones from each of the libraries showed that, despite the low sequence homology to known protein sequences, a substantial proportion of the library members containing alpha-helix modules were indeed helical, possess a defined oligomerization state, and showed cooperative chemical unfolding behavior. On the other hand, proteins composed of mainly beta-strand modules tended to form amyloid-like fibrils and were among the least soluble proteins ever reported. We found that a large fraction of members in non-beta-strand-containing protein libraries that are distant from natural proteins in sequence space possess unexpectedly favorable properties. These results reinforce the efficacy of applying binary patterning to design proteins with native-like properties despite lack of restriction in topology. Because of the intrinsic tendency of beta-strand modules to aggregate, their presence requires precise topologic arrangement to prevent fibril formation.

Variants of 3(10)-helices in proteins

An analysis of the shortest 3(10)-helices, containing three helical residues and two flanking capping residues that participate in two consecutive i + 3 --> i hydrogen bonds, shows that not all helices belong to the classic 3(10)-helix, where the three central residues adopt the right-handed helical conformation (alpha(R)). Three variants identified are: 3L10-helix with all residues in the left-handed helical region (alpha(L)), 3EL10-helix where the first residue is in the extended region followed by two residues in the alpha(L) conformation, and its mirror-image, the 3E'R10-helix. In the context of these helices, as well as the equivalent variants of alpha-helices, the length dependence of the handedness of secondary structures in protein structure is discussed. There are considerable differences in the amino acid preferences at different positions in the various types of 3(10)-helices. Each type of 3(10)-helix can be thought to be made up of an extension of a particular type of beta-turn (made up of residues i to i + 3) such that the (i + 3)th residue assumes the same conformation as the preceding residue. Distinct residue preferences at i and i + 3 positions seem to decide whether a particular stretch of four residues will be a beta-turn or a 3(10)-helix in the folded structure.

Occurrence, conformational features and amino acid propensities for the pi-helix

The most abundant helix type in proteins is the alpha-helix, accounting for about 31% of amino acid secondary structure states, while the 3(10)-helix accounts for about 4%. The pi-helix appears to be extremely rare and is considered to be unstable. Existing secondary structure definition methods find very few within the Protein Data Bank. Using an improved pi-helix definition algorithm to search a non-redundant subset of high-resolution and well-refined protein structures, we found that almost every tenth protein contained a pi-helix. This enabled us to show for the first time that the pi-helix has structural parameters that are different from the hypothesized model values. It also has distinctive amino acid preferences and it is conserved within functionally related proteins. Features that may contribute to the stability of the pi-helical structure have also been identified. In addition to hydrogen bonds, several other factors contribute to the stability of pi-helices. The pi-helix may have some functional advantages over other helical structures. Thus, we describe cases where the side chains of functionally important residues at every fourth position within a pi-helix could be aligned and brought close together in a way that would not be allowed by any other helix type.

Enhancing the T-->R transition of insulin by helix-promoting sequence modifications at the N-terminal B-chain

Structurally, the T-->R transition of insulin mainly consists of a rearrangement of the N-terminal B-chain (residues B1-B8) from extended to helical in one or both of the trimers of the hexamer. The dependence of the transition on the nature of the ligands inducing it, such as inorganic anions or phenolic compounds, as well as of the metal ions complexing the hexamer, has been the subject of extensive investigations. This study explores the effect of helix-enhancing modifications of the N-terminal B-chain sequence where the transition actually occurs, with special emphasis on N-capping. In total 15 different analogues were prepared by semisynthesis. 80% of the hexamers of the most successful analogues with zinc were found to adopt the T3R3 state in the absence of any transforming ligands, as compared to only 4% of wild-type insulin. Transformation with SCN- ions can exceed the T3R3 state where it stops in the case of wild-type insulin. Full transformation to the R6 state can be achieved by only one-tenth the phenol concentration required for wild-type insulin, i.e. almost at the stoichiometric ratio of 6 phenols per hexamer.

Novel protein structural motifs containing two-turn and longer 3(10)-helices

The 3(10)-helix constitutes a small but significant fraction of secondary structural elements in proteins. Protein data base surveys have shown these helices to be present as alpha-helical extensions, in loops and as connectors between beta-strands. The present work focuses on two-turn and longer 3(10)-helices where we establish that two-turn and longer 3(10) helices, unlike the more abundant single-turn 3(10)-helices, frequently occur independent of any other contiguous secondary structural elements. More importantly, a large fraction of these independent two-turn and longer 3(10)-helices, along with alpha-helices and beta-strands, are found to form novel super-secondary structural motifs in several proteins with possible implications for protein folding, local conformational relaxation and biological functions.

A third fibronectin type III domain in the extracellular region of the insulin receptor family

The insulin receptor family consists of the homologous tyrosine kinase receptors, insulin receptor (IR), insulin-like growth factor 1 receptor (IGF1R) and insulin receptor-related receptor. The three-dimensional structures of the tyrosine kinase domain of the IR and the first three extracellular domains (L1, Cys-rich and L2) of the IGF1R are known. Here we present evidence that the connecting domain of the IR family is a member of the fibronectin type II (FnIII) superfamily. Structure-based alignment of FnIII domains reveals several key residues that are also conserved in the sequence of the connecting domain. The alignment of the connecting domain with FnIII domains is in good agreement with secondary structure prediction. A model of the connecting domain shows a hydrophobic core formed by the conserved residues and is consistent with previously known biochemical data. This suggests that the IR family contains three FnIII domains in tandem in the extracellular juxtamembrane region.

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.

Population statistics of protein structures: lessons from structural classifications

Structural classifications aid the interpretation of proteins by describing degrees of structural and evolutionary relatedness. They have also recently revealed strikingly skewed distributions at all levels; for example, a small number of folds are far more common than others, and just a few superfamilies are known to have diverged widely. The classifications also provide an indication of the total number of superfamilies in nature.