Multiple alanine replacements within alpha-helix 126-134 of T4 lysozyme have independent, additive effects on both structure and stability.

In a systematic attempt to identify residues important in the folding and stability of T4 lysozyme, five amino acids within alpha-helix 126-134 were substituted by alanine, either singly or in selected combinations. Together with three alanines already present in the wild-type structure this provided a set of mutant proteins with up to eight alanines in sequence. All the variants behaved normally, suggesting that the majority of residues in the alpha-helix are nonessential for the folding of T4 lysozyme. Of the five individual alanine substitutions it is inferred that four result in slightly increased protein stability and one, the replacement of a buried leucine with alanine, substantially decreased stability. The results support the idea that alanine is a residue of high helix propensity. The change in protein stability observed for each of the multiple mutants is approximately equal to the sum of the energies associated with each of the constituent substitutions. All of the variants could be crystallized isomorphously with wild-type lysozyme, and, with one trivial exception, their structures were determined at high resolution. Substitution of the largely solvent-exposed residues Asp 127, Glu 128, and Val 131 with alanine caused essentially no change in structure except at the immediate site of replacement. Substitutions of the partially buried Asn 132 and the buried Leu 133 with alanine were associated with modest (< or = 0.4 A) structural adjustments. The structural changes seen in the multiple mutants were essentially a combination of those seen in the constituent single replacements. The different replacements therefore act essentially independently not only so far as changes in energy are concerned but also in their effect on structure. The destabilizing replacement Leu 133-->Ala made alpha-helix 126-134 somewhat less regular. Incorporation of additional alanine replacements tended to make the helix more uniform. For the penta-alanine variant a distinct change occurred in a crystal-packing contact, and the "hinge-bending angle" between the amino- and carboxy-terminal domains changed by 3.6 degrees. This tends to confirm that such hinge-bending in T4 lysozyme is a low-energy conformational change.     

Conformational stability and folding mechanisms of dimeric proteins

The folding of multisubunit proteins is of tremendous biological significance since the large majority of proteins exist as protein-protein complexes. Extensive experimental and computational studies have provided fundamental insights into the principles of folding of small monomeric proteins. Recently, important advances have been made in extending folding studies to multisubunit proteins, in particular homodimeric proteins. This review summarizes the equilibrium and kinetic theory and models underlying the quantitative analysis of dimeric protein folding using chemical denaturation, as well as the experimental results that have been obtained. Although various principles identified for monomer folding also apply to the folding of dimeric proteins, the effects of subunit association can manifest in complex ways, and are frequently overlooked. Changes in molecularity typically give rise to very different overall folding behaviour than is observed for monomeric proteins. The results obtained for dimers have provided key insights pertinent to understanding biological assembly and regulation of multisubunit proteins. These advances have set the stage for future advances in folding involving protein-protein interactions for natural multisubunit proteins and unnatural assemblies involved in disease.     

Probing protein folding and stability using disulfide bonds

Disulfide bonds are required to stabilize the folded conformations of many proteins. The rates and equilibria of processes involved in disulfide bond formation and breakage can be manipulated experimentally and can be used to obtain important information about protein folding and stability. A number of experimental procedures for studying these processes, and approaches to interpreting the resulting data, are described here.     

Subdomain interactions as a determinant in the folding and stability of T4 lysozyme.

The folding of large, multidomain proteins involves the hierarchical assembly of individual domains. It remains unclear whether the stability and folding of small, single-domain proteins occurs through a comparable assembly of small, autonomous folding units. We have investigated the relationship between two subdomains of the protein T4 lysozyme. Thermodynamically, T4 lysozyme behaves as a cooperative unit and the unfolding transition fits a two-state model. The structure of the protein, however, resembles a dumbbell with two potential subdomains: an N-terminal subdomain (residues 13-75), and a C-terminal subdomain (residues 76-164 and 1-12). To investigate the effect of uncoupling these two subdomains within the context of the native protein, we created two circular permutations, both at the subdomain interface (residues 13 and 75). Both variants adopt an active wild-type T4 lysozyme fold. The protein starting with residue 13 is 3 kcal/mol less stable than wild type, whereas the protein beginning at residue 75 is 9 kcal/mol less stable, suggesting that the placement of the termini has a major effect on protein stability while minimally affecting the fold. When isolated as protein fragments, the C-terminal subdomain folds into a marginally stable helical structure, whereas the N-terminal subdomain is predominantly unfolded. ANS fluorescence studies indicate that, at low pH, the C-terminal subdomain adopts a loosely packed acid state. An acid state intermediate is also seen for all of the full-length variants. We propose that this acid state is comprised of an unfolded N-terminal subdomain and a loosely folded C-terminal subdomain.     

