INTRODUCCIÓN
La hipertensión arterial (HTA) es una de las enfermedades que más aqueja a la población mundial. Según la Organización Mundial de la Salud, uno de cada cinco adultos la padece, y produce complicaciones derivadas que son la causa de 9,4 millones de defunciones cada año1,2. Atendiendo a esta preocupación, es apremiante proyectar acciones inmediatas en el control de la enfermedad, por lo que todos los esfuerzos destinados a incrementar el conocimiento etiopatogénico y de los mecanismos fisiopatológicos que subyacen en la HTA, constituyen pilares sobre los cuales se asientan las acciones de prevención, diagnóstico y tratamiento de esta enfermedad.
Etapas tempranas de la enfermedad hipertensiva (y algunas instancias de los estados pre-hipertensivos, particularmente en individuos con historia familiar de HTA) se caracterizan por un «estado circulatorio hiperquinético», el cual está mediado por un incremento de la actividad adrenérgica y una reducción de la función parasimpática3. Para la evaluación del desbalance autonómico presente en la HTA se han empleado, por muchos años, diferentes métodos4,5, pero no es hasta las últimas cuatro décadas, con el desarrollo de los métodos computadorizados aplicados a los estudios de la variabilidad de la frecuencia cardíaca (VFC), que se puede observar un avance acelerado en esta temática.
La VFC se define como la variación que ocurre en el intervalo de tiempo, entre latidos cardíacos consecutivos. Actualmente, los métodos más utilizados para evaluarla son los que se basan en los dominios del tiempo, de la frecuencia, de tiempo-frecuencia y los métodos no-lineales6-8. Los primeros tratan a la secuencia de intervalos RR como un grupo de intervalos no ordenados o como intervalos pareados, y emplean diferentes técnicas para expresar la varianza de los datos6-11. El análisis espectral consiste en descomponer el tacograma (registro de frecuencia cardíaca en el tiempo), el cual se asemeja a una onda compleja, de manera que se obtienen los componentes espectrales, y se encuentran componentes de alta frecuencia (HF), relacionados con el tono parasimpático, uno de baja frecuencia (LF), relacionado con la modulación de los sistemas nerviosos simpático y parasimpático; y, aunque se encuentra bajo polémica, se ha propuesto la razón entre estos componentes (LF/HF) como un indicador del balance simpático-vagal8.
Los métodos de tiempo-frecuencia son una extensión de los métodos frecuenciales, que posibilitan analizar los cambios en tiempo presente en el análisis espectral de la VFC. La técnica de la transformada continua de wavelet, función tipo Morlet (CWT-Morlet), es uno de los métodos de tiempo-frecuencia más utilizados en el análisis de la VFC por su adecuada resolución en la descripción de los cambios en el espectro de frecuencias que operan a lo largo del tiempo12. De manera general los indicadores de la VFC se correlacionan con adaptaciones fisiológicas a cambios en el medio interno, externo y a la presencia de enfermedades.
La hiperreactividad cardiovascular y simpática esclarece parte de la etiopatogenia de la HTA y de otras enfermedades cardiovasculares13. Las técnicas para inducir reactividad cardiovascular son varias, las cargas físicas han demostrado ser mucho más sensibles y específicas que todas las restantes14,15. Entre ellas, una de las más difundidas ha sido la prueba de esfuerzo isométrico, que ha mostrado un alto valor predictivo y sensibilidad en el diagnóstico de la HTA16-18. En nuestro medio, se usa una variante de la prueba de esfuerzo isométrico, la prueba del peso sostenido (PPS), que posee gran valor práctico para realizar pesquisas masivas de HTA19,20. Esta consiste en mantener un peso de 500 gramos en la mano izquierda con el brazo extendido en ángulo recto al cuerpo durante 2 minutos. La presión arterial (PA) se toma en el brazo contrario antes del ejercicio y en los últimos 10 segundos de este21. Con la aplicación de la PPS se ha demostrado que los individuos hiperreactivos cardiovasculares tienen un mayor riesgo de padecer HTA que los individuos normorreactivos, de ahí que constituya la hiperreactividad cardiovascular un factor de predicción muy importante para la HTA22.
A pesar de la amplia aplicación de esta prueba en la práctica médica en el contexto nacional e internacional, todos los estudios realizados con la PPS se basan en la respuesta hemodinámica a partir de la PA, por lo que constituye una necesidad conocer cómo es el comportamiento de la dinámica en el tiempo de la regulación autonómica cardíaca durante la PPS en sujetos sanos normotensos y normorreactivos, a partir del análisis de la VFC. Por estas razones, el objetivo de esta investigación ha sido determinar la respuesta autonómica cardiovascular durante la prueba isométrica cubana del peso sostenido.
MÉTODO
Se realizó un estudio no observacional cuasi-experimental (crossover) con 16 sujetos (9 mujeres) sanos, normotensos y normopesos (edad 25,8±7,5 años), que sirvieron como su propio control, 5 minutos en el estado de reposo y otros 5 en el estado de prueba del peso sostenido. La investigación se realizó en el Laboratorio de Ciencias Básicas de la Universidad de Ciencias Médicas de Santiago de Cuba y fue aprobada por el comité de ética de dicho centro.
Criterios de salida
Los que padecieran de afecciones del sistema nervioso (accidente vascular encefálico, afecciones neuromusculares) o generalizadas de la piel, o tuvieran implantados dispositivos electrónicos (marcapasos o desfibriladores automáticos), que presentaran alguno de sus miembros superiores o inferiores amputados o que padecieran de HTA. Además, mujeres en estado de gestación, o con la menstruación. Se excluyeron también aquellos individuos que presentaron algún tipo de arritmia en la inspección visual del registro electrocardiográfico en el período de 10 a 15 minutos de aclimatación al local, así como los sujetos que no desearan ser parte de la investigación y no dieran su consentimiento informado.
Registros fisiológicos
Al comienzo de la sesión de los registros del electrocardiograma, en la mañana (08:30-12:00 am), los sujetos se sentaron en un sillón confortable, situado en una habitación con temperatura controlada entre 24 y 27 grados Celsius, y luz tenue. En estas condiciones se les permitió reposar durante 10-15 minutos hasta lograr una mejor adaptación a las condiciones del local.
Los electrodos de superficie del electrocardiograma se colocaron después de limpiar la piel con alcohol y se ubicaron en la derivación DII extendido a miembros, y se realizó un registro durante 5 minutos (estado de reposo), luego de los cuales se tomó la PA con un esfigmomanómetro y un estetoscopio, certificados. El estado PPS de 5 minutos comprendió los 2 minutos en los que se sostuvo el peso de 500 gramos y los 3 de recuperación.
La señal de electrocardiografía del equipo Powerlab®, con un filtro pasabanda de 0,5-30 Hz, se digitalizó a una frecuencia de muestreo de 1000 muestras/segundo (1kHz), en el paquete de software LabChart® del 2012, ambos de producción australiana por la compañía AD Instruments.
Procesamiento de la señal electrocardiográfica, discriminación de las ondas R y cálculo de los intervalos RR
El procesamiento ulterior de los registros digitalizados incluyó su inspección visual. La discriminación de los picos R de la señal digitalizada y el cálculo de los intervalos RR se realizaron empleando el método de Sabarimalai-Manikandan23. El conjunto de intervalos RR obtenido fue almacenado y constituye la serie de datos a partir de la cual se realizó todo el análisis posterior de la VFC.
Pre-procesamiento de los intervalos RR
El pre-procesamiento de las series de intervalos RR se efectuó con el programa HRVAS (heart rate variability analysis software), Copyright 2015 de John T. Ramshur (https://sourceforge.net/projects/hrvas)24, que permitió realizar el pre-procesamiento de la serie de intervalos RR. Se empleó un filtro de porcentaje con valor de 20% del intervalo previo para detectar los latidos ectópicos, los que se reemplazaron a partir de la interpolación cúbica politómica. El utilitario Wavelet Packet Detrending fue usado para eliminar las tendencias de las muy bajas frecuencias sobre la línea de base.
Análisis de la variabilidad de la frecuencia cardíaca
En el análisis tradicional de la VFC, según lo recomendado por el Consenso Internacional de Expertos sobre la VFC de 19969, se utilizó el programa VFC32 Versión 2.0.3.2 para Windows, de la Facultad de Biología de la Universidad de la Habana25, con el empleo de los siguientes indicadores:
SDNN: Desviación estándar de todos los intervalos RR. Se conoce como Variabilidad Total8,9.
RMSSD : Es la raíz cuadrada del valor medio de la suma de las diferencias al cuadrado de todos los intervalos RR sucesivos. Este parámetro informa de las variaciones a corto plazo de los intervalos RR y se utiliza para observar la influencia del sistema nervioso parasimpático sobre el sistema cardiovascular. Se asocia directamente a la variabilidad a corto plazo8,9.
pNN50: Porcentaje de los intervalos RR consecutivos que discrepan en más de 50 ms entre sí. Un valor elevado de pNN50 nos proporciona valiosa información acerca de variaciones altas espontáneas de la frecuencia cardíaca8,9.
Bajas frecuencias (LF) (n.u.): energía normalizada en el espectro de 0,04 a 0,15 Hz en que se descompone la serie temporal de intervalos RR consecutivos:
Altas frecuencias (HF) (n.u.): energía normalizada en el espectro de 0,15 a 0,4 Hz en que se descompone la serie temporal de intervalos RR consecutivos:
Relación LF/HF: Cociente LF/HF. Aunque bajo polémica, se ha propuesto como un indicador del balance simpático-vagal8,9.
Entropía de Shannon: Se aplica la conocida fórmula de Shannon en el campo de la teoría de la información (Shannon, 1948), para calcular la cantidad de información contenida en la secuencia de intervalos RR:
Cálculo de los indicadores de la VFC en tiempo-frecuencia
Los intervalos RR fueron re-muestreados con una interpolación de 2 Hz (0,5 segundos) para el análisis en tiempo-frecuencia. Se utilizó la CWT-Morlet para el análisis en tiempo-frecuencia de la VFC. Esta técnica utiliza ventanas cortas para las HF y largas para las LF, y puede ser aplicada satisfactoriamente para el procesamiento de señales no-estacionarias, al indicar cuáles frecuencias están presentes en un instante de tiempo, y mostrar buena resolución temporal a las HF y buena resolución espectral a las LF. Teóricamente la función CWT-Morlet está calculada para traslaciones infinitesimalmente pequeñas y para factores de escala. Para una señal x(t) y la función wavelet Ψ ab (t), el coeficiente de transformada continua de wavelet está dado por: , donde Ψ*(t) es la conjugada compleja de la onda madre, α es el parámetro de dilatación y τ es el parámetro de localización. Se utilizó la función CWT-Morlet como una gaussiana, la cual está balanceada en tiempo y frecuencia definida como: , donde ω es una frecuencia adimensional, la cual define el número de ciclos de la función CWT-Morlet, con ω=6 se proveyó una buena calidad en la resolución temporal y frecuencial. La función bivariante W(τ,α) muestra la similaridad de x(t) a una onda escalada por α a un tiempo τ dado.
Para obtener los valores en tiempo-frecuencia se usaron métodos de poder instantáneos, donde el módulo cuadrado del coeficiente de wavelet fue integrado en toda la banda de frecuencia analizada [f1 f2]. El poder instantáneo de una banda de frecuencia [f1 f2] está dado por:
Se utilizaron las bandas de frecuencias tradicionales de la VFC según lo recomendado por el Consenso Internacional de Expertos sobre la VFC de 19969, los cuales fueron para las bajas frecuencias de [LF: 0.04-0.15 Hz]; altas frecuencias de [HF: 0.15-0.4 Hz]. La CWT-Morlet fue exportada en un “archivo.txt”, una matriz de resultados por cada sujeto para registros de 5 minutos (300 s), a partir de la interpolación de 2 Hz (0,5 segundos); quedando 600 valores en las bandas de frecuencias de la VFC (LF, HF). Para su mejor tratamiento estadístico las bandas de frecuencias fueron normalizadas (n.u.), según lo recomendado9.
Procesamiento estadístico de los datos
Con la utilización del Sistema SPSS 22.0, se exponen los valores medios (χ̅) y la desviación estándar de las variables a las cuales se les realizó un análisis estadístico no paramétrico con la prueba de los rangos con signo de Wilcoxon para las muestras relacionadas con un nivel de significación de p<0,05.
Para la elaboración de los gráficos dinámicos en el tiempo (gráficos 1 y 2), se utilizó el programa Matlab R2013b a partir del promedio de las matrices de los valores normalizados de las bandas de frecuencia de la VFC calculadas a partir de la CWT-Morlet, durante la prueba de peso sostenido (300 s) e interpolados a 2 Hz (0,5 s).
RESULTADOS
En la tabla se puede observar que dentro de los parámetros hemodinámicos estudiados solo presentaron diferencias significativas la PA media, sistólica y diastólica, aunque la frecuencia cardíaca tuvo una tendencia al aumento.
Se muestra en la figura 1 la reducción significativa de la entropía de Shannon durante la prueba del peso sostenido en pacientes normotensos.
En la figura 2 se puede observar como hay un incremento del componente simpático de la respuesta autonómica cardíaca cerca del minuto 2, al mismo tiempo que hay un descenso simultáneo del componente parasimpático.
En la figura 3 se puede apreciar que hay un incremento del componente simpático de la respuesta autonómica cardíaca, por lo que predomina su actividad sobre el componente parasimpático alrededor del minuto 2, que luego decrece.
DISCUSIÓN
La presencia de variaciones significativas en cuanto a las presiones arteriales sistólica, diastólica y media que se muestran en la tabla coincide con la bibliografía de referencia26-28. En estudios anteriores en los que se empleó la prueba de esfuerzo isométrico como ya se mencionaba, hubo variaciones en cuanto a los valores de dichas variables producto al efecto de este estímulo estresante28. No se puede dejar de señalar que aunque la frecuencia cardíaca no mostró cambios significativos, sí aumentó. Para comprender este comportamiento es necesario partir de que se mide la frecuencia cardíaca durante 5 minutos, una vez comenzada la PPS, pero esta solo dura 2 minutos por lo que los otros tres restantes responden a un período de recuperación. Todo esto trae consigo que aunque existan variaciones respecto a la frecuencia cardíaca, no serán significativas al promediarla en los 5 minutos correspondientes a la PPS y al período de recuperación, ya que la segunda etapa falseará el promedio real de la primera, por lo que sería recomendable determinar el promedio del período de la PPS y el de recuperación por separado.
Goulopoulou et al29, plantean que ocurren cambios en la frecuencia cardíaca durante los ejercicios estáticos, entre los que se encuentra la PPS. En estos casos los mecanismos neurales postulados son dos: uno que activa los circuitos neuronales centrales que controlan el sistema cardiovascular y el somatomotor, y van a producir cambios en la actividad eferente simpática y parasimpática, lo cual provocará cambios en las respuestas cardiovasculares. El otro mecanismo plantea que los cambios en la actividad eferente autonómica son causa refleja por la estimulación de las fibras musculares de tipo III, las cuales producen un estímulo inhibitorio a nivel del núcleo del tracto solitario que produce una disminución del influjo vagal y un aumento de la frecuencia cardíaca (respuesta temprana mecanorrefleja muscular al ejercicio isométrico)30,31.
Es importante señalar que Lellamo et al32 también informaron que el ejercicio estático produce estimulación de la actividad simpática, lo que trae consigo elevaciones de la PA sistólica, diastólica y la frecuencia cardíaca. Esto reafirma que la excitación simpática mecanorrefleja muscular contribuye a la regulación de la frecuencia cardíaca durante el ejercicio estático y coincide con los resultados obtenidos por Bunsawat y Baynard33en el grupo control integrado por sujetos sanos, en quienes también hubo elevación de los valores de la PA sistólica, diastólica y media durante el ejercicio isométrico.
El caos intrínsecamente aprovecha la riqueza relacionada con su estructura, es por esto que hay beneficios para estos sistemas de la adopción de regímenes caóticos con una amplia gama de conductas posibles34. Así, hay trabajos que muestran que con el envejecimiento se reduce la entropía, al igual que los sistemas patológicos muestran entropías menores que los sanos35,36, que desde la perspectiva de la teoría de la complejidad se pueden interpretar como una pérdida de la adaptabilidad de los sistemas reguladores del ritmo cardíaco. También se ha demostrado cómo la entropía es una medida alternativa no-lineal del balance simpático-vagal, la cual disminuye en pruebas que implican un incremento de la respuesta simpática (prueba de la mesa basculante y estrés mental)36,37, mientras que se produce un aumento de la entropía en condiciones salutogénicas de predominio vagal, como la práctica regular de ejercicios físicos38, la hipnosis39 y la meditación40. Esto es congruente con lo encontrado en este trabajo pues en la figura 1 se muestra como se produce una disminución de la entropía de Shannon en la PPS con respecto al estado de reposo en sujetos normotensos, ya que esta prueba provoca un incremento del tono simpático y de las sustancias presoras, que generan que el sistema se vuelva más rígido, con tendencia hacia un punto de máxima activación.
En la figura 2, con la aplicación de métodos de tiempo-frecuencia, se observa claramente un aumento de la actividad simpática barorrefleja (LF), simultánea a una depresión en la actividad del sistema nervioso parasimpático (HF). Esto se manifiesta en los 2 minutos en los que se realiza la PPS que posteriormente se normaliza en el transcurso del período de recuperación32. Por otra parte, en esa misma figura se aprecia como existe un predominio simpático en el balance simpático vagal (LF/HF) coincidentemente en los dos primeros minutos, durante la PPS, y alcanzando su máximo valor finalizado el minuto 2.
Estos resultados coinciden con los de Tiinanen et al41 que en un estudio realizado en 11 hombres sanos, utilizaron el ejercicio de la prueba de esfuerzo isométrico, y encontraron que en los últimos dos minutos se alcanzó la mayor respuesta simpática, traducida por la elevación de la PA sistólica, LF, así como en la disminución del parámetro HF. Pero los valores más altos registrados en estos parámetros se observaron en el último minuto, por lo que se considera que es precisamente en ese último minuto de registro en los ejercicios isométricos que se puede observar con más intensidad la actividad simpática sobre el sistema cardiovascular.
Con este estudio, al aplicar métodos de tiempo-frecuencia de la VFC, se evidenciaron los cambios dinámicos en el tiempo presentes durante la PPS en el sistema nervioso autónomo, lo que posibilita caracterizar el patrón normal de respuesta durante la PPS en sujetos normotensos y normorreactivos.
CONCLUSIONES
La dinámica de la respuesta autonómica cardíaca durante la prueba de peso sostenido se caracteriza por un incremento en la respuesta simpática cardiovascular, que se hace evidente al finalizar el segundo minuto, acompañado además de elevaciones en la presión arterial, seguida de una recuperación inmediata con activación parasimpática. La respuesta autonómica cardiovascular durante el ejercicio isométrico cubano está asociada a una disminución de la entropía y por ende, de la complejidad del proceso fisiológico, lo cual no es evidente con los métodos lineales tradicionales de la variabilidad de la frecuencia cardíaca.