Modelos de estabilidad acelerada comparados — Arrhenius, ASAP, AccelStab y decisión probabilística

1. Resumen ejecutivo

Cuatro familias de modelos dominan el trabajo de estabilidad acelerada y predictiva en farmacia: la ecuación clásica de Arrhenius, el Arrhenius corregido por humedad empleado por el Accelerated Stability Assessment Program (ASAP), el enfoque cinético de Šesták–Berggren implementado en el paquete de R de código abierto AccelStab, y el enfoque Arrhenius probabilístico + Monte-Carlo que Synlogica Terminus Quality utiliza para la decisión sobre excursiones.

A menudo se presentan como competidores. No lo son. Tres de ellos (Arrhenius, ASAP, AccelStab) responden a una pregunta de desarrollo — «¿qué vida útil puedo declarar a partir de un estudio breve y forzado?» — mientras que Synlogica responde a una pregunta operativa — «dada la cinética ya caracterizada de este producto y este registro real del logger, ¿es liberable este lote y con qué confianza?» Este documento expone con claridad las ecuaciones y los supuestos de cada modelo, los compara y muestra dónde encaja cada uno.

2. La base compartida — Arrhenius

Los cuatro enfoques heredan la dependencia con la temperatura de Arrhenius de la constante de velocidad de degradación k:

Lo que añaden sobre esta ecuación de base — un término de humedad, una forma de degradación no lineal, un modelo de incertidumbre — es precisamente lo que los distingue.

3. Arrhenius clásico y MKT

El enfoque tradicional ajusta k a varias temperaturas elevadas, extrapola a la condición de almacenamiento y reporta una única vida útil. Para exposiciones no isotérmicas (p. ej. una excursión durante la distribución) suele colapsarse en la Temperatura Cinética Media (MKT):

Supuestos: una única reacción de degradación dominante y limitante de la velocidad; E_a constante en todo el rango de temperatura; cinética simple (de orden cero o de primer orden); humedad irrelevante o mantenida constante; la temperatura del sensor igual a la del producto.

Límites: produce una estimación puntual que oculta la incertidumbre en E_a, en la calibración del sensor y en la inhomogeneidad espacial; la MKT descarta la forma de la excursión; sin término de humedad; falla cuando el mecanismo de degradación cambia con la temperatura.

4. ASAP — Arrhenius corregido por humedad

El Accelerated Stability Assessment Program (Waterman et al.) amplía Arrhenius con un término de humedad explícito y log-lineal, y emplea un diseño isoconversional: mantener los productos en varias combinaciones de temperatura/humedad relativa durante días o semanas, medir el tiempo hasta un nivel fijo de degradación y ajustar:

La vida útil y su intervalo de confianza se propagan después mediante simulación de Monte-Carlo sobre los parámetros ajustados. ASAP sustenta las herramientas comerciales de Accelerated Predictive Stability (APS) y se utiliza ampliamente en presentaciones QbD.

Supuestos: un único punto final de degradación isoconversional; dependencia log-lineal con la humedad; equilibrio de humedad con el envase (abierto o definido); comportamiento de Arrhenius en todo el rango forzado; un límite de especificación conocido que define el «fallo».

Límites: requiere un diseño de estudio con humedad controlada; el término log-lineal de RH es una aproximación empírica; el planteamiento isoconversional supone que la vía de degradación no cambia entre condiciones.

5. AccelStab — cinética de Šesták–Berggren

AccelStab es un paquete de R de código abierto (CRAN; AccelStab en GitHub) que combina el término de temperatura de Arrhenius con la función de modelo de reacción de Šesták–Berggren, proporcionando un ajuste no lineal flexible y de forma cerrada a toda la curva de degradación, en lugar de un único punto isoconversional:

Extrapola la degradación tanto en tiempo como en temperatura y cuantifica la incertidumbre con un bootstrap paramétrico (parámetros cinéticos extraídos de una distribución t multivariante) y/o el método delta analítico, obteniendo intervalos de confianza y de predicción sobre la curva predicha.

Supuestos: la forma de Šesták–Berggren captura la verdadera forma de la degradación; dependencia de la temperatura según Arrhenius; las estimaciones de parámetros son aproximadamente t multivariante; puede añadirse un término de Arrhenius modificado por humedad donde esta sea relevante.

Límites: más parámetros (m, n) requieren más datos o más limpios para identificarlos; sigue siendo una herramienta de estudio de desarrollo, no un motor de decisión por envío; al ser de código abierto, el estado validado es responsabilidad del usuario.

6. Synlogica Terminus Quality — decisión probabilística sobre excursiones

Terminus Quality no intenta establecer una vida útil a partir de un estudio forzado. Toma la cinética de Arrhenius ya caracterizada del producto (E_a, A y su incertidumbre, procedentes del estudio de desarrollo ICH Q1A) y la aplica al registro real y no isotérmico del logger de un evento de distribución real para responder a una pregunta de liberación:

Cuando la degradación es no lineal, Terminus ajusta el mismo modelo de reacción de Šesták–Berggren empleado por AccelStab (dα/dt = k·f(α)) mediante mínimos cuadrados no lineales, y reporta intervalos de confianza y de predicción con un bootstrap paramétrico (parámetros cinéticos extraídos de una distribución t multivariante) y el método delta — además de la propagación por Monte-Carlo de E_a, A, la calibración del sensor y la incertidumbre espacial. El motor abarca así tanto el caso simple de Arrhenius integrado como el caso cinético completo de Šesták–Berggren, eligiendo según el producto.

El resultado es una recomendación ponderada por probabilidad — RELEASE / CONDITIONAL / HOLD con un intervalo de confianza declarado — empaquetada junto con cada entrada, versión del modelo y supuesto para auditoría. La persona (QA / QP) toma la decisión vinculante.

Supuestos: los parámetros de Arrhenius de la etapa de desarrollo son válidos y están fijados por versión; el registro del logger, tras la calibración y el cribado de huecos, representa al producto dentro de los límites espaciales declarados; la exposición acumulada (la integral) es la magnitud relevante. La humedad se incluye cuando los datos de estabilidad del producto respaldan un término de humedad.

Lo que no es: no es un sustituto de ASAP/AccelStab en desarrollo, y es consultivo por diseño — no firma la disposición del lote.

7. Comparación lado a lado

8. Supuestos y modos de fallo que comparten

Todos los modelos aquí inducirán a error si se violan sus supuestos. Los modos de fallo comunes:

• Cambio de mecanismo — si la vía de degradación limitante de la velocidad cambia entre las condiciones forzadas y las de almacenamiento, la extrapolación de Arrhenius (y todo lo que se construye sobre ella) deja de ser válida.
• Incertidumbre de E_a ignorada — un error del ±15% en la energía de activación puede desplazar una estimación de vida útil en meses; las estimaciones puntuales lo ocultan.
• Sensor ≠ producto — el error de calibración del logger y la inhomogeneidad espacial dentro de un envío implican que el perfil medido no es exactamente el del producto.
• Humedad ignorada en productos sensibles a la humedad — tanto el Arrhenius clásico como una integración ciega a la humedad sobrestiman la confianza.

9. Qué modelo, cuándo

• Fijar una vida útil registrada con rapidez, cuando la humedad importa: ASAP / APS.
• Modelar la curva de degradación completa, de forma abierta y reproducible: AccelStab (Šesták–Berggren).
• Comprobación interna rápida de plausibilidad: Arrhenius clásico / MKT — conociendo sus puntos ciegos.
• Decidir la liberación de un envío concreto tras una excursión: Synlogica Terminus Quality — consume el modelo de desarrollo y el registro real y devuelve una recomendación auditable y ponderada por probabilidad.

10. Referencias

1. Waterman K.C. et al. — The Application of the Accelerated Stability Assessment Program (ASAP) to Quality by Design (QbD) for Drug Product Stability. AAPS PharmSciTech (PMC3167263).
2. Humidity-corrected Arrhenius equation: the reference condition approach. Int. J. Pharmaceutics, 2016 (S0378517316300473).
3. AccelStab: Accelerated Stability Kinetic Modelling — paquete de CRAN y viñeta; código fuente en github.com/AccelStab/AccelStab (Šesták–Berggren + Arrhenius; intervalos por bootstrap paramétrico / método delta).
4. ICH Q1A(R2) Stability Testing; ICH Q1E Evaluation of Stability Data.

11. Conclusiones

• Los cuatro métodos comparten un núcleo de Arrhenius; difieren en el tratamiento de la humedad, la forma de degradación, el método de incertidumbre y el propósito.
• ASAP y AccelStab son herramientas de desarrollo para la vida útil; Synlogica es una herramienta operativa para la decisión de liberación por lote — son complementarias.
• Sea cual sea el modelo, el resultado defendible es una distribución con supuestos declarados, no un único número — y, para la liberación, un paquete de decisión auditable con la persona como autoridad de decisión.