Identification and analysis of extended strands and β-sheets in globular proteins

A method to identify β-sheets in globular proteins from extended strands, using only α-carbon positions, has been developed. The strands that form β-sheets are picked up by means of simple distance criteria. The method has been tested by applying it to three proteins with accurately known secondary structures. It has also been applied to ten other proteins wherein only α-carbon coordinates are available, and the list of β-sheets obtained. The following points are worth noting: (i) The sheets identified by the algorithm are found to agree satisfactorily with the reported ones based on backbone hydrogen bonding, wherever this information is available. (ii) β-Strands that do not form parts of any sheet are a common feature of protein structures. (iii) Such isolated β-strands tend to be short. (iv) The conformation corresponding to the preferred right-handed twist of the sheet is overwhelmingly observed in both the sheet-forming and isolated β-strands.     

Specific recognition in the tertiary structure of β-sheets of proteins

The frequency of occurrence of nearest neighbour residue pairs on adjacent antiparallel (β_A) and parallel (β_P) strands is obtained from 30 known protein structures. The specificity of interstrand recognition due to such pairing as a factor in the folding of β-sheets is studied by statistical methods. Residues of sufficiently high count for statistical analysis are treated individually while the rest are combined into small groups of similar size, polarity, and/or genetic exchangeability. The hypothesis of specific recognition between individuals and small groups is contrasted with the alternative hypothesis of non-specific recognition between broad classes (hydrophobia, neutral, polar) of residues. A χ² test of pair correlations favours specific recognition against non-specific recognition with a high level of confidence. The largest and most significant correlations are: Ser/Thr (1.9 ± 0.3), Ile/Val (1.7 ± 0.3) and Lys-Arg/Asp-Gln (1.8 ± 0.3) in β_A, and Ile/Leu (1.9 ± 0.4) in β_P. The pair Gly/Gly never occurs in any β-sheet. The specific residue-pair correlations derived here may be useful in statistical prediction methods of protein tertiary structure.     

Analysis and prediction of the packing of α-helices against a β-sheet in the tertiary structure of globular proteins

The packing of α-helices and β-sheets in six $\text{α}\text{β}$ proteins (e.g. flavodoxin) has been analysed. The results provide the basis for a computer algorithm to predict the tertiary structure of an $\text{α}\text{β}$ protein from its amino acid sequence and actual assignment of secondary structure.The packing of an individual α-helix against a β-sheet generally involves two adjacent ± 4 rows of non-polar residues on the α-helix at the positions i, i + 4, i + 8, i + 1, i + 5, i + 9. The pattern of interacting β-sheet residues results from the twisted nature of the sheet surface and the attendant rotation of the side-chains. At a more detailed level, four of the α-helical residues (i + 1, i + 4, i + 5 and i + 8) form a diamond that surrounds one particular β-sheet residue, generally isoleucine, leucine or valine. In general, the α-helix sits 10 Å above the sheet and lies parallel to the strand direction.The prediction follows a combinational approach. First, a list of possible β-sheet structures (10⁶ to 10¹⁴) is constructed by the generation of all β-sheet topologies and β-strand alignments. This list is reduced by constraints on topology and the location of non-polar residues to mediate the sheet/helix packing, and then rank-ordered on the extent of hydrogen bonding. This algorithm was uniformly applied to 16 $\text{α}\text{β}$ domains in 13 proteins. For every structure, one member of the reduced list was close to the crystal structure; the root-mean-square deviation between equivalenced C^α atoms averaged 5.6 Å for 100 residues. For the $\text{α}\text{β}$ proteins with pure parallel β-sheets, the total number of structures comparable to or better than the native in terms of hydrogen bonds was between 1 and 148. For proteins with mixed β-sheets, the worst case is glyceraldehyde-3-phosphate dehydrogenase, where as many as 3800 structures would have to be sampled. The evolutionary significance of these results as well as the potential use of a combinatorial approach to the protein folding problem are discussed.     

Conformational and geometrical properties of β-sheets in proteins: II. Antiparallel and mixed β-sheets

The present work treats the conformational and geometrical properties of antiparallel and mixed β-sheets, following the approach described in the first paper of this series (Salemme & Weatherford, 1981). As shown, antiparallel structures possess conformational degrees of freedom that allow them to assume a greater diversity of spatial configurations than occur in parallel sheets. This configurational flexibility principally owes its origin to the potential variability of the interchain hydrogen bond geometry in antiparallel structures. As a consequence of this inherent elasticity of antiparallel β-sheets, their extended spatial configurations strongly reflect effects arising from local or extended side-chain packing interactions. Although the continuously deformable characteristics of antiparallel sheets frequently result in the attainment of spatial configurations that cannot be quantitatively modeled as hydrogen-bonded arrays of conformationally identical polypeptide chains, several observed sheets are nevertheless shown to be accurately approximated as conformationally regular arrays that typically yield model versus observed fits of less than 1 Å for corresponding α-carbon atoms.     

