1. Introduction
Bioactive peptides (BAPs) are a heterogeneous group of molecules, generally between 5-50 aminoacids with a broad-spectrum of activities, expressed throughout the tree of life, from bacteria to humans [
1]. One interesting source of BAPs if the snake venoms, they are mixtures of a great many distinct proteins and peptides that display an amazing diversity of structure and function [
2]. NA-CATH:ATRA1-ATRA1 (NA) is a cationic (+15) and amphipathic peptide build of 34 amino acids (KRFKKFFKKLKNSVKKRAKKFFKKPKVIGVTFPF). This peptide is a designed synthetic sequence from the modification of the NA-CATH cathelicidin identified in the venom of the elapid snake
Naja atra, and the substitution of the motif ATRA1 (KRFKKFFKKLK) in the structure [
3]. NA has demonstrated to present a potent antimicrobial activity (EC
50 = 0.51 μg/mL or 0.09 μM) against
Staphylococcus aures (
S. aureus) cell cultures and biofilms [
3].
BAPs are mainly known for their direct interactions with microbial membranes, but in the last decade their potential activity against cancer cells has make them interesting pharmaceutical source [
4]. Their most accepted mechanism of action is based on their ability to alter the membranes of target cells through complex molecular interactions between the peptide and the lipid environment of the microbial or cancer membrane [
5]. BAPs mechanism involves several stages, beginning with electrostatic interactions between the positively charged residues of the peptide and the negatively charged lipids of the target membrane. In the case of cancer cells, the presence of phosphatidylserine (PS) in the outer membrane enhances the electrostatic interactions with BAPs [
6,
7], in contrast to normal cell membranes, which are predominantly neutral due to their zwitterionic lipid composition [
8,
9]. Following this, hydrophobic interactions between the phospholipid acyl chains and non-polar residues of the peptide facilitate the integration of the peptide into the membrane bilayer through various modes, such as barrel-stave, carpet detergent, and toroidal pore formations, leading to cell death [
10].
The activity of BAPS makes them a promising avenue for breast cancer treatment, as they can help mitigate resistance mechanisms developed by cancer cells [
11]. While traditional chemotherapeutics must penetrate cancer cells to be effective, often leading to resistance, bioactive peptides (BCPs) have the advantage of acting from outside the membrane—a mechanism that cancer cells struggle to counteract [
12,
13]. Previous results of our group demonstrated that NA presented antitumoral activity against two molecular subtypes of breast cancer, the MCF-7 (IC
50 = 13.4 µM) and the MDA-MB-23 (IC
50 = 6.4 µM) cells [
14]. However, there is no information reported about the biophysical interaction of NA and lipid models representative of cancer cells. The study of the interactions responsible of the biological activity are important to understand and design new therapeutic molecules. For this reason, in this study we investigated the potential antitumoral mechanism of NA against two representative lipid models of MCF-7 and the MDA-MB-23 cell membranes. Additionally, a comparative study was performed with a non-tumoral model system representative of immortalized human keratinocytes (HaCaT) cells. The quantification of the lipid composition of the three cell lines was previously reported by Klaiss-Luna and co-workers [
15]. The lipid systems were employed to investigate the lipid–peptide interaction using differential scanning calorimetry (DSC) and infrared spectroscopy (FT-IR). This study was focused on the thermodynamic properties, structure, and dynamics of both tumor and non-tumor model membranes. Additionally, the study examined the secondary structure–activity relationship of the peptide in specific aqueous and lipid environments. The findings offer valuable insights into how the lipid profile affects membrane–peptide interactions and help elucidate the antitumoral mechanism of action of NA.
3. Discussion
Bioactive peptides are emerging as prospective pharmaceutical sources to treat several types of human diseases. The BAP of interest for this research was the NA-CATH-ATRA1-ATRA1 (NA), KRFKKFFKKLKNSVKKRFKKFFKKLKVIGVTFPF, a synthetic peptide designed and modified from the Chinese cobra
Naja atra host defense peptide NA-CATH, KRFKKFFKKLKNSVKKRAKKFFKKPKVIGVTFPF, by substituting alanine to phenylalanine at position 18 and proline to leucine at position 25, with the aim to increase its hydrophobic character and improve the antimicrobial activity already reported from the parental peptide against
Staphylococcus aureus [
3]. Dean and collaborators reported that NA evidenced greater antimicrobial and antibiofilm potency against the pathogen
S. aureus and associated it to a possible peptide structure factor, especially in the presence of anionic membrane environment [
3]. Since BAPs commonly report a broad spectrum activity, specially a concomitant antimicrobial-antitumoral activity [
17], our research group has been interested in evaluating the antitumoral activity of NA-CATH-ATRA1-ATRA1. Previous results concerning the half maximal inhibitory concentration in breast cancer cell lines MCF-7, MDA-MB-231 evidenced a promising antitumoral activity [
14], however how it exerts its activity is still unknown. The present paper contributes to study the interaction of NA and lipid systems at the membrane level and explores its potential antitumoral molecular mechanism through the biophysical techniques: Differential Scanning Calorimetry and Infrared Spectroscopy using lipid models of biological cell membranes in established conditions.
The outcome of this work is discussed in this section based on the physicochemical properties of the peptide and the lipid composition of model membranes. The latter emphasizes in the presence of phosphatidylserine in the cancer model membranes to resemble the negatively charged environment at membrane level caused by the overexpression of PS in the exoplasmic leaflet of tumoral membranes, and the higher presence of anionic molecules like heparin sulfates, sialylated gangliosides and O-glycosylated mucins [
18,
19]. DSC results indicate that model membranes of breast cancer cells were significantly affected by the presence of NA, the above was expected due to the strong electrostatic interaction between the +15 cationic peptide and the negatively charged vesicles. However, how the interactions occurred seemed to be in different ways considering the effects on the thermodynamic parameters reported in
Section 2.1.1. We suggest that in the vesicles, which are model membranes of MCF-7, NA generates peptide-poor lipid domains that might reflect a partial insertion of the peptide between fatty acids in a concentration-dependent manner, while with MDA-MB-231 model it is proposed that the peptide induced the stabilization of the gel phase due to the increase in T
m. In addition, NA potentially contributes to diminishing the electrostatic repulsion from negatively charged phosphatidylserine in the MDA-MB-231 model, which reduces the distance between lipid molecules and improves the van der Waals interactions between acyl chains [
20]. On the opposite, NA evidenced a fluidizing effect in the HaCaT model membrane with no significant changes in enthalpy, leading to consider a superficial membrane interaction with lipid polar headgroups. Based on the described findings, NA-CATH-ATRA1-ATRA1 is able to alter the physical state of model membranes exerting an effect that is dependent not only on peptide concentration but also on the membrane lipid composition.
Since multicomponent lipid systems exhibited broad chain-melting endotherms, we decided to evaluate the individual interactions with PC, PE, PS and SM vesicles to complement the previous information. Our results indicate that NA strongly interacts with phosphatidylcholine headgroup at the highest evaluated concentrations, perturbing the physical properties of the lipid vesicles. In fact, there are reports in the literature showing the frequent binding of cationic molecules to PC headgroup by electric surface charge that change the direction of the dipole between N
+ and water molecules, affecting the orientation of phosphate group [
21], which explains the abolishment of the ripple phase transition in PC thermogram. The NA effect on PC might predict a cytotoxic effect since phosphatidylcholine is the most abundant phospholipid [
22], indeed, there are reports about the low cytotoxicity from NA in host cells [
3,
23,
24]. The outcome interaction with DPPE vesicles possibly implies a lipid surface aggregation and primary hydrophobic interaction, similar to the report by Willumeit and collaborators with NK-2 peptide [
25].
Although all phospholipids seem to be affected by the presence of the peptide, a closer look at the thermodynamical parameters from
Table 2 reveals a remarkable selective effect they have towards PS and SM lipid systems. Our findings demonstrated a strong and favorable electrostatic interaction between NA and phosphatidylserine that ended up destabilizing the vesicles. The peptide induced lipid microheterogeneity and phase separation to which we suggest NA could be intercalated between fatty acids chains perturbing the lipid packing by reducing the Van der Waals interactions.
Finally, since the phase transition of SM was abolished in the presence of the peptide, we believe it is due to the destabilization of lipid structure following NA binding to the sphingosine moiety through several intra- and intermolecular hydrogen bonds at the interfacial level.
Furthermore, the lipid-peptide interaction assays analyzed by FTIR-ATR are in concordance with DSC evidence, the NA and MCF-7 interaction reflects the strong binding that causes a release of water molecules and at the same time a relief of charge repulsion with a consequent increase in lipid packing [
26]. Spectra portrayed the generation of a more dehydrated and rigid MDA-MB-231 system in presence of NA where the lipid-peptide interaction is more favorable in the liquid-crystalline phase. Taking into consideration that membranes at physiological conditions are in the fluid state, it would be expected a remarkable potential antitumoral activity [
27,
28]. In addition, the peptide incubation with HaCaT SLBs highlights the fluidizing effect caused by their binding and the consecutive release of water molecules from the interfacial region [
29].
Most biological activities of peptides have been related to their structure. Accordingly, Dean and coworkers stated that the increased helical nature of NA generates greater antimicrobial activity [
3]. We observed that negatively charged environments promote the disordered to ordered transition of NA during incubation with PS liposomes until having a proper helical folding of a peptide. Although NA seems to generate a conformational change upon binding, the disordered segment should be considered essential to exert its potential antitumoral activity.
Finally, after assembling the described above we conclude and suggest a peptide localization and interaction, with the breast cancer model membranes MCF-7 and MDA-MB-231 after the initial electrostatic attraction peptide seems to reduce the electrostatic repulsion from PS in the models while it acquires a slightly helical conformation that might promote insertion in the lipid bilayer, having a stronger insertion in the MDA-MB-231 model. On the other hand, NA peptide in the control model HaCaT appears to be located on the lipid surface and exerting an electrostatic interaction with lipid headgroups, thus giving energetically unfavorable conditions for the peptide to fold. These biophysical results are in agreement with previous biological assays and offer insights into the potential antitumoral activity of NA-CATH-ATRA1-ATRA1, however, more studies using complementary techniques and biological assays are required to elucidate its mechanism of action.
4. Materials and Methods
4.1. Peptide and Phospholipids
The Naja atra derived peptide NA-CATH-ATRA1-ATRA1 (KRFKKFFKKLK-NSVKKRFKKFFKKLKVIGVTFPF, Lot. U037QFC180-11/PE6100) was acquired from GenScript (Piscataway Township, NJ, USA), synthesized by the solid-phase method. The 97.5 % purity and molecular weight of the peptide were determined using HPLC and MALDI–TOF mass spectrometry, respectively. 16:0 and 16:0/18:1 acyl chain length phospholipids were purchased from Avanti Polar Lipids (Alabaster, AL, USA). 1,2-Dipalmitoyl-sn-glycero-3-phosphocholine (DPPC, Lot. 160PC-318), 1-palmitoyl-2-oleoyl-sn-glycerol-3-phosphocholine (POPC, Lot. 850457P-500MG-A-211), 1-,2-dipalmitoyl-sn-glycero-3-phosphoethanolamine (DPPE, Lot. 160PE-106), 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine (POPE, Lot. 850757P-500MG-B-151), 1-,2-dipalmitoyl-sn-glycero-3-phospho-L-serine sodium salt (DPPS, Lot. 840037P-500MG-A-078), 1-palmitoyl-2-oleoyl-sn-glycerol-3-phospho-L-serine sodium salt (POPS, Lot. 840034P-25MG-397 A-250), sphingomyelin egg chicken (SM, Lot. 860061P-25MG-A-116). HEPES, EDTA, NaCl, and other analytical grade reagents were purchased from Sigma-Aldrich (St. Louis, MO, USA).
4.2. Phase Transition Measurements by DSC
Individual stock solutions of phospholipids were prepared at 10 mM, DPPC and SM in chloroform, while DPPE and DPPS in chloroform:methanol (70:30 v/v). 1 mM multilamellar vesicles (MLVs) were prepared by adding the exact volume from the individual stock solutions into a test tube to resemble the multi-component phospholipid composition reported for the following breast cancer model membranes, MCF-7, DPPC/DPPE/SM/DPPS 47:35:5:13 (w/w); MDA-MB-231, DPPC/DPPE/SM/DPPS 43:33:4:20 (w/w), and as non-tumoral model membrane HaCaT DPPC/DPPE/SM 59:35:6 (w/w) [
16]. Solvent was evaporated under N
2 stream, and the dried lipid film was hydrated with buffer 10 mM HEPES, 500 mM NaCl, and 1 mM EDTA pH 7.4. MLVs were formed after 6 min of several vortex-hot water immersion cycles at a temperature above the transition temperature of the lipid system and 1 min sonication. To evaluate the peptide effect on the MLVs, the lipid film was hydrated by adding a NA-CATH-ATRA1-ATRA1 solution prepared at 1, 5, and 10 molar% in the same buffer.
Microcalorimeter Nano DSC Series III platinum capillary cell (TA Instruments, New Castle, DE, USA) was used to perform the calorimetric experiments. The reference cell was loaded with 400 µL of buffer, while the sample cell with lipid or lipid-peptide suspension. After sealing, the calorimeter was pressurized at 0.3 MPa and cells were equilibrated for 10 min at the initial temperature. Heating scans were performed at 1 °C min−1 in the temperature intervals 15 °C to 60 °C for DPPC, 40 °C to 80 °C for DPPE, 35°C to 75 °C for DPPS, 15 °C to 50 °C for SM, 15 °C to 75 °C for breast cancer model membranes MCF-7 and MDA-MB-231, and, from 15 °C to 65 °C for the non-tumoral model membrane HaCaT. Main phase transition temperature (Tm) and main phase transition enthalpy (∆H) were analyzed and calculated using the NanoAnalyze 3.12.0 software (TA Instruments, New Castle, DE, USA). The accuracy was ± 0.1 °C for Tm, and ± 1 kJ.mol−1 for ∆H. DSC data was plotted using the software Origin Pro 8.0 (OriginLab Corporation, Northampton, MA, USA).
4.3. Infrared Spectroscopy Experiments
Phospholipid masses were weighed and added into a glass test tube to achieve a 20 mM stock solution of the breast cancer model membranes MCF-7, DPPC/DPPE/SM/DPPS 47:35:5:13 (w/w); MDA-MB-231, DPPC/DPPE/SM/DPPS 43:33:4:20 (w/w), and the non-tumoral model membrane HaCaT, DPPC/DPPE/SM 59:35:6 (w/w). The dried mixture was dissolved in chloroform and stored at -20 °C. For the gel to liquid crystalline phase measurements 20 mM Supported Lipid Bilayers (SLBs) were prepared in situ on the silicon crystal of the BioATR II cell coupled to the Tensor II spectrometer (Bruker Optics, Ettlingen, Germany) with an MCT (Mercury, Cadmium, Tellurium) detector and connected to a circulating water bath Huber Ministat 125 (Huber, Offenburg, Germany) that offers ± 0.1 °C temperature accuracy. SLBs were formed after adding 20 µL of the stock, removing the solvent, and hydrating the lipid film 15 min above the main transition temperature with the same volume of buffer 10 mM HEPES, 500 mM NaCl, and 1 mM EDTA pH 7.4 in the absence or presence of peptide at 1, 5, 10 molar%. Infrared spectra were taken at a temperature range of 43 to 63 °C for MCF-7, 45 to 65 °C for MDA-MB-231, and 40 to 60 °C for HaCaT. The final spectra were the average of 120 consecutive scans per temperature after the buffer subtraction with a resolution of 4 cm−1.
Lipid phase transition temperature was evaluated through the sensitive marker of lipid order, the symmetric stretching vibration of the methylene band (2970 to 2820 cm−1), and the interfacial hydration by the sensitive marker of lipid hydration, the ester carbonyl stretching vibration band (1725 to 1740 cm−1). Data analysis was performed using the OPUS 3D software (Bruker Optics, Ettlingen, Germany); vibrational bands of interest were cut from final spectra with subsequent baseline correction using the Rubberband method with 20% sensitivity. The pick-picking tool indicated the maximum wavenumber for the symmetric extensions per temperature, leading to graph wavenumber as a function of temperature. Tm was determined as the inflection point of the sigmoidal curve that was fitted into a Boltzmann function using the iterative algorithm Levenberg Marquardt by Origin Pro 8.0 software (Origin Lab Corporation, USA).
4.4. Peptide Conformational Analysis
3 mg. mL−1 peptide solution was prepared in buffer (10 mM HEPES, 500 mM NaCl, 1mM EDTA, pH 7.4) in the absence or presence of small unilamellar vesicles (SUVs). 6 mM SUVs were prepared following the composition of the model membranes MCF-7, POPC/POPE/SM/POPS 47:35:5:13 (w/w); and MDA-MB-231, POPC/POPE/SM/POPS 43:33:4:20 (w/w), and HaCaT POPC/POPE/SM 59:35:6 (w/w). Each lipid mass was weighed, mixed, and dissolved in chloroform. The solvent was removed with N2 stream, lipid films were hydrated with buffer and sonicated for 30 min, 50/60 Hz at a temperature above the main phase transition temperature of the lipid model. Peptide and vesicles were mixed to attain a 15 molar% peptide concentration, and its mixture was incubated for 5 min at 37 °C before measurements. The final spectra of the peptide were acquired as the average of 124 scans taken at 37 °C by the Tensor II spectrometer (Bruker Optics, Ettlingen, Germany) using the integrated AquaSpec Transmission Cell. The secondary structure elements alpha-helix, and beta-sheet with ± 4.4 % deviation were predicted by the predefined method BPROT1 from the Confocheck FT-IR system after internal calibrations with the data set of 43 proteins.