Source-linked AI summary
Text-guided multi-property molecular optimization with a diffusion language model
Yida Xiong, Kun Li, Jiameng Chen, Hongzhi Zhang, Di Lin, Yan Che, Wenbin Hu
TL;DR
Existing molecular optimization often relies on property predictors whose errors can accumulate during iterative search. TransDLM guides diffusion with IUPAC-based molecular descriptions and property requirements, and the paper reports stronger benchmark results while maintaining structural similarity.
Problem
Many molecular optimization methods rely on external property predictors, whose approximations can introduce errors that accumulate during iterative search.
Source: §1 · p. 2Method
TransDLM uses a transformer-based diffusion language model that conditions optimization on encoded source molecules and textual descriptions of molecular structure and desired properties.
Source: §1 · p. 2Results
TransDLM outperformed other methods on most structural metrics and all ADMET property criteria, while a case study optimized XAC selectivity from A2AR to A1R.
Source: §1 · p. 2, §7 · p. 10Takeaways & Limitations
The reported results support TransDLM's ability to combine structural retention with property enhancement in benchmark optimization and a ligand–receptor selectivity case study.
Source: §1 · p. 2, §7 · p. 10The multi-stage pipeline adds system complexity and computational overhead, and reliance on external tools such as IUPAC conversion complicates large-scale deployment.
Source: §6 · p. 9
Abstract
from arXiv · showhide
Molecular optimization (MO) is a crucial stage in drug discovery in which task-oriented generated molecules are optimized to meet practical industrial requirements. Existing mainstream MO approaches primarily utilize external property predictors to guide iterative property optimization. However, learning all molecular samples in the vast chemical space is unrealistic for predictors. As a result, errors and noise are inevitably introduced during property prediction due to the nature of approximation. This leads to discrepancy accumulation, generalization reduction and suboptimal molecular candidates. In this paper, we propose a text-guided multi-property molecular optimization method utilizing transformer-based diffusion language model (TransDLM). TransDLM leverages standardized chemical nomenclature as semantic representations of molecules and implicitly embeds property requirements into textual descriptions, thereby mitigating error propagation during diffusion process. By fusing physically and chemically detailed textual semantics with specialized molecular representations, TransDLM effectively integrates diverse information sources to guide precise optimization, which enhances the model's ability to balance structural retention and property enhancement. Additionally, the success of a case study further demonstrates TransDLM's ability to solve practical problems. Experimentally, our approach surpasses state-of-the-art methods in maintaining molecular structural similarity and enhancing chemical properties on the benchmark dataset.
1 Introduction
TransDLM uses molecular descriptions to guide diffusion-based optimization while retaining source structure and avoiding reliance on external property predictors. The paper reports stronger benchmark performance and a case study optimizing XAC selectivity.
1 Introduction
- Industrial molecular optimization requires substantial property gains with only slight changes to molecular architecture, a balance many existing approaches neglect.The introduction also identifies scaffold retention and the scale of chemical space as challenges for manual optimization.Source: §1 · pp. 1–2
- External property predictors can accumulate approximation errors during iterative search, producing inaccurate predictions and suboptimal candidates.The cited discussion also notes that noise may divert the search from optimal regions, especially in novel chemical spaces.Source: §1 · p. 2
- TransDLM starts from encoded source molecules and uses IUPAC-based textual descriptions to guide diffusion, embedding property requirements without an external predictor.The descriptions combine molecular semantics and desired properties, with the stated aim of mitigating error propagation while retaining core scaffolds.Source: §1 · p. 2
- TransDLM outperforms state-of-the-art methods on LogD, Solubility, and Clearance while maintaining structural similarity, and biases XAC selectivity from A2AR to A1R.The reported benchmark result covers three ADMET properties; the selectivity result is a separate case study.Source: §1 · p. 2
2 Related work
The related-work discussion covers molecular optimization strategies and diffusion approaches to language generation.
2.1 Molecular optimization
- Guided search methods explore chemical or latent spaces using property predictors, while molecular-mapping approaches learn transformations from matched molecular pairs.Examples include JTVAE and MolDQN for guided search, and Seq2Seq, DESMILES, and BT4MolGen for molecular mapping.Source: §2.1 · pp. 2–3
- Guided search relies on external property predictors, whose errors and noise can affect iterative optimization.Mapping-based methods instead use molecular transformation rules or learned source-to-target mappings.Source: §2.1 · pp. 2–3
2.2 Language generation with diffusion models
- Language diffusion methods either preserve text's discrete form or embed it in a continuous representation; this work follows the embedding-based strategy for SMILES optimization.The review contrasts discrete diffusion approaches with methods that map text into continuous representations.Source: §2.2 · pp. 3–4
3 Methodology
TransDLM conditions diffusion on molecular structure and textual property descriptions to transform a source molecule into an optimized SMILES output.
3.1 Overall framework
- The framework takes an initial molecule and description as input, then combines embedding, noising, denoising, and rounding to produce an optimized molecule.The denoising process is guided by the textual description and reconstructs the desired molecular representation.Source: §3.1 · p. 4
3.2 Description generator
- RECAP fragments and IUPAC names encode source structure; descriptions combine these details with original property values and desired property changes.The description template also specifies fragment connections and bond types.Source: §3.2 · p. 4
3.3 SMILES tokenizer
- The tokenizer preserves multi-character SMILES semantic groups as tokens and pads each sequence to length n.The paper gives [CH3-] as an example of a group retained as one token.Source: §3.3 · p. 5
3.4 Text-guided molecular optimization
- TransDLM samples from encoded source molecules and uses cross-attention to condition denoising on descriptions before rounding the predicted embedding to SMILES.The model predicts the target embedding from the noisy state and description embedding.Source: §3.4 · p. 5
3.5 TransDLM training
- Training maximizes a variational lower bound while teaching the model to reconstruct the desired data at each denoising step.The paper describes this as learning to reverse diffusion using the data structure and guiding text.Source: §3.5 · p. 5
4 Theoretical analysis
The analysis contrasts iterative predictor-guided optimization, where prediction errors accumulate, with TransDLM’s diffusion-based error suppression. It concludes that TransDLM mitigates error propagation and supports a more stable, reliable multi-property optimization process.
4.1 Error propagation in traditional methods
- Under a Hessian spectral-radius bound H, traditional methods’ cumulative projected error grows linearly with iterations Z, potentially deviating the final optimization direction from its theoretical optimum.The error at each iteration comes from an external predictor’s estimate of the true property.Source: §4.1 · pp. 5–6
4.2 Error suppression mechanism of TransDLM
- Diffusion scheduling attenuates each TransDLM error, and the analysis attributes effective cumulative suppression to convergence of the γ_s series.The process is described as Markovian: each step depends on the current state and the currently predicted x0.Source: §4.2 · p. 6
4.3 Comparative analysis
- For traditional methods relying on external property predictors, total variance grows quadratically with optimization steps Z, whereas TransDLM’s diffusion framework suppresses errors, keeping accumulated errors significantly lower.Source: §4.3 · p. 6
5 Experiments
Experiments on the MMP dataset compare TransDLM with five baselines across structural and property metrics, then assess its design choices and an XAC binding-affinity case study. TransDLM mostly outperformed the baselines, while ablations and the case study illustrate trade-offs and performance beyond the ADMET training setting.
5.1 Dataset
- The study uses 198558 source–target molecule pairs with ADMET properties, split into training, validation, and test sets.The whole dataset was split 90% for training and validation and 10% for testing; the former was further split 90%/10% for training and validation.Source: §5.1 · p. 6
5.2 Metrics
- Eight metrics assess structural construction and six criteria assess property optimization.These include similarity and validity measures, ADMET property-change accuracy, and Pareto-front indicators for QED and SA.Source: §5.2 · p. 6
5.3 Baselines
- Five baselines span autoregressive, deep generative, and VAE models, and all competing approaches use external property predictors.The baselines are MIMOSA, Modof, MolSearch, FRATTVAE, and DyMol.Source: §5.3 · pp. 6–7
5.4 Overall performance and visualization
- TransDLM mostly outperformed baselines across multiple metrics, including a 46.7% increase in the ratio of optimized molecules meeting all ADMET criteria.The reported improvements also include 78% in FCD and 11.8% in Levenshtein distance; the paper reports LogD, Solubility, and Clearance improvements of 57%, 2.4%, and 0.3%, respectively. Validity was 0.993, with extended training slightly improving validity but compromising property optimization and structural similarity.Source: §5.4 · p. 7, Table 1 · p. 7, Table 2 · p. 7
- Representative outputs show TransDLM making targeted substitutions or preserving the source scaffold while meeting ADMET requirements.In the first two examples it makes one-atom or one-functional-group substitutions; competing methods either disrupted structure or missed property requirements.Source: §5.4 · p. 7, Figure 4 · p. 8
5.5 Ablation study
- IUPAC-name descriptions outperform SMILES descriptions, while source-molecule sampling favors structural metrics over sampling from pure noise.Source-molecule sampling gives up some property-criterion advantages; BERT underperforms across metrics, BioBERT and SciBERT are comparable on structural similarity, and SciBERT gives the best property optimization results. The paper relates these differences to their pre-training corpora, including SciBERT’s broader scientific and chemical coverage.Source: §5.5 · p. 9
5.6 Case study analysis
- Optimized XAC had a higher docking score for A1R and a lower score for A2AR than the original, matching the desired preferential binding.Docking scores were evaluated with Schrödinger; lower scores indicate greater binding affinity.Source: §5.6 · p. 9, Figure 5 · p. 9
- The XAC binding-affinity task is out of distribution for a model trained on ADMET data, and its success is presented as evidence of generalization and flexibility.The authors infer that training on broader data including binding-affinity information could extend its use beyond ADMET optimization.Source: §5.6 · p. 9
6 Limitations and future work
TransDLM’s multi-stage pipeline and iterative denoising create complexity and computational costs, while future work targets broader property tasks and improved out-of-distribution robustness.
6 Limitations and future work
- The multi-stage pipeline relies on RECAP decomposition, IUPAC parsing, and pretrained-language-model embeddings, adding external-tool overhead and complicating large-scale deployment.The modular design improves controllability and interpretability, but increases system complexity compared with end-to-end models.Source: §6 · p. 9
- Diffusion-based generation requires iterative denoising and is more expensive than single-step GAN or VAE methods, leaving efficiency a concern.The authors reduced sampling from 2000 to 200 steps with minimal performance loss and propose further acceleration methods.Source: §6 · pp. 9–10
- Future work will extend optimization to QED, SA, and protein-target affinity, and use larger, more diverse libraries to improve robustness on OOD tasks.The authors also aim to streamline the pipeline through deeper integration or end-to-end trainable components.Source: §6 · p. 10
7 Conclusion
TransDLM combines IUPAC-based molecular semantics and detailed textual guidance with source-molecule representations for text-guided optimization. It outperformed other methods on most structural metrics and all ADMET criteria, and successfully optimized XAC.
7 Conclusion
- TransDLM outperformed other methods on most structural metrics and all ADMET property criteria.The conclusion also reports robust structural construction and property enhancement capabilities.Source: §7 · p. 10
- TransDLM uses IUPAC nomenclature as molecular semantic representation and physically and chemically detailed textual descriptions to guide optimization.It samples from encoded source molecules to retain core scaffolds and mitigates error propagation without an external property predictor.Source: §7 · p. 10
- TransDLM successfully optimized XAC, which the authors present as evidence of effectiveness and potential on out-of-distribution problems.The reported achievements did not rely on additional data sources or pretraining.Source: §7 · p. 10
A Details of experimental settings
The experimental settings report TransDLM’s model size, implementation environment, and approximate optimization runtime.
A Details of experimental settings
- TransDLM has around 181 million trainable parameters.The model is implemented with PyTorch 2.0.1 and CUDA 11.7.Source: A Details of experimental settings · p. 12
- Optimizing a molecule takes approximately 1.02 seconds on an AMD EPYC 7763 CPU and a single NVIDIA A6000 GPU.The CPU is specified as 64 cores at 2.450GHz.Source: A Details of experimental settings · p. 12
- Table 5 presents the hyper-parameters used for TransDLM.Source: A Details of experimental settings · p. 12
B Details of RECAP disassembly
RECAP performs synthesis-inspired reverse molecular disassembly into plausible fragments and records the decomposition relationships as a tree.
B Details of RECAP disassembly
- RECAP, built in RDKit, emulates forward laboratory synthesis through reverse conversions and decompositions to produce plausible molecular fragments.Source: B Details of RECAP disassembly · p. 12
- RECAP records disassembly in a tree with the original molecule as root, disassembled molecules as parents, and resulting molecules as children.Source: B Details of RECAP disassembly · p. 12
- Molecules that cannot be further disassembled are designated leaf nodes; related nodes are described as ancestors or descendants.Source: B Details of RECAP disassembly · p. 12
C Statistical Indicator Testing
Repeated-run testing found that most corrected comparisons were statistically significant across five critical metrics. The analysis used paired t-tests with Benjamini-Hochberg correction.
C Statistical Indicator Testing
- Testing covered BLEU, MACCS, FCD, ALL, and HV across five independent runs per method using different random seeds.The metrics cover syntactic similarity, substructure properties, chemical-space similarity, multi-property optimization success, and Pareto-front quality.Source: C Statistical Indicator Testing · p. 12
- Paired t-tests were the primary test, with Benjamini-Hochberg correction used to control false discovery rate across comparisons.The testing assessed whether performance improvements over competing baselines were statistically significant rather than due to randomness.Source: C Statistical Indicator Testing · p. 12
- 24 of 25 corrected p-values were below 0.05, and 22 were below 0.001.Figure 6 marks p-value thresholds at 0.05 and 0.001.Source: Figure 6 · p. 12, C Statistical Indicator Testing · p. 12
D Details of case study analysis
The case study workflow identifies potential binding sites, prepares ligand protonation states, and docks the ligand to rank poses by score.
D Details of case study analysis
- Potential binding sites were identified with SiteMap using a 15×15×15 cubic grid centered on critical receptor functional areas.The configuration followed existing research, and used Schrödinger’s SiteMap tool.Source: D Details of case study analysis · p. 12
- LigPrep generated XAC protonation states at pH 7.0 ± 2.0 using Epik, with up to 32 states per ligand.The passage describes this as preprocessing before docking.Source: D Details of case study analysis · p. 12
- Glide docking in Maestro automatically scored each pose to identify the most favorable binding conformations.Docking followed ligand preparation and binding-site localization.Source: D Details of case study analysis · p. 12
Related papers
This paper cites
- Attention Is All You Need2017
- BERT: Pre-training of Deep Bidirectional Transformers for Language Understanding2018
- PyTorch: An Imperative Style, High-Performance Deep Learning Library2019
- SciBERT: A Pretrained Language Model for Scientific Text2019
- Junction Tree Variational Autoencoder for Molecular Graph Generation2018
- AudioLDM: Text-to-Audio Generation with Latent Diffusion Models2023