Multi-ACPNet: A multi-scale sequence-structure feature fusion framework for anticancer peptide identification and functional prediction

TL;DR

Multi-ACPNet is a new AI model that identifies anticancer peptides and predicts their effectiveness against different cancer types by combining sequence and structural data.

Contribution

The novel integration of sequence and structural features through a multi-scale fusion framework enables dual-task prediction of anticancer peptides and their functional activities.

Findings

Multi-ACPNet achieves high accuracy (0.8140-0.9536) in identifying anticancer peptides across benchmark datasets.

The model outperforms existing methods in functional activity prediction with an AUC of 0.9033 and F1-score of 0.8472.

It enables multi-label prediction of ACP activity across seven cancer cell types and generalizes well to toxicity prediction tasks.

Abstract

Anticancer peptides (ACPs) have emerged as promising therapeutic candidates for cancer treatment due to their high efficacy and low propensity for inducing drug resistance. However, existing ACP identification methods primarily rely on peptide sequence features while neglecting spatial structural characteristics. Moreover, few approaches can simultaneously predict the functional activity of ACPs. To address these limitations, this study proposes Multi-ACPNet, a novel dual-function predictor capable of both ACP identification and activity type classification. This model innovatively integrates sequence and structural features through a multi-stage framework. It employs a hybrid Bidirectional Long Short-Term Memory (BiLSTM) and causal convolutional network to capture both long-range dependencies and local sequence patterns, followed by a multi-scale Graph Convolutional Network (GCN) that dynamically fuses local and long-range structural dependencies using residual connections and adaptive weighting. Experimental results demonstrate that Multi-ACPNet achieves outstanding performance, with Accuracy of 0.8140, 0.9536, and 0.8770 on three benchmark datasets for ACP identification. For functional prediction, it attains an AUC of 0.9033, F1-score of 0.8472, and Hamming loss of 0.1303, significantly outperforming state-of-the-art predictors.

Anticancer peptides (ACPs) have emerged as highly promising therapeutic candidates in cancer treatment due to their high efficacy and low propensity for drug resistance. However, existing prediction methods suffer from two major limitations: first, they predominantly rely on sequence information while neglecting three-dimensional structural features, leading to the loss of crucial spatial interaction information; second, their predictive capability is insufficient—most current models are limited to binary classification tasks for ACP identification and cannot predict their targeting activities against specific cancer cell lines. To address these challenges, we developed Multi-ACPNet, an end-to-end framework that integrates both sequence and structural information for dual-task collaborative prediction. The model innovatively incorporates a multi-scale information fusion mechanism: at the sequence level, it combines Bidirectional Long Short-Term Memory (BiLSTM) with causal convolution to simultaneously capture long-range dependencies and local functional motifs; at the structural level, a multi-hop Graph Convolutional Network is constructed to dynamically learn multi-level spatial interactions within peptide molecules. This architecture not only significantly improves the accuracy of ACP identification but also enables efficient multi-label prediction of functional activities across seven cancer cell types. Furthermore, the model demonstrates excellent generalization capability, achieving competitive performance in extended tasks such as peptide toxicity prediction.

Figures

Overview of the proposed workflow.(A) Multi-ACPNet architecture overview. First, ESM C features are extracted and fed into the (B) Sequence Multi-Scale Network to capture multi-scale sequence features. Then, graphs are constructed using trRosetta, with domain knowledge features and sequence-based learned features as node features. (C) The Graph Multi-Scale Network performs multi-scale structural feature learning on both the full graph and key subgraphs. The Classifier outputs the final prediction results.

ROC curves of the optimal model through 5-fold cross-validation on AntiCP 2.0 training datasets.(A) Cross-validation result on the AntiCP 2.0_Main training dataset. (B) Cross-validation result on the AntiCP 2.0_Alternate training dataset.

Performance comparison between the proposed method and existing methods across seven cancer types.

Performance comparison of sequence-only and structure-only models on the AntiCP 2.0_Main test set.

t-SNE visualization of feature representations from the penultimate layer on AntiCP 2.0_Main test set.(A) ACPNet_no_Graph feature distribution. (B) ACPNet_no_Seq feature distribution. (C) Multi-ACPNet feature distribution.

Agglomerative clustering of feature vectors from the penultimate network layer on ACP-Mixed-80 test set.

Comparative performance on the ACP toxicity test set ((s2).(A) Confusion matrix for ToxGIN. (B) Confusion matrix for Multi-ACPNet.

Equations

Full text

1. Introduction

Cancer remains one of the most severe threats to global human health. Current clinical interventions primarily include surgery, radiotherapy, chemotherapy, and targeted therapy. However, these conventional treatments are significantly limited by tumor heterogeneity, drug resistance, and severe side effects. In recent years, anticancer peptides (ACPs) have emerged as novel therapeutic agents, demonstrating superior efficacy, high selectivity, strong specificity, and excellent tumor penetration capabilities. These advantages allow ACPs to overcome many limitations of traditional therapies, positioning them as promising candidates for cancer treatment [1,2]. Nevertheless, only a limited number of ACPs have been identified, with even fewer approved for clinical use. To date, only 28 ACP-based drugs have been approved by the U.S. Food and Drug Administration (FDA) and the European Medicines Agency (EMA) [3]. Consequently, the development of novel and potent ACPs has become a critical research direction in cancer therapeutics.

ACPs are short peptide molecules that exhibit significant antitumor activity. Traditional ACP discovery relies on experimental screening, which is time-consuming, costly, and inefficient. With the rapid advancement of bioinformatics and artificial intelligence (AI), computational prediction methods have emerged as a powerful tool for accelerating ACP prescreening. These approaches streamline the labor-intensive discovery process [4,5] and significantly improve the efficiency of peptide drug candidate identification.

Current AI-based methods for ACP identification mainly consist of two key steps: feature encoding based on sequence information, followed by classifier training and selection. Traditional machine learning approaches still dominate this field [6]. Tyagi et al. [7] developed the first ACP predictor using Support Vector Machine (SVM) as the classifier. Agrawal et al. [8] extracted multiple input features and implemented six machine learning classifiers on two established datasets. Karim et al. [9] proposed ANNprob-ACPs, where they evaluated 33 machine learning models, selected the six best-performing ones, and integrated their probability scores as input to an Artificial Neural Network (ANN) meta model, achieving significantly improved classification accuracy.

ACP-DL [10] pioneered the first deep learning framework for ACP identification by integrating binary profile features (BPF) and reduced-alphabet k-mer sparse matrices of peptide sequences, utilizing a Long Short-Term Memory (LSTM) network to automatically discriminate between ACPs and non-ACPs, achieving superior performance on benchmark datasets. Liu et al. [11] developed a parallel architecture combining Convolutional Neural Network (CNN) and Recurrent Neural Network (RNN) for learning sequence embeddings through ensemble learning. Zhang et al. [12] proposed an end-to-end ACP identification model utilizing both attention-augmented convolutional neural networks (AAConv) and standard CNNs for sequence feature extraction, with multi-head attention mechanisms further refining high-level features, achieving competitive classification performance. Notable contributions also include iACP-DRLF [13], ACP-OPE [14], ACP-CapsPred [15], CAPTURE [6], GRDF [16], AI4ACP [17], ACP_MS [18], ACP-ML [19], ACPred-BMF [20], ACPred-LAF [21], ME-ACP [22], ACP-CLB [23], and ACP-LSE [24], collectively advancing the field through diverse computational approaches.

Despite significant progress in ACP predictor development, several critical limitations persist. First, current models inadequately capture multi-scale feature interactions, potentially overlooking hierarchical relationships between local and global peptide characteristics. Second, although ACP tertiary structures are relatively simple, the specific residue arrangement may directly influence peptide function. Current predictors predominantly process sequence information alone, largely disregarding these essential structural features. Most notably, most methods are limited to binary classification (ACP or non-ACP), with few capable of predicting the functional types of ACPs. However, ACP research involves not only determining whether peptides have anticancer properties, but also predicting which cancer cell lines they can target, which is an equally important aspect.

To address these limitations, we propose Multi-ACPNet, a graph convolutional-based multidimensional collaborative framework that integrates sequence information, structural features, and multi-scale characteristics for comprehensive ACP prediction. Our primary contributions can be summarized as follows:

To perform sequence analysis, we integrate Bidirectional Long Short-Term Memory (BiLSTM) and causal convolution to comprehensively capture and model both global contextual information and local functional motifs in peptide sequences.We further design a multi-scale residual dynamic Graph Convolutional Network (GCN) that computes adjacency matrices with varying hop distances to adaptively capture critical structural features from both the full graph and key subgraph at varying spatial scales.Our network is capable of dual-functional prediction. Using high-quality benchmark datasets, we first constructed and optimized a binary classifier, and further developed a multilabel classifier to accurately predict peptide sequences’ specific bioactivities against different cancer cell lines, providing critical decision support for personalized therapy.Beyond core anticancer function prediction, we have also successfully validated our model’s generalization capability on peptide toxicity prediction tasks, where it demonstrates competitive and robust performance, thereby expanding its utility for comprehensive safety assessment in therapeutic peptide discovery.

2. Methods

2.1. Datasets

For the first task (predicting whether a peptide is an ACP), we employ three datasets: the ACP-Mixed-80 [21], AntiCP 2.0_Main and AntiCP 2.0_Alternate [8]. The ACP-Mixed-80 dataset was constructed by integrating multiple benchmark datasets [8,25–28], followed by rigorous preprocessing including label verification, label correction or removal, duplicate elimination, and separation of positive and negative samples. Sequences with > 80% similarity were further removed using CD-HIT [29]. The negative samples in the last two datasets consisted of antimicrobial peptides (AMPs) lacking anticancer properties and random peptides generated using SwissProt proteins, respectively. All three datasets were pre-partitioned into training and independent test sets (see Table A in S1 Text for detailed statistics).

For the second task involving multi-label functional prediction, we use the dataset from ACP-MLC [30], which was collected from CancerPPD [31]. The dataset was rigorously processed by retaining only functional categories with over 40 entries, removing duplicate sequences, non-linear peptides, sequences containing non-standard amino acids, and peptides > 100 amino acids in length. Additionally, sequences with over 90% similarity were filtered using CD-HIT. The final dataset comprises seven tissue-specific cancer types: Colon, Breast, Cervix, Skin, Lung, Prostate, and Blood. (detailed in Table B in S1 Text).

2.2. Overview of the model

Our Multi-ACPNet (Fig 1A) integrates peptide sequence and structural information through three key steps: (1) Feature encoding, where we comprehensively encode peptide sequences using ESM C [32] embeddings, BPF, positional encoding, and AAindex [33] features; (2) Sequence-based representation learning, implemented via a Sequence Multi-Scale Network (detailed in Fig 1B) that combines a BiLSTM and parallel CNN module to extract sequence patterns from ESM C embeddings; and (3) Structure-based geometric learning, which represents peptides as graph structures with fused features as node attributes, followed by structural analysis through a Graph Multi-Scale Network (as shown in Fig 1C). The processed features are then aggregated via global pooling and fed into fully connected (FC) layers for prediction.

Overview of the proposed workflow.(A) Multi-ACPNet architecture overview. First, ESM C features are extracted and fed into the (B) Sequence Multi-Scale Network to capture multi-scale sequence features. Then, graphs are constructed using trRosetta, with domain knowledge features and sequence-based learned features as node features. (C) The Graph Multi-Scale Network performs multi-scale structural feature learning on both the full graph and key subgraphs. The Classifier outputs the final prediction results.

2.3. Feature encoding

To encode peptide sequences as numerical vectors and enhance the model’s learning capability through multi-perspective feature representation, we select four different encoding methods: ESM C pretrained embeddings, binary profile features, positional features, and AAindex features. Among these, ESM C provides global deep semantic information, while BPF and position features preserve local sequence patterns, and AAindex supplements physicochemical properties. This hybrid encoding approach improves the model’s peptide sequence representation capability by integrating learned features with domain knowledge features.

In 2024, EvolutionaryScale released the ESM3 [34] series of models along with a parallel series, ESM C, which is specifically designed to capture the intrinsic biological representations of proteins. Based on the Transformer architecture, ESM C automatically learns general protein representations through unsupervised training on large-scale protein sequence data, significantly surpassing its predecessor ESM2 [35] in both predictive performance and computational efficiency. Compared to some very large models, ESM C achieves faster inference speed and higher computational efficiency with fewer parameters (e.g., the 300M and 600M versions). In this study, we used ESMC_600M as the foundational feature extraction model. The model automatically produces residue-level embeddings with dimensions $[eqn]$ , where $[eqn]$ is the peptide sequence length. These embeddings are subsequently processed by the Sequence Multi-Scale Network.

The BPF encoding represents each amino acid residue as a 20-dimensional one-hot vector, where only the position corresponding to the specific residue type is set to 1 while others to 0. For example: Alanine (A) is encoded as $[eqn]$ and Tyrosine (Y) is encoded as $[eqn]$ . This encoding method transforms a peptide into a binary matrix with dimensions $[eqn]$ .

Position encoding is used to introduce the position information of amino acids, generating position encoding using sine and cosine functions of different frequencies, as shown in Formulas 1 and 2.

$$[eqn]$$ $$[eqn]$$

where, pos represents the position of amino acids in the peptide sequence, and i represents the dimension index (( $[eqn]$ ), where $[eqn]$ in our study corresponds to the embedding dimensionality. Based on these formulas, we obtain the position encoding with a dimension of $[eqn]$ .

AAindex is a database storing the physicochemical and biochemical properties of amino acids, containing over 500 characteristics for each amino acid, such as hydrophobicity, charge, volume, etc. In AAindex version 9.2, a total of 566 physicochemical properties were collected, from which we select 553 properties without missing values. In this study, to mitigate potential issues such as redundancy, multicollinearity, and overfitting associated with high-dimensional feature usage, we employ the mRMR (Max-Relevance and Min-Redundancy) [36] algorithm to refine the original AAindex feature set. This method selects a subset of features that maximizes relevance to the target variable while minimizing inter-feature redundancy. Experimental findings in Section 3.6 indicate that the overall performance is optimal after reducing the feature dimension to 60. Beyond this point, further increasing feature dimensionality yields no substantial gains and may even lead to performance decline. This outcome confirms that mRMR-based feature selection effectively reduces noise and redundancy while retaining essential physicochemical information, thereby strengthening the model’s robustness and generalization capacity. Consequently, the final AAindex feature representation has a dimensionality of $[eqn]$ .

2.4. Sequence-based representation learning

ESM C is trained on massive, diverse protein sequences, acquiring universal, protein-level contextual representations. However, these representations remain relatively static and generic. Since our dual tasks, including both ACP identification and functional prediction, require learning specific, subtle functional patterns, we design a multi-scale sequence feature processing network (Sequence Multi-Scale Network; Fig 1B) to comprehensively capture multi-scale features of ACP sequences. This network dynamically transforms the general ESM C embeddings into task-specific, sequence-level representations optimized for these ACP-related tasks.

Peptide sequences and natural language both exhibit context-sensitive characteristics. BiLSTM has successfully modeled long-range dependencies between words in Natural Language Processing (NLP), and similarly can be applied to capture functional relationships between residues. In the Sequence Multi-Scale Network, the pretrained ESM C features are first fed into a BiLSTM module. In peptide sequences, residues at both the N-terminus and C-terminus may jointly determine functionality. Traditional unidirectional LSTM would lose reverse dependency information, whereas BiLSTM employs two opposing LSTMs to process the sequence. This architecture effectively models bidirectional long-range interactions between amino acid residues in the sequence. LSTM selectively remembers or forgets information in a sequence by means of unique gating mechanisms (forget, input, and output gates), thus capturing key residues in the sequence that are far apart but functionally relevant.

The functions of many ACPs rely on short, conserved local motifs, and causal convolution serves as an ideal tool for capturing such features through its local perception and historical dependency. The output of the BiLSTM is then fed into a local information processing module that adopts a dual-branch parallel architecture. The first branch uses causal convolution, LeakyReLU, and Dropout. Causal convolution maintains consistent input and output sequence lengths by zero-padding on the left side of the sequence, while ensuring that the output at each time step relies only on sequence information from the current moment and before. We use a causal convolution with a kernel size $[eqn]$ of 3 and a zero padding of $[eqn]$ , guaranteeing that each output step relies only on the current and two preceding input elements in the sequence. This preserves the causality of local functions while capturing amino acid local patterns. The second branch utilizes standard 1D convolution ( $[eqn]$ ) and LeakyReLU. The two branches are combined through summation, producing an output feature of dimensions $[eqn]$ , which preserves local details while integrating global context.

2.5. Structure-based geometric learning

The function of ACPs may be influenced by structural features at different scales. Therefore, we innovatively propose a framework for learning the structure of ACPs based on Graph Neural Network (GNN), which achieves multi-granularity modeling of spatial interactions in ACPs and improves prediction accuracy by integrating multi-scale neighborhoods with adaptive fusion mechanism.

2.5.1. Graph construction.

We employ trRosetta [37] to predict the distance map of peptides. trRosetta predicts the probability distribution of $[eqn]$ distances between residues directly from multiple sequence alignment (MSA) by end-to-end training. By comparing the experimental performance, we set a distance threshold of $[eqn]$ 10Å to construct the initial adjacency matrix $[eqn]$ :

$$[eqn]$$

where $[eqn]$ denotes the Euclidean distance between the $[eqn]$ atoms of residue $[eqn]$ and residue $[eqn]$ .

The node features of the graph are constructed by concatenating the outputs from the Sequence Multi-Scale Network, BPF encoding features, positional encoding features, and AAindex encoding features, generating an $[eqn]$ dimensional node feature matrix $[eqn]$ .

To construct the multi-scale graph, multi-hop adjacency matrices are generated via matrix power operations:

$$[eqn]$$

where, $[eqn]$ denotes the sign function for matrix binarization. In the k-hop adjacency matrix, a value of 1 indicates the existence of a path of length k between the two nodes.

2.5.2. Graph network architecture.

We develop the Graph Multi-Scale Network (Fig 1C) to perform hierarchical graph convolutions. The network first captures multi-scale structural features across the global peptide graphs. Since the full graph may contain non-functional residues that could interfere with key feature learning, we further retain critical residue nodes to extract multi-scale features from important subgraphs. The network primarily consists of three components: the Multi-GCN Module, adaptive fusion and top-k pooling.

The Multi-GCN Module performs parallel graph convolution operations. After constructing the graphs at three scales (1-hop, 2-hop, 3-hop), graph convolutions are performed in parallel according to Formula 5 to independently process neighborhood information at each scale.

$$[eqn]$$

where $[eqn]$ and $[eqn]$ denote the output feature matrices of the l-th and (l-1)-th graph convolution layers respectively, $[eqn]$ is the initial input, $[eqn]$ is the LeakyReLU activation function, $[eqn]$ is the adjacency matrix with self-loops added, $[eqn]$ is the degree matrix, $[eqn]$ is a trainable weight matrix, and $[eqn]$ represents the graph convolution at 3 scales.

Through a single parallel graph convolution operation, we obtain three feature matrices: $[eqn]$ . To achieve adaptive fusion, we propose a dynamic weighting mechanism, calculated as follows:

$$[eqn]$$

where $[eqn]$ is the learned attention weights, and $[eqn]$ is a row vector with all elements equal to 1. The symbol $[eqn]$ represents outer product. The weights are expanded to an $[eqn]$ -dimensional matrix via the outer product, where each column of the expanded matrix is a copy of $[eqn]$ . Features are scaled by the attention weights through Hadamard product $[eqn]$ . The fused feature $[eqn]$ is obtained by summing all features after attention-weighted scaling.

To obtain the attention weights $[eqn]$ , $[eqn]$ , $[eqn]$ , we first concatenate the features $[eqn]$ and apply a linear transformation:

$$[eqn]$$

where $[eqn]$ represents concatenation, $[eqn]$ and $[eqn]$ are learnable weight and bias, $[eqn]$ is projected into a 3-D space corresponding to the three scales. Then, normalize using softmax to get $[eqn]$ :

$$[eqn]$$

To ensure that critical local features are not completely covered by subsequent processing, we introduce residual connectivity:

$$[eqn]$$

The above describes entire process of the Multi-GCN Module and adaptive fusion over the complete graph. To focus on critical peptide structural regions while suppressing noisy nodes, we subsequently employ top-k pooling, which scores nodes via a learnable projection vector and retains only the top-n nodes. Specifically, we preserve the top 50% highest-scoring nodes to perform multi-scale graph convolution and adaptive fusion over this subgraph.

2.5.3. Output layers.

To enable downstream classification, we employ both global average pooling and global maximum pooling to form the final graph representation vector. The classifier comprises two fully connected layers for nonlinear mapping.

For the ACP recognition task, we use a single output neuron with sigmoid activation to distinguish ACPs from non-ACPs. For functional prediction, samples can simultaneously belong to multiple categories, thus, the output layer contains seven neurons with sigmoid activations, each independently computing the probability for one category.

3. Results and discussion

3.1. Evaluation metrics

In the ACP recognition task, we adopt widely used evaluation criteria including Sensitivity (SE), Specificity (SP), Accuracy (ACC), Precision, F1-score and the Matthews correlation coefficient (MCC), Area Under the ROC Curve (AUC). These metrics are defined as follows:

$$[eqn]$$

where TP (true positives) denotes the number of correctly predicted positive samples, and TN (true negatives) represents the number of correctly predicted negative samples. FP (false positives) and FN (false negatives) indicate misclassified negative and positive samples, respectively.

For the multi-label functional prediction task, we employ label-based multilabel evaluation with two averaging modes: Macro-averaging and Micro-averaging.

Macro-averaging computes metrics separately for each class and then takes the arithmetic mean:

$$[eqn]$$

where L denotes the total number of label categories. Micro-averaging computes global metrics by aggregating TP, FP, TN, and FN for all categories. The calculation process is as follows:

$$[eqn]$$

Additionally, we employ Hamming loss, a key multilabel classification metric that quantifies the discrepancy between predicted and ground-truth labels. It is defined as:

$$[eqn]$$

where $[eqn]$ denotes the number of samples, $[eqn]$ indicates the ground-truth label of sample $[eqn]$ for category $[eqn]$ (1 indicates presence of the function, 0 indicates absence). $[eqn]$ represents the predicted label, and $[eqn]$ is the indicator function that equals 1 for misclassified labels and 0 otherwise.

3.2. Experimental setup

We implement the Multi-ACPNet model based on the PyTorch framework, using Adam as the optimizer with an initial learning rate of 0.0001. The model computations are GPU accelerated on an NVIDIA GeForce RTX 4090, running on a Windows system with an Intel Core i9-14900K processor and 64GB RAM. The hyperparameters are optimized on the AntiCP 2.0 training set through grid search with 5-fold cross-validation (see Table C in S1 Text for details). Notably, we maintain identical hyperparameters configurations for both binary classification and multi-label prediction tasks, achieving task switching solely by modifying the output layer dimensionality. This design significantly enhances the model’s extensibility while preserving its generalizability across different prediction scenarios. To prevent potential data leakage arising from sequence similarity across different datasets, we perform rigorous sequence-level redundancy removal prior to hyperparameter optimization. Specifically, to ensure that no highly similar sequences exist between the training data used for hyperparameter optimization and the independent evaluation sets, we systematically remove sequence redundancy using the CD-HIT tool. Operationally, we first merge the AntiCP 2.0 training set with the AntiCP 2.0 test set, the ACP-Mixed-80 test set, and the multi-label functional prediction dataset. CD-HIT is then executed to perform full-sequence alignment and clustering. During this process, only entries from the training set partition that exhibit >90% similarity with sequences in any evaluation set are removed, while the evaluation sets themselves remain unaltered to preserve their completeness for subsequent fair benchmarking comparisons with other methods. This procedure ensures strict sequence-level independence between the training set and all evaluation sets. Table D in S1 Text presents the sequence similarity analysis between the AntiCP 2.0 training set and each evaluation dataset. It is important to clarify that for the multi-label functional dataset, all sequences are included in the similarity analysis because the entire dataset undergoes 10-fold cross-validation (as described in Section 3.4), meaning every sequence serves as a test sample during the evaluation process. Fig 2 displays the ROC curves of the optimal model evaluated through 5-fold cross-validation on both the AntiCP 2.0_Main and AntiCP 2.0_Alternate training datasets. The corresponding confusion matrices and key evaluation metrics are presented in Fig A in S1 Text.

ROC curves of the optimal model through 5-fold cross-validation on AntiCP 2.0 training datasets.(A) Cross-validation result on the AntiCP 2.0_Main training dataset. (B) Cross-validation result on the AntiCP 2.0_Alternate training dataset.

To conduct the evaluation of each dataset, we train a separate model. We first load the model hyperparameters. Then, we split the training set into training and validation subsets at a 5:1 ratio. Model selection is performed based on validation performance, and the model achieving the highest Accuracy on the validation set is subsequently used for final evaluation on the independent test set. To ensure a fair comparison with existing benchmark methods (which are trained and evaluated on complete public datasets), we likewise utilize the complete original datasets. Concurrently, we systematically analyze the internal sequence similarity within each dataset: the sequence similarity between the training and test sets of ACP-Mixed-80 is below 80%, and the internal similarity of the multi-label functional prediction dataset is below 90%. These thresholds align with the dataset preprocessing standards described in Section 2.1. For AntiCP 2.0_Main and AntiCP 2.0_Alternate, a certain degree of sequence similarity exists between their training and independent test sets (details in Table D in S1 Text), but no identical sequences are present. This analysis ensures effective separation between training and test sets at the sequence-independence level, thereby maintaining evaluation fairness while effectively preventing data leakage.

3.3. Comparison of the proposed method with existing methods for classifying ACPs and non-ACPs

In the task of ACP identification, we compare Multi-ACPNet with existing ACP prediction methods, and the results of all the compared methods are obtained from their original publications. Table 1 presents the performance comparison of the proposed Multi-ACPNet with the existing advanced methods on the ACP-Mixed-80 dataset. Our method achieves superior performance across multiple key metrics: it achieves the highest Accuracy, MCC, Precision and AUC scores. Notably, while maintaining high Sensitivity (0.9016), our method also attains competitive Specificity (0.8525), achieving an excellent balance between these metrics. In contrast, existing methods commonly show performance imbalance issues.

Table 1: Performance comparison between the proposed method and existing methods on the ACP-Mixed-80 dataseta.

For the AntiCP 2.0_Main and AntiCP 2.0_Alternate datasets, we first train models on the complete training sets and evaluate them on the independent test sets, consistent with other comparative methods. As shown in Table 2, our Multi-ACPNet achieves the highest ACC (0.8140 and 0.9536) and MCC (0.6283 and 0.9103) values on both datasets. And, on the AntiCP 2.0_Alternate dataset, our predictor attains a remarkable SP of 0.9948, outperforming all existing methods. To preclude potential performance overestimation caused by redundant sequences, we remove training data with >90% sequence similarity to the test set and evaluate model performance under this condition (this version is denoted as Multi-ACPNet-CD). Experimental results demonstrate that even after removing redundant data, our model still exhibits superior performance compared to existing methods that do not perform redundancy elimination. On the AntiCP 2.0_Main dataset, most existing methods exhibit either high SE with low SP, or conversely low SE with high SP. For instance, AntiCP achieves perfect SE (1) but extremely low SP (0.0116), reflecting complete positive sample detection at the cost of an excessively high false positive rate. Conversely, ACPred-LAF shows the highest SP (0.8895) but relatively low SE (0.6337), suggesting accurate negative sample exclusion but compromised detection rate for true positives. However, our approach maintains a good balance (SE = 0.8214, SP = 0.8021) while achieving competitive rankings, substantially reducing misclassification risks in practical applications.

Table 2: Performance comparison between the proposed method and existing methods on the AntiCP 2.0 independent test dataseta.

Multi-ACPNet achieves a false positive rate (FPR) of 19.79% on AntiCP 2.0_Main (where negative samples are AMPs) and a Specificity of 0.9948 (FPR as low as 0.52%) on AntiCP 2.0_Alternate (where negative samples are random peptides). This indicates that compared to structurally random peptides, the model finds it more difficult to distinguish AMPs from ACPs. However, the two models are independently trained on different data distributions, which could itself contribute to performance differences. Therefore, we conduct an additional experiment: we merge the training sets of both datasets to train a unified model and evaluate it separately on their respective independent test sets. This approach ensures that the model is exposed to both negative sample types during training, providing a more controlled test of the hypothesis that distinguishing ACPs from AMPs is inherently more difficult. The results (Fig B in S1 Text) show that the model still attains a high FPR of 19.77% on AntiCP 2.0_Main, while achieving only 1.55% on AntiCP 2.0_Alternate. These experiments suggest that ACPs and AMPs may possess similar structural characteristics. Consequently, employing AMPs as negative samples substantially elevates the classification challenge, which likely accounts for the notably inferior performance of Multi-ACPNet on AntiCP 2.0_Main compared to AntiCP 2.0_Alternate.

A common issue in the ACP prediction literature is the frequent neglect of model scale signatures, such as the number of parameters, computational cost, and inference time. However, this evaluation is standard practice in the broader machine learning literature to ensure fair and insightful comparisons [48]. To comprehensively evaluate model efficiency, we compare the model scale and inference performance of all deep learning methods listed in Tables 1 and 2. Under unified hardware conditions, we systematically measure the number of parameters (Params), floating-point operations per sample (FLOPs), and inference time per sample for each model. As shown in Table E in S1 Text, although Multi-ACPNet exhibits a higher number of parameters and greater computational cost per sample than all compared deep learning models, its inference time remains competitive without significant delay. In contrast, the inference time for a single sample in models such as AI4ACP, GRCI-Net, and ACPred-BMF exceeds one second, indicating substantially lower efficiency. Furthermore, we evaluate the batch inference efficiency of the model by processing the entire test set as a single batch. The total processing time for the ACP-Mixed-80 dataset (122 samples) is only 35.68 ms, while the AntiCP 2.0_Main (344 samples) and AntiCP 2.0_Alternate (388 samples) datasets require only 110.50 ms and 168.80 ms, respectively. These results demonstrate that Multi-ACPNet possesses strong real-time processing capability. To rule out the possibility that performance improvements stem merely from increased model scale, we conduct in-depth ablation experiments in Section 3.5. The results indicate that the performance advantage primarily arises from the deep integration of sequence and structural information, rather than from simply increasing the number of parameters.

The comparative results demonstrate that our method maintains reasonable inference overhead and achieves an improved balance between Sensitivity and Specificity across all three datasets while attaining optimal performance in other metrics. This is primarily attributed to the introduction of graph convolution in the proposed model, which captures peptide structural features, overcoming the limitation of previous methods that relied solely on sequence information. Our network employs BiLSTM and causal convolution to extract multi-scale sequential patterns (preventing missed detections) while utilizing distance maps to identify spatial interactions (reducing false positives). The synergy between sequence and structural features significantly enhances the recognition capability for complex peptides.

3.4. Performance of functional activity prediction of ACP against various cancers

In the functional activity prediction task, we compare our model with two benchmark methods: the current state-of-the-art DUO-ACP [49] and ACP-MLC. For a fair comparison, we use the same strategy of 10-fold cross-validation as employed in the papers.

Table 3 presents the comparison of three models on seven evaluation metrics, including Macro-averaging scores and Hamming loss. The results reveal that our Multi-ACPNet achieves statistically superior performance on all metrics. Specifically, it achieves an ACC of 0.8697 and MCC of 0.7309, significantly outperforming the other two models. Notably, it attains a Hamming loss of only 0.1303, indicating minimal discrepancy between predicted and ground truth labels. This marked improvement in predictive performance is achieved with only a marginal increase in inference time, underscoring the efficiency of our architecture. Furthermore, in order to deeply analyze the model’s specific performance for each cancer category, we comprehensively evaluate the predictive performance across seven tumor cell categories (Table F in S1 Text). Multi-ACPNet achieves excellent prediction ability on all cancer categories. However, it shows relatively lower Sensitivity in the Blood category, primarily due to the limited samples that hindered effective feature learning. To demonstrate our model’s superior performance, we similarly evaluate all three competing models across cancer types. It should be noted that only four metrics (ACC, AUC, F1-score, and MCC) were reported for DUO-ACP in the original publication. Here we also show the performance of the DUO-ACP model on only four metrics. As illustrated in Fig 3, ACP-MLC shows competitive performance in Colon cancer but shows limited performance in other tissue types. DUO-ACP generally achieves higher scores across its four reported metrics, yet its overall performance remains inferior to our Multi-ACPNet, which attains optimal results for most evaluation metrics across all cancer types.

Table 3: Performance comparison between the proposed method and existing methods for ACP functional activity predictiona.

Performance comparison between the proposed method and existing methods across seven cancer types.

For Micro-averaging, Multi-ACPNet achieves an ACC of 0.8697, Precision of 0.8763, F1-score of 0.8472, and SE of 0.8371, outperforming DUO-ACP by 1.52% in F1-score and 1.31% in SE. These results demonstrate that our proposed model can more reliably handle the complexity of multi-label tasks and exhibits superior performance in ACP activity prediction. This advantage is mainly due to our innovative design of a multi-scale sequence-structure cooperative learning framework, which simultaneously captures short-chain functional motifs and long-range sequence patterns at the sequence level, while resolving multi-scale spatial structures of peptide chains at the structural level. This integrated approach enables the model to more effectively capture discriminative features across different cancer types and specifically recognize the anticancer activity of peptides against various cancer cell lines, thereby maintaining optimal performance under a comprehensive evaluation system.

3.5. Effectiveness analysis of sequence and structure learning

To analyze the synergistic mechanism between sequence and structural information, validate the necessity of sequence learning, and demonstrate the effectiveness of structural information, we design two ablated models while maintaining identical parameter counts: (1) ACPNet_no_Graph, which, following a similar approach to Ortega et al. [50], retains only sequence learning by setting all elements in the adjacency matrix to 1 and all distance values to 1, thereby disabling graph structural learning, and (2) ACPNet_no_Seq, which disables sequence information by setting all input ESM C features to 1. Comprehensive evaluations on the AntiCP 2.0_Main dataset (Fig 4) reveal that ACPNet_no_Graph outperforms the structure-only baseline (ACPNet_no_Seq) across all evaluation metrics. Specifically, it achieves a significantly higher MCC (0.5988). Ultimately, Multi-ACPNet attains optimal performance by integrating both sequence and structural features.

Performance comparison of sequence-only and structure-only models on the AntiCP 2.0_Main test set.

To better demonstrate the effectiveness of sequence-structure synergy and enhance model interpretability, we employ t-distributed stochastic neighbor embedding (t-SNE) [51] for feature visualization. Specifically, we extract features from the penultimate fully connected layer of models and reduce them into a 2D space. Fig 5 presents the visualization results on the AntiCP 2.0_Main independent test set. ACPNet_no_Seq exhibits limited clustering capability, showing satisfactory aggregation only for partial negative samples. ACPNet_no_Graph produces fewer clustering errors. Multi-ACPNet, through sequence-structure fusion, generates two distinct clusters for ACPs and non-ACPs, confirming its superior capability in learning both common and discriminative features. Notably, there remains a substantial number of challenging samples that are difficult to classify correctly, improving the model’s ability to discriminate these difficult samples will be the focus of our future research. Furthermore, Figs C to F in S1 Text present the visualization results on the AntiCP 2.0_Alternate and ACP‑Mixed‑80 datasets. Consistent with the observations on the AntiCP 2.0_Main dataset, ACPNet_no_Graph achieves superior classification performance compared to ACPNet_no_Seq in the ACP identification task. However, as shown in Fig G in S1 Text, a different pattern emerges on the ACPs functional prediction dataset (using Macro-averaging metrics). For this multi‑label functional classification task, ACPNet_no_Seq outperforms ACPNet_no_Graph on most metrics, suggesting that structural information may play a more dominant role in distinguishing specific cancer‑type activities. Nonetheless, the full Multi‑ACPNet model achieves the most optimal classification results across all test sets, further validating the effectiveness and adaptability of our sequence‑structure fusion approach for both binary ACP identification and multi‑label functional prediction tasks.

t-SNE visualization of feature representations from the penultimate layer on AntiCP 2.0_Main test set.(A) ACPNet_no_Graph feature distribution. (B) ACPNet_no_Seq feature distribution. (C) Multi-ACPNet feature distribution.

The above phenomenon indicates that in our proposed model, the Sequence Multi-Scale Network dominates the identification process of ACPs. On the other hand, the Graph Multi-Scale Network further captures residue interactions across different spatial scales, elucidating critical tertiary structural features of peptides, thereby playing a pivotal role in enhancing model performance. The sequence and structural modalities exhibit complementary advantages through synergistic cooperation, collectively establishing a comprehensive mechanism for ACP recognition. These findings demonstrate the macro-effectiveness of the synergistic mechanism between sequence and structural modalities. To further validate the contributions of key internal modules, we conduct ablation studies targeting critical components within the sequence and graph networks. The specific configurations of the ablated models are as follows:

(1) No‑BiLSTM: The BiLSTM layer is removed from the Sequence Multi‑Scale Network, retaining only the causal convolutional module for sequence feature extraction.(2) No‑CausalCNN: The causal convolutional module is removed, with sequence features extracted solely via BiLSTM.(3) No‑TopKPooling: The top‑k pooling layer and the subsequent subgraph convolution module are removed from the Graph Multi‑Scale Network. Global pooling is applied directly to the full‑graph features instead, thereby eliminating the key‑residue selection and hierarchical subgraph learning mechanism.(4) No‑Seq: The entire sequence branch (including both BiLSTM and causal convolution) is removed, and the model relies solely on the multi-scale GNN for prediction.(5) No‑Graph: The graph structure learning branch is removed. The output of the Sequence Multi‑Scale Network is fed directly into the classifier, relying solely on sequence information.

The experimental results are presented in Tables G to J in S1 Text. The dominant role of sequence information is clearly evidenced by the dramatic performance drop observed when the entire sequence branch is removed (No‑Seq). For instance, on the AntiCP 2.0_Main and ACP functional prediction datasets, the MCC scores of No‑Seq decrease by 0.3024 and 0.2327, respectively, compared to the full model. Within the sequence branch, both long‑range dependencies and local patterns are essential. On the functional prediction dataset, removing BiLSTM (No‑BiLSTM) causes a more severe performance degradation than removing causal convolution (No‑CausalCNN), underscoring the greater importance of global context modeling for this task. Conversely, on the three ACP identification datasets, removing causal convolution results in a more significant performance degradation, highlighting the critical role of local motif extraction in distinguishing ACPs from non‑ACPs. In the graph branch, removing top‑k pooling (No‑TopKPooling) consistently reduces performance across all datasets, validating the importance of key‑residue selection and hierarchical subgraph learning.

Critically, performance is not merely a function of model scale. For the first ablation study, where we create ACPNet_no_Seq and ACPNet_no_Graph by holding the number of parameters constant and manipulating the input features to disable sequence or structural learning, respectively, the model performance changes substantially. This confirms that the observed effectiveness depends on the integration of both information modalities, not on scale alone. In the second ablation study involving the removal of key components, no consistent pattern emerges where larger models yield better performance. Comparisons with external algorithms further support this conclusion. For instance, on the ACP-Mixed-80 dataset, the lightweight No-Seq variant (0.44M parameters) outperforms much larger models such as ACPred-LAF (2.46M). Similarly, the No-BiLSTM model (4.11M) achieves better results than the larger MA-PEP (6.41M) and ACP-PDAFF (7.88M). On the AntiCP 2.0_Alternate dataset, No-Seq surpasses ACPred-LAF by 0.0475 in MCC despite having only 17.76% of its parameters. These comparisons clearly indicate that the superiority of Multi-ACPNet stems not from increased model size, but from the effective integration of multi-scale sequence and structural feature extraction.

3.6. Analysis of feature encoding methods

In this section, we conduct an ablation study to evaluate three feature combination schemes (progressively incorporating BPF, position, and AAindex features) across all datasets. Tables K and L in S1 Text present the experimental results. The results demonstrate that starting from BPF, the progressive integration of position and AAindex features not only introduces precise residue position information but also provides a biochemical foundation for structure-function relationships. This integration results in a significant and stable improvement in the model’s predictive performance across both tasks, validating the efficacy of the multi-feature fusion strategy.

To refine the high-dimensional AAindex representation, we apply the mRMR feature selection algorithm. As shown in Tables M and N in S1 Text which present the performance across datasets at different feature dimensionalities, the model achieves its highest ACC and MCC values on the AntiCP 2.0_Main, ACP-Mixed-80, and ACP function prediction datasets when the AAindex features are reduced to 60 dimensions. Furthermore, the highest SP value is also attained on the ACP-Mixed-80 dataset, while performance on the AntiCP 2.0_Alternate dataset remains sufficiently high. Beyond 60 dimensions, further increase in feature number yields no substantial gains and may even lead to performance degradation. Therefore, by selecting 60-dimensional features, we effectively reduce computational cost while retaining the physicochemical properties most relevant to peptide prediction, thereby avoiding the redundancy and interference associated with high-dimensional representations.

In recent years, a wide variety of pre-trained language models (PLMs) have rapidly emerged. In order to demonstrate the superiority of ESMC_600M, we conduct a systematic comparison with several mainstream PLMs, including ESM2, ESM3, ProtT5 [52], and ProteinBERT [53] (detailed in Tables O and P in S1 Text). The experimental results demonstrate that the ESMC_600M model achieves the best performance across multiple key metrics on three ACP classification datasets and attains leading results on all evaluation metrics for the ACP function prediction dataset. Considering its balanced computational efficiency and comprehensive performance superiority, we ultimately select ESMC_600M as our foundational feature extraction pre-trained model.

3.7. Parameter stability analysis

To evaluate the robustness of the Multi-ACPNet model and its sensitivity to key hyperparameters, we analyze the impact of four core hyperparameters: learning rate, distance threshold ( $[eqn]$ ), dropout rate, and TopKPooling ratio. Detailed performance metrics for each parameter across different datasets are provided in Tables Q to X in S1 Text, while Fig H in S1 Text visually illustrates the trend of AUC as the parameters vary. To further quantify the model’s sensitivity to these parameters, we calculate the Coefficient of Variation (CV) of the AUC within the tested ranges for each parameter. The CV is defined as the ratio of the standard deviation to the mean, where a lower value indicates that model performance is less sensitive to fluctuations in that parameter and thus reflects stronger robustness.

The comprehensive analysis shows that Multi-ACPNet maintains relatively stable performance across a wide range of most hyperparameters, though its sensitivity varies significantly depending on the parameter and the dataset. The AntiCP 2.0_Alternate dataset consistently exhibits the highest AUC and the smallest CV values across all four parameter variations, indicating that the prediction task for this dataset is relatively simpler, and the model’s performance is least sensitive to parameter adjustments. In contrast, the ACP-Mixed-80 dataset shows the largest CV values when $[eqn]$ , dropout rate, and TopKPooling ratio are varied, suggesting that the prediction performance for this dataset is more sensitive to adjustments in structure construction and regularization strategies. We also observe that the ACPs functional prediction dataset exhibits substantial fluctuations across different learning rates, with a CV as high as 5.35%. Notably, when the learning rate is set to 1e-5, performance deteriorates drastically, with the MCC dropping to as low as 0.4181. This highlights the critical importance of selecting an appropriate learning rate for ensuring stable and reliable model performance.

3.8. Interpretability analysis

Beyond the t-SNE-based feature visualization in Section 3.5, we further employ agglomerative clustering [54] on the learned feature representations to investigate the decision-making mechanism of Multi-ACPNet and enhance the model’s interpretability. We extract 256-dimensional feature vectors from the penultimate layer of the model and select 20 positive ACPs and 20 non-ACPs from the ACP-Mixed-80 independent test set, resulting in a total of 40 samples for cluster analysis. As shown in Fig 6, the clustering results of the 40 samples in the 256-dimensional feature space demonstrate a clear separation between the two classes. The ACPs (bottom) and non-ACPs (top) form two distinct clusters with a noticeable boundary between them, indicating that the features learned by the model effectively differentiate the two types of peptides. In particular, the ACPs exhibit a more consistent distribution in the feature space, reflecting shared physicochemical properties that are likely critical determinants of their anticancer function. In contrast, the non-ACPs show greater diversity in their feature representation, which aligns with the inherent heterogeneity of this category—encompassing various functional peptides such as antimicrobial peptides, signaling peptides, and random peptides without anticancer activity. Although the high-dimensional features learned by the model are abstract in nature, the distinct clustering patterns strongly suggest that they capture systematic biological differences. These salient feature dimensions not only enhance the credibility of the classification decisions but also provide important clues for subsequent biological interpretation: they may correspond to functional motifs or structural domains that have not been fully characterized, offering guidance for the rational design and functional discovery of novel ACPs.

Agglomerative clustering of feature vectors from the penultimate network layer on ACP-Mixed-80 test set.

Based on the residue importance scores derived from TopKPooling, we select 50 correctly predicted ACP sequences from the ACP-Mixed-80 test set to visualize positional residue importance. In the corresponding heatmap (Fig IA in S1 Text), a redder color indicates a greater contribution of the residue at that position to the correct prediction of ACPs. The results clearly demonstrate that residues in the $[eqn]$ -terminal region, especially the first amino acid position, contribute most significantly. This finding further confirms the critical influence of the $[eqn]$ -terminal sequence on the anticancer activity of peptides, which is consistent with previous research [15].

Fig IB in S1 Text further presents the sequence logos of the first ten amino acid positions at the $[eqn]$ -terminus for the positive samples in the ACP-Mixed-80 test set. Statistical analysis reveals that Glycine (G) occurs most frequently at the first position, while Phenylalanine (F), Leucine (L), and Lysine (K) are significantly enriched within the first ten positions. The crucial functional roles of these amino acids in ACPs have been well-established in prior studies [55,56]. For instance, cationic residues (e.g., K) can penetrate and disrupt cancer cell membranes, thereby inducing cytotoxicity. Concurrently, hydrophobic residues (e.g., F) enhance anticancer efficacy by destabilizing membrane integrity. Importantly, our model successfully captures and prioritizes these functionally critical residues, demonstrating its biological interpretability and alignment with known mechanisms of action.

3.9. Extended application: toxicity prediction

From a translational perspective, the therapeutic potential of ACPs depends not only on their anticancer efficacy but equally on their safety profiles—peptides cytotoxic to normal cells have limited clinical applicability. To comprehensively evaluate peptide candidates, we apply Multi-ACPNet to toxicity prediction, extending its application scope to address this critical safety dimension in peptide therapeutic development.

In the absence of publicly available datasets specifically for ACP toxicity, we construct two dedicated datasets to comprehensively evaluate model performance: a general peptide toxicity dataset S1 for benchmark validation, and an ACP toxicity test set S2 for targeted evaluation. First, to establish the initial data pool, we retrieve 1,783 reviewed toxic peptide sequences from the UniProt [57] database using the search query “KW-0800 AND (length: [10 TO 50])”, and obtain 10,493 non-toxic peptide sequences using the query “NOT KW-0800 AND NOT KW-0020 AND (length: [10 TO 50])”. To specifically assess the model’s ability to predict the toxicity of ACPs, we further screen sequences that overlap with the above‑mentioned data pool from the four ACP datasets employed in this study: ACP-Mixed‑80, AntiCP 2.0_Main, AntiCP 2.0_Alternate, and the ACP functional activity prediction dataset. From this overlap, we construct an independent test set S2, which contains 41 non-toxic ACPs and, due to limited data availability, only 14 toxic ACPs. After excluding the sequences included in S2 from data pool, we reduce redundancy in the remaining peptides using CD‑HIT with a 90% threshold, producing 1,202 unique toxic peptide sequences and 5,811 non-toxic sequences. To ensure class balance, we randomly select an equal number of non‑toxic peptides. We then randomly select approximately 15% of peptides from both the toxic and non-toxic subsets to form the test set, yielding the general peptide toxicity dataset S1 with 2,050 training and 354 test samples.

We first evaluate the model’s performance on the general peptide toxicity dataset S1 and compare it with the recently proposed ToxGIN [58] algorithm, a Graph Isomorphism Network (GIN)-based peptide toxicity predictor. We uniformly employ cross-validation on the training set to determine the model weights and perform the final evaluation on the test set. As shown in Table Y in S1 Text, our model outperforms ToxGIN across all five evaluation metrics with significantly fewer parameters, achieving a Specificity of 0.9096 and a MCC of 0.7748.

We then perform independent external validation of the two trained general models on S2, specifically to examine their generalization capability and practical potential in real-world ACPs toxicity prediction tasks. As shown in Fig 7, ToxGIN demonstrates a slightly higher recognition rate for non-toxic ACPs. However, it exhibits a substantially higher error rate on toxic samples, misclassifying six out of the 14 toxic ACPs as non-toxic. In contrast, our model shows a more balanced and robust performance: it correctly identifies the majority of toxic peptides, while maintaining a comparable recognition rate for non-toxic samples. These results demonstrate that, although our model is not specifically designed for toxicity prediction, it still yields competitive performance in predicting the toxicity of ACPs, highlighting its potential utility in screening for non-toxic ACP candidates.

Comparative performance on the ACP toxicity test set ((s2).(A) Confusion matrix for ToxGIN. (B) Confusion matrix for Multi-ACPNet.

However, the ACP toxicity dataset S2 is limited in size, which to some extent affects the statistical power of model performance evaluation on this subset. Future research will focus on collecting more experimentally validated toxic and non-toxic ACP data to build a more representative benchmark dataset.

4. Limitations

The proposed method possesses clear practical value and significance. From an application perspective, our proposed framework demonstrates robust performance across three critical tasks in peptide drug discovery: anticancer peptide (ACP) identification, functional activity prediction, and toxicity assessment. This capability allows it to serve as an efficient computational screening tool that can accelerate the discovery of novel anticancer peptides by evaluating their efficacy, mechanism, and safety. Methodologically, this work systematically validates the complementarity and synergistic effects between sequence and structural information in peptide prediction, providing a transferable framework for similar tasks. A key future direction is to develop an integrated multi-task model capable of simultaneously identifying ACPs, assessing their toxicity, and predicting their functional activity. Leveraging the powerful multi-scale feature extraction capabilities of Multi-ACPNet and by incorporating a multi-task learning architecture, it is expected to optimize three closely related prediction objectives in parallel based on shared sequence-structure representations. This approach holds promise for providing a powerful computational tool to enable efficient and rational design as well as virtual screening of peptide-based drugs.

Despite these promising aspects and future directions, it is important to acknowledge several current limitations of the study. First, the model’s performance is constrained by the limited availability of high-quality tertiary structural features of peptides. Although we employ advanced tools such as trRosetta for structure prediction, discrepancies inevitably exist between predicted and real experimental structures, which may introduce potential biases during model training and generalization. Future work could explore more robust structure representation learning methods. Furthermore, the current model framework is primarily validated on short peptide sequences. For ultra-long peptides or small proteins, their more complex long-range interactions and folding patterns may require specialized architectural adjustments. Extending the model to handle a broader range of peptide and protein lengths will be a key direction for improving its general applicability. Additionally, the current model architecture is relatively complex. Future work will focus on further optimizing the model architecture to reduce its complexity and parameter size while maintaining predictive performance, thereby shortening inference time and enhancing deployment efficiency. Finally, the imbalanced data distribution across cancer types potentially affects prediction performance for minority classes. As the current model is trained primarily on public datasets and only predicts seven anticancer activity types, future work will focus on enhancing the model’s generalizability and expanding its predictive coverage through integration of more diverse datasets.

5. Conclusion

This study proposes Multi-ACPNet, an innovative dual-function framework capable of simultaneous ACP identification and activity type classification, addressing the limitation of existing methods that rely solely on sequence information. The model achieves multi-level information extraction from sequence features to spatial structural dependencies through deep integration of domain prior knowledge with ESM-2 embeddings, and an original hybrid architecture combining BiLSTM-causal convolution with multi-scale GCN. Experimental results demonstrate that this novel approach significantly improves prediction performance for both tasks, representing a critical advancement for biomedical applications.

Supporting information

S1 TextSupplementary Figures and Tables.(DOCX)