Surface Area Calculations in Simbiology

조회 수: 1 (최근 30일)
Sondra Rahmeh
Sondra Rahmeh 2018년 7월 8일
답변: Arthur Goldsipe 2018년 7월 9일
Is it possible to have non-volumetric compartments in SimBiology? I am trying to replicate a paper with ODEs based on reactions on the surface area of the cell. Whether I run my simulation with concentration or substance, I get a similar error regarding dimensional analysis "SimBiology cannot identify a compartment volume for this conversion when a reaction has reactants in multiple compartments, or when the reactant is null and the reaction has products in multiple compartments."

답변 (1개)

Arthur Goldsipe
Arthur Goldsipe 2018년 7월 9일
Hi,
In SimBiology, compartments can also represent 2-dimensional areas, such as the surface of a cell. I believe the error you are seeing is not really about any 2-dimensional compartments, but it part of a larger error message indicating that you need to write your reaction rate with units of amount per time. The full error message should look something like this:
"The reaction rate for reaction ... must be defined with dimensions of amount per time not concentration per time. SimBiology converts each reaction rate to dimensions of amount per time, even when dimensional analysis is disabled. SimBiology cannot identify a compartment volume for this conversion when a reaction has reactants in multiple compartments, or when the reactant is null and the reaction has products in multiple compartments."
Usually, this means you need to multiply your rate constants by the appropriate volumes or areas. Here's a simple example based on the model in this paper. Consider a reaction for a surface receptor binding with a ligand to form a surface complex: Surface.Receptor + Volume.Ligand <-> Surface.Complex. Regardless of the units used to measure Receptor, Ligand, and Complex, you will need to write the reaction in units of amount/time, in this case molecule/minute. In the paper, Receptor and Complex are measured in molecule, and Ligand is measure in nanomole. The forward reaction in units of molecule/time is k_on*Receptor*Ligand, where k_on has units of 1/nanomolar/minute. The reverse reaction in units of molecule/time is k_off*Complex, where k_off has units of 1/minute. So in SimBiology, you would write this reaction rate as k_on*Surface.Receptor*Volume.Ligand - k_off*Surface.Complex. SimBiology then knows how to convert this to the appropriate rate of change for Ligand (specifically, it knows to divide by Avogadro's number and the cellular volume).
I hope this example clears things up. If not, I recommend making a new MATLAB Answers question that provides the exact reactions and rates that you are having trouble converting. It would also be helpful if you can attach a file that allows us to see how you've currently implemented the model.
-Arthur

커뮤니티

더 많은 답변 보기:  SimBiology Community

카테고리

Help CenterFile Exchange에서 Scan Parameter Ranges에 대해 자세히 알아보기

제품

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by