Next Article in Journal
Homo sapiens May Incorporate Daily Acute Cycles of “Conditioning–Deconditioning” to Maintain Musculoskeletal Integrity: Need to Integrate with Biological Clocks and Circadian Rhythm Mediators
Next Article in Special Issue
Functional Mimicry of Eukaryotic Actin Assembly by Pathogen Effector Proteins
Previous Article in Journal
Targeting the Mitochondrial Protein VDAC1 as a Potential Therapeutic Strategy in ALS
Previous Article in Special Issue
Photorhabdus luminescens TccC3 Toxin Targets the Dynamic Population of F-Actin and Impairs Cell Cortex Integrity
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Clostridium novyi’s Alpha-Toxin Changes Proteome and Phosphoproteome of HEp-2 Cells

1
Institute of Toxicology, Hannover Medical School, Carl-Neuberg-Str. 1, 30625 Hannover, Germany
2
Core Facility Proteomics, Hannover Medical School, Carl-Neuberg-Str. 1, 30625 Hannover, Germany
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2022, 23(17), 9939; https://doi.org/10.3390/ijms23179939
Submission received: 20 July 2022 / Revised: 17 August 2022 / Accepted: 27 August 2022 / Published: 1 September 2022

Abstract

:
C. novyi type A produces the alpha-toxin (TcnA) that belongs to the large clostridial glucosylating toxins (LCGTs) and is able to modify small GTPases by N-acetylglucosamination on conserved threonine residues. In contrast, other LCGTs including Clostridioides difficile toxin A and toxin B (TcdA; TcdB) modify small GTPases by mono-o-glucosylation. Both modifications inactivate the GTPases and cause strong effects on GTPase-dependent signal transduction pathways and the consequent reorganization of the actin cytoskeleton leading to cell rounding and finally cell death. However, the effect of TcnA on target cells is largely unexplored. Therefore, we performed a comprehensive screening approach of TcnA treated HEp-2 cells and analyzed their proteome and their phosphoproteome using LC-MS-based methods. With this data-dependent acquisition (DDA) approach, 5086 proteins and 9427 phosphosites could be identified and quantified. Of these, 35 proteins were found to be significantly altered after toxin treatment, and 1832 phosphosites were responsive to TcnA treatment. By analyzing the TcnA-induced proteomic effects of HEp-2 cells, 23 common signaling pathways were identified to be altered, including Actin Cytoskeleton Signaling, Epithelial Adherens Junction Signaling, and Signaling by Rho Family GTPases. All these pathways are also regulated after application of TcdA or TcdB of C. difficile. After TcnA treatment the regulation on phosphorylation level was much stronger compared to the proteome level, in terms of both strength of regulation and the number of regulated phosphosites. Interestingly, various signaling pathways such as Signaling by Rho Family GTPases or Integrin Signaling were activated on proteome level while being inhibited on phosphorylation level or vice versa as observed for the Role of BRCA1 in DNA Damage Response. ZIP kinase, as well as Calmodulin-dependent protein kinases IV & II, were observed as activated while Aurora-A kinase and CDK kinases tended to be inhibited in cells treated with TcnA based on their substrate regulation pattern.

1. Introduction

Clostridium novyi is an obligate anaerobic pathogen that causes serious disease in both animals and humans. It is the main cause of gas gangrene in humans, while it is involved in infectious necrotic hepatitis in animals [1,2]. The exact identification of C. novyi and its contribution to the course of both diseases is often difficult due to mixed clostridial populations of Clostridium perfringens, C. novyi, Clostridium septicum, Clostridium hystolyticum, or others causing these diseases [1]. Pathogenic C. novyi strains are classified by different toxin productions. Type A strains produce alpha, gamma, delta, and epsilon toxins [3,4,5]. Type B produces alpha, beta, and zeta toxins, and C. haemolyticum, also known as C. novyi type C is capable of producing beta, eta, and theta toxins [3,4,5,6,7]. C. novyi type D does not produce any toxin and is generally referred to as nonvirulent.
Alpha Toxin (TcnA) is one of the major virulence factors of C. novyi and belongs to the family of large clostridial glucosylating toxins (LCGTs) [8]. Other LCGTs are the C. difficile toxins A and B (TcdA, TcdB), the lethal and hemorrhagic toxin of Clostridium sordellii (TcsH, TcsL), and the large cytotoxin of C. perfringens (TpeL) [9]. LCGTs exhibit a four-domain structure. The receptor-binding domain is located at the c-terminus and mediates the initial binding of the toxin to the cell surface. Receptor-mediated endocytosis follows, in which the translocation domain, which is adjacent to the receptor-binding domain, enables toxin translocation through the membrane of the acidified endosome into the cytosol. The N-terminal glucosyltransferase domain is cleaved off by the cysteine protease domain and exerts its enzymatic activity in the cytosol [10,11].
LCGTs inactivate small GTPases by glycosylating conserved threonine residues (Thr-35 in Rac/Thr-37 in RhoA) [12]. Most LCGTs exploit UDP-glucose as the sugar donor for this process, which leads to mono-O-glucosylation. However, TcnA preferably uses UDP-N-acetylglucosamine (UDP-GlcNac), which leads to Mono-O-N-acetylglucosamination of the GTPases [12,13,14]. Crucial for co-substrate specificity of LCGTs are the amino acids at positions 383 and 385 in the catalytic cleft. For UDP-glucose favoring toxins such as TcdB, those positions contain Ile-383 and Glu-385 while TcnA favoring UDP-GlcNac comprises Ser-383 and Ala-385. Isoleucine and glutamic acid contain bulkier side chains than serine and alanine, respectively, thus limiting the available space for the co-substrate [15]. By alteration of GTPase signaling, there is also an impact on the actin cytoskeleton organization [16]. Data on the consequences of TcnA effects on GTPases are hardly available while in contrast, TcdA and TcdB have been extensively studied.
Herein we describe the effects of TcnA on an epithelial cell line called HEp-2. It has been established via HeLa cell contamination and expresses LDLR as well as SLC35B2, which are key factors in initial cellular binding and uptake of TcnA [17]. We performed a global DDA proteome analysis [18] to screen for TcnA responsive proteins and phosphosites after treatment with TcnA. The phosphoproteome responds much faster to altered conditions, e.g., addition of TncA, and should lead to the identification of important regulatory proteins. Phosphorylation is a post-translational modification that proceeds much faster than the synthesis or degradation of mature proteins. With the reproducible enrichment and purification of phosphoproteins and peptides using metal affinity strategies, it has become possible to analyze the phosphoproteome in a comprehensive way [19,20]. Additionally, we compared the impact of TcnA on the proteome and phosphoproteome of target cells with the corresponding effects of TcdA and TcdB [21,22].

2. Results

2.1. Morphological Alterations of HEp-2 Cells Treated with C. novyi’s Alpha Toxin

Morphological changes of HEp-2 cells were analyzed using phase-contrast microscopy (Figure 1). Three different treatments were compared: no toxin (Ctrl), TcnA inactivated with formaldehyde (FA_TcnA) and active TcnA (TcnA). TcnA concentration of 255 ng/mL were used. First morphological changes were visible after 3–4 h (data not shown) and after 24 h most of the cells were completely round (Figure 1), except those treated with no toxin and FA_TcnA. Along with cell rounding, cell adhesion to the flask was reduced and led to cell strands ranging into the medium, indicated by white arrows in Figure 1. Increasing TcnA concentrations or incubation times did not lead to a higher percentage of completely rounded cells without being accompanied by undesirable side effects such as apoptotic or necrotic cells (data not shown). Cell morphology of FA_TcnA looked similar to control cells.

2.2. Generation of a Control Toxin to Elucidate TcnA Effects on Proteome and Phosphoproteome

To detect the effects of TcnA on cellular protein homeostasis using MS-based proteomics, an appropriate control needs to be used. Since no catalytically inactive mutant was available for this purpose (as for TcdA and TcdB), the toxin was inactivated with formaldehyde and subsequently tested on HEp-2 cells. Effects of this inactivated toxin were analyzed on proteome and phosphoproteome level in comparison to the active toxin and non-treated control cells.
For the comparison on proteome level, data from 5086 identified and quantified protein groups were used that have been obtained from the three treatment conditions TcnA, FA_TcnA, and control (proteome data see Supplementary Table S1). In a Benjamini Hochberg (BJH; FDR = 0.05) [23] corrected analysis of variance (ANOVA) test 530 proteins showed significant differences in their abundance across the three treatments. Hierarchical clustering with BJH positive proteins was performed after z-scoring (Figure 2A). FA_TcnA and Ctrl samples showed similar regulated protein clusters, whereas a significant difference was observable in comparison to the TcnA treated cells. Only one protein (Zinc Finger Protein 611, Znf611, see Supplementary Figure S1) was detected as significantly down-regulated between control and FA_TcnA treatment, showing that TcnA inactivation by formaldehyde was successful and the influence of potentially co-enriched other factors during toxin purification on HEp-2 cells should be unlikely. Znf611 was excluded from further proteome analysis of protein regulation between TcnA and FA_TcnA treated replicates. Further analyses of detailed TcnA-induced alterations concerning the proteome are presented below.
To further control the suitability of the FA-inactivated TcnA a comparison of phosphoproteome was performed. Overall, 9427 phosphosites belonging to 2125 phosphoproteins that had been normalized to the corresponding protein abundance determined in the proteome measurements were used. This way, direct comparison of active and inactivated TcnA was possible (phosphoproteome data see Supplementary Table S2).
BJH (FDR = 0.05) corrected ANOVA testing identified 430 phosphosites as significantly regulated across the three treatments. Hierarchical clustering with BJH-positive phosphosites revealed that FA_TcnA-treated cells and control cells shared more similarities than TcnA treated cells (Figure 2B). However, within FA_TcnA treated cells eight phosphosites showed a significant down-regulation while four phosphosites were significantly up-regulated (Supplementary Figure S1). These twelve phosphosites were excluded from further bioinformatic analysis of regulated phosphosites in TcnA versus FA_TcnA treated cells. Overall, most of the alterations induced by FA_TcnA clustered with control conditions, both on proteome and phosphoproteome level.
Thus, the formaldehyde inactivated TcnA has been established as a valuable control to study the effects of TcnA on HEp-2 cells.

2.3. Proteomic Effects Induced by TcnA in HEp-2 Cells

