begin parameters
  A0   1000
  B0   500
  kp1  0.46e-5
  km1  1.0
  p 100
  d 5
  t 10
end parameters

begin species
  A(b,loc~Cyt)   A0    
  B(a,Y~U,loc~Cyt) B0
end species

begin reaction rules
  A(b,loc~Cyt) + B(a,loc~Cyt) <-> A(b!1,loc~Cyt).B(a!1,loc~Cyt) kp1, km1
  A(b!1).B(a!1,Y~U) -> A(b!1).B(a!1,Y~P)  p
  B(Y~P,a,loc~Cyt) -> B(Y~P,a,loc~Nuc)  t
  B(Y~P,loc~Nuc) -> B(Y~U,loc~Nuc)  p
end reaction rules

begin observables
  Molecules A  A()
  Molecules BP  B(a!?,Y~P)
  Molecules BPC  B(a!?,Y~P,loc~Cyt)
  Molecules BPN  B(a!?,Y~P,loc~Nuc)
end observables

generate_network();
writeSBML();
simulate_ode({t_end=>50,n_steps=>20});

#%VC% mergeReversibleReactions
#%VC% speciesRenamePattern("\." , "_") 
#%VC% speciesRenamePattern("[\(,][a-zA-Z]\w*", "") 
#%VC% speciesRenamePattern("~|!\d*", "") 
#%VC% speciesRenamePattern("\(\)", "") 
#%VC% speciesRenamePattern("\)", "")  
#%VC% setUnit("all", "default")  
#%VC% compartmentalizeSpecies("loc~Nuc", "Nucleus", "Cytoplasm") 
#%VC% compartmentalizeSpecies("loc~Cyt", "Cytoplasm", "")