INTRODUCCIÓN
El modelado permite la representación matemática de diferentes tipos de procesos para su análisis en diferentes escalas. En el campo de la medicina, los ingenieros biomédicos han utilizado esta herramienta para el estudio de los sistemas fisiológicos con el fin de lograr una mejor comprensión de estos. Por otro lado, el sistema cardiovascular ha sido objeto de estudio para el desarrollo de diferentes modelos matemáticos donde se han abstraído las principales características, dependiendo de las necesidades y los objetivos de cada estudio. En nuestro caso, se propone el uso de un modelo del sistema cardiovascular y su control nervioso sobre la actividad cardíaca para el análisis de la progresión del daño miocárdico en la ECCH.
La enfermedad de Chagas es una enfermedad tropical endémica producida por el Trypanosoma cruzi. Se ha propagado en varios países según casos notificados desde el sur de los Estados Unidos hasta el sur de Argentina, principalmente en las zonas rurales pobres de México, América Central y América del Sur. Por tal razón. Se ha convertido en un importante problema de salud en América Latina (1).
Al principio, la mayoría de los casos permanecen asintomáticos, pero después de 10 a 20 años, algunos pacientes comienzan a mostrar progresivamente signos electrocardiográficos de daño miocárdico según cada estadio de la enfermedad. Este daño se caracteriza por una considerable ampliación de los ventrículos que puede conducir a una muerte súbita. Sin embargo, un alto porcentaje de sujetos infectados no muestra ningún síntoma ni declaración clínica en el transcurso de su vida, por lo que la enfermedad permanece en una forma indeterminada. Estos sujetos actúan como reservorios del parásito y ayudan a mantener el ciclo de su vida (2).
Aunque esta enfermedad ha sido ampliamente estudiada, no hay una razón clara y definitiva que explique la causa del daño del miocardio. Existen dos hipótesis que lo justifican. La primera se conoce como teoría Neurogénica. Propuesta por Koberle, quien estableció que el daño es producido por la denervación parasimpática y la influencia simpática, que traen como consecuencia el adelgazamiento del ápex del ventrículo izquierdo y la progresiva dilatación de las cámaras cardíacas (3,4). Esta teoría está respaldada por algunas investigaciones realizadas en Brasil (5). Sin embargo, en estudios realizados en pacientes venezolanos con la enfermedad de Chagas, se detectaron alteraciones parasimpáticas justo después de haberse producido un importante avance en el deterioro de las cámaras cardíacas; insinuando así, la segunda hipótesis: neurogénica modificada y unificadora. Esta establece que la denervación parasimpática no representa una causa, sino una consecuencia de la remodelación ventricular (6,7).
Se realizó un estudio in silico utilizando los datos y resultados reportados por Carrasco et al., en una investigación relacionada con la detección precoz del daño miocárdico como consecuencia de la enfermedad de Chagas y los diferentes estadios evolutivos (I, II y III) (6). Específicamente, en este trabajo de investigación, se representaron en un modelo matemático las lesiones cardíacas relacionadas con la contractilidad y distensibilidad del ventrículo izquierdo y se utilizaron los hallazgos hemodinámicos informados en el estudio para la comparación.
El modelo matemático del sistema cardiovascular que se utilizó para este análisis se desarrolló bajo el enfoque compartimental. Este enfoque ha demostrado en diversas investigaciones ser excelente para modelar los procesos fisiológicos del sistema cardiovascular en diferentes condiciones (8-16). El modelo incluye el sistema cardiovascular (SC), corazón y vasos, y el reflejo del barorreceptor como mecanismo de control de la presión arterial (PA). Este reflejo es el principal mecanismo de control del SC. Una de sus principales tareas es mantener la PA en los niveles adecuados para conservar el equilibrio de las funciones cardíacas junto con la acción del sistema nervioso central (SNC) a través del sistema nervioso autónomo (SNA), sistema simpático y parasimpático. Básicamente, este reflejo se inicia por receptores de distensión denominados barorreceptores o presorreceptores, localizados en las paredes de las grandes arterías sistémicas.
La elevación de la PA distiende los barorreceptores y les hace transmitir señales al interior del sistema nervioso central, desde donde vuelven señales de retroacción a la circulación para reducir la presión hasta su valor normal, realizando igualmente el proceso inverso. Para el mecanismo de control el modelo utiliza las ecuaciones matemáticas de los controladores proporcional, integral y derivativo (PID), que funcionan bajo el principio de realimentación de lazo cerrado adaptándose al principal mecanismo homeostático fisiológico.
El estudio consistió en utilizar los parámetros y variables del modelo para simular el avance progresivo de las lesiones que aparecen en la enfermedad de Chagas, tanto en el corazón como en el sistema nervioso. Al realizar las simulaciones de los cambios reportados en cada una de las teorías sobre la patología cardíaca en la enfermedad, se intenta establecer, a través del modelo, que el daño miocárdico no es producto de la denervación parasimpática, respaldando así a la teoría Neurogénica modificada.
METODOLOGÍA
En esta sección se explican las etapas en el desarrollo del estudio. Se detalla en primer lugar el modelo matemático base, se describen los datos experimentales utilizados para las pruebas realizadas en el modelo. Posteriormente, los datos experimentales son simulados en el modelo ajustado con base en las dos teorías que intentan explicar el daño miocárdico y la función de las características de cada estadio de la ECCH.
Modelo matemático
Se utiliza como base para este estudio un modelo matemático del sistema cardiovascular desarrollado por Danielsen (12), modificado por Sánchez y Rojas quienes incluyeron el modelo del sistema de control autónomo como una respuesta del sistema cardiovascular al reflejo del barorreceptor (16). Es un modelo compartimental que incluye el corazón, el sistema circulatorio y el reflejo del barorreceptor como mecanismo de control de la presión arterial.
La Figura 1 muestra el esquema del sistema cardiovascular y del control nervioso utilizado en la investigación.
El corazón está representado por ecuaciones basadas en su comportamiento mecánico. Tanto el ventrículo izquierdo como la aurícula se definen por medio de una elastancia de tiempo variable de acuerdo con las ecuaciones ec.1, ec.2 y ec.3, respectivamente (12,13,16,17).
La elastancia del ventrículo izquierdo, E lv , se define por una suma ponderada de dos señales sinusoidales para el período de contracción ventricular y una constante para el período diastólico de la siguiente manera:
Donde q(t) representa la fase de eyección ventricular:
E min y E max son los valores mínimos diastólicos y máximos sistólicos de la función de elastancia ventricular, respectivamente. T t , es el período de tiempo total del ciclo cardíaco y, t s , es el período de tiempo de la sístole, a y b son constantes que definen la amplitud de la onda.
La elastancia, E la , también se define como variable en el tiempo, de acuerdo con la siguiente ecuación (17):
Donde m y p son constantes que definen la amplitud de onda; n y q constantes que definen el período de tiempo de la onda y, t*, es el tiempo aproximado para la ocurrencia de la contracción auricular.
El flujo ventricular izquierdo (Q lv ), el volumen (V lv ) y la presión (P lv ) se definen por medio de ecuaciones diferenciales simples. Las ecuaciones de la ec.4 a la ec.7 muestran las diferentes relaciones entre presión, volumen y flujo en la cámara cardíaca izquierda (las ecuaciones para el lado derecho son similares) (12,13).
V o,lv corresponden al volumen muerto a presión cero en el ventrículo izquierdo; Q la , es el flujo auricular izquierdo; Q lv , el flujo del ventrículo izquierdo; R la e I la , corresponden a los elementos resistivos e inductivos que simulan las propiedades inerciales y viscosas de la sangre; P la , presión auricular izquierda; E la , elastancia auricular izquierda; V la , volumen auricular izquierdo y, V 0 , volumen muerto en la aurícula izquierda.
La circulación sistémica se modela utilizando la relación de flujo de Poiseuille (18), que permite una acción directa del control nervioso sobre el radio de los vasos considerados.
Δp es la presión diferencial, ρ viscosidad sanguínea, L longitud y r el radio del vaso. Como parte de la validación del modelo matemático utilizado, se muestran las variables hemodinámicas generadas por él encargadas de definir el funcionamiento del sistema cardiovascular (Tabla 1). Se puede apreciar que las variables hemodinámicas calculadas a partir del modelo están dentro de los valores referenciales reportados en la literatura, logrando así, simular el comportamiento en un estado normal del sistema cardiovascular.
Variable | Unidades | Valor generado por el modelo | Valor referencial aproximado (30) |
---|---|---|---|
VFD | ml | 126,89 | 120 |
VFS | ml | 54,32 | 50 |
VL | ml | 72,56 | 70 |
FE | % | 57,16 | 60 |
PFS | mmHg | 106,63 | 110 |
FC | lat/min | 73,74 | 70 |
GC | ml/min | 5.336 | 5.000 |
PAM | mmHg | 90 | 90 |
VFD: volumen final diastólico. VFS: volumen final sistólico. VL: volumen latido. FE: fracción de eyección. PFS: presión final sistólica. FC: frecuencia cardíaca. GC: gasto cardíaco. PAM: presión arterial media. Fuente: (16)
El control parasimpático y simpático del sistema cardiovascular se simula en función de las características fisiológicas típicas. El sistema nervioso parasimpático (SNP) tiene sinapsis de neuronas preganglionares cerca o dentro del órgano efector, lo que da como resultado una respuesta limitada, rápida y localizada (19). En este caso se usa el control proporcional derivativo (PD) cuando se produce un cambio en la señal de error. La señal de error (e(t)) en este caso se define como la diferencia entre la variable de salida del modelo (presión sanguínea media [PSM]) y la señal de referencia (valor fisiológico normal de la PSM normal). El controlador tiene características de predicción, lo que aumenta la velocidad de respuesta, si la derivada de e(t) es cero, la acción de control está limitada en régimen permanente (20). Además, la acción del sistema simpático (SNS), el cual tiene largas fibras posganglionares que terminan en el cuerpo efector final inervando, glándulas, órganos y tejidos, también una extensa ramificación periférica, mantiene la respuesta para lograr un gran equilibrio, mostrando la naturaleza amplia y masiva de la estimulación simpática (19). Este comportamiento está asociado con el control proporcional integral (PI), donde a través de la acción integral, un pequeño error positivo producirá un aumento en la señal de control y un error negativo lo disminuirá, no importa cuán pequeño sea el error (20), es esta la acción que se mantiene, lo que refleja la acción masiva del SNS.
Estas relaciones entre el SNS y el SNP se producen para realizar la acción del control rápido de la presión arterial a través del control del ritmo cardíaco, cambios en la resistencia periférica (diámetro o radio de las arteriolas) y en la compliancia venosa (10-12,21). La estructura del sistema de control propuesto establece un conjunto de controladores con acciones de control PD o PI para cada efector, como se muestra en las ecuaciones ec.9 y ec.10, respectivamente.
La Figura 2 muestra cómo a partir de las variaciones de la frecuencia cardíaca y en la resistencia periférica ejecutadas por las acciones del control del SNS y el SNP, el modelo es capaz de simular el control sobre la presión arterial.
Para un mayor detalle del modelo de control de los barorreceptores por medio de las estructuras matemáticas PD y PI se recomienda leer (16).
Datos experimentales
Los datos utilizados en la investigación se estructuran en dos partes. La primera corresponde a los datos de los pacientes chagásicos en diferentes estadios de la enfermedad (Tabla 2). En la segunda, están los valores generados por el modelo cuando se simulan las condiciones correspondientes a las hipótesis que intentan explicar el daño miocárdico. Para estos últimos, inicialmente, se toma el modelo sin ningún tipo de variación, es decir, simulando un sujeto en condiciones normales y seguidamente se realizan variaciones de los parámetros correspondientes del modelo matemático para representar cada teoría.
Grupo | Variables hemodinámicas | |||||
---|---|---|---|---|---|---|
Pra mmHg | ∆Ppa mmHg | EDP mmHg | Pao mmHg | EDV ml | EF % | |
Control | 3 ± 1 | - | 7 ± 3 | 97 ± 10 | 90 ± 18 | 68 ± 10 |
Ia | 3 ± 1 | -1 ± 1 | 7 ± 3 | 97 ± 10 | 96 ± 31 | 67 ± 10 |
Ib | 4 ± 2 | 2 ± 4 | 8 ± 3 | 94 ± 13 | 136 ± 37 | 60 ±1 |
II | 6 ± 6 | 1 ± 3 | 9 ± 3 | 92 ± 14 | 153 ± 53 | 43 ±15 |
III | 7 ± 5 | 17 ± 4 | 17 ± 9 | 88 ± 19 | 255 ± 71 | 20 ±10 |
Pra: presión auricular derecha. ∆Ppa: valor diferencial de la presión arterial pulmonar (valor base: 13 ± 1). EDP: presión diastólica final. Pao: presión aórtica. EDV: volumen diastólico final. EF: fracción de eyección. Fuente: (16)
Los datos correspondientes a los pacientes chagásicos son parte de los resultados hemodinámicos previamente reportados por Carrasco et al., recopilados de pacientes con enfermedad de Chagas donde se identificaron signos cineangiográficos de daño miocárdico temprano y la evolución de la lesión característica de la ECCH. Los datos incluyen información sobre sujetos normales como grupo de Control y sobre pacientes chagásicos crónicos, divididos estos últimos, en tres grupos con diferentes estadios de evolución de la enfermedad. La Tabla 2 resume los valores hemodinámicos registrados que caracterizan a cada grupo de sujetos (6).
126 pacientes fueron incluidos en el estudio, dieron positivo en el test de fijación del complemento (FC) (reacción de Machado-Guerreiro) sin considerar otra cardiopatía. Otro grupo de 29 sujetos resultaron negativos y se establecieron como el grupo de Control. Este estuvo conformado por 17 sujetos masculinos y 12 femeninos con edades promedio entre 38 ± 13 años. Con respecto a los pacientes chagásicos, el subgrupo Ia correspondía a 21 sujetos de sexo masculino y 9 femeninos con edades entre 20 y 46 años. El subgrupo Ib con edad promedio entre 40 ± 14 años estaba compuesto por 17 sujetos de sexo masculino y 4 femeninos. En el grupo II eran en total 39 sujetos: 24 eran de sexo masculino y 15 femenino, con edades comprendidas entre los 30 y 68 años. Por último, el grupo III estaba compuesto por 36 sujetos con un promedio entre 54 ± 11 años, de los cuales 28 eran de sexo masculino y 8 femeninos (6).
El grupo de Control presentó un electrocardiograma (ECG) sin cambios en los segmentos ST y T, mostró valores en los límites normales en los exámenes hemodinámicos y angiográficos. Los siguientes grupos y pacientes se clasificaron según el estado de daño miocárdico. Grupo I fueron positivos en la prueba de reacción para la enfermedad de Chagas, se dividieron en dos subgrupos: los pacientes en el subgrupo Ia presentaron estudios hemodinámicos normales, radiografía (RX) de tórax y cineangiografía del ventrículo izquierdo. Se les consideró como pacientes con solo infección sin evidencia de daño al miocardio. El subgrupo Ib, eran pacientes con ECG y RX de tórax normales con índice diastólico final del ventrículo izquierdo aumentado y, el índice de distensibilidad del ventrículo izquierdo, disminuido en el estudio hemodinámico. El grupo II presentó un ECG anormal sin evidencia clínica de insuficiencia cardíaca congestiva, valores hemodinámicos alterados, dilatación del ventrículo izquierdo, distensibilidad y fracción de eyección disminuida. Finalmente, el grupo III presentó evidencia clínica de insuficiencia cardíaca congestiva, cardiomegalia y ECG anormal, insuficiencia biventricular, elevación significativa tanto de la presión del ventrículo derecho como de la resistencia pulmonar, además, de un aumento de la presión de fin de diástole del ventrículo izquierdo y una reducción en el índice de contractilidad del índice cardíaco y de la distensibilidad (6).
Estudio comparativo
La fase de análisis de las teorías que explican el daño miocárdico se realizó simulando ambas situaciones en el modelo y observando la tendencia de las variables hemodinámicas obtenidas como resultado de dichas simulaciones.
El daño progresivo al miocardio se simuló variando gradualmente los parámetros que representan el ventrículo izquierdo en la ecuación ec.1; la elastancia máxima para la pérdida de la contractilidad y el parámetro de elastancia mínima para la pérdida de distensibilidad correspondientes (E max y Emin, respectivamente). En este caso, no se realizaron variaciones en el mecanismo de control. Las diferencias de variación según el estado de avance de la enfermedad se pueden establecer en diferentes grupos: Control, Leve, Moderado y Alto. El grupo de Control corresponde al modelo sin ningún tipo de cambio, por lo tanto, simula las variables hemodinámicas reales para el modelo: un sujeto normal. Los grupos sucesivos muestran las variables hemodinámicas de la alteración progresiva de los parámetros seleccionados (hipótesis A: teoría Neurogénica modificada).
Como segunda fase y para representar lo que se denomina la hipótesis B (teoría Neurogénica) se simuló la denervación progresiva del control parasimpático, es decir, se alteraron o disminuyeron en un porcentaje de su valor en estado normal los parámetros del controlador PD (K P y K D ) en la ecuación (9), (0 % de variación equivale a un sujeto en estado normal). Se simuló, entonces, la disminución de la actividad parasimpática en 0 % (Control), 20 % (Leve), 40 % (Moderada) y 60 % (Alta). Para ambas hipótesis se registraron las variables hemodinámicas y los resultados se compararon según la tendencia de dichas variables obtenidas en las respectivas simulaciones.
RESULTADOS
La Tabla 3 resume los valores de las variables hemodinámicas del modelo resultado luego de simular el avance progresivo del daño miocárdico sin alteraciones en la inervación parasimpática (hipótesis A). El primer daño miocárdico reflejado en la enfermedad de Chagas es la pérdida de la distensibilidad, en función de esto, se reduce en un 30 % el parámetro correspondiente a la distensibilidad del ventrículo izquierdo en el modelo denominando de estos resultados como grupo Leve. Se midieron las variables observando un aumento considerable, tanto en la presión de las arterias pulmonar y aórtica como en el volumen y la presión diastólica final, además, de una disminución en la fracción de eyección. Este comportamiento se corresponde, principalmente, con los resultados informados para los sujetos de los grupos Ia y Ib analizados en el estudio de Carrasco et al., (6). Después de esto y como siguiente fase de la enfermedad, se disminuyó el parámetro de contractilidad representando el grupo con daño cardíaco moderado, para el cual se simuló una pérdida del 50 % de distensibilidad y 30 % de pérdida en la contractilidad. Las variables hemodinámicas obtenidas reflejan la misma tendencia de los reportados en el grupo II. Finalmente, se considera un cuarto grupo que corresponde a un estado de mayor evolución de la enfermedad (Alto), en este caso se representa mediante una pérdida del 70 % de distensibilidad y una pérdida del 60 % de contractilidad. Los resultados de este grupo muestran el mismo comportamiento de las variables obtenidas en el grupo III correspondientes al mayor daño del miocardio en el estudio considerado.
Estado de evolución | Variables hemodinámicas | |||||
---|---|---|---|---|---|---|
Pra mmHg | Ppa MmHg | EDP mmHg | Pao mmHg | EDV ml | EF % | |
Control | 2,70 | - | 5,70 | 94,60 | 126,57 | 57 |
Leve | 3,10 | 0,10 | 7,91 | 113,17 | 134,21 | 55 |
Moderado | 3,43 | -0,17 | 10,89 | 11,53 | 158,10 | 47 |
Alto | 3,66 | 7,76 | 12,50 | 98,67 | 169,54 | 38 |
Pra: presión auricular derecha. Ppa: valor diferencial de la presión arterial pulmonar (valor base: 13 ± 1). EDP: presión diastólica final. Pao: presión aórtica. EDV: volumen diastólico final. EF: fracción de eyección. Fuente: creación propia
Seguidamente, se realiza el análisis relacionando la inervación parasimpática con el registro de las variables hemodinámicas estudiadas, simulando así, la hipótesis B. En la Tabla 4 se muestran estos resultados (disminución de la actividad parasimpática [PS dec.]).
PS dec. | Variables hemodinámicas | |||||
---|---|---|---|---|---|---|
Pra mmHg | Ppa mmHg | EDP mmHg | Pao mmHg | EDV ml | EF % | |
Control | 2,70 | - | 5,70 | 94,60 | 126,57 | 0,57 |
Leve | 3,07 | -2,57 | 5,50 | 104,21 | 122,29 | 0,56 |
Moderado | 2,73 | -2,77 | 5,38 | 104,47 | 119,70 | 0,53 |
Alto | 2,52 | -2,89 | 5,18 | 104,37 | 115,73 | 0,53 |
PS dec: disminución de la actividad parasimpática. Pra: presión auricular derecha. Ppa: valor diferencial de la presión arterial pulmonar (valor base: 13 ± 1). EDP: presión diastólica final. Pao: presión aórtica. EDV: volumen diastólico final. EF: fracción de eyección. Fuente: creación propia
Estos resultados indican que la denervación parasimpática en el modelo no ocasiona, por si sola cambios inmediatos en las variables hemodinámicas, mientras que los cambios en la distensibilidad y contractilidad sí pueden ser responsables de dichas variaciones.
DISCUSIÓN
De acuerdo con los resultados obtenidos en función a los cambios en los parámetros del modelo relacionados con la remodelación ventricular y la inervación nerviosa, se observó que solo la denervación parasimpática no causa cambios inmediatos en las variables hemodinámicas, mientras que los cambios en la distensibilidad y la contractilidad pueden ser responsables de los mayores cambios hemodinámicos. Aunque los resultados obtenidos no son concluyentes se puede suponer que el avance en el daño miocárdico en la fase aguda de la enfermedad puede deberse a la remodelación ventricular y no solo a la denervación parasimpática.
Es importante tener en cuenta que todas las fallas cardíacas incluyen alteraciones del sistema nervioso autónomo que causan un desequilibrio con una mayor actividad simpática y una actividad parasimpática disminuida (22). El estudio se centró en la discusión sobre la causa y el avance de la disfunción cardíaca en la enfermedad de Chagas utilizando un modelo matemático que permite, a través de la variación de sus parámetros, simular la relación entre el sistema cardiovascular y el sistema nervioso, que están comprometidos en la patología chagásica según las teorías existentes.
Se han utilizado modelos matemáticos en investigaciones sobre la enfermedad de Chagas con diversos objetivos, tales como en el avance de la enfermedad y el sistema inmune (23); otros referidos a la dinámica o forma de transmisión (24-25). Sin embargo, su utilización en el estudio de la progresión del daño miocárdico, según la literatura revisada, no había sido aplicado, lo que permite que a través de esta investigación, se pueda proporcionar información sobre la controversia entre las teorías que intentan explicar el daño, de forma tal, para encontrar un diagnóstico preciso que conlleve a un mejor tratamiento de la enfermedad y que mejore la calidad de vida de los pacientes, demostrando que el modelado matemático es una herramienta alternativa para apoyar la investigación de la evolución de la patología cardíaca en la enfermedad de Chagas.
Si se tiene en cuenta que la patología chagásica es una de las enfermedades parasitarias más importantes del mundo, todo aporte que contribuya a su diagnóstico permitirá un tratamiento efectivo que mejore la calidad y alargue la vida del paciente, además, estas investigaciones también ayudan a que en los centros de salud los diagnósticos estén mejor orientados al tener mayor información sobre las patologías que pudieran presentarse, mejorando así, el servicio y la atención.