Quasirandom distribution of α and β regions in globular proteins

Distributions of α and β regions in globular proteins among clusters containing different numbers of adjacent α-helices, adjacent β regions, and “overlapping” βαβ units are considered. It is shown that these distributions do not differ greatly from what can be expected for random distributions of α and β regions along protein chains. In particular, it is shown that the amounts of relatively long α, β and βαβ clusters (which provide the basis for the conventional classification of globular proteins or domains into α, β, or α/β types) in random sequences also do not differ very much from those in real globular proteins. It follows that the possibility of structural classification of globular proteins (domains) does not imply the existence of a correlation in protein primary structure. This possibility exists even in random sequences of amino acid residues and therefore may not be the result of biological evolution.     

Physical�Cchemical determinants of coil conformations in globular proteins

We present a method with the potential to generate a library of coil segments from first principles. Proteins are built from α‐helices and/or β‐strands interconnected by these coil segments. Here, we investigate the conformational determinants of short coil segments, with particular emphasis on chain turns. Toward this goal, we extracted a comprehensive set of two‐, three‐, and four‐residue turns from X‐ray–elucidated proteins and classified them by conformation. A remarkably small number of unique conformers account for most of this experimentally determined set, whereas remaining members span a large number of rare conformers, many occurring only once in the entire protein database. Factors determining conformation were identified via Metropolis Monte Carlo simulations devised to test the effectiveness of various energy terms. Simulated structures were validated by comparison to experimental counterparts. After filtering rare conformers, we found that 98% of the remaining experimentally determined turn population could be reproduced by applying a hydrogen bond energy term to an exhaustively generated ensemble of clash‐free conformers in which no backbone polar group lacks a hydrogen‐bond partner. Further, at least 90% of longer coil segments, ranging from 5‐ to 20 residues, were found to be structural composites of these shorter primitives. These results are pertinent to protein structure prediction, where approaches can be divided into either empirical or ab initio methods. Empirical methods use database‐derived information; ab initio methods rely on physical–chemical principles exclusively. Replacing the database‐derived coil library with one generated from first principles would transform any empirically based method into its corresponding ab initio homologue.     

Formation of β-sheets in glutamine and alanine tripeptides

The misfolding and aggregation of proteins is associated with many different diseases including the trinucleotide repeat disorders and Prion diseases. We have studied three residue peptides comprising alanine and glutamine in order to understand the short range interactions affecting the formation of β-rich aggregates. Using infrared spectroscopy, we have found that trialanine and triglutamine form significant amounts of β-sheet, but that tripeptides containing alanine and glutamine are only able to form β-sheet if the glutamine side-chains extend outward on both faces of the sheet. From our data, we conclude that different stabilizing interactions are responsible for β-sheet formation in trialanine and triglutamine.     

How does sucrose stabilize the native state of globular proteins?

It is well known that sucrose stabilizes the native state of globular proteins against both chemical denaturants and temperature. A largely accepted explanation of sucrose-induced stabilization is not yet emerged. It is shown that the same theoretical approach able to rationalize the occurrence of cold denaturation, the contrasting role of GdmCl and Gdm₂SO₄, and the TMAO counteraction of urea denaturing activity [PCCP 12 (2010) 14245; PCCP 13 (2011) 12008; PCCP 13 (2011) 17689] works well also in the case of sucrose. The solvent-excluded volume effect plays the fundamental role because sucrose addition to water causes a marked increase in volume packing density due to the large size of sucrose molecules, that act as crowding agents.     

Evolution of proteins formed by β-sheets: I. Plastocyanin and azurin

Plastocyanin and azurin form a family of small copper-containing proteins, active in the electron transport systems of plants and bacteria, respectively. The crystal structures of two members of this family have been determined: poplar leaf plastocyanin and Pseudomonas aeruginosa azurin. Both proteins contain two β-sheets, packed face-to-face. Using computed superpositions of the structures, we have aligned the sequences, identified homologous positions, and studied how the structures have changed as a result of mutations.The residues in the vicinity of the copper-binding site show minimal amino acid substitution and form almost identical structures. Other portions of these proteins are more variable in sequence and in structure. Buried residues tend to maintain their hydrophobic character, but mutations change their volume. The mean variation in volume of homologous buried residues is 54 ų. The differences in size and shape of these buried residues are accommodated by a 3.8 Å shift in relative position of the packed β-sheets. This shift does not affect the copper binding site, because the residues that form this site are in, or adjacent to, just one of the β-sheets.     

