Validação do modelo por citologia vaginal
O sucesso da OVX foi confirmado em todos os animais ovariectomizados, evidenciado pela perda do ciclo estral dentro de 5 dias após a cirurgia (Figura 1A), enquanto os animais submetidos a cirurgia simulada mantiveram um ciclo de 4 dias (Figura 1B). Após 12 semanas, a DXA mostrou menor DMO femoral no grupo modelo OVX em comparação com os controles simulados (0,214 ± 0,020 vs. 0,281 ± 0,019 mg/cm3; Bonferroni-ajustado p < 0,001), confirmando o estabelecimento do modelo de perda óssea induzida por OVX antes do tratamento.
Densidade mineral óssea
A fórmula completa de ZGP (I2A2) apresentou a maior recuperação de DMO entre as subfórmulas fatoriais e foi estatisticamente comparável a ambos os controles positivos (Tabela 2, Figura 2). Todas as subfórmulas fatoriais aumentaram a DMO em relação ao grupo modelo. O alto nutrição do Yin isoladamente (I3A1: 0,253 mg/cm3) e o alto reforço do Yang isoladamente (I1A3: 0,243 mg/cm3) permaneceram abaixo do ZGP completo (I2A2: 0,307 mg/cm3), o que apoia uma interação em vez de um efeito simples de dose.
Marcadores séricos de remodelação óssea
Em comparação com o grupo modelo, a fórmula completa de ZGP foi associada a um aumento de P1NP (+51%; IC 95%: 210,3–327,5 pg/mL; p <0,0001) e diminuição do CTX-1 (-16%; p <0,01), consistente com um equilíbrio formação-reabsorção melhorado (Figura 3A, B). Os braços com monoterapia de alta dose produziram valores absolutos mais elevados de P1NP em algumas comparações, mas o padrão fatorial apoiou uma atividade combinada otimizada em vez de uma resposta máxima à dose em um único braço. Os dados completos no nível dos grupos estão em Tabela 3.
ANOVA fatorial e interação sinérgica
A ANOVA fatorial de dois fatores nas nove células fatoriais identificou efeitos principais significativos dos fatores de nutrição do Yin e de reforço do Yang, bem como termos de interação I por A significativos em todos os desfechos primários (Tabela 4). O efeito de interação para P1NP (ηp2 = 0,74) sugere que os componentes de reforço do Yang amplificaram a atividade dos componentes de nutrição do Yin nessas condições experimentais.
Densidade óssea (validação do modelo)
A densidade óssea femoral foi menor no grupo Modelo (OVX) em comparação com os grupos Controle e Sham-operado, confirmando a indução bem-sucedida do fenótipo de perda óssea por OVX antes do tratamento, conforme mostrado na Figura 4.
Expressão de EphB4/ephrinB2/RhoA/ALP
As análises de RT-qPCR utilizaram amostras biológicas de cada grupo com reações técnicas em triplicata, e a densitometria de Western blot foi realizada em lisados biológicos independentes (Figura 5 e Figura 6). Em comparação com o grupo modelo, a fórmula completa de ZGP esteve associada a níveis mais altos de EphB4 (0,498 ± 0,009 vs. 0,048 ± 0,002; p < 0,0001), ephrinB2 (0,411 ± 0,036 vs. 0,048 ± 0,003; p < 0,0001) e ALP (0,842 ± 0,019 vs. 0,090 ± 0,003; p < 0,0001), com níveis mais baixos de RhoA (1,658 ± 0,028 vs. 7,332 ± 0,430; p < 0,0001). Essas descobertas moleculares são interpretadas como associações com marcadores de via, e não como prova de mediação causal.
Disponibilidade de Dados
Os arquivos de dados originais utilizados para gerar os resultados do presente estudo são fornecidos como um Arquivo Suplementar de Código acompanhando este artigo

Figura 1. Citologia vaginal representativa utilizada para verificar a ciclicidade estral após ovarioectomia (coloração com hematoxilina e eosina, 200x). (A) Esfregaço vaginal representativo de uma rata ovarioectomizada mostrando citologia persistente não cíclica. (B) Esfregaço vaginal representativo de uma rata submetida a cirurgia simulada demonstrando ciclicidade estral normal. Clique aqui para visualizar uma versão maior desta figura.

Figura 2. DMO femoral em todos os grupos. As barras de erro indicam o desvio padrão (DP). Os símbolos estatísticos indicam comparações ajustadas pelo método de Bonferroni em relação ao modelo ou à fórmula completa de ZGP, conforme definido na seção Análise estatística. Clique aqui para visualizar uma versão maior desta figura.

