ObjectiveTo explore the complex network relationships among pain, kinesiophobia, social participation and knee function in patients after total knee arthroplasty (TKA), and to analyze the moderating effects of different socio-structural factors on the rehabilitation network from an ethical equity perspective. MethodsA convenience sampling method was used to select 291 patients who underwent TKA in Qilu Hospital of Shandong University from May to July, 2023. Pain was assessed using Numerical Rating Scale, kinesiophobia with Chinese short version of the Tampa Scale for Kinesiophobia, social participation with Impact on Participation and Autonomy Questionnaire, and knee function with Hospital for Special Surgery Knee Score. A partial correlation network among pain, kinesiophobia, social participation and knee function was constructed using Graphical Least Absolute Shrinkage and Selection Operator. Key variables were identified through node centrality and bridge centrality analysis. Network Comparison Tests (NCT) were used to analyze network differences among subgroups based on different socio-structural characteristics. ResultsIn the network model, the nodes with the highest strength centrality were indoor participation, activity behavior and activity pain. Bridge centrality analysis indicated that activity pain, knee function, indoor participation and activity cognition were key bridge nodes. NCT revealed no significant differences in overall network structure or global strength among subgroups based on residence, education level or payment method (P > 0.05). However, significant differences in edge weights were found for specific edges such as activity cognition-activity behavior and knee function-indoor participation (P < 0.05). ConclusionThere is a network of interactions among pain, kinesiophobia, social participation and knee function in patients after TKA, with nodes such as indoor participation and activity pain playing key roles in the rehabilitation process. Although the overall rehabilitation network is similar across different socio-structural groups, variations exist in specific relational pathways among patients from rural areas, those with lower education levels, and those with out-of-pocket payment. This suggests that clinical rehabilitation interventions should focus on these core nodes and implement targeted support strategies for socio-structurally disadvantaged groups to promote rehabilitation equity.

ObjectiveTo explore the complex network relationships among pain, kinesiophobia, social participation and knee function in patients after total knee arthroplasty (TKA), and to analyze the moderating effects of different socio-structural factors on the rehabilitation network from an ethical equity perspective. MethodsA convenience sampling method was used to select 291 patients who underwent TKA in Qilu Hospital of Shandong University from May to July, 2023. Pain was assessed using Numerical Rating Scale, kinesiophobia with Chinese short version of the Tampa Scale for Kinesiophobia, social participation with Impact on Participation and Autonomy Questionnaire, and knee function with Hospital for Special Surgery Knee Score. A partial correlation network among pain, kinesiophobia, social participation and knee function was constructed using Graphical Least Absolute Shrinkage and Selection Operator. Key variables were identified through node centrality and bridge centrality analysis. Network Comparison Tests (NCT) were used to analyze network differences among subgroups based on different socio-structural characteristics. ResultsIn the network model, the nodes with the highest strength centrality were indoor participation, activity behavior and activity pain. Bridge centrality analysis indicated that activity pain, knee function, indoor participation and activity cognition were key bridge nodes. NCT revealed no significant differences in overall network structure or global strength among subgroups based on residence, education level or payment method (P > 0.05). However, significant differences in edge weights were found for specific edges such as activity cognition-activity behavior and knee function-indoor participation (P < 0.05). ConclusionThere is a network of interactions among pain, kinesiophobia, social participation and knee function in patients after TKA, with nodes such as indoor participation and activity pain playing key roles in the rehabilitation process. Although the overall rehabilitation network is similar across different socio-structural groups, variations exist in specific relational pathways among patients from rural areas, those with lower education levels, and those with out-of-pocket payment. This suggests that clinical rehabilitation interventions should focus on these core nodes and implement targeted support strategies for socio-structurally disadvantaged groups to promote rehabilitation equity.

ObjectiveTo explore the complex network relationships among pain, kinesiophobia, social participation and knee function in patients after total knee arthroplasty (TKA), and to analyze the moderating effects of different socio-structural factors on the rehabilitation network from an ethical equity perspective. MethodsA convenience sampling method was used to select 291 patients who underwent TKA in Qilu Hospital of Shandong University from May to July, 2023. Pain was assessed using Numerical Rating Scale, kinesiophobia with Chinese short version of the Tampa Scale for Kinesiophobia, social participation with Impact on Participation and Autonomy Questionnaire, and knee function with Hospital for Special Surgery Knee Score. A partial correlation network among pain, kinesiophobia, social participation and knee function was constructed using Graphical Least Absolute Shrinkage and Selection Operator. Key variables were identified through node centrality and bridge centrality analysis. Network Comparison Tests (NCT) were used to analyze network differences among subgroups based on different socio-structural characteristics. ResultsIn the network model, the nodes with the highest strength centrality were indoor participation, activity behavior and activity pain. Bridge centrality analysis indicated that activity pain, knee function, indoor participation and activity cognition were key bridge nodes. NCT revealed no significant differences in overall network structure or global strength among subgroups based on residence, education level or payment method (P > 0.05). However, significant differences in edge weights were found for specific edges such as activity cognition-activity behavior and knee function-indoor participation (P < 0.05). ConclusionThere is a network of interactions among pain, kinesiophobia, social participation and knee function in patients after TKA, with nodes such as indoor participation and activity pain playing key roles in the rehabilitation process. Although the overall rehabilitation network is similar across different socio-structural groups, variations exist in specific relational pathways among patients from rural areas, those with lower education levels, and those with out-of-pocket payment. This suggests that clinical rehabilitation interventions should focus on these core nodes and implement targeted support strategies for socio-structurally disadvantaged groups to promote rehabilitation equity.

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.

This study investigates the genetic diversity and evolutionary relationships of different Artemisia argyi germplasm resources to provide a basis for germplasm identification, variety selection, and resource protection. A total of 192 germplasm resources of A. argyi were studied, and EST-based simple sequence repeat(EST-SSR) primers were designed based on transcriptomic data of A. argyi. Polymerase chain reaction(PCR) amplification was performed on these resources, followed by fluorescence capillary electrophoresis to detect genetic diversity and construct DNA fingerprints. From 197 pairs of primers designed, 28 pairs with polymorphic and clear bands were selected. A total of 278 alleles were detected, with an average of 9.900 0 alleles per primer pair and an average effective number of alleles of 1.407 2. The Shannon's diversity index(I) for the A. argyi germplasm resources ranged from 0.148 1 to 0.418 0, with an average of 0.255 7. The polymorphism information content(PIC) ranged from 0.454 5 to 0.878 0, with an average of 0.766 9, showing high polymorphism. Cluster analysis divided the A. argyi germplasm resources into three major groups: Group Ⅰ contained 136 germplasm samples, Group Ⅱ contained 45, and Group Ⅲ contained 11. Principal component analysis also divided the resources into three groups, which was generally consistent with the clustering results. Mantel test results showed that the genetic variation in A. argyi populations was to some extent influenced by geographic distance, but the effect was minimal. Structure analysis showed that 190 germplasm materials had Q≥ 0.6, indicating that these germplasm materials had a relatively homogeneous genetic origin. Furthermore, 8 core primer pairs were selected from the 28 designed primers, which could distinguish various germplasm types. Using these 8 core primers, DNA fingerprints for the 192 A. argyi germplasm resources were successfully constructed. EST-SSR molecular markers can be used to study the genetic diversity and phylogenetic relationships of A. argyi, providing theoretical support for the identification and molecular-assisted breeding of A. argyi germplasm resources.
Artemisia/classification* ; Microsatellite Repeats ; Genetic Variation ; Expressed Sequence Tags ; DNA Fingerprinting ; Phylogeny ; Polymorphism, Genetic ; DNA, Plant/genetics* ; Genetic Markers

OBJECTIVE: To examine the effectiveness of Chinese medicine (CM) Lianhua Qingwen Granule (LHQW) and Jingyin Gubiao Prescription (JYGB) in asymptomatic or mild patients with Omicron infection in the shelter hospital. METHODS: This single-center retrospective cohort study was conducted in the largest shelter hospital in Shanghai, China, from April 10, 2022 to May 30, 2022. A total of 56,244 asymptomatic and mild Omicron cases were included and divided into 4 groups, i.e., non-administration group (23,702 cases), LHQW group (11,576 cases), JYGB group (12,112 cases), and dual combination of LHQW and JYGB group (8,854 cases). The length of stay (LOS) in the hospital was used to assess the effectiveness of LHQW and JYGB treatment on Omicron infection. RESULTS: Patients aged 41-60 years, with nadir threshold cycle (C_T) value of N gene <25, or those fully vaccinated preferred to receive CM therapy. Before or after propensity score matching (PSM), the multiple linear regression showed that LHQW and JYGB treatment were independent influence factors of LOS (both P<0.001). After PSM, there were significant differences in LOS between the LHQW/JYGB combination and the other groups (P<0.01). The results of factorial design ANOVA proved that the LHQW/JYGB combination therapy synergistically shortened LOS (P=0.032). CONCLUSIONS Patients with a nadir C_T value <25 were more likely to accept CM. The LHQW/JYGB combination therapy could shorten the LOS of Omicron-infected individuals in an isolated environment.
Humans ; Drugs, Chinese Herbal/therapeutic use* ; Male ; Female ; Middle Aged ; Adult ; China/epidemiology* ; Hospitalization ; COVID-19 Drug Treatment ; COVID-19/epidemiology* ; SARS-CoV-2 ; Retrospective Studies ; Treatment Outcome ; Length of Stay ; Young Adult ; Aged

OBJECTIVE: To investigate the effects of astragaloside IV (AS-IV) on podocyte injury of diabetic nephropathy (DN) and reveal its potential mechanism. METHODS: In in vitro experiment, podocytes were divided into 4 groups, normal, high glucose (HG), inositol-requiring enzyme 1 (IRE-1) α activator (HG+thapsigargin 1 µmol/L), and IRE-1α inhibitor (HG+STF-083010, 20 µmol/L) groups. Additionally, podocytes were divided into 4 groups, including normal, HG, AS-IV (HG+AS-IV 20 µmol/L), and IRE-1α inhibitor (HG+STF-083010, 20 µmol/L) groups, respectively. After 24 h treatment, the morphology of podocytes and endoplasmic reticulum (ER) was observed by electron microscopy. The expressions of glucose-regulated protein 78 (GRP78) and IRE-1α were detected by cellular immunofluorescence. In in vivo experiment, DN rat model was established via a consecutive 3-day intraperitoneal streptozotocin (STZ) injections. A total of 40 rats were assigned into the normal, DN, AS-IV [AS-IV 40 mg/(kg·d)], and IRE-1α inhibitor [STF-083010, 10 mg/(kg·d)] groups (n=10), respectively. The general condition, 24-h urine volume, random blood glucose, urinary protein excretion rate (UAER), urea nitrogen (BUN), and serum creatinine (SCr) levels of rats were measured after 8 weeks of intervention. Pathological changes in the renal tissue were observed by hematoxylin and eosin (HE) staining. Quantitative reverse transcription-polymerase chain reaction (RT-PCR) and Western blot were used to detect the expressions of GRP78, IRE-1α, nuclear factor kappa Bp65 (NF-κBp65), interleukin (IL)-1β, NLR family pyrin domain containing 3 (NLRP3), caspase-1, gasdermin D-N (GSDMD-N), and nephrin at the mRNA and protein levels in vivo and in vitro, respectively. RESULTS: Cytoplasmic vacuolation and ER swelling were observed in the HG and IRE-1α activator groups. Podocyte morphology and ER expansion were improved in AS-IV and IRE-1α inhibitor groups compared with HG group. Cellular immunofluorescence showed that compared with the normal group, the fluorescence intensity of GRP78 and IRE-1α in the HG and IRE-1α activator groups were significantly increased whereas decreased in AS-IV and IRE-1α inhibitor groups (P<0.05). Compared with the normal group, the mRNA and protein expressions of GRP78, IRE-1α, NF-κ Bp65, IL-1β, NLRP3, caspase-1 and GSDMD-N in the HG group was increased (P<0.05). Compared with HG group, the expression of above indices was decreased in the AS-IV and IRE-1α inhibitor groups, and the expression in the IRE-1α activator group was increased (P<0.05). The expression of nephrin was decreased in the HG group, and increased in AS-IV and IRE-1α inhibitor groups (P<0.05). The in vivo experiment results revealed that compared to the normal group, the levels of blood glucose, triglyceride, total cholesterol, BUN, blood creatinine and urinary protein in the DN group were higher (P<0.05). Compared with DN group, the above indices in AS-IV and IRE-1α inhibitor groups were decreased (P<0.05). HE staining revealed glomerular hypertrophy, mesangial widening and mesangial cell proliferation in the renal tissue of the DN group. Compared with the DN group, the above pathological changes in renal tissue of AS-IV and IRE-1α inhibitor groups were alleviated. Quantitative RT-PCR and Western blot results of GRP78, IRE-1α, NF-κ Bp65, IL-1β, NLRP3, caspase-1 and GSDMD-N were consistent with immunofluorescence analysis. CONCLUSION AS-IV could reduce ERS and inflammation, improve podocyte pyroptosis, thus exerting a podocyte-protective effect in DN, through regulating IRE-1α/NF-κ B/NLRP3 signaling pathway.
Podocytes/metabolism* ; Animals ; Diabetic Nephropathies/metabolism* ; Saponins/therapeutic use* ; Triterpenes/therapeutic use* ; Signal Transduction/drug effects* ; NF-kappa B/metabolism* ; Protein Serine-Threonine Kinases/metabolism* ; Male ; Rats, Sprague-Dawley ; NLR Family, Pyrin Domain-Containing 3 Protein/metabolism* ; Endoribonucleases/metabolism* ; Endoplasmic Reticulum Chaperone BiP ; Rats ; Diabetes Mellitus, Experimental/complications* ; Endoplasmic Reticulum/metabolism* ; Multienzyme Complexes

OBJECTIVE: To investigate the effect of zedoarondiol on neovascularization of atherosclerotic (AS) plaque by exosomes experiment. METHODS: ApoE^-/- mice were fed with high-fat diet to establish AS model and treated with high- and low-dose (10, 5 mg/kg daily) of zedoarondiol, respectively. After 14 weeks, the expressions of anti-angiogenic protein thrombospondin 1 (THBS-1) and its receptor CD36 in plaques, as well as platelet activation rate and exosome-derived miR-let-7a were detected. Then, zedoarondiol was used to intervene in platelets in vitro, and miR-let-7a was detected in platelet-derived exosomes (Pexo). Finally, human umbilical vein endothelial cells (HUVECs) were transfected with miR-let-7a mimics and treated with Pexo to observe the effect of miR-let-7a in Pexo on tube formation. RESULTS: Animal experiments showed that after treating with zedoarondiol, the neovascularization density in plaques of AS mice was significantly reduced, THBS-1 and CD36 increased, the platelet activation rate was markedly reduced, and the miR-let-7a level in Pexo was reduced (P<0.01). In vitro experiments, the platelet activation rate and miR-let-7a levels in Pexo were significantly reduced after zedoarondiol's intervention. Cell experiments showed that after Pexo's intervention, the tube length increased, and the transfection of miR-let-7a minics further increased the tube length of cells, while reducing the expressions of THBS-1 and CD36. CONCLUSION Zedoarondiol has the effect of inhibiting neovascularization within plaque in AS mice, and its mechanism may be potentially related to inhibiting platelet activation and reducing the Pexo-derived miRNA-let-7a level.
Animals ; MicroRNAs/genetics* ; Exosomes/drug effects* ; Plaque, Atherosclerotic/genetics* ; Neovascularization, Pathologic/genetics* ; Human Umbilical Vein Endothelial Cells/metabolism* ; Humans ; Blood Platelets/drug effects* ; Apolipoproteins E/deficiency* ; Thrombospondin 1/metabolism* ; CD36 Antigens/metabolism* ; Platelet Activation/drug effects* ; Male ; Mice ; Mice, Inbred C57BL

Acute burns and chronic wounds frequently fail to heal owing to various reasons. Most drugs currently used for wound therapy in clinical practice have notable drawbacks, making their application a substantial concern. For instance, anti-inflammatory drugs can exert multisystem toxicity, and cellular therapies are costly and difficult to retain. In recent years, natural functional proteins derived from animals and plants have gained increasing attention owing to their unique biological activities, low cost, and broad application prospects in wound therapy. Herein, we isolated a new protein (JH015Y) from amphibians and demonstrated its excellent wound repair and regeneration properties compared with those of epidermal growth factor, both in vitro and in vivo. JH015 protein increased the proliferative ability of human keratinocytes and skin fibroblasts by 47.73 and 41.40%, respectively. In vivo, the medium-dose (0.5 mg/dose) groups of JH015Y protein demonstrated accelerated wound healing from day 4, with wound healing rates 1.26, 1.27, and 1.14 times that of the blank group in acute wounds, burn wounds, and diabetic ulcer, respectively. Histological analysis of Masson-stained sections indicated that the JH015Y protein contributed to collagen deposition on the wound surface, markedly reduced inflammatory cell infiltration, and exhibited low biological toxicity. Accordingly, the JH015Y protein is a promising biotherapeutic agent for accelerated wound repair and regeneration.