En la célula, las reacciones bioquímicas no solo están organizadas por orgánulos membranosos, sino también por condensados sin membrana que se forman por separación de fases. Aunque las etapas iniciales de la formación de estas gotas están bien estudiadas, su evolución a largo plazo, llamada envejecimiento, sigue siendo poco conocida. El envejecimiento describe la transición de gotas similares a líquidos a estados viscoelásticos, gelatinosos o incluso sólidos en tiempos que van de minutos a meses. Esto provoca una ralentización de la difusión de las moléculas, la pérdida de la capacidad de las gotas para fusionarse y la aparición de heterogeneidad conformacional, como ocurre en los vidrios. Un aumento moderado de la viscosidad puede estabilizar ensamblajes funcionales, pero el endurecimiento excesivo se asocia con trastornos neurodegenerativos y cáncer. El trabajo se sitúa junto a las investigaciones sobre el envejecimiento celular y la organización de las gotas en la célula, pero se centra en los mecanismos moleculares que antes solo se suponían.
Los autores utilizaron un modelo de ARN de grano grueso en el que cada nucleótido está representado por una partícula que conserva la identidad de su base. El modelo incluye enlaces armónicos a lo largo del esqueleto, repulsión de volumen excluido y un potencial multipartícula para los pares de bases Watson-Crick. Las simulaciones se realizaron para 64 cadenas (CAG)31 en una caja periódica a una concentración de 400 µM, usando dinámica de Langevin con un paso de 10 fs a 300 K. Para analizar el envejecimiento se siguieron únicamente las moléculas que permanecían dentro de la gota y se dividió la trayectoria en 27 ventanas. La dinámica se cuantificó mediante la función de solapamiento, el desplazamiento cuadrático medio por pares y los tiempos de relajación extraídos de un ajuste biexponencial. La heterogeneidad espacial se estudió dividiendo la gota en capas concéntricas de 5 nm de grosor. La topología de la red de enmarañamiento se analizó con el método de caminos primitivos. Como comparación, se simuló ARN con una secuencia aleatoria de la misma composición.
Las simulaciones mostraron que el condensado (CAG)31 envejece de forma espontánea: la función de solapamiento no decae hasta cero y la relajación se describe con una ley de potencias, lo que indica una ergodicidad incompleta. El ajuste biexponencial reveló un régimen de relajación rápido (τ1) y otro lento (τ2) que aumentan con la edad de la gota, y τ1 y τ2 difieren en dos órdenes de magnitud. La fracción congelada A0 crece con el tiempo. El análisis por capas mostró que el envejecimiento es espacialmente heterogéneo: en el núcleo la difusión se ralentiza más que en ninguna otra zona (exponente de envejecimiento αw,core = 0.53 ± 0.11), mientras que en la periferia lo hace más despacio (αw,outer = 0.26 ± 0.15). En el núcleo, la fracción de contactos que no se relajan alcanza el 30% y los tiempos de relajación son seis veces mayores que en la superficie. Las cadenas de ARN del núcleo están expandidas (radio de giro de alrededor de 7.0–7.5 nm frente a 5.5–6.0 nm en la periferia), tienen un exponente de Flory ν > 0.5 y están parcialmente orientadas de forma radial. El análisis de caminos primitivos mostró un aumento del número de nudos de enmarañamiento de 5 a 8 durante el tiempo de simulación, y la densidad de nudos en el núcleo es cinco veces mayor que en la corteza. La densidad de pares de bases Watson-Crick en el núcleo cuadruplica la de la periferia y se compone principalmente de pares entre cadenas (19–20 pares por cadena). El tiempo de vida de los pares y de los contactos entre cadenas aumenta con la edad. La comparación con una secuencia aleatoria mostró que esta envejece entre 4 y 5 veces más despacio, tiene una distribución de pares de bases casi uniforme y un menor enmarañamiento, pese a una densidad similar.
El trabajo muestra que el envejecimiento de los condensados de ARN no requiere agentes externos ni un ajuste especial de los parámetros: está codificado en la secuencia y surge de las interacciones entre cadenas. Esto lo distingue de investigaciones cercanas, que a menudo estudiaban el envejecimiento de condensados de proteínas o usaban modelos con parámetros ajustables. Aquí se demuestra por primera vez que la heterogeneidad espacial del envejecimiento surge de un gradiente de apareamiento de bases y no de una estructura predefinida. Los resultados unen las interacciones moleculares con las transiciones materiales mesoscópicas y proponen un principio general: un bucle de retroalimentación en el que el apareamiento de bases expande las cadenas, la expansión aumenta el enmarañamiento, el enmarañamiento ralentiza la relajación y la ralentización refuerza la red.
El tema puede avanzar hacia el estudio de secuencias más largas, la influencia de las modificaciones del ARN, la actividad de las helicasas y los cambios en las condiciones iónicas sobre el envejecimiento. Los enfoques multiescala que combinan el modelado de grano grueso con experimentos ayudarán a comprobar las predicciones en escalas de tiempo biológicamente relevantes. También es importante investigar si las células pueden regular el envejecimiento mediante el diseño de la secuencia o procesos activos, para mantener los condensados en un estado funcional.
Los resultados son importantes para comprender las enfermedades neurodegenerativas asociadas a la expansión de repeticiones CAG, en las que los focos de ARN pueden endurecerse de forma irreversible. También ofrecen una base para el diseño de condensados de ARN sintéticos con propiedades a medida.
Comprobar las predicciones en experimentos con gotas de ARN, por ejemplo mediante FRAP y microrreología. Investigar la influencia de otras secuencias repetidas y de las condiciones del entorno sobre el envejecimiento.
En el campo de los condensados biomoleculares sigue abierta la pregunta de cómo la secuencia determina las propiedades materiales y su evolución. Este trabajo vincula el envejecimiento con interacciones intermoleculares concretas, pero no resuelve el problema de la reversibilidad del envejecimiento en la célula ni el papel de procesos activos como la acción de las helicasas dependientes de ATP.
En el núcleo de la gota los tiempos de relajación son seis veces mayores que en la periferia, y la densidad de nudos de enmarañamiento es cinco veces mayor. Y sin embargo, ambas secuencias tienen la misma composición de nucleótidos, pero distinto orden, y eso basta para que una envejezca cinco veces más rápido que la otra.
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.