def test_transport_reaction_gpr_presence(model):
"""
Expect a small fraction of transport reactions not to have a GPR rule.
As it is hard to identify the exact transport processes within a cell,
transport reactions are often added purely for modeling purposes.
Highlighting where assumptions have been made versus where
there is proof may help direct the efforts to improve transport and
transport energetics of the tested metabolic model.
However, transport reactions without GPR may also be valid:
Diffusion, or known reactions with yet undiscovered genes likely lack GPR.
Implementation:
Check which cobra.Reactions classified as transport reactions have a
non-empty "gene_reaction_rule" attribute.
"""
# TODO: Update threshold with improved insight from meta study.
ann = test_transport_reaction_gpr_presence.annotation
ann["data"] = get_ids(
basic.check_transport_reaction_gpr_presence(model)
)
ann["metric"] = len(ann["data"]) / len(
helpers.find_transport_reactions(model)
)
ann["message"] = wrapper.fill(
"""There are a total of {} transport reactions ({:.2%} of all
transport reactions) without GPR:
{}""".format(len(ann["data"]), ann["metric"], truncate(ann["data"])))
assert ann["metric"] < 0.2, ann["message"]
def test_transport_reaction_gpr_presence(model):
"""
Expect a small fraction of transport reactions not to have a GPR rule.
As it is hard to identify the exact transport processes within a cell,
transport reactions are often added purely for modeling purposes.
Highlighting where assumptions have been made versus where
there is proof may help direct the efforts to improve transport and
transport energetics of the tested metabolic model.
However, transport reactions without GPR may also be valid:
Diffusion, or known reactions with yet undiscovered genes likely lack GPR.
Implementation:
Check which cobra.Reactions classified as transport reactions have a
non-empty "gene_reaction_rule" attribute.
"""
# TODO: Update threshold with improved insight from meta study.
ann = test_transport_reaction_gpr_presence.annotation
ann["data"] = get_ids(basic.check_transport_reaction_gpr_presence(model))
ann["metric"] = len(ann["data"]) / len(
helpers.find_transport_reactions(model))
ann["message"] = wrapper.fill(
"""There are a total of {} transport reactions ({:.2%} of all
transport reactions) without GPR:
{}""".format(len(ann["data"]), ann["metric"], truncate(ann["data"])))
assert ann["metric"] < 0.2, ann["message"]
def test_check_transport_reaction_gpr_presence(model, num):
"""Expect amount of transport reactions without gpr to be identified."""
assert len(basic.check_transport_reaction_gpr_presence(model)) == num