eToxPred

What is eToxPred?

eToxPred is a machine learning tool for predicting the toxicity and synthetic accessibility of small molecules from their chemical structures. Developed at Louisiana State University by Limeng Pu, Michal Brylinski, and colleagues, eToxPred filters out potentially toxic or difficult-to-synthesize compounds early in the drug discovery process.

The tool provides two complementary scores:

• Tox-score indicating the probability that a compound is toxic (0–1)
• SA score measuring synthetic accessibility after transformation to a 0–1 scale, where higher values indicate easier synthesis.

Together, these scores help prioritize which drug candidates are worth pursuing.

As such, we recommend using eToxPred for the following screening applications:

• Virtual screening: Filtering large compound libraries to remove high-risk molecules before expensive docking simulations or experimental testing
• Lead optimization: Evaluating whether structural modifications improve or worsen the toxicity and synthetic accessibility profile
• Hit prioritization: Ranking compounds from high-throughput screens by their likelihood of progressing through development
• Library design: Guiding the selection of compounds for purchase or synthesis based on favorable predicted properties

How to use eToxPred online

ProteinIQ provides a web-based interface for running eToxPred without command-line installation or Python environment configuration. Enter SMILES strings and receive toxicity and synthetic accessibility predictions.

Inputs

Results

The output is a spreadsheet with toxicity and synthetic accessibility predictions for each compound.

Interpreting toxicity scores

The Tox-score represents the probability that a compound exhibits general toxicity based on structural similarity to known toxic and non-toxic compounds.

The optimal discrimination threshold is 0.58, which most effectively separates toxic from non-toxic compounds in validation studies. FDA-approved drugs have a median Tox-score of approximately 0.34, while known toxins from the T3DB database typically score above 0.6.

Interpreting SA scores

The SA score is eToxPred’s transformed synthetic accessibility estimate. Higher values indicate compounds predicted to be easier to synthesize; lower values indicate greater synthetic difficulty. The source does not define universal decision bands for this transformed value, so use it primarily to rank compounds analyzed with the same model and confirm promising candidates with synthesis planning or expert review.

How does eToxPred work?

eToxPred combines toxicity prediction from an Extremely Randomized Trees (Extra Trees) classifier with synthetic accessibility scoring based on RDKit molecular fragments and complexity penalties.

Toxicity prediction

The toxicity model was trained on 4,550 compounds: 1,515 FDA-approved drugs representing the non-toxic class and 3,035 compounds from TOXNET representing the toxic class. Independent validation used 3,682 compounds from KEGG-Drug (non-toxic) and 1,283 compounds from T3DB (toxic).

Molecular representation

Each molecule is converted to a 1024-bit Morgan fingerprint with radius 2 using RDKit. These binary fingerprints encode the presence or absence of structural fragments, capturing the chemical features relevant to toxicity.

Extra Trees classifier

The current trained Extra Trees model uses the published eToxPred model archive. Key hyperparameters reported with the model:

• Number of estimators: 550
• Minimum samples per split: 16
• Minimum samples per leaf: 3
• Features per split: 10

The ensemble votes on classification, with the final Tox-score representing the proportion of trees predicting toxicity. This approach handles noisy biological data well and resists overfitting.

Synthetic accessibility scoring

The SA score combines historical synthetic knowledge with complexity penalties. The current eToxPred implementation calculates an Ertl-Schuffenhauer-style raw score and transforms it with:

$\text{eToxPred SA score} = \exp(1 - \text{raw SA score})$

This transformation means higher displayed values correspond to easier predicted synthesis.

The fragment score compares molecular substructures against fragments frequently found in known synthesized compounds. Common fragments score higher (easier to make); unusual fragments score lower.

The complexity penalty accounts for structural features that complicate synthesis:

• Spiro and fused ring systems
• Multiple stereocenters
• Macrocyclic structures
• Non-standard bridging patterns

Performance metrics

eToxPred was validated on independent test sets not used during training.

General toxicity (KEGG-Drug/T3DB test set)

Specific toxicity endpoints

The model was also evaluated on datasets for specific toxicity types:

Synthetic accessibility

The current SA score calculation follows the Ertl-Schuffenhauer synthetic accessibility approach used by the eToxPred package, then transforms the score so higher 0–1 values indicate easier synthesis.

Comparison with related tools

Several tools on ProteinIQ address overlapping aspects of compound evaluation:

eToxPred provides general toxicity screening with synthetic accessibility in a single analysis. The machine learning model captures patterns across diverse toxic compounds but does not distinguish between specific toxicity mechanisms.

ADMET-AI uses graph neural networks to predict 41 specific ADMET endpoints, including hERG inhibition (cardiotoxicity), hepatotoxicity, CYP interactions, and plasma protein binding. For endpoint-specific toxicity predictions, ADMET-AI offers more detailed information.

Toxicity Prediction uses rule-based structural alerts (PAINS, Brenk filters) rather than machine learning. This approach identifies specific problematic substructures like reactive groups or known interference patterns. The two approaches are complementary—eToxPred captures general toxicity patterns while structural alerts identify specific problematic features.

Lipinski's Rule of 5 evaluates oral bioavailability potential using simple physicochemical rules (molecular weight, LogP, hydrogen bond donors/acceptors). This rule-based approach is interpretable but does not predict toxicity.

QEPPi scores drug-likeness specifically for protein-protein interaction inhibitors, which require different physicochemical properties than conventional drugs.

Example workflow

A typical drug discovery screening workflow incorporating eToxPred:

1. Generate compound library: Start with a virtual library or vendor catalog
2. Filter with Lipinski's rules: Remove compounds unlikely to be orally bioavailable using Lipinski's Rule of 5
3. Structural alert screening: Flag compounds with problematic substructures using Toxicity Prediction
4. Toxicity screening with eToxPred: Remove compounds with Tox-score > 0.6
5. Synthetic accessibility check: Rank compounds by higher SA score, then review synthesis feasibility for finalists
6. Detailed ADMET profiling: Run remaining candidates through ADMET-AI
7. Docking and binding: Proceed to structure-based screening with AutoDock Vina

Example compounds

Try aspirin (acetylsalicylic acid), CC(=O)Oc1ccccc1C(=O)O, and caffeine, Cn1cnc2c1c(=O)n(c(=O)n2C)C, in one job to compare their predictions. Treat the scores as model outputs for prioritization rather than experimental safety or synthesis guarantees.