Protein Folding Dynamics

Many small single-domain proteins fold in a two-state manner, with only the native (N) and unfolded (U) states significantly populated:

Equilibrium Thermodynamics

The equilibrium constant and free energy of folding are:

where $\Delta G = G_U - G_N > 0$ for a stable protein. At the midpoint concentration of denaturant$C_m$, $\Delta G = 0$ and $K = 1$ ($k_f = k_u$).

The fraction unfolded at equilibrium is:

Kinetics: The Chevron Plot

For a two-state protein, relaxation to equilibrium follows single-exponential kinetics with observed rate:

Both rates depend on denaturant concentration [D] according to linear free energy relationships:

where $m_f < 0$ (folding slows with denaturant) and $m_u > 0$ (unfolding accelerates). A plot of $\ln k_{\text{obs}}$ vs [D] produces the characteristic V-shaped chevron plot: the left arm reflects folding and the right arm reflects unfolding.

Tanford $\beta$-Value

The Tanford $\beta$ value quantifies the position of the transition state on the folding reaction coordinate:

where $m_{\text{eq}} = m_f - m_u$ is the equilibrium m-value. A $\beta_T \approx 0.7$ indicates that the transition state has ~70% of the solvent-accessible surface area buried relative to the native state — it is compact and native-like.

Linear Extrapolation Method (LEM)

The free energy of folding in water is obtained by linear extrapolation:

At $C_m$, $\Delta G = 0$, so $\Delta G^{\text{H}_2\text{O}} = m_{\text{eq}} \cdot C_m$. Kinetic and equilibrium measurements must agree: $\Delta G^{\text{H}_2\text{O}} = -RT\ln(k_f^{\text{H}_2\text{O}}/k_u^{\text{H}_2\text{O}})$.

$\Phi$-value analysis, developed by Alan Fersht, is the most powerful experimental method for mapping the structure of the transition state (TS) at residue-level resolution.

Definition

For a mutation that destabilizes the native state by $\Delta\Delta G°$ and raises the folding barrier by $\Delta\Delta G^\ddagger$:

where the individual terms are computed from kinetic data:

Interpretation

Requirements for Reliable $\Phi$-Values

Valid $\Phi$-value analysis requires:

Key Findings from $\Phi$-Value Studies

Decades of $\Phi$-value analysis across many proteins have revealed that:

The Zimm-Bragg model (1959) provides an exact statistical mechanical treatment of the helix-coil transition in polypeptides using a transfer matrix formalism.

The Two Parameters

Transfer Matrix Formulation

Each residue is in state c (coil) or h (helix). The statistical weight of each pair transition defines the transfer matrix:

The rows index the state of residue $i$ (top: c, bottom: h) and the columns index the state of residue $i+1$. The partition function for a chain of $N$ residues is:

where $\mathbf{e}_1$ and $\mathbf{e}_2$ are appropriate boundary vectors.

Eigenvalue Solution

The eigenvalues of $\mathbf{M}$ are:

For large $N$, $Z \approx \lambda_1^N$ (the larger eigenvalue dominates), and the helix fraction is:

Key features of the Zimm-Bragg model:

Temperature Dependence

The propagation parameter $s$ depends on temperature through:

where $\Delta H_{\text{res}} \approx -4\;\text{kJ/mol}$ is the enthalpy change per residue for helix formation and $T_m$ is the melting temperature. At $T = T_m$, $s = 1$ and the helix fraction is 0.5.

The GroEL/GroES System

GroEL is a tetradecameric (14-subunit) chaperonin arranged as two stacked heptameric rings, forming a barrel with a central cavity (~45 Å diameter). GroES is a heptameric co-chaperonin lid.

The energetic cost is 7 ATP per folding cycle. Approximately 10–15% of E. coli proteins are GroEL substrates, primarily those with complex $\alpha/\beta$ topologies (TIM barrels, Rossmann folds).

The Hsp70 System

Hsp70 (DnaK in bacteria) is the most abundant cellular chaperone. It works with the co-chaperone Hsp40 (DnaJ) and nucleotide exchange factor (GrpE/BAG). The mechanism:

Hsp70 recognizes short hydrophobic segments (~5 residues) that are exposed in unfolded or misfolded proteins. By repeatedly binding and releasing, Hsp70 prevents aggregation and gives the substrate multiple opportunities to fold correctly — another form of iterative annealing.

Kinetic Partitioning Model

The competition between productive folding and aggregation can be described by kinetic partitioning:

Chaperones act by reducing $k_{\text{agg}}$ (sequestering aggregation-prone intermediates) and effectively increasing the folding yield. This is particularly important during heat shock, where the concentration of unfolded proteins rises dramatically.

Amyloid Diseases

