🔢 mathematics

A coupled geometric surface and surface reaction-diffusion framework for modelling nuclear mechanics during confined cell migration

Este artículo presenta un marco acoplado de superficie geométrica y de reacción-difusión superficial que modela con éxito cómo la regulación epigenética de H3K9me3 por SETDB1 y SUV39H1 influye en la mecánica nuclear y en los tiempos de entrada de migración en células de cáncer de pulmón confinadas.

Autores originales: Francesca Ballatore, Silvia Comunian, Anotida Madzvamuse, Slimane Ait-Si-Ali, Rachele Allena

Publicado 2026-10-05

📖 6 min de lectura🧠 Análisis profundo

CC BY 4.0

Autores originales: Francesca Ballatore, Silvia Comunian, Anotida Madzvamuse, Slimane Ait-Si-Ali, Rachele Allena

Artículo original bajo licencia CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ Esta es una explicación generada por IA del artículo a continuación. No ha sido escrita ni avalada por los autores. Para mayor precisión técnica, consulte el artículo original. Leer descargo de responsabilidad completo

Las células son los bloques fundamentales de la vida, pero no son ladrillos estáticos; son entidades dinámicas y blandas que deben moverse constantemente para cumplir con sus funciones. En el cuerpo humano, este movimiento es a menudo una cuestión de supervivencia. Las células inmunitarias se cuelan a través de estrechos huecos para cazar infecciones y, durante el desarrollo, las células migran para formar nuevos tejidos. Sin embargo, cuando estas células encuentran espacios más estrechos que su propio ancho, se enfrentan a un obstáculo mecánico significativo. La célula en sí es flexible, pero en su interior se encuentra el núcleo, un orgánulo grande y denso que contiene el material genético de la célula. El núcleo es típicamente la parte más rígida de la célula, actuando como un núcleo rígido que resiste la deformación. Cuando una célula intenta pasar por una abertura estrecha, el núcleo suele ser el cuello de botella, y la facilidad con la que puede ser aplastado determina si la célula puede pasar o se queda atascada.

Los científicos saben desde hace tiempo que la rigidez del núcleo no es fija; puede cambiar según cómo se empaqueta el ADN en su interior. El ADN se enrolla alrededor de proteínas llamadas histonas, y se pueden añadir etiquetas químicas a estas proteínas para apretar o aflojar el empaquetado. Una de estas etiquetas, conocida como H3K9me3, actúa como una abrazadera molecular que empaqueta el ADN de forma apretada, haciendo que el núcleo sea más rígido. Los investigadores han observado que en ciertas células de cáncer de pulmón, el equilibrio de las enzimas que añaden o eliminan estas etiquetas puede alterar drásticamente la rigidez del núcleo. Cuando la enzima SETDB1 está presente en grandes cantidades, el núcleo tiende a ser más blando y la célula se mueve más rápido a través de espacios estrechos. Cuando se elimina esta enzima, el núcleo se vuelve mucho más rígido y la célula tiene dificultades para moverse. La pregunta sigue siendo: ¿cómo se traduce exactamente este cambio molecular a nivel microscópico en el movimiento físico de toda la célula a través de un túnel microscópico?

Para responder a esto, un equipo de investigadores desarrolló un nuevo modelo computacional que simula el viaje de una célula a través de un canal estrecho. En lugar de mirar solo la célula como un todo, construyeron un marco que conecta la actividad molecular en la superficie del núcleo con las fuerzas físicas que actúan sobre la membrana celular. Imaginaron la célula y su núcleo como dos pieles distintas y flexibles que evolucionan y cambian de forma con el tiempo. La piel exterior, la membrana celular, y la piel interior, la envoltura nuclear, están gobernadas por leyes físicas que tienen en cuenta la tensión, la resistencia a la flexión y la fricción. Crucialmente, el modelo incluye un sistema que imita las reacciones químicas que ocurren precisamente en la superficie del núcleo. Este sistema rastrea los niveles de las etiquetas moleculares que controlan el empaquetamiento del ADN. A medida que el modelo se ejecuta, el estado químico del núcleo influye directamente en sus propiedades mecánicas, como cuánto resiste ser comprimido.

Los investigadores utilizaron este modelo para simular dos escenarios específicos basados en células reales de cáncer de pulmón. En el primer escenario, modelaron células de control, que tienen un equilibrio normal de enzimas y un núcleo relativamente blando. En el segundo, modelaron células a las que se les había eliminado la enzima SETDB1. En el mundo real, la eliminación de esta enzima provoca una acumulación de las etiquetas moleculares que endurecen el núcleo en la superficie nuclear. La simulación por computadora reflejó esta realidad biológica. Mientras el modelo se ejecutaba, las células simuladas eran empujadas hacia un canal estrecho, tal como ocurriría en un experimento de laboratorio. Los resultados mostraron una clara diferencia en el comportamiento. Las células de control, con sus núcleos más blandos, pudieron deformarse y deslizarse a través del canal con relativa rapidez. Las células con la enzima ausente, sin embargo, tenían núcleos que se volvieron significativamente más rígidos debido a la acumulación de etiquetas moleculares. Esta mayor rigidez hizo que el núcleo fuera más difícil de comprimir, lo que causó que toda la célula se moviera más lentamente a través de la constricción.

Las simulaciones proporcionaron números específicos que podrían compararse con mediciones de laboratorio reales. En experimentos reales, las células de control tardaron un promedio de 1,44 segundos en entrar en un canal estrecho, mientras que las células con la enzima ausente tardaron 5,01 segundos. El modelo computacional no coincidió perfectamente con estos tiempos exactos; predijo que las células de control tardarían 2,17 segundos y las células modificadas tardarían 3,22 segundos. Aunque el modelo no reprodujo la magnitud completa del retraso observado en el laboratorio, capturó con éxito la dirección del efecto. La simulación confirmó que la mayor rigidez del núcleo por sí sola es suficiente para ralentizar la célula, demostrando que los cambios moleculares en la superficie del núcleo tienen una consecuencia mecánica directa. Sin embargo, el hecho de que el modelo no replicara completamente el retraso extremo observado en los experimentos sugiere que la historia es más compleja. Los investigadores descubrieron que, si bien la rigidez nuclear es un factor importante, no es el único. La eliminación de la enzima probablemente afecta a otras partes de la célula también, contribuyendo a la ralentización general de una manera que el modelo actual, que se centra principalmente en el núcleo, aún no logra capturar por completo.

Este trabajo representa un paso significativo en la comprensión de cómo el mundo invisible de la química molecular impulsa el mundo visible del movimiento celular. Al vincular un sistema de reacción superficial con las leyes geométricas del movimiento, los investigadores crearon una herramienta que puede visualizar cómo un cambio en una sola etiqueta química puede propagarse para cambiar el comportamiento de una célula entera. El modelo demuestra que el núcleo actúa como un guardi de control mecánico, donde su organización interna dicta con qué facilidad una célula puede navegar por los espacios congestionados y estrechos del cuerpo. Aunque la simulación es una representación simplificada, ofrece una prueba de concepto clara: la regulación molecular de la organización de la cromatina es un motor clave de la mecánica nuclear, y estas propiedades mecánicas son esenciales para comprender cómo las células migran en entornos confinados. Los hallazgos sugieren que para comprender plenamente el movimiento celular en enfermedades como el cáncer, los científicos deben mirar no solo la cubierta exterior de la célula, sino también el estado químico de su núcleo.

Resumen Técnico: Un marco acoplado de superficie geométrica y reacción-difusión superficial para modelar la mecánica nuclear durante la migración celular confinada

Planteamiento del problema
La migración celular confinada es un proceso crítico en contextos fisiológicos y patológicos, como la metástasis del cáncer, donde las células deben deformarse para pasar a través de poros más pequeños que su diámetro. El núcleo, siendo el orgánulo más grande y rígido, suele actuar como la principal barrera mecánica para esta migración. Si bien la evidencia experimental vincula las modificaciones epigenéticas —específicamente la trimetilación de la histona H3 en la lisina 9 (H3K9me3)— con la rigidez nuclear y la eficiencia de la migración, aún falta un marco teórico cuantitativo que conecte la regulación molecular en la periferia nuclear con la mecánica celular macroscópica. Específicamente, existe la necesidad de modelar cómo las variables de estado molecular (por ejemplo, el equilibrio entre las enzimas SETDB1 y SUV39H1) modulan dinámicamente las propiedades mecánicas nucleares (tensión superficial, rigidez de flexión, fricción) e influyen, consecuentemente, en la translocación de la célula completa a través de espacios confinados.

Metodología
Los autores proponen un marco matemático totalmente acoplado en dos dimensiones que integra la evolución de la superficie geométrica con la dinámica de reacción-difusión superficial.

Balance de fuerzas de superficie geométrica (GS-PDEs):
- La membrana plasmática celular ( $\Gamma_c$ ) y la envoltura nuclear ( $\Gamma_n$ ) se modelan como curvas cerradas en evolución, gobernadas por leyes de balance de fuerzas por unidad de área.
- Las ecuaciones gobernantes incorporan tensión superficial, elasticidad de flexión (modelo de Helfrich), restricciones de área y perímetro, fuerzas de autorrepulsión (para evitar la autointersección) e interacciones externas (fuerzas de barrera de las paredes del canal, presión motriz y acoplamiento elástico entre el núcleo y la membrana celular).
- El sistema se resuelve utilizando el Método de Elementos Finitos de Superficie Evolutiva (ESFEM), que discretiza los límites móviles como variedades poligonales en evolución, permitiendo grandes deformaciones y cambios geométricos complejos.
Sistema de Reacción-Difusión Superficial (S-RDEs):
- Se plantea un sistema de ecuaciones diferenciales parciales directamente sobre la superficie nuclear evolutiva $\Gamma_n(t)$ para describir la dinámica de las variables de estado molecular adimensionales: $E$ (SETDB1), $U$ (SUV39H1) y $W$ (H3K9me3).
- Los términos de reacción siguen una ley de acción de masas fenomenológica que representa la cascada regulatoria SETDB1–SUV39H1–H3K9me3, donde SETDB1 inhibe a SUV39H1, y SUV39H1 promueve la acumulación de H3K9me3.
- El modelo asume una vía de alimentación hacia adelante (feed-forward) sin retroalimentación de H3K9me3 hacia las enzimas, basándose en la caracterización experimental actual.
Acoplamiento Mecano-Químico:
- Los parámetros mecánicos del núcleo —el coeficiente de tipo fricción ( $\omega_n$ ), la tensión superficial ( $k^s_n$ ) y la rigidez de flexión ( $k^b_n$ )— se definen como funciones explícitas del valor local de H3K9me3 ( $W$ ).
- Estas funciones se calibran utilizando datos experimentales de pinzas ópticas (rigidez nuclear) y ensayos microfluídicos (tiempos de entrada) para células de cáncer de pulmón de control (Ctrl) y con inactivación de SETDB1 (SETDB1 KO).
- En la simulación, los parámetros nucleares dependen del promedio espacial de $W$ , creando un bucle de retroalimentación donde la distribución molecular dicta la resistencia mecánica.
Implementación Numérica:
- El sistema acoplado se resuelve utilizando el marco de código abierto FEniCSx.
- Se emplea un esquema de Euler hacia atrás semi-implícito con iteración de Picard para la discretización temporal.
- Se introducen componentes de velocidad tangencial para mantener la calidad de la malla sin alterar la forma física de las membranas.

Contribuciones Clave

Marco Novedoso: El artículo introduce un marco computacional general que acopla las PDE de superficie geométrica con sistemas de reacción-difusión superficial en dominios en evolución. Esto permite investigar cómo los procesos moleculares espacio-temporales en la periferia nuclear modulan el comportamiento mecánico global.
Modelado Fenomenológico de la Epigenética: El modelo traduce con éxito el mecanismo biológico específico de la regulación SETDB1/SUV39H1 en un sistema matemático que predice estados estacionarios distintos para la distribución de H3K9me3 en las condiciones de control (Ctrl) y de SETDB1 KO.
Calibración Cuantitativa: El marco se calibra contra conjuntos de datos experimentales específicos (pinzas ópticas y tiempos de entrada microfluídicos) de células de cáncer de pulmón de células no pequeñas (NSCLC), uniendo observaciones a escala molecular con métricas de migración a escala celular.

Resultados

Dinámica Molecular: El sistema de reacción-difusión converge a estados estacionarios estables distintos dependiendo de las condiciones iniciales. En el escenario de control (Ctrl), el sistema alcanza un nivel asintótico inferior de $W$ (H3K9me3). En el escenario de SETDB1 KO (donde $E=0$ ), el sistema converge a un nivel asintótico superior de $W$ .
Respuesta Mecánica: Los niveles más altos de $W$ en la condición SETDB1 KO conducen a un aumento de la fricción nuclear, la tensión superficial y la rigidez de flexión, lo que resulta en una respuesta nuclear más rígida en comparación con la condición de control.
Dinámica de Migración:
- Células de Control (Ctrl): El tiempo de entrada simulado en un microcanal de 6 µm es de 2.17 s.
- Células SETDB1 KO: El tiempo de entrada simulado es de 3.22 s.
- Comparación: Aunque la simulación reproduce la dirección de la tendencia experimental (las células SETDB1 KO son más lentas), la magnitud de la diferencia es menor que la observada experimentalmente (Experimental: Ctrl 1.44 s vs. KO 5.01 s).
Heterogeneidad Espacial: El estudio probó distribuciones iniciales espacialmente heterogéneas de las variables moleculares. Los resultados mostraron que la redistribución molecular ocurre en una escala de tiempo significativamente más corta que la entrada de la célula, lo que conduce a respuestas mecánicas efectivamente homogéneas independientemente de la heterogeneidad espacial inicial.

Significancia y Afirmaciones
Los autores afirman que este marco proporciona un primer enfoque descriptivo para vincular cuantitativamente las alteraciones moleculares en la envoltura nuclear con la mecánica nuclear y la migración celular confinada. El modelo reproduce con éxito el comportamiento biológico cualitativo observado experimentalmente: la depleción de SETDB1 conduce al aumento de H3K9me3, lo que incrementa la rigidez nuclear y retrasa la migración.

Sin embargo, el artículo es modesto respecto al acuerdo cuantitativo. Los autores establecen explícitamente que la diferencia residual entre los tiempos de entrada simulados y medidos indica que la respuesta celular es multifactorial. El modelo aísla la contribución de la mecánica nuclear, sugiriendo que si bien el endurecimiento nuclear contribuye al retraso, es insuficiente por sí solo para explicar la diferencia experimental completa. Los autores concluyen que probablemente se requieren cambios en los parámetros mecánicos de toda la célula (más allá de solo el núcleo) para reproducir plenamente el fenotipo experimental.

El marco se presenta como una plataforma flexible que puede adaptarse a otros mecanismos de regulación molecular en las interfaces celulares y extenderse en el futuro para incluir el acoplamiento entre el volumen y la superficie, la retroalimentación mecano-química bidireccional y geometrías tridimensionales.

Más como este