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.

IgA nephropathy (IgAN) is the most common primary glomerulonephritis worldwide and a major cause of chronic kidney disease and kidney failure. IgAN exhibits marked heterogeneity in clinical presentation, histopathology, and pathogenic mechanisms, contributing to variable treatment responses and prognosisamong patients. Precise risk assessment and individualized intervention are therefore of critical importance. This review systematically traces the evolution of IgAN management from traditional risk stratification toward biomarker-driven precision medicine. We first review the clinical utility and limitations of established risk stratification tools, including the KDIGO guidelines, the Oxford MEST-C classification, and the International IgAN Prediction Tool. We then discuss emerging biomarkers closely linked to disease pathogenesis, including galactose-deficient IgA1 (Gd-IgA1), anti-Gd-IgA1 autoantibodies, B cell activating factor (BAFF), a proliferation-inducing ligand (APRIL), and complement components, as well as the targeted therapies they have informed. In addition, urinary biomarkers and multi-omics approaches show promise for dynamic disease monitoring and individualized risk stratification.

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.

ObjectiveType I interferon (IFN-I) signaling is essential for antiviral innate immunity, yet its sustained or excessive activation contributes to the pathogenesis of several autoimmune diseases and interferonopathies, such as systemic lupus erythematosus and Aicardi-Goutières syndrome. Current strategies targeting this pathway, exemplified by JAK inhibitors, act mainly on downstream signal transduction and provide limited direct control over upstream IFN-I production, while also carrying the risk of broad immunosuppression. Phyllanthus emblica L. has long been used in traditional medicine for inflammatory disorders, but the bioactive constituent responsible for its regulation of IFN-I signaling and the underlying molecular mechanism have not been clearly defined. This study aimed to identify the active anti-inflammatory component of P. emblica and to characterize its mechanism of action on the IFN-I pathway in macrophages. MethodsActive components ofP. emblica and their candidate targets were screened by network pharmacology using the TCMSP and DrugBank databases (oral bioavailability≥30%, drug-likeness≥0.18) and intersected with inflammation-related genes retrieved from public databases. The predicted interaction between EGCG and IFN-I pathway proteins (TBK1, IRF3, STAT1) was evaluated by molecular docking, with BX795 and GSK8612 used as reference TBK1 inhibitors. Mechanistic experiments were performed in THP-1-derived macrophages and primary bone marrow-derived macrophages (BMDM). Upstream signaling was activated by transfection of the nucleic acid analogs poly(I∶C) and poly(dA∶dT) or by lipopolysaccharide (LPS) stimulation, whereas downstream signaling was activated by exogenous IFN-β. An siRNA-mediated TREX1 knockdown model was used to mimic endogenous nucleic acid-driven interferonopathy. Expression of IFN-β1 and interferon-stimulated genes (ISGs) was measured by RT-qPCR, protein phosphorylation by Western blot, and IFN-β secretion by ELISA. Cellular thermal shift assay (CETSA) and drug affinity responsive target stability (DARTS) were used to probe the interactionbetween EGCG and IRF3. ResultsNetwork pharmacology identified (-)-epigallocatechin-3-gallate (EGCG) as a candidate IFN-I-suppressive constituent of P. emblica, with predicted binding to TBK1, IRF3, and STAT1. Molecular docking yielded binding energies of -9.2, -7.2, and -8.2 kcal/mol for TBK1, IRF3, and STAT1, respectively, indicating an affinity for TBK1 comparable to that of the reference inhibitors BX795 (-5.7 kcal/mol) and GSK8612 (-6.4 kcal/mol). EGCG suppressed IFN-β1 and ISG mRNA expression under poly (I∶C), poly (dA∶dT), and LPS stimulation in both THP-1 macrophages and BMDM. At the protein level, EGCG reduced the phosphorylation of TBK1 and IRF3 without affecting the levels of the upstream sensors cGAS and RIG-I, and lowered IFN-β secretion in a concentration-dependent manner. CETSA and DARTS showed that EGCG did not enhance the thermal stability or protease resistance of IRF3, indicating that its effect on IRF3 is indirect. Following IFN-β stimulation, prolonged EGCG treatment reduced STAT1 phosphorylation in a time-dependent manner without an apparent change in IRF9, and partially attenuated ISG transcription; this effect was not monotonicly concentration-dependent, and CXCL10 showed the most consistent suppression. In TREX1-knockdown cells, the elevated mRNA levels of ISG15, ISG56, and CXCL10 were reduced by EGCG. ConclusionEGCG suppresses IFN-I responses by concurrently inhibiting TBK1-IRF3-dependent IFN‑β production and JAK-STAT1-mediated downstream transcription. These in vitro findings provide a mechanistic basis for the anti-inflammatory use of P. emblica in traditional medicine and identify EGCG as a candidate for further evaluation in interferon-driven autoimmune disease models.

