ABSTRACT
This study demonstrates the antioxidative stress potential of
Shanyao–Fuling–Yiyiren (SFY) decoction—a Chinese polyherbal formulation
derived from Si Fang decoction—by establishing a systematic framework
that integrates network pharmacology, molecular docking, in vitro
synergy assays, cellular experiments, and 3D printing. Despite its long
traditional use, the molecular and cellular mechanisms underlying its
antioxidative effects remain unclear, and its formulations are based
more on empirical methods than on systematic design. To fill this gap,
a fuzzy mathematical model was used to select the optimal polyherbal
combination. A central composite circumscribed design determined that a
Shanyao:Fuling:Yiyiren ratio of 2:2:1 maximized radical scavenging,
with a strong correlation (R ^2 = 0.9665) between antioxidant activity
and the combination index. Furthermore, network pharmacology, molecular
docking, and cell‐based assays jointly confirm the AKT1/GSK3β/HIF1α
pathway plays a crucial role in preventing the antioxidant effects of
SFY. Finally, the development of 3D printing SFY‐inks with the
optimized shape fidelity suggests promising applications for both
nutraceuticals and hepatocellular carcinoma diagnosis. Overall, the
results prove that 3D printing SFY‐based polyherbal formulation with
promising antioxidant potential and maximum synergism may indeed be a
potential source of preventing oxidant damages in pharmaceutical and
food industries.
Keywords: 3D printing food, AKT1/GSK3β/HIF1α pathway, fuzzy
mathematical model, network pharmacology, Shanyao–Fuling–Yiyiren
(SFY)‐based polyherbal formulation
__________________________________________________________________
SFY‐based polyherbal formulation was selected based on fuzzy
mathematical model. The best combinations with synergistic effect were
optimized and identified by central composite circumscribed design.
Network pharmacology and molecular docking‐based approach revealed that
the potential mechanisms of decoction against oxidant damage via
multicomponent, multitarget, and multipathway. Decoction presented
ideal antioxidative protection activity against H[2]O[2] damage.
Decoction‐based inks have great practical printability and 3D printing
performance for potential application in functional food production.
graphic file with name FSN3-13-e70349-g005.jpg
1. Introduction
About 95% of humans worldwide stay up late at night, but handling jet
lag can disrupt circadian rhythms, which is considered a stress state
related to many diseases, including metabolic disorders and liver
cancer. As a vital organ in charge of detoxification, the liver is very
sensitive to damages from physiological stress, resulting in
life‐threatening conditions (e.g., liver injury) with a high mortality
rate and poor prognosis. Among the many causes of liver injury,
uninterrupted production of reactive oxygen species (ROS) from both
endogenous and exogenous origins has received much attention, including
hydroxyl radical (OH˙), hydrogen peroxide (H[2]O[2]), hydroxyl ion
(OH^−), singlet oxygen (^1O[2]), superoxide anion (O[2]˙^−), and ozone
(O[3]) (Chen et al. [44]2020). Studies have demonstrated a connection
between prolonged exposure to ROS and liver diseases (Prieto and
Monsalve [45]2017). The elevated levels of ROS are intracellular
signaling factors that remarkably enhance hepatic stellate cell
activation and extracellular matrix generation amid liver injuries.
Additionally, ROS activates significant oncogenic pathways leading to
hepatocarcinogenesis, such as extracellular signal‐regulated kinase,
protein kinase B, c‐Jun N‐terminal kinase, hypoxia‐inducible factor
(HIF), microtubule‐related protein kinase, and intensified cellular DNA
mutations (Hammouda et al. [46]2020). Concurrently, concerns about the
adverse effects of synthetic antioxidants in foods have diverted more
attention to novel sources of natural antioxidants from the safety
aspect. Given the critical role of oxidative stress in liver disease
occurrence, the role of phospholipid peroxidation in liver injuries
shall be urgently explored. Against this backdrop, we hypothesize that
phospholipid peroxidation represents a key mechanistic link between
circadian disruption and stress‐related liver injuries.
Emerging evidence has described that traditional medicinal herbs
potentially contribute to preventing oxidative stress‐related chronic
diseases (Allison et al. [47]2025; Ashraf et al. [48]^2024). These
herbs have been applied conventionally for millennia by many cultures
as medicine, flavoring reagents, and even food preservatives, and are
basically cheap and available for poor populations. Researchers have
tested their antioxidant abilities and potential replacements of
synthetic additives in protecting food and cosmetic products from
oxidative damages. Specifically, Shanyao (Rhizoma Dioscoreae), Fuling
(Poria cocos (Schw.) Wolf.), and Yiyiren (Coicis Semen) are widely used
in traditional medicine, owing to their potent antioxidant properties.
Shanyao in the crude form has long been used as a spice, dietary
supplement, and a constituent of many traditional Asian medicines (Luo
et al. [49]2024). Shanyao contains curcumin, which has these
pharmacological abilities owing to its basic beneficial antioxidant,
anti‐inflammatory, antibacterial, and anticancer abilities (Alam
et al. [50]2024). Fuling has diverse pharmacological activities against
rheumatoid arthritis, type II diabetes, multiple sclerosis,
atherosclerosis, Alzheimer's disease, and other chronic diseases (Guo
et al. [51]2025). Free radicals, which are key stimuli for
carcinogenesis, can be inhibited by Fuling from modulating lipid
peroxidation of membranes or oxidative DNA harms. As for the
antioxidant role, Yiyiren is proved to effectively scavenge diverse
risky free radicals, including ROS, NO[2] radicals, O[2]˙^−, and OH˙
(Zhang et al. [52]2024).
Regardless of the abundant research, observations about the antioxidant
synergism in mixtures are still deficient. Moreover, recent studies
have raised significant controversy regarding the standardization of
such evaluations, reflecting a wider gap in understanding the
underlying mechanisms of antioxidant synergism (Cnudde et al. [53]2024;
Eawsakul and Bunluepuech [54]^2024; Shen et al. [55]^2025). This lack
of consensus has led to speculation and explanations in the literature,
many of which lack empirical verification. In view of antioxidant
synergism effects, the mixing of antioxidant extracts may induce a
synergism to generate a better antioxidant effect than the sum produced
by single extracts (additives) (Eawsakul and Bunluepuech [56]2024;
Kurnia et al. [57]^2025). Moreover, the ratio and type of individual
herb extracts are pivotal in deciding the chemo‐preventive potential or
antioxidant capacity in a health‐benefiting herb combination. Thus, the
regulation of antioxidant properties shall be investigated according to
the proportions of herbs in a mixture that can be used to develop
functional foods and pharmaceutical products at different contents. For
instance, the combination of Osmanthus fragrans flower extract with
four types of tea (Longjing, Tieguanyin, black, or Pu'er Tea) can
synergistically scavenge 2,2‐diphenyl‐1‐picrylhydrazyl free radicals
(Mao et al. [58]2017), demonstrating how specific combinations can
enhance antioxidant performance beyond individual components.
To our knowledge, though the antioxidant activities of Shanyao, Fuling,
and Yiyiren have been extensively reported, there is little research on
the bio‐effects of extracts combined with bioactive compounds, which
are believed to improve the antioxidant benefit of free radical
scavenging. Therefore, this study was aimed primarily to explore the
antioxidant interaction among Shanyao, Fuling, and Yiyiren at various
proportions using response surface methodology (RSM) and the
combination index. On the basis of synergism and optimization of the
decoction process, we simultaneously conducted network pharmacology and
molecular docking experiments to elucidate the mechanisms of the
Shanyao, Fuling, and Yiyiren compound (SFY) binding with AKT1, GSK3B,
TP53, HIF1A, and PTGS2‐related targets. Then whether SFY can mitigate
oxidative stress via the AKT1/GSK3β/HIF1α antioxidant system was
evaluated, aiming to provide insight into the intervention of
H[2]O[2]‐induced oxidative injuries. Finally, a decoction‐based broad
spectrum of edible inks was selected and developed for a 3D printing
ink in a balanced diet. Collectively, the present study provides a
molecular basis for using SFY as a promising antioxidant in the future.
2. Materials and Methods
2.1. Reagents and Materials
Shanyao, Fuling, and Yiyiren were obtained from a local market of
Yangzhou. 2,2′‐Azino‐bis (3‐ethylbenzothiazoline‐6‐sulphonic acid)
(ABTS) and 2,4,6‐tris (2‐pyridyl)‐1,3,5‐triazine (TPTZ) were made in
Aldrich‐Sigma (St. Louis, MO, USA). H[2]O[2] (30%), chloroform,
anhydrous ethanol, and isopropanol were bought from Sinopharm Chemical
Reagent Co. Ltd. (China). TRIzol was obtained from Thermo Fisher
Scientific (USA). A real‐time fluorescence quantitative PCR system
(qRT‐PCR) and a reverse transcriptional kit were purchased from
TransGen Biotech (China). Superoxide dismutase (SOD) and catalase (CAT)
activity detection kits were obtained from Beyotime Biotechnology and
Grace Biotechnology (both China), respectively. Dimethyl sulfoxide
(DMSO), ethylene diamine tetraacetic acid (EDTA), phosphate buffer
solution (PBS) phosphate dry powder, and
3‐(4,5‐dimethylthiazol‐2‐yl)‐2,5‐diphenyltetrazolium bromide (MTT) were
offered by Beijing Solarbio (China). All other agents were analytically
pure.
2.2. Optimization and Preparation of Formulations of Decoction
The preparations of medicinal and edible herbals were determined
according to the method from Guan et al. ([59]2021). Based on a
modified classic folk recipe, Sifang Tang, Shanyao, Fuling, and Yiyiren
were selected, cleaned, and air‐dried at 60°C overnight. Then the dried
materials were filtrated through an 80‐mesh sieve and kept in an
airtight container at ambient temperature until used. With the
solid/liquid ratio at 1:20, the three types of herbal powders were
randomly combined and thoroughly mixed in an RK103H ultrasonic bath
(BANDE‐LIN SONOREX, Germany) at 80 kHz. After 4 h of extraction, the
extracts were centrifuged at 4500 rpm for 10 min. The final sensory
test samples were obtained by mixing the three extracts.
2.3. Fuzzy Mathematical‐Based Sensory Evaluation
Sensory evaluation was carried out according to the procedure from Jaya
and Das ([60]2003). Briefly, 10 reportedly healthy, nonsmoking
potential evaluators from the staff and students of the College of
Tourism and Culinary Science, Yangzhou University were selected at a
60% success rate in triangle tests and enrolled to evaluate the color,
odor, taste, and texture of sensory test samples. The testing was
finished in the laboratory as per regulations in ASTM MNL‐26 (1996).
The panelists each sensorily analyzed 10 samples. The tested samples
were categorized into four distinct grades: excellent, good, fair, and
poor, which reflected their distinctiveness. The criteria for sensory
evaluation were listed in Table [61]S1, and the voting outcomes for
each criterion were recorded in Table [62]1.
TABLE 1.
Sensory evaluation indexes of polyherbal formulations.
Num. Sample Color Odor Taste Texture
V[1] V[2] V[3] V[4] V[1] V[2] V[3] V[4] V[1] V[2] V[3] V[4] V[1] V[2]
V[3] V[4]
T [1] Euryale ferox ; Lotus seed; Shanyao 4 2 4 0 3 4 2 1 3 5 1 1 2 6 2
0
T [2] Euryale ferox ; Lotus seed; Yiyiren 3 3 4 0 2 4 4 0 6 6 1 0 1 3 5
1
T [3] Euryale ferox ; Lotus seed; Fuling 0 4 3 3 0 1 4 5 1 5 3 1 5 3 2
0
T [4] Lotus seed; Shanyao; Yiyiren 5 2 2 1 4 2 2 2 2 4 4 0 3 4 2 1
T [5] Lotus seed; Shanyao; Fuling 4 5 1 0 5 3 1 1 5 5 0 0 4 4 2 0
T [6] Shanyao; Yiyiren; Fuling 6 4 0 0 7 2 1 0 7 3 0 0 3 4 1 2
T [7] Euryale ferox ; Shanyao; Yiyiren 4 4 2 0 6 3 1 0 5 4 1 0 2 8 0 0
T [8] Euryale ferox ; Shanyao; Fuling 2 1 4 3 3 4 2 1 3 7 0 0 0 6 4 0
T [9] Euryale ferox ; Yiyiren; Fuling 3 5 1 1 2 8 0 0 7 1 1 1 5 3 1 1
T [10] Lotus seed; Yiyiren; Fuling 1 6 1 2 5 3 2 0 4 4 2 0 2 5 2 1
[63]Open in a new tab
Fuzzy mathematical sensory evaluation was integrated into food sensory
assessment to quantify evaluation factors, substantially mitigate the
influence of personal bias, and produce more accurate and scientific
scoring results (Ranneh et al. [64]2021). In detail, a set of sensory
factors for the prepared solutions was defined as U, and a set of
grades as V. U comprises color (U[1]), odor (U[2]), taste (U[3]), and
texture (U[4]), and V includes excellent (V[1]), good (V[2]), average
(V[3]), and poor (V[4]). The relative weight of each sensory factor in
the overall flavor profile was determined using the fuzzy binary
comparative decision method. Ten sensory evaluators compared these
factors pairwise. A factor deemed important was given 1 point, whereas
those considered less important were assigned 0 point. The total score
of each sensory factor divided by the maximum score of 100 was used as
the weight of this sensory factor (Table [65]S2). Based on the
evaluators' scoring on the importance of color, odor, taste, and
texture in the overall sensory evaluation, the distribution of weights
among these factors is X = {color, odor, taste, texture} = {0.24, 0.28,
0.33, 0.15}.
2.4. Determination of ABTS Radical (ABTS^+) Scavenging Capacity
The ABTS^+ scavenging ability was tested following the method from
Guan, Li, et al. ([66]2024) with modifications. Briefly, an ABTS^+
stock solution was prepared by mixing 5 mL of 7 mM ABTS and 5 mL of
2.45 mM potassium persulfate and put for 12–16 h at room temperature
(RT) without light. Then the mixture was diluted using ultrapure water
until the absorbance at 405 nm was 1.4, forming an ABTS working
solution. Gradient diluted extracts or ascorbic acid (0.5 mL) were
blended with 0.5 mL of the ABTS working solution and stood for 30 min
in the dark at RT. Then the absorbance at 734 nm was recorded. The
percent of ABTS^+ scavenging activity of the extracts was computed
using Equation ([67]1):
[MATH: ABTS%=1−A1−A2/A0×100% :MATH]
(1)
where A [1], A [2], and A [0] are the absorbance of the mixture of
ABTS^+ and the sample solution, the mixture of ABTS^+ and the control
sample, and the mixture of ABTS^+ and deionized water, respectively.
2.5. Experiment Design and Mixture Optimization
Mixture design, a special type of RSM, was used to optimize the
composition of herbal mixtures and test the interactive effect between
components. Based on Section [68]2.3, the polyherbal formulations for
further optimization consisted of Shanyao, Fuling, and Yiyiren. In
total, 16 formulations were formed and the responses were analyzed
using DesignExpert 13 (Minitab Inc., State College, PA, USA). The
layout of the herbal formulations is shown in Table [69]2. The
dependent variable was in vitro ABTS antioxidant ability. The canonical
model of was used for each response after adjustment based on the
testing data. Linear, quadratic, and special cubic models were tested
to determine regression coefficients, which were kept only at the
significant level. Data were refitted to obtain the final model for
each index. The adequacy and goodness of fitting for each model were
statistically analyzed and fitted to a second‐order polynomial
regression model involving the coefficients of linear, quadratic, and
interactive terms. For validation, the optimal formulation was examined
in triplicate and expressed as mean ± standard deviation. Analysis of
variance (ANOVA) was conducted to calculate the model significance and
suitability of the factors and interactions. The formulations with
highly desirable functions were chosen for further analysis.
TABLE 2.
Mixture design experimental arrangement and results.
Standard order Factor A (Fuling, g) Factor B (Shanyao, g) Factor C
(Yiyiren, g) Response
ABTS^+ radical scavenging capacity (%)
1 0.75 0.55 0.40 41.30 ± 0.73
2 0.65 0.75 0.30 56.98 ± 0.41
3 0.65 0.65 0.40 54.35 ± 0.98
4 0.75 0.62 0.33 50.76 ± 0.47
5 0.65 0.65 0.40 61.18 ± 0.63
6 0.75 0.65 0.30 55.61 ± 0.24
7 0.65 0.75 0.30 53.13 ± 0.30
8 0.72 0.62 0.37 61.32 ± 1.09
9 0.65 0.75 0.30 54.63 ± 0.38
10 0.75 0.75 0.20 46.27 ± 0.29
11 0.55 0.75 0.40 61.15 ± 0.66
12 0.62 0.72 0.37 67.47 ± 0.68
13 0.55 0.75 0.40 57.46 ± 0.34
14 0.75 0.68 0.27 53.50 ± 0.64
15 0.68 0.68 0.33 67.60 ± 0.60
16 0.65 0.65 0.40 54.94 ± 0.21
[70]Open in a new tab
2.6. Determination of Antioxidant Synergism
Based on the Chou–Talalay combined drug theory, the synergistic
antioxidant effects of different combinations were explored
(Chou [71]2018). Specifically, antioxidant synergism was determined as
per the sum of IC[50] from seven concentrations based on RSM to
generate a total of 16 solution combinations. The absorbance of the
solution combination was then detected to determine the percent of
ABTS^+ scavenging capacity. To quantify the synergistic, additive, or
antagonistic impact of the combinations of decoction, the testing data
were converted to the combination index (CI):
[MATH: CI=MCa/SCa+MCb/SCb :MATH]
where MCa and MCb are the concentrations of compounds A and B in the
mixture to achieve 50% antioxidant activity respectively; SCa and SCb
are the EC[50] of the single compounds A and B, respectively. CI < 1,
= 1, or > 1 implies a synergistic, additive, or antagonistic effect,
respectively.
2.7. Identification of Core Ingredients, Potential Targets of Decoction, and
Network Construction
The active constituents and related targets in SFY‐based polyherbal
formulations were cited from Chinese Medicine System Pharmacology
Database ([72]https://old.tcmsp‐e.com/tcmsp.php) with pharmacokinetic
parameters (oral bioavailability [OB] ≥ 30% and drug‐like activity
[DL] ≥ 0.18). Then the target names of the corresponding core active
ingredients were turned into a unified format with Homo sapiens via
Unified Protein Database ([73]https://www.uniprot.org/). In addition,
the potential antioxidant associated targets were gathered from
GeneCards ([74]http://www.genecards.org/), Herb 2.0
([75]http://www.disgenet.org/search), and OMIM ([76]https://omim.org/).
Then duplicates were removed after the targets were merged from the
databases above. The targets of antioxidants and SFY‐based polyherbal
formulations and their common targets were acquired via Venny 2.1.0. To
systematically examine the network, the “herb‐compound‐common target”
network was imported and constructed via Cytoscape 3.9.1.
2.8. Protein–Protein Interaction (PPI) Network Construction, Topology
Analysis, and GO/KEGG Enrichment Analysis of Core Targets
To investigate the core targets and their interactions against oxidant
damages, the overlapping targets between SFY‐based polyherbal
formulations and oxidant damages were examined using the STRING network
platform ([77]https://stringdb.org/). Furthermore, a PPI network was
created at a confidence level > 0.4 and with the “ Homo sapiens ”
filter (Ren et al. [78]2024). The network was further analyzed in
Cytoscape 3.9.1 using MCODE to identify the associations between hub
genes and network clusters with specific parameters. Core clusters were
then identified based on degree centrality. GO and KEGG pathway
enrichment analyses of the core targets were done using David 6.8
([79]https://david.ncifcrf.gov/) to explore specific mechanisms and
signaling pathways. Enrichment results with p and q < 0.05 were
visualized using bar and bubble charts on the online tool imageGP.
2.9. Molecular Docking
Molecular docking, an increasingly important computational tool for
exploring the behaviors of biomacromolecule complexes, was conducted as
per the approach of Guan et al. ([80]2023) with modifications. Based on
the degrees from Section [81]2.8, the top targets and core ingredients
from the SFY‐based polyherbal formulations were chosen for affinity
calculations. The 3D structures of AKT1, GSK3B, TP53, HIF1A, and PTGS2
were cited and constructed via AutoDockTools 1.5.6. Polar hydrogens
were added to demand charges after water was removed. The smallest
binding energy was computed on AutoDock Vina. Finally, the binding
details were illustrated on PyMOL and Discovery Studio 2018.
2.10. Cytotoxicity Evaluation and Antioxidant Enzyme Activity Analysis
The BNLCL.2 mouse embryonic liver cell line from Meilun Biotechnology
Co. Ltd. (Dalian, China) was cultured in high‐glucose DMEM added with
10% FBS and 1% P/S at 37°C with 5% CO[2]. Trypsin was used for cell
passage. For H[2]O[2] cytotoxicity assessment, the BNLCL.2 cells
(1 × 10^4 cells/well) were planted in 96‐well plates and processed with
0, 200, 400, 600, 800, or 1000 μmol/L H[2]O[2] for 24 h. To evaluate
the protective effect of SFY, the cells were pretreated with 600 μM
H[2]O[2] for 4 h, and co‐cultured with 0.5, 1, 2, 4, or 8 mg/mL SFY for
24 h. In MTT assays, each well was added with 100 μL of 5 mg/mL MTT,
incubated for 4 h, and then the formazan crystals were solubilized with
DMSO. The absorbance at 490 nm was recorded using a microplate reader
(Flefel et al. [82]2019). Additionally, BNLCL.2 cells were pretreated
with H[2]O[2] or SFY at specific concentrations. Then CAT and SOD
activities were assessed using commercial kits (Cat. no. S0101S,
G0105W48) as instructed by the manufacturer, and were detected with
colorimetry and the xanthine oxidase method, respectively (Ding
et al. [83]2020).
2.11. RNA Extraction and qRT‐PCR Analysis
The BNLCL.2 cells were incubated in 6‐well plates for 24 h, pretreated
with 600 μM H[2]O[2] for 4 h, and exposed to SFY seeds (0, 2, 4, and
8 mg/mL) for 24 h. The treated BNLCL.2 cells were collected. Total RNA
was extracted from the BNLCL.2 cells using a TRIzol reagent, and
reverse transcribed to cDNA using a first‐strand cDNA synthesis
SuperMix kit for qPCR (gDNA digester plus) from TransGen Biotech
(China). Real‐time PCR was conducted with fast SYBR green master mix
(Wang et al. [84]2023). The quantities of transcripts were standardized
to that of GAPDH. The primers (Sangon Biotech, China) were presented in
Table [85]S3.
2.12. 3D Printing Ink Preparation and Printability Assessment
3D printing inks were prepared and assessed in terms of printability by
forming different structures and printing them using an in‐house 3D
printing system. The SFY decoction‐based 3D‐printed gels were prepared
following a modified version of a previous method. Specifically, 2% and
2.7% high acyl gellan gum (HAG) and gelatin (GL) hydrogels were made by
dissolving the powder in a solution at about 50°C under continuous
stirring. These hydrocolloids were then mixed with 1%, 2%, 3%, or 4%
glycerin to form the final testing gels. An extrusion‐based 3D printer
(Luckybot One, Wiiboox Technology Co. Ltd., Nanjing, China) with a
syringe and a retractable plunger was operated at around 26°C. A
0.25 mm inner diameter nozzle was used for printing, and images were
captured from the top side on a food‐grade glass plate using a mobile
phone. Uniform lighting was achieved with a white mini‐studio light box
(Rui Teng Digital, Zhejiang, China) containing 144 LED aerial lights
above the inks.
2.13. Statistical Analysis
All assays were conducted in biological triplicate unless specified
otherwise. The means of two independent groups were compared with an
unpaired two‐tailed Student's t‐test. The normally distributed data
among more than two groups were compared via one‐way ANOVA with Tukey's
post hoc test. Data were expressed as mean ± standard deviation (SD)
and analyzed on GraphPad Prism 9.5.0 (San Diego, CA, USA). The
significance levels were *p < 0.05, **p < 0.01, and ***p < 0.001.
3. Results and Discussion
3.1. Sensory and Fuzzy Comprehensive Evaluation of Optimized Compatibility
Components
To enhance the medicinal value and reduce production costs of Sifang
Tang, sensory evaluation and fuzzy mathematical methods were utilized
to optimize the active ingredients and compositions. The compatibility
of the five components (Qianshi, Shanyao, Fuling, Lianzi, and Yiyiren)
in Sifang Tang was analyzed through sensory evaluation (Table [86]1)
and radar plotting (Figure [87]1A–D). In detail, after the individual
sensory attribute for each polyherbal formulation was calculated, the
highest color was found on sample 6 (V[1] = 6), followed by sample 4
(V[1] = 5), and samples 1, 5, and 7 (V[1] = 4). As for odor, the
highest polyherbal formulation on the response scale was found in
sample 6, which linguistically means “Good.” The highest SI on taste,
flavor, and mouthfeel was found all on the response scale F4, which
linguistically implies “Good.” The SI of color was on the response
scale F3, which linguistically refers to “Satisfactory.” For sample 6,
the frequency of odor under “Excellent,” “Good,” “Average,” and “Poor”
is 7, 2, 1, and 0, and the frequency of odor is 7, 3, 0, and 0
respectively. From the above ranking, in the EFSSE fortified bread
sample, sample 6 (Shanyao, Fuling, and Yiyiren) showed better sensory
evaluation compatibility.
FIGURE 1.
FIGURE 1
[88]Open in a new tab
Optimized compound formula compatibility: sensory evaluation (A–D),
antioxidant effect (E–G), and component quantification analysis (H–J).
Sums of sensory scores for quality attributes of tested samples: T [1]
( Euryale ferox , Lotus seed, Shanyao), T [2] ( Euryale ferox , Lotus
seed, Yiyiren), T [3] ( Euryale ferox , Lotus seed, Fuling), T [4]
(Lotus seed, Shanyao, Yiyiren), T [5] (Lotus seed, Shanyao, Fuling), T
[6] (Shanyao, Yiyiren, Fuling), T [7] ( Euryale ferox , Shanyao,
Yiyiren), T [8] ( Euryale ferox , Shanyao, Fuling), T [9] ( Euryale
ferox , Yiyiren, Fuling), T [10] (Lotus seed, Yiyiren, Fuling).
Traditional sensory evaluation methods are frequently hard to achieve
consensus, owing to personal subjective perceptions, environmental
variations, and psychological fluctuations. In comparison, membership
function theory‐based fuzzy mathematical sensory evaluation can enhance
the reliability of sensory evaluations, and thus has been innovatively
integrated into food sensory assessment (Pallavi [89]2025). The sensory
evaluation data from the 10 evaluators (Table [90]1) were aggregated
into matrix R by dividing the number of votes received for each grade
by 10. The weight set X was then integrated with matrix R to form an
evaluation matrix Y = X × R. Then matrix Y was processed using a
comprehensive scoring matrix T. The set of evaluation grades K = {90,
70, 50, 30} was utilized, and each grade was multiplied by its
corresponding weight and summed to compute the total fuzzy
comprehensive evaluation score for each sample. For instance, when the
colors of medicinal solutions in the first group were evaluated: four
evaluators rated it as excellent, two as good, four4 as average, and
none as poor, resulting in R [color] = (0.4, 0.2, 0.4, 0), R
[odor] = (0.3, 0.4, 0.2, 0.1), R [taste] = (0.3, 0.5, 0.1, 0.1), and R
[texture] = (0.2, 0.6, 0.2, 0). Then we have
[MATH: Y1=R1×X=0.40.30.20.40.4
0.200.1
0.30.2<
/mtable>0.50.60.10.2<
mtd>0.10×0.24,0.28,0.33
mn>,0.15=0.309,0.415,0.21
5,0.061 :MATH]
Similarly, we have Y [2] = (0.341, 0.33, 0.427, 0.015), Y [3] = (0.108,
0.446, 0.313, 0.133), Y [4] = (0.343, 0.296, 0.266, 0.095), Y
[5] = (0.461, 0.429, 0.082, 0.028), Y [6] = (0.616, 0.311, 0.043,
0.03), Y [7] = (0.459, 0.432, 0.109, 0), Y [8] = (0.231, 0.457, 0.212,
0.1), Y [9] = (0.434, 0.422, 0.072, 0.072), and Y [10] = (0.326, 0.435,
0.176, 0.072). Then, we get T
[1] = 0.309 × 90 + 0.415 × 70 + 0.215 × 50 + 0.061 × 30 = 69.44, T
[2] = 75.59, T [3] = 60.58, T [4] = 65.74, T [5] = 72.41, T
[6] = 80.26, T [7] = 77, T [8] = 66.38, T [9] = 74.36, and T
[10] = 70.75.
With sensory evaluation of fuzzy mathematical models considered, the
polyherbal formulations can be ranked as follows: T [6] > T [7] > T
[2] > T [9] > T [5] > T [10] > T [1] > T [8] > T [4] > T [3]. After
quantitative analysis of sensory evaluation using fuzzy mathematics,
the group members/consumers believed sample T [6] (Shanyao, Fuling, and
Yiyiren) had higher overall sensory acceptability than the average.
Moreover, the overall acceptability parameter was mainly dependent on
quantitative analysis of personalized sensory attributes of the
samples, including taste, flavor, and texture. These findings indicate
the sensory characteristics of the SFY decoction are closely linked to
the proportions of Shanyao, Fuling, and Yiyiren, which could be
explained by the different molecular profiles of each herb that
contribute uniquely to taste and texture. However, the sample scope of
the current study was limited to school‐based participants and thus may
not fully reflect broader market preferences. Future work will focus on