PMODiff:Physics-Informed Multi-Objective Optimization Diffusion Model for Protein-Specific 3D Molecule Generation
The code has been tested in the following environment:
conda create -n tagmol python=3.8.17
conda activate tagmol
conda install pytorch=1.13.1 pytorch-cuda=11.6 -c pytorch -c nvidia
conda install pyg=2.2.0 -c pyg
conda install rdkit=2022.03.2 openbabel=3.1.1 tensorboard=2.13.0 pyyaml=6.0 easydict=1.9 python-lmdb=1.4.1 -c conda-forge
# For Vina Docking
pip install meeko==0.1.dev3 scipy pdb2pqr vina==1.2.2
python -m pip install git+https://github.com/Valdes-Tresanco-MS/AutoDockTools_py3The resources can be found here. The data are inside data directory, and the guide checkpoints are inside logs.
python scripts/train_diffusion.py configs/training.ymlpython scripts/train_dock_guide.py configs/training_dock_guide.ymlpython scripts/train_dock_guide.py configs/training_dock_guide_qed.ymlpython scripts/train_dock_guide.py configs/training_dock_guide_sa.ymlNOTE: The outputs are saved in logs/ by default.
python scripts/sample_diffusion.py configs/sampling.yml --data_id {i} # Replace {i} with the index of the data. i should be between 0 and 99 for the testset.We have a bash file that can run the inference for the entire test set in a loop.
bash scripts/batch_sample_diffusion.sh configs/sampling.yml backbonepython scripts/sample_multi_guided_diffusion.py [path-to-config.yml] --data_id {i} # Replace {i} with the index of the data. i should be between 0 and 99 for the testset.To run inference on all 100 targets in the test set:
bash scripts/batch_sample_multi_guided_diffusion.sh [path-to-config.yml] [output-dir-name]The outputs are stored in experiments_multi/[output-dir-name]when run using the bash file. The config files are available in configs/noise_guide_multi.
- Single-objective guidance
- BA:
sampling_guided_ba_1.yml - QED:
sampling_guided_qed_1.yml - SA:
sampling_guided_sa_1.yml
- BA:
- Multi-objective guidance (our main model)
- QED + SA + BA:
sampling_guided_qed_0.33_sa_0.33_ba_0.34.yml
- QED + SA + BA:
For example, to run the multi-objective setting (i.e., our model):
bash scripts/batch_sample_multi_guided_diffusion.sh configs/noise_guide_multi/sampling_guided_qed_0.33_sa_0.33_ba_0.34.yml qed_0.33_sa_0.33_ba_0.34python scripts/eval_dock_guide.py --ckpt_path [path-to-checkpoint.pt]python scripts/evaluate_diffusion.py {OUTPUT_DIR} --docking_mode vina_score --protein_root data/test_setThe docking mode can be chosen from {qvina, vina_score, vina_dock, none}
This codebase was build on top of TargetDiff