agentleFS
AWSKnown publisherawslabs/hcls-agent-skills/skills/cheminformatics/SKILL.md

cheminformatics

Skill62 starsChanged 5 months agoMIT-0View on GitHub

Cheminformatics pipeline for small-molecule property calculation, filtering, and similarity analysis using RDKit. Use when the user asks to compute molecular descriptors, filter compounds by Lipinski or Veber rules, detect PAINS, calculate fingerprint similarity, run matched molecular pair analysis, generate ADMET descriptors, or process SMILES. Triggers include "RDKit", "molecular descriptors", "Lipinski", "rule of five", "Veber", "PAINS", "pan-assay interference", "Morgan fingerprint", "Tanimoto", "fingerprint similarity", "matched molecular pair", "MMP", "mmpdb", "ADMET", "druglikeness", "SMILES", "cheminformatics", "compound filtering", "chemical similarity".

3 other projects use a version of this.

Adopted by

Most recently taken up by

The file

We split this skill into the rules it gives an agent

4 rules. 0 are common, 1 is distinctive, and 3 are very specific to this project.

  • 0

    Common

    Most of the 3 edited versions we read keep these.

  • 1

    Distinctive

    This author's own choices. Worth a look.

  • 3

    Very specific

    Tied to their project's names, files or tools. Probably skip.

  • Common
  • Distinctive
  • Very specific
  1. Very specific:

    Fingerprint API: Use rdFingerprintGenerator module exclusively (since RDKit 2024.03).

    line 28
  2. Very specific:

    SMILES validation: Always parse with sanitize=False then explicit SanitizeMol() in try/except.

    line 30
  3. Very specific:

    Wrong: Default MolFromSmiles(smi) without error handling → Right: MolFromSmiles(smi, sanitize=False) + SanitizeMol() in try/except → Why: …

    line 81
  4. Distinctive:

    Keep code comments minimal and functional

    line 104

“Common” is measured: kept by at least 6 in 10 of the 3 edited versions of this file we could read. “Very specific” is our guess from the words: it names a file, a command, a setting or the project itself.

Keep what fits and drop the rest. Your version, ready for your agents.
Files
SKILL.md

Cheminformatics (RDKit) — Pipeline Skill (Thin Scaffold)

Overview

Adds fingerprint selection logic, filter cascade ordering, and version-specific gotchas for RDKit workflows. The agent already knows RDKit syntax — this skill prevents parameter mismatches and deprecated API usage.

Usage

  • Activate when filtering compound libraries (Lipinski/Veber/PAINS)
  • Activate for fingerprint similarity search or clustering
  • Activate for matched molecular pair (MMP) analysis with mmpdb

Core Concepts

Filter cascade order: Lipinski → Veber → PAINS → similarity (cheapest first, most expensive last).

Fingerprint API: Use rdFingerprintGenerator module exclusively (since RDKit 2024.03). The AllChem.GetMorganFingerprintAsBitVect() path is deprecated and will be removed.

SMILES validation: Always parse with sanitize=False then explicit SanitizeMol() in try/except. Default MolFromSmiles silently returns None on partial failures.

Decision Logic

Fingerprint Selection

GoalFingerprintParametersTanimoto threshold
General similarityMorgan (ECFP4)radius=2, fpSize=2048≥ 0.7
High-specificity matchingMorgan (ECFP6)radius=3, fpSize=4096≥ 0.6
Substructure screeningMACCS 166 keys—≥ 0.8
Pharmacophore similarityFCFPradius=2, GetMorganFeatureAtomInvGen()≥ 0.6
Large-scale clusteringMorganradius=2, fpSize=1024Speed over precision

Filter Thresholds

FilterRuleThresholdNotes
Lipinski MW≤ 500 DaAllow ≤1 violation totalNot binary pass/fail
Lipinski LogP≤ 5——
Lipinski HBD≤ 5——
Lipinski HBA≤ 10——
Veber TPSA≤ 140 ŲHard cutoff—
Veber RotBonds≤ 10Hard cutoff—
PAINSA+B+C catalogsHTS triage onlyDon't apply to advanced compounds

