Reorganização Hierárquica Bifásica e Plasticidade Conectômica Riemanniana: Um Estudo Longitudinal de Mapeamento Funcional de Precisão da Psilocibina

1. Introdução: Dos Artefatos Euclidianos aos Manifolds Riemannianos

A psiquiatria contemporânea vivencia um renascimento nas pesquisas com psicodélicos serotoninérgicos, impulsionada pelo imperativo de tratar transtornos afetivos refratários. Contudo, os arcabouços computacionais que sustentam esse movimento permanecem majoritariamente limitados por metodologias estatísticas euclidianas que pressupõem uma variância de erro plana e homocedástica. Tais premissas sofrem colapso geométrico sob o massivo impacto entrópico da ativação de receptores $5\text{-HT}_{2\text{A}}$, gerando miopia metodológica (Arsigny et al., 2006) e frequentes falsos negativos topológicos (Nichols et al., 2017). Ferramentas lineares clássicas equivalem algebricamente ao uso de uma régua reta para mensurar distâncias sobre a superfície de uma esfera — uma abordagem que inflaciona artificialmente a variância do conectoma e mascara traços subagudos de plasticidade.

O objetivo central deste trabalho é divergir dessa tradição, tratando o conectoma humano como uma topologia matemática rígida residente em um manifold Riemanniano. Utilizando um conjunto de dados denso e longitudinal de mapeamento funcional de precisão (PFM) composto por $N = 101$ sessões de fMRI, busquei mapear com precisão a cinemática da neuroplasticidade induzida por psilocibina em três janelas temporais estritas: homeostase prévia (Linha de Base / Baseline), tempestade entrópica (Fase Aguda) e acomodação estrutural (Fase Persistente).

Para assegurar transparência inferencial e fidelidade anatômica, empreguei o parcelamento funcional Glasser/CABNP. O pipeline computacional implementa uma purga estrutural mandatória da parede medial executada na escala de alta densidade de grayordinates, excluindo 76 vértices não-teciduais do substrato de entrada da superfície crua (reduzindo a grade coordenada local de 91.282 para 91.206 grayordinates) para prevenir completamente a poluição de sinal da linha média. A projeção espacial downstream na configuração funcional CABNP mapeia esses vetores filtrados diretamente em sua partição invariante de 717 nós cortico-subcorticais. Essa rotina viabiliza o rastreamento preciso dos centros corticais ao longo de seus eixos estereotáxicos exatos. Ao impor integridade de manifold Positivo-Definido Simétrico (SPD) via encolhimento de Ledoit-Wolf, este estudo desenterra o sinal biológico puro da remodelação cortical de longo prazo que sobrevive semanas após a completa metabolização do composto.

Os Quatro Pilares da Descoberta Translacional

Ao alinhar rigor computacional com desfechos clínicos, este estudo estabelece um panorama translacional fundamentado em quatro pilares principais:

2. Métodos: Arquitetura de Mapeamento Funcional de Precisão

Para avaliar as modificações topológicas induzidas por psilocibina, a estrutura analítica abordou o problema da alta dimensionalidade em neuroimagem. Ao estimar redes funcionais cerebrais usando parcelamentos espaciais extensos em intervalos temporais discretos, as matrizes de covariância empíricas resultantes podem se tornar mal-condicionadas (Varoquaux & Craddock, 2013).

Garantir a integridade do manifold Positivo-Definido Simétrico (SPD) fornece uma base matematicamente consistente para rastrear distâncias geodésicas de longo prazo entre janelas temporais (Arsigny et al., 2006).

Matriz da Coorte e Arquitetura Longitudinal

Os dados foram obtidos de um framework de amostragem densa de mapeamento funcional de precisão:

• Matriz da Coorte: A amostra consistiu em 7 voluntários adultos saudáveis ($M = 34,1\text{ anos}$, $DP = 9,8$; $n = 3\text{ mulheres}$; 6 participantes identificados como caucasianos).
• Arquitetura Longitudinal: O protocolo de imagem abrangeu 11 fluxos longitudinais de sujeitos capturando três janelas temporais: homeostase basal (Linha de Base), janela de pico agudo (60–90 min pós-ingestão; Fase Aguda) e janela de acomodação subaguda (até 2 semanas pós-exposição; Fase Persistente).
• Protocolo de Replicação: Um subconjunto de 4 indivíduos (P1, P3, P4, P5) foi submetido a um protocolo de replicação após um período de washout superior a 6 meses, totalizando 101 sessões densas de neuroimagem funcional.

3. Resultados: Transformações Topológicas nas Janelas Temporais

Integridade Riemanniana & Descoberta da Reorganização Hierárquica Bifásica

