Single-Component Adsorption Isotherm Models

Key Insight

Adsorption isotherms express the equilibrium concentration of adsorbed protein in the stationary phase relative to its concentration in the mobile phase. This relationship is generally linear at low protein concentrations but becomes non-linear, eventually leveling off to a maximum capacity at higher concentrations. The linear limit depends on accessible binding sites and specific protein affinity, while maximum capacity is typically limited by accessible surface area or binding site concentration. Mobile phase composition, including salt concentration and pH, significantly influences protein adsorption, necessitating constant conditions for isotherm measurements. For example, ion-exchange chromatography isotherms are highly favorable at low salt concentrations but less so as salt increases, while hydrophobic interaction chromatography is promoted by increasing kosmotropic salt concentrations, such as ammonium sulfate.

The Langmuir isotherm model, originally developed for gas adsorption, is commonly employed for proteins, assuming a stoichiometric association of an adsorbate molecule with a surface ligand. It is described by the equation q = q_m K C / (1 + K C), where q is adsorbed concentration, q_m is maximum adsorption capacity, C is solution concentration, and K is the equilibrium constant. At low C, the isotherm approaches a linear limit (q ≈ q_m K C), and at high C, it approaches q_m. A linearized form (Scatchard plot) can determine parameters, though strong initial slopes in protein chromatography can lead to uncertain K-value determination. The separation factor R, defined as 1 / (1 + K C_ref), is a process-dependent parameter that describes the isotherm's curvature, where R approaching 1 indicates a linear isotherm (for dilute feeds) and R approaching 0 indicates a nearly rectangular isotherm (for concentrated feeds).

Beyond Langmuir, various empirical models describe complex protein adsorption behaviors. The Freundlich isotherm (q = a C^(1/b)) implies a heterogeneous surface but lacks a linear limit or maximum capacity. The Temkin isotherm (q = a ln(1 + K C)) is proposed for strongly heterogeneous surfaces. The Toth isotherm (q = q_m C / (1 + (C/K_b)^b)^(1/b)) addresses heterogeneity, linearity at dilute concentrations, and a maximum capacity. The Langmuir-Freundlich or Sips isotherm (q = q_m (K C)^b / (1 + (K C)^b)) can describe sigmoidal curves, suggesting cooperative adsorption when 'b' is greater than 1. The Brunauer, Emmett, and Teller (BET) isotherm models multilayer adsorption, relevant for proteins that aggregate on surfaces.

For ion-exchange systems, the Stoichiometric Displacement (SD) model posits protein adsorption occurs via stoichiometric exchange with counter-ions (e.g., P^z+ + z R^-Na+ <=> R^-P^z+ + z Na+), where 'z' is the effective binding charge. The Steric Mass Action (SMA) model refines this by incorporating a steric hindrance factor 'sigma', accounting for additional ligands shielded by the protein's footprint. Both models predict similar linear behavior at low loadings, but their maximum capacities differ: q_max = q0 / z for SD and q_max = q0 / (z + sigma) for SMA, where q0 is the total concentration of protein-accessible ligands. The effective charge 'z', typically ranging from 3 to 12, is generally smaller than the protein's net charge and can be determined from the slope of a log-log plot relating isotherm slope to counter-ion concentration. For instance, in cytochrome c adsorption, the SMA model can fit data across broad ranges with fewer parameters than Langmuir, which requires different parameters per salt concentration.