The dissimilarity between the Reynolds stresses and the heat fluxes in perturbed turbulent channel and plane Couette flows was studied using direct numerical simulation. The results demonstrate that the majority of the dissimilarity was due to the difference between the wall-normal fluxes, while the difference between the streamwise fluxes was lower. The main causes for the dissimilarity were the production terms, followed by the velocity-pressure interaction terms. Further insights into the importance of the velocity-pressure interaction in the origin of the dissimilarity are provided using two-point correlation. Furthermore, an octant conditional averaged dataset reveals that not only the wall-normal heat flux but also the streamwise heat flux is strongly related to the wall-normal gradient of the mean temperature. A simple Reynolds-averaged Navier–Stokes (RANS) heat flux model is proposed as a function of the Reynolds stresses. A comparison of the direct numerical simulation data with an “a priori” prediction suggests that this simple model performs reasonably well.