misc cleanups to Gene/Reaction relationship

cobrapy can use the _associate_gene and _dissociate_gene new functions
internally for reaction/gene relationships, and the add_gene and
remove_gene functions (which, being unprefixed, are more user-facing)
can be deprecated
add a unit test