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

Sintesi Tecnica: Complessazione di Polielettroliti Modulata dagli Ioni tra DNA e Acido Poliacrilico

Definizione del Problema
La formazione di complessi tra polielettroliti con la stessa carica sfida l'intuizione elettrostatica convenzionale, che prevede la repulsione tra macromolecole con carica simile. Sebbene meccanismi come le correlazioni tra controioni e la "patchiness" (eterogeneità di carica) siano noti per mediare l'attrazione in sistemi specifici (ad esempio, la condensazione del DNA da parte di cationi multivalenti), l'associazione di polimeri con cariche uguali ma chimicamente dissimili rimane scarsamente compresa. Nello specifico, l'organizzazione atomistica degli ioni all'interfaccia tra un biopolimero semirigido (DNA) e un polimero sintetico anionico flessibile (acido poliacrilico, PAA) non è stata esaminata sistematicamente. Questo studio affronta il vuoto nella comprensione di come la valenza dei controioni e la dimensione ionica modulino l'attrazione effettiva e il "bridging" (formazione di ponti) tra due polimeri carichi negativamente in soluzione acquosa.

Metodologia
Gli autori hanno impiegato simulazioni di dinamica molecolare (MD) atomistica per investigare le interazioni DNA–PAA in tre distinti ambienti ionici: NaCl monovalente, e MgCl₂ e CaCl₂ bivalenti.

• Costruzione del Sistema: Le simulazioni hanno utilizzato un duplex di DNA B a 24 coppie di basi (carica netta −46 e⁻) e una singola catena di PAA di 40 monomeri con deprotonazione al 50% (carica netta −20 e⁻). I sistemi sono stati solvatati in acqua TIP3P con una forza ionica fisiologica di 0,15 M.
• Campi di Forza: Il DNA è stato modellato utilizzando il campo di forza AMBER parmBSC1, mentre il PAA ha utilizzato il Generalized Amber Force Field (GAFF). I parametri degli ioni sono stati presi dal set Joung e Cheatham, ottimizzato per l'acqua TIP3P.
• Protocollo di Simulazione: Le corse di produzione sono state condotte nell'ensemble NPT a 300 K e 1 bar, con simulazioni DNA–PAA di 1 μs e simulazioni dei soli sistemi DNA di 500 ns, con tre repliche indipendenti per ogni condizione salina per garantire l'affidabilità statistica.
• Framework di Analisi:
  • Classificazione degli Stati: Uno schema risolto nel tempo ha definito gli stati "legati" basandosi su un cutoff della distanza tra atomi pesanti (≤ 0,6 nm) e la presenza di almeno due legami a idrogeno, richiedendo una persistenza ≥ 100 ns per distinguere la complessazione sostenuta dagli incontri transitori.
  • Metriche Strutturali: RMSD, Raggio di Giraizazione (Rg), momenti principali di inerzia e asfericità sono stati calcolati condizionatamente sugli stati legati e non legati.
  • Analisi delle Interazioni: Le energie di interazione a corto raggio (Lennard-Jones e Coulombiche) sono state decomposte. I contatti mediati dagli ioni sono stati quantificati identificando triplette [DNA–ione–PAA] entro un cutoff di 0,35 nm. Le funzioni di distribuzione radiale (RDF) sono state calcolate per gli ioni attorno ai fosfati del DNA e ai carbossilati del PAA.

Risultati Chiave
Le simulazioni hanno rivelato una chiara gerarchia di forza di associazione e organizzazione strutturale dipendente dalla valenza ionica e dalle proprietà di idratazione:

1. Stabilità del Complesso e Frazione di Legame:

   • CaCl₂: Ha prodotto i complessi più stabili, con una frazione media legata ($f_{bound}$) di 0,89 ± 0,06. Il complesso DNA–PAA è rimasto associato per quasi tutta la durata della simulazione.
   • MgCl₂: Ha mostrato una stabilità intermedia ($f_{bound}$ = 0,51 ± 0,21) con distacchi più dinamici e intermittenti.
   • NaCl: Ha esibito la maggiore variabilità ($f_{bound}$ = 0,53 ± 0,28) e frequenti transizioni tra stati legati e non legati, con alcune repliche che hanno mostrato lunghi periodi di dissociazione.

