Back to Case Studies

Case study · Crystal-structure molecular docking

Imatinib Docking Against ABL1 (PDB 2HYY)

Pipette prepared the ABL1 crystal structure 2HYY, generated ten imatinib poses with AutoDock Vina, and examined residue contacts with ProLIF. The top score was −11.9 kcal/mol. The crystal-pose RMSD check failed, leaving pose recovery unvalidated.

Homo sapiensPDB 2HYY · Chain AImatinibVina + ProLIF

Original 18-page session report · PDF · 1.5 MB

At a glance

A crystal structure, ten docking poses, and an explicit validation gap

Session 3b376bba used the ABL1 kinase crystal structure 2HYY, chain A, with crystallographic imatinib retained as a site reference before receptor cleaning. Vina generated ten poses; ProLIF characterized their contacts. Docking completed, but the heavy-atom RMSD comparison to the crystal ligand did not. These results describe the reported run, not a validated recovery of the experimental pose.

Receptor2HYY · Chain A
Docked poses10
Top Vina score−11.9 kcal/mol
Redocking RMSDNot computed
01 / Scientific context

Can the workflow dock imatinib into a known ABL1 binding site?

The session requested structure retrieval, receptor cleaning, pocket detection, ligand preparation, docking, and a visualization of the resulting pose. The source was the 2HYY crystal structure, reported at 2.4 Å resolution, with chain A selected from four kinase-domain copies.

What poses, Vina scores, and residue contacts does imatinib produce in the cleaned 2HYY chain A receptor, and can the workflow assess pose recovery against crystallographic imatinib?
Session3b376bba; analysis dated September 6, 2026
Receptor sourceRCSB PDB 2HYY; chain A
Ligand inputUser-supplied imatinib SMILES; 37 heavy atoms reported
Site referenceCrystallographic imatinib STI:A:600
Pocket detectionfpocket: 17 sites; pocket 1 selected
Docking and contactsAutoDock Vina 1.1.2; ProLIF 2.2.1 / MDAnalysis 2.10.0
02 / Analysis workflow

From structure retrieval to pose inspection

Retrieve and inspect 2HYYRetrieve the RCSB mmCIF and inspect chains, waters, and co-crystal ligands. The report records a gemmi conversion after an initial conversion dropped HETATM records.
Prepare chain ARetain STI:A:600 as a site reference, then remove the ligand and 51 waters from the docking receptor. The report lists 263 polymer residues and 2,073 non-hydrogen receptor atoms.
Detect and define the siteRun fpocket and compare sites with crystallographic imatinib. Pocket 1 was closest and highest-ranked; the executed Vina box came from the co-crystal ligand coordinates plus 5 Å padding.
Prepare and dock imatinibGenerate a 3D conformer with RDKit ETKDG, optimize with MMFF, and prepare PDBQT with Meeko. Run Vina with exhaustiveness 32, ten modes, seed 42, and four CPUs.
Inspect poses and validationAnalyze all ten poses with ProLIF, export visualizations, and record the failed crystal-pose RMSD check alongside the docking outputs.
Plan versus executionThe proposed plan called for an fpocket-derived box and PLIP contacts. The methods and run summary instead document a co-crystal-derived box and ProLIF analysis. This case study describes the reported execution.
Box center (Å)14.878, 16.836, 17.241
Box dimensions (Å)25.9 × 23.6 × 18.8
Box sourceSTI:A:600 coordinates plus 5 Å padding
Site coverage37/37 reference ligand atoms; 491 receptor atoms enclosed
fpocket pocket 1Score 0.894; reported volume 474 ų; 5.02 Šfrom STI centroid
03 / Docking results

Ten poses, with a top Vina score of −11.9 kcal/mol

Pose scores

Ten returned poses scored from −11.9 to −9.0 kcal/mol in this receptor and search box.

Contact analysis

ProLIF recorded 691 interaction occurrences across ten poses, including 67 for the top pose.

Shared contacts

ASP381, MET290, and PHE382 had reported hydrophobic contacts in all ten poses.

Validation status

The reference-ligand match failed; no heavy-atom redocking RMSD was obtained.

