I want to model the relationship between a gene and its product's molecular function.

graph LR
    gene["AKT1 (HGNC:391)"] -- "has gene product (RO:0002205)" --> protein["RAC-alpha serine/threonine-protein kinase (uniprot:P31749)"] -- "enables (RO:0002327)" --> activity["non-specific serine/threonine protein kinase (EC:2.7.11.1)"]
    gene -- has gene product that enables --> activity

I see this as being a simple property chain of has gene product (RO:0002205) and enables (RO:0002327)

Further discussion about EC/Rhea/GO being molecular functions in #783

Motivation

This is directly motivated by making a more accurate ontologization of HGNC in biopragmatics/pyobo#455