Two observables of one wall: how surface relaxivity can bias the diffusion intra-axonal fraction and the myelin water fraction

Resumo Técnico: O Viés da Relaxividade de Superfície nas Frações de Microestrutura de Difusão e Relaxometria

Definição do Problema
A relaxividade de superfície ($\rho$) e a difusão dependente do tempo são dois observáveis que surgem do mesmo processo físico: colisões de moléculas de água com as paredes de axônios mielinizados. Embora a dMRI (ressonância magnética de difusão) seja frequentemente considerada "cega à relaxividade" devido à normalização por $b=0$, este artigo argumenta que tal suposição só é válida para um único compartimento. Na realidade, a água intra-axonal e a extra-axonal residem contra paredes distintas com razões superfície-volume ($S/V$) distintas. Consequentemente, elas experimentam taxas de relaxação transversal diferentes ($R_2 = 1/T_2 + \rho S/V$). Esse $T_2$ diferencial repondera os compartimentos de sinal no tempo de eco ($TE$) utilizado para a normalização, introduzindo um viés sistemático nas estimativas de microestrutura. O artigo investiga como esse termo geométrico não observado enviesa duas métricas fundamentais: a fração de sinal intra-axonal de difusão ($f_{intra}$) e a fração de água da mielina (MWF).

Metodologia
O estudo combina derivações analíticas com simulações de Monte Carlo (MC) de primeira ordem para validar a física e quantificar o viés.

1. Derivações Analíticas:

   • Taxa de Interior: Derivou a taxa de superfície de Brownstein–Tarr para a água intra-axonal sobre uma população de calibres axonais distribuída via Gamma. A taxa é a média ponderada pela área de $\rho(4/d)$.
   • Taxa de Exterior: Derivou a taxa de superfície de Novikov–Burcaw para a água extra-axonal, que satura em um patamar determinado pela densidade de superfície exterior média ($S_{ext}/V_{ext}$) no limite de estreitamento de movimento (motional-narrowing).
   • Razão de Forma Fechada: Demonstrou que a razão entre o $S/V$ interior e o exterior é independente da distribuição de calibre e depende exclusivamente da fração de volume das fibras ($V_F$) e do $g$-ratio: $(1-V_F)/(g V_F)$.
   • Fórmulas de Viés: Derivou expressões de forma fechada para o viés em $f_{intra}$ (Eq. 9) e o deslocamento na MWF aparente através do cruzamento dos valores de $T_2$ de axônios finos abaixo da janela da mielina.

2. Validação por Monte Carlo:

   • Desenvolveu um simulador de contagem de paredes por MC (usando um estimador de tempo local de fronteira de realized-overshoot) para propagar caminhantes aleatórios através de geometrias explícitas compactadas sem modelos analíticos de relaxividade.
   • Validou que a atenuação de superfície do MC corresponde às taxas de Brownstein–Tarr e Novikov–Burcaw dentro de $\sim0,2$ pontos percentuais.
   • Confirmou que axônios finos na simulação de MC exibem valores de $T_2$ que caem abaixo da janela de 25 ms da mielina, validando o mecanismo para o viés de MWF.

3. Modelagem Direta (Forward Modeling):

   • Difusão: Simulou sinais de dMRI multi-casco com relaxividade de superfície habilitada e os ajustou usando técnicas padrão de média esférica (SMT) para quantificar o viés de $f_{intra}$.
   • Relaxometria: Simulou decaimentos CPMG multi-eco com um pool intra-axonal de calibre distribuído e estimou a MWF usando NNLS (mínimos quadrados não negativos) padrão com uma janela de mielina fixa.
   • Trajetória Longitudinal: Modelou uma trajetória de desmielinização intra-sujeito (preservando o lúmen do axônio, afinando a mielina) para testar se o viés se cancela em mudanças longitudinais.

Resultados Principais

• Viés da Fração de Sinal Intra-axonal ($f_{intra}$):

  • O viés é impulsionado pela diferença de $S/V$ entre as paredes interior e exterior.
  • Lei do Sinal: O sinal do viés é determinado pela fração de volume de fibras ($V_F$). Abaixo de um ponto de cruzamento $V_F^* = 1/(1+g) \approx 0,59$, a água intra-axonal relaxa mais rápido (subestimação). Acima de $V_F^*$, a água extra-axonal relaxa mais rápido, causando uma sobreponderação do sinal intra-axonal.
  • Magnitude: Em substância branca fisiologicamente densa ($V_F \approx 0,65$), a fração intra-axonal é superestimada em aproximadamente 12% em um tempo de eco clínico ($TE=80$ ms) usando um valor de relaxividade citado ($\rho = 1,16$ $\mu$m/s). Este viés escala linearmente com $\rho$ e cresce monotonicamente com o $TE$.
  • Natureza: Este é um erro sistemático de primeira ordem e espacialmente estruturado, não um ruído de média zero, o que significa que ele não se anula entre os voxels, mas desloca as médias regionais e de coorte.

• Viés da Fração de Água da Mielina (MWF):

  • Mecanismo: A relaxividade de superfície encurta o $T_2$ aparente da água intra-axonal. Para os axônios mais finos (diâmetro interno $d < d^* \approx 0,17$ $\mu$m para o $\rho$ citado), o $T_2$ cai abaixo da janela padrão de 25 ms da mielina. Esta água é classificada incorretamente como água da mielina.
  • Direção: A substância branca fina (alto $S/V$) lê-se como mais rica em mielina do que a substância branca grossa, mesmo que o volume real de mielina seja idêntico.
  • Magnitude: No $\rho$ citado, o viés é pequeno (0,33 pontos percentuais em MWF, ou ~0,82 pp em conteúdo de mielina inferido). No entanto, é super-linear em relação a $\rho$; se $\rho$ for um fator de dez maior, o viés torna-se significativo (3,4 pp).
  • Detectabilidade: O viés está oculto sob o ruído de voxel único, mas emerge como um deslocamento sistemático ao realizar a média sobre regiões de interesse (ROIs) com diferentes distribuições de calibre.

• Demiesinização Longitudinal:

  • Em uma trajetória de desmielinização que preserva o lúmen (desmielinização primária), o viés de MWF é definido pelo diâmetro interno do axônio fixo. Assim, o viês é quase constante e se cancela ao calcular a mudança ($\Delta$MWF) entre os pontos temporais.
  • Inversamente, o viés de $f_{intra}$ depende da diferença entre o $S/V$ interior e exterior. À medida que a mielina afina, a parede exterior recua, alterando o $S/V$ exterior enquanto o interior permanece fixo. Consequentemente, o viés de $f_{intra}$ desloca-se significativamente (ex: de -5% para -21%) durante a desmielinização, criando um grande artefato dependente de $TE$ em estudos de difusão longitudinais.

Significância e Alegações
O artigo afirma identificar uma fonte de viés anteriormente não relatada tanto em imagens de microestrutura de difusão quanto de relaxometria, que decorre da inseparabilidade da $T_2$ de bulk e da relaxividade de superfície dependente do tempo em medições de múltiplos ecos.

1. Física Unificada: Estabelece que a relaxividade de superfície e a difusão dependente do tempo são leituras duais do mesmo processo de colisão de parede. O "$T_2$" medido em relaxometria é contaminado por uma taxa dependente da geometria que a dMRI tenta normalizar, mas falha em eliminar totalmente devido às diferenças de $S/V$ dos compartimentos.
2. Erros Sistemáticos: O viés de $f_{intra}$ é caracterizado como um erro sistemático de primeira ordem e determinístico que escala com a densidade de empacotamento e o tempo de eco, podendo confundir comparações entre centros ou entre coortes se os tempos de eco diferirem.
3. Interpretação da MWF: O viés de MWF é enquadrado como um confundimento estrutural onde a "água da mielina" mede efetivamente a "água aprisionada de alto $S/V$". Isso implica que regiões de substância branca fina (ex: genu do corpo caloso) podem parecer ter maior conteúdo de mielina do que regiões de substância branca grossa (ex: corpo médio), independentemente do volume real de mielina.
4. Contraste Longitudinal: Um achado fundamental é a divergência em como esses vieses se comportam longitudinalmente. Enquanto as mudanças de MWF em desmielinização primária são robustas a este viés (pois o deslocamento se cancela), as mudanças de $f_{intra}$ baseadas em difusão são fortemente contaminadas pelo deslocamento da geometria da parede exterior.
5. Limitações: Os autores observam modestamente que a magnitude absoluta do viés depende do valor mal estabelecido de $\rho$ (que pode variar em ordens de magnitude) e que a geometria idealizada (cilindros que não se tocam) pode subestimar os efeitos em pacotes de fibras realistas e em contato. Eles não propõem um novo protocolo de aquisição para resolver isso, observando que a degenerescência é estrutural e requer restrições ortogonais (ex: gradientes oscilantes de alto gradiente) que estão atualmente além do alcance clínico.