Folding, Reasoning, and Scaling with Open-source Drug Discovery Engine
Summary
Introduces OpenDDE, an open-source all-atom biomolecular foundation model that uses co-folding as a shared structural reasoning layer for drug discovery, achieving high accuracy and identifying scaling laws for further improvement.
View Cached Full Text
Cached at: 07/07/26, 04:35 AM
# Folding, Reasoning, and Scaling with Open-source Drug Discovery Engine
Source: [https://arxiv.org/html/2607.03787](https://arxiv.org/html/2607.03787)
AbstractAccurately modeling biomolecular interactions is a central bottleneck in biology and therapeutic discovery\. Here, we introduce Open Drug Discovery Engine \(OpenDDE\), an open\-source, all\-atom biomolecular foundation model that uses co\-folding as the entry point to a scalable AI\-driven drug discovery engine\. Rather than treating structure prediction as an isolated endpoint,OpenDDEis designed as a shared structural reasoning layer for modeling sequence–structure–function relationships across biomolecular complexes, enabling complex structure prediction today while providing a foundation forde novodesign, affinity estimation, structure\-conditioned optimization, and more\.OpenDDEintegrates advances in all\-atom architecture, atomic latent reasoning, inference optimization, and large\-scale data processing to achieve IsoDDE\-level co\-folding accuracy within a reproducible and openly accessible framework\. We also identify two scaling\-law directions for co\-folding models, revealing practical routes for continued improvement through data, model, inference, and training scaling\. By releasing training code, inference pipelines, checkpoints, and benchmarks,OpenDDEaims to democratize access to frontier biomolecular intelligence, accelerate global collaboration, and lay an open foundation for next\-generation drug discovery systems that can move from predicting molecular structures toward designing, scoring, and optimizing therapeutic candidates for human health\.GitHub:[https://github\.com/aurekaresearch/OpenDDE](https://github.com/aurekaresearch/OpenDDE) HuggingFace:[https://huggingface\.co/aurekaresearch/OpenDDE](https://huggingface.co/aurekaresearch/OpenDDE)
![[Uncaptioned image]](https://arxiv.org/html/2607.03787v1/x1.png)## Overview
Key Contributions1\.Atomic latent reasoning over biomolecular tokens\.We incorporate latent reasoning into biomolecular modelling, enabling the model to refine representations of local geometry, chemical context, and cross\-molecular interfaces before all\-atom structure generation\.2\.A folding\-centered foundation for an extensible drug\-discovery engine\.OpenDDEcurrently focuses on structure prediction, but its unified architecture is designed to support futurede novomolecular design, affinity prediction, and other structure\-conditioned modules\.3\.Scaling laws and data distillation\.We study scaling directions along model\-parameter and data axes, and revise distillation strategies for monomeric and multimeric structures to improve robustness, generalization, and training efficiency\.
Characterising biomolecular interactions among proteins, nucleic acids, small\-molecule ligands and other cellular components is fundamental to understanding biological mechanisms and to modulating disease\-relevant processes through therapeutic intervention\[[1](https://arxiv.org/html/2607.03787#bib.bib1),[2](https://arxiv.org/html/2607.03787#bib.bib2)\]\. In silico models that can predict these interactions with experimental\-level accuracy would therefore provide a structural basis for scalable and reliable drug discovery\.
AlphaFold2\[[3](https://arxiv.org/html/2607.03787#bib.bib3)\]transformed protein structure prediction by achieving near\-experimental accuracy for many protein monomers\. AlphaFold3\[[4](https://arxiv.org/html/2607.03787#bib.bib4)\]extended this paradigm to joint, atomic\-level modelling of complexes containing proteins, nucleic acids, small molecules, ions and modified residues\. More broadly, AlphaFold\-style and structure\-conditioned models have become reusable engines for epitope prediction, binder design\[[5](https://arxiv.org/html/2607.03787#bib.bib5),[6](https://arxiv.org/html/2607.03787#bib.bib6),[7](https://arxiv.org/html/2607.03787#bib.bib7)\], virtual screening\[[8](https://arxiv.org/html/2607.03787#bib.bib8),[9](https://arxiv.org/html/2607.03787#bib.bib9)\], conformational sampling\[[10](https://arxiv.org/html/2607.03787#bib.bib10),[11](https://arxiv.org/html/2607.03787#bib.bib11)\]and the acceleration or interpretation of experimental workflows\[[12](https://arxiv.org/html/2607.03787#bib.bib12)\]\. The recent introduction of IsoDDE further shows that co\-folding can serve as the core of an AI drug discovery engine, with reported gains in protein\-ligand generalisation, antibody\-antigen interface prediction, pocket identification and binding affinity estimation\. However, the most capable system remains closed, limiting reproducibility, independent validation and community\-driven extension\. At the same time, incremental improvements in co\-folding accuracy have begun to slow\[[13](https://arxiv.org/html/2607.03787#bib.bib13)\], raising a central question for the field:whether current tokenization, data curation, model architectures, and inference procedures are approaching practical bottlenecks, and where the next scaling directions should come from?
In this manuscript, we present OpenDDE, to our knowledge the first Apache\-licensed all\-atom generative foundation model reported to reach IsoDDE\-level co\-folding accuracy\. Building on recent open co\-folding systems such as Protenix\-v1 and OpenFold3,OpenDDEintroduces three advances that extend the co\-folding paradigm across both pre\-training and post\-training:
In the following sections, we present the model architecture, training, inference, and scaling laws behindOpenDDE, followed by benchmarks against open and closed folding systems across diverse biomolecular complex prediction tasks\.OpenDDEperforms on par with state\-of\-the\-art closed models and substantially improves over existing models on antibody\-antigen generalization benchmarks\.
As an evolving project,OpenDDEand its repository will be regularly updated with contributions from our team and the broader community\. We release checkpoints, code, inference pipelines, and benchmarks under the Apache\-2\.0 license to enable the community to explore, refine, and deploy biological intelligence at scale\. Beyond strong co\-folding performance,OpenDDEprovides an open platform for studying how folding, reasoning, and scaling can advance biomolecular intelligence\.
## 1 Performance
### 1\.1 Structure Prediction Result
Figure 2:Antibody–antigen structure prediction performance across PXMeter\-AB, FoldBench\-AB, and 2026ARK\-AB\.Results are reported as DockQ success rates, decomposed into acceptable\-, medium\-, and high\-quality predictions\.Left:Ranking\-based selection\.Right:Oracle\-based selection\.BenchmarkingOpenDDEon Antibody\-Antigen Dataset \(Figure[2](https://arxiv.org/html/2607.03787#S1.F2)\)\.Antibody\-antigen complexes provide a stringent test for structure\-based discovery, as their therapeutic relevance is matched by the difficulty of accurately modelling flexible and chemically diverse binding interfaces\. Figure[2](https://arxiv.org/html/2607.03787#S1.F2)comparesOpenDDEwith AlphaFold3\[[4](https://arxiv.org/html/2607.03787#bib.bib4)\], Chai\-1\[[14](https://arxiv.org/html/2607.03787#bib.bib14)\], Boltz\-1\[[15](https://arxiv.org/html/2607.03787#bib.bib15)\], OpenFold3\[[16](https://arxiv.org/html/2607.03787#bib.bib16)\], Protenix\-v1\[[17](https://arxiv.org/html/2607.03787#bib.bib17)\], and ESMFold2\[[18](https://arxiv.org/html/2607.03787#bib.bib18)\]across three antibody\-antigen benchmark collections\. Success rates are decomposed by interface quality using DockQ thresholds, separating acceptable, medium\-quality, and high\-quality predictions\.OpenDDEachieves the highest overall success rate on all evaluated benchmarks, reaching 51\.0% on PXMeter\-AB, 70\.0% on FoldBench\-AB, and 66\.4% on the 2026ARK\-AB benchmark\.
These results substantially improve over the strongest baselines, including ESMFold2 on PXMeter\-AB and FoldBench\-AB and ESMFold2/Protenix\-v1 on 2026ARK\-AB\. Importantly, the gains are not limited to low\-threshold recovery:OpenDDEconsistently increases the fraction of medium\- and high\-DockQ interfaces, indicating improved modeling of antibody\-antigen geometry rather than merely producing marginally acceptable complexes\. Together, these results show thatOpenDDEprovides a strong open\-source foundation for antibody\-antigen co\-folding and brings a broader set of therapeutic targets into a practically useful accuracy regime\.
Structure Prediction of New Biological Systems \(Figure[2](https://arxiv.org/html/2607.03787#S1.F2)\)\.We examineOpenDDEon newly released antibody\-antigen structures by curating the 2026ARK\-AB benchmark, following the data protocol of PXMeter\[[19](https://arxiv.org/html/2607.03787#bib.bib19)\]\. This benchmark contains 164 PDB complexes, corresponding to 159 unique antibody\-antigen interface clusters, and is curated to assess prediction performance on recent, low\-homology systems that are not available to earlier released models\. Antibody\-antigen interfaces are defined by residue\-level spatial proximity and annotated using SAbDab\[[20](https://arxiv.org/html/2607.03787#bib.bib20)\]\. We retain high\-quality biological assemblies using standard filters for structure quality, assembly size, chain completeness, and interface geometry\. Protein entities are clustered with MMseqs2 easy\-cluster using a 40% minimum sequence identity threshold and an 80% alignment coverage requirement, and each interface cluster is defined as the unordered pair of the two interacting entity clusters\.
On this newest benchmark,OpenDDEachieves the highest overall success rate, correctly predicting 66\.4% of interfaces at DockQ\>\>0\.23\. It outperforms ESMFold2 \(51\.0%\), Protenix\-v1 \(50\.0%\), OpenFold3 \(21\.9%\), Boltz\-1 \(25\.0%\), and Chai\-1 \(16\.7%\)\.OpenDDEalso shows a clear advantage in the high\-fidelity regime, with a large fraction of complexes reaching DockQ\>\>0\.8\. These results demonstrateOpenDDE’s strong ability to predict newly released antibody\-antigen interfaces and suggest improved generalization to emerging biological systems\.
Figure 3:Cumulative distribution of DockQ across PXMeter\-AB, FoldBench\-AB, and 2026ARK\-AB\.Each curve shows the fraction of targets whose top\-ranked prediction reaches DockQ≥\\geqthe x\-axis cutoff, for Chai\-1, OpenFold3, Protenix\-v1, and OpenDDE; curves closer to the upper\-right are better\.New Model Capacity \(Figure[2](https://arxiv.org/html/2607.03787#S1.F2)\)\.We evaluate model capacity under an oracle selection setting, where the best prediction is selected from each model’s generated candidates\. This setting approximates the maximal achievable performance of a model when ranking or confidence estimation is not the limiting factor, and therefore reflects the upper bound of its sampling capacity\. As shown in Figure[2](https://arxiv.org/html/2607.03787#S1.F2),OpenDDEachieves the highest oracle\-selected success rate across all antibody\-antigen benchmarks, reaching 65\.9% on PXMeter\-AB, 81\.9% on FoldBench\-AB, and 80\.1% on 2026ARK\-AB\. These results substantially exceed the corresponding non\-oracle benchmark performance, indicating thatOpenDDEfrequently generates high\-quality antibody\-antigen conformations even when the top\-ranked prediction is not always optimal\. The large oracle gains suggest thatOpenDDEhas strong latent modeling capacity for antibody\-antigen interfaces, and that further improvements in confidence estimation or candidate selection could translate this capacity into higher realized prediction accuracy\.
DockQ Success\-rate Curve \(Figure[3](https://arxiv.org/html/2607.03787#S1.F3)\)\.The DockQ success\-rate curves provide a threshold\-sweep view of interface prediction quality\. For each DockQ cutoff, the curve reports the fraction of targets whose predicted interface exceeds that quality threshold, so curves that remain higher across the x\-axis indicate better performance over the full range from acceptable to high\-quality predictions\. Across PXMeter\-AB, FoldBench\-AB, and 2026ARK\-AB,OpenDDEconsistently maintains the highest curve, showing that its advantage is not restricted to a single DockQ threshold\. The gap is especially clear in the medium\- and high\-quality regimes, whereOpenDDEretains a larger fraction of successful predictions as the cutoff becomes stricter\. This indicates thatOpenDDEimproves both coarse interface recovery and high\-fidelity antibody\-antigen docking accuracy, producing more predictions with accurate binding geometry rather than only increasing low\-threshold success\.
### 1\.2 Scaling Laws in Biomolecular Foundation Models
Figure 4:Scaling laws in biomolecular foundation models\.Left:Model performance as a function of estimated training tokens, computed as training steps×\\timesnumber of GPUs \(batch size\)×\\timescrop tokens\.Right:Model performance as a function of estimated training cost, computed as estimated training tokens×\\timestrainable model parameters\.We observe a clear scaling trend in biomolecular foundation models\. To examine this trend, we compareOpenDDEwith AlphaFold3\[[4](https://arxiv.org/html/2607.03787#bib.bib4)\], Protenix\-v1\[[17](https://arxiv.org/html/2607.03787#bib.bib17)\], Protenix\-v2\[[21](https://arxiv.org/html/2607.03787#bib.bib21)\], SeedFold\[[22](https://arxiv.org/html/2607.03787#bib.bib22)\], and ESMFold2\[[18](https://arxiv.org/html/2607.03787#bib.bib18)\]on antibody–antigen structure prediction \(Figure[4](https://arxiv.org/html/2607.03787#S1.F4)\)\. We use two estimates of training scale\. The first is the estimated number of training tokens, computed as training steps×\\timesnumber of GPUs \(batch size\)×\\timescrop tokens\.111For OpenDDE, this is estimated as400×\(10k×384\+60k×384\+15k×544\+15k×544\+10k×768\)≈2\.04×1010400\\times\(10\\mathrm\{k\}\\times 384\+60\\mathrm\{k\}\\times 384\+15\\mathrm\{k\}\\times 544\+15\\mathrm\{k\}\\times 544\+10\\mathrm\{k\}\\times 768\)\\approx 2\.04\\times 10^\{10\}\.The second is the estimated training cost, computed as estimated training tokens×\\timestrainable model parameters\.222For OpenDDE, this is estimated as2\.04×1010×655M≈1\.33×10192\.04\\times 10^\{10\}\\times 655\\mathrm\{M\}\\approx 1\.33\\times 10^\{19\}\.
Evidence of Scaling Laws \(Figure[4](https://arxiv.org/html/2607.03787#S1.F4)\)\.The scaling\-law curves show a clear positive relationship between training scale and Ab\-Ag performance\. As the estimated training compute increases, model performance improves from the lower\-scale models such as Protenix\-v1 and AlphaFold3 toward higher\-performing models such as Protenix\-v2 and OpenDDE\. The trend is smooth and sublinear, resembling scaling behavior observed in large language models: early increases in scale produce relatively large gains, while later gains become more gradual as the curve approaches the IsoDDE reference level\.
OpenDDEsits at the high\-scale, high\-performance end of the curve, achieving the best observed Ab\-Ag performance among the compared models and narrowing the gap to IsoDDE\. The remaining gap suggests that further scaling may continue to improve performance, but with diminishing returns\. Deviations among individual models also indicate that scale alone does not fully determine performance; architecture, data quality, and training strategy likely contribute to where each model falls relative to the overall scaling trend\.
Entering the Era of Scaling \(Figure[4](https://arxiv.org/html/2607.03787#S1.F4)\)\.Across these models and performance, the results improve with both data scale and effective compute scale\. This observation suggests that biomolecular structure modeling is entering the same scaling regime that has transformed \(large\) language models:larger effective training corpora, larger\-scale models, and greater compute can systematically translate into stronger biological reasoning and structure prediction capability\.
### 1\.3 Test\-Time Scaling
Figure 5:Test\-time scaling on FoldBench\-AB and 2026ARK\-AB\.Top:Ranking\-based selection of DockQ\>\>0\.8\.Bottom: Oracle\-based selection\.We evaluate whether additional inference\-time sampling improves antibody\-antigen prediction performance \(Figure[5](https://arxiv.org/html/2607.03787#S1.F5)\)\. For each target, we generate multiple stochastic samples with different random seeds and select the final prediction using the model ranking score\. We report two settings: ranking\-based high\-quality success rate with DockQ\>\>0\.8, and oracle success rate, where the best candidate among all generated samples is selected using the ground\-truth DockQ\. The oracle setting measures the maximal performance available in the sampled candidate set and separates sampling capacity from ranking quality\.
As shown in Figure[5](https://arxiv.org/html/2607.03787#S1.F5),OpenDDEexhibits consistent test\-time scaling on both FoldBench\-AB and 2026ARK\-AB\. On FoldBench\-AB, increasing the number of seeds improves the ranked high\-DockQ rate from roughly 28% with one seed to about 34% at large sampling budgets\. Similar behavior is observed on 2026ARK\-AB, whereOpenDDEimproves from approximately 35% to 38% ranked high\-DockQ successThese gains indicate that additional samples expose better antibody\-antigen conformations rather than simply producing redundant predictions\.
The oracle curves show a larger scaling effect\. On FoldBench\-AB,OpenDDE’s oracle success rises from about 66% with one seed to above 90% with hundreds of seeds\. On 2026ARK\-AB, the oracle success similarly increases from about 67% to nearly 90%\. This gap between ranked and oracle performance suggests thatOpenDDEoften generates high\-quality antibody\-antigen structures that are not always selected as the top\-ranked prediction\. Therefore, test\-time scaling reveals strong latent sampling capacity, and further improvements in confidence calibration or candidate ranking could convert more of this oracle performance into realized prediction accuracy\.
Across both metrics,OpenDDEscales more favorably with inference\-time compute than the compared open models, including Boltz\-1, OpenFold3, and Protenix\-v1\. These results show thatOpenDDEbenefits predictably from larger sampling budgets and that its diffusion\-based generation process provides a practical route for improving antibody\-antigen prediction accuracy at test time\.
### 1\.4 Case Study
Figure 6:Case studies of antibody\-antigen structure predictions, from Protenix\-v1 and OpenDDE\.Left:ALB8 VHH binding to human serum albumin \(PDB: 8Z8V, 2024\-09\)\.Middle:Antibody binding to human norepinephrine transporter NET \(PDB: 8ZP2, 2024\-09\)\.Right:Anti\-HER2 Fab H2Mab\-119 binding to HER2 domain I \(PDB: 8JYR, 2024\-03\)\. OpenDDE predictions are shown in purple; Protenix\-v1 predictions are shown in teal, each superimposed on the ground\-truth structure are shown in gray\.We show predicted structures of new biological systems in Figure[6](https://arxiv.org/html/2607.03787#S1.F6)\. OpenDDE predictions, colored in purple, closely recover the experimentally resolved binding geometry shown in gray across all three systems\. In 8Z8V, OpenDDE places the ALB8 VHH on the correct surface of human serum albumin, preserving both the global orientation and the local interface geometry\. In 8ZP2, OpenDDE recovers the nanobody position on the NET transporter complex, despite the added difficulty of a membrane protein target and cryo\-EM\-derived structure\. In 8JYR, OpenDDE accurately aligns the antibody variable region against HER2 domain I, capturing the correct epitope\-facing pose\.
These case studies suggest that OpenDDE improves structure prediction at the level most relevant for biological discovery: the geometry of molecular recognition\. While the examples are qualitative and should be complemented by cutoff\-controlled quantitative benchmarks, their recent release dates make them useful stress tests for assessing OpenDDE on newly solved biological systems rather than only on saturated historical cases\.
## 2 Training Details
### 2\.1 Model Parameters
OpenDDEhas 655M training parameters using NVIDIA computing kernels\. We report model, training, and sampling hyperparameters in Table[1](https://arxiv.org/html/2607.03787#S2.T1)\.
Table 1:Model, training, and sampling hyperparameters of OpenDDE\.\(a\) Model hyperparameters
\(b\) Training and sampling hyperparameters
Table 2:Training stages for OpenDDE\. Stage IV\(a\) inherits the checkpoint from Stage III, but freezes all major modules and trains only the confidence head\. Stage IV\(b\) also inherits the checkpoint from Stage III and increases the crop size from 544 to 768\. Note: The Teddymer dataset was introduced in\[[23](https://arxiv.org/html/2607.03787#bib.bib23)\], and we reproduced the dataset for model training\.ParameterStage warmupStage IStage IIStage IIIStage IV\(a\)Stage IV\(b\)Steps10k60k15k15k10k10kGPUs50\*8 80G50\*8 80G50\*8 80G50\*8 80G50\*8 80G8\*8 141GCrop tokens384384544544544768Max atoms3,5003,5006,0006,0006,0007,500Max MSA depth2,0482,0482,0482,0482,0482,048MSA dropout probability0\.20\.20\.20\.20\.20\.2Diffusion batch size646432323232Smooth LDDT loss✓\\checkmark✓\\checkmark✓\\checkmark✓\\checkmark×\\times✓\\checkmarkShape Complementary loss✓\\checkmark✓\\checkmark✓\\checkmark✓\\checkmark×\\times✓\\checkmarkTrain structure prediction✓\\checkmark✓\\checkmark✓\\checkmark✓\\checkmark×\\times✓\\checkmarkTrainde novodesign×\\times×\\times×\\times✓\\checkmark×\\times✓\\checkmarkDe novodesign probability0000\.100\.2Train confidence head✓\\checkmark✓\\checkmark✓\\checkmark✓\\checkmark✓\\checkmark✓\\checkmarkWeighted PDB\[[17](https://arxiv.org/html/2607.03787#bib.bib17)\]70%50%35%40%50%40%AFDB\-multimer\[[24](https://arxiv.org/html/2607.03787#bib.bib24)\]4%5%11%9%5%9%Teddymer Reproduction4%5%11%9%5%9%MGnify long monomer\[[16](https://arxiv.org/html/2607.03787#bib.bib16)\]10%28%30%25%20%20%MGnify short monomer\[[16](https://arxiv.org/html/2607.03787#bib.bib16)\]4%5%4%4%4%4%Swiss\-Prot\[[3](https://arxiv.org/html/2607.03787#bib.bib3)\]6%5%4%4%5%4%Disordered\[[16](https://arxiv.org/html/2607.03787#bib.bib16)\]2%2%2%1%1%1%SAbDab\[[20](https://arxiv.org/html/2607.03787#bib.bib20)\]0%0%3%8%10%13%
Compared with AlphaFold3, the main architectural difference is thatOpenDDEscales the Pairformer hidden dimension from 128 to 384, which increases the Pairformer parameter count by approximately3×3\\timesand its computational cost by approximately9×9\\times\. Beyond this scaling,OpenDDEintroduces atomic\-level reasoning blocks \(Section[A\.1](https://arxiv.org/html/2607.03787#A1.SS1)\), shape\-complementarity losses for interface fitting \(Section[3\.2](https://arxiv.org/html/2607.03787#S3.SS2)\), and a unified training formulation for structure prediction andde novodesign \(Section[3\.3](https://arxiv.org/html/2607.03787#S3.SS3)\)\.
### 2\.2 Training Data
Our training data follow the data construction protocols used in AlphaFold3\[[4](https://arxiv.org/html/2607.03787#bib.bib4)\], Protenix\[[17](https://arxiv.org/html/2607.03787#bib.bib17)\], and ESMFold2\[[18](https://arxiv.org/html/2607.03787#bib.bib18)\]\. To avoid test\-set leakage, all ground\-truth structures and template databases are filtered using a cutoff date of 2021\-09, as we directly adopted from\[[4](https://arxiv.org/html/2607.03787#bib.bib4),[17](https://arxiv.org/html/2607.03787#bib.bib17),[16](https://arxiv.org/html/2607.03787#bib.bib16),[18](https://arxiv.org/html/2607.03787#bib.bib18)\]\.
We collect training examples from multiple structure and sequence resources, including weighted PDB data from Protenix\[[17](https://arxiv.org/html/2607.03787#bib.bib17)\], AFDB multimer data from Han et al\.\[[24](https://arxiv.org/html/2607.03787#bib.bib24)\], reproduced Teddymer dataset according to DiDi et al\.\[[23](https://arxiv.org/html/2607.03787#bib.bib23)\], MGnify long and short monomers reported in AlphaFold3\[[4](https://arxiv.org/html/2607.03787#bib.bib4)\]and released by OpenFold3\-Preview\[[16](https://arxiv.org/html/2607.03787#bib.bib16)\], Swiss\-Prot entries from AlphaFold DB\[[3](https://arxiv.org/html/2607.03787#bib.bib3)\], disordered motifs reported in AlphaFold3\[[4](https://arxiv.org/html/2607.03787#bib.bib4)\]and released by OpenFold3\-Preview\[[16](https://arxiv.org/html/2607.03787#bib.bib16)\], and SAbDab antibody structures reported in ESMFold2\[[18](https://arxiv.org/html/2607.03787#bib.bib18)\]and collected from SAbDab\[[20](https://arxiv.org/html/2607.03787#bib.bib20)\]\.
### 2\.3 Training Schedule
OpenDDEis trained in four major stages with one additional warmup stage\. For reproducibility, we report the training stages and the corresponding data distributions in Table[2](https://arxiv.org/html/2607.03787#S2.T2)\.
OpenDDEis initially trained on Ampere GPUs with a warmup learning rate of 0\.0012 for 10k steps, followed by a learning\-rate decay factor of 0\.95 every 60k steps\. After the first three stages, we conduct two Stage IV variants: Stage IV\(a\), in which all major modules are frozen and only the confidence module is trained, and Stage IV\(b\), in which the full model is continuously trained on Hopper GPUs for larger crop size\.
Training IntuitionThe intuition behind the stage\-wise data distribution is to trainOpenDDEfrom precision, to breadth, and then back to precision\. In the early stages, we emphasize large\-scale experimentally validated structures so that the model first learns accurate folding patterns and reliable local geometry from experimental data\. In the middle stages, we increase the use of distilled datasets to broaden coverage over sequence and structure space, allowing the model to learn more general fold\-level and geometric priors\. In the later stages, we shift the distribution back toward high\-quality experimental and task\-specific data to improve atomic accuracy, interface fitting, and confidence calibration\. At the same time, the crop size is gradually increased, enabling the model to first learn stable local geometry and then adapt to longer\-range interactions and larger biomolecular complexes\.
### 2\.4 Training Computation
We spent five weeks processing the metadata into trainable examples using10×810\\times 8NVIDIA 80GB GPUs\. This preprocessing stage included data downloading, feature processing, and MSA search\.
OpenDDEwas then trained on50×850\\times 8NVIDIA 80GB Ampere GPUs for six weeks, followed by one additional week on8×88\\times 8NVIDIA 141GB Hopper GPUs\. Additionally, we prepared5×85\\times 8NVIDIA 80GB Ampere GPUs to perform online validation\. In total,OpenDDErequired approximately 414K GPU\-hours for training, that is 54 years in a single computing unit\.
## 3 OpenDDE\-preview
\(a\)OverallOpenDDEarchitecture\.
\(b\)Atom37\-aware structural\-token reasoning\.
\(c\)Atomic shape complementary\.
Figure 7:Overview ofOpenDDEand its unified structure\-reasoning framework\. \(a\)OpenDDEbuilds residue\-token and pair representations, expands them into structural tokens, and generates all\-atom coordinates through diffusion\. \(b\) Atom37\-level geometry is organized into semantic structural tokens, enabling reasoning over backbone, side\-chain, base, and ligand roles before all\-atom refinement\. \(c\) Shape\-complementarity losses guide interface fitting by rewarding compatible contacts while penalizing gaps and steric clashes\.OpenDDEis a coarse\-to\-fine biomolecular structure generation model that unifies structure prediction and conditional design within one diffusion framework \(Details can be found in Appendix[A](https://arxiv.org/html/2607.03787#A1)\)\. Given sequence, atom, MSA, template, and constraint features,OpenDDEfirst builds residue\-token representations with a Pairformer\-style trunk\. These residue\-level states capture global relational context across chains, residues, ligands, and other molecular units\. The model then expands each residue\-token representation into chemically meaningful structural tokens, including protein backbone, side chain, nucleic\-acid backbone, base, and ligand/atom tokens\. A Structural Refiner performs fine\-grained relational reasoning over these structural tokens before a diffusion\-based coordinate generator denoises atom positions conditioned on the refined representations\. Confidence and distance heads provide residue\-level and structure\-level prediction outputs\.
### 3\.1 Reasoning in Biomolecular Foundation Models
We define structure reasoning as latent relational inference over molecular tokens before coordinate generation\. Rather than directly mapping sequence tokens to atom coordinates,OpenDDEfirst constructs single and pair representations that encode hypotheses about molecular geometry\. The residue\-token trunk captures global context, while the structural\-token branch decomposes each residue or molecular unit into finer chemical components and refines their interactions\.
This gives a coarse\-to\-fine reasoning path:
\(Sr,Zr\)→\(Ss,Zs\)→X,\(S^\{r\},Z^\{r\}\)\\rightarrow\(S^\{s\},Z^\{s\}\)\\rightarrow X,where\(Sr,Zr\)\(S^\{r\},Z^\{r\}\)are residue\-level single and pair states,\(Ss,Zs\)\(S^\{s\},Z^\{s\}\)are refined structural\-token states, andXXdenotes all\-atom coordinates\. Pair\-conditioned attention and triangular pair updates allow the model to reason not only about individual token\-token relations, but also about geometric consistency across indirect paths\. In this view, structural tokens form an intermediate layer between residue\-level context and atom\-level coordinates, enabling the model to represent backbone, side\-chain, base, and ligand interactions before denoising atoms\.
### 3\.2 Learning theLock\-and\-KeyBiomolecular Interactions
Accurate molecular structure generation requires more than placing atoms near their native positions\. For complexes, the interacting surfaces must face each other, maintain realistic gaps, avoid steric clashes, and form locally compatible packing\. To encourage this behavior,OpenDDEaugments diffusion training with geometry\-aware objectives that supervise both interface complementarity and all\-atom refinement\.
The shape\-complementarity objective compares predicted and native interface geometry using token\-level surface orientation, cross\-chain spacing, and clash\-aware contact quality\. This guides the model toward binding interfaces that are physically compatible rather than merely close under coordinate distance\. In parallel, local\-frame and torsion supervision improve all\-atom packing by constraining side chains, nucleic\-acid bases, and ligand/atom tokens in chemically meaningful coordinate systems\. Together, these losses teach the model the “lock\-and\-key” geometry of biomolecular interactions: global chains must fit, and local atoms must pack\.
### 3\.3 Structure Prediction andDe NovoDesign in One Framework
OpenDDEtreats structure prediction and de novo design as two instances of the same conditional denoising problem\. The input complex is represented as a partially observed molecular system, where some atoms can be fixed as structural context and the remaining atoms are generated as targets\. Standard structure prediction is the special case in which no atoms are fixed, while conditional design fixes one or more chains, motifs, pockets, or molecular contexts\.
This is implemented through known\-target atom masks:
maknown∈\{0,1\},matarget∈\{0,1\}\.m^\{\\mathrm\{known\}\}\_\{a\}\\in\\\{0,1\\\},\\qquad m^\{\\mathrm\{target\}\}\_\{a\}\\in\\\{0,1\\\}\.During conditional diffusion, known atoms remain fixed as structural context, while target atoms are noised and denoised:
pθ\(Xtarget∣ℐ,Xknown\)\.p\_\{\\theta\}\(X^\{\\mathrm\{target\}\}\\mid\\mathcal\{I\},X^\{\\mathrm\{known\}\}\)\.Whenmknown=0m^\{\\mathrm\{known\}\}=0, the same model performs full structure prediction\. Whenmknown≠0m^\{\\mathrm\{known\}\}\\neq 0, it performs conditional generation or design around the provided structure\.
This unified formulation allows prediction and design to share the same representation learning pipeline: residue\-token reasoning, structural\-token refinement, and all\-atom diffusion\. As a result,OpenDDEdoes not require a separate design\-specific architecture\. Instead, de novo design becomes masked molecular generation under structural conditions, while structure prediction remains the fully unconditioned case\.
## 4 Next Steps
We are actively collaborating with the NVIDIA Fold\-CP team\[[25](https://arxiv.org/html/2607.03787#bib.bib25)\]to improve the computational scalability of OpenDDE\. Our immediate goal is to enable efficient parallelization of a single biomolecular sample across multiple GPUs, while also improving GPU utilization by packing multiple smaller samples onto the same GPU when memory allows\.
Building on these system\-level improvements, we plan to scaleOpenDDEbeyond one billion parameters and train the next\-generation model on NVIDIA B300 GPUs\. This will allow us to further test model and data scaling laws for biomolecular foundation models and to extendOpenDDEtoward larger complexes, higher\-resolution all\-atom refinement, and more challenging structure design tasks\.
## 5 Limitations
Despite its strong empirical performance,OpenDDEhas several limitations\. First, our comparison with IsoDDE is necessarily constrained by the limited public information available about the system\. We do not have access to its full training recipe, data mixture, inference procedure, post\-training strategy, or engineering optimizations\. As a result, althoughOpenDDEreaches the level of accuracy reported for IsoDDE on several co\-folding benchmarks, we cannot yet determine which differences arise from model architecture, data processing, scaling, distillation, or inference\-time procedures\. Further work is needed to systematically dissect these factors and to understand which components are essential for robust all\-atom co\-folding\.
Second, this work focuses primarily on protein\-protein and antibody\-antigen interactions rather than protein\-small molecule complexes\. This choice reflects both our scientific objective and our empirical observations\.OpenDDEis designed around biomolecular co\-folding as a structural core for modelling macromolecular interfaces, and our current training and evaluation emphasis follows this setting\. In preliminary analyses, we observed limited correlation between performance on protein\-protein interfaces and protein\-ligand pose prediction, suggesting that these regimes may depend on partly distinct data, inductive biases, and evaluation criteria\. We therefore do not claim that the present model is optimized for ligand docking, virtual screening, or affinity ranking\. ExtendingOpenDDEto these tasks will likely require ligand\-focused data curation, chemistry\-aware post\-training, and task\-specific evaluation\.
Third,OpenDDEshould be viewed as the folding\-centered foundation of an extensible drug\-discovery engine, not as a complete drug\-discovery system\. The present release establishes the structural modelling core and demonstrates that open co\-folding models can scale competitively with closed systems\. However, downstream capabilities such as molecular design, affinity prediction, conformational ensemble modelling, active learning, and experimental feedback remain future extensions rather than fully validated components of the current system\. Similarly, our scaling analyses identify promising directions within the model and data regimes explored here, but they do not exhaust the possible scaling behavior of all\-atom biomolecular models\. Continued community\-driven development, broader benchmarks, and more systematic ablations will be required to turnOpenDDEfrom a strong open folding model into a mature open drug\-discovery engine\.
## Acknowledgements
We thank the AlphaFold, Protenix, and OpenFold teams for releasing their codebases and, in some cases, providing training data resources\. We also thank the NVIDIA infrastructure and architecture teams for developing computing kernels that accelerate model training and inference\.
## Contributions
Hua, Ma for architecture; Hua, Ma for training; Hua, Ma for inference; Ma for data; Lu, Tao for 2026ARK\-AB; Wang for parallel computing; Hua, Ma, Zheng, Wang, Rehman for writing; Hua, Zheng for supervision\. Zheng for organization and conceptualization\.
## References
- Kuntz \[1992\]Irwin D Kuntz\.Structure\-based strategies for drug design and discovery\.*Science*, 257\(5073\):1078–1082, 1992\.
- Santos et al\. \[2017\]Rita Santos, Oleg Ursu, Anna Gaulton, A Patrícia Bento, Ramesh S Donadi, Cristian G Bologa, Anneli Karlsson, Bissan Al\-Lazikani, Anne Hersey, Tudor I Oprea, et al\.A comprehensive map of molecular drug targets\.*Nature reviews Drug discovery*, 16\(1\):19–34, 2017\.
- Varadi et al\. \[2022\]Mihaly Varadi, Stephen Anyango, Mandar Deshpande, Sreenath Nair, Cindy Natassia, Galabina Yordanova, David Yuan, Oana Stroe, Gemma Wood, Agata Laydon, et al\.Alphafold protein structure database: massively expanding the structural coverage of protein\-sequence space with high\-accuracy models\.*Nucleic acids research*, 50\(D1\):D439–D444, 2022\.
- Abramson et al\. \[2024\]Josh Abramson, Jonas Adler, Jack Dunger, Richard Evans, Tim Green, Alexander Pritzel, Olaf Ronneberger, Lindsay Willmore, Andrew J Ballard, Joshua Bambrick, et al\.Accurate structure prediction of biomolecular interactions with alphafold 3\.*Nature*, 630\(8016\):493–500, 2024\.
- Watson et al\. \[2023\]Joseph L Watson, David Juergens, Nathaniel R Bennett, Brian L Trippe, Jason Yim, Helen E Eisenach, Woody Ahern, Andrew J Borst, Robert J Ragotte, Lukas F Milles, et al\.De novo design of protein structure and function with rfdiffusion\.*Nature*, 620\(7976\):1089–1100, 2023\.
- Stark et al\. \[2025\]Hannes Stark, Felix Faltings, MinGyu Choi, Yuxin Xie, Eunsu Hur, Timothy O’Donnell, Anton Bushuiev, Talip Uçar, Saro Passaro, Weian Mao, et al\.Boltzgen: Toward universal binder design\.*bioRxiv*, pages 2025–11, 2025\.
- Zhang et al\. \[2025\]Odin Zhang, Xujun Zhang, Haitao Lin, Cheng Tan, Qinghan Wang, Yuanle Mo, Qiantai Feng, Gang Du, Yuntao Yu, Zichang Jin, et al\.Odesign: A world model for biomolecular interaction design\.*arXiv preprint arXiv:2510\.22304*, 2025\.
- Lu et al\. \[2022\]Wei Lu, Qifeng Wu, Jixian Zhang, Jiahua Rao, Chengtao Li, and Shuangjia Zheng\.Tankbind: Trigonometry\-aware neural networks for drug\-protein binding structure prediction\.*Advances in neural information processing systems*, 35:7236–7249, 2022\.
- Passaro et al\. \[2025\]Saro Passaro, Gabriele Corso, Jeremy Wohlwend, Mateo Reveiz, Stephan Thaler, Vignesh Ram Somnath, Noah Getz, Tally Portnoi, Julien Roy, Hannes Stark, et al\.Boltz\-2: Towards accurate and efficient binding affinity prediction\.*BioRxiv*, 2025\.
- Lu et al\. \[2024\]Wei Lu, Jixian Zhang, Weifeng Huang, Ziqiao Zhang, Xiangyu Jia, Zhenyu Wang, Leilei Shi, Chengtao Li, Peter G Wolynes, and Shuangjia Zheng\.Dynamicbind: predicting ligand\-specific protein\-ligand complex structure with a deep equivariant generative model\.*Nature Communications*, 15\(1\):1071, 2024\.
- Lewis et al\. \[2025\]Sarah Lewis, Tim Hempel, José Jiménez\-Luna, Michael Gastegger, Yu Xie, Andrew YK Foong, Victor García Satorras, Osama Abdin, Bastiaan S Veeling, Iryna Zaporozhets, et al\.Scalable emulation of protein equilibrium ensembles with generative deep learning\.*Science*, 389\(6761\):eadv9817, 2025\.
- Pacesa et al\. \[2025\]Martin Pacesa, Lennart Nickel, Christian Schellhaas, Joseph Schmidt, Ekaterina Pyatova, Lucas Kissling, Patrick Barendse, Jagrity Choudhury, Srajan Kapoor, Ana Alcaraz\-Serna, et al\.One\-shot design of functional protein binders with bindcraft\.*Nature*, 646\(8084\):483–492, 2025\.
- Škrinjar et al\. \[2025\]Peter Škrinjar, Jérôme Eberhardt, Janani Durairaj, and Torsten Schwede\.Have protein\-ligand co\-folding methods moved beyond memorisation?*BioRxiv*, pages 2025–02, 2025\.
- team et al\. \[2024\]Chai Discovery team, Jacques Boitreaud, Jack Dent, Matthew McPartlon, Joshua Meier, Vinicius Reis, Alex Rogozhonikov, and Kevin Wu\.Chai\-1: Decoding the molecular interactions of life\.*BioRxiv*, pages 2024–10, 2024\.
- Wohlwend et al\. \[2025\]Jeremy Wohlwend, Gabriele Corso, Saro Passaro, Noah Getz, Mateo Reveiz, Ken Leidal, Wojtek Swiderski, Liam Atkinson, Tally Portnoi, Itamar Chinn, et al\.Boltz\-1 democratizing biomolecular interaction modeling\.*BioRxiv*, pages 2024–11, 2025\.
- The OpenFold3 Team \[2025\]The OpenFold3 Team\.Openfold3\-preview, 2025\.URL[https://github\.com/aqlaboratory/openfold\-3](https://github.com/aqlaboratory/openfold-3)\.
- Team et al\. \[2025\]ByteDance AML AI4Science Team, Xinshi Chen, Yuxuan Zhang, Chan Lu, Wenzhi Ma, Jiaqi Guan, Chengyue Gong, Jincai Yang, Hanyu Zhang, Ke Zhang, et al\.Protenix\-advancing structure prediction through a comprehensive alphafold3 reproduction\.*BioRxiv*, pages 2025–01, 2025\.
- Candido et al\. \[2026\]Salvatore Candido, Thomas Hayes, Alexander Derry, Roshan Rao, Zeming Lin, Robert Verkuil, Bryan Z Wu, Jin Sub Lee, Elise S Bruguera, Jehan A Keval, et al\.Language modeling materializes a world model of protein biology\.*bioRxiv*, pages 2026–06, 2026\.
- Ma et al\. \[2025\]Wenzhi Ma, Zhenyu Liu, Jincai Yang, Chan Lu, Hanyu Zhang, and Wenzhi Xiao\.From dataset curation to unified evaluation: Revisiting structure prediction benchmarks with pxmeter\.*bioRxiv*, pages 2025–07, 2025\.
- Dunbar et al\. \[2014\]James Dunbar, Konrad Krawczyk, Jinwoo Leem, Terry Baker, Angelika Fuchs, Guy Georges, Jiye Shi, and Charlotte M Deane\.Sabdab: the structural antibody database\.*Nucleic acids research*, 42\(D1\):D1140–D1146, 2014\.
- Zhang et al\. \[2026\]Yuxuan Zhang, Chengyue Gong, Jinyuan Sun, Jiaqi Guan, Milong Ren, Song Xue, Hanyu Zhang, Wenzhi Ma, Zhenyu Liu, Xinshi Chen, et al\.Protenix\-v2: Broadening the reach of structure prediction and biomolecular design\.*bioRxiv*, pages 2026–04, 2026\.
- Zhou et al\. \[2025\]Yi Zhou, Chan Lu, Yiming Ma, Wei Qu, Fei Ye, Kexin Zhang, Lan Wang, Minrui Gui, and Quanquan Gu\.Seedfold: Scaling biomolecular structure prediction\.*arXiv preprint arXiv:2512\.24354*, 2025\.
- Didi et al\. \[2026\]Kieran Didi, Zuobai Zhang, Guoqing Zhou, Danny Reidenbach, Zhonglin Cao, Sooyoung Cha, Tomas Geffner, Christian Dallago, Jian Tang, Michael M Bronstein, et al\.Scaling atomistic protein binder design with generative pretraining and test\-time compute\.*arXiv preprint arXiv:2603\.27950*, 2026\.
- Han et al\. \[2026\]Yewon Han, Maxim I Tsenkov, Niccolò AE Venanzi, Damian Bertoni, Sooyoung Cha, Alejandro Chacón, Nick Dietrich, Boris Fomitchev, Yonathan Goldtzvik, Darren Hsu, et al\.Alphafold database expands to proteome\-scale quaternary structures\.*bioRxiv*, pages 2026–03, 2026\.
- Lin et al\. \[2026\]Dejun Lin, Simon Chu, Vishanth Iyer, Youhan Lee, John St John, Kevin Boyd, Brian Roland, Xiaowei Ren, Guoqing Zhou, Zhonglin Cao, et al\.Fold\-cp: A context parallelism framework for biomolecular modeling\.*arXiv preprint arXiv:2603\.14806*, 2026\.
## Appendix AModel Architecture
OpenDDEis organized as a coarse\-to\-fine molecular structure generation model\. It first encodes input sequence, atom, MSA, template, and constraint features into a residue\-token representation\. A residue\-level trunk, built on Pairformer\-style single and pair updates, performs global relational reasoning over residues and molecular units\. The resulting residue\-token outputs are then expanded into structural tokens, separating chemically meaningful components such as protein backbone, side chain, nucleic\-acid backbone, base, and ligand/atom tokens\. A Structural Refiner further updates these structural\-token representations to model fine\-grained geometric interactions\. Finally, a diffusion\-based coordinate generator denoises atom positions conditioned on the refined structural\-token features, while confidence and distance heads produce residue\-level or structure\-level prediction outputs\.
### A\.1 Reasoning in Biomolecular Foundation Models
In this work, structure reasoning denotes neural relational inference over molecular tokens: the model builds, conditions, and refines latent pairwise representations before using them as geometric constraints for atom\-level denoising\.
We refer to this process as structure reasoning because the model does not directly map sequence tokens to coordinates\. Instead, it first constructs and updates latent relational states that encode hypotheses about molecular geometry\. The pair representationZZstores token\-token structural relations, while the single representationSSstores token\-wise contextual states\. Through structural\-token expansion and refinement,OpenDDEinfers how residue\-level context should be decomposed into backbone, side\-chain, base, and ligand/atom units, and how these units should interact before coordinate generation\.
#### A\.1\.1 Structure Reasoning over Structural Tokens
OpenDDEperforms structure reasoning by converting residue\-level trunk outputs into a structural\-token graph and then applying relational message passing over this graph\. Each structural token is treated as a node, while the pair representation between two structural tokens is treated as an edge:
𝒢s=\(𝒱s,ℰs\),u∈𝒱s,zuvs∈ℰs\.\\mathcal\{G\}^\{s\}=\(\\mathcal\{V\}^\{s\},\\mathcal\{E\}^\{s\}\),\\qquad u\\in\\mathcal\{V\}^\{s\},\\quad z\_\{uv\}^\{s\}\\in\\mathcal\{E\}^\{s\}\.
The Structural Token Expander initializes this graph from the residue\-token trunk:
su0=sπ\(u\)r\+MLP\(sπ\(u\)r\)\+eρ\(u\),s\_\{u\}^\{0\}=s\_\{\\pi\(u\)\}^\{r\}\+\\mathrm\{MLP\}\(s\_\{\\pi\(u\)\}^\{r\}\)\+e\_\{\\rho\(u\)\},zuv0=zπ\(u\),π\(v\)r\+Pρ\(u\),ρ\(v\)\(zπ\(u\),π\(v\)r\)\+b\(ϕuv\),z\_\{uv\}^\{0\}=z\_\{\\pi\(u\),\\pi\(v\)\}^\{r\}\+P\_\{\\rho\(u\),\\rho\(v\)\}\(z\_\{\\pi\(u\),\\pi\(v\)\}^\{r\}\)\+b\(\\phi\_\{uv\}\),whereπ\(u\)\\pi\(u\)is the parent residue of structural tokenuu,ρ\(u\)\\rho\(u\)is its structural role, andϕuv\\phi\_\{uv\}encodes explicit structural relations such as same\-residue membership, backbone–sidechain/base pairing, and polymer backbone adjacency\.
The Structural Refiner then performs reasoning by repeatedly updating node and edge states\. In each layer, structural tokens attend to other structural tokens under pair\-conditioned attention:
αuv\(l\)=softmaxv\(qu\(l\)kv\(l\)⊤d\+g\(zuv\(l\)\)\+auvs\),\\alpha\_\{uv\}^\{\(l\)\}=\\mathrm\{softmax\}\_\{v\}\\left\(\\frac\{q\_\{u\}^\{\(l\)\}\{k\_\{v\}^\{\(l\)\}\}^\{\\top\}\}\{\\sqrt\{d\}\}\+g\(z\_\{uv\}^\{\(l\)\}\)\+a\_\{uv\}^\{s\}\\right\),su\(l\+1\)=su\(l\)\+∑vαuv\(l\)Vv\(l\)\.s\_\{u\}^\{\(l\+1\)\}=s\_\{u\}^\{\(l\)\}\+\\sum\_\{v\}\\alpha\_\{uv\}^\{\(l\)\}V\_\{v\}^\{\(l\)\}\.
The pair representation is also refined through triangular updates:
zuv\(l\+1\)=zuv\(l\)\+TriUpdate\(zuv\(l\),\{zuk\(l\),zkv\(l\)\}k\)\.z\_\{uv\}^\{\(l\+1\)\}=z\_\{uv\}^\{\(l\)\}\+\\mathrm\{TriUpdate\}\\left\(z\_\{uv\}^\{\(l\)\},\\\{z\_\{uk\}^\{\(l\)\},z\_\{kv\}^\{\(l\)\}\\\}\_\{k\}\\right\)\.
This triangular update is important because it allows the model to reason about geometric consistency through indirect paths\. For example, the relation between two tokensuuandvvcan be updated using their shared relations with a third tokenkk\. Thus, the model does not reason about each pair independently; it enforces consistency across many token\-token\-token configurations\.
Therefore, structure reasoning inOpenDDEcan be summarized as:
\(Sr,Zr\)→\(Ss,0,Zs,0\)→\(Ss,L,Zs,L\)→X\.\(S^\{r\},Z^\{r\}\)\\rightarrow\(S^\{s,0\},Z^\{s,0\}\)\\rightarrow\(S^\{s,L\},Z^\{s,L\}\)\\rightarrow X\.
Here,Ss,LS^\{s,L\}andZs,LZ^\{s,L\}are refined structural\-token representations that encode both global residue\-level context and local chemical geometry\. These representations serve as geometric constraints for the downstream diffusion module, which generates atom\-level coordinates\.
#### A\.1\.2 Intuition Behind Structure Reasoning
The intuition behind structure reasoning inOpenDDEis that molecular structure should not be generated as an immediate coordinate regression problem\. Instead, the model first builds an internal relational representation of the molecule, where tokens encode structural units and pair features encode their geometric and chemical dependencies\.
At the residue\-token level, the trunk captures global context from sequence, MSA, templates, chain identity, and pairwise priors\. This stage answers a coarse structural question: which residues or molecular units are likely to be close, constrained, or functionally coupled? However, residue tokens are too coarse to fully describe local geometry, since a single residue contains chemically distinct components such as backbone atoms, side\-chain atoms, or nucleic\-acid bases\.
Structural token expansion introduces a finer representation that exposes these chemically meaningful sub\-units\. A residue is no longer treated as a single object; instead, it is decomposed into structural roles such as backbone, side chain, base, or ligand/atom tokens\. This allows the model to reason over the parts that directly determine atomic geometry\.
The Structural Refiner then performs relational reasoning over this expanded structural\-token graph\. Tokens exchange information through pair\-conditioned attention and triangular pair updates, allowing local structural components to adjust to one another\. For example, a side chain is conditioned not only on its parent residue but also on neighboring backbones, nearby side chains, ligand atoms, and longer\-range pair constraints\. Thus, the model refines a globally consistent structural hypothesis before atom\-level coordinate generation\.
In this view, structure reasoning is a coarse\-to\-fine process:
sequence/context→residue\-level relations→structural\-token relations→atom\-level coordinates\.\\text\{sequence/context\}\\rightarrow\\text\{residue\-level relations\}\\rightarrow\\text\{structural\-token relations\}\\rightarrow\\text\{atom\-level coordinates\}\.
The key benefit is thatOpenDDEseparates the problem of understanding molecular relationships from the problem of generating coordinates\. The trunk and refiner infer geometric constraints in latent space, while the diffusion module uses these refined constraints to produce physically plausible atomic structures\.
#### A\.1\.3 Between Atom37 Representation and Structure Reasoning
The Atom37 representation describes each protein residue using a fixed set of up to 37 atom slots:
Xi37=\{xi,a\}a=137,mi,a∈\{0,1\},X\_\{i\}^\{37\}=\\\{x\_\{i,a\}\\\}\_\{a=1\}^\{37\},\\qquad m\_\{i,a\}\\in\\\{0,1\\\},wherexi,ax\_\{i,a\}is the coordinate of atom slotaain residueii, andmi,am\_\{i,a\}indicates whether this atom exists for the corresponding amino\-acid type\. Thus, Atom37 provides a residue\-local atom layout: it tells the model which atoms should be generated for each residue and how atom identities are consistently indexed across residues\.
However, Atom37 itself does not perform structure reasoning\. It is mainly an atom\-level coordinate representation\. The reasoning happens in the latent token representations before coordinates are produced\. At the residue level,OpenDDEfirst learns contextual single and pair representations:
Sr=\{sir\}i=1Nr,Zr=\{zijr\}i,j=1Nr,S^\{r\}=\\\{s\_\{i\}^\{r\}\\\}\_\{i=1\}^\{N\_\{r\}\},\\qquad Z^\{r\}=\\\{z\_\{ij\}^\{r\}\\\}\_\{i,j=1\}^\{N\_\{r\}\},whereSrS^\{r\}captures residue\-wise context andZrZ^\{r\}captures residue\-residue geometric relations\.
Structure reasoning can be viewed as the process of using these latent representations to infer constraints that will later determine Atom37 coordinates:
\(Sr,Zr\)→geometric constraints→X37\.\(S^\{r\},Z^\{r\}\)\\rightarrow\\text\{geometric constraints\}\\rightarrow X^\{37\}\.
In OpenDDE, this idea is further refined by introducing structural tokens\. Instead of mapping one residue token directly to all Atom37 atom slots, the model decomposes a residue into chemically meaningful structural units, such as backbone and side chain:
i→\{u:π\(u\)=i\},i\\rightarrow\\\{u:\\pi\(u\)=i\\\},whereuuis a structural token andπ\(u\)\\pi\(u\)is its parent residue\.
Under this view, Atom37 is the final atom\-level coordinate container, while structure reasoning is the latent process that decides how those atom slots should be spatially arranged\. The structural\-token branch acts as an intermediate layer between residue\-level reasoning and Atom37\-like atom\-level geometry:
residue tokens→structural tokens→atom coordinates\.\\text\{residue tokens\}\\rightarrow\\text\{structural tokens\}\\rightarrow\\text\{atom coordinates\}\.
Therefore, Atom37 and structure reasoning are complementary\. Atom37 defines the atom\-level output space, whereas structure reasoning provides the learned geometric constraints that make the generated atom coordinates globally and locally consistent\.
### A\.2 Learning theLock\-and\-KeyBiomolecular Interactions
OpenDDEaugments diffusion\-based coordinate generation with geometry\-aware auxiliary objectives\. After residue\-token reasoning and structural\-token refinement, the model predicts all\-atom coordinates
X=\{xa∈ℝ3\}a=1Natom,X=\\\{x\_\{a\}\\in\\mathbb\{R\}^\{3\}\\\}\_\{a=1\}^\{N\_\{\\mathrm\{atom\}\}\},conditioned on structural\-token representations\(Ss,Zs\)\(S^\{s\},Z^\{s\}\)\. The predicted coordinates are then supervised not only by point\-wise coordinate errors, but also by losses that describe interface fitting and local atomic packing\.
#### A\.2\.1 Shape\-Complementarity Losses Guide Interface Fitting
For molecular complexes, accurate fitting requires more than placing chains near each other\. The contacting surfaces should face each other, maintain an appropriate gap, avoid steric clashes, and exhibit complementary surface normals\.OpenDDEintroduces a differentiable shape\-complementarity objective to encode these geometric requirements\.
Given predicted coordinatesXX, each active tokenuuis assigned a centercuc\_\{u\}\. In residue\-token space, this center is usually selected from a representative atom\. In structural\-token space, side\-chain tokens can use the centroid of supervised atoms, allowing the loss to better describe side\-chain surface geometry\.
A local surface normal is estimated from the same\-chain atomic density around the token center:
gu=∑a∈𝒞\(u\)exp\(−‖cu−xa‖22σ2\)cu−xaσ2,nu=gu‖gu‖\+ϵ,g\_\{u\}=\\sum\_\{a\\in\\mathcal\{C\}\(u\)\}\\exp\\left\(\-\\frac\{\\\|c\_\{u\}\-x\_\{a\}\\\|^\{2\}\}\{2\\sigma^\{2\}\}\\right\)\\frac\{c\_\{u\}\-x\_\{a\}\}\{\\sigma^\{2\}\},\\qquad n\_\{u\}=\\frac\{g\_\{u\}\}\{\\\|g\_\{u\}\\\|\+\\epsilon\},where𝒞\(u\)\\mathcal\{C\}\(u\)denotes atoms from the same chain as tokenuu\. This gradient points away from the local molecular density and approximates a surface normal\.
For a cross\-chain token pair\(u,v\)\(u,v\), let
duv=cv−cu,d^uv=duv‖duv‖\+ϵ\.d\_\{uv\}=c\_\{v\}\-c\_\{u\},\\qquad\\hat\{d\}\_\{uv\}=\\frac\{d\_\{uv\}\}\{\\\|d\_\{uv\}\\\|\+\\epsilon\}\.OpenDDEcomputes several geometric factors:
fface\(u,v\)=ReLU\(nu⊤d^uv\)ReLU\(nv⊤\(−d^uv\)\),f\_\{\\mathrm\{face\}\}\(u,v\)=\\mathrm\{ReLU\}\(n\_\{u\}^\{\\top\}\\hat\{d\}\_\{uv\}\)\\,\\mathrm\{ReLU\}\(n\_\{v\}^\{\\top\}\(\-\\hat\{d\}\_\{uv\}\)\),fopp\(u,v\)=1−nu⊤nv2,f\_\{\\mathrm\{opp\}\}\(u,v\)=\\frac\{1\-n\_\{u\}^\{\\top\}n\_\{v\}\}\{2\},fgap\(u,v\)=exp\(−\(‖duv‖−μgapτgap\)2\),f\_\{\\mathrm\{gap\}\}\(u,v\)=\\exp\\left\(\-\\left\(\\frac\{\\\|d\_\{uv\}\\\|\-\\mu\_\{\\mathrm\{gap\}\}\}\{\\tau\_\{\\mathrm\{gap\}\}\}\\right\)^\{2\}\\right\),fclash\(u,v\)=1−σsigmoid\(dclash−‖duv‖τclash\)\.f\_\{\\mathrm\{clash\}\}\(u,v\)=1\-\\sigma\_\{\\mathrm\{sigmoid\}\}\\left\(\\frac\{d\_\{\\mathrm\{clash\}\}\-\\\|d\_\{uv\}\\\|\}\{\\tau\_\{\\mathrm\{clash\}\}\}\\right\)\.
The pairwise shape\-complementarity score is then defined as
quv=fface\(u,v\)fopp\(u,v\)fgap\(u,v\)fclash\(u,v\),q\_\{uv\}=f\_\{\\mathrm\{face\}\}\(u,v\)\\,f\_\{\\mathrm\{opp\}\}\(u,v\)\\,f\_\{\\mathrm\{gap\}\}\(u,v\)\\,f\_\{\\mathrm\{clash\}\}\(u,v\),computed only for valid cross\-chain token pairs within an interface cutoff\. Token\-level and global scores are obtained by pooling these pair scores:
qu=Poolv\(quv\),qglobal=Meanu\(qu\)\.q\_\{u\}=\\mathrm\{Pool\}\_\{v\}\(q\_\{uv\}\),\\qquad q\_\{\\mathrm\{global\}\}=\\mathrm\{Mean\}\_\{u\}\(q\_\{u\}\)\.
During training, the same fields are computed from the ground\-truth coordinates, yieldingquv∗q^\{\*\}\_\{uv\},qu∗q^\{\*\}\_\{u\}, andqglobal∗q^\{\*\}\_\{\\mathrm\{global\}\}\. The shape\-complementarity objective is:
ℒshape=λpHuber\(quv,quv∗\)\+λtHuber\(qu,qu∗\)\+λgHuber\(qglobal,qglobal∗\)\.\\mathcal\{L\}\_\{\\mathrm\{shape\}\}=\\lambda\_\{p\}\\mathrm\{Huber\}\(q\_\{uv\},q^\{\*\}\_\{uv\}\)\+\\lambda\_\{t\}\\mathrm\{Huber\}\(q\_\{u\},q^\{\*\}\_\{u\}\)\+\\lambda\_\{g\}\\mathrm\{Huber\}\(q\_\{\\mathrm\{global\}\},q^\{\*\}\_\{\\mathrm\{global\}\}\)\.
This loss guides the model toward interface geometries that are not only close in Euclidean distance but also physically compatible in surface orientation, spacing, and anti\-clash behavior\. It is especially useful for protein\-protein, protein\-ligand, and protein\-nucleic\-acid fitting, where correct binding geometry depends on surface complementarity rather than isolated atom placement\.
#### A\.2\.2 All\-Atom Refinement Improves Local Packing
OpenDDEfurther improves local packing by supervising coordinates in all\-atom and local\-frame spaces\. The diffusion module predicts atom coordinates conditioned on refined structural\-token features:
X^0=Dθ\(Xt,t∣Ss,Zs,𝒜\),\\hat\{X\}\_\{0\}=D\_\{\\theta\}\(X\_\{t\},t\\mid S^\{s\},Z^\{s\},\\mathcal\{A\}\),where𝒜\\mathcal\{A\}contains atom\-to\-token and atom\-to\-structural\-token mappings\. This design allows each atom to receive both global pair context and local chemical\-role information\.
A standard global coordinate loss can penalize atom displacement, but it may not fully capture residue\-internal geometry\. Therefore,OpenDDEintroduces local\-frame supervision\. For each residue or structural token with a valid frameTuT\_\{u\}, predicted and ground\-truth atoms are expressed in the same local frame:
x~apred=Tu−1x^a,x~atrue=Tu−1xa∗\.\\tilde\{x\}^\{\\mathrm\{pred\}\}\_\{a\}=T\_\{u\}^\{\-1\}\\hat\{x\}\_\{a\},\\qquad\\tilde\{x\}^\{\\mathrm\{true\}\}\_\{a\}=T\_\{u\}^\{\-1\}x^\{\*\}\_\{a\}\.The local coordinate loss is:
ℒlocal=1\|𝒜u\|∑a∈𝒜uSmoothL1\(x~apred,x~atrue\)\.\\mathcal\{L\}\_\{\\mathrm\{local\}\}=\\frac\{1\}\{\|\\mathcal\{A\}\_\{u\}\|\}\\sum\_\{a\\in\\mathcal\{A\}\_\{u\}\}\\mathrm\{SmoothL1\}\\left\(\\tilde\{x\}^\{\\mathrm\{pred\}\}\_\{a\},\\tilde\{x\}^\{\\mathrm\{true\}\}\_\{a\}\\right\)\.
For protein residues,OpenDDEalso supervises side\-chain torsions:
ℒχ=1−cos\(χ^−χ∗\),\\mathcal\{L\}\_\{\\chi\}=1\-\\cos\(\\hat\{\\chi\}\-\\chi^\{\*\}\),which directly encourages correct rotameric geometry\. For structural tokens, side\-chain and nucleic\-acid base atoms are supervised in frames inherited from their parent backbone token:
ℒsc/base=∑u∈𝒱sc/basewρ\(u\)1\|𝒜u\|∑a∈𝒜uSmoothL1\(Tu−1x^a,Tu−1xa∗\)\.\\mathcal\{L\}\_\{\\mathrm\{sc/base\}\}=\\sum\_\{u\\in\\mathcal\{V\}\_\{\\mathrm\{sc/base\}\}\}w\_\{\\rho\(u\)\}\\frac\{1\}\{\|\\mathcal\{A\}\_\{u\}\|\}\\sum\_\{a\\in\\mathcal\{A\}\_\{u\}\}\\mathrm\{SmoothL1\}\\left\(T\_\{u\}^\{\-1\}\\hat\{x\}\_\{a\},T\_\{u\}^\{\-1\}x^\{\*\}\_\{a\}\\right\)\.
The full geometry objective can be summarized as:
ℒgeom=ℒdiff\+λshapeℒshape\+λlocalℒlocal\+λsc/baseℒsc/base\+λχℒχ\.\\mathcal\{L\}\_\{\\mathrm\{geom\}\}=\\mathcal\{L\}\_\{\\mathrm\{diff\}\}\+\\lambda\_\{\\mathrm\{shape\}\}\\mathcal\{L\}\_\{\\mathrm\{shape\}\}\+\\lambda\_\{\\mathrm\{local\}\}\\mathcal\{L\}\_\{\\mathrm\{local\}\}\+\\lambda\_\{\\mathrm\{sc/base\}\}\\mathcal\{L\}\_\{\\mathrm\{sc/base\}\}\+\\lambda\_\{\\chi\}\\mathcal\{L\}\_\{\\chi\}\.
This all\-atom refinement pathway improves local packing because it supervises atoms in chemically meaningful coordinate systems\. Backbone, side\-chain, base, and ligand/atom tokens are not treated as anonymous points; they are connected through structural roles, parent\-residue frames, atom mappings, and pairwise geometric context\. As a result, the model learns to place side chains, bases, and interface atoms in locally consistent configurations before producing the final all\-atom structure\.
### A\.3 Structure Prediction andDe NovoDesign in One Framework
OpenDDEformulates structure prediction and de novo design as two instances of the same conditional denoising problem\. The key idea is to represent the input complex as a partially observed molecular system, where some atoms may be treated as fixed structural context and the remaining atoms are treated as prediction or design targets\.
For each training example, we define two atom masks:
maknown∈\{0,1\},matarget∈\{0,1\},m^\{\\mathrm\{known\}\}\_\{a\}\\in\\\{0,1\\\},\\qquad m^\{\\mathrm\{target\}\}\_\{a\}\\in\\\{0,1\\\},wheremaknown=1m^\{\\mathrm\{known\}\}\_\{a\}=1indicates atoms provided as structural condition, andmatarget=1m^\{\\mathrm\{target\}\}\_\{a\}=1indicates atoms to be generated or predicted\. In standard structure prediction, no chain is fixed:
mknown=0,mtarget=mresolved\.m^\{\\mathrm\{known\}\}=0,\\qquad m^\{\\mathrm\{target\}\}=m^\{\\mathrm\{resolved\}\}\.In conditional design, one or more chains, motifs, pockets, or structural contexts are fixed:
mknown≠0,mtarget=mresolved−mknown\.m^\{\\mathrm\{known\}\}\\neq 0,\\qquad m^\{\\mathrm\{target\}\}=m^\{\\mathrm\{resolved\}\}\-m^\{\\mathrm\{known\}\}\.
Thus, the same model can be trained on both tasks by changing only the data condition masks, rather than introducing a separate architecture\.
#### A\.3\.1 Known\-Chain Data Conditioning
During training,OpenDDEsamples partial\-complex conditioning cases from multi\-chain structures\. A polymer chain is selected as known structural context, while the remaining resolved atoms are treated as target atoms\. This produces the conditioning fields
𝒞=\{mknown,mtarget,is\_conditioned\}\.\\mathcal\{C\}=\\\{m^\{\\mathrm\{known\}\},m^\{\\mathrm\{target\}\},\\mathrm\{is\\\_conditioned\}\\\}\.
The input feature dictionary therefore contains both ordinary structure prediction features and conditional\-generation features:
ℐ=\{sequence,MSA,template,atomfeatures,tokenfeatures,𝒞\}\.\\mathcal\{I\}=\\\{\\mathrm\{sequence\},\\mathrm\{MSA\},\\mathrm\{template\},\\mathrm\{atom\\ features\},\\mathrm\{token\\ features\},\\mathcal\{C\}\\\}\.
The residue\-token trunk and structural\-token branch encode the full molecular system:
ℐ→\(Sr,Zr\)→\(Ss,Zs\),\\mathcal\{I\}\\rightarrow\(S^\{r\},Z^\{r\}\)\\rightarrow\(S^\{s\},Z^\{s\}\),where\(Sr,Zr\)\(S^\{r\},Z^\{r\}\)are residue\-token representations and\(Ss,Zs\)\(S^\{s\},Z^\{s\}\)are structural\-token representations\. The condition masks do not require a separate design model; they only change how noise, alignment, and loss are applied\.
#### A\.3\.2 Conditional Diffusion Objective
LetX0X\_\{0\}be the ground\-truth all\-atom coordinates\.OpenDDEfirst applies a shared random rigid augmentation:
X¯0=RX0\+t\.\\bar\{X\}\_\{0\}=RX\_\{0\}\+t\.For ordinary structure prediction, Gaussian noise is added to all supervised atoms:
Xt=X¯0\+σtϵ\.X\_\{t\}=\\bar\{X\}\_\{0\}\+\\sigma\_\{t\}\\epsilon\.
For data\-conditioned design, noise is added only to target atoms, while known atoms remain fixed as structural context:
Xt=mknown⊙X¯0\+mtarget⊙\(X¯0\+σtϵ\)\.X\_\{t\}=m^\{\\mathrm\{known\}\}\\odot\\bar\{X\}\_\{0\}\+m^\{\\mathrm\{target\}\}\\odot\(\\bar\{X\}\_\{0\}\+\\sigma\_\{t\}\\epsilon\)\.
The denoising network predicts coordinates conditioned on both the latent representations and the partially observed structure:
X^0=Dθ\(Xt,t∣Ss,Zs,ℐ,mknown\)\.\\hat\{X\}\_\{0\}=D\_\{\\theta\}\\left\(X\_\{t\},t\\mid S^\{s\},Z^\{s\},\\mathcal\{I\},m^\{\\mathrm\{known\}\}\\right\)\.
The training loss is then applied primarily on the target atoms:
ℒdiff=𝔼t,ϵ\[‖mtarget⊙\(X^0−X¯0\)‖2\]\.\\mathcal\{L\}\_\{\\mathrm\{diff\}\}=\\mathbb\{E\}\_\{t,\\epsilon\}\\left\[\\left\\\|m^\{\\mathrm\{target\}\}\\odot\\left\(\\hat\{X\}\_\{0\}\-\\bar\{X\}\_\{0\}\\right\)\\right\\\|^\{2\}\\right\]\.
This objective naturally supports both tasks\. Whenmknown=0m^\{\\mathrm\{known\}\}=0, the model learns standard structure prediction\. Whenmknown≠0m^\{\\mathrm\{known\}\}\\neq 0, the model learns to generate missing or designable parts conditioned on a fixed molecular context\.
#### A\.3\.3 Unified Prediction\-Design Framework
The unified framework can be summarized as:
full complex data→condition sampling→\(mknown,mtarget\)→conditional diffusion training\.\\text\{full complex data\}\\rightarrow\\text\{condition sampling\}\\rightarrow\(m^\{\\mathrm\{known\}\},m^\{\\mathrm\{target\}\}\)\\rightarrow\\text\{conditional diffusion training\}\.
For structure prediction, the model receives sequence and feature information and predicts the full structure:
pθ\(X∣ℐ\)\.p\_\{\\theta\}\(X\\mid\\mathcal\{I\}\)\.
For de novo design, the model receives an additional structural condition and generates the target region:
pθ\(Xtarget∣ℐ,Xknown\)\.p\_\{\\theta\}\\left\(X^\{\\mathrm\{target\}\}\\mid\\mathcal\{I\},X^\{\\mathrm\{known\}\}\\right\)\.
Both distributions are implemented by the same denoising model:
Dθ\(Xt,t∣Ss,Zs,ℐ,mknown\)\.D\_\{\\theta\}\(X\_\{t\},t\\mid S^\{s\},Z^\{s\},\\mathcal\{I\},m^\{\\mathrm\{known\}\}\)\.
This design has three advantages\. First, prediction and design share the same representation learning pipeline, including residue\-token reasoning, structural\-token refinement, and all\-atom diffusion\. Second, the model learns from ordinary structure data without requiring a separate curated design dataset\. Third, conditional design becomes a masked generation problem: fixed chains, motifs, or binding contexts are provided as data conditions, while the model denoises and generates the target molecular components around them\.
In this sense,OpenDDEdoes not treat structure prediction and de novo design as separate tasks\. Instead, both are unified as conditional structure generation under different choices of the known\-target partition\.
## Appendix BExperiment
### B\.1 Overall Performance
Overall,OpenDDEachieves strong and balanced performance across the FoldBench tasks \(Figure[8](https://arxiv.org/html/2607.03787#A2.F8)\), with particularly large gains on interaction\-heavy systems\. It reaches the best result on protein monomer LDDT, reaching 0\.890, and obtains the highest protein\-protein score at 0\.769\. The largest improvement appears on antigen\-antibody complexes, whereOpenDDEreaches 0\.700, substantially outperforming ESMFold2 \(0\.581\), AlphaFold3 \(0\.488\), and Protenix\-v1 \(0\.485\)\.OpenDDEalso remains competitive on RNA monomers \(0\.660\), protein\-ligand complexes \(0\.601\), and protein\-RNA complexes \(0\.735\), where it outperforms or is close to the strongest baseline\. These results suggest thatOpenDDEpreserves high single\-chain structure accuracy while improving the more challenging cross\-molecular reasoning needed for protein\-protein and antibody\-antigen interfaces\.
Figure 8:Structure prediction performance on FoldBench\.In the two monomer panels, the metric is mean LDDT\. In the protein\-protein, antigen\-antibody, and protein\-RNA panels, the metric is DockQ success rate\. In the protein\-ligand panel, the metric is ligand pose success rate under the corresponding FoldBench criterion\.
### B\.2 Test\-Time Scaling
Figure 9:Test\-time scaling on FoldBench\-AB and 2026ARK\-AB\.Top:Ranking\-based selection\.Middle:Ranking\-based selection of DockQ\>\>0\.8\.Bottom: Oracle\-based selection\.We evaluate whether additional inference\-time sampling improves antibody\-antigen prediction performance \(Figure[9](https://arxiv.org/html/2607.03787#A2.F9)\)\. For each target, we generate multiple stochastic samples with different random seeds and select the final prediction using the model ranking score\. We report three settings: ranking\-based success rate, ranking\-based high\-quality success rate with DockQ\>\>0\.8, and oracle success rate, where the best candidate among all generated samples is selected using the ground\-truth DockQ\. The oracle setting measures the maximal performance available in the sampled candidate set and separates sampling capacity from ranking quality\.
As shown in Figure[9](https://arxiv.org/html/2607.03787#A2.F9),OpenDDEexhibits consistent test\-time scaling on both FoldBench\-AB and 2026ARK\-AB\. On FoldBench\-AB, increasing the number of seeds improves the ranked success rate from roughly 62% with one seed to about 73% at large sampling budgets, while the ranked high\-DockQ rate increases from about 28% to 34%\. Similar behavior is observed on 2026ARK\-AB, whereOpenDDEimproves from approximately 63% to 69% ranked success and maintains the strongest high\-DockQ performance across seed counts\. These gains indicate that additional samples expose better antibody\-antigen conformations rather than simply producing redundant predictions\.
The oracle curves show a larger scaling effect\. On FoldBench\-AB,OpenDDE’s oracle success rises from about 66% with one seed to above 90% with hundreds of seeds\. On 2026ARK\-AB, the oracle success similarly increases from about 67% to nearly 90%\. This gap between ranked and oracle performance suggests thatOpenDDEoften generates high\-quality antibody\-antigen structures that are not always selected as the top\-ranked prediction\. Therefore, test\-time scaling reveals strong latent sampling capacity, and further improvements in confidence calibration or candidate ranking could convert more of this oracle performance into realized prediction accuracy\.
Across all three metrics,OpenDDEscales more favorably with inference\-time compute than the compared open models, including Boltz\-1, OpenFold3, and Protenix\-v1\. These results show thatOpenDDEbenefits predictably from larger sampling budgets and that its diffusion\-based generation process provides a practical route for improving antibody\-antigen prediction accuracy at test time\.
## Appendix COpenDDEParallel Computing
OpenDDEintegrates Fold\-CP\[[25](https://arxiv.org/html/2607.03787#bib.bib25)\]to enable memory\-efficient inference for large biomolecular systems\. The dominant memory bottleneck inOpenDDEarises from token\-pair representations, whose size scales quadratically with the number of tokens\. Given a pair representation
Z∈ℝN×N×C,Z\\in\\mathbb\{R\}^\{N\\times N\\times C\},Fold\-CP partitions the two token dimensions across a two\-dimensional context\-parallel mesh\. For a mesh rank\(p,q\)\(p,q\), the local pair block is
Zpq=Z\[ℐp,ℐq,:\],Z\_\{pq\}=Z\[\\mathcal\{I\}\_\{p\},\\mathcal\{I\}\_\{q\},:\],whereℐp\\mathcal\{I\}\_\{p\}andℐq\\mathcal\{I\}\_\{q\}denote the token ranges assigned to the row and column dimensions of that rank\. This decomposition avoids replicating the fullN×NN\\times Npair tensor on every GPU, reducing per\-device memory while preserving the original model computation\.
OpenDDE applies this context\-parallel layout to modules dominated by token\-pair operations, including Pairformer\-style trunk blocks, structural\-token refinement, diffusion conditioning, and confidence prediction\. Pair transitions, triangular updates, and pair\-biased attention are evaluated on local pair blocks, with the required communication performed along the row and column axes of the context\-parallel mesh\. Thus, each device stores and updates only a shard of the pair representation, while collective communication reconstructs the information needed for operations that depend on shared row\-wise, column\-wise, or triangular context\.
For attention with pair bias, the Fold\-CP path follows blockwise exact attention computation, closely related to online softmax and memory\-efficient attention\. Rather than materializing the full attention score matrix,OpenDDEprocesses key, value, and pair\-bias blocks incrementally across context\-parallel ranks\. The softmax statistics are accumulated block by block, so the resulting attention output is equivalent to standard softmax attention, while peak memory is reduced\. This is particularly important for biomolecular systems because pair\-derived attention biases scale with the same quadratic token\-pair dimension\.
The Fold\-CP integration is an execution strategy rather than a different model\. It does not change model weights, architecture, or inference semantics\. Instead, it distributes the sameOpenDDEcomputation across multiple GPUs, enabling larger complexes and longer token contexts to be processed under practical memory limits\.
Table 3:Runtime and memory comparison between single\-GPU inference and CP4 inference\.Forward Time measures the model forward pass only\. End\-to\-End Time measures the full wall\-clock runtime, including setup, input processing, model execution, post\-processing, and output writing\. Peak GPU Memory is the maximum physical GPU memory observed by external monitoring\. Torch Stage Peak Memory is the peak memory recorded by PyTorch stage\-level instrumentation\.Similar Articles
Co-folding model guided by structural proteomics
Introduces AIMS-Fold, an inference-time guided-diffusion framework that integrates cross-linking mass spectrometry (XL-MS) and hydrogen-deuterium exchange (HDX-MS) data to improve protein co-folding predictions for induced proximity drug targets.
ToolMol: Evolutionary Agentic Framework for Multi-objective Drug Discovery
ToolMol is an evolutionary agentic framework that combines a multi-objective genetic algorithm with an LLM-based operator to design small-molecule drugs, achieving state-of-the-art binding affinity and drug-likeness on multiple protein targets.
From Holo Pockets to Electron Density: GPT-style Drug Design with Density
This paper introduces EDMolGPT, an autoregressive framework that generates 3D molecular conformations from low-resolution electron density point clouds, improving structure-based drug design by leveraging physically meaningful density signals.
DrugSAGE:Self-evolving Agent Experience for Efficient State-of-the-Art Drug Discovery
DrugSAGE is a framework that accumulates and reuses cross-task memory to build state-of-the-art drug discovery models efficiently, outperforming baseline agents by 10-30% on held-out tasks.
SPADE: Faster Drug Discovery by Learning from Sparse Data
This paper introduces SPADE, a novel algorithm for drug discovery that efficiently identifies high-quality ligands from sparse data using only ~40 tests. It demonstrates superior sample efficiency and speed compared to deep learning and Bayesian optimization methods.