Characterization of the free energy spectrum of peptostreptococcal protein L

Not Always Sticky: Specificity of Protein Stabilization by Sugars Is Conferred by Protein-Water Hydrogen Bonds

Solutes added to buffered solutions directly impact protein folding. Protein stabilization by cosolutes or crowders has been shown to be largely driven by protein-cosolute volume exclusion complemented by chemical and soft interactions. By contrast to previous studies that indicate the invariably destabilizing role of soft protein-sugar attractions, we show here that soft interactions with sugar cosolutes are protein-specific and can be stabilizing or destabilizing. We experimentally follow the folding of two model miniproteins that are only marginally stable but in the presence of sugars and polyols fold into representative and distinct secondary structures: β-hairpin or α-helix. Our mean-field model reveals that while protein-sugar excluded volume interactions have a similar stabilizing effect on both proteins, the soft interactions add a destabilizing contribution to one miniprotein but further stabilize the other. Using molecular dynamics simulations, we link the soft protein-cosolute interactions to the weakening of direct protein-water hydrogen bonding due to the presence of sugars. Although these weakened hydrogen bonds destabilize both the native and denatured states of the two proteins, the resulting contribution to the folding free energy can be positive or negative depending on the amino acid sequence. This study indicates that the significant variation between proteins in their soft interactions with sugar determines the specific response of different proteins, even to the same sugar.

Effects of Molecular Crowders on Single-Molecule Nucleic Acid Folding: Temperature-Dependent Studies Reveal True Crowding vs Enthalpic Interactions

Biomolecular folding in cells can be strongly influenced by spatial overlap/excluded volume interactions (i.e., "crowding") with intracellular solutes. As a result, traditional in vitro experiments with dilute buffers may not accurately recapitulate biomolecule folding behavior in vivo. In order to account for such ubiquitous excluded volume effects, biologically inert polyethylene glycol (PEG) and polysaccharides (dextran and Ficoll) are often used as in vitro crowding agents to mimic in vivo crowding conditions, with a common observation that high concentrations of these polymers stabilize the more compact biomolecule conformation. However, such an analysis can be distorted by differences in polymer interactions with the folded vs unfolded conformers, requiring temperature-dependent analysis of the thermodynamics to reliably assess competing enthalpic vs entropic contributions and thus the explicit role of excluded volume. In this work, temperature-controlled single-molecule fluorescence resonance energy transfer (smFRET) is used to characterize the thermodynamic interaction between nucleic acids and common polymer crowders PEG, dextran, and Ficoll. The results reveal that PEG promotes secondary and tertiary nucleic acid folding by simultaneously increasing the folding rate while decreasing the unfolding rate, with temperature-dependent studies confirming that the source of PEG stabilization is predominantly entropic and consistent with a true excluded volume crowding mechanism. By way of contrast, neither dextran nor Ficoll induces any significant concentration-dependent change in nucleic acid folding stability at room temperature, but instead, stabilization effects gradually appear with a temperature increase. Such a thermal response indicates that both folding enthalpies and entropies are impacted by dextran and Ficoll. A detailed thermodynamic analysis of the kinetics suggests that, instead of true entropic molecular crowding, dextran and Ficoll associate preferentially with the unfolded vs folded nucleic acid conformer as a result of larger solvent accessible surface area, thereby skewing the free energy landscapes through both significant entropic/enthalpic contributions that compete and fortuitously cancel near room temperature.

The effect of charged residue substitutions on the thermodynamics of protein-surface interactions

The interactions of proteins with surfaces are important in both biological processes and biotechnologies. In contrast to decades of study regarding the biophysics of proteins in bulk solution, however, our mechanistic understanding of the biophysics of proteins interacting with surfaces remains largely qualitative. In response, we have set to explore quantitatively the thermodynamics of protein-surface interactions. In this work, we explore systematically the role of electrostatics in modulating the interaction between proteins and charged surfaces. In particular, we use electrochemistry to explore the extent to which a macroscopic, hydroxyl-coated surface held at a slightly negative potential affects the folding thermodynamics of surface-attached protein variants with different composition of charged amino acids. Doing so, we find that attachment to the surface generally leads to a net stabilization, presumably due to excluded volume effects that reduce the entropy of the unfolded state. The magnitude of this stabilization, however, is strongly correlated with the charged-residue content of the protein. In particular, we find statistically significant correlations with both the net charge of the protein, with greater negative charge leading to less stabilization by the surface, and with the number of arginines, with more arginines leading to greater stabilization. Such findings refine our understanding of protein-surface interactions, providing in turn a guiding rationale to achieve the functional deposition of proteins on artificial surfaces for implementation in, for example, protein-based biotechnologies.

Linking thermodynamics and measurements of protein stability

We review the background, theory and general equations for the analysis of equilibrium protein unfolding experiments, focusing on denaturant and heat-induced unfolding. The primary focus is on the thermodynamics of reversible folding/unfolding transitions and the experimental methods that are available for extracting thermodynamic parameters. We highlight the importance of modelling both how the folding equilibrium depends on a perturbing variable such as temperature or denaturant concentration, and the importance of modelling the baselines in the experimental observables.

Global analysis of protein stability by temperature and chemical denaturation

The stability of a protein is a fundamental property that determines under which conditions, the protein is functional. Equilibrium unfolding with denaturants requires preparation of several samples and only provides the free energy of folding when performed at a single temperature. The typical sample requirement is around 0.5-1 mg of protein. If the stability of many proteins or protein variants needs to be determined, substantial protein production may be needed. Here we have determined the stability of acyl-coenzyme A binding protein at pH 5.3 and chymotrypsin inhibitor 2 at pH 3 and pH 6.25 by combined temperature and denaturant unfolding. We used a setup where tryptophan fluorescence is measured in quartz capillaries where only 10 μl is needed. Temperature unfolding of a series of 15 samples at increasing denaturant concentrations provided accurate and precise thermodynamic parameters. We find that the number of samples may be further reduced and less than 10 μg of protein in total are needed for reliable stability measurements. For assessment of stability of protein purified in small scale e.g. in micro plate format, our method will be highly applicable. The routine for fitting the experimental data is made available as a python notebook.

Dissecting the statistical properties of the linear extrapolation method of determining protein stability

