Identification of trypsin-inhibitory site and structure determination of human SPINK2 serine proteinase inhibitor

Methylation and transcriptome analyses construct a prognostic model and reveal the suppressor role of VMO1 in lung adenocarcinoma

BACKGROUND

DNA methylation is an important epigenetic mechanism of gene regulation. The aberrant DNA methylation has been found to play an important role in the initiation and progression of tumors.

RESULTS

Transcriptome and DNA methylation data of lung adenocarcinoma (LUAD) patients were co-analyzed and 95 methylation-driven genes (MDGs) was found in relation to LUAD. A prognostic model based on 3 MDGs (GMNN, SPINK2 and VMO1) was constructed by Univariate and Multivariate cox regression analyses. The risk score generated from the prognostic model could be used to classify LUAD patients into high and low risk groups. Furthermore, it was found that the risk score was associated with tumor microenvironment (TME) and clinical characteristics (survival status and T stage) of patients. Interestingly, we identified and validated that the patients in the low-risk group responded better to immunotherapy treatment. Then, a nomogram model based on the risk score and clinical characteristics was established which showed significant prediction value. The down-regulation and hypermethylation levels of vitelline membrane outer layer protein 1 homolog (VMO1) were verified in paired LUAD tumor and non-tumor tissues by pyrosequencing assay and RT-qPCR. Furthermore, MTT, migration and wound healing assays were performed with lentivector-mediated ectopic over-expression and 5-Aza-dC demethylation followed by siRNA rescue experiments to investigate the role of VMO1 in LUAD cells. Our results indicated that VMO1 could inhibit proliferation and migration of A549 and NCI-H1299 cells.

CONCLUSIONS

In summary, our experiments constructed a prognostic model with high capacity for risk prediction in LUAD patients. VMO1 had a malignant suppressor role in LUAD cells. The correlation between risk score and TME might elucidate a potential mechanism of oncogenesis and provide an avenue for further therapeutic targets.

RARRES1 inhibits hepatocellular carcinoma progression and increases its sensitivity to lenvatinib through interaction with SPINK2

BACKGROUND

Lenvatinib is an oral small molecule inhibitor approved for treating patients with unresectable hepatocellular carcinoma (HCC) worldwide. Increasing cell sensitivity to lenvatinib would be an effective method of improving therapeutic efficacy.

METHODS

High throughput methods was used to scan the differentially expressed genes (DEGs) related to lenvatinib sensitivity in HCC cells. Gain- and loss-function experiments were used to explore the functions of these DEGs in HCC and lenvatinib sensitivity. CO-IP assay and rescue experiments were utilized to investigate the mechanism.

RESULTS

We identified that RAR responder protein 1 (RARRES1), a podocyte-specific growth arrest gene, was among significantly upregulated DEGs in HCC cells following lenvatinib treatment. Functional analysis showed that ectopic RARRES1 expression decreased HCC progression in vitro and in vivo, as well as improving tumor sensitivity to lenvatinib, while RARRES1 silencing increased HCC cell proliferation and migration. Mechanistically, co-immunoprecipitation assays demonstrated that RARRES1 interacted with serine protease inhibitor Kazal-type 2 (SPINK2) in HCC cells. Further, SPINK2 overexpression suppressed HCC cell proliferation and migration, as well as increasing sensitivity to lenvatinib whereas SPINK2 knockdown promoted cell progression and decreased lenvatinib sensitivity. The mRNA and protein levels of RARRES1 and SPINK2 were low in HCC tissue samples, relative to those in normal liver tissue.

CONCLUSIONS

Our findings highlighted that RARRES1 can inhibit HCC progression and regulate HCC sensitivity to lenvatinib by interacting SPINK2, representing a new tumor suppressor RARRES1/SPINK2 axis in HCC that modulates sensitivity to lenvatinib.

Temporary serine protease inhibition and the role of SPINK2 in human bone marrow

Protease temporary inhibitors are true substrates that bind the catalytic site with high affinity but are slowly degraded, thus acting as inhibitor for a defined time window. Serine peptidase inhibitor Kazal type (SPINK) family is endowed with such functional property whose physiological meaning is poorly explored. High expression of SPINK2 in some hematopoietic malignancies prompted us to investigate its role in adult human bone marrow. We report here the physiological expression of SPINK2 in hematopoietic stem and progenitor cells (HSPCs) and mobilized cluster differentiation 34 (CD34)⁺ cells. We determined the SPINK2 degradation constant and derived a mathematical relationship predicting the zone of inhibited target protease activity surrounding the SPINK2-secreting HSPCs. Analysis of putative target proteases for SPINK2 revealed the expression of PRSS2 and PRSS57 in HSPCs. Our combined results suggest that SPINK2 and its target serine proteases might play a role in the intercellular communication within the hematopoietic stem cell niche.

SPINK2 Protein Expression Is an Independent Adverse Prognostic Marker in AML and Is Potentially Implicated in the Regulation of Ferroptosis and Immune Response

There is an urgent need for the identification as well as clinicopathological and functional characterization of potent prognostic biomarkers and therapeutic targets in acute myeloid leukemia (AML). Using immunohistochemistry and next-generation sequencing, we investigated the protein expression as well as clinicopathological and prognostic associations of serine protease inhibitor Kazal type 2 (SPINK2) in AML and examined its potential biological functions. High SPINK2 protein expression was an independent adverse biomarker for survival and an indicator of elevated therapy resistance and relapse risk. SPINK2 expression was associated with AML with an NPM1 mutation and an intermediate risk by cytogenetics and European LeukemiaNet (ELN) 2022 criteria. Furthermore, SPINK2 expression could refine the ELN2022prognostic stratification. Functionally, an RNA sequencing analysis uncovered a potential link of SPINK2 with ferroptosis and immune response. SPINK2 regulated the expression of certain P53 targets and ferroptosis-related genes, including SLC7A11 and STEAP3, and affected cystine uptake, intracellular iron levels and sensitivity to erastin, a specific ferroptosis inducer. Furthermore, SPINK2 inhibition consistently increased the expression of ALCAM, an immune response enhancer and promoter of T-cell activity. Additionally, we identified a potential small-molecule inhibitor of SPINK2, which requires further characterization. In summary, high SPINK2 protein expression was a potent adverse prognostic marker in AML and might represent a druggable target.

SPINKs in Tumors: Potential Therapeutic Targets

The serine protease inhibitor Kazal type (SPINK) family includes SPINK1-14 and is the largest branch in the serine protease inhibitor family. SPINKs play an important role in pancreatic physiology and disease, sperm maturation and capacitation, Nager syndrome, inflammation and the skin barrier. Evidence shows that the unregulated expression of SPINK1, 2, 4, 5, 6, 7, and 13 is closely related to human tumors. Different SPINKs exhibit various regulatory modes in different tumors and can be used as tumor prognostic markers. This article reviews the role of SPINK1, 2, 4, 5, 6, 7, and 13 in different human cancer processes and helps to identify new cancer treatment targets.

NUP-98 Rearrangements Led to the Identification of Candidate Biomarkers for Primary Induction Failure in Pediatric Acute Myeloid Leukemia

Conventional chemotherapy for acute myeloid leukemia regimens generally encompass an intensive induction phase, in order to achieve a morphological remission in terms of bone marrow blasts (<5%). The majority of cases are classified as Primary Induction Response (PIR); unfortunately, 15% of children do not achieve remission and are defined Primary Induction Failure (PIF). This study aims to characterize the gene expression profile of PIF in children with Acute Myeloid Leukemia (AML), in order to detect molecular pathways dysfunctions and identify potential biomarkers. Given that NUP98-rearrangements are enriched in PIF-AML patients, we investigated the association of NUP98-driven genes in primary chemoresistance. Therefore, 85 expression arrays, deposited on GEO database, and 358 RNAseq AML samples, from TARGET program, were analyzed for “Differentially Expressed Genes” (DEGs) between NUP98+ and NUP98-, identifying 110 highly confident NUP98/PIF-associated DEGs. We confirmed, by qRT-PCR, the overexpression of nine DEGs, selected on the bases of the diagnostic accuracy, in a local cohort of PIF patients: SPINK2, TMA7, SPCS2, CDCP1, CAPZA1, FGFR1OP2, MAN1A2, NT5C3A and SRP54. In conclusion, the integrated analysis of NUP98 mutational analysis and transcriptome profiles allowed the identification of novel putative biomarkers for the prediction of PIF in AML.

Identification of the hub genes and pathways involved in acute myeloid leukemia using bioinformatics analysis

BACKGROUND

We identified the hub genes and pathways dysregulated in acute myeloid leukemia and the potential molecular mechanisms involved.

METHODS

We downloaded the GSE15061 gene expression dataset from the Gene Expression Omnibus database and used weighted gene co-expression network analysis to identify hub genes. Differential expression of the genes was evaluated using the limma package in R software. Subsequently, we built a protein-protein interaction network followed by functional enrichment analysis. Then, the prognostic significance of gene expression was explored in terms of overall survival. Finally, transcription factor-mRNA (ribonucleic acid) and microRNA-mRNA interaction analysis was also explored.

RESULTS

We identified 100 differentially expressed hub genes. Functional enrichment analysis indicated that the genes were principally involved in immune system regulation, host defense, and negative regulation of apoptosis and myeloid cell differentiation. We identified 4 hub genes, the expression of which was significantly correlated with overall survival. Finally, 26 key regulators for hub genes and 38 microRNA-mRNA interactions were identified.

CONCLUSION

We performed a comprehensive bioinformatics analysis of hub genes potentially involved in acute myeloid leukemia development. Further molecular biological experiments are required to confirm the roles played by these genes.

A six-gene leukemic stem cell score identifies high risk pediatric acute myeloid leukemia

Recently, mRNA-expression signature enriched in LSCs was used to create a 17-gene leukemic stem cell (LSC17) score predictive of prognosis in adult AML. By fitting a Cox-LASSO regression model to the clinical outcome and gene-expression levels of LSC enriched genes in 163 pediatric participants of the AML02 multi-center clinical trial (NCT00136084), we developed a 6-gene LSC score of prognostic value in pediatric AML (pLSC6). In the AML02 cohort, the 5-year event-free survival (EFS) of patients within low-pLSC6 group (n=97) was 78.3 (95% CI=70.5–86.9%) as compared to 34.5(95% CI=24.7–48.2 %) in patients within high-pLSC6 group (n=66 subjects), p<0.00001. pLSC6 remained significantly associated with EFS and overall survival (OS) after adjusting for induction 1-MRD status, risk-group, FLT3-status, WBC-count at diagnosis and age. pLSC6 formula developed in the AML02 cohort was validated in the pediatric AML-TARGET project data (n=205), confirming its prognostic value in both single-predictor and multiple-predictor Cox regression models. In both cohorts, pLSC6 predicted outcome of transplant patients, suggesting it as a useful criterion for transplant referrals. Our results suggest that pLSC6 score holds promise in redefining initial risk-stratification and identifying poor risk AML thereby providing guidance for developing novel treatment strategies.

A protein scaffold, engineered SPINK2, for generation of inhibitors with high affinity and specificity against target proteases

Proteases are one of attractive therapeutic targets to play key roles in pharmacological action. There are many protease inhibitors in nature, and most of them structurally have cystine knot motifs. Their structures are favorable for recognition of active pockets of proteases, leading to the potent inhibition. However, they also have drawbacks, such as broad cross-reactivity, on the therapeutic application. To create therapeutic proteins derived from a disulfide-rich scaffold, we selected human serine protease inhibitor Kazal type 2 (SPINK2) through a scaffold screening, as a protein scaffold with requirements for therapeutic proteins. We then constructed a diverse library of the engineered SPINK2 by introducing random mutations into its flexible loop region with the designed method. By phage panning against four serine proteases, we isolated potent inhibitors against each target with picomolar K_D and sub-nanomolar K_i values. Also, they exhibited the desired specificities against target proteases without inhibiting non-target proteases. The crystal structure of kallikrein related peptidase 4 (KLK4)-engineered SPINK2 complex revealed the interface with extensive conformational complementarity. Our study demonstrates that engineered SPINK2 can serve as a scaffold to generate therapeutic molecules against target proteins with groove structures.

Elevated SPINK2 gene expression is a predictor of poor prognosis in acute myeloid leukemia

Acute myeloid leukemia (AML) has a high mortality rate and its clinical management remains challenging. The aim of the present study was to identify the hub genes involved in AML. In order to do so, the gene expression data of the GSE9476 database, including 26 AML and 10 normal samples, were downloaded from the Gene Expression Omnibus database. Differentially expressed genes (DEGs) were then identified via bioinformatics analysis. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway analyses were performed on DEGs. Furthermore, the most upregulated genes were selected for further investigation in the Oncomine, gene expression profiling interactive analysis and UALCAN datasets. In total, 1,744 upregulated and 1,956 downregulated genes were detected. The GO and KEGG results revealed that upregulated genes were enriched in metabolic processes, while downregulated genes were associated with the immune response. Serine protease inhibitor Kazal-type 2 (SPINK2) ranked first among all the upregulated genes and was regarded as a hub gene in the development of AML. The overexpression of SPINK2 was validated in 12 patients with AML from the Linyi Central Hospital and in data from the Oncomine and Gene Expression Profiling Interactive Analysis (GEPIA) databases. Furthermore, the UALCAN and GEPIA datasets demonstrated that patients with high SPINK2 levels had shorter survival times. In conclusion, the results from the present study revealed that the SPINK2 gene was upregulated in patients with AML and that elevated SPINK2 expression was associated with poor outcomes in these patients.

Comparative expression profiling of testis-enriched genes regulated during the development of spermatogonial cells

The testis has been identified as the organ in which a large number of tissue-enriched genes are present. However, a large portion of transcripts related to each stage or cell type in the testis still remains unknown. In this study, databases combined with confirmatory measurements were used to investigate testis-enriched genes, localization in the testis, developmental regulation, gene expression profiles of testicular disease, and signaling pathways. Our comparative analysis of GEO DataSets showed that 24 genes are predominantly expressed in testis. Cellular locations of 15 testis-enriched proteins in human testis have been identified and most of them were located in spermatocytes and round spermatids. Real-time PCR revealed that expressions of these 15 genes are significantly increased during testis development. Also, an analysis of GEO DataSets indicated that expressions of these 15 genes were significantly decreased in teratozoospermic patients and polyubiquitin knockout mice, suggesting their involvement in normal testis development. Pathway analysis revealed that most of those 15 genes are implicated in various sperm-related cell processes and disease conditions. This approach provides effective strategies for discovering novel testis-enriched genes and their expression patterns, paving the way for future characterization of their functions regarding infertility and providing new biomarkers for specific stages of spematogenesis.

The solution structure of the kallikrein-related peptidases inhibitor SPINK6

Kallikrein-related peptidases (KLKs) are crucial for epidermal barrier function and are involved in the proteolytic regulation of the desquamation process. Elevated KLK levels were reported in atopic dermatitis. In skin, the proteolytic activity of KLKs is regulated by specific inhibitors of the serine protease inhibitor of Kazal-type (SPINK) family. SPINK6 was shown to be expressed in human stratum corneum and is able to inhibit several KLKs such as KLK4, -5, -12, -13 and -14. In order to understand the structural traits of the specific inhibition we solved the structure of SPINK6 in solution by NMR-spectroscopy and studied its interaction with KLKs. Thereby, beside the conserved binding mode, we identified an alternate binding mode which has so far not been observed for SPINK inhibitors.

Gene Expression Profiling in Peripheral Blood Cells and Synovial Membranes of Patients with Psoriatic Arthritis

Background

Psoriatic arthritis (PsA) is an inflammatory arthritis whose pathogenesis is poorly understood; it is characterized by bone erosions and new bone formation. The diagnosis of PsA is mainly clinical and diagnostic biomarkers are not yet available. The aim of this work was to clarify some aspects of the disease pathogenesis and to identify specific gene signatures in paired peripheral blood cells (PBC) and synovial biopsies of patients with PsA. Moreover, we tried to identify biomarkers that can be used in clinical practice.

Methods

PBC and synovial biopsies of 10 patients with PsA were used to study gene expression using Affymetrix arrays. The expression values were validated by Q-PCR, FACS analysis and by the detection of soluble mediators.

Results

