Fermionic approach to variational quantum simulation of Kitaev spin models