Adme Property Predictor

ADME Property Predictor

Overview

Comprehensive pharmacokinetic prediction tool that assesses drug-likeness and ADME properties of small molecules using validated cheminformatics models, molecular descriptors, and structure-property relationships.

Key Capabilities:

• Multi-Property Prediction: Absorption, Distribution, Metabolism, Excretion
• Drug-Likeness Scoring: Lipinski's Rule of 5, Veber rules, QED score
• Batch Processing: Analyze compound libraries efficiently
• Structure-Based Insights: Identify liability hotspots and optimization opportunities
• Comparative Analysis: Rank candidates by predicted PK profile

When to Use

✅ Use this skill when:

• Screening compound libraries for drug-like properties in early discovery
• Prioritizing lead compounds for advancement based on predicted PK
• Identifying ADME liabilities requiring structural optimization
• Comparing analogs to select candidates with optimal ADME profiles
• Filtering virtual screening hits before synthesis
• Generating ADME data for regulatory pre-submission packages
• Teaching pharmacokinetics and drug design principles

❌ Do NOT use when:

• Exact PK parameters needed for dosing → Use experimental PK studies
• Biologics (antibodies, proteins) → Use antibody-pk-predictor
• Natural products with complex structures → Models trained on synthetic small molecules
• Prodrugs requiring metabolic activation → Use prodrug-activation-predictor
• Prediction for clinical dosing decisions → CRITICAL: Experimental validation required
• Assessing toxicity or safety → Use toxicity-structure-alert or admetox-predictor

Related Skills:

• 上游: chemical-structure-converter (structure preparation), lipinski-rule-filter (rule-based filtering)
• 下游: drug-candidate-evaluator (integrated scoring), molecular-dynamics-sim (detailed binding)

Integration with Other Skills

Upstream Skills:

• chemical-structure-converter: Convert between SMILES, InChI, MOL formats
• lipinski-rule-filter: Initial rule-based drug-likeness screening
• chemical-structure-converter: Generate 3D conformers for structure-based predictions
• smiles-de-salter: Remove salt counterions before analysis

Downstream Skills:

• drug-candidate-evaluator: Multi-parameter optimization including ADME
• toxicity-structure-alert: Assess safety alongside ADME
• target-novelty-scorer: Evaluate target uniqueness for selected candidates
• biotech-pitch-deck-narrative: Create investor materials with PK data

Complete Workflow:

Chemical Structure Converter (prepare structures) → 
  Lipinski Rule Filter (initial filtering) → 
    ADME Property Predictor (this skill, detailed PK) → 
      Drug Candidate Evaluator (integrated scoring) → 
        Toxicity Structure Alert (safety check)

Core Capabilities

1. Absorption (A) Prediction

Predict intestinal absorption, solubility, and permeability:

from scripts.adme_predictor import ADMEPredictor

predictor = ADMEPredictor()

# Predict absorption properties
absorption = predictor.predict_absorption(
    smiles="CC(=O)Oc1ccccc1C(=O)O",  # Aspirin
    properties=["all"]  # or specific: ["hia", "caco2", "solubility"]
)

print(absorption.summary())

Predicted Properties:

Best Practices:

• ✅ Consider HIA and solubility together (high HIA but low solubility = dissolution-limited)
• ✅ Caco-2 good for oral absorption prediction; poor for BBB penetration
• ✅ Use both rule-based (Lipinski) and ML-based predictions for consensus
• ✅ Check solubility at physiological pH (not just intrinsic)

Common Issues and Solutions:

Issue: Lipinski pass but poor solubility

• Symptom: "Passes Rule of 5 but LogS = -5"
• Solution: Lipinski checks MW and LogP, not solubility directly; use explicit solubility prediction

Issue: Caco-2 predicts high absorption but HIA low

• Symptom: "Caco-2 = 85 (high) but HIA = 60%"
• Solution: Models have different training sets; Caco-2 is in vitro, HIA in vivo; HIA generally more reliable

2. Distribution (D) Prediction

Predict tissue distribution, protein binding, and brain penetration:

# Predict distribution properties
distribution = predictor.predict_distribution(
    smiles="CC(=O)Oc1ccccc1C(=O)O",
    properties=["vd", "ppb", "bbb"]
)

