CODesign: Consistency from Data to Trajectory in All-Atom Protein Binder Co-Design
Organizations: Valhalla Technology · The Chinese University of Hong Kong · University of Washington · Westlake University · Zhejiang University
Abstract
The central challenge in de novo protein design is generating plausible, mutually compatible structures and sequences, such that each designed sequence folds into its intended structure and the structure accommodates that sequence. Compared to typical two-stage design methods, which decouple the modeling of the interdependent modalities, co-design models improve the cross-modal consistency by jointly generating sequences and structures. However, naively generating sequences and structures simultaneously does not ensure their consistency. To address this challenge, we propose CODesign framework. We improve data consistency by generating approximately 105,000 consistency-distilled dimers. We further promote consistency through a multimodal joint flow model that captures the joint distribution of sequences, backbone structures, and local atomic configurations, together with a consistency-aware joint resampling strategy that iteratively refines sequences and side chains. Experiments show that CODesign achieves state-of-the-art performance with the highest in silico success rates on both protein- and ligand-target binder design. Ablation studies also demonstrate our distilled dataset increases performance by 70.9%, which can be further improved by our proposed resampling mechanism with negligible additional computational cost. Code, model weights and the new dataset will be completely open-source.
Figures & tables
| Model | # Unique Successes | scRMSD / Å (%) | Diversity | Novelty | ||||
| Self | MPNN@1 | MPNN@8 | Self | MPNN@1 | MPNN@8 | |||
| two-stage | ||||||||
| RFDiffusion3 | – | 6.46 | 10.96 | – | 7.24 / 36.12 | 4.17 / 53.02 | 25.66 | 0.91 |
| PXDesign | – | 4.59 | 7.01 | – | 5.21 / 51.15 | 2.29 / 74.54 | 15.39 | 0.92 |
| Protpardelle-1c | – | 0.28 | 0.92 | – | 16.44 / 10.26 | 12.61 / 20.40 | 75.15 | 0.88 |
| co-design | ||||||||
| Model | # Unique Successes | Diversity | Novelty | Time (s) | ||||
| SAM | OQO | FAD | IAI | Mean | ||||
| two-stage | ||||||||
| RFdiffusionAA | 2.0 | 3.4 | 1.6 | 5.0 | 3.0 | 135.65 | 0.74 | 73.78 |
| RFDiffusion3 | 9.4 | 14.8 | 20.8 | 11.0 | 14.0 | 52.50 | 0.83 | 14.96 |
| co-design | ||||||||
| Complexa | 3.2 | 3.0 | 5.8 | 8.4 | 5.1 | 147.70 | 0.76 | 10.27 |
Appendix figures & tables18 assets
Supplementary material from the paper’s appendix.
Appendix
| Source | Reference dimers | Codimer pairs |
| AFDB high-confidence heterodimers | 25,837 | 23,823 |
| HumanPPI | 9,757 | 5,177 |
| PINDER | 156,242 | 75,528 |
| Total | 191,836 | 104,528 |
| Dataset | Protein I | Protein II | Ligand adaptation |
| Teddymer | 0.80 | 0.70 | – |
| PDB binder–target examples | 0.19 | 0.10 | – |
| PDB monomers | 0.01 | 0.01 | – |
| Codimer | – | 0.19 | – |
| CDDB monomers | – | – | 0.60 |
| PLINDER | – | – | 0.40 |
| Hyperparameter | Protein-conditioned | Ligand-conditioned |
| Residue trunk | ||
| Representation dimension | 768 | 768 |
| Conditioning dimension | 256 | 256 |
| Pair dimension | 256 | 256 |
| Transformer layers | 14 | 14 |
| Attention heads | 12 | 12 |
| Setting | Protein I | Protein II | Ligand adaptation |
| Number of GPUs | 24 | 64 | 32 |
| Batch size per GPU | 5 | 5 | 3 |
| Global batch size | 120 | 320 | 96 |
| Updates at selected checkpoint | 60,000 | 10,000 | 10,000 |
| Trainable parameters | Full model | Full model | LoRA adapters |
| Optimizer | Adam | Adam | Adam |
| Setting | Backbone | Sequence | Local atoms |
| Active updates | 500 | 500 | 1,000 |
| Time schedule | Exponential | Quadratic | Exponential |
| Noise-variance factor | 0.3 | – | 0.6 |
| Categorical temperature | – | 0.3 | – |
| Discrete stochasticity | – | 5 | – |
| ODE switching time | 0.98 | – | 0.98 |
| Model | Sequence | Unique successes | scRMSD (Å) | Å (%) |
| RFDiffusion3 | MPNN@1 | |||
| RFDiffusion3 | MPNN@8 | |||
| PXDesign | MPNN@1 | |||
| PXDesign | MPNN@8 | |||
| Protpardelle-1c | MPNN@1 | |||
| Protpardelle-1c | MPNN@8 |
| Target | Complexa | CODesign |
| IFNAR2 | ||
| BHRF1 | ||
| BBF14 | ||
| DerF21 | ||
| TrkA | ||
| PD1 |
| Target | RFDiffusion3 | PXDesign | Protpardelle-1c | Complexa | CODesign |
| IFNAR2 | |||||
| BHRF1 | |||||
| BBF14 | |||||
| DerF21 | |||||
| TrkA | |||||
| PD1 |
| Target | RFDiffusion3 | PXDesign | Protpardelle-1c | Complexa | CODesign |
| IFNAR2 | |||||
| BHRF1 | |||||
| BBF14 | |||||
| DerF21 | |||||
| TrkA | |||||
| PD1 |
| Target | Complexa | CODesign |
| IFNAR2 | ||
| BHRF1 | ||
| BBF14 | ||
| DerF21 | ||
| TrkA | ||
| PD1 |
| Target | RFDiffusion3 | PXDesign | Protpardelle-1c | Complexa | CODesign |
| IFNAR2 | |||||
| BHRF1 | |||||
| BBF14 | |||||
| DerF21 | |||||
| TrkA | |||||
| PD1 |
| Target | RFDiffusion3 | PXDesign | Protpardelle-1c | Complexa | CODesign |
| IFNAR2 | |||||
| BHRF1 | |||||
| BBF14 | |||||
| DerF21 | |||||
| TrkA | |||||
| PD1 |
| Target | Complexa | CODesign |
| IFNAR2 | ||
| BHRF1 | ||
| BBF14 | ||
| DerF21 | ||
| TrkA | ||
| PD1 |
| Target | RFDiffusion3 | PXDesign | Protpardelle-1c | Complexa | CODesign |
| IFNAR2 | |||||
| BHRF1 | |||||
| BBF14 | |||||
| DerF21 | |||||
| TrkA | |||||
| PD1 |
| Target | RFDiffusion3 | PXDesign | Protpardelle-1c | Complexa | CODesign |
| IFNAR2 | |||||
| BHRF1 | |||||
| BBF14 | |||||
| DerF21 | |||||
| TrkA | |||||
| PD1 |
| Model | SAM | OQO | FAD | IAI | Mean |
| RFdiffusionAA | |||||
| RFDiffusion3 | |||||
| Complexa | |||||
| Pallatom-Ligand | |||||
| DISCO | |||||
| CODesign |
| Variant | Codimer | CAR | Unique successes | Time (s) |
| Without Codimer , without CAR | No | No | – | |
| Without CAR | Yes | No | 27.96 | |
| CODesign | Yes | Yes | 28.41 |
| Target | Without Codimer without CAR | With Codimer without CAR | CODesign |
| IFNAR2 | |||
| BHRF1 | |||
| BBF14 | |||
| DerF21 | |||
| TrkA | |||
| PD1 |