Фон

Предположим, у вас есть однокамерная модель с видом по имени drug это представляет собой общее количество препарата в организме. Лекарственное средство выводится из организма посредством элиминации первого порядка, представляемой реакцией drug -> null, с постоянной скоростью исключения ke. Другими словами, концентрация лекарственного средства по сравнению с профилем времени следует за моноэкспоненциальным снижением $${}_{\mathrm{Ct}}{=}_{\mathrm{}}{}^{{}_{}}$$C0e-ket, $${\mathrm{где}}_{}$$Ct - концентрация лекарственного средства в момент времени t$${,}_{\mathrm{}}$$C0 - начальная концентрация, и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);