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

Technische Samenvatting: Ion-gemoduleerde polyelektrolytcomplexatie van DNA en polyacrylzuur

Probleemstelling
De vorming van complexen tussen gelijkgeladen polyelektrolyten daagt de conventionele elektrostatische intuïtie uit, die repulsie voorspelt tussen vergelijkbaar geladen macromoleculen. Hoewel mechanismen zoals tegenion-correlaties en lading-patchiness bekend zijn als mediërende factoren voor attractie in specifieke systemen (bijv. DNA-condensatie door multivalente kationen), blijft de associatie van chemisch ongelijksoortige, gelijkgeladen polymeren onvoldoende begrepen. Specifiek is de atomistische organisatie van ionen aan het grensvlak van een semi-rigide biopolymeer (DNA) en een flexibel anionisch synthetisch polymeer (polyacrylzuur, PAA) nog niet systematisch onderzocht. Deze studie adresseert het gat in het begrip van hoe tegenion-valentie en ionische grootte de effectieve attractie en brugvorming tussen twee negatief geladen polymeren in een waterige oplossing moduleren.

Methodologie
De auteurs gebruikten atomistische moleculaire dynamica (MD) simulaties om de DNA–PAA interacties te onderzoeken in drie verschillende ionische omgevingen: monovalent NaCl, en divalent MgCl₂ en CaCl₂.

• Systeemconstructie: De simulaties maakten gebruik van een 24-basenpaar B-DNA duplex (netto lading −46 e⁻) en een enkel PAA-keten van 40 monomeren met 50% deprotonatie (netto lading −20 e⁻). De systemen werden gesolvateerd in TIP3P water met een fysiologische ionsterkte van 0,15 M.
• Force Fields: Voor DNA werd het AMBER parBSC1 force field gebruikt, terwijl PAA het Generalized Amber Force Field (GAFF) gebruikte. Ionparameters werden genomen uit de Joung en Cheatham set, geoptimaliseerd voor TIP3P water.
• Simulatieprotocol: Productieruns werden uitgevoerd in het NPT-ensemble bij 300 K en 1 bar. Alleen DNA-systemen werden gesimuleerd gedurende 500 ns, terwijl DNA–PAA complexen werden gesimuleerd gedurende 1 μs, met drie onafhankelijke replica's voor elke zoutconditie om statistische betrouwbaarheid te waarborgen.
• Analyseframework:
  • Toestandsclassificatie: Een tijdsopgeloste methode definieerde "gebonden" toestanden op basis van een zware-atoom afstandscut-off (≤ 0,6 nm) en de aanwezigheid van ten minste twee waterstofbruggen, waarbij een persistentie van ≥ 100 ns vereist was om duurzame complexatie te onderscheiden van transiënte ontmoetingen.
  • Structurele metrieken: RMSD, straal van gyratie (Rg), de hoofdmomenten van traagheid en asfericiteit werden conditioneel berekend voor gebonden en ongebonden toestanden.
  • Interactieanalyse: Kortbereik interactie-energieën (Lennard-Jones en Coulombic) werden gedecomposeerd. Ion-gemedieerde contacten werden gekwantificeerd door [DNA–ion–PAA] tripletjes te identificeren binnen een 0,35 nm cut-off. Radiale distributiefuncties (RDFs) werden berekend voor ionen rond DNA-fosfaten en PAA-carboxylaten.

Belangrijkste Resultaten
De simulaties onthulden een duidelijke hiërarchie in associatiekracht en structurele organisatie, afhankelijk van ionvalentie en hydratatie-eigenschappen:

1. Complexstabiliteit en Bindingsfractie:

   • CaCl₂: Produceerde de meest stabiele complexen, met een gemiddelde bindingsfractie ($f_{bound}$) van 0,89 ± 0,06. Het DNA–PAA complex bleef gedurende bijna de gehele simulatieduur geassocieerd.
   • MgCl₂: Toonde een intermediaire stabiliteit ($f_{bound}$ = 0,51 ± 0,21) met meer dynamische, intermitterende loskoppelingen.
   • NaCl: Vertoonde de hoogste variabiliteit ($f_{bound}$ = 0,53 ± 0,28) en frequente transities tussen gebonden en ongebonden toestanden, waarbij sommige replica's lange perioden van dissociatie vertoonden.

2. Ioncoördinatie en Brugvormingsmechanismen:

   • Ca²⁺: Fungeerde als een persistente moleculaire cross-linker. RDF's indiceerden sterke inner-sphere coördinatie (0,25–0,30 nm) met zowel DNA-fosfaten als PAA-carboxylaten. Ca²⁺ vormde gemiddeld 2,17 bruggen per gebonden frame en behield ten minste één brug tijdens ~72% van de gebonden tijd.
   • Mg²⁺: Toonde zwakkere, transiënte brugvorming (gemiddeld 1,08 bruggen; ~44% van de gebonden tijd). De strak gebonden hydratatieschil bevoordeelde outer-sphere coördinatie, wat leidde tot meer dynamische, aan de groef geassocieerde interacties in plaats van stabiele ruggengraat-verankering.
   • Na⁺: Zorgde primair voor elektrostatische afscherming met verwaarloosbare brugvorming (~10% van de gebonden tijd). Associatie in NaCl werd gedreven door afscherming en zwakke van der Waals-contacten in plaats van directe ion-brugvorming.

3. Conformationele Respons:

   • DNA: Gebonden DNA in CaCl₂ vertoonde een grotere straal van gyratie (Rg ≈ 3,0 nm) vergeleken met ongebonden segmenten (Rg ≈ 2,5 nm) en DNA alleen, wat wijst op een expansie van het DNA-framework door sterke cross-linking en uitlijning. In contrast hiermee vertoonden MgCl₂ en NaCl minimale toestand-afhankelijke veranderingen in DNA-conformatie.
   • PAA: PAA onderging compactie binnen het complex en nam een meer gecondenseerde conformatie aan.
   • Oriëntatie: Ca²⁺ induceerde een duidelijke directionele uitlijning, waarbij de lange as van DNA en de korte as van PAA een stabiele, laag-hoekige configuratie aannamen (piek in de $\cos\theta$-distributie). Mg²⁺ en Na⁺ vertoonden flexibelere, minder geordende associaties.

4. Bindingsmodi:

   • Ca²⁺: Werd gedomineerd door persistente ruggengraat-associatie, waarbij brugvorming bijna uitsluitend gelokaliseerd was bij de fosfaat-ruggengraat.
   • Mg²⁺: Faciliteerde transiënte interacties binnen de major en minor grooves, waardoor PAA dieper in het DNA-oppervlak kon penetreren dan in Ca²⁺-systemen.
   • Na⁺: Contacten waren schaars en grotendeels beperkt tot de ruggengraat zonder stabiele brugvorming.

Betekenis en Claims
Het artikel claimt een moleculair kader te bieden voor het begrijpen van hoe ion-specifieke eigenschappen (valentie, de labiliteit van de hydratatieschil en mobiliteit) de zelfassemblage van gelijkgeladen polyelektrolyten sturen. De auteurs benadrukken dat multivalente ionen niet uitwisselbaar zijn; Ca²⁺, met zijn labiele hydratatieschil en flexibele coördinatie, fungeert als een echte cross-linker die complexen stabiliseert en de polymeer-topologie reorganiseert, terwijl Mg²⁺ een matige stabilisatie biedt via transiënte interacties, en Na⁺ vertrouwt op afscherming.

De studie benadrukt dat het scheiden van gebonden en ongebonden intervallen essentieel is voor het accuraat karakteriseren van dergelijke dynamische systemen, aangezien eenvoudige tijdgemiddelde metrieken sterke zoutafhankelijke gedragingen kunnen maskeren. Bovendien laten de resultaten zien dat sterke lokale cross-linking (zoals gezien met Ca²⁺) paradoxaal genoeg de globale grootte van de DNA-keten kan vergroten (expansie) terwijl het geassocieerde polymeer compacter wordt, wat de eenvoudige interpretaties van complexatie als louter condensatie uitdaagt. Deze bevindingen bieden een basis voor het begrijpen van ion-responsieve nucleïnezuur–polymeer assemblages die relevant zijn voor genlevering, biosensing en het ontwerp van responsieve biomaterialen in multivalente ionische omgevingen.