The impact of TcnA on target cells was compared to FA_TcnA, which was used as control. The comprehensive proteome analysis was based on 5086 identified and quantified protein groups using a TMT approach (Supplementary Table S1). In TcnA-treated cells, 14 proteins were significantly up-regulated, while 21 proteins were significantly down-regulated (Figure 3A). The log2 expression differences ranged from −2.34 to 2.58 between treated and untreated samples. Among the strongest up-regulated proteins were Legumain (Lgmn) and Immunoglobulin-like and fibronectin type III domain-containing protein 1 (Igfn1), while proteins like Troponin I (Tnni2) as well as ATP-binding cassette sub-family D member 2 (Abcd2) were identified to be strongly down-regulated.
A BJH (FDR = 0.05) corrected t-test identified 1021 positive proteins which were used for Fisher’s exact test to identify TcnA-responsive gene ontology (GO) groups (Figure 3B). GO annotations based on three main classes, cellular component (GOCC), molecular function (GOMF) and biological process (GOBP) were the basis for this Fisher’s exact test. There are various subterms within the three GO main classes, which were annotated to the proteins, based on their localization, known biological or molecular functions. With the Fisher’s exact test, a comparison between the frequency of these subterms within the whole proteome and across the significantly regulated proteins was made. Most of the significantly enriched terms within the subset belonged to GOCC. The only identified GOMF term was cytoskeletal protein binding, while no terms concerning biological processes were enriched. It has to be mentioned that the annotation follows a hierarchical order. Therefore, all proteins with the annotation actin cytoskeleton were also annotated as cytoskeleton. The same applies to the term actin filament bundle which was always accompanied by the terms stress fibers and actomyosin being the least stringent term.
For canonical pathway analysis via IPATM all changed proteins with a significant p-value < 0.05 were used. IPATM uses p-value and z-score as statistical parameters to describe the quality of the predicted events. The p-value includes the amount of genes participating in a pathway (provided by IPATM) and the overlap of those genes within the uploaded data. With a p-value < 0.05 a non-random association between pathway and data is present [24]. The z-score additionally describes the observed regulation status. A z-score below zero indicates inactivation of the pathway. The lower the value, the better the measured protein regulations agree with the inactivation of the pathway. A z-score above 0 shows activation of the pathway and, comprehensibly, the higher the value, the better the measured protein regulations agree with the activation of the pathway [25]. Some of the most significantly regulated pathways are illustrated (Figure 3C). The strongest activated pathways in TcnA treated cells were Signaling by Rho Family GTPases and RHOA Signaling, which contained a big overlap of proteins included in both. Of the significantly regulated proteins 37 were identified as part of RhoA signaling pathway, while 57 proteins were assigned to Rho Family GTPases pathway. An overlap of 29 proteins was present in both. The pathways RHOGDI Signaling, Role of BRCA1 in DNA Damage Response and IGF-1 Signaling were inhibited, while every other was predicted as at least slightly activated according to the z-score.
Comparison of IPATM canonical pathways altered by TcdA [21], TcdB [22], and TcnA in Hep-2 cells was performed in order to identify co- and differentially regulated cellular processes. It has to been mentioned that TcdA and TcdB analyses were done by stable isotope labeling with amino acids in cell culture (SILAC) approaches. Those were characterized by high reproducibility but lower degrees of protein regulation as well as overall lower numbers of quantifiable proteins than the TcnA TMT approach. This is why there were only 169 proteins passing significance level of p < 0.05 for TcdA and 168 proteins for TcdB while the TcnA analysis resulted in 2064 proteins with p < 0.05 used for the IPA analysis. For TcdA 39 regulated (z-score > 0; <0; p < 0.05) pathways could be identified. In average, the pathways tended to be inactivated after TcdA treatment. In contrast, for TcdB 118 regulated pathways were observed, most of which were activated. Following TcnA treatment, 162 regulated pathways were identified. In Figure 4, some of the pathways regulated by every LCGT treatment are depicted (for regulated pathways including z-scores see Table 1, all identified pathways are listed in Supplementary Table S3). The overlap between TcdB and TcnA treated cells is highest in IPA canonical pathway analysis, which is due to the fact that only 39 regulated pathways were identified with TcdA treatment. However, looking directly on the proteome level the overlap between TcdA and TcdB treatment is higher. A comparison of the top 20 regulated proteins per treatment is shown in Table 2. The proteins RhoB and Hmox1 showed pronounced up-regulation within TcdA and TcdB treatment, while Yap1 and Tpm2 are both down regulated. No overlap between any of the strongest regulated proteins after TcnA treatment to treatment with one of the other toxins was observable.

2.4. Phosphoproteomic Effects of TcnA on HEp-2 Cells

Phosphoproteomic analysis of TcnA treated HEp-2 cells revealed overall, 9427 identified phosphosites (Supplementary Table S2) with log2 difference between TcnA and FA_TcnA treatment ranging from −2.92 to 2.97. Out of these 161 phosphosites were identified as significantly down-regulated while 69 phosphosites were significantly up-regulated (Figure 5A). The range and amount of regulated phosphosites both reflect a higher degree of regulation than on protein level. Strongest up-regulated phosphosites were, e.g., Thr-124 of Serine/Threonine Protein Kinase N2 (Pkn2) and the Thr-594 of Mitogen-activated protein kinase 7 (Mapk7). Phosphosites Ser-313 of Plakophilin 3 (Pkp3) and Thr-343 of NK-tumor recognition protein (Nktr) showed pronounced down-regulation.
For phosphorylation analysis via the IPATM software all 1832 significantly different (p < 0.05) phosphosites were uploaded and 889 IPA entries could be matched. Twelve of the pathways with the highest regulation status are shown in Figure 5B. What has to be mentioned in the context of phosphorylation analysis with IPA is that every phosphosite is attributed as ‘activating the pathway’ so in this case the resulting information would not be ‘activated or inhibited pathway’ based on the z-score, but rather pathways which exhibit net-positive or net-negative phosphosite regulations. In comparison to the proteome data (Figure 3C) more pathways showed decreased phosphorylation events which goes along with more phosphosites being observed as downregulated in TcnA treated cells. Actin Cytoskeleton and RHOA Signaling showed upregulation on proteome level as well as on phosphorylation levels. Signaling by Rho Family GTPases, Regulation of Actin-based Motility by Rho, Integrin and Insulin Receptor Signaling showed decreased phosphorylation, while being up-regulated on protein level. Rac-, Paxillin and Epithelial Adherens Junction Signaling, which did not belong to the strongest regulated pathways on proteome level, also showed pronounced down-regulation of phosphosites. Role of BRCA1 in DNA Damage Response was identified as the pathway with the strongest, significant up-regulation of phosphosites (while being inhibited on proteome level). As a subnetwork of Role of BRCA1 in DNA Damage Response, ATM Signaling also exhibited pronounced up-regulation of attributed phosphosites, however with a lower z-score.
A noteworthy protein analyzed in phosphoprotein analysis was Serine/arginine repetitive matrix protein 2 (Srrm2), a nearly 300 kDa protein which is part of the spliceosome and exhibited 71 significant phosphosites.
By application of known kinase–substrate interactions to the collected phosphorylation data, predictions about the regulation of kinase activities are possible. In order to do this a Fisher’s exact test (p-value < 0.05) was performed based on the presence of kinase–substrate motif annotations. The log2 mean regulation of substrates (p < 0.05) were calculated for each kinase in order to predict their activation status. This bioinformatic analysis highlighted ZIP kinase as well as Calmodulin-dependent protein kinase IV and II (CaMK VI; CaMK II) as activated kinases with mean log2 ratios of 0.28; 0.25 and 0.07, respectively (Figure 5). Especially cyclin-dependent kinases showed predicted inactivation with CDK5 and Cdc2 (also known as CDK1) showing slightly more pronounced downregulation of substrates in comparison with the overall CDK kinase motif (log2 ratios −0.39; −0.27, −0.27). The protein kinases with strongest inactivation however are Aurora-A kinase as well as AMP-activated protein kinase (log2 ratio −0.42; −0.41). Overall kinase prediction identified most kinases as inactivated, with only ZIP kinase, Calmodulin-dependent protein kinase IV and II and Chk1 kinase eliciting increased phosphorylated substrates after TcnA treatment.

3. Discussion

