El modelado farmacocinético basado en la fisiología (PBPK) es una técnica de modelado matemático para predecir la absorción, distribución, metabolismo y excreción (ADME) de sustancias químicas sintéticas o naturales en humanos y otras especies animales. El modelado PBPK se utiliza en la investigación farmacéutica y el desarrollo de fármacos, y en la evaluación de riesgos para la salud de los cosméticos o productos químicos en general.
Los modelos PBPK intentan ser mecanicistas al transcribir matemáticamente descripciones anatómicas, fisiológicas, físicas y químicas de los fenómenos involucrados en los complejos procesos ADME. Un alto grado de simplificación residual y empirismo todavía está presente en esos modelos, pero tienen un dominio de aplicabilidad extendido comparado con el de los modelos farmacocinéticos clásicos basados en funciones empíricas . Los modelos PBPK pueden tener usos puramente predictivos, pero otros usos, como la inferencia estadística, han sido posibles gracias al desarrollo de herramientas estadísticas bayesianas capaces de manejar modelos complejos. [1] Esto es cierto tanto para la evaluación del riesgo de toxicidad como para el desarrollo de fármacos terapéuticos.
Los modelos PBPK intentan basarse a priori en la estructura anatómica y fisiológica del cuerpo y, en cierta medida, en la bioquímica. Suelen ser modelos multicompartimentales , con compartimentos correspondientes a órganos o tejidos predefinidos, con interconexiones correspondientes a flujos sanguíneos o linfáticos (más raramente a difusiones). Se puede escribir un sistema de ecuaciones diferenciales para la concentración o cantidad de sustancia en cada compartimento, y sus parámetros representan flujos sanguíneos, tasa de ventilación pulmonar, volúmenes de órganos, etc., para lo cual hay información disponible en publicaciones científicas. De hecho, la descripción que hacen del cuerpo es simplificada y es necesario encontrar un equilibrio entre complejidad y simplicidad. Además de la ventaja de permitir el reclutamiento de información a priori sobre los valores de los parámetros, estos modelos también facilitan las transposiciones interespecies o la extrapolación de un modo de administración a otro ( por ejemplo , inhalación a oral). En la figura de la derecha se da un ejemplo de un modelo PBPK de 7 compartimentos, adecuado para describir el destino de muchos disolventes en el cuerpo de los mamíferos.
El primer modelo farmacocinético descrito en la literatura científica [2] fue, de hecho, un modelo PBPK. Sin embargo, dio lugar a cálculos que en aquel momento eran intratables. El enfoque se desplazó entonces a modelos más simples [3] , para los que se podían obtener soluciones analíticas (tales soluciones eran sumas de términos exponenciales, lo que dio lugar a mayores simplificaciones). La disponibilidad de ordenadores y algoritmos de integración numérica marcó un renovado interés en los modelos fisiológicos a principios de la década de 1970 [4] [5] Para las sustancias con cinéticas complejas, o cuando se requerían extrapolaciones entre especies, los modelos simples eran insuficientes y la investigación continuó con los modelos fisiológicos [6] [7] [8] En 2010, cientos de publicaciones científicas habían descrito y utilizado modelos PBPK, y al menos dos empresas privadas han basado su negocio en su experiencia en esta área.
Las ecuaciones del modelo siguen los principios de transporte de masa, dinámica de fluidos y bioquímica para simular el destino de una sustancia en el cuerpo. [9] Los compartimentos generalmente se definen agrupando órganos o tejidos con una tasa de perfusión sanguínea y un contenido de lípidos similares ( es decir, órganos para los cuales los perfiles de concentración de sustancias químicas frente al tiempo serán similares). Los puertos de entrada (pulmón, piel, tracto intestinal...), los puertos de salida (riñón, hígado...) y los órganos diana para el efecto terapéutico o la toxicidad a menudo se dejan separados. El hueso se puede excluir del modelo si la sustancia de interés no se distribuye en él. Las conexiones entre compartimentos siguen la fisiología ( p. ej. , el flujo sanguíneo en la salida del intestino va al hígado, etc. )
La distribución de un fármaco en un tejido puede verse limitada por la velocidad de perfusión o por la permeabilidad. [10] [11] La cinética limitada por la velocidad de perfusión se aplica cuando las membranas tisulares no presentan ninguna barrera a la difusión. El flujo sanguíneo, suponiendo que el fármaco se transporta principalmente por sangre, como suele ser el caso, es entonces el factor limitante de la distribución en las distintas células del cuerpo. Esto suele ser cierto para fármacos lipofílicos pequeños. En caso de limitación por perfusión, la velocidad instantánea de entrada de la cantidad de fármaco en un compartimento es simplemente igual a la velocidad de flujo volumétrico (sanguíneo) a través del órgano multiplicada por la concentración de sangre entrante. En ese caso, para un compartimento genérico i , la ecuación diferencial para la cantidad Q i de sustancia, que define la velocidad de cambio de esta cantidad, es:
donde F i es el flujo sanguíneo (indicado como Q en la figura anterior), C es la concentración de sangre arterial entrante, P i es el coeficiente de partición entre tejido y sangre y V i es el volumen del compartimento i .
Por lo tanto, un conjunto completo de ecuaciones diferenciales para el modelo de siete compartimentos mostrado arriba podría darse mediante la siguiente tabla:
Tejido | Ecuación diferencial |
---|---|
Intestino | |
Riñón | |
Tejidos mal perfundidos (músculos y piel) | |
Cerebro | |
Corazón y pulmón | |
Páncreas | |
Hígado |
Las ecuaciones anteriores incluyen únicamente términos de transporte y no tienen en cuenta las entradas ni las salidas. Estas pueden modelarse con términos específicos, como se muestra a continuación.
Es necesario modelar las entradas para obtener una descripción significativa de la farmacocinética de una sustancia química. Los siguientes ejemplos muestran cómo escribir las ecuaciones correspondientes.
Cuando se trata de una dosis oral en bolo ( por ejemplo , la ingestión de un comprimido), la absorción de primer orden es una suposición muy común. En ese caso, la ecuación intestinal se amplía con un término de entrada, con una constante de velocidad de absorción K a :
Esto requiere definir una ecuación para la cantidad ingerida y presente en el lumen intestinal:
En ausencia de un compartimento intestinal, la entrada puede realizarse directamente en el hígado. Sin embargo, en ese caso, el metabolismo local en el intestino puede no estar correctamente descrito. El caso de absorción aproximadamente continua ( por ejemplo, a través del agua potable) puede modelarse mediante una tasa de absorción de orden cero (aquí R en unidades de masa a lo largo del tiempo):
Se puede utilizar un modelo de absorción intestinal más sofisticado. En esos modelos, compartimentos adicionales describen las distintas secciones del lumen y el tejido intestinal. Se pueden tener en cuenta el pH intestinal, los tiempos de tránsito y la presencia de transportadores activos. [12]
La absorción de una sustancia química depositada sobre la piel también se puede modelar utilizando términos de primer orden. En ese caso, lo mejor es separar la piel de los demás tejidos, para diferenciar aún más la piel expuesta de la no expuesta y diferenciar la piel viable (dermis y epidermis) del estrato córneo (la capa superior de la piel realmente expuesta). Este es el enfoque adoptado en [Bois F., Diaz Ochoa JG Gajewska M., Kovarich S., Mauch K., Paini A., Péry A., Sala Benito JV, Teng S., Worth A., en prensa, Enfoques de modelado multiescala para evaluar la seguridad de los ingredientes cosméticos, Toxicología. doi: 10.1016/j.tox.2016.05.026]
El estrato córneo no expuesto simplemente intercambia con la piel viable subyacente por difusión:
donde es el coeficiente de partición, es la superficie total de la piel, la fracción de la superficie de la piel expuesta, ...
Para la piel viable no expuesta:
Para el estrato córneo de la piel expuesto:
Para la piel viable expuesta:
dt(QSkin_u) y dt(QSkin_e) se alimentan de sangre arterial y regresan a la sangre venosa.
Se han publicado modelos de difusión más complejos [referencia para agregar].
La inyección intravenosa es una vía de administración clínica común. (a completar)
La inhalación se produce a través del pulmón y es difícilmente disociable de la exhalación (a completar).
Existen varias formas de modelar el metabolismo. En algunos modelos, se prefiere una tasa de excreción lineal, lo que se puede lograr con una ecuación diferencial simple. De lo contrario, una ecuación de Michaelis-Menten , como la que se muestra a continuación, suele ser adecuada para obtener un resultado más preciso.
Los modelos PBPK son modelos compartimentados como muchos otros, pero tienen algunas ventajas sobre los llamados modelos farmacocinéticos "clásicos", que tienen menos fundamento fisiológico. Los modelos PBPK se pueden utilizar, en primer lugar, para abstraer y, en última instancia, conciliar datos dispares (de experimentos fisicoquímicos o bioquímicos, experimentos farmacológicos o toxicológicos in vitro o in vivo , etc. ). También brindan acceso a las concentraciones corporales internas de sustancias químicas o sus metabolitos, y en particular en el lugar de sus efectos, ya sean terapéuticos o tóxicos. Finalmente, también ayudan a la interpolación y extrapolación de conocimientos entre:
Algunas de estas extrapolaciones son "paramétricas": solo se necesitan cambios en los valores de entrada o de los parámetros para lograr la extrapolación (este suele ser el caso de las extrapolaciones de dosis y tiempo). Otras son "no paramétricas" en el sentido de que se necesita un cambio en la estructura del modelo en sí ( por ejemplo , al extrapolar a una mujer embarazada, se deben agregar ecuaciones para el feto).
Debido a la base mecanicista de los modelos PBPK, otro uso potencial de los modelos PBPK es la prueba de hipótesis. Por ejemplo, si un compuesto farmacológico mostró una biodisponibilidad oral menor a la esperada, se pueden evaluar varias estructuras de modelos (es decir, hipótesis) y valores de parámetros para determinar qué modelos y/o parámetros se ajustan mejor a los datos observados. Si la hipótesis de que el metabolismo en los intestinos era responsable de la baja biodisponibilidad arrojaba el mejor ajuste, entonces los resultados del modelo PBPK apoyaban esta hipótesis sobre las otras hipótesis evaluadas.
Como tal, el modelado PBPK se puede utilizar, entre otras cosas , para evaluar la participación del transporte mediado por transportador, la saturación del aclaramiento, la recirculación enterohepática del compuesto original, la eliminación extrahepática/extraintestinal; mayor solubilidad in vivo que la predicha in vitro ; retrasos en el vaciamiento gástrico inducidos por fármacos; pérdida intestinal y variación regional en la absorción intestinal. [15]
Cada tipo de técnica de modelado tiene sus puntos fuertes y sus limitaciones. El modelado PBPK no es una excepción. Una limitación es la posibilidad de utilizar una gran cantidad de parámetros, algunos de los cuales pueden estar correlacionados. Esto puede dar lugar a problemas de identificación y redundancia de los parámetros. Sin embargo, es posible (y se hace habitualmente) modelar explícitamente las correlaciones entre parámetros (por ejemplo, las relaciones no lineales entre la edad, la masa corporal, los volúmenes de los órganos y los flujos sanguíneos).
Después de asignar valores numéricos a cada parámetro del modelo PBPK, se suele utilizar un software informático especializado o general para integrar numéricamente un conjunto de ecuaciones diferenciales ordinarias como las descritas anteriormente, con el fin de calcular el valor numérico de cada compartimento en valores de tiempo específicos (véase Software). Sin embargo, si dichas ecuaciones implican solo funciones lineales de cada valor compartimental, o en condiciones límite (por ejemplo, cuando los valores de entrada siguen siendo muy pequeños) que garantizan que dicha linealidad se aproxime estrechamente, dichas ecuaciones se pueden resolver analíticamente para producir ecuaciones explícitas (o, en esas condiciones límite, aproximaciones muy precisas) para el valor promedio ponderado en el tiempo (TWA) de cada compartimento como una función del valor TWA de cada entrada especificada (véase, por ejemplo, [16] [17] ).
Los modelos PBPK pueden basarse en modelos de predicción de propiedades químicas ( modelos QSAR o modelos de química predictiva) por un lado. Por ejemplo, los modelos QSAR pueden utilizarse para estimar coeficientes de partición. También se extienden a los modelos de biología de sistemas de vías metabólicas, pero no están destinados a reemplazarlos. También son paralelos a los modelos fisiológicos , pero no apuntan a modelar funciones fisiológicas más allá de la circulación de fluidos en detalle. De hecho, los cuatro tipos de modelos anteriores pueden reforzarse entre sí cuando se integran. [18]
Referencias adicionales:
Software dedicado:
Software general: