9S1R nullomer peptide induces mitochondrial pathology, metabolic suppression, and enhanced immune cell infiltration, in triple-negative breast cancer mouse model

Nullomers are the shortest strings of absent amino acid (aa) sequences in a species or group of species. Primes are those nullomers that have not been detected in the genome of any species. 9S1R is a 5-aa peptide prime sequence attached to 5-arginine aa, used to treat triple negative breast cancer (TNBC) in an in vivo mouse model. This unique peptide, administered with a trehalose carrier (9S1R-NulloPT), offers enhanced solubility and exhibits distinct anti-cancer effects against TNBC. In our study, we investigated the effect of 9S1R-NulloPT on tumor growth, metabolism, metastatic burden, tumor immune-microenvironment (TME), and transcriptome of aggressive mouse TNBC tumors. Notably, treated mice had smaller tumors in the initial phase of the treatment, as compared to untreated control, and diminished in vivo and ex vivo bioluminescence at later-stages - indicative of metabolically quiescent, dying tumors. The treatment also caused changes in TME with increased infiltration of immune cells and altered tumor transcriptome, with 365 upregulated genes and 710 downregulated genes. Consistent with in vitro data, downregulated genes were enriched in cellular metabolic processes (179), specifically mitochondrial TCA cycle/oxidative phosphorylation (44), and translation machinery/ribosome biogenesis (45). The upregulated genes were associated with the developmental (13), ECM organization (12) and focal adhesion pathways (7). In conclusion, our study demonstrates that 9S1R-NulloPT effectively reduced tumor growth during its initial phase, altering the TME and tumor transcriptome. The treatment induced mitochondrial pathology which led to a metabolic deceleration in tumors, aligning with in vitro observations.


Background
The chance of a woman being diagnosed with breast cancer (BC) during her lifetime is 1 in 8 [1].In the US, invasive breast cancer is predicted to be diagnosed in an estimated 297,790 women and 2800 men in 2023 [2] and by 2040, with no major changes in prevention or treatment, 1.4 million women will die from BC worldwide [3].There is a well-established heterogeneity in BC subtypes with presence, absence or a combination of estrogen receptor alpha (ERɑ), human epidermal growth factor receptor (HER2) and progesterone receptors (PR), that dictate treatment strategies.Based on SEER (Surveillance,Epidemiology,and End Results) [4] data in the United States 12% of breast tumors are triple-negative, which lack ER-, PR-and HER2-(triple-negative BC -TNBC) [5].This subtype is the most aggressive, with the highest level of invasiveness and worst prognosis, resulting in 25-46% of brain metastasis [4].BC treatment consists principally of surgery, radiation therapy, chemotherapy, hormonal therapy, targeted antibody or small-molecule therapy.Ground breaking approaches such as CAR-T cell therapy show only modest effects on solid tumors [6] and few drugs have been approved by the U.S. FDA for patients with TNBC tumors.PD-1 (programmed cell death protein-1) inhibitor in combination with chemotherapy is one such drug but is limited to patients with tumors expressing PD-L1 (programmed cell death-ligand 1), and those with high-risk, early stage TNBC [7].The only approved targeted therapy for TNBC is a Trop-2 targeted antibody-and topoisomerase I inhibitor conjugate, used as a second-line treatment for unresectable locally advanced or metastatic triplenegative breast cancer (mTNBC) [8,9].Recent cancer immunotherapy drug development is focused on immunological agents that augment the natural immune response of patients against cancer-specific antigens [10].While therapies involving adjuvant immunotherapies, oncolytic viruses, cytokines, antibodies, peptides, and their combinations are in clinical trials for BC treatment, there are certain drawbacks associated with each such as limited efficacy, lack of response, higher cost and poor accessibility for lower income nations.Peptide-based drugs are showing great promise in clinical studies [11][12][13].Peptides offer the advantage of small size, specificity, effects on a broad range of cancers, low toxicity and low manufacturing cost [13].In the current study, we investigated the therapeutic potential of a novel 10 amino acid nullomer-peptide (NulloP) in a TNBC mouse model.
Nullomers are the shortest absent sequences in a species (or group of species).These can be nucleotide strings, or amino acid strings (Nullopeps, or NulloPs).The set of nullomers absent from the entire biome (as represented in available databases) are called primes (DNA) or peptoprimes (amino acids) [14][15][16].Both nullomers and primes are seeds from which other absent sequences can be constructed by adding to either end of the string.The nullomer approach to drug design is based on a simple proposition that the smallest absent sequences can affect reactions when injected into hosts.There are two main hypotheses regarding the effects of injected nullomers on cancer cells, they may be toxic or immune stimulatory.The toxic hypothesis is based on the idea that evolution is the longest running experiment in biology, and selection acts to restrict or amplify the permutations of nucleotide and amino acid sequences that arise by chance.In the most extreme case, sequences that conflict with the basic common elements of cell physiology (ribosomes, mitochondria etc.) may be completely "forbidden" [14,17].The immunogenic hypothesis is based on the idea that the immune system identifies short sequences in a polymer as self or non-self.We could view the immune system as a nullomer search and destroy mission, thus nullomer-based drugs may increase immune responses.The original peptoprime sequences described in 2007 [14] are five amino acids long and have been modified with the addition of five arginine molecules to enhance cell-penetrating properties [17,18].When dissolved in trehalose, several of the peptoprimes (NulloPTs) were shown to be preferentially lethal to breast and prostate cancer cells, as opposed to normal primary cell lines [17][18][19].Trehalose is reported to have an independent effect against BC [20][21][22][23][24] and is an important component in our drug preparation.Nullomers have also been used to design vaccines from epitopes present in the host but absent in a pathogen, and to make DNA "watermarks" for labeling by combining them in series [19].
The original 198 peptoprimes (and scrambled versions of their sequences) were assayed in vitro against the NCI 60 cancer lines [18], and 9S1R, a potent killer of cancer cells was selected for the mouse experiments reported here.9S1R is a scrambled version of an original peptoprime, that is a nullomer (absent) from the biome.9S1R-NulloPT is absent from the human and mouse proteomes and has been shown to dramatically decrease ATP production, inhibit mitochondrial F o F 1 -ATP synthase, reduce mitochondrial membrane potential and increase superoxide free radicals inside the mitochondria of MDA-MB-231 human TNBC cells and the NCI 60 cancer cell panel, ultimately killing them [17,18].This is the first in vivo study of nullomer-based peptides.Here we report the therapeutic effects of 9S1R-NulloPT in the 4T1.2-Luc mouse TNBC model.The 9S1R-NulloPT was initially tested in a bilateral tumor model (pilot study-1, details in Supplementary) with six administrations at two different doses (50 mg/kg and 100 mg/kg).Following these results, we moved to a unilateral tumor model with eight drug administrations at the most effective dose (100 mg/kg).Although our preclinical in vitro studies showed a substantial response against diverse panels of cancer cells, the present in vivo study in mice showed that the drug effectively reduces TNBC tumor size at the initial growing phase of the tumor, but not in the advanced stage.9S1R-NulloPT changes the tumor immune microenvironment and affects tumor energy metabolism specifically targeting mitochondrial and ribosomal genes, providing a segue for future studies targeting tumor mitochondria and ribosome synthesis in TNBC.

Drug preparation-
The 9S1R nullomer peptide (sequence RRRRRWCMNW) was synthesized by Genscript (USA) (lyophilized, HPLC purified, purity > 98%) and stored at − 20 °C.For preparation of 9S1R-NulloPT, briefly, 20 μg/μl stock solution of the 9S1R peptide was prepared in 100 mM Trehalose (Sigma-Aldrich, USA), which was then injected in mice as mg/kg body weight.For example, a 25 g body weight mouse assigned to a dose of 100 mg/kg, received 200 μl of drug formulation containing 100 μg/g bodyweight of 9S1R peptide mixed with 2.5 mM/g bodyweight of trehalose prepared in PBS.Mice in the trehalose only group received 200 μl of 2.5 mM trehalose in PBS per gram bodyweight, and those in the PBS group received 200 μl of PBS.All drug preparations were filtered through a 0.22 μm filter (Millipore, USA) before IP administration.

Single-dose acute toxicity study-
A single dosage of the 9S1R-NulloPT drug at 5, 25, 50, 100 mg/kg and highest equivalent dose of trehalose alone (2.5 mM/g bodyweight) was injected in mice (n = 4) intraperitoneally (IP).Mice were observed every 2 h for clinical signs and symptoms and body weights were taken every 12 h.At the end of the experiment (36 h) mice were euthanized, necropsy was performed and organs preserved in 10% formalin.

Generation of the TNBC mouse model and treatment paradigm-
The model was created by orthotopically injecting 4T1.2-Luc cells in the 4th mammary fat pad of Balb/c mice, with 1 × 10 5 cells per mouse.The cell line was developed by Dr. Cheryl L. Jorcyk [25] and expresses the luciferase gene which serves as a luminescent indicator of gene expression or tumorigenesis, making the tumor cell traceable in vivo.Treatment started after the formation of palpable tumors and was detected by in vivo BLI (IVIS, Perkin Elmer) followed by randomization of the animals per group.For evaluating the effect of the peptides on this model, we carried out two separate studies.The first pilot study involved a set of 4 mice per group with bilateral mammary tumors, which received six IP administration of PBS, 9S1R-NulloPT drug at 50 mg/kg, 100 mg/kg body weight or Trehalose, over a period of 2 weeks.The present study includes a set of 7-9 mice per group with unilateral tumors including the same groups as in the pilot study, but received a total of eight injections per mouse over a period of 2 weeks.The animals were euthanized after 29 days post-tumor cell transplantation, the details are provided in the Supplementary section.The group which received IP injections of 9S1R-NulloPT drug at 100 mg/kg dose is regarded as the treated group, whereas the control group is PBS, unless noted otherwise.

Animal treatment groups-
The 9S1R peptide at dose 100 mg/kg in trehalose is referred to as 9S1R-NulloPT or the treatment group throughout this paper.PBS was used as the control.Results from all groups (PBS, 9S1R-NulloPT at dose 50 mg/kg, 9S1R-NulloPT at 100 mg/kg, and trehalose) are provided in the supplementary materials.

Tumor volume measurement-
Mice were anesthetized by isoflurane followed by measurement of the tumor length and width using a manual Vernier caliper.This was performed before every drug administration and the tumor volume was calculated using the formula (length × width 2 )/2.Body weights were also measured before every dose of the drug.

In vivo bioluminescence imaging (BLI)
Whole-body bioluminescence imaging (BLI) using the In Vivo Imaging System (IVIS ® Spectrum, Perkin Elmer) was performed to detect in vivo tumor burden and metastasis, as evaluated by the bioluminescence signal from the 4T1.2-Luc cells.Mice were injected once with 200 μl of D-Luciferin (150 mg/kg) before in vivo BLI.For confirmation ex vivo BLI was performed from the excised tumor, other organs and secondary metastasis sites, when same dose of D-Luciferin was reinjected once more at the endpoint before necropsy and ex vivo BLI.The images were analyzed by Perkin Elmer software and Aura Version 4.0.7 (Spectral instruments Imaging).

Histology of excised tumors-
Half of the tumor tissue was excised and sent for histopathology analysis by a practicing pathologist at the COBRE-Histopathology Imaging Core at the Boise Veterans Affairs Research Department.The tissues were fixed, paraffin-embedded, sectioned (1 μm) and H&E stained.The slides were analyzed for tumor grade, stage, necrosis, aggressiveness, margin inflammation and immune cell infiltration.Representative brightfield microscope images from the H&E stained slides were captured using ECHO Revolve (Bico, USA) and EVOS M50000 (Invitrogen, USA).Images were captured close to the edge of tumor border and stroma of all tumors.

Half-life of peptide 9S1R in serum-
The half-life of peptide 9S1R was determined by spiking the peptide into fetal bovine serum (FBS) for exposure times of 30 s and 5, 30, 60 or 90 min.During the exposure, FBS containing the spike peptides were incubated at 37 °C.The exposure was halted, and the peptide was extracted from FBS using protein precipitation.Protein precipitation was accomplished using 75% ice cold acetonitrile, followed by incubation at − 20 °C for one hour, and centrifugation at 9000 rpm for 10 min at − 4 °C.The supernatant was removed and saved for analysis.Peptide samples were analyzed by High pressure liquid chromatography (HPLC) mass spectrometry (MS) using an ultra-high-resolution Quadrupole Time of Flight (QTOF) instrument (Bruker maxis, Bruker Corporation, Billerica, MA, USA).HPLC mobile phase consisted of 18 MΩ H 2 O and HPLC grade formic acid and acetonitrile (> 99% purity, Fisher Scientific, Pittsburgh, PA, USA).The electrospray ionization (ESI) source was operated under the following conditions: positive ion mode; nebulizer pressure: 1.2 Bar; flow rate of drying gas (N2): 8 L/min; drying gas temperature: 200 °C; voltage between HV capillary and HV end-plate offset: 3000 V to -500 V; mass range was set from 250 to 2900 m/z; and the quadrupole ion energy was 4.0 eV.Low concentration ESI tuning mix (Agilent Technologies, Santa Clara, CA, USA) was used to calibrate the system in the mass range.HPLC separation was achieved using a Dionex UltiMate ® 3000 RSLCnano system (Dionex Corporation, Sunnyvale, CA, USA) equipped with a Waters XTerra C18 column (4.6 × 100 mm, 3.5 μm) (Waters Corporation, Milford, MA, USA).The mobile phase was 0.1% formic acid in water (Buffer A) and acetonitrile (Buffer B) with a flow rate of 0.2 ml/min.A linear gradient method was used to separate the mixture starting at 5% acetonitrile and ending at 60% acetonitrile over 20 mins.The sample injection volume was 5 μl.Data were analyzed using the Compass Data Analysis software package (Bruker Corporation, Billerica, MA, USA).

RNA sequencing and transcriptomic analysis from excised tumors-
At termination, tumors (n = 3) were excised, snap-frozen and sent to Novogene for RNA sequencing.Paired-end raw sequences (Illumina platform (PE150) were subjected to quality check using FastQC v0.11.9.Sequence reads for all genes were mapped to the Genome Reference Consortium Mouse build 38 (GRCm38) release 99 by HISAT2 v2.1.0using Ensemble ID.Aligned sequence reads were assessed by HTSeq v0.11.3 using the GRCm38 (release 99) ensembl reference genome annotation (gtf format) file.Genes with greater than 10 sequence read counts (for each row considering 3 replicates in a sample) per gene were used for further analysis.Differential gene expression for each comparison (3 replicates for each sample) was performed by default Wald test using Deseq2 v1.38.1 in R. Normalization of sequence reads per gene count for each sample was done using median-ratio-normalization.Volcano plots, PCA plots, heatmap of samples (using Euclidean distances), heatmap of expression of top variable genes were generated using gplots, ggplot2 and RColorBrewer programs (R based).Fold-change for each gene was calculated by comparing counts in all the treatments relative to PBS.Genes with Benjamini-Hochberg (BH) adjusted p-value < 0.05 were considered differentially expressed genes (DEGs).DEGs having log2 fold change value of at least 0.58 or greater were used for investigating significant biological processes, molecular functions, cellular compartments and biological pathways using gProfiler web server.

Interactome and cluster analysis from the RNA sequencing data-
To explore if the DEGs are involved in known and predicted protein-protein interactions, the TNBC-related PPI network was constructed using the online analysis tool STRING (Ver 11.2).The network nodes denote the downregulated and upregulated DEGs and network edges shown in confidence view without any disconnected nodes in the network, and with active interaction sources from experiment, databases, co-expression, neighborhood, gene fusion and co-occurrence [26].The network was clustered to 3 groups following kmeans clustering with hidden edges between clusters for a simplified view.The interaction score was > 0.4 [27].

In vitro Luciferase assay in TNBC cell lines-
4T1.2-Luc cells were plated in white 96 well plates for 48 h without or with 9S1R-NulloPT at concentrations 25, 50 and 100 μM based on previously published IC 50 ranges [17,18].At 48 h post treatment, cells were washed in 1 × PBS, exposed to 150 μg/ml D-Luciferin (Thermo-Fisher Scientific, USA) and incubated in PR free media for 10 mins following previously published protocol [25].Luminescence was measured in a Glomax luminometer (Promega, USA).

Cell viability assay-
Cell viability as a measure of cellular metabolism was examined by MTT assay (Sigma Aldrich, USA) following previously published protocol [28].Briefly, 4T1.2, 4T1.2-Luc and MDA-MB-231 cells were cultured in 96 well plates for 48 h without or with 9S1R-NulloPT at concentrations 25, 50 and 100 μM based on previously published IC 50 ranges [17,18].Post treatment cells were exposed to MTT (Invitrogen, USA) at 5 mg/ml concentration in cell culture media and incubated for 3 h at 37 °C, 5% CO 2 , and 100% humidity.Following media removal, the purple formazan crystals were solubilized in 100 μl DMSO per well and absorbance was taken at 560 nm using a microplate reader (BioTek, Agilent, USA).

In vitro mitochondrial physiology assay-
Measurement of mitochondrial reactive oxygen species (ROS) production was performed using MitoSox dye (Invitrogen, USA), and mitochondrial membrane potential was measured by proton gradient specific TMRM dye (Invitrogen, USA), following previously published protocol [26,29].Briefly, cells were grown in opaque 96 well plates and incubated with 2.5 nM MitoSox for 15 mins, or 20 nM TMRM for 30 mins at 37 °C, 5% CO 2 , and 100% humidity.Fluorescence was measured using a microplate reader (BioTek, Agilent, USA) at an excitation of 488 nm and emission of 510 nm for MitoSox and at an excitation of 540 nm and emission of 570 nm for TMRM.

Confocal imaging of TAMRA tagged 9S1R NulloPT uptake in vitro-
9S1R tagged with Tetramethylrhodamine (TAMRA) at its N-terminus was procured from Genscript (USA) lyophilized, HPLC purified and with > 98% purity.MDA-MB-231 and 4T1.2-Luc cells grown overnight at confocal dishes were treated with nontoxic dose (20 μM) of TAMRA-9S1R NulloPT or trehalose control, following previously published short poly arginine peptide uptake protocol [30,31].The cells with or without peptide were incubated at 37 °C for 2 h in serum free media, followed by washing in PBS.Hoechst 33342 (Invitrogen, USA) at 1x concentration was added to the cells and incubated for 15 min for counterstaining the nuclei.After washing with PBS, live cell images were captured in PR free complete media (GIBCO, USA) using a Leica Stellaris 5 Confocal microscope with attached Uno stage top incubator, which maintains cells at 37 °C and 5% CO 2 .Z-stack images were captured at 40 × magnification using a water immersion objective and maximum intensity projections were created using Fiji software (Ver 1.54 f).

Statistical analysis-
Two-tailed unpaired Student's t-test was used for experimental statistics with p < 0.05 [31]GraphPad Prism (ver 9.5.1) was used for analysis.

Characterization of the peptide 9S1R
The nullomer 9S1R (N-RRRRR-WCMNW-C) has a molecular weight of 1519.81Da and iso-electric point of pH 12.28.It has a net charge of + 5 at pH 7, hydrophobicity + 12.93 Kcal/mol with amphipathic nature and fair solubility in water (Fig. 1A, B).We evaluated the internalization and cellular uptake of TAMRA labeled 9S1R NulloPT after 2 h of incubation in 4T1.2-Luc and MDA-MB-231 cells by confocal microscopy (Fig. 1 C).We found at physiological condition and at a non-toxic dose (20 μM) [30,31] the peptide appears to be internalized within the cells, through endocytosis and primarily located at the cytosol, which is typical for short polyarginine containing peptides [32][33][34].We further confirmed the uptake of the native untagged peptide inside 4T1.2 cells by mass spectrometry (MS) and found that the peptide is available inside the cells at two molecular forms with charges 3 + (m/z 507.2637Da) and 4 + (m/z 380.7002Da), (extracted ion chromatogram (EIC) peaks, Supplementary Fig 1).We determined using MS that the half-life of the peptide in fetal bovine serum (FBS) was 8.7 mins (Fig. 1D).The previous in vitro study by Alileche et al. showed that this peptide targets mitochondrial function [18].This agrees with the above properties of this molecule, in that it is very similar to known mitochondria-targeting peptides (MTPs) [35][36][37][38][39] and, we hypothesize that the peptide ultimately localizes to mitochondria after its endosomal processing in the cytosol.Thus, we determined that the peptide enters cells and has a short half-life in serum.

Single dose acute toxicity study of 9S1R NulloPT in mice
To determine the safety of our peptide formulation in vivo, we performed a single dose toxicity study in female BALB/c mice (Fig. 2 A).The drug doses (administered I.P.) used are 5, 25, 50 and 100 mg/kg, with the control being the highest trehalose dose volume.After 36 h of observation the experiment was terminated, and we found no significant change in the animal behavior, clinically observable signs of acute toxicity or conspicuous change in body weights among any groups (Fig. 2B).There was no significant difference between the drug and trehalose alone group.All the drug doses were well tolerated by the animals, allowing us to move forward with the preclinical in vivo cancer studies.

Treatment with 9S1R-NulloPT inhibits tumor growth in early stage cancer but does not affect tumor volume in late stage growth
We evaluated the effect of 9S1R-NulloPT on female BALB/c mice (8 weeks old), using a syngeneic 4T1.2-Luc orthotropic model of metastatic TNBC.We investigated the effect of the drug on tumor growth over one month in both a bilateral model (initial pilot study, Supplementary Fig 2 .)and a unilateral model (Fig. 3A and Supplementary Fig 3).In both studies, there was a decrease in weight following the first week's treatment, however the mice later returned to near-control level weight (Fig. 3B, Supplementary Fig 2 C).The weight loss was never greater than 15% of the initial body weight, ruling out any adverse effects.We found that although the 9S1R-NulloPT treated mice maintained a reduced tumor volume as compared to untreated controls (Fig. 3C, Supplementary Fig 2B ), the decrease was statistically significant only at the fourth dose (day 18), with a highest inter-group difference (evaluated by Tumgrowth [40]) of 2.4 fold (p < 0.03), followed by the fifth dosage (day 21) with a difference of 1.7 fold (p < 0.2, not significant) (Fig. 3 D, E).There was no statistically significant decrease in any other timepoints, with inter group difference at the third dosage (day 16) being 1.5-fold (p < 0.3), sixth dosage (day 23) 1.7-fold (p < 0.2), seventh dosage (day 25) 1.4-fold (p < 0.3) and eighth dosage (day 25) 1.4 fold (p < 0.4).This indicates a 9S1R-NulloPT therapy responsive window between doses 4 and 5 (Day 18-21) [40].At the endpoint on Day 29, the size and weight of the excised tumors showed no significant difference between the treated and control groups (Fig. 3F, G, Supplementary Fig 3A-D).After harmonizing (aligning the injection days), the combined results of the two studies show a significant decrease in tumor size for the treatment group, with a difference of 2.5 fold (p < 0.015, dose 4), 1.8 fold (p < 0.06, dose 5), and ~2 fold (p < 0.06, dose 6) (Fig. 3H) [35,36,38,39].This suggests that the 9S1R-NulloPT drug is effective in reducing growth during the early phase of tumor development.

9S1R-NulloPT alters tumor metabolism but not metastasis
Monitoring bioluminescence of the tumor cells in vivo revealed that, although the initial tumors (day 8) had similar signals in both control and treatment group, the control group tumors displayed increased in vivo bioluminescence with time.Five of 8 treated animals showed a reduction at the later time points (day 22 and day 28) (Fig. 4 A, B).We found that the excised tumor weights among the groups were not significantly different at the end of the study (Fig. 3); however, the ex vivo BLI of the excised tumors (post mortem) showed a significant decrease in signal with treatment (Fig. 4C, D).This suggests that tumors from the treatment group have a decreased number of metabolically active 4T1.2-Luccells.Firefly luciferase oxidizes luciferin by using ATP, Mg 2 + and O 2 [41], and depletion of cellular ATP results in lower luminescence along with loss of metabolic functioning.It has been previously shown in vitro that 9S1R-NulloPT completely depletes the cellular ATP of BC cell lines within 3 h and are highly lethal to hormone independent and triple negative BC cell lines (MDA-MB-231, BT-549, HS-578 T, and MDA-MB-468), as well as hormone dependent BC cell lines (MCF-7 and T-47D) [17,18].The direct association between decrease of luminescence with changes in cell viability was confirmed in cultured 4T1.2-Luc cells (Fig. 4 E, F).We observed that the 9S1R-NulloPT induced a dose-dependent loss of luminescence and a decrease in cell viability as measured by cellular metabolic activity using MTT assay [42] (Fig. 4E).Additionally, these decreases were found to be significantly positively correlated with each other (Fig. 4

F).
Although there was a trend of reduction, there was no significant change in the bioluminescence of cells metastasizing to lung (Fig. 4 G, H) and no change was seen in the number of metastasis to lungs after treatment (Fig. 4 I).These findings corroborated our previous study (Supplementary Fig 2 A-I), and support the role of 9S1R-NulloPT in altering tumor metabolism yet having no conclusive effect on metastasis [35,36,38,39].

9S1R NulloPT alters tumor immune microenvironment
Histopathological scoring of tumor grade and stage (on a scale of I-IV) from the H&E stained sections was performed (similar to that used for human breast tumors), per NCI's recommendation.All the 4T1.2-Luc tumors were of very advanced stage, and had a necrotic center.Although tumors from all groups were similar in aggressiveness via grade and stage (Fig. 5 A-D), there was ~1.5-fold increase in margin inflammation, with presence of inflammatory cells on tumor margins in the treated group.There was a 3-fold increase in the immune cell infiltration within the treated tumors compared to the control untreated group (Fig. 5E-G).The pathological report also suggests infiltration of plasma cells in the treated tumors, while none were found in the control groups.The results indicate that 9S1R-NulloPT treatment alters the TME and enhances the immune and inflammatory response.This is corroborated by the increased number of plasma cells [43], tumor-infiltrating lymphocytes (TILs) [44,45] and higher immune scoring [46], which are all associated with a positive prognosis in TNBC patients.
Interestingly the Igkv (immunoglobulin kappa variable) family (Igkv 12-44, Igkv 3-4, Igkv 4-58, Igkv 6-15, Igkv 6-23 and Ighe) were among the top 10 downregulated genes (Fig. 6A).We investigated the presence of clusters in the interactome network of the differentially expressed genes (DEGs) from total upregulated and downregulated genes, using STRING software (Ver 11.5) [27].Of the upregulated DEGs, we found clusters specific to cancer, focal adhesion, signal transduction, cell adhesion and nervous system development (Fig. 6C).Network analysis of downregulated genes revealed an overall cluster for metabolic process related genes with specific clusters for mitochondria, ribosome and translation machinery, immune system and myofibril assembly (Fig. 6D).Functional enrichment analysis using the five highest significant gene ontology (GO)-terms (biological process, molecular function and cellular component) and pathways (from Kegg, Reactome and Wikipathways) of upregulated and downregulated DEGs confirmed the results [47,48].The details of the enrichment analysis and list of pathways are provided in Table 1.

9S1R-NulloPT specifically targets metabolic and bioenergetics pathways
We found cancer related pathways, and at least 65 cancer pathway-associated DEGs, were upregulated with treatment.Since, we compared treated TNBC tumors with untreated TNBC tumors, we investigated whether such cancer related pathways and DEGs are contributions of the cancer model or an effect of drug-induced change.To achieve this, we used publicly available RNA-Seq data by Schrors et al. [49], comparing gene expression of 4T1.2-LucTNBC cells (derived from BALB/c mammary gland) versus normal BAlb/C mammary tissue.This comparison identifies the DEGs that result from BC (TNBC) alone.Our 4T1.2-Luc cell line, developed by Jorcyk et al. [25], is a luciferase tagged single cell clone of 4T1.2 cells [50] with similar genomic background and tumor properties to the original 4T1, but with a higher tendency to metastasize.A comparison of the DEGs from 9S1R-NulloPT treated vs untreated tumors (group a), and 4T1 tumor cells vs mammary cells (group b) identified those DEGs that resulted from 9S1R-NulloPT treatment (unique from group a).DEGs resulting from TNBC [51,52] are common to group a and b (Fig. 7).
The comparison between the RNA-Seq results revealed that 688 DEGs were due to TNBC (common to group a and b), with 233 genes upregulated and 455 genes downregulated (Fig. 7A).Cluster analysis of these genes using STRING (Ver 11.5) revealed the TNBC downregulated clusters consist of myofibril assembly, immune related and mitochondrial electron transport chain (ETC)TC related genes.The TNBC upregulated cluster comprises focal adhesion and cancer signaling related genes (Fig. 7 C).Although no noteworthy cluster was found in upregulated genes unique to 9S1R-NulloPT treated cancer, Abl-1 and Shc-1 formed the center node of a small cluster.Remarkably, when analyzed in details (Supplementary Fig 5), among the above-mentioned cancer related pathways we found only 20 associated genes to be uniquely upregulated in 9S1R-NulloPT treated tumors (Supplementary Fig 5 A-C) including pan-cancer related genes such as Shc-1 [53], Mtor [54], Lrp6 and Wnt5a [55].Interestingly, a recent report suggests that Wnt5a [56] and Abl-1 [57] are potent suppressors of TNBC progression and are associated with a better prognosis in BC.We also found 44 DEGs associated with at least 12 cancer pathways that were upregulated in both 9S1R-NulloPT treated and TNBC untreated tumors (Supplementary Fig 5, D-F), signifying a contribution from the TNBC model.A list of cancer pathway related DEGs in TNBC untreated tumors and 9S1R-NulloPT treated tumors is provided in Supplementary figure 6.The unique downregulated clusters, which signify the effect of 9S1R-NulloPT only (and not breast cancer), comprise mitochondrial ATP synthesis coupled proton transport (Atp5e, Atp5j, Atp5k, Atp5h and Atp5j2), mitochondrial electron transport chain (Uqcrb, Uqcrh, Ndufv2 and Ndufs6), Oxidative phosphorylation (Cox6c, Atp5e, Uqcrb, Ndufs6, Atp5j, Cox7a2, Atp5k, Uqcrh, Atp5h, Ndufv2 and Atp5j2), Mitochondrial respiratory chain complex-1(Ndufv2, Ndufs6, Ndufs5, Ndufa5, Ndufb5, Ndufs4 and Ndufb3) and ribosome/ translation related genes (Rpl17, Rpl19, Rps27l, Rpl22l1, Rpl9, Rps14, and Rps27) (Fig. 7 B, E, F).The top 5 unique downregulated genes were from the Igkv family (Igkv 3-4, Igkv 12-44, Igkv 4-58, Ighe and Igkv 6-23) and top 5 unique upregulated genes were Adgrl, Mgat3, Lhx6, Map1a and Spef1 (Fig. 7D).Enrichment analysis of the mitochondria cluster revealed Cristae formation, Complex I biogenesis, Oxidative phosphorylation, Electron transport chain, TCA cycle pathways and the GO-terms such as respirasome, mitochondrial respiratory chain assembly and mitochondrion organization (Fig. 7E).Enrichment analysis of the ribosomal cluster revealed Ribonucleoprotein complex, RNA binding, Ribosomal biogenesis, Translation and Peptide biosynthesis pathways and the GO-terms such as rRNA processing, SRP dependent cotranslational protein targeting to membrane, L13a mediated silencing of Ceruloplasmin expression (Fig. 7F).Interestingly, low ceruloplasmin expression correlates with a favorable prognosis and tumor immune cell infiltration in BC patients [58].We found 149 genes which were upregulated with TNBC (in group-b) were downregulated as an effect of 9S1R-NulloPT treatment (in group-a).These genes comprise the same cluster of Ribosome and Mitochondria genes with similar GO-terms and pathways (Supplementary Fig 6 A-C).Likewise, 99 genes which were downregulated as an effect of TNBC (group-b) were upregulated due to 9S1R-NulloPT treatment (group-a), and these includes Focal adhesion, nervous system development, and ECM organization related genes (Supplementary Fig 6 D-F).In summary, transcriptomic analysis suggests that 9S1R-NulloPT targets metabolic pathway related genes, specifically the Mitochondrial energy metabolism and Ribosome associated pathways such as Ribosome biogenesis and translation.

9S1R-NulloPT reduces mitochondrial membrane potential and increases mitochondrial ROS in triple negative breast cancer cell lines
We confirmed the effect of 9S1R-NulloPT on mitochondrial physiology in in vitro cultures of mouse TNBC cell line 4T1.2-Luc and human TNBC cell line MDA-MB-231.We found that in both cell lines (Fig. 8A, B) treatment induces a dose-dependent decrease in cell viability and cellular metabolism as measured by MTT assay [28,42], and a decrease in mitochondrial membrane potential measured by TMRM fluorescence.We also found an increase in the generation of mitochondrial ROS production with an increased dose of 9S1R-NulloPT as indicated by enhanced MitoSox fluorescence.Additionally, the mitochondrial membrane potential and mitochondrial ROS had a significant negative correlation (r = −0.92,4T1.2-Luc and −0.95, MDA-MB-231), mitochondrial membrane potential and cell viability had a significant positive correlation (r = 0.99, 4T1.2-Luc and 0.81, MDA-MB-231), and mitochondrial ROS and cell viability had a significant negative correlation (r = −0.98,4T1.2-Luc and −0.88, MDA-MB-231) (Fig. 8C, D).These results in two separate TNBC cell lines from two different species along with our previous findings [17,18] demonstrates that 9S1R NulloPT treatment alters cellular metabolism by inducing changes in mitochondrial physiology.

Discussion
This is the first in vivo study to demonstrate the physiological impact of nullomer-based peptide drugs.9S1R-NulloPT (9S1R peptide in trehalose) was very well tolerated in vivo when delivered intraperitoneally in normal mice.Our in vivo RNAseq results support previous work showing that 9S1R-NulloPT targets metabolic pathways, specifically the genes for mitochondrial and ribosomal functions.The effect of this drug on mitochondrial physiology and ATP production has previously been demonstrated in different types of cancer cells [17,18] as well as in several TNBC cell lines.It has been reported that 9S1R-NulloPT reduces ATP formation and mitochondrial membrane potential, and increases mitochondrial reactive oxygen species (ROS) generation.This report, provides further evidence that 9S1R-NulloPT kills TNBC cells by altering cellular metabolism via the mitochondria.A moderate to mild mitochondrial stress and fission is beneficial for TNBC cell growth and aggressiveness, but severe mitochondrial stress and increased fission leads to excessive ROS generation and death of cancer cells [59].The drastic alteration in the mitochondrial pathways in the 9S1R-NulloPT treated tumors shows that the peptides directly or indirectly target mitochondria, which could be exploited as a potential mitochondria-targeted anti-cancer/tumor therapy.Ribosome synthesis in the nucleolus increases in cancer cells to cope with increased demand for protein synthesis.Ribosome biogenesis-targeting 61] and BMH [62]) is still in its infancy, however, it is emerging as an effective cancer therapy used in multiple clinical trials [63].Biogenesis of ribosomes make an interesting target for cancer chemotherapy for many reasons: (1) the inhibition of ribosome biogenesis induces cell cycle arrest in a p53-independent manner [63], (2) these inhibitions don't affect the resting cells, possibly due to the long half-life of cytoplasmic ribosomes [64], and (3) it could lead to apoptosis of neoplastic cells that have a high nucleolar ribosomal biogenesis rate [65].Ribosome synthesis is extremely complex and one of the most energetically demanding cellular activities [66].We hypothesize that the drug 9S1R-NulloPT targets mitochondrial ATP production and ribosome biogenesis leading to a loss of metabolic activity.To ensure more efficient killing of these aggressive TNBC cells in vivo and to enable tumor size reduction, we plan to use our drug synergistically along with other anticancer drugs that act through different cytotoxic pathways [67][68][69][70].
Transcriptomic analysis revealed upregulation of: unique clusters of ECM organization genes (collagen formation and degradation genes such as Col6a3, Col5a3 Plod1, Mmp3, Loxl3 and Plec); focal adhesion genes (Itga5, Lamc1, shc1, Arghap35, and Pxn), and cytoskeletal protein binding genes (Abl1, Map1a, Map1b, Map6, Myo10 and Trak1).Among the highest upregulated genes Lhx6 was noteworthy in relevance to BC, as it is reported to suppress activation of the PI3K/Akt/mTOR signaling, inhibiting the progression of BC [71].The most downregulated genes belonged to the Igkv family, which has been suggested as an identifying biomarker for TNBC cancers [72]; and the Fcmr gene whose knockdown leads to increased phagocytosis, enhanced antigen presentation, and heightened T cell activation.Fcmr is also a promising anti-cancer target [73].
Although treatment did not change tumor size post necropsy, the in vivo as well as ex vivo BLI of the tumors showed reduced bioluminescence in the treated group.There were multiple mice (~50%) in the treatment group that showed a reduction in bioluminescence which implies fewer cells were present in the tumor due to decreased proliferation or increased cell death and loss of metabolic activity of the tumor.This corroborates the results obtained in the pilot study.Ex vivo BLI of the tumors in both the studies showed a significant decrease in bioluminescence and thus a decrease in the metabolic activity of the tumor cells in the treated groups but not in controls.This suggests that although the treatment does not change the size of tumors drastically, it altered the metabolism of the tumor cells rendering them inactive.This makes sense in light of the tumor transcriptomics: metabolism related genes were mostly downregulated in the treated tumors, specifically the mitochondrial and ribosomal genes that are essential for energy production and metabolism.The in vivo BLI from day 22 (24 h after dose 5) clearly shows the reduced bioluminescence in treated tumors during this period.In the pilot study (using a bilateral tumor model) we imaged mice after the 4th 9S1R-NulloPT dose (day 22,Supplementary Fig 4A) and found a similar decrease in in vivo bioluminescence in the treatment groups.
Histopathological evaluation confirmed that treatment and control tumors had the same grade, stage and aggressiveness.However, it was interesting to find that there was a significant change in the tumor immune-microenvironment of the treated tumors.Immune cell infiltration increased significantly in the 9S1R-NulloPT treated tumors, as measured by plasma cell numbers and margin inflammation.A positive correlation of plasma cells with favorable patient outcomes has recently been reported [43].It is also well established that the outcome of immunotherapy treatment, and thus prognosis of BC, is dictated by the tumor microenvironment.Recent reports suggest that the immune score of BC patients could be useful for treatment decisions and survival prediction [46]: activated immune cell infiltration in tumors correlates with better prognosis [44].Increased levels of tumor infiltrating lymphocytes (TILs) have been associated with disease-free status and overall survival rates in TNBC patients with and without any treatment.The presence of TILs in the breast tumor microenvironment can also predict responses to neoadjuvant therapy and adjuvant chemotherapy treatments, and high numbers of TILs correlate with increased pathological complete responses in TNBC [45].We also found that the 9S1R-NulloPT treatment resulted in less externally visible necrosis as compared to untreated tumors or trehalose treated groups (data not shown).There was no conspicuous difference between the PBS control and trehalose groups (Supplementary Fig 4 A-G), suggesting that the immunological changes in the tumor microenvironment (TME) are a specific contribution of the nullomer peptide.
The results of this study should be seen in the light of the fact that the peptides were not stabilized and have a half-life in serum of about 9 min.Small peptides such as 9S1R tend to have shorter half-life in circulation, but when internalized their response could be effective immediately as we have seen in MDA-MB-231 cell lines [18].Cellular internalization studies (Supplementary Fig 1 ) show that the peptide could be identified from crude lysed cytosol fraction (cytosol + mitochondria fraction) in several valence forms (2 +, 3 +, 4 +) from m/z 380.7 ± 0.1 when incubated with 4T1.2 cells for at least two hours.
Biochemical properties of 9S1R are similar to known mitochondria penetrating peptides in terms of charge, hydrophobicity, poly arginine content [74], and the effects on mitochondrial functions [17], all of which suggests a potential for mitochondrial localization.Confocal micrographs of 4T1.2-Luc and MDA-MB-231 cells treated with fluorescently labeled 9S1R NulloPT for two hours under physiological condition (Fig. 1C) shows an initial cytosolic and seemingly endosomal localization of the peptide [30,33].However, further colocalization studies with mitochondria across different time points are warranted to gain insights to the ultimate localization of the peptide.Existing literature supports the idea that short cell penetrating peptides reach into the mitochondria from cytosol after endosomal escape as early as 6 h or even after 24 h in vitro [30,31,33,34].Thus, here we confirmed that the 9S1R-NulloPT are internalized within the cells and are sequestered into other forms, which are more stable than those in serum, however, the mechanism of the 9S1R NulloPT's uptake, metabolism, and final sub-cellular localization is still under investigation.
The Nullomer peptides are readily solubilized in the sugar trehalose, which has known anticancer properties [23,75].In previous work, we used the nullomer-trehalose combination with the peptoprime 9R-and its scrambled version 9S1R to screen the NCI 60 cancer panel and the TNBC model MDA-MB-231 and found them to be effective in vitro.Another peptoprime linked to 5 arginines, did not show anticancer effects, and has been used as a negative control peptoprime [17,18].As noted above, trehalose is reported to exhibit anti-cancer properties, and in this study, we did see some effects (Supplementary Fig 2-5).Trehalose alone had similar potency in tumor size reduction as 9S1R-NulloPT groups in study-2 (Supplementary Fig 3), but a contrasting higher than control tumor size in study-1 (Supplementary Fig 2 ), with no significant change in post mortem tumor size.Similar to the control group, trehalose did not show any effect on the immune cell infiltration or margin inflammation on the tumors (Supplementary Fig 4).The 4T1.2-Luc cells have a specific tendency to metastasize into lungs and bone.We found specific metastasis to the lungs, and although there was a trend of reduced bioluminescence in the treated lungs as compared to controls, it was not statistically significant.There was also no significant change in metastatic foci counts in the lungs.These metastatic findings are still preliminary, primarily because of the large variations in the control groups.
Overall, we show a moderate level of therapeutic potential of 9S1R-NulloPT, in terms of tumor growth.In both of our in vivo studies, the treatment decreased tumor volume moderately with statistical significance during the growth phase of the tumor (day 18) seen by the fourth dose, however at the end of the experiment the control and treatment group had similar tumor weights.This suggests that 9S1R-NulloPT acts during an early therapeutic window in this TNBC model.

1.
The efficacy of the peptide in reducing tumor volume wasn't as large as the metabolic inhibition and this could be due to (a) the short half-life of the peptide in serum, (b) the dramatic rate of proliferation of TNBC cells that may be greater than rate of clearance of the cells from the tumor, making them necrotic and inactive but still large in size.To circumvent these challenges in the future we plan to (a) enhance peptide stability by loading the drug in an LNP or extracellular vesicle carrier (b) resect the tumor after dose 5 and continue drug administration to evaluate the tumor signatures and metastasis by in vivo BLI.

2.
The RNA sequencing results revealed that mitochondrial and ribosomal genes are greatly affected by 9S1R-NulloPT treatment, but this needs to be confirmed by PCR.We plan to identify the targets more directly by the subcellular localization of the 9S1R-NulloPT.

3.
To improve the therapeutic effect of the drug we plan to combine it with Doxorubicin along with 9S1R-NulloPT in in vivo TNBC models.

4.
Transcriptomics analysis should be repeated with mammary tissue from normal mice, those with TNBC tumors, and those with TNBC tumors treated with 9S1R-NulloPT, in at least two different stages (day 18 and day 27).

5.
We did not perform any intravenous administration as we found previously that 9S1R-NulloPT has a low but significant RBC hemolytic activity of 0.5% at 10 μM [17].

Conclusion
Previously NulloPTs showed very promising in vitro results using the NCI 60 panel of cells.The current in vivo mouse study demonstrates that the drug reduces early tumor volume (up to the 5th treatment) with no significant change in terminal tumor mass.

Fig. 1 .
Fig. 1.Characterization of 9S1R nullomer peptide.(A) Projected biochemical structure, blue represents poly-arginine and red is the 5 aa nullomer sequence (WCMNW).(B) physical properties of the peptide as obtained from Pepdraw (©2015).(C) Maximum intensity projection images of confocal sections from 4T1.2-Luc and MDA-MB-231 cells treated with non-toxic dose of 20 μM TAMRA labeled 9S1R NulloPT (red) or unlabeled trehalose control.Nuclei stained with Hoechst 33342 (blue).Arrows (yellow) in DIC merged indicate internalized labeled peptide in cells.Images captured at a magnification of 40x, scale bar represents 20 μm.(D) The extracted ion chromatogram (EIC) peaks of the peptide after LC-MS at different incubation times in serum (0.5-90 min).The X-axis depicts the retention time and Y-axis peak intensity.Inset shows the half-life log 2 peak area under the curve versus treatment time for the EIC peaks.t 1/2 value for the peptide is 8.7 min.

Fig. 2 .
Fig. 2. Toxicity study in mice (A) In vivo safety study timeline.(B) Changes in mean body weight plotted against time among mice treated with trehalose alone or 9S1R-NulloPT drug at 5, 25, 50 and 100 mg/kg; N = 3.

Fig. 3 .
Fig. 3. Effect of 9S1R-NulloPT on tumor volume: (A) Treatment schedule.(B) Body weight over time, in control and 9S1R-NulloPT treated mice.(C) Tumor volume by caliper measurement plotted against days post tumor cell injection with injection timeline of 8 doses, green box highlights drug-responsive therapeutic window with highest inter-group difference.(D) Cross-sectional analysis by Tumor growth [40] showing the highest responsive window on day 18, 4th dose (p < 0.03, n = 6-7) and (E) day 21, 5th dose (p < 0.2, n = 7).(F) Excised tumors at endpoint from control and treatment groups; blue scale in inches.(G) Change in excised tumor weight from control and treatment groups.(H) Tumor volume measurement by harmonizing the data from both current and pilot studies, data aligned by injection timepoints.The green box depicts the responsive window with a statistically significant difference between the control and treated groups by the 4th dose (p < 0.01, n = 10-11).Data presented as Mean+ SEM followed by two-tailed unpaired Student's t-test, *p < 0.05 considered statistically significant.

Fig. 4 .
Fig. 4. Effect of 9S1R-NulloPT treatment on tumor metabolism and lung metastasis: (A) Representative images taken by IVIS Spectrum (Perkin Elmer) on days 8, 22, and 28 post 4T1.2-Luc cell injection, in control and 9S1R-NulloPT treated mice, showing in vivo bioluminescence signal from the tumor cells.Pink star marks no signal in the treatment group.(B) Bioluminescence intensity derived from in vivo BLI at day 22 and 28 showing treated (red) and control (blue) groups with dotted lines denoting a decrease in intensity from individual mice.(C) Ex vivo tumor images on the same luminescent intensity scale from representative control and treated mice.(D) Tumors from the treated group show statistically significant (p < 0.02) loss of bioluminescence represented as photons/second when compared to control tumors.(E) 4T1.2-Luc cells in culture showing dose dependent decrease in cell viability and luciferin induced luminescence, upon treatment with 9S1R-NulloPT.(F) Pearson correlation matrix showing a positive correlation with r = 0.99 (p = 0.007, N = 4-5) between cell viability and luciferin luminescence as found in E. (G) Representative photographs, images from ex vivo bioluminescent imaging (BLI) and fluorescence only signals from lungs, the secondary metastasis sites of control and treated mice.(H) Intensity analysis of total emission from Lungs BLIs showing a trend of reduction (not-significant, p < 0.20) in bioluminescence after treatment.(I) Total count of metastatic foci in lungs with no significant (p < 0.14) change among the groups.Data presented as Mean+ SEM followed by two-tailed unpaired Student's t-test, *p < 0.05 considered statistically significant, N = 5-8.

Fig. 5 .
Fig. 5. Tumor histology: Pathological scoring from the Hematoxylin and Eosin stained sections of mammary tumors from control and treated mice, revealing (A) aggressiveness of the tumor, (B) stage of tumor, (C) tumor grade, (D) score of observed necrosis, (E) marginal inflammation of tumor, and (F) infiltration of immune cells within the tumor.(G) Representative image from the sections at 4x (scale 200 μm), 10x (scale 100 μm) and 40x (scale 25 μm) magnification, showing increased immune cell infiltration in the treated tumor (yellow arrows).Data represented as Mean ± SEM, Two-tailed student's t test with * P ≤ 0.05, N = 5-6.

Fig. 6 .
Fig. 6.RNA sequencing and network analysis: RNA sequencing from advanced stage TNBC tumors (n = 3) showing the effect of 9S1R-NulloPT treatment.(A) Volcano plot showing total DEGs with top 10 upregulated (orange) and downregulated (blue) genes.(B) Top 20 downregulated and upregulated DEGs and their Log2 fold change values (X-axis) and the p-adjusted values (padj) represented as heat map.(C-D) Interactome of the genes following STRING analysis showing DEGs involved in known and predicted protein-protein interactions with network nodes representing upregulated (C) and downregulated genes (D).Network edges (protein-protein interactions) shown in confidence view, the unique enriched GO terms are highlighted to best represent the clusters with counts and GO IDs for upregulated and downregulated interactome.

