Backbones of Folded Proteins Reveal Novel Invariant Amino Acid Neighborhoods

Abstract

Folding of naturally occurring proteins has eluded a universal molecular level explanation till date. Rather, there is an abundance of diverse views on dominant factors governing protein folding. Through rigorous analyses of several thousand crystal structures, we observe that backbones of folded proteins display some remarkable invariant features. Folded proteins are characterized by spatially well-defined, distance dependent, and universal, neighborhoods of amino acids which defy any of the conventionally prevalent views. These findings present a compelling case for a newer view of protein folding which takes into account solvent mediated and amino acid shape and size assisted optimization of the tertiary structure of the polypeptide chain to make a functional protein.     

Single amino acid substitutions globally suppress the folding defects of temperature-sensitive folding mutants of phage P22 coat protein.

The amino acid sequence of a polypeptide defines both the folding pathway and the final three-dimensional structure of a protein. Eighteen amino acid substitutions have been identified in bacteriophage P22 coat protein that are defective in folding and cause their folding intermediates to be substrates for GroEL and GroES. These temperature-sensitive folding (tsf) substitutions identify amino acids that are critical for directing the folding of coat protein. Additional amino acid residues that are critical to the folding process of P22 coat protein were identified by isolating second site suppressors of the tsf coat proteins. Suppressor substitutions isolated from the phage carrying the tsf coat protein substitutions included global suppressors, which are substitutions capable of alleviating the folding defects of numerous tsf coat protein mutants. In addition, potential global and site-specific suppressors were isolated, as well as a group of same site amino acid substitutions that had a less severe phenotype than the tsf parent. The global suppressors were located at positions 163, 166, and 170 in the coat protein sequence and were 8-190 amino acid residues away from the tsf parent. Although the folding of coat proteins with tsf amino acid substitutions was improved by the global suppressor substitutions, GroEL remained necessary for folding. Therefore, we believe that the global suppressor sites identify a region that is critical to the folding of coat protein.     

Energy barriers, cooperativity, and hidden intermediates in the folding of small proteins

Current theoretical views of the folding process of small proteins (< approximately 100 amino acids) postulate that the landscape of potential mean force (PMF) for the formation of the native state has a funnel shape and that the free energy barrier to folding arises from the chain configurational entropy only. However, recent theoretical studies on the formation of hydrophobic clusters with explicit water suggest that a barrier should exist on the PMF of folding, consistent with the fact that protein folding generally involves a large positive activation enthalpy at room temperature. In addition, high-resolution structural studies of the hidden partially unfolded intermediates have revealed the existence of non-native interactions, suggesting that the correction of the non-native interactions during folding should also lead to barriers on PMF. To explore the effect of a PMF barrier on the folding behavior of proteins, we modified Zwanzig's model for protein folding with an uphill landscape of PMF for the formation of transition states. We found that the modified model for short peptide segments can satisfy the thermodynamic and kinetic criteria for an apparently two-state folding. Since the Levinthal paradox can be solved by a stepwise folding of short peptide segments, a landscape of PMF with a locally uphill search for the transition state and cooperative stabilization of folding intermediates/native state is able to explain the available experimental results for small proteins. We speculate that the existence of cooperative hidden folding intermediates in small proteins could be the consequence of the highly specific structures of the native state, which are selected by evolution to perform specific functions and fold in a biologically meaningful time scale.     

Amino acid sequence predicts folding rate for middle-size two-state proteins

The significant correlation between protein folding rates and the sequence-predicted secondary structure suggests that folding rates are largely determined by the amino acid sequence. Here, we present a method for predicting the folding rates of proteins from sequences using the intrinsic properties of amino acids, which does not require any information on secondary structure prediction and structural topology. The contribution of residue to the folding rate is expressed by the residue's Omega value. For a given residue, its Omega depends on the amino acid properties (amino acid rigidity and dislike of amino acid for secondary structures). Our investigation achieves 82% correlation with folding rates determined experimentally for simple, two-state proteins studied until the present, suggesting that the amino acid sequence of a protein is an important determinant of the protein-folding rate and mechanism.     

The relationship between synonymous codon usage and protein structure in Escherichia coli and Homo sapiens

The role of silent position in the codon on the protein structure is an interesting and yet unclear problem. In this paper, 563 Homo sapiens genes and 417 Escherichia coli genes coding for proteins with four different folding types have been analyzed using variance analysis, a multivariate analysis method newly used in codon usage analysis, to find the correlation between amino acid composition, synonymous codon, and protein structure in different organisms. It has been found that in E. coli, both amino acid compositions in differently folded proteins and synonymous codon usage in different gene classes coding for differently folded proteins are significantly different. It was also found that only amino acid composition is different in different protein classes in H. sapiens. There is no universal correlation between synonymous codon usage and protein structure in these two different organisms. Further analysis has shown that GC content on the second codon position can distinguish coding genes for different folded proteins in both organisms.     

The folding type of a protein is relevant to the amino acid composition

The folding types of 135 proteins, the three-dimensional structures of which are known, were analyzed in terms of the amino acid composition. The amino acid composition of a protein was expressed as a point in a multidimensional space spanned with 20 axes, on which the corresponding contents of 20 amino acids in the protein were represented. The distribution pattern of proteins in this composition space was examined in relation to five folding types, alpha, beta, alpha/beta, alpha + beta, and irregular type. The results show that amino acid compositions of the alpha, beta, and alpha/beta types are located in different regions in the composition space, thus allowing distinct separation of proteins depending on the folding types. The points representing proteins of the alpha + beta and irregular types, however, are widely scattered in the space, and the existing regions overlap with those of the other folding types. A simple method of utilizing the "distance" in the space was found to be convenient for classification of proteins into the five folding types. The assignment of the folding type with this method gave an accuracy of 70% in the coincidence with the experimental data.     

The analysis of protein folding kinetic data produced in protein engineering experiments

Over the past decade, the "protein engineering method" has been used to investigate the folding pathways of more than 20 different proteins. This method involves measuring the folding and unfolding rates of mutant proteins with single amino acid substitutions spread across the sequence. Comparison of folding rates of the mutant proteins to that of the wild-type protein allows the calculation of the phi value, which can be used to evaluate the stabilizing contribution of an amino acid side chain to the structure of the folding transition state. Here, we review the methodology for analysing data collected in protein engineering folding kinetics studies. We discuss the calculation of folding rates and kinetic m values, the estimation of errors in folding kinetics experiments, phi value calculation including potential pitfalls of the analysis, Br?nsted plots, detecting Hammond behaviour, and the analysis of curved chevron plots.     

不同类型氨基酸网络参量与蛋白质折叠的关系

理论和实验研究表明,蛋白质天然拓扑结构对其折叠过程具有重要的影响．采用复杂网络的方法分析蛋白质天然结构的拓扑特征,并探索蛋白质结构特征与折叠速率之间的内在联系．分别构建了蛋白质氨基酸网络、疏水网、亲水网、亲水-疏水网以及相应的长程网络,研究了这些网络的匹配系数(assortativity coefficient)和聚集系数(clustering coefficient)的统计特性．结果表明,除了亲水-疏水网,上述各网络的匹配系数均为正值,并且氨基酸网和疏水网的匹配系数与折叠速率表现出明显的线性正相关,揭示了疏水残基间相互作用的协同性有助于蛋白质的快速折叠．同时,研究发现疏水网的聚集系数与折叠速率有明显的线性负相关关系,这表明疏水残基间三角结构(triangle construction)的形成不利于蛋白质快速折叠．还进一步构建了相应的长程网络,发现序列上间距较远的残基接触对的形成将使蛋白质折叠进程变慢．     

Compactness determines protein folding type

We have demonstrated here that protein compactness, which we define as the ratio of the accessible surface area of a protein to that of the ideal sphere of the same volume, is one of the factors determining the mechanism of protein folding. Proteins with multi-state kinetics, on average, are more compact (compactness is 1.49 +/- 0.02 for proteins within the size range of 101-151 amino acid residues) than proteins with two-state kinetics (compactness is 1.59 +/- 0.03 for proteins within the same size range of 101-151 amino acid residues). We have shown that compactness for homologous proteins can explain both the difference in folding rates and the difference in folding mechanisms.     

Protein folding and misfolding: mechanism and principles

Two fundamentally different views of how proteins fold are now being debated. Do proteins fold through multiple unpredictable routes directed only by the energetically downhill nature of the folding landscape or do they fold through specific intermediates in a defined pathway that systematically puts predetermined pieces of the target native protein into place? It has now become possible to determine the structure of protein folding intermediates, evaluate their equilibrium and kinetic parameters, and establish their pathway relationships. Results obtained for many proteins have serendipitously revealed a new dimension of protein structure. Cooperative structural units of the native protein, called foldons, unfold and refold repeatedly even under native conditions. Much evidence obtained by hydrogen exchange and other methods now indicates that cooperative foldon units and not individual amino acids account for the unit steps in protein folding pathways. The formation of foldons and their ordered pathway assembly systematically puts native-like foldon building blocks into place, guided by a sequential stabilization mechanism in which prior native-like structure templates the formation of incoming foldons with complementary structure. Thus the same propensities and interactions that specify the final native state, encoded in the amino-acid sequence of every protein, determine the pathway for getting there. Experimental observations that have been interpreted differently, in terms of multiple independent pathways, appear to be due to chance misfolding errors that cause different population fractions to block at different pathway points, populate different pathway intermediates, and fold at different rates. This paper summarizes the experimental basis for these three determining principles and their consequences. Cooperative native-like foldon units and the sequential stabilization process together generate predetermined stepwise pathways. Optional misfolding errors are responsible for 3-state and heterogeneous kinetic folding.     