Amyloidoses are characterized by the deposition of cross-$\beta$ fibrils — highly ordered, insoluble protein aggregates with a shared structural core: a stack of $\beta$-strands running perpendicular to the fibril axis with inter-strand H-bonds running parallel to the axis (the cross-$\beta$ motif).

The thermodynamic driving force for amyloid formation is that the cross-$\beta$ structure is often the global free energy minimum for polypeptide chains. Native protein structures are kinetically trapped metastable states separated from the amyloid state by large barriers.

Prion Diseases

Prions ($\text{PrP}^{\text{Sc}}$) are infectious misfolded forms of the prion protein ($\text{PrP}^{\text{C}}$). The protein-only hypothesis (Prusiner, Nobel 1997) states that $\text{PrP}^{\text{Sc}}$ propagates by templating the conversion of normal$\text{PrP}^{\text{C}}$ to the misfolded form. The conversion involves a dramatic structural rearrangement:$\text{PrP}^{\text{C}}$ is predominantly $\alpha$-helical, while $\text{PrP}^{\text{Sc}}$ is rich in $\beta$-sheet. The mechanism follows a nucleated polymerization model:

where the first term describes elongation (fast) and the second describes de novo nucleation (extremely slow, accounting for long incubation periods).

Drug Design Targeting Misfolded Proteins

Relative Contact Order (RCO)

Plaxco, Simons, and Baker (1998) discovered that the folding rates of two-state proteins correlate remarkably well with a simple topological property: the relative contact order.

where $L$ is the total number of residues, $N_c$ is the number of native contacts, and$|i - j|$ is the sequence separation between contacting residues $i$ and $j$. The correlation with folding rate is:

with correlation coefficient $r \approx 0.8$. Proteins with predominantly local contacts ($\alpha$-helical, low RCO) fold faster than those with many non-local contacts ($\beta$-sheet rich, high RCO). This supports the idea that topology is the primary determinant of folding rate.

Chain Length Dependence

For two-state folders, the folding rate also depends on chain length:

The exponent $\sim 0.6$ is consistent with polymer theory: the conformational search time scales with the number of effective chain segments raised to a power related to the Flory exponent.

Nucleation-Condensation Mechanism

The dominant mechanism for two-state folding is nucleation-condensation(Fersht, 1995). In contrast to the earlier framework model (secondary structure forms first) and hydrophobic collapse model (collapse occurs first), nucleation-condensation proposes that:

Diffusion-Collision Model

For larger proteins, the diffusion-collision model (Karplus and Weaver) describes folding as a hierarchical process: pre-formed microdomains (secondary structure elements) diffuse and collide to form the native tertiary structure. The rate depends on the diffusion rate of the microdomains and the probability that a collision is productive:

where $P_{\text{productive}}$ is the probability that colliding microdomains are correctly oriented and$P_{\text{correct}}$ accounts for the combinatorics of forming all necessary contacts.

Stopped-Flow Kinetics

The workhorse for measuring folding/unfolding kinetics on the millisecond timescale. Two solutions (protein + denaturant at different concentrations) are rapidly mixed (dead time ~1 ms), and the signal (fluorescence, CD, absorbance) is monitored as a function of time. The observed rate constant is extracted by fitting to single or multi-exponential functions:

Temperature Jump (T-Jump)

Ultrafast heating (nanoseconds) using infrared laser pulses or electrical discharge perturbs the folding equilibrium, enabling the study of folding dynamics on the microsecond timescale. The temperature change is typically 5–15°C. Combined with fluorescence or IR spectroscopy, T-jump can reveal the earliest events in folding: helix formation, hydrophobic collapse, and the formation of the folding nucleus.

Hydrogen/Deuterium Exchange (HDX)

Backbone amide hydrogens exchange with solvent D$_2$O at rates that depend on their structural environment. In a native protein, amides involved in H-bonds or buried in the core exchange slowly (protection factors of $10^3$ to $10^8$). The exchange rate is:

where $k_{\text{int}}$ is the intrinsic (unprotected) exchange rate, $k_{\text{op}}$ and$k_{\text{cl}}$ are the local opening and closing rates, and the protection factor $P_f = k_{\text{cl}}/k_{\text{op}}$. Under EX2 conditions (most physiological):

HDX monitored by NMR provides residue-level information; HDX-MS provides peptide-level resolution for larger proteins and complexes.

Single-Molecule FRET

Fluorescence resonance energy transfer between donor and acceptor dyes attached to specific sites on the protein reports on intramolecular distances in real time. The FRET efficiency is:

where $r$ is the donor-acceptor distance and $R_0$ is the Forster radius (typically 40–60 Å). Single-molecule experiments reveal conformational heterogeneity and rare folding intermediates that are hidden in ensemble-averaged experiments.