The linear extrapolation method to determine protein stability from denaturant-induced unfolding experiments is based on the observation that the free energy of unfolding is often a linear function of the denaturant concentration. The value in the absence of denaturant is then estimated by extrapolation from this linear relationship. Parameters and their confidence intervals are typically estimated by nonlinear least-squares regression. We have compared different methods for calculating confidence intervals and found that a simple method based on linear theory gives accurate results. We have also compared three different parameterizations of the linear extrapolation method and show that the most commonly used form is problematic since the stability and m-value are correlated in the nonlinear least-squares analysis. Parameter correlation can in some cases causes problems in the estimation of confidence intervals and regions and should be avoided when possible. Two alternative parameterizations show much less correlation between parameters.

CalFitter: a web server for analysis of protein thermal denaturation data

Despite significant advances in the understanding of protein structure-function relationships, revealing protein folding pathways still poses a challenge due to a limited number of relevant experimental tools. Widely-used experimental techniques, such as calorimetry or spectroscopy, critically depend on a proper data analysis. Currently, there are only separate data analysis tools available for each type of experiment with a limited model selection. To address this problem, we have developed the CalFitter web server to be a unified platform for comprehensive data fitting and analysis of protein thermal denaturation data. The server allows simultaneous global data fitting using any combination of input data types and offers 12 protein unfolding pathway models for selection, including irreversible transitions often missing from other tools. The data fitting produces optimal parameter values, their confidence intervals, and statistical information to define unfolding pathways. The server provides an interactive and easy-to-use interface that allows users to directly analyse input datasets and simulate modelled output based on the model parameters. CalFitter web server is available free at https://loschmidt.chemi.muni.cz/calfitter/.

Early Folding Events, Local Interactions, and Conservation of Protein Backbone Rigidity

Protein folding is in its early stages largely determined by the protein sequence and complex local interactions between amino acids, resulting in lower energy conformations that provide the context for further folding into the native state. We compiled a comprehensive data set of early folding residues based on pulsed labeling hydrogen deuterium exchange experiments. These early folding residues have corresponding higher backbone rigidity as predicted by DynaMine from sequence, an effect also present when accounting for the secondary structures in the folded protein. We then show that the amino acids involved in early folding events are not more conserved than others, but rather, early folding fragments and the secondary structure elements they are part of show a clear trend toward conserving a rigid backbone. We therefore propose that backbone rigidity is a fundamental physical feature conserved by proteins that can provide important insights into their folding mechanisms and stability.

Designing Fluorinated Proteins

As methods to incorporate noncanonical amino acid residues into proteins have become more powerful, interest in their use to modify the physical and biological properties of proteins and enzymes has increased. This chapter discusses the use of highly fluorinated analogs of hydrophobic amino acids, for example, hexafluoroleucine, in protein design. In particular, fluorinated residues have proven to be generally effective in increasing the thermodynamic stability of proteins. The chapter provides an overview of the different fluorinated amino acids that have been used in protein design and the various methods available for producing fluorinated proteins. It discusses model proteins systems into which highly fluorinated amino acids have been introduced and the reasons why fluorinated residues are generally stabilizing, with particular reference to thermodynamic and structural studies from our laboratory. Lastly, details of the methodology we have developed to measure the thermodynamic stability of oligomeric fluorinated proteins are presented, as this may be generally applicable to many proteins.

Mechanism of Protein Denaturation: Partial Unfolding of the P22 Coat Protein I-Domain by Urea Binding

The I-domain is an insertion domain of the bacteriophage P22 coat protein that drives rapid folding and accounts for over half of the stability of the full-length protein. We sought to determine the role of hydrogen bonds (H-bonds) in the unfolding of the I-domain by examining (3)JNC' couplings transmitted through H-bonds, the temperature and urea-concentration dependence of (1)HN and (15)N chemical shifts, and native-state hydrogen exchange at urea concentrations where the domain is predominantly folded. The native-state hydrogen-exchange data suggest that the six-stranded β-barrel core of the I-domain is more stable against unfolding than a smaller subdomain comprised of a short α-helix and three-stranded β-sheet. H-bonds, separately determined from solvent protection and (3)JNC' H-bond couplings, are identified with an accuracy of 90% by (1)HN temperature coefficients. The accuracy is improved to 95% when (15)N temperature coefficients are also included. In contrast, the urea dependence of (1)HN and (15)N chemical shifts is unrelated to H-bonding. The protein segments with the largest chemical-shift changes in the presence of urea show curved or sigmoidal titration curves suggestive of direct urea binding. Nuclear Overhauser effects to urea for these segments are also consistent with specific urea-binding sites in the I-domain. Taken together, the results support a mechanism of urea unfolding in which denaturant binds to distinct sites in the I-domain. Disordered segments bind urea more readily than regions in stable secondary structure. The locations of the putative urea-binding sites correlate with the lower stability of the structure against solvent exchange, suggesting that partial unfolding of the structure is related to urea accessibility.

Design, synthesis, and study of fluorinated proteins

Highly fluorinated analogs of hydrophobic amino acids have proven to be generally effective in increasing the thermodynamic stability of proteins. These non-proteogenic amino acids can be incorporated into both α-helix and β-sheet structural motifs and generally enhance protein stability towards unfolding by heat and chemical denaturants, and retard their degradation by proteases. Recent detailed structural and thermodynamic studies have demonstrated that the increase in buried hydrophobic surface area that accompanies fluorination is primarily responsible for the stabilizing properties of fluorinated side chains. Fluorination appears to be a particularly useful strategy for increasing protein stability because fluorinated amino acids closely retain the shape of the side chain, and are thus minimally perturbing to protein structure and function. The first part of this chapter discusses some examples of highly fluorinated model proteins designed by our laboratory and protocols for their synthesis. In the second part, methods for determining their thermodynamic stability, along with conditions that have proven to be useful for crystallizing these proteins, are presented.

OneG: a computational tool for predicting cryptic intermediates in the unfolding kinetics of proteins under native conditions

Understanding the relationships between conformations of proteins and their stabilities is one key to address the protein folding paradigm. The free energy change (ΔG) of unfolding reactions of proteins is measured by traditional denaturation methods and native hydrogen-deuterium (H/D) exchange methods. However, the free energy of unfolding (ΔG(U)) and the free energy of exchange (ΔG(HX)) of proteins are not in good agreement, though the experimental conditions of both methods are well matching to each other. The anomaly is due to any one or combinations of the following reasons: (i) effects of cis-trans proline isomerisation under equilibrium unfolding reactions of proteins (ii) inappropriateness in accounting the baselines of melting curves (iii) presence of cryptic intermediates, which may elude the melting curve analysis and (iv) existence of higher energy metastable states in the H/D exchange reactions of proteins. Herein, we have developed a novel computational tool, OneG, which accounts the discrepancy between ΔG(U) and ΔG(HX) of proteins by systematically accounting all the four factors mentioned above. The program is fully automated and requires four inputs: three-dimensional structures of proteins, ΔG(U), ΔG(U)(*) and residue-specific ΔG(HX) determined under EX2-exchange conditions in the absence of denaturants. The robustness of the program has been validated using experimental data available for proteins such as cytochrome c and apocytochrome b(562) and the data analyses revealed that cryptic intermediates of the proteins detected by the experimental methods and the cryptic intermediates predicted by the OneG for those proteins were in good agreement. Furthermore, using OneG, we have shown possible existence of cryptic intermediates and metastable states in the unfolding pathways of cardiotoxin III and cobrotoxin, respectively, which are homologous proteins. The unique application of the program to map the unfolding pathways of proteins under native conditions have been brought into fore and the program is publicly available at http://sblab.sastra.edu/oneg.html.

Mechanically unfolding protein L using a laser-feedback-controlled cantilever

Force spectroscopy using the atomic force microscope (AFM) can yield important information on the strength and lifetimes of the folded states of single proteins and their complexes when they are loaded with force. For example, by mechanically unfolding concatenated proteins at different velocities, a dynamic force spectrum can be built up that allows reconstruction of the energy landscape that the protein traverses during unfolding. To characterize fully the unfolding landscape, however, it is necessary both to explore the entire force spectrum and to characterize each species populated during unfolding. In the conventional AFM apparatus, force is applied to the protein construct through a compliant cantilever. This limits the dynamic range of the force spectrum that can be probed, and the cantilever recoil after unfolding may mask the presence of metastable intermediates. Here, we describe to our knowledge a new technique-constant-deflection AFM-in which the compliance of the AFM cantilever is removed. Using this technique, we show that protein L exhibits a more complex unfolding energy landscape than previously detected using the conventional technique. This technique is also able to detect the presence of a refolding intermediate whose formation is otherwise prevented by cantilever recoil.

Temperature dependence of acetylcholine receptor channels activated by different agonists

The temperature dependence of agonist binding and channel gating were measured for wild-type adult neuromuscular acetylcholine receptors activated by acetylcholine, carbamylcholine, or choline. With acetylcholine, temperature changed the gating rate constants (Q(10) ≈ 3.2) but had almost no effect on the equilibrium constant. The enthalpy change associated with gating was agonist-dependent, but for all three ligands it was approximately equal to the corresponding free-energy change. The equilibrium dissociation constant of the resting conformation (K(d)), the slope of the rate-equilibrium free-energy relationship (Φ), and the acetylcholine association and dissociation rate constants were approximately temperature-independent. In the mutant αG153S, the choline association and dissociation rate constants were temperature-dependent (Q(10) ≈ 7.4) but K(d) was not. By combining two independent mutations, we were able to compensate for the catalytic effect of temperature on the decay time constant of a synaptic current. At mouse body temperature, the channel-opening and -closing rate constants are ∼400 and 16 ms(-1). We hypothesize that the agonist dependence of the gating enthalpy change is associated with differences in ligand binding, specifically to the open-channel conformation of the protein.

Protein stabilization and the Hofmeister effect: the role of hydrophobic solvation

Using the IGg binding domain of protein L from Streptoccocal magnus (ProtL) as a case study, we investigated how the anions of the Hofmeister series affect protein stability. To that end, a suite of lysine-to-glutamine modifications were obtained and structurally and thermodynamically characterized. The changes in stability introduced with the mutation are related to the solvent-accessible area of the side chain, specifically to the solvation of the nonpolar moiety of the residue. The thermostability for the set of ProtL mutants was determined in the presence of varying concentrations (0-1 M) of six sodium salts from the Hofmeister series: sulfate, phosphate, fluoride, nitrate, perchlorate, and thiocyanate. For kosmotropic anions (sulfate, phosphate, and fluoride), the stability changes induced by the cosolute (encoded in m(3)=deltaDeltaG(0)/deltaC(3)) are proportional to the surface changes introduced with the mutation. In contrast, the m(3) values measured for chaotropic anions are much more independent of such surface modifications. Our results are consistent with a model in which the increase in the solution surface tension induced by the anion stabilizes the folded conformation of the protein. This contribution complements the nonspecific and weak interactions between the ions and the protein backbone that shift the equilibrium toward the unfolded state.

How accurate are polymer models in the analysis of Förster resonance energy transfer experiments on proteins?

Single molecule Förster resonance energy transfer (FRET) experiments are used to infer the properties of the denatured state ensemble (DSE) of proteins. From the measured average FRET efficiency, <E>, the distance distribution P(R) is inferred by assuming that the DSE can be described as a polymer. The single parameter in the appropriate polymer model (Gaussian chain, wormlike chain, or self-avoiding walk) for P(R) is determined by equating the calculated and measured <E>. In order to assess the accuracy of this "standard procedure," we consider the generalized Rouse model (GRM), whose properties [<E> and P(R)] can be analytically computed, and the Molecular Transfer Model for protein L for which accurate simulations can be carried out as a function of guanadinium hydrochloride (GdmCl) concentration. Using the precisely computed <E> for the GRM and protein L, we infer P(R) using the standard procedure. We find that the mean end-to-end distance can be accurately inferred (less than 10% relative error) using <E> and polymer models for P(R). However, the value extracted for the radius of gyration (R(g)) and the persistence length (l(p)) are less accurate. For protein L, the errors in the inferred properties increase as the GdmCl concentration increases for all polymer models. The relative error in the inferred R(g) and l(p), with respect to the exact values, can be as large as 25% at the highest GdmCl concentration. We propose a self-consistency test, requiring measurements of <E> by attaching dyes to different residues in the protein, to assess the validity of describing DSE using the Gaussian model. Application of the self-consistency test to the GRM shows that even for this simple model, which exhibits an order-->disorder transition, the Gaussian P(R) is inadequate. Analysis of experimental data of FRET efficiencies with dyes at several locations for the cold shock protein, and simulations results for protein L, for which accurate FRET efficiencies between various locations were computed, shows that at high GdmCl concentrations there are significant deviations in the DSE P(R) from the Gaussian model.

