Identification of isolated compound
Compound 1 is a yellow crystal, showing fluorescence at UV λ254 and 365 nm, and appears as a yellow stain with 10% sulfuric acid spray and heating, shown in S-1 (Supplementary Data) UV (MeOH) λmax: 358 and 208.9 nm. IR (KBr pellet) νmax: 3788, 3249, 2946, 1664, and 1593 cm-1. TQD-MS with negative ionization mode (ESI-) showed m/z 343.93. The 1H-NMR (CDCl₃, 700 MHz) showed signals at δH 12.64 ppm (1H, s, 5-OH), 7.71 ppm (1H, dd, J = 1.98 Hz, H-2′), 7.68 ppm (1H, dd, J = 2.02, 8.4 Hz, H-6′), 6.46 ppm (1H, d, J = 2.2 Hz, H-8), 6.37 ppm (1H, d, J = 2.5 Hz, H-6), 3.98 ppm (3H, s, 3′-OMe), 3.88 ppm (3H, s, 7-OMe), and 3.86 ppm (3H, s, 3-OMe). The spectrum of 13C-NMR (CDCl₃, 175 MHz) displayed 18 chemical shifts at δC 178.9 (C-4), 165.6 (C-7), 162.2 (C-5), 156.9 (C-2), 156.1 (C-9), 148.5 (C-4′), 146.5 (C-3′), 139.0 (C-3), 122.8 (C-6′), 122.6 (C-1′), 114.7 (C-5′), 111.0 (C-2′), 106.2 (C-10), 98.0 (C-6), 92.3 (C-8), 60.3 (3-OMe), 56.3 (3′-OMe), and 56.0 (7-OMe).
The 13C-NMR and DEPT 135° spectra of compound 1 displayed 18 carbon signals, including one carbonyl (δC 178.9 ppm), seven oxygenated aromatic quaternary carbons (δC 162.2-178.9 ppm), one non-oxygenated sp2 quaternary carbon (δC 122.6 ppm), five sp2 methine carbons (δC 92.3-122.8 ppm), and three oxygenated methyl carbons (δC 55.8-60.2 ppm). The 1H-NMR spectrum showed a singlet at δH 12.64 ppm, typical of a chelated hydroxyl group, along with aromatic proton signals consistent with ortho- and meta-substituted protons, and three methoxy singlets at δH 3.98, 3.88, and 3.86 ppm. These features are characteristic of a flavonoid with a flavone backbone, further supported by 2D NMR data. The 1D and 2D NMR spectra were consistent with the data reported in the literature [35, 36], thereby confirming the identity of the compound as pachypodol. The chemical structures of pachypodol (1) are shown in Fig. 1.
Identification of Compound 2
Compound 2 is a colorless oil, showing fluorescence at UV light λ254 nm, and appears as a yellow stain with 10% sulfuric acid spray and heating (See Supplementary Data S-2) UV (MeOH) λmax: 267.9 and 208.9 nm. IR (KBr pellet) νmax: 3436, 3076, 2938, 1686, 1589, 1509, 1288, and 1125 cm-1. TQD-MS with negative ionization mode (ESI-) showed m/z 417.28. The NMR spectrum of compound 2 (CDCl₃; 1H-NMR 700 MHz; 13C-NMR 175 MHz) showed 1H signals at δH 6.28 (1H, s, H-2), 6.26 (1H, s, H-6), 3.33 (1H, d, J = 6.5 Hz, H-7), 1.30 (3H, d, J = 7.1 Hz, H-9), 4.01 (1H, s, H-4′), 6.14 (1H, s, H-6′), 2.77 (2H, d, J = 7.0 Hz, H-7′), 5.90 (1H, ddt, J = 10.2, 17.0 Hz, H-8′), and 5.21–5.30 (2H, m, H-9′). The methoxy protons resonated at δ 3.68 (3H, s, 3-OMe), 3.78 (3H, s, 5-OMe), 3.77 (3H, s, 4-OMe), 3.78 (3H, s, 3′-OMe), and 3.43 (3H, s, 5′-OMe). The 13C-NMR spectrum displayed signals at δC 192.4 (C-2′), 153.1 (C-3, C-5), 137.2 (C-1′, C-4), 134.6 (C-8′), 133.8 (C-1), 125.2 (C-6′), 119.0 (C-9′), 106.2 (C-2, C-6), 95.9 (C-3′, C-5′), 79.1 (C-4′), 59.6 (C-7), 49.3 (C-8), 35.0 (C-7′), and 17.7 (C-9), together with methoxy carbons at δC 61.0 (4-OMe), 56.1 (5-OMe, 3′-OMe), 55.6 (3-OMe), and 54.8 (5′-OMe).
The 13C-NMR and DEPT 135° spectra of compound 2 revealed 23 carbons, including eight quaternary carbons (two sp3 at δC 49.2 and 95.9 ppm, and six sp2 at δC 133.8, 137.2, 152.7, 153.1 (2C), and 192.4 ppm), two methylenes at δC 35.1 and 119.2 ppm, seven methines (δC 49.3, 59.6, 79.1, 106.2 (2C), 125.5, and 134.6 ppm), and six methyl groups comprising one sp3 (δC 17.7 ppm) and five methoxy carbons (δC 54.8, 55.6, 56.1 (2C), and 61.0 ppm). The 1H-NMR spectrum showed resonances corresponding to five methoxy groups (δH 3.43, 3.68, 3.77, 3.78 (2×) ppm), aromatic/olefinic protons (δH 6.14, 6.26, 6.28 ppm), and a methyl doublet at δH 1.30 ppm (J = 7.1 Hz), along with methine and allylic signals at δH 3.33, 4.01, 5.21, 5.30, and 5.90 ppm. Further NOESY analysis revealed stereochemical correlations among protons at C-4, C-16, C-11, and C-9. The 1D and 2D NMR spectra were consistent with the data reported in the literature [37], thereby confirming the identity of the compound as nectamazin B. Its chemical structure is shown in Fig. 1.
Antimicrobial activity assay of compounds
Investigation of the antimicrobial effects of pachypodol (1) and nectamazin B (2) against S. mutans and C. albicans was carried out through Kirby-Bauer agar diffusion, MIC and MBC assay. The MIC value is referred to as strong category if MIC 101-500, moderate if 501-1000, and weak if MIC > 1000 μg/mL [38, 39].
The findings, as displayed in Table 1, indicates that compounds 1 and 2 have strong inhibitory activity against the two microbes. The MBC values exceed the MIC, which clearly indicates that this compound possesses bacteriostatic properties rather than being bactericidal. This suggests that both pachypodol and nectamazin B is capable of inhibiting S. mutans and C. albicans growth without eliminating the bacteria, thereby exhibiting a growth-inhibitory effect.
Table 1 In vitro antibacterial and antifungal activity of pachypodol and nectamazin B from P. crocatum against S. mutans and C. albicans
|
Compound
|
Microorganism
|
Inhibition zone (mm)
|
MIC (µg/mL)
|
MBC/MFC (µg/mL)
|
|
1
|
S.mutans ATCC 25175
|
10.2
|
156.25
|
1250
|
|
C. albicans ATCC 10231
|
9.0
|
156.25
|
625
|
|
2
|
S.mutans ATCC 25175
|
10.67
|
156.25
|
625
|
|
C. albicans ATCC 10231
|
9.61
|
156.25
|
625
|
Molecular docking study
The docking results of pachypodol (1) and nectamazin B (2) against C. albicans virulence proteins (Gtf B and Sap 5) are presented in Table 2. Pachypodol exhibited the strongest binding affinity toward Gtf B with a ΔG of −8.05 kcal/mol and an inhibition constant (Ki) of 1.26 μM, which is superior compared to the native ligand (ΔG −6.85 kcal/mol, Ki 9.51 μM). In contrast, its interaction with Sap 5 was relatively weaker (ΔG −6.67 kcal/mol, Ki 12.89 μM).
Nectamazin B demonstrated moderate affinity toward Gtf B (ΔG −6.70 kcal/mol, Ki 12.2 μM), which is comparable to the native ligand. Interestingly, nectamazin B showed a stronger affinity toward Sap 5 (ΔG −7.56 kcal/mol, Ki 2.88 μM), nearly equivalent to that of the native ligand (ΔG −7.61 kcal/mol, Ki 2.63 μM). These findings suggest that pachypodol may serve as a potent inhibitor of Gtf B, while nectamazin B displays preferential activity toward Sap 5, indicating their potential as selective antifungal agents targeting distinct virulence factors of C. albicans.
Table 2 Binding affinity data and inhibition constants of compounds 1 and 2 against Gtf B and Sap 5 enzymes.
|
Compound
|
Protein Target
|
ΔG (kcal/mol)
|
Ki (µM)
|
|
Pachypodol (1)
|
Gtf B
|
-8.05
|
1.26
|
|
Sap 5
|
-6.67
|
12.89
|
|
Nectamazin B (2)
|
Gtf B
|
-6.70
|
12.2
|
|
Sap 5
|
-7.56
|
2.88
|
|
Ligan Native
|
Gtf B
|
-6.85
|
9.51
|
|
Sap 5
|
-7.61
|
2.63
|
The molecular docking interaction profiles of pachypodol (1) and nectamazin B (2) with enzymes Gtf B and Sap 5 are summarized in Fig. 2 and Fig. 3. For the Gtf B enzyme (Fig. 2), pachypodol formed multiple stabilizing interactions, including conventional hydrogen bonds with ASP A:454, TYR A:404, and HIS A:561, as well as a covalent interaction with ASP A:883, TYR A:890, and GLU A:489. Additionally, several hydrophobic contacts were observed (Pi–alkyl with LEU A:407 and ALA A:452; Pi–Pi T-shaped with TRP A:491; and alkyl with PHE A:881, VAL A:453, ALA A:490). These diverse interactions may explain the strong binding affinity of pachypodol to Gtf B (ΔG −8.05 kcal/mol, Ki 1.26 μM). In comparison, nectamazin B interacted with fewer residues. It mainly formed hydrogen bonds with ASP A:552 and TRP A:491, along with a covalent bond involving ASP A:451 and GLU A:489. The reduced number of hydrophobic and electrostatic contacts relative to pachypodol may underlie its weaker affinity toward Gtf B.
In Sap 5 enzyme (shown in Fig. 3), both compounds exhibited multiple binding interactions. Pachypodol formed conventional hydrogen bonds with TYR B:225 and ARG B:120, as well as covalent interactions with residues such as TYR B:221, ASP B:86, and SER B:118. Furthermore, hydrophobic stabilization was mediated by Pi–alkyl interactions with ILE B:123, ILE B:30, and ALA B:119. Interestingly, nectamazin B displayed a richer binding profile with Sap 5. It established three hydrogen bonds (TYR B:225, ASP B:86, and TYR B:225), as well as multiple covalent and hydrophobic interactions (e.g., LYS B:50, TRP B:51, ARG B:86). Importantly, nectamazin B also engaged in Pi–anion and Pi–hydrogen donor interactions, which may contribute to its strong binding affinity (ΔG −7.56 kcal/mol, Ki 2.88 μM), comparable to the native ligand. Overall, these findings suggest that pachypodol exhibits a higher binding strength toward Gtf B, supported by multiple hydrogen bonds and hydrophobic interactions, whereas nectamazin B preferentially binds Sap 5 through diverse hydrogen bonding, hydrophobic, and Pi interactions. This differential binding behavior highlights their potential as selective inhibitors against distinct virulence proteins of C. albicans.
Molecular dynamics
The stability and binding affinity of the protein–ligand complexes were assessed using molecular dynamics simulations combined with MM/GBSA calculations. The RMSD trajectories (Fig. 4a) indicated that the GtfB–native ligand and GtfB–compound 1 (Pachypodol) complexes remained relatively stable throughout the 100 ns simulation, with mean RMSD values of 2.096 ± 0.325 Å and 1.823 ± 0.312 Å, respectively. In contrast, the GtfB–compound 2 (Nectamazin B) complex exhibited higher structural deviations (2.983 ± 0.410 Å), reflecting reduced stability.
Residue-level flexibility, as shown in the RMSF profiles (Fig. 4b), revealed that both the native ligand and compound 2 complexes followed similar fluctuation patterns across most residues, while compound 1 induced localized flexibility in loop regions. These shifts may reflect conformational adjustments that facilitate stronger binding interactions.
The energetic stability of the complexes was further evaluated using MM/GBSA analysis (Table 5). Compound 1 displayed the most favorable binding free energy (−28.223 ± 5.394 kcal/mol), significantly stronger than both the native ligand (−7.092 ± 4.877 kcal/mol) and compound 2 (−15.042 ± 3.020 kcal/mol). To explore binding stability over time, MM/GBSA binding free energy values were computed using a 10-ns sliding window approach (Fig. 5). Compound 1 complex consistently exhibited the lowest ΔG values, remaining between −30 and −20 kcal/mol for most of the trajectory, highlighting its robust binding affinity. In contrast, compound 2 showed moderate stability with binding energies fluctuating around −15 to −20 kcal/mol, while the native ligand displayed weaker and more variable interactions, with values often closer to −10 kcal/mol.
Taken together, the RMSD, RMSF, and MM/GBSA results consistently demonstrate that compound 1 (Pachypodol) forms the most stable and energetically favorable complex with GtfB, outperforming both the native ligand and compound 2 (Nectamazin B).
Table 5 RMSD and MMGBSA value of complex protein and ligand after 100 ns molecular dynamics simulation
|
Complex
(Protein-Ligand)
|
RMSD (Å)
|
MMGBSA (kcal/mol)
|
|
Mean SD
|
Median MAD
|
Mean SD
|
Median MAD
|
|
Gtf B - Native
|
2.096 0.325
|
2.078 0.314
|
-7.092 4.877
|
-5.832 5.313
|
|
Gtf B - Compound 1
|
1.823 0.312
|
1.841 0.305
|
-28.223 5.394
|
-29.516 3.372
|
|
Gtf B - Compound 2
|
2.983 0.410
|
2.994 0.361
|
-15.042 3.020
|
-13.983 2.373
|
The hydrogen bonding profiles of the complexes were also analyzed to evaluate the stability of ligand–protein interactions over the simulation period (Fig. 6). Compound 1 maintained between one and three hydrogen bonds consistently throughout the trajectory, with intermittent formation of up to four bonds, indicating stable and persistent interactions within the binding pocket. In contrast, compound 2 formed a higher number of hydrogen bonds overall, but these interactions were more fluctuating and less stable, suggesting weaker binding persistence despite frequent contacts. The native ligand showed fewer hydrogen bonds, typically fluctuating between zero and two, with intermittent peaks up to three, reflecting less stable hydrogen-bonding interactions compared to compound 1 and compound 2.
These findings highlight that although compound 2 exhibited a greater number of hydrogen bonds, their dynamic instability reduces their overall contribution to binding affinity. Conversely, compound 1 formed fewer but more stable hydrogen bonds, complementing its favorable MM/GBSA binding energy and lower RMSD, thereby reinforcing its potential as the most stable binder among the tested ligands.
ADMET and drug-likeness
The ADMET predictions of the lead compounds provided further insights into their potential suitability as oral care agents (Table 6). Both compound 1 (Pachypodol) and compound 2 (Nectamazin B) showed good intestinal absorption as indicated by favorable Caco-2 permeability values, suggesting that they could efficiently permeate epithelial membranes if administered systemically. However, for a mouthwash formulation, local retention and safety are more critical than systemic absorption.
In terms of distribution, neither compound was predicted to cross the blood–brain barrier (BBB−), which minimizes potential central nervous system side effects. Both exhibited high plasma protein binding (>85%), although this parameter is less relevant for topical use in oral formulations. With respect to metabolism, the two compounds showed distinct profiles. Pachypodol was predicted to be a substrate and potential inhibitor of CYP1A2, while Nectamazin B was associated with CYP3A4 interaction. Although such interactions are important in systemic drug design, they are of lower concern for a mouthwash, as the formulation is not intended for extensive systemic circulation.
Excretion analysis indicated low clearance values and moderate half-lives for both compounds, suggesting reasonable metabolic stability. Toxicity assessment revealed low LD50 values for both ligands, classifying them as compounds with relatively low acute toxicity risk. Neither compound was predicted to be hepatotoxic or mutagenic in the Ames test. Importantly, Pachypodol was classified as non-sensitizing and non-immunotoxic, making it safer for repeated oral exposure. In contrast, Nectamazin B showed a potential for skin sensitization and immunotoxicity, raising safety concerns for topical use.
Taken together, the ADMET evaluation suggests that Pachypodol exhibits a more favorable pharmacological and safety profile for development as an oral rinse agent compared to Nectamazin B, aligning well with its stronger binding stability and favorable interaction patterns observed in molecular dynamics simulations.
Table 6 Predicted ADMET properties of Pachypodol (1) and Nectamazin B (2) based on in silico analysis, including absorption, distribution, metabolism, excretion, and toxicity parameters
|
Properties
|
Pachypodol (1)
|
Nectamazin B (2)
|
|
Value
|
Interpretation
|
Value
|
Interpretation
|
|
Absorption
|
Caco-2
|
-5,118
|
Good permeability
|
-4,993
|
Good permeability
|
|
Distribution
|
BBB
|
0,001
|
BBB-
|
0,685
|
BBB-
|
|
PPB
|
96,172
|
High plasma protein binding
|
88,605
|
High plasma protein binding
|
|
Metabolism
|
CYP1A2
|
0,999
|
Substrate and potential inhibitor
|
0,003
|
Not a substrate or inhibitor
|
|
CYP2D6
|
0,054
|
Not a substrate or inhibitor
|
0,011
|
Not a substrate or inhibitor
|
|
CYP3A4
|
0,008
|
Not a substrate or inhibitor
|
0,901
|
Substrate and potential inhibitor
|
|
CYP2B6
|
0,026
|
Not a substrate or inhibitor
|
0,0
|
Not a substrate or inhibitor
|
|
Excretion
|
Clearance
|
4,708
|
Low clearance
|
7,197
|
Low clearance
|
|
T1/2
|
1,308
|
Moderate half-life
|
1,428
|
Moderate half-life
|
|
Toxicity
|
LD50
|
0,427
|
Low toxicity
|
0,633
|
Low toxicity
|
|
Hepatotoxicity
|
0,404
|
Non-hepatotoxic
|
0,423
|
Non-hepatotoxic
|
|
Mutagenicity (Ames)
|
0,588
|
Non-mutagenic (Ames test)
|
0,314
|
Non-mutagenic (Ames test)
|
|
Skin sensitization
|
0,54
|
Non-sensitizer
|
0,93
|
Potential sensitizer
|
|
Immunotoxicity
|
0,088
|
Non-immunotoxic
|
0,206
|
Potential immunotoxicity
|