Parkinson's disease （PD） is a complex neurodegenerative disease involving multiple systems and neurotransmitters. Due to the high clinical heterogeneity of PD，it is urgent to establish a comprehensive and long-term traditional Chinese medicine （TCM） management model. In this paper，the conceptual framework of full-cycle management of PD is preliminarily constructed：based on the evolution of the pathophysiological mechanisms of protein deposition and neurotransmitter disorder in PD，the three-stage syndrome characteristics of the prodromal stage （predominant healthy Qi with subtle pathogenic factors），the early clinical stage （declining healthy Qi with growing pathogenic factors） and the middle and late stages （overwhelming pathogenic factors with deficient healthy Qi） are longitudinally described. Through the syndrome differentiation of visceral manifestations，the etiology and pathogenesis of PD motor and non-motor symptoms were comprehensively analyzed，while the matching treatment methods and prescriptions were inferred，and the modular scheme of the combining main symptoms，accompanying symptoms and secondary symptoms was proposed. The conceptual gap of TCM regarding motor complications （'variable syndrome'） and PD-related hyperpyrexia syndrome （'critical syndrome'） was explained. This framework reflects the characteristics of combination of disease and syndrome and overall constant motion，and provides new theories and research ideas for individualized and whole-process management of PD in TCM.

Pulmonary fibrosis（PF） is a progressive and life-threatening disease with limited therapeutic options， highlighting the urgent need for innovative drug discovery strategies. To address this challenge， the authors propose the formula-originated rational intelligent screening&translation（FIRST）， a systematic framework for developing anti-fibrotic monomers derived from classical traditional Chinese medicine（TCM）. The strategy integrates three key dimensions， including tissue-oriented intelligent screening of active compounds， structural optimization based on drug-target spatial interactions and plant biosynthetic pathways， and cross-scale validation of drug. We further highlight its applications in discovering tissue-oriented novel drugs from clinically validated TCM， the development and mechanistic elucidation of anti-fibrotic therapeutics， as well as the clinical translation and secondary development of candidate drugs. This strategy paves the way for first-in-class， formula-derived monomeric drugs with defined structures， clarified mechanisms， and proven safety， offering a transformative avenue to meet the urgent therapeutic needs of PF and setting a new paradigm for TCM-based drug innovation.

Parkinson's disease （PD） is a complex neurodegenerative disease involving multiple systems and neurotransmitters. Due to the high clinical heterogeneity of PD，it is urgent to establish a comprehensive and long-term traditional Chinese medicine （TCM） management model. In this paper，the conceptual framework of full-cycle management of PD is preliminarily constructed：based on the evolution of the pathophysiological mechanisms of protein deposition and neurotransmitter disorder in PD，the three-stage syndrome characteristics of the prodromal stage （predominant healthy Qi with subtle pathogenic factors），the early clinical stage （declining healthy Qi with growing pathogenic factors） and the middle and late stages （overwhelming pathogenic factors with deficient healthy Qi） are longitudinally described. Through the syndrome differentiation of visceral manifestations，the etiology and pathogenesis of PD motor and non-motor symptoms were comprehensively analyzed，while the matching treatment methods and prescriptions were inferred，and the modular scheme of the combining main symptoms，accompanying symptoms and secondary symptoms was proposed. The conceptual gap of TCM regarding motor complications （'variable syndrome'） and PD-related hyperpyrexia syndrome （'critical syndrome'） was explained. This framework reflects the characteristics of combination of disease and syndrome and overall constant motion，and provides new theories and research ideas for individualized and whole-process management of PD in TCM.

Pulmonary fibrosis（PF） is a progressive and life-threatening disease with limited therapeutic options， highlighting the urgent need for innovative drug discovery strategies. To address this challenge， the authors propose the formula-originated rational intelligent screening&translation（FIRST）， a systematic framework for developing anti-fibrotic monomers derived from classical traditional Chinese medicine（TCM）. The strategy integrates three key dimensions， including tissue-oriented intelligent screening of active compounds， structural optimization based on drug-target spatial interactions and plant biosynthetic pathways， and cross-scale validation of drug. We further highlight its applications in discovering tissue-oriented novel drugs from clinically validated TCM， the development and mechanistic elucidation of anti-fibrotic therapeutics， as well as the clinical translation and secondary development of candidate drugs. This strategy paves the way for first-in-class， formula-derived monomeric drugs with defined structures， clarified mechanisms， and proven safety， offering a transformative avenue to meet the urgent therapeutic needs of PF and setting a new paradigm for TCM-based drug innovation.

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.

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.

ObjectiveHepatocellular carcinoma (HCC) represents 90% of all primary liver cancers. The main risk factors associated with HCC include viral hepatitis (B and/or C), alcohol abuse, and metabolic dysfunction-associated steatotic liver disease (MASLD), which progressively advance to liver fibrosis, cirrhosis, and ultimately evolve into HCC. Surgical resection represents the most effective treatment for HCC, while recent advances in immunotherapy, including immune checkpoint inhibitors and adoptive cell therapies, have provided improved treatment prospects for patients with unresectable HCC. However, the complex metabolic heterogeneity of HCC limits the therapeutic efficacy. Metabolic intermediates acyl-CoA not only provide energy and substrates for numerous biochemical reactions but also serve as donors for protein lysine acylation, a major class of post-translational modification (PTM). Therefore, a deeper understanding of the molecular mechanisms underlying protein lysine acylation and hepatocarcinogenesis is urgently needed. MethodsThe levels of protein lysine acylation and silence information regulator 5 (SIRT5) expression levels in clinical HCC samples were analyzed by Western blot. Quantitative malonylome and succinylome of HCC samples were analyzed by antibody-based affinity enrichment coupled with tandem mass spectrometry. The proliferation of HCC cells was analyzed with Cell Counting Kit-8 (CCK-8) assays, the apoptosis was quantified by Annexin V-FITC/propidium iodide (PI) staining coupled with flow cytometry, and the ability of cells to migrate was assayed by Transwell assays. The enzymatic activity of glutathione S-transferase Mu 1 (GSTM1) was quantified. Transgenic mice with hepatic overexpression of SIRT5 were constructed using CRISPR-Cas9, and primary hepatocarcinogenesis was induced by administration of diethylnitrosamine. ResultsWestern blot analysis indicated that the expression level of SIRT5 was elevated in clinical samples from HCC patients, and the levels of lysine malonylation, glutarylation, and succinylation were significantly reduced in HCC tissues. Knockout of SIRT5 in MHCC-97H and MHCC-97L hepatoma cells suppressed cell proliferation, and increased the percentage of apoptotic cells significantly. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses of the differentially malonylome and succinylome of HCC samples revealed significant enrichment in two major classes of biological processes: core energy metabolism (e.g., glycolysis/gluconeogenesis, tricarboxylic acid metabolic process, fatty acid beta oxidation) and detoxification and oxidative stress response (e.g., response to toxic substance, chemical carcinogenesis, reactive oxygen species (ROS)). SIRT5 removes malonylation from lysine residues in GSTM1 and restores its detoxification activity, which is crucial for the survival of hepatocytes under stressed conditions. More importantly, in vivo experiment indicated that hepatic-specific overexpression of SIRT5 in mice accelerated diethylnitrosamine-induced liver fibrosis and hepatocarcinogenesis, indicating the critical role of SIRT5 in HCC progression. ConclusionThis study highlights the previously unrecognized SIRT5-GSTM1 axis as a key regulator in hepatocarcinogenesis, and suggests a potential target for the treatment of patients with HCC.

Objective: To analyze the healthcare-seeking delay and its influencing factors among pulmonary tuberculosis patients in Kaihua County, Zhejiang Province from 2016 to 2025, so as to provide a basis for identifying high-risk groups and optimizing comprehensive prevention and control measures for pulmonary tuberculosis. Methods: Data on pulmonary tuberculosis patients registered in Kaihua County from 2016 to 2025 were collected through the Tuberculosis Information Management System of Chinese Disease Prevention and Control Information System. The data included demographic information, disease and treatment information, and treatment policy phases. Factors affecting delay in healthcare-seeking among pulmonary tuberculosis patients were analyzed using a multivariable logistic regression model. Results: A total of 2 069 pulmonary tuberculosis patients were registered in Kaihua County from 2016 to 2025. There were 1 480 males (71.53%) and 589 females (28.47%). The median age was 64.00 (interquartile range, 20.00) years. There were 972 patients with healthcare-seeking delay, with a healthcare-seeking delay rate of 46.98%. The rates of healthcare-seeking delay were relatively high among tuberculosis patients aged >40-60 years, those seeking treatment due to symptoms, and those in the basic policy guarantee period, at 51.96%, 56.33%, and 49.25%, respectively. The rates of healthcare-seeking delay among tuberculosis patients from 2016 to 2025 were 55.81%, 44.25%, 50.00%, 47.27%, 42.98%, 51.64%, 50.79%, 44.38%, 43.13% and 38.67%, showing a downward trend (P<0.05). Multivariable logistic regression analysis show that, compared with pulmonary tuberculosis patients seeking medical consultation due to symptoms, those identified through physical examinations/active screening (OR=0.024, 95%CI: 0.011-0.047) and referrals/tracking/other sources (OR=0.781, 95%CI: 0.639-0.954) had a lower risk of healthcare-seeking delay. Conclusions From 2016 to 2025, the overall rate of healthcare-seeking delay among pulmonary tuberculosis patients in Kaihua County showed a downward trend from 2016 to 2025, but remained at a moderate level nationwide. The source of patients was an influencing factor for healthcare-seeking delay. Health examination/active screening and referral/tracing could effectively reduce the risk of healthcare-seeking delay.