The folding pathway of a functionally competent C-terminal domain of nucleophosmin: Protein stability and denatured state residual structure

Nucleophosmin (NPM1) is a nucleolar protein implicated in ribosome biogenesis, centrosome duplication and cell cycle control; the NPM1 gene is the most frequent target for mutations in Acute Myeloid Leukemia. Mutations map to the C-terminal domain of the protein and cause its unfolding, loss of DNA binding properties and aberrant cellular localization. Here we investigate the folding pathway and denatured state properties of a NPM1 C-terminal domain construct encompassing the last 70 residues in the reference sequence. This construct is more stable than the previously characterized domain, which consisted of the last 53 residues. Data reveal that, similarly to what was discovered for the shorter construct, also the 70-residue construct of NPM1 displays a detectable residual structure in its denatured state. The higher stability of the latter domain allows us to conclude that the denatured state is robust to changes in solvent composition and that it consists of a discrete state in equilibrium with the expanded fully unfolded conformation. This observation, which might appear as a technicality, is in fact of general importance for the understanding of the folding of proteins. The implications of our results are discussed in the context of previous works on single domain helical proteins.     

Effect of surface electrostatic interactions on the stability and folding of formate dehydrogenase from Candida methylica

NAD⁺-dependent formate dehydrogenase (FDH-EC 1.2.1.2) is an important enzyme to regenerate valuable NADH required by NAD⁺-dependent oxidoreductases in enzyme catalysis. The limitation in the thermostability of FDH enzyme is a crucial problem for development of biotechnological and industrial processes, despite of its advantages. In this study, to investigate the contribution of surface electrostatic interaction to the thermostability of FDH from Candida methylica (cmFDH) N187E, H13E, Q105R, N300E, N147R N300E/N147R, N187E/Q105R, N187E/N147R,Y160R, Y302R, Y160E and Y302E mutants were designed using a homology model of cmFDH based on Candida boidinii (cb) by considering electrostatic interactions on the protein surface. The effects of site-specific engineering on the stability of this molecule was analyzed according to minimal model of folding and assembly reaction and deduced equilibrium properties of the native system with respect to its thermal and denaturant sensitivities. It was observed that mutations did not change the unfolding pattern of native cmFDH and increased numbers of electrostatic interactions can cause either stabilizing or destabilizing effect on the thermostability of this protein. The thermodynamic and kinetic results suggested that except relatively improved mutants, three out of the nine single mutations increased the melting temperature of cmFDH enzyme.     

Equilibrium and kinetic studies on folding of canine milk lysozyme

Thermal and chemical unfolding studies of the calcium-binding canine lysozyme (CL) by fluorescence and circular dichroism spectroscopy show that, upon unfolding in the absence of calcium ions, a very stable equilibrium intermediate state is formed. At room temperature and pH 7.5, for example, a stable molten globule state is attained in 3 M GdnHCl. The existence of such a pure and stable intermediate state allowed us to extend classical stopped-flow fluorescence measurements that describe the transition from the native to the unfolded form, with kinetic experiments that monitor separately the transition from the unfolded to the intermediate state and from the intermediate to the native state, respectively. The overall refolding kinetics of apo-canine lysozyme are characterized by a significant drop in the fluorescence intensity during the dead time, followed by a monoexponential increase of the fluorescence with k = 3.6 s(-1). Furthermore, the results show that, unlike its drastic effect on the stability, Ca(2+)-binding only marginally affects the refolding kinetics. During the refolding process of apo-CL non-native interactions, comparable to those observed in hen egg white lysozyme, are revealed by a substantial quenching of tryptophan fluorescence. The dissection of the refolding process in two distinct steps shows that these non-native interactions only occur in the final stage of the refolding process in which the two domains match to form the native conformation.     