Structural plasticity in self-assembling transmembrane β-sheets

Here we test the hypothesis that membrane-spanning β-sheets can exhibit structural plasticity in membranes due to their ability to shift hydrogen-bonding patterns. Transmembrane β-sheet forming peptides of the sequence AcWL_n, where n = 5, 6, or 7, which range from 21 to 27 Å in maximum length, were incorporated into bilayers made of phosphatidylcholine lipids with saturated acyl chains containing 14, 16, or 18 carbons, which are 36–50 Å in thickness. The effect of the peptide β-sheets on fluid- and gel-phase bilayers were studied with differential scanning calorimetry and circular dichroism spectroscopy. We show that AcWL₅ forms a stable, peptide-rich gel phase in all three lipids. The whole family of AcWL_n peptides appears to form similarly stable, nonmembrane-disrupting β-sheets in all bilayer phases and thicknesses. Bilayers containing up to 20 mol % peptide, which is the maximum concentration tested, formed gel phases with melting temperatures that were equal to, or slightly higher than, the pure lipid transitions. Given the range of peptide lengths and bilayer thicknesses tested, these experiments show that the AcWL_n family of membrane-inserted β-sheets exhibit remarkable structural plasticity in membranes.     

Chromosomal arrangement of the chicken β-type globin genes

We have isolated the chicken β-type globin genes from a library of chicken DNA-λ Charon 4A recombinant bacteriophage. There are four β-type genes within this segment of the genome; we believe this represents all of the β-type genes of the chicken. The recombinant λCβG1 contains the embryonic ?- and adult β-globin genes. The hatching β^H and embryonic p-globin genes are found in the recombinant λCβG2. Although λCβG1 and λCβG2 do not physically overlap, we present evidence that all four genes are closely linked and transcribed from the same DNA strand. These experiments demonstrate that the chromosomal regions represented by λCβG1 and λCβG2 lie approximately 1.6 kb apart in the chicken genome. A third recombinant λCβG3 extends the genomic locus studied in the vicinity of the β-type globin genes to approximately 39 kb. The physical order of the chicken β-type globin genes within this segment of the chromosome is 5′ … ?-β^H-β-? … 3′. This arrangement is unique among the vertebrate β-type globin gene clusters thus far examined, in that embryonic genes are located at the 5′ and 3′ ends of the cluster while the hatching and adult genes occupy central positions.     

Prediction of β-barrel membrane proteins by searching for restricted domains

Background  

The identification of β-barrel membrane proteins out of a genomic/proteomic background is one of the rapidly developing fields in bioinformatics. Our main goal is the prediction of such proteins in genome/proteome wide analyses.     

Prediction of the β-Hairpins in Proteins Using Support Vector Machine

By using of the composite vector with increment of diversity and scoring function to express the information of sequence, a support vector machine (SVM) algorithm for predicting β-hairpin motifs is proposed. The prediction is done on a dataset of 3,088 non homologous proteins containing 6,027 β-hairpins. The overall accuracy of prediction and Matthew’s correlation coefficient are 79.9% and 0.59 for the independent testing dataset. In addition, a higher accuracy of 83.3% and Matthew’s correlation coefficient of 0.67 in the independent testing dataset are obtained on a dataset previously used by Kumar et al. (Nuclic Acid Res 33:154–159). The performance of the method is also evaluated by predicting the β-hairpins of in the CASP6 proteins, and the better results are obtained. Moreover, this method is used to predict four kinds of supersecondary structures. The overall accuracy of prediction is 64.5% for the independent testing dataset.     

Relationship between structure and amino acid sequence of strongly twisted and coiled β-hairpins in globular proteins

β-Hairpins are widespread in proteins, and it is possible to find them both within β-sheets and separately. In this work, a comparative analysis of amino acid sequences of β-strands within strongly twisted β-hairpins from different structural protein subclasses has been conducted. Strongly twisted and coiled β-hairpin generates in the space a right double helix out of β-strands that are connected by a loop region (connections). The frequencies of amino acid residues on the internal (concave) and external (convex) surfaces of strongly twisted β-hairpins have been determined (220 β-hairpins from nonhomologous proteins were studied). The concave surface of these β-hairpins is mainly generated by hydrophobic residues, while the convex surface by hydrophilic residues; accordingly, the alternation of hydrophobic internal and hydrophilic external residues is observed in their amino acid sequences. Amino acid residues of glycine and alanine (especially in places of the largest twisting of the strands) were anomalously frequently found in internal positions of strongly twisted and coiled β-hairpins. It was established that internal positions never contain the proline residues, while external positions in the twisting region contain them in a relatively large amount. It was demonstrated that at least one amino acid residue in α_L- or ε-conformation is required for generation of relatively short (up to 7 amino acid residues) connection. As a rule, these positions are occupied by glycines. Thus, not only the alternation of hydrophobic and hydrophilic amino acid residues, but also the presence of one or two glycine residues in the connection region and the excess of glycines and alanines in the places of the largest strand twisting on the concave surface, as well as the presence of prolines on the convex surface, are required to generate a strongly twisted and coiled β-hairpin.     

