Non-divergence equations structured on Hörmander vector fields: heat kernels and Harnack inequalities