A kinetic model of shear-induced rupture of short dsDNA

Resumo Técnico: Um Modelo Cinético de Ruptura de dsDNA Curto Induzida por Cisalhamento

Definição do Problema
A dissociação induzida por força de DNA de fita dupla (dsDNA) curto é um processo fundamental na biofísica molecular e na nanotecnologia de DNA. Embora a espectroscopia de força de molécula única permita a extração de taxas de dissociação (off-rates) intrínsecas e distâncias de estado de transição, carecia-se de uma descrição cinética fisicamente fundamentada para a ruptura induzida por cisalhamento de construtos de comprimento finito. Os frameworks existentes frequentemente dependem de suposições que não consideram plenamente a interação entre a termodinâmica do duplex, a elasticidade dos segmentos que suportam a força e a geometria específica da aplicação da força. Especificamente, a tradução dos modelos termodinâmicos de vizinho mais próximo e de cinética de nucleação-zíper para a ruptura de duplex curtos dependente de força permanece incompleta.

Metodologia
Os autores desenvolveram um modelo cinético baseado em um framework de equação mestre para descrever a ruptura induzida por força de cisalhamento. O modelo é construído sobre uma via de nucleação-zíper dependente de força envolvendo transições de base única.

• Definição do Sistema: O sistema é modelado como um conjunto de estados rotulados por $(i, j)$, representando o número de pares de bases rompidos das extremidades esquerda e direita, respectivamente. Para um duplex de 11 pares de bases (pb), isso resulta em um sistema de 67 estados.
• Taxas de Transição: O modelo utiliza taxas de abertura ($k_o$), fechamento ($k_c$) e abertura terminal ($k_{ot}$). Estas taxas são derivadas da constante de equilíbrio, que incorpora mudanças de energia livre dependentes de força. A dependência da força surge do trabalho mecânico realizado durante a formação do par de bases, calculado via a mudança na distância extremidade-a-extremidade ($\Delta x$).
• Modelos de Polímero:
  • dsDNA: Modelado como um polímero do tipo bastão (rod-like). Crucialmente, o modelo considera a geometria helicoidal tridimensional para calcular a distância extremidade-a-extremidade, em vez de usar apenas um comprimento de contorno simples.
  • ssDNA: Modelado usando uma Fórmula de Interpolação Global (GIF) do modelo de cadeia tipo verme (worm-like chain - WLC), parametrizada com o comprimento de persistência ($\lambda_{ss}$) e a distância interfosfato ($l_{ss}$) determinados experimentalmente.
• Análise de Dados: O modelo foi aplicado a dados experimentais de Kabtiyal et al., envolvendo um construto de DNA-nanopartícula de ouro (AuNP)-DNA sob força de cisalhamento constante. Os autores reanalisaram os tempos de ruptura experimentais usando Estimativa de Máxima Verossimilhança (MLE) para derivar as taxas de dissociação de reação, abordando as limitações do estudo original onde os dados eram esparsos e agrupados (binned).
• Efeitos de Temperatura: O framework foi estendido para analisar tanto aumentos globais de temperatura quanto as complicações introduzidas pelo aquecimento plasmônico localizado em construtos acoplados a AuNP.

Principais Resultados

1. Reanálise de Dados Experimentais: Usando MLE, os autores determinaram uma distância de estado de transição (TSD) para a ruptura no escuro de $d_{Exp-MLE}^{Dark} = 0,85 \pm 0,20$ nm, que é maior que o anteriormente relatado de $0,60 \pm 0,16$ nm.
2. Precisão do Modelo: O modelo cinético, incorporando a geometria helicoidal, reproduziu com precisão a TSD experimental ($d_{Model}^{Dark} = 0,92$ nm) e as taxas de dissociação de reação através da faixa de força testada (2,5, 5 e 10 pN).
3. Importância da Geometia Helicoidal: Um achado central é que tratar o dsDNA como um bastão simples com comprimento de contorno leva a um desvio significativo na TSD calculada ($0,26$ nm) em comparação com os valores experimentais. Considerar a geometria helicoidal 3D é essencial para definir corretamente a resposta mecânica do dsDNA curto sob cisalhamento.
4. Comportamento do Regime de Força: O modelo revela que a relação entre $\ln k(F)$ e a força é não linear no regime de baixa força ($\le 10$ pN), apresentando uma curvatura pronunciada. No entanto, no regime de força moderada ($7,5 - 15$ pN), o comportamento torna-se linear, resultando em uma TSD constante por par de bases ($\approx 0,12$ nm/pb), consistente com estudos prévios de AFM e pinças ópticas.
5. Robustez aos Parâmetros de ssDNA: A TSD extraída mostrou-se robusta contra variações nos parâmetros do polímero ssDNA ($l_{ss}$ e $\lambda_{ss}$) dentro de faixas experimentalmente relevantes. Além disso, o modelo mostrou que a TSD é insensível ao comportamento de força-extensão de segmentos de ssDNA muito curtos ($<5$ nucleotídeos), sugerindo que modelos molecularmente refinados para esses comprimentos curtos não são estritamente necessários para este sistema.
6. Dependência de Temperatura: O modelo prevê uma correlação negativa entre temperatura e TSD (a TSD diminui conforme a temperatura aumenta). Ao ser aplicado aos dados de ruptura "excitada" (aquecimento plasmônico localizado), os autores descobriram que o aumento aparente na TSD relatado no estudo original era um artefato do método de análise. A reanálise via MLE não mostrou mudança significativa na TSD entre as condições de escuro e excitação. Os autores atribuem a falta de mudança observável na TSD a um aumento de temperatura real modesto (estimado em $\sim 30^\circ$C, ou um aumento de $\sim 1-2^\circ$C sobre o ambiente) e erro experimental, em vez de uma mudança fundamental no panorama cinético.

Significância e Alegações
O artigo estabelece um framework cinético preditivo e fisicamente fundamentado para interpretar experimentos de ruptura induzida por força em dsDNA curto. Os autores alegam que seu modelo fornece uma interpretação unificada de medições anteriores em diferentes regimes de força e esclarece os ingredientes mecânicos necessários para uma modelagem precisa.

• Acoplamento Mecânico-Cinético: O trabalho demonstra que a descrição mecânica do construto (especificamente a geometria helicoidal do dsDNA) é inseparável da análise cinética; negligenciar essa geometria resulta em parâmetros físicos incorretos.
• Interpretação Experimental: O framework oferece um método para distinguir entre efeitos térmicos globais e as complexidades do aquecimento plasmônico localizado, sugerindo que mudanças aparentes nos parâmetros cinéticos em experimentos excitados podem advir de variações no perfil de temperatura e artefatos de análise.
• Aplicação: Os resultados fornecem uma base para o design de nanoestruturas de DNA atuadas por força e temperatura, particularmente aquelas que integram componentes plasmônicos, ao oferecer um método confiável para interpretar dados de ruptura e prever o comportamento do sistema.