We use molecular dynamics simulations with a dissipative particle dynamics (DPD) thermostat to Study the behavior of nanosized inclusions (Colloids) in a polymer brush which is in contact with an explicit solvent in the NPT ensemble. The brush is described by a bead-spring model for flexible polymer chains, grafted oil a solid Substrate, while the polymer-soluble nanoparticles in the solution are taken as hard spheres. By varying the chain length N, the grafting density of the brush, sigma(g), and the size of the nanoparticles b, we determine the equilibrium particle penetration depth delta and the average concentration of nanoinclusions (phi) over bar (nano) in the penetration layer delta at constant pressure. In agreement with a recent theoretical prediction, we demonstrate that for nanoinclusions of size b >= b* the thickness of this layer delta proportional to h(b*/b)(3) where h is brush height and b* proportional to sigma(-2/3)(g) is a typical size below which smaller particles are uniformly distributed in the brush. We also observe that particles, larger than some threshold value b(max) do not mix with the brush. The mean density of nanoinclusions is found to scale as (phi) over bar (nano) proportional to (b*/b)(3) within the whole range of parameter variation. The diffusivity of nanoparticles, embedded in the polymer brush, in direction perpendicular to the grafting plane is found to be up to 20% higher than parallel to the plane. The variation of the respective diffusion coefficients D-perpendicular to(nano) and D-parallel to(nano) changes with growing Volume fraction of the nanoparticles in agreement with theoretical predictions. (D 2009 Elsevier Inc. All rights reserved.