Are the same or different amino acid residues responsible for correct and incorrect protein folding?

It has been shown for 20 proteins that amino acid residues included into the protein folding nucleus, determined experimentally, are often involved in the theoretically determined amyloidogenic fragments. For 18 proteins, Φ-values indicative of the extent of residue involvement into the folding nucleus are on average higher for amino acid residues within amyloidogenic regions. Amyloidogenic fragments were predicted for 20 proteins by two methods chosen from four on the basis of comparison of prediction of amyloidogenic regions known from experimental data. Since theoretical folding nuclei are detected by the protein three-dimensional structure and amyloidogenic regions by the protein chain primary structure, the detected regularity makes possible predictions of folding nucleation sites on the basis of amino acid sequence.     

Native protein sequences are designed to destabilize folding intermediates

Hydrophobic core mutants of sperm whale apomyoglobin were constructed to investigate the amino acid sequence features that determine the folding properties. Replacements of all of the Ile residues with Leu and of all of the Ile and Val residues with Leu decreased the thermodynamic stability of the folded states against the unfolded states but increased the stability of the folding intermediates against the unfolded states, indicating that the amino acid composition of the protein core is important for the protein stability and folding cooperativity. To examine the effect of the arrangement of these hydrophobic residues, mutant proteins were further constructed: 12 sites out of the 18 Leu, 9 Ile, and 8 Val residues of the wild-type myoglobin were randomly replaced with each other so that the amino acid compositions were similar to that of the wild-type protein. Four mutant proteins were obtained without selection of the protein properties. These residue replacements similarly resulted in the stabilization of both the intermediate and folded states against the unfolded states, as compared to the wild-type protein. Thus, the arrangements of the hydrophobic residues in the native amino acid sequence are selected to destabilize the folding intermediate rather than to stabilize the folded state. The present results suggest that the two-state transition of protein folding or the transient formation of the unstable intermediate, which seems to be required for effective production of the functional proteins, has been a major driving force in the molecular evolution of natural globular proteins.     

The protein side of the central dogma: permanence and change

There are two facets to the central dogma proposed by Francis Crick in 1957. One concerns the relation between the sequence of nucleotides and the sequence of amino acids, the second is devoted to the relation between the sequence of amino acids and the native three-dimensional structure of proteins. 'Folding is simply a function of the order of the amino acids,' i.e. no information is required for the proper folding of a protein other than the information contained in its sequence. This protein side of the central dogma was elaborated in a scientific context in which the characteristics and functions of proteins, and the mechanisms of protein folding, were seen very differently. This context, which made the folding problem a simple one, supported the bold proposition of Francis Crick. The protein side of the central dogma was not challenged by the discovery of prions if one adopts the definition of information given by Francis Crick. It might have been challenged by the discovery that regulatory enzymes exist in different conformations, and the evidence for the existence of chaperones assisting protein folding. But it was not, and folding remains what it was for Francis Crick, 'simply a function of the order of amino acids'. But the meaning of 'function' has dramatically changed. It is no longer the result of simple physicochemical laws, but that of a long evolutionary process which has optimized protein folding. Molecular mechanistic explanations have to be allied with evolutionary explanations, in a way characteristic of present biology.     

GroEL binds a late folding intermediate of phage P22 coat protein

GroEL recognizes proteins that are folding improperly or that have aggregation-prone intermediates. Here we have used as substrates for GroEL, wildtype (WT) coat protein of phage P22 and 3 coat proteins that carry single amino acid substitutions leading to a temperature-sensitive folding (tsf) phenotype. In vivo, WT coat protein does not require GroEL for proper folding, whereas GroEL is necessary for the folding of the tsf coat proteins; thus, the single amino acid substitutions cause coat protein to become a substrate for GroEL. The conformation of WT and tsf coat proteins when in a binary complex with GroEL was investigated using tryptophan fluorescence, quenching of fluorescence, and accessibility of the coat proteins to proteolysis. WT coat protein and the tsf coat protein mutants were each found to be in a different conformation when bound to GroEL. As an additional measure of the changes in the bound conformation, the affinity of binding of WT and tsf coat proteins to GroEL was determined using a fluorescence binding assay. The tsf coat proteins were bound more tightly by GroEL than WT coat protein. Therefore, even though the proteins are identical except for a single amino acid substitution, GroEL did not bind these substrate polypeptides in the same conformation within its central cavity. Therefore, GroEL is likely to bind coat protein in a conformation consistent with a late folding intermediate, with substantial secondary and tertiary structure formed.     

Bacterial proteins fold faster than eukaryotic proteins with simple folding kinetics

Protein domain frequency and distribution among kingdoms was statistically analyzed using the SCOP structural database. It appeared that among chosen protein domains with the best resolution, eukaryotic proteins more often belong to α-helical and β-structural proteins, while proteins of bacterial origin belong to α/β structural class. Statistical analysis of folding rates of 73 proteins with known experimental data revealed that bacterial proteins with simple kinetics (23 proteins) exhibit a higher folding rate compared to eukaryotic proteins with simple folding kinetics (27 proteins). Analysis of protein domain amino acid composition showed that the frequency of amino acid residues in proteins of eukaryotic and bacterial origin is different for proteins with simple and complex folding kinetics.     

Cotranslational Protein Folding

The review analyzes the research concerning the folding of proteins in the course of their synthesis on ribosomes. The experimental data obtained for various proteins using various methods give grounds for concluding that a nascent protein largely acquires its spatial structure while still attached to the ribosome, and final folding into the biologically active conformation takes place as soon as the completed protein is released therefrom. Cotranslational folding is characteristic of both bacterial and eukaryotic cells, and appears to be the universal and the most evolutionarily ancient mechanism.     

Computer-based redesign of a protein folding pathway.

A fundamental test of our current understanding of protein folding is to rationally redesign protein folding pathways. We use a computer-based design strategy to switch the folding pathway of protein G, which normally involves formation of the second, but not the first, beta-turn at the rate limiting step in folding. Backbone conformations and amino acid sequences that maximize the interaction density in the first beta-hairpin were identified, and two variants containing 11 amino acid replacements were found to be approximately 4 kcal mol-1 more stable than wild type protein G. Kinetic studies show that the redesigned proteins fold approximately 100 x faster than wild type protein and that the first beta-turn is formed and the second disrupted at the rate limiting step in folding.     

Aggregation and assembly of phage P22 temperature-sensitive coat protein mutants in vitro mimic the in vivo phenotype

Aggregation is a common side reaction in the folding of proteins which is likely due to inappropriate interactions of folding intermediates. In the in vivo folding of phage P22 coat protein, amino acid substitutions that cause a temperature-sensitive-folding (tsf) phenotype lead to the localization of the mutant coat proteins to inclusion bodies. Investigated here is the aggregation of wild-type (WT) coat protein and 3 tsf mutants of coat protein. The tsf coat proteins aggregated when refolded in vitro at high temperature. If the tsf coat proteins were refolded at 4 degrees C, they were able attain an assembly active state. WT coat protein, on the other hand, did not aggregate significantly even when folded at high temperature. The refolded tsf mutants exhibited altered secondary and tertiary structures and had an increased surface hydrophobicity, which may explain the increased propensity of their folding intermediates to aggregate.     

Cotranslational protein folding

The review analyzes the research concerning the folding of proteins in the course of their synthesis on ribosomes. The experimental data obtained for various proteins using various methods give grounds for concluding that a nascent protein largely acquires its spatial structure while still attached to the ribosome, and final folding into the biologically active conformation takes place as soon as the completed protein is released therefrom. Cotranslational folding is characteristic of both bacterial and eukaryotic cells, and appears to be the universal and the most evolutionarily ancient mechanism.     

Rapid unfolding of a domain populates an aggregation-prone intermediate that can be recognized by GroEL

Some amino acid substitutions in phage P22 coat protein cause a temperature-sensitive folding (tsf) phenotype. In vivo, these tsf amino acid substitutions cause coat protein to aggregate and form intracellular inclusion bodies when folded at high temperatures, but at low temperatures the proteins fold properly. Here the effects of tsf amino acid substitutions on folding and unfolding kinetics and the stability of coat protein in vitro have been investigated to determine how the substitutions change the ability of coat protein to fold properly. The equilibrium unfolding transitions of the tsf variants were best fit to a three-state model, N if I if U, where all species concerned were monomeric, a result confirmed by velocity sedimentation analytical ultracentrifugation. The primary effect of the tsf amino acid substitutions on the equilibrium unfolding pathway was to decrease the stability (DeltaG) and the solvent accessibility (m-value) of the N if I transition. The kinetics of folding and unfolding of the tsf coat proteins were investigated using tryptophan fluorescence and circular dichroism (CD) at 222 nm. The tsf amino acid substitutions increased the rate of unfolding by 8-14-fold, with little effect on the rate of folding, when monitored by tryptophan fluorescence. In contrast, when folding or unfolding reactions were monitored by CD, the reactions were too fast to be observed. The tsf coat proteins are natural substrates for the molecular chaperones, GroEL/S. When native tsf coat protein monomers were incubated with GroEL, they bound efficiently, indicating that a folding intermediate was significantly populated even without denaturant. Thus, the tsf coat proteins aggregate in vivo because of an increased propensity to populate this unfolding intermediate.