Virtual Screening

Dock a compound library to find hit candidates

Give a protein structure and a list of compounds, and each compound is docked into the pocket and scored. You get a ranked table you can filter by score, drug-likeness and chemical family, so you know which compounds to buy before you run an experiment.

GNINA 1.3.1 with CNN rescoring · Doc updated Oct 14, 2026

Input
Structure + compounds
Output
Poses, scores, clusters
Engine
GNINA
Cost
1 cr per compound
Runs on
Up to 200 GPUs at once

When to use

  • You have a target with a known structure and want first hit candidates
  • You want to decide which compounds from a purchasable library to buy first
  • You want to shrink a long list before slower steps such as AI-Cofolding or activity prediction

No structure for your target? Predict one with AI-Cofolding first. Docking scores are not measured binding strengths. Use them to rank compounds, not to promise potency.

How it works

StructureCompounds01Prepare receptorBound ligandremovedWater, altLocremoved02Set the pocketDocking box aroundthe reference ligand03Dock and rescoreCNN scorePoses per compound,rescored by CNN04Rank and group#ScoreCluster12345
  1. 01Prepare the receptor

    Keeps the first model and the chains you chose, removes alternate atom positions (altLoc), water and the bound ligand, and keeps only the ions and cofactors you ticked.

  2. 02Set the pocket

    A docking box is drawn around a reference ligand, from the structure or from a file. You can also give the box center and size yourself.

  3. 03Dock and rescore

    Each compound gets 3D coordinates with RDKit, then GNINA searches poses in the box and a CNN rescores each pose for how plausible it is and how strongly it binds.

  4. 04Filter and group

    Compounds are ranked, drug-likeness properties are computed, and similar structures can be grouped into clusters so you compare chemical families.

Input formats

ItemAcceptedNotes
ReceptorPDB ID, or a .pdb, .pdbqt or .cif upload up to 50 MBPick the chains to keep. Water is removed automatically
PocketA reference ligand (from the structure, or an .sdf, .mol2 or .pdb file), or a box center and sizeManual boxes start at 20 Å per side and can be 1–100 Å
CompoundsPre-built libraries, or your own .smi, .txt, .csv or .tsv file up to 100 MBPre-built: Enamine sets (Hit Locator, Discovery Diversity, PPI, CNS, Covalent, molecular glues) and PharmBricks. Combine up to 20; duplicates are removed
CountUp to 100,000 compounds per jobCompounds with more than 70 rotatable bonds are rejected before you are charged

Pre-built libraries can be narrowed before docking by molecular weight, cLogP, TPSA, hydrogen bond donors and acceptors, and QED.

Tutorial: dock known EGFR inhibitors into the erlotinib pocket

PDB 1M17 is the kinase domain of EGFR, a lung cancer target, with the drug erlotinib bound (ligand code AQ4). We use erlotinib's site as the pocket and dock 40 EGFR inhibitors with measured potency into it, to see whether docking ranks the potent ones near the top. Run example at the top of this page fills all of this in for you.

  1. Load the structure

    In Target Protein, enter 1M17 as the PDB ID. The structure appears in the viewer on the right.

    Target Protein step with PDB ID 1M17 entered and the Mol* preview
  2. Pick the chain and the pocket

    Keep chain A. Under binding site, choose Auto Box and pick AQ4 as the reference ligand. The box is drawn around erlotinib.

    Chain A selected, Auto Box as the binding site type, AQ4 chosen as reference ligand
  3. Add the compounds

    In Ligand Library, choose Upload and add egfr-example-ligands.txt: 40 EGFR inhibitors from ChEMBL, one ID and SMILES per line. For a real screen, pick a pre-built library instead. The estimated credits update with the compound count.

    Ligand Library step with the 40-compound example file uploaded and the credit estimate
  4. Check and submit

    In Project Settings, name the project, check the credits in the summary (40 compounds cost the 50 cr minimum) and submit. We email you when it finishes.

  5. Filter and group the results

    On the result page, apply a filter preset such as Lipinski Ro5, run clustering, and open a compound to see its pose in 3D.

    Result page with filters, the ranked table and the clustering panel
  6. Take it further

    Export the compounds you like as CSV or SDF poses, open Enamine compounds in the Enamine store, or fold the best ones with AI-Cofolding.

