ORIGINAL RESEARCH article

Front. Pharmacol., 02 January 2025

Sec. Ethnopharmacology

Volume 15 - 2024 | https://doi.org/10.3389/fphar.2024.1380335

Integrated approach with UHPLC-Q-Exactive-tandem mass spectrometry, network analysis, and molecular docking to determine potential active compounds and mechanisms of Rhizoma Musae decoction in osteoarthritis treatment

  • JZ

    Jian Zhang 1*

  • WS

    Wanyan Shen 2

  • FL

    Fanzhi Liu 1

  • HH

    Hehe He 2

  • SH

    Shuquan Han 1

  • LL

    Lina Luo 1

  • 1. GuiZhou Institute of Subtropical Crops, Guizhou Academy of Agricultural Sciences, Guiyang, China

  • 2. Research and Development Department, Guizhou Weikang Zifan Pharmaceutical Co., Ltd., Guiyang, China

Article metrics

View details

1

Citations

1,6k

Views

515

Downloads

Abstract

Objective:

This study aimed to identify the potential active compounds in Rhizoma Musae decoction and understand their mechanisms of action in osteoarthritis treatment.

Methods:

UHPLC-Q-Exactive-MS/MS technology was used for an in-depth analysis of the chemical compounds present in Rhizoma Musae decoction. A network analysis approach was used to construct a comprehensive network of compounds, targets, and pathways, which provided insights into the molecular mechanisms of Rhizoma Musae decoction in osteoarthritis treatment.

Results:

The integrated analysis revealed the presence of 534 chemical compounds in Rhizoma Musae decoction, with 7beta-hydroxyrutaecarpine, 7,8-dihydroxycoumarin, pinocembrin diacetate, and scopoletin being identified as potential active compounds. Potential targets such as GAPDH, AKT1, TNF, IL6, and SRC were implicated in key pathways including MAPK signaling pathway, lipid and atherosclerosis, PI3K-Akt signaling pathway, and IL-17 signaling pathway. Molecular docking studies showed significant binding affinity between the core targets and key components. In vitro cell experiments have demonstrated that RM decoction can enhance cell proliferation and upregulates the expression of TNFα, IL-6, and SRC, while down-regulating the expression of GAPDH and AKT1.

Conclusion:

The potential active compounds present in Rhizoma Musae decoction influence specific targets and signaling pathways involved in osteoarthritis pathogenesis, providing new insights for the functional development and utilization of RM.

1 Introduction

Rhizoma Musae (RM), which is known as the dried root or rhizome of Musa basjoo Sied.et Zucc. and belongs to the Musa family, has been used as a traditional Chinese herb in Guizhou Province’s Miao region for centuries. The Quality Standards for Traditional Chinese and Ethnic Medicinal Materials in Guizhou Province (Guizhou Medical Products Administration, 2003) indicates that RM has multiple medicinal benefits (Fu et al., 2018), including heat-clearing, detoxifying, thirst-quenching, and diuretic properties. It is widely used to treat conditions such as wind-heat headaches, edema, athlete’s foot, thirst, metrorrhagia, vaginal discharge, erysipelas, and hematospermia. RM is rich in medicinal compounds such as essential oils, phenols, benzophenones, alkaloids, and acenaphthene derivatives, exhibiting a broad spectrum of pharmacological activities, such as anti-inflammatory (Wang et al., 2010), analgesic (Liang et al., 2010), antibacterial (Wei et al., 2010), and antioxidant effects (Liu et al., 2013).

Osteoarthritis (OA), a joint disorder characterized by the progressive degeneration of the articular cartilage (Yao et al., 2023), is on the rise due to factors such as the aging population, increasing obesity rates, and high incidence of traumatic knee injuries. This suggests a forthcoming increase in the occurrence and frequency of OA (Chen et al., 2020). Mechanical damage, inflammation, and aging are significant contributors to the development of osteoarthritis (OA), leading to the gradual deterioration of cartilage (Ouyang et al., 2023). In the pathogenesis of OA, an imbalance of cytokines can trigger the repeated activation of pro-inflammatory cytokines, resulting in damage to internal joint structures such as cartilage (Molnar et al., 2021). Pro-inflammatory cytokines like IL-33, IL-17, IL-6, and IL-22 play crucial roles in various signaling pathways within joints (Englund, 2023). Low-grade inflammation is a notable feature in the early stages of OA (Berenbaum, 2013). The increase in inflammatory factors prompts the breakdown of chondrocyte matrix and contributes to the degeneration of OA chondrocytes. This irreversible process further enhances the release of inflammatory factors, accelerating disease progression.

Network analysis is a multidisciplinary field that integrates systems biology, computer science, and bioinformatics. It employs high-throughput technology to investigate the correlation between diseases and drugs efficiently and economically. The approach aligns with traditional Chinese medicine formulas, known for their attributes targeting multiple levels (Zhao et al., 2023). The Gukang capsule, which primarily contains RM, has shown significant effectiveness in treating osteoporosis, fractures, and injuries (Zhu et al., 2022). The aqueous extract of RM promotes the proliferation, differentiation, and mineralization of bone cells in laboratory settings; upregulates alkaline phosphatase levels; and improves blood rheology, thereby reducing blood viscosity and supporting the fracture healing process (Dong et al., 2019). Despite these findings, a comprehensive understanding of the medicinal components of RM decoction and their effects on OA signaling pathways remains elusive. In this study, UHPLC-Q-Exactive-MS/MS technology is utilized to analyze the chemical constituents of RM decoction. Additionally, network analysis and molecular docking techniques are employed to preliminarily explore the intervention mechanism of RM decoction against OA. These findings establish a robust foundation for further clinical studies and the development of drug formulations involving RM.

2 Materials and methods

2.1 Experimental instruments

Ultra-high-pressure liquid chromatography was performed using a Vanquish system (Thermo Fisher Scientific, United States). Mass spectrometry analysis was conducted with a Q Exactive HFX Mass Spectrometer (Thermo Fisher Scientific, United States). Centrifugation steps were conducted using Centrifuge 5430 R (Eppendorf, Germany). Mixing was achieved with a SCI-VS vortex mixer (Scilogex, United States), and an ultrasonic cleaner SB25-12DTD (Ningbo Xinzhi Biotechnology Co., Ltd., China) was used for sonication.

