BIS-Based Model Reference Adaptive Control of Anesthesia Depth under Patient Variability and Time-Varying Delays

Resumen Técnico: Control Adaptativo de Modelo de Referencia Basado en BIS para la Profundidad de la Anestesia bajo Variabilidad del Paciente y Retrasos Variables en el Tiempo

Planteamiento del Problema
La administración de anestesia general requiere la regulación precisa de fármacos hipnóticos (propofol) y analgésicos (remifentanilo) para mantener la inconsciencia y la analgesia. Un desafío primordial en el control de bucle cerrado automatizado es la significativa variabilidad interpaciente en los parámetros farmacocinéticos/farmacodinámicos (PK/PD) y la presencia de retrasos fisiológicos variables en el tiempo en la distribución de los fármacos. Las estrategias de control convencionales, como el PID de ganancia fija o el Control Predictivo Basado en Modelo (MPC), a menudo dependen de la identificación de parámetros específicos del paciente de forma offline o asumen una distribución de fármacos instantánea. Estas limitaciones reducen su robustez en entornos clínicos reales donde la fisiología del paciente cambia (por ejemplo, variaciones en el gasto cardíaco) y los retrasos fluctúan. Además, la relación no lineal entre las concentraciones de los fármacos y la salida del Índice Bispectral (BIS), gobernada por una función de Hill con términos de interacción de fármacos, complica el diseño de controladores de retroalimentación lineal.

Metodología
Los autores proponen un marco de Control Adaptativo de Modelo de Referencia (MRAC) diseñado para regular la profundidad de la anestesia utilizando una infusión de doble fármaco de propofol y remifentanilo. La metodología integra los siguientes componentes:

1. Modelado de la Planta: El paciente se modela utilizando modelos PK estándar de tres compartimentos (Schnider) para ambos fármacos, extendidos para incluir retrasos variables en el tiempo ($\tau(t)$) en la transferencia intercompartimental. El componente farmacodinámico (PD) vincula las concentraciones en el sitio de efecto con la salida del BIS mediante una función de Hill no lineal que considera la interacción sinérgica entre el propofol y el remifentanilo.
2. Reconstrucción Inversa de Hill: Para abordar la no linealidad de la función de Hill y la ambigüedad del control de múltiples entradas, los autores introducen una estrategia de reconstrucción inversa simplificada. Al asumir que el efecto normalizado del remifentanilo es proporcional al efecto normalizado del propofol ($Y_r \approx \mu Y_p$), la compleja interacción de múltiples fármacos se reduce a una relación de una sola variable. Se aplica entonces una función de Hill inversa ($H^{-1}_\mu$) a la señal BIS medida para estimar la concentración del sitio de efecto de propofol linealizada ($\hat{Y}_p$), que sirve como la variable de retroalimentación.
3. Ley de Control Adaptativo: El esquema central del controlador es un MRAC diseñado para sistemas con grado relativo dos y retrasos desconocidos variables en el tiempo.
   • Modelo de Referencia: Un modelo de referencia de segundo orden ($W_m(s)$) define la respuesta transitoria y en estado estacionario deseada (objetivo de BIS de 50).
   • Estructura de Control: La entrada de control se genera mediante regresores filtrados derivados de la entrada y la salida estimada. Para manejar el grado relativo y asegurar la estabilidad, se introduce un filtro estable de primer orden $Z(s) = s + p_0$ para hacer que la dinámica del error sea estrictamente positiva real (SPR).
   • Adaptación y Robustez: Los parámetros del controlador se actualizan en línea mediante una ley adaptativa derivada de un funcional de Lyapunov-Krasovskii. Se incluye un término de robustez específico que involucra el signo del error de seguimiento y un término integral para compensar los retrasos variables en el tiempo desconocidos y las perturbaciones externas (por ejemplo, estímulos quirúrgicos).
4. Análisis de Estabilidad: La estabilidad teórica se demuestra mediante un funcional de Lyapunov-Krasovskii, demostrando que todas las señales de bucle cerrado permanecen acotadas y que el error de seguimiento converge asintóticamente a cero, siempre que la tasa de cambio del retraso sea menor que la unidad.

Contribuciones Clave
El artículo describe cinco contribuciones específicas:

• Independencia de Parámetros: El enfoque MRAC no requiere la identificación exacta de los parámetros PK/PD específicos del paciente ni una sintonización previa, a diferencia de los controladores PID o de eventos convencionales.
• Manejo Explícito de Retrasos: A diferencia de los marcos adaptativos existentes (por ejemplo, MRAC de orden fraccionario), este método incorpora explícitamente los retrasos fisiológicos variables en el tiempo en el diseño del control y el análisis de estabilidad, en lugar de depender únicamente de efectos de memoria de orden fraccionario.
• Compensación No Lineal: La reconstrucción inversa de Hill simplificada propuesta permite la estimación de una salida lineal a partir de la señal no lineal del BIS, facilitando el control de retroalimentación de salida adaptativo a pesar de las complejas interacciones de múltiples fármacos.
• Robustez: La ley adaptativa actualiza continuamente los parámetros en línea, asegurando la estabilidad contra la variabilidad interpaciente y las perturbaciones intraoperatorias sin necesidad de recalibrar el controlador.
• Seguridad Clínica: El controlador regula las tasas de infusión dentro de los rangos estándar clínicamente aceptados, evitando tanto la sobredosificación como la anestesia insuficiente.

Resultados
El método propuesto fue validado mediante simulaciones en una población heterogénea de 12 pacientes virtuales con diversas características biométricas (edad, peso, altura, género) y parámetros PK/PD generados aleatoriamente.

• Fase de Inducción: El controlador logró llevar el BIS desde la línea base hasta el objetivo de 50 en todos los pacientes. El Tiempo hacia el Objetivo (TT) osciló entre 104 y 193 segundos (media ~145s), y el Tiempo de Asentamiento (ST10) varió entre 141 y 360 segundos (media ~263s). El BIS-NÁDIR (BIS mínimo) se mantuvo por encima de 40 en todos los casos, evitando una profundidad excesiva.
• Fase de Mantenimiento: Bajo perturbaciones quirúrgicas simuladas (entradas de doble escalón), el controlador restauró rápidamente el BIS al rango objetivo. Los tiempos de recuperación ante perturbaciones positivas y negativas fueron cortos (media ~0.6 min), y el BIS se mantuvo dentro del rango 40–60 con un sobreimpulso mínimo (<5%) y desviaciones pequeñas (±2).
• Administración de Fármacos: Los perfiles de infusión para el propofol y el remifentanilo se mantuvieron dentro de los límites clínicos típicos tanto para la fase de inducción como para la de mantenimiento.
• Comparación: El MRAC propuesto superó o igualó a las estrategias existentes (MRAC de orden fraccionario, PID optimizado, MPC no lineal) al lograr un seguimiento competitivo sin requerir parámetros de modelos específicos del paciente o sintonización offline.

Significancia y Reivindicaciones
Los autores afirman que este estudio proporciona una solución prometedora para la gestión de la anestesia de bucle cerrado robusta e individualizada. La principal significancia radica en la capacidad del controlador para operar eficazmente sin la identificación previa de parámetros específicos del paciente, abordando así una barrera importante para el despliegue clínico de los sistemas de anestesia automatizados. Al modelar y compensar explícitamente los retrasos variables en el tiempo y las interacciones no lineales de los fármacos, el marco MRAC propuesto ofrece un enfoque teóricamente fundamentado y escalable que mantiene la seguridad y la estabilidad del paciente a través de diversas condiciones fisiológicas. El trabajo sugiere que tal sistema podría mejorar la seguridad del paciente, reducir el consumo de fármacos y mejorar los tiempos de recuperación en comparación con los métodos manuales o automatizados menos adaptativos.