Le problème du repliement du génome dans la capside virale se situe à l'interface de la physique, de la chimie et de la virologie. Des travaux voisins dans les archives explorent des aspects connexes : comment les gènes bactériens s'assemblent en gouttelettes vivantes, comment les transistors ioniques fonctionnent dans les nanopores, comment les effets quantiques se manifestent dans les systèmes biologiques. Cet article se distingue en incluant pour la première fois de manière systématique, dans une dynamique moléculaire, à la fois des ions explicites, de l'eau explicite, des sites de charge discrets sur la capside et une électrostatique coulombienne complète. Auparavant, de tels calculs utilisaient soit un solvant implicite, soit des modèles de polymère simplifiés, soit des conditions ioniques fixes. Ici, l'écrantage ionique émerge naturellement de la dynamique des ions, ce qui permet de tester les prédictions des théories de champ moyen au niveau moléculaire.
Les auteurs ont construit un modèle à gros grains dans lequel l'ARN est représenté par une chaîne linéaire de 1122 billes de trois types (phosphate, sucre, sucre) basée sur le champ de force MARTINI-2. La capside est modélisée comme une double coque sphérique avec un rayon interne d'environ 12 nm, sur la surface interne de laquelle sont fixés 180 sites chargés positivement imitant les domaines N-terminaux des protéines. L'eau est décrite par un modèle polarisable à trois billes, et les ions sodium et chlore par des particules MARTINI standard. Les interactions électrostatiques sont calculées par la méthode PPPM avec une permittivité diélectrique de 2.5. Les simulations sont effectuées avec le logiciel LAMMPS à une température de 310 K avec un thermostat de Langevin. Pour valider le modèle, la longueur de persistance de la chaîne libre a été mesurée à différentes concentrations de sel (10–100 mM) et comparée aux données expérimentales pour l'ARN simple brin. Ensuite, le polymère a été placé à l'intérieur de la capside et les profils de densité radiaux et angulaires des monomères ont été analysés après l'atteinte de l'équilibre.
Le modèle a reproduit la dépendance expérimentale de la longueur de persistance de l'ARN simple brin en fonction de la concentration en sel : lorsque le sel augmente de 10 à 100 mM, la longueur de persistance diminue d'environ 2.5 à environ 1.0 nm, ce qui est cohérent avec les données de Chen et al. pour le poly(rU). À l'intérieur de la capside, l'augmentation de la concentration en sel de 10 à 100 mM déplace systématiquement la distribution radiale du génome de la paroi vers le centre, comme le prédisent les théories de champ moyen. Cependant, le temps d'établissement de l'équilibre augmente d'un facteur sept : de 64 ns à 10 mM à 470 ns à 100 mM. La variation indépendante de la charge de la capside et de la densité des sites de charge a montré que la liaison électrostatique totale contrôle l'organisation globale, tandis que les détails géométriques modulent le contact local du génome avec la paroi. Les profils angulaires ont révélé un double pic autour de chaque site de charge pour 180 sites, qui disparaît lorsque le nombre de sites augmente à 1470, lorsque la distance entre eux devient inférieure à la longueur de persistance de la chaîne. La représentation des domaines N-terminaux par une seule bille (saillante) divise le maximum unique en deux, ce qui indique que le génome s'enroule autour de la saillie plutôt que de s'y coller directement.
Ce travail se distingue des études voisines : contrairement aux travaux sur les gènes bactériens ou les transistors ioniques, il teste directement les prédictions des théories de champ moyen pour le repliement d'un polyélectrolyte dans une coque sphérique chargée, en tenant compte explicitement des détails moléculaires. Les travaux voisins utilisent soit des modèles continus, soit un solvant implicite, soit des conditions ioniques fixes. Ici, il est montré pour la première fois que les ions discrets et les sites de charge saillants conduisent à des écarts quantitatifs par rapport au champ moyen, en particulier dans la distribution angulaire du génome. Cela signifie que les approximations continues sont valables pour une image radiale grossière, mais ne capturent pas l'organisation locale, importante pour comprendre la spécificité de l'empaquetage.
Le développement naturel est l'inclusion de la structure secondaire de l'ARN (appariement des bases) et d'une topologie ramifiée, ce qui pourrait modifier le rôle relatif de la charge totale et de la distribution des charges. Il est également important d'étendre le modèle aux ions multivalents (Mg2+ et autres), qui influencent plus fortement le repliement de l'ARN. Les longs domaines N-terminaux flexibles, dépassant une ou deux billes, pourraient ajouter de nouveaux effets conformationnels. Enfin, le passage de l'organisation à l'équilibre dans une capside préformée à la cinétique d'assemblage permettrait une comparaison directe avec les expériences d'auto-assemblage des virus.
Les résultats sont importants pour la virologie (compréhension de l'empaquetage du génome), la physique des polyélectrolytes (comportement des chaînes chargées dans un volume confiné) et les nanotechnologies (conception de capsides pour la délivrance d'acides nucléiques thérapeutiques).
Les prochaines étapes sont l'inclusion de la structure secondaire de l'ARN, la simulation de l'assemblage de la capside autour du génome et le test des effets des ions multivalents.
En physique des polyélectrolytes, la question du rôle exact des ions discrets et des corrélations en géométrie confinée reste ouverte. Ce travail montre que ces effets sont notables déjà avec un sel monovalent et des charges modérées, ce qui pose le défi de développer des théories plus précises allant au-delà du champ moyen. Pour la virologie, le problème non résolu clé est le lien entre la séquence de l'ARN, sa structure secondaire et l'efficacité de l'empaquetage dans la capside ; ce modèle crée une plateforme pour l'étudier.
Cette section s'adresse à l'auteur du travail et à ses collègues. Les travaux voisins ont été trouvés par une recherche sémantique dans notre archive, et non par correspondance de mots. Nous n'évaluons pas : nous montrons le champ autour.
Ce que fait ce travail
Les auteurs ont construit un modèle à gros grains de l'ARN simple brin à l'intérieur de la capside virale, en tenant compte explicitement des ions, de l'eau et de l'électrostatique. Le modèle a reproduit la dépendance expérimentale de la longueur de persistance de l'ARN en fonction de la concentration en sel et a montré que l'augmentation du sel déplace le génome vers l'intérieur de la capside, tandis que le temps d'établissement de l'équilibre est multiplié par sept.
Où est sa force
La force de ce travail réside dans l'inclusion systématique des détails moléculaires (charges discrètes, ions explicites, eau polarisable) dans le modèle électrostatique du compactage du génome, ce qui permet de tester les prédictions des théories de champ moyen au niveau moléculaire.
Vers où avancer
À essayer
Pourquoi cela peut compter
Comprendre comment les conditions ioniques et la discrétisation des charges contrôlent le compactage du génome pourrait aider à développer des stratégies visant à perturber l'assemblage des virus ou à stabiliser les capsides pour des applications biotechnologiques.
Travaux voisins dans notre archive
L'archive est vide autour de ce travail : l'analyse la plus proche n'est qu'à 0.51, contre 0,64 d'ordinaire. Il semble qu'on ait peu foré à côté.
Le temps d'établissement de l'équilibre dans le modèle augmente d'un facteur sept lorsque la concentration en sel passe de 10 à 100 mM, ce qui signifie que dans un environnement plus salé, le génome « met plus de temps à trouver » son repliement final. Pour comparaison, à 10 mM, l'équilibre est atteint en 64 nanosecondes, tandis qu'à 100 mM, il faut près d'une demi-microseconde.