Beyond Message Passing: A Symbolic Alternative for Expressive and Interpretable Graph Learning
Authors: Chuqin Geng, Li Zhang, Haolin Ye, Ziyu Zhao, Yuhe Jiang, Tara Saba, Xinyu Wang, Xujie Si
Organizations: School of Computer Science, McGill University, Montreal, Canada · Department of Computer Science, University of Toronto, Toronto, Canada
Graph Neural Networks (GNNs) have become essential in high-stakes domains such as drug discovery, yet their black-box nature remains a significant barrier to trustworthiness. While self-explainable GNNs attempt to bridge this gap, they often rely on standard message-passing backbones that inherit fundamental limitations, including the 1-Weisfeiler-Lehman (1-WL) expressivity barrier and a lack of fine-grained interpretability. To address these challenges, we propose SymGraph, a symbolic framework designed to transcend these constraints. By replacing continuous message passing with discrete structural hashing and topological role-based aggregation, our architecture theoretically surpasses the 1-WL barrier, achieving superior expressiveness without the overhead of differentiable optimization. Extensive empirical evaluations demonstrate that SymGraph achieves state-of-the-art performance, outperforming existing self-explainable GNNs. Notably, SymGraph delivers 10x to 100x speedups in training time using only CPU execution. Furthermore, SymGraph generates rules with superior semantic granularity compared to existing rule-based methods, offering great potential for scientific discovery and explainable AI.
Figures & tables
Figure 1: Visualization of a Local Structural Predicate. The extracted predicate identifies specific node and edge constraints, providing superior semantic granularity compared to canonical subgraph-based explanations. See more details in Appendix H .
Figure 2: An overview of the SymGraph framework.
Ba2Motifs
BaMultiShapes
PROTEINS
Mutagenicity
NCI1
BBBP
Model
Acc ↑
Time ↓
Acc ↑
Time ↓
Acc ↑
Time ↓
Acc ↑
Time ↓
Acc ↑
Time ↓
Acc ↑
Time ↓
GCN
58.3 ± 2.3
0.2 ± 0.0
70.7 ± 9.7
0.5 ± 0.0
71.1 ± 1.5
2.8 ± 0.4
81.6 ± 1.0
14.4 ± 1.6
76.0 ± 1.3
12.6 ± 1.1
86.1 ± 1.0
5.3 ± 1.0
GraphSAGE
48.2 ± 0.8
0.3 ± 0.0
50.6 ± 2.1
0.3 ± 0.0
72.3 ± 2.3
2.2 ± 1.1
81.9 ± 1.6
13.9 ± 3.2
79.6 ± 0.3
12.9 ± 1.4
86.9 ± 1.0
4.6 ± 1.7
GAT
48.2 ± 0.8
0.3 ± 0.0
49.4 ± 2.1
0.3 ± 0.0
70.9 ± 1.0
1.8 ± 1.1
82.2 ± 1.2
18.2 ± 1.9
78.0 ± 2.0
16.4 ± 4.3
88.2 ± 2.0
6.9 ± 1.0
GIN
100.0 ± 0.0
0.3 ± 0.0
95.0 ± 1.5
0.4 ± 0.2
68.5 ± 2.3
3.0 ± 0.4
83.0 ± 1.6
13.1 ± 1.4
81.2 ± 1.3
11.9 ± 2.2
87.7 ± 1.0
5.5 ± 0.6
GNAN
49.5 ± 0.5
51.8 ± 1.2
51.4 ± 1.6
50.4 ± 1.2
59.6 ± 0.0
65.2 ± 1.8
55.4 ± 0.1
163.6 ± 5.7
52.1 ± 1.7
326.7 ± 13.1
23.5 ± 0.0
81.8 ± 3.4
Table 1: Performance Benchmark. Comparison of test accuracy (%) and total runtime ( 102 seconds) against state-of-the-art self-explainable GNNs and standard message-passing GNN baselines. Results are averaged over five random seeds. Bold and underlined indicate the best and second-best performance, respectively. SymGraph -DT and SymGraph -RF denote the decision tree and random forest variants of our framework. Notably, our approach achieves superior accuracy on complex molecular datasets while being orders of magnitude faster than competing self-explainable methods like LogiX-GIN and PiGNN . See more results in Appendix E .
Figure 4: SymGraph identifies expert-validated chemical motifs that provide ground-truth explanations for molecular properties. Our predicates enforce structural constraints on node and edge attributes, which are readily mapped to SMARTS patterns, as illustrated by the red highlights on the chemical compounds. This provides a direct bridge between learned logic and established domain knowledge.
Model Variant
PROTEINS
Mutagenicity
NCI1
BBBP
Ablation 1: Topological Role Encoding via Orbits (Orbits)
DT w/o Orbits
70.0 ± 3.4
78.4 ± 0.9
73.0 ± 1.2
84.4 ± 1.4
RF w/o Orbits
75.3 ± 2.6
81.6 ± 0.7
81.5 ± 1.1
86.7 ± 1.3
Ablation 2: Structural Histograms (Counts)
DT w/o Counts
66.7 ± 4.7
80.1 ± 0.9
74.2 ± 1.0
84.6 ± 1.4
RF w/o Counts
75.2 ± 1.4
83.4 ± 1.7
82.9 ± 1.7
88.5 ± 1.2
Table 2: Ablation Study. Impact of removing key components from SymGraph . Results are reported as mean ± standard deviation over 5 random seeds. The consistent performance degradation across all variants validates the design principles of SymGraph .
Appendix figures & tables10 assets
Supplementary material from the paper’s appendix.
Appendix
Ba2Motifs
BAMultiShapes
Mutagenicity
BBBP
NCI1
PROTEINS
#Graphs
1,000
1,000
4,337
2,050
4,110
1,113
Avg. ∣V∣
25.0
40.0
30.3
23.9
29.9
39.1
Avg. ∣E∣
51.0
87.0
30.8
51.6
32.3
145.6
Appendix
Table 3: Statistics of standard graph classification datasets used in the main paper.
MUTAG
BACE
HIV
IMDB-BINARY
TreeGrid
BaShapes
#Graphs
188
1,513
1,484
1,000
1
1
Avg. ∣V∣
17.9
34.1
26.1
19.8
1,231
700
Avg. ∣E∣
39.6
73.7
55.5
193.1
1,550
2,105
Appendix
Table 4: Statistics of additional graph and node classification datasets used in the supplementary material. Note that TreeGrid and BaShapes are single-graph node classification tasks, so ∣V∣ and ∣E∣ refer to the total size of the single graph.
Dataset
K_HOPS
TREE_MODE
MAX_LEAVES_CAP
POPULATION_SIZE
GENERATION
JUDGE_MAX_LEAVES
JUDGE_CCP_ALPHA
Ba2Motifs
1
’LEAVES’
—
—
—
—
—
BaMultiShapes
1
’LEAVES’
—
—
—
—
—
PROTEINS
1
’LEAVES’
16
100
5
48
0.001
Mutagenicity
2
’LEAVES’
16
100
5
48
0.001
NCI1
2
’LEAVES’
16
100
5
48
0.001
BBBP
2
’LEAVES’
16
100
5
32
0.001
Appendix
Table 5: Hyperparameter configuration for SymGraph models.
Figure 5: Raw subgraph visualization via Node IDs. These figures display partially the raw output from the LogiX-GIN extraction process. The highlighted regions correspond to the internal node indices (IDs) identified in the DNF rules as critical for classification, prior to semantic mapping.
Figure 6: Semantically interpreted visualization with atom identities. The corresponding molecular visualization of the graphs in Figure 5 . By mapping internal IDs to atomic identities and manual filtering, we collect these highlighted subgraphs that are easier to interpret.
MUTAG
BACE
ClinTox
IMDB
TreeGrid
BaShapes
Model
Acc ↑
Time ↓
Acc ↑
Time ↓
Acc ↑
Time ↓
Acc ↑
Time ↓
Acc ↑
Time ↓
Acc ↑
Time ↓
GCN
74.7 ± 6.1
0.5 ± 0.2
80.7 ± 1.7
4.8 ± 1.1
92.3 ± 0.6
8.0 ± 4.6
49.4 ± 2.1
0.3 ± 0.0
71.6 ± 1.1
44.1 ± 0.8
85.5 ± 2.0
34.5 ± 1.7
GIN
88.4 ± 3.0
0.4 ± 0.0
82.4 ± 2.7
3.8 ± 1.0
92.3 ± 0.7
7.8 ± 2.4
69.5 ± 1.9
0.7 ± 0.2
82.6 ± 0.2
34.1 ± 4.1
96.9 ± 1.0
27.6 ± 1.3
GraphSAGE
75.3 ± 4.8
0.4 ± 0.1
82.6 ± 2.2
4.4 ± 0.8
92.2 ± 1.0
8.2 ± 0.5
48.6 ± 1.5
0.2 ± 0.0
41.5 ± 2.4
24.9 ± 1.6
42.9 ± 1.7
21.5 ± 0.2
GAT
75.3 ± 3.5
0.5 ± 0.2
81.7 ± 1.1
6.3 ± 1.0
92.6 ± 0.9
9.4 ± 3.8
50.6 ± 2.1
0.3 ± 0.0
41.5 ± 2.4
32.6 ± 2.0
42.9 ± 1.7
28.7 ± 1.3
KerGNN
69.0 ± 2.2
19.5 ± 0.4
57.0 ± 2.7
91.8 ± 0.1
92.4 ± 1.0
10.3 ± 1.0
51.3 ± 1.6
65.9 ± 0.2
—
—
—
—
Appendix
Table 6: Test accuracy (%) and runtime ( 102 seconds) on additional benchmarks for XAI graph and node classification. Results are averaged over five random seeds. Bold indicates the best performance, and underlined indicates the second-best. SymGraph -DT and SymGraph -RF denote decision tree and random forest variants of our framework, respectively. — indicates that the method is not supported.
Figure 7: Visualization of logic rules induced by local structural predicates. The symbolic rule (Orbit 0: #O=1)∧(Orbit 1: #(C-C)≥2) is grounded to a specific molecular substructure. By decomposing the motif into Orbit 0 (node constraint) and Orbit 1 (edge constraint), SymGraph achieves a level of semantic granularity that exceeds traditional coarse-grained subgraph explanations.
Aspect
SMARTS Rules (Chemical Standard)
SymGraph ’s Rules
Structural Constraints
Atom and bond constraints at specific positions (e.g., [O] , [C]-[O] )
Orbit-wise node and edge constraints (e.g., #O=1 , #(C,C)≥2 )
Topological Roles
Defined by bond positions and connectivity
Orbit indices encoding structural symmetry via hashing
Logical Flexibility
Logical OR over pattern variants
Disjunctive Normal Form (DNF) clauses
Compositionality
Recursive conjunctions of conditions
Conjunctions of structural predicates
Generalization
Matches molecules not present in training data
Covers structural families beyond observed motifs
Implementation
Hand-crafted or mined patterns
Automatically induced via RDKit-compatible predicates
Appendix
Table 7: Conceptual correspondence between chemical SMARTS patterns and SymGraph ’s rules.
Figure 8: Visualization of the structural predicates identified by SymGraph for the BAMultiShapes dataset.
Figure 9: Visualization of the structural predicates identified by SymGraph for the BA2Motifs dataset.
Graph Neural Networks (GNNs) face fundamental limitations in expressivity and capturing structural heterogeneity. Standard message-passing architectures are constrained by the 1-dimensional Weisfeiler-Leman (1-WL) test, unable to distinguish graphs beyond degree sequences, and aggregate information uniformly from neighbors, failing to capture how nodes occupy different structural positions within higher-order patterns. While methods exist to achieve higher expressivity, they incur prohibitive computational costs and lack unified frameworks for flexibly encoding diverse structural properties. To address these limitations, we introduce Invariant-Stratified Propagation (ISP), a framework comprising both a novel WL variant (ISP-WL) and its efficient neural network implementation (ISPGNN). ISP stratifies nodes according to graph invariants, processing them in hierarchical strata that reveal structural distinctions invisible to 1-WL. Through hierarchical structural heterogeneity encoding, ISP quantifies differences in nodes' structural positions within higher-order patterns, distinguishing interactions where participants occupy different roles from those with uniform participation. We provide formal theoretical analysis establishing enhanced expressivity beyond 1-WL, convergence guarantees, and inherent resistance to oversmoothing. Extensive experiments across graph classification, node classification, and influence estimation demonstrate consistent improvements over standard architectures and state-of-the-art expressive baselines.
Asela Hevapathige, Ahad N. Zehmakan, Asiri Wijesinghe +1
Department of Mechanical Engineering, University of Melbourne, Melbourne, Australia · School of Computing, Australian National University, Canberra, Australia · Data61, CSIRO, Canberra, Australia
Graph Neural Networks (GNNs) excel at relational reasoning but face two persistent challenges: the lack of interpretable attribution for heterogeneous node types, and the computational overhead of message passing over large, noisy graphs. We propose the Hierarchical Attention-based Heterogeneous GNN (HA-HeteroGNN), a framework that addresses both issues through a unied explainability-to-pruning pipeline. A two-tier attention mechanism separates sensor-level and context-level computation across 16 node types and 18 edge types, producing per-node relevance scores via an attention-based GNN Explainer without requiring gradient backpropagation. These relevance scores then serve as a principled pruning criterion: removing nodes identied as consistently uninformative yields a 27% reduction in graph edges while simultaneously improving classication accuracy by 2.46.1% across all model variants, challenging the conventional assumption that pruning necessarily trades accuracy for eciency. Experiments on a 50,000-record synthetic dataset spanning 11 report categories demonstrate 97.5% cross-strategy explanation stability and domain consistent sensor attribution, with training-time reductions of up to 43.9% and real-time inference latency of approximately 5860 ms per sample.
Seungwoo Kum
Korea Electronics Technology Institute (KETI) Seongnam, Republic of Korea
We introduce Graph Concept Bottleneck (GCB) as a new paradigm for self-explainable text-attributed graph learning. GCB maps graphs into a subspace, concept bottleneck, where each concept is a meaningful phrase, and predictions are made based on the activation of these concepts. Unlike existing interpretable graph learning methods that primarily rely on subgraphs as explanations, the concept bottleneck provides a new form of interpretation. To refine the concept space, we apply the information bottleneck principle to focus on the most relevant concepts. This not only yields more concise and faithful explanations but also explicitly guides the model to "think" toward the correct decision. We empirically show that GCB achieves intrinsic interpretability with accuracy on par with black-box Graph Neural Networks. Moreover, it delivers better performance under distribution shifts and data perturbations, showing improved robustness and generalizability, benefitting from concept-guided prediction.
Xiaoxue Han, Libo Zhang, Zining Zhu +1
Department of Computer Science, Stevens Institute of Technology, Hoboken, NJ, USA.