Translation‐coupled protein folding assay using a protease to monitor the folding status

Protein folding is an essential prerequisite for proteins to execute nearly all cellular functions. There is a growing demand for a simple and robust method to investigate protein folding on a large‐scale under the same conditions. We previously developed a global folding assay system, in which proteins translated using an Escherichia coli‐based cell‐free translation system are centrifuged to quantitate the supernatant fractions. Although the assay is based on the assumption that the supernatants contain the folded native states, the supernatants also include nonnative unstructured proteins. In general, proteases recognize and degrade unstructured proteins, and thus we used a protease to digest the unstructured regions to monitor the folding status. The addition of Lon protease during the translation of proteins unmasked subfractions, not only in the soluble fractions but also in the aggregation‐prone fractions. We translated ～90 E. coli proteins in the protease‐inclusion assay, in the absence and presence of chaperones. The folding assay, which sheds light on the molecular mechanisms underlying the aggregate formation and the chaperone effects, can be applied to a large‐scale analysis.     

Hydrophobic folding units at protein-protein interfaces: implications to protein folding and to protein-protein association.

A hydrophobic folding unit cutting algorithm, originally developed for dissecting single-chain proteins, has been applied to a dataset of dissimilar two-chain protein-protein interfaces. Rather than consider each individual chain separately, the two-chain complex has been treated as a single chain. The two-chain parsing results presented in this work show hydrophobicity to be a critical attribute of two-state versus three-state protein-protein complexes. The hydrophobic folding units at the interfaces of two-state complexes suggest that the cooperative nature of the two-chain protein folding is the outcome of the hydrophobic effect, similar to its being the driving force in a single-chain folding. In analogy to the protein-folding process, the two-chain, two-state model complex may correspond to the formation of compact, hydrophobic nuclei. On the other hand, the three-state model complex involves binding of already folded monomers, similar to the association of the hydrophobic folding units within a single chain. The similarity between folding entities in protein cores and in two-state protein-protein interfaces, despite the absence of some chain connectivities in the latter, indicates that chain linkage does not necessarily affect the native conformation. This further substantiates the notion that tertiary, non-local interactions play a critical role in protein folding. These compact, hydrophobic, two-chain folding units, derived from structurally dissimilar protein-protein interfaces, provide a rich set of data useful in investigations of the role played by chain connectivity and by tertiary interactions in studies of binding and of folding. Since they are composed of non-contiguous pieces of protein backbones, they may also aid in defining folding nuclei.     

Influence of denatured and intermediate states of folding on protein aggregation

We simulate the aggregation thermodynamics and kinetics of proteins L and G, each of which self-assembles to the same alpha/beta [corrected] topology through distinct folding mechanisms. We find that the aggregation kinetics of both proteins at an experimentally relevant concentration exhibit both fast and slow aggregation pathways, although a greater proportion of protein G aggregation events are slow relative to those of found for protein L. These kinetic differences are correlated with the amount and distribution of intrachain contacts formed in the denatured state ensemble (DSE), or an intermediate state ensemble (ISE) if it exists, as well as the folding timescales of the two proteins. Protein G aggregates more slowly than protein L due to its rapidly formed folding intermediate, which exhibits native intrachain contacts spread across the protein, suggesting that certain early folding intermediates may be selected for by evolution due to their protective role against unwanted aggregation. Protein L shows only localized native structure in the DSE with timescales of folding that are commensurate with the aggregation timescale, leaving it vulnerable to domain swapping or nonnative interactions with other chains that increase the aggregation rate. Folding experiments that characterize the structural signatures of the DSE, ISE, or the transition state ensemble (TSE) under nonaggregating conditions should be able to predict regions where interchain contacts will be made in the aggregate, and to predict slower aggregation rates for proteins with contacts that are dispersed across the fold. Since proteins L and G can both form amyloid fibrils, this work also provides mechanistic and structural insight into the formation of prefibrillar species.     

Reinvestigating the codon and amino acid usage of S. cerevisiae genome: a new insight from protein secondary structure analysis

Biased usage of synonymous codons has been elucidated under the perspective of cellular tRNA abundance for quite a long time now. Taking advantage of publicly available gene expression data for Saccharomyces cerevisiae, a systematic analysis of the codon and amino acid usages in two different coding regions corresponding to the regular (helix and strand) as well as the irregular (coil) protein secondary structures, have been performed. Our analyses suggest that apart from tRNA abundance, mRNA folding stability is another major evolutionary force in shaping the codon and amino acid usage differences between the highly and lowly expressed genes in S. cerevisiae genome and surprisingly it depends on the coding regions corresponding to the secondary structures of the encoded proteins. This is obviously a new paradigm in understanding the codon usage in S. cerevisiae. Differential amino acid usage between highly and lowly expressed genes in the regions coding for the irregular protein secondary structure in S. cerevisiae is expounded by the stability of the mRNA folded structure. Irrespective of the protein secondary structural type, the highly expressed genes always tend to encode cheaper amino acids in order to reduce the overall biosynthetic cost of production of the corresponding protein. This study supports the hypothesis that the tRNA abundance is a consequence of and not a reason for the biased usage of amino acid between highly and lowly expressed genes.     

Surprisingly high correlation between early and late stages in non-two-state protein folding

We collected quantitative kinetic data on early and late stages of folding in non-two-state proteins from the literature, and studied the relationship between the kinetics of the two stages. There was a surprisingly high correlation between the rate constants of these stages. The correlation coefficient of the logarithmic rate constants was as high as 0.97, which could not be caused by chance. We also studied relationships of the logarithmic rate constants of the two stages with native three-dimensional structures represented by the residue-residue contact map. There were again surprisingly high correlations between the logarithmic rate constants and the number of non-local contact clusters obtained from the contact maps. Because the number of non-local contact clusters represents overall arrangement of substructures in a native protein, the results strongly suggested the importance of the arrangement of the substructures for the kinetics of both early and late stages of protein folding.     

The relation between mRNA folding and protein structure

About 200 mRNA sequences of Escherichia coli and human with matching protein secondary structure data were studied. The mRNA folding for each native sequence and for corresponding randomized sequences was calculated through free energy minimization. We have found that the folding energy of mRNA segments in different protein secondary structures is significantly different. The average Z score is more negative for regular secondary structure (alpha-helix and beta-strand) than that for coil. This suggests that the codon choice in native mRNA sequence coding for protein regular structure contributes more to the mRNA folding stability.     

Molecular basis of cooperativity in protein folding IV. Core: A general cooperative folding model

The cooperative nature of the protein folding process is independent of the characteristic fold and the specific secondary structure attributes of a globular protein. A general folding/unfolding model should, therefore, be based upon structural features that transcend the peculiarities of α-helices, β-sheets, and other structural motifs found in proteins. The studies presented in this paper suggest that a single structural characteristic common to all globular proteins is essential for cooperative folding. The formation of a partly folded state from the native state results in the exposure to solvent of two distinct regions: (1) the portions of the protein that are unfolded; and (2) the “complementary surfaces,” located in the regions of the protein that remain folded. The cooperative character of the folding/unfolding transition is determined largely by the energetics of exposing complementary surface regions to the solvent. By definition, complementary regions are present only in partly folded states; they are absent from the native and unfolded states. An unfavorable free energy lowers the probability of partly folded states and increases the cooperativity of the transition. In this paper we present a mathematical formulation of this behavior and develop a general cooperative folding/unfolding model, termed the “complementary region” (CORE) model. This model successfully reproduces the main properties of folding/unfolding transitions without limiting the number of partly folded states accessible to the protein, thereby permitting a systematic examination of the structural and solvent conditions under which intermediates become populated. It is shown that the CORE model predicts two-state folding/unfolding behavior, even though the two-state character is not assumed in the model. © 1993 Wiley-Liss, Inc.     

Protein folding following synthesis in vitro and in vivo: association of newly synthesized protein with 50S subunit of E. coli ribosome

In the accompanying paper, it was shown that a protein, while reverting to native form from the unfolded state in vitro with the help of bacterial 70S ribosome, split the latter into its subunits (50S and 30S) and remains associated with the 50S subunit. Here, we follow the fate of nascent proteins both in case of in vivo and in vitro translation system. The newly synthesised protein was found to associate with the 50S subunit in both the cases.     

Physical principles of protein structure and protein folding

Physical principles determining the protein structure and protein folding are reviewed: (i) the molecular theory of protein secondary structure and the method of its prediction based on this theory; (ii) the existence of a limited set of thermodynamically favourable folding patterns of α- and β-regions in a compact globule which does not depend on the details of the amino acid sequence; (iii) the moderns approaches to the prediction of the folding patterns of α- and β-regions in concrete proteins; (iv) experimental approaches to the mechanism of protein folding. The review reflects theoretical and experimental works of the author and his collaborators as well as those of other groups.     

Polymer principles and protein folding.

This paper surveys the emerging role of statistical mechanics and polymer theory in protein folding. In the polymer perspective, the folding code is more a solvation code than a code of local phipsi propensities. The polymer perspective resolves two classic puzzles: (1) the Blind Watchmaker's Paradox that biological proteins could not have originated from random sequences, and (2) Levinthal's Paradox that the folded state of a protein cannot be found by random search. Both paradoxes are traditionally framed in terms of random unguided searches through vast spaces, and vastness is equated with impossibility. But both processes are partly guided. The searches are more akin to balls rolling down funnels than balls rolling aimlessly on flat surfaces. In both cases, the vastness of the search is largely irrelevant to the search time and success. These ideas are captured by energy and fitness landscapes. Energy landscapes give a language for bridging between microscopics and macroscopics, for relating folding kinetics to equilibrium fluctuations, and for developing new and faster computational search strategies.     

Effect of denaturant and protein concentrations upon protein refolding and aggregation: a simple lattice model.

We present a study of the competition between protein refolding and aggregation for simple lattice model proteins. The effect of solvent conditions (i.e., the denaturant concentration and the protein concentration) on the folding and aggregation behavior of a system of simple, two-dimensional lattice protein molecules has been investigated via (dynamic Monte Carlo simulations. The population profiles and aggregation propensities of the nine most populated intermediate configurations exhibit a complex dependence on the solution conditions that can be understood by considering the competition between intra- and interchain interactions. Some of these configurations are not even seen in isolated chain simulations; they are observed to be highly aggregation prone and are stabilized primarily by the aggregation reaction in multiple-chain systems. Aggregation arises from the association of partially folded intermediates rather than from the association of denatured random-coil states. The aggregation reaction dominates over the folding reaction at high protein concentration and low denaturant concentration, resulting in low refolding yields at those conditions. However, optimum folding conditions exist at which the refolding yield is a maximum, in agreement with some experimental observations.     

Side-chain conformational entropy in protein folding.

An important, but often neglected, contribution to the thermodynamics of protein folding is the loss of entropy that results from restricting the number of accessible side-chain conformers in the native structure. Conformational entropy changes can be found by comparing the number of accessible rotamers in the unfolded and folded states, or by estimating fusion entropies. Comparison of several sets of results using different techniques shows that the mean conformational free energy change (T delta S) is 1 kcal.mol-1 per side chain or 0.5 kcal.mol-1 per bond. Changes in vibrational entropy appear to be negligible compared to the entropy change resulting from the loss of accessible rotamers. Side-chain entropies can help rationalize alpha-helix propensities, predict protein/inhibitor complex structures, and account for the distribution of side chains on the protein surface or interior.     

Causal signals between codon bias,mRNA structure,and the efficiency of translation and elongation

Ribosome profiling data report on the distribution of translating ribosomes, at steady‐state, with codon‐level resolution. We present a robust method to extract codon translation rates and protein synthesis rates from these data, and identify causal features associated with elongation and translation efficiency in physiological conditions in yeast. We show that neither elongation rate nor translational efficiency is improved by experimental manipulation of the abundance or body sequence of the rare AGG tRNA. Deletion of three of the four copies of the heavily used ACA tRNA shows a modest efficiency decrease that could be explained by other rate‐reducing signals at gene start. This suggests that correlation between codon bias and efficiency arises as selection for codons to utilize translation machinery efficiently in highly translated genes. We also show a correlation between efficiency and RNA structure calculated both computationally and from recent structure probing data, as well as the Kozak initiation motif, which may comprise a mechanism to regulate initiation.     

Regulation of phytochrome mRNA abundance in green oat leaves

Abstract. Avena sativa L. (oat) seedings were grown 4 d in continuous white light followed by 3 d in darkness. Probes derived from an oat phytochrome cDNA clone (pAP 3.2) were used in slot blot analyses to measure the abundance of phytochrome mRNA in the distinct etiolated and green portions of the leaves produced by these seedlings. Both the green and etiolated portions accumulated phytochrome mRNA to a level of about 85% of the etiolated seedling level. Subsequent experiments with similar seedlings showed that both the green and etiolated portions were capable of inducing a dramatic decline in phytochrome mRNA abundance in response to a saturating red light pulse. Despite the ability of green portions of oat leaves to accumulate phytochrome mRNA and to down-regulate phytochrome mRNA abundance in response to light, no substantial variation in phytochrome mRNA abundance was observed in green oat seedlings maintained on a 12-h day/12-h night cycle. In the same oat seedlings, the abundance of chlorophyll a/b binding protein mRNA fluctuated dramatically during the day/night cycle.     

Charge centers and formation of the protein folding core

Electrostatic interactions are important for protein folding. At low resolution, the electrostatic field of the whole molecule can be described in terms of charge center(s). To study electrostatic effects, the centers of positive and negative charge were calculated for 20 small proteins of known structure, for which hydrogen exchange cores had been determined experimentally. Two observations seem to be important. First, in all 20 proteins studied 30-100% of the residues forming hydrogen exchange core(s) were clustered around the charge centers. Moreover, in each protein more than half of the core sequences are located near the centers of charge. Second, the general architecture of all-alpha proteins from the set seems to be stabilized by interactions of residues surrounding the charge centers. In most of the alpha-beta proteins, either or both of the centers are located near a pair of consecutive strands, and this is even more characteristic for alpha/Beta and all-beta structures. Consecutive strands are very probable sites of early folding events. These two points lead to the conclusion that charge centers, defined solely from the structure of the folded protein may indicate the location of a protein's hydrogen exchange/folding core. In addition, almost all the proteins contain well-conserved continuous hydrophobic sequences of three or more residues located in the vicinity of the charge centers. These hydrophobic sequences may be primary nucleation sites for protein folding. The results suggest the following scheme for the order of events in folding: local hydrophobic nucleation, electrostatic collapse of the core, global hydrophobic collapse, and slow annealing to the native state. This analysis emphasizes the importance of treating electrostatics during protein-folding simulations.     

Kinetic barriers and the role of topology in protein and RNA folding

This review compares the folding behavior of proteins and RNAs. Topics covered include the role of topology in the determination of folding rates, major folding events including collapse, properties of denatured states, pathway heterogeneity, and the influence of the mode of initiation on the folding pathway.     

Reducing the computational complexity of protein folding via fragment folding and assembly

Understanding, and ultimately predicting, how a 1-D protein chain reaches its native 3-D fold has been one of the most challenging problems during the last few decades. Data increasingly indicate that protein folding is a hierarchical process. Hence, the question arises as to whether we can use the hierarchical concept to reduce the practically intractable computational times. For such a scheme to work, the first step is to cut the protein sequence into fragments that form local minima on the polypeptide chain. The conformations of such fragments in solution are likely to be similar to those when the fragments are embedded in the native fold, although alternate conformations may be favored during the mutual stabilization in the combinatorial assembly process. Two elements are needed for such cutting: (1) a library of (clustered) fragments derived from known protein structures and (2) an assignment algorithm that selects optimal combinations to "cover" the protein sequence. The next two steps in hierarchical folding schemes, not addressed here, are the combinatorial assembly of the fragments and finally, optimization of the obtained conformations. Here, we address the first step in a hierarchical protein-folding scheme. The input is a target protein sequence and a library of fragments created by clustering building blocks that were generated by cutting all protein structures. The output is a set of cutout fragments. We briefly outline a graph theoretic algorithm that automatically assigns building blocks to the target sequence, and we describe a sample of the results we have obtained.     

Water as a conformational editor in protein folding

As molecules approach one another in aqueous solution, desolvation free energy barriers to association are encountered. Experiments suggest these (de)solvation effects contribute to the free energy barriers separating the folded and unfolded states of protein molecules. To explore their influence on the energy landscapes of protein folding reactions, we have incorporated desolvation barriers into a semi-realistic, off-lattice protein model that uses a simplified physico-chemical force-field determined solely by the sequence of amino acids. Monte Carlo sampling techniques were used to study the effects on the thermodynamics and kinetics of folding of a number of systems, diverse in structure and sequence. In each case, desolvation barriers increase the stability of the native conformation and the cooperativity of the major folding/unfolding transition. The folding times of these systems are reduced significantly upon inclusion of desolvation barriers, demonstrating that the particulate nature of the solvent engenders a more defined route to the native fold.     

The competition between protein folding and aggregation: off-lattice minimalist model studies

Protein aggregation has been associated with a number of human diseases, and is a serious problem in the manufacture of recombinant proteins. Of particular interest to the biotechnology industry is deleterious aggregation that occurs during the refolding of proteins from inclusion bodies. As a complement to experimental efforts, computer simulations of multi-chain systems have emerged as a powerful tool to investigate the competition between folding and aggregation. Here we report results from Langevin dynamics simulations of minimalist model proteins. Order parameters are developed to follow both folding and aggregation. By mapping natural units to real units, the simulations are shown to be carried out under experimentally relevant conditions. Data pertaining to the contacts formed during the association process show that multiple mechanisms for aggregation exist, but certain pathways are statistically preferred. Kinetic data show that there are multiple time scales for aggregation, although most association events take place at times much shorter than those required for folding. Last, we discuss results presented here as a basis for future work aimed at rational design of mutations to reduce aggregation propensity, as well as for development of small-molecular weight refolding enhancers.