Diffusion-based Evolutionary Optimization for 3D Multi-Objective Molecular Generation
Organizations: College of Computer Science and Technology, National University of Defense Technology, 410073, Changsha, China · Bioinformatics Center, Academy of Military Medical Sciences (AMMS), 100080, Beijing, China
Abstract
De novo 3D molecular design with multiple properties and fragment-retention constraints is formulated as a constrained multi-objective optimization problem. Existing 3D diffusion models provide strong capabilities for generating novel molecules, but adapting them to new conditions and multi-objective optimization tasks remains difficult. We present Diffusion-based Evolutionary Molecular Optimization (DEMO), a framework that draws solely on the unconditional molecular prior of a frozen pretrained diffusion model to generate chemically valid 3D candidates, while evolutionary search handles conditional generation and constrained multi-objective optimization, without retraining or fine-tuning. DEMO performs crossover and mutation on intermediate diffusion states through an Evolutionary-Guided Diffusion (EGD) operator, combines 2D topology and 3D geometry in Structure-Aware Environmental Selection (SAES) to preserve structural diversity, and uses three cooperating populations to coordinate structural exploration, fragment-constraint satisfaction, and property optimization. Experimental results demonstrate that EGD outperforms task-specific conditional generators in property targeting; SAES preserves structural diversity while retaining competitive hypervolume in multi-objective optimization; and DEMO improves objective convergence and generates diverse feasible molecules under fragment constraints. These results show how evolutionary search can reuse existing diffusion models for diverse molecular design requirements.
Figures & tables
| Method | ||||||
|---|---|---|---|---|---|---|
| cEDM | 2.76 | 655 | 356 | 584 | 1.11 | 1.10 |
| cGeoLDM | 2.37 | 587 | 340 | 522 | 1.11 | 1.03 |
| EEGSDE | 2.50 | 487 | 302 | 447 | 0.78 | 0.94 |
| cGCDM | 1.99 | 595 | 346 | 480 | 0.86 | 0.70 |
| PropMolFlow | 1.36 | 391 | 246 | 312 | 0.63 | 0.64 |
| NExT-Mol | 1.16 | 297 | 205 | 235 | 0.51 | 0.51 |
| Task | Prop. | Baselines | Top-N | EGD & Ablation | |||||
|---|---|---|---|---|---|---|---|---|---|
| MUDM | EEGSDE | Base | (SE) | (SR) | w/o CO | w/o MT | Ours | ||
| Backbone: GeoLDM | |||||||||
| 1.47 | 0.98 | 1.23 | 1.21 | 2.25 | 0.64 | 1.67 | 0.69 | ||
| 0.69 | 0.91 | 1.12 | 0.34 | 0.59 | 0.28 | 0.29 | 0.09 | ||
| 544 | 563 | 664 | 234 | 461 | 78 | 228 | 186 | ||
| 0.58 | 0.87 | 1.13 | 0.31 | 0.54 | 0.25 | 0.27 | 0.48 | ||
| Task | HV D | HV 100K | GeoLDM | EDM | ||||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| Top-N | EGD Ablation & Selection | Top-N | EGD Ablation & Selection | |||||||||||||
| (SR) | (SE) | w/o CO | w/o MT | SPEA2 | SPEA2-SD | SAES | (SR) | (SE) | w/o CO | w/o MT | SPEA2 | SPEA2-SD | SAES | |||
| – | 0.854 | 0.802 | 0.559 (0.017) | 0.590 (0.024) | 0.568 (0.054) | 0.760 (0.015) | 0.758 (0.019) | 0.759 (0.017) | 0.764 (0.016) | 0.535 (0.016) | 0.629 (0.029) | 0.545 (0.042) | 0.761 (0.012) | 0.766 (0.012) | 0.769 (0.016) | 0.772 (0.019) |
| [0.997] | [1.000] | [0.734] | [0.895] | [0.312] | [0.278] | [0.896] | [1.000] | [1.000] | [0.698] | [0.932] | [0.300] | [0.277] | [0.880] | |||
| – | 0.702 | 0.678 | 0.461 (0.020) | 0.543 (0.015) | 0.489 (0.042) | 0.644 (0.019) | 0.635 (0.018) | 0.637 (0.027) | 0.639 (0.018) | 0.475 (0.020) | 0.536 (0.026) | 0.478 (0.036) | 0.641 (0.020) | 0.640 (0.017) | 0.648 (0.024) | 0.640 (0.017) |
| [0.982] | [1.000] | [0.776] | [0.962] | [0.251] | [0.247] | [0.876] | [1.000] | [1.000] | [0.725] | [0.951] | [0.276] | [0.231] | [0.871] | |||
| Target | Top-N (SE) | Evolutionary Framework | |
|---|---|---|---|
| Unconditional | SPEA2 | SAES (Ours) | |
| Backbone: GeoLDM | |||
| 4tqr | 0.718 (.080) - 1.00 | 0.798 (.095) - 0.84 | 0.793 (.081) - 0.97 |
| 1k9t | 0.339 (.038) - 0.60 | 0.338 (.028) - 0.47 | 0.328 (.018) - 0.64 |
| 1ai4 | 0.458 (.033) - 0.95 | 0.481 (.043) - 0.71 | 0.490 (.049) - 0.92 |
| 5mgl | 0.485 (.039) - 0.91 | 0.517 (.043) - 0.65 | 0.519 (.038) - 0.87 |
| Task | Top-N | Traditional CMOEAs | Ablation & DEMO | Task | Top-N | Traditional CMOEAs | Ablation & DEMO | ||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| (SE) | SPEA2 | CCMO | CMOEA-CD | w/o Pop B | Only Pop C | Ours | (SE) | SPEA2 | CCMO | CMOEA-CD | w/o Pop B | Only Pop C | Ours | ||
| Backbone: GeoLDM | Backbone: EDM | ||||||||||||||
| — | 0.035 (.085) | 0.352 (.023) | 0.343 (.028) | 0.338 (.033) | 0.399 (.042) | 0.386 (.062) | 0.394 (.036) | — | 0.104 (.014) | 0.141 (.137) | 0.111 (.133) | 0.093 (.012) | 0.091 (.143) | 0.182 (.136) | 0.143 (.147) |
| [0.04/0.01] | [1.00/0.13] | [1.00/0.17] | [1.00/0.18] | [1.00/0.46] | [0.15/0.11] | [1.00/0.53] | [0.01/0.01] | [0.47/0.05] | [0.31/0.03] | [0.31/0.02] | [0.24/0.03] | [0.03/0.03] | [0.44/0.03] | ||
| - – | 0 (.0) | 0.015 (.069) | 0.011 (.048) | 0.025 (.077) | 0.024 (.081) | 0 (.0) | 0.025 (.067) | - – | 0.251 (.068) | 0.297 (.084) | 0.316 (.038) | 0.297 (.074) | 0.358 (.047) | 0.358 (.044) | 0.356 (.044) |
| [0.00/0.00] | [0.05/0.01] | [0.05/0.01] | [0.10/0.01] | [0.01/0.01] | [0.00/0.00] | [0.10/0.01] | [0.12/0.11] | [0.95/0.11] | [0.96/0.12] | [0.95/0.15] | [1.00/0.34] | [0.11/0.11] | [1.00/0.42] | ||
3D Molecular Generation
3D molecular generation creates molecules by modeling their chemical structures and three-dimensional geometries.
See all 56 papers