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

Resumen Técnico: Complejación de Polielectrolitos Modulada por Iones de ADN y Ácido Poliacrílico

Planteamiento del Problema
La formación de complejos entre polielectrolitos de la misma carga desafía la intuición electrostática convencional, que predice la repulsión entre macromoléculas con la misma carga. Si bien mecanismos como las correlaciones de contraiones y la heterogeneidad de carga (charge patchiness) son conocidos por mediar la atracción en sistemas específicos (por ejemplo, la condensación del ADN por cationes multivalentes), la asociación de polímeros de carga similar pero químicamente disímiles sigue estando insuficientemente comprendida. Específicamente, la organización atomística de los iones en la interfaz de un biopolímero semirrígido (ADN) y un polímero sintético aniónico y flexible (ácido poliacrílico, PAA) no ha sido examinada sistemáticamente. Este estudio aborda la brecha en la comprensión de cómo la valencia de los contraiones y el tamaño iónico modulan la atracción efectiva y el puenteo entre dos polímeros cargados negativamente en solución acuosa.

Metodología
Los autores emplearon simulaciones de dinámica molecular (MD) atomística para investigar las interacciones ADN–PAA en tres entornos iónicos distintos: NaCl monovalente, y MgCl₂ y CaCl₂ divalentes.

• Construcción del Sistema: Las simulaciones utilizaron un dúplex de ADN de 24 pares de bases (carga neta −46 e⁻) y una cadena única de PAA de 40 monómeros con un 50% de desprotonación (carga neta −20 e⁻). Los sistemas se solvataron en agua TIP3P con una fuerza iónica fisiológica de 0.15 M.
• Campos de Fuerza: El ADN se modeló utilizando el campo de fuerza AMBER parmBSC1, mientras que el PAA utilizó el Campo de Fuerza de Amber Generalizado (GAFF). Los parámetros de los iones se tomaron del conjunto de Joung y Cheatham optimizado para agua TIP3P.
• Protocolo de Simulación: Las corridas de producción se llevaron a cabo en el ensamble NPT a 300 K y 1 bar, con tres réplicas independientes para cada condición de sal con el fin de asegurar la fiabilidad estadística, mientras que los sistemas de solo ADN fueron simulados durante 500 ns y los complejos ADN–PAA durante 1 μs.
• Marco de Análisis:
  • Clasificación de Estados: Un esquema basado en el tiempo definió los estados "unidos" basándose en un corte de distancia de átomos pesados (≤ 0.6 nm) y la presencia de al menos dos enlaces de hidrógeno, requiriendo una persistencia de ≥ 100 ns para distinguir la complejación sostenida de los encuentros transitorios.
  • Métricas Estructurales: Se calcularon la Desviación de la Raíz Cuadrática Media (RMSD), el Radio de Giro (Rg), los momentos principales de inercia y la asfericidad, condicionados a los estados unidos y no unidos.
  • Análisis de Interacción: Se descompusieron las energías de interacción de corto alcance (Lennard-Jones y Coulombicas). Los contactos mediados por iones se cuantificaron identificando tríos [ADN–ion–PAA] dentro de un corte de 0.35 nm. Las funciones de distribución radial (RDF) se computaron para los iones alrededor de los fosfatos del ADN y los carboxilatos del PAA.

Resultados Clave
Las simulaciones revelaron una clara jerarquía de fuerza de asociación y organización estructural dependiente de la valencia iónica y las propiedades de hidratación:

1. Estabilidad del Complejo y Fracción de Unión:

   • CaCl₂: Produjo los complejos más estables, con una fracción media de unión ($f_{bound}$) de 0.89 ± 0.06. El complejo ADN–PAA permaneció asociado durante casi toda la duración de la simulación.
   • MgCl₂: Mostró una estabilidad intermedia ($f_{bound}$ = 0.51 ± 0.21) con desprendimientos más dinámicos e intermitentes.
   • NaCl: Exhibió la mayor variabilidad ($f_{bound}$ = 0.53 ± 0.28) y transiciones frecuentes entre los estados unidos y no unidos, con algunas réplicas mostrando largos periodos de disociación.

