Ion-Modulated Polyelectrolyte Complexation of DNA and Polyacrylic Acid from Molecular Dynamics Simulations

Résumé technique : Complexation de polyélectrolytes modulée par les ions entre l'ADN et l'acide polyacrylique

Énoncé du problème
La formation de complexes entre polyélectrolytes de même charge défie l'intuition électrostatique conventionnelle, qui prédit une répulsion entre macromolécules de même charge. Bien que des mécanismes tels que les corrélations de contre-ions et la présence de domaines de charge (charge patchiness) soient connus pour médier l'attraction dans des systèmes spécifiques (par exemple, la condensation de l'ADN par des cations multivalents), l'association de polymères chimiquement dissemblables et de même charge reste mal comprise. Plus précisément, l'organisation atomistique des ions à l'interface d'un biopolymère semi-rigide (ADN) et d'un polymère anionique flexible (acide polyacrylique, PAA) n'a pas été examinée de manière systématique. Cette étude aborde le manque de compréhension concernant la façon dont la valence des contre-ions et la taille ionique modulent l'attraction effective et le pontage entre deux polymères chargés négativement en solution aqueuse.

Méthodologie
Les auteurs ont employé des simulations de dynamique moléculaire (MD) atomistiques pour étudier les interactions ADN–PAA dans trois environnements ioniques distincts : NaCl monovalent, et MgCl₂ et CaCl₂ divalents.

• Construction du système : Les simulations ont utilisé un duplex d'ADN de 24 paires de bases (charge nette −46 e⁻) et une chaîne unique de PAA de 40 monomères avec une déprotonation de 50 % (charge nette −20 e⁻). Les systèmes ont été solvatés dans de l'eau TIP3P à une force ionique physiologique de 0,15 M.
• Champs de force : L'ADN a été modélisé avec le champ de force AMBER parmBSC1, tandis que le PAA a utilisé le Generalized Amber Force Field (GAFF). Les paramètres des ions ont été tirés de l'ensemble de Joung et Cheatham optimisé pour l'eau TIP3P.
• Protocole de simulation : Les cycles de production ont été menés dans l'ensemble NPT à 300 K et 1 bar, avec des simulations de l'ADN seul de 500 ns, tandis que les complexes ADN–PAA ont été simulés pendant 1 μs, avec trois répliques indépendantes pour chaque condition de sel afin d'assurer la fiabilité statistique.
• Cadre d'analyse :
  • Classification des états : Un schéma résolu dans le temps a défini les états « liés » sur la base d'une coupure de distance entre atomes lourds (≤ 0,6 nm) et de la présence d'au moins deux liaisons hydrogène, nécessitant une persistance ≥ 100 ns pour distinguer une complexation soutenue d'une rencontre transitoire.
  • Métriques structurelles : La déviation quadratique moyenne (RMSD), le rayon de giration (Rg), les moments principaux d'inertie et l'asphéricité ont été calculés de manière conditionnelle sur les états liés et non liés.
  • Analyse des interactions : Les énergies d'interaction à courte portée (Lennard-Jones et Coulombiques) ont été décomposées. Les contacts médiés par les ions ont été quantifiés en identifiant des triplets [ADN–ion–PAA] dans une coupure de 0,35 nm. Les fonctions de distribution radiale (RDF) ont été calculées pour les ions autour des phosphates de l'ADN et des carboxylates du PAA.

Résultats clés
Les simulations ont révélé une hiérarchie claire de la force d'association et de l'organisation structurelle dépendant de la valence des ions et des propriétés d'hydratation :

1. Stabilité du complexe et fraction de liaison :

   • CaCl₂ : A produit les complexes les plus stables, avec une fraction liée moyenne ($f_{bound}$) de 0,89 ± 0,06. Le complexe ADN–PAA est resté associé pendant la quasi-totalité de la durée de la simulation.
   • MgCl₂ : A montré une stabilité intermédiaire ($f_{bound}$ = 0,51 ± 0,21) avec des détachements plus dynamiques et intermittents.
   • NaCl : A présenté la plus grande variabilité ($f_{bound}$ = 0,53 ± 0,28) et des transitions fréquentes entre les états liés et non liés, certaines répliques montrant de longues périodes de dissociation.

2. Coordination ionique et mécanismes de pontage :

   • Ca²⁺ : A agi comme un agent de réticulation moléculaire persistant. Les RDF ont indiqué une coordination de sphère interne forte (0,25–0,30 nm) avec les phosphates de l'ADN et les carboxylates du PAA. Le Ca²⁺ a formé une moyenne de 2,17 ponts par image liée et a maintenu au moins un pont pendant environ 72 % du temps lié.
   • Mg²⁺ : A présenté un pontage plus faible et transitoire (moyenne de 1,08 ponts ; ~44 % du temps lié). Sa couche d'hydratation étroitement liée favorisait une coordination de sphère externe, conduisant à des interactions plus dynamiques associées aux sillons plutôt qu'à un verrouillage stable du squelette.
   • Na⁺ : A principalement assuré un écrantage électrostatique avec une formation de pont négligeable (~10 % du temps lié). L'association dans le NaCl était pilotée par l'écrantage et de faibles contacts de van der Waals plutôt que par un pontage ionique direct.

3. Réponses conformationnelles :

   • ADN : Dans le CaCl₂, l'ADN lié présentait un rayon de giration plus grand (Rg ≈ 3,0 nm) par rapport aux segments non liés (Rg ≈ 2,5 nm) et à l'ADN seul, indiquant une expansion de la structure de l'ADN due à un fort réticulation et un alignement. En revanche, le MgCl₂ et le NaCl ont montré des changements de conformation de l'ADN minimaux selon l'état.
   • PAA : Dans le CaCl₂, le PAA a subi une compaction au sein du complexe, adoptant une conformation plus condensée.
   • Orientation : Le Ca²⁺ a induit un alignement directionnel distinct, où le grand axe de l'ADN et le petit axe du PAA ont adopté une configuration stable à faible angle (pic dans la distribution $\cos\theta$). Le Mg²⁺ et le Na⁺ ont montré des associations plus flexibles et moins ordonnées.

4. Modes de liaison :

   • Ca²⁺ : Dominé par une association persistante au squelette, le pontage étant localisé presque exclusivement au niveau du squelette phosphate.
   • Mg²⁺ : A facilité des interactions transitoires dans les sillons majeur et mineur, permettant au PAA de pénétrer plus profondément la surface de l'ADN que dans les systèmes Ca²⁺.
   • Na⁺ : Les contacts étaient rares et largement confinés au squelette sans pontage stable.

Signification et affirmations
L'article affirme fournir un cadre moléculaire pour comprendre comment les propriétés spécifiques des ions (valence, labilité de la couche d'hydratation, mobilité) régissent l'auto-assemblage de polyélectrolytes de même charge. Les auteurs soulignent que les ions multivalents ne sont pas interchangeables ; le Ca²⁺, avec sa couche d'hydratation labile et sa coordination flexible, agit comme un véritable agent de réticulation qui stabilise les complexes et réorganise la topologie des polymères, tandis que le Mg²⁺ offre une stabilisation modérée par des interactions transitoires, et le Na⁺ repose sur l'écrantage.

L'étude souligne que la séparation des intervalles liés et non liés est essentielle pour caractériser avec précision ces systèmes dynamiques, car les simples métriques moyennées dans le temps peuvent masquer des comportements fortement dépendants du sel. De plus, les résultats démontrent qu'un fort réticulation locale (comme observé avec le Ca²⁺) peut paradoxalement augmenter la taille globale de la chaîne d'ADN (expansion) tout en compactant le polymère associé, défiant les interprétations simples de la complexation uniquement comme une condensation. Ces conclusions offrent une base pour comprendre les assemblages de nucléotides et de polymères sensibles aux ions, pertinents pour la délivrance de gènes, la biodétection et la conception de biomatériaux réactifs dans des environnements ioniques multivalents.