Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Missing compartment topology connections in VCell bngl output #246

Open
jrfaeder opened this issue Apr 14, 2022 · 1 comment
Open

Missing compartment topology connections in VCell bngl output #246

jrfaeder opened this issue Apr 14, 2022 · 1 comment

Comments

@jrfaeder
Copy link
Contributor

This issue needs further investigation using VCell itself and with VCell developers, but putting a placeholder here. VCell seems to export compartmental models that don't have parent compartments listed in the compartments block. Not sure if this is all compartmental models exported from VCell or just a subset.

Here's a very simple example with two compartments EC and mem:

begin compartments
    EC	3	1
    mem	2	1 
end compartments

Because mem does not have a parent compartment, compartmental rules involving species in adjacent compartments do not generate reactions:

r0:	@EC:Lig(Site0) + @mem:Receptor(Site0~U) -> @mem:Lig(Site0!1).Receptor(Site0~U!1)		1.0E-4

A complete example model is
lig_receptor.bngl.txt

@jrfaeder
Copy link
Contributor Author

The simple fix for this model is to add a parent compartment to the compartments block:

begin compartments
EC	3	1
mem	2	1 EC
end compartments

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant