Optimizing the stability of single-chain proteins by linker length and composition mutagenesis

Linker length and composition were varied in libraries of single-chain Arc repressor, resulting in proteins with effective concentrations ranging over six orders of magnitude (10 μM–10 M). Linkers of 11 residues or more were required for biological activity. Equilibrium stability varied substantially with linker length, reaching a maximum for glycine-rich linkers containing 19 residues. The effects of linker length on equilibrium stability arise from significant and sometimes opposing changes in folding and unfolding kinetics. By fixing the linker length at 19 residues and varying the ratio of Ala/Gly or Ser/Gly in a 16-residue-randomized region, the effects of linker flexibility were examined. In these libraries, composition rather than sequence appears to determine stability. Maximum stability in the Ala/Gly library was observed for a protein containing 11 alanines and five glycines in the randomized region of the linker. In the Ser/Gly library, the most stable protein had seven serines and nine glycines in this region. Analysis of folding and unfolding rates suggests that alanine acts largely by accelerating folding, whereas serine acts predominantly to slow unfolding. These results demonstrate an important role for linker design in determining the stability and folding kinetics of single-chain proteins and suggest strategies for optimizing these parameters.

Random circular permutation of genes and expressed polypeptide chains: application of the method to the catalytic chains of aspartate transcarbamoylase.

Recent studies on proteins whose N and C termini are in close proximity have demonstrated that folding of polypeptide chains and assembly of oligomers can be accomplished with circularly permuted chains. As yet no methodical study has been conducted to determine how extensively new termini can be introduced and where such termini cannot be tolerated. We have devised a procedure to generate random circular permutations of the catalytic chains of Escherichia coli aspartate transcarbamoylase (ATCase; EC 2.1.3.2) and to select clones that produce active or stable holoenzyme containing permuted chains. A tandem gene construct was made, based on the desired linkage between amino acid residues in the C- and N-terminal regions of the polypeptide chain, and this DNA was treated with a suitable restriction enzyme to yield a fragment containing the rearranged coding sequence for the chain. Circularization achieved with DNA ligase, followed by linearization at random with DNase I, and incorporation of the linearized, repaired, blunt-ended, rearranged genes into a suitable plasmid permitted the expression of randomly permuted polypeptide chains. The plasmid with appropriate stop codons also contained pyrI, the gene encoding the regulatory chain of ATCase. Colonies expressing detectable amounts of ATCase-like molecules containing permuted catalytic chains were identified by an immunoblot technique or by their ability to grow in the absence of pyrimidines in the growth medium. Sequencing of positive clones revealed a variety of novel circular permutations. Some had N and C termini within helices of the wild-type enzyme as well as deletions and insertions. Permutations were concentrated in the C-terminal domain and only few were detected in the N-terminal domain. The technique, which is adaptable generally to proteins whose N and C termini are near each other, can be of value in relating in vivo folding of nascent, growing polypeptide chains to in vitro renaturation of complete chains and determining the role of protein sequence in folding kinetics.

Evidence for nonrandom hydrophobicity structures in protein chains.

The question of whether proteins originate from random sequences of amino acids is addressed. A statistical analysis is performed in terms of blocked and random walk values formed by binary hydrophobic assignments of the amino acids along the protein chains. Theoretical expectations of these variables from random distributions of hydrophobicities are compared with those obtained from functional proteins. The results, which are based upon proteins in the SWISS-PROT data base, convincingly show that the amino acid sequences in proteins differ from what is expected from random sequences in a statistically significant way. By performing Fourier transforms on the random walks, one obtains additional evidence for nonrandomness of the distributions. We have also analyzed results from a synthetic model containing only two amino acid types, hydrophobic and hydrophilic. With reasonable criteria on good folding properties in terms of thermodynamical and kinetic behavior, sequences that fold well are isolated. Performing the same statistical analysis on the sequences that fold well indicates similar deviations from randomness as for the functional proteins. The deviations from randomness can be interpreted as originating from anticorrelations in terms of an Ising spin model for the hydrophobicities. Our results, which differ from some previous investigations using other methods, might have impact on how permissive with respect to sequence specificity protein folding process is-only sequences with nonrandom hydrophobicity distributions fold well. Other distributions give rise to energy landscapes with poor folding properties and hence did not survive the evolution.

Inhibition of immunoglobulin folding and secretion by dominant negative BiP ATPase mutants.

A group of resident ER proteins have been identified that are proposed to function as molecular chaperones. The best characterized of these is BiP/GRP78, an hsp70 homologue that binds peptides containing hydrophobic residues in vitro and unfolded or unassembled proteins in vivo. However, evidence that mammalian BiP plays a direct role in protein folding remains circumstantial. In this study, we examine how BiP interacts with a particular substrate, immunoglobulin light chain (lambda LC), during its folding. Wild-type hamster BiP and several well-characterized BiP ATPase mutants were used in transient expression experiments. We demonstrate that wild-type lambda LCs showed prolonged association with mutant BiP which inhibited their secretion. Both wild-type and mutant BiP bound only to unfolded and partially folded LCs. The wild-type BiP was released from the incompletely folded LCs, allowing them to fold and be secreted, whereas the mutant BiP was not released. As a result, the LCs that were bound to BiP mutants were unable to undergo complete disulfide bond formation and were retained in the ER. Our experiments suggest that LCs undergo both BiP-dependent and BiP-independent folding steps, demonstrating that both ATP binding and hydrolysis activities of BiP are essential for the completion of LC folding in vivo and reveal that BiP must release before disulfide bond formation can occur in that domain.

Kinetic traps in lysozyme folding.

Folding of lysozyme from hen egg white was investigated by using interrupted refolding experiments. This method makes use of a high energy barrier between the native state and transient folding intermediates, and, in contrast to conventional optical techniques, it enables one to specifically monitor the amount of native molecules during protein folding. The results show that under strongly native conditions lysozyme can refold on parallel pathways. The major part of the lysozyme molecules (86%) refold on a slow kinetic pathway with well-populated partially folded states. Additionally, 14% of the molecules fold faster. The rate constant of formation of native molecules on the fast pathway corresponds well to the rate constant expected for folding to occur by a two-state process without any detectable intermediates. The results suggest that formation of the native state for the major fraction of lysozyme molecules is retarded compared with the direct folding process. Partially structured intermediates that transiently populate seem to be kinetically trapped in a conformation that can only slowly reach the native structure.

On constructing folding heteropolymers.

Simplified models of the protein-folding process have led to valuable insights into the generic properties of the folding of heteropolymers. On the basis of theoretical arguments, Shakhnovich and Gutin [(1993) Proc. Natl. Acad. Sci. USA 90, 7195-7199] have proposed a specific method to generate folding sequences for one of these. Here we present a model of folding in heteropolymers that is comparable in simplicity but different in spirit to the one studied by Shakhnovich and Gutin. In our model, the proposed recipe for constructing folding sequence fails. We find that, as a rule, the construction of folding sequences is impossible to achieve by looking at the native conformation only. Rather, competing conformations have to be taken into account too. An evolutionary algorithm that generates folding sequences by optimizing both stability of the native state and folding time is described. Remarkably, this algorithm produces, among others, sequences that fold reproducibly to metastable states.

Gas-phase folding and unfolding of cytochrome c cations.

Water is thought to play a dominant role in protein folding, yet gaseous multiply protonated proteins from which the water has been completely removed show hydrogen/deuterium (H/D) exchange behavior similar to that used to identify conformations in solution. Indicative of the gas-phase accessibility to D2O, multiply-charged (6+ to 17+) cytochrome c cations exchange at six (or more) distinct levels of 64 to 173 out of 198 exchangeable H atoms, with the 132 H level found at charge values 8+ to 17+. Infrared laser heating and fast collisions can apparently induce ions to unfold to exchange at a higher distinct level, while charge-stripping ions to lower charge values yields apparent folding as well as unfolding.

Significance of chaperonin 10-mediated inhibition of ATP hydrolysis by chaperonin 60

Chaperonins are essential for the folding of proteins in bacteria, mitochondria, and chloroplasts. We have functionally characterized the yeast mitochondrial chaperonins hsp60 and hsp10. In the presence of ADP, one molecule of hsp10 binds to hsp60 with an apparent Kd of 0.9 nM and a second molecule of hsp10 binds with a Kd of 24 nM. In the presence of ATP, the purified yeast chaperonins mediate the refolding of mitochondrial malate dehydrogenase. Hsp10 inhibits the ATPase activity of hsp60 by about 40%. Hsp10(P36H) is a point mutant of hsp10 that confers temperature-sensitive growth to yeast. Consistent with the in vivo phenotype, refolding of mitochondrial malate dehydrogenase in the presence of purified hsp10(P36H) and hsp60 is reduced at 25°C and abolished at 30°C. The affinity of hsp10(P36H) to hsp60 as well as to Escherichia coli GroEL is reduced. However, this decrease in affinity does not correlate with the functional defect, because hsp10(P36H) fully assists the GroEL-mediated refolding of malate dehydrogenase at 30°C. Refolding activity, rather, correlates with the ability of hsp10(P36H) to inhibit the ATPase of GroEL but not that of hsp60. Based on our findings, we propose that the inhibition of ATP hydrolysis is mechanistically coupled to chaperonin-mediated protein folding.

Respiratory chain is required to maintain oxidized states of the DsbA-DsbB disulfide bond formation system in aerobically growing Escherichia coli cells

DsbA, the disulfide bond catalyst of Escherichia coli, is a periplasmic protein having a thioredoxin-like Cys-30-Xaa-Xaa-Cys-33 motif. The Cys-30–Cys-33 disulfide is donated to a pair of cysteines on the target proteins. Although DsbA, having high oxidizing potential, is prone to reduction, it is maintained essentially all oxidized in vivo. DsbB, an integral membrane protein having two pairs of essential cysteines, reoxidizes DsbA that has been reduced upon functioning. It is not known, however, what might provide the overall oxidizing power to the DsbA–DsbB disulfide bond formation system. We now report that E. coli mutants defective in the hemA gene or in the ubiA-menA genes markedly accumulate the reduced form of DsbA during growth under the conditions of protoheme deprivation as well as ubiquinone/menaquinone deprivation. Disulfide bond formation of β-lactamase was impaired under these conditions. Intracellular state of DsbB was found to be affected by deprivation of quinones, such that it accumulates first as a reduced form and then as a form of a disulfide-linked complex with DsbA. This is followed by reduction of the bulk of DsbA molecules. These results suggest that the respiratory electron transfer chain participates in the oxidation of DsbA, by acting primarily on DsbB. It is remarkable that a cellular catalyst of protein folding is connected to the respiratory chain.

All cyclophilins and FK506 binding proteins are, individually and collectively, dispensable for viability in Saccharomyces cerevisiae

The cyclophilins and FK506 binding proteins (FKBPs) bind to cyclosporin A, FK506, and rapamycin and mediate their immunosuppressive and toxic effects, but the physiological functions of these proteins are largely unknown. Cyclophilins and FKBPs are ubiquitous and highly conserved enzymes that catalyze peptidyl-prolyl isomerization, a rate-limiting step during in vitro protein folding. We have addressed their functions by a genetic approach in the yeast Saccharomyces cerevisiae. Five cyclophilins and three FKBPs previously were identified in yeast. We identified four additional enzymes: Cpr6 and Cpr7, which are homologs of mammalian cyclophilin 40 that have also recently been independently isolated by others, Cpr8, a homolog of the secretory pathway cyclophilin Cpr4, and Fpr4, a homolog of the nucleolar FKBP, Fpr3. None of the eight cyclophilins or four FKBPs were essential. Surprisingly, yeast mutants lacking all 12 immunophilins were viable, and the phenotype of the dodecuplet mutant resulted from simple addition of the subtle phenotypes of each individual mutation. We conclude that cyclophilins and FKBPs do not play an essential general role in protein folding and find little evidence of functional overlap between the different enzymes. We propose that each cyclophilin and FKBP instead regulates a restricted number of unique partner proteins that remain to be identified.

Chaperonin filaments: The archaeal cytoskeleton?

Chaperonins are high molecular mass double-ring structures composed of 60-kDa protein subunits. In the hyperthermophilic archaeon Sulfolobus shibatae the two chaperonin proteins represent ≈4% of its total protein and have a combined intracellular concentration of >30 mg/ml. At concentrations ≥ 0.5 mg/ml purified chaperonins form filaments in the presence of Mg2+ and nucleotides. Filament formation requires nucleotide binding (not hydrolysis), and occurs at physiological temperatures in biologically relevant buffers, including a buffer made from cell extracts. These observations suggest that chaperonin filaments may exist in vivo and the estimated 4600 chaperonins per cell suggest that such filaments could form an extensive cytostructure. We observed filamentous structures in unfixed, uranyl-acetate-stained S. shibatae cells, which resemble the chaperonin filaments in size and appearance. ImmunoGold (Janssen) labeling using chaperonin antibodies indicated that many chaperonins are associated with insoluble cellular structures and these structures appear to be filamentous in some areas, although they could not be uranyl-acetate-stained. The existence of chaperonin filaments in vivo suggests a mechanism whereby their protein-folding activities can be regulated. More generally, the filaments themselves may play a cytoskeletal role in Archaea.

Hydrophobic sequence minimization of the α-lactalbumin molten globule

The molten globule, a widespread protein-folding intermediate, can attain a native-like backbone topology, even in the apparent absence of rigid side-chain packing. Nonetheless, mutagenesis studies suggest that molten globules are stabilized by some degree of side-chain packing among specific hydrophobic residues. Here we investigate the importance of hydrophobic side-chain diversity in determining the overall fold of the α-lactalbumin molten globule. We have replaced all of the hydrophobic amino acids in the sequence of the helical domain with a representative amino acid, leucine. Remarkably, the minimized molecule forms a molten globule that retains many structural features characteristic of a native α-lactalbumin fold. Thus, nonspecific hydrophobic interactions may be sufficient to determine the global fold of a protein.

Tissue-specific Expression of Dominant Negative Mutant Drosophila HSC70 Causes Developmental Defects and Lethality

