In silico identification of compounds from Piper sarmentosum Roxb leaf fractionated extract inhibit interleukin-6 to prevent rheumatoid arthritis

Objective: Medicinal herbs with a phytonutrient background has been applied globally as major alternatives to ameliorate the continuous increase in rheumatoid arthritis cases worldwide. We herein aimed to critically examine the bioactive components of the medicinal herb Piper sarmentosum Roxb leaf fractionated extract for its potential to inhibit the influx of interleukin-6 (IL-6) in rheumatoid arthritis. Methods: The Schrödinger platform was employed as the main computational acumen for the screening of bioactive compounds identified and reference compounds subjected to molecular simulation (MDS) for analyzing the stability of docked complexes to assess fluctuations and conformational changes during protein–ligand interactions. Results: The values of the simulatory properties and principal component analysis (PCA) revealed the good stability of these phytochemicals in the active pocket of interleukin-6 (IL-6). Discussion: Our findings reveal new strategies in which these phytochemicals are potential inhibitory agents that can be modified and further evaluated to develop more effective agents for the management of rheumatoid arthritis, thereby providing a better understanding and useful model for the reproduction and/or discovery of new drugs for the management of rheumatoid arthritis and its complications.


Introduction
Rheumatoid arthritis (RA) is known as inflammatory arthropathy and is characterized by the presence of long-standing chronic autoimmune inflammatory disease primarily affecting synovial joints, leading to cartilage damage and bone erosion-a disabling and painful condition that can result in significant loss of function and mobility if left untreated (Scott et al., 2010;Amaya-Amaya et al., 2012).The protracted effect of this inflammatory disorder has been reported to ameliorate systemic challenges in several tissues and organs, thereby causing serious effects that eventually result in serious comorbidities such as cardiovascular diseases, osteoporosis, fractures, cognitive dysfunction, and depression (Amaya-Amaya et al., 2013).Research has shown the predominance and disease status of RA on account of geographic region to region, with higher estimates in industrialized countries and urban areas generally being greater.RA can manifest at any age, but its prevalence is greater among middle-aged individuals, particularly females (Finckh et al., 2022).
Furthermore, genetic factors have been identified as one of the major contributors to the pathogenic effect of RA, with higher heritability estimates for seropositive rheumatoid arthritis patients than for seronegative patients (McInnes and Schett, 2011;Aletaha and Smolen, 2018).However, given substantial recent advancements in understanding RA pathogenesis, this condition is incurable.Therefore, early evaluation, diagnosis, and management are paramount for achieving optimal therapeutic outcomes.Interleukin-6 (IL-6) is a pivotal immunomodulatory cytokine that plays a crucial role in the pathogenesis of various diseases, including autoimmune disorders, chronic inflammatory conditions, and cancer (Hirano, 2021).In the context of RA pathogenesis, IL-6 plays a central role.IL-6 interacts with various cell types, including B cells, T cells, osteoclasts, osteoblasts, osteocytes, and synovial fibroblasts, thereby triggering the migration of neutrophils into joints and the transition from acute to chronic inflammation (Srirangan and Choy, 2010;Takeuchi et al., 2021).
Moreover, in the B and T-cell differentiation and activation processes, IL-6 plays an important role by influencing the process, thereby resulting in heightened induction of the adaptive immune system, sustaining chronic inflammation and joint degradation, and significantly contributing to RA-related pathogenesis.These mechanisms also underlie the systemic manifestations of RA, such as pain, fatigue, morning stiffness, anemia, and weight loss, ultimately culminating in bone and cartilage erosion.This change is accompanied by synovial membrane inflammation and thickening, resulting in irreversible joint damage (Takeuchi et al., 2021;Alivernini et al., 2022).Nonetheless, IL-6R blockade has emerged as a crucial target for cytokine blockade in RA treatment, and effective management strategies for treating this disease significantly mitigate the impact of IL-6.Therapeutic approaches for treating RA aim to impede disease progression and prevent joint deformities, and IL-6R blockade represents a targeted treatment strategy.The benefits of IL-6R blockade include improvements in many extra-articular manifestations, including pain, fatigue, and anemia.Furthermore, it may have potentially positive effects on certain comorbidities (Favalli, 2020).Among the innovative treatments for RA and its associated manifestations, within and outside the joint, tocilizumab and sarilumab are biologic drugs approved as humanized IL-6 receptor inhibitors, as recommended in the 2021 treatment guidelines from the American College of Rheumatology (Fraenkel et al., 2021).
Furthermore, herbal medicine research offers promising adjunctive therapies to improve the quality of life of individuals with RA.These include the significant effectiveness of herbal remedies in alleviating RA symptoms, reducing disease manifestations, lowering disease activity scores, minimizing serum acute-phase reactants and RA biomarkers, and decreasing the likelihood of adverse events associated with disease-modifying antirheumatic drugs (DMARDs) (Soeken et al., 2003;Ernst and Posadzki, 2011;Daily et al., 2016;Aryaeian et al., 2019;Li and Zhang, 2020;Zhang et al., 2021;Zeng et al., 2022).A nationwide population-based study conducted in Taiwan revealed a high prevalence of RA patients opting for traditional Chinese medicine (TCM) herbal treatments as complementary therapies.This preference was particularly prominent among patients dealing with anxiety, depression, allergic rhinitis, osteoporosis, menstrual disorders, and menopausal syndrome (Huang et al., 2015).
Considering both the plant's safety profile and its longstanding reputation as a PS for managing inflammation-associated arthralgia and as a potential osteogenic agent.Flavonoids are recognized as key components in traditional herbal medicine for RA treatment, and our study aimed to assess potential flavonoids present in PS leaves through methanol extraction from previous studies (Miean and Mohamed, 2001;Vimala et al., 2003;Ugusman et al., 2012;Purba et al., 2021;Ware et al., 2023).This assessment was performed using computational methodologies, including molecular docking and dynamic simulations.Molecular docking and simulation are important for understanding the interactions between target protein molecules and prospective drug candidates, providing valuable insights into the mechanism of drug activity at the binding site of The ADEMET predictive study of (A) WSP water soluble properties, (B) pharmacokinetics, lipophilicity, (C) drug likeness and medicinal chemistry of the top-scoring compounds from PS leaf fractionated extract.target proteins (Ekins et al., 2007).To the best of our knowledge, there has been no prior research on the molecular docking and dynamic simulation of bioactive compounds from PS.Therefore, our study endeavors to elucidate the molecular-level interaction mechanisms involved and identify potential flavonoid compounds as IL-6 inhibitors of PS.This research contributes to our understanding of novel drug opportunities for impeding the local and systemic progression of rheumatoid arthritis.
2 Materials and methods

