Model validation by vaginal cytology
OVX success was confirmed in all ovariectomized animals, as evidenced by loss of estrous cycling within 5 days after surgery (Figure 1A), whereas sham-operated animals retained a 4-day cycle (Figure 1B). At 12 weeks, DXA showed lower femoral BMD in the OVX model group than in sham controls (0.214 ± 0.020 vs. 0.281 ± 0.019 mg/cm3; Bonferroni-adjusted p < 0.001), supporting the establishment of the OVX-induced bone-loss model prior to treatment.
Bone mineral density
The full ZGP formula (I2A2) showed the highest BMD recovery among factorial subformulations and was statistically comparable to both positive controls (Table 2, Figure 2). All factorial subformulations increased BMD relative to the model group. High Yin-nourishing alone (I3A1: 0.253 mg/cm3) and High Yang-enhancing alone (I1A3: 0.243 mg/cm3) remained lower than full ZGP (I2A2: 0.307 mg/cm3), supporting an interaction rather than a simple dose effect.
Serum bone turnover markers
Compared with the model group, full-formula ZGP was associated with increased P1NP (+51%; 95% CI: 210.3–327.5 pg/mL; p <0.0001) and decreased CTX-1 (-16%; p <0.01), consistent with improved formation-resorption balance (Figure 3A,B). High-dose monotherapy arms produced higher absolute P1NP values in some comparisons, but the factorial pattern supported optimized combined activity rather than maximum single-arm dose response. Complete group-level data are in Table 3.
Factorial ANOVA and synergistic interaction
Two-way factorial ANOVA of the nine factorial cells identified significant main effects of Yin-nourishing and Yang-enhancing factors and significant I by A interaction terms across primary endpoints (Table 4). The interaction effect for P1NP (ηp2 = 0.74) suggests that Yang-enhancing components amplified the activity of Yin-nourishing components under these experimental conditions.
Bone density (model validation)
Femoral bone density was lower in the Model (OVX) group than in the Blank and Sham-operated groups, confirming successful induction of the OVX bone-loss phenotype prior to treatment, as shown in Figure 4.
EphB4/ephrinB2/RhoA/ALP expression
RT-qPCR analyses used biological samples from each group with technical triplicate reactions, and Western blot densitometry was performed on independent biological lysates (Figure 5 and Figure 6). Compared with the model group, full-formula ZGP was associated with higher 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), and ALP (0.842 ± 0.019 vs. 0.090 ± 0.003; p < 0.0001), with lower RhoA (1.658 ± 0.028 vs. 7.332 ± 0.430; p < 0.0001). These molecular findings are interpreted as pathway-marker associations rather than proof of causal mediation.
Data Availability
The original data files used to generate the current study results are provided as a Supplementary Coding File accompanying this article

Figure 1. Representative vaginal cytology used to verify estrous cyclicity following ovariectomy (hematoxylin and eosin staining, 200x). (A) Representative vaginal smear from an ovariectomized rat showing persistent non-cycling cytology. (B) Representative vaginal smear from a sham-operated rat demonstrating normal estrous cyclicity. Please click here to view a larger version of this figure.

Figure 2. Femoral BMD across all groups. Error bars denote SD. Statistical symbols indicate Bonferroni-adjusted comparisons versus model or full-formula ZGP as defined in the Statistical analysis section. Please click here to view a larger version of this figure.

Figure 3. Serum bone turnover markers. (A) P1NP, a bone formation marker. (B) CTX-1, a bone resorption marker. Error bars denote SD. Statistical symbols indicate Bonferroni-adjusted comparisons relative to the model or the full-formula ZGP. Please click here to view a larger version of this figure.

Figure 4. Femoral bone density confirming OVX-induced bone loss before treatment. Bar graph comparing femoral bone density (mg/cm3) among the Blank (non-ovariectomized, untreated), Sham-operated, and Model (OVX) groups at the pre-treatment model-validation time point. Bars represent mean ± SD (n = 10/group). *** = p < 0.001 vs. Blank and vs. Sham (one-way ANOVA with Bonferroni correction). Please click here to view a larger version of this figure.

Figure 5. Relative mRNA gene expression. Expression levels of (A) EphB4 (B) ephrinB2 (C) RhoA (D) ALP using the 2^-ΔΔCt method normalized to GAPDH. qPCR used technical triplicate reactions for each biological sample. Error bars denote SD. Statistical symbols indicate Bonferroni-adjusted comparisons. Please click here to view a larger version of this figure.

Figure 6. Western blot protein expression. Panels show representative immunoblots for EphB4, ephrinB2, RhoA, ALP, and GAPDH loading control, selected from independent biological samples. Densitometry was normalized to GAPDH using ImageJ. Observed molecular weights: EphB4, 110 kDa; ephrinB2, 36 kDa; RhoA, 22 kDa; ALP, 39 kDa; GAPDH, 37 kDa (see Supplementary Table 1). Please click here to view a larger version of this figure.
| Yin-nourishing factor (I) |
| I₁ (0 g) | I₂ (27 g) | I₃ (54 g) |
| Yang-enhancing factor (A) | A₁ (0 g) | I₁A₁ (Model) | 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) |
Table 1: Three by three factorial design and experimental groups. The table presents the 13 experimental groups arranged in a 3-by-3 factorial design of Yin-nourishing (I) and Yang-enhancing (A) dose levels, plus 4 reference groups. †I₂A₂ = full Zuogui Pill formula at clinical-equivalent dose. Dose conversion using the body surface area method (rat : human = 6.25).
| Group | Factorial Cell | BMD (mg/cm³, x̄ ± SD) | vs. Model |
| Blank control | — | 0.279 ± 0.030 | **** |
| Sham-operated | — | 0.281 ± 0.019 | **** |
| Model (OVX) | I₁A₁ | 0.214 ± 0.020 | — |
| Yang-enhancing | I₁A₂ | 0.258 ± 0.011 | ** |
| High Yang-enhancing | I₁A₃ | 0.243 ± 0.021 | * |
| Yin-nourishing | I₂A₁ | 0.262 ± 0.015 | *** |
| High Yin-nourishing | I₃A₁ | 0.253 ± 0.017 | ** |
| Yin + Hi-Yang | I₂A₃ | 0.240 ± 0.014 | * |
| Hi-Yin + Yang | I₃A₂ | 0.248 ± 0.014 | * |
| Full-High | I₃A₃ | 0.257 ± 0.017 | ** |
| ZGP full formula | I₂A₂ | 0.307 ± 0.010 | **** |
| TCM positive ctrl | — | 0.307 ± 0.029 | **** |
| Western medicine ctrl | — | 0.307 ± 0.009 | **** |
Table 2: Femoral BMD across all groups (n=10/group, mean ± SD). * = p <0.05, ** = p <0.01, *** = p<0.001, **** = p<0.0001 vs. model (Bonferroni-corrected).
| Group | Cell | P1NP (pg/mL, x̄ ± SD) | CTX-1 (ng/mL, x̄ ± SD) |
| Blank control | — | 1388.8 ± 29.7****#### | 33.4 ± 6.8****#### |
| Sham-operated | — | 1268.7 ± 124.1****#### | 45.0 ± 6.7****#### |
| Model (OVX) | I₁A₁ | 526.4 ± 78.7 | 105.5 ± 9.9 |
| Yang-enhancing | I₁A₂ | 845.1 ± 76.7**** | 84.2 ± 6.6**** |
| High Yang-enhancing | I₁A₃ | 1013.0 ± 73.7 ****## | 65.5 ± 6.6****#### |
| Yin-nourishing | I₂A₁ | 870.7 ± 63.6**** | 88.2 ± 7.6** |
| High Yin-nourishing | I₃A₁ | 1068.0 ± 48.7****#### | 64.5 ± 4.0****#### |
| Yin + Hi-Yang | I₂A₃ | 951.5 ± 146.5****# | 63.1 ± 6.0****#### |
| Hi-Yin + Yang | I₃A₂ | 990.2 ± 116.8****## | 62.6 ± 4.4****#### |
| Full-High | I₃A₃ | 739.5 ± 93.5** | 97.1 ± 8.9 |
| ZGP full formula | I₂A₂ | 794.8 ± 66.0**** | 88.2 ± 6.0**** |
| TCM positive ctrl | — | 967.7 ± 49.9****# | 66.3 ± 6.1****#### |
| Western medicine ctrl | — | 1026.6 ± 65.1****### | 68.8 ± 7.3****### |
Table 3: Serum P1NP and CTX-1 across all groups (n=10/group, mean ± SD). ** = p <0.01, **** = p <0.0001 vs. model; # = p <0.05, ## = p <0.01, ### = p <0.001, #### = p <0.0001 vs. ZGP (Bonferroni). 95% CI for ZGP P1NP vs. model: 210.3–327.5 pg/mL.
| Outcome | Factor | df | F | p-value | Partial η² |
| 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 |
Table 4: Three by three factorial ANOVA main effects and interaction terms with Bonferroni-corrected post-hoc comparisons. Partial eta-squared is reported as the effect size. df indicates degrees of freedom.
Supplementary Table 1. Western blot antibody information. The table lists the target proteins, observed molecular weights, suppliers, catalog numbers, host species, and working dilutions of the primary and secondary antibodies used in this study.Please click here to download this file.
Supplementary Coding File 1. Raw experimental data supporting the study. This archive contains the original data files used to generate the results presented in the manuscript, including individual BMD measurements, ELISA raw data, RT-qPCR Ct values with amplification and melt curve outputs, original Western blot images and densitometry analysis files, representative histological microscopy images, and instrument-generated output files associated with these experiments.Please click here to download this file.