Figura 3. Marcadores séricos de remodelação óssea. (A) P1NP, um marcador de formação óssea. (B) CTX-1, um marcador de reabsorção óssea. As barras de erro indicam o desvio padrão (DP). Os símbolos estatísticos indicam comparações ajustadas pelo método de Bonferroni em relação ao modelo ou à fórmula completa de ZGP. Clique aqui para visualizar uma versão maior desta figura.

Figura 4. Densidade óssea femoral confirmando a perda óssea induzida pela OVX antes do tratamento. Gráfico de barras comparando a densidade óssea femoral (mg/cm3) entre os grupos Controle (não ovariectomizado, sem tratamento), Sham-operado e Modelo (OVX) no ponto temporal de validação do modelo pré-tratamento. As barras representam média ± DP (n = 10/grupo). *** = p < 0,001 versus Controle e versus Sham (ANOVA unifatorial com correção de Bonferroni). Clique aqui para visualizar uma versão maior desta figura.

Figura 5. Expressão gênica relativa de mRNA. Níveis de expressão de (A) EphB4 (B) ephrinB2 (C) RhoA (D) ALP utilizando o método 2^-ΔΔCt normalizado em relação a GAPDH. A qPCR utilizou reações em triplicata técnica para cada amostra biológica. As barras de erro indicam o desvio padrão (DP). Os símbolos estatísticos indicam comparações ajustadas pelo teste de Bonferroni. Clique aqui para visualizar uma versão maior desta figura.

Figura 6. Expressão proteica por Western blot. Os painéis mostram imunoblotagens representativas para EphB4, ephrinB2, RhoA, ALP e controle de carga GAPDH, selecionadas a partir de amostras biológicas independentes. A densitometria foi normalizada em relação ao GAPDH utilizando o ImageJ. Pesos moleculares observados: EphB4, 110 kDa; ephrinB2, 36 kDa; RhoA, 22 kDa; ALP, 39 kDa; GAPDH, 37 kDa (ver Tabela Suplementar 1). Clique aqui para visualizar uma versão maior desta figura.
| Fator nutritivo do Yin (I) |
| I₁ (0 g) | I₂ (27 g) | I₃ (54 g) |
| Fator de potencialização de Yang (A) | A₁ (0 g) | I₁A₁ (Modelo) | I₂A₁ (Yin) | I₃A₁ (Hi-Yin) |
| A₂ (8 g) | I₁A₂ (Yang) | I₂A₂ (ZGP†) | I₃A₂ (Hi-Yin+Yang) |
| A₃ (16 g) | I₁A₃ (Hi-Yang) | I₂A₃ (Yin+Hi-Yang) | I₃A₃ (Full-Hi) |
Tabela 1: Desenho fatorial 3x3 e grupos experimentais. A tabela apresenta os 13 grupos experimentais organizados em um desenho fatorial 3x3 de níveis de dose de Yin-nutriente (I) e Yang-estimulante (A), além de 4 grupos de referência. †I₂A₂ = fórmula completa da Pílula Zuogui na dose equivalente clínica. Conversão de dose utilizando o método da área da superfície corporal (rato : humano = 6,25).
| Grupo | Célula fatorial | DBM (mg/cm³, x̄ ± DP) | vs. Modelo |
| Controle em branco | — | 0.279 ± 0.030 | **** |
| Sham-operado | — | 0.281 ± 0.019 | **** |
| Modelo (OVX) | I₁A₁ | 0.214 ± 0.020 | — |
| Fortalecimento do Yang | I₁A₂ | 0.258 ± 0.011 | ** |
| Alto fortalecimento do Yang | I₁A₃ | 0.243 ± 0.021 | * |
| Nutrição do Yin | I₂A₁ | 0.262 ± 0.015 | *** |
| Alta nutrição do Yin | I₃A₁ | 0.253 ± 0.017 | ** |
| Yin + Alto-Yang | I₂A₃ | 0.240 ± 0.014 | * |
| Alto-Yin + Yang | I₃A₂ | 0.248 ± 0.014 | * |
| Totalmente alto | I₃A₃ | 0.257 ± 0.017 | ** |
| Fórmula completa ZGP | I₂A₂ | 0.307 ± 0.010 | **** |
| Controle positivo da MTC | — | 0.307 ± 0.029 | **** |
| Controle de medicina ocidental | — | 0.307 ± 0.009 | **** |
Tabela 2: DMO femoral em todos os grupos (n=10/grupo, média ± DP). * = p <0,05, ** = p <0,01, *** = p<0,001, **** = p<0,0001 versus modelo (correção de Bonferroni).
| Grupo | Célula | P1NP (pg/mL, x̄ ± DP) | CTX-1 (ng/mL, x̄ ± DP) |
| Controle em branco | — | 1388,8 ± 29,7****#### | 33,4 ± 6,8****#### |
| Sham-operado | — | 1268,7 ± 124,1****#### | 45,0 ± 6,7****#### |
| Modelo (OVX) | I₁A₁ | 526,4 ± 78,7 | 105,5 ± 9,9 |
| Fortalecedor do Yang | I₁A₂ | 845,1 ± 76,7**** | 84,2 ± 6,6**** |
| Alto fortalecedor do Yang | I₁A₃ | 1013,0 ± 73,7 ****## | 65,5 ± 6,6****#### |
| Nutritivo do Yin | I₂A₁ | 870,7 ± 63,6**** | 88,2 ± 7,6** |
| Alto nutritivo do Yin | I₃A₁ | 1068,0 ± 48,7****#### | 64,5 ± 4,0****#### |
| Yin + Alto Yang | I₂A₃ | 951,5 ± 146,5****# | 63,1 ± 6,0****#### |
| Alto Yin + Yang | I₃A₂ | 990,2 ± 116,8****## | 62,6 ± 4,4****#### |
| Total-Alto | I₃A₃ | 739,5 ± 93,5** | 97,1 ± 8,9 |
| Fórmula completa ZGP | I₂A₂ | 794,8 ± 66,0**** | 88,2 ± 6,0**** |
| Controle positivo MTC | — | 967,7 ± 49,9****# | 66,3 ± 6,1****#### |
| Controle medicina ocidental | — | 1026,6 ± 65,1****### | 68,8 ± 7,3****### |
Tabela 3: P1NP e CTX-1 séricos em todos os grupos (n=10/grupo, média ± DP). ** = p <0,01, **** = p <0,0001 versus modelo; # = p <0,05, ## = p <0,01, ### = p <0,001, #### = p <0,0001 versus ZGP (Bonferroni). IC 95% para ZGP P1NP versus modelo: 210,3–327,5 pg/mL.
| Resultado | Fator | df | F | Valor-p | η² parcial |
| P1NP (pg/mL) | Yin (I) | 2, 81 | 17,84 | <0,001 | 0,31 |
| Yang (A) | 2, 81 | 6,22 | 0,003 | 0,13 |
| I × A | 4, 81 | 57,21 | <0,001 | 0,74 |
| BMD (mg/cm³) | Yin (I) | 2, 81 | 31,90 | <0,0001 | 0,44 |
| Yang (A) | 2, 81 | 28,70 | <0,0001 | 0,41 |
| I × A | 4, 81 | 24,60 | <0,0001 | 0,55 |
| CTX-1 (ng/mL) | Yin (I) | 2, 81 | 14,22 | <0,001 | 0,26 |
| Yang (A) | 2, 81 | 5,88 | 0,004 | 0,13 |
| I × A | 4, 81 | 19,43 | <0,001 | 0,49 |
Tabela 4: ANOVA fatorial 3x3 com efeitos principais e termos de interação, e comparações post hoc corrigidas por Bonferroni. Eta quadrado parcial é apresentado como medida do tamanho do efeito. gl indica graus de liberdade.
Tabela Suplementar 1. Informações sobre os anticorpos utilizados na análise por Western blot. A tabela lista as proteínas-alvo, os pesos moleculares observados, fornecedores, números de catálogo, espécies hospedeiras e diluições de trabalho dos anticorpos primários e secundários utilizados neste estudo.Clique aqui para baixar este arquivo.
Arquivo de Código Suplementar 1. Dados experimentais brutos que apoiam o estudo. Este arquivo contém os arquivos de dados originais utilizados para gerar os resultados apresentados no manuscrito, incluindo medições individuais de DMO, dados brutos de ELISA, valores de Ct de RT-qPCR com amplificação e curvas de dissociação, imagens originais de Western blot e arquivos de análise de densitometria, imagens representativas de microscopia histológica e arquivos de saída gerados pelos instrumentos associados a esses experimentos.Clique aqui para baixar este arquivo.