Lactoperoxidase (LP), which is a high-value minor whey protein, has recently drawn extensive attention from research scientists and industry due to its multi-function and potential therapeutic applications. In this study, the separation and optimization of two similar-sized proteins, LP and lactoferrin (LF) were investigated using strong cation exchange column chromatography. Optimization was started with central composite design based experiments to characterize the importance of different decision variables. The three variables used in the optimization were flow rate, length of gradient and final salt concentration in the linear elution gradient step. The obtained empiric functional model represented the effect of the significant factors on the yield as the objective function. Afterwards, the calibrated mechanistic model was employed to predict accurate optimal set of variables. The optimal operating points were found and the results were compared with validation experiments. Predictions respecting yield confirmed a very good agreement with experimental results while keeping purity, a product quality characteristic, equal or above to a predefined value.