Todos os experimentos foram aprovados pelo Conselho de Revisão Institucional da Universidade de Medicina Chinesa de Nanjing (Código de Protocolo: 202501A018) e foram conduzidos de acordo com o Guia dos Institutos Nacionais de Saúde (NIH) para o Cuidado e Uso de Animais de Laboratório. Todos os procedimentos experimentais seguiram as diretrizes ARRIVE.
Medicamentos experimentais
A solução amostral YQHXD continha 20 g de raiz de Astragalus membranaceus , 20 g de raiz de Angelica sinensis , 15 g de rizoma Ligusticum chuanxiong , 6 g de raiz de Stephania tetrandra , 15 g de fruto de Chaenomeles speciosa e 10 g de sementes maduras de Sinapis alba , com um peso total bruto de 86 g por prescrição. Todos os materiais herbais foram autenticados de acordo com a Farmacopeia Chinesa (edição de 2020) pelo Departamento de Farmácia, Hospital Ruikang, afiliado à Universidade de Medicina Chinesa de Guangxi. A solução amostra YQHXD foi preparada usando um método de decocção em água. Água destilada estéril com 8 vezes o volume total de drogas brutas foi adicionada, e as ervas foram deixadas de molho por 30 minutos até ficarem totalmente submersas. A mistura foi então fervida e decoctada em fogo baixo por 30 minutos antes da filtração. Os resíduos foram extraídos novamente com o mesmo volume de água destilada estéril após o mesmo procedimento. Os filtrados das duas extrações foram combinados e concentrados a 55,5 mL sob pressão reduzida usando evaporação rotativa para obter a solução original YQHXD. Com base no peso total do medicamento bruto, a concentração final de medicamento bruto foi de 1,55 g/mL, com rendimento de extração de aproximadamente 64,5%. A dosagem administrada a ratos experimentais foi convertida a partir da dose clínica adulta usando o método de normalização da área superficial corporal. O coeficiente de conversão entre humanos e ratos era aproximadamente 6,3; Portanto, a dosagem para ratos (mg/kg) foi calculada como a dosagem humana (mg/kg) × 6,3. O grupo de alta dose recebeu 15,48 g/kg/dia como dose de referência, enquanto os grupos de médio e baixo receberam 7,74 g/kg/dia e 3,87 g/kg/dia, respectivamente. Para um volume de gavage de 10 mL/kg, as concentrações de trabalho correspondentes foram 1,548, 0,774 e 0,387 g/mL para os grupos de alta, média e baixa dose, respectivamente. A solução padrão foi diluída de acordo para atingir a dosagem-alvo para cada grupo experimental antes da administração. Após a preparação, as decocções eram resfriadas, aliquotadas em condições estéreis e armazenadas a −40 °C até o uso.
Análise de farmacologia de redes
Previsão das proteínas-alvo do princípio ativo
Em 24 de novembro de 2025, compostos ativos foram recuperados do banco de dados TCMSP usando os termos de busca "Huangqi, Angelica, Chuanxiong, Stephania, Mugua e Jiezi." Os compostos foram avaliados com base na biodisponibilidade oral (OB ≥ 30%) e na semelhança do medicamento (DL ≥ 0,18)9. A justificativa para esses limites é a seguinte: o banco de dados TCMSP integra dados de compostos de fontes autoritativas como PubChem, ChEMBL e ChemSpider e fornece parâmetros relacionados ao ADME, incluindo biodisponibilidade oral. Biodisponibilidade oral refere-se à proporção de uma dose administrada oralmente que chega à circulação sistêmica sem alterações, e um alto índice de obstetricidade é essencial para identificar moléculas bioativas. A semelhança com medicamentos é um indicador usado para avaliar a capacidade de uso de compostos e prever sua solubilidade e estabilidade. Neste estudo, foi aplicado um limiar OB mais rigoroso (≥30%), que exceda o valor recomendado pelo banco de dados (≥20%), para identificar compostos candidatos com absorção oral favorável e perfis farmacocinéticos. Enquanto isso, DL foi usado para estimar o potencial semelhante a fármacos dos compostos, e um corte de DL ≥ 0,18 foi adotado com base no valor médio de DL observado em compostos da medicina tradicional chinesa. Os compostos ativos selecionados foram recuperados no formato SMILES do banco de dados PubChem e importados para o banco de dados SwissTargetPrediction para previsão de alvos, com a espécie definida como Homo sapiens. As proteínas-alvo foram padronizadas usando símbolos genéticos oficiais, e apenas alvos com valores de probabilidade maiores que 0 foram mantidos para melhorar a confiabilidade da previsão. Compostos sem alvos previstos foram excluídos da análise subsequente da rede. Após a fusão de todos os alvos, entradas duplicadas foram removidas usando software de planilha para garantir que cada alvo fosse representado apenas uma vez no conjunto de dados final.
Coleção de alvos-chave
Genes relacionados à doença foram recuperados dos bancos de dadosGeneCards 10, OMIM11 e TTD12 usando a palavra-chave "degeneração do disco lombar". Os alvos gênicos coletados foram fundidos e entradas duplicadas removidas para obter um conjunto final de alvos associados à doença.
Potenciais alvos da YQHXD para o tratamento da LDD
A plataforma Venny 2.1.0 foi usada para determinar os potenciais alvos do YQHXD no LDD. Uma rede foi então gerada para mostrar as associações entre os componentes ativos e seus alvos previstos. Os componentes ativos centrais foram selecionados de acordo com seus valores de grau na rede.
Construção de rede de componentes de medicamentos e alvos de interseção de doenças
Ingredientes ativos, proteínas-alvo e genes relacionados ao LDD da fórmula TCM YQHXD foram importados para o softwareCytoscape 13 para construir um diagrama de rede de componente medicamento-alvo-doença.
Construção da rede PPI e identificação de alvos centrais chave
Alvos que se cruzavam eram submetidos à plataformaSTRING 14 para análise, com o organismo configurado para Homo sapiens. Para construir a rede de interação proteína-proteína (PPI), o limiar de pontuação de interação foi definido em >0,4, nós desconectados foram ocultados e todas as outras configurações foram mantidas como padrão. Os dados resultantes foram importados para o software, onde a ferramenta NetworkAnalyzer foi usada para avaliar propriedades topológicas e calcular os graus dos nós. Os cinco principais nós classificados por grau foram identificados como alvos centrais centrais.
Análise de enriquecimento de alvos centrais
Alvos que se cruzaram foram enviados para o banco de dadosDAVID 15 para anotação funcional de Ontologia Gênica (GO) e análise de enriquecimento de vias da Kyoto Encyclopedia of Genes and Genomes (KEGG). O fundo do organismo foi definido como Homo sapiens, com o corte de enriquecimento definido de acordo com os parâmetros padrão de DAVID, incluindo um limiar de EASE < 0,1 e um limiar de contagem ≥ 2. O universo genético foi baseado no fundo padrão do genoma completo fornecido pelo DAVID. Testes múltiplos foram corrigidos usando o método de Benjamini-Hochberg para controlar a taxa de falsas descobertas (FDR), e os resultados do DAVID foram ajustados com base nisso. As vias enriquecidas foram ordenadas pelo valor de contagem em ordem decrescente. Os plots de barras de enriquecimento GO e KEGG foram produzidos usando uma plataforma online de bioinformática.
Acoplamento molecular
Para prever interações entre os principais componentes ativos do YQHXD e os alvos centrais do LDD, foi realizada uma análise de acoplamento molecular neste estudo. Os cinco principais compostos ativos, classificados de acordo com valores de grau topológico da rede e relevância do alvo, foram selecionados como ligandos, enquanto as respectivas proteínas alvo centrais foram usadas como receptores. Arquivos SDF de ligantes foram recuperados do PubChem, e as estruturas tridimensionais das proteínas alvo foram baixadas do Banco de Dados de Proteínas (PDB). A análise de acoplamento foi realizada na plataforma CB-Dock2 usando seu protocolo totalmente automatizado de acoplamento cego. Após o upload das estruturas proteicas, a plataforma removeu automaticamente os heteroátomos, incluindo pequenas moléculas, íons e moléculas de água, e posteriormente reparou átomos ausentes e átomos de hidrogênio. O modo Auto Blind Docking foi usado. Nesse modo, o CB-Dock2 identifica potenciais bolsões de ligação à superfície da proteína com um algoritmo de detecção de cavidades baseado em curvatura e então realiza cálculos de acoplamento usando o motor integrado AutoDockVina 16. Durante o processo de acoplamento, as proteínas foram tratadas como estruturas rígidas, enquanto os ligantes foram considerados flexíveis, e todos os parâmetros de acoplamento foram definidos para os valores padrão da plataforma. Os resultados do acoplamento foram submetidos a visualização tridimensional e análise de interação usando as ferramentas integradas do CB-Dock217, e as energias de ligação foram calculadas para estimar as afinidades de ligação previstas entre os compostos ativos e as proteínas alvo.
Simulação de dinâmica molecular
Simulações de dinâmica molecular foram realizadas para o complexo proteína-ligante-alvo usando o GROMACS 2022.2. A proteína foi modelada com o campo de força Amber14SB, e o sistema de solventes foi representado usando o modelo de água TIP3P. Cargas atômicas de ligantes foram calculadas usando o método AM1-BCC implementado no Antechamber, seguido pela atribuição dos tipos de átomos GAFF2. Os arquivos de topologia resultantes foram convertidos para um formato compatível com GROMACS usando ACPYPE. Parâmetros iônicos consistentes com o modelo de água TIP3P foram atribuídos com base no conjunto de parâmetros de Joung-Cheatham. O complexo inicial foi colocado em uma caixa de simulação dodecaedrica truncada com pelo menos 1,2 nm entre a superfície da proteína e a borda da caixa. Moléculas de água TIP3P foram adicionadas para solventar o sistema, e íons Na+ e Cl− foram introduzidos para neutralizar a carga do sistema e ajustar a força iônica para 0,15 M. A minimização de energia foi realizada com o algoritmo de descida mais íngreme até que a força máxima ficasse abaixo de 1.000 kJ·mol−1·nm−1. O sistema foi então equilibrado para 100 ps sob o conjunto NVT a 298 K e 1 bar. Após o equilíbrio, uma simulação de produção de 100 ns foi executada sob o conjunto NPT com um passo de tempo de 2 fs. Interações não vinculadas eram tratadas com o esquema de corte Verlet. Interações eletrostáticas foram calculadas usando o método de Ewald de malha de partículas (PME), e um cutoff de 1,2 nm foi usado tanto para as interações de van der Waals quanto de Coulomb. As ligações contendo hidrogênio foram restringidas pelo algoritmo LINCS. A temperatura era controlada com o termostato Nosé-Hoover, e a pressão era controlada com o barostato Parrinello-Rahman, mantendo o sistema em 298 K e 1 bar. Trajetórias eram salvas a cada 10 ps. Dinâmicas estruturais e interações intermoleculares foram analisadas usando ferramentas integradas do GROMACS, VMD e PyMOL. A energia livre de ligação foi estimada com gmx_MMPBSA quando aplicável.
Estudos com animais
Criação de animais
Setenta ratos machos de Sprague-Dawley (SD) entre 3 e 8 semanas foram aleatoriamente atribuídos a sete grupos experimentais (n = 10 por grupo) usando números aleatórios gerados por computador, garantindo que não haja diferenças estatisticamente significativas no peso corporal inicial entre os grupos. Todos os animais foram alojados sob condições controladas (22 ± 1 °C, 55 ± 10% de umidade relativa, 12 horas de ciclo claro/escuro) com livre acesso a comida e água. As posições das gaiolas foram randomizadas a cada mudança de gaiola para minimizar o viés ambiental. Após 1 semana de aclimatação, ratos dos grupos modelo passaram por punção caudal por IVD para induzir degeneração. Anestesia local de infiltração com bupivacaína foi administrada pré-operatório para reduzir a dor. Os animais foram monitorados de perto durante a cirurgia e durante toda a recuperação pós-operatória. As observações pós-operatórias incluíram tempo de recuperação, restauração da atividade espontânea, consumo de alimentos e água, e mudanças no peso corporal. Uma semana após o estabelecimento do modelo, dois ratos de cada grupo foram selecionados aleatoriamente para radiografia digital pós-operatória (DR) e coloração por hematoxilina e eosina (HE). A análise radiográfica mostrou redução da altura do IVD caudal, enquanto a coloração com HE revelou uma redução significativa nas células notocordiais do núcleo pulposo (NP) e a desorganização da estrutura lamelar do anel fibroso18, confirmando o estabelecimento bem-sucedido do modelo. Oito semanas após a cirurgia, ratos foram eutanasiados por inalação de uma overdose de isoflurano, e amostras de tecido caudal e sangue foram coletadas. Após validação preliminar, foram usadas três amostras independentes por grupo que atendessem aos critérios de qualidade pré-definidos para as análises finais de histologia, ELISA, PCR e western blot. As seções histológicas foram avaliadas sob condições de simples cegueira, com lâminas rotuladas apenas por números de identificação dos animais. Para ensaios ELISA, PCR e western blot, a codificação cega foi aplicada durante a alicotação e rotulagem da amostra, e a ordem de teste de todas as amostras foi randomizada antes dos experimentos. Dados brutos, incluindo valores de absorção, valores Ct e intensidades de banda, foram decodificados em identidades de grupo por um estatístico não envolvido em agrupamento experimental ou processamento de amostras e posteriormente analisados por pessoal estatístico. Todas as análises estatísticas foram realizadas por estatísticos que não participaram do agrupamento experimental ou processamento da amostra. Quaisquer eventos adversos, incluindo infecção de ferida, deiscência, sangramento ou letargia severa, foram imediatamente registrados e avaliados. Condições leves foram monitoradas de perto com cuidados de suporte aprimorados. Critérios humanos foram aplicados se ratos apresentassem perda de peso corporal superior a 20% do valor inicial, dor intensa não aliviada por analgésicos, perda de mobilidade, incapacidade de acessar alimentos ou água, infecção, necrose, outras complicações pós-operatórias graves, dificuldade respiratória ou qualquer condição potencialmente fatal. Nesses casos, os animais eram eutanasiados humanamente por inalação de uma overdose de isoflurano.
Modelagem experimental
Os ratos foram anestesiados com isoflurano inalado e posicionados de bruços. Os IVDs de cauda Co6/7, Co7/8 e Co8/9 estavam localizados sob orientação fluoroscópica de raios X. Uma agulha de 21 G era avançada verticalmente para dentro do IVD até alcançar o lado oposto da pele, após o que a ponta da agulha era retraída em direção ao centro da cauda. Uma segunda agulha servia como referência posicional. Para induzir a degeneração sustentada da DIVI, a agulha era girada 360°, mantida no lugar por 30 segundos e depois removida. No grupo Sham, a profundidade da perfuração através da pele e músculo foi de aproximadamente 2 mm, sem danificar o DI. Nos outros grupos, a agulha foi inserida 5 mm em direção ao centro do DIVI, evitando os vasos sanguíneos. Após a inserção, a agulha era girada 360° e deixada no lugar por 30 segundos. O local da perfuração foi desinfetado, marcado e protegido com gaze estéril. Uma semana após a modelagem, a imagem por raio-X da coluna caudal mostrou diminuição na altura do DIVI. A coloração com HE revelou uma redução significativa nas células notocordiais NP e desorganização da estrutura da lamela anular, confirmando o estabelecimento bem-sucedido do modelo. Após modelagem bem-sucedida, os grupos que receberam tratamentos com MTC e o grupo controle positivo (PC) receberam as doses apropriadas, calculadas com base na razão da área superficial corporal entre humanos e ratos.
Agrupamento e intervenção
Um total de 70 ratos machos SD de 3 a 8 semanas foram usados neste experimento, divididos em sete grupos:
(1) Grupo controle (Ctrl): Ratos deste grupo foram designados aleatoriamente ao mesmo tempo e sob as mesmas condições dos grupos experimentais. Eles eram alojados em gaiolas individuais em uma instalação animal de nível SPF, com acesso livre a comida e água potável padrão.
(2) Grupo simulado: Os ratos foram anestesiados por inalação de isoflurano. Após a anestesia ser estabelecida, a pele da cauda foi desinfetada com povidona-iodo, e os ratos foram posicionados pronos. Os IVDs Co6/7 e Co7/8 foram localizados sob orientação fluoroscópica de raios X. A agulha foi avançada aproximadamente 2 mm através da pele e do músculo sem ferir o DI. Água destilada (10 mL/kg) era administrada uma vez ao dia por gavage durante 8 semanas.
(3) Grupo modelo: Os ratos foram anestesiados com isoflurano inalado. Após a anestesia ser estabelecida, a pele da cauda foi desinfetada com povidona-iodo, e os ratos foram posicionados pronos. Sob orientação fluoroscópica de raios X, os IVDs Co6/7, Co7/8 e Co8/9 foram localizados. Uma agulha de 21 G era avançada verticalmente para dentro do DIV até alcançar o lado oposto da pele, e a ponta da agulha era então retirada em direção ao centro da cauda.
Uma segunda agulha era usada como referência posicional. A agulha era girada 360° e mantida no lugar por 30 segundos para induzir degeneração contínua do DIV antes da remoção. Água destilada (10 mL/kg) foi administrada diariamente por gavage durante 8 semanas, e a recuperação dos ratos foi monitorada de perto após a cirurgia.
(4) Grupo PC: O grupo controle positivo foi submetido ao mesmo procedimento de modelagem do grupo modelo e, posteriormente, recebeu gavage oral de celecoxib (20 mg/kg/dia), preparado em suspensão em carboximetilcelulose de sódio 0,5% (CMC-Na), por 8 semanas consecutivas, com um volume de gavage de 10 mL/kg de peso corporal.
(5) YQHXD grupo de baixa dose [YQHXD(L)]: A modelagem foi realizada como no grupo modelo. Os ratos receberam YQHXD a 3,87 g/kg/dia por gavage durante 8 semanas, com volume de gavage de 10 mL/kg e concentração de trabalho de 0,387 g/mL.
(6) YQHXD grupo de dose média [YQHXD(M)]: A modelagem foi realizada como no grupo modelo. Os ratos receberam YQHXD a 7,74 g/kg/dia por gavage durante 8 semanas, com volume de gavage de 10 mL/kg e concentração de trabalho de 0,774 g/mL.
(7) YQHXD grupo de alta dose [YQHXD(H)]: A modelagem foi realizada como no grupo modelo. Os ratos receberam YQHXD a 15,48 g/kg/dia por gavage durante 8 semanas, com volume de gavage de 10 mL/kg e concentração de trabalho de 1,548 g/mL.
Mudanças patológicas
O tecido caudal da DIV do rato foi isolado cirurgicamente, e foram observadas alterações na DIV e nos tecidos musculares ao redor. O tecido coletado foi então processado em seções patológicas, com amostras preparadas em seções longitudinais para observação morfológica. O IVD vertebral caudal, fixado em formaldeído por 96 horas, foi rapidamente descalcificado usando uma solução descalcificante, desidratado por uma série de álcool graduado e eliminado com xileno. Após a incorporação na parafina, foram preparadas seções de 5 μm com parafina. Após a desparafinização com xileno e desidratação por uma série de gradientes de álcool, as seções foram coradas com hematoxilina e eosina. Depois, foram desidratados, limpos e montados. As amostras foram observadas sob um microscópio óptico para examinar mudanças morfológicas na estrutura do DIV caudal entre os grupos. Posteriormente, foi realizada a coloração tricromática de Masson. As seções foram coradas sequencialmente com solução de dicromato de potássio e hematoxilina por 3 minutos cada, seguidas de diferenciação com etanol ácido e depois coloridas com Vermelho de Picrosírio por 10 minutos. Após tratamento com ácido fosfotungstico por 10 minutos, as seções foram tingidas com azul de toluidina por 5 minutos. Após a desidratação, as seções foram limpas com xileno, montadas e observadas ao microscópio para examinar as alterações na DIV caudal.
Detecção ELISA de proteínas-alvo séricas
Oito semanas após a cirurgia, os ratos foram anestesiados com isoflurano inalado, e sangue da aorta abdominal foi coletado. As amostras periféricas de sangue foram mantidas em temperatura ambiente por 1 hora e centrifugadas a 1.400 × g por 10 minutos. O soro era separado e armazenado a −80 °C até o uso. Os níveis séricos de AKT1, STAT3 e Bcl-2 foram quantificados por meio de ensaio imunossorvente ligado a enzimas (ELISA). Os ensaios foram realizados de acordo com os protocolos dos fabricantes. A absorvância foi registrada em 450 nm, e as concentrações foram determinadas a partir de curvas padrão. Os intervalos de detecção para AKT1, STAT3 e BCL-2 foram de 0,625–40 ng/mL, 31,25–1000 pg/mL e 0,16–10 ng/mL, respectivamente, com sensibilidades de 0,1 ng/mL, 10 pg/mL e 0,1 ng/mL. Todas as amostras séricas foram analisadas sem diluição, e todos os resultados se enquadraram nas faixas correspondentes de curvas padrão. Todas as amostras e padrões foram analisados em duplicado, e a normalização foi realizada com base em volumes séricos iguais.
Detecção por PCR dos níveis de expressão gênica
Tecidos intervertebrais do disco foram dissecados das vértebras caudais de rato, e aproximadamente 20 mg foram coletados de cada amostra. Os tecidos foram pulverizados em nitrogênio líquido e lisados em 500 μL de tampão de lise não enzimática por 15 minutos. O RNA total foi então extraído, e sua concentração e pureza foram medidas com um espectrofotômetro. Amostras com valores A260/A 280 de 1,8–2,0 foram aceitas para análise, e as concentrações de RNA ajustadas para 500 ng/μL. A transcrição reversa foi realizada de acordo com o protocolo do fabricante. A PCR quantitativa (qPCR) foi realizada em um sistema de PCR em tempo real para medir a expressão de mRNA de Akt1, Bcl-2 e Stat3 . As sequências de primeiro são fornecidas na Tabela 1. O programa de PCR consistia em desnaturação inicial a 95 °C por 30 s, seguida por 40 ciclos de 95 °C por 5 s e 60 °C por 30 s. A análise da curva de fusão foi realizada após amplificação para verificar a especificidade do produto. O Gapdh foi usado como gene de referência, e todas as amostras foram analisadas em triplicado. As eficiências de amplificação dos primers foram de 90% a 110%. A expressão gênica relativa foi calculada usando o método 2−ΔΔCt .
Western blot para níveis de expressão proteica
Proteína total foi isolada dos tecidos caudais do disco intervertebral de rato de cada grupo. Os tecidos foram lisados em tampão RIPA suplementado com inibidores de protease. A concentração de proteína foi medida usando um kit de ensaio de proteínas BCA e ajustada para 2 μg/μL. Quantidades iguais de proteína (15 μg por amostra) foram separadas por 10% de SDS-PAGE e transferidas para membranas de difluoreto de polivinilideno (PVDF) a uma corrente constante de 300 mA por 90 minutos usando um sistema de transferência úmida. As membranas foram então bloqueadas com leite desnatado a 5% por 1 hora em temperatura ambiente e incubadas durante a noite a 4 °C com anticorpos primários contra AKT1 (1:5.000), STAT3 (1:6.000), Bcl-2 (1:1.500) e β-actina (1:4.000). Após a lavagem com TBST, as membranas foram incubadas com anticorpos secundários conjugados com peroxidase de raiz-forte (HRP) em uma diluição de 1:5.000 por 1 h em temperatura ambiente. As bandas proteicas foram detectadas usando um reagente de quimioluminescência aprimorada (ECL), e imagens foram obtidas com um sistema de imagem em gel sob condições idênticas de exposição. As intensidades das bandas foram medidas usando ImageJ e normalizadas para β-actina para determinar a expressão relativa da proteína. Cada experimento foi realizado independentemente em triplicado, com cada amostra obtida de um animal diferente.
Estatísticas
Os resultados são expressos como média ± desvio padrão (DS). As diferenças entre os grupos foram avaliadas por análise unilateral da variância (ANOVA), seguidas pelo teste de diferença de menor significância (LSD) post hoc. A significância estatística foi definida como P < 0,05.