Identification of amino acids involved in protein structural uniqueness: implication for de novo protein design

Structural uniqueness is characteristic of native proteins and is essential to express their biological functions. The major factors that bring about the uniqueness are specific interactions between hydrophobic residues and their unique packing in the protein core. To find the origin of the uniqueness in their amino acid sequences, we analyzed the distribution of the side chain rotational isomers (rotamers) of hydrophobic amino acids in protein tertiary structures and derived deltaS(contact), the conformational-entropy changes of side chains by residue-residue contacts in each secondary structure. The deltaS(contact) values indicate distinct tendencies of the residue pairs to restrict side chain conformation by inter-residue contacts. Of the hydrophobic residues in alpha-helices, aliphatic residues (Leu, Val, Ile) strongly restrict the side chain conformations of each other. In beta-sheets, Met is most strongly restricted by contact with Ile, whereas Leu, Val and Ile are less affected by other residues in contact than those in alpha-helices. In designed and native protein variants, deltaS(contact) was found to correlate with the folding-unfolding cooperativity. Thus, it can be used as a specificity parameter for designing artificial proteins with a unique structure.     

Amino acid composition of protein termini are biased in different manners

An exhaustive statistical analysis of the amino acid sequences at the carboxyl (C) and amino (N) termini of proteins and of coding nucleic acid sequences at the 5' side of the stop codons was undertaken. At the N ends, Met and Ala residues are over-represented at the first (+1) position whereas at positions 2 and 5 Thr is preferred. These peculiarities at N-termini are most probably related to the mechanism of initiation of translation (for Met) and to the mechanisms governing the life-span of proteins via regulation of their degradation (for Ala and Thr). We assume that the C-terminal bias facilitates fixation of the C ends on the protein globule by a preference for charged and Cys residues. The terminal biases, a novel feature of protein structure, have to be taken into account when molecular evolution, three-dimensional structure, initiation and termination of translation, protein folding and life-span are concerned. In addition, the bias of protein termini composition is an important feature which should be considered in protein engineering experiments.     

Principles of protein folding--a perspective from simple exact models.

General principles of protein structure, stability, and folding kinetics have recently been explored in computer simulations of simple exact lattice models. These models represent protein chains at a rudimentary level, but they involve few parameters, approximations, or implicit biases, and they allow complete explorations of conformational and sequence spaces. Such simulations have resulted in testable predictions that are sometimes unanticipated: The folding code is mainly binary and delocalized throughout the amino acid sequence. The secondary and tertiary structures of a protein are specified mainly by the sequence of polar and nonpolar monomers. More specific interactions may refine the structure, rather than dominate the folding code. Simple exact models can account for the properties that characterize protein folding: two-state cooperativity, secondary and tertiary structures, and multistage folding kinetics--fast hydrophobic collapse followed by slower annealing. These studies suggest the possibility of creating "foldable" chain molecules other than proteins. The encoding of a unique compact chain conformation may not require amino acids; it may require only the ability to synthesize specific monomer sequences in which at least one monomer type is solvent-averse.     

Genetic code redundancy and the evolutionary stability of protein secondary structure

The genetic code has an inherent bias towards some amino acids because of the variable number of synonymous codons per amino acid. The extent to which these biases are expressed in protein secondary structure is described through the analysis of the overall amino acid compositions of the alpha-helix, beta-sheet, beta-turn and random coil segments elucidated by X-ray crystallography. Given the concept of neutral mutation in proteins, the allocation of synonyms in the genetic code appears to protect secondary structures from amino acid changes and discourages the appearance of chemically complex residues. The level of protection is similar for each structural form, despite their clear preferences for certain amino acids. The organization of the code is therefore relevant to the preservation of conformation seen in the evolution of many protein families.     

The behaviour of polyamino acids reveals an inverse side chain effect in amyloid structure formation

Amyloid fibrils and prions are proteinaceous aggregates that are based on a unique form of polypeptide configuration, termed cross-beta structure. Using a group of chemically distinct polyamino acids, we show here that the existence of such a structure does not require the presence of specific side chain interactions or sequence patterns. These observations firmly establish that amyloid formation and protein folding represent two fundamentally different ways of organizing polypeptides into ordered conformations. Protein folding depends critically on the presence of distinctive side chain sequences and produces a unique globular fold. By contrast, the properties of different polyamino acids suggest that amyloid formation arises primarily from main chain interactions that are, in some environments, overruled by specific side chain contacts. This side chain effect can be thought of as the inverse of the one that characterizes protein folding. Conditions including Alzheimer's and Creutzfeldt-Jakob diseases represent, on this basis, pathological cases in which a natural polypeptide chain has aberrantly adopted the conformation that is primarily defined by main chain interactions and not the structure that is determined by specific side chain contacts that depend on the polypeptide sequence.     

Electronic properties of the amino acid side chains contribute to the structural preferences in protein folding.

A database of 118 non-redundant proteins was examined to determine the preferences of amino acids for secondary structures: alpha-helix, beta-strand and coil conformations. To better understand how the physicochemical properties of amino acid side chains might influence protein folding, several new scales have been suggested for quantifying the electronic effects of amino acids. These include the pKa at the amino group, localized effect substituent constants (esigma), and a composite of these two scales (epsilon). Amino acids were also classified into 5 categories on the basis of their electronic properties: O (strong electron donor), U (weak donor), Z (ambivalent), B (weak electron acceptor), and X (strong acceptor). Certain categories of amino acid appeared to be critical for particular conformations, e.g., O and U-type residues for alpha-helix formation. Pairwise analysis of the database according to these categories revealed significant context effects in the structural preferences. In general, the propensity of an amino acid for a particular conformation was related to the electronic features of the side chain. Linear regression analyses revealed that the electronic properties of amino acids contributed about as much to the folding preferences as hydrophobicity, which is a well-established determinant of protein folding. A theoretical model has been proposed to explain how the electronic properties of the side chain groups might influence folding along the peptide backbone.     

Amino acid propensities for the collagen triple-helix

Determination of the tendencies of amino acids to form alpha-helical and beta-sheet structures has been important in clarifying stabilizing interactions, protein design, and the protein folding problem. In this study, we have determined for the first time a complete scale of amino acid propensities for another important protein motif: the collagen triple-helix conformation with its Gly-X-Y repeating sequence. Guest triplets of the form Gly-X-Hyp and Gly-Pro-Y are used to quantitate the conformational propensities of all 20 amino acids for the X and Y positions in the context of a (Gly-Pro-Hyp)(8) host peptide. The rankings for both the X and Y positions show the highly stabilizing nature of imino acids and the destabilizing effects of Gly and aromatic residues. Many residues show differing propensities in the X versus Y position, related to the nonequivalence of these positions in terms of interchain interactions and solvent exposure. The propensity of amino acids to adopt a polyproline II-like conformation plays a role in their triple-helix rankings, as shown by a moderate correlation of triple-helix propensity with frequency of occurrence in polyproline II-like regions. The high propensity of ionizable residues in the X position suggests the importance of interchain hydrogen bonding directly or through water to backbone carbonyls or hydroxyprolines. The low propensity of side chains with branching at the C(delta) in the Y position supports models suggesting these groups block solvent access to backbone C=O groups. These data provide a first step in defining sequence-dependent variations in local triple-helix stability and binding, and are important for a general understanding of side chain interactions in all proteins.     

Predicted Folding of β-Structure in Myelin Basic Protein

Predictions of myelin basic protein secondary structure have not previously considered a major role for beta-structure in the organization of the native molecule because optical rotatory dispersion and circular dichroism studies have provided little, if any, evidence for beta-structure, and because a polycationic protein is generally considered to resist folding into a compact structure. However, the Chou-Fasman, Lim, and Robson algorithms identify a total of five beta-strands in the amino acid sequence. Four of these hydrophobic amino acid sequences (37-45, 87-95, 110-118, and 150-158) could form a hairpin intermediate that initiates folding of a Greek-key-type beta-structure. A second fold on the more hydrophobic side, with the addition of a strand from the N-terminus (residues 13-21), would complete the five-stranded antiparallel beta-sheet. A unique strand alignment can be predicted by phasing the hydrophobic residues. The unusual triproline sequence of myelin basic protein (100-102) is enclosed in the 14-residue hairpin loop. If these prolines are in the trans conformation, models show that a reverse turn could occur at residues 102-105 (Pro-Ser-Gln-Gly). Algorithms do not agree on the prediction of alpha-helices, but each of the two large loops could accommodate an alpha-helix. Myelin basic protein is known to be phosphorylated in vivo on as many as five Ser/Thr residues. Phosphorylation might alter the dynamics of folding if the nascent polypeptide were phosphorylated in the cytoplasm. In particular, phosphorylation of Thr-99 could neutralize cationic residues Lys-106 and Arg-108 within the hairpin loop. In addition, the methylation of Arg-108 might stabilize the hairpin loop structure through hydrophobic interaction with the side chain of Pro-97. The cationic side chains of arginine and lysine residues located on the faces of the beta-sheet (Arg-43, Arg-114, Lys-13, Lys-92, Lys-153, and Lys-156) could provide sites for interaction with phospholipids and other anionic structures on the surface of the myelin lipid bilayer.     

A computer model to dynamically simulate protein folding: studies with crambin

The current work describes a simplified representation of protein structure with uses in the simulation of protein folding. The model assumes that a protein can be represented by a freely rotating rigid chain with a single atom approximating the effect of each side chain. Potentials describing the attraction or repulsion between different types of amino acids are determined directly from the distribution of amino acids in the database of known protein structures. The optimization technique of simulated annealing has been used to dynamically sample the conformations available to this simple model, allowing the protein to evolve from an extended, random coil into a compact globular structure. Many characteristics expected of true proteins, such as the sequence-dependent formation of secondary structure, the partitioning of hydrophobic residues, and specific disulfide pairing, are reproduced by the simulation, suggesting the model may accurately simulate the folding process.     

The DeltaF508 cystic fibrosis mutation impairs domain-domain interactions and arrests post-translational folding of CFTR

Misfolding accounts for the endoplasmic reticulum-associated degradation of mutant cystic fibrosis transmembrane conductance regulators (CFTRs), including deletion of Phe508 (DeltaF508) in the nucleotide-binding domain 1 (NBD1). To study the role of Phe508, the de novo folding and stability of NBD1, NBD2 and CFTR were compared in conjunction with mutagenesis of Phe508. DeltaF508 and amino acid replacements that prevented CFTR folding disrupted the NBD2 fold and its native interaction with NBD1. DeltaF508 caused limited alteration in NBD1 conformation. Whereas nonpolar and some aliphatic residues were permissive, charged residues and glycine compromised the post-translational folding and stability of NBD2 and CFTR. The results suggest that hydrophobic side chain interactions of Phe508 are required for vectorial folding of NBD2 and the domain-domain assembly of CFTR, representing a combined co- and post-translational folding mechanism that may be used by other multidomain membrane proteins.     

Electronic Properties of the Amino Acid Side Chains Contribute to the Structural Preferences in Protein Folding

Abstract

A database of 118 non-redundant proteins was examined to determine the preferences of amino acids for secondary structures: α-helix, β-strand and coil conformations. To better understand how the physicochemical properties of amino acid side chains might influence protein folding, several new scales have been suggested for quantifying the electronic effects of amino acids. These include the pKa at the amino group, localized effect substituent constants (eσ), and a composite of these two scales (ε). Amino acids were also classified into 5 categories on the basis of their electronic properties: O (strong electron donor), U (weak donor), Z (ambivalent), B (weak electron acceptor), and X (strong acceptor). Certain categories of amino acid appeared to be critical for particular conformations, e.g., O and U-type residues for α-helix formation. Pairwise analysis of the database according to these categories revealed significant context effects in the structural preferences. In general, the propensity of an amino acid for a particular conformation was related to the electronic features of the side chain. Linear regression analyses revealed that the electronic properties of amino acids contributed about as much to the folding preferences as hydrophobicity, which is a well-established determinant of protein folding. A theoretical model has been proposed to explain how the electronic properties of the side chain groups might influence folding along the peptide backbone.     

Frequencies of amino acid strings in globular protein sequences indicate suppression of blocks of consecutive hydrophobic residues

Patterns of hydrophobic and hydrophilic residues play a major role in protein folding and function. Long, predominantly hydrophobic strings of 20-22 amino acids each are associated with transmembrane helices and have been used to identify such sequences. Much less attention has been paid to hydrophobic sequences within globular proteins. In prior work on computer simulations of the competition between on-pathway folding and off-pathway aggregate formation, we found that long sequences of consecutive hydrophobic residues promoted aggregation within the model, even controlling for overall hydrophobic content. We report here on an analysis of the frequencies of different lengths of contiguous blocks of hydrophobic residues in a database of amino acid sequences of proteins of known structure. Sequences of three or more consecutive hydrophobic residues are found to be significantly less common in actual globular proteins than would be predicted if residues were selected independently. The result may reflect selection against long blocks of hydrophobic residues within globular proteins relative to what would be expected if residue hydrophobicities were independent of those of nearby residues in the sequence.     

Structural insights into charge pair interactions in triple helical collagen-like proteins

The collagen triple helix is the most abundant protein fold in humans. Despite its deceptively simple structure, very little is understood about its folding and fibrillization energy landscape. In this work, using a combination of x-ray crystallography and nuclear magnetic resonance spectroscopy, we carry out a detailed study of stabilizing pair-wise interactions between the positively charged lysine and the negatively charged amino acids aspartate and glutamate. We find important differences in the side chain conformation of amino acids in the crystalline and solution state. Structures from x-ray crystallography may have similarities to the densely packed triple helices of collagen fibers whereas solution NMR structures reveal the simpler interactions of isolated triple helices. In solution, two distinct types of contacts are observed: axial and lateral. Such register-specific interactions are crucial for the understanding of the registration process of collagens and the overall stability of proteins in this family. However, in the crystalline state, there is a significant rearrangement of the side chain conformation allowing for packing interactions between adjacent helices, which suggests that charged amino acids may play a dual role in collagen stabilization and folding, first at the level of triple helical assembly and second during fibril formation.     

Correlation between amino acid residues converted by RNA editing and functional residues in protein three-dimensional structures in plant organelles

Background

In plant organelles, specific messenger RNAs (mRNAs) are subjected to conversion editing, a process that often converts the first or second nucleotide of a codon and hence the encoded amino acid. No systematic patterns in converted sites were found on mRNAs, and the converted sites rarely encoded residues located at the active sites of proteins. The role and origin of RNA editing in plant organelles remain to be elucidated.

Results

Here we study the relationship between amino acid residues encoded by edited codons and the structural characteristics of these residues within proteins, e.g., in protein-protein interfaces, elements of secondary structure, or protein structural cores. We find that the residues encoded by edited codons are significantly biased toward involvement in helices and protein structural cores. RNA editing can convert codons for hydrophilic to hydrophobic amino acids. Hence, only the edited form of an mRNA can be translated into a polypeptide with helix-preferring and core-forming residues at the appropriate positions, which is often required for a protein to form a functional three-dimensional (3D) structure.

Conclusion

We have performed a novel analysis of the location of residues affected by RNA editing in proteins in plant organelles. This study documents that RNA editing sites are often found in positions important for 3D structure formation. Without RNA editing, protein folding will not occur properly, thus affecting gene expression. We suggest that RNA editing may have conferring evolutionary advantage by acting as a mechanism to reduce susceptibility to DNA damage by allowing the increase in GC content in DNA while maintaining RNA codons essential to encode residues required for protein folding and activity.     

Identification of sites influencing the folding and subunit assembly of the P22 tailspike polypeptide chain using nonsense mutations

Amber mutations have been isolated and mapped to more than 60 sites in gene 9 of P22 encoding the thermostable phage tailspike protein. Gene 9 is the locus of over 30 sites of temperature sensitive folding (tsf) mutations, which affect intermediates in the chain folding and subunit association pathway. The phenotypes of the amber missense proteins produced on tRNA suppressor hosts inserting serine, glutamine, tryosine and leucine have been determined at different temperatures. Thirty-three of the sites are tolerant, producing functional proteins with any of the four amino acids inserted at the sites, independent of temperature. Tolerant sites are concentrated at the N-terminal end of the protein indicating that this region is not critical for conformation or function. Sixteen of the sites yield temperature sensitive missense proteins on at least one nonsense suppressing host. Most of the sites with ts phenotypes map to the central region of the gene which is also the region where most of the tsf mutations map. Mutations at 15 of the sites have a lethal phenotype on at least one tRNA suppressor host. For nine out of ten sites tested with at least one lethal phenotype, the primary defect was in the folding or subunit association of the missense polypeptide chain. This analysis of the tailspike missense proteins distinguishes three classes of amino acid sites in the polypeptide chain; residues whose side chains contribute little to folding, subunit assembly or function; residues critical for maintaining the folding and subunit assembly pathway at the high end of the temperature range of phage growth; and residues critical over the entire temperature range of growth.     

Determinants of protein stability and folding: comparative analysis of beta-lactoglobulins and liver basic fatty acid binding protein

A new energy decomposition approach, aimed at identifying residues playing a folding key role, has been applied here to three homologous proteins, belonging to the calycin superfamily, namely bovine and porcine beta-lactoglobulins and Liver basic fatty acid binding protein, sharing the same beta-barrel fold and different degree of sequence identities. All-atom, explicit solvent molecular dynamics simulations around the native conformation were used to generate, for each of the three proteins, energy maps which were further simplified through eigenvalue decomposition. Analysis of the components of the eigenvector associated with the lowest eigenvalue singled out those residues (hot sites) behaving as strongly interacting and possible nucleation centers. The results fit well with experimental folding data and, especially, with the analysis of side chain-side chain interaction conservation.     

Mapping and characterization of the entomocidal domain of the Bacillus thuringiensis CrylA(b) protoxin

The amino acid sequences necessary for entomocidal activity of the CryIA(b) protoxin of Bacillus thuringiensis were determined. Introduction of stop codons behind codons Arg601, Phe604 or Ala607 showed that amino acid residues C-terminal to Ala607 are not required for insecticidal activity and that activation by midgut proteases takes place distal to Ala607. The two shortest polypeptides, deleted for part of the highly conserved β-strand, were prone to proteolytic degradation, explaining their lack of toxicity. Apparently, this β-strand is essential for folding of the molecule into a stable conformation. Proteolytic activation at the N-terminus was investigated by removing the first 28 codons, resulting in a translation product extending from amino acid 29 to 607. This protein appeared to be toxic not only to susceptible insect larvae such as Manduca sexta and Heliothis virescens, but also to Escherichia coli cells. An additional mutant, encoding only amino acid residues 29–429, encompassing the complete putative pore forming domain, but lacking a large part of the receptor-binding domain, was similarly toxic to E. coli cells. This suggests a role for the N-terminal 28 amino acids in rendering the toxin inactive in Bacillus thuringiensis, and indicates that the cytolytic potential of the pore forming domain is only realized after proteolytic removal of these residues by proteases in the insect gut. In line with this hypothesis are results obtained with a mutant protein in which Arg28 at the cleavage site was replaced by Asp. This substitution prevented the protein from being cleaved by trypsin in vitro, and reduced its toxicity to M. sexta larvae.     

Protein structure and the sequential structure of mRNA: α-Helix and β-sheet signals at the nucleotide level

A direct comparison of experimentally determined protein structures and their corresponding protein coding mRNA sequences has been performed. We examine whether real world data support the hypothesis that clusters of rare codons correlate with the location of structural units in the resulting protein. The degeneracy of the genetic code allows for a biased selection of codons which may control the translational rate of the ribosome, and may thus in vivo have a catalyzing effect on the folding of the polypeptide chain. A complete search for GenBank nucleotide sequences coding for structural entries in the Brookhaven Protein Data Bank produced 719 protein chains with matching mRNA sequence, amino acid sequence, and secondary structure assignment. By neural network analysis, we found strong signals in mRNA sequence regions surrounding helices and sheets. These signals do not originate from the clustering of rare codons, but from the similarity of codons coding for very abundant amino acid residues at the N- and C-termini of helices and sheets. No correlation between the positioning of rare codons and the location of structural units was found. The mRNA signals were also compared with conserved nucleotide features of 16S-like ribosomal RNA sequences and related to mechanisms for maintaining the correct reading frame by the ribosome. © 1996 Wiley-Liss, Inc.     

Exploring synonymous codon usage preferences of disulfide-bonded and non-disulfide bonded cysteines in the E. coli genome

High-quality data about protein structures and their gene sequences are essential to the understanding of the relationship between protein folding and protein coding sequences. Firstly we constructed the EcoPDB database, which is a high-quality database of Escherichia coli genes and their corresponding PDB structures. Based on EcoPDB, we presented a novel approach based on information theory to investigate the correlation between cysteine synonymous codon usages and local amino acids flanking cysteines, the correlation between cysteine synonymous codon usages and synonymous codon usages of local amino acids flanking cysteines, as well as the correlation between cysteine synonymous codon usages and the disulfide bonding states of cysteines in the E. coli genome. The results indicate that the nearest neighboring residues and their synonymous codons of the C-terminus have the greatest influence on the usages of the synonymous codons of cysteines and the usage of the synonymous codons has a specific correlation with the disulfide bond formation of cysteines in proteins. The correlations may result from the regulation mechanism of protein structures at gene sequence level and reflect the biological function restriction that cysteines pair to form disulfide bonds. The results may also be helpful in identifying residues that are important for synonymous codon selection of cysteines to introduce disulfide bridges in protein engineering and molecular biology. The approach presented in this paper can also be utilized as a complementary computational method and be applicable to analyse the synonymous codon usages in other model organisms.     

Comparison of scales of amino acid side chain properties by conservation during evolution of four proteins

As the amino acid sequence of a given protein changes along the phylogenetic tree, enough of the overall folding pattern must be conserved to ensure that the protein still fulfils its biological function. Eighteen published scales which tabulate various side chain properties are compared here by computing the variance of each scale when applied to each of several protein families. The conservation of each scale of side chain properties is examined for the 20,627 residues in 60 mammalian myoglobins, 31 mammalian ribonucleases, insulin A and B chains (29 sequences each), 29 vertebrate and 28 plant cytochrome c's. Those scales which are the most highly conserved through the evolution of each protein family may well be the best predictors of protein folding patterns. The mean-area-buried scale and the optimized matching hydrophobicities scale are more conserved than other scales. An additional result is the relatively poor conservation across evolution of the Chou-Fasman secondary structure predictors.