Module 6: Molecular Docking
Understand how computer algorithms predict ligand binding configurations inside a receptor pocket. Explore conformational searches, scoring functions, and validation strategies.
1. What is Molecular Docking?
Molecular docking is a computational simulation technique that models the interaction between a small molecule (ligand) and a macromolecular target (typically a protein).
The docking process has two primary components:
- Search Algorithm: Explores the conformational space of the ligand inside the active site.
- Scoring Function: Evaluates the binding energy of each generated pose to identify the most stable orientation.
2. Conformational Search Methods
Systematic Search
Exhaustively rotates every rotatable bond by increments (e.g. 60°). It guarantees finding global minima but suffers from a combinatorial explosion.
Stochastic Search
Performs random conformational moves. Prominent algorithms include Monte Carlo and Genetic Algorithms, where coordinates mutate and cross over.
Incremental Construction
Splits the ligand into fragments. The anchor fragment is docked first, and the remaining segments are attached step-by-step.
Interactive Playground: Manual Pose Matcher
Use the translation (X, Y) and rotation sliders to dock the ligand into the active site pocket. Align the positive charge (+) with the negative charge (-) to optimize the binding score.
-8.28kcal/mol
3. Target Coordinate Selection: Navigating the PDB
Before running a docking simulation, researchers must select a protein target coordinate set. The Worldwide Protein Data Bank (PDB) houses hundreds of thousands of structures, but they are not created equal. Poor target coordinates will invalidate all downstream scoring.
Checklist for Systematic PDB Selection:
- Sequence Fidelity via UniProt: Search by UniProt accession ID, not common protein names, to filter out homologous species variants and ensure exact sequence identity.
- Experimental Method & Resolution: X-ray crystal structures are the default choice for docking as they represent frozen energy snapshots with defined coordinate densities. NMR provides structural ensembles that require complex ensemble docking. Resolution should ideally be below 2.5 Å.
- Refinement Metrics:Inspect the structure's R-free (relative to R-work) and Ramachandran outlier percentages.
- Engineered Mutations: Crystallographers frequently introduce mutations to prompt crystallization. Ensure the active site loop residues are wild-type.
- Loop Completeness: Weak electron densities in flexible loops lead crystallographers to omit loops entirely. If gaps occur in the active site, loops must be modeled in using loop refinement tools before docking.
- Biological Assembly vs. Asymmetric Unit: Download the biological assembly (the functional form under physiological conditions), not the asymmetric unit (a mathematical crystallographic convenience). Extracting a monomer from a functional dimer exposes hydrophobic interfaces to solvent, causing artificial coordinate collapses.
The legacy .pdb file format (designed in the 1970s for punch cards) suffers from a hard limit of 99,999 atoms and 62 chains. In the era of cryo-EM macromolecular assemblies, the format breaks down. The PDB has formally deprecated it in favor of PDBx/mmCIF (.cif). By mid-2027, new depositions will only receive 12-character IDs (e.g. pdb_00002v5z) and .pdb formats will no longer be generated. Modern pipelines should migrate to .cif inputs exclusively.
When Experimental Structures Are Unavailable
Docking may begin from a predicted or comparative receptor model, but generating and validating that model is a structural-bioinformatics task. This module therefore focuses on the narrower question: does the resulting pocket support the proposed docking experiment?
Module 13: build and validate the receptor model
Covers structure selection, comparative modeling, local quality, receptor chemistry, PDB/mmCIF interpretation, and preparation provenance.
Sequence supportModule 18: find and assess templates
Covers BLAST, alignment coverage, conserved residues, gaps, profiles, and evolutionary context used to justify a comparative-model template.
Docking-specific acceptance gate
- Inspect confidence and completeness locally at the pocket, not only across the global fold.
- Challenge protonation, side-chain rotamers, cofactors, metals, waters, and alternative pocket states.
- Use an ensemble when plausible receptor conformations materially change the site.
- Where reference ligands exist, test pose recovery and screening enrichment before prospective use.
A well-formed model is not automatically a useful docking receptor. Task-specific retrospective performance is the final acceptance test.
4. Ligand Coordinate Preparation: Electrostatic Charge Selection
To compute electrostatic interactions (like hydrogen bonds, salt bridges, and dipole-dipole contacts), docking scoring functions rely on ligand partial atomic charges. The choice of charge model propagates directly into the final binding scores.
While full ab initio quantum mechanics methods (like RESP) are highly accurate, they are computationally prohibitive for virtual screening libraries. Semi-empirical charge models approximate this quality at a fraction of the cost:
| Charge Model | Electrostatic Quality | Calculation Speed | Recommended Use Case |
|---|---|---|---|
| AM1-BCC | Very High (errors < 0.1 e) | Fast | Default choice for general protein-ligand docking pipelines. |
| PM6 | High (improves H-bonding) | Medium | Challenging systems requiring precise hydrogen-bonding parameters. |
| RESP | Excellent (Gold Standard) | Very Slow | Lead optimization and final pose verification. Too slow for screening. |
| Gasteiger | Poor | Extremely Fast | Quick initial pre-filtering or triage of massive compound libraries. |
The PDBQT Format: Charges and Atom Types
Standard PDB files from the Protein Data Bank contain only element symbols and Cartesian coordinates. For docking engines like AutoDock Vina to parse steric and electrostatic interactions, files must be pre-processed into the PDBQT format:
Partial Charge (Q)
Appends a column representing the partial atomic charge (e.g., in electrons, calculated using AM1-BCC or Gasteiger models) to compute electrostatic potentials.
AutoDock Atom Type (T)
Maps atoms to specific docking categories (e.g., HD for donor hydrogen, A for aromatic carbon, OA for acceptor oxygen) to look up force field parameters.
Torsional Topology
Defines the molecule as a tree with ROOT and BRANCH statements, specifying which bonds are frozen and which can rotate.
5. Types of Scoring Functions
Force-Field-Based Functions
Calculate binding energy from electrostatics (Coulomb) and van der Waals interactions, using mechanical parameters. Examples: DOCK, GOLDScore.
Empirical Scoring Functions
Sum of weighted terms (hydrogen bonds, contacts, rotatable bonds) parameterized using experimental binding affinities. Examples: ChemScore, GlideScore.
Knowledge-Based Functions
Derived from statistical distributions of atom-atom contact pairs in structural databases. Favorable shapes yield lower potentials. Examples: PMF, DrugScore.
Empirical Scoring Function Formulation:
Empirical scoring models (like the one used in AutoDock Vina) estimate the free energy of binding by summing distinct, weighted interaction contributions:
Where each contribution corresponds to a specific physical force:
- vdW (van der Waals): Steric contacts modeled by attractive/repulsive Lennard-Jones terms based on interatomic distance (r_ij).
- hbond (Hydrogen Bonding): Directional and distance-dependent interaction scoring for donor-acceptor atom pairs.
- electro (Electrostatics): Coulombic attraction or repulsion calculated using partial charges and distance-dependent dielectric shielding.
- desolv (Desolvation): Energy penalty for stripping away water molecules from polar groups as they enter hydrophobic pockets.
- w_rotor × N_rot: Entropic penalty proportional to the number of rotatable bonds (N_rot) frozen upon binding, weighted by w_rotor.
Docking Software Comparison
Different docking programs employ distinct algorithms and scoring paradigms. The choice of software depends on the target system, licensing constraints, and screening scale.
| Software | Algorithm | Scoring | License | Best For |
|---|---|---|---|---|
| AutoDock Vina | Iterated Local Search | Empirical hybrid | Open-source | Academic VS, large screens |
| Glide | Hierarchical HTVS/SP/XP | GlideScore | Commercial | Pharma lead optimization |
| GOLD | Genetic Algorithm | ChemScore, GoldScore | Commercial (academic free) | Metalloenzymes, flexible |
| FlexX | Incremental Construction | Bohm empirical | Commercial | Fragment-based |
| rDock | Cavity-based + GA | Empirical + vdW | Open-source | RNA-ligand docking |
Practical Docking Protocol
A systematic end-to-end docking workflow. Each step builds upon the previous one to ensure reliable, reproducible binding predictions.
Target Preparation
Download the PDB structure, remove crystallographic waters (except conserved ones mediating key interactions), assign protonation states using PropKa or H++ at physiological pH, and add missing hydrogens.
Binding Site Definition
Define a grid box centered on the co-crystallized ligand with 0.375 Å spacing and a 10–15 Å buffer around the binding pocket to allow full ligand exploration.
Ligand Preparation
Generate 3D coordinates from SMILES, assign AM1-BCC partial charges, enumerate relevant tautomers at pH 7.4, and generate diverse low-energy conformers.
Docking Run
Set the exhaustiveness parameter to control thoroughness: 8 for fast screening, 32 for thorough sampling, and 128 for exhaustive poses. Generate 9–20 poses per ligand for adequate conformational coverage.
Post-Processing
Rank poses by docking score, filter candidates by interaction fingerprint similarity, visually inspect the top 20 hits for chemically sensible binding modes, and apply consensus scoring across multiple scoring functions to improve reliability.
6. Covalent Docking: Modeling Chemical Warheads
Traditional small-molecule drugs bind via reversible, non-covalent interactions. However, a major class of modern therapeutics consists of covalent inhibitors(e.g. Sotorasib targeting KRAS G12C, or Ibrutinib targeting Bruton's Tyrosine Kinase). These molecules contain a mildly electrophilic "warhead" that forms a permanent chemical bond with a nucleophilic amino acid residue in the binding pocket.
Because standard docking scoring functions are parameterized only for non-covalent forces (electrostatics, van der Waals), simulating covalent binding requires specialized covalent docking protocols:
1. Non-Covalent Fit
The ligand first docks non-covalently. The scoring function evaluates whether the scaffolds fit the pocket topography and positions the reactive warhead (e.g. acrylamide, sulfonyl fluoride) within reaction distance (< 3.5 Å) of the target residue.
2. Chemical Linkage
The algorithm forms a virtual covalent bond between the nucleophile (typically Cys, Ser, or Lys) and the electrophilic carbon. The system conformational search resolves the dihedral angles of the newly formed linkage.
3. Complex Minimization
The bound ligand-protein complex is energy-minimized. Scoring functions evaluate the local strain of the covalent linker alongside the remaining non-covalent interactions of the scaffold to rank candidate covalent hits.
7. Ternary Complex Docking: PROTACs & Molecular Glues
Covalent docking bent the rules by adding a bond. Targeted Protein Degradation (TPD) breaks a deeper assumption: that docking means fitting one ligand into one pocket. Here you must dock a three-body complex.
As Module 2 showed, most disease-driving proteins have no druggable pocket at all. TPD sidesteps the problem entirely: rather than inhibiting the target, it recruits the cell's own Ubiquitin-Proteasome System to destroy it. You no longer need a deep pocket — you only need a grip.
PROTACs vs. Molecular Glues
Cooperativity (α) — why ternary docking is not just docking twice
The two binding events are not independent. Cooperativity relates the ternary dissociation constant (binding of the second protein to the binary complex) to the binary one (binding to the isolated protein):
α > 1 is positive cooperativity: favourable protein-protein contacts induced by the PROTAC make the ternary complex more stable than either binary interaction predicts. α < 1 is negative cooperativity — the two protein surfaces clash. This is why a PROTAC built from two excellent binders can still fail completely: the complex, not the compound, is the drug.
Ternary docking breaks assumptions from three different modules at once. The search problem explodes: you sample two rigid-body placements plus a floppy linker with many rotatable bonds (Module 4). The scoring problem becomes protein-protein, not protein-ligand. And PROTACs sit far outside Lipinski space — routinely 800–1,100 Da — so the Rule of Five you meet in Module 12 simply does not apply. TPD is the clearest demonstration in this course that the standard methods are approximations with edges.
Interactive Playground: PROTAC Ternary Assembly & the Hook Effect
Tune the linker, the cooperativity, and the dose. Then run the experiment that breaks everyone's intuition: set a good linker, then keep raising the concentration. Degradation collapses. More drug makes it work worse — the hook effect.
α > 1 positive · α < 1 negative (surfaces clash)
99.4%
Optimal Ternary Assembly: The linker distance and concentration are balanced, forming a stable Target-PROTAC-E3 complex that actively degrades target proteins.
At low dose there is too little PROTAC to bridge both proteins. At the optimum, each molecule links a target to an E3. But push the concentration high enough and every protein gets its own PROTAC molecule — target-PROTAC binaries and E3-PROTAC binaries — and none are bridged. The productive ternary complex is competed away by the drug itself. This is why a PROTAC dose-response is a bell, not a sigmoid, and why "more drug" is not a valid strategy for a degrader.
8. Redocking Validation (RMSD)
Before trust is placed in docking scores, protocols must be validated. The standard validation method is self-docking (redocking):
Root Mean Square Deviation (RMSD): The measure of spatial distance between atoms of the predicted docking pose and the actual co-crystallized native structure in the crystal. A docking protocol is considered valid if the redocked ligand reproduces the experimental structure with an RMSD ≤ 2.0 Å.