Reading the results

ColumnMeaningHow to read it
TNCOCombined score of CNN affinity, CNN score and ΔG, used to pick each compound's best pose and the default orderHigher is better. It is normalized over one compound's own poses, so also check CNN Aff. and ΔG when comparing compounds
CNNCNN pose score (0–1): how likely the pose is a real binding poseHigher means a more trustworthy pose
CNN Aff. (pKd)CNN-predicted binding strengthHigher means tighter predicted binding. 1 unit is a tenfold change
ΔG (kcal/mol)Physics-based docking scoreLower (more negative) is better
MW, cLogP, QED, SAMolecular weight, lipophilicity, drug-likeness (0–1) and synthetic accessibilityUse the filter presets (Lipinski, lead-like, fragment, drug-like)
ClusterGroup of structurally similar compoundsLook at a few per family. The cluster's top-ranked compound is its representative

Clustering uses Morgan fingerprints (radius 2, 512 bits) and Tanimoto similarity with average-linkage hierarchical clustering over the top 5,000 compounds by default.

Example result: EGFR 1M17

A real job, shown read-only. Erlotinib itself (CHEMBL553) ranks second of the 40 compounds, and the viewer shows its docked pose in the pocket it occupies in the crystal structure.

CompoundTNCOCNNCNN Aff.ΔG
CHEMBL2964070.9300.8535.13-6.5
CHEMBL553 (erlotinib)0.8610.8386.88-7.2
CHEMBL1626220.8570.7446.52-7.5
CHEMBL248280.8530.5486.71-7.8
CHEMBL1266230.8000.7526.33-6.6

Cost and time

CompoundsCreditsPrice
1,0001,000 cr$10
10,00010,000 cr$100
100,000100,000 cr$1,000

One compound costs 1 cr ($0.01), with a minimum of 50 cr per job. Compounds are split into chunks docked in parallel on up to 200 GPUs, so up to 10,000 compounds run in a single wave and take about as long as a few hundred. The EGFR example (40 compounds) took about 12 minutes of docking once a GPU was assigned. The wizard shows the credits before you submit, and the result page shows the remaining time while it runs.

If the job fails or finishes with no results, the credits come back in full. If you stop it yourself, half comes back.

Citations

If you use these results in a report or paper, please cite the methods below.

  • McNutt AT, Francoeur P, Aggarwal R, Masuda T, Meli R, Ragoza M, Sunseri J, Koes DR. GNINA 1.0: molecular docking with deep learning. J Cheminform 13, 43 (2021). doi:10.1186/s13321-021-00522-2Docking and CNN rescoring
  • Ragoza M, Hochuli J, Idrobo E, Sunseri J, Koes DR. Protein-ligand scoring with convolutional neural networks. J Chem Inf Model 57, 942-957 (2017). doi:10.1021/acs.jcim.6b00740CNN scoring model
  • Trott O, Olson AJ. AutoDock Vina: improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading. J Comput Chem 31, 455-461 (2010). doi:10.1002/jcc.21334Vina scoring function behind ΔG
  • RDKit contributors. RDKit: Open-source cheminformatics. Zenodo (2026). doi:10.5281/zenodo.5916373D coordinates, properties and fingerprints
  • Rogers D, Hahn M. Extended-connectivity fingerprints. J Chem Inf Model 50, 742-754 (2010). doi:10.1021/ci100050tMorgan (ECFP) fingerprints for clustering
  • Stamos J, Sliwkowski MX, Eigenbrot C. Structure of the epidermal growth factor receptor kinase domain alone and in complex with a 4-anilinoquinazoline inhibitor. J Biol Chem 277, 46265-46272 (2002). doi:10.1074/jbc.M207135200EGFR–erlotinib structure used in the tutorial (PDB 1M17)

FAQ