Exploring the potential of phytoconstituents from Phaseolus vulgaris L against C-X-C motif chemokine receptor 4 (CXCR4): a bioinformatic and molecular dynamic simulations approach

Introduction The CXCR4 chemokine receptor is a G protein-coupled receptor that plays a role in many physiological processes and diseases, such as cancer metastasis, HIV infection, and immune response. Because of this, it may be possible to target it therapeutically. In addition, the active ingredient of Phaseolus vulgaris L (PVL) has been reported to have anti-inflammatory, antioxidant, and anticancer properties. Novel CXCR4 antagonists from natural resources can be a promising drug development product using a computational approach. This study aims to explore the active compound in PVL that has the responsibility to inhibit CXCR4 using molecular docking and dynamics simulation. Materials and methods Pharmacokinetic analysis were performed using the pkCSM, OSIRIS for toxicity risk analysis, and the PerMM for membrane permeability assessment. Molecular docking was performed using PyRx software to determine the interaction between the CXCR4 target protein from the PDB database and the active component of PVL from the PubChem database. A molecular dynamics (MD) simulation was performed to determine the stability of the interaction using the WEBGRO Macromolecular Simulations online server. The analysis were performed by comparing the results with plerixafor as a control ligand. Results and discussion The pharmacokinetic analysis of quercetin, kaempferol, myricetin, catechin, 3,4-dihydroxy-benzoic acid, and daidzin in PVL showed that they met the drug-like


