ObjectiveTo establish an in vitro culture model of rat corpora cavernous smooth muscle cells （CCSMCs） using a modified tissue block adherence method. MethodsCorpus cavernosum smooth muscle tissue was digested with collagenase type I and subsequently cultured using an adherent method. Cells were purified via differential adhesion and identified through immunofluorescence and Western blotting. ResultsCCSMCs began to emerge from the tissue block after 3 days， increased significantly by day 7， and converged by day 12. Post-passage， CCSMCs exhibited strong proliferation and a “peak-to-valley” phenomenon. After purification， the cells tested positive for α-smooth muscle actin （α-SMA）， confirming the successful establishment of the in vitro culture model. ConclusionThe modified tissue block adherence method is a cost-effective and efficient way to obtain high-purity CCSMCs.

Osteoarthritis (OA), a highly prevalent degenerative joint disease worldwide, is defined by articular cartilage degradation, abnormal bone remodeling, and persistent chronic inflammation. It severely compromises patients’ quality of life, and currently, there is no radical cure. Abnormal mechanical stress is widely regarded as a core driver of OA pathogenesis, and the exploration of mechanical signal perception and transduction mechanisms has become crucial for deciphering OA’s pathophysiological processes. Piezo1, a key mechanosensitive cation channel belonging to the Piezo protein family, has recently gained significant attention due to its pivotal role in mediating cellular responses to mechanical stimuli in joint tissues. This review systematically examines Piezo1’s expression patterns, regulatory mechanisms, and pathological functions in OA, with a particular focus on its dual roles in modulating chondrocyte homeostasis and bone metabolism disorders, while also delving into the underlying molecular signaling pathways and potential therapeutic implications. Piezo1, consisting of approximately 2 500 amino acids and forming a unique trimeric propeller-like structure, is widely expressed in chondrocytes, osteocytes, mesenchymal stem cells, and synovial cells. It exhibits permeability to cations such as Ca²⁺, K⁺, and Na⁺, and directly responds to membrane tension changes induced by mechanical stimuli like fluid shear stress and mechanical overload. In OA patients and animal models, Piezo1 expression is significantly upregulated, especially in cartilage regions subjected to abnormal mechanical stress (e.g., human temporomandibular joint cartilage). This overexpression is closely associated with aggravated cartilage degeneration, increased chondrocyte apoptosis, accelerated cellular senescence, and intensified inflammatory responses. Mechanical overload and pro-inflammatory cytokines (e.g., IL-1β) are key inducers of Piezo1 upregulation: IL-1β activates the PI3K/AKT/mTOR signaling pathway to enhance Piezo1 expression, forming a pathogenic positive feedback loop that inhibits chondrocyte autophagy, promotes apoptosis, and further accelerates joint degeneration. Mechanistically, Piezo1 mediates OA progression through multiple interconnected pathways. When activated by mechanical stress, Piezo1 triggers excessive Ca²⁺ influx, leading to endoplasmic reticulum stress (ERS) and mitochondrial dysfunction, which directly induce chondrocyte apoptosis. This process involves the activation of downstream signaling cascades such as cGAS-STING and YAP-MMP13/ADAMTS5. YAP, a transcriptional regulator, upregulates the expression of matrix metalloproteinase 13 (MMP13) and aggrecanase (ADAMTS5), thereby accelerating cartilage matrix degradation. Additionally, Piezo1-driven Ca²⁺ overload promotes the accumulation of reactive oxygen species (ROS) and upregulates senescence markers (p16 and p21), accelerating chondrocyte senescence via the p38MAPK and NF-κB pathways. Senescent chondrocytes secrete senescence-associated secretory phenotype (SASP) factors (e.g., IL-6, IL-1β), further amplifying joint inflammation. In terms of bone metabolism, Piezo1 maintains joint homeostasis by promoting the differentiation of fibrocartilage stem cells into chondrocytes and balancing bone formation and resorption through regulating the FoxC1/YAP axis and RANKL/OPG ratio. Therapeutically, targeting Piezo1 shows promising potential. Preclinical studies have demonstrated that Piezo1 inhibitors (e.g., GsMTx4) can reduce joint damage and alleviate pain in OA mice. Simultaneously, siRNA-mediated co-silencing of Piezo1 and TRPV4 (another mechanosensitive channel) decreases intracellular Ca²⁺ concentration, inhibits chondrocyte apoptosis, and promotes cartilage repair. Conditional knockout of Piezo1 using Gdf5-Cre transgenic mice alleviates cartilage degeneration in post-traumatic OA models by downregulating MMP13 and ADAMTS5 expression. Despite existing challenges, such as off-target effects of inhibitors, inefficient local drug delivery, and interindividual genetic variability, strategies like developing selective Piezo1 antagonists, optimizing targeted nanocarriers, and combining Piezo1-targeted therapy with physical therapy provide viable avenues for clinical translation. The authors propose that Piezo1 serves as a critical therapeutic target for OA, and future research should focus on deciphering its context-dependent regulatory networks, developing tissue-specific intervention strategies, and validating their efficacy and safety in clinical trials to address the unmet medical needs of OA patients.

ObjectiveTo construct and validate a clinical prediction model for inflammatory remission outcomes in rheumatoid arthritis （RA） patients with liver and kidney deficiency syndrome treated with Bushen Zhiwang Decoction （补肾治尪汤， BZD） based on metabolomics. MethodsA prospective cohort study was conducted， enrol-ling 60 RA patients with liver and kidney deficiency syndrome. All patients were treated with BZD and conventional-dose oral conventional synthetic disease-modifying antirheumatic drugs （csDMARDs） for 12 months. Clinical data were collected， and the change in disease activity score in 28 joints （DAS28） after treatment compared with baseline （△DAS28） was used as the primary outcome and grouping criterion. Peripheral blood samples were collected before treatment to analyze plasma metabolites. Differential analysis and least absolute shrinkage and selection operator （LASSO） regression were used to preliminarily screen differential metabolites， followed by machine learning algorithms to further identify a core metabolite combination. Based on the expression levels of the core metabolite combination， a novel metabolite index， namely the metabolomics-based inflammatory remission score （Met-IRS）， was calculated using standar-dized metabolite values， and its clinical applicability was evaluated. A clinical prediction model was constructed by integrating clinical characteristics and Met-IRS， and the model performance was assessed. ResultsAmong the 60 patients， those with △DAS28 ≥ 0.27 were assigned to the high inflammatory remission group， while those with △DAS28 ＜ 0.27 were assigned to the low inflammatory remission group， with 30 cases in each group. Compared to the low inflammatory remission group， the high inflammatory remission group showed a higher frequency of methotrexate use and a lower positive rate of rheumatoid factor （RF） （P＜0.05）. Seven core metabolites were identified as the optimal combination， including mangiferic acid， fatty acid-hydroxy fatty acid ester 40∶6， fatty acid-hydroxy fatty acid ester 18∶0， fatty acid-hydroxy fatty acid ester 36∶1， glucosylceramide， lysophosphatidylcholine 22∶5， and pregnanetriol ketone. The calculated Met-IRS comprehensively reflected the characteristics of differential metabolites and demonstrated clinical applicability. Met-IRS was significantly higher in the high inflammatory remission group than in the low inflammatory remission group， and was positively correlated with high inflammatory remission outcomes （P＜0.05）. Based on the variables Met-IRS， methotrexate use， leflunomide use， and RF positivity， a clinical prediction model for inflammatory remission in RA treatment （Cj-RTRM） was constructed. Model performance evaluation demonstrated that the model had good clinical predictive ability， with an area under the receiver operating characteristic curve （AUC） of 0.880， sensitivity 0.967， specificity 0.700 and Youden's index 0.667. ConclusionThe clinical prediction model Cj-RTRM constructed based on the metabolomics-based inflammatory remission score Met-IRS can effectively predict clinical inflammatory remission outcomes in RA patients treated with BZD and accurately identify the advantageous population for this treatment. This model provides guiding evidence for dynamic inflammation monitoring， targeted management， and identification of populations with advantages in traditional Chinese medicine.

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.

Andrological diseases have become an important public health problem threatening men's health worldwide， which significantly affects the quality of life of patients and brings a heavy disease burden. Western medicine often faces the dilemma of obvious side effects and limited efficacy. Traditional Chinese medicine has unique advantages in the prevention and treatment of andrological diseases and has accumulated rich clinical experience. Epimedii Folium， as a traditional Chinese medicine for strengthening kidney and Yang， exerts a key therapeutic effect on andrology diseases through multi-component synergy， multi-target regulation， and multi-pathway intervention. Recent studies have found that the main active components of Epimedii Folium， such as icariin， icariside， and icaritin， are the key material basis for the treatment of andrological diseases. The active components of Epimedii Folium can play a role in common andrological diseases such as erectile dysfunction， male infertility， and prostate cancer by regulating the activity of the nitric oxide/cyclic guanosine monophosphate （NO/cGMP） pathway， participating in oxidative stress response， regulating the secretion of hypothalamic-pituitary-gonadal axis hormones， improving spermatogenic dysfunction， and inhibiting the proliferation of cancer cells. However， the systematic action network and molecular mechanisms of the active components of Epimedii Folium have not been fully elucidated， thereby limiting its potential for clinical translation and application. In the future， it is necessary to combine cutting-edge technologies such as metabolomics， single-cell sequencing， and targeted nanoscale drug delivery systems， strengthening the research on the compatibility rules of active components and organ-specific delivery， providing a scientific basis for the development of innovative andrology traditional Chinese medicine formulas with international competitiveness， and promoting the innovation and breakthrough of andrology disease treatment modes.

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.

Objective: To evaluate the feasibility of dynamic contrast-enhanced MRI (DCE-MRI) in differentiating tumor deposits (TDs) from metastatic lymph nodes (MLNs) in rectal cancer. Materials and Methods: A retrospective analysis was conducted on 70 patients with rectal cancer, including 168 lesions (70 TDs and 98 MLNs confirmed by histopathology), who underwent pretreatment MRI and subsequent surgery between March 2019 and December 2022. The morphological characteristics of TDs and MLNs, along with quantitative parameters derived from DCE-MRI (K trans , kep, and v e) and DWI (ADCmin, ADCmax, and ADCmean), were analyzed and compared between the two groups.Multivariable binary logistic regression and receiver operating characteristic (ROC) curve analyses were performed to assess the diagnostic performance of significant individual quantitative parameters and combined parameters in distinguishing TDs from MLNs. Results: All morphological features, including size, shape, border, and signal intensity, as well as all DCE-MRI parameters showed significant differences between TDs and MLNs (all P < 0.05). However, ADC values did not demonstrate significant differences (all P > 0.05). Among the single quantitative parameters, v e had the highest diagnostic accuracy, with an area under the ROC curve (AUC) of 0.772 for distinguishing TDs from MLNs. A multivariable logistic regression model incorporating short axis, border, v e, and ADC mean improved diagnostic performance, achieving an AUC of 0.833 (P = 0.027). Conclusion The combination of morphological features, DCE-MRI parameters, and ADC values can effectively aid in the preoperative differentiation of TDs from MLNs in rectal cancer.

Objective To investigate the expression of GJB4 gene in pancreatic cancer tissue and its correlation with clinicopathology.Methods The expression levels of GJB4 mRNA in pancreatic cancer and adjacent cancer tissues were analyzed using bioinformatics to analyze the Cancer Genome Atlas(TCGA)and Genotype-Tissue Expression(GTEx)RNA sequencing datasets.A total of 120 pancreatic cancer samples and 40 adjacent cancer samples from the Pathology Department of The Second Affiliated Hospital of Kunming Medical University from January 2019 to December 2023 were collected and sorted.Immunohistochemistry staining method was used to detect the expression difference of GJB4 protein between the two groups.RT-qPCR method was used to detect the expression levels of GJB4 in four pancreatic cancer cell lines.Univariate and multivariate Cox regression and Kaplan-Meier curves were used to analyze the clinical pathological data of GJB4 and pancreatic cancer patients.DAVID functional annotation bioinformatics and GSEA enrichment analysis were used to explore the relevant pathways of GJB4 in pancreatic cancer.Results The expression level of GJB4 mRNA in pancreatic cancer was higher than that in adjacent tissues,and the high expression of GJB4 was significantly associated with poor prognosis of patients(P<0.05).Immunohistochemical analysis showed that GJB4 protein was brown-yellow granular in pancreatic cancer tissues,mainly expressed in cytoplasm and cell membrane,and GJB4 protein expression was up-regulated in pancreatic cancer(P<0.05).The RT-qPCR test results showed that out of 4 pancreatic cancer cell lines,3 of them had upregulated expression(P<0.05).COX regression analysis showed that GJB4 gene was an independent risk factor in the prognosis of pancreatic cancer patients.The results of GO enrichment analysis showed that GJB4 was related to the transmembrane transport,ion channel,signal release and membrane potential regulation of pancreatic cancer.GSEA analysis showed that GJB4 was enriched in the Wnt/β-catenin signaling pathway.Conclusion In pancreatic cancer,the high expression level of GJB4 is closely related to the clinicopathological features of the patients,which may predict the poor prognosis of the patients.GJB4 may be involved in regulating the Wnt/β-catenin signaling pathway of pancreatic cancer,and is expected to be one of the potential biomarkers of pancreatic cancer prognosis.

Energy metabolism regulation plays a pivotal role in metabolic engineering. It mainly achieves the balance of material and energy metabolism or maximizes the utilization of materials and energy by regulating the supply intensity and mode of ATP and reducing electron carriers in cells. On the one hand, the production efficiency can be increased by changing the distribution of material metabolic flow. On the other hand, the thermodynamic parameters of enzyme-catalyzed reactions can be altered to affect the reaction balance, and thus the production costs are reduced. Therefore, energy metabolism regulation is expected to become a favorable tool for the modification of microbial cell factories, thereby increasing the production of target metabolites and reducing production costs. This article introduces the commonly used energy metabolism regulation methods and their effects on cell factories, aiming to provide a reference for the efficient construction of microbial cell factories.
Energy Metabolism/physiology* ; Metabolic Engineering/methods* ; Adenosine Triphosphate/metabolism* ; Industrial Microbiology/methods*