Starting from physical insight on the energy transfer phenomena in wall turbulent flows, it is shown how modeling of subgrid stresses in large-eddy simulations can be improved. Each model should aim at reproducing the double feature of energy sink and source of the small scales of wall flows which become relevant when large filter lengths are considered. Here we propose one possible choice where the main ingredient is the coupling of the classical linear formulation of eddy viscosity with the nonlinear anisotropic features of the velocity increments tensor. This approach, which actually presents most of the features of the mixed models, captures the near-wall dynamics for very large filter lengths reproducing the small scales source physics responsible for backward energy transfer. A posteriori tests show excellent agreement with direct numerical simulation of turbulent channel flows even when very coarse grids are considered. The capability of the balance of the filtered second order structure function as a post-processing tool to evaluate the physics of any model is also shown.