Protein preparation
The Protein Data Bank (PDB:1P9M) repository was used to retrieve the three-dimensional (3D) X-ray crystal structure of the interleukin-6 (IL-6) receptor, followed by removal of the cocrystallized ligands.Next, Glide's protein preparation wizard panel (Schrödinger Suite 2022-3) was notified for assigning bond ordering, adding hydrogen, or forming disulfide bonds, while Prime was used to fill in missing side chains and loops.OPLS2005 was used to reduce the structure, after which the structure was optimized using PROPKA after water molecules beyond 3.0 of the heteroatoms were eliminated (Ekins et al., 2007;Umar et al., 2021).Finally, the receptor grid file was created to define the ligand-binding pocket.

Ligand preparation
The LigPrep module from the Schrödinger Suite 2022-3 was used to prepare a total of fifty-five (55) phytocompounds of PS leaf fractionated extract identified from ethnobotanical databases for molecular docking.The 3D structures were created at low energy with the proper chiralities.At a physiological pH of 7.2 ± 0.2, possible ionization states for every ligand structure were generated.The stereo-isomerities of each ligand were calculated by keeping certain chiralities constant while varying each other (Ojo et al., 2021).

Receptor grid generation
For ligand docking, receptor grid generation allows determination of the position and size of the protein's active region.Using the receptor grid construction tool in Schrödinger Maestro 12.5, the scoring grid was supported by the crystal structure of interleukin-6 (IL-6).The van der Waals (vdW) radius scaling factor of the nonpolar receptor atoms was set to 1.0, with a partial charge cutoff of 0.25.Consequently, the top-scoring compound from the molecular docking analysis was further subjected to induced fit docking using Maestro 12.5's induced fit docking panel.

