Computational approach to determine the combination of polyherbs based on the interaction of their metal complexes on the mucoadhesive properties of type II mucin

Mucin is the main mucus component besides water, which comes from the glycoprotein group. Mucoadhesive properties can be analyzed based on the binding interactions between compounds and mucin in mucus. This research is an in-silico approach aiming to obtain bioactive compounds that have a stable interaction with mucin-type II (MUC2) by screening ten compounds used in Indonesian herbal medicines. Molecular docking was used as a procedure for forming the compound-mucin complexes. Furthermore, molecular dynamic simulation with YASARA Dynamics software was used to obtain the stability of compound-mucin complexes. The curcumin-metal ions and kaempferol-metal ions complexes showed the most stable binding affinity and interaction with MUC2 compared to eight other ligands. Based on the molecular dynamic simulation, kaempferol-metal ions complexes are predicted to interact in the D3 chain, while curcumin-metal ions complexes more stably interact with the CysD1 chain. Both complex interactions simultaneously strengthen and synergize interactions with MUC2. Herbs containing one or both bioactive compounds, curcumin, and kaempferol are supposed to be composed into a polyherbal composition to obtain optimal mucoadhesive properties for MUC2


INTRODUCTION
The mucus has several crucial functions. The first function is as a lubricant because it can keep the epithelium hydrated and help with the transport and mobility of food (Ali and Pearson, 2007). In addition, this type of mucus is also a permeable gel layer where gas and nutrient exchange occurs. Mucus is a slippery secretion that coats the epithelial surface in the middle ear, eyes, respiratory tract, digestive tract, and urogenital tract (Linden et al., 2008). The mucus consists of two layers. The first layer, with a depth of about 7 µm, is the periciliary fluid layer (also known as the solid/sol phase) which lies above the epithelial cilia (Rubin, 2010). The second layer covers the slime layer, known as the gel phase. Mucus gels generally have a neutral pH with a thickness varying from 5 µm to 50 µm. Mucus has a short lifespan as it is replaced every 10 to 20 minutes (Lai et al., 2009). Mucus is a complex consisting of ± 97% water (Lu and Zheng, 2013). The remaining 3% consists of enzymes, lysozyme, defensins, cytokines, inorganic salts, proteins, immunoglobulins, glycoproteins (known as mucin), lipids, and cellular debris (Vareille et al., 2011).
Mucins, which are included in 3% of the main components of mucous membranes, are a type of glycoprotein with a molecular weight ranging from 2 to 50 MDa. Mucin comprises 80% of carbohydrates, namely, N-acetylglucosamine, traces of mannose, fucose, N-acetylgalactosamine, galactose, sulfate, and sialic acid (N-acetylneuraminic acid). This carbohydrate component is a signaling mediator between mucins and external components such as pathogens and drugs. These oligosaccharide chains are attached to the protein core, especially serine and threonine hydroxyl side chains (Wagner et al., 2018).
The protein core itself is organized into different regions. A vital characteristic of mucins is a proline-threonine-serine (PTS) rich domain in the central glycosylation region, which contains residues that form glycosidic links with N-Acetylgalactosamine (GalNAc). PTS domain is expressed as a repeating tandem, resulting in a domain that carries large numbers of glycan. Furthermore, the second type region between the PTS repeats has a small number of N-glycosylation and O-glycosylation sites. This region is known as the cysteine-rich domain because the cysteine content in this region is more than 10%. Due to this structure, several interactions are present and define mucin properties, including hydrogen bonding, electrostatic bonding, and hydrophobic interactions (Yang et al., 2012). Molecular docking analysis as an in-silico approach helps explain how interactions are formed and the types of bonds.
Various ethnobotanical and ethnopharmaceutical studies have made a list of medicinal plants used by the Indonesian people and their formulation and use, mainly based on ethnicity and regional origin. Several herbal medicines that have been studied for their efficacy and are composed as polyherbal include asam jawa (Tamarindus indica), tembelekan (Lantana camara), rosella (Hibiscus sabdariffa), kencur (Kaempferia galanga), sirih merah (Piper crocatum), kelor (Moringa oleifera), meniran hijau (Phyllanthus niruri), temulawak (Curcuma xanthorrhiza), jeruk nipis (Citrus aurantifolia), and jahe (Zingiber officinale) (Aristya et al., 2018;Fatimah et al., 2017;Grosvenor et al., 1995;Jadid et al., 2020;Kristianto et al., 2022;Silalahi et al., 2015;Puspitarini et al., 2022;Wijayakusuma, 2008). This study investigated the interaction mechanism between Indonesian herbal medicine and mucin in the digestive tract. Herbs and drugs in solution for the treatment of pathological conditions are most often consumed orally; therefore, herbal preparations should support mucoadhesive properties to interact or initiate the penetration process with mucosal components, especially mucin. Mucoadhesive properties are strong adhesion binding to the mucin component that maintains an attractive interaction with each other for an extended period based on interfacial forces (Smart, 2005).
The pathological mechanisms associated with immune escape have also been widely discussed regarding metalloprotein action. Examples are the enzymes SOD (superoxide dismutase) and CAT (catalase) (Agoro and Mura, 2019;Liao et al., 2013). From the point of view of bioenergetics and electron transfer systems, this pathogenesis mechanism can occur due to uncontrolled sequestration and acquisition of metalloprotein cofactors in the transition metal ions form. The study of roles and mechanisms played by non-carbon (inorganic) elements in a life process (biochemistry) from a multidisciplinary point of view is called bioinorganic (Kaim et al., 2013). Thus, this study constructed bioactive compounds from herbs used as samples bioinorganically in complexes with transition metal ions. This bioinorganic form will be closer to the native condition of bioactive compounds from herbs, commonly called low molecular weight (LMW) (O'Brien and Nunn, 2001). As an antioxidant defense and natural scavenger, especially an indirect-acting and preventive mechanism, LMW needs to interact with transition metals as cofactors (Engwa, 2018).
In comparison, radical scavenging of a single LMW will form new radicals as LMW loses its electrons (Ahmadinejad et al., 2017). As a potential chelator, LMW complex bioinorganic has prospects as a pharmaceutical agent, especially in energy transfer systems Sumitro and Sukmaningsih, 2018). Besides that, transition metals such as Fe, Cr, Zn, Co, Mn, Ni, Cu, Mg, and Ca are also required in small amounts to maintain the homeostasis of the human body, while either a deficiency or an excess of essential metals may result in various diseases (Bencini et al., 2010;Gupta, 2018;Jomova et al., 2022). A transition metal can be defined as any set of metallic elements occupying the central block in groups IVB, VIII, IB, and IIB, or 4-12 in the periodic table. The characteristics of this group exhibit variable valence and a strong tendency to form coordination compounds or complexes; many of their compounds are colored. Examples of the set of transition metals that are also the focus of this study are Fe (iron), Cu (copper), Mn (manganese), and Zn (zinc). The four metal elements are also included in the four elements that have crucial activities and roles in human biochemical metabolic processes and several pathological mechanisms. Metal ions are involved in catalytic mechanisms and play an active role in maintaining stability/complementing metalloproteins' tertiary and quaternary structures. Iron or Fe is the most abundant essential transition metal in living things, with 2-6 g of Fe found in normal adult humans (Alexander and Kowdley, 2009). In the cell, Fe forms in the form of ions (Fe +2 or Fe +3 ), where the higher the oxidation state is, the more it will initiate acceleration or shorten the reaction in the enzyme catalytic cycle (Andreini et al., 2007). Fe binds to metalloproteins as a cofactor and is also found in many cells. Examples of metalloproteins with Fe cofactors are the most widely known heme or iron-sulfur clusters which have functions in oxygen transport. Furthermore, copper (Cu) is commonly found as a cofactor of various proteins in living systems (Andreini et al., 2008). The protein needs to bind Cu ions to reach its active form. However, the free ion concentration of Cu needs to be kept at the essential threshold, which is a maximum of 100 mg, because it is toxic (Cox, 1999). In addition, transition metals tend to catalyze the formation of free radicals that can easily damage cell membranes (Vulpe and Packman, 1995). Meanwhile, manganese (Mn), with average concentrations in the adult human body as much as 12-20 mg, negatively correlates with Fe because they are both transported by DMT1 (Brandt and Schramm, 1986;Tuschl et al., 2013). Manganese metal is also a cofactor of mitochondrial superoxide dismutase (SOD). Zinc metal (Zn) is a constituent of living systems other than Fe, with an average concentration of 1.4-2.3 g (Aggett and Harries, 1979). This metal also dominantly functions as an enzyme cofactor with regulatory, catalytic, and structural activities. An example of a metalloprotein that has a Zn metal cofactor is cytoplasmic SOD together with Cu metal (Ferreira and Gahl, 2017). Therefore, it is crucial to maintain complex bioinorganic forms. Thus, this study used bioactive compounds complex with transition metal ions.
The proper mucin conformation is affected by mucin source, ionic strength, pH, and purity (Balabushevich et al., 2018). Therefore, the target mucin used in this study is gel-forming mucin commonly found in the digestive tract, namely, MUC2 (Corfield, 2018;Johansson et al., 2013). Furthermore, this study also analyzed the stability of the bioactive compound in bioinorganic form with mucin, using molecular dynamics simulation in physiological pH (pH 7.4).

Preparation of MUC2 as a protein target
The three-dimensional crystal structure of Homo sapiens MUC2 was retrieved from Protein Data Bank (PDB ID: 7A5O) (Javitt et al., 2020). The MUC2 3D structure was prepared by removing water molecules and ions using Discovery Studio 2021 Client. Ligand NAG (N-Acetyl Glucosamine) was not removed to validate the ability of herbal bioactive compounds as competitors of NAG to form electrostatic bonds.

QSAR analysis and drug-likeness properties
A quantitative structure-activity relationship or QSAR analysis was carried out using the Pass Online server from Way2Drug (Lagunin et al., 2010;Filimonov and Poroikov, 2008) to predict the potential bioactivity of each ligand. Drug-likeness properties were determined for each ligand with Lipinski's rules of 5. These methods were done to ensure that all ligands meet the requirements of drug-likeness properties to be used as drug molecules. These properties are analyzed using the Protox II database (Banerjee et al., 2018;Drwal et al., 2014) and ADMETLab 2.0 database (Dong et al., 2018;. Input for both databases was the canonical SMILE of each ligand.

Molecular docking
All molecular docking was performed using AutoDock Vina v.1.1.2, integrated into PyRx v.09.8, operating under Microsoft Windows 10. A rigid docking method was used with 50 exhausted and 9 modes parameters. The size of the D3 backbone protein grid box was adjusted to the binding position of the NAG, with the center xyz (312.596; 123.177; 211.069) and dimensions xyz (14.757; 18.817; 20.541), while the CysD1 protein grid box uses maximized values with center (176,539; 182,435; 192,295) and dimensions (53,644; 31,481; 36,168). The program calculated the binding energy using the Vina scoring formula [1]. This formula can be interpreted as a sum of the intermolecular part of the lowestscoring conformation and the intramolecular contributions (Trott and Olson, 2010). Once docking was completed, the results were saved in PDB (.pdb) format. Post-docking results were visualized and interpreted using PyMOL v.2.3.2.1 (Schrödinger and DeLano, 2020) and Discovery Studio 2021 Client (Biovia, 2020). (1)

