BioMed Research International

BioMed Research International / 2015 / Article

Research Article | Open Access

Volume 2015 |Article ID 124082 | 18 pages | https://doi.org/10.1155/2015/124082

Proteomic Study to Survey the CIGB-552 Antitumor Effect

Academic Editor: Zheng Li
Received11 Mar 2015
Accepted26 Aug 2015
Published20 Oct 2015

Abstract

CIGB-552 is a cell-penetrating peptide that exerts in vitro and in vivo antitumor effect on cancer cells. In the present work, the mechanism involved in such anticancer activity was studied using chemical proteomics and expression-based proteomics in culture cancer cell lines. CIGB-552 interacts with at least 55 proteins, as determined by chemical proteomics. A temporal differential proteomics based on iTRAQ quantification method was performed to identify CIGB-552 modulated proteins. The proteomic profile includes 72 differentially expressed proteins in response to CIGB-552 treatment. Proteins related to cell proliferation and apoptosis were identified by both approaches. In line with previous findings, proteomic data revealed that CIGB-552 triggers the inhibition of NF-κB signaling pathway. Furthermore, proteins related to cell invasion were differentially modulated by CIGB-552 treatment suggesting new potentialities of CIGB-552 as anticancer agent. Overall, the current study contributes to a better understanding of the antitumor action mechanism of CIGB-552.

1. Introduction

Peptides, owing to their feasible rational design, high specificity, and low toxicity, have become attractive therapeutic agents to treat different diseases [1]. Ongoing advances in cancer therapy include the discovery of peptides with a potent antineoplastic effect [2]. Peptides inhibiting angiogenesis or blocking protein-protein interactions have already been evaluated as potential anticancer agents [3]. Additionally, proapoptotic peptides or peptides functioning as receptor antagonists have been proposed to restrict tumor progression [3].

CIGB-552 is a second-generation peptide derived from the antitumor peptide L-2. Initially, L-2 was identified by screening a peptide library corresponding to the region 32–51 of Limulus anti-LPS factor (LALF) [4]. The cytotoxic effect of L-2 was demonstrated on different murine and human tumor cell lines [4]. A transcriptomic study on tumor cells suggested that L-2 could induce apoptosis by modulating glycolysis, mitosis, protein biosynthesis, and other cancer related processes [4].

Such biological findings, in combination with the peptide cell-penetrating capacity, made L-2 an attractive lead molecule for further structural optimization. Therefore, the primary sequence of L-2 was modified, including substitution by D-amino acids and N-terminal acetylation [5]. These modifications increased the antitumor effect of the resultant peptide known as CIGB-552 [5, 6].

The CIGB-552 in vitro antineoplastic effect has been documented by using human cancer cells of different histological origins [5]. In vivo, a significant reduction in tumor growth after treatment with CIGB-552 was demonstrated in both syngenic murine tumors and patient-derived xenograft models [6]. Additionally to induce apoptosis in the tumor mass, CIGB-552 administration decreased the microvessels’ density in the human HT-29 xenograft tumor model suggesting its antiangiogenic effect [6].

Significant advances have been made to discover the functional mediators of CIGB-552 biological response. Copper metabolism (Murr1) domain-containing protein 1 (COMMD1) has been identified as a major target of CIGB-552. Such interaction increased the COMMD1 stability and induced ubiquitination of RelA with subsequent inhibition of the antiapoptotic activity regulated by NF-κB [5]. Nevertheless, the proteome regulated by CIGB-552 can provide new insights to support the antitumor action mechanism of such peptide.

In this work, to identify additional CIGB-552 targets and modulated proteins in tumor cells, two proteomics approaches were used: the chemical proteomics and the expression-based proteomics. The proteomics-derived data suggested that CIGB-552 could be considered as a multitarget drug, which exerts its antitumor effect by modulating proteins related to NF-κB activation, cell cycle regulation, and apoptosis. Network analysis of proteomic results indicates the molecular basis by which the CIGB-552 peptide regulates cancer related processes. The analysis provided here is the starting point for further investigations about the role of CIGB-552 as an anticancer drug.

2. Materials and Methods

2.1. Peptide Synthesis

CIGB-552 is a cell-penetrating peptide with sequence Ac-HARIKPTFRRLKWKYKGKFW, where proline and leucine are D-amino acids; and the N-terminal was blocked by acetylation [5]. The CIGB-552 peptide was synthesized on solid phase using the Fmoc strategy; it was purified by reverse phase high performance liquid chromatography (RP-HPLC) to >95% purity on an acetonitrile/H2O-trifluoroacetic acid gradient [7] and confirmed by electrospray mass spectrometry (Micromass, UK). For chemical proteomics, the CIGB-552 peptide biotinylated at the N-terminal end (CIGB-552-B) was synthesized using the same procedure.

2.2. Cell Culture and Treatments

For proteomic studies two CIGB-552 sensitive cancer cell lines were selected, the larynx carcinoma and colon adenocarcinoma cells, Hep-2 and HT-29, respectively, which were obtained from ATCC (Rockville, MD). Both cell lines were cultured at 37°C and 5% CO2 in RPMI 1640 (Life Technologies, USA) supplemented with 10% fetal bovine serum (FBS; PAA, Canada) and 30 μg/mL gentamicin (Sigma, USA).

For expression-based proteomics, 20 × 106 HT-29 cells were seeded in appropriate vessels and incubated for 24 h. Subsequently, the tumor cells were incubated with 150 μmol/L of CIGB-552 during 40 min, 2 h, and 5 h.

2.3. Chemical Proteomics
2.3.1. Isolation of Proteins

For chemical proteomics, 20 × 106 Hep-2 cells were seeded in appropriate vessels and cultured for 24 h. Subsequently, Hep-2 cells were collected by centrifugation, washed twice with cold phosphate buffered saline (PBS), and lysed in hypotonic PBS solution (0.1x), containing 1 mM of DDT (Sigma) and complete protease inhibitor (Roche, USA), by three freeze-thaw (37°C) cycles. Cellular lysate was cleared by centrifugation at 12 000 rpm at 4°C for 15 min.

2.3.2. Affinity Purification

The CIGB-552-B peptide was incubated with 50 μL of streptavidin-sepharose matrix (binding capacity: 300 nmol/mL; GE Healthcare, USA) for 1 h. As a control, streptavidin-sepharose matrix without CIGB-552-B peptide was used. Then, 300 μg of total protein from Hep-2 cell lysate, as determined by the Bradford assay (Bio-Rad, USA), was added to each matrix and incubated for 2 h at 4°C. The streptavidin-sepharose matrix was collected by short spin and extensively washed with PBS containing 1 mM DTT and 0.5% NP40. Retained proteins were eluted by heat denaturing with SDS sample buffer.

2.3.3. Sample Preparation and LC-MS/MS Analysis

Affinity-purified proteins were reduced, alkylated with acrylamide, and separated by sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE). Analytical and preparative gels were silver [8] and Coomassie blue stained, respectively. Each lane of the preparative gel was cut into 12 slices that were faded, washed, and in situ digested with sequencing grade trypsin (Promega, USA) during 18 h at 37°C. The resulting peptide mixtures were extracted and desalted with stage tips (Thermo Scientific, USA).

Purified samples were analyzed in an Agilent 1100 series nano LC system (Agilent, USA) coupled online to a QTof-2 tandem mass spectrometer (Micromass, UK). The capillary and cone voltages of the electrospray ionization source were operated with 1.8 kV and 35 V, respectively. Samples were applied at 20 μL/min to a PepMap C18 Precolumn Cartridge (5 mm × 300 μm i.d.) from LC-Packings (USA) and were extensively desalted for 10 min using 0.1% formic acid. The precolumn was switched back onto a C18 capillary column (15 cm × 75 μm i.d., packed with 5 μm, Zorbax 300 SB) and the tryptic peptides were separated using a mobile phase containing 0.1% formic acid, 5–45% acetonitrile gradient over 90 min at 300 nL/min flow rate. Survey scans were acquired during 1 s and a maximum of 4 concurrent MS/MS acquisitions were only triggered for 2+, 3+ charged precursors ions detected at an intensity above a threshold of 15 counts/s. Each MS/MS acquisition was completed and switched back to MS mode when the total ion current fell below a threshold of 2 counts/s or after a maximum of 6 s of continuous acquisition. Data acquisition and processing were performed using MassLynx v3.5 (Micromass, UK).

2.3.4. Protein Identification

Acquired data were searched against the human proteins in the UniProtKB database using MASCOT (version 2.2, Matrix Science, UK) [9]. Search parameters were set to a mass tolerance of 1.2 Da for the precursor ions and 0.6 Da for the fragment ions. One trypsin missed cleavage site was allowed. Propionamide-cysteine and oxidized methionine were set as fixed and variable modifications, respectively. MS/MS spectra of identified proteins with one or two peptides were manually inspected. The identification of a protein or peptide was considered positive using the consensus of several criteria: the peptide score >20, the assignment of four intense consecutive fragments ions in the MS/MS spectrum, and the most intense signals which must be explained considering the proposed sequence.

2.3.5. Bioinformatics Analysis

