Phytochemical, In Vitro, In Vivo, and In Silico Research on the Extract of Ajuga chamaepitys (L.) Schreb.

The present study focuses on the chemical characterization of a dry extract obtained from the species Ajuga chamaepitys (L.) Schreb, evaluating its antioxidant properties, toxicity, and in silico profile. Quantitative analysis of the dry extract revealed a notable amount of phytochemical compounds: 59.932 ± 21.167 mg rutin equivalents (mg REs)/g dry weight, 45.864 ± 4.434 mg chlorogenic acid equivalents (mg ChAEs)/g dry weight and, respectively, 83.307 ± 3.989 mg tannic acid equivalents (TAEs)/g dry weight. By UHPLC-HRMS/MS, the following were quantified as major compounds: caffeic acid (3253.8 μg/g extract) and kaempherol (3041.5 μg/g extract); more than 11 types of polyphenolic compounds were quantified (genistin 730.2 μg/g extract, naringenin 395 μg/g extract, apigenin 325.7 μg/g extract, galangin 283.3 μg/g extract, ferulic acid 254.3 μg/g extract, p-coumaric acid 198.2 μg/g extract, rutin 110.6 μg/g extract, chrysin 90.22 μg/g extract, syringic acid 84.2 μg/g extract, pinocembrin 32.7 μg/g extract, ellagic acid 18.2 μg/g extract). The antioxidant activity was in accordance with the amount of phytochemical compounds: IC50DPPH = 483.6 ± 41.4 µg/mL, IC50ABTS•+ = 127.4 ± 20.2 µg/mL, and EC50FRAP = 491.6 ± 2 µg/mL. On the larvae of Artemia sp., it was found that the extract has a low cytotoxic action. In silico studies have highlighted the possibility of inhibiting the activity of protein kinases CDK5 and GSK-3b for apigenin, galangin, and kaempferol, with possible utility for treating neurodegenerative pathologies and neuropathic pain. Further studies are warranted to confirm the predicted molecular mechanisms of action and to further investigate the therapeutic potential in animal models of neurological disorders.

Ajuga chamaepitys (L.) Schreb., called Ground Pine, is an annual plant that belongs to the Mediterranean area and is used for several effects in traditional medicine.Infusions prepared from its aerial parts are used to treat diarrhea, hemorrhoids, and various intestinal infections by drinking 3-4 cups of tea a day.Also, the plant has been found to have insecticide effects [3].
Among the aforementioned phytochemicals, neo-clerodane diterpenoids are the characteristic and dominant constituents of the Ajuga species.In addition, iridoids (e.g., ajugol and ajugoside) are chemotaxonomic markers of the genus Ajuga [15].
In Romania, the species Ajuga chamaepitys (L.) Schreb. is dominant in the Dobrogea region (Dobrogea Gorges-44 • 29 ′ 16.5 ′′ N-28 • 27 ′ 09.5 ′′ E), where it is spread spontaneously in many areas, and is traditionally used in this region for analgesic and anti-inflammatory effects.We mention that the species that grows in this region of Romania has not been studied, and our research aims to form the basis of a larger study regarding the scientific foundation of therapeutic effects, correlated with the type of active chemical constituents.
In our previous study [16], the content of phytochemical compounds (flavonoids, phenolic acids, and total polyphenols) was analyzed in the aerial parts, flowers, leaves, roots, and stems collected in four flowering periods (May, June, July, August) using 70% ethanol (v/v).We found that plant materials collected in May contained the highest content in phytochemicals.Further analysis comparing extraction solvents-50% ethanol (v/v), 20% ethanol (v/v), and distilled water-and correlated with the results obtained on 70% ethanol extraction revealed that 50% ethanol (v/v) is the optimal solvent for extracting these phytochemical constituents, and the plant material collected in May had the most abundant content of phytochemicals.Although flowers are a plant product with an increased content of phytoconstituents, a small quantity of raw material is generated.We decided to continue the analyses on the aerial parts since they contain similar compounds, aiming to obtain a dry extract with potential therapeutic value.Extractive solutions in 50% ethanol, 20% ethanol, and water fall under the trend of implementing green technology and the use of green solvents [17,18].
The aim of this research was to evaluate the quality of an extract of Ajuga chamaepitys (L.) Schreb (spectrophotometric and UHPLC-HRMS/MS-ultra-performance liquid chromatography coupled with high-performance mass spectrometry resolution-analysis of polyphenols) obtained from the aerial parts of the species harvested from the spontaneous flora of Dobrogea Gorges, Romania.(44 • 29 ′ 16.5 ′′ N-28 • 27 ′ 09.5 ′′ E).For the projection of the following research, the extract was also subjected to in vitro evaluations (antioxidant assay by three methods: DPPH, ABTS •+ , FRAP), in vivo evaluations (test of cytotoxicity on Artemia salina (L.) and Daphnia species and the teratogenic risk on Daphnia magna embryos), and in silico determinations (molecular docking studies for the potential of inhibiting the activity of protein kinases CDK5 and GSK-3b for apigenin, galangin, and kaempferol, with possible utility for treating neurodegenerative pathologies and neuropathic pain).

Establishing the Identity and Quality of the Plant Raw Material
Ajugae chamaepitysis herba (ACH) (Figure 1), called Ground Pine, is an annual plant that blooms in summer from May to August.It presents a thin, quadrangular, 10-30 cm long stem, its leaves are 1-4 cm long, with flowers which are up to 15 mm long.The corolla ranges from yellow to pale red, and it is 10-15 mm long, whereas its fruit appearance is a tetranucule, and the whole plant emits a pine scent when crushed [7,11,19].
long stem, its leaves are 1-4 cm long, with flowers which are up to 15 mm long ranges from yellow to pale red, and it is 10-15 mm long, whereas its fruit ap tetranucule, and the whole plant emits a pine scent when crushed [7,11,19].The results for flavonoids, phenolic acids, and total polyphenols conten tive solutions are shown in Table 1.Notes.Results are presented as the mean ± SD (standard deviation) ( = 3).Legend: H samples), FC = flavonoids content; PAC = phenolic acids content; TPC = total phenol each sample, corresponding index to the solvent extracted was added, 50% ethanol ethanol (v/v)-b and distilled water-c.
The results obtained guide the type of ethanol used (ethanol 50% (v/v)).

Obtaining and Characterizing the Dry Extract
A quantity of 52 g of dry extract Ajugae chamaepitysis herbae extractum obtained with ethanol 50% (v/v), which is presented in the form of a green powder with a characteristic persistent smell.The extraction yielded 28.45% was kept in a desiccator).

Spectrophotometric Determination of Flavonoids, Phenolic Acids, and phenols
The results of the quantitative analysis of flavonoids, phenolic acids, an phenols for the selective dry extract are shown in Figure 2. The results for flavonoids, phenolic acids, and total polyphenols content for extractive solutions are shown in Table 1.Notes.Results are presented as the mean ± SD (standard deviation) (n = 3).Legend: H (aerial parts samples), FC = flavonoids content; PAC = phenolic acids content; TPC = total phenolic content.For each sample, corresponding index to the solvent extracted was added, 50% ethanol (v/v)-a, 20% ethanol (v/v)-b and distilled water-c.
The results obtained guide the type of ethanol used (ethanol 50% (v/v)).

Obtaining and Characterizing the Dry Extract
A quantity of 52 g of dry extract Ajugae chamaepitysis herbae extractum (ACHE) was obtained with ethanol 50% (v/v), which is presented in the form of a green hygroscopic powder with a characteristic persistent smell.The extraction yielded 28.45% (the extract was kept in a desiccator).

Spectrophotometric Determination of Flavonoids, Phenolic Acids, and Total Polyphenols
The results of the quantitative analysis of flavonoids, phenolic acids, and total polyphenols for the selective dry extract are shown in Figure 2.

Evaluation of the Scavenger Capacity of the ABTS •+ Radical
To evaluate the scavenger capacity of the ABTS •+ radical, the equations of the calibration curves, inhibition (%) vs. concentration (µg/mL) for the dry extract (Figure S6), for a reaction time of 6 min, are frequently used (n = 3).It was observed that at a concentration of 409.6 µg/mL, the inhibition was over 80%.The IC 50 values were calculated at 6 min (I) (Figure S6a), showing a value of 104 µg/mL, then 138.8 µg/mL at 6 min (II) (Figure S6b) and 139.3 µg/mL at 6 min (III) (Figure S6c), having an average of 127.4 ± 21.2 µg/mL.
In Figure 5, antioxidant activity is statistically represented for the three methods used for ACHE.The results obtained by BSLA indicate significant effects of the extract at the tested concentration of 2600 µg/mL (LC50).The first records of mortality were from 557.85 µg/mL (LC10).The statistically analysis (Table 4) indicate statistically significant results and an important concentration-effects correlation (r 2 = 0.9998, Supplementary Materials Figures The results obtained by BSLA indicate significant effects of the extract at the tested concentration of 2600 µg/mL (LC 50 ).The first records of mortality were from 557.85 µg/mL (LC 10 ).The statistically analysis (Table 4) indicate statistically significant results and an important concentration-effects correlation (r 2 = 0.9998, Supplementary Materials Figures S9 and S11).LC 50 mortality recorded at extract concentrations above 1000 µg/mL indicates a low cytotoxic effect.In the control samples, there was no larval mortality during the observation period.Observations allowed the identification of reduced, abnormal movements of larvae exposed to C 7 -C 11 concentrations.These changes may be associated with the interaction of the extract (compounds) with cells that are not yet differentiated but play a role in contraction (myocytes).

Daphnia Species Toxicity Assay
The extract was virtually non-toxic to Daphnia magna throughout the exposure period, with the highest recorded average lethality being 10%.On Daphnia pulex, the extract exhibited toxicity at high concentrations, indicating a greater sensitivity of this species to the extract.While the maximum lethality was 20% at 24 h of exposure, it ranged from 10% at the lowest concentration to 95% at the highest concentration at 48 h.Furthermore, the effect was concentration-dependent (r 2 = 0.86).The LC 50 at 48 h was 584.10 µg/mL, and the 95% confidence interval (CI 95%) was between 470.6 and 724.8 µg/mL, indicating low to moderate toxicity.The lethality curve is presented in Figure 6.The mobility of the crustaceans was significantly reduced at concentrations between 500 and 1000 µg/mL.the extract.While the maximum lethality was 20% at 24 h of exposure, it ranged from 10% at the lowest concentration to 95% at the highest concentration at 48 h.Furthermore, the effect was concentration-dependent (r 2 = 0.86).The LC50 at 48 h was 584.10 µg/mL, and the 95% confidence interval (CI 95%) was between 470.6 and 724.8 µg/mL, indicating low to moderate toxicity.The lethality curve is presented in Figure 6.The mobility of the crustaceans was significantly reduced at concentrations between 500 and 1000 µg/mL.At a concentration of 50 µg/mL, the developmental stage and viability of the embryos were similar to those of the untreated control.Thus, the formation of the compound eye, the antennae, the post-abdominal claw, and the shape of the rostrum and dorsal outline At a concentration of 50 µg/mL, the developmental stage and viability of the embryos were similar to those of the untreated control.Thus, the formation of the compound eye, the antennae, the post-abdominal claw, and the shape of the rostrum and dorsal outline of the carapaces were similar to those observed in the untreated control.At 500 µg/mL, the viability was slightly lower than the control, with no abnormalities recorded in development, but there was a delay in the developmental stage.At the highest tested concentration of 1000 µg/mL, the extract affected the viability of the embryos, with the number of undeveloped embryos being three times higher than the control.The developed larvae exhibited reduced mobility, and the formation of the compound eye and antennae indicated a retardation effect (Figure 7).The toxicity assays on Daphnia species revealed that the extract exhibited low to moderate toxicity, with effects being concentration-dependent.While higher concentrations affected Daphnia pulex significantly and delayed Daphnia magna embryo development, no teratogenic effects were observed, indicating relative safety at lower concentrations.Despite observing some retardation effects and a relatively high rate of toxicity on the Daphnia eggs, these effects were only registered at the highest concentration, which showed toxicity comparable to the control on the adult crustaceans.Therefore, the extract could be considered safe as it showed no evident teratogenic effects on the crustaceans.
The toxicity assays on Daphnia species revealed that the extract exhibited low to moderate toxicity, with effects being concentration-dependent.While higher concentrations affected Daphnia pulex significantly and delayed Daphnia magna embryo development, no teratogenic effects were observed, indicating relative safety at lower concentrations.

Prediction of Molecular Targets
Potential molecular targets were predicted for several phytochemicals (apigenin, galangin, genistein, kaempferol, naringenin, caffeic acid, ferulic acid, and p-coumaric acid), which were found in higher quantities based on UHPLC-HRMS/MS analysis.Following an analysis of prediction results, we selected CDK5 and GSK-3β as potential target candidates for the investigated polyphenols, both kinases being involved in the physiopathology of neurodegenerative diseases and neuropathic pain [20,21].As shown in Table 5, apigenin was predicted as a CDK5 inhibitor with a 100% probability by SwissTargetPrediction, and with a 7.8% probability according to PASS.Moreover, apigenin had a max TC value of 61.54% when paired with known CDK5 inhibitors.Furthermore, apigenin had a 100% predicted probability of showing GSK-3β inhibitory activity according to SwissTargetPrediction algorithm and a 100% max TC value when paired with known GSK-3β inhibitors.After searching the ChEMBL database (10.1093/nar/gky1075), we found that apigenin was shown to inhibit both CDK5/p25 and GSK-3β activities, with similar potencies (IC 50 values of 1.6 µM and 1.9 µM, respectively).Another interesting compound was kaempferol, for which SwissTargetPrediction algorithm predicted a 51.79% probability of inhibiting CDK5 and a 65.80% probability of inhibiting GSK-3β activity, while PASS algorithm yielded a 6.9% probability of acting as a GSK-3β inhibitor.A previous study showed that kaempferol inhibits CDK5/p25 with an IC 50 value of 51 µM, but the activity on GSK-3β was not investigated.Furthermore, galangin was predicted with SwissTargetPrediction to have a 14.97% probability of inhibiting CDK5 and 16.61% of inhibiting GSK-3β, while PASS predicted a 9.80% probability of being active on CDK5.Interestingly, SEA calculated a 45.65% similarity between galangin and CDK5 inhibitors.Although naringenin was predicted as a potential CDK5 and GSK-3β inhibitor only by SwissTargetPrediction, previous studies showed that its glycoside, naringin, inhibits GSK-3β with an IC 50 of 100 µM.

Molecular Docking
Molecular docking studies were further performed to estimate the binding potential between the assessed polyphenols and the protein targets that were selected based on PASS, SwissTargetPrediction, and SEA predictions.Thus, the docking simulations were carried out against CDK5 and GSK-3β, both crystal structures being in active conformations (αCin architecture), since the αC-helix glutamate residues (Glu51 in CDK5 and Glu97 in GSK-3β) The predicted binding energies and ligand efficiencies of the screened phytochemicals are shown in Table 6, and the docked poses of all compounds are illustrated in Figure 8c,d.According to the predictions, apigenin showed the lowest binding energy (−9.315 kcal/mol) after docking against the ATP binding site of CDK5, followed by naringenin (−9.232 kcal/mol) and galangin (−9.039 kcal/mol).After docking with GSK-3β, the lowest binding energies were obtained for galangin (−8.553 kcal/mol), apigenin (−8.551 kcal/mol), and naringenin (−8.537 kcal/mol).We further chose to discuss the predicted molecular interactions between apigenin (as a positive control) and both protein kinases, since the flavonoid was confirmed as a CDK5 and GSK-3β inhibitor according to the data deposited in ChEMBL database.As shown in Figure 9a,b, the hydroxyphenyl sub-structure of apigenin formed a hydrogen bond with β3-strand Lys33 and Asp144 within the conserved kinase DFG motif, while the 7-hydroxy moiety engaged in hydrogen bonding with Cys83 within the hinge and Glu81.Furthermore, the catalytic Asp144 interacted with the phenyl moiety of apigenin also through pi-anion interactions.Moreover, several hydrophobic interactions (pi-sigma, pi-pi T-shaped, pi-alkyl) are formed between the benzopyran-4-one scaffold and residues Ile10, Val18, Ala31, Phe80, and β7-strand residue Leu133, which is the "floor" of the adenine binding pocket.The predicted binding energies and ligand efficiencies of the screened phytochemicals are shown in Table 6, and the docked poses of all compounds are illustrated in Figure 8c,d.According to the predictions, apigenin showed the lowest binding energy (−9.315 kcal/mol) after docking against the ATP binding site of CDK5, followed by naringenin (−9.232 kcal/mol) and galangin (−9.039 kcal/mol).After docking with GSK-3β, the lowest binding energies were obtained for galangin (−8.553 kcal/mol), apigenin (−8.551 kcal/mol), and naringenin (−8.537 kcal/mol).We further chose to discuss the predicted molecular interactions between apigenin (as a positive control) and both protein kinases, since the flavonoid was confirmed as a CDK5 and GSK-3β inhibitor according to the data deposited in ChEMBL database.As shown in Figure 9a,b, the hydroxyphenyl sub-structure of apigenin formed a hydrogen bond with β3-strand Lys33 and Asp144 within the conserved kinase DFG motif, while the 7-hydroxy moiety engaged in hydrogen bonding with Cys83 within the hinge and Glu81.Furthermore, the catalytic Asp144 interacted with the phenyl moiety of apigenin also through pi-anion interactions.Moreover, several hydrophobic interactions (pi-sigma, pi-pi T-shaped, pi-alkyl) are formed between the benzopyran-4-one scaffold and residues Ile10, Val18, Ala31, Phe80, and β7-strand residue Leu133, which is the "floor" of the adenine binding pocket.Based on the results of target prediction approaches and molecular docking simulations, we further chose to discuss the predicted interactions between the potential CDK5 inhibitor galangin and the ATP-binding site (Figure 9c,d).Similar to apigenin, galangin Based on the results of target prediction approaches and molecular docking simulations, we further chose to discuss the predicted interactions between the potential CDK5 inhibitor galangin and the ATP-binding site (Figure 9c,d).Similar to apigenin, galangin engaged in hydrophobic interactions with Ile10, Val18, Ala31, Phe80, and Leu133 and formed hydrogen bonds with Glu81 and Cys83.However, the non-hydroxylated phenyl moiety formed only pi-alkyl interactions with Lys33, while a pi-donor hydrogen bond is formed between Phe80 and the 3-hydroxy radical, which is absent in apigenin.
Next, we analyzed the predicted interactions between the positive control apigenin and the ATP binding site of GSK-3β (Figure 10a,b).The 5-hydroxy radical of apigenin formed hydrogen bonds with adenine pocket residues Val135 and Asp133, while the hydroxy-phenyl radical engaged in hydrogen bonding and pi-anion interactions with catalytic Asp200 within the DFG motif and in pi-alkyl interactions with Val70 and catalytic Lys85.Furthermore, hydrophobic interactions (pi-sigma, pi-pi stacked, pi-alkyl) were formed between the benzopyran-4-one scaffold and Ala83, Tyr134, hinge residue Cys199 and "floor" residue Leu188.Similarly, kaempferol, which was predicted as a potential GSK-3β inhibitor by PASS and SwissTargetPrediction, formed the same interactions as apigenin with Asp200 within DFG motif, catalytic Lys85, "floor" residue Leu188, hinge residue Cys199, and residues Val70m, Ala83, Asp133, and Tyr134.However, no hydrogen bonds were formed with Val135, but a pi-donor hydrogen bond was predicted between the 7-hydroxy radical and Tyr134.Nonetheless, the 3-hydroxy radical that is absent in apigenin showed no interactions with the ATP-binding site.
Plants 2024, 13, x FOR PEER REVIEW engaged in hydrophobic interactions with Ile10, Val18, Ala31, Phe80, and Leu1 formed hydrogen bonds with Glu81 and Cys83.However, the non-hydroxylated moiety formed only pi-alkyl interactions with Lys33, while a pi-donor hydrogen formed between Phe80 and the 3-hydroxy radical, which is absent in apigenin.
Next, we analyzed the predicted interactions between the positive control a and the ATP binding site of GSK-3β (Figure 10a,b).The 5-hydroxy radical of a formed hydrogen bonds with adenine pocket residues Val135 and Asp133, while droxy-phenyl radical engaged in hydrogen bonding and pi-anion interactions wi lytic Asp200 within the DFG motif and in pi-alkyl interactions with Val70 and c Lys85.Furthermore, hydrophobic interactions (pi-sigma, pi-pi stacked, pi-alky formed between the benzopyran-4-one scaffold and Ala83, Tyr134, hinge residue and "floor" residue Leu188.Similarly, kaempferol, which was predicted as a p GSK-3β inhibitor by PASS and SwissTargetPrediction, formed the same interac apigenin with Asp200 within DFG motif, catalytic Lys85, "floor" residue Leu188 residue Cys199, and residues Val70m, Ala83, Asp133, and Tyr134.However, no hy bonds were formed with Val135, but a pi-donor hydrogen bond was predicted b the 7-hydroxy radical and Tyr134.Nonetheless, the 3-hydroxy radical that is ab apigenin showed no interactions with the ATP-binding site.

Discussion
The flora in the Dobrogea Region, Romania, includes approximately 2.000 ta resenting approximately 50% of Romania's flora.The particular flora is impressive

Discussion
The flora in the Dobrogea Region, Romania, includes approximately 2.000 taxa, representing approximately 50% of Romania's flora.The particular flora is impressive from a scientific point of view considering that it includes a large number of species, some of which are rare plants (approximately 200), among which we mention the species protected at the European level: Achillea clypeolata, Centaurea jankae, Campanula romanica, Moehringia jankae, and Paeonia tenuifolia.The vegetation consists of numerous continentals, Balkan, sub-Mediterranean, and Mediterranean species [22].We chose to study the species Ajuga chamaepitys (L.) Schreb, known for its traditional uses in the Dobrogea region but not scientifically documented.
Our research aimed to evaluate the quality of a lyophilized Ajuga chamaepitys (L.) Schreb.extract for its antioxidant, cytotoxic actions and to predict its potential neuropathic pain relief properties.
The analysis of the data obtained by us compared with the existing literature is limited due to the small number of studies on this species and the different methods of extraction and analysis.In previous reports on Ajuga species from Romania, Ajuga laxmanii (L.), the following polyphenolic compounds were isolated in the methanol extract [31] using an optimized HPLC/MS method: rutin 6721.49µg/g dw (dry weight), isoquercitrin 636.1 µg/g dw, apigenin 126.53 µg/g dw, luteolin 88.24 µg/g dw, quercitrin 36.5 µg/g dw, chlorogenic acid (caffeic acid derivative) 19.33 µg/g dw.Analyzing these data, rutin showed low values in our extract, 110.6 µg/g; instead, a higher value of caffeic acid (3253.8µg/g) was quantified.Ghiţă et al. conducted a comparative phytochemical study between the species Ajuga reptans (L.) and Ajuga genevensis (L.), both collected from the spontaneous flora of Romania, where identified compounds such as chlorogenic acid, caffeic acid, apigenol, and luteolin-7-O-glucoside [32] were also present in our extract.Toiu et al. (2017) analyzed the methanolic extract from the flowers of the species Ajuga reptans (L.) collected from the spontaneous flora of Cluj, Romania.The research group identified and quantified the polyphenols by LC-MS, recording the majority component caffeic acid 34.39 µg/g dw (dry weight), followed by luteolin 26.13 µg/g dw, p-coumaric acid 26.02 µg/g dw, apigenin 23.56 µg/g dw, ferulic acid 20.15 µg/g dw, and quercitrin 4.97 µg/g dw [33].Our analysis recorded higher amounts of these compounds, except for quercitrin, which was not detected, and luteolin, whose isomer, kaempferol, was detected in our extract.
Our findings showed that ACHE has significant antioxidant activity, confirming the results of other researchers.Moreover, even though the free radical scavenging effect of ACHE is several folds lower than the reference antioxidant activity (IC 50 DPPH extract = 529.4µg/mL higher than IC 50 ascorbic acid = 16.5 µg/mL), we can observe that the Pearson coefficients (Supplementary Materials Table S3) regarding the DPPH values and the concentrations of chemical compounds that generate the effect reached the highest significant level (DPPH vs. FC: r = 0.9959, p < 0.05; DPPH vs. PAC: r = 0.9991, p < 0.001; DPPH vs. TPC: r = 0.9704, p < 0.05).These results highlighted the precision of the correlation between the variables as well as the high degree of accuracy of the DPPH method in contrast to the other antioxidant techniques.Turkoglu et al. (2010) concluded that the antioxidant effect of an aqueous extract of Ajuga chamaepitys (L.) ssp.euphratica was superior to extracts obtained with methanol and chloroform [35].The aerial parts of Ajuga chamaepitys (L.) Schreb.from Iran, extracted with methanol of various concentrations, had the best inhibition on testing the DPPH free radical scavenging activity [4].Similar results were obtained for the acetone extract of Ajuga chamaepitys ssp.chamaepitys collected from Serbia [6].In another research on the polar extracts (in ethanol and in water) and the essential oil of Ajuga chamaepitys (from herba), collected from central Italy, it was found that the antioxidant activity of the polar extracts was higher than that of the essential oil, and moderate compared with that of Trolox [7].Thus, our results are consistent with those reported in the literature, demonstrating that ACHE has comparable antioxidant properties with extracts obtained from Ajuga chamaepitys (L.), collected from different geographical areas.
Testing of Artemia species revealed no effects on cell division.Evaluation at 24-30 h after exposure showed that all larvae passed into larval stage II (microscopic observations).This passage is made by molting, a phenomenon that requires repeated divisions; therefore, the cellular processes specific to the division were not disturbed by the composition of the extract.The LC 50 value, which expresses the mortality of 50% of the subjects, was higher than 1000 µg/mL, indicating low cytotoxic action.
The toxicity evaluation on Daphnia magna and Daphnia pulex revealed different sensitivities to the Ajuga extract.D. magna demonstrated significant resistance, with negligible lethality caused by the extract.In contrast, D. pulex exhibited high sensitivity at elevated concentrations, particularly after a 48h exposure period.Although other studies showed that overall, the differences in responses between D. magna and species considered more sensitive, such as D. pulex or Ceriodaphnia dubia, were not substantial; the variability of each population could be an important factor that influences the response of the crustaceans [36,37].Our research found the use of both species advantageous because each yielded unique sensitivity responses, thus serving as a valuable tool for identifying potential toxic effects induced by plant extracts.However, toxicity was evident only at concentrations significantly higher than those typically used in phytotherapy.
The embryo test on D. magna embryos can be used to identify potential teratogenic effects [38,39].In our research, the D. magna embryo developmental assay further elucidated the impact of the extract.Thus, between 50 and 500 µg/mL, the developmental stages and viability were slightly lower or comparable with the untreated control, whereas at 1000 µg/mL, a very high concentration, the viability was reduced, the development of the larvae was delayed, and abnormal modifications were observed.
Ajuga chamaepitys (L.) Schreb from Dobrogea, Romania, is known in traditional medicine for its anti-inflammatory and analgesic effects, and therefore, we focused our research on predicting its potential use in neuropathic pain.Jaffal et al. (2019) showed that the methanolic extract of Ajuga chamaepitys from Jordan mediated the antinociceptive activity in experimental models on mice.Mediation of the analgesic effect by the involvement of an opioid receptor was evaluated by administration of an opioid receptor antagonist (naloxone), which resulted in antagonism of the analgesic effect of Ajuga chamaepitys extract.The results showed that the methanolic extract induced significant, dose-dependent analgesic effects [5].Similar results were obtained by Khanavi et al. (2014) using extracts from Ajuga chamaecistus ssp.Tomentella on mice, showing a significant analgesic effect in the chronic phase, especially the extracts from water, hexane, and ethyl ether [40].The mechanism underlying this phenomenon is thought to be the inhibition of the biosynthesis of proinflammatory prostaglandins, and the main compound responsible for this action would be 8-acetylharpagide [40], a compound also found in Ajuga chamaepitys (L.) Schreb.Although other mechanisms can be responsible for the analgesic and anti-inflammatory effects, these observations support the traditional use of Ajuga species for pain and other inflammatory diseases.Thus, we predicted potential molecular targets for several polyphenols, which we identified in the ACHE.Following the predictions, we identified serine/threonine protein kinases CDK5 and GSK-3β as potential therapeutic targets for the flavonoids apigenin, galangin, and kaempferol.Molecular docking studies further supported the interaction potential between apigenin and the ATP-binding site of both kinases, between galangin and CDK5, and between kaempferol and GSK-3β, respectively.According to bioassay data deposited in the ChEMBL database, apigenin inhibits CDK5 and GSK-3β at micromolar concentrations, while kaempferol inhibits CDK5.Since both protein kinases are dysregulated in several neurodegenerative diseases (e.g., tauopathies) and in neuropathic pain [20,21], the investigated extract could show therapeutic effects in such disorders through dual inhibition of CDK5 and GSK-3β.One limitation of the present study is the lack of experimental validation of the computational results.Further studies are warranted to support the present findings and to evaluate the therapeutic efficacy of Ajuga chamaepitys (L.) Schreb extract in animal models of neuropathic pain and other neurological conditions.

Plant Material
The aerial parts of the species Ajuga chamaepitys (L.) Schreb.were collected in May 2022 from the spontaneous flora of the Dobrogea region, Romania (44 • 29 ′ 16.5 ′′ N-28 • 27 ′ 09.5 ′′ E), in the morning, in dry and sunny weather.The species was identified by Prof. PhD Laura Bucur from the Department of Pharmacognosy, Phytotherapy and Phytochemistry, and Assoc.Prof. PhD Verginica Schröder from the Department of Celullar Biology and Genetics, Faculty of Pharmacy, "Ovidius" University Constanţa, and several specimens are in the collection of the Pharmacognosy Laboratory.The solvents and experimental conditions used in this study are presented at each step.Also, we mention that this species' spread in the Dobrogea area does not indicate problems of endangering the plant raw material; the species is adapted to the pedoclimatic conditions of this region of the country, and it spreads spontaneously and is very abundant.
Our research was carried out in several stages, as follows: (1) establishing the quality of the plant raw material; (2) obtaining and chemical characterization of the dry extract by spectrophotometric and chromatographic methods (UHPLC-HRMS/MS); (3) evaluation of the antioxidant capacity; (4) determination of cytotoxicity in vivo on the larvae of Artemia sp.(BSLA) and Daphnia; (5) determination of the teratogenic potential by the embryo test on Daphnia magna embryos; (6) molecular docking predictions; and (7) potential interactions between the quantified compounds and types of enzymes.All the results were statistically processed.

Establishing the Quality of the Plant Material
In this stage of work, we verified the quality of the aerial parts of the plant (dried in laboratory conditions) by spectrophotometric determination of three types of phytochemical compounds (flavonoids, phenolic acids, and total polyphenols).Extractions were performed separately for each type of sample in 50% ethanol, 20% ethanol, and water.Extracts were prepared by refluxing for 30 min.One g of herbal product (herba-noted with H) was used for each type of extractive solution.After cooling, the solutions were filtered in a 25 mL volumetric flask and filled to the mark with the same solvent.The extracts were stored at 4 • C until analysis.For each sample, an index corresponding to the solvent with which it was extracted was added, 50% ethanol (v/v)-indexed by a, 20% ethanol (v/v)-indexed by b, and dis+9/tilled water-indexed by c. Ethanol was purchased from Sigma-Aldrich (Darmstadt, Germany).

Quantitative Analysis of Flavonoids (FC)
In a series of 10 mL volumetric flasks, different volumes of the extractive solution were poured: 0.2, 0.4, 0.6, 0.8, and 1 mL; then, 2 mL of sodium acetate 100 g/L and 1 mL of aluminum chloride solution 25 g/L were added and filled to the mark with the same solvent used for extraction.In parallel, the control samples were prepared by measuring 0.2, 0.4, 0.6, 0.8, and 1 mL of the extractive solution and filling up to the mark with the same solvent.After 45 min, the absorbance reading was performed at λ = 427 nm, calibrating each time with the blank sample.The absorbance value must be between 0.1511 and 1.1300 (the range respects the Lambert-Beer law).FC of the extract was expressed in rutin equivalents (mg REs)/g dry herbal product.All spectrophotometric measurements were performed using a JASCO V-530 UV-VIS (Tokyo, Japan) spectrophotometer.

Quantitative Analysis of Phenolic Acids (PAC)
In a series of 10 mL volumetric flasks, different volumes of the extractive solution were poured: 0.2, 0.4, 0.6, 0.8, and 1 mL.Next, 2 mL of 0.5 M hydrochloric acid, 2 mL of Arnow reagent, and 2 mL of 85 g/L sodium hydroxide were added and filled to the mark with distilled water.In parallel, the control samples were prepared under the same conditions but without the Arnow reagent.After 5 min, the absorbance reading was performed at λ = 525 nm, calibrating each time with the blank sample.The absorbance value must be between 0.1923 and 0.9546 (the range respects the Lambert-Beer law).PAC of the extract was expressed in chlorogenic acid equivalents (mg ChAEs)/g dry herbal product.

Quantitative Analysis of Total Polyphenols (TPC)
In a series of 10 mL volumetric flasks, different volumes of the extractive solution (stock solution = 1:10) were poured: 0.4, 0.6, 0.8, and adjusted to 1 mL by adding distilled water.Then, after a few minutes, 1 mL of Folin-Ciocâlteu reagent was added and filled to the mark with 200 g/L sodium carbonate solution.In parallel, the control samples were prepared under the same conditions but without the analyzed solution.After 40 min in the absence of light, the absorbance reading was performed at λ = 763 nm, calibrating each time with the blank sample.The absorbance value must be between 0.2314 and 0.9752 (the range respects the Lambert-Beer law).TPC of the extract was expressed in tannic acid equivalents (mg TAEs/g) dry herbal product.

Obtaining and Characterizing the Dry Extract
Hundred g of dried vegetable product (after soaking the vegetable product in approximately 50-100 mL of 50% ethanol (v/v)) was refluxed with 1.5 L of 50% (v/v) ethanol for 30 min.A new amount of 50% (v/v) ethanol 1 L was added and refluxed for another 30 min.Two reflux extractions were performed.After cooling and filtering, the extractive solution obtained was homogenized and subjected to concentration up to 200 mL on a rotary evaporator (Buchi R-215 with vacuum controller V-850, Marshall Scientific, Hampton, VA, USA).The dry extract was obtained by lyophilization (lyophilizer Christ Alpha 1-2 B (Braun Biotech International GmbH, Melsungen, Germany) equipped with vacuum pump RZ 2.5).

Spectrophotometric Determination of Flavonoids, Phenolic Acids, and Total Polyphenols
For spectrophotometric analyses, 0.1 g of dry extract was solubilized in 25 mL of 50% (v/v) ethanol.The contents of flavonoids, phenolic acids, and total polyphenols were determined using spectrophotometric methods under the same conditions as previously described.
Flavonoids content of the extract was expressed as mg rutin equivalents (mg REs)/g (dw).Phenolic acid content of the extract was expressed as mg chlorogenic acid equivalents (mg ChAEs)/g (dw).
Total polyphenol content of the extract was expressed as mg tannic acid equivalents (TAEs)/g (dw).

Identification and Quantification of Polyphenol Content by UHPLC-HRMS/MS
Materials used: analytical standards (Sigma-Aldrich, Germany), methanol, ethanol, formic acid (98-100%), and ultrapure water (Merck, Bucharest, Romania).Stock solutions of 1000 mg/L, prepared by dissolving the standard substance in a suitable solvent (acethonitrile, methanol, dimethylformamide), were kept at −20 • C. Individual stock solutions were used to prepare a mixed stock solution of 4 µg/mL.By diluting the volume corresponding to the stock solution mixture in methanol, the desired concentrations were obtained in the range of 0.05-1 µg/mL for plotting the standard curve.The working solutions were stored at −20 • C for a period of one month.
Although the solvent used for the extraction was ethanol, for the analysis of these samples, methanol was used as the solvent; methanol does not significantly influence the results from a qualitative/quantitative point of view.
The parameters specific to the analysis, LC parameters, and MS parameters, are presented in the Supplementary Materials.
For compounds for which analytical standards were not available but were expected to be found in the extracts based on the literature, a list of masses was compiled, and these were identified in total ion current (TIC) by manual search.However, for unknown compounds, the most reasonable molecular formula with a lower mass error was searched in the Chemspider database (www.chemspider.comaccessed on 1 October 2023).Considering that flavones, isoflavones, and phenolic acids have many structural similarities, the ion fragments resulting from MS-MS analysis were used to confirm the chemical structure of the identified compounds using the NORMANMassBank spectral databases (https://massbank.eu/MassBankaccessed on 1 October 2023), mzCloudeTM, Advanced-Mass (https://www.mzcloud.org/accessed on 1 October 2023), and PubChem (https://pubchem.ncbi.nlm.nih.gov/accessed on 5 November 2023).ACDLabs MS Fragmenter 2019.2.1 software (https://www.acdlabs.com/products/adh/ms/index.phpaccessed on 1 October 2023) was used to generate a fragmentation model of the identified compounds for comparative analysis.
Quantification method validation, quantitative method validation, and optimization of the UHPLC and MS conditions are presented in the Supplementary Materials (Tables S1 and S2).

Evaluation of Antioxidant Activity
The antioxidant activity of the extract was determined on the basis of the scavenging capacity of the free radical ABTS •+ (2,2'-azinobis-3-ethylbenzothiazoline-6-sulfonic acid), the DPPH radical (2,2-diphenyl-2-picryl-hydrazyl), and the ferric-reducing capacity (FRAP).The antioxidant capacity was evaluated based on the value of IC 50 (µg/mL) (the concentration of the extract that inhibits free radical activity (DPPH and ABTS) by 50%) and EC 50 , the concentration at which the absorbance has a value of 0.5 for the ferric reduction capacity.The lower the IC 50 or EC 50 value, the more pronounced the antioxidant effect.

Evaluation of the Scavenger Capacity of the DPPH Radical
The method involved the reduction of the 2,2-diphenyl-2-picryl-hydrazyl radical, a stable, violet nitric radical, with the formation of the corresponding yellow hydrazine, whose concentration is proportional to the concentration of the antioxidant.The DPPH method was used according to Ohnishi M et al. with slight modifications [43].DPPH and ethanol were purchased from Sigma-Aldrich, Germany.For the experimental determinations, 0.1 mM DPPH solution was prepared daily and stored in the dark.
The absorbance value of the DPPH solution ranged between 0.9 and 1 read at λ = 515 nm.A 0.1033 g sample of extract was solubilized in 25 mL ethanol 50% (v/v) filled to the mark in the volumetric flask.Volumes of 0.1, 0.2, 0.3, 0.4, 0.7, 1, 1.5, 2, 2.5 mL of the previously obtained solution were poured into volumetric flasks of 10 mL and filled to the mark with 50% (v/v) ethanol.Each 0.5 mL of the solution corresponding to the extracts was treated with 3 mL of 0.1 mM DPPH ethanolic solution, stirred, and stored in the dark for 30 min (respectively, 60 and 90 min).The absorbances were measured at λ = 515 nm, against ethanol, which was used as a blank.In order to control the accuracy, ascorbic acid (Sigma-Aldrich, Germany) was used as a reference for the calibration curve in the range of concentration between 2 and 22 µg/mL (Supplementary Materials, Figure S4).
The formula used to calculate the inhibition (%) of DPPH radical activity is shown in the Supplementary Materials.
To evaluate the scavenger capacity of the DPPH radical, the equations of the calibration lines, inhibition (%) vs. concentration (µg/mL), for the dry extract at 30, 60, and 90 min and the IC 50 values were determined.

Evaluation of the Scavenger Capacity of the ABTS •+ Radical
The ABTS •+ radical, in the presence of an antioxidant substrate, loses its blue coloration, and the color changes are accompanied by a decrease in absorbance [44].For our determination, we chose a reaction time of 6 min, which is frequently used.The antioxidant activity was evaluated according to the method described by Re R. et al. (1999) [45].ABTS •+ ammonium salt and ethanol were purchased from Sigma-Aldrich, Germany, and potassium persulfate was purchased from Merck, Germany.The ABTS •+ radical was obtained following the reaction between the ammonium salt of 2,2-azino-bis(3-ethyl-benzothiazoline-6-sulfonic acid) (7.4 mM solution) and potassium persulfate (2.6 mM solution mM) [45], the mixture being stored in the absence of light and in contact for 16 h.A volume of 2 mL of the obtained solution was diluted in a 100 mL volumetric flask with absolute ethanol, so that the absorbance at λ = 734 nm was 0.700 ± 0.02.An amount of 0.1024 g of the extract was solubilized in 25 mL ethanol 50% (v/v) filled to the mark in the volumetric flask.Volumes of 0.05, 0.2, 0.4, 0.6, 0.8, 1 mL of the previously obtained solution were poured into a 10 mL volumetric flask and filled to the mark with 50% ethanol (v/v).
For the actual determination, 0.5 mL of the solutions were treated with 3 mL of ABTS •+ ethanolic solution, stirred, and stored in the dark for 6 min.The absorbance of the samples was measured at the maximum absorbance of the ABTS •+ solution (λ = 734 nm) compared with absolute ethanol.The reduction in the absorbance values represents the inhibition of the ABTS •+ solution and is calculated according to the equation presented in the Supplementary Materials.
To evaluate the scavenger capacity of the ABTS •+ radical, the equations of the calibration lines, inhibition (%) vs. concentration (µg/mL), for the dry extract at 6 min, with three repetitions, were used to determine the IC 50 values.

Assessment of Antioxidant Capacity Based on Ferric-Reducing Capacity (FRAP)
The ferric reduction capacity was determined using the modified FRAP [46] method.Potassium ferricyanide (Fe 3+ ) is reduced to potassium ferrocyanide (Fe 2+ ) under the action of antioxidant compounds and reacts with ferric chloride to form a complex, whose absorbance can be determined at an absorption maximum of 700 nm, with a color shift from orange to blue.A 0.1020 g extract was solubilized in 25 mL ethanol 50% (v/v) filled to the mark in the volumetric flask.Volumes of 0.1, 0.4, 0.5, 0.7, 1, 1.5, 2, 2.5 mL of the previously obtained solution were poured into 10 mL volumetric flasks and filled to the mark with 50% ethanol (v/v).From each dilution, 2.5 mL was taken and placed in test tubes, 2.5 mL phosphate buffer pH 6.6 (Sigma-Aldrich, Germany) and 2.5 mL K 3 (FeCN) 6 1% (Sigma-Aldrich, Germany) were added, and the mixture was kept in a water bath at 50 • C for 20 min.Next, 2.5 mL of 10% trichloroacetic acid solution (Sigma-Aldrich, Germany) was added to each sample.Then, 2.5 mL of solution from each test tube was placed in other test tubes, over which 2.5 mL of water and 0.5 mL of FeCl 3 0.1% (Sigma-Aldrich, Germany) were added.After 10 min., the absorbance of the samples was measured at λ = 700 nm against a blank obtained from 5 mL of water and 0.5 mL of 0.1% FeCl 3 .
The antioxidant capacity was determined using the EC 50 value (µg/mL), which represents the concentration of the solutions to be analyzed at which the absorbance has a value of 0.5.The EC 50 value (µg/mL) was determined from the equation of the regression line, concentration (µg/mL) vs. absorbance, at 30, 60, and 90 min.The BSLA test for the assessment of cytotoxicity in Artemia salina larvae is one of the fastest and most effective methods because of the sensitivity of the larvae to a variety of chemical substances [47][48][49].The test is considered a useful tool for evaluating the preliminary toxicity of plant or animal extracts [50] and is an efficient, cheap, and relatively fast way to detect toxic compounds, requiring only small amounts of sample, e.g., less than 20 mg [51,52].
The evaluation of the cytotoxic activity of the extract was assessed using in vivo assay, the BSLA test; it was quantified by the mortality rate of the Artemia salina (L.) larvae.The BSLA test was performed according to the EBPI protocol (Environmental Bio-detection Products) with modifications regarding the control sample and according to the toxicity protocol (Artoxkit M) obtaining larvae by incubating Artemia salina eggs (JBL ARTEMIO cysts) in artificial seawater.We would also like to specify that the larvae from this test offered the possibility of studying the organism with high growth with the undifferentiated cell characteristics.Also, in the first hours of life, the tested substance effects were correlated with cytotoxicity, without interferences with digestive enzymes, because the larvae did not require feeding [53].
The conditions of the experiment: the testing was carried out in a static system in microplates, adapted according to the ARTOXKIT protocol, with a volume of 1 mL, at t = 21 • C, pH = 7-7.2,saline water 5‰.
Obtaining Artemia larvae: Artemia salina larvae were obtained by introducing the cysts into 35‰ saline solution for 24 h under conditions of continuous lighting and aeration.After hatching, larvae in larval stage I (instar I) were separated and placed in experimental vessels, Plexiglas plates with wells (total volume of 1 mL), in saline solutions of 5‰ (Figure S10 in the Supplementary Materials).
The evaluation was carried out with a stereomicroscope to record the effects, and for the analysis of details such as checking the stages or highlighting the development of the larvae, microscopic preparations were made.Organisms exposed to different concentrations of dry extract were evaluated; the movements and passage of larvae in stages II and III were followed, and mortality was recorded after 24 h of exposure, which was the response quantification parameter.

Daphnia Species Toxicity Assay
The bioassays were performed according to the described protocol, with some modifications.Briefly, daphnids from the species Daphnia magna and Daphnia pulex were selected from parthenogenetic cultures prior to the assay.The bioassay was conducted using 12-well tissue culture plates (Greiner Bio-One), with each well housing 10 organisms.The extract was tested at six different concentrations, ranging from 50 to 1000 µg/mL, with each test being duplicated, and an untreated control was used as a negative control.The samples Plants 2024, 13, 1192 20 of 24 were maintained at constant temperature and humidity (25 • C, 75% RH) in a climatic chamber (MLR-351H; Sanyo).Lethality was observed at the 24 and 48 h marks, and the LC 50 and 95% confidence interval (CI95%) of the LC 50 values for each sample were calculated using the least square fit method, facilitated by GraphPad Prism v 5.1 software.

Daphnia magna Embryo Developmental Assay
The embryo test was performed according to the described protocol in the literature [54] in 12-well tissue culture plates (Greiner Bio-One, Kremsmünster, Austria), with each well containing five embryos.The extract was tested in three concentrations: 50, 500, and 1000 µg/mL.The concentrations were selected following the results obtained in the toxicity assay, and the test was performed against the untreated control.For each sample, two replicates were performed, and the results were compared with an untreated control in the conditions described in the Toxicity Assay section.Every 24 h, the embryos were subjected to microscopic examination (Euromex bScope, Arnhem, The Netherlands) in order to identify the developmental stages and abnormalities.

Prediction of Biological Activity for Selected Phytochemicals
The biological activities of the phytoconstituents that were identified in higher quantities were predicted using three web-servers: PASS (Prediction of Activity Spectra for Substances), SwissTargetPrediction, and SEA (Similarity Ensemble Approach).PASS is a predictive algorithm based on Level 2 Multilevel Neighborhoods of Atom descriptors that uses a Bayesian approach to estimate for input ligands the probability of being active (Pa) or inactive (Pi) on a specific target or the probability of showing a certain therapeutic effect [55].SwissTargetPrediction is a tool that predicts the probability of interacting with target proteins for small-molecule ligands based on a combination of 2D and 3D similarity measures between query molecules and known active ligands [56].SEA is an approach that relates proteins based on the chemical similarity of known ligands and calculates Tanimoto similarity scores (Tc) for each query ligand based on 2D topological Daylight fingerprints [57].SMILES codes of the selected polyphenolic compounds were used as input variables for each approach.The results were manually filtered in order to select key therapeutic targets relevant to neurological disorders.

Molecular Docking Studies
Molecular docking studies were performed to further investigate the potential interactions between the selected polyphenols and previously predicted protein targets.Based on PASS, SwissTargetPrediction and SEA results, we selected CDK5 (cyclin-dependent kinase 5) and GSK-3β (glycogen synthase kinase-3β) for molecular docking experiments.Crystal structures of human CDK5 in complex with activator p25 and a naphthyridine inhibitor (PDB ID: 7VDR, 2.55 Å resolution [58]) and human GSK-3β in complex with a pyrazine inhibitor (PDB ID: 6HK3, 2.35 Å resolution [59]) were retrieved from the RCSB PDB database.Prior to docking, macromolecule structures were prepared using the YASARA structure [60] as follows: water, solvent molecules, and ligands were removed, structural errors were corrected, side-chains were protonated according to the physiological pH (7.4), the hydrogen-bonding network was optimized, and the conformation was optimized by energy minimization with a YASARA2 forcefield.
Ligands were processed by generating 3D coordinates with DataWarrior 5.2.1 [61], followed by energy minimization using the MMFF94s+ forcefield, and protonated according to the physiological pH.Docking simulations were performed using AutoDock Vina v1.1.2algorithm [62].The docking search space was set to include the ATP-binding site of both kinases.The docking protocol was validated by docking the co-crystallized inhibitors into the binding site, followed by the determination of RMSD (root mean square deviation) values after superposition of the predicted poses on the experimental coordinates.The docking results were retrieved as predicted binding poses, binding energies (∆G, kcal/mol), and ligand efficiencies (LE, ∆G/no. of heavy atoms).The predicted interactions between the ligands and target proteins were assessed using UCSF ChimeraX v1.6.1 [63] and BIOVIA Discovery Studio Visualizer (BIOVIA, Discovery Studio Visualizer, Version 17.2.0,Dassault Systèmes, 2016, San Diego, CA, USA).

Statistical Analysis
First, descriptive statistics were performed, and all data are expressed as mean ± standard deviation.Before applying any statistical tests, key conditions were evaluated for each set of experimental data, such as data normality and variance homogeneity.Normality of the data was assessed using the Shapiro-Wilk test.For non-normally distributed data, we transformed them (using base the 10 logarithm) to achieve normality before conducting statistical tests.
To identify significant differences between data groups, a one-way ANOVA was employed with a Tukey post hoc test for multiple comparisons with statistical significance at a 95% confidence level (p < 0.05).Principal component analysis (PCA) with Pearson's correlation coefficient was performed to identify the relationships between phytochemical compounds and antioxidant capacity.
Clustering analysis was conducted to create groups with similar behavior within and differences between groups.Hierarchical cluster analysis (HCA) was used to identify homogeneous compound groups within each extract based on the exact mass, adduct ion (m/z)/monitored negative ion, and retention time.The "silhouette" method was used to determine the optimal number of groups.
All statistical analyses of the data obtained were performed using GraphPad Software 9, Inc., San Diego, CA, USA, and XLSTAT for Excel 2021 version (Addinsoft, New York, NY, USA).

Conclusions
In conclusion, according to the performed analysis, it has been shown that Ajuga chamaepitys (L.) Schreb., which is widespread in the spontaneous flora of Dobrogea, Romania, can be an important source of polyphenols with potential inhibitory activity on specific targets in neurodegenerative diseases and neuropathic pain.Further studies are warranted to confirm the predicted molecular mechanisms of action and to further investigate the therapeutic potential in animal models of neurological disorders.

Figure 2 .
Figure 2. Results of the spectrophotometric determinations of the dry extract of ACHE.The data are represented as mean ± SD and analyzed using ANOVA with a Tukey post hoc test.* p < 0.05, ns: p > 0.05 no statistical differences.

Figure 5 . 1 .
Figure 5. Antioxidant activity of ACHE.The data are represented as mean ± SD and analyzed using ANOVA with a Tukey post hoc test.*** p < 0.001, ns: p > 0.05 no statistical differences.

Figure 6 .
Figure 6.The lethality curves obtained after 48 h of exposure of Daphnia pulex to ACHE; error bars represent the SE of two replicates.2.4.3.Daphnia Magna Embryo Developmental Assay

Figure 6 . 3 .
Figure 6.The lethality curves obtained after 48 h of exposure of Daphnia pulex to ACHE; error bars represent the SE of two replicates.
Figure S9.Probit values correlated with the tested concentrations.
Figure S10.The test system in microplates with wells (1mL test volume).
Figure S11.Regression calculated for the evaluated effects (mortality).Table

Table 1 .
Results of the spectrophotometric determinations for extractive solutions.

Table 1 .
Results of the spectrophotometric determinations for extractive solutions.

Table 4 .
Concentration values at which lethal effects were recorded (at the 95% confidence level).

Table 5 .
Predicted activity on CDK5 and GSK-3β using three different algorithms.

Table 6 .
Predicted binding energies and ligand efficiencies following docking in the ATP-binding site of CDK5 and GSK-3β.

Table 6 .
Predicted binding energies and ligand efficiencies following docking in the ATP-binding site of CDK5 and GSK-3β.