All materials, including databases, software, and experimental reagents, used in this study are listed in the Table of Materials. All experimental procedures were approved by the Animal Ethics Committee of the Suzhou Dushu Lake Hospital (permit number: 2410008).
Screening of active components and target acquisition of DCQD
The active components of the constituent herbs rhubarb, mirabilite, magnolia bark, and immature bitter orange in Da-Cheng-Qi Decoction (DCQD) were screened from the traditional Chinese medicine systems pharmacology (TCMSP) database using the following criteria: oral bioavailability (OB) ≥30% and drug-likeness (DL) ≥0.18. The corresponding targets of the identified active components were then collected. These thresholds are widely adopted in network pharmacology studies to filter for compounds with favorable pharmacokinetic properties and higher potential for being drug-like molecules34,35. To improve the reliability of compound selection and reduce bias from single-database screening, cross-validation was performed using the BATMAN-TCM database. BATMAN-TCM is an integrative database that stores known and predicted interactions between TCM ingredients and target proteins, aiding in the exploration of TCM pharmacological mechanisms and drug discovery. The updated version 2.0 offers a significantly expanded TTI dataset and enhanced features, making it a valuable resource for understanding TCM's molecular mechanisms and developing new treatments for complex diseases. Compounds without predicted targets or with BATMAN-TCM target scores lower than 20 were excluded. In addition, compounds with relatively low OB values but with predicted targets in BATMAN-TCM and well-documented pharmacological activities (e.g., anthraquinones and flavonoids) were retained through literature curation. To further reduce potential false-positive components, endogenous hormones and neurotransmitter-like substances (e.g., progesterone and serotonin) were manually excluded. The final compound set was used for subsequent target prediction and network construction. The detailed workflow is illustrated in Supplementary Figure 1. The final list of active components is provided in Supplementary Table 1.
The SMILES identifiers (a linear notation method that uses ASCII strings to describe the chemical structure of molecules. By employing specific symbols and rules, it encodes atoms, chemical bonds, and molecular topology into a continuous string of characters without spaces, enabling concise, unambiguous, and computer-readable storage and transmission of molecular information) of each active component were retrieved from the PubChem database and potential targets were predicted using the swiss target prediction, with the organism parameter set to 'Homo sapiens' and a probability score cut-off of >0 applied. The target sets obtained from both databases were merged, and duplicate entries were removed to generate the final candidate target set for DCQD.
It is important to note that Mirabilitum (Mangxiao, Na₂SO₄·10H₂O), a mineral drug in DCQD, is not included in the TCMSP or similar compound-target databases. Consequently, due to this database limitation, the component-target network constructed in this study does not encompass the potential contributions of Mirabilitum. The network is therefore based on the phytochemical constituents of the other three herbal components: Rheum palmatum, Magnolia officinalis,and Citrus aurantium.
Enrichment analysis of DCQD targets
DCQD targets were imported into the metascape database for enrichment analysis, with the species restricted to "Homo sapiens". Gene ontology (GO) and Kyoto encyclopedia of genes and genomes (KEGG) enrichment analyses were performed, with the former encompassing biological processes (BP), cellular components (CC), and molecular functions (MF). The top 6 of GO enrichment results and the top 20 of KEGG enrichment results were selected. The KEGG and GO bar plots were generated using the online tool from bioinformatics. Data was imported via the DOSE package in R to construct the disease ontology (DO) bar plot.
Construction of the "Herbal Medicine-Component-Target" network
The relevant components of DCQD and their corresponding targets were imported into Cytoscape 3.10.3 to construct a network diagram illustrating the "herbal medicine-component-target" relationships. The CytoNCA plugin was employed to calculate Degree values, and the top 5 active components were screened based on ranking.
Retrieval of ARDS targets from public databases
Using "acute respiratory distress syndrome" as the keyword, ARDS-related targets were retrieved from the OMIM and GeneCards databases. In Gene Cards, targets with a relevance score >5 were selected. After merging and removing duplicates from both databases, a candidate target set for ARDS was obtained.
Bioinformatics analysis
The sepsis-ARDS-related dataset GSE32707 (GPL570 platform, n = 45), comprising RNA-seq data from peripheral blood mononuclear cells of 30 sepsis-induced ARDS patients and 15 healthy controls, was downloaded from the NCBI GEO database. Data were imported using the GEOquery package in R, and normalization was performed using the limma package to minimize inter-sample variability. Differentially expressed genes (DEGs) were identified using the thresholds |log2FC| > 1 and p < 0.05. Visualization was conducted via volcano plots and heatmaps. GO and KEGG enrichment analyses were performed for pathway analysis, with results presented as bubble plots and bar plots.
Prediction of potential targets for DCQD in ARDS treatment
A Venn diagram was generated using the bioinformatics online tool to compare the DEGs, ARDS-related targets from public databases, and DCQD drug targets, thereby identifying potential therapeutic targets of DCQD for ARDS.
Construction of the protein-protein interaction (PPI) network and screening of core targets
The overlapping targets were imported into the STRING database with the species set to "Homo sapiens" and an interaction score threshold of 0.4. The PPI network data were exported in TSV format and visualized using Cytoscape 3.10.3. The CytoHubba plugin was applied to screen the top 10 targets based on MCC, MNC, degree, EPC, betweenness, closeness, radiality, and stress scores. The intersection of these targets yielded 5 core targets.
Construction of the "Herbal Medicine-Component-Target-Disease" network
Based on the overlapping targets, corresponding DCQD components were identified. The relevant components, along with their associated targets and disease linkages, were imported into Cytoscape 3.10.3 to construct a network diagram illustrating the "herbal medicine-component-target-disease" relationships.
Enrichment analysis of core targets and construction of the "herbal medicine-component-target-disease-pathway" network
To explore the potential biological functions and key signaling pathways of DCQD in ARDS treatment, the core targets were subjected to GO and KEGG enrichment analyses in Metascape (species: "Homo sapiens"). The top 10 of GO and the top 13 of KEGG enrichment results were selected. Visualization was performed using Bioinformatics tools to generate KEGG bubble plots, GO bar plots, and Sankey bubble plots. Additionally, based on the functional enrichment results, a network diagram was constructed in Cytoscape 3.10.3 to illustrate the "herbal medicine-component-target-disease-pathway" relationships.
Molecular docking and visualization
The secondary molecular structures of compounds were retrieved from the PubChem database, while the 3D structures of core targets were obtained from the PDB database (http://www.rcsb.org/). Water molecules and amino acid residues were removed using PyMOL. Target proteins and compounds were imported into AutoDock for molecular docking to calculate binding energies. Visualization was performed using PyMOL. A binding energy heatmap was generated using the pheatmap package in R, with core targets on the x-axis, compounds on the y-axis, and a color gradient representing binding energies. Two-dimensional interaction diagrams of docking results were generated using Discovery Studio Visualizer 4.5. Molecular docking was performed to evaluate the binding feasibility between compounds and targets, rather than to predict actual inhibitory efficacy.
Construction of ARDS mouse model
This study was approved by the medical ethics committee of Suzhou Dushu Lake Hospital (2410008). Healthy male C57BL/6J mice (n = 30; 20 ± 5 g) were obtained. The mice were divided into a control group (n = 10), an LPS group (n = 10, LPS 10 mg/kg), and a DCQD group (n = 10, receiving DCQD 0.9 g/kg via oral gavage combined with LPS 10 mg/kg). Mice in the LPS group received an intraperitoneal injection of LPS. Mice in the DCQD group were administered DCQD 0.9 g/kg via oral gavage immediately following the intraperitoneal injection of LPS. Measurements were taken immediately before and after lung extraction from all mice. Subsequently, the lung tissues were stored at -80 °C.
Western blot analysis (WB)
Proteins were isolated from mouse lung tissue lysates, and their concentrations were measured using a BCA protein quantification kit. Subsequently, the proteins were denatured by boiling at 95 °C for 10 min. The denatured proteins were separated on a 10% SDS-PAGE gel prepared with a gel preparation kit, then transferred to a PVDF membrane. To prevent nonspecific binding, the PVDF membranes were blocked at room temperature for 5 min using a protein-free rapid blocking buffer. After washing with PBST, the membranes were incubated overnight at 4 °C with primary antibodies against MMP9, HSP90, MYC, and ARG1. Following another wash with PBST, the membranes were incubated for 2 h at room temperature with HRP-conjugated goat anti-rabbit IgG(H+L). Finally, the ultra-sensitive chemiluminescence detection kit substrate was applied, and chemiluminescent protein bands were detected using a BIO-RAD chemiluminescence imaging system. Grayscale analysis was performed with ImageJ. The antibodies used in this experiment are listed in Supplementary Table 2.
Enzyme-linked immunosorbent assay analysis (ELISA)
First, bronchoalveolar lavage fluid (BALF) was collected by performing lavage three times with 2 mL of saline. The resulting lavage fluids were pooled and centrifuged at 10,000 g for 10 min to remove cells. The supernatant was then aspirated and stored frozen in liquid nitrogen. Subsequently, the concentrations of target cytokines were determined following the manufacturer's instructions using the Mouse IL-6 ELISA Kit, Mouse IL-18 ELISA Kit, and Mouse TNF-alpha ELISA Kit. Next, the reaction plate was placed in a microplate reader, and the optical density (OD) value of each well was measured at a wavelength of 450 nm, with the blank control well used for zero adjustment. Finally, a standard curve was plotted based on the concentrations and corresponding OD values of the standards. The concentration of the target factor in each sample was calculated according to the standard curve equation.
Statistical analysis
Data were analyzed using GraphPad Prism 10.1.2. The significance of the differences between groups was statistically evaluated through one-way analysis of variance or Student's t-test. If the P value is less than 0.05, the difference is considered significant; if the P value is less than 0.01, it is regarded as highly significant.