The rational selection of optimal protein purification sequences, as well as mathematical models that simulate and allow optimization of chromatographic protein purification processes have been developed for purification procedures such as ion-exchange, hydrophobic interaction and gel filtration chromatography. This paper investigates the extension of such analysis to affinity chromatography both in the selection of chromatographic processes and in the use of the rate model for mathematical modelling and simulation. Two affinity systems were used: Blue Sepharose and Protein A. The extension of the theory developed previously for ion-exchange and HIC chromatography to affinity separations is analyzed in this paper. For the selection of operations two algorithms are used. In the first, the value of η, which corresponds to the efficiency (resolution) of the actual chromatography and, Σ, which determines the amount of a particular contaminant eliminated after each separation step, which determines the purity, have to be determined. It was found that the value of both these parameters is not generic for affinity separations but will depend on the type of affinity system used and will have to be determined on a case by case basis. With Blue Sepharose a salt gradient was used and with Protein A, a pH gradient. Parameters were determined with individual proteins and simulations of the protein mixtures were done. This approach allows investigation of chromatographic protein purification in a holistic manner that includes ion-exchange, HIC, gel filtration and affinity separations for the first time.