Quantitative protein stability measurement in vivo.

The equilibrium between the native and denatured states of a protein can be key to its function and regulation. Traditionally, the folding equilibrium constant has been measured in vitro using purified protein and simple buffers. However, the biological environment of proteins can differ from these in vitro conditions in ways that could significantly perturb stability. Here, we present the first quantitative comparison between the stability of a protein in vitro and in the cytoplasm of Escherichia coli using amide hydrogen exchange detected by MALDI mass spectrometry (SUPREX). The results indicate that the thermodynamic stability of monomeric lambda repressor within the cell is the same as its stability measured in a simple buffer in vitro. However, when the E. coli are placed in a hyperosmotic environment, the in vivo stability is greatly enhanced. The in vivo SUPREX method provides a general and quantitative way to measure protein stabilities in the cell and will be useful for applications where intracellular stability information provides important biological insights.

Phosphorylation of RNA polymerase II CTD fragments results in tight binding to the WW domain from the yeast prolyl isomerase Ess1.

The yeast prolyl isomerase, Ess1, has recently been shown to interact via its WW domain with the hyperphosphorylated form of the RNA polymerase II C-terminal domain (CTD). We have investigated folding of the Ess1 WW domain and its binding to peptides representing the CTD by circular dichroism and fluorescence. Ess1 WW folds and unfolds reversibly, but in the absence of ligand is only marginally stable with a melting temperature of 19 degrees C. The WW domain is stabilized by the addition of anionic ligands, namely, chloride, inorganic phosphate, phosphoserine, and phosphorylated CTD peptides. Dissociation constants were measured to be 70--100 microM for CTD peptides phosphorylated at one serine, and 16--21 microM for peptides with two or more phosphorylated serines. Weaker or no affinity was observed for nonphosphorylated CTD peptides. There is surprisingly little difference in the affinity for peptides phosphorylated at Ser 2 or Ser 5 of the consensus repeat, or for peptides with different patterns of multiple phosphorylation. The binding of Ess1 to phosphorylated CTD peptides is consistent with a model wherein the WW domain positions Ess1 to catalyze isomerization of the many pSer--Pro peptide bonds in the phosphorylated CTD. We suggest that cis/trans isomerization of prolyl peptide bonds plays a crucial role in CTD function during eukaryotic transcription.

Preorganized secondary structure as an important determinant of fast protein folding.

The folding and unfolding kinetics of the B-domain of staphylococcal protein A, a small three-helix bundle protein, were probed by NMR. The lineshape of a single histidine resonance was fit as a function of denaturant to give folding and unfolding rate constants. The B-domain folds extremely rapidly in a two-state manner, with a folding rate constant of 120,000 s-1, making it one of the fastest-folding proteins known. Diffusion-collision theory predicts folding and unfolding rate constants that are in good agreement with the experimental values. The apparent rate constant as a function of denaturant ('chevron plot') is predicted within an order of magnitude. Our results are consistent with a model whereby fast-folding proteins utilize a diffusion-collision mechanism, with the preorganization of one or more elements of secondary structure in the unfolded protein.

Linked folding and anion binding of the Bacillus subtilis ribonuclease P protein.

Ribonuclease P (RNase P) is the endoribonuclease responsible for the 5'-maturation of precursor tRNA transcripts. In bacteria, RNase P is composed of a catalytic RNA subunit and an associated protein subunit that enhances the substrate specificity of the holoenzyme. We have initiated a study of the biophysical properties of the protein subunit from Bacillus subtilis RNase P (P protein) toward the goal of understanding the thermodynamics of RNase P holoenzyme assembly. The P protein is predominantly unfolded in 10 mM sodium cacodylate at neutral pH based on circular dichroism and NMR studies and therefore has several characteristics typical of "intrinsically unstructured" proteins. Furthermore, the P protein folds to its native alpha/beta structure upon addition of various small molecule anions. Anion-induced folding is best attributed to the binding of these anions to the folded state of the protein, and a model is presented which describes the observed tightly coupled folding and binding phenomena. The P protein also undergoes a cooperative folding transition upon addition of the osmolyte trimethylamine N-oxide (TMAO). The equilibrium constant of folding (K(fold)) at 37 degrees C for the P protein was determined to be 0.0071 +/- 0.0005 using a two-state folding model to describe the TMAO titration data. Thus, the folding and binding equilibria observed in the anion-induced folding of the P protein can be uncoupled to determine the intrinsic binding affinities (K(a)'s) of the anionic ligands. Evidence that the osmolyte-induced and the ligand-induced folded conformations of the P protein are structurally similar is also presented.

A quantitative, high-throughput screen for protein stability.

In proteomic research, it is often necessary to screen a large number of polypeptides for the presence of stable structure. Described here is a technique (referred to as SUPREX, stability of unpurified proteins from rates of H/D exchange) for measuring the stability of proteins in a rapid, high-throughput fashion. The method uses hydrogen exchange to estimate the stability of microgram quantities of unpurified protein extracts by using matrix-assisted laser desorption/ionization MS. The stabilities of maltose binding protein and monomeric lambda repressor variants determined by SUPREX agree well with stability data obtained from conventional CD denaturation of purified protein. The method also can detect the change in stability caused by the binding of maltose to maltose binding protein. The results demonstrate the precision of the method over a wide range of stabilities.

Novel disulfide engineering in human carbonic anhydrase II using the PAIRWISE side-chain geometry database.

An analysis of the pairwise side-chain packing geometries of cysteine residues observed in high-resolution protein crystal structures indicates that cysteine pairs have pronounced orientational preferences due to the geometric constraints of disulfide bond formation. A potential function was generated from these observations and used to evaluate models for novel disulfide bonds in human carbonic anhydrase II (HCAII). Three double-cysteine variants of HCAII were purified and the effective concentrations of their thiol groups were determined by titrations with glutathione and dithiothreitol. The effects of the cysteine mutations on the native state structure and stability were characterized by circular dichroism, enzymatic activity, sulfonamide binding, and guanidine hydrochloride titration. These analyses indicate that the PAIRWISE potential is a good predictor of the strength of the disulfide bond itself, but the overall structural and thermodynamic effects on the protein are complicated by additional factors. In particular, the effects of cysteine substitutions on the native state and the stabilization of compact nonnative states by the disulfide can override any stabilizing effect of the cross-link.

Reinterpretation of GCN4-p1 folding kinetics: partial helix formation precedes dimerization in coiled coil folding.

The folding of coiled coil peptides has traditionally been interpreted in terms of native dimer and unfolded monomers. Calculations using AGADIR and experimental studies of fragments suggest that the monomers of the coiled coil peptide, GCN4-p1, contain significant residual helical structure. A simple model based on diffusion-collision theory predicts not only the measured folding rate within an order of magnitude, but also predicts remarkably well the effect of alanine to glyXcine mutations. We suggest that intrinsic helix stability is a major determinant of the folding rate of the GCN4 coiled coil.

Contribution of a buried hydrogen bond to lambda repressor folding kinetics.

A hydrogen bond between the buried residues Asp 14 and Ser 77 in monomeric lambda repressor has been removed by mutation of these residues to alanine. Double mutant cycles show that the interaction stabilizes the native state of the protein by 1.5 kcal/mol. Removal of the interaction affects mainly the unfolding rates and not the folding rates, suggesting that this hydrogen bond is not substantially formed in the rate-limiting steps in the folding pathways of the protein. Mutations in two versions of lambda6-85, wild type and the faster folding G46A/G48A (WT), show similar effects. Diffusion-collision correctly predicts the behavior of WT but not of wild type. Our analysis suggests that folding of helix 3 is a crucial slow step along the various folding pathways and generally occurs before the formation of the 14-77 hydrogen bond. Experiments removing tertiary interactions, combined with experiments altering helical stability and diffusion-collision calculations, provide a strategy to unravel the folding mechanisms of small helical proteins.

The structural distribution of cooperative interactions in proteins: analysis of the native state ensemble.

Cooperative interactions link the behavior of different amino acid residues within a protein molecule. As a result, the effects of chemical or physical perturbations to any given residue are propagated to other residues by an intricate network of interactions. Very often, amino acids "sense" the effects of perturbations occurring at very distant locations in the protein molecule. In these studies, we have investigated by computer simulation the structural distribution of those interactions. We show here that cooperative interactions are not intrinsically bi-directional and that different residues play different roles within the intricate network of interactions existing in a protein. The effect of a perturbation to residue j on residue k is not necessarily equal to the effect of the same perturbation to residue k on residue j. In this paper, we introduce a computer algorithm aimed at mapping the network of cooperative interactions within a protein. This algorithm exhaustively performs single site thermodynamic mutations to each residue in the protein and examines the effects of those mutations on the distribution of conformational states. The algorithm has been applied to three different proteins (lambda repressor fragment 6-85, chymotrypsin inhibitor 2, and barnase). This algorithm accounts well for the observed behavior of these proteins.

Folding kinetics of a fluorescent variant of monomeric lambda repressor.

A tryptophan-containing variant of monomeric lambda repressor has been made, and its folding kinetics were analyzed at 20 degreesC using fluorescence stopped-flow and dynamic NMR. Equilibrium denaturation curves obtained by circular dichroism, fluorescence, and NMR are superimposable. Stopped-flow analysis indicates that in the absence of denaturants the folding reaction is complete within the dead-time of the experiment. Within higher denaturant conditions, where the folding rate is slower, NMR and stopped-flow agree on the folding and unfolding rates of the protein. In 3.4 M urea and 1.8 M GdmCl, we show that the variant folds within 2 ms. Extrapolation indicates that the folding time is 20 micro(s) in the absence of denaturants. All folding and unfolding reactions displayed monoexponential kinetics, and no burst-phases were observed. In addition, the thermodynamic parameters Delta G and meq obtained from the kinetic analysis are consistent with the equilibrium experiments. The results support a two-state Dleft and right arrow N folding model.

Protein folding dynamics: quantitative comparison between theory and experiment.

The development of a quantitative kinetic scheme is a central goal in mechanistic studies of biological phenomena. For fast-folding proteins, which lack experimentally observable kinetic intermediates, a quantitative kinetic scheme describing the order and rate of events during folding has yet to be developed. In the present study, the folding mechanism of monomeric lambda repressor is described using the diffusion-collision model and estimates of intrinsic alpha-helix propensities. The model accurately predicts the folding rates of the wild-type protein and five of eight previously studied Ala left and right arrow Gly variants and suggests that the folding mechanism is distributed among multiple pathways that are highly sensitive to the amino acid sequence. For example, the model predicts that the wild-type protein folds through a small number of pathways with a folding time of 260 micros. However, the folding of a variant (G46A/G48A) is predicted to fold through a large number of pathways with a folding time of 12 micros. Both folding times quantitatively agree with the experimental values at 37 degrees C extrapolated to 0 M denaturant. The quantitative nature of the diffusion-collision model allows for rigorous experimental tests of the theory.

Determinants of backbone dynamics in native BPTI: cooperative influence of the 14-38 disulfide and the Tyr35 side-chain.

15Nitrogen relaxation experiments were used to characterize the backbone dynamics of two modified forms of bovine pancreatic trypsin inhibitor (BPTI). In one form, the disulfide between Cys14 and Cys38 in the wild-type protein was selectively reduced and methylated to generate an analog of the final intermediate in the disulfide-coupled folding pathway. The second form was generated by similarly modifying a mutant protein in which Tyr35 was replaced with Gly (Y35G). For both selectively reduced proteins, the overall conformation of native BPTI was retained, and the relaxation data for these proteins were compared with those obtained previously with the native wild-type and Y35G proteins. Removing the disulfide from either protein had only small effects on the observed longitudinal relaxation rates (R1) or heteronuclear cross relaxation rates (nuclear Overhauser effect), suggesting that the 14-38 disulfide has little influence on the fast (ps to ns) backbone dynamics of either protein. In the wild-type protein, the pattern of residues undergoing slower (micros to ms) internal motions, reflected in unusually large transverse relaxation rates (R2), was also largely unaffected by the removal of this disulfide. It thus appears that the large R2 rates previously observed in native wild-type protein are not a direct consequence of isomerization of the 14-38 disulfide. In contrast with the wild-type protein, reducing the disulfide in Y35G BPTI significantly decreased the number of backbone amides displaying large R2 rates. In addition, the frequencies of the backbone motions in the modified protein, estimated from R2 values measured at multiple refocusing delays, appear to span a wider range than those seen in native Y35G BPTI. Together, these observations suggest that the slow internal motions in Y35G BPTI are more independent in the absence of the 14-38 disulfide and that formation of this bond may lead to a substantial loss of conformational entropy. These effects may account for the previous observation that the Y35G substitution greatly destabilizes the disulfide. The results also demonstrate that the disulfide and the buried side-chain influence the dynamics of the folded protein in a highly cooperative fashion, with the effects of removing either being much greater in the absence of the other.

Enhanced protein flexibility caused by a destabilizing amino acid replacement in BPTI.

A genetically engineered variant of bovine pancreatic trypsin inhibitor (Y35G BPTI) has been shown previously by X-ray crystallography to have a three-dimensional structure dramatically different from that of the wild-type protein, particularly in the protease-binding region of the molecule. Yet, the Y35G variant is a potent trypsin inhibitor. Described here are 15N NMR relaxation studies to compare the backbone dynamics of Y35G BPTI to those of the wild-type protein. The Tyr35 --> Gly substitution increased the transverse relaxation rates of more than one third of all backbone amide groups, but had little effect on the longitudinal relaxation rates, indicating that the substitution facilitates relatively slow backbone motions, estimated to be on the microsecond time-scale. The results indicate that the residues making up the trypsin-binding site undergo large and relatively slow conformational changes in solution, estimated to be on the 5 to 20 micros time-scale. It is thus likely that the crystal structure represents only one of multiple interconverting conformations in solution, only a fraction of which may be competent for binding trypsin. The large thermodynamic destabilization associated with this substitution may arise, in part, from a loss in cooperativity among the multiple stabilizing interactions that are normally favored by the highly ordered structure of the wild-type protein. These results suggest that fully understanding the effects of amino acid replacements on the functional and thermodynamic properties of proteins may often require analysis of the dynamic, as well as the structural, properties of altered proteins.

The energy landscape of a fast-folding protein mapped by Ala-->Gly substitutions.

A moderately stable protein with typical folding kinetics unfolds and refolds many times during its cellular lifetime. In monomeric lambda repressor this process is extremely rapid, with an average folded state lifetime of only 30 milliseconds. A thermostable variant of this protein (G46A/G48A) unfolds with the wild-type rate, but it folds in approximately 20 microseconds making it the fastest-folding protein yet observed. The effects of alanine to glycine substitutions on the folding and unfolding rate constants of the G46A/G48A variant, measured by dynamic NMR spectroscopy, indicate that the transition state is an ensemble comprised of a disperse range of conformations. This structural diversity in the transition state is consistent with the idea that folding chains are directed towards the native state by a smooth funnel-like conformational energy landscape. The kinetic data for the folding of monomeric lambda repressor can be understood by merging the new energy landscape view of folding with traditional models. This hybrid model incorporates the conformational diversity of denatured and transition state ensembles, a transition state activation energy, and the importance of intrinsic helical stabilities.

High-level misincorporation of lysine for arginine at AGA codons in a fusion protein expressed in Escherichia coli.

The expression of eukaryotic genes in Escherichia coli is one of the most frequently used tools of modern science. The arginine codon AGA is a common codon in eukaryotic genes but is particularly rare in E. coli. We report here 36 to 42% misincorporation of lysine at three AGA codons in a well-expressed protein. This misincorporation yields a protein whose electrospray mass spectrum (ESMS) shows peaks at the expected mass (M), M-28, M-56 and M-84 with intensities representing 34.5(+/-0.7), 37.5(+/-1.1), 21.2(+/-1.7) and 6.6(+/-0.5) % of the total intensity, respectively. Replacement of either all three AGA codons or the two closest to the 3' end of the gene by the more common CGC arginine codon gave a protein with a single ESMS peak. Misincorporation could also be eliminated by the co-expression of the tRNA(UCL)Arg gene, argU. These studies demonstrate that misincorporation of amino acids at rare codons of recombinant proteins can be far higher than previously thought.

Microsecond protein folding through a compact transition state.

Dynamic NMR methods have been employed to measure the folding and unfolding rate constants of two extremely fast-folding proteins. lambda 6-85, a truncated, monomeric form of the N-terminal domain of lambda repressor, refolds with a lifetime of approximately 250 microseconds. These methods have also been applied to a thermostable lambda 6-85 variant with alanine substituted for glycine residues 46 and 48 in the third helix (G46A/G48A). Both proteins exhibit linear ln (kf,u) versus [urea] plots, consistent with two-state folding for both proteins. When extrapolated to 0M urea, the data indicate that G46A/G48A folds with a lifetime of less than 20 microseconds. The slopes of the ln (kf,u) versus [urea] curves (mu and mf) indicate that the modest Gly-->Ala double mutation dramatically changes the transition state solvent accessibility. The transition state for lambda 6-85 has a fractional accessibility (mu/(mu-mf)) of 0.61, whereas the transition state for G46A/G48A is much more native-like, with a fractional accessibility of 0.16. The extraordinary change in the folding pathway that these mutations induce suggests that the intrinsic stability of helix 3 is an important determinant of the folding mechanism.

Heat and cold denatured states of monomeric lambda repressor are thermodynamically and conformationally equivalent.

Although the denaturation of proteins by low temperatures is a well-documented phenomenon, little is known about the molecular details of the process. In this study, the parameters describing the denaturation thermodynamics of residues 6-85 of the N-terminal domain of lambda repressor have been determined by fitting the three-dimensional thermal-urea denaturation surface obtained by circular dichroism. The shape of the surface shows cold denaturation at low temperatures and urea concentrations above 2 M, which allows accurate determination of the apparent heat capacity of denaturation (delta Cp). Denaturation curves based on aromatic 1H NMR spectra give identical denaturation curves, confirming purely twostate folding under all conditions studies. The denaturation surface can be fit with constant delta Cp and delta In KD/delta[urea] (KD is the equilibrium constant for denaturation), consistent with a thermodynamically invariant denatured state. In addition, the aromatic 1H NMR spectrum of the cold denatured state at 0 degree C in 3 M uea is essentially identical to the spectrum at 70 degree C in 3 M urea. These observations indicate that the structures of the cold and heat denatured states, in the presence of 3 M urea, are thermodynamically and conformationally equivalent.

Submillisecond folding of monomeric lambda repressor.

The folding kinetics of a truncated form of the N-terminal domain of phage lambda repressor [lambda 6-85] has been investigated by using the technique of dynamic NMR. lambda 6-85 has been shown previously to fold in a purely two-state fashion. This allows the determination of folding and unfolding rates from simulation of the exchange-broadened aromatic resonances of Tyr-22. The folding kinetics were determined over a range of 1.35 to 3.14 M urea. The urea dependence of both folding and unfolding rate constants is exponential, suggesting that the rate-determining step is invariant at the urea concentrations studied. The folding and unfolding rates extrapolated to 0 M urea at 37 degrees C are 3600 +/- 400 s-1 and 27 +/- 6 s-1, respectively. The observed lambda 6-85 folding rate constant exceeds that of other fast-folding globular proteins by a factor of 14-54. The urea dependence of the folding and unfolding rate constants suggests that the transition state of the rate-determining step is considerably more exposed to solvent than previously studied protein-folding transition states. The surprising rapidity of lambda 6-85 folding and unfolding may be the consequence of its all-helical secondary structure. These kinetic results clearly demonstrate that all of the fundamental events of protein folding can occur on the submillisecond time scale.

Structure and stability of monomeric lambda repressor: NMR evidence for two-state folding.

The absence of equilibrium intermediates in protein folding reactions (i.e., two-state folding) simplifies thermodynamic and kinetic analyses but is difficult to prove rigorously. We demonstrate a sensitive method for detecting partially folded species based on using proton chemical shifts as local probes of structure. The coincidence of denaturation curves for probes throughout the molecule is a particularly stringent test for two-state folding. In this study we investigate a new form of the N-terminal domain of bacteriophage lambda repressor consisting of residues 6-85 (lambda 6-85) using nuclear magnetic resonance (NMR) and circular dichroism (CD). This truncated version lacks the residues required for dimerization and is monomeric under the conditions used for NMR. Heteronuclear NMR was used to assign the 1H, 15N, and backbone 13C resonances. The secondary and tertiary structure of lambda 6-85 is very similar to that reported for the crystal structure of the DNA-bound 1-92 fragment [Beamer, L. J., and Pabo, C. O. (1992) J. Mol. Biol. 227, 177-196], as judged by analysis of chemical shifts, amide hydrogen exchange, amide-alpha coupling constants, and nuclear Overhauser enhancements. Thermal and urea denaturation studies were conducted using the chemical shifts of the four aromatic side chains as local probes and the CD signal at 222 nm as a global probe. Plots of the fraction denatured versus denaturant concentration obtained from these studies are identical for all probes under all conditions studied. This observation provides strong evidence for two-state folding, indicating that there are no populated intermediates in the folding of lambda 6-85.

Propensity for a leucine zipper-like domain of human immunodeficiency virus type 1 gp41 to form oligomers correlates with a role in virus-induced fusion rather than assembly of the glycoprotein complex.

For a number of viruses, oligomerization is a critical component of envelope processing and surface expression. Previously, we reported that a synthetic peptide (DP-107) corresponding to the putative leucine zipper region (aa 553-590) of the transmembrane protein (gp41) of human immunodeficiency virus type 1 (HIV-1) exhibited alpha-helical secondary structure and self-associated as a coiled coil. In view of the tendency of this type of structure to mediate protein association, we speculated that this region of gp41 might play a role in HIV-1 envelope oligomerization. However, later it was shown that mutations which should disrupt the structural elements of this region of gp41 did not affect envelope processing, transport, or surface expression (assembly oligomerization). In this report we compare the effects of amino acid substitutions within this coiled-coil region on structure and function of both viral envelope proteins and the corresponding synthetic peptides. Our results establish a correlation between the destabilizing effects of amino acid substitutions on coiled-coil structure in the peptide model and phenotype of virus entry. These biological and physical biochemical studies do not support a role for the coiled-coil structure in mediating the assembly oligomerization of HIV-1 envelope but do imply that this region of gp41 plays a key role in the sequence of events associated with viral entry. We propose a functional role for the coiled-coil domain of HIV-1 gp41.