Objective:To investigate the association between serum IgG concentration and renal prognosis in patients with IgA nephropathy (IgAN).Methods:It was a multi-center retrospective cohort study, patients with biopsy proven primary IgAN who were recorded in the Chinese IgA Nephropathy Information Registration System between April 1996 and September 2018 were included. Exclusion criteria were: (1) age <18 years; (2) <8 glomeruli in biopsy specimens; (3) estimated glomerular filtration rate (eGFR) <15 ml·min -1·(1.73 m 2) -1 at biopsy; (4) missing baseline serum IgG values; (5) incomplete follow-up data; (6) follow-up duration <12 months. Enrolled patients were divided into 3 groups according to the baseline tertiles of serum IgG: ≤9.50 g/L (G1 group), 9.51-11.99 g/L (G2 group), and ≥12.00 g/L (G3 group). Clinical, and pathological parameters were compared across groups. The endpoint events were defined as doubled serum creatinine level from baseline, or end-stage renal disease (ESRD). Results:A total of 1 976 IgAN patients were included in this study, 631 were in G1 group, 664 in G2 group, and 681 in G3 group. The comparison of baseline clinical data showed that there were statistically significant differences among the three groups in terms of gender, age, microscopic hematuria, edema, body mass index, systolic blood pressure, diastolic blood pressure, hemoglobin, serum creatinine, eGFR, 24-hour urine protein quantity, blood uric acid, blood albumin, serum IgA, serum IgM, the proportion of using immunosuppressants, and the proportion of using glucocorticoids (all P<0.05). In terms of pathology, the higher the serum IgG concentration, the relatively less severe the overall renal pathological damage. The results of univariate Cox regression analysis showed that gender, systolic blood pressure, diastolic blood pressure, hemoglobin, serum creatinine, eGFR, 24-hour urine protein quantity, total protein, serum albumin, globulin, serum IgG, Oxford renal pathological classification, glomerular sclerosis ratio, and glomerular IgM deposition were all associated with the occurrence of renal endpoint events (all P<0.05). Based on clinical practice and previous studies, after adjusting for gender, age, systolic blood pressure, diastolic blood pressure, eGFR, 24-hour urine protein quantity, body mass index, Oxford renal pathological classification, glomerular sclerosis ratio, and the use of renin-angiotensin-aldosterone system inhibitors, glucocorticoids, and immunosuppressants, multivariate Cox regression analysis showed that as a continuous variable, the baseline serum IgG level ( HR=0.91, 95% CI 0.87-0.96) was independently associated with the risk of renal endpoint events in IgAN patients; as a categorical variable, with serum IgG ≤ 9.50 g/L as the reference, serum IgG 9.51-11.99 g/L and serum IgG ≥ 12.00 g/L were independent factors for the occurrence of renal endpoint events in IgAN patients ( HR=0.69, 95% CI 0.49-0.96, P=0.027; HR=0.50, 95% CI 0.34-0.74, P<0.001). During a median follow-up of 33(21, 53) months started from the date of renal biopsy and continued until December 31, 2019, the median follow-up duration was 33 (21, 53) months, and a total of 232 patients (11.74%) reached the composite endpoint. Kaplan-Meier survival analysis showed that the higher the serum IgG concentration in patients with IgAN, the higher their cumulative renal survival rate (Log-rank test, χ2=47.176, P<0.001). Conclusion:The higher level of serum IgG at diagnosis is associated with better clinicopathologic features and renal outcomes, and may portend better renal survival in IgAN patients.