INTRODUCCIÓN
El VIH o Virus de Inmunodeficiencia Humana es un lentivirus que ataca al sistema inmunológico de las personas, debilitándolo y haciéndoles vulnerables ante una serie de infecciones que aparecen debido a que las defensas inmunitarias son insuficientes. El sistema inmunológico se considera deficiente cuando la cantidad de linfocitos T CD4+ está por debajo de 500 células/mm 3 y se vuelve incapaz de cumplir su función de lucha contra infecciones y enfermedades oportunistas [1].
Un individuo infectado de VIH es conocido como una persona seropositiva; gracias a los adelantos recientes en el acceso al tratamiento con antirretrovirales, las personas seropositivas pueden vivir más tiempo y en mejor estado de salud. Además, la Organización Mundial de la Salud (OMS) ha confirmado que el tratamiento con antirretrovirales ayuda a disminuir la probabilidad de transmisión del VIH. Si una persona con VIH no recibe tratamiento a tiempo, su condición pasará a ser denominada SIDA (Síndrome de Inmunodeficiencia Humana) que es el estado de la infección por el VIH y se identifica principalmente por bajos niveles de defensas [2].
En el mundo existen alrededor de 34.2 millones de personas infectadas por VIH, de las cuales la mayoría vive en zonas de bajos recursos, desinformadas, con baja o nula sanidad y carecen de atención médica de calidad. De acuerdo con las estadísticas del Centro Nacional para la Prevención y el Control del VIH y SIDA (CENSIDA), en México existen 180,000 personas con VIH, de las cuales el 50% no tenían noción de ser portadoras del virus sino hasta llegar a una etapa avanzada de la enfermedad. Dentro del país, Baja California ocupa el cuarto lugar nacional con mayor número de población infectada del VIH, que de acuerdo con el CENSIDA reporta una tasa de 272.2 personas infectadas por cada 100,000 habitantes [3].
Cabe destacar que se podría evitar un número significativo de muertes si se recurriera al uso constante de antirretrovirales; en el país tan solo 57,073 del total de personas infectadas toma regularmente el medicamento, que según la Secretaria de Salud tiene un costo anual de $44,997 pesos por persona y se proporciona por algunas instituciones de salud pública de manera gratuita. A finales del año 2014, 14.9 millones de personas recibieron terapia antirretrovírica en todo el mundo, lo que representa el 40% de los 36.9 millones de personas que viven con VIH [3].
Con el propósito de comprender la compleja dinámica del VIH, los biólogos y matemáticos han recurrido al modelizado matemático como una herramienta para obtener información acerca de la relación entre la respuesta del sistema inmunológico, el virus libre en el organismo y la aplicación de tratamientos. Algunos modelos se diseñan con datos experimentales y resultados clínicos, en la literatura se encuentran diversas publicaciones, por ejemplo [4,5,6,7].
La Localización de Conjuntos Compactos Invariantes (LCCI) es un método que permite analizar la dinámica global de modelos matemáticos de ecuaciones diferenciales ordinarias de primer orden. El método fue propuesto por Krishchenko [8] y optimizado por Krishchenko y Starkov [9]. Recientemente se ha utilizado para el análisis de sistemas biológicos [10,11,12] y en esta investigación se utiliza como base para modelizar un sistema biológico que describe la evolución del VIH bajo tratamiento antirretroviral.
El trabajo se organiza de la siguiente manera, en la siguiente sección se presenta la teoría correspondiente al método de LCCI. Después, se muestra el modelo matemático desarrollado para describir la dinámica del VIH bajo tratamiento antirretroviral. En la sección que le sigue se aplica el método de LCCI para determinar los límites superiores de cada variable del modelo matemático. En la sección subsecuente se ilustran los resultados al graficar las soluciones del modelo matemático propuesto y los límites de localización. Finalmente, se muestran las conclusiones del trabajo desarrollado.
METODOLOGÍA
Método de localización de conjuntos compactos invariantes
El método de LCCI se utiliza para determinar un dominio en Rn en el cual se localizan todos los conjuntos compactos invariantes que se presentan bajo ciertas condiciones en un sistema específico, estos conjuntos pueden ser: órbitas periódicas, homoclínicas y heteroclínicas, ciclos límite, puntos de equilibrio y atractores caóticos. La importancia del método radica en que el análisis es útil para conocer la dinámica del sistema en el largo plazo. Su característica principal consiste en que es un método estrictamente analítico, lo que implica la solución del problema sin la necesidad de realizar la integración numérica del sistema de ecuaciones diferenciales.
A continuación, se describirán los teoremas, notaciones y definiciones básicas utilizadas. Considere un sistema no lineal de la forma:
donde ƒ es una función vectorial continua para un C∞ y x ∈ Rn es el vector de estados. Sea h(x): Rn → R, la cual es llamada función localizadora y no es la primera integral de (1), entonces, por h|B se denota la restricción de h a un conjunto B ⊂ Rn. Por S(h) se denota el conjunto {x ∈ Rn | Lƒh(x) = 0}, donde Lƒh es la derivada Lie de (1) y está dada por: Lƒh = (∂h/∂x)ƒ(x). Además, se define
A continuación, se definirán el Teorema General de LCCI y el Teorema Iterativo.
Teorema 1. Vea [8,9]. Cada conjunto compacto invariante Γ de (1) está contenido en el conjunto de localización
Un refinamiento del conjunto de localización K(h) puede realizarse con el uso del teorema iterativo que dice:
Teorema 2. Vea [8,9]. Sea hm(x),m = 0,1,2,… una secuencia de funciones de clase infinitamente diferenciable. Los conjuntos
Con
contienen cualquier conjunto compacto invariante del sistema (1) y
Este método ya se ha aplicado en con otras teorías de estabilidad para analizar modelos matemáticos de sistemas biológicos que describen la evolución del cáncer, la respuesta del sistema inmunológico y el efecto de tratamientos como la quimioterapia e inmunoterapia [10,11,12].
Modelo matemático de VIH con acción de antirretrovirales
En esta sección se describe el modelo matemático del VIH con acción de antirretrovirales. Las Ecuaciones (2)-(4) describen el proceso de replicación del VIH en células sanas de dos formas. La primera de manera directa por parte del VIH libre en el organismo y la segunda de manera indirecta, mediante la interacción entre células infectadas y células sanas [6], dicha dinámica se describe mediante las siguientes tres ecuaciones diferenciales ordinarias de primer orden:
Al considerar lo anterior, se propone una modificación al modelo que consiste en la adición de una ecuación que representa el efecto del tratamiento antirretroviral en el proceso de replicación del VIH [3]. El nuevo modelo se define mediante las siguientes cuatro ecuaciones:
La Ecuación (5) modeliza la concentración de células sanas en el sistema; el parámetro s indica la producción de manera natural de linfocitos T CD4+ por el sistema inmunológico. El crecimiento de la población de estos linfocitos se representa mediante la ley de crecimiento logístico, con una tasa de crecimiento r 2 y una carga máxima de m, estos tienen una muerte natural dada por μT y su población disminuye por la interacción con el VIH y las células infectadas a una tasa dada por k 2 y k 3 , respectivamente.
La Ecuación (6) describe la evolución de las células infectadas por el VIH y su interacción con células sanas a una tasa dada por k 2 y k 3 respectivamente. La replicación del virus por las células infectadas produce una lisis celular con una tasa dada por μI.
La Ecuación (7) representa la concentración de VIH que circula libre en el torrente sanguíneo donde q representa el número total de partículas de virus producidas por una célula infectada durante su tiempo de vida y el parámetro δX define la capacidad de los antirretrovirales en conjunto con el sistema inmune para eliminar las partículas de VIH que se encuentran libres en el sistema, el valor de este parámetro se determina con un aumento porcentual debido a la acción del tratamiento en la eliminación de las partículas del virus.
La Ecuación (8) describe la dinámica del tratamiento antirretroviral [3]; el parámetro γ determina la duración del efecto antirretroviral en el sistema, asumiendo que la ingesta del mismo se realiza diariamente [13]. La concentración de medicamento ingerida por el paciente se define mediante el parámetro A X cuyo valor depende de datos personales como peso, altura, tiempo de la infección, edad y nivel de células no infectadas, entre los más importantes.
En la Tabla 1 se presenta la descripción, valores y unidades de cada parámetro del sistema (5)-(8). Estos valores son presentados por Lou et. al. [6].
Adicionalmente, cabe destacar que la dinámica del modelo matemático de VIH con acción de antirretrovirales (5)-(8) se localiza en el ortante no negativo [14] definido como sigue:
Localización de conjuntos compactos invariantes para el modelo matemático (5)-(8)
En esta sección se muestra el análisis matemático necesario para obtener un dominio compacto en el espacio R4+,0. El dominio de localización está definido por los límites superiores e inferiores de cada variable de estado del modelo de VIH en conjunto con una terapia antirretroviral, dichos límites se encuentran descritos por medio de desigualdades en función de los parámetros del sistema. Biológicamente, estos límites representan los valores máximos que pueden alcanzar las poblaciones de células sanas, de células infectadas, el total de partículas libres de VIH y la concentración del antirretroviral en el torrente sanguíneo.
A continuación, se presentan las operaciones matemáticas necesarias para obtener los límites inferior y superior de la concentración del tratamiento antirretroviral presente en el torrente sanguíneo. Para lo cual se resuelve la ecuación diferencial
mediante separación de variables como se muestra a continuación
obteniendo el siguiente resultado para la familia uniparamétrica de soluciones
al considerar la condición inicial
Entonces, al determinar el siguiente límite
se obtiene la cota superior
El límite inferior se considera como el valor de la condición inicial, es decir,
el tratamiento antirretroviral. De los resultados anteriores se puede inferir
Para determinar el límite superior de la población de células sanas se propone la función localizadora lineal
a la cual se calcula su derivada de Lie como se muestra a continuación
posteriormente, se obtiene el conjunto
y entonces se determina el subconjunto
del cual se obtiene el valor máximo de la población de células sanas de la siguiente forma
Ahora, se presenta el análisis matemático para obtener el límite superior de la población de células infectadas que se encuentran en el organismo. En este caso se propone una segunda función localizadora dada por
cuya derivada de Lie se muestra a continuación
posteriormente, se obtiene el conjunto
y entonces se determina el subconjunto
con lo cual se determina el conjunto
en donde se despeja I para establecer Isup como se muestra
Finalmente se propone una tercera función localizadora para encontrar el límite superior de las partículas de VIH que se encuentran libres en el organismo, la función se define como
cuya derivada está dada por
La cual se iguala a cero para obtener
y se obtiene el valor máximo de la función al aplicar el Teorema Iterativo como se indica a continuación
por lo tanto, se concluye el siguiente conjunto
Teorema 3. Todos los conjuntos compactos invariantes del sistema de VIH con acción de antirretrovirales (5)-(8) se encuentran localizados dentro del dominio definido por
RESULTADOS Y DISCUSIÓN
Simulaciones numéricas
En esta sección se presentan simulaciones numéricas con el propósito de ilustrar los resultados en el plano temporal, para esto se grafican las soluciones del sistema y los límites superiores obtenidos, estos permiten observar el valor máximo que pueden alcanzar las poblaciones de células y la concentración del VIH en el torrente sanguíneo.
En la Figura 1 se ilustran las soluciones del modelo matemático (2)-(4), en las cuales se observa que existe una fluctuación considerable tanto de células sanas como de células infectadas a lo largo del tiempo, lo cual se debe al mecanismo de acción del virus debido a que al infectar al cuerpo humano comienza a presentarse como enfermedades oportunistas que usualmente se confunden con malestares pasajeros como lo son un cuadro de gripe o una infección estomacal, es por eso que durante las primeras etapas de la enfermedad es difícil detectar si el virus está presente en el organismo, además de que el recuento de linfocitos se mantiene en valores normales (más de 500 células/mm 3 ) por lo que las células sanas siguen presentes en el sistema; con el paso del tiempo este nivel comienza a decaer considerablemente a medida que el VIH libre continua infectando más células sanas. Por otra parte, el VIH libre actúa con una dinámica muy similar a las células infectadas, solo que con una concentración más alta. El virus continuará viviendo como un agente libre o en la mitocondria de las células replicándose hasta destruir por completo el sistema inmunológico del individuo.
La Figura 2 ilustra la dinámica del proceso de replicación del VIH cuando el paciente se encuentra bajo un tratamiento antirretroviral (5)-(8). A diferencia del modelo (2)-(4) en el sistema biológico propuesto las células sanas y las infectadas no presentan niveles tan variados, esto es debido a la aplicación del tratamiento antirretroviral. Las células sanas mantienen un nivel en el que el individuo se considera “saludable” hasta que el VIH comienza a atacar a los linfocitos y los niveles decaen, pero al estar bajo el efecto de la terapia antirretroviral las células infectadas solo tienen un periodo muy pequeño de alza, justamente cuando entran al cuerpo y conforme surge el efecto de los antirretrovirales baja su nivel a un estado “semi-saludable” para el portador del virus.
En la Figura 2 también se ilustra el valor máximo de cada población de células y del tratamiento, se observa que ninguna de las gráficas supera este punto. Las células sanas alcanzan un valor muy cercano al máximo en un punto de su solución. Sin embargo, a medida que el tiempo transcurre convergen a un punto de equilibrio localizado por debajo del límite superior.
El tratamiento evita que el virus continúe replicándose más no lo elimina por completo del organismo, es visible el hecho de que tanto el VIH libre en el sistema como las células infectadas siguen presentes en el individuo, sin embargo, las tres poblaciones de células alcanzan un nivel en el cual se considera que el paciente está clínicamente sano, es decir, con el virus bajo control. La concentración de la terapia antirretroviral se mantiene en un nivel constante durante todo el tiempo que sea ingerido por el paciente, asumiendo que estará bajo este tratamiento para toda su vida.
Se observa que la terapia antirretroviral cuenta con una concentración máxima y una mínima que estará presente en el torrente sanguíneo; el valor máximo representa la dosis del medicamento antirretroviral consumida por cada individuo que será absorbida paulatinamente en cierto lapso de tiempo, y el valor mínimo es la concentración que debe consumir del medicamento.
CONCLUSIONES
Los resultados obtenidos en este trabajo demuestran que mediante el uso del método de LCCI y la realización de simulaciones numéricas se logró analizar la dinámica del sistema propuesto de VIH cuando se aplica una terapia antirretroviral.
Con el método de LCCI se estimaron los límites máximos de las variables del sistema biológico, los cuales permiten encontrar un dominio acotado dentro del cual se encuentran todos sus conjuntos compactos invariantes, en particular puntos de equilibrio.
Los límites máximos de cada una de las ecuaciones del sistema tienen un significado biológico los cuales se escriben mediante desigualdades en función de los parámetros del sistema. En primera instancia la variable T(t) que representa la población de células sanas tiene un valor máximo el cual representa el valor de linfocitos T CD4+ que tiene la persona en el periodo de infección cuando el sistema inmunológico aún no ha sufrido los efectos de la infección del VIH, pero una vez que el virus comienza a infectar una mayor cantidad de células los niveles decaen. Con ayuda de la terapia antirretroviral se mantienen niveles saludables que no ponen en riesgo la vida del paciente. Cabe destacar que la terapia no elimina por completo al virus, sino que lo mantiene bajo control, esto implica que nos es posible alcanzar el valor máximo de células sanas.
En el caso en que se utiliza la terapia antirretroviral se observa que la variable I(t), que define la concentración de células infectadas en el torrente sanguíneo, converge a un valor por debajo del 50% del caso cuando no se utiliza el medicamento. La población de células infectadas tiene la capacidad de crecer hasta una máxima expresión que por lo regular será igual o menor al nivel máximo de células sanas. Esto se debe a que un sistema enfermo debe contar con un respaldo de células sanas que ayuden a la conservación del mismo, de otro modo el sistema colapsaría produciendo la muerte del paciente.
Las partículas de VIH son representadas por la variable V(t) y tiene un modo de acción muy similar a las células infectadas, estas partículas actúan de la misma forma estando libres en el torrente o en la mitocondria de la célula, solo que en este caso se manejan valores mucho más altos considerando que cada célula infectada produce cientos de partículas de VIH.
La concentración en el torrente sanguíneo del tratamiento antirretroviral X(t) tiene un valor máximo y un mínimo establecidos por la dosis del medicamento que consume el paciente; estos valores fueron obtenidos de la literatura considerando aquellos medicamentos distribuidos por el Sector Salud y son ingeridos diariamente. Como valor máximo se establece la concentración que existe en el cuerpo una vez consumido el medicamento y ha transcurrido un periodo suficiente de tiempo. Como valor mínimo se tiene la dosis de medicamento que es administrada diariamente para conservar su efectividad.
La contribución principal de esta investigación es el establecimiento de una base teórica para analizar la dinámica de un modelo matemático invariante en el tiempo de un sistema biológico compuesto por ecuaciones diferenciales ordinarias de primer orden, el cual describe la acción del VIH cuando un paciente está bajo terapia antirretroviral.