Exemplo n.º 1
0
from pysb import Model, Observable

from mEncoder.mechanistic_model import (add_monomer_synth_deg,
                                        generate_pathway, add_activation,
                                        add_observables, add_inhibitor,
                                        add_gf_bolus, cleanup_unused)

model = Model('EGFR_MAPK')

rtkfs = ['EGF']

# EGFR
for rtkf_name in rtkfs:
    add_monomer_synth_deg(rtkf_name)
    add_gf_bolus(model, rtkf_name, [rtkf_name])

erbb_cascade = [
    ('EGFR', {
        'Y1173': ['EGF']
    }),
    ('ERBB2', {
        'Y1248': ['EGFR__Y1173_p']
    }),
]
generate_pathway(model, erbb_cascade)

active_rtks = ['EGFR__Y1173_p', 'ERBB2__Y1248_p']

mapk_cascade = [
    ('RAF1', {
        'S338': active_rtks + ['EGF']
Exemplo n.º 2
0
from pysb import Model

from mEncoder.mechanistic_model import (
    add_monomer_synth_deg, generate_pathway, add_activation,
    add_abundance_observables, add_phospho_observables
)

model = Model('FLT3_MAPK')

# FLT3
for rtkf_name in ['FL']:
    add_monomer_synth_deg(rtkf_name)


rtk_cascade = [
    ('FLT3',  {'Y843': ['FL']}),
]
generate_pathway(model, rtk_cascade)

# ERK
for ras_name in ['HRAS', 'KRAS', 'NRAS']:
    add_monomer_synth_deg(ras_name, nsites=['N'])
    add_activation(
        model, ras_name, 'N', 'nucleotide_exchange',
        ['FLT3__Y843_p']
    )

mapk_cascade = [
    ('RAF1',   {'S338':      ['KRAS__N_gtp', 'HRAS__N_gtp', 'NRAS__N_gtp']}),
    ('BRAF',   {'S447':      ['KRAS__N_gtp', 'HRAS__N_gtp', 'NRAS__N_gtp']}),
    ('MAP2K1', {'S218_S222': ['RAF1__S338_p', 'BRAF__S447_p']}),
from pysb import Model, Observable

from mEncoder.mechanistic_model import (add_monomer_synth_deg,
                                        generate_pathway, add_activation,
                                        add_observables, add_inhibitor,
                                        add_gf_bolus)

model = Model('EGFR_MAPK_AKT_STAT')

rtkfs = ['EGF']

# EGFR
for rtkf_name in rtkfs:
    add_monomer_synth_deg(rtkf_name)
    add_gf_bolus(model, rtkf_name, [rtkf_name])

erbb_cascade = [
    ('EGFR', {
        'Y1173': ['EGF']
    }),
    ('ERBB2', {
        'Y1248': ['EGF']
    }),
]
generate_pathway(model, erbb_cascade)

active_rtks = ['EGFR__Y1173_p', 'ERBB2__Y1248_p']
stat_rtks = ['EGFR__Y1173_p', 'ERBB2__Y1248_p']

# MAPK
for ras_name in ['HRAS', 'KRAS', 'NRAS']: