This paper addresses the problem of heat transfer in the near-wall region of a turbulent channel flow. We investigate the role of coherent eddy structures on the thermal exchange between the solid and fluid materials. In the present approach, the instantaneous velocity field is modelled by means of a low-order dynamical system based on proper orthogonal decomposition (POD) modes while the temperature field is accounted for by using a stochastic filtered density function (FDF) method. A good comparison of results with reference DNS data on mean temperature and on temperature r.m.s. at various Prandtl numbers and for different wall boundary conditions (isotemperature, isoflux, and conjugate heat transfer) is achieved.