AIM To explore the differential metabolites of different aged Citrus reticulate'Chachi'and their processed cakes.METHODS Non-targeted metabolomics technology of GC-TOF-MS was used to analyze the chemical constituents.The data was processed by principal component analysis and orthogonal partial least squares discriminant analysis,and the differential metabolites were identified.RESULTS A total of 74 differential metabolites were identified,including 16 glycosides,14 organic acids and their derivatives,11 amino acids and their derivatives,and 4 flavonoids.Comparative analysis revealed 40 and 30 differential metabolites between fresh C.reticulate'℃hachi'and 3-year or 5-year aged samples,respectively.Furthermore,27 and 34 differential metabolites were identified between the 3-year or 5-year aged samples and their corresponding processed cakes,respectively.Differential metabolites among fresh,aged C.reticulate'Chachi',and processed cakes were predominantly enriched in 6 metabolic pathways,including the biosynthesis of secondary metabolites.Specifically,differential metabolites between 3-year aged C.reticulate'Chachi'and its processed cake were significantly enriched in 4 pathways,such as ABC transporters.Differential metabolites between 5-year aged C.reticulate'Chachi'and its processed cake were mainly enriched in 5 pathways,including carbon metabolism.CONCLUSION Non-targeted metabolomics technology can elucidate the chemical compositional differences among fresh/aged and processed cakes of C.reticulate'Chachi',laying a foundation for the research into C.reticulate'Chachi'aging processing techniques and the development of processed products.

AIM To establish the preparation process and quality standard for Zhenggu Pills.METHODS With decoction time,decoction frequency and water addition as influencing factors,comprehensive score for extract yield and transfer rates of epicatechin and naringin as an evaluation index,the decoction process was optimized by orthogonal test.With sugarless paste relative density,medicinal powder fineness,sugarless paste-corn starch ratio,drying temperature and drying time as influencing factors,soft material traits,pill formability,moisture and disintegration time limit as evaluation indices,the formability process was optimized by single factor test.TLC was adopted in the qualitative identification of Dipsaci Radix,salt-processed Psoraleae Fructus,cooked Rhei Radix et Rhizoma and Notoginseng Radix et Rhizoma.HPLC was used for the content determination of paeoniflorin and naringin.RESULTS The optimal decoction process was determined to be 0.5 h for decoction time,two times for decoction frequency,and 10 times for water addition,the comprehensive score was 0.93.The optimal formability process was determined to be 1.21-1.22 for sugarless paste relative density,80 mesh for medicinal powder fineness,1∶0.17-1∶0.18 for sugarless paste-corn starch ratio,70 ℃ for drying temperature,and 24 h for drying time,good soft material traits and pill formability were observable,and moisture and disintegration time limit accored with 2020 edition of Chinese Pharmacopoeia requirements.The TLC spots were clear without negative interference.Two constituents showed good linear relationships within 61.30-490.41 μg/mL(r=0.999 8)and 3.27-26.18 μg/mL(r=0.999 8),whose average recoveries were 100.15％and 98.15％with the RSDs of 0.55％and 2.30％,respectively.CONCLUSION This stable,reliable and specific method can be used for the production and quality evaluation of Zhenggu Pills.

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.

In the process of strontium(Sr)isotope analysis,the separation and purification of Sr are crucial for obtaining high-precision Sr isotope ratios.In this work,ammonia solution was used as a precipitating agent to initially separate and purify Sr from geological samples through coprecipitation,followed by the combination with cation-exchange resin to quickly and effectively obtain high-purity Sr components.The experimental conditions of ammonia coprecipitation reaction were optimized using typical high Rb/Sr ratio(mass ratio)standard sample JR-2.The results showed that using 3 mL of ammonia(1 mol/L)as the precipitating agent,at 120℃for 1 h,could effectively remove more than 90%of Rb,Na and K,with Sr recovery exceeding 95%.To verify the stability and applicability of this separation process,geological standard samples(JR-2,BCR-2,AGV-2,BHVO-2,GSP-2,W-2a)with different rock types and Rb/Sr ratios were practically validated.The results demonstrated that the coprecipitation process with ammonia solution removed more than 90%of Rb,Na and K,with Sr recovery greater than 95%.To further assess the applicability of this method for samples with high Rb/Sr ratios,Rb was added to the rhyolite standard material JR-2,adjusting its Rb/Sr ratio to 500.The results showed that this separation process could still effectively achieve the separation of Rb and Sr,and high-precision Sr isotope analysis results were successfully obtained by combining thermal ionization mass spectrometer(TIMS).Ammonia solution,as a precipitating agent,offered advantages such as low cost,low procedural blanks,simplicity,and a wide range of applicability.Therefore,the combination of ammonia solution coprecipitation and cation-exchange resin enabled rapid Sr purification from geological samples,providing a novel separation approach for high-precision Sr isotope analysis of samples with high Rb/Sr ratios,with broad application prospects.

Near infrared spectroscopy(NIRS)analysis technology has become an important process analysis tool in industrial and agricultural production,and has been widely used for qualitative and quantitative analysis in the fields of tobacco,agriculture,and pharmaceuticals.To address issues such as poor generalization ability and low prediction accuracy in NIRS modeling,a two-dimensional convolutional neural network(2DCNN)quantitative analysis model based on competitive adaptive reweighted sampling(CARS)and Gramian angular difference field(GADF)(CARS-GADF-2DCNN)was proposed.CARS-GADF-2DCNN used the CARS method to select an optimal wavelength set from the full spectrum,then employed GADF to encode the selection results into two-dimensional images,and finally used 2DCNN for prediction analysis.The 2DCNN model consisted of convolutional layers,parallel convolution modules,flattening layer,and fully connected layers.Simulation experiments were conducted on three public near-infrared(NIR)spectral datasets encompassing soil,tablet,and grain datasets to evaluate the CARS-GADF-2DCNN model.The results demonstrated that,compared to the one-dimensional convolutional neural network(1DCNN),the GADF-2DCNN model achieved 16.74％,23.40％,and 7.13％improvement in prediction accuracy for the soil,tablet,and grain datasets,respectively.Compared to GADF-2DCNN,VCPA-GADF-2DCNN,and IRIV-GADF-2DCNN models,the CARS-GADF-2DCNN model further improved prediction accuracy.For the soil dataset,prediction accuracy improved by 39.00％,30.78％and 4.13％;for the tablet dataset,the improvements were 9.52％,6.94％and 2.56％;for the grain dataset,the improvements were 20.57％,9.85％and 15.66％.In conclusion,CARS-GADF-2DCNN effectively selected the optimal wavelength subset from near infrared spectra,and revealed the latent features between different wavelengths.CARS-GADF-2DCNN addresses the issues of high complexity in prediction models and low prediction accuracy in near infrared spectral modeling,and could be effectively applied to near infrared spectral prediction analysis of different substances.