Update README.md
Browse files
README.md
CHANGED
|
@@ -1,70 +1,77 @@
|
|
| 1 |
-
|
| 2 |
-
|
| 3 |
-
|
| 4 |
-
|
| 5 |
-
|
| 6 |
-
|
| 7 |
-
|
| 8 |
-
|
| 9 |
-
|
| 10 |
-
|
| 11 |
-
|
| 12 |
-
|
| 13 |
-
|
| 14 |
-
|
| 15 |
-
|
| 16 |
-
|
| 17 |
-
|
| 18 |
-
|
| 19 |
-
|
| 20 |
-
|
| 21 |
-
|
| 22 |
-
|
| 23 |
-
|
| 24 |
-
|
| 25 |
-
|
| 26 |
-
|
| 27 |
-
|
| 28 |
-
|
| 29 |
-
|
| 30 |
-
|
| 31 |
-
|
| 32 |
-
|
| 33 |
-
|
| 34 |
-
|
| 35 |
-
|
| 36 |
-
|
| 37 |
-
|
| 38 |
-
|
| 39 |
-
|
| 40 |
-
|
| 41 |
-
|
| 42 |
-
|
| 43 |
-
|
| 44 |
-
|
| 45 |
-
|
| 46 |
-
|
| 47 |
-
|
| 48 |
-
|
| 49 |
-
|
| 50 |
-
|
| 51 |
-
|
| 52 |
-
|
| 53 |
-
|
| 54 |
-
|
| 55 |
-
)
|
| 56 |
-
|
| 57 |
-
|
| 58 |
-
|
| 59 |
-
|
| 60 |
-
|
| 61 |
-
|
| 62 |
-
|
| 63 |
-
|
| 64 |
-
|
| 65 |
-
|
| 66 |
-
|
| 67 |
-
|
| 68 |
-
|
| 69 |
-
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| 70 |
```
|
|
|
|
| 1 |
+
---
|
| 2 |
+
license: mit
|
| 3 |
+
language:
|
| 4 |
+
- en
|
| 5 |
+
tags:
|
| 6 |
+
- PyTorch
|
| 7 |
+
---
|
| 8 |
+
|
| 9 |
+
# De Novo Drug Generator - RNN-VAE
|
| 10 |
+
|
| 11 |
+
De Novo Drug Generator - RNN-VAE is a deep learning model designed for generating novel drug molecules.
|
| 12 |
+
Training data from ChemBL library
|
| 13 |
+
|
| 14 |
+
Full project file at https://github.com/teohyc/drug_agent
|
| 15 |
+
|
| 16 |
+
|
| 17 |
+
## Usage
|
| 18 |
+
|
| 19 |
+
```python
|
| 20 |
+
from rdkit import Chem
|
| 21 |
+
from rdkit.Chem import Draw, Descriptors
|
| 22 |
+
from tree_rnn_vae_infer import generate_candidate_mol
|
| 23 |
+
from tree_rnn_vae_model import TreeEncoder, LatentHead, TreeVAE, TreeDecoder
|
| 24 |
+
|
| 25 |
+
|
| 26 |
+
def compute_molecule_props(mol):
|
| 27 |
+
return {
|
| 28 |
+
"MW": Descriptors.MolWt(mol),
|
| 29 |
+
"logP": Descriptors.MolLogP(mol),
|
| 30 |
+
"HBD": Descriptors.NumHDonors(mol),
|
| 31 |
+
"HBA": Descriptors.NumHAcceptors(mol),
|
| 32 |
+
}
|
| 33 |
+
|
| 34 |
+
|
| 35 |
+
# display molecule
|
| 36 |
+
def render_molecule_grid(selected):
|
| 37 |
+
if not selected:
|
| 38 |
+
return
|
| 39 |
+
|
| 40 |
+
mols, legends = [], []
|
| 41 |
+
|
| 42 |
+
if isinstance(selected, dict):
|
| 43 |
+
iterable = selected.items()
|
| 44 |
+
else:
|
| 45 |
+
iterable = enumerate(selected, 1)
|
| 46 |
+
|
| 47 |
+
for i, item in iterable:
|
| 48 |
+
if isinstance(selected, dict):
|
| 49 |
+
smi, props = i, item
|
| 50 |
+
else:
|
| 51 |
+
smi, props = item, None
|
| 52 |
+
|
| 53 |
+
mol = Chem.MolFromSmiles(smi)
|
| 54 |
+
if mol:
|
| 55 |
+
mols.append(mol)
|
| 56 |
+
if props is None:
|
| 57 |
+
props = compute_molecule_props(mol)
|
| 58 |
+
legends.append(
|
| 59 |
+
f"M{i}
|
| 60 |
+
MW={props['MW']:.0f}, logP={props['logP']:.2f}, "
|
| 61 |
+
f"HBD={props['HBD']}, HBA={props['HBA']}"
|
| 62 |
+
)
|
| 63 |
+
|
| 64 |
+
img = Draw.MolsToGridImage(
|
| 65 |
+
mols,
|
| 66 |
+
molsPerRow=3,
|
| 67 |
+
subImgSize=(400, 400),
|
| 68 |
+
legends=legends,
|
| 69 |
+
useSVG=False,
|
| 70 |
+
)
|
| 71 |
+
return img
|
| 72 |
+
|
| 73 |
+
|
| 74 |
+
molecules = generate_candidate_mol(num_samples=6, max_len=20) #change to your desired molecule size and number
|
| 75 |
+
img = render_molecule_grid(molecules)
|
| 76 |
+
img.show()
|
| 77 |
```
|