2.2 Plant materials and chemicals

The RM samples were provided by Guizhou Weikang Zifan Pharmaceutical Co., Ltd., China, with the batch number YC01-003-20230801. Solvents including methanol, acetonitrile, formic acid, isopropanol (all of analytical grade from Anpu), and 95% ethanol (of analytical grade from Shanghai Shenggong) were used. Ultrapure water was utilized throughout all procedures.

2.3 Analysis of the chemical composition of RM decoction

2.3.1 Sample preparation

Referred to Gao’s report (Gao et al., 2019), RM decoction was prepared as follows: 10 g of RM samples were powdered and added to a container with 100 mL of ultrapure water. The mixture was then brought to a boil and simmered at low heat, maintaining a slight boil, for 1 h. After cooling for 1 h, 1 mL of the supernatant was mixed with 2 mL of a methanol-acetonitrile solution (1:1, v/v), vortexed for 60 s, and then sonicated at low temperature for 30 min. The mixture was centrifuged at 12,000 rpm for 10 min at 4°C, and the supernatant was collected and chilled at −20°C for 1 h to precipitate the proteins. This process was followed by another centrifugation under the same conditions, after which the supernatant was freeze-dried, resuspended in 100 μL of 50% acetonitrile, vortexed, and centrifuged again at 12,000 rpm for 10 min at 4°C to collect the final supernatant for analysis.

2.3.2 Detection conditions

  • (1) Chromatographic analysis was performed using a UPLC-Orbitrap-MS system. The UPLC conditions were set with a Waters HSS T3 column (100 × 2.1 mm, 1.8 μm), a column temperature of 40°C, a flow rate of 0.3 mL/min, and an injection volume of 2 μL. The solvent system consisted of water with 0.1% formic acid and acetonitrile with 0.1% formic acid. The gradient program started at 0% B for 1 min, increased to 95% B over 12 min, held at 95% B for 1 min, then returned to 0% B by 17 min.

  • (2) Mass spectrometric analysis included the use of a Q Exactive HFX Hybrid Quadrupole Orbitrap mass spectrometer with a heated ESI source utilizing the Full-ms-ddMS2 MS acquisition methods. Parameters for the ESI source were set to a spray voltage of −2.8 kV/3.0 kV, sheath gas pressure at 40 arb, aux gas pressure at 10 arb, sweep gas pressure at 0 arb, capillary temperature at 320°C, and aux gas heater temperature at 350°C. The scan range was set from 70 to 1050 Da, with a resolution of 70,000 at the first level and 17,500 at the second level. The compounds were characterized using Sanshu Biotech’s proprietary database, which is specialized in Traditional Chinese Medicine (TCM), along with their unique secondary mass spectrometry fragmentation pattern matching technique. The relative contents of each compound in RM decoction were determined using the area normalization method.

2.4 Network analysis research

2.4.1 Obtaining targets related to RM decoction