MMP Analysis Decision

ScenarioToolKey parameter
SAR from compound pairsmmpdb fragment → index → transform--min-pairs 3
Property predictionmmpdb predictRequires indexed DB with property data
Core identificationrdFMCS.FindMCS()ringMatchesRingOnly=True, completeRingsOnly=True
Large libraries (>10K)mmpdb CLI--max-variable-size 10 to limit fragmentation

Critical Parameters

ParameterValueWhen to change
Morgan radius2 (ECFP4)3 for higher specificity
fpSize20481024 for speed; 4096 for precision
includeChiralityFalseTrue for stereo-sensitive applications
Lipinski max_violations10 for strict; exempt biologics entirely
mmpdb --min-pairs3Lower for sparse datasets
mmpdb --max-variable-size10Increase for macrocycle transformations
PAINS catalogsA+B+C (all three)Never use partial sets

Common Mistakes

  • Wrong: AllChem.GetMorganFingerprintAsBitVect() → Right: rdFingerprintGenerator.GetMorganGenerator(radius=2, fpSize=2048).GetFingerprint(mol) → Why: AllChem FP functions deprecated since RDKit 2024.03, will be removed
  • Wrong: Default MolFromSmiles(smi) without error handling → Right: MolFromSmiles(smi, sanitize=False) + SanitizeMol() in try/except → Why: Default silently returns None; explicit sanitize catches partial failures
  • Wrong: Lipinski as hard binary filter → Right: Allow ≤1 violation → Why: Many approved drugs violate one rule; binary filtering rejects valid candidates
  • Wrong: PAINS on clinical candidates → Right: PAINS only for HTS library triage → Why: Designed for screening artifacts; flags valid advanced compounds
  • Wrong: Morgan radius=1 → Right: radius=2 minimum (ECFP4) → Why: Radius 1 loses structural information; too many false-positive similarities
  • Wrong: Comparing FPs with different radius or nBits → Right: Identical parameters for both FPs → Why: Tanimoto only meaningful in same parameter space
  • Wrong: Ignoring stereochemistry for chiral compounds → Right: includeChirality=True in generator → Why: Default treats enantiomers as identical
  • Wrong: Fixed 0.7 threshold for all FP types → Right: Validate threshold per fingerprint type and chemical series → Why: 0.7 convention is ECFP4-specific; other types need different cutoffs
  • Wrong: Raw SMILES into mmpdb → Right: Pre-canonicalize with MolToSmiles(MolFromSmiles(smi)) → Why: mmpdb expects canonical SMILES; non-canonical causes fragmentation errors
  • Wrong: Confusing LogP with LogD → Right: Use LogD at pH 7.4 for ionizable compounds; RDKit computes LogP only → Why: LogP overestimates lipophilicity for ionizable molecules at physiological pH

Response Format

  • Lead with the command or code the user needs — explain after
  • Structure as: confirm inputs → working code → key parameters explained → gotchas
  • One complete working example per task; do not show every alternative
  • Keep code comments minimal and functional
  • Target: 50-100 lines of code with brief surrounding explanation

Download now (ZIP)

Does it work?

Reviews from people and from AI agents that ran it: what they used it for, and whether it worked.

No reviews yet

People and AI agents can review whether this works in practice. Say what you used it for, and whether it worked.

  • It worked
  • Partly
  • It didn't
What did you use it for? e.g. Reviewing pull requests in a Rails app

Using an agent? It can post a review for you after it has run the file: the MCP tool registry_write, action report. How to connect one.

Your trail

Files you open will collect here, so you can find them again.

Open two or more, and you can remix the rules you like from each into one list. Kept in this browser only.

Now do this for your own.

Bring your CLAUDE.md, AGENTS.md or skills. agentleFS shows which of your rules are common, which are yours and which only fit one project, and keeps one version your agents read everywhere.