Русская Википедия:Физиологически обоснованное фармакокинетическое моделирование (РВРК)
Физиологически обоснованное фармакокинетическое моделирование (Physiologically based pharmacokinetic modelling, PBPK) — это метод математического моделирования позволяющий прогнозировать всасывание, распределение, метаболизм и выведение (ADME) синтетических или природных химических веществ в организме человека и животных. PBPK-моделирование используют в фармацевтических исследованиях и разработке лекарственных препаратов, а также для оценки риска для здоровья применительно к косметическим средствам или химической продукции общего назначения.
Механистичность PBPK-моделей обеспечивается за счет математического преобразования анатомических, физиологических, физических и химических характеристик явлений, задействованных в сложных процессах ADME. Несмотря на значительную степень остаточного упрощения и эмпиризма, эти модели обладают более широкой областью применимости, чем классические фармакокинетические модели, основанные на эмпирических функциях. Модели PBPK можно использовать исключительно в прогностических целях, однако с появлением методов байесовской статистики для работы со сложными моделями стало возможно их применение и в других областях (например, в статистическом анализе) [1] Это применимо как к оценке риска токсичности, так и к разработке лекарственных препаратов.
В PBPK-моделях изначально учитывается анатомическая и физиологическая структура организма, а также в определенной степени биохимические параметры. Обычно это мультикомпартментные модели, где компартменты соответствуют заранее определенным органам или тканям, а взаимосвязи отражают циркуляцию крови или лимфы (реже — диффузию). Можно разработать систему дифференциальных уравнений для концентрации или количества вещества в каждом компартменте модели. Параметры этой системы будут соответствовать циркуляции крови, скорости вентиляции легких, объемам органов и прочим показателям, для которых имеются данные в научных публикациях. Полученные характеристики организма будут упрощенными, поэтому необходимо найти баланс между сложностью и простотой. Такие модели не только позволяют использовать априорную информацию о значениях параметров, но и упрощают межвидовое преобразование или экстраполяцию данных с одного способа применения на другой (например, с ингаляционного введения на прием внутрь). На рисунке справа представлен пример семикомпартментной PBPK-модели, позволяющей охарактеризовать трансформацию различных растворителей в организме млекопитающего.
История
Первая фармакокинетическая модель, описанная в научной литературе [2], фактически представляла собой PBPK-модель. Однако необходимые вычисления в то время было сложно реализовать. Затем акцент сместился на более простые модели [3], для которых можно было получить аналитические решения (такие решения представляли собой суммы экспоненциальных членов, что привело к дальнейшему упрощению). В начале 1970-х годов с появлением компьютеров и алгоритмов численного интегрирования интерес к физиологическим моделям возобновился [4][5] Для веществ со сложной кинетикой или в случаях, когда требовались межвидовые экстраполяции, простых моделей было недостаточно. Это обусловило необходимость дальнейших исследований физиологических моделей [6][7][8]. К 2010 году PBPK-модели и их применение были описаны в сотнях научных публикаций. В настоящее время деятельность по меньшей мере двух частных компаний основана на опыте в данной области.
Построение модели PBPK
Уравнения модели основаны на принципах переноса массы, динамики жидкостей и биохимии. Это позволяет смоделировать трансформацию вещества в организме [9]. Компартменты, или камеры, модели обычно определяют путем группировки органов или тканей со схожими показателями скорости перфузии крови и содержания липидов (то есть органов со схожими профилями зависимости концентрации химического вещества от времени). Каналы входа (легкие, кожа, кишечный тракт и прочие), каналы выхода (почки, печень и прочие) и органы-мишени при изучении терапевтического эффекта или токсичности часто оценивают отдельно. Костную ткань можно исключить из модели, если оцениваемое вещество в ней не распределяется. Взаимосвязи между компартментами модели соответствуют физиологии (например, кровь на выходе из кишечника поступает в печень и т. д.).
Основные уравнения переноса
Скорость распределения лекарственного препарата в ткани может быть ограничена перфузией или проницаемостью [10][11]. Кинетика, ограниченная скоростью перфузии, применяется, когда тканевые мембраны не являются барьером для диффузии. При этом кровоток служит ограничивающим фактором для распределения лекарственного препарата в различных клетках организма (исходя из того, что препарат переносится в основном кровью, что наблюдается во многих случаях). Это, как правило, применимо к низкомолекулярным липофильным лекарственным средствам. При ограничении, обусловленном перфузией, мгновенная скорость поступления определенного количества препарата в компартмент равна объемной скорости потока крови, проходящей через орган, умноженной на концентрацию препарата в поступающей крови. В этом случае для типичного компартмента i дифференциальное уравнение для количества вещества Qi, определяющее скорость изменения этого количества, будет следующим:
<math>{dQ_i \over dt} = F_i (C_{art} - {{Q_i} \over {P_i V_i}})</math>,
где Fi — кровоток (обозначенный как Q на рисунке выше); Cart — концентрация в поступающей артериальной крови; Pi — коэффициент распределения «ткань — кровь» и Vi — объем компартмента i.
Таким образом, полный набор дифференциальных уравнений для семикомпартментной модели, рассмотренной выше, можно представить в виде следующей таблицы.
| Ткань | Дифференциальное уравнение |
|---|---|
| Кишечник | <math>{dQ_g \over dt} = F_g (C_{art} - {{Q_g} \over {P_g V_g}})</math> |
| Почка | <math>{dQ_k \over dt} = F_k (C_{art} - {{Q_k} \over {P_k V_k}})</math> |
| Ткани с низкой перфузией (мышцы и кожа) | <math>{dQ_p \over dt} = F_p (C_{art} - {{Q_p} \over {P_p V_p}})</math> |
| Головной мозг | <math>{dQ_b \over dt} = F_b (C_{art} - {{Q_b} \over {P_b V_b}})</math> |
| Сердце и легкие | <math>{dQ_h \over dt} = F_h (C_{art} - {{Q_h} \over {P_h V_h}})</math> |
| Поджелудочная железа | <math>{dQ_{pn} \over dt} = F_{pn} (C_{art} - {{Q_{pn}} \over {P_{pn} V_{pn}}})</math> |
| Печень | <math>{dQ_l \over dt} = F_a C_{art}
+ F_g ({{Q_g} \over {P_g V_g}})
+ F_{pn} ({{Q_{pn}} \over {P_{pn} V_{pn}}})
- (F_a + F_g + F_{pn}) ({{Q_l} \over {P_l V_l}}) </math>
|
Эти уравнения включают только условия переноса, без учета входных и выходных параметров. Их можно смоделировать с использованием специальных членов уравнения, как показано ниже.
Входные данные для моделирования
Входные данные при моделировании необходимы, чтобы получить информативное описание фармакокинетики химического вещества. Далее приведены примеры составления соответствующих уравнений.
Прием внутрь
При оценке болюсной дозы препарата для приема внутрь (например, при приеме таблетки внутрь) всасывание первого порядка — широко распространенное допущение. В этом случае уравнение для кишечника дополняют входным параметром с константой скорости всасывания Ka:
<math>{dQ_g \over dt} = F_g (C_{art} - {{Q_g} \over {P_g V_g}}) + K_a Q_{ing}</math>.
Для этого требуется определить уравнение для количества препарата, принятого внутрь и присутствующего в просвете кишечника:
<math>{dQ_{ing} \over dt} = - K_a Q_{ing}</math>.
При отсутствии компартмента модели для кишечника в качестве входного параметра можно использовать печень. Однако в этом случае возможна некорректная оценка локального метаболизма в кишечнике. Для моделирования аппроксимативно непрерывного всасывания (например, с помощью питьевой воды) можно использовать скорость всасывания нулевого порядка (где значение Ring представлено в единицах массы в зависимости от времени):
<math>{dQ_g \over dt} = F_g (C_{art} - {{Q_g} \over {P_g V_g}}) + R_{ing}</math>.
Можно использовать и более сложную модель всасывания в кишечнике. В этих моделях дополнительные компартменты соответствуют различным участкам просвета и тканей кишечника. Можно также учитывать такие показатели, как рН кишечника, время транзита и наличие активных переносчиков [12].
Кожное депо
Всасывание химического вещества, нанесенного на кожу, также можно смоделировать с использованием членов первого порядка. В этом случае рекомендуется оценивать кожу отдельно от других тканей, чтобы достоверно различить участки кожи, подверженные и не подверженные воздействию, а также отличить жизнеспособную кожу (дерма и эпидермис) от рогового слоя (верхний слой кожи, подвергающийся воздействию). Данный подход используется в работе Bois F., Diaz Ochoa J.G., Gajewska M., Kovarich S., Mauch K., Paini A., Péry A., Sala Benito J.V., Teng S., Worth A. Multiscale modelling approaches for assessing cosmetic ingredients safety. Toxicology (в печати; doi: 10.1016/j.tox.2016.05.026).
Между роговым слоем, не подвергающимся воздействию, и находящейся под ним жизнеспособной кожей происходит обмен путем диффузии:
<math>{dQ_{{sc}_{u}} \over dt} = K_p \times S_s \times (1 - f_{S_{e}}) \times ({Q_{s_u} \over {P_{sc} V_{{sc}_{u}}}} - C_{{sc}_u})</math>,
где Kp — коэффициент распределения; Ss — общая площадь поверхности кожи; fSe — доля площади поверхности кожи, подвергшейся воздействию. Для жизнеспособной кожи, не подвергающейся воздействию:
<math> {dQ_{s_u} \over dt} = F_s (1 - f_{S_{e}}) (C_{art} - {{Q_{s_u}} \over {P_s V_{s_u}}}) - {dQ_{{sc}_{u}} \over dt} </math>.
Для рогового слоя кожи, подвергающегося воздействию:
<math>{dQ_{{sc}_{e}} \over dt} = K_p \times S_s \times f_{S_{e}} \times ({Q_{s_e} \over {P_{sc} V_{{sc}_e}}} - C_{{sc}_e})</math>.
Для жизнеспособной кожи, подвергающейся воздействию:
<math> {dQ_{s_e} \over dt} = F_s f_{S_{e}} (C_{art} - {{Q_{s_e}} \over {P_s V_{s_e}}}) - {dQ_{{sc}_{e}} \over dt} </math>.
Поступление dt(QSkin_u) и dt(QSkin_e) из артериальной крови и обратно в венозную кровь.
Внутривенное введение
Внутривенное введение — распространенный способ применения лекарственных препаратов в клинической практике.
Ингаляционное введение
Ингаляционное введение осуществляется через легкие, и этот процесс практически неотделим от выдыхания.
Моделирование метаболизма
Существует несколько способов моделирования метаболизма. В некоторых моделях предпочтительно использовать линейную скорость выведения. Для этого можно воспользоваться простым дифференциальным уравнением. Для получения более точных результатов обычно применяют уравнение Михаэлиса — Ментен:
- <math> v = \frac{d [P]}{d t} = \frac{ V_\max {[S]}}{K_m + [S]} </math>.
Области применения PBPK-моделирования
PBPK-модели по аналогии со многими другими моделями являются компартментными. Тем не менее они обладают рядом преимуществ по сравнению с так называемыми классическими фармакокинетическими моделями, которые в меньшей степени основаны на физиологии. PBPK-модели можно первоначально использовать для абстрактного представления и последующего согласования разрозненных данных (полученных в физико-химических или биохимических исследованиях, фармакологических или токсикологических исследованиях iin vitro или in vivo и проч.). Они также позволяют оценить внутренние концентрации химических веществ или их метаболитов в организме, в частности возможна оценка терапевтических или токсических концентраций таких химических веществ или метаболитов в местах их воздействия. Наконец такие модели способствуют интерполяции и экстраполяции данных для следующих показателей:
- дозы: например с высоких концентраций, обычно используемых в лабораторных экспериментах, на концентрации, наблюдаемые в окружающей среде;
- продолжительность воздействия: например с непрерывного на прерывистое или с однократного на многократное воздействие;
- способы применения: например с ингаляционного введения на прием внутрь;
- биологический вид: например преобразование данных для грызунов в данные для человека, до первого применения лекарственного препарата в клиническом исследовании или в случаях, когда эксперименты с участием человека считаются неэтичными (например, когда соединение токсично и не оказывает терапевтического эффекта);
- группы пациентов: например мужчины и женщины, взрослые пациенты и дети, небеременные и беременные женщины и проч.;
- с условий iin vitro на условиях in vivo [13].
Некоторые из этих экстраполяций считаются параметрическими: для экстраполяции требуются только изменения входных данных или значений параметров (обычно применительно к экстраполяции значений дозы и времени). Другие экстраполяции считаются непараметрическими, то есть требуют изменения самой структуры модели (например, при экстраполяции данных на беременную женщину следует добавить уравнения для плода).
В связи с механистической основой PBPK-моделей, еще одна потенциальная область применения PBPK-моделирования — это проверка гипотез. Например, если пероральная биодоступность лекарственного соединения оказывается ниже ожидаемой, возможна оценка различных структур (т. е. гипотез) и значений параметров модели, чтобы определить, какие модели и (или) параметры обеспечивают наилучшее соответствие наблюдаемым данным. Если наилучшее соответствие наблюдается для гипотезы, согласно которой низкая биодоступность обусловлена метаболизмом в кишечнике, то результаты PBPK‑моделирования лучше подтверждают эту гипотезу по сравнению с другими оцениваемыми гипотезами.
Таким образом, PBPK-моделирование можно использовать, помимо прочего, для оценки роли опосредованного переносчиками переноса, насыщения клиренса, кишечно-печеночной рециркуляции исходного соединения, внепеченочной/внекишечной элиминации, растворимости in vivo, превышающей расчетное значение in vitro, лекарственно-индуцированного замедления эвакуации содержимого желудка, снижения всасывания в кишечнике и местных изменений всасывания в кишечнике [14].
Ограничения и дополнительные возможности применения PBPK-моделирования
Каждая методика моделирования имеет свои сильные и слабые стороны. PBPK‑моделирование не исключение. Одно из ограничений — возможность использования большого количества параметров, некоторые из которых могут быть взаимосвязанными. Это может отрицательно повлиять на идентифицируемость параметров или привести к их избыточности. Тем не менее существует возможное (и широко применяемое) решение: прямое моделирование взаимосвязей между параметрами (например, нелинейные отношения между возрастом, массой тела, объемами органов и циркуляцией крови).
После присвоения численных значений каждому параметру PBPK-модели обычно используют специализированное или общее компьютерное программное обеспечение для численного интегрирования набора обычных дифференциальных уравнений, подобных описанным выше, чтобы рассчитать численное значение для каждого компартмента модели при заданных значениях времени. Иногда такие уравнения включают только линейные функции каждого значения для компартмента или к ним применяются ограничивающие условия (например, когда входные значения очень низкие), гарантирующие точное приближение такой линейности. В таких случаях эти уравнения можно решать аналитически для получения уравнений в явном виде (или, в случае ограничивающих условий, очень точных приближений) для средневзвешенного по времени (TWA) значения для каждого компартмента как функции значения TWA каждого указанного входного значения (см., например, [15][16]).
С одной стороны, PBPK-модели могут опираться на модели прогнозирования химических свойств (модели количественного соотношения «структура — активность» (QSAR) или прогностические модели химических свойств). Например, модели QSAR можно использовать для оценки коэффициентов распределения. Они также применимы к моделям системной биологии метаболических путей, но не заменяют их. Кроме того, они схожи с моделями физиома, но не предназначены для детального моделирования физиологических функций за пределами циркуляции жидкостей. В целом четыре типа моделей, упомянутые выше, могут взаимно дополнять друг друга при интеграции [17].
Рекомендации
- ↑ Шаблон:Cite journal
- ↑ Шаблон:Cite journal
- ↑ Шаблон:Cite book
- ↑ Шаблон:Cite journal
- ↑ Шаблон:Cite journal
- ↑ Шаблон:Cite journal
- ↑ Шаблон:Cite journal
- ↑ Шаблон:Cite journal
- ↑ Шаблон:Cite book
- ↑ Шаблон:Cite journal
- ↑ Шаблон:Cite journal
- ↑ Шаблон:Cite journal
- ↑ Шаблон:Cite journal
- ↑ Шаблон:Cite journal
- ↑ Шаблон:Cite journal
- ↑ Шаблон:Cite journal
- ↑ Шаблон:Cite journal
Дополнительные ссылки:
Форумы
- Ecotoxmodels — это веб-сайт, посвященный математическим моделям в экотоксикологии.
Программное обеспечение
Специальное программное обеспечение:
- БиоДМЕТ
- ГастроПлюс
- Максим2
- ПК-Сим
- ПКквест
- PSE: gCOAS
- Симцип Симулятор
- Верстак ADME
Общее программное обеспечение:
- АДАПТ 5
- Беркли Мадонна
- COPASI: Симулятор биохимической системы
- Эколего
- Бесплатное программное обеспечение для моделирования : GNU MCSIM
- Октава ГНУ
- Matlab PottersWheel
- Модельмейкер
- Физиолаб
- Пакет R deSolve
- СААМ II
- Phoenix WinNonlin/NLME/IVIVC/Trial Simulator