Molecular dynamics simulation
Molecular dynamics simulations were performed using YASARA Dynamic v.14.12.2 (Krieger et al., 2002), operating under Microsoft Windows 10. YASARA program makes preparation, minimization, equilibration, and production processes in one parameter file, namely, the md_run macro. Variables from the preparation step include the water solvent model TIP3P (transferable intermolecular potential) with a density of 0.997 g/l in cubic simulation cell (Sgrignani et al., 2018); periodic type of cell boundaries; water box cell size 20 Å larger than the protein; 0.9% ions of NaCl; forcefield AMBER03 12-6-4 (Duan et al., 2003;Panteva and Giambasu, 2015). The variables for the minimization steps include the number of steepest gradient descent minimization steps of 50 steps; pH and temperature were edited to physiological (pH 7.4 and 310K), and without any constrain (Cojocaru and Clima, 2019;Krieger et al., 2006). The variables of the equilibration step are carried out with an ensemble in the form of NVT and a protonation scheme where pKas are only predicted for the side-chains of Asp, Glu, His, and Lys residues, but neither for Tyr, Cys, and Arg residues nor for the N-and C-termini. Variables from the production step include a 2.5 fs timestep automatically saved as a .sim file every 25 ps and a running time of 20,000 ps (20 ns). After production, the step is an analysis run using two-parameter files, namely, the macro md_ analyze and md_analyzeres. The analysis of root means square deviation (RMSD), the RG, and potential energy was obtained by running macro md_analyze. Furthermore, macro md_analyzeres was run to analyze root means square fluctuation (RMSF).

RESULTS AND DISCUSSION
Analysis was carried out between single compounds selected based on their tendencies as fingerprint compounds and their ability to bind transition metal ions. This chelating process i<j can be determined by the presence or the absence of functional groups ROH (hydroxyl) and RCOOH (carboxyl) (Fig. 1). In addition, mucoadhesive properties were predicted by interpreting the interaction of the ligand with the two D3 and CysD1 chains of MUC2. The D3 chain is a part of the backbone protein rich in PTS (Pro-Thr-Ser) domains and interacts with carbohydrate groups, especially NAG. The CysD1 chain is a structure at the end of each mucin peptide containing many residues of the amino acid cysteine (Cys).
Lipinski's rule of 5 helps to distinguish molecules of drugs from non-drugs. Lipinski can predict the probability of a compound's success or failure metabolism based on its similarity to the drug. The criteria for compounds that meet Lipinski's rules have at least 2 of the seven properties (Lipinski et al., 1997). The molecular mass of the compound is less than 500 Dalton. The compound is highly lipophilic (Log p value of less than 5). Then it has not more than ten hydrogen bond acceptors type, not more than five hydrogen bond donor type, and less than ten rotatable bonds. TPSA should also be between 0 and 140. At the same time, molar refractivity should be between 40 and 130. All variables used as Lipinski rule relate to the tendency of a drug/herbal compound to interact with macromolecules such as target proteins in cells as lipophilic environments, compared to its tendency to interact with solvents (Pollastri, 2010). Therefore, it is closely related to its absorption, distribution, metabolism, excretion, and toxicity (ADMET). All ligands meet Lipinski's rule of 5 (Fig. 2a) and are included as safe compound because it has a toxicity class ≥3 (Fig.  2b). The results of the ADMET analysis showed that the bioactive metal complex would potentially increase bioavailability, decrease Lipinski's rule violations, increase the L50 value, and reduce the level of toxicity of a compound compared to its single form.
The results of the QSAR analysis showed that almost all herbal fingerprint compounds, except moringyne, were predicted to have potential as mucomembranous protectors. However, carvacrol, phyllanthin, citric acid, and gingerol were predicted to have significant activity as mucositis treatment, while geraniol, carvacrol, moringyne, and gingerol compounds are also predicted to have activity as membrane permeability enhancers. Bioactivities such as steroid n-acetylglucosaminyl transferase inhibitor, TNF expression inhibitor, free radical scavenger, antioxidant, antiinflammatory, lipid peroxidase inhibitor, peroxidase inhibitor, and MMP9 expression inhibitor are also likely to be possessed by almost all compounds. However, some bioactivity is predicted to show optimal values only for two bioactive compounds. These bioactivities are intestinal-specific anti-inflammatory (curcumin and gingerol), immunosuppressant, and immunostimulant (geraniol and moringyne).
Meanwhile, when viewed from the overall average bioactivity, compounds with potential Pa (probability to be active) values above the threshold include geraniol, hibiscetin, kaempferol, carvacrol, curcumin, epicatechin, and gingerol (Fig.  3). The value of Pa describes the potential bioactivity of a sample being tested. Pa value analyzed above 0.7 indicates that the compound has significant potential and high similarity with the compound in the database that has been proven as the treatment in the wet-lab experiment. We use a score of 0.5 as the cut-off score for the minimum potential for bioactivity. The Pa value means the accuracy of the bioactivity prediction. The better the level of accuracy displayed is, the higher the Pa value of the bioactivity will be (Lagunin et al., 2010;Filimonov and Poroikov, 2008).

Molecular docking analysis
Molecular docking analysis was performed on two essential domains of MUC2, the D3 domain composed of PTS and NAG glycosylation and the CysD1 cysteine-rich domain (Yang et al., 2012). The molecular docking results show that almost all compounds that have formed complexes by binding to transition metal ions have a much lower binding affinity than a single compound when interacting with the two mucin chains. In addition, NAG used as a control compound required a relatively large binding affinity, videlicet −6.2 kcal/mol for the D3 chain and −4.6 kcal/mol for the CysD1 chain. The compound hibiscetin is predicted to have a stable interaction with the CysD1 chain but is less stable in interacting with D3. Meanwhile, kaempferol, moringyne, curcumin, and epicatechin showed the opposite results, with more stable binding to D3 than CysD1 based on the value of binding affinity (Fig. 4). This prediction is based on a minimum binding affinity value of −7 kcal/mol as the threshold (Trott and Olson, 2010).

