Dans la cellule, les réactions biochimiques sont organisées non seulement par les organites membranaires, mais aussi par des condensats sans membrane, qui se forment par séparation de phases. Si les étapes initiales de la formation de ces gouttelettes sont bien connues, leur évolution à long terme, appelée vieillissement, reste mal comprise. Le vieillissement décrit le passage de gouttelettes liquides à des états viscoélastiques, gélatineux, voire solides, sur des durées allant de quelques minutes à plusieurs mois. Il en résulte un ralentissement de la diffusion des molécules, une perte de la capacité des gouttelettes à fusionner et l'apparition d'une hétérogénéité conformationnelle, comme dans les verres. Une augmentation modérée de la viscosité peut stabiliser des assemblages fonctionnels, mais un durcissement excessif est associé aux troubles neurodégénératifs et au cancer. Ce travail se situe dans le prolongement des recherches sur le vieillissement cellulaire et l'organisation des gouttelettes dans la cellule, mais il se concentre sur les mécanismes moléculaires qui n'étaient jusqu'ici que supposés.
Les auteurs ont utilisé un modèle d'ARN à gros grains, dans lequel chaque nucléotide est représenté par une seule particule conservant l'identité de sa base. Le modèle comprend des liaisons harmoniques le long du squelette, une répulsion de volume exclu et un potentiel multiparticulaire pour les paires de bases Watson-Crick. Les simulations ont été menées sur 64 chaînes (CAG)31 dans une boîte périodique à une concentration de 400 µM, en dynamique de Langevin avec un pas de 10 fs à 300 K. Pour analyser le vieillissement, seules les molécules restant en permanence à l'intérieur de la gouttelette ont été suivies, et la trajectoire a été découpée en 27 fenêtres. La dynamique a été quantifiée au moyen de la fonction de recouvrement, du déplacement quadratique moyen par paires et des temps de relaxation extraits d'un ajustement biexponentiel. L'hétérogénéité spatiale a été étudiée en divisant la gouttelette en couronnes concentriques de 5 nm d'épaisseur. La topologie du réseau d'enchevêtrement a été analysée par la méthode des chemins primitifs. À titre de comparaison, un ARN de séquence aléatoire de même composition a été simulé.
Les simulations ont montré que le condensat (CAG)31 vieillit spontanément : la fonction de recouvrement ne retombe pas à zéro, et la relaxation suit une loi de puissance, ce qui indique une ergodicité incomplète. L'ajustement biexponentiel a révélé un régime rapide (τ1) et un régime lent (τ2) de relaxation, qui augmentent avec l'âge de la gouttelette, τ1 et τ2 différant de deux ordres de grandeur. La fraction gelée A0 croît avec le temps. L'analyse par couronnes a montré que le vieillissement est spatialement hétérogène : dans le cœur, la diffusion ralentit le plus vite (exposant de vieillissement αw,core = 0,53 ± 0,11), et en périphérie plus lentement (αw,outer = 0,26 ± 0,15). Dans le cœur, la fraction de contacts non relaxants atteint 30 %, et les temps de relaxation y sont six fois plus grands qu'en surface. Les chaînes d'ARN du cœur sont étendues (rayon de giration d'environ 7,0 à 7,5 nm contre 5,5 à 6,0 nm en périphérie), présentent un exposant de Flory ν > 0,5 et sont partiellement orientées radialement. L'analyse des chemins primitifs a montré une croissance du nombre de nœuds d'enchevêtrement de 5 à 8 au cours de la simulation, la densité de nœuds dans le cœur étant cinq fois plus élevée que dans la couronne. La densité de paires de bases Watson-Crick dans le cœur est quatre fois plus élevée qu'en périphérie, et elle est constituée principalement de paires interbrins (19 à 20 paires par chaîne). La durée de vie des paires et des contacts entre chaînes augmente avec l'âge. La comparaison avec une séquence aléatoire a montré que celle-ci vieillit 4 à 5 fois plus lentement, présente une distribution des paires de bases presque uniforme et un enchevêtrement moindre, malgré une densité similaire.
Ce travail montre que le vieillissement des condensats d'ARN ne nécessite ni stimulus externe ni réglage particulier des paramètres : il est encodé dans la séquence et résulte des interactions interbrins. Cela le distingue des études voisines, qui portaient souvent sur le vieillissement de condensats protéiques ou utilisaient des modèles à paramètres ajustables. Ici, on montre pour la première fois que l'hétérogénéité spatiale du vieillissement découle d'un gradient d'appariement des bases, et non d'une structure prédéfinie. Les résultats relient les interactions moléculaires aux transitions matérielles mésoscopiques et proposent un principe général : une boucle de rétroaction dans laquelle l'appariement des bases étend les chaînes, l'extension accroît l'enchevêtrement, l'enchevêtrement ralentit la relaxation, et ce ralentissement renforce le réseau.
Le sujet pourrait s'orienter vers l'étude de séquences plus longues, de l'influence des modifications de l'ARN, de l'activité des hélicases et de l'évolution des conditions ioniques sur le vieillissement. Les approches multi-échelles, combinant modélisation à gros grains et expériences, aideront à tester les prédictions sur des échelles de temps biologiquement pertinentes. Il importe aussi de déterminer si les cellules peuvent réguler le vieillissement par la conception de leurs séquences ou par des processus actifs, afin de maintenir les condensats dans un état fonctionnel.
Ces résultats sont importants pour comprendre les maladies neurodégénératives liées à l'expansion des répétitions CAG, où les foyers d'ARN peuvent durcir de façon irréversible. Ils fournissent aussi une base pour concevoir des condensats d'ARN synthétiques aux propriétés choisies.
Vérifier les prédictions par des expériences sur des gouttelettes d'ARN, par exemple au moyen de la FRAP et de la microrhéologie. Étudier l'influence d'autres séquences répétées et des conditions du milieu sur le vieillissement.
Dans le domaine des condensats biomoléculaires, la question de savoir comment la séquence détermine les propriétés matérielles et leur évolution reste ouverte. Ce travail relie le vieillissement à des interactions intermoléculaires précises, mais il ne résout pas le problème de la réversibilité du vieillissement dans la cellule ni celui du rôle des processus actifs, comme l'action des hélicases dépendantes de l'ATP.
Dans le cœur de la gouttelette, les temps de relaxation sont six fois plus grands qu'en périphérie, et la densité de nœuds d'enchevêtrement y est cinq fois plus élevée. Pourtant, les deux séquences ont la même composition en nucléotides, mais un ordre différent, et cela suffit pour que l'une vieillisse cinq fois plus vite que l'autre.
Non publié sur arXiv. Voici sa page, le PDF et la formalisation : vous pouvez vérifier vous-même, pas seulement lire notre analyse.