# Access specific predictions
vd = distribution.volume_of_distribution
bbb = distribution.blood_brain_barrier
ppb = distribution.plasma_protein_binding

Predicted Properties:

Best Practices:

• ✅ High PPB (>90%) may require higher doses but longer half-life
• ✅ Low Vd (<0.3) = mainly in plasma; high Vd (>3) = extensive tissue distribution
• ✅ BBB penetration critical for CNS drugs; avoid for peripherally-acting drugs
• ✅ fu (free fraction) drives pharmacological activity, not total concentration

Common Issues and Solutions:

Issue: BBB predictions unreliable for certain chemotypes

• Symptom: "BBB model gives conflicting predictions for peptides"
• Solution: Models trained on small molecules; use specialized BBB predictors for peptides, macrocycles

Issue: PPB overestimated for acidic drugs

• Symptom: "PPB predicted 95% but experimental is 70%"
• Solution: Some models biased toward neutral/basic compounds; check model training set overlap

3. Metabolism (M) Prediction

Predict metabolic stability, CYP interactions, and liability sites:

# Predict metabolism properties
metabolism = predictor.predict_metabolism(
    smiles="CC(=O)Oc1ccccc1C(=O)O",
    include_site_prediction=True
)

# Check CYP interactions
cyp_profile = metabolism.cyp_profile
stability = metabolism.metabolic_stability

Predicted Properties:

Best Practices:

• ✅ Screen for CYP3A4 inhibition early (most common DDI)
• ✅ Check if compound is CYP substrate (for polymorphism concerns)
• ✅ Identify metabolic hotspots for structural blocking
• ✅ Consider species differences (human vs rodent metabolism)

Common Issues and Solutions:

Issue: False negatives for time-dependent inhibition (TDI)

• Symptom: "No CYP inhibition predicted but TDI observed experimentally"
• Solution: Standard models predict reversible inhibition; use specialized TDI predictors

Issue: Metabolic site prediction shows multiple hotspots

• Symptom: "5 different atoms flagged as metabolic liabilities"
• Solution: Prioritize by reactivity score; consider blocking highest-risk site first

4. Excretion (E) Prediction

Predict clearance routes and elimination kinetics:

# Predict excretion properties
excretion = predictor.predict_excretion(
    smiles="CC(=O)Oc1ccccc1C(=O)O",
    properties=["clearance", "half_life", "route"]
)

# Access predictions
clearance = excretion.clearance_ml_min_kg
t12 = excretion.half_life_hours
route = excretion.primary_route

Predicted Properties:

Best Practices:

• ✅ Half-life determines dosing frequency (T1/2 × 5 = time to steady state)
• ✅ Renal clearance predictable for polar compounds; hepatic less predictable
• ✅ High clearance (>15) may require high doses or prodrug approach
• ✅ Very long T1/2 (>24h) good for adherence but risk accumulation

Common Issues and Solutions:

Issue: Clearance predictions highly variable

• Symptom: "Same compound, different models give CL = 5 vs 20 mL/min/kg"
• Solution: Allometry-based methods unreliable for novel scaffolds; use average of multiple models

Issue: Route prediction contradicts structure

• Symptom: "Highly polar compound predicted biliary, expected renal"
• Solution: Check LogP/LogD; polar compounds (<0) usually renal; neutral/lipophilic (>1) usually hepatic

5. Integrated Drug-Likeness Scoring

Overall assessment combining all ADME properties:

# Generate comprehensive drug-likeness score
druglikeness = predictor.calculate_druglikeness(
    smiles="CC(=O)Oc1ccccc1C(=O)O",
    methods=["qed", "muegge", "golden_triangle"]
)

# Multi-parameter optimization
mpo_score = predictor.mpo_score(
    smiles="CC(=O)Oc1ccccc1C(=O)O",
    target_profile={"hia": >80, "bbb": <0.3, "t12": "2-8h"}
)

Scoring Methods:

Best Practices:

• ✅ Use QED as quick overall metric; MPO for property-weighted scoring
• ✅ Don't rely solely on drug-likeness; efficacy and safety equally important
• ✅ Compare to marketed drugs in same class for context
• ✅ Track drug-likeness trends during optimization (should improve)

Common Issues and Solutions:

Issue: Drug-likeness score conflicts with project needs

• Symptom: "CNS drug has low QED (0.5) because high LogP needed for BBB"
• Solution: Drug-likeness rules biased toward oral drugs; use category-specific models (CNS, oncology, etc.)

6. Batch Processing and Library Screening

Analyze compound libraries efficiently:

# Batch process library
results = predictor.batch_predict(
    input_file="library.smi",  # SMILES file
    properties=["all"],
    output_format="csv",
    n_workers=4  # Parallel processing
)

# Filter by criteria
filtered = results.filter(
    lipinski_pass=True,
    hia__gt=80,
    t12__between=(2, 8)
)

# Rank by multi-parameter score
ranked = results.rank(by="mpo_score", ascending=False)

Best Practices:

• ✅ Process in batches of 1000-10000 for memory efficiency
• ✅ Save intermediate results (crash recovery)
• ✅ Apply filters sequentially (Lipinski first, then detailed ADME)
• ✅ Check property distributions to identify outliers

Common Issues and Solutions:

Issue: Batch processing runs out of memory

• Symptom: "Killed: Out of memory" with 50K compounds
• Solution: Process in chunks; use generators instead of loading all into RAM

Issue: Some compounds fail prediction

• Symptom: "30% of library returns NaN"
• Solution: Check for invalid SMILES, unusual atoms, or molecules outside training set domain

Complete Workflow Example

From SMILES to prioritized candidates:

# Step 1: Predict ADME for single compound
python scripts/main.py \
  --smiles "CC(=O)Oc1ccccc1C(=O)O" \
  --properties all \
  --output aspirin_adme.json

# Step 2: Batch process compound library
python scripts/main.py \
  --input library.smi \
  --properties absorption,distribution \
  --format csv \
  --output library_adme.csv

# Step 3: Filter and rank
python scripts/main.py \
  --input library_adme.csv \
  --filter "lipinski_pass=True,hia>80" \
  --rank-by qed \
  --top-n 100 \
  --output top_candidates.csv

Python API Usage:

from scripts.adme_predictor import ADMEPredictor
from scripts.batch_processor import BatchProcessor

# Initialize
predictor = ADMEPredictor()
batch = BatchProcessor()

# Single compound analysis
aspirin = predictor.predict_all("CC(=O)Oc1ccccc1C(=O)O")
print(f"HIA: {aspirin.absorption.hia}%")
print(f"Half-life: {aspirin.excretion.t12} hours")

# Batch screening
results = batch.process(
    input_file="library.smi",
    predictor=predictor,
    properties=["absorption", "distribution"],
    n_workers=4
)

# Filter good candidates
good_candidates = results[
    (results.lipinski_pass == True) &
    (results.hia > 80) &
    (results.bbb < 0.3) &
    (results.t12.between(2, 8))
]

Expected Output Files:

output/
├── aspirin_adme.json           # Single compound detailed results
├── library_adme.csv            # Batch screening results
├── top_candidates.csv          # Filtered and ranked candidates

Quality Checklist

Pre-Prediction Checks:

• SMILES string is valid and canonical
• Salt forms removed (if analyzing parent compound)
• Tautomeric state appropriate for physiological pH
• Stereochemistry specified (if relevant for activity)

During Prediction:

• Compound within model applicability domain (check similarity to training set)
• No unusual atoms or functional groups (models trained on typical drug-like space)
• MW in range 100-800 Da (outside range predictions less reliable)
• Predictions complete (no missing values for critical properties)

Post-Prediction Verification:

• Drug-likeness scores in reasonable range (sanity check)
• Individual properties internally consistent (e.g., high LogP predicts low solubility)
• CRITICAL: Comparison to experimental data if available (validate model for chemotype)
• Rankings align with medicinal chemistry intuition

Before Making Decisions:

• CRITICAL: Predictions are NOT experimental data; use for prioritization only
• Multiple orthogonal models give consistent results
• Structural alerts checked (toxicity, reactivity)
• Top candidates selected for experimental validation
• Documentation of model versions and confidence intervals

For Regulatory Submissions:

• Model validation documented (training set, test set performance)
• Applicability domain clearly defined
• Prediction uncertainty quantified
• Experimental confirmation for key predictions