Functional classification of identified proteins was based on the information annotated in the Gene Ontology (GO) database (http://www.geneontology.org). The analysis was performed using the functional enrichment tool GeneCodis (version 3.0) (http://genecodis.cnb.csic.es/) [10]. To identify significant enriched biological processes ( values lower than 0.05), the hypergeometric distribution and the false discovery rate (FDR) correction method were computed by GeneCodis, as statistical analysis. Protein complexes associated with the CIGB-552 target profile were identified by using the CORUM database (http://mips.helmholtz-muenchen.de/genre/proj/corum) [11].

A target deconvolution strategy was applied to filter the CIGB-552 target profile. Nonspecific proteins, retained in unloaded streptavidin-sepharose matrix (without CIGB-552-B), were ruled out. Further, proteins reported by Burkard et al. [12] as part of the central proteome were subtracted from the analysis. To identify contaminants or background proteins the potential CIGB-552 target profile was queried against the CRAPome database (http://www.crapome.org). For each protein, the average spectral count was retrieved as a measure of its abundance in affinity purification followed by mass spectrometry (AP-MS) experiments [13].

Functional subnetworks perturbed by CIGB-552 were identified by using drugDisruptNet (http://bioinformatics.cemm.oeaw.ac.at/index.php/downloads-left/87-disruption-of-functional-networks). The impact of CIGB-552 on functional subnetworks was estimated by the score , which was calculated as described by Burkard et al. [14] but with some modifications. In this regard, the affinity of CIGB-552 to its targets (affinity score) was computed irrespective of the protein abundance in a competitive pull-down (fixed to 1 for all potential CIGB-552 targets). Hence, the affinity score was set to be only proportional to the amount of protein pulled down by CIGB-552-B. As a measure of protein abundance, the exponentially modified Protein Abundance Index (emPAI) of each identified protein was retrieved from MASCOT results. Briefly, the emPAI value is calculated as (PAI = observed peptides/observable peptides), being the number of different observed peptides that cover the protein sequence, a rough estimate of the protein amount in mass spectrometric analysis [15].

2.4. Expression-Based Proteomics
2.4.1. Isolation of Cytosolic Proteins

HT-29 cells were collected by trypsinization. After washing with PBS, the cells were suspended in 700 μL of the isotonic buffer containing 10 mM Tris-HCl adjusted to pH 7.5, 0.25 M sucrose, 1 mM EGTA, and protease inhibitors. For plasma membrane solubilization Triton X-100 at final concentration of 0.25% was added. After 15 min at 4°C, the cell lysate was centrifuged for 15 min at 12 000 rpm and 4°C. The supernatant containing cytoplasmic proteins was kept at −70°C until subsequent analysis.

2.4.2. Tryptic Digestion and Isobaric Labeling (iTRAQ)

The cytoplasmic protein extracts (120 μg) of CIGB-552-treated and untreated (control) HT-29 cells were precipitated with acetone/TCA and independently dissolved in 20 μL of buffer containing 2 M guanidinium hydrochloride (GuCl) and 500 mM tetraethylammonium bromide (TEAB), pH 8.5. Proteins were reduced using 49 mM tris-2-carboxyethyl phosphine (TCEP) at 60°C for 1 h and alkylated with 84 mM iodoacetamide for 30 min at ambient temperature in the dark. The pool of reduced and S-alkylated proteins was diluted until 50 μL by adding 500 mM TEAB buffer. Proteins were digested with sequencing grade trypsin (Promega, USA) at an enzyme-to-substrate mass ratio of 1 : 10 for 18 h at 37°C.

Tryptic peptides were labeled with iTRAQ reagent according to the manufacturer’s protocol (Applied Biosystems, USA). Briefly, each vial of iTRAQ reagent (114–117 tags) was dissolved in 70 μL ethanol by vortexing for 1 min at room temperature. Equal amounts of tryptic peptides derived from different samples were labeled by adding iTRAQ reagent 114 (untreated HT-29 cells, control), iTRAQ reagent 115 (HT-29 cells treated with CIGB-552 for 40 min), iTRAQ reagent 116 (HT-29 cells treated with CIGB-552 for 2 h), and iTRAQ reagent 117 (HT-29 cells treated with CIGB-552 for 5 h). The reaction for iTRAQ labeling was incubated at room temperature for 1 h. Labeled peptides were mixed and dried in a centrifugal evaporator.

2.4.3. LC-MS/MS Analysis

A high pH reversed-phase chromatography step was used to separate the complex mixture of peptides prior to LC-MS/MS analysis. The labeled peptides were resuspended in 500 μL buffer A (0.1% NH4OH, pH 10.5) and separated in 24 fractions using a 4.6 mm × 10 cm RP column packed with POROS R2 resin (Applied Biosystems, USA). Peptides were eluted at a flow rate of 0.8 mL/min using a stepwise buffer B (0.1% NH4OH, 60% acetonitrile, pH 10.5) gradient. Collected fractions were acidified by adding 200 μL of 5% formic acid, lyophilized, and further dissolved in 200 μL of 0.2% formic acid.

For each fraction, three replicate aliquots of 40 μL were separated in independent LC-MS/MS experiments by reverse phase chromatography. Peptides were separated and analyzed using an Agilent 1100 series nano-LC system (Agilent, USA) coupled online to a QTof-2 orthogonal hybrid tandem mass spectrometer (Micromass, UK) operated as described above. Peptides were eluted at a flow rate of 300 nL/min using a 60 min gradient starting with 5% acetonitrile to 45% acetonitrile with a two-buffer system (Buffer A: 0.2% formic acid; Buffer B: 0.2% formic acid, 80% acetonitrile). Data dependent acquisition MS/MS spectra of the eluted peptides were acquired in three m/z ranges (400–600, 590–900, or 890–2000) using the gas-phase fractionation approach [16]. Subsequent analysis proceeded as described above.

2.4.4. Protein Identification

Raw files were processed using MASCOT Distiller software (version 2.3). Peptides were assigned to MS/MS spectra using MASCOT search engine (version 2.2) against the human proteins in the UniProtKB database. The following search parameters were selected: 1.2 Da precursor mass tolerance, 0.6 Da daughter ions mass tolerance, and tryptic search with up to one missed cleavage site. Variable modifications including deamidation of glutamine and asparagine, methionine sulfoxide, and the side reaction of iTRAQ labeling on tyrosine residues were taken into account. Carbamidomethyl cysteine and iTRAQ label on lysine and at the peptide N-terminus were specified as fixed modifications. To accept a peptide hit as positive we considered a false discovery rate (FDR) of 3% based on the target-decoy strategy [17].

Quantification was achieved using ISOTOPICA software [18, 19]. The software enables the relative peptide quantification based on the detailed analysis of the observed isotopic ion distribution. The software calculated the best ratio of the reporter ions (tags 114, 115, 116, and 117) to obtain an in silico isotopic ion distribution best matching with the isotopic ion distribution observed experimentally. To evaluate the quality of this adjustment, the software calculates the difference between the area of both in silico and experimental isotopic ion distributions, and it is expressed as a GOF coefficient (goodness of fitting). The relative quantifications of peptides with GOF below 0.8 were not considered for further analysis. Additionally all relative quantifications were manually inspected. The relative peptide expression ratios (fold changes) were determined with respect to the control sample (tag 114). The fold changes of all peptides corresponding to the same protein were averaged. The resulting protein fold changes were normalized; population median and standard deviation set the threshold ratio to consider, with a value < 0.1, a protein as differentially expressed in each condition (HT-29 cells treated with CIGB-552 for 40 min, 2 h, or 5 h) with respect to control (untreated HT-29 cells).

2.4.5. Bioinformatics Analysis

The functional classification of differentially expressed proteins and enrichment analysis were performed as described above. For biological network analysis, the Cytoscape software (version 2.8) [20] and accessory applications were used. Protein-protein interaction network was generated using the BisoGenet application (version 1.4) [21] which retrieves information from interaction databases including DIP, BioGrid, HPRD, and BIND. The fold changes of differentially expressed proteins were visualized in the network context using the MultiColoredNodes application [22].

3. Results and Discussion

CIGB-552 is a peptide-based drug with antitumor effect and cell-penetrating capacity [4]. The first evidence about CIGB-552 mechanism of action has been already reported [5]. CIGB-552 interacts with COMMD1 protein, increasing its abundance levels. Such evidences were corroborated by Western blot and immunofluorescence detection of COMMD1 in human cancer cells of different histological origins [5]. CIGB-552 induces cellular cytotoxicity in a variety of tumor cell lines [4, 5]. Among them, the Hep-2 larynx carcinoma cell line is highly sensitive to CIGB-552 cytotoxic effect.

To identify other proteins that interact with CIGB-552 a chemical proteomic approach was then conducted (Figure 1(a)). PBS-soluble proteins derived from Hep-2 cells were incubated with CIGB-552-B previously bound to the streptavidin-sepharose matrix (Figure 1(a)). As a negative control, the same pool of proteins was also incubated with the unloaded streptavidin-sepharose matrix. CIGB-552 interacting proteins were eluted and analyzed by SDS-PAGE (Figure 1(b)). The electrophoretic pattern of the whole cell extract and CIGB-552-matrix eluted proteins showed appreciable differences in terms of band intensities indicating the potential CIGB-552 target profile. The proteins identified in the negative control correspond to proteins that interact with the matrix, known as either nonspecific binders or sticky proteins.

A total of 265 proteins were identified by mass spectrometry, of which 104 proteins were also identified as nonspecific binders (see Supplementary Table Sa in the Supplementary Material available online at http://dx.doi.org/10.1155/2015/124082). Therefore, 161 proteins constitute the potential CIGB-552 target profile identified by chemical proteomics. Biological processes related to carbohydrate metabolism, protein modification, and cell cycle are significantly represented on this dataset (Figure 2). Interestingly, such biological processes are also represented in the transcriptomic profile regulated by L-2 peptide in Hep-2 tumor cells [4].

The potential CIGB-552 target profile identified in vitro includes five biological complexes (Figure 3). Components of the minichromosome maintenance (MCM) complex were identified as potential CIGB-552 targets. The MCM complex is related to DNA replication and cell cycle regulation. Therefore, the MCM proteins are frequently upregulated in different cancer types including meningioma, lung cancer, and laryngeal carcinoma [2326]. Four components of the 60S ribosomal large subunit were also identified. Accordingly, translation is a biological process overrepresented in the CIGB-552 target profile (Figure 2). Related to protein modification process, five components of the ubiquitin E3 ligase complex were found. In addition, nine subunits of the proteasome complex were included in the CIGB-552 target profile.

Chemical proteomics do not distinguish between drugs direct and indirect binders. In addition, a high background is typically identified in these experiments [27]. Therefore, the 161 proteins found in the chemical proteomic profile might include nonspecific binders of CIGB-552. Potential solutions to overcome the specificity problem of chemical proteomics are (i) to include unrelated drugs or the matrix itself as negative controls, (ii) to identify and subtract from the analysis the abundant proteins or core proteome, and (iii) to perform a competitive pull-down to rank the most relevant targets [28].

In the present study, besides identifying nonspecific binders, the core proteome was subtracted from the potential CIGB-552 target profile. A set of 1124 proteins previously identified as part of the human central proteome was used [12]. The potential CIGB-552 target profile (161 proteins) includes 106 proteins of the central proteome. The majority of these proteins (91/106) were also reported with an average spectral count greater than two in the CRAPome database (Supplementary Table Sa-b). The spectral count is used as a measure of protein abundance in affinity purification (AP) experiments [13]. Therefore, the CRAPome repository gives information about abundant proteins or common contaminants of AP-MS experiments.

By subtracting the central proteome, 55 proteins were identified as the most probable targets of CIGB-552 (Table 1). This result suggests the potentialities of CIGB-552 as a multitarget drug. In contrast to previously published results [5], the COMMD1 protein was not identified in the chemical proteomic profile, probably due to its low expression levels in tumor cells [29].


UniProt ACCDescriptionGene symbolScoreemPAI

O95218Zinc finger Ran-binding domain-containing protein 2ZRANB2390.21
Q9NX65Zinc finger and SCAN domain-containing protein 32ZSCAN32260.05
Q9UBQ0Vacuolar protein sorting-associated protein 29VPS29560.66
Q14376UDP-glucose 4-epimerase GALE380.1
Q12888Tumor suppressor p53-binding protein 1TP53BP1280.02
O14773Tripeptidyl-peptidase 1TPP1760.06
Q13428Treacle proteinTCOF1270.02
P04183Thymidine kinase, cytosolicTK1240.15
Q9NXG2THUMP domain-containing protein 1THUMPD1580.09
Q8WW59SPRY domain-containing protein 4SPRYD4390.16
P63151Serine/threonine-protein phosphatase 2A 55 kDa regulatory subunit B alpha isoformPPP2R2A870.07
Q9NVA2Septin-11SEPT11290.08
O00764Pyridoxal kinasePDXK590.11
Q16740ATP-dependent Clp protease proteolytic subunit, mitochondrialCLPP670.26
Q9Y4X5E3 ubiquitin-protein ligase ARIH1ARIH1290.06
O15460Prolyl 4-hydroxylase subunit alpha-2P4HA2410.12
P40261Nicotinamide N-methyltransferaseNNMT460.27
P19105Myosin regulatory light chain 12AMYL12A2002.45
Q16539Mitogen-activated protein kinase 14MAPK14820.19
Q9UNF1Melanoma-associated antigen D2MAGED22380.18
P50213Isocitrate dehydrogenase [NAD] subunit alpha, mitochondrialIDH3A640.2
P05161Ubiquitin-like protein ISG15ISG15340.21
O00629Importin subunit alpha-3KPNA4490.06
Q8WTS6Histone-lysine N-methyltransferase SETD7SETD7420.09
P06737Glycogen phosphorylase, liver formPYGL400.04
P30712Glutathione S-transferase theta-2GSTT2960.14
P47755F-actin-capping protein subunit alpha-2CAPZA21090.38
P52907F-actin-capping protein subunit alpha-1CAPZA11180.54
P41091Eukaryotic translation initiation factor 2 subunit 3EIF2S3430.07
O95571Persulfide dioxygenase ETHE1, mitochondrialETHE1350.13
Q9BVJ7Dual specificity protein phosphatase 23DUSP23500.23
P52564Dual specificity mitogen-activated protein kinase kinase 6MAP2K6290.1
Q14691DNA replication complex GINS protein PSF1GINS1240.16
Q3LXA3Bifunctional ATP-dependent dihydroxyacetone kinase/FAD-AMP lyase (cyclizing)DAK450.06
P21291Cysteine and glycine-rich protein 1CSRP1470.39
Q99627COP9 signalosome complex subunit 8COPS8690.35
Q9UBW8COP9 signalosome complex subunit 7aCOPS7A460.12
P61201COP9 signalosome complex subunit 2COPS21360.23
Q07021Complement component 1 Q subcomponent-binding protein, mitochondrialC1QBP6481.74
Q9NT62Ubiquitin-like-conjugating enzyme ATG3ATG3580.1
Q9BZZ5Apoptosis inhibitor 5API51730.13
P12814Alpha-actinin-1ACTN113451.76
P11766Alcohol dehydrogenase class 3ADH5780.19
P36404ADP-ribosylation factor-like protein 2ARL2380.4
P40616ADP-ribosylation factor-like protein 1ARL1390.19
P84077ADP-ribosylation factor 1ARF1990.66
P00813Adenosine deaminaseADA350.09
Q9BTT0Acidic leucine-rich nuclear phosphoprotein 32 family member EANP32E1690.41
P8409860S ribosomal protein L19RPL19380.16
P61586Transforming protein RhoARHOA320.17
Q07960Rho GTPase-activating protein 1ARHGAP1630.07
Q14657EKC/KEOPS complex subunit LAGE3LAGE3340.26
Q96IU4Alpha/beta hydrolase domain-containing protein 14BABHD14B670.17
P39687Acidic leucine-rich nuclear phosphoprotein 32 family member AANP32A2861.38
O75607Nucleoplasmin-3NPM3330.2

UniProtKB/Swiss-Prot entry accession number (http://www.uniprot.org/).
Recommended gene name (official gene symbol) as provided by UniProtKB/Swiss-Prot.
MASCOT protein score.
emPAI: exponentially modified Protein Abundance Index.

The central proteome is mainly enriched in cell vital processes [12]. Consequently, after the central proteome subtraction, biological processes such as translation, transport, response to stress, and cell death are no longer significantly represented in the CIGB-552 chemical proteomic profile. However, proteins related to cell cycle, carbohydrate metabolic process, and signal transduction are corroborated as mediators of the CIGB-552 antitumor activity (Supplementary Table Sa).

Target deconvolution based on the subtractions of sticky and core proteome proteins might remove a real drug target [28]. Nevertheless, this strategy has been previously used to identify the targets of BCR-ABL kinase inhibitor INNO-406 in chronic myeloid leukemia [30]. To overcome the potential disadvantage of removing real targets, a computational strategy was applied for identifying the functional subnetworks perturbed by CIGB-552. Such computational approach, named drugDisruptNet, was presented by Burkard et al. to predict the mechanism and potential side effects of drugs [14]. A functional subnetwork is defined as a connected fraction of the interactome in which all the proteins share the same function [14]. A drug can impact functional subnetworks directly (the drug target is part of the subnetwork) or indirectly (the drug target interacts with the protein subnetwork) [14]. Therefore, even if a CIGB-552 target is removed during the target deconvolution, the remaining potential targets might indicate the functional subnetworks perturbed by the drug.

To identify the functional subnetworks perturbed by CIGB-552, the filtered target profile (55 proteins) was used. As other studies have corroborated the role of COMMD1 in CIGB-552 cytotoxic effect [5], this protein was included in the network analysis. According to drugDisruptNet results (Table 2), different subnetworks related to carcinogenesis were perturbed by CIGB-552, “antiapoptosis” and “negative regulation of cell cycle” among them. These results corroborate the antitumor effect of CIGB-552.


Perturbed functional subnetworksGO id.NodesTargetsScore

Insulin receptor signaling pathwayGO:000828657150.175
Positive regulation of biosynthetic processGO:000989173190.173
Positive regulation of NF-κB transcription factor activityGO:005109240140.165
Protein autoprocessingGO:001654067140.165
Extracellular structure organizationGO:004306266130.143
Mitochondrion organizationGO:000700568170.136
AntiapoptosisGO:0006916178260.112
Interleukin-6 productionGO:00326352180.105
Response to hypoxiaGO:000166668140.103
Negative regulation of cell cycleGO:00457863780.076

GO id.: Gene Ontology identifier.
Number of proteins related to the functional subnetwork and identified as potential CIGB-552 targets.
: score calculated by drugDisruptNet, predicting the impact of CIGB-552 over the functional subnetwork.

Additionally, “extracellular structure organization” and “response to hypoxia” were modulated by the CIGB-552 target profile (Table 2). Both the composition and organization of the extracellular matrix and hypoxia are microenvironment signals that contribute to metastatic spread of cancer cells [31]. Therefore, CIGB-552 could inhibit metastasis in treated cells.

The “positive regulation of NF-κB transcription factor activity” is disrupted by CIGB-552 (Table 2). Although the target profile is not annotated in this signaling pathway, it acts at the periphery of the functional subnetwork (Figure 4). Previous results demonstrated that CIGB-552 treatment upregulates COMMD1 levels [5]. Such event increases RelA ubiquitination and consequently inhibits the NF-κB signaling pathway [5]. In the present study, COMMD1 was not identified in the target profile. However, the “positive regulation of NF-κB transcription factor activity” is disrupted by the CIGB-552 target profile essentially at two network nodes: RELA and TRAF6 (Figure 4).

NF-κB transcription factor is a homo- or heterodimeric complex composed of REL proteins (RelA/p65, RelB, c-Rel, p50/p105, and p52/p100). RelA/p50 appears to be the most common heterodimer. After activation of I-kappa-B kinase complex (IKK), the active NF-κB complex is translocated into the nucleus to regulate gene transcription [32]. NF-κB target genes promote tumor cell proliferation, inhibition of apoptosis, migration, inflammation, and angiogenesis [33]. Therefore, upregulated NF-κB activity has been reported in different tumor types including laryngeal and pancreatic cancers [3436].

The TNF receptor-associated factor 6 (TRAF6) activates NF-κB [37]. This mechanism depends on TRAF6-mediated TGF-beta-activated kinase 1 (MAP3K7/TAK1) activation, which in turn phosphorylates and activates the I-kappa-B kinase complex (IKK) [38]. MAP3K7/TAK1 also phosphorylates the dual specificity mitogen-activated protein kinase kinase 6 (MAP2K6), which in turn activates the p38 MAPK signaling pathway [37]. In the present study, MAP2K6 and MAPK14 (mitogen-activated protein kinase 14) were identified as potential CIGB-552 targets (Figure 4). MAPK14 is a member of the p38 MAPK family. The p38 MAPKs phosphorylate different substrates, including the nuclear mitogen- and stress-activated protein kinase 1 (RPS6KA5/MSK1). MSK1 phosphorylates RelA and increases the transcriptional activity of NF-κB [39]. Therefore, CIGB-552 could inhibit the NF-κB signaling pathway in tumor cells by blocking the MAP2K6/MAPK14 functions.

Furthermore, persulfide dioxygenase ETHE1 (ETHE1) and lysine N-methyltransferase 7 (SETD7) were identified as potential CIGB-552 targets (Figure 4). Other studies demonstrate that both proteins regulate RelA transcriptional activity. ETHE1 protein promotes accumulation of RelA in the cytoplasm and consequently inhibits the NF-κB transcriptional activity [40]. SETD7 protein monomethylates RelA subunit; this event triggers RelA degradation and downregulates NF-κB target gene expression [41, 42]. The identification of ETHE1 and SETD7 in the target profile suggests that CIGB-552 could exert, by interacting with these proteins, different mechanisms to regulate RelA function.

Interleukin-6 (IL-6) expression is activated by NF-κB transcription factor [43, 44]. The “interleukin-6 production” subnetwork was perturbed by CIGB-552 target profile (Table 2). Therefore, disruption of this subnetwork could be a consequence of inhibiting NF-κB transcriptional activity. Importantly, IL-6 downregulation could be validated as a subrogate marker of CIGB-552 antitumor effect.

CIGB-552 target interactions trigger cellular signaling events that finally accomplish the antitumor response. To investigate these signaling events at the protein level, an expression-based proteomics study was conducted. The majority of proteins identified by chemical proteomics as potential CIGB-552 targets (109/161) are located in the cytosol. Cytoplasmic proteins are also enriched after target deconvolution (36/55) (Supplementary Table S2b-c). Therefore, the cytoplasmic proteome regulated in the presence of 150 μmol/L of CIGB-552 was studied using the HT-29 cell line as target cells. Such peptide dose represents the inhibitory concentration 80 (IC80) for CIGB-552 in HT-29 cells (unpublished results). Confocal microscopy demonstrates that after 10 min of incubation CIGB-552 is able to penetrate the cells [4]. In the present work, the proteomic profile modulated by CIGB-552 was investigated at three incubation times: 40 min, 2 h, and 5 h.

In both proteomic experiments, different tumor cell lines were used to evaluate consistence among results. HT-29 cells, as well as Hep-2 cells, have shown to be sensitive to CIGB-552 anticancer effect [4]. In addition, CIGB-552 was able to elicit significant antitumor activity in both murine CT-26 (colon carcinoma cells) and human HT-29 implanted tumors [6]. Consequently, the differentially expressed proteins in treated HT-29 cells could support the findings of the CIGB-552 target profile identified in total cell extracts from Hep-2 cells.

As a result, 658 protein hits were identified (Supplementary Table S3), of which 72 proteins were differentially modulated in at least one experimental condition (CIGB-552: 40 min, 2 h, and 5 h) with respect to control (Table 3). A total of 68 proteins of the proteomic profile (658 protein hits), including eight differentially modulated proteins in response to CIGB-552 treatment, were also identified by chemical proteomics as potential CIGB-552 targets (Supplementary Table S3).


UniProt ACCDescriptionGene symbolScoreFold change
40 min (2.6)2 h (2.6)5 h (2.4)

Downregulated proteins
Q14566DNA replication licensing factor MCM6 MCM676.01−4.2−1.1−2
Q14764Major vault protein MVP88.9−1.2−1.9−2.6
P35579Myosin-9 MYH9166.85−1.7−2.9−1.6
P61160Actin-related protein 2 ACTR260.97−2.6−1.4−1.8
Q02338D-beta-hydroxybutyrate dehydrogenase, mitochondrial BDH171.17−1.7−1.4−2.6
P20290Transcription factor BTF3 BTF362.091.2−2.9−1.5
Q7L1Q6Basic leucine zipper and W2 domain-containing protein 1BZW195.07−1.3−3.8−2.3
P17655Calpain-2 catalytic subunit CAPN2101.09−1.3−3.1−1.9
A8K010Putative uncharacterized protein C6orf176 C6orf17633.72−1.1−3.3−1.6
P53618Coatomer subunit beta COPB164.4−2.9−1.5−1.4
O60888Protein CutA CUTA62.421.2−2.2−2.9
O14976Cyclin-G-associated kinase GAK47.67−3.4−3.8−2.9
O95395Beta-1,3-galactosyl-O-glycosyl-glycoprotein beta-1,6-N-acetylglucosaminyltransferase 3 GCNT359.62−3.2−1.2−1.6
P52306Rap1 GTPase-GDP dissociation stimulator 1 RAP1GDS166.74−1.2−3−1.1
O15228Dihydroxyacetone phosphate acyltransferase GNPAT46.58−4.5−3.1−1.9
Q86Y56HEAT repeat-containing protein 2 HEATR263.07−1.5−5−1.4
Q6YP21 Kynurenine-oxoglutarate transaminase 3 CCBL239.91.7OFF1.5
P43243Matrin-3 MATR3269.161−1.2−4.5
O75439Mitochondrial-processing peptidase subunit beta PMPCB99.73−2.91.2−1.1
O43847Nardilysin NRD141.16−1.3−2.91.6
O15270Serine palmitoyltransferase 2SPTLC262.71−3.2−2−2.5
P29122Proprotein convertase subtilisin/kexin type 6 PCSK630.52−1.3−1.5−16.7
P00491Purine nucleoside phosphorylase PNP91.96−2−5.61.7
P25787Proteasome subunit alpha type 2 PSMA275.421.6−2.91
O14818Proteasome subunit alpha type 7 PSMA7103.36−1.3−2.8−1.1
Q13813Spectrin alpha chain, nonerythrocytic 1SPTAN1112.76−1.6−3.1−1.6
Q9UL25Ras-related protein Rab-21 RAB21129.33−2.61.1−1.1
Q9NV58E3 ubiquitin-protein ligase RNF19A RNF19A31.17−1.4−1.8−4.8
P6208140S ribosomal protein S7 RPS765.3−4.2−3.8−2.6
Q96FQ6Protein S100-A16 S100A1678−2.6−1.41
P63208S-phase kinase-associated protein 1 SKP172.23−7.1−3.11.2
O00193Small acidic protein SMAP74.25−1.5−2.6−1.3
Q96JX3Protein SERAC1 SERAC132.13−1.11.3−2.5
Q9Y5M8Signal recognition particle receptor subunit beta SRPRB54.95−1.7−1.1−2.9
Q96I99Succinyl-CoA ligase [GDP-forming] subunit beta, mitochondrial SUCLG285.7−1.5−2.4−2.4

Upregulated proteins
O75223Gamma-glutamylcyclotransferaseGGCT37.331.123.8
P61758Prefoldin subunit 3VBP1152.642.12.22.5
Q9BWD1Acetyl-CoA acetyltransferase, cytosolicACAT269.22.523.2
Q86TX2Acyl-coenzyme A thioesterase 1ACOT173.322.81.62.4
P02647Apolipoprotein A-IAPOA176.161.92.7−1.1
O00192Armadillo repeat protein deleted in velocardiofacial syndromeARVCF46ONONON
Q9NYF8Bcl-2-associated transcription factor 1BCLAF186.181ONON
O14523C2 domain-containing protein 2-likeC2CD2L31.32.91.72
P60953Cell division control protein 42 homologCDC4252.582.31.53.1
Q9UJS0Calcium-binding mitochondrial carrier protein Aralar2SLC25A1357.952.43.41.9
O43809Cleavage and polyadenylation specificity factor subunit 5NUDT2170.351.2−1.72.5
Q96C86Scavenger mRNA-decapping enzyme DcpSDCPS70.591.52.12.4
Q5TDH0Protein DDI1 homolog 2DDI231.551ONON
P18858DNA ligase 1LIG149.532.32.22.6
Q5JZY3Ephrin type A receptor 10EPHA1037.22ONONON
O76003Glutaredoxin-3GLRX350.313.21.82.4
Q14103Heterogeneous nuclear ribonucleoprotein D0HNRNPD138.391.431.4
Q96P70Importin-9IPO988.043.92.92.9
P33176Kinesin-1 heavy chainKIF5B102.254.33.65.1
P49006MARCKS-related proteinMARCKSL1112.06ONON1
Q9H3R2Mucin-13MUC1368.361.84.21.9
Q7Z406Myosin-14MYH1490.59−1.19.7−1.9
Q9NX40OCIA domain-containing protein 1OCIAD187.82.62.92.3
Q15102Platelet-activating factor acetylhydrolase IB subunit gammaPAFAH1B384.52.462.41.6
O43252Bifunctional 3′-phosphoadenosine 5′-phosphosulfate synthetase 1PAPSS151.911ON
P36871Phosphoglucomutase-1PGM1105.243.3−2.22.2
P35232ProhibitinPHB90.528.44.73.2
Q9UHJ6SedoheptulokinaseSHPK59.091 1ON
P52788Spermine synthaseSMS60.56−1.22.13.2
Q01082Spectrin beta chain, brain 1SPTBN163.84ONONON
Q9Y3F4Serine-threonine kinase receptor-associated proteinSTRAP80.32.93.74.5
P26640Valyl-tRNA synthetaseVARS54.771.733.5
Q86VP1Tax1-binding protein 1TAX1BP131.292.311.71.2
Q9HD45Transmembrane 9 superfamily member 3TM9SF371.541.53.1−1.8
P55327Tumor protein D52TPD5264.695.91.81.8
P38606V-type proton ATPase catalytic subunit AATP6V1A66.31−1.24.42.8
P45880Voltage-dependent anion-selective channel protein 2VDAC259.142.142

UniProtKB/Swiss-Prot entry accession number (http://www.uniprot.org/). The UniProt accession numbers highlighted in bold indicate those proteins (8) that were also identified by chemical proteomics.
Recommended gene name (official gene symbol) as provided by UniProtKB/Swiss-Prot.
MASCOT protein score.
Fold change of differentially regulated proteins after 40 min, 2 h, and 5 h of treatment with CIGB-552 [(−) downregulated]. For each incubation time with CIGB-552, is given the fold change threshold to consider a protein significantly up- or downregulated. The fold changes that have an absolute value higher than the threshold are highlighted in boldface.

Bioinformatics analysis of comparative proteomic data demonstrated the enrichment of biological processes such as gene expression, proteolysis, and response to drug (Figure 5). Importantly, carbohydrate and nitrogen metabolic processes, protein transport, cell cycle, and regulation of apoptosis were found to be significantly represented in chemical and comparative proteomic profiles (Figures 2 and 5). These results demonstrated that even in different cell lines (HT-29, Hep-2) the CIGB-552 antitumor effect is exerted by modulating similar biological processes.

Differentially expressed proteins related to apoptosis and cell cycle regulation support the CIGB-552 antitumor effect. Among them, the DNA replication licensing factor MCM6 (MCM6), the cyclin-G-associated kinase (GAK), and the S-phase kinase-associated protein 1 (SKP1), which function as positive regulators of cell cycle [25, 45, 46], were significantly downregulated after 40 min of CIGB-552 treatment.

As expected, on CIGB-552-treated cells, the apoptosis was a much later event than the negative regulation of cell cycle (Figure 6). Only after 2 h, the CIGB-552 treatment increased the expression of proapoptotic proteins. For instance, the voltage-dependent anion-selective channel protein 2 (VDAC2) was significantly upregulated after 2 h of CIGB-552 treatment. This protein is a component of the mitochondrial permeability transition pore complex, which releases apoptogenic factors, such as cytochrome C, during apoptosis [47, 48]. In line with this evidence, the expression of the cytochrome C-releasing factor 21 (GGCT) was increased sequentially on CIGB-552-treated cells, the upregulation being significant after 5 h of treatment. Furthermore, the expression of Bcl-2-associated transcription factor 1 (BCLAF1), a transcriptional repressor that promotes apoptosis [49, 50], was increased after 2 h of treatment with CIGB-552. Concomitant with upregulation of proapoptotic proteins, the expression of the major vault protein (MVP) was significantly decreased on HT-29 cells after 5 h of CIGB-552 treatment (Figure 5). MVP functions as a multidrug resistance factor [5153]. Therefore, the downregulation of MVP could facilitate the proapoptotic effect in response to CIGB-552 treatment and overcome the chemotherapeutic resistance usually developed by cancer cells.

Most of the proteins related to proteolysis were downregulated by CIGB-552 treatment (Figure 5). Two peptidase M16 family members, which have metalloendopeptidase activity, were identified: the mitochondrial-processing peptidase subunit beta (PMPCB) and nardilysin (NRD1). Particularly, NRD1 protein promotes cell growth and invasion of cancer cells [54, 55]. In the present work, Calpain-2 (CAPN2) and proprotein convertase subtilisin/kexin type 6 (PCSK6) were also identified. Both proteins are proteases related to tumor progression. CAPN2 upregulation increases the invasive potential of different tumor cells; such effect has been associated with secretion of matrix metalloproteinases (MMP-2 and MMP-9) [5658]. Besides, in the tumor microenvironment, the inhibition of Calpain blocks angiogenesis [59]. PCSK6 substrates include precursors of cancer related proteins such as matrix metalloproteinases and the vascular endothelial growth factor (VEGF) [60]. Therefore, upregulation of PCSK6 increases the invasiveness of carcinoma cells [6163]. Indeed inhibition of PCSK6 has been proposed as an anticancer therapeutic strategy [64, 65]. Altogether, the concomitant downregulation of NRD1, CAPN2, and PCSK6 suggests that CIGB-552 treatment could inhibit cancer cell invasion, exerting a potential antimetastatic effect. As previously mentioned, the antimetastatic effect of CIGB-552 is also supported by chemical proteomics results. To validate such hypothesis other studies have been conducted.

To analyze chemical proteomic data in combination with differentially expressed proteins, an interaction network was represented (Figure 7). According to interaction databases, 18 differentially modulated proteins interact with the potential CIGB-552 target profile (18 proteins). Prohibitin (PHB) interacts with three proteins included in the CIGB-552 target profile (ACTN1, SET, and C1QBP) (Figure 7). The expression of PHB, which is a negative regulator of cell proliferation [66], was increased on CIGB-552-treated cells. Besides, prefoldin subunit 3 (PFDN3), upregulated by CIGB-552 treatment, interacts with the acidic leucine-rich nuclear phosphoprotein 32 family member A (ANP32A) (Figure 7). It is a tumor suppressor protein [67] that was identified as a potential CIGB-552 target. Additionally, PFDN3 binds to von Hippel-Lindau protein (VHL), which regulates ubiquitination and proteasomal degradation of the hypoxia-inducible factor 1-alpha (HIF-α) [68]. Functional subnetworks of “negative regulation of cell cycle” and “response to hypoxia” were perturbed by CIGB-552 target profile (Table 2). Accordingly, such effects are supported by the upregulation of PHB and PFDN3 on HT-29 cells treated with CIGB-552.

The intersection between chemical proteomic data and differentially expressed proteins indicates potential mediators of CIGB-552 antitumor effect. The Rho GTPase-activating protein 1 (ARHGAP1) was identified as a potential CIGB-552 target (Figure 7). This protein inactivates signal transduction mediated by Rho-family GTPases, such as RhoA, Rac1, and CDC42 [69, 70]. Accordingly, the functional category of “small GTPase mediated signal transduction” was significantly represented in the proteomic profile modulated by CIGB-552 on HT-29 cells (Figure 5). In addition, ARHGAP1 interacts with Tax1-binding protein 1 (TAX1BP1) (Figure 7) [71]. The expression of TAX1BP1 was increased at 2 h of CIGB-552 treatment (Table 3). TAX1BP1, in conjugation with the zinc finger protein A20, the E3 ubiquitin-protein ligase Itch, and RING finger protein 11, constitutes the A20 ubiquitin-editing protein complex [72]. Downstream to Toll-like receptor 4 (TLR4) and tumor necrosis factor receptor 1 (TNFR1) the A20 complex inactivates TRAF6 and RIP1, repressing the NF-κB signaling pathway [73, 74]. Furthermore, in the present work, several potential CIGB-552 targets interact with TRAF6 (Figure 4). Altogether, these results suggest that CIGB-552 could inhibit the TRAF6/NF-κB axis.

Other proteins related to NF-κB signaling pathway were also modulated on HT-29 cells by CIGB-552. The S-phase kinase-associated protein 1 (SKP1) was downregulated at 40 min and 2 h after CIGB-552 treatment (Figure 7). This protein is an essential component of the SCF (SKP1-CUL1-F-box protein) E3 ubiquitin ligase complex, which mediates the ubiquitination of proteins involved in cell cycle progression, signal transduction, and transcription [75]. As part of SCF complex, the F-box/WD repeat-containing protein 1A (BTRC) ubiquitinates IκB leading to its proteasomal degradation and consequently NF-κB activation [76]. In the present work, the F-box protein BTRC was not identified. However, by decreasing the expression of SKP1, CIGB-552 could compromise the function of SCF complex. SKP1 is a proteasome interacting protein (Figure 7) [77]. The expression of two proteasome subunits (PSMA2, PSMA7) was decreased at 2 h of CIGB-552 treatment (Figure 7). Besides, several components of the proteasome complex, including the subunit PSMA2, were identified by chemical proteomics as potential CIGB-552 targets (Figure 3). Interestingly, modulation of SKP1 and proteasome subunits is concomitant with TAX1BP1 upregulation (Figure 6). These results suggest that CIGB-552 potentiates, at 2 h following treatment, the inhibition of NF-κB signaling pathway. Such effect is also supported by chemical proteomic results, as the “positive regulation of NF-κB transcription factor activity” is a functional subnetwork disrupted by the CIGB-552 target profile (Figure 4).

A potential mechanism of NF-κB inhibition by CIGB-552 was proposed based on both chemical and comparative proteomic profiles (Figure 8). Previous results demonstrated that CIGB-552 interacts with COMMD1 triggering the RelA ubiquitination [5]. Other potential CIGB-552 interactors include the proteins ETHE1 and SETD7. ETHE1 and SETD7, as well as COMMD1, regulate the function of RelA [4042]. Therefore, CIGB-552 by interacting with ETHE1 and SETD7 could repress the transcription of NF-κB target genes (Figure 8). As the MAPK cascade increases the transcriptional activity of RelA [39], interaction between CIGB-552 and MAPKs proteins (MAP2K6, MAPK14) could also inhibit the NF-κB signaling pathway (Figure 8). Finally, the upregulation of TAX1BP1 and the downregulation of SKP1 induced by CIGB-552 are of great relevance. From a mechanistic view, the modulation of both proteins by CIGB-552 supports the inhibition of NF-κB signaling pathway (Figure 8).

4. Conclusions

This study explored at the proteomic level the CIGB-552 antitumor mechanism of action. The CIGB-552 target profile was investigated in Hep-2 cells using a chemical proteomic approach. In total, 161 proteins were identified as potential CIGB-552 targets. This target profile was reduced to 55 proteins by a target deconvolution strategy. According to chemical proteomic results, CIGB-552 could be a multitarget drug. Downstream to drug-target interactions proteins that mediate cellular response to CIGB-552 treatment were identified by a comparative proteomic experiment including time series (40 min, 2 h, and 5 h). The proteomic profile modulated by CIGB-552 on HT-29 cells includes 72 differentially expressed proteins. The results of both experimental approaches, chemical and expression-based proteomics, were highly consistent. Proteins related to NF-κB signal transduction were identified as potential CIGB-552 targets and were significantly modulated by CIGB-552 treatment. Such findings are in line with previous results demonstrating that CIGB-552 blocks NF-κB signaling by upregulation of COMMD1 and consequently ubiquitination of RelA. The proteomic data revealed new mediators for NF-κB inhibition in response to CIGB-552 treatment. Proteins related to cell proliferation and apoptosis were identified by chemical and expression-based proteomics. These differentially expressed proteins could represent subrogate biomarkers of the CIGB-552 effect on tumor cells. The negative regulation of cell cycle and promotion of apoptosis seem to be early and late events, respectively, triggered by treatment with CIGB-552. Proteins related to cell invasion were differentially modulated on CIGB-552-treated cells. This biological process was also represented in the functional subnetworks perturbed by the potential CIGB-552 targets. Therefore, an antimetastatic effect could be expected as a consequence of CIGB-552 treatment. Further experiments are required to corroborate the role of identified proteins on mediating the anticancer effect of CIGB-552. Overall, the current study contributes to a better understanding of the CIGB-552 potentialities for cancer therapy.

Conflict of Interests

No potential conflict of interests was disclosed.

Acknowledgments

This work was supported by Heber Biotec SA. The Proteomic and Bioinformatics Departments at the Center for Genetic Engineering and Biotechnology would like to thank INSPUR Company (China) for their kind donation of a computer cluster used for data storage and processing.

Supplementary Materials

Supplementary Table S1a contains the proteins purified on the CIGB-552 streptavidin sepharose matrix from the Hep-2 cell lysate.

Supplementary Table S1b contains the CRAPome information about the central proteome proteins which were identified as part of the CIGB-552 target profile.

Supplementary Table S2a contains biological processes enriched in the CIGB-552 target profile after target deconvolution.

Supplementary Table S2b contains cellular components enriched in the CIGB-552 target profile.

Supplementary Table S2c contains cellular components enriched in the CIGB-552 target profile after target deconvolution.

Supplementary Table S3 contains the proteomic profile modulated by CIGB-552 in HT-29 cells.

  1. Supplementary Material

References

  1. F. Albericio and H. G. Kruger, “Therapeutic peptides,” Future Medicinal Chemistry, vol. 4, no. 12, pp. 1527–1531, 2012. View at: Publisher Site | Google Scholar
  2. D. Wu, Y. Gao, Y. Qi, L. Chen, Y. Ma, and Y. Li, “Peptide-based cancer therapy: opportunity and challenge,” Cancer Letters, vol. 351, no. 1, pp. 13–22, 2014. View at: Publisher Site | Google Scholar
  3. J. Thundimadathil, “Cancer treatment using peptides: current therapies and future prospects,” Journal of Amino Acids, vol. 2012, Article ID 967347, 13 pages, 2012. View at: Publisher Site | Google Scholar
  4. M. G. Vallespi, J. R. Fernandez, I. Torrens et al., “Identification of a novel antitumor peptide based on the screening of an Ala-library derived from the LALF32-51 region,” Journal of Peptide Science, vol. 16, no. 1, pp. 40–47, 2010. View at: Publisher Site | Google Scholar
  5. J. R. F. Massó, B. P. Argüelles, Y. Tejeda et al., “The antitumor peptide CIGB-552 increases COMMD1 and inhibits growth of human lung cancer cells,” Journal of Amino Acids, vol. 2013, Article ID 251398, 13 pages, 2013. View at: Publisher Site | Google Scholar
  6. M. G. Vallespí, G. Pimentel, A. Cabrales-Rico et al., “Antitumor efficacy, pharmacokinetic and biodistribution studies of the anticancer peptide CIGB-552 in mouse models,” Journal of Peptide Science, vol. 20, no. 11, pp. 850–859, 2014. View at: Publisher Site | Google Scholar
  7. G. B. Fields and R. L. Noble, “Solid phase peptide synthesis utilizing 9-fluorenylmethoxycarbonyl amino acids,” International Journal of Peptide and Protein Research, vol. 35, no. 3, pp. 161–214, 1990. View at: Google Scholar
  8. J. Heukeshoven and R. Dernick, “Characterization of a solvent system for separation of water-insoluble poliovirus proteins by reversed-phase high-performance liquid chromatography,” Journal of Chromatography A, vol. 326, pp. 91–101, 1985. View at: Publisher Site | Google Scholar
  9. D. N. Perkins, D. J. C. Pappin, D. M. Creasy, and J. S. Cottrell, “Probability-based protein identification by searching sequence databases using mass spectrometry data,” Electrophoresis, vol. 20, no. 18, pp. 3551–3567, 1999. View at: Publisher Site | Google Scholar
  10. D. Tabas-Madrid, R. Nogales-Cadenas, and A. Pascual-Montano, “GeneCodis3: a non-redundant and modular enrichment analysis tool for functional genomics,” Nucleic Acids Research, vol. 40, no. 1, pp. W478–W483, 2012. View at: Publisher Site | Google Scholar
  11. A. Ruepp, B. Brauner, I. Dunger-Kaltenbach et al., “CORUM: the comprehensive resource of mammalian protein complexes,” Nucleic Acids Research, vol. 36, no. 1, pp. D646–D650, 2008. View at: Publisher Site | Google Scholar
  12. T. R. Burkard, M. Planyavsky, I. Kaupe et al., “Initial characterization of the human central proteome,” BMC Systems Biology, vol. 5, article 17, 2011. View at: Publisher Site | Google Scholar
  13. D. Mellacheruvu, Z. Wright, A. L. Couzens et al., “The CRAPome: a contaminant repository for affinity purification-mass spectrometry data,” Nature Methods, vol. 10, no. 8, pp. 730–736, 2013. View at: Publisher Site | Google Scholar
  14. T. R. Burkard, U. Rix, F. P. Breitwieser, G. Superti-Furga, and J. Colinge, “A computational approach to analyze the mechanism of action of the kinase inhibitor bafetinib,” PLoS Computational Biology, vol. 6, no. 11, Article ID e1001001, 2010. View at: Publisher Site | Google Scholar
  15. Y. Ishihama, Y. Oda, T. Tabata et al., “Exponentially modified protein abundance index (emPAI) for estimation of absolute protein amount in proteomics by the number of sequenced peptides per protein,” Molecular and Cellular Proteomics, vol. 4, no. 9, pp. 1265–1272, 2005. View at: Publisher Site | Google Scholar
  16. C. S. Spahr, M. T. Davis, M. D. McGinley et al., “Towards defining the urinary proteome using liquid chromatography-tandem mass spectrometry. I. Profiling an unfractionated tryptic digest,” Proteomics, vol. 1, no. 1, pp. 93–107, 2001. View at: Publisher Site | Google Scholar
  17. J. E. Elias and S. P. Gygi, “Target-decoy search strategy for increased confidence in large-scale protein identifications by mass spectrometry,” Nature Methods, vol. 4, no. 3, pp. 207–214, 2007. View at: Publisher Site | Google Scholar
  18. J. Fernandez-de-Cossio, L. J. Gonzalez, Y. Satomi et al., “Isotopica: a tool for the calculation and viewing of complex isotopic envelopes,” Nucleic Acids Research, vol. 32, supplement 2, pp. W674–W678, 2004. View at: Publisher Site | Google Scholar
  19. J. Fernández-de-Cossio, L. J. Gonzalez, Y. Satomi et al., “Automated interpretation of mass spectra of complex mixtures by matching of isotope peak distributions,” Rapid Communications in Mass Spectrometry, vol. 18, no. 20, pp. 2465–2472, 2004. View at: Publisher Site | Google Scholar
  20. P. Shannon, A. Markiel, O. Ozier et al., “Cytoscape: a software environment for integrated models of biomolecular interaction networks,” Genome Research, vol. 13, no. 11, pp. 2498–2504, 2003. View at: Publisher Site | Google Scholar
  21. A. Martin, M. E. Ochagavia, L. C. Rabasa, J. Miranda, J. Fernandez-de-Cossio, and R. Bringas, “BisoGenet: a new tool for gene network building, visualization and analysis,” BMC Bioinformatics, vol. 11, article 91, 2010. View at: Publisher Site | Google Scholar
  22. G. Warsow, B. Greber, S. S. I. Falk et al., “ExprEssence—revealing the essence of differential experimental data in the context of an interaction/regulation net-work,” BMC Systems Biology, vol. 4, article 164, 2010. View at: Publisher Site | Google Scholar
  23. P. Chatrath, I. S. Scott, L. S. Morris et al., “Aberrant expression of minichromosome maintenance protein-2 and Ki67 in laryngeal squamous epithelial lesions,” British Journal of Cancer, vol. 89, no. 6, pp. 1048–1054, 2003. View at: Publisher Site | Google Scholar
  24. J. Kikuchi, I. Kinoshita, Y. Shimizu et al., “Minichromosome maintenance (MCM) protein 4 as a marker for proliferation and its clinical and clinicopathological significance in non-small cell lung cancer,” Lung Cancer, vol. 72, no. 2, pp. 229–237, 2011. View at: Publisher Site | Google Scholar
  25. M. Lei, “The MCM complex: its role in DNA replication and implications for cancer therapy,” Current Cancer Drug Targets, vol. 5, no. 5, pp. 365–380, 2005. View at: Publisher Site | Google Scholar
  26. O. Saydam, O. Senol, T. B. M. Schaaij-Visser et al., “Comparative protein profiling reveals minichromosome maintenance (MCM) proteins as novel potential tumor markers for meningiomas,” Journal of Proteome Research, vol. 9, no. 1, pp. 485–494, 2010. View at: Publisher Site | Google Scholar
  27. U. Rix and G. Superti-Furga, “Target profiling of small molecules by chemical proteomics,” Nature Chemical Biology, vol. 5, no. 9, pp. 616–624, 2009. View at: Publisher Site | Google Scholar
  28. M. Raida, “Drug target deconvolution by chemical proteomics,” Current Opinion in Chemical Biology, vol. 15, no. 4, pp. 570–575, 2011. View at: Publisher Site | Google Scholar
  29. B. van de Sluis, X. Mao, Y. Zhai et al., “COMMD1 disrupts HIF-1α/β dimerization and inhibits human tumor cell invasion,” Journal of Clinical Investigation, vol. 120, no. 6, pp. 2119–2130, 2010. View at: Publisher Site | Google Scholar
  30. U. Rix, L. L. R. Rix, A. S. Terker et al., “A comprehensive target selectivity survey of the BCR-ABL kinase inhibitor INNO-406 by kinase profiling and chemical proteomics in chronic myeloid leukemia cells,” Leukemia, vol. 24, no. 1, pp. 44–50, 2010. View at: Publisher Site | Google Scholar
  31. D. M. Gilkes, G. L. Semenza, and D. Wirtz, “Hypoxia and the extracellular matrix: drivers of tumour metastasis,” Nature Reviews Cancer, vol. 14, no. 6, pp. 430–439, 2014. View at: Publisher Site | Google Scholar
  32. N. D. Perkins and T. D. Gilmore, “Good cop, bad cop: the different faces of NF-κB,” Cell Death and Differentiation, vol. 13, no. 5, pp. 759–772, 2006. View at: Publisher Site | Google Scholar
  33. A. S. Baldwin, “Control of oncogenesis and cancer therapy resistance by the transcription factor NF-kappaB,” The Journal of Clinical Investigation, vol. 107, no. 3, pp. 241–246, 2001. View at: Publisher Site | Google Scholar
  34. B. Rayet and C. Gélinas, “Aberrant rel/nfkb genes and activity in human cancer,” Oncogene, vol. 18, no. 49, pp. 6938–6947, 1999. View at: Publisher Site | Google Scholar
  35. S. Liptay, C. K. Weber, L. Ludwig, M. Wagner, G. Adler, and R. M. Schmid, “Mitogenic and antiapoptotic role of constitutive NF-κB/Rel activity in pancreatic cancer,” International Journal of Cancer, vol. 105, no. 6, pp. 735–746, 2003. View at: Publisher Site | Google Scholar
  36. C. Huang, K. Huang, C. Wang et al., “Overexpression of mitogen-activated protein kinase kinase 4 and nuclear factor-kappaB in laryngeal squamous cell carcinoma: a potential indicator for poor prognosis,” Oncology Reports, vol. 22, no. 1, pp. 89–95, 2009. View at: Publisher Site | Google Scholar
  37. B. Lamothe, A. Besse, A. D. Campos, W. K. Webster, H. Wu, and B. G. Darnay, “Site-specific Lys-63-linked tumor necrosis factor receptor-associated factor 6 auto-ubiquitination is a critical determinant of IκB kinase activation,” The Journal of Biological Chemistry, vol. 282, no. 6, pp. 4102–4112, 2007. View at: Publisher Site | Google Scholar
  38. M. Windheim, M. Stafford, M. Peggie, and P. Cohen, “Interleukin-1 (IL-1) induces the Lys63-linked polyubiquitination of IL-1 receptor-associated kinase 1 to facilitate NEMO binding and the activation of IκBα kinase,” Molecular and Cellular Biology, vol. 28, no. 5, pp. 1783–1791, 2008. View at: Publisher Site | Google Scholar
  39. L. Vermeulen, G. De Wilde, P. Van Damme, W. V. Berghe, and G. Haegeman, “Transcriptional activation of the NF-κB p65 subunit by mitogen- and stress-activated protein kinase-1 (MSK1),” The EMBO Journal, vol. 22, no. 6, pp. 1313–1324, 2003. View at: Publisher Site | Google Scholar
  40. H. Higashitsuji, H. Higashitsuji, T. Nagao et al., “A novel protein overexpressed in hepatoma accelerates export of NF-kappa B from the nucleus and inhibits p53-dependent apoptosis,” Cancer Cell, vol. 2, no. 4, pp. 335–346, 2002. View at: Publisher Site | Google Scholar
  41. Y. Li, M. A. Reddy, F. Miao et al., “Role of the histone H3 lysine 4 methyltransferase, SET7/9, in the regulation of NF-κB-dependent inflammatory genes. Relevance to diabetes and inflammation,” Journal of Biological Chemistry, vol. 283, no. 39, pp. 26771–26781, 2008. View at: Publisher Site | Google Scholar
  42. X.-D. Yang, B. Huang, M. Li, A. Lamb, N. L. Kelleher, and L.-F. Chen, “Negative regulation of NF-κB action by Set9-mediated lysine methylation of the relA subunit,” The EMBO Journal, vol. 28, no. 8, pp. 1055–1066, 2009. View at: Publisher Site | Google Scholar
  43. T. A. Libermann and D. Baltimore, “Activation of interleukin-6 gene expression through the NF-κB transcription factor,” Molecular and Cellular Biology, vol. 10, no. 5, pp. 2327–2334, 1990. View at: Publisher Site | Google Scholar
  44. B. Hoesel and J. A. Schmid, “The complexity of NF-kappaB signaling in inflammation and cancer,” Molecular Cancer, vol. 12, article 86, 2013. View at: Publisher Site | Google Scholar
  45. M. E. Tanenbaum, T. Vallenius, E. F. Geers, L. Greene, T. P. Mäkelä, and R. H. Medema, “Cyclin G-associated kinase promotes microtubule outgrowth from chromosomes during spindle assembly,” Chromosoma, vol. 119, no. 4, pp. 415–424, 2010. View at: Publisher Site | Google Scholar
  46. D. Frescas and M. Pagano, “Deregulated proteolysis by the F-box proteins SKP2 and β-TrCP: tipping the scales of cancer,” Nature Reviews Cancer, vol. 8, no. 6, pp. 438–449, 2008. View at: Publisher Site | Google Scholar
  47. S. Shimizu, M. Narita, and Y. Tsujimoto, “Bcl-2 family proteins regulate the release of apoptogenic cytochrome c by the mitochondrial channel VDAC,” Nature, vol. 399, no. 6735, pp. 483–487, 1999. View at: Publisher Site | Google Scholar
  48. Y. Tsujimoto and S. Shimizu, “The voltage-dependent anion channel: an essential player in apoptosis,” Biochimie, vol. 84, no. 2-3, pp. 187–193, 2002. View at: Publisher Site | Google Scholar
  49. G. M. Kasof, L. Goyal, and E. White, “Btf, a novel death-promoting transcriptional repressor that interacts with Bcl-2-related proteins,” Molecular and Cellular Biology, vol. 19, no. 6, pp. 4390–4404, 1999. View at: Google Scholar
  50. H. Liu, Z.-G. Lu, Y. Miki, and K. Yoshida, “Protein kinase C delta induces transcription of the TP53 tumor suppressor gene by controlling death-promoting factor Btf in the apoptotic response to DNA damage,” Molecular and Cellular Biology, vol. 27, no. 24, pp. 8480–8491, 2007. View at: Publisher Site | Google Scholar
  51. C. M. Laurençot, G. L. Scheffer, R. J. Scheper, and R. H. Shoemaker, “Increased LRP mRNA expression is associated with the MDR phenotype in intrinsically resistant human cancer cell lines,” International Journal of Cancer, vol. 72, no. 6, pp. 1021–1026, 1997. View at: Publisher Site | Google Scholar
  52. S. Tang, C. Bai, P. Yang, and X. Chen, “14-3-3ε boosts bleomycin-induced DNA damage response by inhibiting the drug-resistant activity of MVP,” Journal of Proteome Research, vol. 12, no. 6, pp. 2511–2524, 2013. View at: Publisher Site | Google Scholar
  53. A. Losert, D. Lötsch, A. Lackner et al., “The major vault protein mediates resistance to epidermal growth factor receptor inhibition in human hepatoma cells,” Cancer Letters, vol. 319, no. 2, pp. 164–172, 2012. View at: Publisher Site | Google Scholar
  54. K. Kanda, H. Komekado, T. Sawabu et al., “Nardilysin and ADAM proteases promote gastric cancer cell growth by activating intrinsic cytokine signalling via enhanced ectodomain shedding of TNF-α,” EMBO Molecular Medicine, vol. 4, no. 5, pp. 396–411, 2012. View at: Publisher Site | Google Scholar
  55. N. Uraoka, N. Oue, N. Sakamoto et al., “NRD1, which encodes nardilysin protein, promotes esophageal cancer cell invasion through induction of MMP2 and MMP3 expression,” Cancer Science, vol. 105, no. 1, pp. 134–140, 2014. View at: Publisher Site | Google Scholar
  56. H. Roumes, L. Leloup, E. Dargelos, J.-J. Brustis, L. Daury, and P. Cottin, “Calpains: markers of tumor aggressiveness?” Experimental Cell Research, vol. 316, no. 9, pp. 1587–1599, 2010. View at: Publisher Site | Google Scholar
  57. H. S. Jang, S. Lal, and J. A. Greenwood, “Calpain 2 is required for glioblastoma cell invasion: rRegulation of matrix metalloproteinase 2,” Neurochemical Research, vol. 35, no. 11, pp. 1796–1804, 2010. View at: Publisher Site | Google Scholar
  58. B. Chen, J. Tang, Y.-S. Guo, Y. Li, Z.-N. Chen, and J.-L. Jiang, “Calpains are required for invasive and metastatic potentials of human HCC cells,” Cell Biology International, vol. 37, no. 7, pp. 643–652, 2013. View at: Publisher Site | Google Scholar
  59. Y. Su, Z. Cui, Z. Li, and E. R. Block, “Calpain-2 regulation of VEGF-mediated angiogenesis,” The FASEB Journal, vol. 20, no. 9, pp. 1443–1451, 2006. View at: Publisher Site | Google Scholar
  60. D. E. Bassi, J. Fu, R. L. de Cicco, and A. J. P. Klein-Szanto, “Proprotein convertases: “Master switches” in the regulation of tumor growth and progression,” Molecular Carcinogenesis, vol. 44, no. 3, pp. 151–161, 2005. View at: Publisher Site | Google Scholar
  61. M. Lapierre, G. Siegfried, N. Scamuffa et al., “Opposing function of the proprotein convertases furin and PACE4 on breast cancer cells' malignant phenotypes: role of tissue inhibitors of metalloproteinase-1,” Cancer Research, vol. 67, no. 19, pp. 9030–9034, 2007. View at: Publisher Site | Google Scholar
  62. S. Delic, N. Lottmann, K. Jetschke, G. Reifenberger, and M. J. Riemenschneider, “Identification and functional validation of CDH11, PCSK6 and SH3GL3 as novel glioma invasion-associated candidate genes,” Neuropathology and Applied Neurobiology, vol. 38, no. 2, pp. 201–212, 2012. View at: Publisher Site | Google Scholar
  63. D. E. Bassi, J. Cenna, J. Zhang, E. Cukierman, and A. J. Klein-Szanto, “Enhanced aggressiveness of benzopyrene-induced squamous carcinomas in transgenic mice overexpressing the proprotein convertase PACE4 (PCSK6),” Molecular Carcinogenesis, vol. 54, no. 10, pp. 1122–1131, 2015. View at: Publisher Site | Google Scholar
  64. F. Couture, C. Levesque, V. Dumulon-Perreault et al., “PACE4-based molecular targeting of prostate cancer using an engineered 64Cu-radiolabeled peptide inhibitor,” Neoplasia, vol. 16, no. 8, pp. 634–643, 2014. View at: Publisher Site | Google Scholar
  65. D. E. Bassi, J. Zhang, J. Cenna, S. Litwin, E. Cukierman, and A. J. P. Klein-Szanto, “Proprotein convertase inhibition results in decreased skin cell proliferation, tumorigenesis, and metastasis,” Neoplasia, vol. 12, no. 7, pp. 516–526, 2010. View at: Publisher Site | Google Scholar
  66. S. Wang, B. Zhang, and D. V. Faller, “Prohibitin requires Brg-1 and Brm for the repression of E2F and cell growth,” The EMBO Journal, vol. 21, no. 12, pp. 3019–3028, 2002. View at: Publisher Site | Google Scholar
  67. W. Pan, L. S. de Graca, Y. Shao, Q. Yin, H. Wu, and X. Jiang, “PHAPI/pp32 suppresses tumorigenesis by stimulating apoptosis,” The Journal of Biological Chemistry, vol. 284, no. 11, pp. 6946–6954, 2009. View at: Publisher Site | Google Scholar
  68. H. Tsuchiya, T. Iseda, and O. Hino, “Identification of a novel protein (VBP-1) binding to the von Hippel-Lindau (VHL) tumor suppressor gene product,” Cancer Research, vol. 56, no. 13, pp. 2881–2885, 1996. View at: Google Scholar
  69. C. A. Lancaster, P. M. Taylor-Harris, A. J. Self, S. Brill, H. E. Van Erp, and A. Hall, “Characterization of rhoGAP: a GTPase-activating protein for rho-related small GTPases,” The Journal of Biological Chemistry, vol. 269, no. 2, pp. 1137–1142, 1994. View at: Google Scholar
  70. M. Lazarini, F. Traina, J. A. Machado-Neto et al., “ARHGAP21 is a RhoGAP for RhoA and RhoC with a role in proliferation and migration of prostate adenocarcinoma cells,” Biochimica et Biophysica Acta - Molecular Basis of Disease, vol. 1832, no. 2, pp. 365–374, 2013. View at: Publisher Site | Google Scholar
  71. G. M. Nagaraja and R. P. Kandpal, “Chromosome 13q12 encoded Rho GTPase activating protein suppresses growth of breast carcinoma cells, and yeast two-hybrid screen shows its interaction with several proteins,” Biochemical and Biophysical Research Communications, vol. 313, no. 3, pp. 654–665, 2004. View at: Publisher Site | Google Scholar
  72. E. Jacque and S. C. Ley, “RNF11, a new piece in the A20 puzzle,” The EMBO Journal, vol. 28, no. 5, pp. 455–456, 2009. View at: Publisher Site | Google Scholar
  73. N. Shembade, A. Ma, and E. W. Harhaj, “Inhibition of NF-κB signaling by A20 through disruption of ubiquitin enzyme complexes,” Science, vol. 327, no. 5969, pp. 1135–1139, 2010. View at: Publisher Site | Google Scholar
  74. I. E. Wartz, K. M. O'Rourke, H. Zhou et al., “De-ubiquitination and ubiquitin ligase domains of A20 downregulate NF-κB signalling,” Nature, vol. 430, no. 7000, pp. 694–699, 2004. View at: Publisher Site | Google Scholar
  75. Y. Yoshida, A. Murakami, and K. Tanaka, “Skp1 stabilizes the conformation of F-box proteins,” Biochemical and Biophysical Research Communications, vol. 410, no. 1, pp. 24–28, 2011. View at: Publisher Site | Google Scholar
  76. M. Shirane, S. Hatakeyama, K. Hattori, K. Nakayama, and K. Nakayama, “Common pathway for the ubiquitination of IkappaBalpha, IkappaBbeta, and IkappaBepsilon mediated by the F-box protein FWD1,” The Journal of Biological Chemistry, vol. 274, no. 40, pp. 28169–28174, 1999. View at: Google Scholar
  77. M.-P. Bousquet-Dubouch, E. Baudelet, F. Guérin et al., “Affinity purification strategy to capture human endogenous proteasome complexes diversity and to identify proteasome-interacting proteins,” Molecular & Cellular Proteomics, vol. 8, no. 5, pp. 1150–1164, 2009. View at: Publisher Site | Google Scholar

Copyright © 2015 Arielis Rodríguez-Ulloa et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


More related articles

997 Views | 386 Downloads | 3 Citations
 PDF  Download Citation  Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

We are committed to sharing findings related to COVID-19 as quickly and safely as possible. Any author submitting a COVID-19 paper should notify us at help@hindawi.com to ensure their research is fast-tracked and made available on a preprint server as soon as possible. We will be providing unlimited waivers of publication charges for accepted articles related to COVID-19. Sign up here as a reviewer to help fast-track new submissions.