A theoretical framework for the ion-exchange behaviour of bioactive substances in non-linear ion-exchange chromatogaphy is described. The aim of the study was the creation of a model basis to support a process design for production-scale ion-exchange chromatography. The theory can be applied to a whole variety of biological substances, such as amino acids, polysaccharides, peptides and proteins and either isocratic or gradient elution can be carried out. The influence of the eluent concentration on the ion-exchange as well as on the characteristic charge was considered. Experimental measurements showed a strong non-linear ion-exchange equilibrium with a transition from a Langmuir-type to a sigmoidal isotherm at higher eluent concentrations. Hereby, the compound binds to the surface though it is not ionic. Therefore, the model considered the possibility of ion-exchange as well as adsorption. A simplified distribution of the counter-ions based on the Gouy-Chapman theory with a discrete distribution of the counter-ions was used. The theory was extended by a selectivity in the double layer to allow specific adsorption. Calculations of adsorption-elution cycles showed, in agreement with the experimental observations, the development of non-linear elution profiles with a desorption fronting. As a result, the column loading and the eluent concentration were varied. The effect of contaminants, in this case sodium ions, was investigated and included in the model. Finally, the model was extended to multicomponent systems to investigate the effect of side components on the retention behaviour. The development of the characteristic elution profiles and the effect of the column loading on the separation are discussed. Calculated concentration profiles along the column at discrete time steps were used to reveal the influence of side components and the underlying separation mechanism. The simulations provided a new insight into the phenomena involved in biochromatography and make convenient design concepts at least doubtful as the separation is in this case mainly determined by the loading step and not by the choice of the elution gradient.