Prediction of β-Turns

Kohonen's self-organization model, a neural network model, is applied to predict the -turns in proteins. There are 455 -turn tetrapeptides and 3807 non--turn tetrapeptides in the training database. The rates of correct prediction for the 110 -turn tetrapeptides and 30,229 non--turn tetrapeptides in the testing database are 81.8% and 90.7%, respectively. The high quality of prediction of neural network model implies that the residue-coupled effect along a polypeptide chain is important for the formation of reversal turns, such as -turns, during the process of protein folding.     

Single amino acid substitutions producing instability of globular proteins. Calculation of their frequencies in the entire mutational spectra of the α- and β-subunits of human hemoglobin

Summary The frequencies of substitutions resulting in protein instability were calculated by a method estimating changes in stability produced by amino acid substitutions. The method takes into account the accessibility of an amino acid position to a solvent and changes in the specificity of amino acid interactions. When tested on human mutant hemoglobins, the method yielded predictions with a preciseness of 80%. The consideration of the evolutionary homologous proteins in the analysis allowed us to estimate the evolutionary constraints imposed on stability of their spatial structure. With these limitations, approximately 50% of amino acid substitutions in the entire mutational spectra of the - and -subunits of human hemoglobin were found to damage the spatial structure of the globular proteins.     

Prediction of β-Turns

Kohonen's self-organization model, a neural network model, is applied to predict the β-turns in proteins. There are 455 β-turn tetrapeptides and 3807 non-β-turn tetrapeptides in the training database. The rates of correct prediction for the 110 β-turn tetrapeptides and 30,229 non-β-turn tetrapeptides in the testing database are 81.8% and 90.7%, respectively. The high quality of prediction of neural network model implies that the residue-coupled effect along a polypeptide chain is important for the formation of reversal turns, such as β-turns, during the process of protein folding.     

Conformational analysis and design of cross-strand disulfides in antiparallel β-sheets

Cross‐strand disulfides bridge two cysteines in a registered pair of antiparallel β‐strands. A nonredundant data set comprising 5025 polypeptides containing 2311 disulfides was used to study cross‐strand disulfides. Seventy‐six cross‐strand disulfides were found of which 75 and 1 occurred at non‐hydrogen‐bonded (NHB) and hydrogen‐bonded (HB) registered pairs, respectively. Conformational analysis and modeling studies demonstrated that disulfide formation at HB pairs necessarily requires an extremely rare and positive χ¹ value for at least one of the cysteine residues. Disulfides at HB positions also have more unfavorable steric repulsion with the main chain. Thirteen pairs of disulfides were introduced in NHB and HB pairs in four model proteins: leucine binding protein (LBP), leucine, isoleucine, valine binding protein (LIVBP), maltose binding protein (MBP), and Top7. All mutants LIVBP T247C V331C showed disulfide formation either on purification, or on treatment with oxidants. Protein stability in both oxidized and reduced states of all mutants was measured. Relative to wild type, LBP and MBP mutants were destabilized with respect to chemical denaturation, although the sole exposed NHB LBP mutant showed an increase of 3.1°C in T _m . All Top7 mutants were characterized for stability through guanidinium thiocyanate chemical denaturation. Both exposed and two of the three buried NHB mutants were appreciably stabilized. All four HB Top7 mutants were destabilized (ΔΔG ⁰ = ?3.3 to ?6.7 kcal/mol). The data demonstrate that introduction of cross‐strand disulfides at exposed NHB pairs is a robust method of improving protein stability. All four exposed Top7 disulfide mutants showed mild redox activity. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

An easy NMR method to study the formation of parallel β-sheets in peptide aggregates

Summary There is currently great interest in the study of peptide aggregation by β-sheet formation because of its relevance in pathological states or in the design of self-assembling systems of technological interest. NMR studies of β-sheet aggregates are difficult because of their long correlation times and spectral degeneracy. In this communication we demonstrate the combination of a semiselective TOCSY-NOESY experiment with partial deuterium exchange of labile protons to assign inter-molecular NOE cross peaks and prove the presence of a soluble parallel β-sheet in fast exchange with monomeric Ac-ASTTNYT-NH₂ (Ac-T-NH₂) in solution.