Computer-aided methods have been broadly used in pharmaceutical research to identify potential ligands and design effective therapeutics. Most of the approaches rely on the binding affinity prediction and approximate thermodynamic properties of the system. Our alternative approach focuses on structural stability, provided by native protein-ligand interactions, in particular hydrogen bonds. Based on this idea, we designed new fast computational method, called dynamic undocking (DUck), that evaluates stability by calculating the work necessary to break the most important native contact in a ligand-receptor complex. This property is effective in distinguishing true ligands from decoys and is orthogonal to currently existing docking methods, thus making it exceptionally useful in virtual screening. Here, we present a protocol suitable for DUck's application in drug design strategy, as well as notes that will help to solve common problems addressed by users.
Keywords: Drug discovery; Hydrogen bonds; Molecular docking; Molecular dynamics; Structural stability; Structure-based drug design; Virtual screening.