Protein-ligand docking
The outcome of the receptor grid file was utilized to perform molecular docking investigations with the Glide tool of Schrödinger Maestro 12.5.Standard precision (SP) was employed to dock the protein interleukin-6 (IL-6), and therefore, ligands were prepared from P. sarmentosum Roxb along with the quality activator, keeping the ligand sampling set to flexible and ligand sampling set to none (refine only).For ligand atoms, the vdW radius scaling factor was scaled at 0.80 with a partial charge cutoff of 0.15.

Multiple-ligand and receptor-ligand complex modeling
The chemistry of the ligands was properly standardized and extrapolated, and the structure files (SDF) of the test compounds sourced from the PubChem database were prepared using the LigPrep panel of the Schrödinger suite (Schrödinger 2022-3, Comprehensive 2D and 3D molecular interaction data for 5281672 with the receptor (interleukin-6 (IL-6)) obtained by docking analysis.
LLC, New York, NY, United States) and were used for modeling via PHASE.PHASE automatically aligns the ligands and supports the optimum arrangement and mutual characteristics of the ligands.PHASE was employed to create a receptor-ligand complex model employing the highest nine compounds with the best binding affinity against the target protein regarding the quality.The hypothesis was set with a maximum number of features to be created as 7.00, a minimum feature-feature distance of 2.00, and a minimum feature-feature distance of 4.00 for components of the identical type and donors as vectors.Aside from the donor and negative ionic features, which were set to 1, the hypothesis difference criteria were retained, and the acceptor and negative features were made equivalent.

Binding free energy calculation
The Prime Molecular Mechanics-Generalized Born area MM-GBSA tool (Schrödinger suite 2022-3) was used to determine the stability of protein-ligand complexes according to their binding free energy.The ligands were prepared beforehand using LigPrep, and therefore, the relevant proteins were prepared using the protein preparation wizard, as described previously.Sitemap anticipated the active sites of the proteins.Glide standard precision (SP) docking was then performed to dock the chemicals with proteins.The MM-GBSA technology offered with Prime was utilized to determine the binding free energy for ligand-protein complexes utilizing the Prime MM-GBSA panel.The OPLS3 physical phenomenon was chosen; therefore, the continuum solvent model was VSGB.The default settings for the opposite options were selected (Ojo et al., 2021).The concentration of each of the ligands containing ΔGbind with the LASV nucleoprotein was calculated via the following equation: ΔGbind = ΔE + ΔGsolv + ΔGSA (1) ΔE = E complex -E protein -E ligand, where E complex , Eprotein and E ligand are the minimized energies of the protein-inhibitor complex, protein, and inhibitor, respectively.Furthermore, ΔGsolv = Gsolv ( complex ) -Gsolv (protein) -Gsolv ( ligand ), where Gsolv ( complex ), Gsolv ( protein ), and Gsolv ( ligand ) contain the solvation free energies of the complex, protein, and inhibitor, respectively.ΔGSA = GSA ( complex ) -GSA ( protein ) -GSA ( ligand ), where GSA ( complex ), GSA ( protein ), and GSA ( ligand ) serve as the surface area energies for the complex, protein, and inhibitor, respectively (Ojo et al., 2021).

Pharmacological parameters
The absorption, distribution, metabolism, excretion, and toxicity (ADMET) of the test drugs were assessed using in silico integrative model predictions from the SwissADME and PROTOX-II servers.The construction of the MD preparation and trajectory analysis methods followed previous similar methods (Adekiya et al., 2022).The docked complexes were subjected to molecular simulation following the Desmond module of the Schrödinger software with an OPLS 2005 force field.The protein-ligand complex was bounded with a predefined transferable intermolecular potential with a three-point water model in an orthorhombic box.The neutralization of the overall charge was minimized by the addition of sodium and chloride ions to mimic physiological conditions.Both the temperature and pressure were kept constant at 310 °C and 1.01325 bar, respectively, by using a Nose-Hoover thermostat and a Martyna-Tobias-Klein barostat made from the United States.The simulation relaxation was undertaken by using an NPT ensemble after considering the number of atoms, the pressure, and the timescale.During the MD simulation, the long-range electrostatic interactions were calculated by using the particle mesh Ewald method.Furthermore, MD simulation analysis was carried out for 100 ns, and trajectory sampling was set at an interval of 100 ps with 1000 frames.The simulation outputs were analyzed and visualized by a simulation interaction diagram and an MS-MD trajectory analysis.The MD analysis was performed in replicate to avoid variation.The data was plotted by using OriginPro version 9.

Results
In molecular docking analysis, the significance of the structure of each molecule identified from a medicinal plant is crucial when compared to established molecules due to their molecular interactions, which may be influenced by the solute state of the water molecule in the aqueous phase.Therefore, ADEMET prediction studies serve as vital bioinformatics tools for in silico analysis, where toxicity predictions for the chosen compounds indicate their relative safety, as illustrated in Figure 1.Moreover, the potential of all the identified agents as small molecules when acknowledging their druggable potency can be easily observed in the active transport system of cells; however, the standards for drug likeness do not apply to natural bioactive substances.Therefore, the identified compounds exhibited therapeutic potential via protein binding and the rift domain in the State one protein complex, which revealed a small number of ligands with outstanding results, as shown in Table 1, with no breaches of Lipinski's rule of five in the categories of compounds identified, as shown in Figures 1A, B. The in silico analysis of the binding free energy potential revealed that the Comprehensive 2D and 3D molecular interaction packages of 5280443 with the receptor (Interleukin-6 (IL-6)) from docking analysis.
predicted molecules had greater free binding energies than did the reference molecules, as shown in Table 1.Amazingly, based on the results obtained with the Qikprop Schrödinger suite version from Figure 1C, the ligands of interest did not inhibit Pgp, suggesting that they favorably compete with reference small molecules predicted to be inhibitors of IL-6 in Pgp.Moreover, in comparison to reference small molecules, the ligands derived from the fractionated extract of PS leaves demonstrated gastrointestinal absorption.Consequently, a majority of these ligands function as ATP-dependent drug efflux pumps for ADME compounds with significant potency due to their substrate specificity, as indicated in Table 2.
Furthermore, our in silico analysis of drug likeness and ADMET properties revealed that all the compounds (except vitexin, myricetin, and rutin) identified from the fractionated extract of P. sarmentosum Roxb leaf are drug-like candidates based on the Lipinski rules of drug-like small molecules.In addition, the synthetic accessibility of any drug-like compound is very important in many areas of drug discovery.Moreover, historical knowledge of synthetic compounds analyses the information obtained from millions of previously synthesized chemicals and considers the complexity of the molecules.Therefore, the synthetic products obtained from our Insilco platform compete favorably with the reference molecule, as shown in Figure 1C.
Furthermore, one of the compounds identified, 5-hydroxy-7,4dimethoxyflavone, has the potential to compete with reference molecules with high gastrointestinal absorption and the ability to cross the blood-brain barrier (BBB) (Table 2).

