RoBep: a region-oriented deep learning model for B-cell epitope prediction

TL;DR

RoBep is a new deep learning model that improves B-cell epitope prediction by considering spatial clustering of residues, enhancing accuracy and biological relevance.

Contribution

RoBep introduces a region constraint mechanism to model spatial clustering of epitope residues, improving prediction accuracy and biological plausibility.

Findings

RoBep outperforms existing methods with improvements of up to 45% in key performance metrics.

The model ensures predicted epitope residues are spatially compact, enhancing biological plausibility.

RoBep provides both residue-level predictions and antibody–antigen binding regions.

Abstract

Accurate in silico identification of B-cell epitope residues is crucial for antibody design and structure-guided vaccine development. Although recent protein language models and structure-aware methods can capture spatial information of tertiary structure when generating residue embeddings, most existing epitope predictors use these embeddings to perform classification for individual residues one by one, without enforcing spatial continuity for reported epitope residues. Such methods often result in biologically implausible predictions because B-cell epitope residues always cluster together on the antigen surface.

We present RoBep, a region-oriented B-cell epitope predictor that explicitly models the spatial clustering of epitope residues. RoBep introduces a novel region constraint mechanism and combines the advanced protein language model ESM-Cambrian with an equivariant graph neural network. Our method outperforms existing structure-based methods on the benchmark dataset, demonstrating improvements of 26%, 45%, 13%, and 43% in F1, Matthews correlation coefficient, area under the precision–recall curve, and AUROC0.1, respectively. In addition to residue-level predictions, RoBep can also provide antibody–antigen binding regions. Importantly, the predicted epitope residues are ensured to be spatially compact, enhancing biological plausibility and practical relevance for immunotherapeutic design.

A user-friendly website for using RoBep is provided at https://huggingface.co/spaces/NielTT/RoBep. All datasets, source code used in this work, and implementation instructions of the website are publicly available at https://github.com/YitaoXU/RoBep.

Figures

Overview of the RoBep framework for B-cell epitope prediction on a single antigen. The model identifies candidate surface regions using a rolling-sphere strategy, encodes each region as a residue-level graph with both sequence-based and structure-based features, and applies an E(n)-equivariant GNN to extract residue representations. High-confidence regions are scored and selected for downstream residue-level probability prediction, which integrates both local and regional information. The final epitope probability for each residue is obtained by averaging its predictions across all selected regions.

Precision–recall (AUPR) curves on the external test set. RoBep equipped with either ESM-2 or ESM-C, consistently achieves higher precision across all recall levels.

Comparison of cMPD and NoC among predicted epitopes and true epitopes on the independent test dataset. Lower cMPD values indicate tighter clustering. Asterisks denote statistical significance from one-sided Dunnett’s test against RoBep. The horizontal dashed line indicates the mean cMPD of true epitopes across all test samples.

Performance comparison on the test dataset when removing the Region Constraint or other individual components.

Mean AUPR of RoBep across species–length strata in the test set, with counts of non-redundant antigens per stratum.

Predicted antigenicity of antigen and non-structural protein. The red dashed line indicates the threshold of 0.525.

Visual comparison of predicted B-cell epitopes on antigen 5I9Q (chain A) by six methods: (a) RoBep, (b) RoBep without Region Constraint, (c) GraphBepi, (d) SEMA-2.0, (e) DiscoTope-3.0, and (f) CALIBER. The grid sphere for (a) indicates the region with the highest predicted likelihood of being an antibody–antigen interface, as identified by RoBep.

Tables

Table 1. Performance comparison of our model with state-of-the-art methods with variant PLMs on the independent test dataset. Bold indicates the highest value for each metric.

Model	PLM	F1	MCC	AgIoU	AUPR	AUROC0.1	Threshold
Seppa-3.0	–	0.194	0.131	0.107	0.118	0.125	0.15
SEMA-2.0	SaProt	0.284	0.196	0.166	0.236	0.201	0.75
GraphiBepi	ESM-2	0.276	0.189	0.156	0.215	0.172	0.15
DiscoTope-3.0	ESM-IF	0.300	0.234	0.177	0.237	0.231	0.30
CALIBER	ESM-2	0.241	0.164	0.137	0.158	0.156	0.10
RoBep	ESM-2	0.346	0.295	0.209	0.262	0.294	0.15
RoBep	ESM-C	0.379	0.340	0.234	0.268	0.330	0.35

Equations

Full text

1 Introduction

B cells are central to the humoral immune response by producing antibodies that bind to and neutralize potentially harmful antigens (Delves et al. 2017). The specific residues of an antigen recognized by an antibody are termed B-cell epitopes (BCEs). BCEs can be categorized as either linear, contiguous in the primary sequence, or conformational, composed of discontinuous residues brought into proximity by the protein’s 3D structure, with the latter accounting for over 90% of all BCEs (Barlow et al. 1986). Accurate identification of these epitopes is essential for understanding immune mechanisms (Davidson 2000), as well as advancing applications such as mRNA vaccine development (May 2021) and antibody design (Boitreaud et al. 2025). For example, the de novo antibody design model Chai-2 requires antigen target residues as input, underscoring the importance of precise conformational BCE prediction. Therefore, reliable localization of such epitopes is crucial for both basic immunology research and relevant applications.

Current experimental methods, such as X-ray crystallography and nuclear magnetic resonance, offer reliable epitope mapping but are costly and time-consuming, limiting scalability (Gershoni et al. 2007). This has spurred the development of in silico approaches for B-cell epitope prediction. Early in silico methods typically relied on handcrafted features derived from structural and physicochemical properties, fed into machine learning models like support vector machines. For example, the DiscoTope series (Kringelum et al. 2012) used features such as solvent accessibility and contact number, while the SEPPA series (Zhou et al. 2019) added biological factors like glycosylation. However, their overall performance remained unsatisfactory (Cia et al. 2023), likely due to data scarcity and the complexity of immune recognition, which handcrafted descriptors fail to effectively capture.

The recent emergence of large protein language models (PLMs) (Hsu et al. 2022, Lin et al. 2023, Su et al. 2023), based on the transformer architecture (Vaswani et al. 2017), has markedly enriched the sequence- and structure-level representations available for downstream tasks, including epitope prediction. Leveraging these models, many recent BCE prediction approaches have achieved substantial advances. Several methods adopted the sequence-based PLM ESM-2 (Lin et al. 2023) to generate residue-level embeddings, which were then passed to classifiers for epitope prediction (Clifford et al. 2022, Israeli and Louzoun 2024). GraphBepi (Zeng et al. 2023) further incorporated structural context through graph neural networks, enhancing predictive accuracy by modeling residue–residue spatial relationships. Meanwhile, methods such as SEMA-2.0 (Ivanisenko et al. 2024) and DiscoTope-3.0 (Høie et al. 2024), both of which leverage structure-aware PLMs, have demonstrated impressive performance, reaffirming the importance of structural information in predicting conformational BCEs. Nonetheless, earlier work remained largely sequence-centric and focused on linear epitopes, since antigen structures were often unavailable at prediction time in real applications, and experimentally determined antigen–antibody complexes providing true conformational BCE for training were scarce.

These limitations are being alleviated by advances on two fronts: AlphaFold (Jumper et al. 2021), which enables accurate structure prediction directly from sequence and thus makes structure-based models widely accessible at inference time, and single-particle cryo-EM (Cheng 2018), which has accelerated the accumulation of experimentally resolved antigen–antibody complexes to expand structure-mapped epitopes. Consistent with these trends, recent studies (Zeng et al. 2023, Høie et al. 2024) report comparable performance for conformational BCE prediction when using either experimental or AlphaFold-predicted structures, further motivating structure-guided predictors.

Despite these advances, most existing methods incorporate spatial patterns by extracting informative features for each residue, but feed all features equally into the classifier during the final prediction step. This overlooks a key biological property of BCEs, their natural tendency to form spatially clustered surface regions rather than being dispersed across the whole antigen surface. Reis et al. (2022) observed that, in most antigens, the epitope surface typically forms a single connected patch, which can often be covered by an elliptical plane (Kringelum et al. 2013). These findings suggest that spatial compactness is a defining characteristic of BCEs. Ignoring this constraint often leads to fragmented and unrealistic epitope predictions, limiting the utility of current tools in downstream applications such as antibody design. Motivated by this limitation, we aim to develop a biologically consistent epitope predictor that enforces spatial coherence and harnesses modern deep learning techniques to enhance predictive accuracy and structural awareness.

