Phytoconstituents of Butterbur (P. japonicus), their metabolic pathway and ability to modulate bone morphogenic protein (BMP) signaling

ABSTRACT A library of natural sesquiterpene and phenolic compounds from Petasites japonicus are being investigated through different computational techniques to study their ability to target BMP. Lipinski rule, ADMET, molecular docking studies and metabolism were used to reach promising candidates with proposed activity against BMP. Four sesquiterpenes (kablicin, petisinol, bakkenolide D and bakkenolide IIIa) and four phenolic compounds exhibited drug-like properties (caffeic acid, petasiphenol, petasitesin A and petasitesin B), so they deserve further clinical exploration as bone loss modulators. The phenolic compounds specially fukinolic acid and petasiphenol showed lower binding energy with both BMPRIA and BMPRII than Icariin agonist and sesquiterpenes. Bakkenolide IIIa showed dual potential on both BMPRIA and BMPRII with binding energies equal – 7.82 and – 9.9 Kcal/mol respectively, which is more better score than Betulinic acid agonist. This research is focusing on plant-human interactions and exploring the ability of plant constituents to modulate a human protein such as BMP.


Introduction
Bone fractures are a major public health concern and a significant financial burden, particularly for those who have osteoporosis.Fractures are a significant burden on individuals, communities and healthcare systems because they can result in loss of productivity, incapacity, diminished quality of life, health decline and high health-care expenditures (Wu et al. 2021).Osteoporosis is one of the most prevalent metabolic skeletal illnesses; it is characterized by reduced bone strength, limited osteogenic activity, and increased osteoclast resorption, which increases the risk of bone fractures (Ferrari et al. 2018).Fractures caused by osteoporosis result in a considerable decrease in quality of life, as well as an increase in morbidity, mortality, and disability (Varacallo and Fox 2014).Numerous variables contribute to osteoporosis and fractures.Major predisposing factors are food and physical activity.Fracture history in the family and personal fracture history are both risk factors for osteoporosis.More than 8.9 million fractures occur annually due to osteoporosis, resulting in a fracture every three seconds (Johnell and Kanis 2006).Consequently, there is an urgent need to develop novel treatments for this disease.Osteoblasts facilitate the process of bone formation by secreting bone matrix and boosting calcium deposition.Nonetheless, osteoclasts mediate the remodeling process by resorbing bone matrix (Guo et al. 2021).Bone morphogenetic proteins (BMPs) are traditional growth factors that play crucial roles in the regulation of proper bone formation and bone resorption (Baek et al. 2021).BMPs regulate anabolic and catabolic mechanisms by controlling the development and activity of osteoblasts (OBs) and osteoclasts (OCs), the principal cell types responsible for bone production and resorption, respectively.BMPs are released from the hydrolyzed bone and transported by OC-like cells to stimulate osteogenesis, which is triggered by the recruitment of mesenchymal OB progenitors.Next, differentiating chondrocytes and activated OBs produce the organic bone matrix, which is then matured by the addition of inorganic salts.In mature calcified bone, OBs undergo apoptosis or develop into quiescent osteocytes (Segredo-Morales et al. 2018;Wang et al. 2021).More than fifteen known BMPs are structurally related and can be further classified based on amino acid or nucleotide similarity.Due to their pervasive expression and role as bodywide regulators, a loss in BMP synthesis or function typically results in pronounced abnormalities or severe diseases (Zhang et al. 2021).Usually, alterations in TGF-b/BMP activity are associated with a vast array of clinical consequences, including skeletal and extra-skeletal abnormalities, autoimmune disorders, cancer, and cardiovascular diseases (Tang et al. 2021) Some BMPs can act as bone formation inhibitors.BMP3, for instance, is a negative regulator of bone density, whereas BMP13 is a potent inhibitor of bone development.Based on their powerful bone-inducing activities, BMP2, 4, 6, 7, and 9 are usually referred to as osteogenic BMPs (Lademann et al. 2020).BMP-2 is one of the most osteogenic BMPs capable of considerably promoting bone formation.Although BMP2 is unnecessary for bone production, it is necessary for the onset of fracture healing and induced both activity of osteoblast and osteoclast (Lowery et al. 2011).BMP-2 deficiency can slow the transition of osteoprogenitor cells into osteoblasts (Wang et al. 2014).The Food and Drug Administration (FDA) has approved BMP2 for the treatment of long bone fractures and spinal fusion operations and recombinant human BMP2 (rhBMP-2) for human usage (James et al. 2016;Chao et al. 2021).BMP-2 can strongly promote bone growth in osteoblasts via the SMAD signaling pathway (Fitzpatrick et al. 2017).BMP2 can bind to the type I receptor (BMPRIa) and type II receptor, which are involved for beginning the canonical and noncanonical signaling cascades, respectively (BMPRII) (Figure 1).Signaling caused by BMP2 begins when the protein either recruits the dimerization of two receptors or binds to receptors that are already dimerized (Okamoto et al. 2006;Bragdon et al. 2011;Durbano et al. 2020).The type II receptor will phosphorylate the type I receptor's glycine-serine-rich (GS) box homeodomain once all three components have formed a complex (Mace et al. 2006).Consequently, SMAD-independent anddependent signaling may be stimulated (Yang et al. 2014;Wang et al. 2020).
Petasites japonicas (butterbur), is known as giant butterbur (Figure 2) and sweet-coltsfoot, the plant family is Asteraceae, it has medicinal significance as it has been used in treatment and prophylaxis of migraine, headache and tension in China, Japan and Korea, It has antiallergic effects, anti-asthma, anti-cancer, antioxidant activity, neuroprotective, Improvement of metabolism and anti-inflammatory properties (Lee et al. 2019;Hiemoriikondo 2020;Borlak et al. 2022).
Moreover, Petasites japonicus, a perennial member of the Compositae family that grows wild across much of Japan from Hokkaido to Kyushu, is known as the Japanese butterbur (Iwamoto 2009).Butterbur (Petasites japonicus) consists of herbaceous perennial plants with huge, hat-shaped leaves and thick, creeping underground rhizomes (Figure 2).
Sesquiterpene analogs have played a good role in bone biology as they had proved to play an important role towards the development of new and effective approaches to prevent and treat osteolytic bone diseases (Qin et al. 2015).Additionally, it has been reported that some sesquiterpenes can inhibit osteoclastogenesis (Qin et al. 2015).As instance, sesquiterpen lactone compound I has been shown to inhibit bone morphogenic protein (BMP) with an IC 50 value of 51 μM (Ohte et al. 2021).Natural Sesquiterpenes II-IV were also reported to Inhibit bone resorption and osteoclast formation (Qin et al. 2015).The natural phenolic compounds oleuropein and resveratrol (Figure 3) have important roles in bone formation.Puel et al. found that oleuropein reduced bone loss (Puel et al. 2006).Resveratrol has shown ability to inhibit the MAP kinase pathways in osteoblasts.Furthermore, resveratrol inhibited the expression of bone morphogenic protein-4 (BMP-4) (Nicolin et al. 2019).
Also, Petasites japonicus (Siebold & Zucc.) is one of Petasites species which is known with its content of sesquiterpenes and phenolic acids constituents (Kulinowski et al. 2022), so we decided in this research work to investigate the possible interactions of P. Japonicous sesquiterpenes and polyphenols with bone morphogenetic protein.
These findings encouraged us to investigate the natural sesquiterpenes and phenolic constituents from Butterbur as an alternative strategy to enhance BMP signaling by using these naturally extracted products.
Icariin is a flavonoid glucoside, and various studies have demonstrated that treatment with Icariin increases osteogenic differentiation and decreases bone loss in vivo and in vitro (Figure 3).Icariin augments BMP2-mediated osteoblastic development of myoblast cells (Chen et al. 2019).Betulinic acid is a pentacyclic triterpenoid that occurs naturally in a variety of fruits and vegetables.Betulinic acid can boost the osteogenic potentials of BMP2 (Figure 3) in vivo, possibly through activating the Smad 1/5/8 and p38 pathways, and the combination of both drugs can be regarded as a treatment for bone disorders (Choi et al. 2016).
In this research work, we will try to find new modulating agents for BMP-2 that may be useful in repairing bones.We decided to investigate some sesquiterpenes and phenolic compounds (Bradford et al. 2019) for their potentials to act as BMP-2 up-regulators.The selected library (Figure 4 and Figure 5) was subjected to different computational techniques as Lipinski rule, ADMET, toxicity and docking  studies to get new promising candidates with expected activity against BMP-2.

Drug-likeness properties and admet analysis
In this work, a library of 37 sesquiterpene compounds and 12 phenolic constituents of Butterbur (Petasites japonicus) compounds (Figures 4 and 5) were subjected to various computational techniques to investigate the effect of chemical structure changes on the ability of these two scaffolds to target BMP-2 receptors, leading to promising candidates with proposed agonist activity on BMP-2 receptors.

Lipinski's and veber's rules
As per the guidelines of Veber's and Lipinski's rule of five, a computer analysis was performed in this study to ascertain the number of rotatable bonds, topological polar surface area (TPSA), and other physicochemical features of the examined candidates.
2.1.1.1.Sesquiterpene.Table 1S (Supp Data) showed no violation with either Lipinski's or Veber's rules, indicating that they all have a potential for pharmaceutical use.Hydrogen bond acceptors (HBA) in the sesquiterpene compounds ranged from 0 to 6, which is below the threshold of 10 to meet the first parameter of Lipinski's rule.Furthermore, concerning the number of hydrogen bond donors (HBD), all investigated sesquiterpene compounds have HBD ranging from 0-2 (i.e. less than 5) which meets the second parameter of Lipinski's rule.The third parameter of Lipinski's rule is supported by the fact that all of these molecules have a molecular weight of less than 500.The logP values of these compounds were also quite varied, going from 1.70-4.63.The results are consistent with Lipinski's rule's fourth parameter, which states that logP values below 5 indicates promising formulated candidates.The molecular flexibility and oral bioavailability of the drug candidates were evaluated by determining the number of rotatable bonds using Veber's principles.The study showed that sesquiterpene molecules had between 1 and 7 rotatable bonds, which is within the range required by Veber's rules for such compounds.Furthermore, the number of TPSA a physicochemical characteristic defining the polarity of such molecules is less than 140 Å², which is within the allowed values (range from 0 to 104.2 Å²) (Table 2S, Supp Data).

