Chemical Safety & Toxicology Assessment
Comprehensive chemical safety and toxicology analysis integrating predictive AI models, curated toxicogenomics databases, regulatory safety data, and chemical-biological interaction networks. Generates structured risk assessment reports with evidence grading.
When to Use This Skill
Triggers:
- "Is this chemical toxic?" / "What are the toxicity endpoints for [compound]?"
- "Assess the safety profile of [drug/chemical]"
- "What are the ADMET properties of [SMILES]?"
- "What genes does [chemical] interact with?"
- "What diseases are linked to [chemical] exposure?"
- "Predict toxicity for these molecules"
- "Drug safety assessment for [drug name]"
- "Environmental health risk of [chemical]"
- "Chemical hazard profiling"
- "Toxicogenomic analysis of [compound]"
Use Cases:
- Predictive Toxicology: AI-predicted toxicity endpoints (AMES mutagenicity, DILI, LD50, carcinogenicity, skin reactions) for novel compounds via SMILES
- ADMET Profiling: Full absorption, distribution, metabolism, excretion, toxicity characterization
- Toxicogenomics: Chemical-gene interaction mapping, gene-disease associations from CTD
- Regulatory Safety: FDA label warnings, boxed warnings, contraindications, adverse reactions
- Drug Safety Assessment: Combined DrugBank safety + FDA labels + adverse event data
- Chemical-Protein Interactions: STITCH-based chemical-protein binding and interaction networks
- Environmental Toxicology: Chemical-disease associations for environmental contaminants
KEY PRINCIPLES
- Report-first approach - Create report file FIRST, then populate progressively
- Tool parameter verification - Verify params via
get_tool_info before calling unfamiliar tools
- Evidence grading - Grade all safety claims by evidence strength (T1-T4)
- Citation requirements - Every toxicity finding must have inline source attribution
- Mandatory completeness - All sections must exist with data minimums or explicit "No data" notes
- Disambiguation first - Resolve compound identity (name -> SMILES, CID, ChEMBL ID) before analysis
- Negative results documented - "No toxicity signals found" is data; empty sections are failures
- Conservative risk assessment - When evidence is ambiguous, flag as "requires further investigation"
- English-first queries - Always use English chemical/drug names in tool calls
Evidence Grading System (MANDATORY)
Grade every toxicity claim by evidence strength:
| Tier |
Symbol |
Criteria |
Examples |
| T1 |
[T1] |
Direct human evidence, regulatory finding |
FDA boxed warning, clinical trial toxicity, human case reports |
| T2 |
[T2] |
Animal studies, validated in vitro |
Nonclinical toxicology, AMES positive, animal LD50 |
| T3 |
[T3] |
Computational prediction, association data |
ADMET-AI prediction, CTD association, QSAR model |
| T4 |
[T4] |
Database annotation, text-mined |
Literature mention, database entry without validation |
Required Evidence Grading Locations
Evidence grades MUST appear in:
- Executive Summary - Key toxicity findings graded
- Toxicity Predictions - Every ADMET-AI endpoint with confidence note
- Regulatory Safety - FDA findings marked [T1]
- Chemical-Gene Interactions - CTD data marked by curation status
- Risk Assessment - Final risk classification with supporting evidence tiers
Core Strategy: 8 Research Dimensions
Chemical/Drug Query
|
+-- PHASE 0: Compound Disambiguation (ALWAYS FIRST)
| +-- Resolve name -> SMILES, PubChem CID, ChEMBL ID
| +-- Get molecular formula, weight, canonical structure
|
+-- PHASE 1: Predictive Toxicology (ADMET-AI)
| +-- Mutagenicity (AMES)
| +-- Hepatotoxicity (DILI, ClinTox)
| +-- Carcinogenicity
| +-- Acute toxicity (LD50)
| +-- Skin reactions
| +-- Stress response pathways
| +-- Nuclear receptor activity
|
+-- PHASE 2: ADMET Properties
| +-- Absorption: BBB penetrance, bioavailability
| +-- Distribution: clearance, volume of distribution
| +-- Metabolism: CYP interactions (1A2, 2C9, 2C19, 2D6, 3A4)
| +-- Physicochemical: solubility, lipophilicity, pKa
|
+-- PHASE 3: Toxicogenomics (CTD)
| +-- Chemical-gene interactions
| +-- Chemical-disease associations
| +-- Affected biological pathways
|
+-- PHASE 4: Regulatory Safety (FDA Labels)
| +-- Boxed warnings (Black Box)
| +-- Contraindications
| +-- Adverse reactions
| +-- Warnings and precautions
| +-- Nonclinical toxicology
|
+-- PHASE 5: Drug Safety Profile (DrugBank)
| +-- Toxicity data
| +-- Contraindications
| +-- Drug interactions affecting safety
|
+-- PHASE 6: Chemical-Protein Interactions (STITCH)
| +-- Direct chemical-protein binding
| +-- Interaction confidence scores
| +-- Off-target effects
|
+-- PHASE 7: Structural Alerts (ChEMBL)
| +-- Known toxic substructures (PAINS, Brenk)
| +-- Structural alert flags
|
+-- SYNTHESIS: Integrated Risk Assessment
+-- Aggregate all evidence tiers
+-- Risk classification (Low/Medium/High/Critical)
+-- Data gaps and recommendations
Phase 0: Compound Disambiguation (ALWAYS FIRST)
CRITICAL: Resolve compound identity before any analysis.
Input Types Handled
| Input Format |
Resolution Strategy |
| Drug name (e.g., "Aspirin") |
PubChem_get_CID_by_compound_name -> get SMILES from properties |
| SMILES string |
Use directly for ADMET-AI; resolve to CID for other tools |
| PubChem CID |
PubChem_get_compound_properties_by_CID -> get SMILES + name |
| ChEMBL ID |
ChEMBL_get_molecule -> get SMILES + properties |
Resolution Steps
- Input detection: Determine if input is name, SMILES, CID, or ChEMBL ID
- SMILES: contains typical SMILES characters (=, #, [, ], (, ), c, n, o and no spaces in middle)
- CID: numeric only
- ChEMBL: starts with "CHEMBL"
- Otherwise: treat as compound name
- Name to CID:
PubChem_get_CID_by_compound_name(name=<compound_name>)
- CID to properties:
PubChem_get_compound_properties_by_CID(cid=<cid>)
- Extract SMILES: Get SMILES from PubChem properties (field:
ConnectivitySMILES, CanonicalSMILES, or IsomericSMILES depending on response format)
- Store resolved IDs: Maintain dict with
name, smiles, cid, formula, weight, inchi
Disambiguation Output
## Compound Identity
| Property | Value |
|----------|-------|
| **Name** | Acetaminophen |
| **PubChem CID** | 1983 |
| **SMILES** | CC(=O)Nc1ccc(O)cc1 |
| **Formula** | C8H9NO2 |
| **Molecular Weight** | 151.16 |
| **InChI** | InChI=1S/C8H9NO2/... |
Phase 1: Predictive Toxicology (ADMET-AI)
When: SMILES is available (from Phase 0 or provided directly)
Objective: Run comprehensive AI-predicted toxicity endpoints
Tools Used
All ADMET-AI tools take the same parameter format:
| Tool |
Predicted Endpoints |
Parameter |
ADMETAI_predict_toxicity |
AMES, Carcinogens_Lagunin, ClinTox, DILI, LD50_Zhu, Skin_Reaction, hERG |
smiles: list[str] |
ADMETAI_predict_stress_response |
Stress response pathway activation (ARE, ATAD5, HSE, MMP, p53) |
smiles: list[str] |
ADMETAI_predict_nuclear_receptor_activity |
AhR, AR, ER, PPARg, Aromatase nuclear receptor activity |
smiles: list[str] |
Workflow
- Call
ADMETAI_predict_toxicity(smiles=[resolved_smiles])
- Call
ADMETAI_predict_stress_response(smiles=[resolved_smiles])
- Call
ADMETAI_predict_nuclear_receptor_activity(smiles=[resolved_smiles])
- For each endpoint, interpret prediction:
- Classification endpoints: Active (1) = toxic signal, Inactive (0) = no signal
- Regression endpoints (LD50): Report numerical value with context
- All predictions graded [T3] (computational prediction)
Decision Logic
- Multiple SMILES: Can batch up to ~10 SMILES in single call
- Failed prediction: If ADMET-AI fails, note "prediction unavailable" (don't fail entire report)
- Confidence: Note that AI predictions are [T3] evidence, not definitive
- hERG flag: If hERG = Active, flag prominently (cardiac safety risk)
- AMES flag: If AMES = Active, flag prominently (mutagenicity concern)
- DILI flag: If DILI = Active, flag prominently (liver toxicity concern)
Output Table
### Toxicity Predictions [T3]
| Endpoint | Prediction | Interpretation | Concern Level |
|----------|-----------|---------------|---------------|
| AMES Mutagenicity | Inactive | No mutagenic signal | Low |
| Carcinogenicity | Inactive | No carcinogenic signal | Low |
| ClinTox | Active | Clinical toxicity signal | HIGH |
| DILI | Active | Drug-induced liver injury risk | HIGH |
| LD50 (Zhu) | 2.45 log(mg/kg) | ~282 mg/kg (moderate) | Medium |
| Skin Reaction | Inactive | No skin sensitization signal | Low |
| hERG Inhibition | Active | Cardiac arrhythmia risk | HIGH |
*All predictions from ADMET-AI. Evidence tier: [T3] (computational prediction)*
Phase 2: ADMET Properties
When: SMILES is available
Objective: Full ADMET characterization beyond toxicity
Tools Used
| Tool |
Properties Predicted |
Parameter |
ADMETAI_predict_BBB_penetrance |
Blood-brain barrier crossing probability |
smiles: list[str] |
ADMETAI_predict_bioavailability |
Oral bioavailability (F20%, F30%) |
smiles: list[str] |
ADMETAI_predict_clearance_distribution |
Clearance, VDss, half-life, PPB |
smiles: list[str] |
ADMETAI_predict_CYP_interactions |
CYP1A2, 2C9, 2C19, 2D6, 3A4 inhibition/substrate |
smiles: list[str] |
ADMETAI_predict_physicochemical_properties |
LogP, LogD, LogS, MW, pKa |
smiles: list[str] |
ADMETAI_predict_solubility_lipophilicity_hydration |
Aqueous solubility, lipophilicity, hydration free energy |
smiles: list[str] |
Workflow
- Call all 6 ADMET tools in parallel (independent calls)
- Compile results into Absorption / Distribution / Metabolism / Excretion sections
- Assess Lipinski Rule of 5 compliance from physicochemical properties
- Flag drug-drug interaction risks from CYP inhibition profiles
Decision Logic
- BBB penetrant + toxicity: If BBB = Yes and any CNS toxicity endpoint active, flag as neurotoxicity risk
- Low bioavailability: If F20% = Low, note absorption concerns
- CYP inhibitor: If CYP3A4 inhibitor = Yes, flag high DDI risk
- Lipinski violations: Count violations and report drug-likeness assessment
Output Format
### ADMET Profile [T3]
#### Absorption
| Property | Value | Interpretation |
|----------|-------|----------------|
| BBB Penetrance | Yes | Crosses blood-brain barrier |
| Bioavailability (F20%) | 85% | Good oral absorption |
#### Distribution
| Property | Value | Interpretation |
|----------|-------|----------------|
| VDss | 1.2 L/kg | Moderate tissue distribution |
| PPB | 92% | Highly protein bound |
#### Metabolism
| CYP Enzyme | Substrate | Inhibitor |
|------------|-----------|-----------|
| CYP1A2 | No | No |
| CYP2C9 | Yes | No |
| CYP2C19 | No | No |
| CYP2D6 | No | No |
| CYP3A4 | Yes | Yes (DDI risk) |
#### Excretion
| Property | Value | Interpretation |
|----------|-------|----------------|
| Clearance | 8.5 mL/min/kg | Moderate clearance |
| Half-life | 6.2 h | Moderate half-life |
Phase 3: Toxicogenomics (CTD)
When: Compound name is resolved
Objective: Map chemical-gene-disease relationships from curated CTD data
Tools Used
| Tool |
Function |
Parameter |
CTD_get_chemical_gene_interactions |
Genes affected by chemical |
input_terms: str (chemical name) |
CTD_get_chemical_diseases |
Diseases linked to chemical exposure |
input_terms: str (chemical name) |
Workflow
- Call
CTD_get_chemical_gene_interactions(input_terms=compound_name)
- Call
CTD_get_chemical_diseases(input_terms=compound_name)
- Parse gene interactions: extract gene symbols, interaction types (increases/decreases expression, binding, etc.)
- Parse disease associations: extract disease names, evidence types (marker/mechanism/therapeutic)
- Identify most affected biological processes from gene list
Decision Logic
- Direct evidence vs inferred: CTD separates curated direct evidence from inferred associations
- Therapeutic vs toxic: Disease associations can be therapeutic (drug treats disease) or adverse (chemical causes disease)
- Gene interaction types: Distinguish between expression changes, binding, and activity modulation
- Prioritize marker/mechanism: These indicate stronger causal evidence than simple associations
- Grade curated as [T2]: Direct curated CTD evidence from literature
- Grade inferred as [T3]: Computationally inferred associations
Output Format
### Toxicogenomics (CTD) [T2/T3]
#### Chemical-Gene Interactions (Top 20)
| Gene | Interaction | Type | Evidence |
|------|------------|------|----------|
| CYP1A2 | increases expression | mRNA | [T2] curated |
| TP53 | affects activity | protein | [T2] curated |
| ... | ... | ... | ... |
**Total interactions found**: 156
**Top affected pathways**: Xenobiotic metabolism, Apoptosis, DNA damage response
#### Chemical-Disease Associations (Top 10)
| Disease | Association Type | Evidence |
|---------|-----------------|----------|
| Liver Neoplasms | marker/mechanism | [T2] curated |
| Contact Dermatitis | therapeutic | [T2] curated |
| ... | ... | ... |
Phase 4: Regulatory Safety (FDA Labels)
When: Compound has an approved drug name
Objective: Extract regulatory safety information from FDA drug labels
Tools Used
| Tool |
Information Retrieved |
Parameter |
FDA_get_boxed_warning_info_by_drug_name |
Black box warnings (most serious) |
drug_name: str |
FDA_get_contraindications_by_drug_name |
Absolute contraindications |
drug_name: str |
FDA_get_adverse_reactions_by_drug_name |
Known adverse reactions |
drug_name: str |
FDA_get_warnings_by_drug_name |
Warnings and precautions |
drug_name: str |
FDA_get_nonclinical_toxicology_info_by_drug_name |
Animal toxicology data |
drug_name: str |
FDA_get_carcinogenic_mutagenic_fertility_by_drug_name |
Carcinogenicity/mutagenicity/fertility data |
drug_name: str |
Workflow
- Call all 6 FDA tools in parallel (independent queries by drug name)
- Parse and structure each response
- Prioritize: Boxed Warnings > Contraindications > Warnings > Adverse Reactions
- All FDA label data is [T1] evidence (regulatory finding based on human/animal data)
Decision Logic
- Boxed warning present: Flag as CRITICAL safety concern in executive summary
- No FDA data: Chemical may not be an approved drug; note "Not an FDA-approved drug" and continue with other phases
- Multiple warnings: Categorize by organ system (hepatic, cardiac, renal, CNS, etc.)
- Nonclinical toxicology: Grade as [T2] (animal data supporting human risk)
Output Format
### Regulatory Safety (FDA) [T1]
#### Boxed Warning
**PRESENT** - Hepatotoxicity risk with doses >4g/day. Liver failure reported. [T1]
#### Contraindications
- Severe hepatic impairment [T1]
- Known hypersensitivity [T1]
#### Adverse Reactions (by frequency)
| Reaction | Frequency | Severity |
|----------|-----------|----------|
| Nausea | Common (>1%) | Mild |
| Hepatotoxicity | Rare (<0.1%) | Severe |
| ... | ... | ... |
#### Nonclinical Toxicology [T2]
- **Carcinogenicity**: No carcinogenic potential in 2-year rat/mouse studies
- **Mutagenicity**: Negative in Ames assay and in vivo micronucleus test
- **Fertility**: No effects on fertility at doses up to 10x human dose
Phase 5: Drug Safety Profile (DrugBank)
When: Compound is a known drug
Objective: Retrieve curated drug safety data from DrugBank
Tools Used
| Tool |
Information |
Parameters |
drugbank_get_safety_by_drug_name_or_drugbank_id |
Toxicity, contraindications |
query: str, case_sensitive: bool, exact_match: bool, limit: int |
Workflow
- Call
drugbank_get_safety_by_drug_name_or_drugbank_id(query=drug_name, case_sensitive=False, exact_match=False, limit=5)
- Parse toxicity information, overdose data, contraindications
- Cross-reference with FDA data from Phase 4
Decision Logic
- Toxicity field: Contains LD50 values, overdose symptoms, organ toxicity data
- DrugBank ID: Note if found for cross-referencing
- Conflict with FDA: If DrugBank and FDA disagree, note discrepancy and defer to FDA [T1]
- Not found: Chemical may not be in DrugBank; continue with other phases
Phase 6: Chemical-Protein Interactions (STITCH)
When: Compound can be identified by name or SMILES
Objective: Map chemical-protein interaction network for off-target assessment
Tools Used
| Tool |
Function |
Parameters |
STITCH_resolve_identifier |
Resolve chemical name to STITCH ID |
identifier: str, species: int (9606=human) |
STITCH_get_chemical_protein_interactions |
Get chemical-protein interactions |
identifiers: list[str], species: int, required_score: int |
STITCH_get_interaction_partners |
Get interaction network |
identifiers: list[str], species: int, limit: int |
Workflow
- Resolve compound:
STITCH_resolve_identifier(identifier=compound_name, species=9606)
- Get interactions:
STITCH_get_chemical_protein_interactions(identifiers=[stitch_id], species=9606, required_score=700)
- Identify off-target proteins (not the intended drug target)
- Flag safety-relevant targets: hERG (cardiac), CYP enzymes (metabolism), nuclear receptors (endocrine)
Decision Logic
- High confidence (>900): Well-established interaction [T2]
- Medium confidence (700-900): Probable interaction [T3]
- Low confidence (400-700): Possible interaction, needs validation [T4]
- Safety-relevant targets: Flag interactions with known safety targets
- No STITCH data: Chemical may be too novel; note and continue
Phase 7: Structural Alerts (ChEMBL)
When: ChEMBL molecule ID is available (from Phase 0)
Objective: Check for known toxic substructures
Tools Used
| Tool |
Function |
Parameters |
ChEMBL_search_compound_structural_alerts |
Find structural alert matches |
molecule_chembl_id: str, limit: int |
Workflow
- If ChEMBL ID available:
ChEMBL_search_compound_structural_alerts(molecule_chembl_id=chembl_id, limit=20)
- Parse alert types: PAINS (pan-assay interference), Brenk (medicinal chemistry), Glaxo (GSK structural alerts)
- Categorize severity: Some alerts are informational, others indicate likely toxicity
Decision Logic
- PAINS alerts: May cause false positives in screening; note for medicinal chemistry
- Brenk alerts: Known problematic substructures; flag if present
- No alerts: Good sign but not definitive proof of safety
- No ChEMBL ID: Skip this phase gracefully; note "structural alert analysis not available"
Synthesis: Integrated Risk Assessment (MANDATORY)
Always the final section. Integrates all evidence into actionable risk classification.
Risk Classification Matrix
| Risk Level |
Criteria |
| CRITICAL |
FDA boxed warning present OR multiple [T1] toxicity findings OR active DILI + active hERG |
| HIGH |
FDA warnings present OR [T2] animal toxicity OR multiple active ADMET endpoints |
| MEDIUM |
Some [T3] predictions positive OR CTD disease associations OR structural alerts |
| LOW |
All ADMET endpoints negative AND no FDA/DrugBank safety flags AND no CTD concerns |
| INSUFFICIENT DATA |
Fewer than 3 phases returned data; cannot make confident assessment |
Synthesis Template
## Integrated Risk Assessment
### Overall Risk Classification: [HIGH]
### Evidence Summary
| Dimension | Finding | Evidence Tier | Concern |
|-----------|---------|--------------|---------|
| ADMET Toxicity | DILI active, hERG active | [T3] | HIGH |
| FDA Label | Boxed warning for hepatotoxicity | [T1] | CRITICAL |
| CTD Toxicogenomics | 156 gene interactions, liver neoplasms | [T2] | HIGH |
| DrugBank | Known hepatotoxicity at high doses | [T2] | HIGH |
| STITCH | Binds CYP3A4, hERG | [T3] | MEDIUM |
| Structural Alerts | 2 Brenk alerts | [T3] | MEDIUM |
### Key Safety Concerns
1. **Hepatotoxicity** [T1]: FDA boxed warning + ADMET-AI DILI prediction + CTD liver disease associations
2. **Cardiac Risk** [T3]: ADMET-AI hERG prediction + STITCH hERG interaction
3. **Drug Interactions** [T3]: CYP3A4 substrate/inhibitor, potential DDI risk
### Data Gaps
- [ ] No in vivo genotoxicity data available
- [ ] STITCH interaction scores moderate (700-900)
- [ ] No environmental exposure data
### Recommendations
1. Avoid doses >4g/day (hepatotoxicity threshold) [T1]
2. Monitor liver function in chronic use [T1]
3. Screen for CYP3A4 interactions before co-administration [T3]
4. Consider cardiac monitoring for at-risk patients [T3]
Mandatory Completeness Checklist
Before finalizing any report, verify:
Tool Parameter Reference
Critical Parameter Notes (verified from source code):
| Tool |
Parameter Name |
Type |
Notes |
| All ADMETAI tools |
smiles |
list[str] |
Always a list, even for single compound |
| All CTD tools |
input_terms |
str |
Chemical name, MeSH name, CAS RN, or MeSH ID |
| All FDA tools |
drug_name |
str |
Brand or generic drug name |
| drugbank_get_safety_* |
query, case_sensitive, exact_match, limit |
str, bool, bool, int |
All 4 required |
| STITCH_resolve_identifier |
identifier, species |
str, int |
species=9606 for human |
| STITCH_get_chemical_protein_interactions |
identifiers, species, required_score |
list[str], int, int |
required_score=400 default |
| PubChem_get_CID_by_compound_name |
name |
str |
Compound name (not SMILES) |
| PubChem_get_compound_properties_by_CID |
cid |
int |
Numeric CID |
| ChEMBL_search_compound_structural_alerts |
molecule_chembl_id |
str |
ChEMBL ID (e.g., "CHEMBL112") |
Response Format Notes
- ADMET-AI: Returns
{status: "success", data: {...}} with prediction values
- CTD: Returns list of interaction/association objects
- FDA: Returns
{status, data} with label text
- DrugBank: Returns
{data: [...]} with drug records
- STITCH: Returns list of interaction objects with scores
- PubChem CID lookup: Returns
{IdentifierList: {CID: [...]}} (may or may not have data wrapper)
- PubChem properties: Returns dict with
CID, MolecularWeight, ConnectivitySMILES, IUPACName
Fallback Strategies
Compound Resolution
- Primary: PubChem by name -> CID -> properties -> SMILES
- Fallback 1: ChEMBL search by name -> molecule -> SMILES
- Fallback 2: If SMILES provided directly, skip name resolution
Toxicity Prediction
- Primary: All 9 ADMET-AI endpoints
- Fallback: If ADMET-AI fails for a compound, note "prediction failed" and continue with database evidence
- Note: ADMET-AI may fail for very large or unusual SMILES
Regulatory Data
- Primary: FDA labels by drug name
- Fallback: If FDA returns no data, try alternative drug names (brand vs generic)
- Note: Non-drug chemicals (pesticides, industrial) will not have FDA labels
CTD Data
- Primary: Search by common chemical name
- Fallback: Try MeSH name if common name fails
- Note: Novel compounds may not be in CTD
Common Use Patterns
Pattern 1: Novel Compound Assessment
Input: SMILES string for new molecule
Workflow: Phase 0 (SMILES->CID) -> Phase 1 (toxicity) -> Phase 2 (ADMET) -> Phase 7 (structural alerts) -> Synthesis
Output: Predictive safety profile for novel compound
Pattern 2: Approved Drug Safety Review
Input: Drug name (e.g., "Acetaminophen")
Workflow: All phases (0-7 + Synthesis)
Output: Complete safety dossier with regulatory + predictive + database evidence
Pattern 3: Environmental Chemical Risk
Input: Chemical name (e.g., "Bisphenol A")
Workflow: Phase 0 -> Phase 1 -> Phase 2 -> Phase 3 (CTD, key for env chemicals) -> Phase 6 -> Synthesis
Output: Environmental health risk assessment focused on gene-disease associations
Pattern 4: Batch Toxicity Screening
Input: Multiple SMILES strings
Workflow: Phase 0 -> Phase 1 (batch) -> Phase 2 (batch) -> Comparative table -> Synthesis
Output: Comparative toxicity table ranking compounds by safety
Pattern 5: Toxicogenomic Deep-Dive
Input: Chemical name + specific gene or disease interest
Workflow: Phase 0 -> Phase 3 (CTD expanded) -> Literature search -> Synthesis
Output: Detailed chemical-gene-disease mechanistic analysis
Output Report Structure
All analyses generate a structured markdown report with progressive sections:
# Chemical Safety & Toxicology Report: [Compound Name]
**Generated**: YYYY-MM-DD HH:MM
**Compound**: [Name] | SMILES: [SMILES] | CID: [CID]
## Executive Summary
[2-3 sentence overview with risk classification and key findings, all graded]
## 1. Compound Identity
[Phase 0 results - disambiguation table]
## 2. Predictive Toxicology
[Phase 1 results - ADMET-AI toxicity endpoints]
## 3. ADMET Profile
[Phase 2 results - absorption, distribution, metabolism, excretion]
## 4. Toxicogenomics
[Phase 3 results - CTD chemical-gene-disease relationships]
## 5. Regulatory Safety
[Phase 4 results - FDA label information]
## 6. Drug Safety Profile
[Phase 5 results - DrugBank data]
## 7. Chemical-Protein Interactions
[Phase 6 results - STITCH network]
## 8. Structural Alerts
[Phase 7 results - ChEMBL alerts]
## 9. Integrated Risk Assessment
[Synthesis - risk classification, evidence summary, data gaps, recommendations]
## Appendix: Methods and Data Sources
[Tool versions, databases queried, date of access]
Limitations & Known Issues
Tool-Specific
- ADMET-AI: Predictions are computational [T3]; should not replace experimental testing
- CTD: Curated but may lag behind latest literature by 6-12 months
- FDA: Only covers FDA-approved drugs; not applicable to environmental chemicals or supplements
- DrugBank: Primarily drugs; limited coverage of industrial chemicals
- STITCH: Score thresholds affect sensitivity; lower scores increase false positives
- ChEMBL: Structural alerts require ChEMBL ID; not all compounds have one
Analysis
- Novel compounds: May only have ADMET-AI predictions (no database evidence)
- Environmental chemicals: FDA/DrugBank phases will be empty; rely on CTD and ADMET-AI
- Batch mode: ADMET-AI can handle batches; other tools require individual queries
- Species specificity: Most data is human-centric; animal data noted where applicable
Technical
- SMILES validity: Invalid SMILES will cause ADMET-AI failures
- Name ambiguity: Chemical names can be ambiguous; always verify with CID
- Rate limits: Some FDA endpoints may rate-limit for rapid queries
Summary
Chemical Safety & Toxicology Assessment Skill provides comprehensive safety evaluation by integrating:
- Predictive toxicology (ADMET-AI) - 9 tools covering toxicity, ADMET, physicochemical properties
- Toxicogenomics (CTD) - Chemical-gene-disease relationship mapping
- Regulatory safety (FDA) - 6 tools for label-based safety extraction
- Drug safety (DrugBank) - Curated toxicity and contraindication data
- Chemical interactions (STITCH) - Chemical-protein interaction networks
- Structural alerts (ChEMBL) - Known toxic substructure detection
Outputs: Structured markdown report with risk classification, evidence grading, and actionable recommendations
Best for: Drug safety assessment, chemical hazard profiling, environmental toxicology, ADMET characterization, toxicogenomic analysis
Total tools integrated: 25+ tools across 6 databases
1---2name: tooluniverse-chemical-safety3description: Comprehensive chemical safety and toxicology assessment integrating ADMET-AI predictions, CTD toxicogenomics, FDA label safety data, DrugBank safety profiles, and STITCH chemical-protein interactions. Performs predictive toxicology (AMES, DILI, LD50, carcinogenicity), organ/system toxicity profiling, chemical-gene-disease relationship mapping, regulatory safety extraction, and environmental hazard assessment. Use when asked about chemical toxicity, drug safety profiling, ADMET properties, environmental health risks, chemical hazard assessment, or toxicogenomic analysis.4---5
6# Chemical Safety & Toxicology Assessment
7
8Comprehensive chemical safety and toxicology analysis integrating predictive AI models, curated toxicogenomics databases, regulatory safety data, and chemical-biological interaction networks. Generates structured risk assessment reports with evidence grading.
9
10## When to Use This Skill
11
12**Triggers**:
13- "Is this chemical toxic?" / "What are the toxicity endpoints for [compound]?"
14- "Assess the safety profile of [drug/chemical]"
15- "What are the ADMET properties of [SMILES]?"
16- "What genes does [chemical] interact with?"
17- "What diseases are linked to [chemical] exposure?"
18- "Predict toxicity for these molecules"
19- "Drug safety assessment for [drug name]"
20- "Environmental health risk of [chemical]"
21- "Chemical hazard profiling"
22- "Toxicogenomic analysis of [compound]"
23
24**Use Cases**:
251. **Predictive Toxicology**: AI-predicted toxicity endpoints (AMES mutagenicity, DILI, LD50, carcinogenicity, skin reactions) for novel compounds via SMILES
262. **ADMET Profiling**: Full absorption, distribution, metabolism, excretion, toxicity characterization
273. **Toxicogenomics**: Chemical-gene interaction mapping, gene-disease associations from CTD
284. **Regulatory Safety**: FDA label warnings, boxed warnings, contraindications, adverse reactions
295. **Drug Safety Assessment**: Combined DrugBank safety + FDA labels + adverse event data
306. **Chemical-Protein Interactions**: STITCH-based chemical-protein binding and interaction networks
317. **Environmental Toxicology**: Chemical-disease associations for environmental contaminants
32
33---
34
35## KEY PRINCIPLES
36
371. **Report-first approach** - Create report file FIRST, then populate progressively
382. **Tool parameter verification** - Verify params via `get_tool_info` before calling unfamiliar tools
393. **Evidence grading** - Grade all safety claims by evidence strength (T1-T4)
404. **Citation requirements** - Every toxicity finding must have inline source attribution
415. **Mandatory completeness** - All sections must exist with data minimums or explicit "No data" notes
426. **Disambiguation first** - Resolve compound identity (name -> SMILES, CID, ChEMBL ID) before analysis
437. **Negative results documented** - "No toxicity signals found" is data; empty sections are failures
448. **Conservative risk assessment** - When evidence is ambiguous, flag as "requires further investigation"
459. **English-first queries** - Always use English chemical/drug names in tool calls
46
47---
48
49## Evidence Grading System (MANDATORY)
50
51Grade every toxicity claim by evidence strength:
52
53| Tier | Symbol | Criteria | Examples |
54|------|--------|----------|----------|
55| **T1** | [T1] | Direct human evidence, regulatory finding | FDA boxed warning, clinical trial toxicity, human case reports |
56| **T2** | [T2] | Animal studies, validated in vitro | Nonclinical toxicology, AMES positive, animal LD50 |
57| **T3** | [T3] | Computational prediction, association data | ADMET-AI prediction, CTD association, QSAR model |
58| **T4** | [T4] | Database annotation, text-mined | Literature mention, database entry without validation |
59
60### Required Evidence Grading Locations
61
62Evidence grades MUST appear in:
631. **Executive Summary** - Key toxicity findings graded
642. **Toxicity Predictions** - Every ADMET-AI endpoint with confidence note
653. **Regulatory Safety** - FDA findings marked [T1]
664. **Chemical-Gene Interactions** - CTD data marked by curation status
675. **Risk Assessment** - Final risk classification with supporting evidence tiers
68
69---
70
71## Core Strategy: 8 Research Dimensions
72
73```
74Chemical/Drug Query
75|
76+-- PHASE 0: Compound Disambiguation (ALWAYS FIRST)
77| +-- Resolve name -> SMILES, PubChem CID, ChEMBL ID
78| +-- Get molecular formula, weight, canonical structure
79|
80+-- PHASE 1: Predictive Toxicology (ADMET-AI)
81| +-- Mutagenicity (AMES)
82| +-- Hepatotoxicity (DILI, ClinTox)
83| +-- Carcinogenicity
84| +-- Acute toxicity (LD50)
85| +-- Skin reactions
86| +-- Stress response pathways
87| +-- Nuclear receptor activity
88|
89+-- PHASE 2: ADMET Properties
90| +-- Absorption: BBB penetrance, bioavailability
91| +-- Distribution: clearance, volume of distribution
92| +-- Metabolism: CYP interactions (1A2, 2C9, 2C19, 2D6, 3A4)
93| +-- Physicochemical: solubility, lipophilicity, pKa
94|
95+-- PHASE 3: Toxicogenomics (CTD)
96| +-- Chemical-gene interactions
97| +-- Chemical-disease associations
98| +-- Affected biological pathways
99|
100+-- PHASE 4: Regulatory Safety (FDA Labels)
101| +-- Boxed warnings (Black Box)
102| +-- Contraindications
103| +-- Adverse reactions
104| +-- Warnings and precautions
105| +-- Nonclinical toxicology
106|
107+-- PHASE 5: Drug Safety Profile (DrugBank)
108| +-- Toxicity data
109| +-- Contraindications
110| +-- Drug interactions affecting safety
111|
112+-- PHASE 6: Chemical-Protein Interactions (STITCH)
113| +-- Direct chemical-protein binding
114| +-- Interaction confidence scores
115| +-- Off-target effects
116|
117+-- PHASE 7: Structural Alerts (ChEMBL)
118| +-- Known toxic substructures (PAINS, Brenk)
119| +-- Structural alert flags
120|
121+-- SYNTHESIS: Integrated Risk Assessment
122 +-- Aggregate all evidence tiers
123 +-- Risk classification (Low/Medium/High/Critical)
124 +-- Data gaps and recommendations
125```
126
127---
128
129## Phase 0: Compound Disambiguation (ALWAYS FIRST)
130
131**CRITICAL**: Resolve compound identity before any analysis.
132
133### Input Types Handled
134
135| Input Format | Resolution Strategy |
136|-------------|---------------------|
137| Drug name (e.g., "Aspirin") | PubChem_get_CID_by_compound_name -> get SMILES from properties |
138| SMILES string | Use directly for ADMET-AI; resolve to CID for other tools |
139| PubChem CID | PubChem_get_compound_properties_by_CID -> get SMILES + name |
140| ChEMBL ID | ChEMBL_get_molecule -> get SMILES + properties |
141
142### Resolution Steps
143
1441. **Input detection**: Determine if input is name, SMILES, CID, or ChEMBL ID
145 - SMILES: contains typical SMILES characters (=, #, [, ], (, ), c, n, o and no spaces in middle)
146 - CID: numeric only
147 - ChEMBL: starts with "CHEMBL"
148 - Otherwise: treat as compound name
1492. **Name to CID**: `PubChem_get_CID_by_compound_name(name=<compound_name>)`
1503. **CID to properties**: `PubChem_get_compound_properties_by_CID(cid=<cid>)`
1514. **Extract SMILES**: Get SMILES from PubChem properties (field: `ConnectivitySMILES`, `CanonicalSMILES`, or `IsomericSMILES` depending on response format)
1525. **Store resolved IDs**: Maintain dict with `name`, `smiles`, `cid`, `formula`, `weight`, `inchi`
153
154### Disambiguation Output
155
156```markdown
157## Compound Identity
158
159| Property | Value |
160|----------|-------|
161| **Name** | Acetaminophen |
162| **PubChem CID** | 1983 |
163| **SMILES** | CC(=O)Nc1ccc(O)cc1 |
164| **Formula** | C8H9NO2 |
165| **Molecular Weight** | 151.16 |
166| **InChI** | InChI=1S/C8H9NO2/... |
167```
168
169---
170
171## Phase 1: Predictive Toxicology (ADMET-AI)
172
173**When**: SMILES is available (from Phase 0 or provided directly)
174
175**Objective**: Run comprehensive AI-predicted toxicity endpoints
176
177### Tools Used
178
179All ADMET-AI tools take the same parameter format:
180
181| Tool | Predicted Endpoints | Parameter |
182|------|---------------------|-----------|
183| `ADMETAI_predict_toxicity` | AMES, Carcinogens_Lagunin, ClinTox, DILI, LD50_Zhu, Skin_Reaction, hERG | `smiles`: list[str] |
184| `ADMETAI_predict_stress_response` | Stress response pathway activation (ARE, ATAD5, HSE, MMP, p53) | `smiles`: list[str] |
185| `ADMETAI_predict_nuclear_receptor_activity` | AhR, AR, ER, PPARg, Aromatase nuclear receptor activity | `smiles`: list[str] |
186
187### Workflow
188
1891. Call `ADMETAI_predict_toxicity(smiles=[resolved_smiles])`
1902. Call `ADMETAI_predict_stress_response(smiles=[resolved_smiles])`
1913. Call `ADMETAI_predict_nuclear_receptor_activity(smiles=[resolved_smiles])`
1924. For each endpoint, interpret prediction:
193 - Classification endpoints: Active (1) = toxic signal, Inactive (0) = no signal
194 - Regression endpoints (LD50): Report numerical value with context
195 - All predictions graded [T3] (computational prediction)
196
197### Decision Logic
198
199- **Multiple SMILES**: Can batch up to ~10 SMILES in single call
200- **Failed prediction**: If ADMET-AI fails, note "prediction unavailable" (don't fail entire report)
201- **Confidence**: Note that AI predictions are [T3] evidence, not definitive
202- **hERG flag**: If hERG = Active, flag prominently (cardiac safety risk)
203- **AMES flag**: If AMES = Active, flag prominently (mutagenicity concern)
204- **DILI flag**: If DILI = Active, flag prominently (liver toxicity concern)
205
206### Output Table
207
208```markdown
209### Toxicity Predictions [T3]
210
211| Endpoint | Prediction | Interpretation | Concern Level |
212|----------|-----------|---------------|---------------|
213| AMES Mutagenicity | Inactive | No mutagenic signal | Low |
214| Carcinogenicity | Inactive | No carcinogenic signal | Low |
215| ClinTox | Active | Clinical toxicity signal | HIGH |
216| DILI | Active | Drug-induced liver injury risk | HIGH |
217| LD50 (Zhu) | 2.45 log(mg/kg) | ~282 mg/kg (moderate) | Medium |
218| Skin Reaction | Inactive | No skin sensitization signal | Low |
219| hERG Inhibition | Active | Cardiac arrhythmia risk | HIGH |
220
221*All predictions from ADMET-AI. Evidence tier: [T3] (computational prediction)*
222```
223
224---
225
226## Phase 2: ADMET Properties
227
228**When**: SMILES is available
229
230**Objective**: Full ADMET characterization beyond toxicity
231
232### Tools Used
233
234| Tool | Properties Predicted | Parameter |
235|------|---------------------|-----------|
236| `ADMETAI_predict_BBB_penetrance` | Blood-brain barrier crossing probability | `smiles`: list[str] |
237| `ADMETAI_predict_bioavailability` | Oral bioavailability (F20%, F30%) | `smiles`: list[str] |
238| `ADMETAI_predict_clearance_distribution` | Clearance, VDss, half-life, PPB | `smiles`: list[str] |
239| `ADMETAI_predict_CYP_interactions` | CYP1A2, 2C9, 2C19, 2D6, 3A4 inhibition/substrate | `smiles`: list[str] |
240| `ADMETAI_predict_physicochemical_properties` | LogP, LogD, LogS, MW, pKa | `smiles`: list[str] |
241| `ADMETAI_predict_solubility_lipophilicity_hydration` | Aqueous solubility, lipophilicity, hydration free energy | `smiles`: list[str] |
242
243### Workflow
244
2451. Call all 6 ADMET tools in parallel (independent calls)
2462. Compile results into Absorption / Distribution / Metabolism / Excretion sections
2473. Assess Lipinski Rule of 5 compliance from physicochemical properties
2484. Flag drug-drug interaction risks from CYP inhibition profiles
249
250### Decision Logic
251
252- **BBB penetrant + toxicity**: If BBB = Yes and any CNS toxicity endpoint active, flag as neurotoxicity risk
253- **Low bioavailability**: If F20% = Low, note absorption concerns
254- **CYP inhibitor**: If CYP3A4 inhibitor = Yes, flag high DDI risk
255- **Lipinski violations**: Count violations and report drug-likeness assessment
256
257### Output Format
258
259```markdown
260### ADMET Profile [T3]
261
262#### Absorption
263| Property | Value | Interpretation |
264|----------|-------|----------------|
265| BBB Penetrance | Yes | Crosses blood-brain barrier |
266| Bioavailability (F20%) | 85% | Good oral absorption |
267
268#### Distribution
269| Property | Value | Interpretation |
270|----------|-------|----------------|
271| VDss | 1.2 L/kg | Moderate tissue distribution |
272| PPB | 92% | Highly protein bound |
273
274#### Metabolism
275| CYP Enzyme | Substrate | Inhibitor |
276|------------|-----------|-----------|
277| CYP1A2 | No | No |
278| CYP2C9 | Yes | No |
279| CYP2C19 | No | No |
280| CYP2D6 | No | No |
281| CYP3A4 | Yes | Yes (DDI risk) |
282
283#### Excretion
284| Property | Value | Interpretation |
285|----------|-------|----------------|
286| Clearance | 8.5 mL/min/kg | Moderate clearance |
287| Half-life | 6.2 h | Moderate half-life |
288```
289
290---
291
292## Phase 3: Toxicogenomics (CTD)
293
294**When**: Compound name is resolved
295
296**Objective**: Map chemical-gene-disease relationships from curated CTD data
297
298### Tools Used
299
300| Tool | Function | Parameter |
301|------|----------|-----------|
302| `CTD_get_chemical_gene_interactions` | Genes affected by chemical | `input_terms`: str (chemical name) |
303| `CTD_get_chemical_diseases` | Diseases linked to chemical exposure | `input_terms`: str (chemical name) |
304
305### Workflow
306
3071. Call `CTD_get_chemical_gene_interactions(input_terms=compound_name)`
3082. Call `CTD_get_chemical_diseases(input_terms=compound_name)`
3093. Parse gene interactions: extract gene symbols, interaction types (increases/decreases expression, binding, etc.)
3104. Parse disease associations: extract disease names, evidence types (marker/mechanism/therapeutic)
3115. Identify most affected biological processes from gene list
312
313### Decision Logic
314
315- **Direct evidence vs inferred**: CTD separates curated direct evidence from inferred associations
316- **Therapeutic vs toxic**: Disease associations can be therapeutic (drug treats disease) or adverse (chemical causes disease)
317- **Gene interaction types**: Distinguish between expression changes, binding, and activity modulation
318- **Prioritize marker/mechanism**: These indicate stronger causal evidence than simple associations
319- **Grade curated as [T2]**: Direct curated CTD evidence from literature
320- **Grade inferred as [T3]**: Computationally inferred associations
321
322### Output Format
323
324```markdown
325### Toxicogenomics (CTD) [T2/T3]
326
327#### Chemical-Gene Interactions (Top 20)
328| Gene | Interaction | Type | Evidence |
329|------|------------|------|----------|
330| CYP1A2 | increases expression | mRNA | [T2] curated |
331| TP53 | affects activity | protein | [T2] curated |
332| ... | ... | ... | ... |
333
334**Total interactions found**: 156
335**Top affected pathways**: Xenobiotic metabolism, Apoptosis, DNA damage response
336
337#### Chemical-Disease Associations (Top 10)
338| Disease | Association Type | Evidence |
339|---------|-----------------|----------|
340| Liver Neoplasms | marker/mechanism | [T2] curated |
341| Contact Dermatitis | therapeutic | [T2] curated |
342| ... | ... | ... |
343```
344
345---
346
347## Phase 4: Regulatory Safety (FDA Labels)
348
349**When**: Compound has an approved drug name
350
351**Objective**: Extract regulatory safety information from FDA drug labels
352
353### Tools Used
354
355| Tool | Information Retrieved | Parameter |
356|------|---------------------|-----------|
357| `FDA_get_boxed_warning_info_by_drug_name` | Black box warnings (most serious) | `drug_name`: str |
358| `FDA_get_contraindications_by_drug_name` | Absolute contraindications | `drug_name`: str |
359| `FDA_get_adverse_reactions_by_drug_name` | Known adverse reactions | `drug_name`: str |
360| `FDA_get_warnings_by_drug_name` | Warnings and precautions | `drug_name`: str |
361| `FDA_get_nonclinical_toxicology_info_by_drug_name` | Animal toxicology data | `drug_name`: str |
362| `FDA_get_carcinogenic_mutagenic_fertility_by_drug_name` | Carcinogenicity/mutagenicity/fertility data | `drug_name`: str |
363
364### Workflow
365
3661. Call all 6 FDA tools in parallel (independent queries by drug name)
3672. Parse and structure each response
3683. Prioritize: Boxed Warnings > Contraindications > Warnings > Adverse Reactions
3694. All FDA label data is [T1] evidence (regulatory finding based on human/animal data)
370
371### Decision Logic
372
373- **Boxed warning present**: Flag as CRITICAL safety concern in executive summary
374- **No FDA data**: Chemical may not be an approved drug; note "Not an FDA-approved drug" and continue with other phases
375- **Multiple warnings**: Categorize by organ system (hepatic, cardiac, renal, CNS, etc.)
376- **Nonclinical toxicology**: Grade as [T2] (animal data supporting human risk)
377
378### Output Format
379
380```markdown
381### Regulatory Safety (FDA) [T1]
382
383#### Boxed Warning
384**PRESENT** - Hepatotoxicity risk with doses >4g/day. Liver failure reported. [T1]
385
386#### Contraindications
387- Severe hepatic impairment [T1]
388- Known hypersensitivity [T1]
389
390#### Adverse Reactions (by frequency)
391| Reaction | Frequency | Severity |
392|----------|-----------|----------|
393| Nausea | Common (>1%) | Mild |
394| Hepatotoxicity | Rare (<0.1%) | Severe |
395| ... | ... | ... |
396
397#### Nonclinical Toxicology [T2]
398- **Carcinogenicity**: No carcinogenic potential in 2-year rat/mouse studies
399- **Mutagenicity**: Negative in Ames assay and in vivo micronucleus test
400- **Fertility**: No effects on fertility at doses up to 10x human dose
401```
402
403---
404
405## Phase 5: Drug Safety Profile (DrugBank)
406
407**When**: Compound is a known drug
408
409**Objective**: Retrieve curated drug safety data from DrugBank
410
411### Tools Used
412
413| Tool | Information | Parameters |
414|------|------------|------------|
415| `drugbank_get_safety_by_drug_name_or_drugbank_id` | Toxicity, contraindications | `query`: str, `case_sensitive`: bool, `exact_match`: bool, `limit`: int |
416
417### Workflow
418
4191. Call `drugbank_get_safety_by_drug_name_or_drugbank_id(query=drug_name, case_sensitive=False, exact_match=False, limit=5)`
4202. Parse toxicity information, overdose data, contraindications
4213. Cross-reference with FDA data from Phase 4
422
423### Decision Logic
424
425- **Toxicity field**: Contains LD50 values, overdose symptoms, organ toxicity data
426- **DrugBank ID**: Note if found for cross-referencing
427- **Conflict with FDA**: If DrugBank and FDA disagree, note discrepancy and defer to FDA [T1]
428- **Not found**: Chemical may not be in DrugBank; continue with other phases
429
430---
431
432## Phase 6: Chemical-Protein Interactions (STITCH)
433
434**When**: Compound can be identified by name or SMILES
435
436**Objective**: Map chemical-protein interaction network for off-target assessment
437
438### Tools Used
439
440| Tool | Function | Parameters |
441|------|----------|------------|
442| `STITCH_resolve_identifier` | Resolve chemical name to STITCH ID | `identifier`: str, `species`: int (9606=human) |
443| `STITCH_get_chemical_protein_interactions` | Get chemical-protein interactions | `identifiers`: list[str], `species`: int, `required_score`: int |
444| `STITCH_get_interaction_partners` | Get interaction network | `identifiers`: list[str], `species`: int, `limit`: int |
445
446### Workflow
447
4481. Resolve compound: `STITCH_resolve_identifier(identifier=compound_name, species=9606)`
4492. Get interactions: `STITCH_get_chemical_protein_interactions(identifiers=[stitch_id], species=9606, required_score=700)`
4503. Identify off-target proteins (not the intended drug target)
4514. Flag safety-relevant targets: hERG (cardiac), CYP enzymes (metabolism), nuclear receptors (endocrine)
452
453### Decision Logic
454
455- **High confidence (>900)**: Well-established interaction [T2]
456- **Medium confidence (700-900)**: Probable interaction [T3]
457- **Low confidence (400-700)**: Possible interaction, needs validation [T4]
458- **Safety-relevant targets**: Flag interactions with known safety targets
459- **No STITCH data**: Chemical may be too novel; note and continue
460
461---
462
463## Phase 7: Structural Alerts (ChEMBL)
464
465**When**: ChEMBL molecule ID is available (from Phase 0)
466
467**Objective**: Check for known toxic substructures
468
469### Tools Used
470
471| Tool | Function | Parameters |
472|------|----------|------------|
473| `ChEMBL_search_compound_structural_alerts` | Find structural alert matches | `molecule_chembl_id`: str, `limit`: int |
474
475### Workflow
476
4771. If ChEMBL ID available: `ChEMBL_search_compound_structural_alerts(molecule_chembl_id=chembl_id, limit=20)`
4782. Parse alert types: PAINS (pan-assay interference), Brenk (medicinal chemistry), Glaxo (GSK structural alerts)
4793. Categorize severity: Some alerts are informational, others indicate likely toxicity
480
481### Decision Logic
482
483- **PAINS alerts**: May cause false positives in screening; note for medicinal chemistry
484- **Brenk alerts**: Known problematic substructures; flag if present
485- **No alerts**: Good sign but not definitive proof of safety
486- **No ChEMBL ID**: Skip this phase gracefully; note "structural alert analysis not available"
487
488---
489
490## Synthesis: Integrated Risk Assessment (MANDATORY)
491
492**Always the final section**. Integrates all evidence into actionable risk classification.
493
494### Risk Classification Matrix
495
496| Risk Level | Criteria |
497|-----------|----------|
498| **CRITICAL** | FDA boxed warning present OR multiple [T1] toxicity findings OR active DILI + active hERG |
499| **HIGH** | FDA warnings present OR [T2] animal toxicity OR multiple active ADMET endpoints |
500| **MEDIUM** | Some [T3] predictions positive OR CTD disease associations OR structural alerts |
501| **LOW** | All ADMET endpoints negative AND no FDA/DrugBank safety flags AND no CTD concerns |
502| **INSUFFICIENT DATA** | Fewer than 3 phases returned data; cannot make confident assessment |
503
504### Synthesis Template
505
506```markdown
507## Integrated Risk Assessment
508
509### Overall Risk Classification: [HIGH]
510
511### Evidence Summary
512| Dimension | Finding | Evidence Tier | Concern |
513|-----------|---------|--------------|---------|
514| ADMET Toxicity | DILI active, hERG active | [T3] | HIGH |
515| FDA Label | Boxed warning for hepatotoxicity | [T1] | CRITICAL |
516| CTD Toxicogenomics | 156 gene interactions, liver neoplasms | [T2] | HIGH |
517| DrugBank | Known hepatotoxicity at high doses | [T2] | HIGH |
518| STITCH | Binds CYP3A4, hERG | [T3] | MEDIUM |
519| Structural Alerts | 2 Brenk alerts | [T3] | MEDIUM |
520
521### Key Safety Concerns
5221. **Hepatotoxicity** [T1]: FDA boxed warning + ADMET-AI DILI prediction + CTD liver disease associations
5232. **Cardiac Risk** [T3]: ADMET-AI hERG prediction + STITCH hERG interaction
5243. **Drug Interactions** [T3]: CYP3A4 substrate/inhibitor, potential DDI risk
525
526### Data Gaps
527- [ ] No in vivo genotoxicity data available
528- [ ] STITCH interaction scores moderate (700-900)
529- [ ] No environmental exposure data
530
531### Recommendations
5321. Avoid doses >4g/day (hepatotoxicity threshold) [T1]
5332. Monitor liver function in chronic use [T1]
5343. Screen for CYP3A4 interactions before co-administration [T3]
5354. Consider cardiac monitoring for at-risk patients [T3]
536```
537
538---
539
540## Mandatory Completeness Checklist
541
542Before finalizing any report, verify:
543
544- [ ] **Phase 0**: Compound fully disambiguated (SMILES + CID at minimum)
545- [ ] **Phase 1**: At least 5 toxicity endpoints reported or "prediction unavailable" noted
546- [ ] **Phase 2**: ADMET profile with A/D/M/E sections or "not available" noted
547- [ ] **Phase 3**: CTD queried; gene interactions and disease associations reported or "no data in CTD"
548- [ ] **Phase 4**: FDA labels queried; results or "not an FDA-approved drug" noted
549- [ ] **Phase 5**: DrugBank queried; results or "not found in DrugBank" noted
550- [ ] **Phase 6**: STITCH queried; results or "no STITCH data available" noted
551- [ ] **Phase 7**: Structural alerts checked or "ChEMBL ID not available" noted
552- [ ] **Synthesis**: Risk classification provided with evidence summary
553- [ ] **Evidence Grading**: All findings have [T1]-[T4] annotations
554- [ ] **Data Gaps**: Explicitly listed in synthesis section
555
556---
557
558## Tool Parameter Reference
559
560**Critical Parameter Notes** (verified from source code):
561
562| Tool | Parameter Name | Type | Notes |
563|------|---------------|------|-------|
564| All ADMETAI tools | `smiles` | `list[str]` | Always a list, even for single compound |
565| All CTD tools | `input_terms` | `str` | Chemical name, MeSH name, CAS RN, or MeSH ID |
566| All FDA tools | `drug_name` | `str` | Brand or generic drug name |
567| drugbank_get_safety_* | `query`, `case_sensitive`, `exact_match`, `limit` | str, bool, bool, int | All 4 required |
568| STITCH_resolve_identifier | `identifier`, `species` | str, int | species=9606 for human |
569| STITCH_get_chemical_protein_interactions | `identifiers`, `species`, `required_score` | list[str], int, int | required_score=400 default |
570| PubChem_get_CID_by_compound_name | `name` | `str` | Compound name (not SMILES) |
571| PubChem_get_compound_properties_by_CID | `cid` | `int` | Numeric CID |
572| ChEMBL_search_compound_structural_alerts | `molecule_chembl_id` | `str` | ChEMBL ID (e.g., "CHEMBL112") |
573
574### Response Format Notes
575
576- **ADMET-AI**: Returns `{status: "success", data: {...}}` with prediction values
577- **CTD**: Returns list of interaction/association objects
578- **FDA**: Returns `{status, data}` with label text
579- **DrugBank**: Returns `{data: [...]}` with drug records
580- **STITCH**: Returns list of interaction objects with scores
581- **PubChem CID lookup**: Returns `{IdentifierList: {CID: [...]}}` (may or may not have `data` wrapper)
582- **PubChem properties**: Returns dict with `CID`, `MolecularWeight`, `ConnectivitySMILES`, `IUPACName`
583
584---
585
586## Fallback Strategies
587
588### Compound Resolution
589- **Primary**: PubChem by name -> CID -> properties -> SMILES
590- **Fallback 1**: ChEMBL search by name -> molecule -> SMILES
591- **Fallback 2**: If SMILES provided directly, skip name resolution
592
593### Toxicity Prediction
594- **Primary**: All 9 ADMET-AI endpoints
595- **Fallback**: If ADMET-AI fails for a compound, note "prediction failed" and continue with database evidence
596- **Note**: ADMET-AI may fail for very large or unusual SMILES
597
598### Regulatory Data
599- **Primary**: FDA labels by drug name
600- **Fallback**: If FDA returns no data, try alternative drug names (brand vs generic)
601- **Note**: Non-drug chemicals (pesticides, industrial) will not have FDA labels
602
603### CTD Data
604- **Primary**: Search by common chemical name
605- **Fallback**: Try MeSH name if common name fails
606- **Note**: Novel compounds may not be in CTD
607
608---
609
610## Common Use Patterns
611
612### Pattern 1: Novel Compound Assessment
613```
614Input: SMILES string for new molecule
615Workflow: Phase 0 (SMILES->CID) -> Phase 1 (toxicity) -> Phase 2 (ADMET) -> Phase 7 (structural alerts) -> Synthesis
616Output: Predictive safety profile for novel compound
617```
618
619### Pattern 2: Approved Drug Safety Review
620```
621Input: Drug name (e.g., "Acetaminophen")
622Workflow: All phases (0-7 + Synthesis)
623Output: Complete safety dossier with regulatory + predictive + database evidence
624```
625
626### Pattern 3: Environmental Chemical Risk
627```
628Input: Chemical name (e.g., "Bisphenol A")
629Workflow: Phase 0 -> Phase 1 -> Phase 2 -> Phase 3 (CTD, key for env chemicals) -> Phase 6 -> Synthesis
630Output: Environmental health risk assessment focused on gene-disease associations
631```
632
633### Pattern 4: Batch Toxicity Screening
634```
635Input: Multiple SMILES strings
636Workflow: Phase 0 -> Phase 1 (batch) -> Phase 2 (batch) -> Comparative table -> Synthesis
637Output: Comparative toxicity table ranking compounds by safety
638```
639
640### Pattern 5: Toxicogenomic Deep-Dive
641```
642Input: Chemical name + specific gene or disease interest
643Workflow: Phase 0 -> Phase 3 (CTD expanded) -> Literature search -> Synthesis
644Output: Detailed chemical-gene-disease mechanistic analysis
645```
646
647---
648
649## Output Report Structure
650
651All analyses generate a structured markdown report with progressive sections:
652
653```markdown
654# Chemical Safety & Toxicology Report: [Compound Name]
655
656**Generated**: YYYY-MM-DD HH:MM
657**Compound**: [Name] | SMILES: [SMILES] | CID: [CID]
658
659## Executive Summary
660[2-3 sentence overview with risk classification and key findings, all graded]
661
662## 1. Compound Identity
663[Phase 0 results - disambiguation table]
664
665## 2. Predictive Toxicology
666[Phase 1 results - ADMET-AI toxicity endpoints]
667
668## 3. ADMET Profile
669[Phase 2 results - absorption, distribution, metabolism, excretion]
670
671## 4. Toxicogenomics
672[Phase 3 results - CTD chemical-gene-disease relationships]
673
674## 5. Regulatory Safety
675[Phase 4 results - FDA label information]
676
677## 6. Drug Safety Profile
678[Phase 5 results - DrugBank data]
679
680## 7. Chemical-Protein Interactions
681[Phase 6 results - STITCH network]
682
683## 8. Structural Alerts
684[Phase 7 results - ChEMBL alerts]
685
686## 9. Integrated Risk Assessment
687[Synthesis - risk classification, evidence summary, data gaps, recommendations]
688
689## Appendix: Methods and Data Sources
690[Tool versions, databases queried, date of access]
691```
692
693---
694
695## Limitations & Known Issues
696
697### Tool-Specific
698- **ADMET-AI**: Predictions are computational [T3]; should not replace experimental testing
699- **CTD**: Curated but may lag behind latest literature by 6-12 months
700- **FDA**: Only covers FDA-approved drugs; not applicable to environmental chemicals or supplements
701- **DrugBank**: Primarily drugs; limited coverage of industrial chemicals
702- **STITCH**: Score thresholds affect sensitivity; lower scores increase false positives
703- **ChEMBL**: Structural alerts require ChEMBL ID; not all compounds have one
704
705### Analysis
706- **Novel compounds**: May only have ADMET-AI predictions (no database evidence)
707- **Environmental chemicals**: FDA/DrugBank phases will be empty; rely on CTD and ADMET-AI
708- **Batch mode**: ADMET-AI can handle batches; other tools require individual queries
709- **Species specificity**: Most data is human-centric; animal data noted where applicable
710
711### Technical
712- **SMILES validity**: Invalid SMILES will cause ADMET-AI failures
713- **Name ambiguity**: Chemical names can be ambiguous; always verify with CID
714- **Rate limits**: Some FDA endpoints may rate-limit for rapid queries
715
716---
717
718## Summary
719
720**Chemical Safety & Toxicology Assessment Skill** provides comprehensive safety evaluation by integrating:
721
7221. **Predictive toxicology** (ADMET-AI) - 9 tools covering toxicity, ADMET, physicochemical properties
7232. **Toxicogenomics** (CTD) - Chemical-gene-disease relationship mapping
7243. **Regulatory safety** (FDA) - 6 tools for label-based safety extraction
7254. **Drug safety** (DrugBank) - Curated toxicity and contraindication data
7265. **Chemical interactions** (STITCH) - Chemical-protein interaction networks
7276. **Structural alerts** (ChEMBL) - Known toxic substructure detection
728
729**Outputs**: Structured markdown report with risk classification, evidence grading, and actionable recommendations
730
731**Best for**: Drug safety assessment, chemical hazard profiling, environmental toxicology, ADMET characterization, toxicogenomic analysis
732
733**Total tools integrated**: 25+ tools across 6 databases