In this paper we examine the mutual interactions between microbubbles and turbulence in vertical channel flow. An Eulerian-Lagrangian approach based on pseudo-spectral direct numerical simulation is used: bubbles are momentum coupled with the fluid and are treated as pointwise spheres subject to gravity, drag, added mass, pressure gradient, Basset and lift forces. Two different flow configurations (upward and downward channel flow of water at shear Reynolds number Re-tau = 150) and four different bubble diameters are considered, assuming that bubbles are non-deformable (i.e. small Eotvos number) and contaminated by surfactants (i.e. no-slip condition applies at bubble surface). Confirming previous knowledge, we find macroscopically differen...