Dynamic behavior of a two-dimensional discrete-time HIV infection model: stability and bifurcation analysis

El virus de la inmunodeficiencia humana, conocido como VIH, es un patógeno implacable que ataca el sistema inmunológico del cuerpo, dirigiéndose específicamente a las células responsables de coordinar la defensa contra las enfermedades. Sin tratamiento, esta infección puede progresar a una etapa en la que el cuerpo ya no puede combatir enfermedades comunes. Durante décadas, los científicos han utilizado modelos matemáticos para comprender cómo se comporta el virus dentro del cuerpo humano. Tradicionalamente, estos modelos se han construido sobre la idea del tiempo continuo, tratando la infección como un flujo suave e ininterrumpido donde los cambios ocurren en cada instante. Sin embargo, en el mundo real, los datos médicos rara vez se recolectan en un flujo continuo. Los médicos e investigadores toman mediciones en intervalos específicos —días, semanas o meses de diferencia—. Esta brecha entre la teoría suave del tiempo continuo y la naturaleza paso a paso de los datos del mundo real ha llevado a un equipo de investigadores de la Universidad de Alejandría a explorar un enfoque diferente: un modelo que se mueve en pasos discretos, reflejando la forma en que realmente observamos la enfermedad.

Los investigadores, Sarah Ramadan, Sanaa Salman y Ahmed El-Sayed, se propusieron construir una simulación computacional simplificada pero poderosa de la infección por VIH. Comenzaron con un modelo continuo conocido que describe la interacción entre las células inmunitarias no infectadas y el virus, luego lo tradujeron a un sistema paso a paso. En este nuevo marco, el estado de la infección se actualiza en saltos distintos en lugar de un deslizamiento suave. El equipo se centró en dos poblaciones principales: las células inmunitarias sanas que el virus ataca y las células infectadas que producen activamente nuevas partículas virales. Al ajustar la tasa a la que el virus infecta las células sanas, observaron cómo respondía el sistema a lo largo del tiempo. Su objetivo no era solo ver si el virus desaparecería o persistiría, sino descubrir los patrones ocultos que surgen cuando la tasa de infección cambia, buscando signos de un comportamiento complejo e impredecible que los modelos más simples podrían pasar por alto.

Lo que encontraron fue un mundo de sorprendente complejidad oculto dentro de las reglas simples de la infección. Cuando la tasa de infección es baja, el sistema se comporta de manera predecible, asentándose en un estado estacionario donde el virus es eliminado o mantenido en un nivel constante y manejable. Sin embargo, a medida que los investigadores aumentaban la tasa de infección en sus simulaciones, el sistema comenzaba a cambiar de maneras dramáticas. En lugar de estabilizarse, los números que representan al virus y a las células inmunitarias comenzaban a oscilar. No solo fluctuaban aleatoriamente; entraban en ciclos rítmicos, subiendo y bajando en un patrón regular. En algunos escenarios, estos ciclos duplicaban su longitud, creando un patrón donde el virus alcanza su pico dos veces más seguido antes de que el ciclo se repita. En otros casos, el comportamiento se volvía aún más intrincado, con el sistema trazando un bucle suave y cerrado que nunca repite exactamente el mismo camino, un estado conocido como movimiento cuasiperiódico.

Quizás lo más sorprendente fue que el equipo descubrió que, bajo ciertas condiciones, el sistema podía caer en el caos. En este estado, la carga viral y los recuentos de células inmunitarias fluctúan de una manera que parece aleatoria y es imposible de predecir a largo plazo, a pesar de que las reglas que gobiernan el sistema son fijas y deterministas. Esto significa que las oscilaciones erráticas en los niveles virales que se ven ocasionalmente en pacientes podrían no ser causadas por ruido externo o azar, sino que podrían ser una característica inherente de la propia interacción biológica. Los investigadores utilizaron una herramienta llamada exponente de Lyapunov para medir este comportamiento, un método que actúa como un termómetro para el caos. Cuando la lectura era negativa, el sistema era estable o rítmico; cuando se volvía positiva, el sistema había entrado en un régimen caótico. Sus simulaciones mostraron que la transición del orden al caos ocurría a través de una ruta específica, donde el sistema primero comenzaba a oscilar, luego las oscilaciones se duplicaban y, finalmente, el patrón se rompía en el caos.

El estudio también identificó umbrales específicos donde ocurren estos cambios dramáticos. Al variar el parámetro de la tasa de infección, el equipo localizó los momentos exactos en los que el sistema pierde su estabilidad. Encontraron que el modelo experimenta dos tipos distintos de transiciones. La primera es una bifurcación de salto (flip bifurcation), donde el sistema pasa de un estado estacionario a un ciclo que se repite cada dos pasos, y luego a ciclos de cuatro, ocho, etcétera, hasta volverse caótico. La segunda es una bifurcación de Neimark-Sacker, donde el sistema se mueve directamente desde un estado estacionario hacia una oscilación suave y en bucle que nunca se repite exactamente. Estos hallazgos sugieren que el comportamiento complejo y a veces errático del VIH en el cuerpo no es necesariamente un signo de un tratamiento fallido o de un entorno impredecible, sino una consecuencia natural de la retroalimentación no lineal entre el virus y el sistema inmunológico.

Las implicaciones de este trabajo se extienden más allá de las matemáticas. Los investigadores sugieren que los "picos virales" (viral blips) —aumentos repentinos y temporales en la carga viral que ocurren incluso en pacientes que de otro modo estarían estables con la medicación— podrían ser un resultado determinista de estas dinámicas internas en lugar de errores aleatorios. Si el sistema es naturalmente propenso a estas oscilaciones y estados caóticos, entonces el tiempo de las oleadas virales podría ser una propiedad intrínseca de la propia infección. Esta perspectiva desplaza el enfoque de intentar eliminar todas las fluctuaciones hacia la comprensión de los ritmos subyacentes de la enfermedad. El estudio no pretende haber resuelto el problema del VIH ni ofrecer una nueva cura, pero proporciona un marco teórico robusto para comprender por qué el virus se comporta de la manera en que lo hace. Al demostrar que los modelos de tiempo discreto pueden capturar comportamientos ricos y complejos como la periodicidad y el caos, el trabajo ofrece una nueva lente a través de la cual ver la persistencia del virus y las posibles estrategias para controlarlo. Los resultados, confirmados mediante extensas simulaciones computacionales, resaltan que la batalla entre el virus y el sistema inmunológico es una danza numérica dinámica y siempre cambiante que puede derivar en la imprevisibilidad con el más mínimo cambio en las condiciones.