The initiation of adaptive immune responses relies on the precise recognition and interpretation of antigenic information. In this process, the specific binding of T cell receptors (TCRs) to peptide-major histocompatibility complex (pMHC) molecules represents one of the key molecular events in the initiation of adaptive immune responses. Accordingly, the structural features of TCR-pMHC complexes provide a fundamental basis for dissecting antigen recognition mechanisms and support rational vaccine design, therapeutic target discovery in TCR-based immunotherapy, and TCR identification and optimization. However, experimental determination of TCR-pMHC structures remains costly, time-consuming, and limited in coverage, making computational approaches essential for rapidly obtaining reliable structural information. Computational methods for predicting the structures of TCR-pMHC complexes have advanced rapidly in recent years, driven by progress in deep learning-based modeling frameworks and the increasing availability of structural and sequence resources. Despite these developments, most existing tools do not adequately distinguish the key structural and biophysical differences between MHC class I (MHC-I) and MHC class II (MHC-II) complexes during model construction. As a consequence, their predictive performance differs substantially between class I and class II complexes. In general, structural predictions for class I complexes outperform those for class II complexes. This discrepancy may be related to several fundamental differences between the two systems, including the architecture of the peptide-binding groove, the distribution of peptide lengths, and the properties of peptide flanking residues (PFRs). Compared with MHC-I molecules, MHC-II molecules usually bind longer antigenic peptides, which typically range from 13 to 25 amino acids in length. PFRs at both termini of these peptides participate in regulating the overall conformation of TCR-pMHC class II complexes and exert a pronounced effect on the geometric and physicochemical characteristics of the TCR-pMHC binding interface. Furthermore, within the TCR recognition interface, the complementarity-determining regions (CDRs) consist of segments that differ markedly in conformational behavior. They commonly include regions that are relatively rigid and structurally stable, together with highly flexible segments exhibiting substantial conformational plasticity. These rigidity-flexibility features constitute an essential structural basis enabling TCRs to recognize diverse peptide-MHC ligands and to accommodate conformational heterogeneity at the interface. However, many current modeling tools, in an effort to enforce global conformational stability or reduce structural noise, tend to over-constrain intrinsically flexible regions. Such oversimplification may lead to inappropriate rigidification of flexible CDR loops, resulting in local structural distortions, compromised interface geometry, or even complete modeling failure for specific complexes. Against this background, the review approaches the field from the perspective of computational differences between MHC-I and MHC-II complexes. We first systematically organize and summarize available resources related to TCRs and pMHCs, including structural datasets, sequence databases, prediction tools, and benchmarking studies. We then focus on five representative tools capable of predicting both class I and class II complexes—AlphaFold2, AlphaFold3, TCRmodel2, tFold-TCR, and TCR-pHLA_ModellerS. After excluding structures present in the training sets of these tools, we constructed a benchmark dataset comprising 25 class I and 10 class II TCR-pMHC complexes in the bound state and conducted a systematic evaluation using this dataset. We first employ widely used general evaluation metrics, including All-Atom Root Mean Square Deviation (All-Atom RMSD), Backbone RMSD, Template Modeling score (TM-score), and DockQ, to assess the global conformational accuracy and interface modeling quality of class I and class II complexes. For class II complexes, we propose for the first time a peptide flanking residue deviation index, including the PFRs-Deviation Index (PFRs-DI), N-PFR-Deviation Index (N-PFR-DI), and C-PFR-Deviation Index (C-PFR-DI), to quantitatively characterize conformational deviations in PFRs. In addition, we propose the CDR conformational consistency index (CCC) designed to qualitatively evaluate the ability of prediction tools to capture TCR CDR conformational flexibility. These metrics collectively assess a tool’s ability to model both overall conformation and critical functional regions, thereby addressing the limitations of existing evaluation criteria that overemphasize global structure while inadequately capturing modeling quality in key functional areas. This establishes a unified analytical framework for MHC-I and MHC-II complexes to guide data resource selection, modeling strategy formulation, and evaluation system development. The framework further advances computational modeling and provides crucial support for multi-scale analysis of TCR-pMHC recognition mechanisms and their biological functions.

The initiation of adaptive immune responses relies on the precise recognition and interpretation of antigenic information. In this process, the specific binding of T cell receptors (TCRs) to peptide-major histocompatibility complex (pMHC) molecules represents one of the key molecular events in the initiation of adaptive immune responses. Accordingly, the structural features of TCR-pMHC complexes provide a fundamental basis for dissecting antigen recognition mechanisms and support rational vaccine design, therapeutic target discovery in TCR-based immunotherapy, and TCR identification and optimization. However, experimental determination of TCR-pMHC structures remains costly, time-consuming, and limited in coverage, making computational approaches essential for rapidly obtaining reliable structural information. Computational methods for predicting the structures of TCR-pMHC complexes have advanced rapidly in recent years, driven by progress in deep learning-based modeling frameworks and the increasing availability of structural and sequence resources. Despite these developments, most existing tools do not adequately distinguish the key structural and biophysical differences between MHC class I (MHC-I) and MHC class II (MHC-II) complexes during model construction. As a consequence, their predictive performance differs substantially between class I and class II complexes. In general, structural predictions for class I complexes outperform those for class II complexes. This discrepancy may be related to several fundamental differences between the two systems, including the architecture of the peptide-binding groove, the distribution of peptide lengths, and the properties of peptide flanking residues (PFRs). Compared with MHC-I molecules, MHC-II molecules usually bind longer antigenic peptides, which typically range from 13 to 25 amino acids in length. PFRs at both termini of these peptides participate in regulating the overall conformation of TCR-pMHC class II complexes and exert a pronounced effect on the geometric and physicochemical characteristics of the TCR-pMHC binding interface. Furthermore, within the TCR recognition interface, the complementarity-determining regions (CDRs) consist of segments that differ markedly in conformational behavior. They commonly include regions that are relatively rigid and structurally stable, together with highly flexible segments exhibiting substantial conformational plasticity. These rigidity-flexibility features constitute an essential structural basis enabling TCRs to recognize diverse peptide-MHC ligands and to accommodate conformational heterogeneity at the interface. However, many current modeling tools, in an effort to enforce global conformational stability or reduce structural noise, tend to over-constrain intrinsically flexible regions. Such oversimplification may lead to inappropriate rigidification of flexible CDR loops, resulting in local structural distortions, compromised interface geometry, or even complete modeling failure for specific complexes. Against this background, the review approaches the field from the perspective of computational differences between MHC-I and MHC-II complexes. We first systematically organize and summarize available resources related to TCRs and pMHCs, including structural datasets, sequence databases, prediction tools, and benchmarking studies. We then focus on five representative tools capable of predicting both class I and class II complexes—AlphaFold2, AlphaFold3, TCRmodel2, tFold-TCR, and TCR-pHLA_ModellerS. After excluding structures present in the training sets of these tools, we constructed a benchmark dataset comprising 25 class I and 10 class II TCR-pMHC complexes in the bound state and conducted a systematic evaluation using this dataset. We first employ widely used general evaluation metrics, including All-Atom Root Mean Square Deviation (All-Atom RMSD), Backbone RMSD, Template Modeling score (TM-score), and DockQ, to assess the global conformational accuracy and interface modeling quality of class I and class II complexes. For class II complexes, we propose for the first time a peptide flanking residue deviation index, including the PFRs-Deviation Index (PFRs-DI), N-PFR-Deviation Index (N-PFR-DI), and C-PFR-Deviation Index (C-PFR-DI), to quantitatively characterize conformational deviations in PFRs. In addition, we propose the CDR conformational consistency index (CCC) designed to qualitatively evaluate the ability of prediction tools to capture TCR CDR conformational flexibility. These metrics collectively assess a tool’s ability to model both overall conformation and critical functional regions, thereby addressing the limitations of existing evaluation criteria that overemphasize global structure while inadequately capturing modeling quality in key functional areas. This establishes a unified analytical framework for MHC-I and MHC-II complexes to guide data resource selection, modeling strategy formulation, and evaluation system development. The framework further advances computational modeling and provides crucial support for multi-scale analysis of TCR-pMHC recognition mechanisms and their biological functions.

Objective This study reports on the obligatory parasitism of water mites Arrenurus sp. on Anopheles messeae at the Manzhouli Port, Inner Mongolia, China. Methods Duing July 2024, a survey on the mosquito diversity was conducted at the Manzhouli Port. Captured mosquitoes and their ectoparasites were identified to species level. Results A total of 1840 adult mosquitoes were collected, representing species from three genera: Culex(Cx. modestus, Cx. pipiens pallens), Aedes(Ae. dorsalis, Ae. flavidorsalis, Ae. flavescens), and Anopheles (An. messeae). Among all the mosqutioes specimens,3 out of 150 captured An. messeae were found to carry ectoparasitic mites, with number of 2,4,27 mites separately. Morphological and molecular identification reached the same result as water mites(Hydrachnidiae, Hydracrina). COI gene sequence showed 94% similarity with the closest species Arrenurus truncatellus. Conlusions Literature review suggests water mites are host-specific parasitism of mosquito species and herein with the first record of Arrenurus sp. parasiting on An. Messeae in the most high-latitude region globally.

ObjectiveTo establish a model that can quickly identify the aroma components in different parts of Nardostachys jatamansi， so as to provide a quality control basis for the market circulation and clinical use of N. jatamansi. MethodsHeadspace solid-phase microextraction-gas chromatography-mass spectrometry（HS-SPME-GC-MS） combined with Smart aroma database and National Institute of Standards and Technology（NIST） database were used to characterize the aroma components in different parts of N. jatamansi， and the aroma components were quantified according to relative response factor（RRF） and three internal standards， and the markers of aroma differences in different parts of N. jatamansi were identified by orthogonal partial least squares-discriminant analysis（OPLS-DA） and cluster thermal analysis based on variable importance in the projection（VIP） value >1 and P<0.01. The odor data of different parts of N. jatamansi were collected by Heracles Ⅱ Neo ultra-fast gas phase electronic nose， and the correlation between compound types of aroma components collected by the ultra-fast gas phase electronic nose and the detection results of HS-SPME-GC-MS was investigated by drawing odor fingerprints and odor response radargrams. Chromatographic peak information with distinguishing ability≥0.700 and peak area≥200 was selected as sensor data， and the rapid identification model of different parts of N. jatamansi was established by principal component analysis（PCA）， discriminant factor alysis（DFA）， soft independent modeling of class analogies（SIMCA） and statistical quality control analysis（SQCA）. ResultsThe HS-SPME-GC-MS results showed that there were 28 common components in the underground and aboveground parts of N. jatamansi， of which 22 could be quantified and 12 significantly different components were screened out. Among these 12 components， the contents of five components（ethyl isovalerate， 2-pentylfuran， benzyl alcohol， nonanal and glacial acetic acid，） in the aboveground part of N. jatamansi were significantly higher than those in the underground part（P<0.01）， the contents of β-ionone， patchouli alcohol， α-caryophyllene， linalyl butyrate， valencene， 1，8-cineole and p-cymene in the underground part of N. jatamansi were significantly higher than those in the aboveground part（P<0.01）. Heracles Ⅱ Neo electronic nose results showed that the PCA discrimination index of the underground and aboveground parts of N. jatamansi was 82， and the contribution rates of the principal component factors were 99.94% and 99.89% when 2 and 3 principal components were extracted， respectively. The contribution rate of the discriminant factor 1 of the DFA model constructed on the basis of PCA was 100%， the validation score of the SIMCA model for discrimination of the two parts was 99， and SQCA could clearly distinguish different parts of N. jatamansi. ConclusionHS-SPME-GC-MS can clarify the differential markers of underground and aboveground parts of N. jatamansi. The four analytical models provided by Heracles Ⅱ Neo electronic nose（PCA， DFA， SIMCA and SQCA） can realize the rapid identification of different parts of N. jatamansi. Combining the two results， it is speculated that terpenes and carboxylic acids may be the main factors contributing to the difference in aroma between the underground and aboveground parts of N. jatamansi.

In the central nervous system (CNS), the myelin sheath, a specialized membrane structure that wraps around axons, is formed by oligodendrocytes through a highly coordinated spatiotemporal developmental program. The process begins with the directed differentiation of neural precursor cells into oligodendrocyte precursor cells (OPCs), followed by their migration, proliferation, differentiation, and maturation, ultimately leading to the formation of a multi-segmental myelin sheath structure. Recent single-cell sequencing research has revealed that this process involves the temporal regulation of over 200 key genes, with a regulatory network composed of transcription factors such as Sox10 and Olig2 playing a central role. The primary function of the myelin sheath is to accelerate nerve signal transmission and protect nerve fibers from damage. Its insulating properties not only increase nerve conduction speed by 50-100 times but also ensure the long-term functional integrity of the nervous system by maintaining axonal metabolic homeostasis and providing mechanical protection. The pathological effects of myelin sheath injury exhibit a cascade amplification pattern: acute demyelination leads to action potential conduction block, while chronic lesions may cause axonal damage and neuronal death in severe or long-term cases, ultimately resulting in irreversible neurological dysfunction with neurodegenerative characteristics. Multiple sclerosis (MS) is a neurodegenerative disease characterized by chronic inflammatory demyelination of the CNS. Clinically, the distribution of lesions in MS exhibits spatial heterogeneity, which is closely related to differences in the regenerative capacity of oligodendrocytes within the local microenvironment. Emerging evidence suggests that astrocytes form a dynamic “neural-immune-metabolic interface” and play a multidimensional regulatory role in myelin development and regeneration by forming heterogeneous populations composed of different subtypes. During embryonic development, astrocytes induce the targeted differentiation of OPCs in the ventricular region through the Wnt/β-catenin pathway. In the mature stage, they secrete platelet-derived growth factor AA (PDGF-AA) to establish a chemical gradient that guides the precise migration of OPCs along axonal bundles. Notably, astrocytes also provide crucial metabolic support by supplying energy substrates for high-energy myelin formation through the lactate shuttle mechanism. In addition, astrocytes play a dual role in myelin regulation. During the acute injury phase, reactive astrocytes establish a triple defense system within 72 h: upregulating glial fibrillary acidic protein (GFAP) to form scars that isolate lesions, activating the JAK-STAT3 regeneration pathway in oligodendrocytes via leukemia inhibitory factor (LIF), and releasing tumor necrosis factor-stimulated gene-6 (TSG-6) to inhibit excessive microglial activation. However, in chronic neurodegenerative diseases, the phenotypic transformation of astrocytes contributes to microenvironmental deterioration. The secretion of chondroitin sulfate proteoglycans (CSPGs) inhibits OPC migration via the RhoA/ROCK pathway, while the persistent release of reactive oxygen species (ROS) leads to mitochondrial dysfunction and the upregulation of complement C3-mediated synaptic pruning. This article reviews the mechanisms by which astrocytes regulate the development and regeneration of myelin sheaths in the CNS, with a focus on analyzing the multifaceted roles of astrocytes in this process. It emphasizes that astrocytes serve as central hubs in maintaining myelin homeostasis by establishing a metabolic microenvironment and signaling network, aiming to provide new therapeutic strategies for neurodegenerative diseases such as multiple sclerosis.

Objective To establish a stable,reliable,and clinically relevant porcine model of endotoxin-induced acute respiratory distress syndrome(ARDS).Methods Ten 8-month-old male Bama minipigs were deeply sedated,followed by invasive mechanical ventilation and electrocardiographic monitoring.Lipopolysaccharide(LPS)was intravenously pumped at 600 μg/(kg·h)for 3 hours,then maintained at 15 μg/(kg·h)thereafter.Dynamic monitoring was performed at five time points after LPS injection(LPS 0,1,3,5,and 8 h),including arterial blood gas analysis and chest computed tomography(CT)scans.Pathological examination of lung tissues obtained via bronchoscopic biopsy(HE staining and transmission electron microscopy)was conducted.These indicators were comprehensively used to evaluate the success of the animal model.Results At 5 hours after LPS administration,8 minipigs developed symptoms such as skin cyanosis,elevated body temperature,and respiratory distress.The oxygenation index decreased to<300 mmHg.Chest CT scans showed diffuse pulmonary infiltrates.Histopathology revealed alveolar edema and hyaline membrane formation.Transmission electron microscopy demonstrated disruption of pulmonary blood-air barrier,depletion of lamellar bodies in type Ⅱ pneumocytes,inflammatory cell infiltration,and exudation of plasma proteins and fibrin.Compared with LPS 0 h,at LPS 8 h,the oxygenation index and arterial blood pH were significantly decreased(P<0.001),while blood lactic acid and serum potassium were significantly increased(P<0.05);serum calcium and base excess were significantly decreased(P<0.05),and the lung injury score based on HE-stained lung sections was significantly increased(P<0.01).Conclusion The porcine ARDS model established by continuous LPS injection can dynamically simulate the pathophysiological characteristics and typical pathological manifestations of clinical septic ARDS,making it an effective tool to study the pathogenesis,prevention,and treatment strategies of septic ARDS.

Objective:To apply the Timed Up and Go Test(TUGT)to investigate the correlation between motor function, emotional state, cognitive function, and sleep quality among elderly individuals in the Chinese community.Methods:A cross-sectional study was conducted, involving 739 subjects aged 60 to 90 years, who were randomly recruited from December 2021 to August 2023 across Beijing, Tianjin, Zhejiang, Guangdong, and Hainan Provinces in China.Basic demographic information was collected, and the TUGT was utilized to assess motor function.Based on the TUGT time(t), the subjects were divided into three groups: normal motor function group, mild motor abnormality group, and significant motor abnormality group.Cognitive function was evaluated using the Chinese Revised Mini-Mental State Examination(MMSE), while the Patient Health Questionnaire Depression Scale(PHQ-9)was employed to measure the degree of depression.Additionally, the Epworth Sleepiness Scale(ESS)was used to assess excessive daytime sleepiness.The correlation between subjects' motor function and their cognitive abilities, mood, and sleep was subsequently analyzed.Results:Systolic blood pressure, heart rate, PHQ-9, MMSE, and ESS scores were identified as significant factors influencing TUGT time.Specifically, TUGT time was positively correlated with PHQ-9 and ESS scores, while exhibiting negative correlations with systolic blood pressure, heart rate, and MMSE scores.Additionally, TUGT time was negatively correlated with the MMSE subcomponents of orientation, immediate memory, and verbal ability.All observed differences were statistically significant(all P<0.05).Logistic regression analysis indicated that an increase in the PHQ-9 score was associated with an odds ratio( OR)of 1.099(95% CI: 1.045-1.155, P<0.001)(mild motor abnormality group)and 1.150(95% CI: 1.066-1.242, P<0.001)(Significant motor abnormality group).Additionally, a reduction in the MMSE score was observed, with an OR of 0.939(95% CI: 0.886-0.995, P<0.001)(mild motor abnormality group)and 0.793(95% CI: 0.729-0.862, P<0.001)(Significant motor abnormality group).Furthermore, an increase in the ESS score was noted, with ORs of 1.139(95% CI: 1.094-1.186, P<0.001)(mild motor abnormality group)and 1.203(95% CI: 1.132-1.279, P<0.001)(Significant motor abnormality group).These findings suggest that these variables are independently related to decreased motor function. Conclusions:Depression, cognitive impairment, and excessive daytime sleepiness are independent risk factors for motor dysfunction among elderly individuals in community settings.The Timed Up and Go Test TUGT can be utilized for the early screening of motor function decline in this population.

AIM To investigate the variations in mercury dissolution from HgS-containing traditional medicines in three kinds of simulated gastrointestinal dissolution media.METHODS 39 batches of 15 types of HgS-containing traditional medicines were collected,total mercury content and dissolved mercury concentrations in simulated gastric fluid,simulated intestinal fluid,and L-cysteine-containing simulated intestinal fluid were measured.The maximum daily intake of total mercury and soluble mercury was calculated based on the maximum daily clinical dosage.RESULTS Among the 15 types of medicines,the maximum daily intake of total mercury varied by 156 times,the daily intake of soluble mercury varied by 3 502 times in simulated gastric fluid,313 times in simulated intestinal fluid,and 10 663 times in L-cysteine-containing simulated intestinal fluid,approximately.CONCLUSION For the 15 types of HgS-containing traditional medicines,the daily maximum intake of soluble mercury showed greater variations than that of total mercury.Soluble mercury concentration is more closely correlated with intestinal absorption of mercury and thus represents a more rational quality control indicator for HgS-containing traditional medicines.

