Objective To investigate the causal relationship between gut microbiota and idiopathic pulmonary fibrosis (IPF). Methods Genome-wide association studies (GWAS) data of gut microbiota and IPF were obtained from MiBioGen and IEU OpenGWAS, respectively. Instrumental variables were screened by means of significance, linkage disequilibrium, weak instrumental variable screening, and removal of confounding factors (genetics, smoking, host characteristics). Inverse variance weighted (IVW) was used as the main Mendelian randomization (MR) analysis method, and the weighted median, simple mode, MR-Egger, and weighted mode were used to perform MR to reveal the causal effect of gut microbiota and IPF. The Cochrane's Q, leave-one-out, MR-Egger-intercept, and Mendelian randomization pleiotropy residual sum and outlier (MR-PRESSO) and Steiger tests were used to analyze the heterogeneity, horizontal pleiotropy, outliers, and directionality, respectively. Results IVW analysis results showed that Actinobacteria [OR=1.773, 95%CI (1.323, 2.377), P<0.001], Erysipelatoclostridium [OR=2.077, 95%CI (1.107, 3.896), P=0.023], and Streptococcus [OR=1.35, 95%CI (1.100, 1.657), P=0.004] could increase the risk of IPF. Bifidobacterium [OR=0.668, 95%CI (0.620, 0.720), P<0.001], Ruminococcus [OR=0.434, 95%CI (0.222, 0.848), P=0.015], and Tyzzerella [OR=0.479, 95%CI (0.304, 0.755), P=0.001] could reduce the risk of IPF. No significant heterogeneity, horizontal pleiotropy, outliers, and reverse causality were found. Conclusion Actinobacteria, Erysipelatoclostridium and Streptococcus may increase the risk of IPF, while Bifidobacterium, Ruminococcus and Tyzzerella may reduce the risk of IPF. Regulation of the above gut microbiota may become a new direction in the study of the pathogenesis of IPF.

ObjectiveTo evaluate the antitumor activity of Periplaneta americana extract（PAE） against human-derived lung adenocarcinoma organoids（LUAD-PDOs） and to elucidate its potential mechanism based on transcriptomics. MethodsFresh tumor and adjacent normal tissues from patients with LUAD were collected to construct LUAD-PDOs and normal lung organoid（Nor-PDOs） models using 3D organoid culture technology. The effective intervention concentration of PAE was determined using the cell counting kit-8（CCK-8） assay. Experimental groups included the model group（LUAD-PDOs）， normal group， model administration group（LUAD-PDOs+PAE）， and normal administration group（Nor-PDOs+PAE）. Hematoxylin-eosin（HE） staining was used to observe the pathological structures of PDOs， immunohistochemistry（IHC） was performed to detect the expressions of the proliferation marker Ki-67 and lung adenocarcinoma differentiation markers cytokeratin-7（CK-7） and Napsin A， TUNEL staining was applied to detect cell apoptosis. RNA sequencing（RNA-Seq） was conducted to identify differentially expressed genes（DEGs）， followed by Gene Ontology（GO）， Kyoto Encyclopedia of Genes and Genomes（KEGG）， and Gene Set Enrichment Analysis（GSEA）， alongside protein-protein interaction（PPI） network analysis to screen core mechanisms. Finally， key targets were validated by integrating external database analysis with immunofluorescence（IF）. ResultsNor-PDOs and LUAD-PDOs that highly recapitulated the pathological characteristics of the primary tissues were successfully established. The CCK-8 assay determined that the effective intervention concentration of PAE was 16 g·L^-1. Morphological observation showed that Nor-PDOs exhibited lumen-forming structures， whereas LUAD-PDOs displayed dense， solid structures. CCK-8 and TUNEL assays revealed that， compared with the model group， PAE intervention inhibited the proliferation of LUAD-PDOs and promoted apoptosis in LUAD cells， while showing no significant effect on the viability of Nor-PDOs. Transcriptomic analysis identified 719 DEGs that were significantly reversed after PAE intervention（347 up-regulated and 372 down-regulated）（P<0.05）. GO enrichment analysis indicated that DEGs in the model administration group were significantly enriched in biological processes related to cell cycle regulation compared to the model group. KEGG pathway analysis revealed that PAE affected pathways related to proliferation and metabolism， including pathways in cancer and the p53 signaling pathway. GSEA further confirmed that PAE significantly enhanced the activity of the p53 signaling pathway（P<0.05）. PPI network analysis indicated that breast cancer type 1 susceptibility protein（BRCA1） and checkpoint kinase 1（CHEK1） were the core down-regulated targets in the p53 pathway. IF verified the high expression of BRCA1 and CHEK1 in LUAD-PDOs and their significant downregulation after PAE intervention（P<0.05）. Furthermore， survival analysis based on The Cancer Genome Atlas（TCGA） database indicated that low expression of BRCA1 and CHEK1 was significantly associated with prolonged overall survival in patients with LUAD（P<0.05）. ConclusionPAE effectively inhibits proliferation of LUAD-PDOs and promotes their apoptosis， its anti-tumor mechanism is potentially associated with the activation of the p53 signaling pathway， with BRCA1 and CHEK1 genes likely serving as key downstream targets for the effects of PAE.

ObjectiveCleft palate (CP) is a common congenital deformity often associated with velopharyngeal insufficiency (VPI), which disrupts the physiological coupling between respiration and speech. Conventional clinical assessments, such as nasometry and spirometry, provide limited static data and fail to visualize the dynamic spatiotemporal distribution of lung ventilation during phonation. This study introduces spatiotemporal electrical impedance tomography (ST-EIT) to evaluate speech-respiratory functional features in CP patients compared to normal controls (NC). The aim is to characterize multi-domain respiratory patterns and to validate an interpretable machine learning framework for providing objective, quantitative evidence for clinical assessment. MethodsSeventy-five participants were enrolled in this study, comprising 37 patients with surgically repaired CP and 38 healthy volunteers matched for age, gender, and body mass index (BMI). All subjects performed standardized sustained phonation tasks while undergoing synchronous monitoring with a 16-electrode EIT system and a pneumotachograph. A comprehensive feature engineering pipeline was developed to extract physiological parameters across 3 complementary domains. (1) Temporal domain: including inspiratory/expiratory phase duration (tPhase), time constants (Tau), and inspiratory-to-expiratory time ratios (TI/TE); (2) airflow domain: comprising mean flow, peak flow, and instantaneous flow at 25%, 50%, and 75% of tidal volume; and (3) spatial domain: quantifying global and regional tidal impedance variation (TIV), global inhomogeneity (GI), and center of ventilation (CoV). Extreme Gradient Boosting (XGBoost) classifiers were trained using 5 distinct data sources (Spirometry, Nasometry, Inspiratory-EIT, Expiratory-EIT, and fused ST-EIT). Model performance was rigorously evaluated via stratified 5-fold cross-validation, and Shapley additive explanations (SHAP) were employed to quantify global and local feature contributions. ResultsThe CP group exhibited a distinct respiratory phenotype compared to controls. In the temporal domain, CP patients showed significantly shorter inspiratory (1.60 s vs.1.85 s, P<0.001) and expiratory phase durations (2.45 s vs. 3.95 s, P<0.001), indicating a rapid, shallow breathing rhythm. In the airflow domain, while inspiratory flows were comparable, the CP group demonstrated significantly elevated mean and peak flows during the expiratory phase (P<0.001), reflecting compensatory respiratory effort. Spatially, CP patients presented significant ventilation redistribution, characterized by higher regional TIV in the right-anterior (ROI1) and left-posterior (ROI4) quadrants, but lower TIV in the left-anterior (ROI2) quadrant. In terms of diagnostic accuracy, the multi-modal ST-EIT model achieved the highest performance (AUC: 0.915±0.012, Accuracy: 0.843±0.019, F1-score: 0.872±0.017), substantially outperforming models based on spirometry (AUC: 0.721) or nasometry (AUC: 0.625) alone. Interpretability analysis revealed that spatial domain features were the most critical, contributing 53.4% to the model’s decision-making, followed by temporal (25.0%) and airflow (21.6%) features. ConclusionST-EIT successfully captures the temporal, airflow, and spatial deviations in CP speech respiration that are undetectable by conventional methods—specifically, rapid phase transitions, hyperdynamic expiratory airflow, and regional ventilation heterogeneity. This study validates ST-EIT as a robust, non-invasive, and radiation-free tool for characterizing speech-respiratory dysfunction, offering high clinical value for bedside screening, rehabilitation planning, and longitudinal monitoring of patients with cleft palate.

