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

Resumo Técnico: Complexação de Polieletrólitos Modulada por Íons de DNA e Ácido Poliacrílico

Definição do Problema
A formação de complexos entre polieletrólitos de mesma carga desafia a intuição eletrostática convencional, que prevê a repulsão entre macromoléculas de cargas semelhantes. Embora mecanismos como correlações de contraíons e heterogeneidade de carga (charge patchiness) sejam conhecidos por mediar a atração em sistemas específicos (ex: condensação de DNA por cátions multivalentes), a associação de polímeros quimicamente dissimilares e de mesma carga permanece pouco compreendida. Especificamente, a organização atomística de íons na interface entre um biopolímero semirrígido (DNA) e um polímero aniônico flexível (ácido poliacrílico, PAA) não foi examinada sistematicamente. Este estudo aborda a lacuna na compreensão de como a valência e o tamanho iônico modulam a atração efetiva e a ponte (bridging) entre dois polímeros carregados negativamente em solução aquosa.

Metodologia
Os autores empregaram simulações de dinâmica molecular (MD) atomística para investigar as interações DNA–PAA em três ambientes iônicos distintos: NaCl monovalente, e MgCl₂ e CaCl₂ divalentes.

• Construção do Sistema: As simulações utilizaram um duplex de DNA B de 24 pares de bases (carga líquida −46 e⁻) e uma única cadeia de PAA de 40 monômeros com 50% de desprotonação (carga líquida −20 e⁻). Os sistemas foram solvatados em água TIP3P com uma força iônica fisiológica de 0,15 M.
• Campos de Força: O DNA foi modelado utilizando o campo de força AMBER parmBSC1, enquanto o PAA utilizou o Generalized Amber Force Field (GAFF). Os parâmetros dos íons foram retirados do conjunto de Joung e Cheatham, otimizado para água TIP3P.
• Protocolo de Simulação: As corridas de produção foram conduzidas no ensemble NPT a 300 K e 1 bar, com simulações de DNA isolado por 500 ns, enquanto os complexos DNA–PAA foram simulados por 1 μs, com três réplicas independentes realizadas para cada condição de sal para garantir a confiabilidade estatística.
• Estrutura de Análise:
  • Classificação de Estados: Um esquema baseado em tempo definiu estados "ligados" baseando-se em um corte de distância de átomos pesados (≤ 0,6 nm) e na presença de pelo menos duas ligações de hidrogênio, exigindo uma persistência de ≥ 100 ns para distinguir a complexação sustentada de encontros transitórios.
  • Métricas Estruturais: O Desvio Quadrático Médio (RMSD), Raio de Giro (Rg), momentos principais de inércia e assfericidade foram calculados condicionalmente sobre os estados ligados e não ligados.
  • Análise de Interação: As energias de interação de curto alcance (Lennard-Jones e Coulomb) foram decompostas. Os contatos mediados por íons foram quantificados pela identificação de tripletos [DNA–íon–PAA] dentro de um corte de 0,35 nm. Funções de distribuição radial (RDFs) foram computadas para íons ao redor dos fosfatos de DNA e carboxilatos de PAA.

Resultados Principais
As simulações revelaram uma clara hierarquia de força de associação e organização estrutural dependente da valência iônica e das propriedades de hidratação:

1. Estabilidade do Complexo e Fração de Ligação:

   • CaCl₂: Produziu os complexos mais estáveis, com uma fração média ligada ($f_{bound}$) de 0,89 ± 0,06. O complexo DNA–PAA permaneceu associado durante quase toda a duração da simulação.
   • MgCl₂: Mostrou estabilidade intermediária ($f_{bound}$ = 0,51 ± 0,21) com desprendimentos mais dinâmicos e intermitentes.
   • NaCl: Exibiu a maior variabilidade ($f_{bound}$ = 0,53 ± 0,28) e transições frequentes entre estados ligados e não ligados, com algumas réplicas mostrando longos períodos de dissociação.

2. Coordenação Iônica e Mecanismos de Ponte (Bridging):

   • Ca²⁺: Atuou como um reticulador (cross-linker) molecular persistente. As RDFs indicaram coordenação de esfera interna forte (0,25–0,30 nm) tanto com fosfatos de DNA quanto com carboxilatos de PAA. O Ca²⁺ formou uma média de 2,17 pontes por quadro ligado e manteve pelo menos uma ponte durante ~72% do tempo de ligação.
   • Mg²⁺: Exibiu uma ponte mais fraca e transitória (média de 1,08 pontes; ~44% do tempo ligado). Sua camada de hidratação fortemente ligada favoreceu a coordenação de esfera externa, levando a interações mais dinâmicas associadas aos sulcos (grooves) em vez de um travamento estável do esqueleto (backbone).
   • Na⁺: Forneceu principalmente a blindagem eletrostática (screening) com negligível formação de pontes (~10% do tempo ligado). A associação em NaCl foi impulsionada pela blindagem e contatos de Van der Waals fracos, em vez de ponte direta por íons.

3. Respostas Conformacionais:

   • DNA: Em CaCl₂, o DNA ligado exibiu um Raio de Giro (Rg ≈ 3,0 nm) maior em comparação com segmentos não ligados (Rg ≈ 2,5 nm) e o DNA isolado, indicando uma expansão da estrutura do DNA devido à forte reticulação e alinhamento. Em contraste, MgCl₂ e NaCl mostraram mudanças mínimas de conformação dependentes do estado.
   • PAA: O PAA sofreu compactação no complexo, adotando uma conformação mais condensada.
   • Orientação: O Ca²⁺ induziu um alinhamento direcional distinto, onde o eixo longo do DNA e o eixo curto do PAA adotaram uma configuração estável de baixo ângulo (pico na distribuição de $\cos\theta$). Mg²⁺ e Na⁺ mostraram associações mais flexíveis e menos ordenadas.

4. Modos de Ligação:

   • Ca²⁺: Dominado pela associação persistente ao esqueleto, com a ponte localizada quase exclusivamente no esqueleto de fosfato.
   • Mg²⁺: Facilitou interações transitórias dentro dos sulcos maior e menor, permitindo que o PAA penetrasse mais profundamente na superfície do DNA do que nos sistemas de Ca²⁺.
   • Na⁺: Os contatos foram esparsos e amplamente confinados ao esqueleto, sem formação de pontes estáveis.

Significância e Alegações
O artigo afirma fornecer um arcabouço em nível molecular para compreender como as propriedades específicas dos íons (valência, labilidade da camada de hidratação e mobilidade) governam a auto-montagem de polieletrólitos de mesma carga. Os autores enfatizam que os íons multivalentes não são intercambiáveis; o Ca²⁺, com sua camada de hidratação lábil e coordenação flexível, atua como um verdadeiro reticulador que estabiliza complexos e reorganiza a topologia dos polímeros, enquanto o Mg²⁺ oferece estabilização moderada através de interações transitórias, e o Na⁺ depende da blindagem.

O estudo destaca que separar os intervalos ligados e não ligados é essencial para caracterizar com precisão tais sistemas dinâmicos, pois métricas simples de tempo médio podem ocultar comportamentos fortemente dependentes de sal. Além disso, os resultados demonstram que a forte reticulação local (como visto em Ca²⁺) pode, paradoxalmente, aumentar o tamanho global da cadeia de DNA (expansão) enquanto compacta o polímero associado, desafiando interpretações simples da complexação apenas como condensação. Essas descobertas oferecem uma base para entender montagens de polímero–ácido nucleico responsivas a íons, relevantes para a entrega de genes, biosensoriamento e o design de biomateriais responsivos em ambientes iônicos multivalentes.