A procedure to determine effective diffusion coefficients of proteins in chromatographic gels, required as model input parameters, is presented. An experimental methodology based on frontal liquid chromatography was combined with a numerical methodology based on a detailed mathematical model describing the chromatographic process including the extra-column dispersion, the dispersion due to the packed bed, the external mass transfer from the bulk phase to the stationary phase, and the diffusive transport within the stationary phase. The procedure has several advantages compared to previously reported methods to determine diffusion coefficients in that no other equipment than an HPLC is required, any class of stationary phases can be investigated as long as the experiments are performed under non-binding conditions, and no modification, e.g., molding of slabs or membranes, to the stationary phase is required. To show the applicability of the procedure, the effective diffusion coefficients of lysozyme, bovine serum albumin, and immunoglobulin γ in Sepharose™ CL-4B were determined and shown to be comparable to those determined by other methods. Molecularly imprinted polymers (MIPs) are man-made polymeric materials with molecular recognition abilities. They mimic the molecular recognition of naturally occurring molecular recognition elements such as receptors and antibodies by binding target molecules by either non-covalent, covalent, or metal-coordinating interactions. Traditionally, MIPs are synthesized in the form of monolithic polymers which are subsequently crushed, ground, and sieved to an appropriate size range. In this thesis, a suspension polymerization method to prepare MIPs in the shape of spherical beads is presented. The method involves suspending a prepolymerization solution in mineral oil, used as the continuous phase. The droplets are transformed into solid spherical beads by free-radical polymerization. The beads have been shown to compare well to the traditional irregularly shaped particles prepared from monoliths. The advantages of the method compared to previously reported methods are the low cost and commercial availability of the continuous phase and the absence of the need for stabilizers for the formation of droplets of pre-polymerization solution in the mineral oil. When compared to the method to prepare particles from monolithic polymers, this method is advantageous due to the spherical shape of the resulting beads and the reduction in time needed to prepare a MIP. When a new MIP is designed, the traditional approach is to use eithera previously reported protocol or rules of thumb based on previous knowledge. This results in non optimized MIPs. The number of possible combinations of monomers, cross-linkers, solvents, and initiators are huge. A full optimization of a MIP formulation therefore requires a large number of experiments. To facilitate the efforts, chemometrics was applied to the work described in this thesis. Three factors (i.e., the amount of monomer, the amount of crosslinker, and the amount of porogen) were chosen as the factors in the model. Multivariate data analysis of the binding to the MIPs indicated how the factors influenced the binding and an optimized MIP composition was identified. The combined use of the suspension polymerization method to produce spherical beads with the application of chemometrics was shown in this thesis to drastically reduce the number of experiments and the time needed to design and optimize a new MIP.