In this work, we present a region-oriented B-cell epitope predictor (RoBep), which introduces a novel region constraint mechanism to promote spatially coherent predictions. RoBep integrates multiple deep learning components, including the sequence-based PLM ESM-Cambrian (ESM-C) (Hayes et al. 2025) and enhanced Equivariant Graph Neural Networks (EGNN) (Satorras et al. 2021), to effectively capture both sequence and structural features of antigens. Beyond providing residue-level epitope predictions, RoBep is also capable of identifying high-likelihood antibody–antigen interface regions, offering richer biological insights. Moreover, we demonstrate through comprehensive evaluations that RoBep exhibits high specificity on antigenic proteins, further supporting its practical relevance. Our model not only ensures more realistic spatial distributions of predicted epitope residues but also achieves the best performance on curated benchmark datasets, underscoring its utility for both accurate prediction and downstream immunological applications.

2 Materials and methods

2.1 Dataset

We constructed a high-quality antigen dataset for structure-based B-cell epitope prediction following the strategy of Shashkova et al. (2022). All protein complexes from the PDB (Berman et al. 2000) with a resolution smaller than 4.0 Å (as of 4 November 2024) were retrieved. Antibody–antigen (Ab–Ag) complexes containing heavy or light chains were identified using ANARCI (Dunbar and Deane 2016), and only those with both chains present were retained. In each complex, the remaining protein chains were designated as candidate antigen chains and processed independently, including multi-chain antigens. All chains with sequence lengths outside 12–2046 residues were removed. For each retained antigen, a residue was labeled as an epitope if any of its heavy atoms were within 4.0 Å of a heavy atom from an antibody, as defined by Kringelum et al. (2012). In addition, antigens with fewer than five epitope residues were excluded. To reduce redundancy, antigen sequences were then clustered at 95% identity using CD-HIT (Li and Godzik 2006), and within each cluster, sequences were aligned using MAFFT (Katoh et al. 2002). The longest sequence in each cluster was selected as the representative, and epitope annotations from other members were propagated based on aligned positions. This merging step prevents highly similar antigen sequences from carrying discrepant epitope labels into the training data, thereby avoiding harmful supervision conflicts. The curated dataset consisted of 756 antigens comprising 176 846 residues, of which 11 848 (6.7%) were labeled as epitopes, including both linear and conformational BCEs. Finally, the dataset was split into training and test sets by re-running CD-HIT on the non-redundant set at a 70% identity threshold and assigning entire clusters to the two splits to approximate an 8:2 ratio by sequence count, with all member sequences of each assigned cluster included in that split. We additionally double-checked the split with CD-HIT-2D at 70% identity to confirm that no cross-split pairs reach $[eqn]$ identity. This procedure yielded 612 antigens for training (140 106 residues, 6.6% epitopes) and 144 for testing (36 740 residues, 7.0% epitopes), ensuring that no antigen in the training set shares $[eqn]$ sequence identity with any antigen in the test set, thus avoiding information leakage.

Meanwhile, for the protein-level antigenicity evaluation in Section 3.5, we constructed an auxiliary antigenicity dataset comprising 144 antigens (positives) from the test set and 85 non-structural proteins from viruses (negatives). This choice follows the fact that non-structural proteins (nsps) generally exhibit low antigenicity and are rarely recognized by antibodies Islam et al. (2023). This dataset is used solely to assess RoBep’s ability to distinguish antigenic proteins from weakly antigenic nsps. Further details are provided in Supplementary Materials S1, available as supplementary data at Bioinformatics online.

2.2 RoBep framework

The overall architecture of RoBep is illustrated in Fig. 1. Given an antigen structure, RoBep first enumerates candidate Ab-Ag binding regions using a rolling-sphere strategy. Each region is modeled as a residue-level graph, with node features combining ESM-C embeddings and structural descriptors, and edge features capturing spatial and sequential relationships. Two layers of E(n)-EGNN encode these graphs to produce structure-aware residue representations. Region-level scores are then computed via attention pooling, and the top-ranked regions are selected. Residue-level epitope prediction is performed by integrating residue and region embeddings through a multilayer perceptron module (MLP), with the final probability for each residue obtained by averaging its predictions across all regions in which it appears. All predicted epitope residues lie entirely within the selected regions, a design referred to as the region constraint mechanism, which encourages spatial coherence. Detailed descriptions of each component are provided in the following subsections.

Overview of the RoBep framework for B-cell epitope prediction on a single antigen. The model identifies candidate surface regions using a rolling-sphere strategy, encodes each region as a residue-level graph with both sequence-based and structure-based features, and applies an E(n)-equivariant GNN to extract residue representations. High-confidence regions are scored and selected for downstream residue-level probability prediction, which integrates both local and regional information. The final epitope probability for each residue is obtained by averaging its predictions across all selected regions.

2.2.1 Region enumeration

To support the region constraint mechanism, we first enumerate candidate surface regions that may serve as potential Ab–Ag binding regions. We employ a sphere-based scanning strategy, a method conceptually related to the classical rolling-sphere view of solvent exposure (Lee and Richards 1971): the $[eqn]$ atoms of all solvent-accessible residues with relative solvent accessibility (RSA) ≥ 20% are selected as sphere centers, and spheres with a fixed radius of $[eqn]$  Å are constructed around each residue.

Each candidate region is defined as the set of residues whose $[eqn]$ atoms fall within the intersection of the sphere and the antigen surface. The radius is chosen based on Ag-Ab interface statistics (mean area $[eqn]$ Å^2^; Reis et al. 2022), ensuring adequate coverage for most epitopes. All these candidate regions are subsequently used for scoring and residue-level prediction. To enforce spatial coherence, final residue predictions are constrained to lie within selected high-scoring regions, forming the basis of our region constraint mechanism.

2.2.2 Geometric graph encoding

Each candidate region is represented as a fully connected residue-level graph $[eqn]$ , where nodes correspond to residues and edges connect all residue pairs within the region. Node features are composed of (i) projected ESM-C embeddings, derived from the 2560-dimensional pretrained output of ESM-C (Hayes et al. 2025), mapped to the hidden dimension via an MLP, and (ii) structural descriptors, including RSA and backbone dihedral angles. Edge features include (i) spatial distances, transformed from the Euclidean distance between $[eqn]$ atoms using 16 Gaussian radial basis functions, and (ii) sequence offsets, encoded as 16-dimensional sinusoidal positional embeddings following Vaswani et al. (2017). For a given graph, we denote the initial node feature of residue i as $[eqn]$ , and the edge feature between residue i and j as $[eqn]$ . Details of node and edge feature construction are provided in the Supplementary Materials S2, available as supplementary data at Bioinformatics online.

For simplicity, we use $[eqn]$ to denote an MLP, where each instance (e.g. $[eqn]$ ) has its own learnable parameters. Unless otherwise specified, MLPs are composed of linear layers with activation functions, followed by dropout and normalization for regularization and stability.

To capture the geometric context of each region, RoBep applies L layers of EGNNs (Satorras et al. 2021), enhanced with architectural improvements inspired by Luo et al. (2025). At each layer l, each residue i is associated with a hidden representation $[eqn]$ and a 3D coordinate $[eqn]$ . Let $[eqn]$ denote the set of neighboring residues within the region (all other nodes). For every residue pair $[eqn]$ , the message is computed based on the current node features, coordinates, and edge attributes:

$$[eqn]$$

where $[eqn]$ is an MLP. The 3D coordinate of residue i is then updated via a weighted sum of relative displacements:

$$[eqn]$$

To update node features, the aggregated messages are fused with the current representation through two MLPs:

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

where $[eqn]$ is the node embeddings of residue i at layer l and all $[eqn]$ denote MLPs with separate parameters. The node embeddings $[eqn]$ produced by the last EGNN layer are used for both region- and residue-level prediction tasks in subsequent steps.

2.2.3 Region scoring and ranking

To predict the proportion of epitopes $[eqn]$ in each candidate region j, we apply attention-based pooling over the set of residue embeddings $[eqn]$ in this region to obtain a global region representation $[eqn]$ :

$$[eqn]$$

where $[eqn]$ are attention weights and $[eqn]$ is a learnable vector. The pooled vector is then transformed into a region-level proportion:

$$[eqn]$$

All candidate regions are ranked by $[eqn]$ , and the top-7 regions, $[eqn]$ , are retained for residue-level prediction. This step enforces the Region Constraint by restricting residue predictions to these selected regions, ensuring spatial coherence and avoiding scattered outputs across the antigen surface. The selection of the value of k (number of top-ranked regions) is detailed in Supplementary Materials S3, available as supplementary data at Bioinformatics online.