Phenolic compounds.
There are four molecules with drug-like molecular structures (Table 1).Caffeic acid, Petasiphenol, Petasitesin A, and Petasitesin B all had sufficient numbers of rotatable bonds, HBA, HBD, logP, and TPSA values (Table 1, 3S and 4S, Supp Data).As demonstrated in Tables 3S and 4S, the other phenolic compounds showed some violations of Lipinski's rule or Veber's rule because either their molecular weight exceeds 500 or the number of HBA and HBD exceeds appropriate levels (Supp.Data).

Pharmacokinetic properties parameters
Blood-brain barrier (BBB) penetrability, intestinal absorption, and plasma protein binding (PPB) characteristics were determined for all compounds.In details represented in (Table 2, 1S, and 3S), 32 sesquiterpene compounds revealed ability to permeate BBB.Eleven Phenolic compounds, on the other hand, did not traverse BBB and were therefore predicted to be safe for the CNS.In general, all sesquiterpene compounds except eremophilene showed high intestinal absorption, and caffeic acid, petasiphenol, petasite-sinA, and petasitesin B exhibited high intestinal absorption.Furthermore, absorption percent (%ABS) was predicted using the equation %ABS = 109 − (0.345 × TPSA) (Zhao et al. 2002), 37 sesquiterpenes and 12 phenolic compounds had a range between 36.10% to 109%, indicating that the necessary cell membrane permeability of sesquiterpene compounds (73.05-109%) and caffeic acid (82.17%) was greater than that of other phenolic compounds (36.10-72%) (Table 2S and 4S, Supp.Data).Bioavailability, which is an indicator of the amount of drug content in plasma, is the most significant factor influencing absorption.Interestingly, the bioavailability of all sesquiterpene compounds and just four phenolic molecules was determined to be high.
Due to the fact that none of the sesquiterpenes compounds (with the exception of furanofukinol) are Pgp protein substrates, they have a low risk of effluxing out of the cell, resulting in a maximal effect.Five phenolic compounds out of twelve are not Pgp protein substrates.These include chlorogenic acid and caffeic acid, fukinolic acid, fukiic acid, and petasiphenol.
On the basis of the results obtained from previous physicochemical studies, the Pharmacokinetic properties of compounds that obey Lipinski's and do not pass the BBB, have high bioavailability scores and high gastrointestinal absorption were evaluated.Four sesquiterpenes compounds (kablicin, petisinol, bakkenolide D, and bakkenolide IIIa) and four phenolic compounds exhibit drug-like properties (caffeic acid, petasiphenol, petasitesin A, and petasitesin B).Using the Molsoft programme to predict the solubility value of the compounds, it was determined that the derivatives met the solubility criteria with values ranging from 19.91- 6782.68 mg/l, which is greater than 0.0001 mg/l.The Pre-ADME website was also utilized to assess oral drug delivery, central nervous system disposition, and plasma-protein binding (PPB) model efficacy (Thangavelu et al. 2017).Human colon adenocarcinoma-derived Caco2 cells were examined for their permeability.In the Caco2 cell model, the medium cell permeability scores for the tested compounds ranged from 18.07-35.51nm/s, as shown in Table 3.The most permeable molecule to MDCK (Mandin Darby Canine Kidney) cells in vitro was petasitesin A, while the least permeable compound was bakkenolide IIIa.Moreover, intestinal absorption (HIA) values for humans ranged from 67.60% to 99.22%, indicating that substances are highly absorbed.According to the PPB abilities investigation, Kablicin, Bakkenolide D, Bakkenolide IIIa, and petasitesin A revealed strong protein binding (above 90%), but petasinol caffeic acid, petasiphenol, and Petasitesin B showed low protein binding (less than 90%).The DataWarrior programme was used to investigate the mutagenicity and tumorigenicity profiles of eight compounds, none of which displayed mutagenicity or carcinogenicity with the exception of caffeic acid, which exhibited strong tumorigenicity (Table 3).

Docking studies
Docking studies serve to validate the model's quality and comprehend the binding manner of drugs to their biological targets in terms of affinity.Consequently, docking research was conducted to elucidate the potential binding manner of sesquiterpene and phenolic compounds with distinct activity profiles (Belal et al. 2022).Human BMP-2, in its mature form, is an acidic glycoprotein of about 32 kDa.The BMP2 ligand, a homodimer of 114 amino acid residues, has been found to include two receptor-binding motifs, rendering it as a multifunctional growth factor.A high-affinity binding site for BMPRIA is found in the wrist epitope, which contains residues from both BMP2 subunits, whereas a binding site for BMPRII is found in the knuckle epitope but is weaker (Rahman et al. 2015).In BMP-2, three intra-chain disulfide bonds are produced by six conserved cysteine residues.There are four glycosylation sites in the human BMP-2 proprotein.The mature form of BMP-2 contains a single glycosylation site in its polypeptide chain.Upregulation of the BMP-2 gene and the activation of multiple downstream signaling pathways result from BMP-2 binding to BMP type I and type II serine/threonine kinase receptors (Figure 6) (Miyazono et al. 2010).After binding to the BMPRs, BMP-2 triggers the development of adipocytes, chondrocytes, and osteoblasts via BMPRIA phosphorylation by BMPRII and apoptosis/cell death via BMPRIB phosphorylation (Ebara and Nakayama 2002).All compounds docked in the active sites of both BMPRIA and BMPRII.Table 2 showed docking scores of all compounds, Icariin and Betulinic acid.Figures 6-17 and Table 4, Table 5 and Table 6 displayed the binding mode of the highest docking scores of sesquiterpene and phenolic compounds.

BMPRIA receptor docking simulation
During embryonic development, the Bmp receptor activinlike kinase 3 (Alk3; or BMPRIA) is expressed in a wide variety of tissues.However, in the postnatal period, it is mostly expressed in osteo-lineage cells and bone marrow cells.Throughout osteogenesis and chondrogenesis, ALK3 has been repeatedly demonstrated to be a critical receptor for regulating BMP signaling (Mueller and Nickel 2012).ALK3 is one of the key factors for regulating the specification of pre-chondrogenic mesenchyme as well as chondro-lineage differentiation and maturation, postnatal chondrogenesis and the maintenance of articular cartilage (Lin et al. 2016).Two β-sheets and an α-helix compensate the extracellular domain of BMPRIA.Additional residues from monomer B and helix 1 make up the rest of the wrist epitope.The interaction between Phe 85 amino acid and BMP-2 ligand is a hallmark of type I receptors.Pro45 and His43 found in loop 1 of BMPRIA whose backbone flexibility is restricted among type I receptors (Lin et al. 2016).Phe 60, Met 78, and Ile 99 are highly conserved and functionally important residues in BMPRIA receptor.Asn 59, along with Ile 62, Val 63, and Leu 66 of BMP-2 ligand, provides an almost entirely hydrophobic surface that makes contact with receptor residues Glu 81, Phe 85, Gln 86, and Asp 89 (Keller et al. 2004;Mace et al. 2006).The docking results for BMPRIA revealed that phenolic compounds are superior to sesquiterpene by more than -8 kcal/mol, while the energy of binding of sesquiterpene is greater than -7.5 kcal/mol.Consequently, then we can conclude that both are promising regulators, however, sesquiterpene has a lower affinity for the BMPRIA receptor when compared with phenolic compounds.The best phenolic compounds were 3,4,5-tricaffeoylquinic acid, Fukinolic acid, and Petasiphenol because their docking scores were greater than Icariin's (BMP-2 agonist), see Table 5.The binding energy was used to calculate the binding affinity constant (Ki value) through the equation (ki = 10 [Binding Energy/ 1.366] ) (Farghaly et al. 2021) and the most promising target compounds were chosen based on their binding energy and interaction within the binding site, Table 5.The hypothesized binding mechanism of Icariin and Betulinic acid with the BMPIA receptor was depicted in Figure 7 and Figure 8. Furthermore, the binding energy of icariin was lower than that of Betulinic acid, -10.62 versus -7.49kcal/mole, indicating a greater affinity of Icariin for the BMPIA receptor than Betulinic acid.Icariin is anchored by multiple forms of hydrogen bonds with Asp89, Lys92, Ser90, His43, Phe85, and Gln86.Icariin also forms four hydrophobic bonds with Arg97 and Asp89 (Figure 6).Betulinic acid has one hydrogen bond and four hydrophobic interactions (Figure 7).
Bakkenolide IIIA, which has a favorable ADME profile and the highest (sesquiterpene compound) docking score of -7.82 kcal/mol at the BMPIA receptor, is held in place by hydrogen bonding connections with Ser90 and four hydrophobic contacts with the essential amino acids Phe85, Lys88, and Pro91 (Figure 9).
In the postulated binding mode of the phenolic compound 3,4,5-tricaffeoylquinic acid, two catechol moieties are anchored by hydrogen bonding interactions with Gln86, Asp89, and Asp84, as well as the formation of carbon hydrogen bonds and Pi-doner-hydrogen interactions with Asp89, and Phe85 (Figure 10).
The final results for Fukinolic acid and Petsiphenol with the BMPIA receptor are consistent with the general pattern observed for Icariin with high docking scores (−14.07, 13.90, and -10.62 kcal/mol respectively).The hydrogen bonding interactions continue to exist as before.Figure 11 shows that the di hydroxy phenyl moieties of Fukinolic acid interacted favorably with Gln86, Asp89, Pro91, Phe85, and Lys88.Also, in Petasiphenol, the two catechol rings may be important for the compound's accommodation at the binding site (Figure 12).

BMPRII receptor docking simulation:
BMPRII is a serine threonine kinase receptor of type II TGFb that is essential to the BMP signaling pathway.Crystals of the BMPRII ectodomain were produced in two distinct forms, and both forms diffracted with a high degree of resolution (Miyazono et al. 2010).The tetragonal form exhibited some disorder, whereas the orthorhombic form included the full polypeptide.They exhibit various conformations of the A-loop at the periphery of the hypothesized ligand-binding interface, as well as a disulfide bridge rearrangement within the loop.This specific disulfide (Cys94 -Cys117) is exclusive to BMPRII and activin receptors, suggesting that it is essential for their shared manner of binding.Cys99 and Cys116 are connected by an extra disulfide bridge in BMPRII.His87 has a crucial part in ligand recognition (Guo et al. 2022).Based on the docking scores and interactions with the binding site, the most promising natural compounds were selected for binding affinity predictions.Equation (ki = 10 [Binding Energy/ 1.366] ) was used to get the binding affinity (Ki value) using the binding energy scores (Farghaly et al. 2021).
Icariin binds to the BMPRII receptor via multiple interactions, particularly those with Ser119, Cys94, Pro91, and His87 that are of critical importance (Figure 13).Betulinic acid interacts via three hydrogen bonds and four alkyl interactions, similar to Icariin (Figure 14), but with a higher binding energy.The docking scores (S) for Icariin and Betulinic acid were -11.60 and -7.50 kcal/mol, respectively.
As shown in Figure 15 and Figure 16, the sesquiterpene compounds Furanopetasin and Bakkenolide IIIa formed multiple interactions with the binding site of the BMPII receptor.Furanopetasin's octahydronaphtho[2,3-b]furan formed five interactions with the essential residue Cys94, as well as interactions with Gly89, Ile88 and Ser119.In addition, each carbonyl group served as an acceptor for a hydrogen bond with the Asp90 amino acid.Finally, a Pi-Sigma bond was formed between the methyl group of the side chain and the crucial residue His87.The docking score for Bakkenolide IIIa was predicted to be favourable (−9.90 kcal/mole).Bakkenolide IIIa's furan-3,2'-inden-2one formed a hydrogen bond with Thr65, a sulfur bond with Cys94, and an alkyl bond with Pro91, whereas one of the carbonyl groups formed a hydrogen bond with Ser119.
The majority of phenolic compounds demonstrated superior docking scores compared to sesquiterpene compounds.The two highest compounds were Fukinolic acid and Petasiphenol, whose catechol groups significantly contributed to the binding interactions at the active site by forming a hydrogen bond with the His87 residue.In addition, Fukinolic acid formed five hydrogen bonds with the residues Thr65, Asp90, Pro91, and Ser119 (Figure 17).Likewise, it binds via Pi-cation to His 87 and Pi-Anion to Glu93.Four additional hydrogen bonds were observed between the hydroxyl of Petasiphenol and Gln92, Cys94, and His95, in regards to the Pi-Alkyl interaction with the essential amino acid Cys94 (Figure 18).

