G-10615-TiS Template (Full Article) - 36658-1-2-20250509
G-10615-TiS Template (Full Article) - 36658-1-2-20250509
https://doi.org/
Highlights
1. Anhuienoside E was successfully isolated from the methanol extract of N. sativa seeds.
2. Anhuienoside E has antibacterial activity against S. mutans, S. sanguinis, E. faecalis based on MIC and MBC
tests.
3. Anhuienoside E has the potential to inhibit Gtf and MurA enzymes based on In silico molecular docking
study.
Graphical abstract
Trends Sci. 202x; xx(xx): xxxxx 2 of 22
Abstract
Dental caries is a chronic disease suffered by almost the entire population in the world. The main bacterium
that causes dental caries is Streptococcus mutans, which has the enzyme glucosyltransferase as a virulence factor.
Other fatogenic bacteria that play a role in exacerbating the biofilm form of dental caries such as Streptococccus
sanguinis and Enterococcus faecalis have cell wall defenses catalyzed by the MurA enzyme. Chlorhexidine has
been reported as a treatment for dental caries but has developed resistance over time. Nigella sativa L. seeds have
been widely recognized to have many health benefits such as antibacterial, antioxidant and antifungal. The caries-
causing antibacterial activity of N. sativa seeds has not been widely reported. The aim of this study was to isolate
antibacterial compounds from N. sativa against S. mutans, S. sanguinis, E. faecalis, predict the mechanism of
compounds and their derivatives by in silico molecular docking and pharmacokinetic analysis by ADMET and
drug-likeness methods. Isolation of compounds was carried out by column chromatography with bioassay
guidance, inhibitory mechanism and sugar substituent effects were predicted by in silico molecular docking,
pharmacokinetics and drug-likeness analysis were predicted by ADMET and Lipinski rules. Anhuienoside E was
successfully isolated from N. sativa extract and demonstrated antibacterial activity, which had moderate MIC and
MBC values of 625 μg/mL and 1000 μg/mL against S. mutans, S. sanguinis, E. faecalis. Anhuienoside E had
moderate binding affinity value compared to its derivatives with ∆G values of −6.84 and −8.67 Kcal/mol against
MurA and Gtf. ADMET pharmacokinetic analysis and drug-likeness evaluation suggest that Anhuienoside E and
its derivatives may serve as non-toxic, non-oral drug candidates. In conclusion, Anhuienoside E has potential
antibacterial activity to be developed as an alternative to chlorhexidine. The amount and type of sugar substituents
greatly affect the antibacterial activity of Anhuienoside E.
Introduction
Dental caries is a chronic disease that affects nearly the entire global population. According to Indonesia's
2018 Basic Health Research (Riskesdas) data, 95% of children aged 5 to 6 years suffer from dental caries [1]. The
main bacteria that causes dental caries is Streptococcus mutans. S. mutans secretes extracellular polysaccharides
by glucosyltransferase enzymes and gives the signals to other bacteria to form colonies on the biofilm layer [2–
3]. Streptococcus sanguinis and Enterococcus faecalis receiving signals from S. mutans work together to make a
stronger biofilm layer into caries [4–5]. Both of these bacteria have cell walls composed of peptidoglycans, which
are synthesized by the MurA enzyme. Various treatments of dental caries have been reported in previous studies.
Chlorhexidine is one of antiseptics that has been used commercially. Chlorhexidine targets the destruction of
bacterial cell membranes. The ability of bacteria to regulate the DltA operon causes the increasing of
hydrophobicity on the surface of cell membranes so chlorhexidine is becoming resistance [6]. Treatment of dental
caries with more specific targets needs to be explored, one of which is by targeting the glucosyltransferase (Gtf)
and Muramidase A (MurA) enzymes.
Gtf is an enzyme that catalyzes the synthesis of extracellular polysaccharides on the tooth surface. It is
secreted by Streptococcus mutans as a virulence factor and breaks the glycosidic bond in sucrose, producing
glucose and fructose. Glucose monomers are polymerized and act as a signal for other bacteria to form colonies
on the teeth surface (biofilm) which then form caries [7].
MurA is the first enzyme that catalysed the formation of peptidoglycans to build bacterial cell wall. Uridine
diphosphate N-acetylglucosamine is converted into UDP-MurNac by MurA enzyme. In the next stage, it will be
attached to sugar molecules to form peptidoglycans. Therefore, inhibition of this enzyme becomes very potent so
bacteria will grow abnormalities or die [8–9]. Exploration of natural compounds that have potential as caries
antibacterials by in silico molecular docking of inhibitory mechanisms against MurA and Gtf enzymes needs to
be done. One of the natural sources which has antibacterial agent is the seeds of Black Cumin (Nigella sativa
Linn.).
N. sativa Linn. is an herbal plant and has been known for a long time among people in all parts of the world.
The high benefits of this plant, especially in the seeds, make it often used as an alternative to conventional
medicine. N. sativa seeds are widely used as an herbal medicine in Indonesia called “Jamu”. The method of use
is pounded then boiled and filtered water. The efficacy of this “Jamu” is believed to treat various diseases such as
inflammation, diabetes and cold medicine. Several researchers revealed the benefits of extracts from N. sativa
seeds as antibacterial, antioxidant, immunomodulatory and antiviral [10–11].
The content of secondary metabolite compounds found in N. sativa seeds such as flavonoids, triterpenoids,
saponins, alkaloids and essential oils provide diverse bioactivities [12]. Natural product compounds have diverse
functional groups such as hydroxy, carbonyl, amine and ester. The role of these functional groups involves redox
reactions in prokaryotic cells. So that the presence of functional groups from compounds in N. sativa will increase
antibacterial activity [13–16]. The main compound isolated from the essential oil of N.sativa seeds is
thymoquinone, which has potential therapeutic effects against atherosclerosis, cancer, and diabetes [17–19].
Trends Sci. 202x; xx(xx): xxxxx 3 of 22
The treatment of dental caries with N. sativa seed compounds has not been widely reported. The study of
the mechanism of inhibition of Gtf and MurA enzymes by compounds from N. sativa through molecular docking
approach is very useful to know the treatment action of compounds. Therefore, compounds from N. sativa have
the potential to be explored as natural caries drug candidates with a review of mechanism of action studies. This
study aimed to identify antibacterial compounds from N. sativa seed extract effective against S. mutans, S.
sanguinis, and E. faecalis. to determine their mechanisms of action targeting MurA and Gtf enzymes, and to
predict their pharmacokinetic properties through ADMET and drug-likeness analyse.
N. sativa seeds were obtained from Indonesian local market at Jl. Pasar Barat No.44, Kebun Jeruk, Andir,
Bandung in July 2023. Dried N. sativa seeds were ground then macerated, extracted and purified using n-hexane,
ethyl acetate, methanol and water solvents from Sigma Aldrich Co. Ltd. (St. Louis, MO, USA) and Merck Co.
Ltd. A total of 2 kg of dried N. sativa seeds were ground and macerated using 4 liters of methanol for 3 days, then
filtered and concentrated using a rotatory evaporator. The methanol extract was tested for antibacterial activity
against S. mutans ATCC 25175, S. sanguinis ATCC 10556 and E. faecalis ATCC 29212 at concentrations of 1,
2, 3, 4 and 5%. The extract that was known to have oral antibacterial activity was partitioned as much as 30 g
using a normal phase chromatography column Silica G 60 F254 as much as 300 g. The extract was eluted by n-
hexane, ethyl acetate, methanol and water of 300 mL for 3 times, respectively until get eleven fractions. The
fraction that had good antibacterial activity with the large amount of mass was fraction 7 (F-7) (methanol fraction,
14.52 g) which is continued to the purification stage. Analysis of spot patterns in F-7 used TLC (Thin Layer
Chromatography) ODS (Octadecylsilane) RP-18 F254S and Silica G 60 F254 (Merck, Darmstadt, Germany). It was
observed under UV light at wavelengths of 254 and 365 nm and sprayed with 10% H2SO4 solution in ethanol
(v/v). 1 g of F-7 was performed column chromatography using ODS RP-18 F254S and elution was performed using
30 mL of a methanol–water mixture (9:1, v/v) with a 2% v/v gradient of increasing polarity, resulting in 32
subfractions.[20] Column results were analyzed using TLC and observed the stain pattern under UV light 254 and
365 nm, then sprayed with H2SO4 10% stain reagent in ethanol (v/v) [21].
3. Antibacterial Testing
The methanol extract of N. sativa seeds was tested for antibacterial activity against S. mutans ATCC 25175,
S. sanguinis ATCC 10556 and E. faecalis ATCC 29212 at concentrations of 1, 2, 3, 4 and 5% by disc diffusion
method. Positive control (chlorhexidine 2%), negative control methanol and extract as much as 20 μL were placed
on paper disks and then placed on Muller Hinton agar media (purchased from Merck Co. Ltd. and Sigma Aldrich)
containing 100 μL bacteria ½ McFarland. Agar media was incubated for 24 hours at 37 ℃. Determination of the
active fraction used the same procedure as the N. sativa seed methanol extract assay.
MIC testing was carried out on compound 1 by microdilution method at a concentration of 2% using
microplate-96 wells. All wells of microplate-96 wells were filled with Muller Hinton liquid media, then 100 μL
of methanol solvent was added as control. Compound 1 was added to the wells without methanol as much as 100
μL. Dilutions were performed by the microdilution method on Microplate-96 wells. Bacteria with a turbidity level
of 0.5 McFarland as much as 5 μL were added to the wells. The microplate-96 wells were incubated for 24 hours
at 37°C. After incubation, the absorbance was observed on the Biochrom EZ Read 400 ELISA tool with a
wavelength of 620 nm and the MIC value was obtained. MBC testing is done by spreading clear media microplate-
96 wells containing compound 1 and bacteria on Muller Hinton agar media then incubated at 37°C for 24 hours
so that the MBC value can be determined [24]. Antibacterial testing of zone of inhibition was performed in
triplicate, MIC and MBC were performed in duplicate.
Trends Sci. 202x; xx(xx): xxxxx 4 of 22
The enzymes used in molecular docking were MurA and Gtfs which have PDB ID: 8FKL[25] and 1UAE
[26] whose 3D structures were downloaded on the website www.rcsb.org. Conformation validation of MurA and
Gtf enzymes using Ramachandran Plot [27]. The two Ramachandran-conforming enzyme PDB codes 8FKL and
1UAE were separated with bile and ligand natives using the Biovia discovery studio. Coordinates of the enzyme
active site were analyzed on the web server http://sts.bioe.uic.edu/castp/index [28]. Enzyme preparation using
Biovia Discovery Studio 2021 software. The structures of compound 1 and its derivatives were drawn in
ChemDraw 15.0 and energy minimization was performed in Chem3D Ultra 8.0. Chlorohexidine and apigenin as
positive controls were downloaded on the page https://pubchem.ncbi.nlm.nih.gov with CID: 9552079 and
5280443 codes.
Molecular docking was run using AutoDockTools-1.5.7 software.[29] Receptors (Gtf and MurA) were
pretreated by adding Kollman charge while ligands compound 1, derivates and positive control were added to
Compute Gasteiger. Docking was performed on receptor active site regions obtained from calculations on the
CASTp webserver [28]. Docking parameters using Genetic algorithms with 100 runs. Visualization of docking
results (type of interaction and number of interactions) was observed in Biovia discovery studio. Adding grid box
area is better
ADMET analysis of compound 1 and derivatives using web server https://biosig.lab.uq.edu.au/pkcsm. The
chemical structure of compound 1 and its derivatives is made in smile format which is drawn on the web
http://www.swissadme.ch. Determination of drug-likeness using web server https://tox.charite.de/protox3 [30–
31][32].
Maceration of 2 kg of N. sativa seeds using 4L methanol was concentrated to obtain 46.5 g of extract. The
methanol extract of N. sativa seeds was tested for antibacterial activity against S. mutans, S. sanguinis and E.
faecalis through inhibition zone test. The test results of the extract showed active as antibacterial against the three
bacteria (test results are in Table 1) which indicated the compound components potentially as antibacterial agent
[33]. The methanol extract was further fractionated into eleven fractions and tested for antibacterial activity. Based
on the results of antibacterial testing of the eleven fractions, it showed that fraction seven had good oral
antibacterial activity, so it was continued to the purification stage (test results are in Table 2). The isolation of
compounds from the active fraction aims to obtain the most active isolate. Some isolates showed more active
antibacterial activity than the fraction and others were less active depending on their nature in the antagonistic or
synergistic fraction [34].
Fraction F-7 from the initial partition was purified using ODS RP-18 F254S chromatography column with
2% gradient H2O-methanol eluent starting from 9:1 solvent ratio until 32 fractions were obtained. The purification
results were analyzed using TLC. Based on the KLT analysis on the ODS RP-18 plate, F-7-30 (100 mg) has a
single stain that does not fluoresce under UV lights 254 and 365 nm and is brown after being sprayed with H2SO4.
This fraction is thought to be a pure isolate which is called compound 1 [35–36].
The results of measurements using IR spectrophotometer indicated that compound 1 has an OH group at
maks 3367 cm-1, C-H sp3 stretch at maks 2940 cm-1, C=O stretch vibrations at maks 1634 cm-1, aliphatic C=C
stretch vibrations at maks 1455 cm-1, gem dimethyl C(CH3)2 stretch vibrations at maks 1385 cm-1, C-O stretch
vibrations at maks 1058 cm-1, and alkene C-H bending at maks 635 cm-1 [37].
Trends Sci. 202x; xx(xx): xxxxx 5 of 22
Chemical shift of carbon compound (1) (brown solids, 100 mg) showed 13C-NMR (CD3OD, 700 MHz) δ
13.8 ppm (C-1), 16.6 ppm (C-2), 17.8 ppm (C-3), 17,9 ppm (C-4), 18.0 ppm (C-5), 18.8 ppm (C-6), 23.9 ppm (C-
7), 24.2 ppm (C-8), 24.6 ppm (C-9), 26,4 ppm (C-10), 26.6 ppm (C-11), 28.9 ppm (C-12), 31.6 ppm (C-13), 33.3
ppm (C-14), 33.3 ppm (C-15), 33.5 ppm (C-16), 34.8 ppm (C-17), 37.6 ppm (C-18), 39.7 ppm (C-19), 40.6 ppm
(C-20), 42.5 ppm (C-21), 43.9 ppm (C-22), 44.0 ppm (C-23), 47.1 ppm (C-24), 48.0 ppm (C-25), 48.1 ppm (C-
26), 49.9 ppm (C-27 Rha), 61.8-82.4 ppm (C-(28-52 sugar), 95.7, 101.3, 104.2, 104.6 and 106.5 ppm (C-(53-
57anomeric), 123.7 ppm (C-58), 144.0 ppm (C-59), 178.0 ppm (C-60); 1H-NMR (CD3OD, 700 MHz) δ 0.72 ppm
(3H, s; H-1), 0.82 ppm (3H, s; H-4), 0.94 ppm (3H, s; H-16), 0.97 ppm (3H, s; H-8), 1.00 ppm (3H, s; H-2), 1.19
ppm (3H, s; H-10), 1.25 ppm (2H, t J = 6, H-14), 1.28 ppm (3H, s; H-3), 1.29 ppm (3H, s; H-5), 1.39 ppm (2H, t;
J = 12, H-17), 1.51 ppm (2H, t; J = 12, H-6), 1.63 ppm (2H, t; J = 6, H-7), 1.64 ppm (3H, s; H-28), 1.67 ppm (4H,
t; J = 6, H-9&15), 1.71 ppm (4H, t; J = 6, H-11&12), 1.78 ppm (2H, d; J = 6, H-24), 1.79 ppm (2H, q; J = 6, H-
24), 1.92 ppm (1H, t; J = 6, H-27), 1.93 ppm (1H, t; J = 6, H-25), 2.89 ppm (1H, q; J = 6, H-21), 3.87 ppm (1H,
t; J= 6; H-52), 3.89 ppm (1H, t; J= 6; H-58), 3.90-4.1 ppm (21H; H-sugar), 4.42 ppm (1H, d; J= 7.8 H-55
anomeric), 4.51 ppm (1H, d; J = 7,5 H-57 anomeric), 4.55 ppm (1H, d; J = 6.2 H-56 anomeric), 5.26 ppm (1H, d;
J = 16.3 H-54 anomeric), 5.36 ppm (1H, d; J = 8.1 H-53 anomeric). 1H, 13C- and 2D-NMR chemical shift data of
compound 1 are shown in Table 3-4.
Based on the 13C-NMR data, compound 1 contained sixty types of carbon. DEPT, 1H-NMR spectrum and
HSQC correlation showed compound 1 contained eight quaternary carbons, five methine carbons, twenty-four
oxidized methine carbons, one sp2 methine, nine methyl, and thirteen methylene. HMBC spectrum shows proton
to-carbon correlation as far as two to five bonds and the COSY spectrum shows proton to proton correlation as
far as three bonds are suspected of compound 1 being a burnt triterpenoid group [38]. Carbon sixty and carbon
fifty-two experienced the addition of sugar substituents as many as five pieces shown by the presence of carbon
and anomeric protons as many as five. Three sugars namely glucose’- glucose”-rhamnose’ were bound to carbon
sixty and two glc-rha sugars at carbon fifty-two. The placement of the number and type of sugar was based on the
HMBC correlation of the anomeric proton to carbon sixty and fifty-two proton to anomeric carbon.[39]
Configuration of sugar group in compound 1 was β which was indicated by the shear of the anomeric carbon ≥99
and α at ≤95 [40–41]. NMR spectrum showed that compound 1 had the molecular formula C60H98O6 with a degree
of unsaturation of twelve. Measurement of the molecular mass of compound 1 using a Water Acquit UPLC type
triquadrupole [M+Na]+ at m/z 1257 confirmed the alleged structure of compound 1 (Shown in Figure 1) [42].
Based on the literature, compound 1 has proton and carbon chemical shifts similar to Anhuienoside E isolated
from N. sativa seeds. In addition, based on MS measurements, it has a similar molecular weight, thus confirming
that compound 1 is Anhuienoside E (all spectra in Figure S1-S7) [12–43].
Antibacterial testing of Anhuienoside E (1) against S. mutans, S. sanguinis and E. faecalis was carried out
by microdilution method. The test results are in Table 5 which showed moderate activity with MIC values of 625
μg⁄mL and MBC of 1000 μg⁄mL [44]. Anhuienoside E (1) showed that did not give activity for killing E. faecalis
but only gave inhibition activity. Anhuienoside E had five sugar groups and other functional groups such as
alkenes and carbonyls that could inhibit bacteria through the destruction of cell membranes and inhibit enzymes
that have important catalytic activities in bacteria [45]. The moderate MIC values were very relevant to MBC
values that were in the inactive range. The study of the antibacterial mechanism of Anhuienoside E (1) and the
effect of sugar groups on antibacterial activity was predicted by molecular docking against key pathogenic
enzymes of oral bacteria.
Trends Sci. 202x; xx(xx): xxxxx 6 of 22
Figure 1: Structure of Anhuienoside E (1) with thick line COSY spectrum correlation and red line HMBC correlation.
Figure 3: Chemical structure of Anhuienoside D (derivate 2), derivate 8, Anhuienoside F (derivate 4), and Oleanolic acid (derivate 9)
Trends Sci. 202x; xx(xx): xxxxx 8 of 22
Figure 4: Chemical structure of Anhuienoside E (1), Flaccisdoside II (derivate 6), Flaccisdoside III
(derivate 7), chlorhexidine and Apigenin (positive control)
Trends Sci. 202x; xx(xx): xxxxx 9 of 22
4. Molecular docking
Molecular docking was performed on Anhuienoside E (1) and derivates on two key enzymes from caries-
causing bacteria (the structures are shown in Figure 2-4). MurA enzyme was an enzyme that catalysed
peptidoglycan at the initiation stage representing E. faecalis and S. sanguinis bacteria. Gtf enzyme with PDBID
8FKL isolated from S. mutans catalysed the formation of extracellular polysaccharides starting from the cleavage
of sucrose glycosides [46]. Both enzymes were analysed for conformational similarity of amino acid bile before
characterization and after visualization in the protein data bank using the Ramachandran plot. Ramachandaran
plot describes the conformation of α-helical and β-sheet bonds in protein peptides with two coordinate torsions ϕ
and ψ (ψ 180, ϕ>0). Based on the Karplus equation, the Ramachandran plot was divided into four regions to
determine whether the conformation of amino acids was still stable during X-ray and does not damage the α-
helical and β-sheet bonds. One [A, B, L] region was highly favoured, additional allowed regions [a, b, l, p],
generously allowed regions [~a, ~b, ~l, ~p] and disallowed regions as shown in Figure 5. Amino acid residues
from proteins that will be used for molecular docking should not occupy disallowed regions of more than 2% to
ensure conformational conformity with the actual protein structure [47–49]. Both MurA and Gtf have a disallowed
region value of 0%, indicating that their conformations are consistent with the original structures and can be used
for in silico molecular docking.
(1) (2)
Figure 5: Ramachandran plot of MurA enzyme (1) and Gtf enzyme (2)
Anhuienoside E (1), Anhuienoside C (3), Anhuienoside D (2), Anhuienoside F (4), Cussonoside B (5),
Flaccidoside II (6), and Flaccidoside III (7), Derivate (8), Oleanolic acid (9), chlorhexidine, apigenin were docked
with MurA and Gtf enzymes. Compounds 1–9 were subjected to molecular docking to compare their inhibitory
activity against MurA and Gtf enzymes, with a focus on variations in the number and type of sugar substituents.
Sugar groups had many hydroxyl functional groups that can interact by hydrogen bond with enzymes. Molecular
docking results demonstrate that sugar groups significantly influence the strength of interactions between ligands
(compounds) and receptors (enzymes), as characterized by varying delta G (∆G) values, as shown in Figures 5
and 6. ∆G is a parameter representing the strength of ligand-receptor binding, measured in kcal/mol, and is derived
from the following equation:
∆G= -RT ln K
The ∆G value decreases as the K value (the equilibrium constant between protein and ligand) increases [50]. The
greater the value of K the reaction goes towards the product or the stronger the protein-ligand association. Docking
is performed at the receptor (enzyme) active site, meaning the ligand (compound) acts as a competitive inhibitor
against the native ligand or cofactor. Ligand binding in the active site region can alter enzyme stability or inhibit
its activity [51].
Anhuienoside E (1) had a relatively moderate ∆G value compared to other compounds and positive controls
against the MurA enzyme. The number of glucose units significantly affects ligand–protein binding affinity, as
evidenced by Anhuienoside F (4), which exhibited the lowest binding affinity (the highest ∆G value) toward
Trends Sci. 202x; xx(xx): xxxxx 10 of 22
MurA among the tested compounds, due to the presence of six sugar groups (Figure 6). The results of docking
visualization in Figure 8 showed that almost all compounds interact with ligands through sugar groups excepted
compounds four and seven with carbonyl esters. Sugar groups covered other functional groups in terpenoid
structure because the interactions formed (hydrogen bonds) were stronger. MurA enzyme structure was also
relatively small, so bulk compounds such as Anhuienoside E were difficult to interact with the active site area.
MurA facilitates glycosidic bond formation in peptidoglycan biosynthesis, and the presence of sugar groups may
disrupt this activity, preventing bacterial cell wall formation [12]. The type of sugar affected the value of ∆G.
Almost all compounds (1-8) interacted by hydrogen bond with the enzyme through the sugar rhamnose. The
hydroxy position at number four on rhamnose tended to be the most electronegative because there was an EDG
(electron donating group) group, namely either of C-anomeric and hydroxy group at carbon five which pushed
electrons through the induction effect [52]. Compound 9 did not show the lowest ∆G value, likely due to the
formation of a single hydrogen bond, highlighting the influence of the type and number of sugar groups on MurA
enzyme inhibition. In addition, the positive control, chlorhexidine, did not show a higher ∆G value than some of
the tested compounds (1-9), suggesting that Anhuienoside E, with proper derivatization, has potential as an
alternative to chlorhexidine.
The docking results of compounds 1-9 with Gtf enzyme showed various ∆G values (shown in Figure 7).
Flaccisdoside III (7) exhibited the lowest ∆G value through its interaction with the Gtf enzyme, specifically at the
Glc'-Glc'' glycosidic bond (Figure 9). The Gtf enzyme catalyzes the synthesis of extracellular polysaccharides,
initiating the process by cleaving the glycosidic bond of sucrose (Glc–fructose) [7]. The presence of Flaccisdoside
III (7) in the active site region decreased the catalytic activity of the enzyme through hydrogen bonds formed from
aspartic amino acids with Glc'-acid protons and Glc'-oxygen. In addition, in other compounds 1-8 there were
hydrogen bonds (Shown in Figure 9) with Gtf which were very dominant from Glc to Gtf amino acid residues in
the active site region. Anhuienoside E had inhibitory activity similar to the positive control Apigenin as shown
by the ∆G values in Figure 7. This finding is supported by the number of hydrogen bond interactions formed by
the glucose moieties of Anhuienoside E. All sugar-containing derivatives demonstrated stronger inhibitory
activity against the Gtf enzyme compared to the positive control.
ADMET analysis results of compounds 1-9 are show in Table 6. Analysis of Absorption (A), Distribution
(D), Metabolism (M), Excretion (E) and Toxicity (T) was conducted to investigate the feasibility of these as
antibacterial drugs. Absorption parameters (Intestinal Absorption and Water Solubility) indicated the drug's
ability to be absorbed by the body based on its water solubility. The range of water solubility was log 95% of
drugs: -6.0 to 0.5 [53]. Only compounds five, eight and nine were absorbed by the body and the rest were carried
by the blood to all body tissues. This data showed that Anhuienoside E, can reach the target of pathogenic bacteria.
Distribution parameters (VDss, CNS and BBB): all compounds are poor distributed in the blood or in the
bloodstream. A good VDss range was 0.5 to 3 log L/Kg. The drug candidate should not affect the central nervous
system (CNS) and BBB (Blood-Brain Barrier). The drug concentration should be focused on the disease target
[54–55]. BBB and CNS >2.0 indicated strong absorption, whereas 0.1-2.0 indicated moderate absorption and <0.1
indicated low absorption. Analysis of the effect of compounds 1-9 on central nervous system metabolism of
CYP1A2, CYP2D6, CYP2C9, and CYP2C19 showed that all compounds did not inhibit the four enzymes. These
enzymes were CYP isoenzymes that catalysed the metabolism of drugs so they could be digested or absorbed by
the body [56]. Good drugs did not inhibit the activity of these enzymes. Compounds 1-9 exhibited substantial total
clearance values, indicating how easily these compounds, particularly thymoquinone, can be secreted and their
overall benefit to the body. Acute oral toxicity analysis with the LD50 parameter in compounds 1-9 showed the
possibility of being harmful if ingested. This parameter was used to measure how toxic or non-toxic the oral drug
is when tested on rats. LD50 values in the range of 50 mg/kg ≤ LD50 >5 mg/kg indicated the drug fatal if swallowed,
toxic in the range of 300 mg/kg ≤ LD50 >50 mg/kg, oral drugs could be harmful if swallowed in the range of 2000
mg/Kg-5000 mg/Kg, and non-toxic in the range of LD50 >5000 mg/kg [56].
The potential of a compound as an oral drug candidate can also be assessed through its physicochemical
properties, particularly its compliance with Lipinski's Rule of Five (RO5) [57]. RO5 rules must have a specific
gravity of less than 500 Daltons, the number of H bond acceptors was less than 10 bonds, H bond donors were
less than 5, lipophilicity was less than 5 bonds and molar refractivity was 40-130. A compound was allowed to
be oral drug-like if there had only a maximum one deviation from the RO5 rule. RO5 analysis of compounds 1-9
was contained in Table 7. The compound that complied with RO5 parameters was only compound 9. Compounds
1-8 do not qualify as oral drugs, as they have a Violation score of 4. However, they have the potential to be used
as mouthwash as an alternative to chlorhexidine [58].
Trends Sci. 202x; xx(xx): xxxxx 11 of 22
-5
-5.38 -4.75
-10 -6.86 -6.45 -7.31 -8.13
-9.18 -9.39 -8.94
-11.66 -10.96
-15
1 2 3 4 5 6 7 8 9 A B
Figure 6: Graph of binding affinity (∆G) values between compounds (1-9), apigenin (A), chlorhexidine (B) with
MurA enzyme
-5
Figure 7: Graph of binding affinity (∆G) values between compounds (1-9), apigenin (A), chlorhexidine (B) with
Gtf enzyme
Trends Sci. 202x; xx(xx): xxxxx 12 of 22
1. 2. 5. VAL 4. ARG
A:163 A:120 PRO
A:121
VAL
THR GLU
A:122
A:326 A:329
THR
A:326 ARG
A:120 PRO
A:298
VAL
VAL LYS
A:327
VAL VAL A:161 A:160
LYS A:161
A:327
A:160 THR
ALA VAL
PRO A119 A:362
A:327
A:122
3.
A. PHE
ASN
B. ARG SER
PRO A:332 THR
A:303 A:67 A:120 A:162
VAL A:326 LEU
PRO A:161 A:124
VAL A:298 GLU
VAK A:326 LYS
A:163
A:327 A:160
ARG
A:227 TRP ARG
GLU HIS
A:298 A:91
A:329 ASN A:125
A:253
6. VAL
7. 9. 8.
A:161
THR
A:326
LYS
A:160 ARG
ALA A:91
VAL A:119 VAL
VAL PRO
A:122 PRO A:327
A:161 A:121
A:298
ARG ALA LYS
ASN ASP GLU
A:120 A:119 A:160
A:184 A:159 A:329
C.
Interactions:
Hydrogen bond
π -Alkyl
π -sigma
π -donor Hydrogen bond
Figure 8: 2D visualization of docking results of Anhuienoside E (1), Anhuienoside D (2), Anhuienoside C (3), Anhuienoside F (4), Cussonoside B (5), Flaccidoside II (6),
Flaccidoside III (7), their derivatives (9), Oleanolic acid (9), Chlorhexidine (A), Apigenin (B) docked with MurA and 3D visualization (C)
Trends Sci. 202x; xx(xx): xxxxx 13 of 22
3.
A. PRO
PHE
ASP
B. ARG SER
A:332
A:303 THR A:688 A:120 A:162
VAL A:326 LEU
A:162 GLU
VAL PRO A:124
A:325 SER LYS
A:163 A:298 VAL A:160
GLY A:690
A:327
A:689 PHE
GLU A:396 TRP HIS ARG
A:329 A:95 A:125 A:91
LYS PHE TYR
A:712 A:619 A:394
6. 7. 9. 8.
GLY PHE
ASP A:689 A:650
LYS A:688
A:712
LYS
ASN A:664 ASP
A:655 ALA
A:688
A:649
LYS GLN
A:457 SER A:648
LYS
A:395 PRO GLY
A:712
A:713 A:451
C.
Interactions:
Hydrogen bond
π -Alkyl
π -sigma
π -donor Hydrogen bond
Figure 9: 2D visualization of docking results of Anhuienoside E (1), Anhuienoside D (2), Anhuienoside C (3), Anhuienoside F (4), Cussonoside B (5), Flaccidoside II (6),
Flaccidoside III (7), their derivatives (8), Oleanolic acid (9), Chlorhexidine (A), Apigenin (B) docked with Gtf and 3D visualization (C).
Trends Sci. 202x; xx(xx): xxxxx 14 of 22
Conclusions
The pure compound Anhuienoside E (1) was successfully isolated from the methanol extract of N. sativa seeds and
demonstrated moderate antibacterial activity against S. mutans, S. sanguinis, and E. faecalis. In silico molecular docking
analysis indicated that derivatization of Anhuienoside E (1) is necessary to improve its antibacterial potency by targeting
the inhibition of MurA and Gtf enzymes. ADMET and drug-likeness analyses revealed that Anhuienoside E (1) is not
suitable for oral administration. However, it shows potential as a mouthwash agent and could serve as a promising
alternative to chlorhexidine, which has shown increasing resistance in the treatment of dental caries. In the present study,
derivatization of Anhuienoside E (1) was not performed. Further research is required to enhance its antibacterial efficacy,
particularly through sugar moiety modifications aimed at MurA and Gtf enzyme inhibition.
Acknowledgements
The authors are grateful to Academic Leadership Grant (ALG) Prof. Dikdik Kurnia, M.Sc., Ph.D. Universitas
Padjadjaran Indonesia for all research facilities
References
[1] R Amalia, F Chairunisa, MF Alfian, and A Supartinah. Indonesia: Epidemiological profiles of early childhood
caries. Frontiers in Public Health 2019; 7,. https://doi.org/10.3389/fpubh.2019.00210.
[2] JA Lemos, SR Palmer, L Zeng, ZT Wen, JK Kajfasz, IA Freires, J Abranches, and LJ Brady. The Biology of
Streptococcus mutans . Microbiology Spectrum 2019; 7,. https://doi.org/10.1128/microbiolspec.gpp3-0051-2018.
[3] S Javed, M Zakirulla, RU Baig, SM Asif, and AB Meer. Development of artificial neural network model for
prediction of post-streptococcus mutans in dental caries. Computer Methods and Programs in Biomedicine 2020;
186,. https://doi.org/10.1016/j.cmpb.2019.105198.
[4] B Zhu, LC Macleod, T Kitten, and P Xu. Streptococcus sanguinis biofilm formation & interaction with oral
pathogens. Future Microbiology 2018; 13, 915–932. https://doi.org/10.2217/fmb-2018-0043.
[5] L Xin, Z Hu, R Han, X Xu, C Wang, D Li, Y Guo, and F Hu. Asp50Glu mutation in MurA results in fosfomycin
resistance in Enterococcus faecium. Journal of Global Antimicrobial Resistance 2022; 30, 50–55.
https://doi.org/10.1016/j.jgar.2022.05.026.
[6] S Huang, M Wu, Y Li, J Du, S Chen, S Jiang, X Huang, and L Zhan. The dlt operon contributes to the resistance
to chlorhexidine in Streptococcus mutans. International Journal of Antimicrobial Agents 2022; 59, 106540.
https://doi.org/https://doi.org/10.1016/j.ijantimicag.2022.106540.
[7] Y Lin, J Chen, X Zhou, and Y Li. Inhibition of Streptococcus mutans biofilm formation by strategies targeting
the metabolism of exopolysaccharides. Critical Reviews in Microbiology 2021; 47, 667–677.
https://doi.org/10.1080/1040841X.2021.1915959.
[8] MM Miyachiro, D Granato, DM Trindade, C Ebel, AF Paes Leme, and A Dessen. Complex Formation between
Mur Enzymes from Streptococcus pneumoniae. Biochemistry 2019; 58, 3314–3324.
https://doi.org/10.1021/acs.biochem.9b00277.
[9] LD Garbinski, BP Rosen, and M Yoshinaga. Organoarsenicals inhibit bacterial peptidoglycan biosynthesis by
targeting the essential enzyme MurA. Chemosphere 2020; 254,.
https://doi.org/10.1016/j.chemosphere.2020.126911.
[10] E Sutrisna, S Wahyuni, and A Fitriani. Antibacterial Effect of Nigella sativa L. Seed from Indonesia.
Pharmacognosy Journal 2022; 14, 1029–1032. https://doi.org/10.5530/pj.2022.14.206.
[11] M Ikhsan, N Hiedayati, K Maeyama, and F Nurwidya. Nigella sativa as an anti-inflammatory agent in asthma.
BMC Research Notes 2018; 11,. https://doi.org/10.1186/s13104-018-3858-8.
[12] D Kurnia, R Padilah, E Apriyanti, and HDA Dharsono. Phytochemical Analysis and Anti-Biofilm Potential That
Cause Dental Caries from Black Cumin Seeds (Nigella sativa Linn.). Drug design, development and therapy
2024; 18, 1917–1932. https://doi.org/10.2147/DDDT.S454217.
[13] A Shah, S Akhtar, F Mahmood, S Urooj, AB Siddique, MI Irfan, M Naeem-ul-Hassan, M Sher, A Alhoshani, A
Rauf, HMA Amin, and A Abbas. Fagonia arabica extract-stabilized gold nanoparticles as a highly selective
colorimetric nanoprobe for Cd2+ detection and as a potential photocatalytic and antibacterial agent. Surfaces and
Interfaces 2024; 51, 104556. https://doi.org/10.1016/j.surfin.2024.104556.
[14] S Ullah, R Khalid, MF Rehman, MI Irfan, A Abbas, A Alhoshani, F Anwar, and HMA Amin. Biosynthesis of
phyto-functionalized silver nanoparticles using olive fruit extract and evaluation of their antibacterial and
antioxidant properties. Frontiers in Chemistry 2023; 11,. https://doi.org/10.3389/fchem.2023.1202252.
Trends Sci. 202x; xx(xx): xxxxx 15 of 22
[15] A Jabbar, A Abbas, N Assad, M Naeem-ul-Hassan, HA Alhazmi, A Najmi, K Zoghebi, M Al Bratty, A Hanbashi,
and HMA Amin. A highly selective Hg 2+ colorimetric sensor and antimicrobial agent based on green synthesized
silver nanoparticles using Equisetum diffusum extract. RSC Advances 2023; 13, 28666–28675.
https://doi.org/10.1039/D3RA05070J.
[16] AB Siddique, MA Shaheen, A Abbas, Y Zaman, MA Bratty, A Najmi, A Hanbashi, M Mustaqeem, HA Alhazmi,
Z ur Rehman, K Zoghebi, and HMA Amin. Thermodynamic and kinetic insights into azo dyes photocatalytic
degradation on biogenically synthesized ZnO nanoparticles and their antibacterial potential. Heliyon 2024; 10,
e40679. https://doi.org/10.1016/j.heliyon.2024.e40679.
[17] AF Majdalawieh, SM Yousef, and IA Abu-Yousef. Thymoquinone, a major constituent in Nigella sativa seeds,
is a potential preventative and treatment option for atherosclerosis. European Journal of Pharmacology 2021;
909,. https://doi.org/10.1016/j.ejphar.2021.174420.
[18] MF Mahomoodally, MZ Aumeeruddy, LJ Legoabe, D Montesano, and G Zengin. Nigella sativa L. and Its Active
Compound Thymoquinone in the Clinical Management of Diabetes: A Systematic Review. International Journal
of Molecular Sciences 2022; 23,. https://doi.org/10.3390/ijms232012111.
[19] P Gnanasekaran, A Roy, N Sirpu Natesh, V Raman, P Ganapathy, and MK Arumugam. Removal of microbial
pathogens and anticancer activity of synthesized nano-thymoquinone from Nigella sativa seeds. Environmental
Technology & Innovation 2021; 24, 102068. https://doi.org/https://doi.org/10.1016/j.eti.2021.102068.
[20] W Huang. Open tubular ion chromatography: A state-of-the-Art review. Analytica Chimica Acta 2021; 1143,
210–224. https://doi.org/10.1016/j.aca.2020.08.012.
[21] S Šegan, D Opsenica, and D Milojković-Opsenica. Thin-layer chromatography in medicinal chemistry. Journal
of Liquid Chromatography and Related Technologies 2019; 42, 238–248.
https://doi.org/10.1080/10826076.2019.1585615.
[22] IA Evangelina, Y Herdiyati, A Laviana, R Rikmasari, C Zubaedah, . Anisah, and D Kurnia. Bio-Mechanism
Inhibitory Prediction of β-Sitosterol from Kemangi (Ocimum basilicum L.) as an Inhibitor of MurA Enzyme of
Oral Bacteria: In vitro and in silico Study. Advances and Applications in Bioinformatics and Chemistry 2021;
Volume 14, 103–115. https://doi.org/10.2147/AABC.S301488.
[23] JR Hanson. A hundred years in the elucidation of the structures of natural products. Science Progress 2017; 100,
63–79. https://doi.org/10.3184/003685017X14876775256129.
[24] P Parvekar, J Palaskar, S Metgud, R Maria, and S Dutta. The minimum inhibitory concentration (MIC) and
minimum bactericidal concentration (MBC) of silver nanoparticles against Staphylococcus aureus . Biomaterial
Investigations in Dentistry 2020; 7, 105–109. https://doi.org/10.1080/26415275.2020.1796674.
[25] N Schormann, M Patel, L Thannickal, S Purushotham, R Wu, JL Mieher, H Wu, and C Deivanayagam. The
catalytic domains of Streptococcus mutans glucosyltransferases: a structural analysis. Acta Crystallographica
Section F: Structural Biology Communications 2023; 79, 119–127.
https://doi.org/10.1107/S2053230X23003199.
[26] D Kurnia, GS Hutabarat, D Windaryanti, T Herlina, Y Herdiyati, and MH Satari. Potential allylpyrocatechol
derivatives as antibacterial agent against oral pathogen of S. Sanguinis ATCC 10,556 and as inhibitor of MurA
Enzymes: In vitro and in silico study. Drug Design, Development and Therapy 2020; 14, 2977–2985.
https://doi.org/10.2147/DDDT.S255269.
[27] O V. Sobolev, P V. Afonine, NW Moriarty, ML Hekkelman, RP Joosten, A Perrakis, and PD Adams. A Global
Ramachandran Score Identifies Protein Structures with Unlikely Stereochemistry. Structure 2020; 28, 1249-
1258.e2. https://doi.org/10.1016/j.str.2020.08.005.
[28] W Tian, C Chen, X Lei, J Zhao, and J Liang. CASTp 3.0: Computed atlas of surface topography of proteins.
Nucleic Acids Research 2018; 46, W363–W367. https://doi.org/10.1093/nar/gky473.
[29] MS Valdés-Tresanco, ME Valdés-Tresanco, PA Valiente, and E Moreno. AMDock: a versatile graphical tool for
assisting molecular docking with Autodock Vina and Autodock4. Biology Direct 2020; 15,.
https://doi.org/10.1186/s13062-020-00267-2.
[30] LLG Ferreira and AD Andricopulo. ADMET modeling approaches in drug discovery. Drug Discovery Today
2019; 24, 1157–1165. https://doi.org/10.1016/j.drudis.2019.03.015.
[31] A Daina, O Michielin, and V Zoete. SwissADME: A free web tool to evaluate pharmacokinetics, drug-likeness
and medicinal chemistry friendliness of small molecules. Scientific Reports 2017; 7,.
https://doi.org/10.1038/srep42717.
[32] CY Jia, JY Li, GF Hao, and GF Yang. A drug-likeness toolbox facilitates ADMET study in drug discovery. Drug
Discovery Today 2020; 25, 248–258. https://doi.org/10.1016/j.drudis.2019.10.014.
[33] Y Bibi, S Nisa, FM Chaudhary, and M Zia. Antibacterial activity of some selected medicinal plants of Pakistan.
BMC Complementary and Alternative Medicine 2011; 11,. https://doi.org/10.1186/1472-6882-11-52.
[34] MD Awouafack, LJ McGaw, S Gottfried, R Mbouangouere, P Tane, M Spiteller, and JN Eloff. Antimicrobial
activity and cytotoxicity of the ethanol extract, fractions and eight compounds isolated from Eriosema robustum
Trends Sci. 202x; xx(xx): xxxxx 16 of 22
[55] L Zhu, J Zhao, Y Zhang, W Zhou, L Yin, Y Wang, Y Fan, Y Chen, and H Liu. ADME properties evaluation in
drug discovery: in silico prediction of blood–brain partitioning. Molecular Diversity 2018; 22, 979–990.
https://doi.org/10.1007/s11030-018-9866-8.
[56] T Siswina, MM Rustama, D Sumiarsa, E Apriyanti, H Dohi, and D Kurnia. Antifungal Constituents of Piper
crocatum and Their Activities as Ergosterol Biosynthesis Inhibitors Discovered via In Silico Study Using
ADMET and Drug-Likeness Analysis. Molecules 2023; 28,. https://doi.org/10.3390/molecules28237705.
[57] CA Lipinski. Lead- and drug-like compounds: The rule-of-five revolution. Drug Discovery Today: Technologies
2004; 1, 337–341. https://doi.org/10.1016/j.ddtec.2004.11.007.
[58] MP Pollastri. Overview on the rule of five. Current Protocols in Pharmacology 2010;
https://doi.org/10.1002/0471141755.ph0912s49.
Table 2. Inhibition zone test results of fractions 1-11 of N. sativa seed methanol extract against three oral pathogenic bacteria
Inhibition Zone (mm) at 5% Concentration
Chlorhexidine
No. Bacteria 1 2 3 4 5 6 7 8 9 10 11
(2%)
M±SD M±SD M±SD M±SD M±SD M±SD M±SD M±SD M±SD M±SD M±SD
M±SD
S. mutans
1 nd 7.5±0.1 8.9±0.1 8.3±0.05 7.3±0.1 7.7±0.1 7.2±0.1 - - - - 17,8±0.0
ATCC 25175
S. sanguinis
2 nd 9.5±0.1 9.0±0.15 9.0±0.1 8.9±0.05 9.6±0.1 9.7±0.1 8.6±0.15 8.4±0.05 8.7±0.1 8.8±0.15 21,9±0.0
ATCC 10556
E. faecalis
3 nd nd 13.4± 14.4±0.1 nd nd 8.6±0.1 nd nd nd 10.9±0.15 17,25±0.0
ATCC 29212
Abbreviation: nd (Note Determine), mm (millimeters), - (no inhibition), M (mean), SD (Standard Deviation)
Trends Sci. 202x; xx(xx): xxxxx 20 of 22
Table 5. MIC and MBC test results of Anhuienoside E (1) against three oral pathogenic bacteria
MIC (μg⁄mL) MBC( μg⁄mL)
No. Bacteria
M±SD M
1 S. mutans ATCC 25175 625±0.05 1000
2 S. sanguinis ATCC 10556 625±0.05 1000
3 E. faecalis ATCC 29212 625±0.1 -
Abbreviation: MIC (minimum inhibition concentration), MBC (minimum bactericidal concentration), - (inactive) M (mean), SD (Standard Deviation)