For the first time a global proteome and phosphoproteome analysis of TcnA treated target cells was performed. HEp-2 cells were used as target cells since they have proven before to be sensitive to various LCGTs and express required key proteins for cellular uptake of TcnA [17,24,25]. In this way, it was possible to identify cellular pathways regulated through toxin administration and deepen the understanding of particular effects of TcnA and general effects of LCGTs on target cells.
First, we tried to generate a suitable control for TcnA-induced cellular effects. An inactive toxin variant was not available as it is the case for Clostridioides difficile toxins TcdA and TcdB, for which genetically engineered mutants had been generated with impaired glucosyltransferase activity [26,27]. Thus, we inactivated the purified TcnA by treatment with formaldehyde and named it FA_TcnA. The quality of this inactivated version was analyzed by comparing it to non-treated control cells on three different levels: (1) general morphology of treated cells observable via light-microscopy; (2) alterations on proteome level; and (3) by alteration on phosphoproteome level. Cell morphology was similar for cells treated with FA_TcnA and non-treated control cells. On proteome level only the zinc-finger protein Znf611 was affected by FA_TcnA treated cells and from 9500 quantified phosphosites only 12 were regulated between FA_TcnA treatment and control. This comparison has been performed in order to exclude effects of high or low molecular weight compounds that might have been co-purified during toxin purification. Overall, we can say that all three tests proved FA_TcnA to be a very well suited control to characterize TcnA effects on HEp-2 cells.
Following, the focus was set on observed proteomic alterations induced by TcnA. A total of 5086 proteins were identified and quantified in every replicate (Supplementary Tables S1 and S2). Of those proteins, 14 expressed significant up-regulation and 21 significant down-regulation in response to TcnA (Figure 3). Two of those strongly regulated proteins (Tnni2 and Lgmn) and the interpretation of the data on single protein level will be discussed hereinafter.
The protein Legumain (Lgmn) showed up-regulation of factor 2.4 with strongest significance (−log2 p-value: 4.5). It is also referred to as asparagine endopeptidase (AEP) and primarily hydrolyzes asparaginyl bonds as well as aspartyl bonds at low pH, thus contributing to protein degradation [28,29,30]. It is mainly located within lysosomes but has also been detected extracellularly, in the nucleus or in the cytosol [31,32,33,34,35]. Each compartment with different pH influences its enzymatic functions [36]. In lysosomes, it takes part in antigen processing, while a carboxypeptidase function has been described at neutral pH and ligase activities have also been observed [37,38,39]. Since whole cell lysates were generated in our proteome analysis, it is not possible to tell whether Lgmn was enriched in a certain compartment, making hypothesis according to its function in cellular responses to TcnA tenuous. Probably Lgmn might be induced due to an increased protein turnover rate. However, to further specify the role of Lgmn in TcnA treated cells future experiments need to be carried out.
The protein Troponin I, fast skeletal muscle (Tnni2) is one of the strongest down-regulated proteins among the TcnA responsive proteins with a regulation factor of −4.7 and a p-value of 0.03. It is part of the inhibitory unit of the troponin complex which is mostly associated with striated muscles (skeletal or cardiac) and regulates muscle contraction [40]. The troponin complex is a combination of the Ca2+ binding protein TnC that connects the inhibitory subunit (TnI) and the tropomyosin binding subunit (TnT) [41]. By binding of Ca2+ a conformational change is induced which leads to exposition of myosin binding sides on the actin strands [42]. However, while the presence of Tnni2 in various epithelial cells has been reported, its exact function remains poorly explored. Sawaki et al. observed Tnni2 to be enriched in gastric cancer prone to metastasis and proposed a potential biomarker function [43]. Recently, Fu et al. reported a similar impact of Tnni2 in pancreatic cancer [44]. Along with the orphan nuclear receptor ERRα, Tnni2 was identified as a crucial protein contributing to increased cell proliferation, migration and invasion of pancreatic cancer by up-regulation of Sirt1 and downstream Syt8 [44]. Within the acquired data ERRα or Syt8 could not be detected. However, Sirt1 was detected in TcnA treated HEp-2 cells with a log2 difference of −0.1, and a p-value of 0.07. TcnA treatment thus seems to cause a reduction of proliferation and migration which fits very well with the morphological changes and is also consistent with the effects of the other LCGTs. Modulating effects of TcnA on cellular proliferation would also be supported by the predicted kinase regulations (see below). As described by the examples of Lgmn and Tnni2 strong regulation of certain single proteins are not always easily explainable. However, regulation of pathways identified with IPATM seems to correlate well to known effects of TcnA.
Since TcnA inactivates Rho GTPases in a similar process to TcdA and TcdB of C. difficile, comparable downstream effects should be expectable [12]. Due to previous investigations from our laboratory on the cellular response after TcdA or TcdB treatment [21,22], we were able to compare the effects of TcnA with those of other LCGTs. On proteome level, activation of Signaling by Rho Family GTPases, Actin Cytoskeleton Signaling, Epithelial Adherens Junction and Regulation of Actin-based Motility by Rho were expectedly observed. In contrast, RHOGDI Signaling was inactivated (Supplementary Table S3). An overlap of 23 regulated pathways with both C. difficile toxins TcdA and TcdB was identified (Table 1). However, predicted pathway regulation between TcdB and TcnA was closer to each other than to TcdA regulations. Overall TcdA was more prone to inactivation of pathways, while many activated pathways were observed for TcnA and TcdB treatment (Figure 4). When having a deeper look into the pathway analysis (Supplementary Table S3) various proteins are assigned in different pathways, such as RhoA which obviously takes part in RhoA signaling, but also in ILK and Integrin Signaling, Regulation of Actin-based Motility by Rho, RhoGDI Signaling and of course Signaling by Rho Family GTPases. In order to draw reliable conclusions, not only the activation status of pathways but also the individual regulation of participating proteins is important, since a single pathway can also contribute to various cellular outcomes (data not shown). In order to identify these outcomes in more detail, functional analysis needs to be performed. Looking at the top-regulated single proteins there was a stronger correlation among C. difficile toxins instead of TcnA compared to TcdA/TcdB. TcdA and TcdB share RhoB and Heme oxygenase-1 (Hmox1) as some of the strongest up-regulated proteins. RhoB is a direct target of inactivation by C. difficile toxins and up-regulation of RhoB followed by toxin administration is a commonly observed effect of TcdA and TcdB [21,22,45,46]. Hmox1 however has not yet been discussed in the course of C. difficile toxin administration. It is a heat shock protein and catalyzes the conversion of heme to biliverdin, ferrous ion, and carbon monoxide. Activation is mainly induced as a response to stress especially when cells are challenged with oxidative stress which would go along with the observed increase of reactive oxygen species (ROS) production especially after TcdB administration [47,48,49,50]. Additionally, cytoprotective effects by preventing TNF-α induced apoptosis are mediated by Hmox1 [51]. In TcnA treated cells, only Hmox2 could be detected, but exhibited insignificant regulation. This goes along with the statement that Hmox1 and 2 catalyze the same reaction, however Hmox1 is inducible, while Hmox2 is constitutively expressed [52]. Hmox1 has been identified as a promising target for inflammatory bowel disease treatment due to its modulating effects on the immune response. This causes reduction of inflammation in the gut in a variety of animal models and promotes bacterial clearance by interfering in the bacterial electron transport chain [53,54,55]. Since no observation of Hmox1 regulation could be observed for TcnA treated cells it is tempting to postulate that TcnA treatment is not accompanied by excessively increased ROS production. This hypothesis would also be supported by superoxide dismutase [Cu-Zn] (SOD1), catalase (CAT) and thioredoxin (TXN) levels observed in the proteome. All of those proteins participate in cellular oxidative stress defense [56,57,58]. However, only SOD1 was observed as significantly regulated with a log2 difference of only 0.18 while the other two proteins were only slightly down-regulated and did not show significance.
Another eye-catching protein family being down regulated by treatment with LCGTs are tropomyosins (Tpm1-3, Table 2). Apart from the same isoform of Tpm2 being strongly down regulated in both C. difficile toxin treatments, TcdB treatment leads to down regulation of Tpm2 as well as Tpm3. In TcdA treatment, Tpm1 was identified as one of the top-3 down-regulated proteins. All three tropomyosins share binding to actin filament organization as a biological process [59]. This activity is present in muscle cells but is also observed in non-muscle cells. With Tpm’s binding to actin-filaments the cellular actin-cytoskeleton is stabilized. Vice versa down regulation of Tpms could be regarded as a proteomic marker for actin-cytoskeleton reorganization. In TcnA treatment seven isoforms of Tpm1, Tpm2 and Tpm4 were measured. Of those proteins six were significantly altered between TcnA treated cells and control cells, but regulation factors were below the threshold of factor 2 and only one isoform of Tpm2 showed a minimal down regulation. However, for Tpm1, interaction with the troponin complex has been reported and Tnni2 as mentioned above is one of the strongest down-regulated proteins in TcnA treatment. Putatively, tropomyosin regulation might occur in shorter incubation time and thus, the chosen 24 h proteomic snapshot might be too long to detect stronger alterations.
In contrast to the proteome, more pronounced effects could be observed on phosphoproteome level. In every replicate 9427 phosphosites could be identified and quantified underlining the high reproducibility of the analysis. The regulatory effects of TcnA became particularly evident in the 230 responsive phosphosites which were characterized by a significant p-value (p < 0.05) and strongly altered abundance (regulation factor > 2; <−2). Some of the strongest regulated phosphosites will be discussed hereinafter.
Serine/threonine-protein kinase N2 (Pkn2) has been identified as a protein worth discussing due to its regulated phosphosite Thr-124. TcnA treatment led to significantly (−log10 p-value: 3) increased phosphorylation at this site with a log2 regulation factor of 2.7 (Figure 5). Pkn2 is a Rho/Rac effector protein which takes part in various cellular processes such as signal transduction, cell cycle progression, actin cytoskeleton assembly, cell adhesion and others [60,61,62,63]. Pkn2 is known for binding RhoA, subsequently promoting a variety of RhoA-regulated processes. The regulated phosphosite at position 124 is located within the binding domain of RhoA [60]. Phosphorylation at this site leads to increased binding of Pkn2 to RhoA, thus promoting RhoA signaling. Degradation of RhoA via Smurf1 is hampered when Pkn2 is bound [60], which would go along with the observed slight increase in RhoA protein abundance (log2 difference of 0.36), and the observed activation of RhoA signaling pathway. What remains unclear was the effect of the N-acetylglucosamination of RhoA, which would result in an inactivation of RhoA downstream effects. Additional analysis would be helpful to identify whether up-regulation of RhoA could be another side effect of TcnA treatment similar to the inactivation and up-regulation of RhoB in TcdA treatment [45,46].
Another Protein, which possessed a phosphosite sensible to TcnA treatment, was Plakophilin 3 (Pkp3). The phosphosite Ser-313 of Pkp3 has been identified as significantly (−log10 p-value: 5) down regulated with a log2 regulation factor of −2.9 (Figure 5). Pkp3 is localized in desmosomes, structures at the cell membrane that are extracellularly connected to cadherins of surrounding cells and intracellularly bound to intermediate filaments [64,65]. Ser-313 phosphosite has been identified as responsible for a decreased membrane localization of Pkp3 leading to reduced tricellular contacts [66]. The observed decrease in phosphorylation after TcnA treatment should lead to increased cell–cell adhesion, which contradicts the observed phenotype of cell rounding accompanied by loss of cell–cell contacts. However, in former analysis a loss of Pkp3 did not exhibit pronounced adhesion defects, because plakophilin 1 (Pkp1), another protein of the plakophilin family, can rescue Pkp3 knockout concerning adherens junctions [67]. Pkp1 could not be detected in any of our treatment conditions. Aside from its functions concerning organization of the cytoskeleton, additional effects within cellular stress response or control of protein biosynthesis are known, although not associated with the regulated phosphosite [68,69].
The focal adhesion adapter protein Paxillin is another protein that was identified to have TcnA responsive phosphosites. In TcnA measurements four phosphosites of Paxillin showed downregulation, the strongest at Serine-258 with a log2 difference of −2.2 (Supplementary Table S2). Additionally, downregulation of the Paxillin pathway based on phosphorylation data could be identified (Figure 5). Paxillin serves as a scaffold protein in reactions mainly associated with cell movement and migration depending on its phosphorylation status on specific Tyrosine and Serine residues [70]. Anti-apoptotic effects of Paxillin, depending on its phosphorylation status on various phosphosites, have been reported as well [71]. Interestingly, Paxillin dephosphorylation has been reported to be sensitive to treatment with other LCGTs, namely TcdA and TcsL, as well [14,72]. In TcdA experiments, Paxillin dephosphorylation could be identified at the Tyrosine in position 118 [72,73]. This phosphosite is primarily regulated by the tyrosine-protein kinase Src [74]. TcdA is able to directly bind to the catalytic Src domain, thereby reducing the phosphorylation of Tyrosine 118 of Paxillin independently of Rho glucosylation [72]. After mutating the Tyrosine at position 118 to Phenylalanine the anti-apoptotic effects of Paxillin were suppressed [71]. TcsL dephosphorylation of Paxillin does not result from inactivation of Src but has been contributed to earlier cellular response events [14]. The phosphosite Ser-258 which was identified as down-regulated within our experiment has been identified to be regulated by Pak1 [75]. Against the background of HIV infection, its phosphorylation is associated with inhibition of the TNFα converting enzyme ADAM 17 binding to Paxillin, leading to an inhibition of lipid raft transfer [75]. It remains to be elucidated whether TcsL and TcnA treated cells address dephosphorylation of Paxillin in a similar manner or if there is yet another possibility for LCGTs to target Paxillin signaling.
Finally, based on phosphosite regulation predictions about the activity of associated kinases are possible. In our analysis activation for ZIP kinase, CaMK II & IV and Chk1 kinase was predicted. Shared downstream effects result in alterations concerning cell cycle progression (CaMK IV & II, Chk1), apoptosis (Zip, CaMK II, Chk1), as well as cell morphology (CaMK II &IV) [76,77,78,79,80]. On the other hand, Amp-activated protein Kinase has been predicted as down-regulated based on substrate regulation. It is a key protein in keeping cellular ATP homeostasis and, upon its activation, especially energy-consuming processes such as CDK-mediated cell-cycle progression are suppressed [81,82,83]. However, there have to be additional regulatory effects surpassing the effect of Amp-activated protein kinase for CDKs after TcnA treatment since substrates expressing the general CDK kinase substrate motif tended to be down-regulated as well as CDK1 and CDK5 substrates, respectively. A possible explanation would be the regulation of CDKs via cyclins, regulatory proteins such as p38 or inhibiting enzymes such as Wee1 or Myt1 [84,85]. Aurora-A kinase was the kinase with strongest inactivation predicted and has already been observed as inhibited after treatment of Hela cells with Toxin B from C. difficile serotype F strain 1470 (TcdBF) [86]. This takes part in cell proliferation and is usually activated during G2 phase to M phase transition. Degradation of Aurora A kinase has been related to enhanced apoptosis while its up-regulation is frequently observed in breast cancer [87]. All in all, kinase prediction leads to the conclusion that TcnA treatment decreases cellular proliferation in a complex pattern.
In this study we started with establishing a well-suited control concerning the analysis of TcnA effects on HEp-2 cells by using inactivated FA_TcnA. We were able to observe distinct regulations of TcnA on the proteome and phosphoproteome of HEp-2 cells, which could be illustrated not only via direct regulation of single proteins or phosphosites but also by observing differentially affected pathways on phospho- or protein level. Various pathways seemed to be affected similarly to TcdB-treated HEp-2 cells such as Integrin Signaling and Paxillin Signaling while other pathways such as Actin Cytoskeleton Signaling or Signaling by Rho Family GTPases showed activation across treatments with TcnA, TcdA and TcdB. A similar impact on the cellular proteome in comparison to TcdB would also fit to the close substrate GTPase profiles of TcdB and TcnA [73]. With the phosphosite regulation of Paxillin, a common target protein not only of TcnA and TcdA but also TcsL was identified. It might be promising to have a closer look, especially concerning the time scale when these phosphorylation events take place in order to match the identified regulations with cellular functions. Interestingly, compared with TcdA and B, there has been no evidence of increased ROS production by TcnA under the chosen conditions. All in all, a broad first insight into (phospho-) proteomic regulations of TcnA treated HEp-2 cells could be illustrated. However, remodeling time scaling and experimental set up should enable conclusions that are even more precise and, of course, functional analysis of regulated phosphosites/proteins would help break down the detailed effects caused in target cells. Apart from this compartmental enrichment for e.g., the case of legumain would be helpful in order to identify the functional effects underlying its regulation status.