Fig. 7 .
Fig. 7. Comparative analysis of the effect of cancer vs effect of 9S1R-NulloPT treatment, in mice mammary tissue: DEGs from a published study [49] comparing 4T1.2-LucTNBC cells originating from BALB/c mammary gland vs normal BALB/c mammary tissue (group-b).These were compared with the DEGs obtained from the present study of 4T1.2-Luc tumors treated with 9S1R-NulloPT vs untreated tumors (group-a) to differentiate the drug-induced changes vs cancer-induced changes on the mammary tissue.(A) 132 upregulated and 255 downregulated genes were unique to group-a as an effect of 9S1R-NulloPT treatment, 455 genes were commonly downregulated and 233 commonly upregulated in the two groups, and thus attributed to cancer-related effects.(B) Cluster analysis of the unique DEGs from our study (group a) emphasizes the drug-induced changes.No definite cluster of upregulated genes were found; however, the downregulated genes show two prominent clusters of Mitochondria related genes and Ribosome assembly/Translational machinery related genes.(C) Results from cluster analysis of the common DEGs emphasizing cancer induced gene expression, with focal adhesion related cluster upregulated and at least three distinct downregulated clusters comprising Mitochondrial Electron Transport Chain genes, Immune system-related genes, and Myofibril assembly genes.(D) Analysis of unique clusters from group-a (effect of 9S1R-NulloPT) revealed a list of top 10 up-and downregulated genes with Log2 fold change values (Y-axis) and padj scores inside the columns; a list of unique

Fig. 8 .
Fig. 8. 9S1R-NulloPT shows a dose-dependent alteration on mitochondrial physiology and cell viability, in mouse and human TNBC cell lines: 4T1.2-Luc a mouse TNBC cell line (A) and MDA-MB-231 human TNBC cells (B) showing best fit curve indicating a dose dependent reduction in mitochondrial membrane potential, cell viability, and an increase in mitochondrial ROS production, in in vitro cultures.Pearson's r correlation matrix between mitochondrial membrane potential, cell viability, and mitochondrial ROS, showing positive correlations (blue) and negative correlations (red), with r and p values depicted inside the matrices, in 4T1.2-Luc (C) and MDA-MB-231 (D) cell lines.Heat map for the value of r closer to + 1 (blue) or -1 (red) shows the degree of positive or negative correlation respectively.Data presented as Mean ± SEM (N = 4-6).
The treatment modifies the cancer cell microenvironment by rendering it metabolically inactive, as shown by the reduced tumor bioluminescence, and previous in vitro studies and by the downregulation of metabolic pathway-related genes involved with mitochondrial function, ATP production, and ribosome assembly.These results corroborate our previous findings from the NCI 60 panel, where we have shown that the 9S1R-NulloPT leads to a drastic decrease in mitochondrial membrane potential and ATP generation, and a noticeable increase in mitochondrial ROS generation.

Table 1 -
Functional enrichment analysis from STRING interactome: List of top five (A) GO terms and (B) pathways, based on significance of enrichment or false discovery rate (FDR) shown as a measure of p-values corrected using Benjamini-Hochberg procedure for the upregulated and downregulated DEGs. (A)

Functional enrichments in network: Top 5 GO terms based on significance FDR Upregulated genes GO-term Description Count in network Strength FDR Biological process (Gene ontology)
Biomed Pharmacother.Author manuscript; available in PMC 2024 February 16.