Non-divergence equations structured on H"{o}rmander vector fields: heat kernels and Harnack inequalities