Animal model research on spleen and stomach diseases in traditional Chinese medicine （TCM） is of great significance for elucidating the nature of diseases and syndromes and for revealing the mechanisms of action of Chinese herbal medicinals. At present， studies on classical TCM syndrome models of spleen and stomach diseases mainly focus on spleen deficiency syndrome， liver constraint syndrome， and damp-heat syndrome. Model construction is mostly based on the etiological and pathophysiological characteristics of syndrome， and model evaluation primarily involves macroscopic manifestations and physicochemical indicators. This paper summarizes the current research status of animal models integrating disease and syndrome for seven common spleen and stomach diseases， including chronic gastritis and gastric precancerous lesions， gastroesophageal reflux disease， functional dyspepsia， inflammatory bowel disease， irritable bowel syndrome， functional constipation， and functional diarrhea. The modeling methods and characteristics of disease-syndrome combined animal models for each disease are analyzed. It is proposed that future research on disease-syndrome integration in spleen and stomach diseases should move toward syste-matic， precise， and integrative development， and that interdisciplinary and cross-disciplinary research approaches should be adopted to enhance the predictive value and application efficiency of disease-syndrome combined animal models.

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

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

Objective: This prospective study evaluated the feasibility of chemical exchange saturation transfer (CEST) MRI for assessing disease activity in thyroid-associated ophthalmopathy (TAO). Materials and Methods: A total of 88 patients with active TAO, 76 with inactive TAO, and 30 healthy controls were enrolled. CEST MRI-derived magnetization transfer ratio (MTR) and MTR asymmetry (MTRasym) at 1 ppm, 2 ppm, and 3.5 ppm were calculated. Clinical data, MTR, and MTRasym values for the extraocular muscles (one representative muscle per eye, yielding two measurements per participant) were compared among the groups. Spearman’s correlation was used to examine associations between imaging parameters and the clinical activity score (CAS) in patients with TAO. Logistic regression analysis was used to identify independent associations between imaging parameters and disease activity in patients with TAO (active vs. inactive). Receiver operating characteristic (ROC) analysis was conducted to evaluate the diagnostic performance for discriminating active from inactive TAO. Results: Patients with active TAO showed lower MTR values (P < 0.001) and higher MTRasym (1 ppm), MTRasym (2 ppm), and MTRasym (3.5 ppm) (all P < 0.001) compared with those with inactive TAO. MTR was negatively correlated with CAS (r = -0.402; P < 0.001), while MTRasym (1 ppm), MTRasym (2 ppm), and MTRasym (3.5 ppm) were positively correlated with CAS (r = 0.369, 0.350, and 0.349, respectively;all P < 0.001). MTR and MTRasym (1 ppm) were independently associated with TAO activity. The areas under the ROC curve (AUCs) for MTR and MTRasym (1 ppm) in discriminating active from inactive TAO were 0.772 and 0.730, respectively. Combining MTR with MTRasym (1 ppm) significantly improved diagnostic performance compared with either parameter alone, achieving an AUC of 0.805 (P = 0.029 and 0.001). Conclusion MTR and MTRasym (1 ppm) were independently associated with TAO activity. Their combination further enhanced diagnostic performance in distinguishing active from inactive TAO, suggesting their potential as quantitative imaging biomarkers to guide treatment in patients with TAO.

ObjectiveTo investigate the effect of honokiol on proliferation， apoptosis， migration， and invasion of gastric cancer cells and its underlying mechanistic. MethodsHuman gastric cancer cell lines HGC-27 and AGS were treated with Honokiol at concentrations of 0， 15， and 25 μmol/L. CCK-8 assays were conducted to determine the half maximal inhibitory concentration （IC₅₀） for both cell lines. Cell viability， proliferation， migration， and invasion capabilities were assessed using CCK-8， colony formation， wound healing， Transwell migration and Transwell invasion assays. Apoptosis rates were measured via flow cytometry. Western blot analysis examined proteins related to proliferation， apoptosis， migration， invasion， and the endoplasmic reticulum stress pathway ATF4-CHOP-TRIB3. ResultsCompared with the control group， treatment with 15 and 25 μmol/L Honokiol significantly reduced the proliferation， colony formation， migration， and invasion capabilities of the two gastric cancer cell lines， while significantly increasing the apoptosis rate （P<0.05）. Additionally， compared to the control group， the protein expression levels of neural cadherin（N-cadherin）， Vimentin， proliferating cell nuclear antigen（PCNA）， and B-cell lymphoma/leukemia-2 protein（Bcl-2）decreased in the two gastric cancer cell lines after treatment with 15 and 25 μmol/L Honokiol， while the protein expression levels of epithelial cadherin（E-cadherin）， Bcl-2-associated X protein（Bax）， activating transcription factor 4（ATF4）， endoplasmic reticulum stress-related protein（CHOP）， and tribbles homolog 3（TRIB3）increased（P<0.05）. ConclusionHonokiol promotes apoptosis and inhibits proliferation， migration， and invasion of gastric cancer HGC-27 and AGS cells by regulating the ERS signaling pathway ATF4/CHOP/TRIB3.

Gastric cancer remains a highly prevalent malignancy worldwide, with surgical resection currently constituting the cornerstone of treatment aimed at improving long-term patient survival. Owing to their notable advantages, including reduced surgical trauma and accelerated postoperative recovery, minimally invasive procedures are progressively supplanting conventional open surgery and have become the mainstream approach in gastric cancer management. Concurrently, the rapid advancement of artificial intelligence (AI) technologies has enabled real-time intraoperative monitoring of surgical scenes, thereby furnishing novel technical support for adjunctive decision-making, surgical navigation, and skill assessment during gastrectomy. This article provides a systematic review of the current status of AI applications in minimally invasive gastric cancer surgery, with a particular focus on research progress pertaining to instrument recognition, surgical phase identification, delineation of normal anatomical structures, detection of metastatic foci, and early warning of intraoperative adverse events. Furthermore, we discuss the potential value of AI in enhancing surgical efficiency, ensuring patient safety, and optimizing surgical education. On this basis, we further analyze the principal challenges and inherent risks confronting current AI systems, with the aim of informing future technological innovation and facilitating clinical translation.

ObjectiveTo investigate the mechanism by which Feixin decoction treats hypoxic pulmonary hypertension （HPH） by regulating the peroxisome proliferator-activated receptor gamma （PPARγ）/nuclear factor-kappa B （NF-κB）/NOD-like receptor pyrin domain containing 3 （NLRP3） signaling pathway. MethodsForty-eight male SD rats were randomly allocated into normal， hypoxia， and low-， medium- and high-dose （5.85， 11.7， 23.4 g·kg^-1， respectively） Feixin decoction groups， with 8 rats in each group. Except the normal group， the remaining five groups were placed in a hypoxia chamber with an oxygen concentration of （10.0±0.5）% for 8 h per day， 28 days， and administrated with corresponding drugs during the modeling process. After 4 weeks of treatment， echocardiographic parameters ［pulmonary artery acceleration time （PAT）， pulmonary artery ejection time （PET）， right ventricular anterior wall thickness （RVAWd）， and tricuspid annular plane systolic excursion （TAPSE）］ were measured for each group. The right ventricular systolic pressure （RVSP） was measured by the right heart catheterization method， and the right ventricular hypertrophy index （RVHI） was calculated by weighing the heart. The pathological changes in pulmonary arterioles were observed by hematoxylin-eosin staining. The co-localization of α-smooth muscle actin （α-SMA） with NLRP3， N-terminal gasdermin D （N-GSDMD）， and cysteinyl aspartate-specific proteinase-1 （Caspase-1） in pulmonary arteries was detected by immunofluorescence. The protein levels of PPARγ， NF-κB， NLRP3， apoptosis-associated speck-like protein containing a CARD （ASC）， N-GSDMD， interleukin-1β （IL-1β）， interleukin-18（IL-18）， and cleaved Caspase-1 in the lung tissue was determined by Western blot. The ultrastructural changes in pulmonary artery smooth muscle cells （PASMCs） were observed by transmission electron microscopy. ResultsCompared with the normal group， the hypoxia group showed increased RVSP and RVHI （P<0.01）， decreased right heart function （P<0.01）， increased pulmonary vascular remodeling （P<0.01）， increased co-localization of α-SMA with NLRP3， N-GSDMD， and Caspase-1 in pulmonary arterioles （P<0.01）， up-regulated protein levels of NF-κB， NLRP3， ASC， N-GSDMD， IL-1β， IL-18， and cleaved Caspase-1 in the lung tissue （P<0.05， P<0.01）， a down-regulated protein level of PPARγ （P<0.05， P<0.01）， and pyroptosis in PASMCs. Compared with the hypoxia group， Feixin decoction reduced RVSP and RVHI， improved the right heart function and ameliorated pulmonary vascular remodeling （P<0.05， P<0.01）， decreased the co-localization of α-SMA with NLRP3， N-GSDMD， and Caspase-1 （P<0.05， P<0.01）， down-regulated the protein levels of NF-κB， NLRP3， ASC， N-GSDMD， IL-1β， IL-18， and cleaved Caspase-1 in the lung tissue （P<0.05， P<0.01）， up-regulated the protein level of PPARγ （P<0.05， P<0.01）， and alleviated pyroptosis in PASMCs. ConclusionFeixin decoction can ameliorate pulmonary vascular remodeling and right heart dysfunction in chronically induced HPH rats by regulating pyroptosis in PASMCs through the PPARγ/NF-κB/NLRP3 pathway.

ObjectiveTraditional Chinese medicine (TCM) constitutes a valuable cultural heritage and an important source of antitumor compounds. Poria (Poria cocos (Schw.) Wolf), the dried sclerotium of a polyporaceae fungus, was first documented in Shennong’s Classic of Materia Medica and has been used therapeutically and dietarily in China for millennia. Traditionally recognized for its diuretic, spleen-tonifying, and sedative properties, modern pharmacological studies confirm that Poria exhibits antioxidant, anti-inflammatory, antibacterial, and antitumor activities. Pachymic acid (PA; a triterpenoid with the chemical structure 3β-acetyloxy-16α-hydroxy-lanosta-8,24(31)-dien-21-oic acid), isolated from Poria, is a principal bioactive constituent. Emerging evidence indicates PA exerts antitumor effects through multiple mechanisms, though these remain incompletely characterized. Neuroblastoma (NB), a highly malignant pediatric extracranial solid tumor accounting for 15% of childhood cancer deaths, urgently requires safer therapeutics due to the limitations of current treatments. Although PA shows multi-mechanistic antitumor potential, its efficacy against NB remains uncharacterized. This study systematically investigated the potential molecular targets and mechanisms underlying the anti-NB effects of PA by integrating network pharmacology-based target prediction with experimental validation of multi-target interactions through molecular docking, dynamic simulations, and in vitro assays, aimed to establish a novel perspective on PA’s antitumor activity and explore its potential clinical implications for NB treatment by integrating computational predictions with biological assays. MethodsThis study employed network pharmacology to identify potential targets of PA in NB, followed by validation using molecular docking, molecular dynamics (MD) simulations, MM/PBSA free energy analysis, RT-qPCR and Western blot experiments. Network pharmacology analysis included target screening via TCMSP, GeneCards, DisGeNET, SwissTargetPrediction, SuperPred, and PharmMapper. Subsequently, potential targets were predicted by intersecting the results from these databases via Venn analysis. Following target prediction, topological analysis was performed to identify key targets using Cytoscape software. Molecular docking was conducted using AutoDock Vina, with the binding pocket defined based on crystal structures. MD simulations were performed for 100 ns using GROMACS, and RMSD, RMSF, SASA, and hydrogen bonding dynamics were analyzed. MM/PBSA calculations were carried out to estimate the binding free energy of each protein-ligand complex. In vitro validation included RT-qPCR and Western blot, with GAPDH used as an internal control. ResultsThe CCK-8 assay demonstrated a concentration-dependent inhibitory effect of PA on NB cell viability. GO analysis suggested that the anti-NB activity of PA might involve cellular response to chemical stress, vesicle lumen, and protein tyrosine kinase activity. KEGG pathway enrichment analysis suggested that the anti-NB activity of PA might involve the PI3K/AKT, MAPK, and Ras signaling pathways. Molecular docking and MD simulations revealed stable binding interactions between PA and the core target proteins AKT1, EGFR, SRC, and HSP90AA1. RT-qPCR and Western blot analyses further confirmed that PA treatment significantly decreased the mRNA and protein expression of AKT1, EGFR, and SRC while increasing the HSP90AA1 mRNA and protein levels. ConclusionIt was suggested that PA may exert its anti-NB effects by inhibiting AKT1, EGFR, and SRC expression, potentially modulating the PI3K/AKT signaling pathway. These findings provide crucial evidence supporting PA’s development as a therapeutic candidate for NB.

