To use all functions of this page, please activate cookies in your browser.
my.chemeurope.com
With an accout for my.chemeurope.com you can always see everything at a glance – and you can configure your own website and individual newsletter.
 My watch list
 My saved searches
 My saved topics
 My newsletter
Equilibrium unfolding
In biochemistry, equilibrium unfolding is the process of unfolding a protein or RNA molecule by gradually changing its solution conditions, i.e., its environment. Since equilibrium is maintained at all steps, the process is reversible (equilibrium folding). Equilibrium unfolding is used to determine the conformational stability of the molecule. Additional recommended knowledge
Theoretical backgroundIn its simplest form, equilibrium unfolding assumes that the molecule may belong to only two thermodynamic states, the folded state (typically denoted N for "native" state) and the unfolded state (typically denoted U). This "allornone" model of protein folding was first proposed by Tim Anson (1945), but is believed to hold only for small, single structural domains of proteins (Jackson, 1998); larger domains and multidomain proteins often exhibit intermediate states. As usual in statistical mechanics, these states correspond to ensembles of molecular conformations, not just one conformation. The molecule may transition between the native and unfolded states according to a simple kinetic model with rate constants k_{f} and k_{u} for the folding () and unfolding () reactions, respectively. The dimensionless equilibrium constant can be used to determine the conformational stability ΔG by the equation
where R is the gas constant and T is the absolute temperature in kelvins. Thus, ΔG is positive if the unfolded state is less stable (i.e., disfavored) relative to the native state. The most direct way to measure the conformational stability ΔG of a molecule with twostate folding is to measure its kinetic rate constants k_{f} and k_{u} under the solution conditions of interest. However, since protein folding is typically completed in milliseconds, such measurements can be difficult to perform, usually requiring expensive stoppedflow or (more recently) continuousflow mixers to provoke folding with a high time resolution. Chemical denaturationIn the less expensive technique of equilibrium unfolding, the fractions of folded and unfolded molecules (denoted as p_{N} and p_{U}, respectively) are measured as the solution conditions are gradually changed from those favoring the native state to those favoring the unfolded state, e.g., by adding a denaturant such as guanidinium hydrochloride or urea. (In equilibrium folding, the reverse process is carried out.) Given that the fractions must sum to one and their ratio must be given by the Boltzmann factor, we have Protein stabilities are typically found to vary linearly with the denaturant concentration. A number of models have been proposed to explain this observation prominent among them being the denaturant binding model, solventexchange model (both by Schellman, JA^{[1]}) and the Linear Energy Model (LEM; Pace, NC^{[2]}). All of the models assume that only two thermodynamic states are populated/depopulated upon denaturation. They could be extended to interpret more complicated reaction schemes. The denaturant binding model assumes that there are specific but independent sites on the protein molecule (folded or unfolded) to which the denaturant binds with an effective (average) binding constant k. The equilibrium shifts towards the unfolded state at high denaturant concentrations as it has more binding sites for the denaturant relative to the folded state (Δn). In other words, the increased number of potential sites exposed in the unfolded state is seen as the reason for denaturation transitions. An elementary treatment results in the following functional form: where ΔG_{w} is the stability of the protein in water and [D] is the denaturant concentration. Thus the analysis of denaturation data with this model requires 7 parameters: ΔG_{w},Δn, k, and the slopes and intercepts of the folded and unfolded state baselines. The solvent exchange model (also called the ‘weak binding model’ or ‘selective solvation’) of Schellman invokes the idea of an equilibrium between the water molecules bound to independent sites on protein and the denaturant molecules in solution. It has the form: where K is the equilibrium constant for the exchange reaction and X_{d} is the molefraction of the denaturant in solution. This model tries to answer the question of whether the denaturant molecules actually bind to the protein or they seem to be bound just because denaturants occupy about 2030 % of the total solution volume at high concentrations used in experiments, i.e. nonspecific effects – and hence the term ‘weak binding’. As in the denaturantbinding model, fitting to this model also requires 7 parameters. One common theme obtained from both these models is that the binding constants (in the molar scale) for urea and guanidinium hydrochloride are small: ~ 0.2 M^{ − 1} for urea and 0.6 M^{ − 1} for GuHCl. Intuitively, the difference in the number of binding sites between the folded and unfolded states is directly proportional to the differences in the accessible surface area. This forms the basis for the LEM which assumes a simple linear dependence of stability on the denaturant concentration. The resulting slope of the plot of stability versus the denaturant concentration is called the mvalue. In pure mathematical terms, mvalue is the derivative of the change in stabilization free energy upon the addition of denaturant. However, a strong correlation between the accessible surface area (ASA) exposed upon unfolding, i.e. difference in the ASA between the unfolded and folded state of the studied protein (dASA), and the mvalue has been documented by Pace and coworkers. In view of this observation, the mvalues are typically interpreted as being proportional to the dASA. There is no physical basis for the LEM and is purely empirical, though it is widely used in interpreting solventdenaturation data. It has the general form: where the slope m is called the "mvalue"(> 0 for the above definition) and (also called C_{m}) represents the denaturant concentration at which 50% of the molecules are folded (the denaturation midpoint of the transition, where p_{N} = p_{U} = 1 / 2). In practice, the observed experimental data at different denaturant concentrations are fit to a twostate model with this functional form for ΔG, together with linear baselines for the folded and unfolded states. The m and are two fitting parameters, along with four others for the linear baselines (slope and intercept for each line); in some cases, the slopes are assumed to be zero, giving four fitting parameters in total. The conformational stability ΔG can be calculated for any denaturant concentration (including the stability at zero denaturant) from the fitted parameters m and . When combined with kinetic data on folding, the mvalue can be used to roughly estimate the amount of buried hydrophobic surface in the folding transition state. Structural probesUnfortunately, the probabilities p_{N} and p_{U} cannot be measured directly. Instead, we assay the relative population of folded molecules using various structural probes, e.g., absorbance at 287 nm (which reports on the solvent exposure of tryptophan and tyrosine), farultraviolet circular dichroism (180250 nm, which reports on the secondary structure of the protein backbone) and nearultraviolet fluorescence (which reports on changes in the environment of tryptophan and tyrosine). However, nearly any probe of folded structure will work; since the measurement is taken at equilibrium, there is no need for high time resolution. Thus, measurements can be made of NMR chemical shifts, intrinsic viscosity, solvent exposure (chemical reactivity) of side chains such as cysteine, backbone exposure to proteases, and various hydrodynamic measurements. To convert these observations into the probabilities p_{N} and p_{U}, one generally assumes that the observable A adopts one of two values, A_{N} or A_{U}, corresponding to the native or unfolded state, respectively. Hence, the observed value equals the linear sum
By fitting the observations of A under various solution conditions to this functional form, one can estimate A_{N} and A_{U}, as well as the parameters of ΔG. The fitting variables A_{N} and A_{U} are sometimes allowed to vary linearly with the solution conditions, e.g., temperature or denaturant concentration, when the asymptotes of A are observed to vary linearly under strongly folding or strongly unfolding conditions. Thermal denaturationAssuming a two state denaturation as stated above, one can derive the fundamental thermodynamic parameters namely, ΔH, ΔS and ΔG provided one has knowledge on the ΔC_{p} of the system under investigation. The thermodynamic observables of denaturation can be described by the following equations:
→
→
→ → where , and indicate the enthalpy,entropy and Gibbs free energy of unfolding under a constant pH and pressure. The temperature, is varied to probe the thermal stability of the system and is the temperature at which half of the molecules in the system are unfolded. The last equation is known as the GibbsHelmholtz equation. Determining the heat capacity of proteinsIn principle one can calculate all the above thermodynamic observables above from a single differential scanning calorimetry thermogram of the system assuming that the is independent of the temperature. However, it is difficult to obtain accurate values for this way. More accurately, the can be derived from a the variations in vs. which can be achieved from measurements with slight variations in or protein concentration. The slope of the linear fit is equal to the . Note that any nonlinearity of the datapoints indicates that is probably not independent of the temperature. Alternatively, the can be estimated very accurately from the calculation of the accessible solvent area (ASA) of a protein prior and after thermal denaturation as follows:
For proteins that have a known 3d structure, the can be calculated through computer programs such as Deepview (also known as swiss PDB viewer. The ASA_{unfolded} can be calculated from tabulated values of each amino acid through the semiempirical equation:
where the subscripts polar, nonpolar and aromatic indicate the parts of the 20 naturally occurring amino acids. Finally for proteins there is a linear correlation between and through the following equation^{[3]}:
Assessing twostate unfoldingFurthermore, one can assess whether the folding proceeds according to a twostate unfolding as described above. The can be done with differential scanning calorimetry by comparing the calorimetric enthalpy of denaturation i.e. the area under the peak, to the van 't Hoff enthalpy described as follows:
at the can be described as:
When a twostate unfolding is observed the . Other forms of denaturationAnalogous functional forms are possible for denaturation by pressure, pH, etc. References


This article is licensed under the GNU Free Documentation License. It uses material from the Wikipedia article "Equilibrium_unfolding". A list of authors is available in Wikipedia. 