Фон

Предположим, что у вас есть модель с одним отделением с видом под названием 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 автоматически создается и добавляется в отсек. The 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);