2.2.4 Epitope prediction

For each residue i within a selected region j, we concatenate its embedding $[eqn]$ with the region-level representation $[eqn]$ , and compute the epitope probability as:

$$[eqn]$$

This design allows residue-level predictions to incorporate both local features and global region context. For residues appearing in multiple selected regions, we average their predictions:

$$[eqn]$$

where $[eqn]$ is the number of selected regions that include residue i and $[eqn]$ is the set of residues in $[eqn]$ . Residues not included in any selected region are assigned zero.

2.3 Objective function

To effectively train both region- and residue-level predictions, we adopt a joint loss function with dynamic weights:

$$[eqn]$$

where $[eqn]$ and $[eqn]$ are loss weights dynamically adjusted via GradNorm (Chen et al. 2018), which is detailed in Supplementary Materials S4, available as supplementary data at Bioinformatics online. $[eqn]$ controls the strength of the consistency term. Each loss component is defined below.

Region-level loss: For each candidate region j, we supervise the predicted proportion $[eqn]$ of epitopes in the region using a target-weighted mean squared error (MSE):

$$[eqn]$$

where $[eqn]$ is the fraction of true-epitope residues in region j, and $[eqn]$ is a weight that increases with $[eqn]$ encouraging accurate prediction of highly enriched regions.

Residue-level loss: For each residue i in region j, the predicted epitope probability $[eqn]$ is supervised with focal loss:

$$[eqn]$$

where $[eqn]$ denotes the set of residues in region j, and $[eqn]$ if $[eqn]$ , otherwise $[eqn]$ . This formulation encourages the model to focus on sparse and difficult-to-identify epitopes.

Cross-level consistency: To align region- and residue-level predictions, we include a consistency term:

$$[eqn]$$

which aligns each region’s proportion with its residue-level scores, encouraging mutual consistency. This multi-level loss provides complementary supervision across levels, and reinforces region–residue consistency, improving prediction accuracy.

2.4 Antigenicity prediction

For any given protein k, we define its antigenicity score as the average predicted epitope probability across residues within its top-7 predicted regions. Formally, the score is computed as:

$$[eqn]$$

where $[eqn]$ denotes the union of residues in the top-7 $[eqn]$ selected in Section 2.2.1 and $[eqn]$ denotes the predicted probability of residue i being an epitope in Section 2.2.4. This score predicts the antigenicity of a given protein.

2.5 Evaluation metrics and experiment setting

To comprehensively evaluate our method against baselines, we report five metrics: F1, Matthews correlation coefficient (MCC), area under the precision–recall curve (AUPR), area under the ROC curve at low false-positive rates ( $[eqn]$ ) (Richardson et al. 2024), and an epitope–overlap metric called the antigen-level intersection over union (AgIoU) (You et al. 2025). All these metrics can capture model performance well in scenarios with significant class imbalance, given that epitope residues constitute only 6.7% of the entire dataset.

To further evaluate the spatial coherence of predicted epitopes, we calculate the clustered mean pairwise distance (cMPD) between $[eqn]$ atoms of predicted residues for each antigen, a metric modified from Colavin et al. (2022) that provides a fair and accurate assessment of residues’ compactness, including multi-region cases. Lower cMPD values indicate tighter clustering of predicted residues, reflecting higher biologically plausible prediction results. We also compare the predicted cluster count (NoC), the number of distinct epitope regions identified by each model, to the ground-truth NoC to assess regional correspondence. Details of metrics are provided in Supplementary Materials S5, available as supplementary data at Bioinformatics online.

Models were trained with a batch size of 64 for up to 50 epochs, using early stopping with a patience of 10. We used the AdamW optimizer with a maximum learning rate of $[eqn]$ , a 10% linear warm-up, and cosine annealing with restarts. All experiments were run on two NVIDIA Quadro GV100 GPUs. More training implementation details are provided in Supplementary Materials S6, available as supplementary data at Bioinformatics online.

3 Results

3.1 Benchmark comparison

We compared our method, RoBep, with popular state-of-the-art methods for structure-based prediction of conformational B-cell epitopes on the test set. These methods leverage structural information through various strategies, including handcrafted physicochemical features (e.g., Seppa-3.0), structure-based PLNs (e.g. SEMA-2.0, DiscoTope-3.0), and GNN (e.g. GraphiBepi). For CALIBER, we used the best-performing reproducible configuration (ESM-2 + BiLSTM) reported in the original paper of Israeli and Louzoun (2024). To ensure a fair comparison, we select the decision threshold for each method by sweeping $[eqn]$ (step 0.05) to maximize its F1, and then fix that $[eqn]$ for all thresholded metrics.

Since some existing methods, such as GraphiBepi and CALIBER, adopt earlier-generation sequence-level PLMs like ESM-2 as an encoder, we also evaluated our model using the same encoder to ensure a fair comparison. As summarized in Table 1, our model still consistently outperformed all baselines under this setting, with improvements of 10.5% in AUPRC, 15.3% in F1 score, 18.1% in AgIoU, and 27.3% in $[eqn]$ over the strongest competitor (DiscoTope-3.0). Figure 2 further illustrates this advantage by showing that our method achieves the highest precision across almost the entire recall range. These results highlight the effectiveness of our region-constrained design, which not only enhances biological plausibility but also enhances epitope localization accuracy.

Precision–recall (AUPR) curves on the external test set. RoBep equipped with either ESM-2 or ESM-C, consistently achieves higher precision across all recall levels.

Moreover, replacing the PLM with ESM-C, a model with more parameters and a more advanced architecture, led to further gains over the ESM-2-based version of RoBep: AUPRC increased by 0.9%, F1 score by 9.5%, AgIoU by 12.0%, and $[eqn]$ by 12.2%. These improvements are consistent with ESM-C’s enhanced sequence modeling capacity and further demonstrate that our model can effectively leverage high-quality embeddings to boost predictive performance.

Finally, the notably weaker performance of Seppa-3.0, which does not utilize any PLM, underscores the importance of modern PLMs in capturing informative sequence or structure-derived representations.

3.2 Spatial compactness of predicted epitope residues

In this section, we evaluate the spatial compactness of epitopes predicted by each model, an important characteristic reflecting the biological plausibility and practical significance of predicted epitopes. To quantify compactness, we computed the clustered mean pairwise distance (cMPD) of predicted residues across each antigen in the independent test dataset, as detailed in Section 2.5. As shown in Fig. 3, RoBep demonstrated significantly tighter spatial clustering compared to all baseline methods ( $[eqn]$ , one-sided Dunnett’s test). Specifically, RoBep achieves an averaged cMPD of 18.49 Å with a standard deviation of 3.38 Å, closest to the distribution of true-epitope cMPD (13.62 Å $[eqn]$ Å; red dashed line). In comparison, the averaged cMPDs for the baseline models are: CALIBER (21.91 Å), GraphiBepi (25.40 Å), SEMA-2.0 (25.91 Å), and DiscoTope-3.0 (28.02 Å).

Comparison of cMPD and NoC among predicted epitopes and true epitopes on the independent test dataset. Lower cMPD values indicate tighter clustering. Asterisks denote statistical significance from one-sided Dunnett’s test against RoBep. The horizontal dashed line indicates the mean cMPD of true epitopes across all test samples.

Beyond compactness, we also compare the number of clusters (NoC) implied by each model’s predictions to that of true epitopes. It can be found in Fig. 3 that RoBep’s NoC is the closest to the ground truth, indicating preservation of epitope regional organization (avoiding both over-fragmentation into many small patches and over-merging into a single cluster). By contrast, SEMA-2.0 tends to produce one large connected cluster (NoC $[eqn]$ ) with high cMPD, whereas other baselines tend to over-fragment epitopes. Together with the lower cMPD, these results demonstrate the effectiveness of our approach in producing spatially coherent predictions that align with biologically realistic epitope distributions.

3.3 Effect of region constraint and other components

To assess the contribution of region constraint and other modules in our architecture to BCE prediction, we conducted a series of ablation experiments. First, to evaluate the impact of region constraint, we removed the region enumeration and region ranking module to perform epitope prediction over the entire antigen, as in most existing methods. In this setting, the graph construction was changed from a region-based complete graph to a radius-based graph (18 Å cutoff) to maintain compatibility with EGNN. As visualized in Fig. 4, removing the region constraint led to a substantial performance drop: F1 score fell from 0.3794 to 0.3247, MCC from 0.3398 to 0.2357, and AUPRC from 0.2676 to 0.2520. These decreases underscore that focusing predictions on a few high-likelihood surface patches both preserves biological consistency and significantly benefits BCE prediction.