2. Coordinación Iónica y Mecanismos de Puenteo:

   • Ca²⁺: Actuó como un reticulador molecular persistente. Las RDF indicaron una fuerte coordinación de esfera interna (0.25–0.30 nm) tanto con los fosfatos del ADN como con los carboxilatos del PAA. El Ca²⁺ formó un promedio de 2.17 puentes por cuadro unido y mantuvo al menos un puente durante aproximadamente el 72% del tiempo de unión.
   • Mg²⁺: Mostró un puenteo más débil y transitorio (promedio de 1.08 puentes; ~44% del tiempo de unión). Su capa de hidratación fuertemente unida favoreció la coordinación de esfera externa, lo que condujo a interacciones más dinámicas asociadas a los surcos en lugar de un bloqueo estable del esqueleto.
   • Na⁺: Proporcionó principalmente el apantallamiento electrostático con una formación de puentes insignificante (~10% del tiempo de unión). La asociación en NaCl fue impulsada por el apantallamiento y contactos débiles de van der Waals en lugar de un puenteo iónico directo.

3. Respuestas Conformacionales:

   • ADN: En presencia de CaCl₂, el ADN unido exhibió un Radio de Giro mayor (Rg ≈ 3.0 nm) en comparación con los segmentos no unidos (Rg ≈ 2.5 nm) y el ADN solo, lo que indica una expansión del armazón del ADN debido al fuerte reticulado y alineación. En contraste, el MgCl₂ y el NaCl mostraron cambios mínimos en la conformación del ADN dependientes del estado.
   • PAA: El PAA experimentó una compactación dentro del complejo, adoptando una conformación más condensada.
   • Orientación: El Ca²⁺ indujo una alineación direccional distinta, donde el eje largo del ADN y el eje corto del PAA adoptaron una configuración estable de ángulo bajo (pico en la distribución de $\cos\theta$). El Mg²⁺ y el Na⁺ mostraron asociaciones más flexibles y menos ordenadas.

4. Modos de Unión:

   • Ca²⁺: Dominado por una asociación persistente al esqueleto, con el puenteo localizado casi exclusivamente en el esqueleto de fosfato.
   • Mg²⁺: Facilitó interacciones transitorias dentro de los surcos mayor y menor, permitiendo que el PAA penetrara la superficie del ADN más profundamente que en los sistemas de Ca²⁺.
   • Na⁺: Los contactos fueron escasos y se limitaron en gran medida al esqueleto sin un puenteo estable.

Significancia y Reivindicaciones
El artículo afirma proporcionar un marco a nivel molecular para comprender cómo las propiedades específicas de los iones (valencia, labilidad de la capa de hidratación y movilidad) gobiernan el autoensamblaje de polielectrolitos de la misma carga. Los autores enfatizan que los iones multivalentes no son intercambiables; el Ca²⁺, con su capa de hidratación lábil y coordinación flexible, actúa como un verdadero reticulador que estabiliza los complejos y reorganiza la topología de los polímeros, mientras que el Mg²⁺ ofrece una estabilización moderada mediante interacciones transitorias, y el Na⁺ depende del apantallamiento.

El estudio destaca que separar los intervalos unidos de los no unidos es esencial para caracterizar con precisión tales sistemas dinámicos, ya que las métricas simples promediadas en el tiempo pueden ocultar comportos fuertes dependientes de la sal. Además, los resultados demuestran que el fuerte reticulado local (como se observa con el Ca²⁺) puede paradójicamente aumentar el tamaño global de la cadena de ADN (expansión) mientras compacta el polímero asociado, desafiando las interpretaciones simples de la complejación basada únicamente en la condensación. Estos hallazgos ofrecen una base para comprender los ensamblajes de nucleoproteínas sensibles a los iones relevantes para la entrega de genes, la biosensórica y el diseño de biomateriales responsivos en entornos iónicos multivalentes.