2. Coordinazione Ionica e Meccanismi di Bridging:

   • Ca²⁺: Ha agito come un cross-linker molecolare persistente. Le RDF hanno indicato una forte coordinazione a sfera interna (0,25–0,30 nm) sia con i fosfati del DNA che con i carbossilati del PAA. Il Ca²⁺ ha formato una media di 2,17 ponti per frame legato e ha mantenuto almeno un ponte durante circa il 72% del tempo di legame.
   • Mg²⁺: Ha mostrato un bridging più debole e transitorio (media di 1,08 ponti; ~44% del tempo di legame). Il suo guscio di idratazione strettamente legato ha favorito la coordinazione a sfera esterna, portando a interazioni più dinamiche associate ai solchi piuttosto che a un bloccaggio stabile dello scheletro.
   • Na⁺: Ha fornito principalmente uno screening elettrostatico con un trascurabile formazione di ponti (~10% del tempo di legame). L'associazione in NaCl è stata guidata dallo screening e da deboli contatti di van der Waals piuttosto che dal diretto bridging ionico.

3. Risposte Conformazionali:

   • DNA: In CaCl₂, il DNA legato ha esibito un Raggio di Giraizazione (Rg ≈ 3,0 nm) maggiore rispetto ai segmenti non legati (Rg ≈ 2,5 nm) e al DNA da solo, indicando un'espansione del framework del DNA dovuta al forte cross-linking e all'allineamento. Al contrario, MgCl₂ e NaCl hanno mostrato cambiamenti minimi della conformazione del DNA in base allo stato.
   • PAA: In CaCl₂, il PAA ha subito una compattazione all'interno del complesso, adottando una conformazione più condensata.
   • Orientamento: Il Ca²⁺ ha indotto un allineamento direzionale distinto, in cui l'asse lungo del DNA e l'asse corto del PAA hanno adottato una configurazione stabile a basso angolo (picco nella distribuzione $\cos\theta$). Mg²⁺ e Na⁺ hanno mostrato associazioni più flessibili e meno ordinate.

4. Modalità di Legame:

   • Ca²⁺: Dominato da una persistente associazione allo scheletro, con il bridging localizzato quasi esclusivamente sul backbone fosfato.
   • Mg²⁺: Ha facilitato interazioni transitorie all'interno dei solchi maggiore e minore, permettendo al PAA di penetrare la superficie del DNA più profondamente rispetto ai sistemi con Ca²⁺.
   • Na⁺: I contatti sono stati sparsi e ampiamente confinati allo scheletro senza un bridging stabile.

Significato e Rivendicazioni
L'articolo sostiene di fornire un quadro a livello molecolare per comprendere come le proprietà specifiche degli ioni (valenza, labilità del guscio di idratazione e mobilità) governino l'auto-assemblaggio di polielettroliti con la stessa carica. Gli autori sottolineano che gli ioni multivalenti non sono intercambiabili: il Ca²⁺, con il suo guscio di idratazione labile e la sua coordinazione flessibile, agisce come un vero cross-linker che stabilizza i complessi e riorganizza la topologia dei polimeri, mentre il Mg²⁺ offre una stabilizzazione moderata attraverso interazioni transitorie, e il Na⁺ si affida allo screening.

Lo studio evidenzia che separare gli intervalli legati da quelli non legati è essenziale per caratterizzare accuratamente tali sistemi dinamici, poiché le semplici metriche mediate nel tempo possono oscurare comportamenti fortemente dipendenti dai sali. Inoltre, i risultati dimostrano che il forte cross-linking locale (come visto in CaCl₂) può paradossalmente aumentare la dimensione globale della catena di DNA (espansione) pur compattando il polimero associato, sfidando le semplici interpretazioni della complessazione intesa solo come condensazione. Questi risultati offrono una base per comprendere gli assemblaggi acido nucleico-polimero responsivi agli ioni, rilevanti per la consegna di geni, il biosensing e la progettazione di biomateriali responsivi in ambienti ionici multivalenti.