Mechanical characterization of protein L in the low-force regime by electromagnetic tweezers/evanescent nanometry

Mechanical manipulation at the single molecule level of proteins exhibiting mechanical stability poses a technical challenge that has been almost exclusively approached by atomic force microscopy (AFM) techniques. However, due to mechanical drift limitations, AFM techniques are restricted to experimental recordings that last less than a minute in the high-force regime. Here we demonstrate a novel combination of electromagnetic tweezers and evanescent nanometry that readily captures the forced unfolding trajectories of protein L at pulling forces as low as 10-15 pN. Using this approach, we monitor unfolding and refolding cycles of the same polyprotein for a period of time longer than 30 min. From such long-lasting recordings, we obtain ensemble averages of unfolding step sizes and rates that are consistent with single-molecule AFM data obtained at higher stretching forces. The unfolding kinetics of protein L at low stretching forces confirms and extends the observations that the mechanical unfolding rate is exponentially dependent on the pulling force within a wide range of stretching forces spanning from 13 pN up to 120 pN. Our experiments demonstrate a novel approach for the mechanical manipulation of single proteins for extended periods of time in the low-force regime.

Molecular origin of constant m-values, denatured state collapse, and residue-dependent transition midpoints in globular proteins

Experiments show that for many two-state folders the free energy of the native state, DeltaG(ND)([C]), changes linearly as the denaturant concentration, [C], is varied. The slope {m = [dDeltaG(ND)([C])]/(d[C])}, is nearly constant. According to the transfer model, the m-value is associated with the difference in the surface area between the native (N) and denatured (D) state, which should be a function of DeltaR(g)(2), the difference in the square of the radius of gyration between the D and N states. Single-molecule experiments show that the R(g) of the structurally heterogeneous denatured state undergoes an equilibrium collapse transition as [C] decreases, which implies m also should be [C]-dependent. We resolve the conundrum between constant m-values and [C]-dependent changes in R(g) using molecular simulations of a coarse-grained representation of protein L, and the molecular transfer model, for which the equilibrium folding can be accurately calculated as a function of denaturant (urea) concentration. In agreement with experiment, we find that over a large range of denaturant concentration (>3 M) the m-value is a constant, whereas under strongly renaturing conditions (<3 M), it depends on [C]. The m-value is a constant above [C] > 3 M because the [C]-dependent changes in the surface area of the backbone groups, which make the largest contribution to m, are relatively narrow in the denatured state. The burial of the backbone and hydrophobic side chains gives rise to substantial surface area changes below [C] < 3 M, leading to collapse in the denatured state of protein L. Dissection of the contribution of various amino acids to the total surface area change with [C] shows that both the sequence context and residual structure are important. There are [C]-dependent variations in the surface area for chemically identical groups such as the backbone or Ala. Consequently, the midpoints of transition of individual residues vary significantly (which we call the Holtzer effect) even though global folding can be described as an all-or-none transition. The collapse is specific in nature, resulting in the formation of compact structures with appreciable populations of nativelike secondary structural elements. The collapse transition is driven by the loss of favorable residue-solvent interactions and a concomitant increase in the strength of intrapeptide interactions with a decreasing [C]. The strength of these interactions is nonuniformly distributed throughout the structure of protein L. Certain secondary structure elements have stronger [C]-dependent interactions than others in the denatured state.

Ruggedness in the folding landscape of protein L

By exploring the folding pathways of the B1 domain of protein L with a series of equilibrium and rapid kinetic experiments, we have found its unfolded state to be more complex than suggested by two-state folding models. Using an ultrarapid mixer to initiate protein folding within approximately 2-4 microseconds, we observe folding kinetics by intrinsic tryptophan fluorescence and fluorescence resonance energy transfer. We detect at least two processes faster than 100 mus that would be hidden within the burst phase of a stopped-flow instrument measuring tryptophan fluorescence. Previously reported measurements of slow intramolecular diffusion are commensurate with the slower of the two observed fast phases. These results suggest that a multidimensional energy landscape is necessary to describe the folding of protein L, and that the dynamics of the unfolded state is dominated by multiple small energy barriers.

Effects of denaturants and osmolytes on proteins are accurately predicted by the molecular transfer model

Interactions between denaturants and proteins are commonly used to probe the structures of the denatured state ensemble and their stabilities. Osmolytes, a class of small intracellular organic molecules found in all taxa, also profoundly affect the equilibrium properties of proteins. We introduce the molecular transfer model, which combines simulations in the absence of denaturants or osmolytes, and Tanford's transfer model to predict the dependence of equilibrium properties of proteins at finite concentration of osmolytes. The calculated changes in the thermodynamic quantities (probability of being in the native basin of attraction, m values, FRET efficiency, and structures of the denatured state ensemble) with GdmCl concentration [C] for the protein L and cold shock protein CspTm compare well with experiments. The radii of gyration of the subpopulation of unfolded molecules for both proteins decrease (i.e., they undergo a collapse transition) as [C] decreases. Although global folding is cooperative, residual secondary structures persist at high denaturant concentrations. The temperature dependence of the specific heat shows that the folding temperature (T(F)) changes linearly as urea and trimethylamine N-oxide (TMAO) concentrations increase. The increase in T(F) in TMAO can be as large as 20 degrees C, whereas urea decreases T(F) by as much as 35 degrees C. The stabilities of protein L and CspTm also increase linearly with the concentration of osmolytes (proline, sorbitol, sucrose, TMAO, and sarcosine).

Anion modulation of the 1H/2H exchange rates in backbone amide protons monitored by NMR spectroscopy

The ability of three anionic cosolutes (sulfate, thiocyanate, and chloride) in modulating the (1)H/(2)H exchange rates for backbone amide protons has been investigated using nuclear magnetic resonance (NMR) for two different proteins: the IGg-binding domain of protein L (ProtL) and the glucose-galactose-binding protein (GGBP). Our results show that moderate anion concentrations (0.2 M-1 M) regulate the exchange rate following the Hofmeister series: Addition of thiocyanate increases the exchange rates for both proteins, while sulfate and chloride (to a less extent) slow down the exchange reaction. In the presence of the salt, no alteration of the protein structure and minimal variations in the number of measurable peaks are observed. Experiments with model compounds revealed that the unfolded state is modulated in an equivalent way by these cosolutes. For ProtL, the estimated values for the local free energy change upon salt addition (m (3,DeltaG )) are consistent with the previously reported free energy contribution from the cosolute's preferential interaction/exclusion term indicating that nonspecific weak interactions between the anion and the amide groups constitute the dominant mechanism for the exchange-rate modulation. The same trend is also found for GGBP in the presence of thiocyanate, underlining the generality of the exchange-rate modulation mechanism, complementary to more investigated effects like the electrostatic interactions or specific anion binding to protein sites.

Partially folded states of staphylococcal nuclease highlight the conserved structural hierarchy of OB-fold proteins

We have been interested in whether three proteins that share a five-stranded beta-barrel "OB-fold" structural motif but no detectable sequence homology fold by similar mechanisms. Here we describe native-state hydrogen exchange experiments as a function of urea for SN (staphylococcal nuclease), a protein with an OB-fold motif and additional nonconserved elements of structure. The regions of structure with the largest stability and unfolding cooperativity are contained within the conserved OB-fold portion of SN, consistent with previous results for CspA (cold shock protein A) and LysN (anticodon binding domain of lysyl tRNA synthetase). The OB-fold also has the subset of residues with the slowest unfolding rates in the three proteins, as determined by hydrogen exchange experiments in the EX1 limit. Although the protein folding hierarchy is maintained at the level of supersecondary structure, it is not evident for individual residues as might be expected if folding depended on obligatory nucleation sites. Rather, the site-specific stability profiles appear to be linked to sequence hydrophobicity and to the density of long-range contacts at each site in the three-dimensional structures of the proteins. We discuss the implications of the correlation between stability to unfolding and conservation of structure for mechanisms of protein structure evolution.

Native state energetics of the Src SH2 domain: evidence for a partially structured state in the denatured ensemble

We have defined the free-energy profile of the Src SH2 domain using a variety of biophysical techniques. Equilibrium and kinetic experiments monitored by tryptophan fluorescence show that Src SH2 is quite stable and folds rapidly by a two-state mechanism, without populating any intermediates. Native state hydrogen-deuterium exchange confirms this two-state behavior; we detect no cooperative partially unfolded forms in equilibrium with the native conformation under any conditions. Interestingly, the apparent stability of the protein from hydrogen exchange is 2 kcal/mol greater than the stability determined by both equilibrium and kinetic studies followed by fluorescence. Native-state proteolysis demonstrates that this "super protection" does not result from a deviation from the linear extrapolation model used to fit the fluorescence data. Instead, it likely arises from a notable compaction in the unfolded state under native conditions, resulting in an ensemble of conformations with substantial solvent exposure of side chains and flexible regions sensitive to proteolysis, but backbone amides that exchange with solvent approximately 30-fold slower than would be expected for a random coil. The apparently simple behavior of Src SH2 in traditional unfolding studies masks the significant complexity present in the denatured-state ensemble.

An Experimental Investigation of Conformational Fluctuations in Proteins G and L

The B1 domains of streptococcal proteins G and L are structurally similar, but they have different sequences and they fold differently. We have measured their NMR spectra at variable temperature using a range of concentrations of denaturant. Many residues have curved amide proton temperature dependence, indicating that they significantly populate alternative, locally unfolded conformations. The results, therefore, provide a view of the locations of low-lying, locally unfolded conformations. They indicate approximately 4-6 local minima for each protein, all within ca. 2.5 kcal/mol of the native state, implying a locally rough energy landscape. Comparison with folding data for these proteins shows that folding involves most molecules traversing a similar path, once a transition state containing a beta hairpin has been formed, thereby defining a well-populated pathway down the folding funnel. The hairpin that directs the folding pathway differs for the two proteins and remains the most stable part of the folded protein.

Mechanically unfolding the small, topologically simple protein L

beta-sheet proteins are generally more able to resist mechanical deformation than alpha-helical proteins. Experiments measuring the mechanical resistance of beta-sheet proteins extended by their termini led to the hypothesis that parallel, directly hydrogen-bonded terminal beta-strands provide the greatest mechanical strength. Here we test this hypothesis by measuring the mechanical properties of protein L, a domain with a topology predicted to be mechanically strong, but with no known mechanical function. A pentamer of this small, topologically simple protein is resistant to mechanical deformation over a wide range of extension rates. Molecular dynamics simulations show the energy landscape for protein L is highly restricted for mechanical unfolding and that this protein unfolds by the shearing apart of two structural units in a mechanism similar to that proposed for ubiquitin, which belongs to the same structural class as protein L, but unfolds at a significantly higher force. These data suggest that the mechanism of mechanical unfolding is conserved in proteins within the same fold family and demonstrate that although the topology and presence of a hydrogen-bonded clamp are of central importance in determining mechanical strength, hydrophobic interactions also play an important role in modulating the mechanical resistance of these similar proteins.

A large scale test of computational protein design: folding and stability of nine completely redesigned globular proteins

A previously developed computer program for protein design, RosettaDesign, was used to predict low free energy sequences for nine naturally occurring protein backbones. RosettaDesign had no knowledge of the naturally occurring sequences and on average 65% of the residues in the designed sequences differ from wild-type. Synthetic genes for ten completely redesigned proteins were generated, and the proteins were expressed, purified, and then characterized using circular dichroism, chemical and temperature denaturation and NMR experiments. Although high-resolution structures have not yet been determined, eight of these proteins appear to be folded and their circular dichroism spectra are similar to those of their wild-type counterparts. Six of the proteins have stabilities equal to or up to 7kcal/mol greater than their wild-type counterparts, and four of the proteins have NMR spectra consistent with a well-packed, rigid structure. These encouraging results indicate that the computational protein design methods can, with significant reliability, identify amino acid sequences compatible with a target protein backbone.

Folding subdomains of thioredoxin characterized by native-state hydrogen exchange

Native-state hydrogen exchange (HX) studies, used in conjunction with NMR spectroscopy, have been carried out on Escherichia coli thioredoxin (Trx) for characterizing two folding subdomains of the protein. The backbone amide protons of only the slowest-exchanging 24 amino acid residues, of a total of 108 amino acid residues, could be followed at pH 7. The free energy of the opening event that results in an amide hydrogen exchanging with solvent (DeltaG(op)) was determined at each of the 24 amide hydrogen sites. The values of DeltaG(op) for the amide hydrogens belonging to residues in the helices alpha(1), alpha(2), and alpha(4) are consistent with them exchanging with the solvent only when the fully unfolded state is sampled transiently under native conditions. The denaturant-dependences of the values of DeltaG(op) provide very little evidence that the protein samples partially unfolded forms, lower in energy than the unfolded state. The amide hydrogens belonging to the residues in the beta strands, which form the core of the protein, appear to have higher values of DeltaG(op) than amide hydrogens belonging to residues in the helices, suggesting that they might be more stable to exchange. This apparently higher stability to HX of the beta strands might be either because they exchange out their amide hydrogens in a high energy intermediate preceding the globally unfolded state, or, more likely, because they form residual structure in the globally unfolded state. In either case, the central beta strands-beta(3,) beta(2), and beta(4)-would appear to form a cooperatively folding subunit of the protein. The native-state HX methodology has made it possible to characterize the free energy landscape that Trx can sample under equilibrium native conditions.

The origins of asymmetry in the folding transition states of protein L and protein G

Topology has been shown to be an important determinant of many features of protein folding; however, the delineation of sequence effects on folding remains obscure. Furthermore, differentiation between the two influences proves difficult due to their intimate relationship. To investigate the effect of sequence in the absence of significant topological differences, we examined the folding mechanisms of segment B1 peptostreptococcal protein L and segment B1 of streptococcal protein G. These proteins share the same highly symmetrical topology. Despite this symmetry, neither protein folds through a symmetrical transition state. We analyzed the origins of this difference using theoretical models. We found that the strength of the interactions present in the N-terminal hairpin of protein L causes this hairpin to form ahead of the C-terminal hairpin. The difference in chain entropy associated with the formation of the hairpins of protein G proves sufficient to beget initiation of folding at the shorter C-terminal hairpin. Our findings suggest that the mechanism of folding may be understood by examination of the free energy associated with the formation of partially folded microstates.

Thermodynamic stability measurements on multimeric proteins using a new H/D exchange- and matrix-assisted laser desorption/ionization (MALDI) mass spectrometry-based method

We recently reported on a new H/D exchange- and matrix-assisted laser desorption/ionization (MALDI) mass spectrometry-based technique, termed SUPREX, that removes several important limitations associated with measuring the thermodynamic stability of proteins. In contrast to conventional spectroscopy-based techniques for characterizing the equilibrium unfolding behavior of proteins, SUPREX is amenable to the thermodynamic analysis of both purified and unpurified proteins using mg to ng quantities of material. Here we report on the application of SUPREX to the analysis of multimeric protein systems. Included in this work are the SUPREX results we obtained in studies on six model multimeric proteins including the GCN4p1 dimer, the coil-V(a)L(d) trimer, the 4-oxalocrotonate tautomerase (4-OT) hexamer, the Trp repressor (TrpR) dimer, the Arc repressor (ArcR) dimer, and an ArcR mutant (the (DOA20)ArcR) dimer which contained two destabilizing mutations including an Asp to Ala mutation at position 20 and an amide to ester bond mutation between amino acid (aa) residues 19 and 20. As part of the work described here, we present a new method for the analysis of SUPREX data that is generally applicable to both monomeric and multimeric protein systems. Our results on the model proteins in this study indicate that this new method can be used to determine folding free energies for proteins with the accuracy and precision of conventional spectroscopy-based methods.

Absence of stable intermediates on the folding pathway of barnase

Barnase is one of the few protein models that has been studied extensively for protein folding. Previous studies led to the conclusion that barnase folds through a very stable submillisecond intermediate ( approximately 3 kcal/mol). The structure of this intermediate was characterized intensively by using a protein engineering approach. This intermediate has now been reexamined with three direct and independent methods. (i) Hydrogen exchange experiments show very small protection factors ( approximately 2) for the putative intermediate, indicating a stability of approximately 0.0 kcal/mol. (ii) Denaturant-dependent unfolding of the putative intermediate is noncooperative and indicates a stability less than 0.0 kcal/mol. (iii) The logarithm of the unfolding rate constant of native barnase vs. denaturant concentrations is not linear. Together with the measured rate ("I" to N), this nonlinear behavior accounts for almost all of the protein stability, leaving only about 0.3 kcal/mol that could be attributed to the rapidly formed intermediate. Other observations previously interpreted to support the presence of an intermediate are now known to have alternative explanations. These results cast doubts on the previous conclusions on the nature of the early folding state in barnase and therefore should have important implications in understanding the early folding events of barnase and other proteins in general.

Ensemble modulation as an origin of denaturant-independent hydrogen exchange in proteins

Native state hydrogen exchange (HX) has become a powerful tool for the analysis of conformational states that exist under native conditions. However, the interpretation of HX data in terms of conformational fluctuations is still controversial. In particular, it has been shown that many residues display exchange behavior that is independent of denaturant concentration. It has been postulated that this lack of denaturant dependence results from local fluctuations that do not expose appreciable amounts of buried surface area. Here, we use a general thermodynamic description of HX to explore the different possibilities for this behavior. We find that the denaturant dependence seen in HX experiments under native conditions is not a de facto indication of the amount of surface area exposure required for exchange. Instead, this behavior results from the relatively homogenous character of the conformational ensemble that exists under native conditions and the non-specific nature of denaturant effects. Furthermore, a comparison of the HX behavior from a stabilized mutant of Staphylococcal nuclease (SNase) with that predicted for the wild-type SNase from the COREX algorithm suggests that denaturant-independent exchange of many residues is consistent with significant (approximately 10 %) surface area exposure for this protein.

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.

NMR characterization of residual structure in the denatured state of protein L

Triple-resonance NMR experiments were used to assign the (13)C(alpha), (13)C(beta), (15)N and NH resonances for all the residues in the denatured state of a destabilized protein L variant in 2 M guanidine. The chemical shifts of most resonances were very close to their random coil values. Significant deviations were observed for G22, L38 and K39; increasing the denaturant concentration shifted the chemical shifts of these residues towards theory random coil values. Medium-range nuclear Overhauser enhancements were detected in segments corresponding to the turn between the first two strands, the end of the second strand through the turn between the second strand and the helix, and the turn between the helix and the third strand in 3D H(1), N(15)-HSQC-NOESY-HSQC experiments on perdeuterated samples. Longer-range interactions were probed by measuring the paramagnetic relaxation enhancement produced by nitroxide spin labels introduced via cysteine residues at five sites around the molecule. Damped oscillations in the magnitude of the paramagnetic relaxation enhancement as a function of distance along the sequence suggested native-like chain reversals in the same three turn regions. The more extensive interactions within the region corresponding to the first beta-turn than in the region corresponding to the second beta-turn suggests that the asymmetry in the folding reaction evident in previous studies of the protein L folding transition state is already established in the denatured state.

Microscopic stability of cold shock protein A examined by NMR native state hydrogen exchange as a function of urea and trimethylamine N-oxide

Native state hydrogen exchange of cold shock protein A (CspA) has been characterized as a function of the denaturant urea and of the stabilizing agent trimethylamine N-oxide (TMAO). The structure of CspA has five strands of beta-sheet. Strands beta1-beta4 have strongly protected amide protons that, based on experiments as a function of urea, exchange through a simple all-or-none global unfolding mechanism. By contrast, the protection of amide protons from strand beta5 is too weak to measure in water. Strand beta5 is hydrogen bonded to strands beta3 and beta4, both of which afford strong protection from solvent exchange. Gaussian network model (GNM) simulations, which assume that the degree of protection depends on tertiary contact density in the native structure, accurately predict the strong protection observed in strands beta1-beta4 but fail to account for the weak protection in strand beta5. The most conspicuous feature of strand beta5 is its low sequence hydrophobicity. In the presence of TMAO, there is an increase in the protection of strands beta1-beta4, and protection extends to amide protons in more hydrophilic segments of the protein, including strand beta5 and the loops connecting the beta-strands. TMAO stabilizes proteins by raising the free energy of the denatured state, due to highly unfavorable interactions between TMAO and the exposed peptide backbone. As such, the stabilizing effects of TMAO are expected to be relatively independent of sequence hydrophobicity. The present results suggest that the magnitude of solvent exchange protection depends more on solvent accessibility in the ensemble of exchange susceptible conformations than on the strength of hydrogen-bonding interactions in the native structure.

The hydrogen exchange core and protein folding

A database of hydrogen-deuterium exchange results has been compiled for proteins for which there are published rates of out-exchange in the native state, protection against exchange during folding, and out-exchange in partially folded forms. The question of whether the slow exchange core is the folding core (Woodward C, 1993, Trends Biochem Sci 18:359-360) is reexamined in a detailed comparison of the specific amide protons (NHs) and the elements of secondary structure on which they are located. For each pulsed exchange or competition experiment, probe NHs are shown explicitly; the large number and broad distribution of probe NHs support the validity of comparing out-exchange with pulsed-exchange/competition experiments. There is a strong tendency for the same elements of secondary structure to carry NHs most protected in the native state, NHs first protected during folding, and NHs most protected in partially folded species. There is not a one-to-one correspondence of individual NHs. Proteins for which there are published data for native state out-exchange and theta values are also reviewed. The elements of secondary structure containing the slowest exchanging NHs in native proteins tend to contain side chains with high theta values or be connected to a turn/loop with high theta values. A definition for a protein core is proposed, and the implications for protein folding are discussed. Apparently, during folding and in the native state, nonlocal interactions between core sequences are favored more than other possible nonlocal interactions. Other studies of partially folded bovine pancreatic trypsin inhibitor (Barbar E, Barany G, Woodward C, 1995, Biochemistry 34:11423-11434; Barber E, Hare M, Daragan V, Barany G, Woodward C, 1998, Biochemistry 37:7822-7833), suggest that developing cores have site-specific energy barriers between microstates, one disordered, and the other(s) more ordered.

NMR hydrogen exchange of the OB-fold protein LysN as a function of denaturant: the most conserved elements of structure are the most stable to unfolding

The structure of LysN contains an OB-fold motif composed of a structurally conserved five-stranded beta-barrel capped by a poorly conserved alpha-helix between strands beta3 and beta4. Two additional alpha-helices, unique to the LysN structure, flank the N terminus of the OB-fold. The stability of LysN to unfolding has been investigated with NMR native state hydrogen exchange measurements as a function of guanidinium hydrochloride concentration, and equilibrium unfolding transitions monitored by ellipticity at 222 nm and fluorescence at 350 nm. The spectrophotometric measurements suggest an apparent two-state unfolding transition with DeltaGu(0) approximately 6 kcal/mol and m approximately 3 kcal/(molM). By contrast, NMR hydrogen exchange measurements manifest a distribution of DeltaGu(0) and m values which indicate that the protein can undergo subglobal unfolding. The largest DeltaGu(0) values from hydrogen exchange are for residues in the beta-sheet of the protein. These values, which reflect complete unfolding of the protein, are between 3 and 4 kcal/mol higher than those obtained from circular dichroism or fluorescence. This discrepancy may be due to the comparison of NMR hydrogen exchange parameters measured at residue-level resolution, with spectrophotometric parameters that reflect an unresolved super position of unfolding transitions of the alpha-helices and beta-strands. The largest DeltaGu(0) values obtained from hydrogen exchange for the subset of residues in the alpha-helices of the protein, agree with the DeltaGu(0) values obtained from circular dichroism or fluorescence. Based on the hydrogen exchange data, however, the three alpha-helices of LysN are on average 3 kcal/mol less stable than the beta-sheet. Consistent with the subglobal unfolding of LysN evinced by hydrogen exchange, a deletion mutant that lacks the first alpha-helix of the protein retains a cooperatively folded structure. Taken together with previous results on the OB-fold proteins SN and CspA, the present results for LysN suggest that the most conserved elements of structure in the OB-fold motif are the most resistant to denaturation. In all three proteins, stability to denaturation correlates with sequence hydrophobicity.

