molecule-evolution-agent

$npx mdskill add swaruplab/operon/molecule-evolution-agent

Evolves molecules to optimize binding for a protein target.

  • Improves lead compounds for drug discovery and potency optimization.
  • Uses SMILES format and an LLM to suggest chemical modifications.
  • Iteratively modifies structure based on LLM suggestions and mock scoring.
  • Outputs top candidate SMILES and evolution history.

SKILL.md

.github/skills/molecule-evolution-agentView on GitHub ↗
---
name: 'molecule-evolution-agent'
description: 'Evolve Molecules'
measurable_outcome: Execute skill workflow successfully with valid output within 15 minutes.
allowed-tools:
  - read_file
  - run_shell_command
---

# Molecule Evolution Agent

The **Molecule Evolution Agent** acts as an autonomous medicinal chemist. It takes a starting molecule (or uses a default like Aspirin) and iteratively modifies its structure to optimize binding for a specific protein target.

## When to Use This Skill

*   **Lead Optimization**: When you have a hit molecule and want to improve its potency.
*   **De Novo Design**: To explore chemical space around a target protein.
*   **Idea Generation**: To get creative structural modifications suggested by an LLM.

## Core Capabilities

1.  **SMILES Manipulation**: Reads and writes chemical structures in SMILES format.
2.  **LLM Chemist**: Uses an LLM to suggest chemically valid modifications (e.g., "Add a fluorine group to the ring").
3.  **Mock Scoring**: (Currently) Uses a mock scoring function to simulate docking affinity.

## Workflow

1.  **Input**: Target Protein Name (e.g., "GPRC5D").
2.  **Process**: 
    *   Start with a seed molecule.
    *   Loop for *N* generations.
    *   Ask LLM for a modification.
    *   Score the new molecule.
    *   Keep the best candidate.
3.  **Output**: Top candidate SMILES and the evolution history.

## Example Usage

**User**: "Design a better binder for GPRC5D."

**Agent Action**:
```bash
python3 Skills/Drug_Discovery/Molecule_Design/evolution_agent.py
# (Note: The script currently defaults to GPRC5D, but can be extended for arguments)
```

More from swaruplab/operon

SkillDescription
agentd-drug-discoveryUse the AgentD workflow to mine evidence, design molecules, and rank candidates with SAR plus ADMET annotations for early drug discovery tasks.
alphafold-databaseAccess AlphaFold 200M+ AI-predicted protein structures. Retrieve structures by UniProt ID, download PDB/mmCIF files, analyze confidence metrics (pLDDT, PAE), for drug discovery and structural biology.
antibody-design-agentAn advanced agent for de novo antibody design and optimization using state-of-the-art protein language models (MAGE, RFdiffusion).
archrSingle-cell ATAC-seq analysis with ArchR (R). The mature R-based scATAC pipeline — Arrow files, doublet inference, iterative LSI + Harmony, clustering, gene scores, MACS2 peak calling, motif enrichment, chromVAR deviations, footprinting, scRNA-seq integration (label transfer), trajectory analysis (built-in + Monocle3 + Slingshot), and ArchR's interactive genome browser. Sister protocol to snapatac2 — pick this for established R workflows.
automated-bioinformatics-pipelinesPipeline maestro
autonomous-biomedical-agentAutonomous biomedical AI agent framework for executing complex research tasks across genomics, drug discovery, molecular biology, and clinical analysis. Use this skill when conducting multi-step biomedical research including CRISPR screening design, single-cell RNA-seq analysis, ADMET prediction, GWAS interpretation, rare disease diagnosis, or lab protocol optimization. Leverages LLM reasoning with code execution and integrated biomedical databases.
autonomous-oncology-agentPrecision Oncology
basecallingConvert raw Nanopore signal data (FAST5/POD5) to nucleotide sequences using Dorado basecaller. Covers model selection, GPU acceleration, modified base detection, and quality filtering. Use when processing raw Nanopore data before alignment. Note: Guppy is deprecated; use Dorado for all new analyses.
batch-processingProcess multiple sequence files in batch using Biopython. Use when working with many files, merging/splitting sequences, or automating file operations across directories.
bindingdb-databaseQuery BindingDB for measured drug-target binding affinities (Ki, Kd, IC50, EC50). Search by target (UniProt ID), compound (SMILES/name), or pathogen. Essential for drug discovery, lead optimization, polypharmacology analysis, and structure-activity relationship (SAR) studies.