Perturbations of the denatured state ensemble: modeling their effects on protein stability and folding kinetics.

By considering the denatured state of a protein as an ensemble of conformations with varying numbers of sequence-specific interactions, the effects on stability, folding kinetics, and aggregation of perturbing these interactions can be predicted from changes in the molecular partition function. From general considerations, the following conclusions are drawn: (1) A perturbation that enhances a native interaction in denatured state conformations always increases the stability of the native state. (2) A perturbation that promotes a non-native interaction in the denatured state always decreases the stability of the native state. (3) A change in the denatured state ensemble can alter the kinetics of aggregation and folding. (4) The loss (or increase) in stability accompanying two mutations, each of which lowers (or raises) the free energy of the denatured state, will be less than the sum of the effects of the single mutations, except in cases where both mutations affect the same set of partially folded conformations. By modeling the denatured state as the ensemble of all non-native conformations of hydrophobic-polar (HP) chains configured on a square lattice, it can be shown that the stabilization obtained from enhancement of native interactions derives in large measure from the avoidance of non-native interactions in the D state. In addition, the kinetic effects of fixing single native contacts in the denatured state or imposing linear gradients in the HH contact probabilities are found, for some sequences, to significantly enhance the efficiency of folding by a simple hydrophobic zippering algorithm. Again, the dominant mechanism appears to be avoidance of non-native interactions. These results suggest stabilization of native interactions and imposition of gradients in the stability of local structure are two plausible mechanisms involving the denatured state that could play a role in the evolution of protein folding and stability.     

Effects of macromolecular crowding agents on protein folding in vitro and in silico

Proteins fold and function inside cells which are environments very different from that of dilute buffer solutions most often used in traditional experiments. The crowded milieu results in excluded-volume effects, increased bulk viscosity and amplified chances for inter-molecular interactions. These environmental factors have not been accounted for in most mechanistic studies of protein folding executed during the last decades. The question thus arises as to how these effects—present when polypeptides normally fold in vivo—modulate protein biophysics. To address excluded volume effects, we use synthetic macromolecular crowding agents, which take up significant volume but do not interact with proteins, in combination with strategically selected proteins and a range of equilibrium and time-resolved biophysical (spectroscopic and computational) methods. In this review, we describe key observations on macromolecular crowding effects on protein stability, folding and structure drawn from combined in vitro and in silico studies. As expected based on Minton’s early predictions, many proteins (apoflavodoxin, VlsE, cytochrome c, and S16) became more thermodynamically stable (magnitude depends inversely on protein stability in buffer) and, unexpectedly, for apoflavodoxin and VlsE, the folded states changed both secondary structure content and, for VlsE, overall shape in the presence of macromolecular crowding. For apoflavodoxin and cytochrome c, which have complex kinetic folding mechanisms, excluded volume effects made the folding energy landscapes smoother (i.e., less misfolding and/or kinetic heterogeneity) than in buffer.     

Mutational analysis of the BPTI folding pathway: II. Effects of aromatic-->leucine substitutions on folding kinetics and thermodynamics.

The rates of the individual steps in the disulfide-coupled folding and unfolding of eight BPTI variants, each containing a single aromatic to leucine amino acid replacement, were measured. From this analysis, the contributions of the four phenylalanine and four tyrosine residues to the stabilities of the native protein and the disulfide-bonded folding intermediates were determined. While the substitutions were found to destabilize the native protein by 2 to 7 kcal/mol, they had significantly smaller effects on the intermediates that represent the earlier stages of folding, even when the site of the substitution was located within the ordered regions of the intermediates. These results suggest that stabilizing interactions contribute less to conformational stability in the context of a partially folded intermediate than in a fully folded native protein, perhaps because of decreased cooperativity among the individual interactions. The kinetic analysis also provides new information about the transition states associated with the slowest steps in folding and unfolding, supporting previous suggestions that these transition states are extensively unfolded. Although the substitutions caused large changes in the distribution of folding intermediates and in the rates of some steps in the folding pathway, the kinetically-preferred pathway for all of the variants involved intramolecular disulfide rearrangements, as observed previously for the wild-type protein. These results suggest that the predominance of the rearrangement mechanism reflects conformational constraints present relatively early in the folding pathway.     