Synovial biopsies of patients showed a modulation of approximately 200 genes when compared to the biopsies of healthy donors. Among the differentially expressed genes we observed the upregulation of Th17 related genes and of type I interferon (IFN) inducible genes. FACS analysis confirmed the Th17 polarization. Moreover, the synovial trascriptome shows gene clusters (bone remodeling, angiogenesis and inflammation) involved in the pathogenesis of PsA. Interestingly 90 genes are modulated in both compartments (PBC and synovium) suggesting that signature pathways in PBC mirror those of the inflamed synovium. Finally the osteoactivin gene was upregulared in both PBC and synovial biopsies and this finding was confirmed by the detection of high levels of osteoactivin in PsA sera but not in other inflammatory arthritides.

Conclusions

We describe the first analysis of the trancriptome in paired synovial tissue and PBC of patients with PsA. This study strengthens the hypothesis that PsA is of autoimmune origin since the coactivity of IFN and Th17 pathways is typical of autoimmunity. Finally these findings have allowed the identification of a possible disease biomarker, osteoactivin, easily detectable in PsA serum.

Spink2 modulates apoptotic susceptibility and is a candidate gene in the Rgcs1 QTL that affects retinal ganglion cell death after optic nerve damage

The Rgcs1 quantitative trait locus, on mouse chromosome 5, influences susceptibility of retinal ganglion cells to acute damage of the optic nerve. Normally resistant mice (DBA/2J) congenic for the susceptible allele from BALB/cByJ mice exhibit susceptibility to ganglion cells, not only in acute optic nerve crush, but also to chronic inherited glaucoma that is characteristic of the DBA/2J strain as they age. SNP mapping of this QTL has narrowed the region of interest to 1 Mb. In this region, a single gene (Spink2) is the most likely candidate for this effect. Spink2 is expressed in retinal ganglion cells and is increased after optic nerve damage. This gene is also polymorphic between resistant and susceptible strains, containing a single conserved amino acid change (threonine to serine) and a 220 bp deletion in intron 1 that may quantitatively alter endogenous expression levels between strains. Overexpression of the different variants of Spink2 in D407 tissue culture cells also increases their susceptibility to the apoptosis-inducing agent staurosporine in a manner consistent with the differential susceptibility between the DBA/2J and BALB/cByJ strains.

Bioinformatic analysis of chemokine (C-C motif) ligand 21 and SPARC-like protein 1 revealing their associations with drug resistance in ovarian cancer

Chemokine (C-C motif) ligand 21 (CCL21) and SPARC-like protein 1 (SPARCL1/MAST9/hevin/SC-1) are associated with various biological behavior in the development of cancers. Although the expression of CCL21 and SPARCL1 is downregulated in many solid tumors, their roles in ovarian cancer and their associations with drug resistance have rarely been studied. We performed a comprehensive bioinformatic analysis consisting of motif analysis, literature co-occurrence, gene/protein-gene/protein interaction network, protein-small molecule interaction network, and microRNAs enrichments which revealed that CCL21 and SPARCL1 directly or indirectly interact with a number of genes, proteins, small molecules and pathways associated with drug resistance in ovarian and other cancers. These results suggested that CCL21 and SPARCL1 may contribute to drug resistance in ovarian cancer. This study provided important information for further investigation of drug resistance-related functions of CCL21 and SPARCL1 in ovarian cancer.

Detection of protease activities by flash chronopotentiometry using a reversible polycation-sensitive polymeric membrane electrode

A novel electrochemical method, termed flash chronopotentiometry (FCP), is used to develop a rapid and sensitive method for detecting protease activities. In this method, an appropriate current pulse is applied across a polycation-selective polymer membrane to induce a strong flux of the polycationic peptides from the sample phase into the organic membrane of the electrode. During this current pulse, the cell potential (EMF) is monitored continuously, and is a function of the polypeptide concentration. The imposed current causes a local depletion of the polypeptide at the sample/membrane interface, which yields a drastic potential change in the observed chronopotentiogram at a characteristic time, called the transition time (τ). For a given magnitude of current, the square root of τ is directly proportional to the concentration of the polypeptide. Proteases cleave polypeptides into smaller fragments that are not favorably extracted into the membrane of the sensor. Therefore, a decrease in the transition time is observed during the proteolysis process. The degree of change in the transition time can be correlated to protease activity. To demonstrate this approach, the activities of trypsin and α-chymotrypsin are detected using protamine and synthetic polycationic oligopeptides that possess specific cleavage sites that are recognized by these proteases.