Depletion interactions in polymer-colloid mixtures in the presence of the attraction between polymers and particle surfaces are computed from the simulation data on polymer partitioning. The gap between two large colloid particles is approximated by the slit plates immersed in a dilute polymer solution. The difference between the intra-slit and bulk pressures of polymers Delta p (net pressure) and the effective pair potential of plates W is calculated as a function of the plate separation D and the attraction strength epsilon of plates to polymer segments. Weak attraction strengths epsilon are considered, including the region around the compensation point epsilon(c) (the depletion/adsorption threshold) where the partition coefficient of polymers K = 1. A fit of the simulation data for the free energy confinement Delta A = -k(B)T In K served as a base of computations of the above depletion characteristics. The analytical functions deduced for Delta p(D,epsilon) and W(D,epsilon) predict a reduction of the range and depth of depletion interaction by an increase in the attraction strength vertical bar epsilon vertical bar. Variation of the thickness delta of depletion layers near the plates with attraction e was also determined. Consequences of reduced depletion interaction for the effective second virial coefficient B-2 and for the stability of colloidal suspensions were discussed. At the compensation point epsilon(c) the depletion effect disappears and the stability of polymer-colloid mixtures in dilute solution is not affected by the polymer concentration. (c) 2005 Elsevier Ltd. All rights reserved.