4. Materials and Methods

4.1. Purification and Inactivation of Clostridium Novyi α-Toxin

Purification of α-Toxin (TcnA) from Clostridium novyi strain ATCC 19,402 was performed as previously described [88]. Briefly, the bacteria were cultivated (BHI, Difco, Sparks, MD, USA) in a dialysis bag for 72 h at 37 °C under microaerophilic conditions. Culture supernatant was ammonium sulfate precipitated and the protein pellet was extracted with buffer A (20 mM Tris, 15 mM NaCl pH 8.0) and dialyzed against the same buffer. This fraction was separated by anion exchange chromatoghraphy (HiPrepTM Q Fast Flow 16/10 column, Sigma-Aldrich, St. Louis, MO, USA). Proteins were eluted by a linear gradient from 15 mM NaCl to 1 M NaCl in buffer A. Fractions were tested for their toxicity on MEF-cells. Fractions that showed rounding of cells after 4 h of incubation were combined and concentrated with Vivaspin 2 centrifugal concentrators (100 MWCO, Sigma-Aldrich). Seventy-five percent Glycerin was added in a ratio of 1:2 (v/v) for stabilization purposes. Aliquots were stored at −80 °C. Protein content of the TcnA concentrate was calculated using DC protein assay (Bio-Rad, Hercules, CA, USA). For control experiments, part of the purified TcnA was inactivated by incubation with 3.7% formaldehyde (Sigma-Aldrich, MO, USA) at 37 °C for 3 h. Inactivated toxin was referred to as FA_TcnA. Activity of TcnA was evaluated by cell rounding of HEp-2 cells using phase-contrast microscopy.

4.2. Cell Culture

Cell Culture was proceeded as previously described [89]. HEp-2 cells were maintained in a 75 cm2 flask in humidified atmosphere at 37 °C and 5% CO2. Cells were cultured in minimal essential medium (GibcoTM, Paisley, UK) supplemented with 10% fetal bovine serum, 100 U/mL penicillin and streptomycin each (Sigma-Aldrich, St. Louis, MO, USA). Cells were split, depending on confluency, in order to maintain vitality.

4.3. Toxin Treatment of HEp-2 Cells and Sample Preparation

At a confluency of 75% medium was replaced with 20 mL of fresh medium. For treatment 255 ng/mL of TcnA or FA_TcnA were added to 20 mL medium. Control cells were grown in medium without addition of toxin. Morphological changes of cells were documented by phase-contrast microscopy after 24 h. These conditions were chosen in order to reach maximal percentage of rounded cells without any undesired side effects as necrosis. Cells were washed twice with ice-cold PBS and stored at −80 °C overnight. The cells were then lysed by scraping them in 300 μL lysis buffer (8 M Urea, 50 mM ammonium bicarbonate (pH 8.0), 1 mM sodium ortho-vanadate, complete EDTA-free protease inhibitor cocktail (Roche, Mannheim, Germany), and phosphoSTOP phosphatase inhibitor cocktail (Roche, Mannheim, Germany)). Lysates were sonicated on ice for 15 s at 30% energy output. Cell debris was removed by centrifugation for 15 min at 13,000× g at 4 °C. The supernatant was used directly or stored at −80 °C.

4.4. Protein Digestion

Tryptic digestion of proteins proceeded as previously described [89]. Protein concentrations were calculated using DC protein assay (Bio-Rad, Hercules, CA, USA). Proteins were reduced with 5 mM DTT at 37 °C for 1 h and afterwards alkylated with 10 mM iodoacetamide at RT for 30 min. Alkylation was quenched by adding DTT to a final concentration of 5 mM. For digestion, lysates were diluted 1:5 with 50 mM ammonium bicarbonate (ABC) buffer. One mg protein was digested per condition at 37 °C for 4 h using 10 μg Lys-C (Wako, Osaka, Japan) followed by overnight digestion with 10 μg trypsin (Promega, Madison, WI, USA). Digestion was stopped by adding TFA to a final concentration of 1%. Peptide solutions were desalted with Sep Pak C18 1cc cartridges (Waters Corporation, Milford, MA, USA).

4.5. Tandem Mass Tags Labeling and Phosphopeptide Enrichment

From each sample 1 mg peptides were labeled with TMT10plex™ Isobaric Label Reagent (Thermo ScientificTM, Rockford, IL, USA) according to the manufacturer’s protocol. Labeled samples were pooled and 1 mg equivalent was used for peptide measurements (proteome data) while 5 mg equivalent was used for phosphopeptide enrichment. Both samples were again desalted with Sep Pak C18 1cc cartridges (Waters Corporation, Milford, MA, USA).
Phosphopeptides were enriched using High-Select™ TiO2 Phosphopeptide Enrichment Kit (Thermo ScientificTM, Rockford, IL, USA) according to the manufacturer’s protocol. All flow-throughs after application of the sample were collected and dried by vacuum centrifugation. The residues were then applied to High-Select™ Fe-NTA Phosphopeptide Enrichment Kit (Thermo ScientificTM, Rockford, IL, USA) enabling reproducible enrichment of phosphopeptides with a higher degree of phosphorylation. Eluates of both enrichment procedures (TiO2; Fe-NTA) were pooled and eight fractions generated by high pH reversed-phase peptide fractionation (Thermo ScientificTM, Rockford, IL, USA).

4.6. LC-MS Analysis

Dried peptides were dissolved in 0.1% TFA/2% ACN and analyzed in an Orbitrap Fusion Lumos mass spectrometer (Thermo Fisher Scientific) equipped with a nano-electrospray source and connected to an Ultimate 3000 RSLC nanoflow system (Thermo Fisher Scientific). Peptides were loaded onto a C18 PepMap100 (5 µm, 100 Å) µ-Precolumn (Thermo Fisher Scientific) and separated by a 50 cm µPACTM analytical column (Pharma fluidics) at 35 °C column temperature. A binary gradient with solvent A consisting of 0.1% formic acid (Merck, Darmstadt, Germany) and solvent B consisting of 100% ACN/0.1% formic acid was used to separate eluting peptides. Columns were equilibrated with 3.4% solvent B for 5 min. Separation of peptides followed by increasing B to 21% within 60 min and up to 42% in the following 32 min. Within 2 min B was increased to 75.6% which was held for 3 min and then decreased back to 3.4%. Re-equilibration of the column was reached by keeping B at 3.4% for 16 min. The spray voltage was set to 2 kV. Measurements were done with a DDA approach, with a cycle time of 3 s and Top N setting. Dynamic exclusion was set to 60 s, AGC target at 4 × 105 and a maximum injection time of 50 ms was applied. Orbitrap resolution for the MS1 scan was set to 120,000 and higher-energy collisional dissociation (HCD) fragmentation at 38% was used followed by MS2 measurements. The first mass was set to 100 m/z in order to be able to quantify the split reporter ions ranging from 126–131 m/z. MS2 maximum injection time was set to 110 ms and corresponding isolation width to 0.8 m/z. Orbitrap resolution for MS2 scans was set to 60,000.

4.7. Data Processing

