The challenge that motivated the work is deceptively simple to state. Effector proteins are extraordinarily diverse. They often share little recognizable sequence similarity with one another, even when they perform related functions, because bacteria have evolved them repeatedly from different starting points. Traditional bioinformatic tools that scan for conserved sequence motifs therefore miss many genuine effectors while occasionally flagging innocent proteins. In recent years, pretrained protein language models, which learn statistical representations of amino acid sequences from vast databases, have substantially improved effector prediction by capturing contextual patterns that simpler methods overlook. Yet even these powerful models share a fundamental blind spot: they see proteins as one-dimensional strings of letters, while the biological activity of an effector depends on how its chain folds in three-dimensional space.

Structure, in principle, holds the key. The virulence functions of effectors, such as binding host proteins, mimicking host enzymes, or disrupting cellular compartments, are encoded in conformational motifs: recurring arrangements of residues in space that may appear at widely separated positions along the linear sequence. The problem has always been data. Experimentally resolved structures, obtained through crystallography or cryo-electron microscopy, exist for only a small fraction of bacterial proteins, and effectors are particularly underrepresented because they are difficult to purify and crystallize. Conventional structure-based prediction methods, however elegant, have been starved of the structural information they need to work at scale.

GeoEPred, developed by Shouzhen Song, Hua Shi, Hongfeng Wu, Dachen Liu, Yihang Lin, Nor Ashidi Mat Isa, Quan Zou, Leyi Wei and colleagues, threads this needle with a clever compromise. Instead of waiting for experimental structures, the framework uses ESMFold, a deep learning system that predicts three-dimensional protein structures directly from sequence, to generate structural models for every protein it examines. These predicted structures are then processed by geometric vector perceptrons, neural network components that operate on the local geometry of the protein: the orientations between neighboring residues, the distances separating them, and the topology of the spatial neighborhood each residue occupies. In effect, the model learns what local conformational neighborhoods look like in proteins that function as effectors, capturing spatial signatures that no sequence-only model can perceive.

But structure alone is not enough either, and this is where the multimodal design of GeoEPred becomes distinctive. The framework runs two parallel streams of representation. The first draws on a pretrained protein language model to produce sequence-contextual embeddings, rich numerical descriptions of each amino acid in its linear context. The second builds the geometric representations from the predicted structure. A dedicated feature projection network then refines the fine-grained sequence signals that are most relevant to effector function, filtering the language model’s general-purpose knowledge toward the specific patterns that matter for secretion and virulence.

The hardest part of any multimodal system is making the two modalities talk to each other. A sequence embedding and a structural embedding describe the same protein, but they live in different mathematical spaces with different statistics, and naive concatenation often fails to exploit their complementary strengths. The GeoEPred team addressed this with a cross-modal alignment and feature-tokenized self-attention module. The alignment component uses contrastive learning, a technique that pulls matching representations from the two modalities closer together in the learned space while pushing mismatched ones apart, thereby enforcing consistency between sequence semantics and structural geometry. The tokenized self-attention component then models associations at the level of individual residues, allowing the network to link linear functional motifs, short stretches of sequence with known or learned significance, to the spatial conformational patterns those stretches adopt in the folded protein.

The result is a model that can recognize an effector not merely because its sequence resembles known effectors, but because its residues arrange themselves in space the way effector residues tend to. That distinction matters most at the frontier of discovery, where new effectors share no detectable ancestry with characterized ones. The authors evaluated GeoEPred on multiple benchmark datasets covering the three major secretion systems of Gram-negative bacteria, known as type III, type IV, and type VI secretion systems, each of which delivers effectors in its own characteristic way. Across all three prediction tasks, abbreviated T3SE, T4SE, and T6SE in the field’s shorthand, GeoEPred outperformed existing leading models, demonstrating that the structural signal adds genuine predictive power rather than redundant information.

Perhaps more important than the headline benchmarks is the model’s behavior in remote homolog recognition scenarios, where the test proteins are evolutionarily distant from anything seen during training. This is the regime that determines whether a predictor is useful for real discovery or merely good at memorizing known families. GeoEPred maintained stable performance under these demanding conditions, suggesting that the geometric representations it learns capture functional constraints that persist across deep evolutionary time even as sequences drift beyond recognition. The modular architecture also proved extensible, and the authors report strong generalization ability with substantial application potential for genome-scale effector discovery, the scenario in which a newly sequenced bacterial genome can be scanned end to end for candidate virulence factors.

The practical implications reach well beyond the benchmark leaderboards. Rapid identification of effectors accelerates the elucidation of bacterial pathogenic mechanisms, giving microbiologists a prioritized shortlist of proteins to characterize experimentally instead of a genome full of untested hypotheses. It also informs the development of precise anti-infective strategies: effectors and the machinery that delivers them are attractive drug targets, and knowing the full effector repertoire of a pathogen reveals which host pathways it depends on. In agriculture, where Gram-negative pathogens devastate crops, the same logic applies to engineering resistance. As predicted structures from tools like ESMFold become routine inputs to biological machine learning, frameworks like GeoEPred illustrate a broader shift in computational biology, from reading genomes as text to understanding proteins as physical objects whose shapes, not just their sequences, carry the instructions for what they do.

None of this replaces experimental validation, and the authors are careful to frame GeoEPred as a discovery engine whose predictions guide, rather than conclude, the scientific process. Predicted structures, however accurate on average, carry uncertainties that propagate into any downstream model, and effector biology remains full of surprises that no current dataset fully captures. But the trajectory is clear. By fusing the semantic richness of protein language models with the physical realism of geometric deep learning, and by solving the fusion problem with principled contrastive alignment rather than brute-force concatenation, GeoEPred offers a template for how the next generation of protein function predictors will be built. For the hidden weapons of Gram-negative bacteria, many of which have evaded detection precisely because they look like nothing we have sequenced before, the era of shape-aware prediction may finally be opening the door.

Subject of Research: Multimodal geometric deep learning for predicting Gram-negative bacterial secreted effector proteins

Article Title: GeoEPred: A multimodal structure-aware geometric deep learning framework for Gram-negative bacterial secreted effector prediction with sequence semantics

Article References: GeoEPred: A multimodal structure-aware geometric deep learning framework for Gram-negative bacterial secreted effector prediction with sequence semantics. (n.d.). https://doi.org/10.1371/journal.pcbi.1014344

Image Credits: AI Generated

DOI: 10.1371/journal.pcbi.1014344

Keywords: GeoEPred, Gram-negative bacteria, effector proteins, deep learning, protein language models, ESMFold, protein structure prediction, type III secretion system, type IV secretion system, type VI secretion system, virulence, computational biology

Why does this matter? Proteins almost never work alone. They form intricate interaction networks that govern everything from immune responses to reproduction, growth, and product quality. Understanding these networks at a systems level is essential for improving traits such as disease resistance, fertility, and meat or egg quality—the cornerstones of food security and sustainable farming. Yet existing interaction databases like STRING, BioGRID, IntAct, and MINT offer only sparse coverage of agricultural species, and they typically lack the structural annotations needed to identify binding interfaces, functional hotspots, or the molecular consequences of mutations. Without such detail, linking a genetic variant to a measurable trait—an essential step in molecular breeding—remains a matter of guesswork.

The new database was built through an ambitious integrative pipeline that combines three complementary prediction strategies. The first, interolog mapping, exploits evolutionary conservation: if two proteins in a model organism are known to interact, their homologs in a farm animal are likely to interact as well. The team assembled a rigorously filtered library of experimentally supported interactions from seven model organisms, including human, mouse, fly, worm, yeast, E. coli, and Arabidopsis, retaining only high-confidence interactions scored by the HIPPIE scheme. The second strategy, domain–domain interaction inference, predicts interactions based on conserved protein domain pairings drawn from the 3did catalog and expanded with an expectation–maximization algorithm applied to curated interaction data.

The third and most modern component is a deep learning model. The researchers retrained PPITrans, a transformer-based architecture built on the ProtT5 protein language model, using a curated set of 44,035 high-confidence human interactions as positive samples and 430,370 negative pairs. Training on an NVIDIA A30 GPU with the AdamW optimizer, the model learned sequence-derived interaction patterns that generalize across species. The scores from all three methods were then fused through a logistic regression model trained on a mouse benchmark dataset, producing a single interaction probability for every protein pair. On independent testing, the integrated model achieved an area under the ROC curve of 0.914 and an area under the precision–recall curve of 0.679, outperforming each individual component. Interactions were classified as high-confidence only above a stringent threshold calibrated at a 0.1 percent false positive rate.

The resulting species-specific networks are striking in scale. Cattle lead the pack with 173,331 predicted interactions covering 34.0 percent of the proteome, followed by pig with 90,987 interactions, chicken with 61,593, duck with 43,886, and sheep with 40,624. In total, these interactomes span proteomes at coverage rates ranging from 19.5 to 34.0 percent—a dramatic expansion of molecular data for organisms that have historically been interaction-data deserts. Cross-validation against independent resources lent credibility to the predictions: the cattle network shared 8,127 interactions with a D-SCRIPT-based dataset and captured 207 of 632 experimentally solved cattle complex structures in Interactome3D, both overlaps highly significant by Fisher’s exact test.

What truly distinguishes AASIA, however, is its structural dimension. The team used ESMFold, a cutting-edge structure prediction system powered by the ESM-2 protein language model, to generate three-dimensional models for every protein monomer across the five proteomes and, in multimer mode, for all 410,000-plus interacting pairs. Unlike earlier approaches requiring multiple sequence alignments or structural templates, ESMFold predicts structures directly from amino acid sequences, enabling proteome-scale modeling at unprecedented speed. Most monomer models showed high confidence, with pLDDT scores predominantly above 70. Each complex was then annotated with interface residues defined by geometric contact and hotspot residues identified through computational alanine scanning with the FoldX force field, which flags positions whose mutation would destabilize the interface by at least 0.5 kilocalories per mole.

To demonstrate the database’s power, the researchers dove deep into the chicken interactome. Its network topology follows classic small-world, scale-free architecture, with an average of 24 interactions per protein, a mean shortest path length of 3.3, and a clustering coefficient of 0.63—far higher than randomized networks of the same size. Fifty predicted chicken interactions overlapped significantly with experimentally supported ones, and predicted pairs showed Gene Ontology semantic similarity comparable to known interactions and far above random pairs. A flagship case study examined the actin-capping protein assembly: the predicted CAPZA1–CAPZB heterodimer aligned with the experimentally solved crystal structure with a root-mean-square deviation of just 1.6 angstroms, and the predicted binding interfaces for both subunits with ACTIN matched mutagenesis-verified actin-binding motifs at the C-termini of the CapZ subunits.

Perhaps the most compelling illustration of AASIA’s practical value involves a predicted interaction between XRCC5, a DNA repair factor, and H2AC39, a histone variant, in chicken. The interaction was independently supported by all three prediction methods and scored 0.94 overall. Remarkably, the predicted interface overlaps a known quantitative trait locus associated with eggshell strength. A residue at that locus, arginine 10 of XRCC5, sits directly at the interaction interface and was flagged as a putative hotspot: FoldX predicted that mutating it to alanine would destabilize the complex by 1.09 kilocalories per mole, and the database’s variant effect module predicted a corresponding drop in the interaction score. This offers a testable, structure-based hypothesis that genetic variation at this position could alter chromatin dynamics during eggshell mineralization—precisely the kind of genotype-to-phenotype bridge molecular breeders have long sought.

The platform itself is designed for broad usability, built on a Vue.js frontend, Express.js backend, and MySQL metadata management, with Cytoscape.js for network visualization and Mol* and NGL Viewer for interactive 3D structure exploration. Users can search by UniProt, Ensembl, or Entrez identifiers, gene names, or keywords; browse protein and interaction pages with functional annotations from Gene Ontology and KEGG; and, crucially, submit their own protein sequences to the Predict module for de novo interaction inference or quantitative assessment of mutation effects on interaction scores. This last capability turns AASIA from a static repository into an active hypothesis-generating engine.

The authors acknowledge limitations: public databases mix direct binary interactions with association-level evidence, and ESMFold does not explicitly model intrinsically disordered regions that mediate many interactions. Still, performance on a stricter binary-only test subset was even stronger, with an AUROC of 0.928. Looking forward, the team plans to expand AASIA with host–pathogen interaction data—potentially illuminating conserved interfaces exploited by avian influenza and African swine fever viruses—and to integrate transcriptomic, epigenomic, and genome-wide variant data. By transforming raw omics information into spatially resolved, residue-level insight, AASIA positions itself as a foundational tool for the next generation of livestock genomics and biotech breeding.

Subject of Research: A structure-informed protein–protein interactome database for five agricultural animal species

Article Title: AASIA: A comprehensive protein structural interactome database for agricultural animals

Article References: Li, J., Yuan, M., Jiang, L., Li, D., Yin, Z., Zhu, F., Shi, W., Hou, Z., & Zhang, Z. (2026). AASIA: A comprehensive protein structural interactome database for agricultural animals. Journal of Advanced Research. https://doi.org/10.1016/j.jare.2026.09.015

Image Credits: AI Generated

DOI: 10.1016/j.jare.2026.09.015

Keywords: AASIA, protein-protein interactions, agricultural animals, structural biology, ESMFold, deep learning, interolog mapping, molecular breeding, livestock genomics, variant effect prediction, disease resistance, food security