🔬 condensed matter

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

Atomistic molecular dynamics simulations reveal that ion valency critically modulates the stability and structure of DNA-polyacrylic acid complexes, with Ca $^{2+}$ driving strong, persistent inner-sphere bridging, Mg $^{2+}$ facilitating weaker transient interactions, and Na $^+$ providing only electrostatic screening.

Original authors: Sisem Ektirici, Vagelis Harmandaris, Christos N. Likos, Terpsichori S. Alexiou

Published 2026-01-29

📖 4 min read☕ Coffee break read

CC BY 4.0

Original authors: Sisem Ektirici, Vagelis Harmandaris, Christos N. Likos, Terpsichori S. Alexiou

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine two very stubborn, negatively charged magnets. In the world of physics, these two magnets should push each other away violently, never touching. This is exactly what happens with DNA (the molecule that holds our genetic code) and Polyacrylic Acid (PAA) (a common synthetic polymer). Both are covered in negative electrical charges, so standard logic says they should repel each other.

However, this paper asks a fascinating question: What happens if we throw different types of "glue" (ions) into the mix?

The researchers used powerful computer simulations to watch how these two negative molecules behave when surrounded by three different types of salt solutions:

Sodium (Na⁺): Like the common salt in your kitchen.
Magnesium (Mg²⁺): A stronger, double-charged ion.
Calcium (Ca²⁺): Another double-charged ion, but slightly larger and more flexible.

Here is what they discovered, explained through simple analogies:

1. The Sodium Salt (NaCl): The "Weak Shield"

When the researchers added Sodium, the two negative molecules mostly stayed apart.

The Analogy: Imagine two people wearing heavy, static-charged wool coats trying to hug. Sodium acts like a thin, flimsy blanket thrown over them. It reduces the static shock a little bit, allowing them to get close occasionally, but they can't stick together for long. They bump into each other, drift apart, and bump again.
The Result: The DNA and PAA formed a "complex" (stuck together) about 50% of the time, but it was a shaky, unstable relationship. They were mostly just floating near each other without a strong bond.

2. The Magnesium Salt (MgCl₂): The "Fickle Matchmaker"

Magnesium is a double-charged ion, so it's a stronger glue than Sodium.

The Analogy: Think of Magnesium as a matchmaker who is very eager but a bit clumsy. It tries to hold the two negative molecules together, but it holds on too tightly to the water around it (its "hydration shell"). It can't quite let go of the water to grab the molecules directly.
The Result: The molecules got closer than with Sodium, but the connection was still temporary. The Magnesium ions acted like a bridge that kept collapsing and rebuilding. The DNA and PAA would stick together for a while, then drift apart, then stick again. It was a "transient" relationship—stronger than Sodium, but not permanent.

3. The Calcium Salt (CaCl₂): The "Super Glue"

Calcium was the game-changer.

The Analogy: Calcium is like a flexible, super-strong clamp. Unlike Magnesium, it doesn't hold onto its water tightly. It can easily let go of the water and grab directly onto the DNA and the PAA at the same time. It acts as a literal bridge, locking the two negative chains together.
The Result: Once Calcium grabbed them, they stayed stuck. The complex was stable and persistent. In the simulations, the DNA and PAA remained bound together for nearly the entire time they were watched.

How the Shapes Changed

The paper also looked at how the shapes of these molecules changed when they got stuck together:

With Calcium: The DNA stretched out and became longer, while the PAA curled up and became tighter. It's as if the Calcium clamp pulled the DNA straight like a taut rope while squeezing the PAA into a tight ball. They lined up neatly, like a ruler and a coiled spring pressed together.
With Magnesium: The DNA stayed mostly the same shape, but the PAA would occasionally poke into the "grooves" (the little valleys) of the DNA helix. It was a looser, more chaotic dance.
With Sodium: The molecules barely changed their shapes because they weren't really holding on to each other.

The Big Takeaway

The main lesson from this study is that not all salts are created equal.

Even though DNA and PAA are both negatively charged and should repel each other, the specific type of ion in the water decides their fate:

Sodium just lets them float near each other.
Magnesium lets them hug briefly.
Calcium acts as a molecular cross-linker, permanently welding them together and forcing them into a specific, organized shape.

This research helps us understand the invisible rules of how charged molecules organize themselves in water, showing that the "personality" of the ion (its size and how it holds water) is just as important as the electrical charge itself.

Technical Summary: Ion-Modulated Polyelectrolyte Complexation of DNA and Polyacrylic Acid

Problem Statement
The formation of complexes between like-charged polyelectrolytes challenges conventional electrostatic intuition, which predicts repulsion between similarly charged macromolecules. While mechanisms such as counterion correlations and charge patchiness are known to mediate attraction in specific systems (e.g., DNA condensation by multivalent cations), the association of chemically dissimilar, like-charged polymers remains poorly understood. Specifically, the atomistic organization of ions at the interface of a semi-rigid biopolymer (DNA) and a flexible anionic synthetic polymer (polyacrylic acid, PAA) has not been systematically examined. This study addresses the gap in understanding how counterion valency and ionic size modulate the effective attraction and bridging between two negatively charged polymers in aqueous solution.

Methodology
The authors employed atomistic molecular dynamics (MD) simulations to investigate DNA–PAA interactions in three distinct ionic environments: monovalent NaCl, and divalent MgCl₂ and CaCl₂.

System Construction: Simulations utilized a 24-base-pair B-DNA duplex (net charge −46 e⁻) and a single PAA chain of 40 monomers with 50% deprotonation (net charge −20 e⁻). Systems were solvated in TIP3P water with a physiological ionic strength of 0.15 M.
Force Fields: DNA was modeled using the AMBER parmBSC1 force field, while PAA utilized the Generalized Amber Force Field (GAFF). Ion parameters were taken from the Joung and Cheatham set optimized for TIP3P water.
Simulation Protocol: Production runs were conducted in the NPT ensemble at 300 K and 1 bar. DNA-only systems were simulated for 500 ns, while DNA–PAA complexes were simulated for 1 μs, with three independent replicas performed for each salt condition to ensure statistical reliability.
Analysis Framework:
- State Classification: A time-resolved scheme defined "bound" states based on a heavy-atom distance cutoff (≤ 0.6 nm) and the presence of at least two hydrogen bonds, requiring a persistence of ≥ 100 ns to distinguish sustained complexation from transient encounters.
- Structural Metrics: Root Mean Square Deviation (RMSD), Radius of Gyration (Rg), principal moments of inertia, and asphericity were calculated conditionally on bound and unbound states.
- Interaction Analysis: Short-range interaction energies (Lennard-Jones and Coulombic) were decomposed. Ion-mediated contacts were quantified by identifying [DNA–ion–PAA] triplets within a 0.35 nm cutoff. Radial distribution functions (RDFs) were computed for ions around DNA phosphates and PAA carboxylates.

Key Results
The simulations revealed a clear hierarchy of association strength and structural organization dependent on ion valency and hydration properties:

Complex Stability and Binding Fraction:
- CaCl₂: Produced the most stable complexes, with a mean bound fraction ( $f_{bound}$ ) of 0.89 ± 0.06. The DNA–PAA complex remained associated for nearly the entire simulation duration.
- MgCl₂: Showed intermediate stability ( $f_{bound}$ = 0.51 ± 0.21) with more dynamic, intermittent detachment.
- NaCl: Exhibited the highest variability ( $f_{bound}$ = 0.53 ± 0.28) and frequent transitions between bound and unbound states, with some replicas showing long periods of dissociation.
Ion Coordination and Bridging Mechanisms:
- Ca²⁺: Acted as a persistent molecular cross-linker. RDFs indicated strong inner-sphere coordination (0.25–0.30 nm) with both DNA phosphates and PAA carboxylates. Ca²⁺ formed an average of 2.17 bridges per bound frame and maintained at least one bridge during ~72% of the bound time.
- Mg²⁺: Displayed weaker, transient bridging (average 1.08 bridges; ~44% of bound time). Its tightly bound hydration shell favored outer-sphere coordination, leading to more dynamic, groove-associated interactions rather than stable backbone locking.
- Na⁺: Primarily provided electrostatic screening with negligible bridge formation (~10% of bound time). Association in NaCl was driven by screening and weak van der Waals contacts rather than direct ion bridging.
Conformational Responses:
- DNA: In CaCl₂, bound DNA exhibited a larger Radius of Gyration (Rg ≈ 3.0 nm) compared to unbound segments (Rg ≈ 2.5 nm) and DNA alone, indicating an expansion of the DNA framework due to strong cross-linking and alignment. In contrast, MgCl₂ and NaCl showed minimal state-dependent changes in DNA conformation.
- PAA: In CaCl₂, PAA underwent compaction within the complex, adopting a more condensed conformation.
- Orientation: Ca²⁺ induced a distinct directional alignment, where the long axis of DNA and the short axis of PAA adopted a stable, low-angle configuration (peak in $\cos\theta$ distribution). Mg²⁺ and Na⁺ showed more flexible, less ordered associations.
Binding Modes:
- Ca²⁺: Dominated by persistent backbone association, with bridging localized almost exclusively at the phosphate backbone.
- Mg²⁺: Facilitated transient interactions within the major and minor grooves, allowing PAA to penetrate the DNA surface more deeply than in Ca²⁺ systems.
- Na⁺: Contacts were sparse and largely confined to the backbone without stable bridging.

Significance and Claims
The paper claims to provide a molecular-level framework for understanding how ion-specific properties (valency, hydration shell lability, and mobility) govern the self-assembly of like-charged polyelectrolytes. The authors emphasize that multivalent ions are not interchangeable; Ca²⁺, with its labile hydration shell and flexible coordination, acts as a true cross-linker that stabilizes complexes and reorganizes polymer topology, whereas Mg²⁺ offers moderate stabilization through transient interactions, and Na⁺ relies on screening.

The study highlights that separating bound and unbound intervals is essential for accurately characterizing such dynamic systems, as simple time-averaged metrics can obscure strong salt-dependent behaviors. Furthermore, the results demonstrate that strong local cross-linking (as seen with Ca²⁺) can paradoxically increase the global size of the DNA chain (expansion) while compacting the associated polymer, challenging simple interpretations of complexation solely as condensation. These findings offer a basis for understanding ion-responsive nucleic acid–polymer assemblies relevant to gene delivery, biosensing, and the design of responsive biomaterials in multivalent ionic environments.

1. The Sodium Salt (NaCl): The "Weak Shield"

2. The Magnesium Salt (MgCl₂): The "Fickle Matchmaker"

3. The Calcium Salt (CaCl₂): The "Super Glue"

How the Shapes Changed

The Big Takeaway

More like this