Imatinib top pose shown as green sticks in the ABL1 2HYY chain A receptor
Figure 1. Docked pose in the crystal-structure receptor. PyMOL view from report page 12: imatinib in green, receptor in slate, and selected pocket residues as sticks. This view shows the computed pose; it does not establish agreement with the crystal ligand.
AutoDock Vina scores for ten imatinib poses in ABL1 2HYY chain A, from minus 11.9 to minus 9.0 kcal per mol
Figure 2. Scores for all ten poses. Extracted from report page 15. Dashed lines show the user-provided expected range, not an experimental validation threshold. Vina scores rank poses within this protocol and are not measured affinities.
Pose rankVina score (kcal/mol)
1−11.9
2−11.1
3−10.4
4−10.0
5−10.0
6−9.9
7−9.9
8−9.3
9−9.2
10−9.0
04 / Protein–ligand contacts

Which contacts recur across the docked poses?

ProLIF contact frequencies summarize how many of the ten generated poses contain each residue–interaction pair. They describe this docking ensemble; they are not time occupancies from molecular dynamics or probabilities of binding.

ProLIF residue contact frequencies across ten imatinib docking poses in ABL1 2HYY chain A
Figure 3. Contact frequencies across ten Vina poses. Original interaction-frequency chart from report page 14. The chart distinguishes hydrophobic, van der Waals, hydrogen-bond, and cation–π contact categories; frequencies apply across the pose ensemble.
Residue (chain A)Reported interactionPoses with contact
ASP381Hydrophobic10/10
MET290Hydrophobic10/10
PHE382Hydrophobic10/10
GLU286Hydrophobic8/10
THR315VdWContact8/10
GLU286HBDonor5/10
TYR253CationPi5/10
Reading the reportThe interpretation paragraph groups GLU286 with contacts described as 100% consistent, but its hydrophobic-contact frequency is 80% in the results table and chart. The values above follow those results. A van der Waals contact is not evidence of a hydrogen bond.
05 / Limitations

The docking run completed; pose recovery remains unvalidated

  • The heavy-atom RMSD calculation failed with “No sub-structure match found.” The report attributes this to a mismatch between the prepared protonated ligand and the crystal reference; the actual pose error remains unknown.
  • The report flags the protocol as docking_protocol_unvalidated. A favorable score or a visually plausible pose does not replace a successful comparison with the experimental ligand coordinates.
  • The receptor is rigid, with crystal waters removed. Side-chain relaxation, backbone movement, and water-mediated contacts were not assessed.
  • Protonation is uncertain. The report describes a doubly protonated piperazine and an Open Babel warning about simultaneous -p/-h flags, so the preparation should not be treated as a verified pH-specific state.
  • One receptor chain and one random seed were used. Sensitivity to other chains, seeds, or preparation choices was not assessed.
  • Vina scores are exploratory ranking values. This run does not establish binding affinity, potency, selectivity, or resistance effects.

A next validation step is to reconcile ligand identity, protonation, and atom mapping before repeating the crystal-pose RMSD comparison. Preparation and seed sensitivity can then be examined with the same receptor and scoring protocol.

06 / Reproducibility and outputs

The session report and its recorded outputs

The downloadable PDF preserves session 3b376bba, including the request, plan, methods, scores, figures, limitations, and an exported reasoning excerpt. It was generated on September 10, 2026. The files below are listed in that report; this page provides the PDF and extracted figures, not a separately downloadable raw analysis bundle.

2hyy.cif / 2HYY_chainA.pdbRetrieved structure and selected receptor chain
imatinib.csvLigand SMILES input
fpocket_analysis.jsonPocket rankings and reference-site distances
docking_scores.csv / docking_box.jsonPose scores and executed search-box geometry
docking_summary.json / qc_metrics.jsonDocking settings, summary, and quality flags
redock_validation.jsonFailed RMSD validation attempt
interaction_details.csv / residue_frequency.csvPer-pose contacts and contact frequencies
pose_0000.pdbqt / pose_0000.sdfDocked pose files listed in the report
2HYY_imatinib_docked.psePyMOL session listed with the report figures
Software recorded in the reportAutoDock Vina 1.1.2; RDKit 2025.9.6; ProLIF 2.2.1; MDAnalysis 2.10.0; gemmi 0.6. The report also names fpocket, Meeko, Open Babel, and PyMOL without a complete version record for every tool.
Original Pipette session report18-page PDF · 1.5 MB · generated September 10, 2026. The HTML article is the primary case study; the PDF is the preserved session record.

Run a structure-based docking analysis

Bring a receptor and ligand question. Review the preparation, site definition, pose scores, and validation checks alongside the results.

Start in Pipette