Exploring medication rules and mechanism of Chinese medicine for children with cough variant asthma based on data mining, network pharmacology, and molecular docking

Cough variant asthma (CVA) is a common disease with high incidence among children. Cough is the main clinical symptom and Chinese medicine (CM) has an exact effect on CVA. However, the rules of herb formulation, the pharmacodynamic substances, and the mechanism remained unclear. Therefore, we conducted this article to explore medication rules and molecular mechanism of CM against CVA in children using data mining, network pharmacology, molecular docking, and molecular dynamics simulation. Relevant literatures were collected from China National Knowledge Infrastructure, Chinese Scientific and Technical Journals database, Wanfang database, Pubmed, and Web of science. Excel 2016 was used to extract related data and establish the database for Chinese medical frequency, properties, tastes, and meridian analysis. Association rules were analyzed based on Apriori algorithm using IBM SPSS Modeler 18.0 software and core herb combination was identified. The active ingredients and targets of the core herb combination were acquired form the Traditional Chinese Medicine Systems Pharmacology Database and Analysis Platform database. The main targets of CVA were obtained from the GeneCards and Online Mendelian Inheritance in Man database. Core targets were selected by using STRING platform and Cytoscape 3.7.2 software. Metascape platform was utilized to perform gene ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analysis. The results were verified by molecular docking. SwissADME and pkCSM website were used to analyze the pharmacokinetic profiles and toxicity of key components of the core herb combination. Molecular dynamics simulation was utilized to evaluate the stable of protein-ligand complex. Two hundred seventy-five literatures containing 202 herbs were finally collected. Statistics indicated that these herbs possessed bitter, pungent taste, and warm properties, and belonged to lung meridian. Glycyrrhizae radix et rhizome, Ephedrae herba, and Armeniacae semen amarum were the most frequently used herbs. "Glycyrrhizae radix et rhizoma-ephedrae herba-Armeniacae semen amarum" was the core herb combination with highest support and confidence. Network pharmacology predicted that the main active ingredients, like quercetin, kaempferol, luteolin, etc, might target on RAC-alpha serine/threonine-protein kinase, tumor necrosis factor, interleukin-6, vascular endothelial growth factor A, transcription factor AP-1, interleukin-1 beta, matrix metalloproteinase-9, etc. They played a pivotal role in regulating multiple signaling pathways, such as tumor necrosis factor signaling pathway, IL-17 signaling pathway, and PI3K-Akt signaling pathway. Molecular docking revealed that the key active ingredients were well docked with core targets. The absorption, distribution, metabolism, excretion, and toxicity analysis showed that formononetin, luteolin, naringenin, and quercetin have high gastrointestinal absorption, no AMES toxicity, hepatotoxicity, and skin sensitization. Molecular dynamics simulation revealed that the formononetin-matrix metalloproteinase-9 complex was relatively stable. This article revealed that CM against CVA in children focused on dispelling wind and reducing phlegm, warming lung, and relieving cough. The mechanism of the core herb combination of CM for CVA through muti-components, muti-targets, and muti-pathways.