How to set boundary conditions for a reaction to occur in SimBiology?

Question

Tessa Warmink 2021-9-30

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1463844-how-to-set-boundary-conditions-for-a-reaction-to-occur-in-simbiology

评论： Arthur Goldsipe 2021-10-3

Hi there!

I am building a model in SimBiology where I want to model the production gas vesicle proteins and subsequent formation of gas vesicles. Ideally, I want to model this such that there is a boundary condition where an X amount of gas vesicle proteins (or above the essential concentration) will lead to an assembly reaction.

I tried this by adding the chemical reaction: X [Gas vesicle protein] -> [Gas vesicle]

However, this leads to a negative values for the gas vesicle species. Is there a way in which I can set a condition for this reaction to occur (for example require that enough gas vesicle proteins are present in order for the reaction to be active)? Or could there maybe be a better way to solve this?

The model I'm building is a stiff model, and I am using the sundials solver.

Thank you in advance!

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Arthur Goldsipe 2021-10-1

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1463844-how-to-set-boundary-conditions-for-a-reaction-to-occur-in-simbiology#answer_799134

Hi,

Take a look at this answer I shared. It's from several years ago, but the information is still accurate. Some additional tips are also listed here. As mentioned there, you might want to look at absolute tolerance scaling. If your model changes too quickly at the beginning of the simulation, you may need to turn off this feature or reduce the absolute tolerance step size.

If you can share your model, we might be able to offer more specific suggestions. For example, we can help decide whether the negative concentrations represent a modeling error, an acceptable value you can ignore, or a numerical instability that requires modifications to your model.

-Arthur (he/him)

2 个评论
显示无隐藏无

Tessa Warmink 2021-10-3

Thank you so much for your helpful answer! I tried using the tips from the links you added, but the negative value that occurs in my model is larger than the absolute tolerance that I have set, so I think that means the problem is in my model.

I cannot upload my model here as the file is too large, but could the problem be with the assembly reaction (where multiple proteins congregate into one large complex)? The reaction looks like this:

23 GvpA --> Gas vesicle

with kf = 1E-8 (1/(second*molarity^22))

Arthur Goldsipe 2021-10-3

Possibly. What's the reaction rate? We should see what happens to the rate if a concentration goes slightly negative. Sometimes this causes instabilities. Such problems can sometimes be fixed by replacing species X with max(0,X) in the rate. What happens if you do that in all rates involving any X that goes negative?

请先登录，再进行评论。