Molecular dynamic analysis
Molecular dynamics simulation shows the dynamic stability of single and complex ion compounds when interacting with mucin in the human digestive tract. This molecular dynamic simulation was chosen as one of the drug discovery methods for predicting protein-ligand interactions by including the dynamic properties and flexibility of the target protein as the measured variables (Liu et al., 2018). The AMBER03 forcefield type used for simulating biological materials (Duan et al., 2003) is suitable for a secondary structure peptide, especially intrapeptide bonds (Smith et al., 2015), while the AMBER03 12-6-4 version is used because it has been optimized as a force field for monovalent, divalent, trivalent, and tetravalent ions, especially Mn 2+ , Zn 2+ , Cu 2+ , and Fe 3+ (Panteva and Giambasu, 2015). The samples subjected to molecular dynamics analysis were the two most potent complexes, curcumin-metal ions, and kaempferol-metal ions. This simulation also uses 2.5 fs of timestep instead of 5 fs because when the timestep for nonbonded forces is larger, hydrogen angle vibrations become critical at 2.5 fs (Krieger and Vriend, 2015).
The total potential energy analysis is obtained after running the md_analyze macro. Furthermore, it generally shows that potential energy will experience a very significant increase from 0 to 0.6 ns running time. This pattern indicates an energy initiation process to achieve energy stability. However, fluctuations began at 0.25 ns of running time, indicating that both samples changed the molecular bond energy. Even though it has reached a stable potential energy range (equilibrium phase), this fluctuation continues. Up-fluctuations can be interpreted as a strengthening of  bonds in molecules, while on the contrary a decrease in potential energy can be interpreted as a result of the relaxation of molecular bonds (Krieger et al., 2006).
In addition to the potential energy, the complex's stability can also be analyzed based on the RMSD value. RMSD is a score that provides information on conformational changes in a macromolecule that acts as a receptor after interacting with a particular ligand. RMSD can also be used as a standard deviation of conformational changes, with standards < 2 Å and > 2 Å generally being applied to docking results (Trott and Olson, 2010). The dynamic stability in question is the absence of significant conformational changes, better known as the unfolding process. The RMSD standard for a simulated protein receptor is 3 Å that the protein has undergone a conformational change (Smith et al., 2015). The pH treatment does not change the fluctuation pattern, and it only increases the average value of the distance in Å. The MUC2 chain D3 complex showed stability better if the control in the form of NAG binds to the PTS domain. Furthermore, the addition of curcumin-metal ions to the MUC2-NAG complex will significantly reduce the total RMSD value. However, the simulation results of the kaempferol compound bound with mucin showed the most unstable complex with average and maximum RMSD values of 5.413 Å and 9.976 Å (Fig. 5a). Meanwhile, the MUC2 chain CysD1 complex showed better stability if there was no interaction with curcumin or kaempferol although kaempferol showed a lower mean total RMSD value (Fig. 5b). However, all samples had a maximum total RMSD value not exceeding 3.5 Å. Therefore, it could be predicted that bonding between the mucin and the tethered compound would not be released.
The comparison of curcuminutes-metal ions complexes shows the order from the most stable to the most likely to be apart from the complex form with D3 chain after 20 ns simulation are three samples (single curcuminutes, CcMn, and CcFe) have almost the same stability, CcZn, and CcCu (Fig. 6a). Curcumin in a complex form with Fe and Mn has the best binding affinity (−7.9 kcal/mol) compared to curcumin without metal ions and other curcumin-metal complexes. Meanwhile, comparing the kaempferol-metal ions complex showed that the order KmCu, KmFe, KmZn, KmMn, and single kaempferol became the least stable complex (Fig. 6b). These results indicate that the complex kaempferol metals are stronger and more stable than kaempferol per se. Fluctuations values above 4 Å in the compound-ion complex are predicted to release the ions in the complex to find the most stable conformation. Next, the comparison of curcumin-metal ion complexes shows the order from the most stable to the most likely to be apart from the complex form with CysD1 after 20 ns of simulation are CcFe, CcMn, CcZn, single curcumin, and CcCu (Fig. 6c). Fluctuations in the curcumin-Cu complex that binds to CysD1 show an increase in the average number of RMSD close to 4 Å, indicating that Cu ions are no longer bound to curcumin. Meanwhile, comparing the kaempferol-metal ions complex showed the order KmCu, KmZn, KmFe, single kaempferol, and finally KmMn (Fig. 6d). These results show a consistent pattern of complex stability similar to the interaction between the kaempferol-metal ion and the CysD1 mucin chain.
The YASARA program's RMSD value can also be used to see the RMSD of the ligand configuration, which in this study consisted of two single compounds (curcumin and kaempferol) and each complex structure with four types of transition metal ions. Ligand configuration RMSD was calculated by superposing the ligand with its reference structure as a function of the simulation period. The lower value of ligand configuration RMSD usually indicates proximity to the reference and hence better ligand binding affinity (Odhar et al., 2021). The RMSD ligand value above 4 Å indicates a release of ligand, especially the metal ion portion chelated by the compound (Smith et al., 2015). The RMSD of the curcumin ligand configuration was in line with the total RMSD values.
Meanwhile, in the kaempferol ligand configuration, different results between the RMSD ligand configuration and the total RMSD were seen in the KmCu (chain D3) and KmZn (CysD1 chain) complexes. This result indicates that the transition metal ion Cu tends to make the complex interactions of Cu compounds with the D3 MUC2 chain of the electrostatic type more unstable (Fig. 7a), while in the CysD1 chain, it can be predicted that the release of Zn ions from the kaempferol complex after molecular dynamic simulations was due to the RMSD value of the ligand configuration, which exceeded the value of 4 Å (Fig. 7b).
A comparison of the results before and after the simulation for 20 ns was performed by superimposing the complex at 0 ns (colored dark blue-red ligand) with the sample at 20 ns (colored dark green-white ligand). The single curcumin compound bonded to MUC2 chain D3 did not experience significant changes in position and conformation (Fig. 8a). The curcumin-Cu compound that bonds to the MUC2 chain D3 makes the NAG ligand shift position so that it successfully forms electrostatic bonds, which are also needed in the signaling process for exogenous compounds from the diet to the intestinal flora. The compound curcumin complexed with Fe, Mn, and Zn binding to the MUC2 chain D3 made the NAG ligand shift position but failed to form an electrostatic bond. Meanwhile, the single kaempferol compound that binds to the MUC2 chain D3 experienced a significant change in position and conformation, thus supporting the average total RMSD results (Fig. 8b). The kaempferol-metal ions compound attached to the MUC2 chain D3 causes the NAG ligand to shift its position and form an electrostatic bond (Fig. 8c-g).    electrostatic bonds of Pi-cation, and GLU will form Pi-anion. Cationic environments in the Pi-cation bond are provided by the protonation of primary arginine, lysine, and histidine groups. At the same time, deprotonation of the carboxylic acid side chain will occur in GLU or aspartic acid residues to form an anionic environment (Schaeffer, 2008).
The metal ion successfully initiates the formation of an electrostatic bond, which shows that the bioinorganic form of the bioactive compound is a functional form that naturally has more potent activity than the single form without metal ions. This natural form is a LMWC. LMW can be divided into primary and secondary metabolites based on molecular weight characteristics below 900 Da (Macielag, 2011). Secondary metabolites of LWM are a group of bioactive compounds commonly studied as medicinal compounds. Transition metals such as Mn, Cu, and Fe can catalyze the electron reduction process of molecular oxygen. Hence, LMWCs usually act as organic ligands in coordination complexes formed by various transition metals as central atoms (Hadacek and Bachmann, 2015).
Furthermore, the interaction between LMW and inorganic transition metals as central ions can initiate chargetransfer complex formation. The formation of charge-transfer complexes occurs when LMW and protein mucins' adjacent molecular orbitals as targets overlap (Szent-Györgyi, 1960;Szent-Györgyi, 1968). Charge-transfer complex formation in protein chemistry can be grouped as a variant of the dipole-dipole interactions (Silverman, 2002). Dipole-dipole interaction is a form of the electrostatic force between two permanent molecular dipoles (Pal, 2020). Hydrogen bonding is a dipole-dipole interaction formed between the proton of a group electronegative atom (X) and atom hydrogen (H), while there are also other electronegative atoms (Y) that contain a pair of nonbonded electrons (Silverman, 2004). Another example of an electrostatic interaction is Pi-anion and Pi-cation. Pi-Anion interactions are favorable non-covalent between an anion and an electron-deficient (Pi-acidic) aromatic system (Schottel, 2008). Opposite to Pi-anion, Pi-cation also forms favorable interaction between the cation and an electron-rich aromatic system. Many studies have analyzed Pi interactions, especially Pi-cation, which can enhance binding energies by 2-5 kcal/mol. Hence, hydrogen bonds and attractive charge interaction (ion pairs) are crucial factors in drug-receptor and protein-protein interactions (Dougherty, 2013). Pi interactions also play an essential role in nature as enzyme catalysis, protein structure, and molecular recognition. Molecular recognition also includes the chemical signaling of exogenous compounds and a NAG assignment on the MUC2 PTS domain. Based on these interpretations, it can be predicted that the kaempferol complex is more stable and has better mucoadhesive properties than the D3 chain because it has succeeded in forming electrostatic bonds for all types of ligand complexes than the curcumin complex. However, due to the tendency of metal ions to change the binding affinity of the compound, several types of unfavorable bonds were also formed from the ligand interaction analysis of the KmCu complex (unfavorable bump and unfavorable donor-donor), KmMn complex (unfavorable donor-donor), and KmZn complex (unfavorable donor-donor and unfavorable positive-positive). This type of bond can reduce the stability of the interaction between the ligand and macromolecules according to the results reflected in the total RMSD average value and the ligand configuration RMSD.
Furthermore, comparisons on the CysD1 chain focused on the hydrophobicity of the positions of the amino acid residues interacting with the ligands. Interpretation results show that adding transition metal ions Zn, Mn, and Cu will increase the tendency of curcumin to interact with hydrophobic amino acid residues, while Fe ions tend to interact with hydrophilic amino acid residues lysine (Fig. 9b). However, Cu (Fig. 9a) and Zn (Fig. 9c) ions are unstable and tend to be apart from the complex after simulation. The addition of transition metal ions Mn, Cu, and especially Zn will increase the tendency of kaempferol compounds to interact with mainly hydrophilic amino acid residues. In contrast, Fe ions interact with hydrophobic amino acid residues such as leucine (Fig. 10b). In addition, similar to the curcumin complex, the Zn ion is unstable and tends to escape after simulation (Fig. 10c). However, Cu ions form covalent bonds with kaempferol, which align with the total RMSD and RMSD ligand configurations (Fig. 10a). Thus, curcumin generally forms more bonds with hydrophobic amino acids in the CysD1 chain than kaempferol to form hydrophobic pockets. This hydrophobic pocket makes the compound environment more hydrophobic and protects the compound to maintain its bioavailability.
The prediction that the kaempferol-metal ions have a higher binding affinity for the D3 domain was also confirmed by the results of RMSF (Fig. 11a) and RG (Fig. 11c). The RMSF showed that kaempferol binding could significantly maintain the bound amino acid residue (in the dotted box), especially in the results of the kaempferol-metal ions complex, which was used as a ligand. This stability can be seen in the average RMSF, which is still below 5 Å as the cut-off value (Krieger et al., 2006). These results align with RG results. The kaempferol-metal ions complex significantly maintains the consistency of the RG pattern of MUC2, which is much more stable than the MUC2-kaempferol complex. RMSF is a score that provides information on conformational changes in more detail because it is associated with fluctuations in the level of amino acid residues. High RMSF values indicate ligand mobility, while low RMSF suggests the interacted ligand stability with MUC2 (Dash et al., 2017). The value of RG and RMSF analysis helps explain that the formation of an unfavorable bond in the kaempferol-metal ions complex with MUC2 only affects the instability of the ligand bond, according to the results of RMSD and RMSD ligand configurations. However, interactions from hydrogen, electrostatic, hydrophobic, and van der Waals help maintain the complex's stability, especially the amino acid residues of MUC2 that bind the ligands.
Meanwhile, curcumin, predicted to have a better binding affinity for the CysD1 domain, maintained conformational stability with an RMSF mean value below 4 Å and a consistent RG pattern (Fig. 10b and d). Furthermore, the RG is a parameter that describes the conformation equilibrium for all the various simulations performed. The RG value, which can also be explained as the radius of rotation of the dynamic movement of a complex, both proteins and protein compounds to the solvent, thus becomes one way of predicting the simulation of the solubility of the sample in a liquid solution or solvent (Lobanov et al., 2008). The lowest value indicates the folded protein condition, while the highest value indicates the protein conformational condition when unfolded (Yamamoto et al., 2021). Therefore, RG is commonly used to predict the dynamics pattern in a solvent environment.  Based on the overall results, it is predicted that kaempferol and curcumin tend to interact with MUC2 in the D3 and the CysD1 chain, respectively. The interaction of both will strengthen and synergize polyherbal interactions with MUC2. In addition, the bioactive compound kaempferol tends to interact with the D3 chain of MUC2 compared to the CysD1 chain, reducing the possibility of crosslinking (Menchicchi et al., 2015). Chain D3 also has a more hydrophilic surface than chain CysD1 to maintain the viscosity of MUC2. Besides having high affinity, both bioactive compounds can significantly increase or maintain conformational stability and dynamic pattern in the bioinorganic form. The conclusion of this study was obtained from only one fingerprint compound, which represented the complex compounds contained in herbs. Therefore, it is necessary to do an in-vitro study of the mucoadhesive properties of Indonesian herbs.
Kaempferol is widely distributed in different genera like Berberis, Camellia, Allium, Citrus, Delphinium, Brassica, and Malus (Devi et al., 2015). Kaempferol is naturally bonded to different glycoside moieties (Montano et al., 2011). This study shows that kaempferol has potential as a mucoadhesive component aided by its ability to form interaction with glucan groups such as NAG. Kaempferol is also commonly known as a flavonolderived and fingerprint compound from Kaempferia galanga L. rhizome (Khairullah et al., 2020). Kencur (Kaempferia galanga L.) is one of the Indonesian medical plants from the family Zingiberaceae, distributed throughout the world (Mitra et al., 2007;Ling, 2009;Techaprasan et al., 2010;Sukenti et al., 2016;Purba et al., 2018). Kencur also has been mentioned to have many bioactivities such as a stimulant, analgesic, expectorant, diuretic, antipyretic, antioxidant, anti-inflammatory, anticholelithiatic, and antimicrobial (Umar et al., 2011). This provides the results of the QSAR analysis. Compounds with choleretic effects, which are also equipped with potential anti-inflammatory, antioxidant, and antimicrobial activities, will interact with gel-forming mucin to treat intestinal wounds (Nazar et al., 2008). However, no research still evaluates the mucoadhesive properties of kaempferol in compound form and in complex form as kencur in more detail. Thus, conducting a wet lab study is necessary to explore the mucoadhesive properties of kaempferol both as pure extracts and herbal components. Meanwhile, curcumin, a fingerprint compound of curcuminoid from Curcuma longa extract, has been investigated for its mucoadhesive properties (Ortiz et al., 2019). However, the binding mechanism of curcumin to the mucin structure is not precise.
Several studies have also proven that plant herbal medicine contains transition metals as one of the micronutrients absorbed with groundwater. Jahe (Zingiber officinale) and meniran (Phyllanthus niruri L.) are proven to contain 6-gingerol and shogaol as phenolic compounds that bond with Fe and Mn (Kristianto et al., 2022). Zaitun (Olea europaea) is also proven to contain oleuropein-Cu complex (Capo et al., 2017). Brotowali (Tinospora crispa) contains various metals such as Fe, Cu, and Mn, significantly predicted in choline-Fe and DL-Carnitine-Fe complex molecules (Widodo et al., 2021). Duwet (Syzygium cumini) is proven to contain flavonoid-Fe complexes, especially in the cyanidin-Fe complex molecule (Sukmaningsih et al., 2018). Flavonoid complexes with transition metals Fe and Cu were of more vigorous radical scavenging activity than their single form (Jabeen et al., 2017;Sukmaningsih et al., 2018).