The response of T4 lysozyme to large-to-small substitutions within the core and its relation to the hydrophobic effect.

To further examine the structural and thermodynamic basis of hydrophobic stabilization in proteins, all of the bulky non-polar residues that are buried or largely buried within the core of T4 lysozyme were substituted with alanine. In 25 cases, including eight reported previously, it was possible to determine the crystal structures of the variants. The structures of four variants with double substitutions were also determined. In the majority of cases the "large-to-small" substitutions lead to internal cavities. In other cases declivities or channels open to the surface were formed. In some cases the structural changes were minimal (mainchain shifts < or = 0.3 A); in other cases mainchain atoms moved up to 2 A. In the case of Ile 29 --> Ala the structure collapsed to such a degree that the volume of the putative cavity was zero. Crystallographic analysis suggests that the occupancy of the engineered cavities by solvent is usually low. The mutants Val 149 --> Ala (V149A) and Met 6 --> Ala (M6A), however, are exceptions and have, respectively, one and two well-ordered water molecules within the cavity. The Val 149 --> Ala substitution allows the solvent molecule to hydrogen bond to polar atoms that are occluded in the wild-type molecule. Similarly, the replacement of Met 6 with alanine allows the two solvent molecules to hydrogen bond to each other and to polar atoms on the protein. Except for Val 149 --> Ala the loss of stability of all the cavity mutants can be rationalized as a combination of two terms. The first is a constant for a given class of substitution (e.g., -2.1 kcal/mol for all Leu --> Ala substitutions) and can be considered as the difference between the free energy of transfer of leucine and alanine from solvent to the core of the protein. The second term can be considered as the energy cost of forming the cavity and is consistent with a numerical value of 22 cal mol(-1) A(-3). Physically, this term is due to the loss of van der Waal''s interactions between the bulky sidechain that is removed and the atoms that form the wall of the cavity. The overall results are consistent with the prior rationalization of Leu --> Ala mutants in T4 lysozyme by Eriksson et al. (Eriksson et al., 1992, Science 255:178-183).     

Electrostatic interactions in the denatured state ensemble: their effect upon protein folding and protein stability

It is now recognized that the denatured state ensemble (DSE) of proteins can contain significant amounts of structure, particularly under native conditions. Well-studied examples include small units of hydrogen bonded secondary structure, particularly helices or turns as well as hydrophobic clusters. Other types of interactions are less well characterized and it has often been assumed that electrostatic interactions play at most a minor role in the DSE. However, recent studies have shown that both favorable and unfavorable electrostatic interactions can be formed in the DSE. These can include surprisingly specific non-native interactions that can even persist in the transition state for protein folding. DSE electrostatic interactions can be energetically significant and their modulation either by mutation or by varying solution conditions can have a major impact upon protein stability. pH dependent stability studies have shown that electrostatic interactions can contribute up to 4 kcal mol^-1 to the stability of the DSE.     

Flash dietary methionine supply over growth requirements in pigs: Multi-facetted effects on skeletal muscle metabolism

Dietary methionine affects protein metabolism, lean gain and growth performance and acts in the control of oxidative stress. When supplied in large excess relative to growth requirements in diets for pigs, positive effects on pork quality traits have been recently reported. This study aimed to decipher the molecular and biochemical mechanisms affected by a dietary methionine supply above growth requirements in the loin muscle of finishing pigs. During the last 14 days before slaughter, crossbred female pigs (n = 15 pigs/diet) were fed a diet supplemented with hydroxy-methionine (Met5; 1.1% of methionine) or not (CONT, 0.22% of methionine). Blood was sampled at slaughter to assess key metabolites. At the same time, free amino acid concentrations and expression or activity levels of genes involved in protein or energy metabolism were measured in the longissimus lumborum muscle (LM). The Met5 pigs exhibited a greater activity of creatine kinase in plasma when compared with CONT pigs. The concentrations of free methionine, alpha-aminobutyric acid, anserine, 3-methyl-histidine, lysine, and proline were greater in the LM of Met5 pigs than in CONT pigs. Expression levels of genes involved in protein synthesis, protein breakdown or autophagy were only scarcely affected by the diet. Among ubiquitin ligases, MURF1, a gene known to target creatine kinase and muscle contractile proteins, and OTUD1 coding for a deubiquitinase protease, were up-regulated in the LM of Met5 pigs. A lower activity of citrate synthase, a reduced expression level of ME1 acting in lipogenesis but a higher expression of PPARD regulating energy metabolism, were also observed in the LM of Met5 pigs compared with CONT pigs. Principal component analysis revealed that expression levels of many studied genes involved in protein and energy metabolism were correlated with meat quality traits across dietary treatments, suggesting that subtle modifications in expression of those genes had cumulative effects on the regulation of processes leading to the muscle transformation into meat. In conclusion, dietary methionine supplementation beyond nutritional requirements in pigs during the last days before slaughter modified the free amino acid profile in muscle and its redox capacities, and slightly affected molecular pathways related to protein breakdown and energy metabolism. These modifications were associated with benefits on pork quality traits.     