Raw data were processed with MaxQuant software (version 1.6.14.0, Martinsried, Germany) [90] using the Andromeda search engine [91]. Spectra were searched against the Swiss-Prot reviewed UniprotKB Homo sapiens database (version 01/2020, 42,352 entries) [92]. Carbamidomethylation of cysteine was set as fixed modification while oxidation of methionine, N-terminal acetylation, deamidation of glutamine and asparagine were set as variable modifications. N-acetylglucosaminylation (C8H13NO5; 203.079 Da) was considered as a variable modification on Threonine residues due to TcnA treatment. For Phosphopeptide-enriched samples the variable modification phosphorylation at serine, threonine, and tyrosine residues was considered as well. Reporter ion MS2 for TMT 10-plex was set with the associated isotopic distributions. Second peptide search was activated and match between runs (0.7 min match time window; 20 min alignment time window) enabled. False discovery rate was set to 0.01 and a maximum of two missed cleavages was allowed.
Only proteins that were quantified in every replicate were used for further data processing. Phosphosites additionally needed to have a localization probability of at least 75%. Both proteome and phosphoproteome were normalized by subtracting the mean intensity of each sample. Phosphosite analysis was solely performed if the corresponding protein could be quantified within the proteome measurements. Phosphosite abundancies were normalized to the corresponding protein abundancy. Thus, identified differences trace back to the degree of phosphorylation instead of protein expression.
Data evaluation and analysis were done using Perseus (1.6.14.0) [93] and Microsoft Excel (Microsoft Corporation, 2018. Microsoft Excel, Available online: https://office.microsoft.com/excel, 26 August 2022). Data visualization was done with GraphPad Prism (GraphPad Prism 7.00 version for Windows, GraphPad Software, San Diego, CA, USA, www.graphpad.com). Upstream and canonical pathway analyses were performed using ingenuity pathway analysis (IPA; QIAGEN Inc., Hilden, Germany, https://digitalinsights.qiagen.com/IPA, 26 August 2022) [94]. This software predicts regulatory proteins by comparing up- and down-regulated proteins (or phosphosites) to a manually curated database. Even if a protein could not be identified itself, conclusions can be made due to known relations, e.g., based on phosphorylation data predictions about the activity of certain kinases can be made. For IPA analysis all significantly different phosphosites and proteins were used (p-value > 0.05).

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms23179939/s1.

Author Contributions

Conceptualization, T.S. and A.P.; methodology, T.S., A.P. and H.G.; formal analysis, T.S., A.P. and H.G.; investigation, T.S.; resources: A.P.; data curation, T.S., A.P. and H.G.; writing—original draft preparation, T.S.; writing—review and editing, T.S., A.P. and H.G.; visualization, T.S.; supervision, A.P.; project administration, T.S. and A.P.; funding acquisition, A.P. All authors have read and agreed to the published version of the manuscript.

Funding

This work was partially funded by the Federal State of Lower Saxony, Niedersächsisches Vorab (VWZN3380, Project B1, B2).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium via the PRIDE [95] partner repository with the dataset identifier PXD036078 and 10.6019/PXD036078.

Acknowledgments

We are grateful to Ralf Gerhard for providing the C. novyi strain, Ilona Schelle for preparing TcnA as well as Ulrike Schrameck and Karsten Heidrich for excellent technical assistance.

Conflicts of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

  1. Buboltz, J.B.; Murphy-Lavoie, H.M. Gas gangrene. In StatPearls; StatPearls Publishing: Tampa, FL, USA, 2022. [Google Scholar]
  2. Navarro, M.A.; Uzal, F.A. Pathobiology and diagnosis of clostridial hepatitis in animals. J. Vet. Diagn. Investig. 2019, 32, 192–202. [Google Scholar] [CrossRef]
  3. Oakley, C.L.; Warrack, G.H.; Clarke, P.H. The Toxins of Clostridium oedematiens (Cl. novyi). J. Gen. Microbiol. 1947, 1, 91–107. [Google Scholar] [CrossRef] [PubMed]
  4. Eklund, M.W.; Poysky, F.T.; Peterson, M.E.; Meyers, J.A. Relationship of bacteriophages to alpha toxin production in Clostridium novyi types A and B. Infect. Immun. 1976, 14, 793–803. [Google Scholar] [CrossRef]
  5. Skarin, H.; Segerman, B. Plasmidome Interchange between Clostridium botulinum, Clostridium novyi and Clostridium haemolyticum Converts Strains of Independent Lineages into Distinctly Different Pathogens. PLoS ONE 2014, 9, e107777. [Google Scholar] [CrossRef] [PubMed]
  6. Nakamura, S.; Kimura, I.; Yamakawa, K.; Nishida, S. Taxonomic Relationships among Clostridium novyi Types A and B, Clostridium haemolyticum and Clostridium botulinum Type C. Microbiology 1983, 129, 1473–1479. [Google Scholar] [CrossRef] [PubMed]
  7. Hunter, L.C.; Poxton, I.R. Clostridium botulinum types C and D and the closely related Clostridium novyi. Rev. Med. Microbiol. 2002, 13, 75–90. [Google Scholar] [CrossRef]
  8. Bryant, A.E.; Stevens, D.L. Clostridial toxins in the pathogenesis of gas gangrene. In The Comprehensive Sourcebook of Bacterial Protein Toxins; Elsevier Ltd.: Amsterdam, The Netherlands, 2015. [Google Scholar] [CrossRef]
  9. Orrell, K.E.; Melnyk, R.A. Large Clostridial Toxins: Mechanisms and Roles in Disease. Microbiol. Mol. Biol. Rev. 2021, 85, e0006421. [Google Scholar] [CrossRef] [PubMed]
  10. Richard, J.F.; Petit, L.; Gibert, M.; Marvaud, J.C.; Bouchaud, C.; Popoff, M.R. Bacterial toxins modifying the actin cytoskeleton. Int. Microbiol. 1999, 2, 185–194. [Google Scholar] [CrossRef]
  11. Jank, T.; Aktories, K. Structure and mode of action of clostridial glucosylating toxins: The ABCD model. Trends Microbiol. 2008, 16, 222–229. [Google Scholar] [CrossRef]
  12. Selzer, J.; Hofmann, F.; Rex, G.; Wilm, M.; Mann, M.; Just, I.; Aktories, K. Clostridium novyi α-Toxin-catalyzed Incorporation of GlcNAc into Rho Subfamily Proteins. J. Biol. Chem. 1996, 271, 25173–25177. [Google Scholar] [CrossRef]
  13. Busch, C.; Schömig, K.; Hofmann, F.; Aktories, K. Characterization of the Catalytic Domain of Clostridium novyi Alpha-Toxin. Infect. Immun. 2000, 68, 6378–6383. [Google Scholar] [CrossRef] [PubMed]
  14. Popoff, M.R.; Geny, B. Rho/Ras-GTPase-dependent and -independent activity of clostridial glucosylating toxins. J. Med. Microbiol. 2011, 60, 1057–1069. [Google Scholar] [CrossRef] [PubMed]
  15. Belyi, Y.; Aktories, K. Bacterial toxin and effector glycosyltransferases. Biochim. Biophys. Acta (BBA) Gen. Subj. 2010, 1800, 134–143. [Google Scholar] [CrossRef] [PubMed]
  16. Oksche, A.; Nakov, R.; Habermann, E. Morphological and biochemical study of cytoskeletal changes in cultured cells after extracellular application of Clostridium novyi alpha-toxin. Infect. Immun. 1992, 60, 3002–3006. [Google Scholar] [CrossRef]
  17. Zhou, Y.; Li, D.; Luo, J.; Chen, A.; Li, X.; Pan, Z.; Wan, L.; He, L.; Li, D.; Li, Y.; et al. Sulfated glycosaminoglycans and low-density lipoprotein receptor mediate the cellular entry of Clostridium novyi alpha-toxin. Cell Res. 2021, 31, 935–938. [Google Scholar] [CrossRef]
  18. Bateman, N.W.; Goulding, S.P.; Shulman, N.J.; Gadok, A.K.; Szumlinski, K.K.; MacCoss, M.J.; Wu, C.C. Maximizing Peptide Identification Events in Proteomic Workflows Using Data-Dependent Acquisition (DDA). Mol. Cell. Proteom. 2014, 13, 329–338. [Google Scholar] [CrossRef]
  19. Engholm-Keller, K.; Larsen, M.R. Technologies and challenges in large-scale phosphoproteomics. Proteomics 2013, 13, 910–931. [Google Scholar] [CrossRef]
  20. Nilsson, C.L. Advances in Quantitative Phosphoproteomics. Anal. Chem. 2011, 84, 735–746. [Google Scholar] [CrossRef]
  21. Junemann, J.; Birgin, G.; Erdmann, J.; Schröder, A.; Just, I.; Gerhard, R.; Pich, A. Toxin A of the nosocomial pathogen Clostridium difficile induces primary effects in the proteome of HEp-2 cells. Proteom. Clin. Appl. 2016, 11, 1600031. [Google Scholar] [CrossRef]
  22. Erdmann, J.; Junemann, J.; Schröder, A.; Just, I.; Gerhard, R.; Pich, A. Glucosyltransferase-dependent and -independent effects of TcdB on the proteome of HEp-2 cells. Proteomics 2017, 17, 1600435. [Google Scholar] [CrossRef]
  23. Hochberg, Y.; Benjamini, Y. More powerful procedures for multiple significance testing. Stat. Med. 1990, 9, 811–818. [Google Scholar] [CrossRef] [PubMed]
  24. IPATM of QIAGEN’s Ingenuity Pathway Analysis. Calculating and Interpreting the p-Values for Functions, Pathways and Lists in IPA. 2016. Available online: https://www.ingenuity.com/wp-content/themes/ingenuity-qiagen/pdf/ipa/functions-pathways-pval-whitepaper.pdf (accessed on 28 February 2022).
  25. Ingenuity® Systems. Ingenuity Downstream Effects Analysis in IPA®. 2012. Available online: http://pages.ingenuity.com/rs/ingenuity/images/0812%20downstream_effects_analysis_whitepaper.pdf (accessed on 28 February 2022).
  26. Gerhard, R.; Nottrott, S.; Schoentaube, J.; Tatge, H.; Olling, A.; Just, I. Glucosylation of Rho GTPases by Clostridium difficile toxin A triggers apoptosis in intestinal epithelial cells. J. Med. Microbiol. 2008, 57, 765–770. [Google Scholar] [CrossRef] [PubMed]
  27. Wohlan, K.; Goy, S.; Olling, A.; Srivaratharajan, S.; Tatge, H.; Genth, H.; Gerhard, R. Pyknotic cell death induced by Clostridium difficile TcdB: Chromatin condensation and nuclear blister are induced independently of the glucosyltransferase activity. Cell. Microbiol. 2014, 16, 1678–1692. [Google Scholar] [CrossRef] [PubMed]
  28. Ishii, S. Asparaginylendopeptidase: An enzyme probably responsible to post-translational proteolysis and transpeptidation of proconcanavalin A. Seikagaku J. Jpn. Biochem. Soc. 1993, 65, 185–189. [Google Scholar]
  29. Dall, E.; Brandstetter, H. Structure and function of legumain in health and disease. Biochimie 2016, 122, 126–150. [Google Scholar] [CrossRef]
  30. Rotari, V.I.; Dando, P.M.; Barrett, A.J. Legumain Forms from Plants and Animals Differ in Their Specificity. Biol. Chem. 2001, 382, 953–959. [Google Scholar] [CrossRef] [PubMed]
  31. Chen, J.-M.; Dando, P.M.; Stevens, R.A.E.; Fortunato, M.; Barrett, A. Cloning and expression of mouse legumain, a lysosomal endopeptidase. Biochem. J. 1998, 335, 111–117. [Google Scholar] [CrossRef] [PubMed]
  32. Haugen, M.H.; Johansen, H.T.; Pettersen, S.J.; Solberg, R.; Brix, K.; Flatmark, K.; Maelandsmo, G.M. Nuclear Legumain Activity in Colorectal Cancer. PLoS ONE 2013, 8, e52980. [Google Scholar] [CrossRef]
  33. Lin, Y.; Qiu, Y.; Xu, C.; Liu, Q.; Peng, B.; Kaufmann, G.F.; Chen, X.; Lan, B.; Wei, C.; Lu, D.; et al. Functional Role of Asparaginyl Endopeptidase Ubiquitination by TRAF6 in Tumor Invasion and Metastasis. JNCI J. Natl. Cancer Inst. 2014, 106, dju012. [Google Scholar] [CrossRef]
  34. Smith, R.; Johansen, H.T.; Nilsen, H.; Haugen, M.H.; Pettersen, S.J.; Mælandsmo, G.M.; Abrahamson, M.; Solberg, R. Intra- and extracellular regulation of activity and processing of legumain by cystatin E/M. Biochimie 2012, 94, 2590–2599. [Google Scholar] [CrossRef] [Green Version]
  35. Solberg, R.; Smith, R.; Almlöf, M.; Tewolde, E.; Nilsen, H.; Johansen, H.T. Legumain expression, activity and secretion are increased during monocyte-to-macrophage differentiation and inhibited by atorvastatin. Biol. Chem. 2014, 396, 71–80. [Google Scholar] [CrossRef] [PubMed]
  36. Vidmar, R.; Vizovisek, M.; Turk, D.; Turk, B.; Fonović, M. Protease cleavage site fingerprinting by label-free in-gel degradomics reveals pH -dependent specificity switch of legumain. EMBO J. 2017, 36, 2455–2465. [Google Scholar] [CrossRef] [PubMed]
  37. Dall, E.; Brandstetter, H. Mechanistic and structural studies on legumain explain its zymogenicity, distinct activation pathways, and regulation. Proc. Natl. Acad. Sci. USA 2013, 110, 10940–10945. [Google Scholar] [CrossRef] [PubMed]
  38. Dall, E.; Stanojlovic, V.; Demir, F.; Briza, P.; Dahms, S.O.; Huesgen, P.F.; Cabrele, C.; Brandstetter, H. The Peptide Ligase Activity of Human Legumain Depends on Fold Stabilization and Balanced Substrate Affinities. ACS Catal. 2021, 11, 11885–11896. [Google Scholar] [CrossRef]
  39. Manoury, B.; Hewitt, E.; Morrice, N.; Dando, P.M.; Barrett, A.J.; Watts, C. An asparaginyl endopeptidase processes a microbial antigen for class II MHC presentation. Nature 1998, 396, 695–699. [Google Scholar] [CrossRef] [PubMed]
  40. Mullen, A.J.; Barton, P.J. Structural characterization of the human fast skeletal muscle troponin I gene (TNNI2). Gene 2000, 242, 313–320. [Google Scholar] [CrossRef]
  41. Farah, C.S.; Reinach, F.C. The troponin complex and regulation of muscle contraction. FASEB J. 1995, 9, 755–767. [Google Scholar] [CrossRef]
  42. Sawaki, K.; Kanda, M.; Miwa, T.; Umeda, S.; Tanaka, H.; Tanaka, C.; Kobayashi, D.; Suenaga, M.; Hattori, N.; Hayashi, M.; et al. Troponin I2 as a Specific Biomarker for Prediction of Peritoneal Metastasis in Gastric Cancer. Ann. Surg. Oncol. 2018, 25, 2083–2090. [Google Scholar] [CrossRef]
  43. Herzberg, O.; Moult, J.; James, M.N. A model for the Ca2+-induced conformational transition of troponin C. A trigger for muscle contraction. J. Biol. Chem. 1986, 261, 2638–2644. [Google Scholar] [CrossRef]
  44. Fu, Z.; Liang, X.; Shi, L.; Tang, L.; Chen, D.; Liu, A.; Shao, C. SYT8 promotes pancreatic cancer progression via the TNNI2/ERRα/SIRT1 signaling pathway. Cell Death Discov. 2021, 7, 390. [Google Scholar] [CrossRef]
  45. Gerhard, R.; Tatge, H.; Genth, H.; Thum, T.; Borlak, J.; Fritz, G.; Just, I. Clostridium difficile Toxin A Induces Expression of the Stress-induced Early Gene Product RhoB. J. Biol. Chem. 2005, 280, 1499–1505. [Google Scholar] [CrossRef] [PubMed]
  46. Genth, H.; Dreger, S.C.; Huelsenbeck, J.; Just, I. Clostridium difficile toxins: More than mere inhibitors of Rho proteins. Int. J. Biochem. Cell Biol. 2008, 40, 592–597. [Google Scholar] [CrossRef] [PubMed]
  47. Parsell, D.A.; Lindquist, S. The Function of Heat-Shock Proteins in Stress Tolerance: Degradation and Reactivation of Damaged Proteins. Annu. Rev. Genet. 1993, 27, 437–496. [Google Scholar] [CrossRef] [PubMed]
  48. Poss, K.D.; Tonegawa, S. Reduced stress defense in heme oxygenase 1-deficient cells. Proc. Natl. Acad. Sci. USA 1997, 94, 10925–10930. [Google Scholar] [CrossRef] [PubMed]
  49. Kamalvand, G.; Pinard, G.; Ali-Khan, Z. Heme-oxygenase-1 response, a marker of oxidative stress, in a mouse model of AA amyloidosis. Amyloid 2003, 10, 151–159. [Google Scholar] [CrossRef]
  50. Frädrich, C.; Beer, L.-A.; Gerhard, R. Reactive Oxygen Species as Additional Determinants for Cytotoxicity of Clostridium difficile Toxins A and B. Toxins 2016, 8, 25. [Google Scholar] [CrossRef]
  51. Physiol, A.J.; Cell, L.; Physiol, M.; Petrache, I.; Otterbein, L.E.; Alam, J.; Wiegand, G.W.; Augustine, M.K.; Petrache, I.; Otterbein LE, O.E.; et al. Heme oxygenase-1 inhibits TNF-α-induced apoptosis in cultured fibroblasts. Am. J. Physiol. Lung Cell. Mol. Physiol. 2000, 278, L312–L319. [Google Scholar]
  52. Cruse, I.; Maines, M.D. Evidence suggesting that the two forms of heme oxygenase are products of different genes. J. Biol. Chem. 1988, 263, 3348–3353. [Google Scholar] [CrossRef]
  53. Onyiah, J.C.; Sheikh, S.Z.; Maharshak, N.; Steinbach, E.C.; Russo, S.M.; Kobayashi, T.; Mackey, L.C.; Hansen, J.J.; Moeser, A.; Rawls, J.; et al. Carbon Monoxide and Heme Oxygenase-1 Prevent Intestinal Inflammation in Mice by Promoting Bacterial Clearance. Gastroenterology 2013, 144, 789–798. [Google Scholar] [CrossRef]
  54. Rana, N.; McLean, S.; Mann, B.E.; Poole, R.K. Interaction of the carbon monoxide-releasing molecule Ru(CO)3Cl(glycinate) (CORM-3) with Salmonella enterica serovar Typhimurium: In situ measurements of carbon monoxide binding by integrating cavity dual-beam spectrophotometry. Microbiology 2014, 160, 2771–2779. [Google Scholar] [CrossRef]
  55. Sebastián, V.P.; Salazar, G.; Coronado-Arrázola, I.; Schultz, B.M.; Vallejos, O.P.; Berkowitz, L.; Álvarez-Lobos, M.M.; Riedel, C.; Kalergis, A.; Bueno, S.M. Heme Oxygenase-1 as a Modulator of Intestinal Inflammation Development and Progression. Front. Immunol. 2018, 9, 1956. [Google Scholar] [CrossRef] [Green Version]
  56. Hwang, J.; Jin, J.; Jeon, S.; Moon, S.H.; Park, M.Y.; Yum, D.-Y.; Kim, J.H.; Kang, J.-E.; Park, M.H.; Kim, E.-J.; et al. SOD1 suppresses pro-inflammatory immune responses by protecting against oxidative stress in colitis. Redox Biol. 2020, 37, 101760. [Google Scholar] [CrossRef] [PubMed]
  57. Nishikawa, M.; Hashida, M.; Takakura, Y. Catalase delivery for inhibiting ROS-mediated tissue injury and tumor metastasis. Adv. Drug Deliv. Rev. 2009, 61, 319–326. [Google Scholar] [CrossRef] [PubMed]
  58. Nordberg, J.; Arnér, E.S. Reactive oxygen species, antioxidants, and the mammalian thioredoxin system. Free Radic. Biol. Med. 2001, 31, 1287–1312. [Google Scholar] [CrossRef]
  59. Janco, M.; Bonello, T.T.; Byun, A.; Coster, A.C.F.; Lebhar, H.; Dedova, I.; Gunning, P.W.; Böcking, T. The impact of tropomyosins on actin filament assembly is isoform specific. BioArchitecture 2016, 6, 61–75. [Google Scholar] [CrossRef] [PubMed]
  60. Guen, V.J.; Gamble, C.; Perez, D.E.; Bourassa, S.; Zappel, H.; Gärtner, J.; Lees, J.A.; Colas, P. STAR syndrome-associated CDK10/Cyclin M regulates actin network architecture and ciliogenesis. Cell Cycle 2016, 15, 678–688. [Google Scholar] [CrossRef]
  61. Kubouchi, K.; Mukai, H. PKN2 is involved in aggregation and spheroid formation of fibroblasts in suspension culture by regulating cell motility and N-cadherin expression. Biochem. Biophys. Rep. 2021, 25, 100895. [Google Scholar] [CrossRef]
  62. Ruby, M.A.; Riedl, I.; Massart, J.; Åhlin, M.; Zierath, J.R. Protein kinase N2 regulates AMP kinase signaling and insulin responsiveness of glucose metabolism in skeletal muscle. Am. J. Physiol. Metab. 2017, 313, E483–E491. [Google Scholar] [CrossRef]
  63. Schmidt, A.; Durgan, J.; Magalhães, A.; Hall, A. Rho GTPases regulate PRK2/PKN2 to control entry into mitosis and exit from cytokinesis. EMBO J. 2007, 26, 1624–1636. [Google Scholar] [CrossRef]
  64. Bonné, S.; Gilbert, B.; Hatzfeld, M.; Chen, X.; Green, K.J.; Van Roy, F. Defining desmosomal plakophilin-3 interactions. J. Cell Biol. 2003, 161, 403–416. [Google Scholar] [CrossRef]
  65. Gosavi, P.; Kundu, S.T.; Khapare, N.; Sehgal, L.; Karkhanis, M.S.; Dalal, S.N. E-cadherin and plakoglobin recruit plakophilin3 to the cell border to initiate desmosome assembly. Experientia 2010, 68, 1439–1454. [Google Scholar] [CrossRef] [PubMed]
  66. Neuholz, M.E. Der Einfluss der Protein-Kinase-C auf Die Regulation der Desmosomalen Proteine in Mauskeratinozyten. Ph.D. Thesis, Martin-Luther-University Halle-Wittenberg, Halle (Saale), Germany, 2021. [Google Scholar] [CrossRef]
  67. Keil, R.; Rietscher, K.; Hatzfeld, M. Antagonistic Regulation of Intercellular Cohesion by Plakophilins 1 and 3. J. Investig. Dermatol. 2016, 136, 2022–2029. [Google Scholar] [CrossRef] [PubMed]
  68. Fischer-Kešo, R.; Breuninger, S.; Hofmann, S.; Henn, M.; Röhrig, T.; Ströbel, P.; Stoecklin, G.; Hofmann, I. Plakophilins 1 and 3 Bind to FXR1 and Thereby Influence the mRNA Stability of Desmosomal Proteins. Mol. Cell. Biol. 2014, 34, 4244–4256. [Google Scholar] [CrossRef]
  69. Neuber, S.; Jäger, S.; Meyer, M.; Wischmann, V.; Koch, P.J.; Moll, R.; Schmidt, A. c-Src mediated tyrosine phosphorylation of plakophilin 3 as a new mechanism to control desmosome composition in cells exposed to oxidative stress. Cell Tissue Res. 2014, 359, 799–816. [Google Scholar] [CrossRef] [PubMed]
  70. López-Colomé, A.M.; Lee-Rivera, I.; Benavides-Hidalgo, R.; López, E. Paxillin: A crossroad in pathological cell migration. J. Hematol. Oncol. 2017, 10, 50. [Google Scholar] [CrossRef]
  71. Nah, A.S.; Chay, K.O. Roles of paxillin phosphorylation in IL-3 withdrawal-induced Ba/F3 cell apoptosis. Genes Genom. 2019, 41, 241–248. [Google Scholar] [CrossRef] [PubMed]
  72. Kim, H.; Rhee, S.H.; Pothoulakis, C.; LaMont, J.T. Clostridium difficile toxin A binds colonocyte Src causing dephosphorylation of focal adhesion kinase and paxillin. Exp. Cell Res. 2009, 315, 3336–3344. [Google Scholar] [CrossRef] [PubMed]
  73. Genth, H.; Pauillac, S.; Schelle, I.; Bouvet, P.; Bouchier, C.; Varela-Chavez, C.; Just, I.; Popoff, M.R. Haemorrhagic toxin and lethal toxin from Clostridium sordellii strain vpi9048: Molecular characterization and comparative analysis of substrate specificity of the large clostridial glucosylating toxins. Cell. Microbiol. 2014, 16, 1706–1721. [Google Scholar] [CrossRef]
  74. Brown, M.C.; Cary, L.A.; Jamieson, J.S.; Cooper, J.A.; Turner, C.E. Src and FAK Kinases Cooperate to Phosphorylate Paxillin Kinase Linker, Stimulate Its Focal Adhesion Localization, and Regulate Cell Spreading and Protrusiveness. Mol. Biol. Cell 2005, 16, 4316–4328. [Google Scholar] [CrossRef]
  75. Lee, J.-H.; Wittki, S.; Bräu, T.; Dreyer, F.S.; Krätzel, K.; Dindorf, J.; Johnston, I.C.; Gross, S.; Kremmer, E.; Zeidler, R.; et al. HIV Nef, Paxillin, and Pak1/2 Regulate Activation and Secretion of TACE/ADAM10 Proteases. Mol. Cell 2013, 49, 668–679. [Google Scholar] [CrossRef] [Green Version]
  76. Bartek, J.; Lukas, J. Chk1 and Chk2 kinases in checkpoint control and cancer. Cancer Cell 2003, 3, 421–429. [Google Scholar] [CrossRef]
  77. Colbran, R.J. Targeting of calcium/calmodulin-dependent protein kinase II. Biochem. J. 2004, 378, 1–16. [Google Scholar] [CrossRef] [PubMed]
  78. Kögel, D.; Reimertz, C.; Mech, P.; Poppe, M.; Frühwald, M.C.; Engemann, H.; Scheidtmann, K.H.; Prehn, J.H.M. Dlk/ZIP kinase-induced apoptosis in human medulloblastoma cells: Requirement of the mitochondrial apoptosis pathway. Br. J. Cancer 2001, 85, 1801–1808. [Google Scholar] [CrossRef] [PubMed]
  79. Najar, M.A.; Modi, P.K.; Ramesh, P.; Sidransky, D.; Gowda, H.; Prasad, T.S.K.; Chatterjee, A. Molecular Profiling Associated with Calcium/Calmodulin-Dependent Protein Kinase Kinase 2 (CAMKK2)-Mediated Carcinogenesis in Gastric Cancer. J. Proteome Res. 2021, 20, 2687–2703. [Google Scholar] [CrossRef]
  80. Naz, H.; Islam, A.; Ahmad, F.; Hassan, M.I. Calcium/calmodulin-dependent protein kinase IV: A multifunctional enzyme and potential therapeutic target. Prog. Biophys. Mol. Biol. 2016, 121, 54–65. [Google Scholar] [CrossRef] [PubMed]
  81. Lacher, M.D.; Pincheira, R.; Zhu, Z.; Camoretti-Mercado, B.; Matli, M.; Warren, R.S.; Castro, A.F. Rheb activates AMPK and reduces p27Kip1 levels in Tsc2-null cells via mTORC1-independent mechanisms: Implications for cell proliferation and tumorigenesis. Oncogene 2010, 29, 6543–6556. [Google Scholar] [CrossRef]
  82. Motoshima, H.; Goldstein, B.J.; Igata, M.; Araki, E. AMPK and cell proliferation—AMPK as a therapeutic target for atherosclerosis and cancer. J. Physiol. 2006, 574, 63–71. [Google Scholar] [CrossRef]
  83. Shen, Y.; Sherman, J.W.; Chen, X.; Wang, R. Phosphorylation of CDC25C by AMP-activated protein kinase mediates a metabolic checkpoint during cell-cycle G2/M-phase transition. J. Biol. Chem. 2018, 293, 5185–5199. [Google Scholar] [CrossRef]
  84. Gubern, A.; Joaquin, M.; Marquès, M.; Maseres, P.; Garcia-Garcia, J.; Amat, R.; González-Nuñez, D.; Oliva, B.; Real, F.X.; de Nadal, E.; et al. The N-Terminal Phosphorylation of RB by p38 Bypasses Its Inactivation by CDKs and Prevents Proliferation in Cancer Cells. Mol. Cell 2016, 64, 25–36. [Google Scholar] [CrossRef]
  85. Obaya, A.; Sedivy, J. Regulation of cyclin-Cdk activity in mammalian cells. Experientia 2002, 59, 126–142. [Google Scholar] [CrossRef]
  86. May, M.; Schelle, I.; Brakebusch, C.; Rottner, K.; Genth, H. Rac1-dependent recruitment of PAK2 to G2phase centrosomes and their roles in the regulation of mitotic entry. Cell Cycle 2014, 13, 2210–2220. [Google Scholar] [CrossRef]
  87. Adhikari, B.; Bozilovic, J.; Diebold, M.; Schwarz, J.D.; Hofstetter, J.; Schröder, M.; Wanior, M.; Narain, A.; Vogt, M.; Stankovic, N.D.; et al. PROTAC-mediated degradation reveals a non-catalytic function of AURORA-A kinase. Nat. Chem. Biol. 2020, 16, 1179–1188. [Google Scholar] [CrossRef] [PubMed]
  88. Genth, H.; Hofmann, F.; Selzer, J.; Rex, G.; Aktories, K.; Just, I. Difference in Protein Substrate Specificity between Hemorrhagic Toxin and Lethal Toxin from Clostridium sordellii. Biochem. Biophys. Res. Commun. 1996, 229, 370–374. [Google Scholar] [CrossRef] [PubMed]
  89. Stieglitz, F.; Gerhard, R.; Pich, A. The Binary Toxin of Clostridioides difficile Alters the Proteome and Phosphoproteome of HEp-2 Cells. Front. Microbiol. 2021, 12, 725612. [Google Scholar] [CrossRef] [PubMed]
  90. Cox, J.; Mann, M. MaxQuant enables high peptide identification rates, individualized p.p.b.-range mass accuracies and proteome-wide protein quantification. Nat. Biotechnol. 2008, 26, 1367–1372. [Google Scholar] [CrossRef] [PubMed]
  91. Cox, J.; Michalski, A.; Mann, M. Software Lock Mass by Two-Dimensional Minimization of Peptide Mass Errors. J. Am. Soc. Mass Spectrom. 2011, 22, 1373–1380. [Google Scholar] [CrossRef]
  92. The UniProt Consortium. UniProt: The universal protein knowledgebase in 2021. Nucleic Acids Res. 2021, 49, D480–D489. [Google Scholar] [CrossRef]
  93. Tyanova, S.; Temu, T.; Sinitcyn, P.; Carlson, A.; Hein, M.Y.; Geiger, T.; Mann, M.; Cox, J. The Perseus computational platform for comprehensive analysis of (prote)omics data. Nat. Methods 2016, 13, 731–740. [Google Scholar] [CrossRef]
  94. Krämer, A.; Green, J.; Pollard, J., Jr.; Tugendreich, S. Causal analysis approaches in Ingenuity Pathway Analysis. Bioinformatics 2014, 30, 523–530. [Google Scholar] [CrossRef]
  95. Perez-Riverol, Y.; Bai, J.; Bandla, C.; García-Seisdedos, D.; Hewapathirana, S.; Kamatchinathan, S.; Kundu, D.J.; Prakash, A.; Frericks-Zipper, A.; Eisenacher, M.; et al. The PRIDE database resources in 2022: A hub for mass spectrometry-based proteomics evidences. Nucleic Acids Res. 2021, 50, D543–D552. [Google Scholar] [CrossRef]
Figure 1. Morphological alterations of HEp-2 cells after 24 h of incubation with TcnA (Magnification 100). Ctrl: untreated cells, FA_TcnA: cells treated with inactivated TcnA 255 ng/mL, TcnA: cells treated with 255 ng/mL TcnA. White arrows indicating cell strands ranging into the medium but still attached to cells adhering to the flask.
Figure 1. Morphological alterations of HEp-2 cells after 24 h of incubation with TcnA (Magnification 100). Ctrl: untreated cells, FA_TcnA: cells treated with inactivated TcnA 255 ng/mL, TcnA: cells treated with 255 ng/mL TcnA. White arrows indicating cell strands ranging into the medium but still attached to cells adhering to the flask.
Ijms 23 09939 g001
Figure 2. Statistical analyses and quality control of effects of TcnA on the proteome of HEp-2 cells. (A) Heat map of significantly changed proteins based on BJH-corrected and FDR-based ANOVA testing of the protein abundancies. (B) Heat map of significantly changed phosphosites based on BJH-corrected and FDR-based ANOVA testing of the protein abundancies. Dendrograms show similarity of samples/treatments.
Figure 2. Statistical analyses and quality control of effects of TcnA on the proteome of HEp-2 cells. (A) Heat map of significantly changed proteins based on BJH-corrected and FDR-based ANOVA testing of the protein abundancies. (B) Heat map of significantly changed phosphosites based on BJH-corrected and FDR-based ANOVA testing of the protein abundancies. Dendrograms show similarity of samples/treatments.
Ijms 23 09939 g002
Figure 3. Effects of TcnA on the proteome of HEp-2 cells. (A) Volcano plot of proteins from HEp-2 cells with TcnA and FA_TcnA. Dashed lines give the minimal regulation factor (log2 difference > 1; <−1) and the significance level −log10 (0.05). (B) Fisher’s exact test of significantly regulated proteins after BJH corrected t-testing. Enrichment factor > 1 indicates enrichment of the term, while an enrichment factor < 1 marks a reduction of this term within the significant proteins. (C) IPATM Pathway analysis of significant proteins (p-value < 0.05). A positive z-scores indicates predicted activation of the pathway, while a negative z-score indicates predicted inhibition of the pathway. Size of the points correlates to the −log10 (p-value).
Figure 3. Effects of TcnA on the proteome of HEp-2 cells. (A) Volcano plot of proteins from HEp-2 cells with TcnA and FA_TcnA. Dashed lines give the minimal regulation factor (log2 difference > 1; <−1) and the significance level −log10 (0.05). (B) Fisher’s exact test of significantly regulated proteins after BJH corrected t-testing. Enrichment factor > 1 indicates enrichment of the term, while an enrichment factor < 1 marks a reduction of this term within the significant proteins. (C) IPATM Pathway analysis of significant proteins (p-value < 0.05). A positive z-scores indicates predicted activation of the pathway, while a negative z-score indicates predicted inhibition of the pathway. Size of the points correlates to the −log10 (p-value).
Ijms 23 09939 g003
Figure 4. Regulated pathways (p-value < 0.05) in TcdA [21], TcdB [22] and TcnA treated HEp-2 cells. Positive z-score indicates activation (z-score > 2 significant activation), while a negative z-score indicates inhibition of the pathway (z-score < −2 significant inhibition).
Figure 4. Regulated pathways (p-value < 0.05) in TcdA [21], TcdB [22] and TcnA treated HEp-2 cells. Positive z-score indicates activation (z-score > 2 significant activation), while a negative z-score indicates inhibition of the pathway (z-score < −2 significant inhibition).
Ijms 23 09939 g004
Figure 5. Effects of TcnA on the phosphoproteome of HEp-2 cells. (A) Volcano plot of proteins from HEp-2 cells with TcnA and FA_TcnA. Dashed lines give the minimal regulation factor (log2 difference > 1; <−1) and the significance level −log10 (0.05). (B) IPATM Pathway analysis of significantly altered phosphosites (p-value < 0.05). Size of the points correlates to the −log10 (p-value). (C) Kinase activation prediction. Top 10 regulated kinases identified by Fisher’s exact test (p-value < 0.05) and activation status prediction based on phosphosite (p < 0.05) regulation.
Figure 5. Effects of TcnA on the phosphoproteome of HEp-2 cells. (A) Volcano plot of proteins from HEp-2 cells with TcnA and FA_TcnA. Dashed lines give the minimal regulation factor (log2 difference > 1; <−1) and the significance level −log10 (0.05). (B) IPATM Pathway analysis of significantly altered phosphosites (p-value < 0.05). Size of the points correlates to the −log10 (p-value). (C) Kinase activation prediction. Top 10 regulated kinases identified by Fisher’s exact test (p-value < 0.05) and activation status prediction based on phosphosite (p < 0.05) regulation.
Ijms 23 09939 g005
Table 1. Regulated pathways (p < 0.05) identified in TcdA, TcdB and TcnA treated HEp-2 cells. Pathways sorted according to regulation following TcnA treatment.
Table 1. Regulated pathways (p < 0.05) identified in TcdA, TcdB and TcnA treated HEp-2 cells. Pathways sorted according to regulation following TcnA treatment.
Pathwayz-Score
TcdATcdBTcnA
Signaling by Rho Family GTPases0.821.073.29
Epithelial Adherens Junction Signaling0.381.263.67
IL-8 Signaling0.822.112.12
Protein Kinase A Signaling0.821.92.11
Actin Cytoskeleton Signaling0.331.91.9
Regulation of Actin-based Motility by Rho0.8211.71
Ephrin Receptor Signaling−1.131.671.46
Leukocyte Extravasation Signaling−1.131.671.3
mTOR Signaling0.451.411.22
Integrin Signaling−0.711.941.15
Hepatic Fibrosis Signaling Pathway−0.631.291.07
ERK/MAPK Signaling−12.241.03
Cholecystokinin/Gastrin-mediated Signaling−0.821.670.73
Cardiac Hypertrophy Signaling0.821.670.67
Actin Nucleation by ARP-WASP Complex0.4510.53
Paxillin Signaling−1.340.820.45
PAK Signaling−11.630.22
HER-2 Signaling in Breast Cancer−11.63−0.18
Glioma Invasiveness Signaling0.450.82−0.24
HGF Signaling−1.341.41−0.24
ERBB Signaling−21.34−0.77
HIPPO signaling−0.82−1−1.15
RHOGDI Signaling−1.41−1−2.12
Table 2. Top 10 up- and down-regulated proteins after TcdA, TcdB and TcnA treatment. Bold: proteins that are regulated in more than one treatment condition.
Table 2. Top 10 up- and down-regulated proteins after TcdA, TcdB and TcnA treatment. Bold: proteins that are regulated in more than one treatment condition.
ExperimentUp-RegulatedDown-Regulated
Gene NameProtein IDp-ValueLog2 Fold-ChangeGene NameProtein IDp-ValueLog2 Fold-Change
TcnAIGFN1Q86VF2-59.9 × 10−42.58C7orf25Q9BPX74.8 × 10−3−2.34
TcnANKTRP304146.4 × 10−42.37TNNI2P48788-22.7 × 10−2−2.24
TcnAFBRSL1Q9HCM71.2 × 10−32.36ABCD2Q9UBJ21.7 × 10−2−2.05
TcnANPY4RP0DQD54.3 × 10−42.33SMU1Q2TAY7-21.5 × 10−2−1.53
TcnAINTS1Q8N2016.8 × 10−42.22C1QTNF6Q9BXI9-12.8 × 10−2−1.50
TcnAMTCL1Q9Y4B5-22.7 × 10−32.11MIFP141743.5 × 10−3−1.37
TcnATAOK1Q7L7X35.0 × 10−42.05CCDC129Q6ZRS4-22.7 × 10−2−1.33
TcnAMYLPFQ96A321.3 × 10−31.81SCAF4O95104-29.7 × 10−3−1.33
TcnAC12orf29Q8N999-34.8 × 10−31.45NUDT9Q9BW91-22.5 × 10−2−1.25
TcnASLC6A12P480658.8 × 10−41.42HSP90AB2PQ58FF82.2 × 10−2−1.25
TcdBRHOBP627451.4 × 10−22.25CDC42EP1Q005872.1 × 10−3−2.48
TcdBCTSLP077111.1 × 10−41.94RND3P615874.5 × 10−3−1.63
TcdBJUNP054125.6 × 10−41.73YAP1P469377.6 × 10−6−1.55
TcdBHMOX1P096012.0 × 10−20.83CDC42EP2O146131.2 × 10−2−1.30
TcdBMYO10Q9HD671.8 × 10−20.77CYR61O006227.2 × 10−3−1.23
TcdBSDCBPO005603.2 × 10−40.77TPM2Q5TCU31.3 × 10−2−1.05
TcdBSEC14L1Q925036.7 × 10−30.76TPM2P07951-21.6 × 10−2−0.86
TcdBSDCBPO00560-22.2 × 10−20.71TPM3Q5VU638.1 × 10−3−0.74
TcdBSTAU1O957932.5 × 10−20.63F3P137263.3 × 10−2−0.69
TcdBRASSF5Q8WWW03.3 × 10−20.61SLC1A3P430031.9 × 10−2−0.67
TcdARHOBP627451.8 × 10−42.72YAP1P469374.5 × 10−4−1.71
TcdAHMOX1P096014.4 × 10−20.93TPM2P07951-21.2 × 10−5−1.14
TcdAAPPE9PEV01.0 × 10−30.67TPM1P09493-37.6 × 10−4−1.07
TcdATIMP3B1AJV73.4 × 10−20.61CLDN1O958325.2 × 10−3−0.78
TcdATM4SF1F8WF272.5 × 10−30.61SLC1A3P43003-23.6 × 10−2−0.74
TcdATUBB2AQ138853.3 × 10−20.55AFAP1Q8N5562.6 × 10−3−0.73
TcdACD63F8VNT96.7 × 10−30.51FHL2J3KNW42.1 × 10−3−0.63
TcdARRM2P313507.2 × 10−30.51MPRIPQ6WCQ11.9 × 10−4−0.51
TcdAS100A16Q96FQ61.4 × 10−30.49STAT3K7ENL34.2 × 10−3−0.46
TcdAITGA5P086485.0 × 10−30.46PRKCIP417433.1 × 10−2−0.46
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Schweitzer, T.; Genth, H.; Pich, A. Clostridium novyi’s Alpha-Toxin Changes Proteome and Phosphoproteome of HEp-2 Cells. Int. J. Mol. Sci. 2022, 23, 9939. https://doi.org/10.3390/ijms23179939

AMA Style

Schweitzer T, Genth H, Pich A. Clostridium novyi’s Alpha-Toxin Changes Proteome and Phosphoproteome of HEp-2 Cells. International Journal of Molecular Sciences. 2022; 23(17):9939. https://doi.org/10.3390/ijms23179939

Chicago/Turabian Style

Schweitzer, Theresa, Harald Genth, and Andreas Pich. 2022. "Clostridium novyi’s Alpha-Toxin Changes Proteome and Phosphoproteome of HEp-2 Cells" International Journal of Molecular Sciences 23, no. 17: 9939. https://doi.org/10.3390/ijms23179939

APA Style

Schweitzer, T., Genth, H., & Pich, A. (2022). Clostridium novyi’s Alpha-Toxin Changes Proteome and Phosphoproteome of HEp-2 Cells. International Journal of Molecular Sciences, 23(17), 9939. https://doi.org/10.3390/ijms23179939

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop