Next Article in Journal
Escherichia coli Alpha-Hemolysin HlyA Induces Host Cell Polarity Changes, Epithelial Barrier Dysfunction and Cell Detachment in Human Colon Carcinoma Caco-2 Cell Model via PTEN-Dependent Dysregulation of Cell Junctions
Next Article in Special Issue
Systemic Effects of Hemorrhagic Snake Venom Metalloproteinases: Untargeted Peptidomics to Explore the Pathodegradome of Plasma Proteins
Previous Article in Journal
Seasonal Variability of the Airborne Eukaryotic Community Structure at a Coastal Site of the Central Mediterranean
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Toxins
Volume 13
Issue 8
10.3390/toxins13080519
toxins-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Bothrops Jararaca Snake Venom Modulates Key Cancer-Related Proteins in Breast Tumor Cell Lines

by
Carolina Yukiko Kisaki
1,
Stephanie Santos Suehiro Arcos
1,
Fabio Montoni
1,
Wellington da Silva Santos
1,
Hamida Macêdo Calacina
1,
Ismael Feitosa Lima
1,
Daniela Cajado-Carvalho
1,
Emer Suavinho Ferro
2,
Milton Yutaka Nishiyama-Jr
1,* and
Leo Kei Iwai
1,*
1
Laboratory of Applied Toxinology (LETA) and Center of Toxins, Immune-Response and Cell Signaling (CeTICS), Butantan Institute, São Paulo 05503-900, Brazil
2
Department of Pharmacology, Biomedical Sciences Institute (ICB), University of São Paulo (USP), São Paulo 05508-000, Brazil
*
Authors to whom correspondence should be addressed.
Toxins 2021, 13(8), 519; https://doi.org/10.3390/toxins13080519
Submission received: 5 July 2021 / Accepted: 19 July 2021 / Published: 25 July 2021
(This article belongs to the Special Issue Celebrating 120 Years of Butantan Institute Contributions for Toxinology)
Browse Figures
Versions Notes

Abstract

:
Cancer is characterized by the development of abnormal cells that divide in an uncontrolled way and may spread into other tissues where they may infiltrate and destroy normal body tissue. Several previous reports have described biochemical anti-tumorigenic properties of crude snake venom or its components, including their capability of inhibiting cell proliferation and promoting cell death. However, to the best of our knowledge, there is no work describing cancer cell proteomic changes following treatment with snake venoms. In this work we describe the quantitative changes in proteomics of MCF7 and MDA-MB-231 breast tumor cell lines following treatment with Bothrops jararaca snake venom, as well as the functional implications of the proteomic changes. Cell lines were treated with sub-toxic doses at either 0.63 μg/mL (low) or 2.5 μg/mL (high) of B. jararaca venom for 24 h, conditions that cause no cell death per se. Proteomics analysis was conducted on a nano-scale liquid chromatography coupled on-line with mass spectrometry (nLC-MS/MS). More than 1000 proteins were identified and evaluated from each cell line treated with either the low or high dose of the snake venom. Protein profiling upon venom treatment showed differential expression of several proteins related to cancer cell metabolism, immune response, and inflammation. Among the identified proteins we highlight histone H3, SNX3, HEL-S-156an, MTCH2, RPS, MCC2, IGF2BP1, and GSTM3. These data suggest that sub-toxic doses of B. jararaca venom have potential to modulate cancer-development related protein targets in cancer cells. This work illustrates a novel biochemical strategy to identify therapeutic targets against cancer cell growth and survival.
Keywords:
mass spectrometry; proteome; snake venom; Bothrops jararaca; breast cancer
Key Contribution: We describe that B. jararaca snake venom modulates specific protein pathways related to cancer cell growth and invasion, which could be a useful approach to identify novel therapeutic targets for cancer treatment.

Graphical Abstract

1. Introduction

Ophidian accidents constitute a serious public health problem in Brazil, with an average of 29,000 cases and 125 deaths reported every year (Brazilian Ministry of Health, 2019) [1]. Approximately 80% of ophidian accidents are caused by snakes of the Viperidae family, more specifically of the Bothrops genus [2]. Among them, about 25% lead to death or sequels capable of generating temporary or permanent incapacity for work and customary activities. Venom from the Bothrops jararaca (B. jararaca) snake is a complex mixture composed of proteins, peptides, amino acids, nucleotides, lipids, and carbohydrates that present a range of different actions when they are isolated or together [3,4,5,6,7], leading to hemotoxic, cardiotoxic, cytotoxic, or neurotoxic effects [8,9]. Several previous reports have defined the proteomics composition of Bothrops venoms [5,10,11,12,13,14,15,16]. These studies have shown that Bothrops venoms are composed of various classes of toxin, including metalloproteinases, serine proteinases, phospholipases A2, and C-type lectins, the most abundant components participating in the local and systemic envenomation effects.
The venom of B. jararaca engenders three main activities: proteolytic, coagulant, and hemolytic. The proteolytic activity causes degradation of extracellular matrix proteins, plasma, and cell surface [17], which represents an important factor for the clinical characterization of a bothropic accident [2,18]. In addition, venom can cause local tissue lesion, myonecrosis, edema, cardiovascular alterations, hypovolemic shock, coagulation alteration and renal alterations, resulting from the combined action of the enzymatic and toxic activity of the venom [19]. About 90 to 95% of the dry weight of the B. jararaca venom is composed of a complex mixture of proteins, mainly metalloproteinases, serine proteinases, phospholipases (PLA2), and L-amino acid oxidases. The metalloproteinases comprise most of the venom composition [20,21]. They are proteolytic enzymes associated with fibrinolysis and coagulation, and they are involved in cell migration and tissue repair, besides being related to pathological effects such as cancer [22,23]. In terms of therapeutic interventions, protease inhibitors have been shown to inhibit homeostasis and thrombosis by acting on the coagulation cascade [24]. The third major component of the venom, PLA2, is an enzyme capable of hydrolyzing the ester bonds at the sn2 position of glycerolphospholipids, releasing arachidonic acid, important for the biosynthesis of many mediators involved in inflammation, such as prostaglandins, thromboxanes, and leukotrienes [25]. Finally, the L-amino acid oxidases (LAAOs), which make up about 1 to 9% of the venom composition [26], are flavoenzymes belonging to the class of oxidoreductases, which produce alpha-keto acid, hydrogen peroxide, and ammonia [27,28]. However, when there is a high production of hydrogen peroxide, it has been found that L-amino acid oxidases can induce apoptosis in mammalian endothelial cells [29].
Snake venom constituents have been isolated and studied for their therapeutic potential in the treatment of various diseases. One example is Eptifibatide, marketed as Integrilin, derived from the Echis carinatus snake venom and produced by Millennium Pharmaceuticals and Schering-Plow. It is used as an antiplatelet drug [30]. Another example is the angiotensin I converting enzyme inhibitor Captopril produced by Bristol-Myers Squibb whose active component was derived from B. jararaca venom. It is used for the treatment against hypertension and renal insufficiency [31,32]. In addition to the potential use of the derivatives of snake venom toxins in the treatment of non-malignant diseases [33], several studies have described anti-tumorigenic characteristics of snake venom, stating that snake venom may be capable of inhibiting cell proliferation and promoting cell death by different means: inducing apoptosis in cancer cells by increasing the influx of Ca2+, inducing the release of cytochrome C, decreasing or increasing the expression of proteins that control the cell cycle, and causing damage to cell membranes [26,34,35,36]. With the goal of searching novel therapy against cancer, studies have characterized the proteins, peptides or enzymes derived from snake venom to identify components that are capable of interfering with the transport of substances or signal transduction across the membrane or disrupting the cell membrane [35,37].
With the rapid advances of nano-scale liquid chromatography (nLC) and mass spectrometry (MS) technologies in the last two decades, nLC-MS/MS-based proteomics analysis has been widely applied as a powerful tool for biomarker discovery to improve cancer therapy [38,39]. A number of studies have characterized the biochemical and physiological action of venom or isolated venom derivatives on cell lines or tissues [40,41,42,43]. In addition, several works have shown proteomic changes of cancer cell lines upon drug treatment that suggest molecular mechanisms of drug action, including diverse effects on proteasome regulation, metabolic processes, and oxidative stress [44,45].
However, to the best of our knowledge, there has been no report that describes the effects of B. jararaca snake venom treatment on breast cancer-related cell proteome. In the present study, nLC-MS/MS was used to characterize the effects of sub-toxic doses of B. jararaca snake venom on two different breast cancer cell lines MCF7 and MDA-MB-231. MCF7 and MDA-MB-231 are non-metastatic and metastatic tumor cell lines, respectively. They are characterized by a high degree of glycolytic efficiency that promotes the interaction between the tumor cell and the extracellular matrix [46]. Although both cell lines are from breast origin, they are molecularly distinct. MCF7 are estrogen and progesterone receptors positive and HER2 negative, while MDA-MB-231 are triple negative (estrogen receptor, progesterone receptor, and HER2 negatives) and prone to cytotoxic agents because of their impaired DNA repairing capability which is in part due to mutation in the p53 gene [47,48]. Proteomic changes observed herein upon treatment with B. jararaca snake venom in these cell lines highlight proteins and cell pathways that could be targeted in cancer therapy.

2. Results

2.1. The Cytotoxicity of B. jararaca Snake Venom in MCF7 and MDA-MB-231 Cells

The B. jararaca venom cytotoxicity assay on MCF7 and MDA-MB-231 cell lines was monitored using the WST-1 reagent. This analysis showed that although cell viability was similar between both MCF7 and MDA-MB-231 cell lines, they had different venom resistance profiles where MDA-MB-231 cells showed to be more resistant to the venom when compared to the MCF7 cells. Although both cell lines started to die at doses higher than 2.5 μg/mL, at the 5.0 μg/mL of venom, only about half of the MDA-MB-231 cells have died while most of the MCF7 have died at this venom concentration (Figure 1). Lethal concentration 50 (LC50) was determined as 4.50 μg/mL for MCF7 and 4.76 μg/mL for MDA-MB-231 cell line. Interestingly, the treatment of cells with the low dose at 0.5 μg/mL of venom also killed more MCF7 cells compared to the MDA-MB-231 cells when compared to the venom doses of 0.63 μg/mL and 1.25 μg/mL (Figure 1).

2.2. Optical Microscopy Analysis of MCF7 and MDA-MB-231 Cells under B. jararaca Venom Treatment

Optical microscopy analysis at 10× magnification of MCF7 and MDA-MB-231 cell lines treated with concentrations higher than 2.5 µg/mL of B. jararaca snake venom showed cellular morphological changes such as cell shrinkage and cell birefringence change (Figure S1). At the 20 µg/mL of venom treatment all MCF7 cells detached from the plate, whereas the MDA-MB-231 cell line continued to show morphological death-like changes, but the cells did not detach from the plate (Figure S1).
Based on the cytotoxicity assays and visualization of cell morphology changes through the microscope images, two working concentrations, representing a low dose of 0.63 µg/mL and a high sub-toxic dose of 2.5 µg/mL of venom, were selected for further proteomics analysis.

2.3. Mass Spectrometry-Based Proteomics of MCF7 and MDA-MB-231 Cells Treated with B. jararaca Venom

Protein identification was performed analyzing the raw data using the MaxQuant software against the Homo sapiens database downloaded from Uniprot. Analysis of the MCF7 cell line treated with different venom concentrations allowed the identification of 789 proteins from which 704 proteins were identified in the control group (non-venom treatment), 725 proteins were identified in cells treated with 0.63 µg/mL of venom, and 713 proteins were identified in the cells treated with 2.5 µg/mL of venom. Comparative analysis of the identified proteins showed that among all of the 789 proteins identified, 637 proteins were in common to all conditions, 657 proteins were in common between the control and the 0.63 µg/mL venom treatment, 656 proteins were in common between control and the 2.5 µg/mL venom treatment, and 677 proteins were in common between the 0.63 µg/mL and 2.5 µg/mL venom treatments (Figure 2a, Table S1a). In addition, we also identified exclusive proteins in each condition: 28 proteins in the control group, 28 proteins in 0.63 µg/mL venom treatment, and 17 proteins exclusive in the 2.5 µg/mL venom treatment (Figure 2a, Table S1a).
In the MDA-MB-231 cell line, mass spectrometry-based proteomics analysis allowed us to identify a total of 1093 proteins from which 1006 proteins in the control group, 992 proteins in the cells treated with 0.63 µg/mL of venom, and 978 proteins in the cells treated with 2.5 µg/mL of venom. Comparative analysis of the identified proteins showed 893 proteins identified in common to all conditions, 949 proteins in common between the control group and 0.63 µg/mL venom treatment, 912 proteins were in common between control group and 2.5 µg/mL venom treatment, and 915 proteins were in common between the 0.63 µg/mL and 2.5 µg/mL venom treatments (Figure 2a, Table S1a). We also observed 38 exclusive proteins in the control group, 21 proteins in 0.63 µg/mL venom treatment, and 44 proteins exclusive in the 2.5 µg/mL venom treatment (Figure 2b, Table S1b).

2.4. Semi-Quantitative Proteomics Analysis: MCF7 and MDA-MB-231 Cell Line Protein Abundance Variation

In general, we observed a higher number of proteins whose abundance had changed more than 1.5× or less than 0.67× in MCF7 cell lines when compared to the MDA-MB-231 cell lines. The semi-quantitative analysis of the MCF7 cell line (Table S2a,) treated with 2.5 µg/mL venom allowed us to identify 137 proteins, whose abundance changed over 1.5× (fold change FC ≥ 1.5, marked in light red) when compared to the control group, from which 55 proteins presented FC ≥ 2.0 (marked in red). We highlight 12 highly abundant proteins with FC ≥ 3.0, marked in dark red: Sorting nexin-3 (SNX3), Purine nucleoside phosphorylase (HEL-S-156an), Peroxisome proliferator activated receptor interacting complex protein (PRIC295), Small nuclear ribonucleoprotein component (SNRP116), Eukaryotic translation initiation factor 4B (EIF4B), Methylcrotonoyl-CoA carboxylase beta chain (MCCC2), 26S proteasome non-ATPase regulatory subunit 5 (PSMD5), Heterogeneous nuclear ribonucleoprotein R (HNRNPR), Full-length cDNA clone CS0DJ015YJ12 of T cells (PSME2), Isoleucyl-tRNA synthetase (IARS), Large proline-rich protein BAG6 (BAG6), and Glutathione S-transferase (GSTM3). In addition, we identified 23 proteins with FC ≤ 0.67 (marked in light green) from which five proteins presented FC ≤ 0.5 (marked in dark green): Histone H4 (HIST1H4J), ATP synthase subunit d, mitochondrial (ATP5PD), Voltage-dependent anion-selective channel protein 2 (VDAC2), 4a-hydroxytetrahydrobiopterin dehydratase (PCBD), and Histone H3 (H3F3B). At the low 0.63 µg/mL venom treatment, we identified 25 proteins with FC ≥ 1.5 (marked in light red) from which only two proteins presented FC ≥ 2.0 (marked in red): Anterior gradient 2 homolog (AGR2) and Leucine-rich PPR-motif containing protein (LRPPRC); and 19 proteins with FC ≤ 0.67 (marked in light green) from which two proteins with FC ≤ 0.5 (marked in dark green): PCBD and 40S ribosomal protein S29 (RPS29). The description of the highlighted proteins is shown in Table 1.
Semi-quantitative proteomic analysis of the MDA-MB-231 cell line (Table S2b) treated with 2.5 µg/mL allowed the identification of 34 proteins whose abundance changed (FC) ≥1.5 over the control (marked in light red), nine proteins with FC ≥ 2 (marked in red) from which we highlight three proteins with FC ≥ 3 (marked in dark red): Histone H3.2 (H3C15/HIST2H3), 14 kDa phosphohistidine phosphatase (HEL-S-132P), and Mitochondrial carrier homolog 2 (MTCH2). Moreover, 41 proteins presented FC ≤ 0.67 (light green) from which we highlight four proteins with FC ≤ 0.5 (marked in dark green): DnaJ homolog subfamily A member 1 (DNAJA1), Insulin-like growth factor 2 mRNA-binding protein 1 (IGF2BP1), Cysteine-rich angiogenic inducer 61 (CYR61), and Thrombospondin-1 (THBS1). At the lower 0.63 µg/mL venom treatment, 16 proteins presented FC ≥ 1.5 (marked in light red) from which we highlight H3C15/HIST2H3 with FC = 3.8 and MTCH2 with FC = 2.3 (marked in red); and 28 proteins with FC ≤ 0.67 (marked in light green) from which 12 proteins with FC ≤ 0.5 (marked in dark green): 60S ribosomal protein L37 (RPL37), D-3-phosphoglycerate dehydrogenase (HEL-S-113), ATPase inhibitor, mitochondrial (ATP5IF1), Non-histone chromosomal protein HMG-14 (HMGN1), RCC2 protein (RCC2), Serine/threonine-protein phosphatase PP1-gamma catalytic subunit (PPP1CC), D-dopachrome decarboxylase (DDT), Ran GTPase-activating protein 1 (RANGAP1), dCTP pyrophosphatase 1 (DCTPP1), IGF2BP1, THBS1, and CYR61 (Table 1).

2.5. Hierarchical Clustering Analysis

Hierarchical clustering analysis of differentially expressed proteins with FC ≥ 1.5 identified in both MCF7 and MDA-MB-231 cell lines characterized those proteins in seven major clusters (Figure 3, Table S3). Clusters 1–4 had proteins with fold change (FC) increased at 2.5 µg/mL of venom treatment in MCF7 cells and presented no change in MDA-MB-231 cell lines. Proteins in these clusters include HEL-S-156 and PRIC295 in cluster 1 and PSMD5 and PSME2 in cluster 2. Cluster 5 identified proteins that decreased in FC at 2.5 µg/mL of venom treatment in both cell lines including HIST1H4J, VDAC1, and VDAC2. Cluster 6 identified proteins that increased FC at 2.5 µg/mL of venom treatment in MDA-MB-231 cell line including H3F3B, LAP3, and KRT1, and cluster 7 identified proteins that did not change when cells were treated with low and high venom treatment, but they presented higher FC change in MDA-MB-231 when compared to MCF7 cell line such as RPS29 (Table 1).

2.6. Principal Component Analysis

PCA was applied to the differentially expressed proteins identified from both MCF7 and MDA-MB-231 cell lines based on the log2 FC of cells treated with low and high B. jararaca venom compared to the PBS treatment control group (Figure 4). The projection into the component space shows a distinct coordinated activity of proteins between the cell lines conditions. Orthogonal vectors show highly positive correlation between both venom concentrations in MDA-MB-231 cell line which may represent similar cell line responses to the different venom concentrations, but they present a highly negative correlation to both venom concentrations in MCF7 and respective set of expressed proteins. In the MCF7 cell line, however, we observe a negative correlation between the 0.63 μg/mL and 2.5 μg/mL venom treatment, indicating a differential response upon low and high dose venom treatment. PCA also shows differential correlation between both cell lines where the first two components showed 40% variation in the PC1 and 26.7% variation in the PC2 between the MCF7 and MDA-MB-231 cell lines. In addition, we observed clusters of proteins positively correlating with both low and high venom treatment in MDA-MB-231 cell line such as LAP3, H3F3B, and KRT1. On the other hand, we also observed proteins such as H2AC20 and TUFM correlating with low venom treatment, and proteins such as HEL-S-156an correlating with high venom treatment on MCF7 cell line (Table 1).

2.7. Gene Ontology Functional Analysis

The most enriched protein families and functional categories were analyzed based on highly abundant proteins with FC ≥ 1.5 for each cell line and treatment at low and high B. jararaca venom conditions (Figures S2–S5). The functional ontology classification analysis of these sets of proteins showed that both MCF7 and MDA-MB-231 venom-treated cell lines showed similar enriched categories. In addition, the most prominent enrichment was identified for treated cells with the sub-toxic dose of 2.5 µg/mL of venom. The molecular function enrichment analysis in both the MCF7 and MDA-MB-231 cell lineages showed an enrichment of proteins related to binding, structural molecule activity, and catalytic activity. In addition, the MCF7 cells had enriched, albeit in a lower amount, proteins related to function and transcriptional regulatory activity and carrier activity (Figure S2). The functional classification analysis related to biological processes showed enriched proteins related to the metabolic process and the cellular component organization or biogenesis. Moreover, the analysis of proteins identified in the MCF7 cell line presented proteins related to the cellular process, localization, biological regulation, stimulus response, developmental process, multicellular organismal process, and the immune system process (Figure S3). The analysis of protein distribution by cellular components showed an enrichment related to the “cell”, protein complex, and organelle (Figure S4), and the enrichment analysis of protein family classification showed an enrichment of cytoskeleton proteins, ligase, nucleic acid binding, signaling molecule, modulating enzyme, calcium binding protein, and hydrase (Figure S5).

2.8. Protein–Protein Interaction Analysis

STRING protein–protein interaction network analysis tool was used to evaluate protein–protein interactions identified among proteins with FC ≥ 1.5 from each cell line treated with B. jararaca venom. According to Doncheva and colleagues [114], STRING indicates interactions according to co-expression analyzes and evolutionary signals in all genomes based on data described in the literature between genes or proteins using functional classification systems such as Gene Ontology, KEGG (Kyoto Encyclopedia of Genes and Genomes), and Reactome.
Analysis of protein–protein interactions (PPIs) of proteins whose abundance increased more than 1.5× in the MCF7 cell line treated with 2.5 μg/mL of venom showed a high interconnection of proteins related to metabolic process (in red) and metabolism pathways (in blue) (Figure 5a). Interestingly, among the proteins identified in the MCF7, we observed clusters of highly connected proteins related to proteasome pathway (in green) and mRNA splicing (in yellow). We highlight the proteins PSMD2 and PSMD11 (26S proteasome non-ATPase regulatory subunit 2 and 11, respectively), PSME1 and PSME2 (Proteasome activator complex subunit 1 and 2, respectively), PSMC1, PSMC4, and PSMC5 (26S proteasome regulatory subunit 4, 6B, and 8, respectively), TXNL (Thioredoxin-like protein 1), and USP14 (Ubiquitin carboxyl-terminal hydrolase 14), which are all members of the proteasome pathway (in green) and proteasome complex, and, together with MST4 (Serine/threonine-protein kinase 26) and PSMD5 (26S proteasome non-ATPase regulatory subunit 5) they are all related to apoptosis (Figure 5a). PPI analysis of proteins identified in the MDA-MB-231 in the same conditions showed less interaction among the proteins identified with FC ≥ 1.5 (Figure 5b).
PPI analysis of the MCF7 cell lines with FC ≥ 1.5 treated with 0.63 μg/mL and FC ≤ 0.67 when cells were treated with 0.63 μg/mL and 2.5 μg/mL of venom are shown in Figure S6. Additionally, PPI analysis of proteins identified in the MDA-MB-231 at 0.63 μg/mL and 2.5 μg/mL venom treatment presenting FC ≤ 0.67 and FC ≥ 2.5 are shown in Figure S7.

2.9. Exclusive Proteins

We further analyzed proteins exclusively identified in one or two of the three conditions in both MCF7 cell line (Table S4a, Figure S8) and MDA-MB-231 cell line (Table S4b, Figure S8). Proteins identified only on PBS treated cells (control group), suggesting that the expression of the proteins was inhibited with the addition of venom, allowed us to identify 28 proteins in the MCF7 cell line from which we highlight proteins: eukaryotic peptide chain release factor subunit 1 (ETF1), serine/arginine repetitive matrix protein 2 (SRRM2), PHD finger-like domain containing protein 5A (PHF5A), and lamina-associated polypeptide 2 (TMPO). The analysis of the MDA-MB-231 cell line allowed us to identify 41 proteins exclusively expressed at the control group from which we highlight Plasminogen activator inhibitor 1 (SERPINE1), MHC class I antigen (HLA-C), Transcription factor BTF3 (BTF3L4), Cytochrome c oxidase subunit 6C (COX6C), Myb-binding protein 1A (MYBBP1A), and MYCBP protein (MYCBP) (Figure S8a, Table S4).
Of the proteins identified exclusively when cells were treated with 0.63 μg/mL of venom there were 28 proteins from MCF7 cell line, from which we highlight Cytochrome b-c1 complex subunit 1, mitochondrial (UQCRC1), Latexin (LXN), Vesicle-associated membrane protein-associated protein B/C (VAPB), Cysteine and glycine-rich protein 1 (CSRP1), and DDB1- and CUL4-associated factor 7 (DCAF7). From MDA-MB-231 cell line, 20 proteins were identified from which we highlight Eukaryotic translation initiation factor 2 subunit 2 (EIF2S2), low molecular weight phosphotyrosine protein phosphatase (ACP1), Transcription factor A, mitochondrial (TFAM), and Rae1 protein homolog (RAE1) (Table 1 and Table S4, Figure S8b).
Proteins identified exclusively at the high dose treatment of 2.5 μg/mL of venom included 17 proteins detected in the MCF7 cell line, from which we highlight 26S proteasome non-ATPase regulatory subunit 6 (PSMD6) and 60S ribosomal protein L36 (RPL36). Similarly, we identified 44 proteins in MDA-MB-231 cell line from which we highlight protein Ubiquitin carboxyl-terminal hydrolase (HEL-117) (Figure S8c, Table S4).
We also identified proteins that were expressed exclusively when cells were treated with low and high doses of B. jararaca venom. From the MCF7 cell line, we identified 40 proteins from which we highlight: Serine/threonine-protein kinase ATR (ATR) and Jupiter microtubule-associated homolog 1 (JPT1). From the MDA-MB-231 cell line, we identified 22 proteins from which we highlight Angiomotin (AMOT) and Small nuclear ribonucleoprotein G (SNRPG) (Table 1, Tables S3 and S4, Figure S8d).
Among the proteins that were identified in the control group and at the 2.5 μg/mL venom treatment, we highlight: Small nuclear ribonucleoprotein Sm D1 (SNRPD1) in MCF7, and Epididymis secretory protein Li 71 (Hel-S-71) and Isoform Far upstream element-binding protein 3 (FUBP3) in MDA-MB-231 cell line (Figure S8e, Table S4). Of proteins that were identified in the control group and low 0.63 μg/mL dose of venom, we highlight Serine/arginine-rich splicing factor 10 (SRSF10) in MCF7 cell line, and Histone H2B type 2-E (HIST2H2BE), 60S ribosomal protein L34 and L37a (RPL34, RPL37A) in MDA-MB-231 cell line (Figure S8f, Tables S3 and S4).

3. Discussion