Metabolic study
Bakkenolide IIIa and Petasiphenol were selected for metabolic investigations due to their satisfactory ADME profiles and high docking scores.Using Biotransformer, the metabolic pathways and predicted metabolites of both compounds were estimated by (http://biotransformer.ca) (Djoumbou-Feunang et al. 2019).The chosen substances were uploaded to the server.Then, metabolic changes during phases I and II were selected.The prediction of phase I metabolites for both compounds did not show any results.The resultant data was arranged into metabolic processes that described the expected metabolites and the enzymes of phase II transformation that may perform this role in humans.The outcomes are depicted in Figure 19 and Figure 20.In the metabolic investigations of substances Bakkenolide IIIa and Petsiphenol, respectively, three likely metabolic pathways and eight predicted metabolites were identified.Figure 19 depicts the metabolic steps involving thioether Smethylation and GSH-conjugation of alpha, beta-unsaturated carbonyl in Bakkenolide IIIa.GSH-conjugation of alpha, beta-unsaturated carbonyl by Glutathione transferase accounted for the majority of the anticipated phase II metabolites of Bakkenolide IIIa.In contrast, Petasiphenol's metabolic pathways included aromatic OH-glucuronidation by UDP-glucuronosyltransferase and catechol O-methylation by catechol O-methyltransferase (Figure 19).

Conclusion
In this research work which represent a journey through natural phyto-constituents to investigate their ability to act as modulators for bone morphogenic protein, several computational tools were used to assess the pharmacokinetic profile of the selected library.Molecular docking studies reveled more promising natural compounds than the reported BMP agonists, which indicate their high promiscuity to act against bone damage, Bakkenolide IIIa showed dual potential on both BMPRIA and BMPRII with docking scores -7.82 and -9.9 Kcal/mol respectively, this is more better score than Betulinic acid agonist.
The phenolic compounds Fukinolic acid and Petasiphenol showed more better docking scores with both BMPRIA and BMPRII than Icariin agonist.All these findings support the need for further preclinical investigations for these hopefully bone regulating candidates to treat bone loss and damage.

Computational admet studies
Physicochemical characteristics and ADME were evaluated using the SwissADME web server (http://www.swissadme.ch/index.php),the Molsoft server (https://www.molsoft.com), and the PreADME server (https://preadmet.bmdrc.kr), as well as the DataWarrior programme.Group Inc., Montreal, Quebec, Canada) to assess the free energy and interaction modes of molecules with respect to BMP-2 receptors.The crystal structure was initially prepared by eliminating water molecules while preserving the crucial chain.The protein structure was then protonated, the missing hydrogens were inserted to provide the desired ionization states, and the hydrogen atoms were concealed.
Then, the amount of energy was reduced.With the default settings, the MOE Alpha Site Finder was used to identify the active site.The obtained alpha spheres were used to construct the inactive site's dummy atoms.Preparation of ligands: Using ChemBioDraw Ultra 16.0, the 2D structures of the investigated drugs were sketched and saved in MDL-SD format.After opening the saved files with MOE, 3D structures were protonated.Next, energy reduction was implemented.
Molecular docking calculations: Molecular docking of the tested drugs against the target receptor was done using the default technique.Using genetic algorithm searches, docked structures were created for each instance.The forcefield's refinement code was MMFF94X.
Selection of best pose: The poses of the compounds with the highest ratings are chosen.Visualizations: We utilized the Biovia discovery-studio 2020 visualizer to examine the representation of protein-ligand interactions in the active site of the complex.In conclusion, the H-bonds, Piinteractions, and hydrophobic are regarded as the most significant interactions discovered.

Metabolic study
(http://biotransformer.ca) estimated the metabolic pathways and predicted metabolites of both compounds using Biotransformer.Firstly, the smile of compounds Bakkenolide IIIa and Petasiphenol were obtained from chemdraw program and were used in http://biotransformer.ca to predicate phase I and II metabolism.

Figure 1 .
Figure 1.BMP2 binds to its receptors and activate the SMAD signaling pathway.

Figure 3 .
Figure3.Reported compounds that can affect bone formation.

Figure 4 .
Figure 4.The chemical structures of the sesquiterpene compounds from P. Japonicous.

Figure 4
Figure 4 Continued

Figure 5 .
Figure 5.The chemical structures of the phenolic compounds from P. Japonicous.

Figure 7 .
Figure 7. 2D and 3D representation of predicted binding mode of Icariin with BMPIA receptor.

Figure 8 .
Figure 8. 2D and 3D representation of predicted binding mode of Batulinic acid with BMPIA receptor.

Figure 9 .
Figure 9. 2D and 3D representation of predicted binding mode of Bakkenolide IIIA with BMPIA receptor.

Figure 11 .
Figure 11.2D and 3D representation of predicted binding mode of Fukinolic acid with BMPIA receptor.

Figure 12 .
Figure 12. 2D and 3D representation of predicted binding mode of Petasiphenol with BMPIA receptor.

Figure 13 .
Figure 13.2D and 3D representation of predicted binding mode of Icariin with BMPII receptor.

Figure 14 .
Figure 14.2D and 3D representation of predicted binding mode of Betulinic acid with BMPII receptor.

Figure 15 .
Figure 15.2D and 3D representation of predicted binding mode of Furanopetasin with BMPII receptor.

Figure 16 .
Figure 16.2D and 3D representation of predicted binding mode of Bakkenolide IIIa with BMPII receptor.

Figure 17 .
Figure 17.2D and 3D representation of predicted binding mode of Fukinolic acid with BMPII receptor.

Figure 18 .
Figure 18.2D and 3D representation of predicted binding mode of Petasiphenol with BMPII receptor.

Figure 19 .
Figure19.Predicted phase II metabolic pathways and metabolites of Bakkenolide IIIa in humans.

Table 1 .
Physicochemical properties of compounds that follow Lipinskís and veber rules.
a Caco2: Low values were less than 4 nm/s, medium values ranged from 4 to 70 nm/s, and high values were greater than 70 nm/s.bHIA: From 0% to 20% poor, 20% to 70%, and 70% to 100% well absorbed.JOURNAL OF PLANT INTERACTIONS

Table 4 .
The binding energy values (kcal/mol) for all compounds.