In [1]:

import actonet
import biolqm

This notebook has been executed using the docker image colomoto/colomoto-docker:2020-07-01

In [2]:

model = biolqm.load("https://raw.githubusercontent.com/jgtz/StableMotifs/master/TLGLNetwork.txt", "booleannet")

Downloading https://raw.githubusercontent.com/jgtz/StableMotifs/master/TLGLNetwork.txt

In [3]:

bn = biolqm.to_minibn(model, ensure_boolean=True)
inputs = bn.constants()
inputs

Out[3]:

{'CD45': False,
 'IL15': True,
 'PDGF': False,
 'Stimuli': True,
 'Stimuli2': False,
 'TAX': False}

In [4]:

for n in inputs:
    bn[n] = n

In [5]:

act = actonet.ActoNet(bn, inputs)

In [6]:

r = act.reprogramming_fixpoints(Apoptosis=1, maxsize=2)
r

Out[6]:

[FromCondition('input', PermanentPerturbation(Apoptosis=1, TCR=0)),
 FromCondition('input', PermanentPerturbation(CTLA4=0, Caspase=1)),
 FromCondition('input', PermanentPerturbation(Caspase=1, TCR=0)),
 FromCondition('input', PermanentPerturbation(CTLA4=0, DISC=1)),
 FromCondition('input', PermanentPerturbation(Apoptosis=1, CTLA4=0)),
 FromCondition('input', PermanentPerturbation(DISC=1, TCR=0)),
 FromCondition('input', PermanentPerturbation(CTLA4=1, Caspase=1)),
 FromCondition('input', PermanentPerturbation(CTLA4=1, DISC=1)),
 FromCondition('input', PermanentPerturbation(Apoptosis=1, CTLA4=1)),
 FromCondition('input', PermanentPerturbation(Caspase=1, TCR=1)),
 FromCondition('input', PermanentPerturbation(DISC=1, TCR=1)),
 FromCondition('input', PermanentPerturbation(Apoptosis=1, TCR=1))]

In [7]:

r.as_table()

Out[7]:

	Apoptosis	CTLA4	Caspase	DISC	TCR
0				1	0
1				1	1
2			1		0
3			1		1
4		0		1
5		0	1
6		1		1
7		1	1
8	1				0
9	1				1
10	1	0
11	1	1

In [8]:

r = act.reprogramming_fixpoints(Apoptosis=0, maxsize=2)
r

Out[8]:

[FromCondition('input', PermanentPerturbation(Apoptosis=0, CTLA4=0)),
 FromCondition('input', PermanentPerturbation(Apoptosis=0, TCR=0)),
 FromCondition('input', PermanentPerturbation(CTLA4=0, Caspase=0)),
 FromCondition('input', PermanentPerturbation(CTLA4=0, TBET=0)),
 FromCondition('input', PermanentPerturbation(Caspase=0, TCR=0)),
 FromCondition('input', PermanentPerturbation(TBET=0, TCR=0)),
 FromCondition('input', PermanentPerturbation(S1P=1, TCR=0)),
 FromCondition('input', PermanentPerturbation(CTLA4=0, S1P=1)),
 FromCondition('input', PermanentPerturbation(Apoptosis=0, CTLA4=1)),
 FromCondition('input', PermanentPerturbation(CTLA4=1, Caspase=0)),
 FromCondition('input', PermanentPerturbation(CTLA4=1, TBET=0)),
 FromCondition('input', PermanentPerturbation(CTLA4=1, S1P=1)),
 FromCondition('input', PermanentPerturbation(Apoptosis=0, TCR=1)),
 FromCondition('input', PermanentPerturbation(Caspase=0, TCR=1)),
 FromCondition('input', PermanentPerturbation(TBET=0, TCR=1)),
 FromCondition('input', PermanentPerturbation(S1P=1, TCR=1))]

In [9]:

r.as_table()

Out[9]:

	Apoptosis	CTLA4	Caspase	S1P	TBET	TCR
0					0	0
1					0	1
2				1		0
3				1		1
4			0			0
5			0			1
6		0			0
7		0		1
8		0	0
9		1			0
10		1		1
11		1	0
12	0					0
13	0					1
14	0	0
15	0	1

In [ ]: