# Antiproliferative Effects of Naja anchietae and Naja senegalensis Venom Peptides on Glioblastoma Cell Lines

Yasmine Boughanmi, Caroline Berenguer-Daizé, Marielle Balzano, Hend Mosrati, Maxime Moulard, Pascal Mansuelle, Patrick Fourquet, Franck Torre, Harold de Pomyers, Didier Gigmes, et al.

► To cite this version:

Yasmine Boughanmi, Caroline Berenguer-Daizé, Marielle Balzano, Hend Mosrati, Maxime Moulard, et al. Antiproliferative Effects of Naja anchietae and Naja senegalensis Venom Peptides on Glioblastoma Cell Lines. Toxins, 2024, 16 (10), pp.433. 10.3390/toxins16100433. hal-04747350

HAL Id: hal-04747350

https://hal.science/hal-04747350v1

Submitted on 22 Oct 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Article

# Antiproliferative Effects of *Naja anchietae* and *Naja senegalensis* Venom Peptides on Glioblastoma Cell Lines

Yasmine Boughanmi<sup>1,2</sup>, Caroline Berenguer-Daizé<sup>3</sup>, Marielle Balzano<sup>4</sup>, Hend Mosrati<sup>1</sup>, Maxime Moulard<sup>4</sup>, Pascal Mansuelle<sup>5 ID</sup>, Patrick Fourquet<sup>6 ID</sup>, Franck Torre<sup>7</sup>, Harold de Pomyers<sup>2</sup>, Didier Gigmes<sup>1</sup>, Lhoucine Ouafik<sup>3 ID</sup> and Kamel Mabrouk<sup>1,* ID</sup>

<figure><img src="image_3.png" /></figure>

check for
updates

**Citation:** Boughanmi, Y.; Berenguer-Daizé, C.; Balzano, M.; Mosrati, H.; Moulard, M.; Mansuelle, P.; Fourquet, P.; Torre, F.; de Pomyers, H.; Gigmes, D.; et al. Antiproliferative Effects of *Naja anchietae* and *Naja senegalensis* Venom Peptides on Glioblastoma Cell Lines. *Toxins* **2024**, *16*, 433. https://doi.org/10.3390/toxins16100433

Received: 31 July 2024

Revised: 3 October 2024

Accepted: 6 October 2024

Published: 10 October 2024

<figure><img src="image_4.png" /></figure>

**Copyright:** © 2024 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

**Abstract:** This study explores the potential of natural bioactive peptides from animal venoms as targeted anti-cancer agents with reduced toxicity. Initially, we screened a broad collection of animal venoms for their antiproliferative activity against cancer cell lines. From this collection, we selected venoms from *Naja anchietae* and *Naja senegalensis* due to their promising activity. Utilizing reverse-phase high-performance liquid chromatography (RP HPLC), mass spectrometry (MALDI-TOF MS and MALDI-TOF TOF MSMS), and Edman degradation sequencing, we isolated and characterized three peptides named CTNanc1, CTNanc2, and CTNanc3 from *Naja anchietae*, and three others named CTNsen1, CTNsen2, and CTNsen3 from *Naja senegalensis*, each with a molecular weight of around 7 kDa. These purified peptides demonstrated inhibition of U87 glioblastoma cell proliferation, but not of U251 and T98G cells, in cell viability assays. To assess the impact of these treatments on cell viability, apoptosis, and necrosis, flow cytometry assays were conducted on U87 cells at 72 h. The results showed a decrease in cell viability and an increase in dead cells, suggesting that the treatments not only promote apoptosis, but may also lead to increased necrosis or late-stage apoptosis as the exposure time increases. These findings suggest that these peptides could be developed as leads for cancer therapy.

**Keywords:** venom; *Naja senegalensis*; *Naja anchietae*; peptide; cobra cytotoxin; antiproliferative activity; apoptosis; necrosis

**Key Contribution:** Six novel cytotoxins isolated from *Naja anchietae* and *Naja senegalensis* venoms showed antiproliferative effects on glioblastoma cell lines. Our results suggest that several of these peptides exhibit selective cytotoxicity against glioblastoma cells while demonstrating reduced toxicity toward normal endothelial cells.

## 1. Introduction

Over the past three decades, cancer cases have continued to rise. The five-year relative survival rate for all types of cancer has also increased, partly due to targeted

<footnote n="1">1 Aix Marseille University, Institut de Chimie Radicalaire UMR 7273, 13397 Marseille, France; yasmine.boughanmi@etu.univ-amu.fr (Y.B.); mosrati_hend@yahoo.fr (H.M.); didier.gigmes@univ-amu.fr (D.G.)</footnote>

<footnote n="2">2 Latoxan, 26800 Portes-lès-Valence, France; harold.pomyers@latoxan.com</footnote>

<footnote n="3">3 Aix-Marseille University, INP—Institute of Neuropathophysiology—UMR 7051 CNRS, 13005 Marseille, France; caroline.berenguer@univ-amu.fr (C.B.-D.); lhoucine.ouafik@univ-amu.fr (L.O.)</footnote>

<footnote n="4">4 BioCytex 140 Chemin de l'Armée d'Afrique, 13010 Marseille, France; marielle.balzano@biocytex.fr (M.B.); maxime.moulard@biocytex.fr (M.M.)</footnote>

<footnote n="5">5 Proteomics Platform, Marseille Proteomics (MaP), Institut de Microbiologie de la Méditerranée (IMM), FR 3479, CNRS, 13009 Marseille, France; pmansuelle@imm.cnrs.fr</footnote>

<footnote n="6">6 Aix-Marseille University, INSERM, CNRS, Institut Paoli-Calmettes, CRCM, Marseille Proteomics, 13009 Marseille, France; patrick.fourquet@inserm.fr</footnote>

<footnote n="7">7 Aix-Marseille University, Institut Méditerranéen de Biodiversité et d'Ecologie (IMBE), Avignon University, CNRS, IRD, 13397 Marseille, France; franck.torre@univ-amu.fr</footnote>

<footnote>* Correspondence: kamel.mabrouk@univ-amu.fr</footnote>
therapies. However, these targeted therapy drugs still have limitations. Therefore, it is crucial to continue developing new biological therapeutic agents that could be used alone or in synergy with FDA-approved drugs to combat cancer. Venoms are a treasure trove of bioactive molecules that have shown potential in treating several diseases, including cancer.

Cytotoxins (CTX), major components of cobra venom, account for approximately 40 to 60% of the total protein content within these venoms. These molecules are typically composed of 59 to 62 amino acids and have a molecular weight ranging between 6 and 9 kDa. The amino acid sequences of these cytotoxins generally share over 60% homology yet exhibit a diverse array of activities. This suggests that even minor differences in amino acid sequences can significantly influence toxin interactions, highlighting the complex nature of these biomolecules [1,2]. The UniProt protein database records over a hundred sequences of CTXs, primarily identified in venoms from snakes of the *Naja* genus, commonly known as cobras. These sequences were discovered predominantly in cobra venom, emphasizing the biological significance and diversity of these toxins [3].

Structurally, CTXs are characterized by their three-dimensional arrangement consisting of five antiparallel β-sheets, which form three functional hydrophobic loops. This distinct structure categorizes them within the three-finger toxin (3FTx) family. The stability of this structure is largely due to the presence of four disulfide bridges that form a robust hydrophobic core, essential for the structural integrity of the proteins [4,5].

The ends of loops I to III in cytotoxins are predominantly composed of hydrophobic residues, which confer amphiphilic properties to these molecules. Around this hydrophobic zone lies a belt of charged residues, dominated by lysine amino acids. This unique structural feature of the molecules, combining hydrophobic and cationic regions, plays a crucial role in their interaction with cellular membranes [6,7]. Notably, loop III exhibits more pronounced structural stability compared to the other loops, starting with the residues located between positions 40 and 45, corresponding to the standard sequence of a CTX composed of 60 amino acids. The residues from 46 to 49, located at the end of loop III, form a type I turn in the β-structure of the protein. Continuing the loop structure, residues from 49 to 54 associate with those from 20 to 26 to establish an antiparallel β-strand. This configuration is of major interest as it plays a key role in the stabilization and functionality of the protein.

Moving to the cellular impact of cytotoxins, the mechanisms through which they induce cellular responses are diverse and complex. Notably, CTXs have been shown to trigger apoptosis through multiple pathways. For example, CTX3 from *Naja atra* venom activates apoptosis in human myeloid leukemia U937 cells through a cascade involving Ca²⁺, protein phosphatase 2A (PP2A), and AMP-activated protein kinase (AMPK), demonstrating the intricacy of intracellular signaling pathways affected by these toxins. Moreover, CTX1 from the same venom species further illustrates this complexity by inducing the transcription of genes associated with the Fas cell death pathway via a sequence of interactions starting from Ca²⁺ influx, followed by activation of NOX4, and resulting in increased ROS production that activates the p38 MAPK pathway. This detailed cascade of events underscores the cytotoxins' potential as targets for therapeutic interventions, particularly in the context of diseases characterized by deregulated apoptosis mechanisms [4,8]. Currently, there are no cytotoxins from cobra venom that have progressed to clinical trials or are available on the market specifically for cancer treatment. However, studies continue to explore the anticancer potential of these cytotoxins.

Research on the application of cytotoxins for glioblastoma treatment is notably sparse, with few studies available in the current literature. Given the highly aggressive nature of glioblastoma, there is a critical need for novel therapeutic approaches. To address this gap, we isolated six distinct peptides from the venoms of two *Naja* species and systematically investigated their antiproliferative activities. Our preliminary results provide initial insights into the potential efficacy of these venom-derived peptides against glioblastoma.
## 2. Results and Discussion

We screened nearly 40 different animal venoms to evaluate their potential anticancer properties. The venoms were solubilized at a concentration of 1.5 mg/mL in water. The antiproliferative activity of these venoms was tested at two concentrations, 3 µg and 15 µg, against MCF-7 cells (Appendix A, Figure A1). From this group, two highly active venoms, previously undescribed in the literature, *Naja anchietae* and *Naja senegalensis*, were selected and subjected to further fractionation and subfractionation using RP-HPLC to isolate and identify the active components responsible for the observed biological activities.

### 2.1. Physicochemical Characterization

#### 2.1.1. Purification of CTNs

The venoms of *Naja anchietae* and *Naja senegalensis* were subjected to purification via reverse-phase high-performance liquid chromatography (RP-HPLC). The initial gradient applied was from 0% to 70% solvent B over 70 min, equivalent to a gradient slope of 1% solvent B per minute, utilizing a C18 column with a 100 Å pore size (Figure 1).

<figure><img src="image_5.png" /><figcaption>**Figure 1.** Preparative HPLC chromatogram of *Naja anchietae* venom using gradient 0–70% for 70 min at 214 nm a C18 100 Å column. Fractions were collected at intervals of every 0.1 min, corresponding to a flow rate of 6 mL per minute. The region containing the active peptides is highlighted with a red square in the figure. These peptides were eluted between 40 and 52 min, which corresponds to an elution percentage of approximately 30%.</figcaption></figure>

Three peptides were isolated from this region (Figure 1), designated as CTNanc1, CTNanc2, and CTNanc3. These peptides were characterized using analytical RP-HPLC under a gradient of 20–50% solvent B over 60 min at a flow rate of 1 mL per minute. Detection was carried out at wavelengths of 214 nm and 280 nm. The profile of the homogeneous products is presented in the Figure 2.

As illustrated in Figure 2, the peptides CTNanc2 and CTNanc3 exhibit very similar hydrophobic properties, with retention times differing by only 0.5 min. This similarity has adversely affected the yield of these two peptides. Interestingly, despite the peptide CTNanc3 appearing highly homogeneous in the analytical RP-HPLC profile as shown in the figure, it presented a complex and contaminated mass spectrum (Figure 3).

Consequently, the purification protocol was adjusted to employ a more gradual gradient of 20% to 50% solvent B over 120 min, corresponding to a gradient slope of 0.25% solvent B per minute. This adjustment was implemented on a Phenomenex Jupiter 300 Å,
5 µm C18 column (250 mm × 21.2 mm), featuring larger pore sizes (Figure 4). This modification led to significantly improved resolution, thereby enhancing the homogeneity and yield of all peptides isolated from these venoms. Indeed, it is believed that materials with wider pores more effectively facilitate the entry of larger proteins into the porous matrix and enhance the partitioning between phases [9].

This new strategy enabled the successful purification of the three peptides, with their purity confirmed by mass spectrometry, as demonstrated in the spectra presented in the following section.

Similarly, the venom of *Naja senegalensis* was purified according to this protocol, as depicted in the Figures 5 and 6.

<figure><img src="image_6.png" /><figcaption>**Figure 2.** Analytical HPLC chromatogram of the isolated peptides from *Naja anchietae* at 214 nm on a C18 100 Å column.</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>**Figure 3.** Mass spectrum of CTNanc3 using MALDI-ToF.</figcaption></figure>
<figure><img src="image_8.png" /><figcaption>**Figure 4.** Preparative HPLC chromatogram of *Naja anchietae* venom using gradient 20–60% for 120 min at 214 nm on a C18 300 Å column. Fractions were collected at intervals of every 0.1 min, corresponding to a flow rate of 6 mL per minute. The region containing the active peptides is highlighted with a red square in the figure. These peptides were eluted between 80 and 120 min, which corresponds to an elution percentage of approximately 37%.</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>**Figure 5.** Preparative HPLC chromatogram of *Naja senegalensis* venom using gradient 20–60% for 120 min at 214 nm on a C18 300 Å column. Fractions were collected at intervals of every 0.1 min, corresponding to a flow rate of 6 mL per minute. The region containing the active peptides is highlighted with a red square in the figure. These peptides were eluted between 90 and 120 min, which corresponds to an elution percentage of approximately 40%.</figcaption></figure>
<figure><img src="image_10.png" /><figcaption>**Figure 6.** Analytical HPLC chromatogram of the isolated peptides from *Naja senegalensis* at 214 nm on a C18 100 Å column.</figcaption></figure>

### 2.1.2. MALDI-TOF Mass Spectrometry

Active peaks isolated and purified from the venoms of *Naja anchietae* and *Naja senegalensis* were characterized using analytical high-performance liquid chromatography (HPLC). Homogeneous molecules identified via HPLC were subsequently analyzed by matrix-assisted laser desorption/ionization time-of-flight (MALDI-TOF) mass spectrometry. This analysis revealed that the peaks designated as CTNsen1 (Aa), CTNsen2 (Ab), CTNsen3 (Ac), CTNanc1 (Bd), CTNanc2 (Be), and CTNanc3 (Bf) had molecular weights of [M + H] 6766.4 Da, 6849.6 Da, 6840.8 Da, 6688.7 Da, 6850.4 Da, and 6862.9 Da, respectively, as shown in Figure 7.

<figure><img src="image_11.png" /><figcaption>**Figure 7.** *Cont.*</figcaption></figure>
<figure><img src="image_12.png" /><figcaption>**Figure 7.** Mass spectrum of peptides CTNsen1 (Aa), CTNsen2 (Ab), CTNsen3 (Ac), CTNanc1 (Bd), CTNanc2 (Be), and CTNanc3 (Bf) using MALDI-ToF.</figcaption></figure>

### 2.1.3. N-Terminal Edman Sequencing

Peptides from *Naja anchietae* and *Naja senegalensis* were subjected to reduction with dithiothreitol (DTT) and alkylation with 4-vinylpyridine (4VP), as detailed in the Materials and Methods section. Subsequently, these peptides underwent N-terminal Edman sequencing. To identify homologous sequences, a BLAST search was conducted to determine the best alignment scores. The results of the sequencing and the alignments are summarized in the table below.

All overall repetitive yields exceed 90%, confirming the efficiency and precision of the sequencing performed. Beyond the 50th cycle, the detected PTH amounts decreased due to carryover effects.

The missing portion of CTNsen1 was determined using MALDI TOF-TOF MS/MS following digestion with the AspN enzyme on the reduced and alkylated peptide. This analysis successfully identified the following C-terminal sequence: 'DV(PEC)PKNSALVKYV(PEC) (PEC)NT(D)', as shown in Table 1.

CTNsen1, CTNanc1, and CTNanc2 demonstrated a 100% sequence identity with previously described cytotoxins, indicating that these are well-conserved within the toxin literature. However, CTNsen2, CTNsen3, and CTNanc3 exhibited 98.33% similarity to known cytotoxins, as determined by BLAST analysis. The minor sequence variations identified suggest that these cytotoxins may represent novel variants. These distinctions could potentially influence their biological functions and interactions.

Interestingly, BLAST alignment revealed that the sequences of CTNsen3 and CTNanc3 are both highly similar to the sequence of Cytotoxin 2 from *Naja nivea*. These sequences were compared in Figure 8 to highlight the subtle yet potentially significant variations that distinguish them from the known Cytotoxin 2 sequence.

This cytotoxin family is distinguished by four conserved disulfide bridges that enhance structural stability. The locations of these crucial disulfide linkages are consistent across the family, specifically at positions 3–21, 14–38, 42–53, and 54–59 in the amino acid sequence, underscoring the highly conserved nature of these cytotoxic proteins.
<table><caption>**Table 1.** Determined sequences by EDMAN sequencing and their homologous sequences by BLAST.</caption><thead><tr><th>Peptide</th><th>Determined Sequence</th><th>Homologous Sequence<br>Identified Using<br>BLAST/UniProtKB/Swiss-<br>Prot-% of Identity</th></tr></thead><tbody><tr><td>CTNsen1</td><td>LKCHQLVPPFWKTCPEGKNLCYKMYMVSSSTVPVKRGCIDVCPKNS<br>ALVKYVCCNTDKCN</td><td>P01457.1<br>Cytotoxin 5<br>Naja haje haje 100%</td></tr><tr><td>CTNsen2</td><td>LKCHKLVPPFWKTCPEGKNLCYKMYMVATPMIPVKRGCIDVCPKNS<br>ALVKYVCCNTNKCN</td><td>P01464.1<br>Cytotoxin 5 Naja haje<br>annulifera 98.33%</td></tr><tr><td>CTNsen3</td><td>LKCHQLVPPFWKTCPEGKNLCYKMYMVATPMIPVKRGCIDVCPKN<br>SALVKYMCCNTDKCN</td><td>P01463.1 Cytotoxin 2 Naja<br>nivea 98.33%</td></tr><tr><td>CTNanc1</td><td>LKCHKLVPPVWKTCPEGKNLCYKMFMVSTSTVPVKRGCIDVCPKDS<br>ALVKYVCCSTDKN</td><td>P01456.1 Cytotoxin 1 Naja<br>nivea 100.00%</td></tr><tr><td>CTNanc2</td><td>LKCHKLVPPFWKTCPEGKNLCYKMYMVATPMLPVKRGCIDVCPKD<br>SALVKYMCCNTDK</td><td>P01462.1<br>Cytotoxin 2 Naja annulifera<br>100.00%</td></tr><tr><td>CTNanc3</td><td>LKCHKLIPPFWKTCPEGKNLCYKMYMVATPMIPVKRGCIDVCPKDS<br>ALVKYMCCNTDKCN</td><td>P01463.1 Cytotoxin 2 Naja<br>nivea 98.33%</td></tr></tbody></table>