CONCLUSION
Herbs containing one or both bioactive compounds, namely, curcumin and kaempferol, are suggested to be composed into a polyherbal composition to obtain optimal mucoadhesive properties for MUC2. One of the suggested compositions based on this research is polyherbal composition formed from Kencur (Kaempferia galanga) with Temulawak (Curcuma xanthorrhiza). In addition, the preparation or manufacture of herbal medicine that maintains the availability of transition metal content is recommended because it is predicted to have much better bioactivity. It is crucial to preserve the natural form of LMW compounds in complexes formed with inorganic transition metals as the central ion. The complex form of these metal-ionic compounds can initiate charge transfers, which are crucial in signaling processes, functional activities, mucoadhesive properties, and other metabolic processes.

AUTHOR CONTRIBUTIONS
All authors made substantial contributions to conception and design, acquisition of data, or analysis and interpretation of data; took part in drafting the article or revising it critically for important intellectual content; agreed to submit to the current journal; gave final approval of the version to be published; and agree to be accountable for all aspects of the work. All the authors are eligible to be an author as per the international committee of medical journal editors (ICMJE) requirements/guidelines.

CONFLICTS OF INTEREST
The authors report no financial or any other conflicts of interest in this work.

ETHICAL APPROVALS
This study does not involve experiments on animals or human subjects.

DATA AVAILABILITY
All data generated and analyzed are included in this research article.

PUBLISHER'S NOTE
This journal remains neutral with regard to jurisdictional claims in published institutional affiliation.