The present study successfully identified cancer-related proteins that undergo significant changes upon B. jararaca venom treatment of MCF7 and MDA-MB-231 cells, including SNX3, HEL-S-156an, UQCRC1, RPL36, and ATR identified in MCF7 cell line and H3C15/HIST2H3, HEL-S-132P, MTCH2, TFAM, KCTD12, RPL34, and RPL37A identified in MDA-MB-231 cell line, and histone H3F3B and LAP3 in both cell lines.
With the discovery of rattlesnake venom’s antitumor activity in 1931 by Essex and Priestley [115] and later the angiotensin-converting enzyme inhibitor drug Captopril in 1981, developed originally from B. jararaca snake venom [116,117,118], several groups have been focusing on to the therapeutic potentials of bioactive compounds in snake venoms by isolating and characterizing components of the venom and analyzing their pharmacological potential that may lead to the development of more efficacious drugs [119,120].
Snake venoms have a complex mixture of proteins that can account for 95% of the total dried weight. Metalloproteases, serine proteases, LAAOs, and PLA2s are some of the most abundant proteins present in bothropic venoms, and the toxic action of these proteins are responsible for myotoxic effects, disruption in the coagulation cascade through thrombin-like proteins (serine proteases), imbalance in the blood homeostasis by the action of some metalloproteinases, apoptosis induction, and changes in the protein expression in oxidative stress and energy metabolism by the increase of Ca2+ influx. Several studies have revealed that these enzymes are critical regulators of cancer pathologies in several types of cancers including breast cancer [121,122,123,124,125,126,127]. Moreover, it has been previously shown that B. jararaca venom may have an important antitumor effect on Ehrilich ascites tumor cells in vivo and in vitro [36].
Aiming at the drug development, several groups have used proteomics approach to observe up- or down-regulated effects of proteins resulting from disease activity or side effects of the treatments [44,45,128,129,130,131]. Moreover, the quest for the identification of biological markers or biomarkers using venom that may help in the early detection of pathologies such as cancer, and the quest for products able to evaluate the metastatic potential and propose new forms of treatment to tumors, remains an important goal of research groups and pharmaceutical companies.
Proteomics studies using tumor cell lines has been performed to systematically characterize protein complexes that may be important in disease or drug development, comparing the protein levels of cells in normal versus pathological situations [132]. Studies have shown that venoms from several snakes such as Calloselasma rhodostoma, Macrovipera lebetina and Bothrops mattogrossensis significantly inhibit cell viability, either with crude or fractionated venoms and in different tumor cell lines such as LS174T (colorectal adenocarcinoma), HCT116 (colorectal carcinoma), HT29 (colorectal adenocarcinoma), HEL92.1.7 (erythroleukemia), and SK-BR-3 (breast adenocarcinoma) cell lines [133,134].
To evaluate the cytotoxicity of B. jararaca venom, we tested two different invasive breast cancer cell lines with different phenotypic and genotypic differences, MCF7 and MDA-MB-231 cell lines. They were treated with different venom concentrations ranging from 0.32 µg/mL to 20 µg/mL. When treated with over 10 µg/mL of venom both MCF7 and MDA-MB-231 cell lines showed rounded and irregular morphologies and the MCF7 cell line detached from the plate. These results are similar to the results observed by Bernardes-Oliveira et al. (2016) when testing 50 µg/mL of venom from B. jararaca and Bothrops erythromelas for 48 h on SiHa HPV-16 and HeLa cells, both derived from cervical tumors. They observed that the cells became more rounded, with gradual size reduction and detachment from the cell monolayer [135]. In addition to the morphological analyses, we evaluated cell viability through mitochondrial cell activity using WST-1 assay for cell proliferation and viability. We observed that cell death of both MCF7 and MDA-MB-231 cell lines started with 2.5 µg/mL of B. jararaca venom treatment. Based on both morphological and cell viability analyses, we selected the concentrations of 0.63 µg/mL (low sub-toxic dose) and 2.5 µg/mL (high sub-toxic dose) for further proteomic analyses.
Semi-quantitative proteomics analysis of the MDA-MB-231 cell line treated with the low 0.63 μg/mL and the sub-lethal dose of 2.5 μg/mL of B. jararaca venom presented less proteins with FC ≥ 1.5 when compared to the MCF7 cell line. We identified 137 proteins with FC ≥ 1.5 at 2.5 μg/mL venom treatment in MCF7 while 34 proteins in MDA-MB-231 at 2.5 μg/mL venom treatment. At the low 0.63 μg/mL venom treatment, 25 proteins in MCF7 presented FC ≥ 1.5 and 16 proteins in MDA-MB-231. However, MDA-MB-231 showed more proteins with FC ≤ 0.67 compared to MCF7 cell lines. MDA-MB-231 showed 28 proteins presenting FC ≤ 0.67 at 0.63 μg/mL venom treatment and 19 proteins in MCF7; and 41 proteins with FC ≤ 0.5 at 2.5 μg/mL venom treatment in MDA-MB-231 cell lines compared to 23 proteins in MCF7 cell line. Among these proteins, our data showed that B. jararaca venom was able to modulate, up or down, the abundance of several proteins playing essential roles in tumorigenicity. For example, the venom treatment increased the abundance of SNX3, HEL-S-156, MCCC2, and GSTM3 in the MCF7 cell line, and H3C15/HIST2H3 and MTCH2 in the MDA-MB-231 cell line whose overexpression has been related to increase in tumorigenicity (Table 1). Furthermore, among the proteins that decreased the abundance presenting FC ≤ 0.67, we highlight PCBD, PSMD5, RPS29, H3F3B, VDAC1 and VDAC2, ATP5PD, and HIST1H4J in MCF7 cell line and IGF2BP1, THBS1, and CYR61 in MDA-MB-231 cell line. These proteins have also been described to be related to several mechanisms and pathways in the onset of different types of cancer (Table 1). Additionally, cytotoxicity analysis performed in our laboratory (data not shown) evaluated the effects of B. jararaca venom on non-tumor HEK293 and HUVEC cell lines, which showed relatively less sensitivity to cell-death induction by snake venoms compared to MCF7 and MDA-MB-231 cancer cells. These data suggest that cancer cells could be more sensitive to snake venom components than non-cancer cells in certain situations.
Proteomic analysis also revealed exclusive proteins, i.e., proteins that were only identified in one or two of the three conditions analyzed. Among these proteins, we highlight the UQCRC1, identified in the MCF7 cell line only when these cells were treated with 0.63 µg/mL venom. High expression of this protein was observed in 74% of cases of breast cancer and 34% of ovarian cancer [104] (Table 1). Among proteins identified exclusively in the MCF7 cell line we highlight ATR (serine/threonine-protein kinase ATR) which plays important roles for cell survival and is considered a major mediator of DNA response in human cells, preventing cells with damaged or incompletely replicated DNA from entering mitosis when cells are damaged by radiotherapy or chemotherapy during cancer treatment [53]. In the MDA-MB-231 cell line, we highlight AMOT (angiomotin) protein identified only when cells were treated with either low or high concentration venom treatment and is known to play a critical role in angiogenesis, proliferation, and migration of cancer cells and is also known to promote proliferation and invasion of several types of tumor cells, including breast, prostate, colon, cervical, and liver cancer (Table 1).
Semi-quantitative proteomics analysis of proteins with FC ≥ 1.5 after addition of venom in both MCF7 and MDA-MB-231 cell lines also identified enrichment of proteins related to protein metabolism pathways as well as cellular stress response pathway. The protein metabolism pathway ranges from protein synthesis to post-translational protein modification and degradation [136]. It is likely that the set of modifications undergone in the protein metabolism pathway undermine the tumorigenic activity of the cell line MCF7 and MDA-MB-231. Therefore, the activation of the apoptotic pathways using snake venom may be a potential treatment approach aiming at tumor degradation preferentially over the normal cells. The ribosome pathway enrichment was also observed in both cell lines. The lack of ribosomes in the cell impairs cell growth even under optimal conditions for it to happen [137]. Protein enrichment in ribosome pathway suggests increased protein synthesis, which may be related to damage repair mechanisms and cellular responses to external stimuli. Alterations in ribosomal pathways are related to cellular alterations and susceptibility to cancer. However, there is still no consensus if deregulation of ribosomes alteration is a cancer consequence or cause [138]. In addition, downregulation of cell maintenance-related proteins may also be related to tumorigenicity, thus opening the door for debate on the role of ribosomes in tumor cell lines [139].
The increase of terms involved in cellular stress response pathway, the metabolic pathway, the proteasome pathway, and the spliceosome pathway in MCF7 cell line suggests that the cells activated the metabolic pathways in response to stress, mRNA transcription and spliceosome activation, and ultimately activation of the proteasome. The proteasome is a cell apparatus that has several cellular functions, such as cell cycle regulation, differentiation, signal transduction pathway, antigen processing for immune responses, stress signaling, inflammatory responses and apoptosis, the latter being mediated by ubiquitination [140,141,142]. Since the inhibition of proteasome is related to cancer onset (i.e., p53 stability), proteolytic activity of the proteasome and its role in both biology and cancer treatment suggest an effective way to treat cancer [143,144].
In the MCF7 cell line treated with B. jararaca snake venom, a higher number of proteins changed in abundance when compared to the MDA-MB-231 cell line. Some of our observations showed that B. jararaca venom was able to up or down modulate the abundance of proteins that plays essential roles in tumorigenicity, conversely to what has been frequently described in the literature to decrease or increase cell tumorigenicity [113,145]. For example, the venom treatment was able to increase abundance of SNX3, HEL-S-156, MCCC2, and GSTM3 in the MCF7 cell line and H3C15/HIST2H3 and MTCH2 in the MDA-MB-231 cell line whose overexpression is related to increase in tumorigenicity. Although we cannot explain the mechanism on how the increase in abundance of these already overexpressed proteins is contributing to cell death of MCF7 and MDA-MB-231 cells, we have clues based on a recent published manuscript by Dias and colleagues (2019). The authors demonstrated that overloading the cell stress pathway in Y1 adrenocortical mouse tumor cell line under the effect of FGF-2 growth factor disrupted the cell homeostasis and sensitized the Y1 tumor cells to stress-oriented therapeutic inhibitors. Dias concludes that further stimulation of the same signaling pathways may further increase mobilization and dependence on stress response pathways in tumor cells, improving the efficacy and selectivity of therapeutic interventions [146]. This data suggests that B. jararaca snake venom or some of its content may be used to down-modulate already known overexpressed proteins in cancer or overexpress oncogenic proteins that lead cells to stress, by disrupting the tumor homeostasis leading cells to death. This is particularly interesting for cancer types that are more resistant to the chemotherapeutic agents, such as the MCF7 cells. In this case, B. jararaca snake venom would work as an adjuvant treatment of cancer, making cells more sensitive to the cytotoxic agents.
Moreover, due to the complexity of B. jararaca venom and the complex response of different cell lines, this snake venom is also a promising candidate in the prognostic aid of different tumors, assisting in the assessment of tumor level, as highlighted for the several overexpressed proteins. For example, targeted-knockdown of some of these proteins may result in tumor growth inhibition. In addition, B. jararaca venom has induced overexpression of proteins related to pathways that impairs antitumor activity such as the protein metabolism pathways, proteasome, and cellular stress response pathways.
It is possible that many responses observed here were due to specific or multiple venom components. Nevertheless, the current approach of applying nLC-MS/MS proteomics to cancer cells treated with crude snake venom is a promising strategy for the identification of proteins with potential application in cancer cell therapy. Overall, this study was able to identify several cancer-related proteins that undergo significant changes upon venom treatment of MCF7 and MDA-MB-231 cells. Future studies should address specific mechanisms by which the snake venom and some of its content may contribute to MCF7 and MDA-MB-231 cell death and survival.

4. Conclusions

Quantitative proteomic analysis of breast cancer cell lines allowed us to identify several proteins whose abundance (FC) increased more than 1.5×, and proteins that abundance decreased less than 0.67×, after 24 h treatment with B. jararaca at either low sub-toxic dose of 0.63 μg/mL or high sub-toxic dose of 2.5 μg/mL. Most of these proteins identified suggest that the treatment with the venom may activate mitochondrial apoptotic pathways leading cells to death. In addition, several of the identified proteins play important roles related to cell proliferation, invasion, metastasis, apoptosis, and stress response. Therefore, these data show that B. jararaca venom or some of its toxin or components can inhibit tumor cell proliferation and survival and can potentially be used to identify novel targets for cancer therapy.

5. Materials and Methods

5.1. Cell Culture and Maintenance

MCF7 and MDA-MB-231 breast cancer cell lines were previously acquired from ATCC (Manassas, VA, USA) and maintained at our laboratory cell bank. Cells were thawed and cultured at 37 °C in 5% CO2 in RPMI 1640 culture medium (Gibco, Life Technologies, Grand Island, NY, USA) supplemented with 10% inactivated fetal bovine serum (Cultilab, Campinas, Brazil) and 25 mg/mL ampicillin and 100 mg/mL streptomycin antibiotics. Cells were split when the confluence reached 80%.

5.2. Bothrops jararaca Venom

The venom of B. jararaca (lot 01/09-2) used in this study was a pool made from the extraction of 697 snakes collected at various Brazilian locations (in the states of São Paulo, Paraná, and Santa Catarina). Extracted venom were pooled and lyophilized at the department of herpetology at Butantan Institute under the coordination of Dr. Marisa Maria Teixeira da Rocha, and assigned for use and testing at the Laboratory of Applied Toxinology at Butantan Institute. The methods and use of venom in this work were approved by Butantan Institute Ethics Committee under the certification CEUAIB #9766150719 (2019).

5.3. Cell Viability Assay

The cytotoxicity assays were carried out by treating the cells with different concentrations of venom in the range from 0.1 μM to 20 μM. The cytotoxicity tests were performed using the colorimetric method of cell viability analysis using the WST-1 cell proliferation reagent kit (Roche, Mannheim, Germany) according to manufacturer’s instructions. Briefly, 104 cells were plated in each well of a flat bottom 96-well microtiter plate. Cells were cultured at 37 °C in 5% CO2 and on the day prior to analysis (70–80% confluency), culture medium from each well was exchanged to 100 μL of medium in the absence or presence of venom in different concentrations. After 24 h, 10 μL of the WST-1 reagent previously dissolved into an Electro Coupling Solution (ECS) were added to each well and incubated for 4 h at 37 °C in 5% CO2. The absorbance at 450 nm of each well was measured using a microplate reader FlexStation 3 spectrophotometer (Molecular Devices, San Jose, CA, USA). The determination of the absorbance and the quantification of viable cells were calculated using SoftMax Pro (version 5.1, Molecular Devices) and Microsoft Excel software.

5.4. Cell Treatment with B. jararaca Venom

In a six-well plate, 2 × 105 cells were plated and after the cells had reached a confluence of 70–80% (2–3 days), different concentrations of B. jararaca venom were added to the culture medium, starting from the cell cytotoxicity threshold up to a 1000-fold dilution for 24 h. For mass spectrometry analysis, cells were washed with ice-cold PBS and lysed with 1 mL ice-cold 8 M urea supplemented with cOmpleteTM protease and phosphatase inhibitors cocktail (Merck Millipore, Burlington, MA, USA). Fifty microliters of the cell lysate were separated for protein quantification using BCA (Bicinchoninic acid assay, Thermo Pierce, Walthan, MA, USA) and the lysate were stocked at −80 °C until further preparation of the sample for mass spectrometry analysis. Both cytotoxicity and proteomics experiments were performed in three biological replicates.

5.5. Sample Preparation for Proteomic Analysis

To the cell lysates, four microliters of 10 mM dithiothreitol (DTT) were added to reduce the disulfide bonds and incubated at 56 °C for 1 h. After incubation with DTT, 40 μL of iodoacetamide (IAA) was added for the alkylation of the cysteines for 1 h in the dark at room temperature. After incubation, 2.5 mL of ammonium acetate and 20 μL of trypsin (Sigma Aldrich, St. Louis, MO, USA) were added enough to have an enzyme:protein ratio of 1:50 and incubated at 37 °C overnight. The reaction was stopped by adding 30 μL of 100% acetic acid. Samples volume were reduced in a SpeedVac (Hetovac VR-1, Heto Lab Equipment, Allerød, Denmark) until the volume has lowered to 50 μL and were further desalted on in-house manufactured stop-and-go extraction tips (Stage-Tip) with three SDB-XC (styrene-divinylbenzene, Empore, 3M, Royersford, PA, USA) membranes as previously described [147]. The stage-tips were initially conditioned with 100% methanol followed by 0.1% formic acid. Samples were applied to the stage-tips and the peptides bound to the membranes were eluted with 50% acetonitrile, 0.1% formic acid. The eluate was lyophilized for further analysis in the mass spectrometer.

5.6. Mass Spectrometry Analysis

Dried samples were resuspended in 0.1% formic acid and analyzed on LTQ-Orbitrap Velos mass spectrometer (Thermo Scientific, Bremen, Germany) coupled to an EASY II nano liquid chromatographer (Thermo Scientific) using the shotgun approach. The spectrometer was equipped with a nanospray source connected to an in-house prepared analytical column (10 cm × ID 75 μm × OD 360 μm) packaged with 7 cm of 5 μm C18 resin (Jupiter, Phenomenex, Torrance, CA, USA). Precolumn (7 cm × ID 75 μm × OD 360 μm) was also prepared in-house packed with 5 cm of 10 μm C18 resin (Acqua, Phenomenex). The LC-MS/MS analyses were carried out by injecting 5 μg of the peptide extract and the peptides eluted from the column with a gradient of 5–40% in 100 min of solvent B (acetonitrile 90%, formic acid 0.1%) at a flow rate of 200 nL/min. The nanospray source was operated at 1.8 kV. The peptide mixture was analyzed by the acquisition of spectra in the full scan mode at a resolution of 30,000 for the determination of molecular masses (MS) of up to 10,000 Da. The 10 most intense peaks were automatically selected via data dependent acquisition (DDA) for the subsequent acquisition of spectra of the ions product to MS/MS for amino acid sequence determination at a resolution of 7500, a maximum injection time of 30 ms, a range of 200 to 2000 m/z, and a dynamic exclusion of 70 s. Protein identification was performed using Mascot (Matrix Science, version 2.4.0, Boston, MA, USA), MaxQuant (version 1.4.1.2, www.maxquant.net accessed on 1 October 2019), and Peaks Studio (version 10, Bioinformatic Solutions Inc, Toronto, Canada) against the Homo sapiens database downloaded from the Uniprot in March 2019 (www.uniprot.org accessed on 1 October 2019). As search parameters, oxidation of methionine was set as a variable modification and carbamidomethylation of cysteine was set as a fixed modification. Searches were performed with mass error tolerance of 10 ppm for MS and 0.3 Da for MS/MS, and trypsin was selected as the proteolytic enzyme used in the digestion of proteins, and up to two miscleavages were allowed. Mass spectrometry-based proteomic analysis of MCF7 and MDA-MB-231 cell lysates, treated or not treated with low or high doses of B. jararaca venom were performed in triplicate. All raw data files for these analyses were uploaded and are available at: http://massive.ucsd.edu/MSV000084138/ accessed on 5 July 2021.