A arquitetura de regularização assegurou com sucesso o status de posto completo e a integridade de manifold Positivo-Definido Simétrico (SPD) em todas as 101 sessões ($\lambda_{\min} > 0$). As matrizes de consenso em nível de grupo convergiram com estabilidade: Linha de Base ($\lambda_{\min} = 1,042 \times 10^{-4}$), Fase Aguda ($\lambda_{\min} = 1,037 \times 10^{-4}$) e Fase Persistente ($\lambda_{\min} = 1,042 \times 10^{-4}$).

A distância entre os polos transmodal e unimodal, que se situava em uma média de 0,002707 na Linha de Base, contraiu-se para 0,000841 durante a Fase Aguda ($\Delta = -68,95\%$), refletindo uma compressão transitória da hierarquia funcional. Na Fase Persistente pós-aguda, a rede exibiu expansão subaguda, atingindo média polar de 0,002947 ($\Delta = +250,66\%$ em relação ao estado agudo).

Desintegração Global e a Avalanche Conectiva Aguda

Durante o pico agudo, o índice de segregação local da DMN sofreu redução expressiva, passando de uma média basal de 0,4608 para 0,2547. Concomitantemente, o coeficiente de participação médio global elevou-se de 0,1910 para 0,2816. A deformação de rede global (RiemannianPlasticity) saltou de 34,0016 para um pico agudo de 53,0813.

Na janela de acompanhamento persistente, a segregação da DMN recuperou-se para 0,4800, enquanto a deformação global acomodou-se em 43,7956, com uma forte redução na variância entre sessões ($DP = 4,0266$, comparado a $DP = 10,5015$ na linha de base).

O Espaço de Interação Fenomenológica (MEQ-30)

Avaliei a relação entre escores místicos subjetivos e a remodelação conectômica objetiva. A coorte exibiu ampla variância subjetiva no MEQ-30 ($M = 3,0697$, $DP = 1,2623$, intervalo 1,1333 – 4,4000).

A modelagem longitudinal via GEE revelou que o efeito principal de longo prazo do escore MEQ-30 centrado na média foi estatisticamente não-significativo ($b = 0,0315$, $p = 0,513$). Entretanto, emergiu uma interação traço-por-estado altamente significativa durante o pico agudo: o deslocamento agudo de GradientShift foi fortemente previsto pela interação entre o estado temporal e a psicometria centrada ($b = 0,1550$, $p = 0,008$). Na janela persistente, esse valor preditivo dissipou-se ($b = 0,0341$, $p = 0,137$).

Modelos de Regressão Alternativos & Robustez Heterocedástica

Para tratar a variância não-constante sem descarte arbitrário de dados, comparei o modelo OLS padrão (Modelo D.2) contra um Modelo Linear Generalizado Gamma com função de ligação log e estimador de covariância sanduíche robusto de Huber-White HC3 (Modelo D.1).

Ambos os modelos convergiram no mesmo comportamento: a perturbação aguda foi altamente significativa ($p < 0,001$), enquanto o estimador da fase persistente cruzou a linha nula ($p = 0,963$ no OLS, $p = 0,981$ no Gamma GLM).

Trajetórias de Sub-redes & Assimetria Transmodal (Depressão vs. Ansiedade)

Avaliando 12 nós de controle transmodal (DMN/FPN, regulação do humor) contra 12 nós de saliência e paralímbicos (complexos insula/amígdala, hipervigilância), observou-se acentuada divergência:

• O polo cortical transmodal Default-66_L-Ctx exibiu um vetor de deslocamento médio de 0,2133 ($\Delta = 0,2963$ Baseline → Agudo).
• O polo de saliência Cingulo-Opercular-42_L-Ctx atingiu apenas 0,1255 ($\Delta = 0,1901$ Agudo → Persistente).

O teste não-paramétrico de Mann-Whitney U rejeitou enfaticamente a distribuição espacial homogênea: $U = 139,00$, $p < 0,001$, correlação bisserial por postos $r = -0,9306$. Circuitos transmodais de depressão absorvem significativamente mais maleabilidade que as vias de ansiedade.

4. Discussão: Da Geometria Conectômica à Biologia Translacional

Estabelecidas as bases topológicas e estatísticas do deslocamento neuronal por psilocibina, o imperativo translacional consiste em traduzir tensores e manifolds em biologia médica. Decodificar a compressão aguda e o subsequente rebote hiper-expansivo exige acoplar essa deformação arquitetônica à farmacodinâmica primária dos receptores $5\text{-HT}_{2\text{A}}$ e definir seu impacto no prognóstico clínico.

Deslocamento Neuronal e Validação do Modelo REBUS

Os vetores de deslocamento de coordenadas ao longo do gradiente principal indicam ajustes precisos nas hierarquias corticais. A contração observada na distância polar (0,002707 → 0,000841, $\Delta = -68,95\%$) sinaliza uma redução temporária na diferenciação hierárquica macroescalar.

Em repouso típico, o córtex mantém fronteiras topológicas estritas separando vias sensoriais primárias de redes transmodais de alta ordem (Margulies et al., 2016). Em transtornos afetivos, essas redes tornam-se hiper-rígidas, perpetuando ruminações depressivas. Os dados agudos confirmam o relaxamento dessas barreiras: a segregação da DMN despenca de 0,4608 para 0,2547, enquanto a participação média global sobe de 0,1910 para 0,2816. Isso fornece corroboração empírica direta ao modelo REBUS (RElaxed Beliefs Under pSychedelics) (Carhart-Harris & Friston, 2019): o agonismo $5\text{-HT}_{2\text{A}}$ suspende momentaneamente os priors rígidos de topo, tornando os circuitos receptivos à reestruturação cognitiva.

O Dogma Místico: Reset Fisiológico vs. Fenomenologia Aguda

Um modelo amplamente divulgado sugere que a duração das alterações terapêuticas depende diretamente da intensidade subjetiva da experiência mística (Carhart-Harris et al., 2017). Os dados deste estudo refinam essa proposição:

1. A remodelação conectômica pós-aguda de longo prazo é estatisticamente ortogonal aos escores de MEQ-30 centrados ($p = 0,513$), apontando para um mecanismo de reinicialização neurobiológica autônomo.
2. Por outro lado, o deslocamento agudo de gradiente é fortemente previsto pela interação traço-por-estado ($p = 0,008$).

A intensidade consciente espelha a magnitude momentânea da desestabilização topológica, mas a retenção de longo prazo é regida por fatores biológicos de plasticidade e suporte ambiental (Letheby, 2021).

Suscetibilidade Regional & A Assimetria de Sub-redes

A identificação não-paramétrica da assimetria transmodal ($U = 139,00$, $p < 0,001$, $r = -0,9306$) desmistifica a hipótese de um reset cerebral indiferenciado. Conjuntos transmodais de alta ordem (DMN/FPN; ex.: Default-66_L-Ctx, $M = 0,2133$) exibem maleabilidade e deriva posicional significativamente maiores que estruturas de saliência ou límbicas (Cingulo-Opercular-42_L-Ctx, $M = 0,1255$).

Essa suscetibilidade diferencial espelha o gradiente de densidade de receptores $5\text{-HT}_{2\text{A}}$, justificando por que a psilocibina atua de maneira tão contundente em esquemas ruminativos de depressão, demandando protocolos psicoterapêuticos mais focados para quadros de ansiedade generalizada.

5. Conclusão e Fronteiras Futuras

Em síntese, o rastreamento regularizado em 101 sessões de mapeamento funcional de precisão define três fases cinéticas fundamentais:

1. Fase Basal: Redes funcionais operam em modularidade típica, preservando fronteiras estritas entre vias sensoriais e redes de controle transmodal.
2. Fase Aguda: Compressão marcante do gradiente hierárquico ($\Delta = -68,95\%$) e dessegregação da DMN (0,4608 → 0,2547), promovendo comunicação de alta entropia e relaxando restrições de topo.
3. Fase Persistente: Expansão de gradiente subaguda ($\Delta = +250,66\%$) e canalização em um atrator de variância estreitada ($CV = 9,19\%$), retendo maleabilidade individualizada.

Limitações Metodológicas e Recomendações

• Amostra vs. Precisão: Embora 101 sessões densas ofereçam poder estatístico intra-sujeito incomparável, a amostra ancora-se em $N = 7$ voluntários únicos. Validações em coortes clínicas mais amplas são recomendadas.
• Integração Multimodal hdEEG-fMRI: A combinação de fMRI denso com EEG de alta densidade permitirá rastrear o colapso do gradiente em resolução de sub-milissegundos.
• Especificidade Farmacológica: Contrastar psilocibina contra cetamina, LSD e MDMA determinará se a assimetria transmodal é exclusiva do agonismo $5\text{-HT}_{2\text{A}}$ ou representa uma propriedade geral de psicoplastógenos rápidos.

💻 Disponibilidade de Dados e Código

• Conjunto de Dados Brutos de Neuroimagem: OpenNeuro ID de Acesso ds006072 (Psilocybin Precision Functional Mapping database).
• Código-Fonte do Pipeline Computacional: O pipeline algorítmico completo e versionado (projeção CABNP, shrinkage de Ledoit-Wolf, tensores Riemannianos AIRM) está disponível publicamente no GitLab:
  https://gitlab.com/lunalab/macroscale-cortical-gradient-disintegration-and-riemannian-connectomics-pipeline/-/tree/v.1.0.0?ref_type=tags
• Arquivo Persistente Zenodo: DOI: 10.5281/zenodo.22144006

📚 Referências Bibliográficas