An all-in-one workflow for pose prediction and scoring of covalently bound ligands
The Advantages of Covalent Docking
With the recent resurgence in covalent drug research, computational insight into covalent docking is becoming key to understanding how covalent inhibitors can be used to address selectivity and potency challenges.
Covalent inhibitors derive their activity not only from the formation of a covalent bond between the target and the ligand but also from stabilizing non-covalent forces in the binding pocket. CovDock selects the top covalent complexes using the extensively validated Prime energy model, and calculates an apparent affinity score that captures these essential elements of a successful covalent docking process:
1. The pre-reactive ligand form occupies the binding pocket with enough residency time to facilitate the reaction of the ligand warhead with the reactive protein residue; and
2. unfavorable steric clashes and poor electrostatic contacts are prevented as the reaction proceeds.
CovDock begins with Glide docking to a receptor with the reactive residue trimmed to alanine. The receptor reactive residue is then added and sampled to form a covalent bond with the ligand in different poses. Covalent complexes are minimized using the Prime VSGB2.0 energy model to score the top covalent complexes. An apparent affinity score, based on the Glide score of pre-reactive and post-reactive poses, is also calculated to estimate binding energies for use in virtual screening.
Accurate binding mode prediction:
CovDock is built upon a foundation of the time-tested Glide docking algorithm and Prime structure refinement methodology for accurate prediction of non-covalently docked poses. Glide quickly samples a large pool of initial poses for the pre-reactive species and Prime simultaneously optimizes the ligand pose and attachment residue to produce a sound physical chemistry. The resultant accuracy outperforms other docking programs in achieving lower RMS deviations from native co-crystallized structures.
CovDock performs a series of automated steps based on a simple setup from the Maestro graphical interface or from the command line. First, CovDock docks the pre-reactive ligand to determine viable poses that bring the reactive group into close proximity with the reactive receptor residue. Then the covalent bond is formed for the top scoring complex structures, the covalently attached ligand is sampled, and the complexes are scored using all-atom molecular mechanics with the OPLS force field and VSGB2.0 implicit solvent model.
Intuitive graphical interface:
Schrödinger’s intuitive graphical user interface, Maestro, provides easy-to-use panels for straightforward set-up of experiments, easy visualization, and efficient analysis of CovDock results.
Citations and Acknowledgements
ö Zhu, K.; Borrelli, K.W.; Greenwood, J.R.; Day, T.; Abel, R.; Farid, R.S.; Harder, E., "Docking covalent inhibitors: A parameter free approach to pose prediction and scoring," J. Chem. Inf. Model., 2014, 54, 1932−1940
ö "A structure-based virtual screening approach for discovery of covalently bound ligands"Toledo Warshaviak, D.; Golan, G.; Borrelli, K.W.; Zhu, K.; Kalid, O., J. Chem. Inf. Model, 2014, 54(7), 1941–1950
ö "Docking covalent inhibitors: A parameter free approach to pose prediction and scoring"Zhu, K.; Borrelli, K.W.; Greenwood, J.R.; Day, T.; Abel, R.; Farid, R.S.; Harder, E., J. Chem. Inf. Model., 2014, 54, 1932−1940
Schrödinger has made available several custom reactions that can be used in CovDock studies. To import any of the following files into Maestro, simply open the Covalent Docking Panel, click on the Reaction Type tab, select Custom from the drop-down menu, and then use the Browse option to load the reaction of interest.
|Custom Covalent Reactions|
|Index||Reaction||Ligand SMARTS||Receptor SMARTS||CDOCK File|
The following reaction types come pre-installed with CovDock and are shown here for reference.
|Reaction||Ligand SMARTS||Receptor Residue|
|Michael Addition||[C,c]=[C,c]-[C,c,S,s]=[O]||Cys, Ser|
|Nucleophilic Addition to a Double Bond||[C,c]=[O,S]||Cys, Ser|
|Nucleophilic Addition to a Triple Bond||[C]#[N]||Cys, Ser, Tyr|
|Nucleophilic Substitution||[F,Cl,Br,I]||Cys, Ser|
|Boronic Acid Addition||[B]([O])[O]||Cys, Ser|
|Epoxide Opening||[C;r3][O;r3][C;r3]||Cys, Ser, His|
|Phosphonate Addition||[P]-[O;H1,-1]||Cys, Ser|
|Beta Lactam Addition||[O-0X1]=[C]1[C][C][N]1||Cys, Ser|
|Conjugate Addition to Alkene (nitrile activated)||[C,c]=[C,c]-[C,c]#[N,n]||Cys, Ser|
|Conjugate Addition to Alkyne (carbonyl activated)||[C-0X2]#[C-0X2][C-0X3]=[O-0X1]||Cys, Ser|
|Conjugate Addition to Alkyne (aryl activated)||[C]#[C]-[c]||Cys, Ser|