The single helix in protein L is largely disrupted at the rate-limiting step in folding

To investigate the role of helix formation in the folding of protein L, a 62 residue alpha/beta protein, we studied the consequences of both single and multiple mutations in the helix on the kinetics of folding. A triple mutant with 11 additional carbon atoms in core residues in the amino-terminal portion of the helix folded substantially faster than wild type, suggesting that hydrophobic association with residues elsewhere in the protein occurs at the rate-limiting step in folding. However, helix-destabilizing mutations had little effect on the rate of folding; in particular, a triple glycine substitution on the solvent-exposed side of the helix increased the unfolding rate 56-fold while reducing the folding rate less than threefold. Thus, in contrast to the predictions of models of folding involving the coalescence of well-formed secondary structure elements, the single helix in protein L appears to be largely disrupted at the rate-limiting step in folding and unfolding.

Limited internal friction in the rate-limiting step of a two-state protein folding reaction

Small, single-domain proteins typically fold via a compact transition-state ensemble in a process well fitted by a simple, two-state model. To characterize the rate-limiting conformational changes that underlie two-state folding, we have investigated experimentally the effects of changing solvent viscosity on the refolding of the IgG binding domain of protein L. In conjunction with numerical simulations, our results indicate that the rate-limiting conformational changes of the folding of this domain are strongly coupled to solvent viscosity and lack any significant "internal friction" arising from intrachain collisions. When compared with the previously determined solvent viscosity dependencies of other, more restricted conformational changes, our results suggest that the rate-limiting folding transition involves conformational fluctuations that displace considerable amounts of solvent. Reconciling evidence that the folding transition state ensemble is comprised of highly collapsed species with these and similar, previously reported results should provide a significant constraint for theoretical models of the folding process.

Global analysis of the thermal and chemical denaturation of the N-terminal domain of the ribosomal protein L9 in H2O and D2O. Determination of the thermodynamic parameters, deltaH(o), deltaS(o), and deltaC(o)p and evaluation of solvent isotope effects

The stability of the N-terminal domain of the ribosomal protein L9, NTL9, from Bacillus stearothermophilus has been monitored by circular dichroism at various temperatures and chemical denaturant concentrations in H2O and D2O. The basic thermodynamic parameters for the unfolding reaction, deltaH(o), deltaS(o), and deltaC(o)p, were determined by global analysis of temperature and denaturant effects on stability. The data were well fit by a model that assumes stability varies linearly with denaturant concentration and that uses the Gibbs-Helmholtz equation to model changes in stability with temperature. The results obtained from the global analysis are consistent with information obtained from individual thermal and chemical denaturations. NTL9 has a maximum stability of 3.78 +/- 0.25 kcal mol(-1) at 14 degrees C. DeltaH(o)(25 degrees C) for protein unfolding equals 9.9 +/- 0.7 kcal mol(-1) and TdeltaS(o)++(25 degrees C) equals 6.2 +/- 0.6 kcal mol(-1). DeltaC(o)p equals 0.53 +/- 0.06 kcal mol(-1) deg(-1). There is a small increase in stability when D2O is substituted for H2O. Based on the results from global analysis, NTL9 is 1.06 +/- 0.60 kcal mol(-1) more stable in D2O at 25 degrees C and Tm is increased by 5.8 +/- 3.6 degrees C in D2O. Based on the results from individual denaturation experiments, NTL9 is 0.68 +/- 0.68 kcal mol(-1) more stable in D2O at 25 degrees C and Tm is increased by 3.5 +/- 2.1 degrees C in D2O. Within experimental error there are no changes in deltaH(o) (25 degrees C) when D2O is substituted for H2O.

Contact order, transition state placement and the refolding rates of single domain proteins

Theoretical studies have suggested relationships between the size, stability and topology of a protein fold and the rate and mechanisms by which it is achieved. The recent characterization of the refolding of a number of simple, single domain proteins has provided a means of testing these assertions. Our investigations have revealed statistically significant correlations between the average sequence separation between contacting residues in the native state and the rate and transition state placement of folding for a non-homologous set of simple, single domain proteins. These indicate that proteins featuring primarily sequence-local contacts tend to fold more rapidly and exhibit less compact folding transition states than those characterized by more non-local interactions. No significant relationship is apparent between protein length and folding rates, but a weak correlation is observed between length and the fraction of solvent-exposed surface area buried in the transition state. Anticipated strong relationships between equilibrium folding free energy and folding kinetics, or between chemical denaturant and temperature dependence-derived measures of transition state placement, are not apparent. The observed correlations are consistent with a model of protein folding in which the size and stability of the polypeptide segments organized in the transition state are largely independent of protein length, but are related to the topological complexity of the native state. The correlation between topological complexity and folding rates may reflect chain entropy contributions to the folding barrier.

Hydrogen exchange and protein folding

Amide hydrogen-deuterium exchange is a sensitive probe of the structure, stability and dynamics of proteins. The significant increase in the number of small, model proteins that have been studied has allowed a better understanding of the structural fluctuations that lead to hydrogen exchange. Recent technical advances enable the methodology to be applied to the study of protein-protein interactions in much larger, more complex systems.

Contrasting roles for symmetrically disposed beta-turns in the folding of a small protein

To investigate the role of turns in protein folding, we have characterized the effects of combinatorial and site-directed mutations in the two beta-turns of peptostreptococcal protein L on folding thermodynamics and kinetics. Sequences of folded variants recovered from combinatorial libraries using a phase display selection method were considerably more variable in the second turn than in the first turn. These combinatorial mutants as well as strategically placed point mutants in the two turns had a similar range of thermodynamic stabilities, but strikingly different folding kinetics. A glycine to alanine substitution in the second beta-turn increased the rate of unfolding more than tenfold but had little effect on the rate of folding, while mutation of a symmetrically disposed glycine residue in the first turn had little effect on unfolding but slowed the rate of folding nearly tenfold. These results demonstrate that the role of beta-turns in protein folding is strongly context-dependent, and suggests that the first turn is formed and the second turn disrupted in the folding transition state.