The study aims to examine the effects and potential mechanisms of disulfiram (DSF) on cardiac hypertrophic injury, focusing on the role of transforming growth factor-β-activated kinase 1 (TAK1)-mediated pan-apoptosis (PANoptosis). H9C2 cardiomyocytes were treated with angiotensin II (Ang II, 1 µmol/L) to establish an in vitro model of myocardial hypertrophy. DSF (40 µmol/L) was used to treat cardiomyocyte hypertrophic injury models, either along or in combination with the TAK1 inhibitor, 5z-7-oxozeaenol (5z-7, 0.1 µmol/L). We assessed cell damage using propidium iodide (PI) staining, measured cell viability with CCK8 assay, quantified inflammatory factor levels in cell culture media via ELISA, detected TAK1 and RIPK1 binding rates using immunoprecipitation, and analyzed the protein expression levels of key proteins in the TAK1-mediated PANoptosis pathway using Western blot. In addition, the surface area of cardiomyocytes was measured with Phalloidin staining. The results showed that Ang II significantly reduced the cellular viability of H9C2 cardiomyocytes and the binding rate of TAK1 and RIPK1, significantly increased the surface area of H9C2 cardiomyocytes, PI staining positive rate, levels of inflammatory factors [interleukin-1β (IL-1β), IL-18, and tumor necrosis factor α (TNF-α)] in cell culture media and p-TAK1/TAK1 ratio, and significantly up-regulated key proteins in the PANoptosis pathway [pyroptosis-related proteins NLRP3, Caspase-1 (p20), and GSDMD-N (p30), apoptosis-related proteins Caspase-3 (p17), Caspase-7 (p20), and Caspase-8 (p18), as well as necroptosis-related proteins p-MLKL, RIPK1, and RIPK3]. DSF significantly reversed the above changes induced by Ang II. Both 5z-7 and exogenous IL-1β weakened these cardioprotective effects of DSF. These results suggest that DSF may alleviate cardiac hypertrophic injury by inhibiting TAK1-mediated PANoptosis.
Animals ; MAP Kinase Kinase Kinases/physiology* ; Rats ; Myocytes, Cardiac/pathology* ; Disulfiram/pharmacology* ; Cardiomegaly ; Apoptosis/drug effects* ; Cell Line ; Angiotensin II ; Necroptosis/drug effects* ; Interleukin-1beta/metabolism* ; Receptor-Interacting Protein Serine-Threonine Kinases/metabolism* ; Lactones ; Resorcinols ; Zearalenone/administration & dosage*

This article aims to explore the effect and mechanism of Liujunzi Pills on the intestinal microbiota of rats with spleen Qi deficiency syndrome. The raw Rhei Radix et Rhizoma water extract(1 g·mL~(-1)) was used to prepare spleen Qi deficiency rat models. A total of 44 SD male rats were randomly divided into a control group, a model group, Liujunzi Pills groups at high(3.24 g·kg~(-1)), medium(1.62 g·kg~(-1)), low(0.81 g·kg~(-1)) doses, and Shenling Baizhu San(2.50 g·kg~(-1)) group. The drug effect was evaluated by observing the following aspects: spleen index, fecal water content, body weight, and intestinal propulsion index. Gut microbiota analysis and 16S rRNA gene sequencing were conducted on feces. Enzyme-linked immunosorbent assay(ELISA) and UV spectrophotometry were used to detect interleukin-1β(IL-1β) and adenosine triphosphate(ATP) levels in small intestine tissues. Hematoxylin-eosin staining and transmission electron microscopy were employed to observe changes in intestinal pathology and microstructure. The results show that, compared with the control group, fecal moisture content is significantly increased while spleen index, body weight, and intestinal propulsion index are significantly reduced in rats of the model group, indicating the successful establishment of the model. The above symptoms can be improved by both Shenling Baizhu San and Liujunzi Pills. Compared with the control group, in the model group, the gut microbiota abundance is changed with an unbalanced development: the abundance of beneficial bacteria within the Bacteroidetes phylum is reduced, accompanied by a significantly decreased Shannon index, and reduced signal levels of nicotinamide adenine dinucleotide phosphate(NADPH)-related enzymes relevant to mitochondria. However, Liujunzi Pills and Shenling Baizhu San can significantly improve the Bacteroidetes phylum abundance in gut microbiota, microbial diversity, and NADPH activity in the model group. Additionally, compared with the control group, the ATP level is decreased and the IL-1β level is increased in small intestinal tissues of the model group, with shorter small intestinal epithelial villi and decreased mitochondrial number. The above symptoms can be improved by Liujunzi Pills and Shenling Baizhu San. In conclusion, Liujunzi Pills can treat spleen Qi deficiency syndrome by enhancing mitochondrial function to regulate gut microbiota balance and diversity.
Animals ; Gastrointestinal Microbiome/drug effects* ; Drugs, Chinese Herbal/pharmacology* ; Male ; Rats, Sprague-Dawley ; Rats ; Qi ; Spleen/metabolism* ; Splenic Diseases/metabolism* ; Humans ; Interleukin-1beta/genetics* ; Bacteria/drug effects* ; Feces/microbiology* ; Adenosine Triphosphate/metabolism*