Preferential exclusion of sucrose from recombinant interleukin-1 receptor antagonist: Role in restricted conformational mobility and compaction of native state

Understanding the mechanism for sucrose-induced protein stabilization is important in many diverse fields, ranging from biochemistry and environmental physiology to pharmaceutical science. Timasheff and Lee [Lee, J. C. & Timasheff, S. N. (1981) J. Biol. Chem. 256, 7193–7201] have established that thermodynamic stabilization of proteins by sucrose is due to preferential exclusion of the sugar from the protein’s surface, which increases protein chemical potential. The current study measures the preferential exclusion of 1 M sucrose from a protein drug, recombinant interleukin 1 receptor antagonist (rhIL-1ra). It is proposed that the degree of preferential exclusion and increase in chemical potential are directly proportional to the protein surface area and that, hence, the system will favor the protein state with the smallest surface area. This mechanism explains the observed sucrose-induced restriction of rhIL-1ra conformational fluctuations, which were studied by hydrogen–deuterium exchange and cysteine reactivity measurements. Furthermore, infrared spectroscopy of rhlL-1ra suggested that a more ordered native conformation is induced by sucrose. Electron paramagnetic resonance spectroscopy demonstrated that in the presence of sucrose, spin-labeled cysteine 116 becomes more buried in the protein’s interior and that the hydrodynamic diameter of the protein is reduced. The preferential exclusion of sucrose from the protein and the resulting shift in the equilibrium between protein states toward the most compact conformation account for sucrose-induced effects on rhIL-1ra.

A transient expansion of the native state precedes aggregation of recombinant human interferon-γ

Aggregation of proteins, even under conditions favoring the native state, is a ubiquitous problem in biotechnology and biomedical engineering. Providing a mechanistic basis for the pathways that lead to aggregation should allow development of rational approaches for its prevention. We have chosen recombinant human interferon-γ (rhIFN-γ) as a model protein for a mechanistic study of aggregation. In the presence of 0.9 M guanidinium hydrochloride, rhIFN-γ aggregates with first order kinetics, a process that is inhibited by addition of sucrose. We describe a pathway that accounts for both the observed first-order aggregation of rhIFN-γ and the effect of sucrose. In this pathway, aggregation proceeds through a transient expansion of the native state. Sucrose shifts the equilibrium within the ensemble of rhIFN-γ native conformations to favor the most compact native species over more expanded ones, thus stabilizing rhIFN-γ against aggregation. This phenomenon is attributed to the preferential exclusion of sucrose from the protein surface. In addition, kinetic analysis combined with solution thermodynamics shows that only a small (9%) expansion surface area is needed to form the transient native state that precedes aggregation. The approaches used here link thermodynamics and aggregation kinetics to provide a powerful tool for understanding both the pathway of protein aggregation and the rational use of excipients to inhibit the process.

Denaturant-induced movement of the transition state of protein folding revealed by high-pressure stopped-flow measurements

The small all-β protein tendamistat folds and unfolds with two-state kinetics. We determined the volume changes associated with the folding process by performing kinetic and equilibrium measurements at variable pressure between 0.1 and 100 MPa (1 to 1,000 bar). GdmCl-induced equilibrium unfolding transitions reveal that the volume of the native state is increased by 41.4 ± 2.0 cm3/mol relative to the unfolded state. This value is virtually independent of denaturant concentration. The use of a high-pressure stopped-flow instrument enabled us to measure the activation volumes for the refolding (ΔVf0‡) and unfolding reaction (ΔVu0‡) over a broad range of GdmCl concentrations. The volume of the transition state is 60% native-like (ΔVf0‡ = 25.0 ± 1.2 cm3/mol) in the absence of denaturant, indicating partial solvent accessibility of the core residues. The volume of the transition state increases linearly with denaturant concentration and exceeds the volume of the native state above 6 M GdmCl. This result argues for a largely desolvated transition state with packing deficiencies at high denaturant concentrations and shows that the structure of the transition state depends strongly on the experimental conditions.

Characterization of residual structure in the thermally denatured state of barnase by simulation and experiment: Description of the folding pathway

Residual structure in the denatured state of a protein may contain clues about the early events in folding. We have simulated by molecular dynamics the denatured state of barnase, which has been studied by NMR spectroscopy. An ensemble of 104 structures was generated after 2 ns of unfolding and following for a further 2 ns. The ensemble was heterogeneous, but there was nonrandom, residual structure with persistent interactions. Helical structure in the C-terminal portion of helix α1 (residues 13–17) and in helix α2 as well as a turn and nonnative hydrophobic clustering between β3 and β4 were observed, consistent with NMR data. In addition, there were tertiary contacts between residues in α1 and the C-terminal portion of the β-sheet. The simulated structures allow the rudimentary NMR data to be fleshed out. The consistency between simulation and experiment inspires confidence in the methods. A description of the folding pathway of barnase from the denatured to the native state can be constructed by combining the simulation with experimental data from φ value analysis and NMR.

The pKa of His-24 in the folding transition state of apomyoglobin

In native apomyoglobin, His-24 cannot be protonated, although at pH 4 the native protein forms a molten globule folding intermediate in which the histidine residues are readily protonated. The inability to protonate His-24 in the native protein dramatically affects the unfolding/refolding kinetics, as demonstrated by simulations for a simple model. Kinetic data for wild type and for a mutant lacking His-24 are analyzed. The pKa values of histidine residues in native apomyoglobin are known from earlier studies, and the average histidine pKa in the molten globule is determined from the pH dependence of the equilibrium between the native and molten globule forms. Analysis of the pH-dependent unfolding/refolding kinetics reveals that the average pKa of the histidine residues, including His-24, is closely similar in the folding transition state to the value found in the molten globule intermediate. Consequently, protonation of His-24 is not a barrier to refolding of the molten globule to the native protein. Instead, the normal pKa of His-24 in the transition state, coupled with its inaccessibility in the native state, promotes fast unfolding at low pH. The analysis of the wild-type results is confirmed and extended by using the wild-type parameters to fit the unfolding kinetics of a mutant lacking His-24.

Partial molar volume, surface area, and hydration changes for equilibrium unfolding and formation of aggregation transition state: High-pressure and cosolute studies on recombinant human IFN-γ

The equilibrium dissociation of recombinant human IFN-γ was monitored as a function of pressure and sucrose concentration. The partial molar volume change for dissociation was −209 ± 13 ml/mol of dimer. The specific molar surface area change for dissociation was 12.7 ± 1.6 nm2/molecule of dimer. The first-order aggregation rate of recombinant human IFN-γ in 0.45 M guanidine hydrochloride was studied as a function of sucrose concentration and pressure. Aggregation proceeded through a transition-state species, N*. Sucrose reduced aggregation rate by shifting the equilibrium between native state (N) and N* toward the more compact N. Pressure increased aggregation rate through increased solvation of the protein, which exposes more surface area, thus shifting the equilibrium away from N toward N*. The changes in partial molar volume and specific molar surface area between the N* and N were −41 ± 9 ml/mol of dimer and 3.5 ± 0.2 nm2/molecule, respectively. Thus, the structural change required for the formation of the transition state for aggregation is small relative to the difference between N and the dissociated state. Changes in waters of hydration were estimated from both specific molar surface area and partial molar volume data. From partial molar volume data, estimates were 25 and 128 mol H2O/mol dimer for formation of the aggregation transition state and for dissociation, respectively. From surface area data, estimates were 27 and 98 mol H2O/mol dimer. Osmotic stress theory yielded values ≈4-fold larger for both transitions.

The folding mechanism of larger model proteins: role of native structure.