The present study aimed to investigate the effects of eccentric treadmill exercise on adaptive hypertrophy of skeletal muscle in rats. Thirty-two 3-month-old Sprague Dawley (SD) rats were selected and randomly assigned to one of the four groups based on their body weights: 2-week quiet control group (2C), 2-week downhill running exercise group (2E), 4-week quiet control group (4C), and 4-week downhill running exercise group (4E). The downhill running protocol for rats in the exercise groups involved slope of -16°, running speed of 16 m/min, training duration of 90 min, and 5 training sessions per week. Twenty-four hours after the final session of training, all the four groups of rats underwent an exhaustion treadmill exercise. After resting for 48 h, all the rats were euthanized and their gastrocnemius muscles were harvested for analysis. HE staining was used to measure the cross-sectional area (CSA) and diameter of muscle fibers. Transmission electron microscope was used to observe the ultrastructural changes in muscle fibers. Purithromycin surface labeling translation method was used to measure protein synthesis rate. Immunofluorescence double labeling was used to detect the colocalization levels of lysosomal-associated membrane protein 2 (Lamp2)-leucyl-tRNA synthetase (LARS) and Lamp2-mammalian target of rapamycin (mTOR). Western blot was used to measure the protein expression levels of myosin heavy chain (MHC) IIb and LARS, as well as the phosphorylation levels of mTOR, p70 ribosomal protein S6 kinase (p70S6K), and eukaryotic translation initiation factor 4E binding protein 1 (4E-BP1). The results showed that, compared with the 2C group rats, the 2E group rats showed significant increases in wet weight of gastrocnemius muscle, wet weight/body weight ratio, running distance, running time, pre- and post-exercise blood lactate levels, myofibrillar protein content, colocalization levels of Lamp2-LARS and Lamp2-mTOR, and LARS protein expression. Besides these above changes, compared with the 4C group, the 4E group further exhibited significantly increased fiber CSA, fiber diameter, protein synthesis rate, and phosphorylation levels of mTOR, p70S6K, and 4E-BP1. Compared with the quiet control groups, the exercise groups exhibited ultrastructural damage of rat gastrocnemius muscle, which was more pronounced in the 4E group. These findings suggest that eccentric treadmill exercise may promote mTOR translocation to lysosomal membrane, activating mTOR signaling via up-regulating LARS expression. This, in turn, increases protein synthesis rate through the mTOR-p70S6K-4E-BP1 signaling pathway, promoting protein deposition and inducing adaptive skeletal muscle hypertrophy. Although the ultrastructural changes of skeletal muscle are more pronounced, the relatively long training cycles during short-term exercise periods have a more significant effect on promoting gastrocnemius muscle protein synthesis and adaptive hypertrophy.
Animals ; Rats, Sprague-Dawley ; Physical Conditioning, Animal/physiology* ; Rats ; Muscle, Skeletal/metabolism* ; TOR Serine-Threonine Kinases/metabolism* ; Male ; Hypertrophy ; Adaptation, Physiological/physiology* ; Adaptor Proteins, Signal Transducing/metabolism* ; Ribosomal Protein S6 Kinases, 70-kDa/metabolism* ; Intracellular Signaling Peptides and Proteins

Recently, male reproductive health has attracted extensive attention, with the adverse effects of circadian disruption on male fertility gradually gaining recognition. However, the mechanism by which circadian disruption leads to damage to male reproductive system remains unclear. In this review, we first summarized the dual regulatory roles of circadian clock genes on the male reproductive system: (1) circadian regulation of testosterone synthesis via the hypothalamic-pituitary-testicular (HPT) and hypothalamic-pituitary-adrenal (HPA) axes; (2) non-circadian regulation of spermatogenesis. Next, we further listed the possible mechanisms by which circadian disruption impairs male fertility, including interference with the oscillatory function of the reproductive system, i.e., synchronization of the HPT axis, crosstalk between the HPT axis and the HPA axis, as well as direct damage to germ cells by disturbing the non-oscillatory function of the reproductive system. Future research using spatiotemporal omics, epigenomic assays, and neural circuit mapping in studying the male reproductive system may provide new clues to systematically unravel the mechanisms by which circadian disruption affects male reproductive system through circadian clock genes.
Male ; Humans ; Animals ; Circadian Clocks/physiology* ; Hypothalamo-Hypophyseal System/physiology* ; Circadian Rhythm/genetics* ; Spermatogenesis/physiology* ; Pituitary-Adrenal System/physiology* ; Testis/physiology* ; Testosterone/biosynthesis* ; CLOCK Proteins ; Infertility, Male/physiopathology*