← back

Zero-shot design of drug-binding proteins via neural iterative selection−expansion

Benjamin Fry (et. al) · Paper · June 27, 2026

Main Takeaways

  1. A self-consistent folding → inverse-folding loop can generate strong binders
    1. Ligand self-consistency is the stronger screen than backbone self-consistency
  2. Pretty easy to add ligand conditioning to inverse-folding models, not sure how I never thought about this lol
  3. Neural proofreading is really good: LASErMPNN can suggest mutations, to improve binding, which yielded some super good results.

Model Approach

Neural Iterative Selection-Expansion (NISE): coordinate ascent over inverse-folding (LASErMPNN) and folding (Boltz-2, RFAA) models:

\[P(\text{sequence} \mid \text{structure, ligand}) \to P(\text{structure, ligand} \mid \text{sequence}) \\ \to P(\text{sequence} \mid \text{structure, ligand}) \to \ldots\]

Self-consistency hypothesis: complexes that are robust to this iteration are structurally sound → good binders. I’ve always felt this idea is pretty natural and nice.

LASErMPNN is inverse-folding with a few main changes:

  1. Condition on protein backbone + ligand atoms
    1. pretrained ligand encoder for ligand atoms:
    The ligand node embeddings (e1, e2, etc) are taken from a pretrained ligand encoder (box labelled pretrained ligand encoder) and locked during training. The ligand encoder is trained to predict the atomic properties of ligands (such as the partial atomic charge) conditioned on atomic coordinates and elements.
  2. Predicts amino-acid identity + sidechain dihedrals jointly
    1. latter matters for bound complexes because overpacking is more relevant (as opposed to protein-only) and ligands mainly interact with sidechains anyway
    Autoregressive sampling (AS) is performed in turn for each L residue of the protein. Each new residue prediction is conditioned on previously decoded residue identities and their associated dihedral angles

Disorganized Notes

Iteratively selecting Rosetta-minimized structures by ligand energy and expanding them using sequence design with LASErMPNN did not reduce the NLL of designed sequences nor increase the pLDDT of ligands in predicted co-structures. This implies that these structures are not sufficiently altered to more strongly encode optimal sequences. These results show that both neural networks in NISE are important for optimization.

This is a good point; use cofolding as basin-escaping behavior, pretty much.

NISE is agnostic to the specific networks used; as these models improve, so will NISE. Indeed, a RFAA-based assessment would have led us to discard our apixaban binders (Supplementary Fig. 36); the use of Boltz-2 was important to select binders for experimental testing.