The folding mechanism of a 125-bead heteropolymer model for proteins is investigated with Monte Carlo simulations on a cubic lattice. Sequences that do and do not fold in a reasonable time are compared. The overall folding behavior is found to be more complex than that of models for smaller proteins. Folding begins with a rapid collapse followed by a slow search through the semi-compact globule for a sequence-dependent stable core with about 30 out of 176 native contacts which serves as the transition state for folding to a near-native structure. Efficient search for the core is dependent on structural features of the native state. Sequences that fold have large amounts of stable, cooperative structure that is accessible through short-range initiation sites, such as those in anti-parallel sheets connected by turns. Before folding is completed, the system can encounter a second bottleneck, involving the condensation and rearrangement of surface residues. Overly stable local structure of the surface residues slows this stage of the folding process. The relation of the results from the 125-mer model studies to the folding of real proteins is discussed.

Definition of amide protection factors for early kinetic intermediates in protein folding

Hydrogen–deuterium exchange experiments have been used previously to investigate the structures of well defined states of a given protein. These include the native state, the unfolded state, and any intermediates that can be stably populated at equilibrium. More recently, the hydrogen–deuterium exchange technique has been applied in kinetic labeling experiments to probe the structures of transiently formed intermediates on the kinetic folding pathway of a given protein. From these equilibrium and nonequilibrium studies, protection factors are usually obtained. These protection factors are defined as the ratio of the rate of exchange of a given backbone amide when it is in a fully solvent-exposed state (usually obtained from model peptides) to the rate of exchange of that amide in some state of the protein or in some intermediate on the folding pathway of the protein. This definition is straightforward for the case of equilibrium studies; however, it is less clear-cut for the case of transient kinetic intermediates. To clarify the concept for the case of burst-phase intermediates, we have introduced and mathematically defined two different types of protection factors: one is Pstruc, which is more related to the structure of the intermediate, and the other is Papp, which is more related to the stability of the intermediate. Kinetic hydrogen–deuterium exchange data from disulfide-intact ribonuclease A and from cytochrome c are discussed to explain the use and implications of these two definitions.

Exploring the folding free energy surface of a three-helix bundle protein

The multidimensional free energy surface for a small fast folding helical protein is explored based on first-principle calculations. The model represents the 46-residue segment from fragment B of staphylococcal protein A. The relationship between collapse and tertiary structure formation, and the order of collapse and secondary structure formation, are investigated. We find that the initial collapse process gives rise to a transition state with about 30% of the native tertiary structure and 50–70% of the native helix content. We also observe two distinct distributions of native helix in this collapsed state (Rg ≈ 12 Å), one with about 20% of the native helical hydrogen bonds, the other with near 70%. The former corresponds to a local minimum. The barrier from this metastable state to the native state is about 2 kBT. In the latter case, folding is essentially a downhill process involving topological assembly. In addition, the order of formation of secondary structure among the three helices is examined. We observe cooperative formation of the secondary structure in helix I and helix II. Secondary structure in helix III starts to form following the formation of certain secondary structure in both helix I and helix II. Comparisons of our results with those from theory and experiment are made.

Folding in vivo of a newly translated yeast cytosolic enzyme is mediated by the SSA class of cytosolic yeast Hsp70 proteins

The nature of chaperone action in the eukaryotic cytosol that assists newly translated cytosolic proteins to reach the native state has remained poorly defined. Actin, tubulin, and Gα transducin are assisted by the cytosolic chaperonin, CCT, but many other proteins, for example, ornithine transcarbamoylase (OTC), a cytosolic homotrimeric enzyme of yeast, do not require CCT action. Here, we observe that yeast cytosolic OTC is assisted to its native state by the SSA class of yeast cytosolic Hsp70 proteins. In vitro, refolding of OTC diluted from denaturant was assisted by crude yeast cytosol and ATP and found to be directed by SSA1/2. In vivo, when OTC was induced in a temperature-sensitive SSA-deficient strain, it exhibited reduced specific activity, and nonnative subunits were detected in the soluble fraction. These findings indicate that, in vivo, the Hsp70 system assists in folding at least some newly translated cytosolic enzymes, most likely functioning in a posttranslational manner.

Detection of residue contacts in a protein folding intermediate

Protein folding can be described in terms of the development of specific contacts between residues as a highly disordered polypeptide chain converts into the native state. Here we describe an NMR based strategy designed to detect such contacts by observation of nuclear Overhauser effects (NOEs). Experiments with α-lactalbumin reveal the existence of extensive NOEs between aromatic and aliphatic protons in the archetypal molten globule formed by this protein at low pH. Analysis of their time development provides direct evidence for near-native compactness of this state. Through a rapid refolding procedure the NOE intensity can be transferred efficiently into the resolved and assigned spectrum of the native state. This demonstrates the viability of using this approach to map out time-averaged interactions between residues in a partially folded protein.

Structural distribution of stability in a thermophilic enzyme

Stability parameters for individual residues in Thermus thermophilus cysteine-free RNase H were determined by native state hydrogen exchange, thus providing a unique comparison of regional thermodynamics between thermophilic and mesophilic homologues. The general distribution of stability in the thermophilic protein is similar to that of its mesophilic homologue, with a proportional increase in stability for almost all residues. As a consequence, the residue-specific stabilities of the two proteins are remarkably similar under conditions where their global stabilities are the same. These results indicate that T. thermophilus RNase H is stabilized in a delocalized fashion, preserving a finely tuned balance of stabilizing interactions throughout the structure. Therefore, although protein stability can be altered by single amino acid substitution, evolution for optimal function may require more subtle and delocalized mechanisms.

Folding thermodynamics of a model three-helix-bundle protein

The calculated folding thermodynamics of a simple off-lattice three-helix-bundle protein model under equilibrium conditions shows the experimentally observed protein transitions: a collapse transition, a disordered-to-ordered globule transition, a globule to native-state transition, and the transition from the active native state to a frozen inactive state. The cooperativity and physical origin of the various transitions are explored with a single “optimization” parameter and characterized with the Lindemann criterion for liquid versus solid-state dynamics. Below the folding temperature, the model has a simple free energy surface with a single basin near the native state; the surface is similar to that calculated from a simulation of the same three-helix-bundle protein with an all-atom representation [Boczko, E. M. & Brooks III, C. L. (1995) Science 269, 393–396].

Amphipols: Polymers that keep membrane proteins soluble in aqueous solutions

Amphipols are a new class of surfactants that make it possible to handle membrane proteins in detergent-free aqueous solution as though they were soluble proteins. The strongly hydrophilic backbone of these polymers is grafted with hydrophobic chains, making them amphiphilic. Amphipols are able to stabilize in aqueous solution under their native state four well-characterized integral membrane proteins: (i) bacteriorhodopsin, (ii) a bacterial photosynthetic reaction center, (iii) cytochrome b6f, and (iv) matrix porin.

Staphylococcal protein A: Unfolding pathways, unfolded states, and differences between the B and E domains

We present multiple native and denaturation simulations of the B and E domains of the three-helix bundle protein A, totaling 60 ns. The C-terminal helix (H3) consistently denatures later than either of the other two helices and contains residual helical structure in the denatured state. These results are consistent with experiments suggesting that the isolated H3 fragment is more stable than H1 and H2 and that H3 forms early in folding. Interestingly, the denatured state of the B domain is much more compact than that of the E domain. This sequence-dependent effect on the dimensions of the denatured state and the lack of correlation with structure suggest that the radius of gyration can be a misleading reaction coordinate for unfolding/folding. Various unfolding and refolding events are observed in the denaturation simulations. In some cases, the transitions are facilitated through interactions with other portions of the protein—contact-assisted helix formation. In the native simulations, the E domain is very stable: after 6 ns, the Cα root-mean-square deviation from the starting structure is less than 1.4 Å. In contrast, the native state of the B domain deviates more and its inter-helical angles fluctuate. In apparent contrast, we note that the B domain is thermodynamically more stable than the E domain. The simulations suggest that the increased stability of the B domain may be due to heightened mobility, and therefore entropy, in the native state and decreased mobility/entropy in the more compact denatured state.