Exploring subdomain cooperativity in T4 lysozyme II: uncovering the C-terminal subdomain as a hidden intermediate in the kinetic folding pathway

Intermediates along a protein's folding pathway can play an important role in its biology. Previous kinetics studies have revealed an early folding intermediate for T4 lysozyme, a small, well-characterized protein composed of an N-terminal and a C-terminal subdomain. Pulse-labeling hydrogen exchange studies suggest that residues from both subdomains contribute to the structure of this intermediate. On the other hand, equilibrium native state hydrogen experiments have revealed a high-energy, partially unfolded form of the protein that has an unstructured N-terminal subdomain and a structured C-terminal subdomain. To resolve this discrepancy between kinetics and equilibrium data, we performed detailed kinetics analyses of the folding and unfolding pathways of T4 lysozyme, as well as several point mutants and large-scale variants. The data support the argument for the presence of two distinct intermediates, one present on each side of the rate-limiting transition state barrier. The effects of circular permutation and site-specific mutations in the wild-type and circular permutant background, as well as a fragment containing just the C-terminal subdomain, support a model for the unfolding intermediate with an unfolded N-terminal and a folded C-terminal subdomain. Our results suggest that the partially unfolded form identified by native state hydrogen exchange resides on the folded side of the rate-limiting transition state and is, therefore, under most conditions, a "hidden" intermediate.     

Role of copper in folding and stability of cupredoxin-like copper-carrier protein CopC

CopC is a periplasmic copper carrier that, in contrast to cytoplasmic copper chaperones, has a beta-barrel fold and two metal-binding sites distinct for Cu(II) and Cu(I). The copper sites are located in each end of the molecule: the Cu(I) site involves His and Met coordination whereas the Cu(II) site consists of charged residues. To reveal biophysical properties of this protein, we have explored the effects of the cofactors on CopC unfolding in vitro. We demonstrate that Cu(II) coordination affects both protein stability and unfolding pathway, whereas Cu(I) has only a small effect on stability. Apo-CopC unfolds in a two-state reaction between pH 4 and 7.5 with maximal stability at pH 6. In contrast, Cu(II)-CopC unfolds in a three-state reaction at pH6 that involves a partly folded intermediate that retains Cu(II). This intermediate exhibits high thermal and chemical stability. Unique energetic and structural properties of different metalated CopC forms may help facilitate metal transport to many partners in vivo.     

Inter-residue interactions in protein folding and stability

During the process of protein folding, the amino acid residues along the polypeptide chain interact with each other in a cooperative manner to form the stable native structure. The knowledge about inter-residue interactions in protein structures is very helpful to understand the mechanism of protein folding and stability. In this review, we introduce the classification of inter-residue interactions into short, medium and long range based on a simple geometric approach. The features of these interactions in different structural classes of globular and membrane proteins, and in various folds have been delineated. The development of contact potentials and the application of inter-residue contacts for predicting the structural class and secondary structures of globular proteins, solvent accessibility, fold recognition and ab initio tertiary structure prediction have been evaluated. Further, the relationship between inter-residue contacts and protein-folding rates has been highlighted. Moreover, the importance of inter-residue interactions in protein-folding kinetics and for understanding the stability of proteins has been discussed. In essence, the information gained from the studies on inter-residue interactions provides valuable insights for understanding protein folding and de novo protein design.     