The Drosophila melanogaster HSC3 and HSC4 genes encode Hsc70 proteins homologous to the mammalian endoplasmic reticulum (ER) protein BiP and the cytoplasmic clathrin uncoating ATPase, respectively. These proteins possess ATP binding/hydrolysis activities that mediate their ability to aid in protein folding by coordinating the sequential binding and release of misfolded proteins. To investigate the roles of HSC3 (Hsc3p) and HSC4 (Hsc4p) proteins during development, GAL4-targeted gene expression was used to analyze the effects of producing dominant negatively acting Hsc3p (D231S, K97S) and Hsc4p (D206S, K71S) proteins, containing single amino acid substitutions in their ATP-binding domains, in specific tissues of Drosophila throughout development. We show that the production of each mutant protein results in lethality over a range of developmental stages, depending on the levels of protein produced and which tissues are targeted. We demonstrate that the functions of both Hsc3p and Hsc4p are required for proper tissue establishment and maintenance. Production of mutant Hsc4p, but not Hsc3p, results in induction of the stress-inducible Hsp70 at normal temperatures. Evidence is presented that lethality is caused by tissue-specific defects that result from a global accumulation of misfolded protein caused by lack of functional Hsc70. We show that both mutant Hsc3ps are defective in ATP-induced substrate release, although Hsc3p(D231S) does undergo an ATP-induced conformational change. We believe that the amino acid substitutions in Hsc3p interfere with the structural coupling of ATP binding to substrate release, and this defect is the basis for the mutant proteins’ dominant negative effects in vivo.

Functions of FKBP12 and Mitochondrial Cyclophilin Active Site Residues In Vitro and In Vivo in Saccharomyces cerevisiae

Cyclophilin and FK506 binding protein (FKBP) accelerate cis–trans peptidyl-prolyl isomerization and bind to and mediate the effects of the immunosuppressants cyclosporin A and FK506. The normal cellular functions of these proteins, however, are unknown. We altered the active sites of FKBP12 and mitochondrial cyclophilin from the yeast Saccharomyces cerevisiae by introducing mutations previously reported to inactivate these enzymes. Surprisingly, most of these mutant enzymes were biologically active in vivo. In accord with previous reports, all of the mutant enzymes had little or no detectable prolyl isomerase activity in the standard peptide substrate-chymotrypsin coupled in vitro assay. However, in a variation of this assay in which the protease is omitted, the mutant enzymes exhibited substantial levels of prolyl isomerase activity (5–20% of wild-type), revealing that these mutations confer sensitivity to protease digestion and that the classic in vitro assay for prolyl isomerase activity may be misleading. In addition, the mutant enzymes exhibited near wild-type activity with two protein substrates, dihydrofolate reductase and ribonuclease T1, whose folding is accelerated by prolyl isomerases. Thus, a number of cyclophilin and FKBP12 “active-site” mutants previously identified are largely active but protease sensitive, in accord with our findings that these mutants display wild-type functions in vivo. One mitochondrial cyclophilin mutant (R73A), and also the wild-type human FKBP12 enzyme, catalyze protein folding in vitro but lack biological activity in vivo in yeast. Our findings provide evidence that both prolyl isomerase activity and other structural features are linked to FKBP and cyclophilin in vivo functions and suggest caution in the use of these active-site mutations to study FKBP and cyclophilin functions.

Equilibrium constants and free energies in unfolding of proteins in urea solutions

A novel thermodynamic approach to the reversible unfolding of proteins in aqueous urea solutions has been developed based on the premise that urea ligands are bound cooperatively to the macromolecule. When successive stoichiometric binding constants have values larger than expected from statistical effects, an equation for moles of bound urea can be derived that contains imaginary terms. For a very steep unfolding curve, one can then show that the fraction of protein unfolded, B̄, depends on the square of the urea concentration, U, and is given by \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} \begin{equation*}\bar {B}=\frac{{\mathit{A}}^{{\mathit{2}}}_{{\mathit{1}}}{\mathit{e}}^{{\mathrm{{\lambda}}}n\bar {B}}{\mathit{U}}^{{\mathit{2}}}}{{\mathrm{1\hspace{.167em}+\hspace{.167em}}}{\mathit{A}}^{{\mathrm{2}}}_{{\mathrm{1}}}{\mathit{e}}^{{\mathrm{{\lambda}}}\bar {B}}{\mathit{U}}^{{\mathrm{2}}}}{\mathrm{.}}\end{equation*}\end{document} A12 is the binding constant as B̄→ 0, and λ is a parameter that reflects the augmentation in affinities of protein for urea as the moles bound increases to the saturation number, n. This equation provides an analytic expression that reproduces the unfolding curve with good precision, suggests a simple linear graphical procedure for evaluating A12 and λ, and leads to the appropriate standard free energy changes. The calculated ΔG° values reflect the coupling of urea binding with unfolding of the protein. Some possible implications of this analysis to protein folding in vivo are described.