Discussion
Despite the surge of several diseases worldwide, the use of bioactive molecules from medicinal plants for the treatment and management of diseases is still important and has received increased attention in recent decades (Udenze et al., 2014;Onikanni et al., 2022).PS plants, including leaf parts, are known for their medicinal value in the treatment and management of human diseases, apart from their traditional famous role in Chinese medicine, as indicated in the Supplementary Appendix of the Chinese pharmacopoeia because of their medicinal importance and edible value (Sun et al., 2020;Xu et al., 2021).In addition to its nutritional value, research has revealed its higher level of consumption as food due to its high proportion of proteins, minerals and fatty acids, which help yin and regulate endocrine functions (Yeoh and Wong, 1993).Moreover, the rise in the prevalence of various debilitating illnesses globally necessitates further investigation.Cardiovascular disease and inflammatory autoimmune conditions, such as rheumatoid arthritis, are prime examples.Rheumatoid arthritis is a significant condition that triggers joint inflammation, leading to exhaustion, anemia, and osteoporosis in the body (Dayer and Choy, 2010).
Moreover, as one of the major important receptors, IL-6 plays a role as a complex whole by combining with the host defense against any attack agent of infection and damage to tissues via the ability to induce acute phase reactions and blood cell production, which also has devastating effects due to its uncontrolled increase in IL-6 Comprehensive 2D and 3D molecular interaction packages of 439246 with the receptor (interleukin-6 (IL-6)) from docking analysis.
production, thereby resulting in different immune-mediated developmental inflammation and diseases (McInnes and Schett, 2007;Yoshida and Tanaka, 2014).However, additional information needs to be obtained on the anti-rheumatoid arthritis action of individual phytocompounds of PS leaf for further development of novel drugs for the treatment and management of rheumatoid arthritis.Utilizing a target-based computational approach, ligands extracted from PS leaves were screened for potential antirheumatoid arthritis activity by modeling the interaction between a receptor or protein target and a drug candidate (Sim et al., 2009;Sun et al., 2020).This approach is also useful for assessing the binding affinities of ligands for proteins, thereby enabling the depiction of the behavior of a drug candidate within the binding cavity of a receptor and hence providing insight into the biological function of the drug.Interestingly, our molecular docking of the identified compounds from the fractionated extract of P. sarmentosum Roxb leaf revealed therapeutic potential via protein binding and the rift domain in the State one protein complex, with a small number of ligands showing outstanding results, as shown in Figure 1.
Furthermore, the BBB is a limiting factor in the therapeutic effects of most medications for several disorders (Onikanni et al., 2022).Interestingly, our druggable candidate exhibited high bioavailability (0.85:85%) and was permeable to the BBB, suggesting its suitability for the treatment of neuronal impairment.Therefore, our study provided preclinical evidence that the fractionated extract of P. sarmentosum Roxb leaf attenuates rheumatoid arthritis complications via the inhibition of IL-6 influx.In addition, based on the molecular docking analysis from the in silico study, the binding energies of the five top-scoring phytocompounds were compared favorably with that of the reference ligand given that one of the promising compounds (5hydroxy-7,4-dimethoxyflavone) is a potent inhibitor of IL-6 with a binding potential energy of −9.46 compared with the reference ligand (−8.77), as shown in Table 1.Hence, this finding provided a better indication that this fractionated extract of P. sarmentosum Roxb leaf compound could inhibit the activity of IL-6 as a reference ligand.
Moreover, historical information on synthetic compounds has been obtained from millions of already synthesized chemicals, and the complexity of the molecules has been considered.Therefore, the synthetic products obtained from our in silico compete favorably with the reference molecule, as shown in Figure 1.Therefore, one of the compounds identified from the fractionated extract of P. sarmentosum Roxb leaf (5-hydroxy-7,4-dimethoxy flavone) strongly competes with the reference molecule because of its high gastrointestinal absorption and ability to cross the blood-brain barrier (BBB), as shown in Table 2.According to previous studies, the polypeptide chain is composed of 212 amino acid residues that are folded into different domains of unequal size with a 28-amino acid signal peptide, and its gene maps to chromosome 7p21.The Molecular Mechanics-Generalized Born Comprehensive 2D and 3D molecular interaction packages of 5280863 with the receptor (interleukin-6 (IL-6)) from the docking analysis.Comprehensive 2D and 3D molecular interaction packages of standard with the receptor (Interleukin-6 (IL-6)) from the docking analysis.Plot of free binding energy against docking score from in silico analysis of the ligands identified from the fractionated Piper sarmentosum Roxb leaf extract.
Expanse (MM-GBSA) was also used to validate the results of the molecular docking analysis.Furthermore, this analysis revealed how the binding affinity score can predict the specificity and binding strength of a ligand in the binding pocket of a receptor (Coleman, 2011;Genheden and Ryde, 2015;Zainab et al., 2020).Therefore, the results of our docking analysis of the complex hexameric human IL-6 of the alpha receptor with 5-hydroxy-7,4-dimethoxy flavone strongly correlated with the prime molecular mechanicsgeneralized band area −38.74.Compared with the reference ligand of −37.03, which produced unique hydrogen bonds, ASN60 was found to reside against the reference ligands with salt bridges containing ARG168 and ASN66 of the identified compound against the reference ligand of the ARG179 amino acid residue, as shown in Table 1.Furthermore, the structural interaction revealed in Figures 3-10 observed similar event with Apigenin when docked with hexameric human IL-6 of alpha receptor having displayed -Prime Molecular Mechanics-Generalized Born area −10.66 and Molecular Mechanics-Generalized Born area −42.62, respectively shown in Table 1 with hydrogen bonds of LEU68 and LYS66 compared with against the reference ligand of ARG179 amino acid residue.However, compared with those of the reference ligands, the  intermolecular interactions with the residues revealed strong interactions, as shown in Table 1, and the structural interactions are displayed in Figures 2-9.In addition, the in silico investigation also unveiled the covalent priming angle of each ligand by plotting free binding energy against the docking score from the in silico analysis of the ligands identified from the PS leaf fractionated extract.This highlighted the potential binding pocket of the promising compound in comparison to the reference ligand, as illustrated in Figure 10, showcasing the interresidue contacts within the binding pocket.
Furthermore, an important aspect of drug discovery requires molecular simulation, which could play a vital role in any Frontiers in Pharmacology frontiersin.org13 experiment (Durrant and McCammon, 2011;Borhani and Shaw, 2012).Many key drug discovery targets have been revealed by structural biology, allowing us to exploit the power of structurebased drug design, which requires the consideration of the dynamic properties of proteins.Moreover, one of the valuable steps in deciphering the functional mechanism of proteins and other biomolecules in revealing the structural basis of several diseases is molecular simulation, which is crucial in the design and optimization of proteins, small molecules and peptides.Therefore, our in silico study captured potential immune protein inhibitors in full atomic detail and at very clear and better temporal resolution.The mean square deviation (RMSD) of the backbone of the target was used to evaluate the stability of the system, thereby measuring the amount of protein change with respect to the initial structure during the simulation period, i.e., the structural distance between the coordinates.The results of the RSMF analysis of the IL-6 targets were carefully studied using 400 ns trajectories, as shown in Figure 11.The importance of protein structure prediction in a competitive assessment based RSMD is a critical point for revealing how the submitted structure fits into the target structure, thereby indicating that a lower RMSD is a better model than is the case for the target structure.The RMSD of the complex in the major system did not deviate drastically over the simulation testing period studied.The RSMD of the system observed mostly in the present study was within the range of 0.3-1.4Å, indicating that the conformation changed with respect to IL-6.Furthermore, the simulation converged and showed some stability (approximately 50-100 ns) in most of the systems, which revealed some system equilibration and simulatory time that was sufficient for rigorous analysis.The results are also presented as the mean ± SEM of the deviation plot, as shown in Table 3.
Furthermore, the dynamic behavior of the protein residues was studied by examining the RSMF pattern, as shown in Figure 11.The fluctuations of the complexed molecules are relatively similar and were observed within 1.5-8.5 Å within the stabilized system up to 100 ns of the simulation.The analysis of the mean ± SEM of the plots in the RFs of receptor 5280863 has the highest RSMF value (4.14 ± 1.82), and that of the receptor 5280443 has the highest RMSF value (3.21 ± 1.15), as shown in Table 3.
Importantly, MolSA fluctuations are mostly associated with the rearrangement of amino acid residues from either the accessible region or the buried region, thereby causing a change in protein-ligand structural modeling.Therefore, our study revealed that the reorganization of amino acid residues from either the available point or the buried point from the ligands to the receptor of interest was relatively similar throughout the simulation period.Both the dynamic behavior of the protein residues and pattern and the differences in the mean deviation were 239.75 ± 1.12, 245.15 ± 1.28, and 253.75 ± 1.37, respectively, as shown in Table 3.Furthermore, the pressure swing adsorption (PSA) analysis of the mean ± SEM of the plots showed similar results for most of the ligands, with deviation plots and a mean ± SEM for the receptor binding complex, as shown in Table 3 and Figure 12, with no noticeable decrease in pattern from the original pattern position until the end of the simulation period.With the strong interactions between IL-6 and its atoms connected with amino acid residues in the binding pockets of said receptor, changes in conformation could be induced, thereby affecting the binding interactions that could lead to receptor inhibition.The stability of the SASA values remains constant when the receptor of interest is folded into a particular conformation.However, the unfolding of this receptor could eventually lead to a larger surface area being exposed to the solvent, thereby indicating that the SASA was folded or unfolded at the position of the receptor of interest.Moreover, the crystal structure of hexameric human IL-6 receptor which examined the interactions between the protein surface and surrounding solvent molecules revealed the event in solvent-accessibility surface area from 50 ns until the end of the simulation period with equilibrium till the end of the simulation pattern as shown in Figure 13.In addition, the compactness of the receptor, perturbations, and folding status in gyration (rGyr) values reveal unfolding events in the structure, revealing higher and lower (rGyr) representations.The radius of gyration was analyzed to determine changes in compactness.The mean plot shows a strongly stabilized system with the complex from 10 ns of the simulation period time throughout the 100 ns simulation period, as shown in Figure 13, with the simulation converging and showing some stability at approximately 240-310 ns in the system.The deviation plots and mean ± SEM of the receptor binding complex are shown in Table 3.Moreover, the primary component, which examines the covariance matrix of the positional variances of the simulated receptors, is a key aspect studied to enhance the understanding of the structural stability of the residues towards a group of linearly related and unrelated variables.Consequently, as depicted in Figure 14, the scatter plot generated for the receptors and ligands clearly demonstrates significant differences in the ligands.Furthermore, through the diagonalization of the covariance matrix and utilization of the Origin-pro interface, as shown in Figure 15, the eigenvalues are revealed, providing essential information about the interconnected movements among the receptors.In summary, our results indicate that the ligands isolated from the fractionated extract of P. sarmentosum Roxb leaf display minor deviations in pattern and fluctuations compared to those observed in the reference.Furthermore, the dynamics of biomolecules often undergo high levels of structural changes during their biological functions, thereby resulting in distinct conformal states.In addition, both Gaussian and Boltzmann fitting accelerated molecular dynamics approaches can be used to simultaneously improve both the sampling and free energy calculations of these biomolecules.Therefore, our data, as shown in Figure 16, revealed that both Gaussian and Boltzmann fitting distributions were useful for further exploring the  enhancement of sampling and free energy calculations of the biomolecules.However, the findings of the present study suggest that the selected ligands might be favorable for total cure of rheumatoid arthritis.Therefore, additional clinical studies are needed to determine the underlying mechanism involved and possibly to validate whether any of the selected ligands could be effective in the treatment of rheumatoid arthritis (Table 4).

Conclusion
Our current study used a steadfast approach to identify potential immune protein inhibitors from a natural source (Piper sarmentosum Roxb leaf fractionated extract) as potential therapeutic agents for rheumatoid arthritis treatment.An in silico study conducted as part of this research has identified phyto-compounds 5280343 (quercetin), 5280443 (apigenin), 5280863 (kaempferol), and 5281601 (5-hydroxy-7,4-dimethoxyflavone) as promising inhibitors of these immune proteins.Furthermore, these phytocompounds that adhere to the Lipinski ROF and other significant parameters need to be considered drug candidates.However, several limitations and defects exist; for example, determining the accuracy and speed of docking calculations is challenging when exploring this approach could serve as a solution for most docking changes.Furthermore, the limited knowledge in academic computational models could be one of the challenges; however, further steps are needed to address these challenges via the integration of ligands and target of interest-, phenotype-, and biological network-based approaches with deeper reinforcement learning techniques, which could increase the predictive power of the models.However, in vivo and/or in vitro assays could further affirm the potency of these compounds as anti-inflammatory drugs for inhibiting IL-6 production.

TABLE 1 MM
-GBSA binding free energy (ΔG) calculations and binding sites of the top-scoring proteins from the Piper sarmentosum Roxb leaf fractionated extract.

TABLE 3
Molecular dynamics simulation properties of the native protein and ligands of interest.
Note: Values are represented as the mean ± SEM, measured in the Armstrong unit (Å).RMSD: root mean square deviation; RMSF: root mean square fluctuation; rGyr: radius of gyration; MolSA: molecular surface area; SASA: solvent-accessibility surface area; PSA: pressure swing adsorption.

TABLE 4
Principal component analysis of the native protein and ligands of interest.
Receptor (M)Ligand Eigen value cumulative Percentage of variance Coefficient of PC2 Coefficient of PC1