A split-ubiquitin-based assay detects the influence of mutations on the conformational stability of the p53 DNA binding domain in vivo

Many mutations in the human tumor suppressor p53 affect the function of the protein by destabilizing the structure of its DNA binding domain. To monitor the effects of those mutations in vivo the stability of the DNA binding domain of p53 and some of its mutants was investigated with the split-ubiquitin (split-Ub) method. The split-Ub-derived in vivo data on the relative stability of the mutants roughly correlate with the quantitative data from in vitro denaturation experiments as reported in the literature. A variation of this assay allows visualizing the difference in stability between the wild-type p53 core and the mutant p53(V143A) by a simple growth assay.     

Analysis of the effectiveness of proline substitutions and glycine replacements in increasing the stability of phage T4 lysozyme.

It was previously shown that the two replacements Gly 77-->Ala (G77A) and Ala 82-->Pro (A82P) increase the thermostability of phage T4 lysozyme at pH 6.5. Such replacements are presumed to restrict the degrees of freedom of the unfolded protein and so decrease the entropy of unfolding [B. W. Matthews, H. Nicholson, and W. J. Becktel (1987) Proceedings of the National Academy of Science USA Vol. 84, pp. 6663-6667]. To further test this approach, three additional replacements--G113A, K60P and A93P--have been constructed. On the basis of model building, each of these three replacements was judged to be less than optimal because it would tend to introduce unfavorable van der Waals contacts with neighboring parts of the protein. The presence of such contacts was verified for G113A and K60P by conformational adjustments seen in the crystal structures of these mutant proteins. In the case of G113A there are backbone conformational changes of 0.5-1.0 A in the short alpha-helix, 108-113, that includes the site of substitution. In the case of K60P the pyrrolidine ring shows evidence of strain. The thermal stability of each of the three variants at both pH 2.0 and pH 6.5 was found to be very close to that of wild-type lysozyme. The results suggest that the procedure used to predict sites for both Xaa-->Pro and Gly-->Ala is, in principle, correct. At the same time, the increase in stability expected from substitutions of this type is modest, and can easily be offset by strain associated with introduction of the alanine or proline. This means that the criteria used to select substitutions that will increase thermostability have to be stringent at least. In the case of T4 lysozyme this severely limits the number of sites. The analysis reveals a significant discrepancy between the conformational energy surface predicted for the residue preceding a proline and the conformations observed in crystal structures.     

The introduction of strain and its effects on the structure and stability of T4 lysozyme

In order to try to better understand the role played by strain in the structure and stability of a protein a series of "small-to-large" mutations was made within the core of T4 lysozyme. Three different alanine residues, one involved in backbone contacts, one in side-chain contacts, and the third adjacent to a small cavity, were each replaced with subsets of the larger residues, Val, Leu, Ile, Met, Phe and Trp. As expected, the protein is progressively destabilized as the size of the introduced side-chain becomes larger. There does, however, seem to be a limit to the destabilization, suggesting that a protein of a given size may be capable of maintaining only a certain amount of strain. The changes in stability vary greatly from site to site. Substitution of larger residues for both Ala42 and Ala98 substantially destabilize the protein, even though the primary contacts in one case are predominantly with side-chain atoms and in the other with backbone. The results suggest that it is neither practical nor meaningful to try to separate the effects of introduced strain on side-chains from the effects on the backbone. Substitutions at Ala129 are much less destabilizing than at sites 42 or 98. This is most easily understood in terms of the pre-existing cavity, which provides partial space to accommodate the introduced side-chains. Crystal structures were obtained for a number of the mutants. These show that the changes in structure to accommodate the introduced side-chains usually consist of essentially rigid-body displacements of groups of linked atoms, achieved through relatively small changes in torsion angles. On rare occasions, a side-chain close to the site of substitution may change to a different rotamer. When such rotomer changes occur, they permit the structure to dissipate strain by a response that is plastic rather than elastic. In one case, a surface loop moves 1.2 A, not in direct response to a mutation, but in an interaction mediated via an intermolecular contact. It illustrates how the structure of a protein can be modified by crystal contacts.