Фон

Предположим, что у вас есть модель с одним отсеком с разновидностью под названием drug, который представляет общую сумму препарата в теле. Препарат удален из тела через устранение первого порядка, представленное реакцией drug -> null с константой скорости устранения ke. Другими словами, концентрация препарата по сравнению с профилем времени следует за моноэкспоненциальным снижением $$C_t=C_0{e}^{-k_{e}t}$$, где $$C_t$$ концентрация препарата во время t, $$C_0$$ начальная концентрация, и ke является константой скорости устранения. Этот пример показывает, как настроить такую модель с одним отсеком и добавить дозу вливания на постоянном уровне 10 мг/час для суммы суммарной дозы 250 мг.

Создайте модель с одним отсеком

Создайте модель SimBiology под названием onecomp.

m1 = sbiomodel('onecomp');

Задайте устранение препарата от системы путем добавления реакции drug -> null в модель.

r1 = addreaction(m1,'drug -> null');

Разновидность drug автоматически создается и добавляется к отсеку. Разновидность null является зарезервированной разновидностью, которая действует как приемник в этой реакции.

Добавьте массовую акцию кинетический закон в реакцию. Этот кинетический закон задает устранение препарата, чтобы следовать за кинетикой первого порядка.

k1 = addkineticlaw(r1,'MassAction');

Задайте параметр уровня устранения ke и добавьте его в кинетический закон.

p1 = addparameter(k1,'ke','Value',1.0,'ValueUnits','1/hour');

Задайте параметр уровня ke как параметр форвардного курса реакции путем установки свойства ParameterVariableNames кинетического объекта k1 закона. Это позволяет SimBiology определять скорость реакции для реакции drug -> null.

k1.ParameterVariableNames = 'ke';

Настройте Дозу Вливания

Добавьте объект дозы в модель с помощью метода adddose. Задайте сумму дозы (Amount), цель дозы (TargetName) и уровень вливания (Rate). Также необходимо установить свойство Active объекта дозы к true так, чтобы доза была применена к модели во время симуляции.

d1 = adddose(m1,'InfusionDose');
d1.Amount = 250;
d1.TargetName = 'drug';
d1.Rate = 10;
d1.RateUnits = 'milligram/hour';
d1.Active = true;

Моделируйте модель

Измените время остановки симуляции на 48 часов, чтобы видеть полный курс времени.

cs = getconfigset(m1);
cs.StopTime = 48;
cs.TimeUnits = 'hour';
sd = sbiosimulate(m1);

Постройте результаты

Постройте концентрацию по сравнению с профилем времени препарата в системе.

sbioplot(sd);