A BGK relaxation model for polyatomic gas mixtures