ObjectiveCleft palate (CP) is a common congenital deformity often associated with velopharyngeal insufficiency (VPI), which disrupts the physiological coupling between respiration and speech. Conventional clinical assessments, such as nasometry and spirometry, provide limited static data and fail to visualize the dynamic spatiotemporal distribution of lung ventilation during phonation. This study introduces spatiotemporal electrical impedance tomography (ST-EIT) to evaluate speech-respiratory functional features in CP patients compared to normal controls (NC). The aim is to characterize multi-domain respiratory patterns and to validate an interpretable machine learning framework for providing objective, quantitative evidence for clinical assessment. MethodsSeventy-five participants were enrolled in this study, comprising 37 patients with surgically repaired CP and 38 healthy volunteers matched for age, gender, and body mass index (BMI). All subjects performed standardized sustained phonation tasks while undergoing synchronous monitoring with a 16-electrode EIT system and a pneumotachograph. A comprehensive feature engineering pipeline was developed to extract physiological parameters across 3 complementary domains. (1) Temporal domain: including inspiratory/expiratory phase duration (tPhase), time constants (Tau), and inspiratory-to-expiratory time ratios (TI/TE); (2) airflow domain: comprising mean flow, peak flow, and instantaneous flow at 25%, 50%, and 75% of tidal volume; and (3) spatial domain: quantifying global and regional tidal impedance variation (TIV), global inhomogeneity (GI), and center of ventilation (CoV). Extreme Gradient Boosting (XGBoost) classifiers were trained using 5 distinct data sources (Spirometry, Nasometry, Inspiratory-EIT, Expiratory-EIT, and fused ST-EIT). Model performance was rigorously evaluated via stratified 5-fold cross-validation, and Shapley additive explanations (SHAP) were employed to quantify global and local feature contributions. ResultsThe CP group exhibited a distinct respiratory phenotype compared to controls. In the temporal domain, CP patients showed significantly shorter inspiratory (1.60 s vs.1.85 s, P<0.001) and expiratory phase durations (2.45 s vs. 3.95 s, P<0.001), indicating a rapid, shallow breathing rhythm. In the airflow domain, while inspiratory flows were comparable, the CP group demonstrated significantly elevated mean and peak flows during the expiratory phase (P<0.001), reflecting compensatory respiratory effort. Spatially, CP patients presented significant ventilation redistribution, characterized by higher regional TIV in the right-anterior (ROI1) and left-posterior (ROI4) quadrants, but lower TIV in the left-anterior (ROI2) quadrant. In terms of diagnostic accuracy, the multi-modal ST-EIT model achieved the highest performance (AUC: 0.915±0.012, Accuracy: 0.843±0.019, F1-score: 0.872±0.017), substantially outperforming models based on spirometry (AUC: 0.721) or nasometry (AUC: 0.625) alone. Interpretability analysis revealed that spatial domain features were the most critical, contributing 53.4% to the model’s decision-making, followed by temporal (25.0%) and airflow (21.6%) features. ConclusionST-EIT successfully captures the temporal, airflow, and spatial deviations in CP speech respiration that are undetectable by conventional methods—specifically, rapid phase transitions, hyperdynamic expiratory airflow, and regional ventilation heterogeneity. This study validates ST-EIT as a robust, non-invasive, and radiation-free tool for characterizing speech-respiratory dysfunction, offering high clinical value for bedside screening, rehabilitation planning, and longitudinal monitoring of patients with cleft palate.

Metabolic dysfunction-associated steatotic liver disease (MASLD) has become the most prevalent chronic liver disease worldwide, affecting approximately 32%-38% of the adult population and posing a growing public health burden. MASLD represents a continuous disease spectrum ranging from simple steatosis to metabolic dysfunction-associated steatohepatitis (MASH), progressive hepatic fibrosis, cirrhosis, and ultimately hepatocellular carcinoma (HCC). The pathological core of MASLD lies in disruption of hepatic lipid metabolic homeostasis, characterized by an imbalance among de novo lipogenesis, fatty acid β-oxidation, and very-low-density lipoprotein (VLDL)-mediated lipid export. This metabolic disequilibrium subsequently drives inflammatory injury and fibrotic progression. Among the multiple regulatory pathways involved, thyroid hormone (TH) signaling has emerged as a central regulator of hepatic metabolic homeostasis. The liver is a major peripheral target organ of TH action, where TH predominantly exerts its metabolic effects through thyroid hormone receptor β (TRβ). Large-scale epidemiological studies and meta-analyses have demonstrated that hypothyroidism is significantly associated with increased MASLD prevalence, more severe histological injury, and advanced hepatic fibrosis, suggesting that dysregulation of TH signaling may participate throughout the entire MASLD disease spectrum. At the molecular level, TH regulates hepatic lipid metabolism by coordinating suppression of lipogenesis, enhancement of mitochondrial fatty acid oxidation, and promotion of VLDL assembly and secretion through integrated genomic actions of the T3-TRβ axis and non-genomic signaling pathways. Across different stages of MASLD, TH signaling exerts stage-dependent protective effects. In the steatosis stage, TH improves metabolic flexibility by modulating insulin sensitivity, glucose metabolism, and lipid droplet clearance, thereby alleviating early lipotoxic stress. During progression to MASH, TH attenuates inflammatory amplification by improving mitochondrial homeostasis, suppressing activation of the NOD-like receptor family pyrin domain containing 3 (NLRP3) inflammasome, and modulating the gut-liver axis microenvironment. In advanced stages, TH signaling influences hepatic stellate cell activation and extracellular matrix deposition, partly through interaction with the transforming growth factor-β (TGF-β)/SMAD pathway, while alterations in intrahepatic TH availability, mediated by dynamic changes in iodothyronine deiodinase 1 (DIO1), contribute to fibrosis progression and hepatocellular dedifferentiation. In hepatocellular carcinoma, coordinated downregulation of TRβ and DIO1 establishes a tumor-associated hypothyroid state that promotes metabolic reprogramming and tumor progression. The clinical relevance of TH signaling in MASLD has been underscored by the recent approval of Resmetirom, a liver-targeted TRβ‑selective agonist, for the treatment of non-cirrhotic MASH with moderate-to-severe fibrosis (F2-F3). This approval represents a landmark transition from mechanistic understanding to metabolism-centered precision therapy in MASLD. Clinical trials have demonstrated that Resmetirom not only improves key histological endpoints, including MASH resolution and fibrosis regression, but also favorably modulates atherogenic lipid profiles, highlighting the therapeutic potential of selectively targeting hepatic TH pathways. This review systematically summarizes the multidimensional regulatory roles of TH across the MASLD disease spectrum and discusses emerging diagnostic and therapeutic implications of TH-based interventions, aiming to inform future mechanistic research and optimize clinical management strategies.

ObjectiveTo investigate the correlation between neutrophil-to-lymphocyte ratio （NLR）， platelet-to-lymphocyte ratio （PLR）， systemic immune-inflammation index （SII） and the severity of intrauterine adhesions （IUA）. MethodsThe retrospective study included 380 patients who underwent transcervical resection of adhesions （TCRA） from December 2019 to March 2025. Based on the American Fertility Society （AFS） classification， patients were divided into mild （n=61）， moderate （n=225）， and severe （n=94） groups. NLR， PLR， and SII were calculated from preoperative blood tests. Statistical analyses included Kruskal-Wallis test and ordinal Logistic regression. ResultsNLR， PLR， and SII were significantly higher in the severe IUA group compared to the mild group （P<0.05）， with SII showing the strongest predictive ability （OR=1.004， P=0.001）. The number of intrauterine procedures was an independent risk factor （OR=1.27/level， P=0.016）. The predictive model ［Logit（P）=-0.676+0.241×operation times+0.004×SII］ effectively identified severe IUA cases. ConclusionInflammatory markers （particularly SII） are correlated with IUA severity and may serve as non-invasive tools for clinical assessment.

BACKGROUND:Choline kinase alpha is a key enzyme in phospholipid metabolism,involved in the synthesis of phosphatidylcholine,and plays an important role in maintaining cell membrane integrity and signal transduction.Research has shown that choline kinase alpha is highly expressed in various tumors and is closely related to cell proliferation,metabolic reprogramming,and tumor progression.As a potential therapeutic target,the role of choline kinase alpha in tumor metabolism and mitochondrial function still needs further exploration.OBJECTIVE:To evaluate the effects and the underlying mechanisms of choline kinase alpha on the proliferation and apoptosis of glioma U87MG and U251 cells.METHODS:Short hairpin RNA of choline kinase alpha and its empty vector control were transfected into U87MG and U251 glioma cells.Mitochondrial morphology was observed by transmission electron microscopy.Mitochondrial structure and functional protein levels were assessed by western blot assay.Reactive oxygen species levels in cells were measured using a reactive oxygen species fluorescent probe.Mitochondrial membrane potential was assessed with a JC-1 assay.Intracellular adenosine triphosphate levels were measured by chemiluminescence.Cell proliferation was evaluated using a CCK-8 assay.Apoptosis levels were analyzed by flow cytometry.The mitochondrial fission inhibitor Mdivi-1 was used to protect the mitochondrial function of the choline kinase α-silenced lentiviral cells.Finally,U87MG cells were subcutaneously injected to construct a subcutaneous tumor model in nude mice.The tumor growth in nude mice was observed before and after choline kinase alpha silencing and after the use of the mitochondrial fission inhibitor Mdivi-1.RESULTS AND CONCLUSION:(1)Compared with the empty control group,the mitochondria of U87MG and U251 cells in the choline kinase alpha silencing lentivirus group exhibited significant structural abnormalities in mitochondria,such as vacuolization and cristae disruption.The expressions of mitochondrial structure and function-related proteins TOM20,ACO2,and ATP5A were significantly decreased(P＜0.01,P＜0.001),the expression of SOD2 was significantly increased(P＜0.01,P＜0.000 1),the fluorescence intensity of reactive oxygen species was significantly increased(P＜0.01),the mitochondrial membrane potential and adenosine triphosphate level were significantly decreased(P＜0.01,P＜0.001),the cell proliferation ability was reduced(P＜0.01),and the apoptosis level was increased(P＜0.001).(2)Following Mdivi-1 treatment,the fluorescence intensity of reactive oxygen species in U87MG and U251 cells decreased(P＜0.05,P＜0.01),mitochondrial membrane potential and adenosine triphosphate levels were significantly restored(P＜0.05,P＜0.01,P＜0.001),cell proliferation ability was improved(P＜0.05,P＜0.01),and apoptosis level was decreased(P＜0.05).(3)In addition,the in vitro subcutaneous tumor formation experiment of nude mice showed that compared with the empty control group,the mass and growth rate of subcutaneous tumors formed by U87MG cells in the choline kinase alpha silencing lentivirus group were significantly reduced(P＜0.000 1).After Mdivi-1 treatment,the mass and growth rate of tumors were significantly increased(P＜0.000 1).(4)The results show that choline kinase alpha silencing affects the proliferation and apoptosis of glioma cells by inducing mitochondrial dysfunction.

BACKGROUND:Choline kinase alpha is a key enzyme in phospholipid metabolism,involved in the synthesis of phosphatidylcholine,and plays an important role in maintaining cell membrane integrity and signal transduction.Research has shown that choline kinase alpha is highly expressed in various tumors and is closely related to cell proliferation,metabolic reprogramming,and tumor progression.As a potential therapeutic target,the role of choline kinase alpha in tumor metabolism and mitochondrial function still needs further exploration.OBJECTIVE:To evaluate the effects and the underlying mechanisms of choline kinase alpha on the proliferation and apoptosis of glioma U87MG and U251 cells.METHODS:Short hairpin RNA of choline kinase alpha and its empty vector control were transfected into U87MG and U251 glioma cells.Mitochondrial morphology was observed by transmission electron microscopy.Mitochondrial structure and functional protein levels were assessed by western blot assay.Reactive oxygen species levels in cells were measured using a reactive oxygen species fluorescent probe.Mitochondrial membrane potential was assessed with a JC-1 assay.Intracellular adenosine triphosphate levels were measured by chemiluminescence.Cell proliferation was evaluated using a CCK-8 assay.Apoptosis levels were analyzed by flow cytometry.The mitochondrial fission inhibitor Mdivi-1 was used to protect the mitochondrial function of the choline kinase α-silenced lentiviral cells.Finally,U87MG cells were subcutaneously injected to construct a subcutaneous tumor model in nude mice.The tumor growth in nude mice was observed before and after choline kinase alpha silencing and after the use of the mitochondrial fission inhibitor Mdivi-1.RESULTS AND CONCLUSION:(1)Compared with the empty control group,the mitochondria of U87MG and U251 cells in the choline kinase alpha silencing lentivirus group exhibited significant structural abnormalities in mitochondria,such as vacuolization and cristae disruption.The expressions of mitochondrial structure and function-related proteins TOM20,ACO2,and ATP5A were significantly decreased(P＜0.01,P＜0.001),the expression of SOD2 was significantly increased(P＜0.01,P＜0.000 1),the fluorescence intensity of reactive oxygen species was significantly increased(P＜0.01),the mitochondrial membrane potential and adenosine triphosphate level were significantly decreased(P＜0.01,P＜0.001),the cell proliferation ability was reduced(P＜0.01),and the apoptosis level was increased(P＜0.001).(2)Following Mdivi-1 treatment,the fluorescence intensity of reactive oxygen species in U87MG and U251 cells decreased(P＜0.05,P＜0.01),mitochondrial membrane potential and adenosine triphosphate levels were significantly restored(P＜0.05,P＜0.01,P＜0.001),cell proliferation ability was improved(P＜0.05,P＜0.01),and apoptosis level was decreased(P＜0.05).(3)In addition,the in vitro subcutaneous tumor formation experiment of nude mice showed that compared with the empty control group,the mass and growth rate of subcutaneous tumors formed by U87MG cells in the choline kinase alpha silencing lentivirus group were significantly reduced(P＜0.000 1).After Mdivi-1 treatment,the mass and growth rate of tumors were significantly increased(P＜0.000 1).(4)The results show that choline kinase alpha silencing affects the proliferation and apoptosis of glioma cells by inducing mitochondrial dysfunction.

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.