You could use something like this
rules := { \gamma^{\star} -> 1/ 24 \epsilon_{m n p q} \Gamma^{m n p q} };
and then compute, say $(\gamma^5)^2$ as follows
gs_gs := \gamma^{\star} \gamma^{\star};
substitute(gs_gs, rules, repeat=True);
epsilon_to_delta(gs_gs);
expand_delta(gs_gs)
distribute(_)
eliminate_kronecker(_)
canonicalise(_)
collect_terms(_)
join_gamma(_);
Cheers.