Drug Screening
COMPUTATIONALThis pipeline computationally filters thousands of ChEMBL compounds down to the best candidates for SMA drug discovery. The process runs in six steps: (1) ChEMBL query — compounds bioactive against top-scored SMA targets are fetched with their SMILES strings; (2) RDKit descriptor calculation — molecular weight, LogP, rotatable bonds, H-bond donors/acceptors, TPSA, and QED are computed from SMILES; (3) Lipinski Rule of 5— MW < 500, LogP < 5, HBD ≤ 5, HBA ≤ 10; (4) BBB permeability estimate— TPSA < 90 Ų and MW < 450; (5) CNS MPO score — 0–6 composite tuned for CNS drug development; (6) PAINS filter — substructure alerts for pan-assay interference compounds.
Why BBB penetration matters for SMA:SMA motor neurons reside in the anterior horn behind the blood-brain barrier. Risdiplam succeeds partly because of its BBB-permeable profile. Compounds with TPSA > 90 Ų or MW > 500 Da are unlikely to achieve meaningful CNS exposure.
Score glossary: Lipinski — binary pass/fail. BBB — heuristic estimate (TPSA + MW). CNS MPO — 0–6; ≥ 4 is CNS-optimized. QED — 0–1 drug-likeness; ≥ 0.5 is high quality. PAINS — substructure alert for reactive scaffolds.
▶How does Computational Drug Screening work?
The screening library starts from ChEMBL (25,000+ compounds with known bioactivity against SMA-relevant targets). Compounds pass through sequential filters:
- Lipinski Ro5 — MW ≤500, LogP ≤5, HBD ≤5, HBA ≤10
- BBB heuristic— TPSA < 90 Ų, MW < 450
- CNS MPO ≥ 4 — composite CNS optimization score
- QED ≥ 0.5 — drug-likeness estimate
- PAINS-free — no pan-assay interference alerts
* QED, CNS MPO, and PAINS columns are estimated values (heuristic from pChEMBL, MW/LogP thresholds). Not RDKit-computed. Use /screen/smiles for exact values.
Note: Drug-likeness predictions use rule-based heuristics (Lipinski Rule of 5, TPSA-based BBB estimate, QED score). These are filtering tools, not validated PK/tox models.
Top Candidates
| ChEMBL ID | Target | MW | LogP | QED* | CNS MPO* | BBB* | Lipinski | PAINS* | pChEMBL |
|---|---|---|---|---|---|---|---|---|---|
| CHEMBL305177 | MTOR_PATHWAY | 438.4 | 1.45 | 1.10 | 4 | YES | PASS | CLEAN | 11.0 |
| CHEMBL305177 | NCALD | 438.4 | 1.45 | 1.10 | 4 | YES | PASS | CLEAN | 11.0 |
| CHEMBL305177 | NEDD4L | 438.4 | 1.45 | 1.10 | 4 | YES | PASS | CLEAN | 11.0 |
| CHEMBL305177 | NMJ_MATURATION | 438.4 | 1.45 | 1.10 | 4 | YES | PASS | CLEAN | 11.0 |
| CHEMBL305177 | PLS3 | 438.4 | 1.45 | 1.10 | 4 | YES | PASS | CLEAN | 11.0 |
| CHEMBL305177 | SPATA18 | 438.4 | 1.45 | 1.10 | 4 | YES | PASS | CLEAN | 11.0 |
| CHEMBL305177 | SULF1 | 438.4 | 1.45 | 1.10 | 4 | YES | PASS | CLEAN | 11.0 |
| CHEMBL305177 | LIMK2 | 438.4 | 1.45 | 1.10 | 4 | YES | PASS | CLEAN | 11.0 |
| CHEMBL305177 | ROCK2 | 438.4 | 1.45 | 1.10 | 4 | YES | PASS | CLEAN | 11.0 |
| CHEMBL305177 | ANK3 | 438.4 | 1.45 | 1.10 | 4 | YES | PASS | CLEAN | 11.0 |
| CHEMBL305177 | CAST | 438.4 | 1.45 | 1.10 | 4 | YES | PASS | CLEAN | 11.0 |
| CHEMBL410215 | MTOR_PATHWAY | 496.0 | 3.97 | 1.05 | 2 | NO | PASS | CLEAN | 10.5 |
| CHEMBL410215 | NCALD | 496.0 | 3.97 | 1.05 | 2 | NO | PASS | CLEAN | 10.5 |
| CHEMBL410215 | NEDD4L | 496.0 | 3.97 | 1.05 | 2 | NO | PASS | CLEAN | 10.5 |
| CHEMBL410215 | NMJ_MATURATION | 496.0 | 3.97 | 1.05 | 2 | NO | PASS | CLEAN | 10.5 |
| CHEMBL410215 | PLS3 | 496.0 | 3.97 | 1.05 | 2 | NO | PASS | CLEAN | 10.5 |
| CHEMBL410215 | SPATA18 | 496.0 | 3.97 | 1.05 | 2 | NO | PASS | CLEAN | 10.5 |
| CHEMBL410215 | SULF1 | 496.0 | 3.97 | 1.05 | 2 | NO | PASS | CLEAN | 10.5 |
| CHEMBL410215 | LIMK2 | 496.0 | 3.97 | 1.05 | 2 | NO | PASS | CLEAN | 10.5 |
| CHEMBL410215 | ROCK2 | 496.0 | 3.97 | 1.05 | 2 | NO | PASS | CLEAN | 10.5 |
| CHEMBL410215 | ANK3 | 496.0 | 3.97 | 1.05 | 2 | NO | PASS | CLEAN | 10.5 |
| CHEMBL410215 | CAST | 496.0 | 3.97 | 1.05 | 2 | NO | PASS | CLEAN | 10.5 |
| CHEMBL38735 | MTOR_PATHWAY | 381.9 | 4.44 | 1.03 | 2 | NO | PASS | CLEAN | 10.3 |
| CHEMBL38735 | NCALD | 381.9 | 4.44 | 1.03 | 2 | NO | PASS | CLEAN | 10.3 |
| CHEMBL38735 | NEDD4L | 381.9 | 4.44 | 1.03 | 2 | NO | PASS | CLEAN | 10.3 |