Methods of theoretical and experimental evaluation of irreversible adsorption of particles, e.g., colloids and globular proteins at heterogeneous surfaces were reviewed. The theoretical models were based on the generalized random sequential adsorption (RSA) approach. Within the scope of these models, localized adsorption of particles occurring as a result of short-ranged attractive interactions with discrete adsorption sites was analyzed. Monte-Carlo type simulations performed according to this model enabled one to determine the initial flux, adsorption kinetics, jamming coverage and the structure of the particle monolayer as a function of the site coverage and the particle/site size ratio, denoted by lambda. It was revealed that the initial flux increased significantly with the site coverage theta(s) and the lambda parameter. This behavior was quantitatively interpreted in terms of the scaled particle theory. It also was demonstrated that particle adsorption kinetics and the jamming coverage increased significantly, at fixed site coverage, when the lambda parameter increased. Practically, for alpha = lambda2theta(s) > 1 the jamming coverage at the heterogeneous surfaces attained the value pertinent to continuous surfaces. The results obtained prove unequivocally that spherically shaped sites were more efficient in binding particles in comparison with disk-shaped sites. It also was predicted that for particle size ratio lambda < 4 the site multiplicity effect plays a dominant role, affecting significantly the structure of particle monolayers and the jamming coverage. Experimental results validating main aspects of these theoretical predictions also have been reviewed. These results were derived by using monodisperse latex particles adsorbing on substrates produced by covering uniform surface by adsorption sites of a desired size, coverage and surface charge. Particle deposition occurred under diffusion-controlled transport conditions and their coverage was evaluated by direct particle counting using the optical and electron microscopy. Adsorption kinetics was quantitatively interpreted in terms of numerical solutions of the governing diffusion equation with the non-linear boundary condition derived from Monte-Carlo simulations. It was proven that for site coverage as low as a few percent the initial flux at heterogeneous surfaces attained the maximum value pertinent to homogeneous surfaces. It also was demonstrated that the structure of larger particle monolayers, characterized in terms of the pair correlation function, showed much more short-range ordering than predicted for homogeneous surface monolayers at the same coverage. The last part of this review was devoted to detection of polyelectrolyte multilayers on various substrates via particle deposition experiments.