Introduction
The C-X-C motif chemokine receptor 4 (CXCR4) is a chemokine receptor that couples with G-proteins and is primarily expressed on endothelial cells and pericytes [1].The activation is caused by the homeostatic chemokines stromal cell-derived factor-1 (SDF1, CXCL12), which initiates several biological processes [2,3], including mesenchymal stem cell (MSC) mobilization to inflamed tissues [3,4].The interplay between CXCL12 and CXCR4 is complex and dependent on physiological and pathological conditions [2,3].Modulation of the CXCL12/CXCR4 axis contributes to disease progression, including cancer, autoimmune disease, and certain cardiac and neurological diseases [2].
CXCL12/CXCR4 is also involved in inflammation through several different pathways, such as PI3K/ AKT, NF-κB, and JAK/STAT [1,3,5,6].The CXCL12/CXCR4 pathway is crucial for the migration of various immune cells, including leukocyte trafficking, T and B cell migration and stimulation to their subtypes, inflammatory cytokine production [1], hematopoietic stem cell homing, survival, and maintenance in bone marrow [3].The maintenance of activated immune cells by the CXCL12/ CXCR4 pathway contributes to chronic inflammation in autoimmune inflammatory arthritis.The NF-kB activation and ERK phosphorylation pathways simulate IL-6, which subsequently promotes osteoclastogenesis [1].Additionally, CXCL12 enhances apoptotic resistance and attracts the precursors of osteoclasts, contributing to an increase in bone-resorbing activity and cartilage destruction [7].Based on the existing evidence regarding the role of CXCL12/CXCR4 in different diseases, there is a potential need to develop inhibition agents for this axis.
Several studies have revealed that the extracellular region of CXCR4 contains a significant negative potential, including crucial negatively charged residues (Asp187, Asp97, Asp262, and Glu288) that make up the essential binding region for its native ligand, CXCL12 [8,9].AMD3100 (Plerixafor), a synthetic antagonist for CXCR4, has significant potential and can strongly bind to specific residues for CXCR4 interaction.A recent NMR study showed that the CXCR4 antagonist AMD3100 could displace the CXCL12 N-terminus from the receptor without displacing the chemokine core domain [10,11].Previous studies also indicated that AMD3100 can mobilize stem cells to the peripheral blood [12].Moreover, various studies have shown that the administration of AMD3100 can stimulate stem cell mobilization, inhibit CXCL12-induced angiogenesis, and reduce inflamed vessels in rheumatoid arthritis, all of which can contribute to the regeneration of damaged bone tissue [2,[12][13][14].However, in addition to its high-cost [15], AMD3100 is not widely available in many developing countries, including Indonesia.Intravenous injection administration poses significant challenges and has numerous immediate side effects.Therefore, studying new antagonists for CXCR4 from herbal medicine is one of the strategies to reduce cost and can be widely distributed worldwide.
Phaseolus vulgaris L., also known as black beans, is a popular dietary source that is widely distributed in Indonesia and belongs to the legume family.P. vulgaris L. (PVL) beans are known for their high flavonoid content, which includes quercetin, kaempferol, myricetin, genistein, and tannic acid [16,17].Several studies have demonstrated that PVL exhibits anti-inflammatory, antioxidant, and anticancer properties, as well as lipid and glucose-lowering properties.PVL extracts in water and acetone can reduce inflammation by blocking TNF-α and CXCR4 [17,18].Moreover, PVL extract facilitates MSC homing and differentiation in the pancreatic tissue [19].This study aims to explore the active constituent from PVL that is responsible for inhibiting CXCR4 using molecular docking and dynamics simulation.

Data mining of active compounds in P. vulgaris L.
We obtained the bioactive compounds of PVL from the literature review [20].Additionally, PubChem ID, Molecular formula, and 2D structure were gathered for each compound.
The membrane permeability was assessed using PerMM (https:// permm.phar.umich.edu/ permm_ server_ cgopm).The PerMM web server and database allow quantitative investigation and visualization of the passive translocation of bioactive chemicals across lipid membranes.A groundbreaking physics-based web tool, the server calculates membrane binding energies and permeability coefficients for a wide range of compounds on the phospholipid layer, PAMPA-DS, blood-brain barrier, and Caco-2/MDCK cell membranes.This study also visualizes transmembrane translocation pathways, showing the successive translational and rotational positions of a permeant as it passes through the lipid layer.This study also illustrates the change in solvation energy resulting from this process.The server predicts permeability coefficients for compounds with many chemical scaffolds.Selection and optimization help identify and improve drug leads [23].

Protein-protein analysis of CXCR4 and CXCL12
We performed protein-protein interaction prediction between CXCR4 and CXCL12 because there is no experimental structure interaction between CXCR4 and CXCL12 in a database.Previous research by Xu et al. [10] predicted the interaction of CXCR4 and CXCL12 using ZDock.In this research, we generate an AI (artificial intelligence) prediction model using AlphaFold2 Multimer.We used CXCR4 and CXCL12 from the protein database (CXCR4: 3ODU; CXCL12: 4UAI), then input the amino acid sequences into AlphaFold2 Multimer and set them as default settings.The AlphaFold2 system will predict the structure and save it as a PDB file [24,25].

Molecular docking analysis
Molecular docking was performed using PyRx 0.95 software on a computer with the Windows 10 operating system, 8 GB of RAM, a 500 GB NVME SSD, and an AMD Athlon 3150U CPU.The three-dimensional structures of every compound were obtained through the PubChem website (https:// pubch em.ncbi.nlm.nih.gov/) and saved in a (.sdf) format file.The CCR4 as the protein target used in this study was obtained from the Protein Data Bank database (PDB) with protein ID 3ODU (https:// www.rcsb.org/) [26,27].The three-dimensional structure of the protein was downloaded in PDB format.The downloaded proteins were stabilized by removing water, ligands, and hydrogen using PyMOL 2.0 software (https:// pymol.org/2/) [28,29].The binding affinity interaction prediction was calculated using the Vina Wizard integrated with PyRx 9.5, using a specific grid box with the center coordinates and dimensions shown in Table 1.Since plerixafor is an inhibitor of CXCR4, we used plerixafor as a control in this study.The coordinates were chosen based on the location of a binding pocket and the residue interaction identified in the literature review [30,31].The protein-ligand interaction was visualized using UCSF ChimeraX (https:// www.cgl.ucsf.edu/ chime rax/) and Discovery Studio 2021 (https:// disco ver.3ds.com/ disco very-studio-visua lizer-downl oad) to discover the binding site and type of molecular interactions [32,33].

Molecular dynamic (MD) simulation
Molecular dynamics simulations are an effective technique for understanding the structure-function relationships of macromolecules.In this research, molecular dynamics simulations were performed using the WEB-GRO Macromolecular Simulations online server (https:// simlab.uams.edu/ Prote inWit hLiga nd/ index.html).This web is available free of charge to all researchers worldwide for academic purposes.Molecular dynamics simulations were performed for phytoconstituents that had a higher interaction affinity than the control.The first step for the MD simulations was to generate a topology of catechin, quercetin, kaempferol, myricetin, and hydroxybenzoic acid using the GlycoBioChemPRODRG2 server (http:// davap.c1.bioch.dundee.ac.uk/ cgi-bin/ prodrg).The molecular dynamics simulation used GROMOS96 43a1 as the force field parameter, a temperature of 310 K, the addition of 0.15 M NaCl, 5000 steps of energy minimization, and a simulation time of 20 ns [34].Parameters measured in this simulation include the root mean square distance (RMSD), the root mean square fluctuation (RMSF), and the number of hydrogen bonds.

Druglikeness and ADME prediction
The drug similarity of the PVL compounds was assessed using Lipinski parameters and other parameters listed in Fig. 1A.Pharmacokinetic properties were evaluated by ADME screening and are shown in Table 3.

Toxicity prediction
The toxicity prediction of the PVL compounds was assessed using OSIRIS software, as seen in Fig. 1B.

Membrane permeability prediction
The membrane permeability prediction of the PVL compounds was assessed using the PerMM web server, which can be seen in Fig. 1C for visualization and Fig. 1D for quantification.

Identification of CXCR4 and CXCL12 interaction and plerixafor inhibition site
To determine the binding mode interaction between CXCL12 and CXCR4, we generated an AlphaFold 2 multimer prediction structure as shown in Fig. 2a.Plerixafor binds to the same binding pocket with CLXCL12 as shown in Fig. 2B.
The binding affinity interaction between the PVL chemical constituent and CXCR4 is shown in Fig. 3A.The binding pocket interaction between the active constituents of PVL and CXCR4 was analyzed and compared with the amino acid residues formed with plerixafor.Quercetin, myricetin, kaempferol, catechin, and 3,4-dihydroxybenzoic acid can bind in the same binding pocket compared to the control in the presence of Asp187, Asp97, and Glu288 as critical amino acid residues, as shown in Fig. 3B.To know the specific interactions that occur in the interaction between the active constituent of PVL and CXCR4, we visualize the interaction using Discovery Studio.The interaction of the active constituent of PVL against CXCR4 is shown in Fig. 5A-E.

Molecular dynamics simulations
To determine the stability of the interaction, we performed a molecular dynamics simulation using the WEB-GRO web server for as long as 20 ns.We visualization of post-MD interaction as shown in Fig. 5F-J.The root mean square distance (RMSD), the root mean square     fluctuation (RMSF), and the number of hydrogen bond results can be seen in Fig. 6.

Phytochemical constituent of P. vulgaris L
The phytochemical constituents of PVL have been explored in some varieties of PVL influenced by distributed area of origin, environmental conditions, or method of extractions [17,20,35].Different studies reported the phytochemical content of PVL, such as saponin, anthocyanins, flavonols, phenolic acids, proanthocyanidins, and other bioactive compounds, as well as its evidence for anti-inflammatory, anticancer, and immunomodulatory agent [17,20,36,37] as summarized in Table 2.The compound with the highest binding affinity will be selected for molecular dynamic simulations analysis as shown in Fig. 5.

Druglikeness, and ADME prediction toxicity prediction of phytochemical constituents of P. vulgaris L
To know the prediction of pharmacokinetics for each chemical component of PVL, we analyzed the Rule of Five (RoF) score.PVL active compounds were retrieved in sdf file format from the PubChem database (http:// pubch em.ncbi.nlm.nih.gov/).Our study revealed that quercetin, kaempferol, myricetin, catechin, 3,4-dihydroxybenzoic acid, and daidzin contained in PVL met the RoF criteria.Besides, other compounds did not meet the RoF criteria due to some factors such as MW > 500, HBA > 10, and LogP > 5, as shown in Fig. 1A and Table 3. RoF is a criterion for evaluating the similarity of a compound to a drug.The requirements for RoF are that the hydrogen bond acceptor (HBA) must be less than 10, the hydrogen bond donor (HBD) must be less than 5, the molecular weight (MW) must be less than 500 g/mol, the H2O partition coefficient (LogP) must be less than 5, and the molar refractivity must be between 40 and 130.
The compounds that meet the RoF are predicted to have drug-like properties [38].Pharmacokinetic properties, including absorption, distribution, metabolism, and excretion, were analyzed among the active compounds of PVL, as shown in Table 3.In absorption, solubility and permeability are essential drug-specific physicochemical properties related to the ability of the drug across the membrane to reach the desired concentration in the systemic circulation [39].Solubility is critical in developing oral drugs, as low solubility can affect intestinal absorption through the portal circulation.Notably, all PVL compounds have low water solubility (< 0.4 mol/L), indicating poor water solubility related to lipophilicity since they have LogP less than 5 [40].Low water solubility can lead to a poor dissolution rate and potentially limit intestinal absorption, whereas high lipophilicity can enhance the intestinal absorption of a drug.Other factors like permeability and formulation play significant roles in determining the overall bioavailability of a drug.Strategies such as the use of surfactants, lipids, permeation enhancers, micronization, salt formation, nanoparticles-drug delivery system, and solid dispersions are employed to overcome issues related to poor aqueous solubility and enhance the bioavailability of lipophilic drugs [41][42][43].Thus, PVL compounds can be potentially used as oral administration drugs.However, further studies are needed to optimize the absorption capacity for future drug development.
It is noteworthy that all PVL compounds have low water solubility (< 0.4 mol/L), indicating poor water solubility related to lipophilicity.However, all compounds have high intestinal absorption in humans, indicating the percentage of drug absorption in the small intestine.Thus, further studies to optimize the absorption capacity are needed for future drug development.
The distribution property shows the distribution of drugs within different body compartments.All PVL drugs have a low blood-brain barrier (BBB) (< 0.1 log BB) and central nervous system (CNS) permeability (< − 3 log PS), indicating that they cannot penetrate the BBB.Metabolic property is another essential factor that affects a drug's pharmacokinetics, pharmacodynamics, and safety profile.Cytochrome P450 is a primary component in drug metabolism located in the liver and intestine and is either induced or inhibited by various substances [40].Most of the compounds contained in PVL do not inhibit the cytochromes, except for CYP1A2 and CYP3A4 in several substances as shown in Table 3.

Toxicity prediction of phytochemical constituents of P. vulgaris L
Toxicity analysis is an essential parameter in determining the safety of compounds.To know the safety of each compound, we predicted the possibility of toxicity using the OSIRIS software, which used some parameters such as mutagenic, tumorigenic, irritant, and reproductive effects.Mutagenic and tumorigenic are parameters to predict the effect of the compound on becoming mutagenic and causing tumors.As a result, most PVL Fig. 3 A Binding affinity from molecular docking result of CXCR4 and the ligand compounds compared to ligand control; B Catechin (Red), Kaempferol (blue), Myricetin (yellow), Quercetin (cyan) and 3,4-Dihydroxybenzoic acid (orange) interact with CXCR4 in the same binding pocket compounds are safe for toxicity analysis instead of quercetin, kaempferol, myricetin, 3,4-dihydroxybenzoic acid, and daidzin.These compounds were predicted to have medium-risk mutagenic and tumorigenic effects, except for daidzin, which has high-risk tumorigenic and reproductive effects.In addition, myricetin-3-glycoside, quercetin-3-glycoside, kaempferol-3-glycoside, catechin, and tannic acid have no risk for toxicity analysis, as shown in Fig. 1B.Nevertheless, wet lab studies are required to determine the optimal dosage, and a significant risk depends on the amount (16).

Membrane permeability prediction phytochemical constituents of P. vulgaris L
The ability of compounds to penetrate the membrane revealed that the compound with the lowest transfer energy along the translocation pathways identified by its calculated transfer energy profile is more permeable across the bilayer membrane [23].To know the membrane penetration ability of each compound, we predicted using the PerMM web server.Figure 1C is a visualization of the conformational change of the PVL active compounds as they penetrated the cell membrane.Each molecule continuously adjusted its position to match the hydrophilic and hydrophobic characteristics of the plasma membrane.The analysis showed that plerixafor had the lowest energy across the membrane.Among the compounds contained in PVL, 3,4-dihydroxybenzoic acid has the lowest energy to cross the membrane, followed by kaempferol, daidzin, catechin, tannic acid, quercetin, myricetin, kaempferol-3-glucoside, quercetin-3-glucoside, and myricetin-3-glucoside, as shown in Fig. 1D.
Hydrophobic (non-polar) compounds can easily penetrate the lipid bilayer.Hydrophilic (polar) molecules, on the other hand, typically rely on transport proteins to cross the membrane because they cannot easily penetrate the hydrophobic interior of the bilayer [23].Typically, the ligand for CXCR4 does not enter the membrane.Instead, it attaches to the extracellular region of the receptor, inducing a structural change that transmits a signal throughout the cell.

CXCR4 and CXCL12 binding mode interaction and plerixafor as inhibitor
The protein-protein interaction results were then analyzed for quality using Ramachandran plots via the PRO-CHECK web server.A good quality model was confirmed based on residues in the most favored regions (> 90%), as shown in Additional file 1: Fig. S1.The interaction appears to be similar to previous studies in that CXCL12 bound to the extracellular region of CXCR4 contains a significant negative potential, including critical negatively charged residues at residues Asp187, Asp97, Asp262, and Glu288 [10,11].
As we mentioned in the introduction plerixafor is a CXCR4 inhibitor.The visualization interaction between plerixafor and CXCR4 shows that it can bind to the extracellular region of CXCR4 (Fig. 2B) at Asp187, Asp97, and Glu288 residue, but slightly far from Asp262 residue as shown in Figs.2C.The mechanism by which SDF-1 binds to CXCR4 involves the interaction of the N-terminus domain of SDF-1 with the extracellular domain CXCR4, followed by further interactions that stabilize the complex and lead to signal transduction.Superimposed interaction shows that plerixafor has a binding mode similar to CXLC12.This binding mode curled up in the CXCR4 extracellular binding pocket of the CXCr4 [44] as shown in Fig. 2D.2D interaction shows that plerixafor can bin to CXCR4 by forming hydrogen bond at Glu288, Tyr45, Asp97, Leu41, Ile185, Arg30, His281, Ala98, Ser285, Cys186, Trp102, Val112, Asp187, Arg188; hydrophobic bond at His113 and Trp94 as shown in Fig. 2E.As we mention before that Glu288, Asp187 and Asp97 is important residue in binding interaction CXCR4/CXCL12, indicating the reasonable of plerixafor as an established drug for CXCR4 antagonist.Moreover, Plerixafor is the only CXCR4 antagonist approved by the FDA and has been commercially distributed.Several studies have developed potential candidates for synthetic CXCR4 antagonists, such as POL6326 (Balixafortide), LY2510924, TN14003, and MSX-122, which exhibit good safety and tolerability profiles yet are still in preclinical and clinical studies phases [45][46][47].So, we used plerixafor as the control ligand in this study.

Molecular docking before and after dynamics simulation of P. vulgaris L active constituent against CXCR4
A molecular docking study is a research model used in drug discovery to determine the binding interaction between ligand and protein [48].The binding affinity measures the strength of the interaction between two molecules [49].Our molecular docking analysis identified that quercetin (− 6.6 kcal/mol), myricetin (− 6.6 kcal/ mol), kaempferol (− 6.3 kcal/mol), catechin (-6.5 kcal/ mol), and 3,4-dihydroxybenzoic acid (− 5.4 kcal/mol) bind to CXCR4 with the highest affinity compared to plerixafor (− 5.0 kcal/mol) as the control ligand.The remaining compounds with the lowest affinity compared to the control ligand are myricetin-3-glucoside (− 2.5 kcal/mol), quercetin-3-glucoside (− 0.6 kcal/ mol), kaempferol-3-glucoside (− 2.5 kcal/mol), daidzin (− 2.8 kcal/mol), and tannic acid (− 4.2 kcal/mol), with lower affinity compared to plerixafor.The molecule with the lowest binding energy will have a constant temperature and pressure, called a stable molecule.The amino acid residues influenced the binding domain of the target protein and the type of chemical interactions in the binding region [50].A lower energy value of binding affinity denotes a more stable and favorable binding relationship between the target protein and the ligand [49].However, further experimental studies are necessary to confirm the actual protein-ligand interaction [51].
In addition, based on MD simulations, kaempferol-CXCR4 interacts by forming a hydrogen bond with Phe13, Glu15, Ser16, Trp57, Leu55, and Val18; and a hydrophobic bond with Phe14 (Fig. 5I), while quercetin post-MD simulation formed a hydrogen bond at Ser28, Leu41, Tyr45, Val112, Tyr116, and Arg188; a hydrophobic bond at Trp94 and Ala98; and a pi-anion bond with Asp97 and His113 (Fig. 5J).Post-MD interaction analysis showed that quercetin, myricetin, and 3,4-dihydroxybenzoic acid still can bind to critical amino acid residue at Asp97, Asp187, and Glu288 as shown in Table 4 and Fig. 5.The stability of the CXCR4 protein refers to the collective pressures determining whether the protein will maintain its folded shape or adopt non-native aggregating configurations.Comparing the structure of protein-ligand complexes at different levels of stimulation up to 20 ns provides valuable structural insights that help to understand the potential changes in ligand position that occur as a result.
The protein-ligand interactions play a vital role in structural biology by providing insight into the mechanisms of these interactions at the molecular level.The similarity of the binding to the control indicates the same function.A similar interaction active constituent of PVL can be seen in Table 4. Identification of molecular interactions and binding orientations on the docked protein-ligand complex revealed that PVL compounds form non-covalent interactions with all target proteins through hydrophobic, pi, and hydrogen bonds.These interactions lead to the development of the proteinligand complex and trigger the initiation of an activity response, including enhancement and inhibition of the target protein [22,48].Together with another additional factor, such as binding affinity (binding affinity analysis), the similarity of those interaction types between a natural compound and the control ligand (interaction analysis) could indicate the accuracy of the docking method and support the potential of a ligand as a drug candidate [52,53].The critical molecular interaction, hydrophobic, pi, and hydrogen bonds can help to stabilize the ligand within the protein's active site, leading to a stronger binding affinity and potentially inhibitory activity [53,54].The similarity interaction between the PVL compound and the control ligand (Plerixafor) suggests the potential inhibitory activity as a CXCR4 antagonist candidate.Since molecular docking solely is insufficient to conclude the promising activity, further methods such as molecular dynamic simulation will provide information about the stability and dynamics of the protein-ligand complex by optimizing the structures of the final complexes from docking, calculating detailed interaction energies, and providing information about the ligand binding mechanism [55], as we conducted in this study.5 Amino acid residues result from the interaction between ligands and CXCR4.Panels A to E show interaction before molecular dynamics, while panels F to J show residues after molecular dynamics simulations indicating that the interaction of myricetin with CXCR4 is unstable.3,4-dihydroxybenzoic acid fluctuates from 0 to 6 ns and is stable from 6 to 20 with an RMSD score of 0.7 nm.Kaemferol was stable from 3 to 12 ns with an RMSD score of 0.4 nm, then fluctuated to 16 ns and stabilized at 20 ns.Interestingly, quercetin has a similar stability interaction with an average RMSD score of 0.3 nm from 0 to 20 ns, indicating that quercetin has a similar stability interaction compared to plerixafor, as shown in Fig. 6A.

The molecular dynamics simulation
The PVL compound causes significant fluctuations in amino acid residues in certain regions of the protein, as shown by the RMSF curve in Fig. 6B.Analysis of residuespecific root mean square fluctuation (RMSF) provides insight into how ligand binding affects the structural flexibility of a protein at the single amino acid level [57].A higher RMSF value indicates greater flexibility of the complex.Leu68, Asp97, Leu146, and Thr178 have the largest backbone variation of 5 nm.The significant variations in amino acid residues suggest the fluctuating interactions between the compound and these residues, possibly involving hydrogen bonding, ionic interactions, or van der Waals forces.The PVL compounds appear to have a remarkable effect on the structure and dynamics of the protein, potentially serving as a modulator or regulator of important protein functions.
To explain the conformational stability of the interaction, we examine the total number of intermolecular hydrogen bonds in the ligand-protein complexes.Hydrogen bonds play an important role in the stability of the protein structure [55].Interestingly, the number of hydrogen bonds in catechin, hydroxybenzoic acid, myricetin, and quercetin compounds has the same value as plerixafor, with an average of 380.but the number of hydrogen bonds in kaempferol compounds has a lower average of 320, as shown in Fig. 6C.

Future potential active constituent of PVL as CXCR4 inhibitor
PVL is rich in a variety of phytochemicals, including flavonoids, phenolic acids, and other antioxidants [17].Some of these compounds have shown biological activity in various contexts, although their specific interaction with CXCR4 would require targeted research.

Compounds Molecular Interaction Molecular interaction after MD
The CXCR4/CXCL12 pathway has been proposed as a target for stem cell mobilization-based therapy in various rheumatic diseases.Several studies have reported that upregulation of CXCR4 and CXCL12 is associated with joint erosion, synovial inflammation, synovial hyperplasia, and synovial angiogenesis [58].Mesenchymal stem cells have been shown to stimulate chondrocyte regeneration and differentiation into cartilage [59].Osteoarthritis, stem cell-based treatment has shown promising clinical effects, including improved joint function, pain threshold, and quality of life [60].In comparison, stem cell-based treatment can modify and restore the balance of inflammatory T cells in rheumatoid arthritis [61].
CXCL12/CXCR4 axis also has a potential therapeutic target for several inflammatory diseases, not only by affecting cell migration but also by altering the immune response.Only one antagonist targeting the CXCR4 ligand binding region, plerixafor, has shown therapeutic relevance (1).In addition, the role of chemokines in immune modulation in autoimmune diseases remains to be explored.Chronic diseases, especially autoimmune diseases such as rheumatoid arthritis (RA), inflammatory bowel disease (IBD), and systemic lupus erythematosus (SLE), are associated with an abnormal inflammatory response because the immune system recognizes the protein as an antigen and attacks itself [62][63][64].
This study shows that several candidate active constituents of PVL can potentially become CXCR4 inhibitors.Computational tools and molecular modeling can predict the interaction between potential phytochemicals in PVL against the CXCR4 receptor.Such in silico studies can be a cost-effective first step before experimental studies, but this finding needs to be clarified with more precise and advanced studies [65].Any potential CXCR4 inhibitors identified would need to undergo rigorous preclinical and clinical testing to establish their safety and efficacy.
The efficacy of natural product derived-CXCR4 inhibitor flavonoids, isoflavones, bioketones, and isoprenoidyl has been observed through CXCR12/ CXCR4 axis.Promising inhibitory activity of flavonoid compounds such as Quercetin [66,67] and myricetin has been reported in downregulating CXCL12 and CXCR4 expression in prostate cancer [68].However, further safety assessment is necessary to validate those efficacious activities of the novel CXCR4 antagonist.Research in CXCR4 antagonist development from natural products is escalating despite the compound's pharmacokinetic, pharmacodynamic, and toxicity testing not being comprehensively studied, which is insufficient to replace the established agent Plerixafor.
We understand that there are limitations to this study.In our study, we only predicted the potential of the active component of PVL against CXCR4 using computational study.Therefore, we cannot confirm the effect of PVL on CXCR4 in experimental studies using cells or animals.At least, this research can support and serve as a basis for further research in vitro and in vivo.Further research is needed to clarify our predictive findings through experimental studies, such as examination of binding interaction and visualization using X-rays or cryo-EM.Experimental toxicity studies are also needed to confirm PVL as a candidate for CXCR4 inhibitor.It is important to note that while natural products are a promising source for drug discovery, the path from identifying a potential lead compound to developing a clinically approved drug is long, complex, and challenging.Any findings would need to be substantiated through rigorous scientific research and clinical trials.At this time, any potential CXCR4 inhibitory compounds in PVL remain speculative and would require significant research to validate.

Conclusion
In sum, virtual screenings of active constituents of PVL using molecular docking and dynamic simulation revealed that quercetin, myricetin, and 3,4-dihydroxybenzoic acid have the potential to become CXCR4 agonists.These three active constituents have good oral bioavailability and toxicity.Further research is required to fully understand and clarify PVL as a CXCR4 inhibitor.

N 8 -Fig. 1 A
Fig.1A Lipinski RoF criteria analysis showed by yellow and black code indicating fulfilled or not fulfilled the criteria; > 1 parameters fail indicated not pass the criteria.B Toxicity analysis of the compounds against mutagenic, tumorigenic, irritant, and reproductive effects parameter categorized by 0 (high risk) (magenta), 0.5 (moderate risk) (light green), and 1 (safe) (red); C Anti-inflammatory probability activity analysis of all compounds; D) Permeability analysis showing the transfer energy for the compounds across the bilayer membrane

Fig. 2 A
Fig.2A Interaction of CXCR4 (green) and CXCR12 (blue) as a native ligand; B Interaction of CXCR4 (green) and Plerixafor (magenta) as a control ligand; C CXCR4-Plerixafor interaction through three critical residues (red))(Asp97, Asp187, Asp262, and Glu288); D Superimposed interactions between CXCR4, CXCL12, and Plerixafor; E 2D visualization interaction of plerixafor against CXCR4 using Discovery studio 2D plots were generated to illustrate the varying behavior of the docked complex at different time intervals during the MD simulation runs.The plots play a critical role in the statistical analysis of the MD simulations.They provide valuable insight into the stability and flexibility of the residues at different time points during the simulations[56].The RMSD results show a stable interaction between plerixafor and CXCR4.The value was approximately 0.3 nm from 0 to 20 ns.The active constituent of PVL, catechin RMSD, shows stability from 2 to 11 ns with an average RMSD value of 1.1 nm, then unstable fluctuation up to 20 ns.Myricetin shows instability from 0 to 20 ns with an average RMSD value of 1 nm,

Fig. 6
Fig. 6 The molecular dynamic of PVL against CXCR4.A Root Mean Square (RMSD), B Root Mean Square Fluctuation (RMSF), C Number of Hydrogen bond

Table 1
Grid box docking dimension of study

Table 2
List of PVL compounds Comp.

Table 3
Pharmacokinetic of PVL prediction

Table 4
Molecular interaction active compound of P. vulgaris L. with CXCR4