The rule rxn1 ~= 0 <=> not gene1 does not compile correctly. It creates two >= inequalities, not > inequalities.