SimBiology

关注

I'm trying to run this model with SimBiology but see an error message. BIOMD0000000268_url.xml SBML L2V4 representation of Reed2008_Glutathione_Metabolism Reed2008_Glutathione_Metabolism | BioModels (ebi.ac.uk) "piecewise" is a function in matlab but it seems SimBiology model bulider recognizes it as a variable... Any idea? --------- RULE Aminoacid_input = piecewise(breakfast,and(le(7,daytime),le(daytime,10)),lunch,and(le(12,daytime),le(daytime,15)),dinner,and(le(18,daytime),le(daytime,21)),fasting) error message: "Rule include varaibles that are unidentified." with red higllight for "piecewise" ---------------------------------

Hi, all I want to calculate the intestinal solubility of drug according to intestinal cyclodextrin concentration in a SimBiology model, so I compile a function like: function [IntestinalAGSolubility] = calculateSolubilityCD(CDConcentration,CDMW,AGMW) % This is a function to calculate solubility of a drug,AG, according to concentration % of cyclodextrin (CD). Equation is transformed from phase solubility test. % AGMW: AG molecular weight (g/mol); CDMW: CD molecular weight (g/mol); % CDConcentration (mg/mL); IntestinalAGSolubility (mg/mL) % 0.4653: slope; 0.4409: intercept (mmol/L) if (0.4653*CDConcentration*1000/CDMW+0.4409)*AGMW>2990.1 IntestinalAGSolubility=2990.1*1000; elseif (0.4653*CDConcentration*1000/CDMW+0.4409)*AGMW<=2990.1 IntestinalAGSolubility=(0.4653*CDConcentration*1000/CDMW+0.4409)*AGMW*1000; end This function is due to the fact that the solubility of AG and CD concentraion is consistent with a equation: S (mmol/L) =0.4653*CD (mmol/L)+0.4409, and maximum solubility is 2990.1 mg/L. Variables in function with their dimension are used in my SimBiology model, which is modified from "generic PBPK model". I have turned on the UnitConversion in program and I call this function in SimBiology APP like ColonSolubility=calculateSolubilityCD(Colon.CD,CDMW,AGMW). The line: IntestinalAGSolubility=2990.1*1000; and IntestinalAGSolubility=(0.4653*CDConcentration*1000/CDMW+0.4409)*AGMW*1000 are coded like this because only by this way I can get the exepexted result that presented in figure. However, value of 2990.1*1000 is not what I want to use in AG solubility (2.990 mg/mL). I am confused why I can not use IntestinalAGSolubility=2.990 in function (which will result in 1000-fold decrease in simulated solubility). Can anyone explain why this happens and do I make it right in other part of this function? Thanks a lot.