La presión del pulso en las arterias cerebrales genera una carga cíclica sobre las paredes y el tejido circundante. Las mediciones por resonancia magnética muestran que el cerebro se desplaza fracciones de milímetro con cada latido, pero las imágenes por sí solas no bastan para separar la contribución de la geometría del vaso, las propiedades de la pared y las características del tejido. La mayoría de los trabajos computacionales abordan o bien el flujo dentro del vaso, o bien la mecánica del cerebro en su conjunto, y hasta ahora no se había construido un modelo totalmente acoplado con anatomía realista. Este trabajo llena ese vacío utilizando geometría de resonancia magnética de alto campo y un cálculo acoplado.
A partir de imágenes de resonancia magnética de 5 teslas de una mujer sana de unos 30 años se segmentaron cuatro regiones: la arteria cerebral media (M1), su bifurcación, la arteria basilar y la arteria carótida interna. Para cada una se construyó un modelo tridimensional de la luz, la pared de grosor variable y el tejido circundante en forma de cubo con tres capas concéntricas alrededor del vaso. La sangre se describió con el modelo de Carreau para el adelgazamiento por cizallamiento, la pared con el modelo hiperelástico de Ogden y el tejido con el modelo hiperviscoelástico de Storakers con varios tiempos de relajación. La interacción sangre-pared-tejido se resolvió por el método de elementos finitos en COMSOL con acoplamiento bidireccional y malla móvil. En la entrada se aplicaron perfiles fisiológicos pulsátiles de caudal y presión, y los resultados se compararon con datos de resonancia magnética de la literatura sobre la deformación volumétrica de la sustancia gris.
Las simulaciones reprodujeron la evolución temporal de la deformación volumétrica: aumento en sístole, pico alrededor del 30% del ciclo y descenso en diástole. La deformación es máxima en la pared y se atenúa rápidamente al alejarse del vaso, manteniéndose dentro de los valores medidos experimentalmente para la sustancia gris. La arteria carótida interna dio la mayor deformación debido a la combinación de alta velocidad de flujo y marcada tortuosidad, mientras que la arteria basilar mostró un patrón más simétrico y local. El grosor de la pared actúa de forma inversamente proporcional: con el grosor mínimo, la deformación volumétrica en la primera capa de tejido alcanzó 7,75×10⁻³, y con el máximo, 1,68×10⁻³. Los parámetros hemodinámicos junto a la pared (TAWSS, OSI, RRT) resultaron heterogéneos: los valores elevados aparecían en bifurcaciones, curvaturas y estrechamientos, pero en conjunto el flujo se mantuvo unidireccional.
El trabajo reúne por primera vez geometría realista de arterias cerebrales, reología sanguínea no lineal, paredes hiperelásticas y tejido viscoelástico en un único modelo acoplado con validación mediante datos de resonancia magnética. Esto permite separar la contribución de la forma del vaso de las propiedades del tejido y de la pared, algo que no ofrecen ni los modelos puramente hemodinámicos ni los puramente tisulares. Los resultados subrayan que la deformación y las perturbaciones hemodinámicas están gobernadas por la geometría de maneras distintas, por lo que evaluar el estado de un vaso exige analizar conjuntamente los índices mecánicos y los parietales.
El modelo puede ampliarse a una anatomía realista con líquido cefalorraquídeo, meninges y espacios perivasculares, así como a geometrías específicas de cada paciente en casos de aneurismas, estenosis o cambios de rigidez asociados a la edad. La incorporación de procesos de transporte permitirá estudiar cómo la deformación pulsátil modula el transporte de fármacos en el espacio intercelular y perivascular.
Los resultados serán útiles en neurocirugía y neurorradiología para evaluar el riesgo de aneurismas y estenosis, y también en farmacocinética para predecir la distribución de fármacos en el cerebro.
El siguiente paso es la validación en un mayor número de sujetos y la comparación con datos de DENSE MRI en forma espacialmente coregistrada. Después se puede pasar a modelar alteraciones patológicas de la geometría y de las propiedades del tejido.
En la mecánica del tejido cerebral sigue abierta la pregunta de cómo se reparten las cargas pulsátiles entre la red vascular, el líquido cefalorraquídeo y el parénquima, y cómo esto afecta al transporte y al daño en patologías. Este modelo acerca la respuesta, pero no tiene en cuenta el líquido cefalorraquídeo ni la anisotropía del tejido.
El cerebro se desplaza aproximadamente entre 0,1 y 0,3 mm con cada latido, una magnitud comparable al grosor de un cabello humano, pero esa deformación basta para influir en el transporte de sustancias en los espacios perivasculares.
No publicado en arXiv. Aquí están su página, el PDF y la formalización: puede comprobarlo usted mismo, no solo leer nuestro análisis.