The Zonal Rate Model (ZRM) has previously been shown to accurately account for contributions to elution band broadening, including external flow nonidealities and radial concentration gradients, in ion-exchange membrane (IEXM) chromatography systems operated under nonbinding conditions. Here, we extend the ZRM to analyze and model the behavior of retained proteins by introducing terms for intra-column mass transfer resistances and intrinsic binding kinetics. Breakthrough curve (BTC) data from a scaled-down anion-exchange membrane chromatography module using ovalbumin as a model protein were collected at flow rates ranging from 1.5 to 20?mL?min-1. Through its careful accounting of transport nonidealities within and external to the membrane stack, the ZRM is shown to provide a useful framework for characterizing putative protein binding mechanisms and models, for predicting BTCs and complex elution behavior, including the common observation that the dynamic binding capacity can increase with linear velocity in IEXM systems, and for simulating and scaling separations using IEXM chromatography. Global fitting of model parameters is used to evaluate the performance of the Langmuir, bi-Langmuir, steric mass action (SMA), and spreading-type protein binding models in either correlating or fundamentally describing BTC data. When combined with the ZRM, the bi-Langmuir, and SMA models match the chromatography data, but require physically unrealistic regressed model parameters to do so. In contrast, for this system a spreading-type model is shown to accurately predict column performance while also providing a realistic fundamental explanation for observed trends, including an observed increase in dynamic binding capacity with flow rate. Biotechnol. Bioeng. 2012; 109:615629. (C) 2011 Wiley Periodicals, Inc.