The effect of solute concentration on the equilibrium partitioning of sphere-like, colloidal solutes in stiff polymer hydrogels is examined theoretically and experimentally. The theoretical development is a statistical mechanics approach, and allows quantitative calculations to be performed to determine the concentration-dependent partition coefficient correct to first order in solute concentration at specific surface charge densities. The theory predicts that repulsive steric and/or electrostatic solute-fiber interactions exclude solute from the gel phase, but that repulsive solute-solute interactions cause partitioning into the gel to increase with increasing solute concentration. These trends are enhanced for larger solutes, increased fiber volume fractions, or stronger electrostatic repulsion. Partition coefficients have also been measured for two proteins, bovine serum albumin (BSA) and alpha -lactalbumin (ALA), in a system consisting of a salt solution and cubes of agarose hydrogel. To investigate the effect of electrostatic interactions, the experiments were performed at 0.15 M KCI and 0.01 M KCI. The theory underpredicts the strong electrostatic repulsion between BSA macromolecules at the lower ionic strength. The experimental results for ALA show the influence of an attractive interaction between the protein macromolecules, in addition to hard-sphere repulsive and electrostatic interactions.