Performance comparison on the test dataset when removing the Region Constraint or other individual components.

Next, omitting the Geometric Graph Encoding module caused notable declines in all metrics (F1 –8.2%, MCC –10.7%, AUPRC –13.2%), demonstrating its critical role in capturing spatial patterns. Finally, removing structural features, RSA, and backbone torsion angles resulted in a modest F1 drop of 2.0%, suggesting these descriptors provide complementary information, although some of their signal may already be captured by pretrained embeddings and geometric graph encoders.

3.4 Impact of antigen species and length

Given the considerable diversity arising from antigen species and their variant lengths, we investigated how RoBep’s predictive performance varies across different antigen species and sequence lengths. Specifically, antigens were classified into three groups: (i) human (autoantigens, tumor antigens, etc.), (ii) viral (all viral antigens), and (iii) other (including artificial proteins, Arrestin, etc.). We further divided them into five distinct length bins: 0–100, 100–200, 200–400, 400–700, and 700+ aa, for a more granular categorization based on antigen lengths. We computed mean AUPR per stratum as well as the number of non-redundant antigens in each stratum to reflect investigation density.

As shown in Fig. 5, RoBep exhibits clear trends with respect to antigen length and species. First, although the 0–100 aa stratum contains fewer antigens than longer strata, it achieves the best performance (viral: AUPR = 0.569; human: AUPR = 0.491), and the mean AUPR decreases as length increases. This suggests that epitopes on shorter antigens are intrinsically easier to localize and may also reflect that representations learned on longer sequences transfer to shorter ones. Meanwhile, viral antigens consistently yielded higher AUPR values compared to human and other categories, even though human antigens are more numerous overall than viral (421 versus 220 in our non-redundant set). It is likely attributed to the fact that several viral antigens [e.g., influenza A (Grandea et al. 2010) and coronaviruses (Kapingidza et al. 2023)] contain relatively conserved epitope regions, whereas human antigens including autoantigens (Vanderlugt and Miller 1996) and tumor antigens (Blass and Ott, 2021) are more heterogeneous across individuals and disease contexts. In addition to biological conservation, species-specific performance differences may also be influenced by the degree of antibody characterization per antigen, as viral antigens tend to be recognized by a more diverse set of antibodies than human antigens (see Supplementary Materials S8, available as supplementary data at Bioinformatics online). Notably, an intriguing performance rebound is observed for viral antigens in the 700+ aa bin (AUPR = 0.491), likely driven by a higher concentration of well-studied viral antigens (e.g. spike-like proteins) that dominate this stratum. The overall trend indicates that shorter and viral antigen sequences may inherently contain epitope information that is easier to predict. Further discussion on the factors contributing to species-specific performance differences is provided in Supplementary Materials S8, available as supplementary data at Bioinformatics online.

Mean AUPR of RoBep across species–length strata in the test set, with counts of non-redundant antigens per stratum.

3.5 Model specificity on antigenic proteins

Identifying whether a given protein contains antigenic epitopes typically requires resource-intensive in vivo assays. Therefore, in silico methods capable of specifically recognizing antigenic proteins can greatly facilitate the initial screening of potential antigenic proteins. To evaluate RoBep’s specificity toward antigenic proteins, we examined its ability to distinguish known antigenic proteins from those relatively lacking antigenicity. Specifically, we conducted a binary classification on our antigenicity dataset (positives: antigens; negatives: viral nsps), using the protein-level antigenicity score defined in Section 2.4. The decision threshold (0.525) was selected on the training set (positives: antigens; negatives: their corresponding antibodies) to avoid data leakage and then fixed for evaluation on the antigenicity dataset introduced in Section 2.1.

As shown in Fig. 6, RoBep assigned significantly lower antigenicity scores to most non-structural proteins, confirming its capacity to suppress false positives on weakly antigenic proteins. Our model also achieved strong discriminative performance, with AUPR = 0.9930, AUROC = 0.9905, F1 = 0.7500, and Recall =0.6042. It can be found that the model’s ability to identify antigenic proteins is weaker than that for non-antigenic proteins, yielding a lower recall, consistent with the intrinsic difficulty of generalizing across diverse antigen structures under limited data availability. Overall, these results confirm that RoBep effectively assigns lower antigenicity scores to proteins without antigenic properties, underscoring its practical value in specifically identifying genuine antigenic targets. More evaluation results on the Ab–Ag dataset and comparison with baseline models are provided in Supplementary Materials S9, available as supplementary data at Bioinformatics online.

Predicted antigenicity of antigen and non-structural protein. The red dashed line indicates the threshold of 0.525.

3.6 Case study

In this section, we present a representative case study on an antigen from the test set (PDB ID: 5i9q, Chain A) to further demonstrate the strength of our RoBep framework and the contribution of the region constraint mechanism. RoBep is compared with its variant without region constraint, as well as with four representative baselines: GraphBepi, SEMA-2.0, DiscoTope-3.0, and CALIBER.

As shown in Fig. 7, RoBep achieves the highest performance across all metrics, with an F1 score of 0.88, AUPR of 0.75, and precision of 0.86. Notably, even without the region constraint mechanism, our model maintains competitive performance, confirming the robustness of its underlying architecture. However, this unconstrained variant exhibits a notable increase in false positives (red), underscoring the value of the region constraint in suppressing noisy predictions and enforcing spatial focus.

Visual comparison of predicted B-cell epitopes on antigen 5I9Q (chain A) by six methods: (a) RoBep, (b) RoBep without Region Constraint, (c) GraphBepi, (d) SEMA-2.0, (e) DiscoTope-3.0, and (f) CALIBER. The grid sphere for (a) indicates the region with the highest predicted likelihood of being an antibody–antigen interface, as identified by RoBep.

In comparison, baseline methods such as SEMA-2.0, DiscoTope-3.0, and CALIBER generate predictions that are more scattered across the antigen surface, despite identifying a comparable number of true-epitope residues. This results in substantially lower precision (0.38, 0.53, and 0.36, respectively) and weaker spatial coherence. By contrast, RoBep produces highly spatially clustered predictions, closely aligning with how B-cell epitopes are typically present in nature, as compact regions. This spatial compactness not only improves biological plausibility but also provides more actionable guidance for downstream tasks such as vaccine design and antibody development. More examples with different performance can be found in the Supplementary Materials S10, available as supplementary data at Bioinformatics online.

4 Discussion

In this study, we introduced RoBep, a region-constrained framework for B-cell epitope prediction. By explicitly modeling the spatial clustering of epitopes through a region constraint mechanism, RoBep generates biologically coherent predictions and effectively avoids the overly scattered outputs commonly observed in existing residue-level models. The framework integrates ESM-C and EGNN to jointly capture informative sequence and structural representations. In addition to providing residue-level epitope probabilities, RoBep also predicts high-likelihood Ab–Ag binding regions, making it a practical and interpretable tool for structure-guided vaccine design and antibody discovery. Notably, RoBep demonstrates a degree of specificity on antigenic proteins and achieves consistently strong performance in benchmark comparison. Note that RoBep can also predict linear BCEs since linear B-cell epitope residues are also clustered on an antigen surface region.

Despite its advantages, RoBep has certain limitations. First, its region enumeration strategy is based on a fixed-radius spherical scanning approach, which may not well adapt to antigens of varying size, shape, or domain architecture. Future work may explore more adaptive or data-driven region selection strategies to improve coverage and generalizability. Second, while the region constraint enforces spatial coherence, it introduces a high-risk, high-reward trade-off: if all top-ranked regions miss the true-epitope location, the model has limited capacity to recover correct predictions. This limitation suggests the possibilities of incorporating hierarchical or soft region constraints to improve robustness without sacrificing spatial focus, allowing the model to make more flexible and resilient predictions.

In summary, RoBep is the first biologically grounded, structure-aware framework to incorporate explicit spatial constraints into B-cell epitope prediction. By combining advanced deep learning models with a region constraint mechanism, RoBep not only achieves accurate and biologically plausible residue-level outputs but also introduces region-level BCE identification for the first time. Beyond epitope prediction, we believe the region constraint strategy may serve as a general design principle for other structure-based localization tasks where spatial coherence and interpretability are essential.

Supplementary Material

btag006_Supplementary_Data