<figure><img src="image_13.png" /><figcaption>**Figure 8.** Comparative analysis of CTNsen3, CTNanc3, and Cytotoxin 2 from *Naja nivea* sequences. The differences are in bold and cysteines forming disulfide bridges are highlighted in yellow.</figcaption></figure>

In-depth studies, such as those by Menez et al., have pinpointed key structural features of CTXs, including three critical disulfide bridges named B1 (C1–C3), B2 (C2–C4), and B4 (C7–C8). These bridges, along with adjacent residues, create essential atomic contacts that are crucial for the protein's functionality. The substitution of cysteine residues with alanines in proteins, resulting in the removal of a specific disulfide bridge, has shown that the elimination of bridges B3 (C5–C6) or B4 drastically alters the three-finger structure of the protein. Conversely, the removal of bridges B1 or B2 does not significantly impact the biological function of the proteins. This discovery underscores the variable importance of disulfide bridges in maintaining the three-dimensional structure of three-finger toxins. While bridges B3 and B4 are essential for the structural integrity required for cytotoxin function, bridges B1 and B2 play a less critical role. Their absence does not prevent the cytotoxin from maintaining a functional conformation [7,10].

Furthermore, the positioning of the Asn60 residue near the terminal cysteine (C8), contributing to a CN motif, is crucial for the conformation of the toxins. This residue, consistently found in this key position across all CTXs, plays a vital role as its side chain forms three hydrogen bonds, thus stabilizing the hydrophobic center of the molecule. Additionally, the transformation of aspartate 57 to asparagine (mutation D57N) has been identified as a disruptive factor for CTX structures in neutral environments. This change weakens interactions within the localized hydrophobic cluster at the C- and N-termini of the molecule, destabilizing the salt bridge between Asp57 and Lys2 [7].

Based on our studies results, CTNsen1 and CTNanc1 are classified as S-type cytotoxins, characterized by their specific structural and functional properties. Conversely, CTNsen2,
CTNsen3, CTNanc2, and CTNanc3 are categorized as Type P cytotoxins. Type S cytotoxins present a contrast in their mode of action compared to Type P cytotoxins. While Type P cytotoxins enhance their binding and disruptive capabilities on cellular membranes through an additional binding site associated with Pro 31 (or 30), Type S cytotoxins, characterized by the Ser 29 (or 28) residue, tend to engage in more hydrophilic interactions. This fundamental structural difference influences how these toxins affect cellular membranes: Type P cytotoxins are adept at deeper penetration into the lipid bilayer, enhancing their membrane destabilization effects as observed in the studies conducted by Suzuki-Matsubara et al. and Dubovskii et al., leading to necrosis or apoptosis [11,12]. In contrast, Type S cytotoxins primarily interact through electrostatic attractions facilitated by their charged lysine rich loops and are often studied for their ability to induce programmed cell death (apoptosis) rather than immediate lysis. This results in a milder impact on membrane structures, with interactions largely confined to the surface rather than extensive penetration. This mode of action does not lead to the extensive disruption typical of Type P, reflecting a more nuanced approach to cellular interaction that could have different biological implications [5,13]. Some of the examples of S and P type cytotoxins from *Naja* species are listed in the Table 2 below.

<table><caption>**Table 2.** Some examples of studies of Type S and Type P CTXs derived from *Naja* species against cancer.</caption><thead><tr><th>Cytotoxin</th><th>Type</th><th>Active Against/Dose</th><th>Venom Source</th><th>Reference</th></tr></thead><tbody><tr><td>CTX-2N</td><td>Type S</td><td>0.8 µM for A549 non-small cell lung adenocarcinoma</td><td><i>Naja nigricollis</i></td><td>Conlon et al., 2020 [14]</td></tr><tr><td>NK-CTX/ Cytotoxin 3</td><td>Type S</td><td>Lung cancer (A549) A549 = 1.22 µg/mL<br>Prostate cancer PC-3 = 4.46 µg/m<br>Breast cancer MCF-7 = 12.23 µg/mL</td><td><i>Naja kaouthia</i></td><td>Chong et al., 2020 [2]</td></tr><tr><td>CTX-II</td><td>Type P</td><td>Breast cancer (MCF-7)</td><td><i>Naja oxiana</i></td><td>Ebrahim et al., 2014 [15]</td></tr><tr><td>Cytotoxin 2a/ NS-CTX</td><td>Type P</td><td>Lung cancer (A549) 0.88 µg/mL<br>Prostate cancer PC-3 = 3.13 µg/m<br>Breast cancer MCF-7 = 9.10 µg/mL</td><td><i>Naja sumatrana</i></td><td>Chong et al., 2020 [2]</td></tr><tr><td>Cytotoxin 1</td><td>Type S</td><td>Acute myeloid leukemia (KG-1a) 3.31 µg/mL</td><td><i>Naja atra</i></td><td>Liu et al. (2019) [16]</td></tr><tr><td>Cytotoxin 3/ NN-32</td><td>Type S</td><td>Breast cancer MCF-7 = 2.5 µg/mL</td><td><i>Naja naja</i></td><td>Attarde and Pandit (2017) [17]</td></tr><tr><td>CT1No</td><td>Type S</td><td>Leukemia WEHI-3 = 340 nM<br>Myelogenous leukemia K562 = 650 nM<br>Acute promyelocytic leukemia HL-60 = 1400 nM</td><td><i>Naja oxiana</i></td><td>Feofanov et al. (2004) [18]</td></tr><tr><td>CT2Nh</td><td>Type P</td><td>Leukemia WEHI-3 = 1190 nM<br>Acute promyelocytic leukemia HL-60 = 3700 nM</td><td><i>Naja haje</i></td><td>Feofanov et al. (2004) [18]</td></tr><tr><td>Cytotoxin 3</td><td>Type P</td><td>Human neuroblastoma SK-N-SH = 0.8 µM<br>Breast adenocarcinoma MDA-MB-231 = 0.09–0.15 mM</td><td><i>Naja atra</i></td><td>Chen et al. (2008) [19], Lin et al. (2010) [20]</td></tr><tr><td>Cytotoxin 4</td><td>Type S</td><td>Human neuroblastoma SK-N-SH = 1.2 µM</td><td><i>Naja atra</i></td><td>Chen et al. (2008) [19]</td></tr></tbody></table>

In the scientific literature, numerous cytotoxins derived from snake venoms have been described for their efficacy against various types of cancer and malignancies, such as leukemia, breast cancers, and lung cancers. However, the exploration of snake venom cytotoxins in the context of glioblastoma treatment remains scarce or even non-existent.

An alternative classification of cytotoxins introduced in 2024 by Dubovskii and Utkin provides a new approach to understanding their interactions with cellular membranes. This new classification system categorizes cytotoxins based on the presence of key amino
acids that significantly influence their membrane activity. According to this scheme, cytotoxins such as CTNsen1, CTNsen2, CTNanc1, CTNanc2, and CTNanc3, which contain the specific residues Pro9, Ser28, and Asn/Asp45, are assigned to Group 1. This group is characterized by a lower membrane disruptive activity, indicating a nuanced mechanism of action where these cytotoxins bind to cell membranes but are less likely to form disruptive pores (1). This result does not align well with our flow cytometry (FACS) findings for the peptides CTNanc1 and CTNanc2 which demonstrated a substantial increase in Propidium iodide (PI)-positive cells. This increase is indicative of late apoptosis or necrosis, suggesting that apoptosis in these cases progresses to a later stage where there is a loss of membrane integrity.

## 2.2. Pharmacological Characterization

Previous research has indicated strong antiproliferative activity of *Naja* venoms in other cancer types, such as MCF-7 [5], which motivated our study in the specific context of glioblastoma. Moreover, cytotoxicity tests to evaluate antiproliferative effects are generally simpler and more straightforward than those required for other tests. Furthermore, in glioblastoma and other aggressive cancers, the rapid proliferation of tumor cells is a major characteristic that contributes to the severity of the disease. Thus, targeting proliferation can have an immediate and measurable therapeutic impact, crucial for developing effective treatments.

### 2.2.1. In Vitro Cytotoxicity Assay

The purified peptides, named CTNsen1, CTNsen2, CTNsen3, CTNanc1, CTNanc2, and CTNanc3 from *Naja senegalensis* and *Naja anchietae* respectively, were subsequently tested for their antiproliferative effects on several glioblastoma cell lines (U87, U251, and T98G) as well as on primary human umbilical vein endothelial cells (HUVECs). A series of increasing concentrations (1 × 10⁻⁷, 0.5 × 10⁻⁶, 1 × 10⁻⁶, 0.5 × 10⁻⁵, 1 × 10⁻⁵, 0.5 × 10⁻⁴, and 1 × 10⁻⁴ Molar) were used. The dose-response curves for each peptide are depicted in Figures 9–12. Antiproliferative assays conducted on U87, U251, and T98G glioblastoma cell lines revealed distinct activities for the peptides derived from *Naja senegalensis* and *Naja anchietae*.

<figure><img src="image_14.png" /><figcaption>**Figure 9.** Antiproliferative activity of CTNsen1 (a), CTNsen2 (b), CTNsen3 (c), CTNanc1 (d), CTNanc2 (e), and CTNanc3 (f) from *Naja senegalensis* and *Naja anchietae*, respectively. These results were obtained using an MTT assay on U87 cells after 72 h of treatment with a series of increasing concentrations (1 × 10⁻⁷, 0.5 × 10⁻⁶, 1 × 10⁻⁶, 0.5 × 10⁻⁵, 1 × 10⁻⁵, 0.5 × 10⁻⁴, and 1 × 10⁻⁴ Molar) of peptides. The OD was measured at 570 nm.</figcaption></figure>
<figure><img src="image_15.png" /><figcaption>**Figure 10.** Antiproliferative activity of CTNsen1 (a), CTNsen2 (b), CTNsen3 (c), CTNanc1 (d), CTNanc2 (e), and CTNanc3 (f) from *Naja senegalensis* and *Naja anchietae* respectively. These results were obtained using an MTT assay on U251 cells after 72 h of treatment with a series of increasing concentrations (1 × 10⁻⁷, 0.5 × 10⁻⁶, 1 × 10⁻⁵, 0.5 × 10⁻⁵, 1 × 10⁻⁵, 0.5 × 10⁻⁴, and 1 × 10⁻⁴ Molar) of peptides. The OD was measured at 570 nm.</figcaption></figure>

<figure><img src="image_16.png" /><figcaption>**Figure 11.** Antiproliferative activity of CTNsen1 (a), CTNsen2 (b), CTNsen3 (c), CTNanc1 (d), CTNanc2 (e), and CTNanc3 (f) from *Naja senegalensis* and *Naja anchietae*, respectively. These results were obtained using an MTT assay on T98G cells after 72 h of treatment with a series of increasing concentrations (1 × 10⁻⁷, 0.5 × 10⁻⁶, 1 × 10⁻⁵, 0.5 × 10⁻⁵, 1 × 10⁻⁵, 0.5 × 10⁻⁴, and 1 × 10⁻⁴ Molar) of peptides. The OD was measured at 570 nm.</figcaption></figure>

To determine whether the six identified peptides exhibit cytotoxicity on non-tumor cells, we chose to test their effects on HUVEC cells, as illustrated in Figure 12.

The IC50 values obtained from these results were calculated and are presented in Table 3.
<figure><img src="image_17.png" /><figcaption>**Figure 12.** Cytotoxic activity of CTNsen1 (a), CTNsen2 (b), CTNsen3 (c), CTNanc1 (d), CTNanc2 (e), and CTNanc3 (f) from *Naja senegalensis* and *Naja anchietae*, respectively. These results were obtained using an MTT assay on HUVEC cells after 72 h of treatment with a series of increasing concentrations (1 × 10⁻⁷, 0.5 × 10⁻⁶, 1 × 10⁻⁵, 0.5 × 10⁻⁵, 1 × 10⁻⁵, 0.5 × 10⁻⁴, and 1 × 10⁻⁴ Molar) of peptides. The OD was measured at 570 nm.</figcaption></figure>

Our results revealed notable antiproliferative activity of CTNanc3 against the U87 glioblastoma cell line, exhibiting an IC50 of 5.33 µg/mL. However, this peptide also demonstrated substantial toxicity towards HUVEC endothelial cells with a comparable IC50, suggesting that its effects might largely be attributed to general cytotoxicity rather than selective anticancer activity. CTNanc3 displayed minimal activity against U251 and T98G glioblastoma cell lines, with IC50 values of 416.37 µg/mL and 420.01 µg/mL, respectively, indicating poor efficacy in these models. In contrast, CTNanc1 emerged as the most potent peptide against U87 cells, with an IC50 of 36.41 µg/mL and no cytotoxicity against HUVEC endothelial cells up to a concentration of 1 × 10⁻⁴ M (686 µg). Interestingly, this activity appears to be U87-selective since CTNanc1 showed little or no activity towards U251 and T98G cell lines even at high concentrations. CTNanc2 displayed considerable activity against the U87 cell line with an IC50 of 38.34 µg/mL. However, it also proved cytotoxic towards HUVEC cells, suggesting a toxicity profile that could limit its therapeutic application despite its efficacy (Figures 9–12, Table 3).

Concerning the CTNsen peptides on U87 cell lines: (i) CTNsen1 was active with an IC50 of 54.69 µg/mL, but showed notable cytotoxicity against HUVEC cells; (ii) CTNsen2 showed no antiproliferative activity on U87 cells even at high concentrations, nor cytotoxicity on normal cells; and (iii) CTNsen3 displayed interesting antiproliferative activity against U87, with an IC50 of 45.38 µg/mL, and importantly, no cytotoxicity was observed
on HUVEC cells at any concentration tested, underscoring its potential safety for non-cancerous tissues. Lastly, considering the cytotoxicity of all peptides on normal cells, no notable antiproliferative activity on the U251 and T98G lines was observed (Figures 9–12, Table 3). Overall, CTNanc1 and CTNsen3 stand out as the most effective and selective peptides against U87 cells; they exhibit good activity without toxicity to non-cancerous cells even at high concentrations.

<table><caption>**Table 3.** Differential cytotoxicity of CTN's in cancerous and non-cancerous cell lines. The half maximal inhibitory concentration (IC50) was estimated using a three-parameter log-logistic fit along with 95% confidence intervals limits (lower, higher). * Non-cytotoxic activity at $1 \times 10^{-4}$ Molar.</caption><thead><tr><th rowspan="2">Tissue Origin</th><th rowspan="2">Cell Line</th><th colspan="6">IC50 (µg/mL)</th></tr><tr><th>CTNanc1</th><th>CTNanc2</th><th>CTNanc3</th><th>CTNsen1</th><th>CTNsen2</th><th>CTNsen3</th></tr></thead><tbody><tr><td rowspan="2">Brain Glioblastoma<br>95% confidence interval</td><td>U87</td><td>36.41<br>32.87–39.96</td><td>38.34<br>34.63–42.04</td><td>5.33<br>4.65–6.02</td><td>54.69<br>46.11–63.28</td><td>Inactive</td><td>45.38<br>36.80–53.97</td></tr><tr><td>U251</td><td>Inactive</td><td>475.8<br>395.04–556.56</td><td>416.37<br>338.15–494.59</td><td>Inactive</td><td>Inactive</td><td>Inactive</td></tr><tr><td rowspan="2">Brain Glioblastoma<br>95% confidence interval</td><td>T98G</td><td>Inactive</td><td>389.59<br>296–483.18</td><td>420.01<br>384.96–455.06</td><td>Inactive</td><td>462.18<br>386.35–538.01</td><td>Inactive</td></tr><tr><td>HUVEC</td><td>Inactive*</td><td>467.21<br>423.44–510.97</td><td>3.58<br>3.42–3.74</td><td>225.65<br>167.03–284.28</td><td>Inactive*</td><td>Inactive*</td></tr></tbody></table>

The variable sensitivity of different glioblastoma cell lines to these peptides suggests that the molecular environment of each cell line, including the expression of specific receptors, membrane composition, and intracellular pathways, plays a critical role in determining their cytotoxic efficacy. Proteomic analysis of U87 and U251 glioblastoma cell lines reveals distinct differences in their protein expression profiles, reflecting their unique cellular characteristics and responses to environmental conditions. U87 cells, known for robust growth and strong adherence in culture, exhibit higher levels of proteins associated with cellular metabolism, cytoskeletal structure, and stress response pathways. This suggests that U87 cells may have adapted to in vitro conditions by optimizing pathways that support their rapid proliferation and structural integrity. In contrast, U251 and T98G cells display a proteomic profile characterized by elevated expression of proteins involved in DNA repair, transcription regulation, and apoptosis, indicating these cells may prioritize genomic stability and controlled cell death mechanisms to maintain their survival in culture. Additionally, U251 cells show a reduction in certain proteins abundantly expressed in U87 cells, highlighting different regulatory mechanisms and potential vulnerabilities [21,22].

Cytotoxins that showed a high percentage of identity with ours via BLAST display varied anticancer activities according to the literature. Cytotoxin 5 from *Naja haje haje*, as reported by Lafnoune et al. in 2021, comprises about 23% of the proteins in the F7 fraction of the venom and exhibits potent anticancer effects. It significantly reduces the size of hepatocellular carcinoma multicellular tumor spheroids, thus impacting cancer cell proliferation directly. This study primarily characterized the anticancer activities of this fraction without further separation into individual components [23]. Conversely, no anticancer activity has been described for Cytotoxin 5 from *Naja haje annulifera* or for Cytotoxin 1 and Cytotoxin 2 from *Naja annulifera*. However, Cytotoxin 2 from *Naja nivea* shows promising results. According to Ebrahim et al. in 2014, this peptide induces apoptosis in the MCF-7 breast adenocarcinoma cell line, demonstrating a significant decrease in cell viability with an IC50 of approximately 4.18 µg/mL, far lower than the standard drug cisplatin. The peptide initiates a cascade of apoptosis-related intracellular events, including the generation of reactive oxygen species (ROS), disruption of mitochondrial membrane potential, and activation of caspase-9 and caspase-3. This indicates a strong potential for therapeutic application, although Cytotoxin 2 was also tested against the non-tumorigenic MCF 10A cell line with an IC50 of 18.12 ± 2.56 µg/mL (around $26.4 \times 10^{-4}$ M), suggesting
its broader cytotoxicity profile [15]. Importantly, no studies have described any positive anticancer activity of these peptides against glioblastoma cells in the literature.

To mitigate the harmful effects of cytotoxins on healthy tissue while still effectively targeting glioblastoma cells, several strategies can be employed. One approach involves the use of specialized delivery systems such as liposomes or nanoparticles, which are designed to direct treatment specifically to cancer cells, thereby reducing unwanted effects on other cells. These carriers can be functionalized with antibodies or ligands that exclusively bind to markers overexpressed on glioblastoma cells, ensuring that the cytotoxin is delivered primarily to the tumor site. For instance, lipid crystal nanoparticles and proliposomes encapsulating temozolomide, a well-established medication for glioblastoma treatment, have been engineered to enhance penetration of the blood–brain barrier, facilitating targeted delivery of this chemotherapeutic agent directly to brain tumor cells, as described by Waghule et al. in 2021 [24]. Another strategy is the modification of cytotoxins through PEGylation or conjugation. For example, conjugated toxins like ANG1005, which consists of paclitaxel conjugated with angiopep, target the LRP-1 receptor that is overexpressed in various cancers, including glioblastoma. This conjugation enhances the drug’s ability to traverse the blood–brain barrier and more effectively target tumor cells, as detailed by Kumthekar et al. in 2020 [25]. Lastly, combination therapy serves as a third strategy. Combining cytotoxins with targeted inhibitors or chemotherapy drugs can synergistically increase cancer cell death and mitigate toxicity to normal cells by lowering the required doses of each agent. This strategy not only amplifies the therapeutic impact but also minimizes the chances of developing resistance to single-agent treatments. For instance, combining bevacizumab, an anti-VEGF antibody, with lomustine is employed in treating recurrent glioblastoma to disrupt the tumor’s blood supply concurrently with a direct cytotoxic attack, enhancing overall tumor control [26].

Finally, to determine the influence of disulfide bridges on the antiproliferative activity against U87 cells, the peptides CTNanc1 and CTNsen1 were reduced and alkylated before being tested on U87 cells (Figure 13). The disruption of disulfide bonds resulted in a loss of activity compared to the native oxidized peptides, underscoring the importance of these bonds for their antiproliferative activity. Additionally, a comparison of the primary structures of the studied cytotoxins reveals very slight variations in amino acids. Notably, CTNanc1, the most active peptide against the proliferation of U87 cell lines, contains a valine at position 10, unlike the phenylalanine found in other cytotoxins. It also replaces a tyrosine with a phenylalanine at position 25 and features an STSTV sequence from positions 28 to 32. These combined observations suggest that antiproliferative activity depends not only on the three-dimensional structure stabilized by disulfide bonds but also on the potential exposure of certain amino acids critical to this activity.

<figure><img src="image_18.png" /><figcaption>**Figure 13.** Antiproliferative activity of reduced and alkylated CTNsen1 (a) and CTNanc1 (b). These results were obtained using an MTT assay on U87 cells after 72 h of treatment with a series of increasing concentrations (1 × 10⁻⁷, 0.5 × 10⁻⁶, 1 × 10⁻⁵, 0.5 × 10⁻⁵, 1 × 10⁻⁵, 0.5 × 10⁻⁴, and 1 × 10⁻⁴ Molar) of peptides. The OD was measured at 570 nm.</figcaption></figure>
2.2.2. Flow Cytometry Analysis Results

a. Gating strategy

The first part of the gating strategy focuses on identifying and isolating viable cells from the total cell population (Figure 14A). This process begins by excluding doublets and aggregates through setting a gate on forward scatter height (FSC-H) versus forward scatter width (FSC-W), ensuring that only individual cells are analyzed. Subsequently, a second gate is applied based on forward scatter area (FSC-A) against side scatter area (SSC-A) to differentiate actual cells from debris and other non-cellular elements. Within this refined population, cells are then analyzed using two gating strategies, as shown in the representative examples for the control in Figure 14B,C.

<figure><img src="image_19.png" /><figcaption>**Figure 14.** Flow cytometry gating strategy: representative example using untreated U87 cells. Doublets and debris are excluded (A); the first gating strategy identifies living cells (PI-negative and Annexin V-negative), apoptotic cells (Annexin V-positive and PI-negative), and dead cells (Annexin V-positive and PI-positive) (B); the second gating assesses the intensity of Annexin V fluorescence and the percentage of Annexin V-positive cells, Annexin V-negative and positive cell subsets are colored in pink and orange, respectively (C).</figcaption></figure>

Cells are analyzed for Annexin V FITC/PI (Figure 14A), which defines the living cells (PI negative and Annexin V negative), the apoptotic cells (Annexin V positive and PI negative), and the dead cells (Annexin V positive and PI positive). The percentages of each cell subset from this analysis are reported in Table 4. The main objective of this gating strategy was to analyze U87 cell death after peptide treatment ([Dead cells] % gated, Table 4). To refine the analysis of the Annexin V-positive population among living cells, the Mean Fluorescence Intensity (MFI) of Annexin V-FITC was recorded (Figure 14B). PI-negative cells are indicative of membrane integrity and, therefore, viability (Figure 14B). Annexin V staining binds to phosphatidylserine on the external membrane of apoptotic cells. The extent of apoptosis within the viable cell subset correlates with the Annexin V fluorescence intensity (MFI). In the histogram from the control, a threshold of positivity was arbitrarily set at 50% to separate the Annexin V-negative and positive cell subsets. The percentage of Annexin V-positive gated cells and the Annexin V fluorescence intensity from the PI-negative (viable) cells are reported in Table 5. The comparison of the MFI and percentage of gated Annexin V-positive cells to the control is correlated with an increase in the apoptosis rate.
<table><caption>**Table 4.** Data from the initial gating strategy show cellular responses after 72 h of treatment with cytotoxin peptides at a concentration of 1 × 10⁻⁴ Molar. The percentages of living and dead cells reflect the effects of various cytotoxin peptides on U87 cell samples.</caption><thead><tr><th>Sample Tested</th><th>% Living Cells Gated</th><th>% Apoptotic Cells Gated</th><th>% Dead Cells Gated</th></tr></thead><tbody><tr><td>Control U87 cells</td><td>81</td><td>2</td><td>14</td></tr><tr><td>CTNanc1-U87</td><td>7</td><td>2</td><td>90</td></tr><tr><td>CTNanc2-U87</td><td>11</td><td>5</td><td>84</td></tr><tr><td>CTNanc3-U87</td><td>83</td><td>3</td><td>12</td></tr><tr><td>CTNsen2-U87</td><td>33</td><td>17</td><td>43</td></tr><tr><td>CTNsen3-U87</td><td>85</td><td>2</td><td>10</td></tr></tbody></table>

<table><caption>**Table 5.** Data from the second part of the gating strategy focus on the analysis of apoptotic markers in PI-negative U87 cells treated with cytotoxin peptides for 72 h at a concentration of 1 × 10⁻⁴ Molar. The Mean Fluorescence Intensity (MFI) of Annexin V-FITC and the percentage of Annexin V-positive cells, derived from histograms, indicate the extent of apoptosis in each treatment group.</caption><thead><tr><th>Sample Tested</th><th>[Cells PI Negative Cells] Number</th><th>[Cells PI Negative Cells] Annexin FITC MFI</th><th>[AnnexinV<sup>+</sup>] % Gated from Histograms</th></tr></thead><tbody><tr><td>Control U87 cells</td><td>71,033</td><td>29,764</td><td>50</td></tr><tr><td>CTNanc1-U87</td><td>312</td><td>59,862</td><td>76</td></tr><tr><td>CTNanc2-U87</td><td>2737</td><td>164,227</td><td>93</td></tr><tr><td>CTNanc3-U87</td><td>68,626</td><td>46,465</td><td>69</td></tr><tr><td>CTNsen2-U87</td><td>3689</td><td>200,997</td><td>96</td></tr><tr><td>CTNsen3-U87</td><td>65,486</td><td>29,041</td><td>49</td></tr></tbody></table>

Cells exhibiting Annexin V positivity are identified as early apoptotic since they are undergoing apoptosis but have not yet lost their membrane integrity to the extent that would permit PI staining.

b. Cytotoxicity analysis by flow cytometry

The Annexin V-FITC/PI staining analysis (Figure 15) revealed that the peptides CT-Nanc1 and CTNanc2 exhibited high cytotoxicity, notably reducing the viability of U87 cells. Specifically, cells treated with CTNanc1 showed a viability of only 7%, while those treated with CTNanc2 had a viability of 11% (Table 4). In addition, compared to the control, a higher percentage of Annexin V positive cells and Mean Fluorescence Intensity (MFI) are shown for CTNanc1 and CTNanc2 in Table 5, suggesting apoptosis. In contrast, the peptides CTNanc3 and CTNsen3 demonstrated much lower cytotoxicity, with cell viability percentages of 83% and 85%, respectively, indicating that they are less toxic to U87 cells. While the MTT assay indicated activity against U87 cells at the tested concentration (1 × 10⁻⁴ Molar), flow cytometry did not show pronounced cytotoxic effects in terms of apoptosis or necrosis. This discrepancy could be due to several factors: (i) differential sensitivity: the MTT assay measures metabolic activity, specifically mitochondrial activity, which can decrease due to factors other than cell death, such as metabolic reprogramming, cell cycle arrest, or temporary metabolic suppression. (ii) Flow cytometry directly measures cell death by quantifying apoptotic and necrotic markers. CTNanc3 and CTNsen3 might induce a form of cell death that was not detectable by the flow cytometry assay used, such as a slower or more regulated process that does not manifest within the same timeframe as the more aggressive cytotoxic actions observed with other CTNs. These toxins might also
exert effects that modulate cell behavior or viability without leading to outright cell death, impacting processes like migration, invasion, or proliferation.

<figure><img src="image_20.png" /><figcaption>**Figure 15.** Dot plots analysis of U87 cell responses to cytotoxin peptides over 72 h by flow cytometry. Cells were labeled with Annexin-V and Propidium iodide (PI). Cells that did not undergo treatment served as the negative control, Annexin V-negative and positive cell subsets are colored in pink and orange, respectively. Annexin V-negative and positive cell subsets are colored in pink and orange, respectively on histograms.</figcaption></figure>
Our hypothesis is that: (i) these CTNanc3 and CTNsen3 might be triggering sub-lethal stress responses that temporarily impair mitochondrial function or overall cellular metabolism without necessarily leading to cell death. Such stress could activate survival pathways, allowing cells to maintain viability even in the face of cytotoxic challenges. (ii) The peptides might induce autophagy, a process that can either be a survival mechanism or lead to cell death. Autophagy can reduce mitochondrial activity as part of a broader cellular recycling process, which could explain the reduced MTT signal.

Notably, CTNsen2 also exhibited considerable cytotoxicity, reducing cell viability to 33%, though it was not as extreme as CTNanc1 and CTNanc2 when compared with the control untreated cells (Table 4 and Figure 14). However, treatment of U87 cells with CTNsen2 showed the highest Annexin V Mean Fluorescence Intensity (MFI) and 96% of the cells were Annexin V-positive (Table 5), compared to untreated cells. This suggests that CTNsen2 induces cytotoxicity leading to the apoptotic stage after 72 h.

Upon further analysis, the rise in Annexin V-positive cells within the PI-negative population (Annexin V+/PI−) suggests that early apoptosis is the primary mechanism of cell death induced by these cytotoxins, particularly for CTNanc1 and CTNanc2, which showed 76% and 93% early apoptotic cells, respectively (Table 5). These percentages reflect early apoptosis, as they indicate cells that have externalized phosphatidylserine (Annexin V positivity) but still retain membrane integrity (PI-negative).

The rise in PI-positive cells, indicative of late-stage apoptosis or necrosis, was not directly measured in the specific columns of Table 5, but the overall increase in cytotoxicity suggests that apoptosis progresses to later stages in treatments with CTNanc1 and CTNanc2.

Apoptosis and necrosis represent two distinct modalities of cellular demise [27]. Historically, necrosis was perceived as a spontaneous and unregulated form of cell death, characterized by the swelling of organelles and fragmentation of the nucleus. However, recent research has increasingly demonstrated that necrosis can also occur through a controlled process, despite being independent of caspase activation [28]. This regulated form of necrosis, often referred to as programmed necrosis or necroptosis, involves specific molecular pathways that orchestrate cellular disintegration.

In 2019, Lui et al. reported that CTX1 from *Naja atra* induces necroptosis in leukemia cells. The team initially observed a reduction in cell viability following CTX1 treatment, which was not mitigated by Z-VAD-fmk, a caspase inhibitor, suggesting a caspase-independent pathway. Crucially, pretreatment with Necrostatin-1 (Nec-1), a specific inhibitor of necroptosis, significantly restored cell viability, indicating that the mode of cell death induced by CTX1 was necroptosis. Further supporting these findings was the absence of classic apoptotic features such as chromatin condensation and nuclear fragmentation. Instead, the researchers documented early loss of plasma membrane integrity and cellular swelling, which are characteristic signs of necroptosis [16]. CTX1 (P60304 · 3SA1_NAJAT) is classified as an S-type cytotoxin, similar to CTNanc1, suggesting that the significant increase in PI-positive U87 cells may be attributed to necroptosis. To further explore this hypothesis, planning experiments using Necrostatin-1 (Nec-1), a specific inhibitor of necroptosis, would be advisable.

**3. Conclusions**

This study represents preliminary research into the cytotoxic effects of novel peptides isolated from the venoms of *Naja anchietae* and *Naja senegalensis*, marking a significant advancement in our understanding of venom cytotoxins. We successfully isolated six cytotoxins that have not previously been described in these species. Each peptide was rigorously characterized both physicochemically and pharmacologically to ascertain their potential as therapeutic agents. Two cytotoxin peptides, CTNanc1 and CTNsen3, demonstrated antiproliferative activity on U87 cell lines without exhibiting notable cytotoxicity against human umbilical vein endothelial cells (HUVECs). Our findings are particularly groundbreaking as they introduce cytotoxins as potential therapeutic agents against glioblastoma, a context in which venom-derived cytotoxins have not been extensively studied. To our knowledge,
this is among the first studies to explore and document the efficacy of snake cytotoxins
from these specific snakes against glioblastoma cell lines.

**4. Materials and Methods**

**4.1. RP-HPLC Purification**

Solvents were purchased from Sigma-Aldrich (Saint-Quentin-Fallavier, France) and crude venoms of *Naja senegalensis* and *Naja anchietae* were supplied by Latoxan (Portes lès Valence, France). To purify our molecules, venoms solubilized in water at a concentration of 50 mg/mL were filtered through 0.45-micron filters. The peptides were purified from the filtrates using RP HPLC Agilent 1260 Infinity (Agilent Technologies, Santa Clara, CA, USA) on RP column Phenomenex Jupiter 300A 5 µm C18 250 mm × 21.2 mm (Phenomenex Inc., Torrance, CA, USA). The gradient used for this purification is 20% to 50% over 120 min of solvent B, which consisted of 0.1% TFA in 90% acetonitrile, and solvent A, containing 0.1% TFA in water. The fractions collected every 0.1 min were analyzed using a Shimadzu LC-2010 HPLC system (Shimadzu, Kyoto, Japan) on RP column MACHEREY-NAGEL NUCLEOSIL 100-5 C18, 5 µm, 250 × 4.6 mm (Düren, Germany).

**4.2. Mass Spectrometry Analysis**

The isolated *Naja* venom peptides were characterized by spectrometry MALDI-ToF using an α-cyanohydrocinnamic acid (HCCA) matrix (Bruker Daltonics, Wissembourg, France) (0.7–1 µL equal volumes of saturated solution at 10 µg/µL in MeCN or in 50% MeCN/0.3% TFA/H₂O). Spectra recorded in linear positive method LP12 kDa (Bruker microflex) or HMASS 5–20 kDa were externally calibrated with suitable standards as Prot cal I of Brucker (Bruker Daltonics, Wissembourg, France) with a laser (λ = 337 nm, W = 75.2 µJ, 500 shots of ν = 10 Hz), and were analyzed by a Bruker Daltonics flex analysis software Flex Control 1.2 SP1 version 3.0 (2008) (Bruker, Palaiseau, France).

**4.3. Reduction and Alkylation of Peptides**

In total, 20 nanomoles of each peptide sample were dissolved in 160 uL of alkylation buffer containing 6 M guanidine hydrochloride (GuHCl), 0.5 M Tris/HCl, 2 mM EDTA (pH 7.5), and 4.8 µM dithiothreitol. Following 20 h of reduction in the dark, at 40 degrees, 6 µM of 4-vinylpyridine (4-VP) was added to initiate the alkylation process. The proteins that were S-alkylated were then purified using a reverse-phase C8 column (5 µm; 4.6 × 250 mm) on a RP HPLC Agilent 1260 Infinity system (Agilent Technologies, USA).

**4.4. N-Terminal Sequencing**

N-terminal sequencing of S-alkylated peptides was conducted through Edman degradation. Samples were deposited onto a glass fiber membrane GFD that had been pre-treated with 10 µL of biobrene and conditioned using the PRETREATGFD program. A calibration was performed using a solution of phenylthiohydantoin (PTH)-amino acid standard at a concentration of 10 pmol. The membrane was then inserted into the cartridge of the Shimadzu PPSQ 31B sequencer (Kyoto, Japan), and analysis was performed. The obtained sequences were subsequently cross-referenced with entries in the Uniprot database.

**4.5. Peptide Mass Fingerprintings**

Peptides were digested by AspN enzyme (Promega, Madison, WI, USA) according to manufacturer recommendations. In total, 10 µL of the sample at 7 mg/mL in H₂O (1 µg/µL) was treated with 10 µL of 0.1 M DTT in Ammonium Bicarbonate 0.1 M for 20 min at 56 °C, then alkylated by 10 µL Iodoacetamide 10 mg/mL in Ammonium Bicarbonate 0.1 M. Then, 2 µg AspN was resuspended in 50 µL H₂O, and 5 µL (200 ng) (1/50 *w*/*w*) was added to the reaction middle and incubated overnight at 37 °C. Peptides were then resuspended in a α-cyano-4-hydroxycinnamic acid matrix solution (saturated solution in 50% acetonitrile/0.3% TFA), and 1 µL was spotted on the Maldi metal target. Mass analyses were performed on a MALDI-TOF-TOF Ultraflex III spectrometer (Bruker Dal-
tonics, Wissembourg, France) controlled by the Flexcontrol 3.0 package (Build 51). This instrument was used at a maximum accelerating potential of 25 kV and was operated in reflector mode and the *m/z* range from 600 to 3500. Six external standards (Peptide Calibration Standard II, Bruker Daltonics) were used to calibrate each spectrum to a mass accuracy within 10 ppm with a minimal resolution of 10,000 for the angiotensin II peak (Monoisotopic mass = 1,046,542 Da). Peak picking was performed with Flexanalysis 3.0 software (Bruker) with an adapted analysis method. Parameters used were as follows: SNAP peak detection algorithm, S/N threshold fixed to 6, and a quality factor threshold of 30.

### 4.6. MS/MS Analysis

The sequence of P1V214 peptide was confirmed by MS/MS analyses on the Maldi TOF-TOF Ultraflex III spectrometer (Bruker Daltonics, Wissembourg, France) using the LIFT_2017 method after calibration. LIFT mass spectra were acquired on a Bruker Ultraflex TOF/TOF mass spectrometer (Coventry, UK) operated in the positive ion mode. Metastable fragmentation was induced by a nitrogen laser (337 nm) without the further use of collision gas. Precursor ions were accelerated to 8 kV and selected in a timed ion gate. In the LIFT-cell the fragments were further accelerated to 19.07 kV. The reflector potential was 29.6 kV [29].

### 4.7. Cell Culture

U87 glioma cells and human umbilical vein endothelial cells (HUVECs) were cultured under controlled laboratory conditions. The U87, T98G and U251 cells were grown in Dulbecco's Modified Eagle Medium (DMEM) supplied by Gibco™ (Sigma, St. Louis, MO, USA), supplemented with 10% fetal bovine serum (FBS) and 1% L-glutamine. HUVECs, procured from Lonza (Paris, France), were maintained in EBM-2 medium from the same supplier. This growth medium was fortified with hydrocortisone (1 µg/mL), bovine brain extract (12 µg/mL), epidermal growth factor (10 ng/mL), and 2% FBS (Life Technologies, Carlsbad, CA, USA). These cells were kept at 37 °C in a humidified environment containing air and 5% CO₂. HUVECs were regularly tested and confirmed to be free from mycoplasma contamination and cultured up to the fifth passage to ensure viability and consistency in experiments.

### 4.8. Cell Viability Assay

U87 and T98G were plated at a density of 4000 cells per well and U251 cells at a density of 2500 cells per well in microtiter plates and allowed to adhere overnight at 37 °C in a 5% CO₂ humidified environment. The following day, the medium was replaced with fresh medium containing our venom peptides under investigation. After a 72-h incubation period, cells were treated with MTT (Promega, Lyon, France) at a concentration of 0.5 mg/mL. The MTT crystals, reduced by the cells' mitochondrial dehydrogenases, were then solubilized using DMSO. Viable cells were then quantified by measuring the absorbance at 570 nm. HUVEC cells were plated at a density of 10,000 cells per well in microtiter plates and tested under the same conditions. Concerning the statistical analysis, modeling the differential cytotoxicity of CTNs in cancerous and non-cancerous cell lines was achieved using a log-logistic fit. Dose-response analyses were performed using the drc package [30], available from R software (R Core Team, 2023) [31]. The best results were obtained through three-parameter log-logistic fitting models after a Box-Cox transformation for each peptide/cell line combination. Model calibration was assessed using the drm() function. The half-maximal inhibitory concentration (IC50) and 95% confidence interval limits were obtained from the ED() function.

### 4.9. Flow Cytometry Analysis

U87 cells underwent flow cytometric analysis following staining with Annexin V-FITC and Propidium Iodide (PI). These cells were evaluated using a CytoFLEX S (Beckman
Coulter, Brea, CA, USA) flow cytometer. For the experimental setup, cells were seeded into 6-well plates and treated with the different peptides, each at a concentration of 10⁻⁴ M. After a 72-h incubation period, cells were trypsinized and harvested, ensuring collection of both suspended and adherent cells. For the staining process, 5 µL of Annexin-V-fluorescein isothiocyanate (FITC) and 10 µL of propidium iodide (PI) were added to approximately 1 × 10⁵ cells. The mixture was incubated for 15 min at room temperature, shielded from light. Subsequently, the cells were immediately subjected to flow cytometric analysis. The categorization of cells was based on their staining patterns: viable (Annexin V−/PI−), early apoptotic (Annexin V+/PI−), late apoptotic or secondary necrosis (Annexin V+/PI+), and necrotic (Annexin V−/PI+).

**Author Contributions:** Y.B.: investigation, methodology, writing, review, editing, and purification of CTX, pharmacological testing, and physicochemical characterizations; C.B.-D.: technical supervision for pharmacological tests; M.B.: flow cytometry experiments; H.M.: participation in venom screening on MCF7; M.M.: discussions on flow cytometry experiments, reading and correction of the manuscript; P.M.: mass spectrometry and N-terminal sequencing; P.F.: peptide mass fingerprinting and MS/MS analysis; F.T.: Statistical analysis and determination of appropriate confidence intervals for all data; H.d.P., D.G. and L.O.: conceptualization; K.M.: conceptualization and design of the experiments, writing, review, editing, and final approval. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research received no external funding.

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** The results of this study are fully documented and presented within the confines of this article.

**Acknowledgments:** We sincerely thank Eric Nebel (Agilent and Advion Interchim Scientific) for his valuable advice and continuous support in HPLC separations and column selection.

**Conflicts of Interest:** The authors declare no conflicts of interest.

## Appendix A

<figure><img src="image_21.png" /><figcaption>**Figure A1.** Antiproliferative activity of different animal venoms. These results were obtained using an MTT assay on MCF7 cells after 72 h of treatment with 3 and 15 µg/mL. The OD was measured at 570 nm after 72 h of treatment. The abbreviations include Na: *Naja anchietae*; Vg: *Viridovipera gumprechti*; Vk: *Vipera kaznakovi*; Tv: *Trimeresurus vogeli*; Tt: *Trimeresurus trigonocephalus*; Ach: *Atheris chlorechis*; Ath ni: *Atheris nitschei nitschei*; Nox: *Naja oxiana*; Np: *Naja philippinensis*; At sq: *Atheris squamigera*; Nhl: *Naja haje legionis*; Por Oph: *Porthidium ophryomegas*; Ph: *Parias hageni*; Ns: *Naja senegalensis*.</figcaption></figure>

Through meticulous HPLC fractionation, we were able to isolate six distinct peptides, named CTNsen1, CTNsen2, CTNsen3, CTNanc1, CTNanc2, and CTNanc3, from *Naja senegalensis* and *Naja anchietae*, respectively. Their activity was confirmed on MCF7 cells as shown in Figure A2.
<figure><img src="image_22.png" /><figcaption>**Figure A2.** Antiproliferative activity of CTNsen1, CTNsen2, CTNsen3, CTNanc1, CTNanc2, and CTNanc3 from *Naja senegalensis* (b) and *Naja anchietae* (a), respectively. These results were obtained using an MTT assay on MCF7 cells after 72 h of treatment at 4 µg and 20 µg of peptide. Na and Ns here were used as a positive control at two concentrations (3 µg and 15 µg). The OD was measured at 570 nm.</figcaption></figure>

## References

1. Dubovskii, P.V.; Utkin, Y.N. Specific Amino Acid Residues in the Three Loops of Snake Cytotoxins Determine Their Membrane Activity and Provide a Rationale for a New Classification of These Toxins. *Toxins* **2024**, *16*, 262. [CrossRef] [PubMed]

2. Chong, H.P.; Tan, K.Y.; Tan, C.H. Cytotoxicity of snake venoms and cytotoxins from two Southeast Asian cobras (*Naja sumatrana*, *Naja kaouthia*): Exploration of anticancer potential, selectivity, and cell death mechanism. *Front. Mol. Biosci.* **2020**, *7*, 58358. [CrossRef]

3. Kalita, B.; Utkin, Y.N.; Mukherjee, A.K. Current insights in the mechanisms of cobra venom cytotoxins and their complexes in inducing toxicity: Implications in antivenom therapy. *Toxins* **2022**, *14*, 839. [CrossRef] [PubMed]

4. Hiu, J.J.; Yap, M.K.K. The myth of cobra venom cytotoxin: More than just direct cytolytic actions. *Toxicon* **X** **2022**, *14*, 100123. [CrossRef] [PubMed]

5. Dubovskii, P.V.; Utkin, Y.N. Antiproliferative activity of cobra venom cytotoxins. *Curr. Top. Med. Chem.* **2015**, *15*, 638–648. [CrossRef]

6. Kumar, T.K.S.; Jayaraman, G.; Lee, C.S.; Arunkumar, A.I.; Sivaraman, T.; Samuel, D.; Yu, C. Snake venom cardiotoxins-structure, dynamics, function and folding. *J. Biomol. Struct. Dyn.* **1997**, *15*, 431–463. [CrossRef]

7. GKonshina, A.; VDubovskii, P.; GEfremov, R. Structure and dynamics of cardiotoxins. *Curr. Protein Pept. Sci.* **2012**, *13*, 570–584. [CrossRef]

8. Chiou, J.T.; Wang, L.J.; Lee, Y.C.; Chang, L.S. *Naja atra* cardiotoxin 1 induces the FasL/Fas death pathway in human leukemia cells. *Cells* **2021**, *10*, 2073. [CrossRef] [PubMed]

9. Boughanmi, Y.; Aboudou, S.; Boyadjian, K.; Charouandi, A.; Bouzid, S.; Barnetche, T.; Mabrouk, K. One-Step Chromatographic Approach for Purifying Peptides and Proteins from Venoms. *Separations* **2024**, *11*, 179. [CrossRef]

10. Galat, A.; Gross, G.; Drevet, P.; Sato, A.; Ménez, A. Conserved structural determinants in three-fingered protein domains. *FEBS J.* **2008**, *275*, 3207–3225. [CrossRef]

11. Suzuki-Matsubara, M.; Athauda, S.B.; Suzuki, Y.; Matsubara, K.; Moriyama, A. Comparison of the primary structures, cytotoxici- ties, and affinities to phospholipids of five kinds of cytotoxins from the venom of Indian cobra, *Naja naja*. *Comp. Biochem. Physiol. Part C Toxicol. Pharmacol.* **2016**, *179*, 158–164. [CrossRef] [PubMed]

12. Dubovskii, P.V.; Lesovoy, D.M.; Dubinnyi, M.A.; Konshina, A.G.; Utkin, Y.N.; Efremov, R.G.; Arseniev, A.S. Interaction of three-finger toxins with phospholipid membranes: Comparison of S-and P-type cytotoxins. *Biochem. J.* **2005**, *387*, 807–815. [CrossRef]

13. Efremov, R.G.; Volynsky, P.E.; Nolde, D.E.; Dubovskii, P.V.; Arseniev, A.S. Interaction of cardiotoxins with membranes: A molecular modeling study. *Biophys. J.* **2002**, *83*, 144–153. [CrossRef] [PubMed]
14. Conlon, J.M.; Attoub, S.; Musale, V.; Leprince, J.; Casewell, N.R.; Sanz, L.; Calvete, J.J. Isolation and characterization of cytotoxic and insulin-releasing components from the venom of the black-necked spitting cobra *Naja nigricollis* (Elapidae). *Toxicon X* **2020**, *6*, 100030. [CrossRef] [PubMed]

15. Ebrahim, K.; Shirazi, F.H.; Vatanpour, H.; Kobarfard, F.; Rabiei, H. Anticancer activity of cobra venom polypeptide, cytotoxin-II, against human breast adenocarcinoma cell line (MCF-7) via the induction of apoptosis. *J. Breast Cancer* **2014**, *17*, 314–322. [CrossRef]

16. Liu, Y.; Ming, W.; Wang, Y.; Liu, S.; Qiu, Y.; Xiang, Y.; Guo, Q. Cytotoxin 1 from *Naja atra* Cantor venom induced necroptosis of leukemia cells. *Toxicon* **2019**, *165*, 110–115. [CrossRef]

17. Attarde, S.S.; Pandit, S.V. Cytotoxic activity of NN-32 toxin from Indian spectacled cobra venom on human breast cancer cell lines. *BMC Complement. Altern. Med.* **2017**, *17*, 1–9. [CrossRef]

18. Feofanov, A.V.; Sharonov, G.V.; Dubinnyi, M.A.; Astapova, M.V.; Kudelina, I.A.; Dubovskii, P.V.; Rodionov, D.I.; Utkin, Y.N.; Arseniev, A.S. Comparative study of structure and activity of cytotoxins from venom of the cobras *Naja oxiana*, *Naja kaouthia*, and *Naja haje*. *Biochemistry* **2004**, *69*, 1148–1157. [CrossRef]

19. Chen, K.C.; Chiou, Y.L.; Kao, P.H.; Lin, S.R.; Chang, L.S. Taiwan cobra cardiotoxins induce apoptotic death of human neuroblastoma SK-N-SH cells mediated by reactive oxygen species generation and mitochondrial depolarization. *Toxicon* **2008**, *51*, 624–634. [CrossRef]

20. Lin, K.L.; Su, J.C.; Chien, C.M.; Chuang, P.W.; Chang, L.S.; Lin, S.R. Down-regulation of the JAK2/PI3K-mediated signaling activation is involved in Taiwan cobra cardiotoxin III-induced apoptosis of human breast MDA-MB-231 cancer cells. *Toxicon* **2010**, *55*, 1263–1273. [CrossRef]

21. Vogel, T.W.; Zhuang, Z.; Li, J.; Okamoto, H.; Furuta, M.; Lee, Y.S.; Weil, R.J. Proteins and protein pattern differences between glioma cell lines and glioblastoma multiforme. *Clin. Cancer Res.* **2005**, *11*, 3624–3632. [CrossRef] [PubMed]

22. Patil, V.; Pal, J.; Somasundaram, K. Elucidating the cancer-specific genetic alteration spectrum of glioblastoma derived cell lines from whole exome and RNA sequencing. *Oncotarget* **2015**, *6*, 43452. [CrossRef] [PubMed]

23. Lafnoune, A.; Lee, S.Y.; Heo, J.Y.; Gourja, I.; Darkaoui, B.; Abdelkafi-Koubaa, Z.; Oukkache, N. Anti-cancer effect of moroccan cobra *Naja haje* venom and its fractions against hepatocellular carcinoma in 3d cell culture. *Toxins* **2021**, *13*, 402. [CrossRef] [PubMed]

24. Waghule, T.; Rapalli, V.K.; Singhvi, G.; Gorantla, S.; Khosa, A.; Dubey, S.K.; Saha, R.N. Design of temozolomide-loaded proliposomes and lipid crystal nanoparticles with industrial feasible approaches: Comparative assessment of drug loading, entrapment efficiency, and stability at plasma pH. *J. Liposome Res.* **2020**, *31*, 158–168. [CrossRef]

25. Kumthekar, P.; Tang, S.C.; Brenner, A.J.; Kesari, S.; Piccioni, D.E.; Anders, C.; Ibrahim, N. ANG1005, a brain-penetrating peptide–drug conjugate, shows activity in patients with breast cancer with leptomeningeal carcinomatosis and recurrent brain metastases. *Clin. Cancer Res.* **2020**, *26*, 2789–2799. [CrossRef]

26. Taal, W.; Oosterkamp, H.M.; Walenkamp, A.M.; Dubbink, H.J.; Beerepoot, L.V.; Hanse, M.C.; van den Bent, M.J. Single-agent bevacizumab or lomustine versus a combination of bevacizumab plus lomustine in patients with recurrent glioblastoma (BELOB trial): A randomised controlled phase 2 trial. *Lancet Oncol.* **2014**, *15*, 943–953. [CrossRef]

27. Chen, D.; Yu, J.; Zhang, L. Necroptosis: An alternative cell death program defending against cancer. *Biochim. Biophys. Acta (BBA)-Rev. Cancer* **2016**, *1865*, 228–236. [CrossRef]

28. Yuan, J.; Najafov, A.; Py, B.F. Roles of caspases in necrotic cell death. *Cell* **2016**, *167*, 1693–1704. [CrossRef]

29. Suckau, D.; Resemann, A.; Schuerenberg, M.; Hufnagel, P.; Franzen, J.; Holle, A. A novel MALDI LIFT-TOF/TOF mass spectrometer for proteomics. *Anal. Bioanal. Chem.* **2003**, *376*, 952–965. [CrossRef]

30. Ritz, C.; Baty, F.; Streibig, J.C.; Gerhard, D. Dose-Response Analysis Using R. *PLoS ONE* **2015**, *10*, e0146021. [CrossRef]

31. R Core Team. *R: A Language and Environment for Statistical Computing R*; Foundation for Statistical Computing: Vienna, Austria, 2023. Available online: https://www.R-project.org/ (accessed on 31 August 2024).

**Disclaimer/Publisher's Note:** The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.