5.7. Proteome Functional and Enrichment Analysis

The entire proteome was analyzed to classify protein profiles by biological processes, molecular function, cellular components, and cellular pathways using Gene Ontology (geneontology.org accessed on 1 October 2019) and PantherDB (pantherdb.org accessed on 1 October 2019). The proteins presenting different profiles in the venom conditions were also analyzed using KEGG (www.genome.jp/kegg accessed on 1 October 2019) canonical pathways to determine the protein group or pathways that have undergone changes with different concentrations of the venom treatment. The PANTHER Classification System (Protein ANalysis THrough Evolutionary Relationships) is a research tool that covers evolutionary and functional information on protein family genes and has been widely used to understand protein evolution and its functional classification [148,149]. We conducted the functional classification for the molecular function, cell component, biological process, and protein family classification profiles [150,151] of the different proteins in each cell line and treatment at low and high B. jararaca venom that presented FC ≥ 1.5 in order to identify an enrichment pattern of protein families and pathways affected by the venom in the cell lines.

5.8. Semi-Quantitative Proteomics Analysis

Semi-quantitative protein analysis was performed using the Label Free Quantification MaxLFQ algorithm from MaxQuant software with an FDR rate of ≤1% to compare relative abundance of proteins in each of the cell lines [152]. Data generated by MaxQuant were entered into the Perseus program to further perform statistical and bioinformatics analyzes [153]. Proteins identified in the contaminant database and the decoy database were removed. As a protein identification criterion, it was considered that only peptides identified with the posterior error probability (PEP) ≤ 0.01 in at least one biological replicate, the minimum identification of eight ions belonging to the b and y ion series in the MS/MS spectra, and the occurrence of at least one unique peptide. We considered the intensity values of the LFQ that are normalized by the Maxquant software based on the sum of the intensity of all peptides of all identified proteins. LFQ data was considered for calculation, when the intensity data were present in at least two out of three replicates. Protein abundance or fold change (FC) analysis was performed using Microsoft Excel software. In addition, proteins that had a zero value in two of the three conditions were analyzed separately. The breast cancer cell lines treated with the whole B. jararaca venom for two different doses were compared according to the differentially expressed proteins with FC ≥ 1.5 by hierarchical clustering using as variables the average log2 fold change for the three replicates for each treatment, normalized with mean-centering. The Clustering analysis was performed using R statistical software version 3.6.3 (http://www.R-project.org accessed on 1 October 2019). The set of protein dissimilarities were computed using the “Euclidean” distance with the function “dist” to the hierarchical clustering based on the package and function “hclust”. There was employed the agglomerative method with “ward.D2”. The fold change and the protein–protein interactions provided by information from the String database [154] was used to explore the biological interactions of proteins identified as differentially abundant between control and cells treated with B. jararaca venom using the Homo sapiens reference genome.

5.9. Principal Component Analysis

Principal component analysis (PCA) was applied to find which linear combinations of the differentially expressed proteins with ≥1.5 would explain most of the variability for the different cell lineage conditions. The PCA analysis was performed using the R-statistics packages FactoMineR (accessed on 1 November 2020) [155] and Factoextra (http://www.sthda.com/english/rpkgs/factoextra accessed on 1 November 2020) for graphical visualization.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/toxins13080519/s1, Figure S1: Optical microscopy analysis of breast cancer cell lines treated with different concentrations of B. jararaca snake venom. MCF7 (a) and MDA-MB-231 (b) cell lines (10× magnification), Figure S2: Functional classification of the proteins presenting fold change (FC) ≥ 1.5 according to Molecular Function GO enrichment analysis. (a) MCF7 cell line. (b) MDA-MB-231 cell line, Figure S3: Functional classification of the proteins presenting FC ≥ 1.5 according to the Biological Process GO enrichment analysis. (a) MCF7 cell line. (b) MDA-MB-231 cell line, Figure S4: Functional classification of the proteins presenting FC ≥ 1.5 according to the Cellular Component GO enrichment analysis. (a) MCF7 cell line. (b) MDA-MB-231 cell line, Figure S5: Functional classification of the proteins presenting FC ≥ 1.5 according to Protein Family Classification GO enrichment analysis. (a) MCF7 cell line. (b) MDA-MB-231 cell line, Figure S6a: Protein-protein interaction of proteins identified in MCF7 cell line presenting FC ≤ 0.67 in cells treated with 2.5 μg/mL of B. jararaca venom compared to control, Figure S6b: Protein-protein interaction of proteins identified in MCF7 cell line presenting FC ≥ 1.5 in cells treated with 0.67 μg/mL of B. jararaca venom compared to control, Figure S6c: Protein-protein interaction of proteins identified in MCF7 cell line presenting FC ≤ 0.67 when in cells treated with 0.63 μg/mL of B. jararaca venom compared to control, Supplemental Figure S7a: Protein-protein interaction of proteins identified in MDA-MB-231 cell line presenting FC ≤ 0.67 in cells treated with 2.5 μg/mL of B. jararaca venom compared to control, Figure S7b: Protein-protein interaction of proteins identified in MDA-MB-231 cell line presenting FC ≤ 0.67 in cells treated with 0.63 μg/mL of B. jararaca venom compared to control, Figure S7c: Protein-protein interaction of proteins identified in MDA-MB-231 cell line presenting FC ≥ 2.5 in cells treated with 0.63 μg/mL of B. jararaca venom compared to control, Figure S8: Intensity plots of proteins identified exclusively in one or two of the conditions tested (Control, 0.63 μg/mL or 2.5 μg/mL of B. jararaca venom) in both MCF7 and MDA-MB-231 cell lines (Values presented in Table S4). (a) Proteins identified exclusively in the control group. (b) Proteins identified exclusively in the 0.63 μg/mL of venom treatment. (c) Proteins identified exclusively in the 2.5 μg/mL of venom treatment. (d) Proteins identified exclusively in the 0.63 and 2.5 μg/mL of venom treatments. (e) Proteins identified exclusively in the control group and 2.5 μg/mL of venom treatments. (f) Proteins identified exclusively in the control group and 0.63 μg/mL of venom treatment, Table S1a: Mass spectrometry-based proteomics analysis of MCF7 cell line treated with PBS (CTRL), 0.63 μg/mL and 2.5 μg/mL of B. jararaca snake venom, Table S2a: Mass spectrometry-based proteomics analysis of MCF7 cell lines treated with PBS (CTRL), 0.63 mg/mL and 2.5 μg/mL of B. jararaca snake venom. Fold change based on the control LFQ intensity, Table S1b: Mass spectrometry-based proteomics analysis of MDA-MB-231 cell line treated with PBS (CTRL), 0.63 μg/mL, and 2.5 μg/mL of B. jararaca snake venom, Table S2b: Mass spectrometry-based proteomics analysis of MDA-MB-231 cell lines treated with PBS (CTRL), 0.63 μg/mL and 2.5 μg/mL of B. jararaca snake venom. Fold change based on the control LFQ intensity, Table S3: Protein list of each of the seven clusters corresponding to Figure 3, Table S4a: MCF7 proteins identified exclusively at one or two of the conditions tested: PBS (CTRL), 0.63 μg/mL and 2.5 μg/mL of B. jararaca snake venom, Table S4b: MDA-MB-231 proteins identified exclusively at one or two of the conditions tested: PBS (CTRL), 0.63 μg/mL and 2.5 μg/mL of B. jararaca snake venom.

Author Contributions

Conceptualization, L.K.I.; methodology, L.K.I.; validation, C.Y.K. and S.S.S.A.; formal analysis, C.Y.K., S.S.S.A., F.M., W.d.S.S., H.M.C., D.C.-C., E.S.F., M.Y.N.-J. and L.K.I., investigation, C.Y.K., S.S.S.A., I.F.L., and D.C.-C.; resources, L.K.I.; data curation, M.Y.N.-J.; writing—original draft preparation, C.Y.K., S.S.S.A., E.S.F., M.Y.N.-J. and L.K.I.; writing—review and editing, C.Y.K., E.S.F., M.Y.N.-J. and L.K.I.; supervision, E.S.F., M.Y.N.-J. and L.K.I.; project administration, L.K.I.; funding acquisition, E.S.F. and L.K.I. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by grants 2013/07467-1, 2016/04000-3, and 2017/17943-6 from the São Paulo Research Foundation (FAPESP). C.Y.K. was supported by FAPESP master’s degree program fellowship 2017/06496-9. E.S.F. was supported by Erna and Jakob Michael “Visiting Professorship”, Department of Biological Regulation, Weizmann Institute of Science, Rehovot, Israel and grant 302809/2016-3 from the Brazilian National Council for Scientific and Technological Development, CNPq, Brazil. F.M., W.d.S.S. and S.S.S.A. were supported by the Coordination for the Improvement of Higher Education Personnel (CAPES, Brazil) and CNPq institutional fellowship #88882.442313/2019-01 (Butantan Institute, Brazil) for W.d.S.S. and # 131408/2019-4 for S.S.S.A. (University of São Paulo, Brazil).

Institutional Review Board Statement

The methods and use of venom in this work were approved by Butantan Institute Ethics Committee under the certification CEUAIB #9766150719 (2019).

Informed Consent Statement

Not Applicable.

Data Availability Statement

All generated raw files for these analyses were uploaded and are available at: http://massive.ucsd.edu/MSV000084138/ (acccessed on 5 July 2021).

Acknowledgments

We thank Solange M.T. Serrano, Hugo A. Armelin and Zhibin Chen for useful discussion.

Conflicts of Interest

The authors declare that there is no conflict of interest.

References

  1. Brazilian Ministry of Health. Accidents at Work by Venomous Animals among Workers at the Countryside, Forest, and Water in Brazil from 2007 to 2017 (Port). 2019; Volume 50. Available online: https://portalarquivos2.saude.gov.br/images/pdf/2019/marco/29/2018-059.pdf (accessed on 1 March 2021).
  2. Ribeiro, L.A.; Jorge, M.T. Bites by snakes in the genus Bothrops: A series of 3139 cases. Rev. Soc. Bras. Med. Trop. 1997, 30, 475–480. [Google Scholar] [CrossRef] [Green Version]
  3. Tanjoni, I.; Weinlich, R.; Della-Casa, M.S.; Clissa, P.B.; Saldanha-Gama, R.F.; de Freitas, M.S.; Barja-Fidalgo, C.; Amarante-Mendes, G.P.; Moura-da-Silva, A.M. Jararhagin, a snake venom metalloproteinase, induces a specialized form of apoptosis (anoikis) selective to endothelial cells. Apoptosis 2005, 10, 851–861. [Google Scholar] [CrossRef] [PubMed]
  4. Calvete, J.J.; Juarez, P.; Sanz, L. Snake venomics. Strategy and applications. J. Mass Spectrom. 2007, 42, 1405–1414. [Google Scholar] [CrossRef] [PubMed]
  5. Fox, J.W.; Serrano, S.M. Exploring snake venom proteomes: Multifaceted analyses for complex toxin mixtures. Proteomics 2008, 8, 909–920. [Google Scholar] [CrossRef] [PubMed]
  6. Doley, R.; Kini, R.M. Protein complexes in snake venom. Cell. Mol. Life Sci. 2009, 66, 2851–2871. [Google Scholar] [CrossRef]
  7. Serrano, S.M.; Shannon, J.D.; Wang, D.; Camargo, A.C.; Fox, J.W. A multifaceted analysis of viperid snake venoms by two-dimensional gel electrophoresis: An approach to understanding venom proteomics. Proteomics 2005, 5, 501–510. [Google Scholar] [CrossRef]
  8. Jin, H.; Varner, J. Integrins: Roles in cancer development and as treatment targets. Br. J. Cancer 2004, 90, 561–565. [Google Scholar] [CrossRef]
  9. White, J. Snake venoms and coagulopathy. Toxicon 2005, 45, 951–967. [Google Scholar] [CrossRef]
  10. Gay, C.; Sanz, L.; Calvete, J.J.; Pla, D. Snake Venomics and Antivenomics of Bothrops diporus, a Medically Important Pitviper in Northeastern Argentina. Toxins 2015, 8, 9. [Google Scholar] [CrossRef] [Green Version]
  11. Sanz, L.; Perez, A.; Quesada-Bernat, S.; Diniz-Sousa, R.; Calderon, L.A.; Soares, A.M.; Calvete, J.J.; Caldeira, C.A.S. Venomics and antivenomics of the poorly studied Brazil’s lancehead, Bothrops brazili (Hoge, 1954), from the Brazilian State of Para. J. Venom. Anim. Toxins Incl. Trop. Dis. 2020, 26, e20190103. [Google Scholar] [CrossRef]
  12. Tashima, A.K.; Zelanis, A.; Kitano, E.S.; Ianzer, D.; Melo, R.L.; Rioli, V.; Sant’anna, S.S.; Schenberg, A.C.; Camargo, A.C.; Serrano, S.M. Peptidomics of three Bothrops snake venoms: Insights into the molecular diversification of proteomes and peptidomes. Mol. Cell. Proteom. 2012, 11, 1245–1262. [Google Scholar] [CrossRef] [Green Version]
  13. Gren, E.C.K.; Kitano, E.S.; Andrade-Silva, D.; Iwai, L.K.; Reis, M.S.; Menezes, M.C.; Serrano, S.M.T. Comparative analysis of the high molecular mass subproteomes of eight Bothrops snake venoms. Comp. Biochem. Physiol. Part D Genom. Proteom. 2019, 30, 113–121. [Google Scholar] [CrossRef]
  14. Zelanis, A.; Tashima, A.K.; Pinto, A.F.; Paes Leme, A.F.; Stuginski, D.R.; Furtado, M.F.; Sherman, N.E.; Ho, P.L.; Fox, J.W.; Serrano, S.M. Bothrops jararaca venom proteome rearrangement upon neonate to adult transition. Proteomics 2011, 11, 4218–4228. [Google Scholar] [CrossRef]
  15. Farias, I.B.; Morais-Zani, K.; Serino-Silva, C.; Sant’Anna, S.S.; Rocha, M.; Grego, K.F.; Andrade-Silva, D.; Serrano, S.M.T.; Tanaka-Azevedo, A.M. Functional and proteomic comparison of Bothrops jararaca venom from captive specimens and the Brazilian Bothropic Reference Venom. J. Proteom. 2018, 174, 36–46. [Google Scholar] [CrossRef]
  16. Augusto-de-Oliveira, C.; Stuginski, D.R.; Kitano, E.S.; Andrade-Silva, D.; Liberato, T.; Fukushima, I.; Serrano, S.M.; Zelanis, A. Dynamic Rearrangement in Snake Venom Gland Proteome: Insights into Bothrops jararaca Intraspecific Venom Variation. J. Proteome Res. 2016, 15, 3752–3762. [Google Scholar] [CrossRef]
  17. Laing, G.D.; Clissa, P.B.; Theakston, R.D.; Moura-da-Silva, A.M.; Taylor, M.J. Inflammatory pathogenesis of snake venom metalloproteinase-induced skin necrosis. Eur. J. Immunol. 2003, 33, 3458–3463. [Google Scholar] [CrossRef] [Green Version]
  18. Azevedo-Marques, M.M.; Cupo, P.; Hering, S.E. Acidentes Por Animais PeÇonhentos: Serpentes PeÇonhentas. Medicina 2003, 36. [Google Scholar] [CrossRef] [Green Version]
  19. Bjarnason, J.B.; Fox, J.W. Hemorrhagic metalloproteinases from snake venoms. Pharmacol. Ther. 1994, 62, 325–372. [Google Scholar] [CrossRef]
  20. Cidade, D.A.; Simao, T.A.; Davila, A.M.; Wagner, G.; Junqueira-de-Azevedo, I.L.; Ho, P.L.; Bon, C.; Zingali, R.B.; Albano, R.M. Bothrops jararaca venom gland transcriptome: Analysis of the gene expression pattern. Toxicon 2006, 48, 437–461. [Google Scholar] [CrossRef]
  21. Pereira, L.M.; Messias, E.A.; Sorroche, B.P.; Oliveira, A.D.N.; Arantes, L.; de Carvalho, A.C.; Tanaka-Azevedo, A.M.; Grego, K.F.; Carvalho, A.L.; Melendez, M.E. In-depth transcriptome reveals the potential biotechnological application of Bothrops jararaca venom gland. J. Venom. Anim. Toxins Incl. Trop. Dis. 2020, 26, e20190058. [Google Scholar] [CrossRef]
  22. Serrano, S.M.; Maroun, R.C. Snake venom serine proteinases: Sequence homology vs. substrate specificity, a paradox to be solved. Toxicon 2005, 45, 1115–1132. [Google Scholar] [CrossRef]
  23. Fox, J.W.; Serrano, S.M. Insights into and speculations about snake venom metalloproteinase (SVMP) synthesis, folding and disulfide bond formation and their contribution to venom complexity. FEBS J. 2008, 275, 3016–3030. [Google Scholar] [CrossRef]
  24. Matsui, T.; Fujimura, Y.; Titani, K. Snake venom proteases affecting hemostasis and thrombosis. Biochim. Biophys. Acta (BBA)-Protein Struct. Mol. Enzym. 2000, 1477, 146–156. [Google Scholar] [CrossRef]
  25. Da Silva Cunha, K.C.; Fuly, A.L.; de Araujo, E.G. A phospholipase A(2) isolated from Lachesis muta snake venom increases the survival of retinal ganglion cells in vitro. Toxicon 2011, 57, 580–585. [Google Scholar] [CrossRef]
  26. Jain, D.; Kumar, S. Snake venom: A potent anticancer agent. Asian Pac. J. Cancer Prev. 2012, 13, 4855–4860. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Tan, N.; Fung, S.Y. Snake Venom L-Amino Acid Oxidases potential biomedical applications. Mal. J. Biochem. Mol. Biol. 2008, 16, 1–10. [Google Scholar]
  28. Tan, K.K.; Bay, B.H.; Gopalakrishnakone, P. L-amino acid oxidase from snake venom and its anticancer potential. Toxicon 2018, 144, 7–13. [Google Scholar] [CrossRef]
  29. Pawelek, P.D.; Cheah, J.; Coulombe, R.; Macheroux, P.; Ghisla, S.; Vrielink, A. The structure of L-amino acid oxidase reveals the substrate trajectory into an enantiomerically conserved active site. EMBO J. 2000, 19, 4204–4215. [Google Scholar] [CrossRef] [Green Version]
  30. Phillips, D.R.; Scarborough, R.M. Clinical Pharmacology of Eptifibatide. Am. J. Cardiol. 1997, 80, 11B–20B. [Google Scholar] [CrossRef]
  31. Ferreira, S.H. A Bradykinin-Potentiating Factor (Bpf) Present in the Venom of Bothrops Jararca. Br. J. Pharmacol. Chemother. 1965, 24, 163–169. [Google Scholar] [CrossRef] [Green Version]
  32. Camargo, A.C.; Ianzer, D.; Guerreiro, J.R.; Serrano, S.M. Bradykinin-potentiating peptides: Beyond captopril. Toxicon 2012, 59, 516–523. [Google Scholar] [CrossRef]
  33. Fox, J.W.; Serrano, S.M. Approaching the golden age of natural product pharmaceuticals from venom libraries: An overview of toxins and toxin-derivatives currently involved in therapeutic or diagnostic applications. Curr. Pharm. Des. 2007, 13, 2927–2934. [Google Scholar] [CrossRef] [Green Version]
  34. Marsh, N.; Williams, V. Practical applications of snake venom toxins in haemostasis. Toxicon 2005, 45, 1171–1181. [Google Scholar] [CrossRef]
  35. Vyas, V.K.; Brahmbhatt, K.; Bhatt, H.; Parmar, U. Therapeutic potential of snake venom in cancer therapy: Current perspectives. Asian Pac. J. Trop. Biomed. 2013, 3, 156–162. [Google Scholar] [CrossRef] [Green Version]
  36. da Silva, R.J.; da Silva, M.G.; Vilela, L.C.; Fecchio, D. Antitumor effect of Bothrops jararaca venom. Mediat. Inflamm. 2002, 11, 99–104. [Google Scholar] [CrossRef]
  37. Yamazaki, Y.; Morita, T. Snake venom components affecting blood coagulation and the vascular system: Structural similarities and marked diversity. Curr. Pharm. Des. 2007, 13, 2872–2886. [Google Scholar] [CrossRef] [PubMed]
  38. Dias, M.H.; Kitano, E.S.; Zelanis, A.; Iwai, L.K. Proteomics and drug discovery in cancer. Drug Discov. Today 2016, 21, 264–277. [Google Scholar] [CrossRef] [PubMed]
  39. Hanahan, D.; Weinberg, R.A. Hallmarks of cancer: The next generation. Cell 2011, 144, 646–674. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Mora, R.; Valverde, B.; Diaz, C.; Lomonte, B.; Gutierrez, J.M. A Lys49 phospholipase A(2) homologue from Bothrops asper snake venom induces proliferation, apoptosis and necrosis in a lymphoblastoid cell line. Toxicon 2005, 45, 651–660. [Google Scholar] [CrossRef] [PubMed]
  41. Bateman, E.; Venning, M.; Mirtschin, P.; Woods, A. The effects of selected Australian snake venoms on tumour-associated microvascular endothelial cells (TAMECs) in vitro. J. Venom. Res. 2013, 4, 21–30. [Google Scholar]
  42. Moreira, V.; Lomonte, B.; Vinolo, M.A.; Curi, R.; Gutierrez, J.M.; Teixeira, C. An Asp49 phospholipase A2 from snake venom induces cyclooxygenase-2 expression and prostaglandin E2 production via activation of NF-kappaB, p38MAPK, and PKC in macrophages. Mediat. Inflamm. 2014, 2014, 105879. [Google Scholar] [CrossRef] [Green Version]
  43. Leiguez, E.; Giannotti, K.C.; Moreira, V.; Matsubara, M.H.; Gutierrez, J.M.; Lomonte, B.; Rodriguez, J.P.; Balsinde, J.; Teixeira, C. Critical role of TLR2 and MyD88 for functional response of macrophages to a group IIA-secreted phospholipase A2 from snake venom. PLoS ONE 2014, 9, e93741. [Google Scholar] [CrossRef] [Green Version]
  44. Frejno, M.; Meng, C.; Ruprecht, B.; Oellerich, T.; Scheich, S.; Kleigrewe, K.; Drecoll, E.; Samaras, P.; Hogrebe, A.; Helm, D.; et al. Proteome activity landscapes of tumor cell lines determine drug responses. Nat. Commun. 2020, 11, 3639. [Google Scholar] [CrossRef] [PubMed]
  45. Ruprecht, B.; Di Bernardo, J.; Wang, Z.; Mo, X.; Ursu, O.; Christopher, M.; Fernandez, R.B.; Zheng, L.; Dill, B.D.; Wang, H.; et al. A mass spectrometry-based proteome map of drug action in lung cancer cell lines. Nat. Chem. Biol. 2020, 16, 1111–1119. [Google Scholar] [CrossRef] [PubMed]
  46. Furtado, C.M.; Marcondes, M.C.; Sola-Penna, M.; de Souza, M.L.; Zancan, P. Clotrimazole preferentially inhibits human breast cancer cell proliferation, viability and glycolysis. PLoS ONE 2012, 7, e30462. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Hui, L.; Zheng, Y.; Yan, Y.; Bargonetti, J.; Foster, D.A. Mutant p53 in MDA-MB-231 breast cancer cells is stabilized by elevated phospholipase D activity and contributes to survival signals generated by phospholipase D. Oncogene 2006, 25, 7305–7310. [Google Scholar] [CrossRef] [Green Version]
  48. Neve, R.M.; Chin, K.; Fridlyand, J.; Yeh, J.; Baehner, F.L.; Fevr, T.; Clark, L.; Bayani, N.; Coppe, J.P.; Tong, F.; et al. A collection of breast cancer cell lines for the study of functionally distinct cancer subtypes. Cancer Cell 2006, 10, 515–527. [Google Scholar] [CrossRef] [Green Version]
  49. Lv, M.; Shen, Y.; Yang, J.; Li, S.; Wang, B.; Chen, Z.; Li, P.; Liu, P.; Yang, J. Angiomotin Family Members: Oncogenes or Tumor Suppressors? Int. J. Biol. Sci. 2017, 13, 772–781. [Google Scholar] [CrossRef]
  50. Zhang, H.; Fan, Q. MicroRNA-205 inhibits the proliferation and invasion of breast cancer by regulating AMOT expression. Oncol. Rep. 2015, 34, 2163–2170. [Google Scholar] [CrossRef]
  51. Couderc, C.; Boin, A.; Fuhrmann, L.; Vincent-Salomon, A.; Mandati, V.; Kieffer, Y.; Mechta-Grigoriou, F.; Del Maestro, L.; Chavrier, P.; Vallerand, D.; et al. AMOTL1 Promotes Breast Cancer Progression and Is Antagonized by Merlin. Neoplasia 2016, 18, 10–24. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Song, K.H.; Kim, J.H.; Lee, Y.H.; Bae, H.C.; Lee, H.J.; Woo, S.R.; Oh, S.J.; Lee, K.M.; Yee, C.; Kim, B.W.; et al. Mitochondrial reprogramming via ATP5H loss promotes multimodal cancer therapy resistance. J. Clin. Investig. 2018, 128, 4098–4114. [Google Scholar] [CrossRef] [PubMed]
  53. Qiu, Z.; Oleinick, N.L.; Zhang, J. ATR/CHK1 inhibitors and cancer therapy. Radiother. Oncol. 2018, 126, 450–464. [Google Scholar] [CrossRef] [PubMed]
  54. Karnitz, L.M.; Zou, L. Molecular Pathways: Targeting ATR in Cancer Therapy. Clin. Cancer Res. 2015, 21, 4780–4785. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. O’Kelly, J.; Chung, A.; Lemp, N.; Chumakova, K.; Yin, D.; Wang, H.J.; Said, J.; Gui, D.; Miller, C.W.; Karlan, B.Y.; et al. Functional domains of CCN1 (Cyr61) regulate breast cancer progression. Int. J. Oncol. 2008, 33, 59–67. [Google Scholar] [CrossRef] [PubMed]
  56. Huang, Y.T.; Lan, Q.; Lorusso, G.; Duffey, N.; Ruegg, C. The matricellular protein CYR61 promotes breast cancer lung metastasis by facilitating tumor cell extravasation and suppressing anoikis. Oncotarget 2017, 8, 9200–9215. [Google Scholar] [CrossRef] [Green Version]
  57. Hayes, J.D.; Flanagan, J.U.; Jowsey, I.R. Glutathione transferases. Annu. Rev. Pharmacol. Toxicol. 2005, 45, 51–88. [Google Scholar] [CrossRef] [PubMed]
  58. Meding, S.; Balluff, B.; Elsner, M.; Schone, C.; Rauser, S.; Nitsche, U.; Maak, M.; Schafer, A.; Hauck, S.M.; Ueffing, M.; et al. Tissue-based proteomics reveals FXYD3, S100A11 and GSTM3 as novel markers for regional lymph node metastasis in colon cancer. J. Pathol. 2012, 228, 459–470. [Google Scholar] [CrossRef] [PubMed]
  59. Saucedo-Cuevas, L.P.; Ruppen, I.; Ximenez-Embun, P.; Domingo, S.; Gayarre, J.; Munoz, J.; Silva, J.M.; Garcia, M.J.; Benitez, J. CUL4A contributes to the biology of basal-like breast tumors through modulation of cell growth and antitumor immune response. Oncotarget 2014, 5, 2330–2343. [Google Scholar] [CrossRef] [Green Version]
  60. Lan, F.; Shi, Y. Histone H3.3 and cancer: A potential reader connection. Proc. Natl. Acad. Sci. USA 2015, 112, 6814–6819. [Google Scholar] [CrossRef] [Green Version]
  61. Ayoubi, H.A.; Mahjoubi, F.; Mirzaei, R. Investigation of the human H3.3B (H3F3B) gene expression as a novel marker in patients with colorectal cancer. J. Gastrointest. Oncol. 2017, 8, 64–69. [Google Scholar] [CrossRef] [Green Version]
  62. Roberts, E.L.; Newton, R.P.; Axford, A.T. Plasma purine nucleoside phosphorylase in cancer patients. Clin. Chim. Acta 2004, 344, 109–114. [Google Scholar] [CrossRef]
  63. Vareed, S.K.; Bhat, V.B.; Thompson, C.; Vasu, V.T.; Fermin, D.; Choi, H.; Creighton, C.J.; Gayatri, S.; Lan, L.; Putluri, N.; et al. Metabolites of purine nucleoside phosphorylase (NP) in serum have the potential to delineate pancreatic adenocarcinoma. PLoS ONE 2011, 6, e17177. [Google Scholar] [CrossRef]
  64. Dummer, R.; Duvic, M.; Scarisbrick, J.; Olsen, E.A.; Rozati, S.; Eggmann, N.; Goldinger, S.M.; Hutchinson, K.; Geskin, L.; Illidge, T.M.; et al. Final results of a multicenter phase II study of the purine nucleoside phosphorylase (PNP) inhibitor forodesine in patients with advanced cutaneous T-cell lymphomas (CTCL) (Mycosis fungoides and Sezary syndrome). Ann. Oncol. 2014, 25, 1807–1812. [Google Scholar] [CrossRef]
  65. Tan, E.; Besant, P.G.; Zu, X.L.; Turck, C.W.; Bogoyevitch, M.A.; Lim, S.G.; Attwood, P.V.; Yeoh, G.C. Histone H4 histidine kinase displays the expression pattern of a liver oncodevelopmental marker. Carcinogenesis 2004, 25, 2083–2088. [Google Scholar] [CrossRef] [Green Version]
  66. Besant, P.G.; Attwood, P.V. Histone H4 histidine phosphorylation: Kinases, phosphatases, liver regeneration and cancer. Biochem. Soc. Trans. 2012, 40, 290–293. [Google Scholar] [CrossRef] [Green Version]
  67. Long, M.; Sun, X.; Shi, W.; Yanru, A.; Leung, S.T.C.; Ding, D.; Cheema, M.S.; MacPherson, N.; Nelson, C.J.; Ausio, J.; et al. A novel histone H4 variant H4G regulates rDNA transcription in breast cancer. Nucleic. Acids Res. 2019, 47, 8399–8409. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  68. Nacev, B.A.; Feng, L.; Bagert, J.D.; Lemiesz, A.E.; Gao, J.; Soshnev, A.A.; Kundra, R.; Schultz, N.; Muir, T.W.; Allis, C.D. The expanding landscape of ‘oncohistone’ mutations in human cancers. Nature 2019, 567, 473–478. [Google Scholar] [CrossRef] [PubMed]
  69. Bell, J.L.; Wachter, K.; Muhleck, B.; Pazaitis, N.; Kohn, M.; Lederer, M.; Huttelmaier, S. Insulin-like growth factor 2 mRNA-binding proteins (IGF2BPs): Post-transcriptional drivers of cancer progression? Cell. Mol. Life Sci. 2013, 70, 2657–2675. [Google Scholar] [CrossRef] [Green Version]
  70. Fakhraldeen, S.A.; Clark, R.J.; Roopra, A.; Chin, E.N.; Huang, W.; Castorino, J.; Wisinski, K.B.; Kim, T.; Spiegelman, V.S.; Alexander, C.M. Two Isoforms of the RNA Binding Protein, Coding Region Determinant-binding Protein (CRD-BP/IGF2BP1), Are Expressed in Breast Epithelium and Support Clonogenic Growth of Breast Tumor Cells. J. Biol. Chem. 2015, 290, 13386–13400. [Google Scholar] [CrossRef] [Green Version]
  71. Kim, T.; Havighurst, T.; Kim, K.; Albertini, M.; Xu, Y.G.; Spiegelman, V.S. Targeting insulin-like growth factor 2 mRNA-binding protein 1 (IGF2BP1) in metastatic melanoma to increase efficacy of BRAF(V600E) inhibitors. Mol. Carcinog. 2018, 57, 678–683. [Google Scholar] [CrossRef] [PubMed]
  72. Han, W.; Hu, C.; Fan, Z.J.; Shen, G.L. Transcript levels of keratin 1/5/6/14/15/16/17 as potential prognostic indicators in melanoma patients. Sci. Rep. 2021, 11, 1023. [Google Scholar] [CrossRef]
  73. Palko, E.; Poliska, S.; Sziklai, I.; Penyige, A. Analysis of KRT1, KRT10, KRT19, TP53 and MMP9 expression in pediatric and adult cholesteatoma. PLoS ONE 2018, 13, e0200840. [Google Scholar] [CrossRef]
  74. Hitzerd, S.M.; Verbrugge, S.E.; Ossenkoppele, G.; Jansen, G.; Peters, G.J. Positioning of aminopeptidase inhibitors in next generation cancer therapy. Amino Acids 2014, 46, 793–808. [Google Scholar] [CrossRef]
  75. Tian, S.Y.; Chen, S.H.; Shao, B.F.; Cai, H.Y.; Zhou, Y.; Zhou, Y.L.; Xu, A.B. Expression of leucine aminopeptidase 3 (LAP3) correlates with prognosis and malignant development of human hepatocellular carcinoma (HCC). Int. J. Clin. Exp. Pathol. 2014, 7, 3752–3762. [Google Scholar]
  76. He, X.; Huang, Q.; Qiu, X.; Liu, X.; Sun, G.; Guo, J.; Ding, Z.; Yang, L.; Ban, N.; Tao, T.; et al. LAP3 promotes glioma progression by regulating proliferation, migration and invasion of glioma cells. Int. J. Biol. Macromol. 2015, 72, 1081–1089. [Google Scholar] [CrossRef] [PubMed]
  77. Zhang, S.; Yang, X.; Shi, H.; Li, M.; Xue, Q.; Ren, H.; Yao, L.; Chen, X.; Zhang, J.; Wang, H. Overexpression of leucine aminopeptidase 3 contributes to malignant development of human esophageal squamous cell carcinoma. J. Mol. Histol. 2014, 45, 283–292. [Google Scholar] [CrossRef]
  78. Liu, Y.; Yuan, Z.; Song, C. Methylcrotonoyl-CoA carboxylase 2 overexpression predicts an unfavorable prognosis and promotes cell proliferation in breast cancer. Biomark. Med. 2019, 13, 427–436. [Google Scholar] [CrossRef] [PubMed]
  79. Arco, A.D.; Satrustegui, J. New mitochondrial carriers: An overview. Cell. Mol. Life Sci. 2005, 62, 2204–2227. [Google Scholar] [CrossRef] [PubMed]
  80. Hao, J.M.; Chen, J.Z.; Sui, H.M.; Si-Ma, X.Q.; Li, G.Q.; Liu, C.; Li, J.L.; Ding, Y.Q.; Li, J.M. A five-gene signature as a potential predictor of metastasis and survival in colorectal cancer. J. Pathol. 2010, 220, 475–489. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  81. Arigoni, M.; Barutello, G.; Riccardo, F.; Ercole, E.; Cantarella, D.; Orso, F.; Conti, L.; Lanzardo, S.; Taverna, D.; Merighi, I.; et al. miR-135b coordinates progression of ErbB2-driven mammary carcinomas through suppression of MID1 and MTCH2. Am. J. Pathol. 2013, 182, 2058–2070. [Google Scholar] [CrossRef] [Green Version]
  82. De Palma, G.; Sallustio, F.; Curci, C.; Galleggiante, V.; Rutigliano, M.; Serino, G.; Ditonno, P.; Battaglia, M.; Schena, F.P. The Three-Gene Signature in Urinary Extracellular Vesicles from Patients with Clear Cell Renal Cell Carcinoma. J. Cancer 2016, 7, 1960–1967. [Google Scholar] [CrossRef] [Green Version]
  83. Nunes, F.D.; de Almeida, F.C.; Tucci, R.; de Sousa, S.C. Homeobox genes: A molecular link between development and cancer. Pesquisa Odontol. Bras. 2003, 17, 94–98. [Google Scholar] [CrossRef] [Green Version]
  84. Samuel, S.; Naora, H. Homeobox gene expression in cancer: Insights from developmental regulation and deregulation. Eur. J. Cancer 2005, 41, 2428–2437. [Google Scholar] [CrossRef]
  85. Pyper, S.R.; Viswakarma, N.; Jia, Y.; Zhu, Y.J.; Fondell, J.D.; Reddy, J.K. PRIC295, a Nuclear Receptor Coactivator, Identified from PPARalpha-Interacting Cofactor Complex. PPAR Res. 2010, 2010, 173907. [Google Scholar] [CrossRef] [Green Version]
  86. Viswakarma, N.; Jia, Y.; Bai, L.; Gao, Q.; Lin, B.; Zhang, X.; Misra, P.; Rana, A.; Jain, S.; Gonzalez, F.J.; et al. The Med1 subunit of the mediator complex induces liver cell proliferation and is phosphorylated by AMP kinase. J. Biol. Chem. 2013, 288, 27898–27911. [Google Scholar] [CrossRef] [Green Version]
  87. Ahmed, N.; Escalona, R.; Leung, D.; Chan, E.; Kannourakis, G. Tumour microenvironment and metabolic plasticity in cancer and cancer stem cells: Perspectives on metabolic and immune regulatory signatures in chemoresistant ovarian cancer stem cells. Semin. Cancer Biol. 2018, 53, 265–281. [Google Scholar] [CrossRef]
  88. Levin, A.; Minis, A.; Lalazar, G.; Rodriguez, J.; Steller, H. PSMD5 Inactivation Promotes 26S Proteasome Assembly during Colorectal Tumor Progression. Cancer Res. 2018, 78, 3458–3468. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  89. Chai, F.L.Y.; Bi, J.; Chen, L.Z.F.; Cui, Y.; Bian, X.; Jiang, J. High expression of REGγ is associated with metastasis and poor prognosis of patients with breast cancer. Int. J. Clin. Exp. Pathol. 2014, 7, 7834. [Google Scholar] [PubMed]
  90. Chai, F.; Liang, Y.; Bi, J.; Chen, L.; Zhang, F.; Cui, Y.; Jiang, J. REGgamma regulates ERalpha degradation via ubiquitin-proteasome pathway in breast cancer. Biochem. Biophys. Res. Commun. 2015, 456, 534–540. [Google Scholar] [CrossRef]
  91. Wang, X.; Tu, S.; Tan, J.; Tian, T.; Ran, L.; Rodier, J.F.; Ren, G. REG gamma: A potential marker in breast cancer and effect on cell cycle and proliferation of breast cancer cell. Med. Oncol. 2011, 28, 31–41. [Google Scholar] [CrossRef] [PubMed]
  92. Coppock, D.; Kopman, C.; Gudas, J.; Cina-Poppe, D.A. Regulation of the quiescence-induced genes: Quiescin Q6, decorin, and ribosomal protein S29. Biochem. Biophys. Res. Commun. 2000, 269, 604–610. [Google Scholar] [CrossRef]
  93. Notterman, D.A.; Alon, U.; Sierk, A.J.; Levine, A.J. Transcriptional gene expression profiles of colorectal adenoma, adenocarcinoma, and normal tissue examined by oligonucleotide arrays. Cancer Res. 2001, 61, 7. [Google Scholar]
  94. Lv, Z.; Wang, Z.; Luo, L.; Chen, Y.; Han, G.; Wang, R.; Xiao, H.; Li, X.; Hou, C.; Feng, J.; et al. Spliceosome protein Eftud2 promotes colitis-associated tumorigenesis by modulating inflammatory response of macrophage. Mucosal Immunol. 2019, 12, 1164–1173. [Google Scholar] [CrossRef]
  95. Sato, N.; Maeda, M.; Sugiyama, M.; Ito, S.; Hyodo, T.; Masuda, A.; Tsunoda, N.; Kokuryo, T.; Hamaguchi, M.; Nagino, M.; et al. Inhibition of SNW1 association with spliceosomal proteins promotes apoptosis in breast cancer cells. Cancer Med. 2015, 4, 268–277. [Google Scholar] [CrossRef] [PubMed]
  96. Pons, V.; Luyet, P.P.; Morel, E.; Abrami, L.; van der Goot, F.G.; Parton, R.G.; Gruenberg, J. Hrs and SNX3 functions in sorting and membrane invagination within multivesicular bodies. PLoS Biol. 2008, 6, e214. [Google Scholar] [CrossRef]
  97. Mendelsohn, J. The epidermal growth factor receptor as a target for cancer therapy. Endocr. Relat. Cancer 2001, 8, 3–9. [Google Scholar] [CrossRef] [Green Version]
  98. John, A.S.; Rothman, V.L.; Tuszynski, G.P. Thrombospondin-1 (TSP-1) Stimulates Expression of Integrin alpha6 in Human Breast Carcinoma Cells: A Downstream Modulator of TSP-1-Induced Cellular Adhesion. J. Oncol. 2010, 2010, 645376. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  99. Jayachandran, A.; Anaka, M.; Prithviraj, P.; Hudson, C.; McKeown, S.J.; Lo, P.H.; Vella, L.J.; Goding, C.R.; Cebon, J.; Behren, A. Thrombospondin 1 promotes an aggressive phenotype through epithelial-to-mesenchymal transition in human melanoma. Oncotarget 2014, 5, 5782–5797. [Google Scholar] [CrossRef] [Green Version]
  100. Tzeng, H.T.; Tsai, C.H.; Yen, Y.T.; Cheng, H.C.; Chen, Y.C.; Pu, S.W.; Wang, Y.S.; Shan, Y.S.; Tseng, Y.L.; Su, W.C.; et al. Dysregulation of Rab37-Mediated Cross-talk between Cancer Cells and Endothelial Cells via Thrombospondin-1 Promotes Tumor Neovasculature and Metastasis. Clin. Cancer Res. 2017, 23, 2335–2345. [Google Scholar] [CrossRef] [Green Version]
  101. Shi, H.; Hayes, M.; Kirana, C.; Miller, R.; Keating, J.; Macartney-Coxson, D.; Stubbs, R. TUFM is a potential new prognostic indicator for colorectal carcinoma. Pathology 2012, 44, 506–512. [Google Scholar] [CrossRef]
  102. Weng, X.; Zheng, S.; Shui, H.; Lin, G.; Zhou, Y. TUFM-knockdown inhibits the migration and proliferation of gastrointestinal stromal tumor cells. Oncol. Lett. 2020, 20, 250. [Google Scholar] [CrossRef]
  103. He, K.; Guo, X.; Liu, Y.; Li, J.; Hu, Y.; Wang, D.; Song, J. TUFM downregulation induces epithelial-mesenchymal transition and invasion in lung cancer cells via a mechanism involving AMPK-GSK3beta signaling. Cell. Mol. Life Sci. 2016, 73, 2105–2121. [Google Scholar] [CrossRef] [PubMed]
  104. Kulawiec, M.; Arnouk, H.; Desouki, M.M.; Kazim, L.; Still, I.; Singh, K.K. Proteomic analysis of mitochondria-to-nucleus retrograde response in human cancer. Cancer Biol. Ther. 2006, 5, 967–975. [Google Scholar] [CrossRef] [PubMed]
  105. Li, W.; Wubulikasimu, G.; Zhao, X.; Wang, C.; Liu, R.; Wang, L.; Zhu, X.; Chen, Z. UQCRC1 downregulation is correlated with lymph node metastasis and poor prognosis in CRC. Eur. J. Surg. Oncol. 2019, 45, 1005–1010. [Google Scholar] [CrossRef]
  106. Wang, Q.; Li, M.; Gan, Y.; Jiang, S.; Qiao, J.; Zhang, W.; Fan, Y.; Shen, Y.; Song, Y.; Meng, Z.; et al. Mitochondrial Protein UQCRC1 is Oncogenic and a Potential Therapeutic Target for Pancreatic Cancer. Theranostics 2020, 10, 2141–2157. [Google Scholar] [CrossRef] [PubMed]
  107. Mazure, N.M. VDAC in cancer. Biochim. Biophys. Acta Bioenerg. 2017, 1858, 665–673. [Google Scholar] [CrossRef]
  108. Sotgia, F.; Fiorillo, M.; Lisanti, M.P. Mitochondrial markers predict recurrence, metastasis and tamoxifen-resistance in breast cancer patients: Early detection of treatment failure with companion diagnostics. Oncotarget 2017, 8, 68730–68745. [Google Scholar] [CrossRef] [Green Version]
  109. Camara, A.K.S.; Zhou, Y.; Wen, P.C.; Tajkhorshid, E.; Kwok, W.M. Mitochondrial VDAC1: A Key Gatekeeper as Potential Therapeutic Target. Front. Physiol. 2017, 8, 460. [Google Scholar] [CrossRef] [Green Version]
  110. Magri, A.; Reina, S.; De Pinto, V. VDAC1 as Pharmacological Target in Cancer and Neurodegeneration: Focus on Its Role in Apoptosis. Front. Chem. 2018, 6, 108. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  111. Thul, P.J.; Akesson, L.; Wiking, M.; Mahdessian, D.; Geladaki, A.; Ait Blal, H.; Alm, T.; Asplund, A.; Bjork, L.; Breckels, L.M.; et al. A subcellular map of the human proteome. Science 2017, 356. [Google Scholar] [CrossRef]
  112. Uhlen, M.; Fagerberg, L.; Hallstrom, B.M.; Lindskog, C.; Oksvold, P.; Mardinoglu, A.; Sivertsson, A.; Kampf, C.; Sjostedt, E.; Asplund, A.; et al. Tissue-based map of the human proteome. Science 2015, 347, 1260419. [Google Scholar] [CrossRef] [PubMed]
  113. Uhlen, M.; Zhang, C.; Lee, S.; Sjostedt, E.; Fagerberg, L.; Bidkhori, G.; Benfeitas, R.; Arif, M.; Liu, Z.; Edfors, F.; et al. A pathology atlas of the human cancer transcriptome. Science 2017, 357. [Google Scholar] [CrossRef] [Green Version]
  114. Doncheva, N.T.; Morris, J.H.; Gorodkin, J.; Jensen, L.J. Cytoscape StringApp: Network Analysis and Visualization of Proteomics Data. J. Proteome Res. 2019, 18, 623–632. [Google Scholar] [CrossRef]
  115. Essex, H.E.; Priestley, J.T. Effect of Rattlesnake Venom on Flexner-Jobling’s Carcinoma in the White Rat (Mus Norvegicus Albinus.). Exp. Biol. Med. 1931, 28, 550–551. [Google Scholar] [CrossRef]
  116. Ferreira, S.H.; Greene, L.H.; Alabaster, V.A.; Bakhle, Y.S.; Vane, J.R. Activity of various fractions of bradykinin potentiating factor against angiotensin I converting enzyme. Nature 1970, 225, 379–380. [Google Scholar] [CrossRef] [PubMed]
  117. Cushman, D.W.; Ondetti, M.A. History of the design of captopril and related inhibitors of angiotensin converting enzyme. Hypertension 1991, 17, 589–592. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  118. Smith, C.G.; Vane, J.R. The discovery of captopril. FASEB J. 2003, 17, 788–789. [Google Scholar] [CrossRef] [PubMed]
  119. Ma, R.; Mahadevappa, R.; Kwok, H.F. Venom-based peptide therapy: Insights into anti-cancer mechanism. Oncotarget 2017, 8, 100908–100930. [Google Scholar] [CrossRef] [Green Version]
  120. Li, L.; Huang, J.; Lin, Y. Snake Venoms in Cancer Therapy: Past, Present and Future. Toxins 2018, 10, 346. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  121. Sukocheva, O.; Menschikowski, M.; Hagelgans, A.; Yarla, N.S.; Siegert, G.; Reddanna, P.; Bishayee, A. Current insights into functions of phospholipase A2 receptor in normal and cancer cells: More questions than answers. Semin. Cancer Biol. 2019, 56, 116–127. [Google Scholar] [CrossRef]
  122. Peng, Z.; Chang, Y.; Fan, J.; Ji, W.; Su, C. Phospholipase A2 superfamily in cancer. Cancer Lett. 2021, 497, 165–177. [Google Scholar] [CrossRef]
  123. Cathcart, J.; Pulkoski-Gross, A.; Cao, J. Targeting Matrix Metalloproteinases in Cancer: Bringing New Life to Old Ideas. Genes Dis. 2015, 2, 26–34. [Google Scholar] [CrossRef] [Green Version]
  124. Di Cara, G.; Marabeti, M.R.; Musso, R.; Riili, I.; Cancemi, P.; Pucci Minafra, I. New Insights into the Occurrence of Matrix Metalloproteases-2 and -9 in a Cohort of Breast Cancer Patients and Proteomic Correlations. Cells 2018, 7, 89. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  125. Quintero-Fabian, S.; Arreola, R.; Becerril-Villanueva, E.; Torres-Romero, J.C.; Arana-Argaez, V.; Lara-Riegos, J.; Ramirez-Camacho, M.A.; Alvarez-Sanchez, M.E. Role of Matrix Metalloproteinases in Angiogenesis and Cancer. Front. Oncol. 2019, 9, 1370. [Google Scholar] [CrossRef] [Green Version]
  126. Bezerra, P.H.A.; Ferreira, I.M.; Franceschi, B.T.; Bianchini, F.; Ambrosio, L.; Cintra, A.C.O.; Sampaio, S.V.; de Castro, F.A.; Torqueti, M.R. BthTX-I from Bothrops jararacussu induces apoptosis in human breast cancer cell lines and decreases cancer stem cell subpopulation. J. Venom. Anim. Toxins Incl. Trop. Dis. 2019, 25, e20190010. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  127. Murray, A.S.; Hyland, T.E.; Sala-Hamrick, K.E.; Mackinder, J.R.; Martin, C.E.; Tanabe, L.M.; Varela, F.A.; List, K. The cell-surface anchored serine protease TMPRSS13 promotes breast cancer progression and resistance to chemotherapy. Oncogene 2020, 39, 6421–6436. [Google Scholar] [CrossRef] [PubMed]
  128. Krasny, L.; Huang, P.H. Data-independent acquisition mass spectrometry (DIA-MS) for proteomic applications in oncology. Mol. Omics 2021, 17, 29–42. [Google Scholar] [CrossRef] [PubMed]
  129. Zhao, W.; Li, J.; Chen, M.M.; Luo, Y.; Ju, Z.; Nesser, N.K.; Johnson-Camacho, K.; Boniface, C.T.; Lawrence, Y.; Pande, N.T.; et al. Large-Scale Characterization of Drug Responses of Clinically Relevant Proteins in Cancer Cell Lines. Cancer Cell 2020, 38, 829–843.e4. [Google Scholar] [CrossRef]
  130. Klaeger, S.; Heinzlmeir, S.; Wilhelm, M.; Polzer, H.; Vick, B.; Koenig, P.A.; Reinecke, M.; Ruprecht, B.; Petzoldt, S.; Meng, C.; et al. The target landscape of clinical kinase drugs. Science 2017, 358. [Google Scholar] [CrossRef] [Green Version]
  131. Qi, Z.; Zhang, L.; Chen, Y.; Ma, X.; Gao, X.; Du, J.; Zhang, F.; Cheng, X.; Cui, W. Biological variations of seven tumor markers. Clin. Chim. Acta 2015, 450, 233–236. [Google Scholar] [CrossRef]
  132. Lamond, A.I.; Mann, M. Cell biology and the genome projects a concerted strategy for characterizing multiprotein complexes by using mass spectrometry. Trends Cell Biol. 1997, 7, 139–142. [Google Scholar] [CrossRef]
  133. Tavares, C.; Maciel, T.; Burin, S.; Ambrosio, L.; Ghisla, S.; Sampaio, S.; Castro, F. l-Amino acid oxidase isolated from Calloselasma rhodostoma snake venom induces cytotoxicity and apoptosis in JAK2V617F-positive cell lines. Rev. Bras. Hematol. Hemoter. 2016, 38, 128–134. [Google Scholar] [CrossRef] [Green Version]
  134. Zakraoui, O.; Marcinkiewicz, C.; Aloui, Z.; Othman, H.; Grepin, R.; Haoues, M.; Essafi, M.; Srairi-Abid, N.; Gasmi, A.; Karoui, H.; et al. Lebein, a snake venom disintegrin, suppresses human colon cancer cells proliferation and tumor-induced angiogenesis through cell cycle arrest, apoptosis induction and inhibition of VEGF expression. Mol. Carcinog. 2017, 56, 18–35. [Google Scholar] [CrossRef] [PubMed]
  135. Bernardes-Oliveira, E.; Gomes, D.L.; Martelli Palomino, G.; Juvenal Silva Farias, K.; da Silva, W.D.; Rocha, H.A.; Goncalves, A.K.; Fernandes-Pedrosa, M.F.; Crispim, J.C. Bothrops jararaca and Bothrops erythromelas Snake Venoms Promote Cell Cycle Arrest and Induce Apoptosis via the Mitochondrial Depolarization of Cervical Cancer Cells. Evid. Based Complement. Altern. Med. 2016, 2016, 1574971. [Google Scholar] [CrossRef]
  136. Knorre, D.G.; Kudryashova, N.V.; Godovikova, T.S. Chemical and functional aspects of posttranslational modification of proteins. Acta Nat. 2009, 1, 23. [Google Scholar]
  137. Freed, E.F.; Bleichert, F.; Dutca, L.M.; Baserga, S.J. When ribosomes go bad: Diseases of ribosome biogenesis. Mol. Biosyst 2010, 6, 481–493. [Google Scholar] [CrossRef]
  138. Ruggero, D.; Pandolfi, P.P. Does the ribosome translate cancer? Nat. Rev. Cancer 2003, 3, 179–192. [Google Scholar] [CrossRef]
  139. Bilanges, B.; Stokoe, D. Mechanisms of translational deregulation in human tumors and therapeutic intervention strategies. Oncogene 2007, 26, 5973–5990. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  140. Hershko, A. Ubiquitin: Roles in protein modification and breakdown. Cell 1983, 34, 11–12. [Google Scholar] [CrossRef]
  141. Hershko, A.; Ciechanover, A. The ubiquitin system. Annu. Rev. Biochem. 1998, 67, 425–479. [Google Scholar] [CrossRef]
  142. Jankowska, E.; Stoj, J.; Karpowicz, P.; Osmulski, P.A.; Gaczynska, M. The proteasome in health and disease. Curr. Pharm. Des. 2013, 19, 19. [Google Scholar]
  143. Pajonk, F.; McBride, W.H. The Proteasome in Cancer Biology and Treatment. Radiat. Res. 2001, 156, 447–459. [Google Scholar] [CrossRef]
  144. Teicher, B.A.; Ara, G.; Herbst, R.; Palombella, V.J.; Adams, J. The proteasome inhibitor PS-341 in cancer therapy. Clin. Cancer Res. 1999, 5, 2638–2645. [Google Scholar] [PubMed]
  145. Thul, P.J.; Lindskog, C. The human protein atlas: A spatial map of the human proteome. Protein Sci. 2018, 27, 233–244. [Google Scholar] [CrossRef] [Green Version]
  146. Dias, M.H.; Fonseca, C.S.; Zeidler, J.D.; Albuquerque, L.L.; da Silva, M.S.; Cararo-Lopes, E.; Reis, M.S.; Noel, V.; Dos Santos, E.O.; Prior, I.A.; et al. Fibroblast Growth Factor 2 lethally sensitizes cancer cells to stress-targeted therapeutic inhibitors. Mol. Oncol. 2019, 13, 290–306. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  147. Montoni, F.; Andreotti, D.Z.; Eichler, R.; Santos, W.D.S.; Kisaki, C.Y.; Arcos, S.S.S.; Lima, I.F.; Soares, M.A.M.; Nishiyama-Jr, M.Y.; Nava-Rodrigues, D.; et al. The impact of rattlesnake venom on mice cerebellum proteomics points to synaptic inhibition and tissue damage. J. Proteom. 2020, 221, 103779. [Google Scholar] [CrossRef]
  148. Mi, H.; Muruganujan, A.; Huang, X.; Ebert, D.; Mills, C.; Guo, X.; Thomas, P.D. Protocol Update for large-scale genome and gene function analysis with the PANTHER classification system (v.14.0). Nat. Protoc. 2019, 14, 703–721. [Google Scholar] [CrossRef] [PubMed]
  149. Mi, H.; Poudel, S.; Muruganujan, A.; Casagrande, J.T.; Thomas, P.D. PANTHER version 10: Expanded protein families and functions, and analysis tools. Nucleic Acids Res. 2016, 44, D336–D342. [Google Scholar] [CrossRef] [Green Version]
  150. Thomas, P.D.; Campbell, M.J.; Kejariwal, A.; Mi, H.; Karlak, B.; Daverman, R.; Diemer, K.; Muruganujan, A.; Narechania, A. PANTHER: A library of protein families and subfamilies indexed by function. Genome Res. 2003, 13, 2129–2141. [Google Scholar] [CrossRef] [Green Version]
  151. Thomas, P.D.; Kejariwal, A.; Campbell, M.J.; Mi, H.; Diemer, K.; Guo, N.; Ladunga, I.; Ulitsky-Lazareva, B.; Muruganujan, A.; Rabkin, S.; et al. PANTHER: A browsable database of gene products organized by biological function, using curated protein family and subfamily classification. Nucleic Acids Res. 2003, 31, 334–341. [Google Scholar] [CrossRef]
  152. Cox, J.; Hein, M.Y.; Luber, C.A.; Paron, I.; Nagaraj, N.; Mann, M. Accurate proteome-wide label-free quantification by delayed normalization and maximal peptide ratio extraction, termed MaxLFQ. Mol. Cell Proteom. 2014, 13, 2513–2526. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  153. 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] [PubMed]
  154. Szklarczyk, D.; Franceschini, A.; Wyder, S.; Forslund, K.; Heller, D.; Huerta-Cepas, J.; Simonovic, M.; Roth, A.; Santos, A.; Tsafou, K.P.; et al. STRING v10: Protein-protein interaction networks, integrated over the tree of life. Nucleic Acids Res. 2015, 43, D447–D452. [Google Scholar] [CrossRef]
  155. Lê, S.; Josse, J.; Husson, F. FactoMineR: AnRPackage for Multivariate Analysis. J. Stat. Softw. 2008, 25. [Google Scholar] [CrossRef] [Green Version]
Toxins 13 00519 g001 550
Figure 1. Cytotoxicity assay of (a) MCF7 and (b) MDA-MB-231 cell lines treated with B. jararaca snake venom ranging from 0 to 20 μg/mL for 24 h. Experiment was performed using the WST-1 reagent kit.
Figure 1. Cytotoxicity assay of (a) MCF7 and (b) MDA-MB-231 cell lines treated with B. jararaca snake venom ranging from 0 to 20 μg/mL for 24 h. Experiment was performed using the WST-1 reagent kit.
Toxins 13 00519 g001
Toxins 13 00519 g002 550
Figure 2. Diagrammatic representation of the comparative analysis of identified proteins under the three conditions: no venom control group, low venom dose at 0.63 µg/mL and high venom sub-toxic dose of 2.5 µg/mL of B. jararaca venom. (a) MCF7 and (b) MDA-MB-231.
Figure 2. Diagrammatic representation of the comparative analysis of identified proteins under the three conditions: no venom control group, low venom dose at 0.63 µg/mL and high venom sub-toxic dose of 2.5 µg/mL of B. jararaca venom. (a) MCF7 and (b) MDA-MB-231.
Toxins 13 00519 g002
Toxins 13 00519 g003 550
Figure 3. Hierarchical clustering of differentially expressed proteins detected in both MCF7 and MDA-MB-231 cells treated with low (0.63 µg/mL) and high (2.5 µg/mL) B. jararaca venom for 24 h. (a) Heatmap representation of the hierarchical clustering of proteins detected in both cell lines with quantification in at least two replicates showing the changes in protein abundance. The protein fold change is log2 transformed and normalized with mean-centering scale. (b) Protein Clusters extracted from the hierarchical clustering. X axis: Cell types treated with different B. jararaca venom concentrations (MCF7 0.63 μg/mL; MCF7 2.5 μg/mL, MDA-MB-231 0.63 μg/mL, MDA-MB-231 2.5 μg/mL); Y axis: mean-centered log2 Fold Change. Grey lines: individual proteins; Black line: average expression values per cluster.
Figure 3. Hierarchical clustering of differentially expressed proteins detected in both MCF7 and MDA-MB-231 cells treated with low (0.63 µg/mL) and high (2.5 µg/mL) B. jararaca venom for 24 h. (a) Heatmap representation of the hierarchical clustering of proteins detected in both cell lines with quantification in at least two replicates showing the changes in protein abundance. The protein fold change is log2 transformed and normalized with mean-centering scale. (b) Protein Clusters extracted from the hierarchical clustering. X axis: Cell types treated with different B. jararaca venom concentrations (MCF7 0.63 μg/mL; MCF7 2.5 μg/mL, MDA-MB-231 0.63 μg/mL, MDA-MB-231 2.5 μg/mL); Y axis: mean-centered log2 Fold Change. Grey lines: individual proteins; Black line: average expression values per cluster.
Toxins 13 00519 g003
Toxins 13 00519 g004 550
Figure 4. Comparison of protein log2 Fold Change profiles across treated cell lines. Principal component analysis in a 2D graph represented by the first two components PC1 and PC2 explains 67.7% of the protein variability among the different conditions. Vectors that are closer are highly correlated. Vectors representing the conditions which are orthogonal or well-spaced in terms of the observed proteome indicate that those proteins can be closely related to each specific cell line condition.
Figure 4. Comparison of protein log2 Fold Change profiles across treated cell lines. Principal component analysis in a 2D graph represented by the first two components PC1 and PC2 explains 67.7% of the protein variability among the different conditions. Vectors that are closer are highly correlated. Vectors representing the conditions which are orthogonal or well-spaced in terms of the observed proteome indicate that those proteins can be closely related to each specific cell line condition.
Toxins 13 00519 g004
Toxins 13 00519 g005a 550Toxins 13 00519 g005b 550
Figure 5. Protein–protein interaction of proteins identified in (a) MCF7 and (b) MDA-MB-231 cell lines presenting FC ≥ 1.5 at 2.5 g/mL B. jararaca venom treatment.
Figure 5. Protein–protein interaction of proteins identified in (a) MCF7 and (b) MDA-MB-231 cell lines presenting FC ≥ 1.5 at 2.5 g/mL B. jararaca venom treatment.
Toxins 13 00519 g005aToxins 13 00519 g005b
Table
Table 1. Description of the highlighted proteins and their association with cancer.
Table 1. Description of the highlighted proteins and their association with cancer.
ProteinProtein NameProtein Description Related to CancerCancer Type AssociationReferences
AMOTAngiomotinPlays a central role in tight junction maintenance. Appears to regulate endothelial cell migration and tube formation. May also play a role in the assembly of endothelial cell-cell junctions. Plays a critical role in angiogenesis, proliferation and migration and invasion of cancer cellsBL, BR, CE, CL, CR, EN, HN, KD, LE, LI, LA, LS, OV, PR, ST[49,50,51]
ATP5PDATP synthase peripheral stalk subunit DMitochondrial ATP synthase catalyzes ATP synthesis, utilizing an electrochemical gradient of protons across the inner membrane during oxidative phosphorylation. Linked to failure of therapy, disease progression, and poor survival in patients with cancer. High expression of ATP5PD has been observed in several types of cancerBR, CE, CL, EN, GL, HN, LI, LU, LY, ME, OV, PA, PR, RE, SK, TE, TY[52]
ATRSerine/threonine-protein kinasePlays important roles for cell survival and is considered a major mediator of DNA response in human cells, preventing cells with damaged or incompletely replicated DNA from entering mitosis when cells are damaged by radiotherapy or chemotherapy during cancer treatmentBL, BR, CE, CR, EN, GB, HN, KD, LU, LA, LS, OV, ST, TY[53,54]
CYR61Cysteine-rich heparin-binding protein 61Plays an important role in cell proliferation, survival, chemotaxis, angiogenesis, adhesion, and migration of different types of cells. Participate in key different cellular events during vascular development, angiogenesis, wound healing and the development and progression of various types of cancersBN, BR, CR, EN, GA, GB, GL, LI, LU, OV, PA, PR, ST, TE, UR[55,56]
GSTM3The glutathione S-transferase Mu 3Part of the GSTs enzymes that have functions such as immunological system evasion and inhibition of apoptosis. Involved in prostaglandin and leukotriene synthesis and metabolization of both endogenous compounds and xenobiotics such as chemotherapeutic drugs, insecticides, carcinogens, and oxidative stress byproductsBL, BR, CR, EN, LE, LA, LS, OV, PA, ST, TY, UR[57,58]
H3F3B/H3C15H3 histone family member 3BCore component of nucleosome. Histones play a central role in transcription regulation, DNA repair, DNA replication and chromosomal stability and are related to different types of cancer. H3F3B mutation has been described to lead to some human cancersBL, BN, BR, CE, CH, CR, EN, GB, HN, LU, OV, UR[59,60,61]
HEL-S-156an PNPPurine nucleoside phosphorylaseCatalyze the phosphorolysis of purine nucleosides. Mutations which result in nucleoside phosphorylase deficiency result in defective T-cell (cell-mediated) immunity but can also affect B-cell immunity and antibody responses. High expression of PNP has been observed in several types of cancerBR, CL, CR, GA, GL, KD, LI, LU, LA, LY, ME, OV, PR, TY[62,63,64]
HIST1H4JHistone H4Core component of nucleosome. Histones play a central role in transcription regulation, DNA repair, DNA replication and chromosomal stability. Post-translational alterations of histones have been shown to affect the activation and repression of oncogenes and tumor suppressor genesBR, CE, CR, GL, HN, LI, LU, ME, OV, PA, PR, SK, ST, TE, TY[65,66,67,68]
IGF2BP1Insulin like growth factor 2 MRNA binding protein 1RNA-binding factor that recruits target transcripts to cytoplasmic protein–RNA complexes (mRNPs). IGF2BP1 has an oncogenic role, characterized by changes in actin dynamics, migration, invasion, proliferation, and self-renewal. Play a role in resistance to drugsBR, CE, CR, EN, GB, HN, LI, LU, LA, LS, ME, OV, PR, ST, TE, UR[69,70,71]
KRT1KeratinMay regulate the activity of kinases such as PKC and SRC via binding to integrin beta-1 (ITB1) and to the receptor of activated protein C kinase 1 (RACK1). High expression of KRT1 protein has been observed in several types of cancer and is correlated with advanced melanoma tumor stage and infiltration of immune cellsBR, CE, CR, EN, GB, HN, KD, LI, LS, LA, OV, SK, ST, UR[72,73]
LAP3Leucine aminopeptidase 3Cytosolic metallopeptidase that catalyzes the removal of unsubstituted N-terminal hydrophobic amino acids from various peptides. Involved in the metabolism of glutathione and in the degradation of glutathione S-conjugates, which may play a role in the control of the cell redox status. Related to protein renewal. Have a potential for determining the prognosis for breast cancerBL, BR, CR, EN, HN, KD, LI, LA, LS, PA, ST, TY[74,75,76,77]
MCCC2Methylcrotonoyl-CoA carboxylase beta chain, mitochondrialEnzyme that catalyzes the conversion of 3-methylcrotonyl-CoA to 3-methylglutaconyl-CoA, a critical step for leucine and isovaleric acid catabolism. Overexpression of MCCC2 is associated with tumor stage, node, metastasis, lymph node metastasis and predicts unfavorable prognosis. Additionally, involved in the development and formation of some tumors, such as breast cancerBL, BR, CR, EN, HN, KD, LI, LA, OV, PA, PR, ST, TY[78]
MTCH2Mitochondrial carrier homolog 2Member of the SLC25 family of nuclear-encoded transporters that are localized in the inner mitochondrial membrane. Members of this superfamily are involved in many metabolic pathways and cell functions. Associated with metastasis and tumor cell survival. Indirect involvement in the expression of miR-135b mRNA, which is one of the proteins responsible for tumorigenicityBL, BR, CR, EN, HN, KD, LI, LA, ME, OV, PA, PR, ST, TY[79,80,81]
PCBD4a-hydroxytetrahydrobiopterin dehydrataseInvolved in tetrahydrobiopterin biosynthesis. Regulates various aspects of cell morphogenesis and differentiation as a cofactor for the homeobox transcription factor. Several types of cancer show expression or alteration in the homeobox genes. PCBD degradation increases cell survival and proliferation, and inhibits tumor cell differentiationBR, CL, CR, EN, LE, LI, LU, OV, PA, PR, RE, SK[82,83,84]
PRIC295Peroxisome proliferator-activated receptor-a (PPARα)-interacting cofactorFunctions as a transcriptional coactivator for nuclear receptors. Enhances the activation of PPARα and PPARγ and plays a key role in lipid metabolism and energy combustion regulating the genes for fatty acid oxidation. Observed to be significantly enhanced in chemotherapy recurrence when compared to chemotherapy treatment in ovarian cancer patientsBL, BR, CR, EN, GB, HN, KD, LE, LA, LS, ME, OV, PA, PR, ST, TY[85,86,87]
PSMD5The 26S proteasome non-ATPase regulatory subunit 5Acts as a chaperone during the assembly of the 26S proteasome. Expression reduced in several types of cancer including intestinal and colorectal tumorsBL, BR, CR, EN, HN, PR, ST, TY[88]
PSME2The proteasome activator complex subunit 2Implicated in immunoproteasome assembly and required for efficient antigen processing. Member of the PSME family that regulates proteasome function. Elevated expression of PSME have also been associated with several types of cancerBR, CR, EN, HN, LC, LA, LS, ME, PR, ST[89,90,91]
RPS29Ribosomal protein S29Belongs to the universal ribosomal protein uS14 family. Related to have tumor suppressor activity for ras-transformed NIH3T3 cells. High expression RPS29 mRNA levels observed in adenomasBR, CE, CR, EN, HN, KD, LA, OV, ST[92,93]
SNRP116/EFTUD2Small nuclear ribonucleoprotein componentRequired for pre-mRNA splicing as a component of the spliceosome, including pre-catalytic, catalytic, and post-catalytic spliceosomal complexes. Knockout of EFTUD2 suppressed the development and tumor progression due to impaired activation of NF-kB signaling in macrophagesBL, BR, CR, EN, GB, HN, KD, LA, LS, ME, OV, PR, ST, TE, TY, UR[94,95]
SNX3The sorting nexin 3Phosphoinositide-binding protein required for multivesicular body formation. Plays a role in protein transport between cellular compartments. The knockdown of SNX3 is associated with degradation of the EGF receptor which is related to resistance to chemotherapy and radiotherapyBR, CR, EN, GL, HN, KD, LU, OV, RE, TY[96,97]
THBS1Thrombospondin 1Adhesive glycoprotein that mediates cell-to-cell and cell-to-matrix interactions. Influences angiogenesis modulation by regulating adhesion, invasion, metastasis, migration, proliferation, and apoptosis and has been implicated in numerous types of cancersBR, CR, EN, ES, HN, KD, LA, LS, LY, ME, OV, PA, PR, SK, ST, TE, TY[98,99,100]
TUFMTu translation elongation factor, mitochondrialPromotes the GTP-dependent binding of aminoacyl-tRNA to the A-site of ribosomes during protein biosynthesis. Plays important roles in the regulation of autophagy and innate immunity. TUFM is highly expressed in several types of cancersBR, CR, EN, ES, GA, GS, HN, LI, LU, LA, ME, PA, PR, RE, SK, ST, TE[101,102,103]
UQCRC1The ubiquinol-cytochrome C reductase core protein 1Component of the ubiquinol-cytochrome c reductase complex, which is part of the mitochondrial respiratory chain. High expression was observed in several types of cancer. Negative expression correlated significantly with clinical and pathological parameters including tumor stage, vascular invasion, and lymph node metastasis, suggesting that the reduction of this protein is associated with tumor progressionBR, CR, EN, GB, GL, HN, LI, LA, LY, ME, OV, PA, PR, RE, ST, TE, TY, UR[104,105,106]
VDAC1 and VDAC2The voltage-dependent anion selective channel 1 and 2Forms a channel through the mitochondrial outer membrane and plasma membrane allowing diffusion of small hydrophilic molecules. In the plasma membrane it is involved in cell volume regulation and apoptosis. The abnormal expression or mal functioning of VDACs has been reported in multiple tumors and it has been considered as a biomarker capable of predicting treatment failure and breast cancer recurrenceBL, BR, CR, EN, GB, HN, LI, LA, LS, ME, OV, PA, PR, RE, ST, TY, UR[107,108,109,110]
Note: This description may include information from UniProtKB, PhosphoSitePlus (v.6.5.9.3), GeneCards (Human Gene Database), and Human Protein Atlas [111,112,113]. Bladder: BL, bone (osteosarcoma): BN, breast: BR, cervical: CE, chondroblastoma: CH colon: CL, colorectal: CR, endometrial: EN, esophageal: ES, gastric: GA, gastrointestinal stromal tumor: GS glioblastoma: GB, glioma: GL, head and neck: HN, kidney: KD, laryngeal carcinoma: LC, leukemia: LE, liver: LI, lung: LU, lung adenocarcinoma: LA, lung squamous: LS, lymphoma: LY, melanoma: ME, ovarian: OV, pancreatic: PA, prostate: PR, renal: RE, skin: SK, stomach: ST, testis: TE, thyroid: TY, urothelial: UR.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kisaki, C.Y.; Arcos, S.S.S.; Montoni, F.; da Silva Santos, W.; Calacina, H.M.; Lima, I.F.; Cajado-Carvalho, D.; Ferro, E.S.; Nishiyama-Jr, M.Y.; Iwai, L.K. Bothrops Jararaca Snake Venom Modulates Key Cancer-Related Proteins in Breast Tumor Cell Lines. Toxins 2021, 13, 519. https://doi.org/10.3390/toxins13080519

AMA Style

Kisaki CY, Arcos SSS, Montoni F, da Silva Santos W, Calacina HM, Lima IF, Cajado-Carvalho D, Ferro ES, Nishiyama-Jr MY, Iwai LK. Bothrops Jararaca Snake Venom Modulates Key Cancer-Related Proteins in Breast Tumor Cell Lines. Toxins. 2021; 13(8):519. https://doi.org/10.3390/toxins13080519

Chicago/Turabian Style

Kisaki, Carolina Yukiko, Stephanie Santos Suehiro Arcos, Fabio Montoni, Wellington da Silva Santos, Hamida Macêdo Calacina, Ismael Feitosa Lima, Daniela Cajado-Carvalho, Emer Suavinho Ferro, Milton Yutaka Nishiyama-Jr, and Leo Kei Iwai. 2021. "Bothrops Jararaca Snake Venom Modulates Key Cancer-Related Proteins in Breast Tumor Cell Lines" Toxins 13, no. 8: 519. https://doi.org/10.3390/toxins13080519

APA Style

Kisaki, C. Y., Arcos, S. S. S., Montoni, F., da Silva Santos, W., Calacina, H. M., Lima, I. F., Cajado-Carvalho, D., Ferro, E. S., Nishiyama-Jr, M. Y., & Iwai, L. K. (2021). Bothrops Jararaca Snake Venom Modulates Key Cancer-Related Proteins in Breast Tumor Cell Lines. Toxins, 13(8), 519. https://doi.org/10.3390/toxins13080519

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