A search was performed in the NCBI PubChem database (https://pubchem.ncbi.nlm.nih.gov/) to identify compounds and extract their structural information related to RM decoction. Isomeric SMILE structural formulas of compounds associated with RM decoction were inputted into the Swiss Target Prediction Database (http://www.swisstargetprediction.ch/) to identify their targets. To identify the target protein of the selected compound, a screening principle based on a predicted target probability greater than 0 was employed. The UniProt database’s ID mapping tool (https://www.uniprot.org/) was then used to convert UniProt IDs into gene symbols, and after removing duplicates, the targets of the potential active compound were successfully identified.

2.4.2 Identification of OA-related target genes

The keyword “osteoarthritis” was used to identify target genes by extracting relevant data from three disease gene databases: GeneCards (https://www.genecards.org/), Comparative Toxicogenomics Database (CTD) database (https://ctdbase.org/), and DisGeNET database (https://www.disgenet.org/). Genes associated with OA were retrieved from these databases, and the collected information was organized to identify disease targets.

2.4.3 Common targets screening and PPI network construction

The online software program Venny (https://bioinfogp.cnb.csic.es/tools/venny/) was used to predict the common targets between RM decoction and OA. Subsequently, the cross-genes were analyzed using the STRING data analysis platform (http://string-db.org/) to create a protein-protein interaction (PPI) network of common targets, applying a minimum interaction score threshold of ≥0.4. The PPI network diagram and data were saved in TSV format. Cytoscape software (version 3.9.0) was then utilized to visualize the PPI network and construct a multidimensional network depicting the relationship between RM decoction compounds and OA.

2.4.4 Analysis of gene ontology and Kyoto Encyclopedia of Genes and Genomes pathways

In-depth analyses were performed using R software to explore the potential biological functions and key signaling pathways of RM decoction in treating OA. Our comprehensive enrichment analyses included Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, with a significance criterion set at a q-value of <0.05. The results were organized by descending p-values, highlighting significant differences in enrichment.

2.4.5 Molecular docking

The core target-related protein structures (in PDB format) for intersection target screening were obtained from the RCSB database (https://www.rcsb.org/), and the core structure of the potential active compound was sourced from the PubChem database. AutoDock software was used to preprocess the target proteins and small molecule compounds and convert their formats to analyze the protein binding sites. The compounds were then docked using AutoDock Vina 1.1.2 software, with defined docking site coordinates and docking verification. A binding energy value below zero indicates spontaneous binding affinity, with values below −5 kcal·mol−1 signifying strong binding activity and those below −7 kcal·mol−1 indicating intense binding activity (Hsin et al., 2013).

2.5 In vitro cell experiments

2.5.1 Cell culture

After resuscitating human osteoblast hFOB1.19 cells (Shanghai Cell Bank, Chinese Academy of Sciences, China), the cells were passaged when they reached 80%–95% confluence. To resuspend the cells, 0.25% Trypsin-EDTA medium (Thermo Fisher, United States), preheated to 37°C, was added. The resuspended cells were then transferred to a new culture dish at a density of 5 × 105 cells/mL and placed in a cell culture incubator set to 37°C with 5% CO2.

2.5.2 Cell viability assay

Collect hFOB1.19 cells that have reached the logarithmic phase of growth and treat them with RM decoction at concentrations of 0, 60, 80, 100, 120, 140, and 160 μg/mL, while adding blank culture medium to the control group. The cells should be seeded in a 96-well plate, with three replicates for each group, and each well should contain 3 × 10³ cells. After culturing for 48 h, add 10 μL of CCK-8 reagent to each well and incubate for an additional 1.5 h. Subsequently, use a microplate reader to measure the absorbance (A) at 450 nm, and calculate the cell proliferation rate for each experimental group.

2.5.3 Quantitative real-time polymerase chain reaction(qPCR)

Total RNA was extracted from cells using an TRIzol Reagent (Life Technologies, United States), followed by the synthesis of cDNA (cDNA).The qPCR experiments were conducted using CFX384 Multiple real-time fluorescent quantitative PCR apparatus (Bio-Rad, United States). The PCR amplification conditions were set at 95°C for 30 s, followed by 40 cycles of 95°C for 5 s and 60°C for 30 s. The relative expression level of mRNA was calculated using β-actin as the internal reference gene and the 2−ΔΔCT method. The qPCR primer design was performed using Primer Premier 6.0 and Beacon Designer 7.8 software. The primers were synthesized by Sangon Bioengineering (Shanghai), with the sequences were listed in Table 1.

TABLE 1

Gene nameForward primerReverse primer
Human β-actinGAT​GAC​CCA​GAT​CAT​GTT​TGA​GACGGA​GTC​CAT​CAC​GAT​GCC​AGT
Human GAPDHCCA​TGA​CAA​CTT​TGG​TAT​CGT​GGA​AGGC​CAT​CAC​GCC​ACA​GTT​TC
Human AKT1GCC​CCA​CTT​CCC​CCA​GTT​CTCCG​CCT​CTC​CAT​CCC​TCC​AA
Human TNFαCCA​TGT​TGT​AGC​AAA​CCC​TCA​AGC​TCCT​TGA​AGA​GGA​CCT​GGG​AGT​AGA​T
Human IL-6CCT​TCG​GTC​CAG​TTG​CCT​TCTGTG​TGG​GGC​GGC​TAC​ATC​TTT
Human SRCCTG​CTT​TGG​CGA​GGT​GTG​GAT​GCCA​CAG​CAT​ACA​ACT​GCA​CCA​G

Primers sequences.

2.5.4 Statistical analysis

Data were analyzed using GraphPad Prism version 10.1.2. All results are presented as mean ± SEM. An independent unpaired two-tailed t-test was employed to compare the two groups.

3 Results

3.1 Identification of compounds from RM decoction

Overall, 534 compounds were identified using Sanshu Biotech’s proprietary database, specialized in TCM, and their unique secondary mass spectrometry fragmentation pattern matching technique (Supplementary Table S1). Of these, 299 were detected in the positive ion mode (Figure 1A) and 235 in the negative ion mode (Figure 1B). These compounds were categorized into various groups, including flavonoids (66), prenol lipids (54), organooxygen compounds (50), carboxylic acids and derivatives (Rose-John, 2012), coumarins and derivatives (Lin et al., 2022), fatty acyls (Jiménez-Orozco et al., 2020), benzene and substituted derivatives (Wu et al., 1995), isoflavonoids (Li, 2004), cinnamic acids and derivatives (Gao et al., 2019), phenols (Berenbaum, 2013), benzopyrans (Englund, 2023), anthracenes (Chen et al., 2020), diarylheptanoids (Yao et al., 2023), purine nucleosides (Yao et al., 2023), steroids and steroid derivatives (Yao et al., 2023), and 172 other compounds.

FIGURE 1

3.2 Identification and PPI network analysis of common targets against OA

The top 60 compounds from the RM decoction were selected based on their relative peak areas (Table 2). A total of 737 compound targets were accessed from the Swiss Target Prediction database, and 1174, 6508, and 368 OA-related targets were identified from the GeneCards, CTD, and DisGeNET databases, respectively. Through target gene intersection, the present study identified 507 potential targets for OA treatment (Figure 2).

TABLE 2

No.CompoundFormulaRt (min)m/zIon modePubChem ID
1BetaineC5H11NO20.94118.0864POS247
23-O-Caffeoylquinic acidC16H18O95.75353.0874NEG1794427
3Cholinesulfuric acidC5H13NO4S0.93184.0636POS485
4Eleutherazine BC22H36N4O67.22453.2698POS20839739
5ADENOSINEC10H13N5O44.06268.1035POS60961
67beta-HydroxyrutaecarpineC18H12N3O212.87303.1010POS15225951
7Ginsenoside Rg1C42H72O1410.54823.4803POS441923
84-Hydroxybenzoic acidC7H6O36.70137.0233NEG135
9TrigonellineC7H7NO21.01138.0549POS5570
10Vanillic acidC8H8O47.24167.0341NEG8468
11D-mannitolC6H14O60.93181.0709NEG6251
12CarnitineC7H15NO30.94162.1123POS288
132,16-Kauranediol 2-O-beta-D-allopyranosideC26H44O712.43491.2970POS73554066
14ScopoletinC10H8O48.09237.0400NEG5280460
154-Hydroxyphenylpyruvic acidC9H8O45.85163.0389POS979
16Caffeic acidC9H8O46.80163.0389POS689043
17Suavioside AC26H44O811.54507.2922POS73821014
18D-LEUCINEC6H13NO22.01132.1019POS439524
19VanillinC8H8O37.76151.0390NEG1,183
20D-Glutamic acidC5H9NO40.93148.0603POS23327
21Myo-InositolC6H12O61.06179.0552NEG892
22GuanosineC10H13N5O52.05284.0983POS135398635
23D-prolineC5H9NO21.03116.0708POS8988
24SugerosideC26H42O811.74505.2764POS3082543
25OpuntiolC7H8O41.06174.0760POS10034839
26XylitolC5H12O50.94151.0601NEG6912
27URIDINEC9H12N2O62.79243.0618NEG6029
28D-ValineC5H11NO20.94235.1648POS71563
29Gallic acidC7H6O53.11169.0133NEG370
307,8-DihydroxycoumarinC9H6O47.13223.0243NEG5280569
31Quinic acidC7H12O63.84173.0446NEG6508
32LucidoneC15H12O49.51301.0714NEG11253859
335,7-DihydroxyphthalideC8H6O46.15211.0242NEG11062751
343,4-Dihydroxyphenylacetic acidC8H8O45.09167.0341NEG547
35D-ArabinoseC5H10O50.95133.0495POS854
36Alpha-IsowighteoneC20H18O512.08303.1009POS91885205
37AconineC25H41NO912.44517.3130POS20054813
38Syringolin AC24H39N5O612.18535.3238POS42601513
39N-Isobutyl-2,4,12-octadecatrienamideC22H39NO12.91378.2750POS25221579
40StemonidineC19H29NO513.29332.1863NEG24721470
41[(3,8,12-trihydroxy-24-oxocholan-24-yl)amino]acetateC26H42NO610.54509.2716POS11834768
42Giffonin RC19H16O310.74257.0956POS134715258
434-AMINOBUTYRIC ACIDC4H9NO20.94104.0710POS119
44Linoleic acidC18H32O214.89313.2729POS5280450
45Paeonilactone AC10H14O46.75216.1228POS10081437
463-Hexen-1-ol O-b-D-glucopyranosideC12H22O69.25261.1340NEG5318045
47D-GlucosamineC6H13NO50.96180.0864POS439213
48Pierisformoside BC26H42O811.09505.2763POS155978780
493′-HydroxyxanthyletinC14H12O410.30227.0699POS129069558
50AdenineC5H5N51.30136.0618POS190
51PreisocalamendiolC15H24O10.54203.1793POS12305706
52GamabufotalinC24H34O512.24403.2450POS259803
53LongicaulenoneC12H18O47.91249.1108POS25750965
547-DeoxyechinosporinC10H9NO46.96190.0497POS11745823
55Notoginsenoside R1C47H80O1810.22955.5221POS131752529
564-EpialyxialactoneC10H16O410.08199.0969NEG14194344
57Pinocembrin diacetateC19H16O69.42385.0925NEG6546286
58Dehydroacerogenin CC19H18O311.52259.1113POS154790969
59Excavatin MC19H20O77.81359.1133NEG15871351
60L (+)-Ascorbic acidC6H8O61.23221.0297NEG54670067

Relative content of the top 60 compounds in Rhizoma Musae decoction.

FIGURE 2

These targets were analyzed using the String database, yielding 507 proteins and 10,739 interaction lines. A drug-disease PPI network diagram was created (Figure 3), highlighting the complex protein relationships. Using Cytoscape software, 103 core targets (Supplementary Table S2) were identified based on the specific criteria of Degree≥42.82936508, betweenness centrality ≥0.002366375, and closeness centrality ≥0.464818506. Among the core targets, GAPDH, AKT1, TNF, IL6, SRC, EGFR, STAT3, HSP90AA1, BCL2, and JUN (Table 3) had higher degrees, suggesting that the compounds present in RM decoction may collectively act on these core targets, leading to its pharmacological effects on OA.

FIGURE 3

TABLE 3

TargetsDegreeBetweenness centralityCloseness centrality
GAPDH2690.0735240.680650
AKT12550.0492350.666225
TNF2440.0450670.655802
IL62370.0422810.651554
SRC2150.0558170.624845
EGFR2140.0295350.628750
STAT32110.0233260.627182
HSP90AA11850.0218140.608222
BCL21820.0112130.602395
JUN1810.0183990.602395

Topological parameter analysis of core targets (top 10).

3.3 Elaboration of the “RM decoction-OA” network diagram

A network diagram of the RM decoction against OA was then constructed to further explore the potential association between the core components of RM decoction (Figure 4). The network featured 40 active compound nodes, 103 target nodes, and one disease node. The topological characteristics of this network were analyzed using the Network Analyzer plug-in of Cytoscape software 3.9.0. The potential active compounds with the highest degree values were identified, including 7beta-Hydroxyrutaecarpine (degree = 22), 7,8-dihydroxycoumarin (degree = 21), pinocembrin diacetate (degree = 19), Scopoletin (degree = 19), and Eleutherazine B (degree = 18), among others, detailed in Table 4.

FIGURE 4

TABLE 4

No.CompoundsDegreeBetweenness centralityCloseness centralityRelative content (%)
17beta-Hydroxyrutaecarpine220.0233800.4005602.54 ± 1.10
27,8-Dihydroxycoumarin210.0153590.3833780.68 ± 0.04
3Pinocembrin diacetate190.0163060.3896460.29 ± 0.06
4Scopoletin190.0128800.3773091.87 ± 0.20
5Eleutherazine B180.0110660.3733683.11 ± 0.27
6Excavatin M180.0119910.3773090.29 ± 0.01
7linoleic acid180.0148200.3833780.44 ± 0.08
8Syringolin A170.0129440.3833780.50 ± 0.04
9Stemonidine160.0085400.3657290.48 ± 0.19
103-Hexen-1-ol O-b-D-glucopyranoside150.0091150.3753280.38 ± 0.04

The potential active compounds of Rhizoma Musae decoction against osteoarthritis (top 10).

3.4 GO analysis and KEGG pathway enrichment analysis

The “Cluster Profiler” software package for the R language facilitated GO and KEGG enrichment analysis on intersection targets, yielding 2821 entries from the GO analysis (Supplementary Table S3). These included 2536 entries in the Biological Process (BP) category, 182 in the molecular function (MF) category, and 103 in the cellular component (CC) category. Figure 5 presents the top 5 BPs, identified based on the lowest p-values: neuron death, cellular response to chemical stress, response to oxidative stress, regulation of neuron death, and cellular response to oxidative stress. The leading MFs were transcription coregulator binding, RNA polymerase II-specific DNA-binding transcription factor binding, ubiquitin protein ligase binding, DNA-binding transcription factor binding, and ubiquitin-like protein ligase binding. The primary CCts included membrane raft, membrane microdomain, neuronal cell body, ficolin-1-rich granule lumen, and ficolin-1-rich granule. KEGG pathway analysis identified 187 signaling pathways (Supplementary Table S4), such as lipid and atherosclerosis, PI3K-Akt signaling pathway, hepatitis B, proteoglycans in cancer, and chemical carcinogenesis–receptor activation, as detailed in Figure 6. An RM decoction component-target-pathway network model was systematically constructed, as displayed in Figure 7.

FIGURE 5

FIGURE 6

FIGURE 7

3.5 Molecular docking results

Utilizing Diacerein and corresponding target inhibitors as positive controls, molecular docking was performed between the initial five core components and the first five identified targets, followed by the calculation of binding energy. A binding energy below zero indicates spontaneous binding affinity, with more negative values suggesting greater binding stability. The docking results, illustrated in Table 5, demonstrated spontaneous binding interactions between the core components and the five identified targets. Among the 15 groups, those with binding energy less than −5 kcal·mol−1 demonstrated strong binding activity, and two groups with less than −7.0 kcal·mol−1 exhibited intense binding activity. Eleutherazine B demonstrates a greater binding affinity towards GAPDH, AKT1, and SRC compared to Diacerein and their respective target inhibitors. Additionally, 7beta-Hydroxyrutaecarpine shows a higher binding capacity to SRC than Diacerein, while Pinocembrin diacetate exhibits a stronger binding affinity to AKT1 than Diacerein.

TABLE 5

TargetsCompounds
NameBinding energy (kcal·mol−1)
GAPDH7,8-Dihydroxycoumarin−4.84
7beta-Hydroxyrutaecarpine−5.79
Eleutherazine B−8.36
Pinocembrin diacetate−6.64
Scopoletin−5.20
Diacereina−7.16
heptelidic acidb−6.44
AKT17,8-Dihydroxycoumarin−4.17
7beta-Hydroxyrutaecarpine−4.97
Eleutherazine B−6.57
Pinocembrin diacetate−5.34
Scopoletin−4.16
Diacereina−5.09
BAY 1125976b−5.81
TNF7,8-Dihydroxycoumarin−4.60
7beta-Hydroxyrutaecarpine−5.44
Eleutherazine B−7.04
Pinocembrin diacetate−5.85
Scopoletin−4.61
Diacereina−7.27
Pomalidomideb−6.00
IL67,8-Dihydroxycoumarin−4.86
7beta-Hydroxyrutaecarpine−5.16
Eleutherazine B−6.68
Pinocembrin diacetate−5.69
Scopoletin−5.13
Diacereina−6.22
Curcuminb−6.83
SRC7,8-Dihydroxycoumarin−4.40
7beta-Hydroxyrutaecarpine−5.20
Eleutherazine B−6.69
Pinocembrin diacetate−5.08
Scopoletin−4.08
Diacereina−5.13
Saracatinibb−6.05

Predicted binding energies of potential active compounds to core targets combined with positive controls.

a

Positive control drug.

b

Target Inhibitors.

Figure 8 methodically illustrates the visualization of docking conformations with the highest binding energies. 7,8-Dihydroxycoumarin forms a hydrogen bond with IL6 at specific sites, resulting in a binding energy of −4.86 kcal·mol−1, indicating robust interaction strength (Figure 8A). The interaction between 7beta-Hydroxyrutaecarpine and GAPDH primarily involves hydrophobic forces, with a binding energy of −5.79 kcal·mol−1 (Figure 8B). Eleutherazine B interacts with GAPDH through hydrogen bonds at ARG P:80 and ASP P:35, exhibiting a binding energy of −8.36 kcal·mol−1, denoting an intense interaction (Figure 8C). Pinocembrin diacetate forms significant hydrogen bonds with GAPDH at ALA P:183 and ASN P:316 reflected in a binding energy of −6.64 kcal·mol−1, indicative of a potent interaction (Figure 8D). Scopoletin engages in hydrogen bonding with GAPDH at ALA P:183 and ASN P:316, displaying a binding energy of −5.20 kcal·mol−1, signifying a strong interaction (Figure 8E).

FIGURE 8

3.6 RM decoction promotes the proliferation of hFOB1.19 cells in vitro

In Figure 9, treatment with RM decoction for 48 h resulted in changes to the proliferation rate of hFOB1.19 cells, which varied with the concentration of the decoction. Notably, at a concentration of 140 μg/mL, the proliferation rate of human osteoblasts was the highest, significantly exceeding that of the control group (p < 0.05).

FIGURE 9

3.7 The effects of RM decoction on the expression of related genes in hFOB1.19 cells

Figure 10 illustrates the alterations in gene expression following the treatment of hFOB1.19 cells with 140 μg/mL RM decoction for a duration of 48 h. Notably, RM decoction significantly upregulates the expression of TNFα, IL-6, and SRC, while down-regulating the expression of GAPDH and AKT1.

FIGURE 10

4 Discussion

RM, documented in the Compendium of Materia Medica (Li, 2004), has been utilized in Miao medicine for centuries. This study employs UHPLC-Q-Exactive-MS/MS technology to identify and analyze 534 chemical components present in RM decoction, with a majority being flavonoids, prenol lipids, organooxygen compounds, and others. Furthermore, a component-target correlation analysis on the top 60 chemical components of RM decoction aids in identifying key components for OA treatment.

The decoction method, which involves boiling medicinal materials in water to extract their juices, is a simple form of leaching (Yang and Gao, 2021). This method, considered one of the simplest forms of leaching, is best suited for medicinal materials with active compounds that are soluble in water and can withstand exposure to moisture and heat (Jia et al., 2022). Through mass spectrometry and network analysis, the key components of RM decoction were identified as 7beta-hydroxyrutaecarpine, 7,8-dihydroxycoumarin, Pinocembrin diacetate, and Scopoletin. 7beta-Hydroxyrutaecarpine, derived from evodiacarpine, is a hydroxylation product and a novel indopyridine-quinazoline alkaloid initially discovered in Phellodendron callus tissue (Wu et al., 1995). Experimentally verified, 7beta-Hydroxyrutaecarpine exhibits a range of pharmacological effects, such as analgesic, anti-inflammatory, anti-ulcer, antiemetic, and antidiarrheal properties. (Hu et al., 2018). Another component, 7,8-dihydroxycoumarin, a natural coumarin found in various edible plants of the Thymelaeaceae family, such as Daphne Korean Nakai, Daphne gnidium, Daphne odora, and Daphne oleoides exhibits diverse biological activities, particularly on anti-arthritic (Jiménez-Orozco et al., 2020). Pinocembrin diacetate, a prominent flavonoid, is widely studied for its antimicrobial, anti-inflammatory, antioxidant, and anticancer properties (Rasul et al., 2013). Scopoletin, a naturally occurring coumarin, falls into the category of secondary metabolites and possesses potential anti-inflammatory and anti-tumorigenesis properties (Sakthivel et al., 2022). Diacerein is a slow-acting drug used for osteoarthritis treatment, with anti-inflammatory, antioxidant, and anti-apoptotic properties (Almezgagi et al., 2020). It works by inhibiting the IL-1b system and its downstream signaling pathways (Panova and Jones, 2015). Molecular docking analysis showed potential active compounds in RM decoction, such as 7beta-hydroxyrutaecarpine, 7,8-dihydroxycoumarin, Pinocembrin diacetate, Eleutherazine B, and Scopoletin, effectively binding to the target protein. Organic extraction is a common method for preparing medicinal materials (Lin et al., 2022). Studies have confirmed the anti-inflammatory and analgesic properties of lupenone and β-sitosterol in the ethyl acetate fraction of Rhizoma Musae (Xu et al., 2014). Further research is necessary to explore the potential effectiveness of RM organic extract in osteoarthritis treatment.

The analysis of the PPI network indicates that key proteins such as GAPDH, AKT1, TNF, IL6, and SRC, which have high node Degree values and multiple connections to potential active components, are likely central targets for RM decoction in OA treatment. GAPDH, integral to glycolysis, significantly influences cell proliferation (Sirover, 2021) and apoptosis (Sen et al., 2009). AKT1, a subset of the AKT serine/threonine kinase family encoded by the PKB gene, plays a pivotal role in the PI3K pathway, influencing various downstream effectors. This PI3K/AKT signaling pathway is involved in the abnormal proliferation of bone cells, synovial inflammation, and the formation and differentiation of osteoclasts (Hayer et al., 2009), leading to bone and articular cartilage damage and joint deformity. TNF, which encodes a multifunctional proinflammatory cytokine in the TNF superfamily, is linked to diseases such as autoimmune disorders, insulin resistance, psoriasis, and rheumatoid arthritis (Nie et al., 2013). IL-6, a versatile cytokine, plays roles in immune and nervous system regulation and is associated with antimicrobial molecule production and cytokine activity, as evidenced by research (Rose-John, 2012). SRC, vital in bone resorption by osteoclasts and inhibiting bone formation by osteoblasts (Matsubara et al., 2022), emerges as a potential therapeutic target in osteoporosis treatment. Comparing our results with previous studies suggests that RM decoction could potentially offer anti-inflammatory benefits and promote bone cell proliferation by targeting these specific mechanisms.

The GO analysis covers a range of biological processes, such as neuron death, cellular response to chemical stress, positive regulation of kinase activity, and the positive regulation of the MAPK cascade, among others. The MAPK signaling pathway, critical for various cellular functions including apoptosis, differentiation, and cell proliferation, operates through the activation of nuclear transcription factors (Zheng et al., 2019). The MAPK family, consisting of members such as JNK, ERK, and p38, serves diverse roles. It has been established that ERK1/2 promotes cell proliferation (Chen et al., 2014), while p38 supports cell differentiation. KEGG enrichment analysis has identified the involvement of core targets across 20 signaling pathways, with particular relevance to diseases found in lipid and atherosclerosis, PI3K-Akt signaling pathway, and IL-17 signaling pathway. The lipid and atherosclerosis pathway, for example, involves several inflammation-related targets such as AKT1, TNF, and IL6, underscoring chronic inflammation’s critical role in atherosclerosis progression (Poznyak et al., 2020). This is consistent with literature highlighting atherosclerosis as primarily a chronic inflammatory condition. The PI3K/AKT signaling pathway is pivotal for cell proliferation, differentiation, invasion, and apoptosis (Ma et al., 2014), with 31 core gene targets identified in this study, including AKT1, IL6, and EGFR, among others. These targets are essential components closely associated with the development of inflammatory diseases and cancer progression. IL-17, a highly adaptable cytokine known for its pro-inflammatory capabilities, is linked to genes encoding antimicrobial agents and other cytokines (Tsai et al., 2013). In this research, the IL-17 pathway involves 16 genes, such as TNF, HSP90AA1, and MMP13, significant for host defense, tissue regeneration, and the onset of inflammatory conditions.

5 Conclusion

This study aimed to investigate the chemical components and mechanisms of RM decoction in treating OA. The potential components of RM decoction were identified as 7beta-hydroxyrutaecarpine, 7,8-dihydroxycoumarin, Pinocembrin diacetate, and Scopoletin. Network analysis indicates that GAPDH, AKT1, TNF, IL6, and SRC may be important target proteins. These potential active compounds may regulate pathways like the MAPK signaling pathway, Lipid and Atherosclerosis, PI3K-Akt, and IL-17 signaling pathway, potentially offering anti-inflammatory effects and promoting osteocyte proliferation and differentiation. The development of diseases and the pharmacological effects of bioactive substances are complex and constantly changing processes. It is crucial to validate these findings with further experimental and clinical research. Despite the limitations of our study, we have gathered valuable data on the chemical components and mechanism of action of RM in treating OA, offering new perspectives for the functional advancement and utilization of RM.

Statements

Data availability statement

The original contributions presented in the study are included in the article/Supplementary material, further inquiries can be directed to the corresponding author.

Author contributions

JZ: Funding acquisition, Resources, Writing–original draft, Writing–review and editing, Formal Analysis. WS: Data curation, Formal Analysis, Investigation, Writing–review and editing. FL: Formal Analysis, Funding acquisition, Writing–review and editing. HH: Funding acquisition, Supervision, Writing–review and editing. SH: Formal Analysis, Writing–review and editing. LL: Formal Analysis, Writing–review and editing.

Funding

The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. The Science and Technology Planning Project of Guizhou Province ([2024]YB091). General fund of Guizhou Academy of Agricultural Sciences ([2024]10). China Agriculture Research System of MOF and MARA (CARS-31). Guizhou Academy of Agricultural Sciences, Research on Innovation and Efcient Key Technologies of Characteristic Crop Germplasm in Guizhou Hot Zone, (Qiannongke Germplasm Resources [2024] No. 08). Guizhou Provincial Programme to Support High-Quality Creation and Application of Intellectual Property Rights ([2023]10). The Science and Technology Support Project of Guizhou Province ([2020]1Y024).

Acknowledgments

We thank Bullet Edits Limited for the linguistic editing and proofreading of the manuscript.

Conflict of interest

Authors WS and HH were employed by Guizhou Weikang Zifan Pharmaceutical Co., Ltd.

The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphar.2024.1380335/full#supplementary-material

References

  • 1

    AlmezgagiM.ZhangY.HezamK.ShamsanE.GamahM.Al-ShaebiF.et al (2020). Diacerein: Recent insight into pharmacological activities and molecular pathways. Biomed. Pharmacother.131, 110594. 10.1016/j.biopha.2020.110594

  • 2

    BerenbaumF. (2013). Osteoarthritis as an inflammatory disease (osteoarthritis is not osteoarthrosis!). Osteoarthr. Cartil.21, 1621. 10.1016/j.joca.2012.11.012

  • 3

    ChenJ.CaoJ.LuoY.XieL.SongJ.XueW.et al (2014). Expression of ERK and p-ERK proteins of ERK signaling pathway in the kidneys of fluoride-exposed carp (Cyprinus carpio). Acta histochem.116 (8), 13371341. 10.1016/j.acthis.2014.08.005

  • 4

    ChenL.YaoF.WangT.LiG.ChenP.BulsaraM.et al (2020). Horizontal fissuring at the osteochondral interface: a novel and unique pathological feature in patients with obesity-related osteoarthritis. Ann. Rheum. Dis.79 (6), 811818. 10.1136/annrheumdis-2020-216942

  • 5

    DongL.XuM.ChenJ. H.GuoT.DongY. X. (2019). Screening of effective components of Banana root promoting osteoblast proliferation and differentiation. China & Foreign Med. Treat.38 (18), 1316. 10.16662/j.cnki.1674-0742.2019.18.013

  • 6

    EnglundM. (2023). Osteoarthritis, part of life or a curable disease? A bird's-eye view. J. Intern Med.293 (6), 681693. 10.1111/joim.13634

  • 7

    FuS. H.PengX.JiangL.LiuT.DongL.XlX.et al (2018). Chemical constituents from rhizomes of Musa basjoo. J. Chin. Med. Mater41 (3), 98102. 10.13863/j.issn1001-4454.2018.03.019

  • 8

    GaoS.LinJ.WangT.ShenY.LiY.YangW.et al (2019). Qingxin kaiqiao fang ameliorates memory impairment and inhibits apoptosis in APP/PS1 double transgenic mice through the MAPK pathway. Drug Des. Dev. Ther.13, 459475. 10.2147/DDDT.S188505

  • 9

    Guizhou Medical Products Administration (2003). Quality standards of traditional Chinese medicine and ethnic medicine in Guizhou Province,Guiyang, China. Guiyang, China: Guizhou Science and Technology Press, 366.

  • 10

    HayerS.PundtN.PetersM. A.WunrauC.KühnelI.NeugebauerK.et al (2009). PI3Kgamma regulates cartilage damage in chronic inflammatory arthritis. FASEB J.23 (12), 42884298. 10.1096/fj.09-135160

  • 11

    HsinK. Y.GhoshS.KitanoH. (2013). Combining machine learning systems and multiple docking simulation packages to improve docking prediction reliability for network pharmacology. PLoS One8 (12), e83922. 10.1371/journal.pone.0083922

  • 12

    HuX.LiD.ChuC.LiX.WangX.JiaY.et al (2018). Antiproliferative effects of alkaloid Evodiamine and its derivatives. Int. J. Mol. Sci.19 (11), 3403. 10.3390/ijms19113403

  • 13

    JiaD.LiuC.ZhuZ.GaoY.WenW.HongZ.et al (2022). Novel transketolase inhibitor oroxylin A suppresses the non-oxidative pentose phosphate pathway and hepatocellular carcinoma tumour growth in mice and patient-derived organoids. Clin. Transl. Med.12 (11), e1095. 10.1002/ctm2.1095

  • 14

    Jiménez-OrozcoF. A.RanđelovićI.HegedüsZ.Vega-LopezA.Martínez-FloresF.TóvaríJ. (2020). In vitro anti-proliferative effect and in vivo antitumor action of daphnetin in different tumor cells. Cir. Cir.88 (6), 765771. 10.24875/CIRU.20000197

  • 15

    LiS. Z. (2004) “Compendium of Materia Medica,” in Chifeng: Inner Mongolia. Chifeng, China: Science and Technology Press, 30.

  • 16

    LiangY. Q.YangL. B.LouF. Y.WangX. P.LouF. Y.WuH. M. (2017). Comparative studies of ultraviolet spectrum group and anti-inflammatory analgesic effect in rhizome, pseudostem and leaf of Musa basjoo. Lishizhen Med. Mater Med. Res.28, 545546. 10.3969/j.issn.1008-0805.2017.03.012

  • 17

    LinS.MengX.TanC.TongY.WanM.WangM.et al (2022). Composition and antioxidant activity of anthocyanins from Aronia melanocarpa extracted using an ultrasonic-microwave-assisted natural deep eutectic solvent extraction method. Ultrason. Sonochem89, 106102. 10.1016/j.ultsonch.2022.106102

  • 18

    LiuY.LuQ.ZhouZ. Y.ChenJ. C. (2013). Advances in studies on Musa basjoo Sieb.Et Zucc. Guangdong Med. Coll.29, 675677. 10.3969/j.issn.1006-8783.2013.06.024

  • 19

    MaP.GuB.XiongW.TanB.GengW.LiJ.et al (2014). Glimepiride promotes osteogenic differentiation in rat osteoblasts via the PI3K/Akt/eNOS pathway in a high glucose microenvironment. PLoS One9 (11), e112243. 10.1371/journal.pone.0112243

  • 20

    MatsubaraT.YasudaK.MizutaK.KawaueH.KokabuS. (2022). Tyrosine kinase Src is a regulatory factor of bone Homeostasis. Int. J. Mol. Sci.23 (10), 5508. 10.3390/ijms23105508

  • 21

    MolnarV.MatišićV.KodvanjI.BjelicaR.JelečŽ.HudetzD.et al (2021). Cytokines and Chemokines involved in osteoarthritis pathogenesis. Int. J. Mol. Sci.22 (17), 9208. Published 2021 Aug 26. 10.3390/ijms22179208

  • 22

    NieH.ZhengY.LiR.GuoT. B.HeD.FangL.et al (2013). Phosphorylation of FOXP3 controls regulatory T cell function and is inhibited by TNF-α in rheumatoid arthritis. Nat. Med.19 (3), 322328. 10.1038/nm.3085

  • 23

    OuyangZ.DongL.YaoF.WangK.ChenY.LiS.et al (2023). Cartilage-related Collagens in osteoarthritis and rheumatoid arthritis: from pathogenesis to therapeutics. Int. J. Mol. Sci.24 (12), 9841. 10.3390/ijms24129841

  • 24

    PanovaE.JonesG. (2015). Benefit-risk assessment of diacerein in the treatment of osteoarthritis. Drug Saf.38 (3), 245252. 10.1007/s40264-015-0266-z

  • 25

    PoznyakA.GrechkoA. V.PoggioP.MyasoedovaV. A.AlfieriV.OrekhovA. N. (2020). The Diabetes Mellitus-atherosclerosis connection: the role of lipid and glucose metabolism and chronic inflammation. Int. J. Mol. Sci.21 (5), 1835. 10.3390/ijms21051835

  • 26

    RasulA.MillimounoF. M.Ali EltaybW.AliM.LiJ.LiX. (2013). Pinocembrin: a novel natural compound with versatile pharmacological and biological activities. Biomed. Res. Int.2013, 379850. 10.1155/2013/379850

  • 27

    Rose-JohnS. (2012). IL-6 trans-signaling via the soluble IL-6 receptor: importance for the pro-inflammatory activities of IL-6. Int. J. Biol. Sci.8 (9), 12371247. 10.7150/ijbs.4989

  • 28

    SakthivelK. M.VishnupriyaS.Priya DharshiniL. C.RasmiR. R.RameshB. (2022). Modulation of multiple cellular signalling pathways as targets for anti-inflammatory and anti-tumorigenesis action of Scopoletin. J. Pharm. Pharmacol.74 (2), 147161. 10.1093/jpp/rgab047

  • 29

    SenN.HaraM. R.AhmadA. S.CascioM. B.KamiyaA.EhmsenJ. T.et al (2009). GOSPEL: a neuroprotective protein that binds to GAPDH upon S-nitrosylation. Neuron63 (1), 8191. 10.1016/j.neuron.2009.05.024

  • 30

    SiroverM. A. (2021). The role of posttranslational modification in moonlighting glyceraldehyde-3-phosphate dehydrogenase structure and function. Amino Acids53 (4), 507515. 10.1007/s00726-021-02959-z

  • 31

    TsaiH. C.VelichkoS.HungL. Y.WuR. (2013). IL-17A and Th17 cells in lung inflammation: an update on the role of Th17 cell differentiation and IL-17R signaling in host defense against infection. Clin. Dev. Immunol.2013, 267971. 10.1155/2013/267971

  • 32

    WangX. P.SunY. C.YanF. Y.LouF. Y.WuH. M.XuS. N. (2010). The constituents analysis of amino acids in rhizoma Musae. Lishizhen Med. Mater Med. Res.21, 2428. 10.3969/j.issn.1008-0805.2010.10.003

  • 33

    WeiJ. F.ZhangQ.ZhaoL.KangW. Y. (2010). Antimicrobial activity of Musa basjoo in vitro. Chin. J. Exp. Tradit. Med. Form.16, 6971. 10.13422/j.cnki.syfjx.2010.17.014

  • 34

    WuT.YehJ. H.WuP.ChenK. T.LinL. C.ChenC. (1995). 7-Hydroxyrutaecarpine from tetradium glabrifolium and tetradium ruticarpum. Heterocycles41, 10711076. 10.3987/com-95-7050

  • 35

    XuF.WuH.WangX.YangY.WangY.QianH.et al (2014). RP-HPLC Characterization of lupenone and β-sitosterol in rhizoma Musae and evaluation of the anti-Diabetic activity of lupenone in Diabetic Sprague-Dawley Rats. Molecules19 (9), 1411414127. 10.3390/molecules190914114

  • 36

    YangD.GaoX. (2021). Research progress on the antioxidant biological activity of beer and strategy for applications. Trends Food Sci. Technol.110 (0), 754764. 10.1016/j.tifs.2021.02.048

  • 37

    YaoQ.WuX.TaoC.GongW.ChenM.QuM.et al (2023). Osteoarthritis: pathogenic signaling pathways and therapeutic targets. Signal Transduct. Target Ther.8 (1), 56. 10.1038/s41392-023-01330-w

  • 38

    ZhaoL.ZhangH.LiN.ChenJ.XuH.WangY.et al (2023). Network pharmacology, a promising approach to reveal the pharmacology mechanism of Chinese medicine formula. J. Ethnopharmacol.309, 116306. 10.1016/j.jep.2023.116306

  • 39

    ZhengW.GuX.HuD.HaoY. (2019). Co-culture with synovial tissue in patients with rheumatoid arthritis suppress cell proliferation by regulating MAPK pathway in osteoblasts. Am. J. Transl. Res.11 (6), 33173327.

  • 40

    ZhuL. G.YuJ.WangS. Q. (2022). Expert Consensus on the clinical Application of Gukang capsule in treating osteoporosis. Chin. Archives Traditional Chin. Med.40 (01), 252258. 10.13193/j.issn.1673-7717.2022.01.059

Summary

Keywords

Rhizoma Musae decoction, UHPLC-Q-Exactive-MS/MS, network analysis, osteoarthritis, signaling pathways

Citation

Zhang J, Shen W, Liu F, He H, Han S and Luo L (2025) Integrated approach with UHPLC-Q-Exactive-tandem mass spectrometry, network analysis, and molecular docking to determine potential active compounds and mechanisms of Rhizoma Musae decoction in osteoarthritis treatment. Front. Pharmacol. 15:1380335. doi: 10.3389/fphar.2024.1380335

Received

01 February 2024

Accepted

03 September 2024

Published

02 January 2025

Volume

15 - 2024

Edited by

Daniela Rigano, University of Naples Federico II, Italy

Reviewed by

Karl Tsim, Hong Kong University of Science and Technology, Hong Kong SAR, China

Cheorl-Ho Kim, Sungkyunkwan University, Republic of Korea

Updates

Copyright

*Correspondence: Jian Zhang,

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics