{"id":6099,"date":"2020-10-06T17:32:37","date_gmt":"2020-10-06T16:32:37","guid":{"rendered":"https:\/\/www.blopig.com\/blog\/?p=6099"},"modified":"2020-10-14T17:37:31","modified_gmt":"2020-10-14T16:37:31","slug":"constrained-docking-for-bump-and-hole-methodology","status":"publish","type":"post","link":"https:\/\/www.blopig.com\/blog\/2020\/10\/constrained-docking-for-bump-and-hole-methodology\/","title":{"rendered":"Constrained docking for bump and hole methodology"},"content":{"rendered":"\n<p>Selectivity is an important trait to consider when designing small molecule probes for chemical biology. If you wish to use a small molecule to study a particular protein, but that small molecule is fairly promiscuous in its binding habits, there are risks that any effects you observe may be due to it binding other proteins with similarly shaped binding pockets, instead of your protein of interest.<\/p>\n\n\n\n<!--more-->\n\n\n\n<p>In this situation, you may decide to investigate so called \u2018bump and hole\u2019 methodology. This technique, which has been successfully employed in a variety of circumstances (e.g. <em>Science<\/em>. 2014 October 31; 346(6209): 638\u2013641. doi:10.1126\/science.1249830.), involves redesigning a small molecule probe to include extra steric bulk &#8211; a \u2018bump\u2019 &#8211; and conversely introducing a mutation into the protein of interest that swaps a bulkier residue for a smaller one &#8211; a \u2018hole\u2019. The new \u2018bump\u2019 in the small molecule will prevent it from binding properly with its original targets, but the mutated &#8211; \u2018holed\u2019 &#8211; protein will still accommodate it.<\/p>\n\n\n\n<p>Under the right circumstances, simple human inspection of a crystal structure containing the original protein and ligand combination will yield plenty of potential bump\/hole combination ideas to try. If there is a sterically bulky side chain cosying up next to portion of your molecule that looks like it could be relatively easy to chemically modify, you don\u2019t necessarily need anything more empirical before starting the relevant synthesis and mutation in the lab, then testing in vitro.<\/p>\n\n\n\n<p>But if no immediately obvious routes present themselves, taking a more methodical approach in silico in the first instance has the potential to shave many dead ends off any proposed lab work, and even to present routes that may not have been obvious by inspection alone.<\/p>\n\n\n\n<p>Assuming there is a crystal structure of the small molecule bound to the protein of interest that can be used as a basis, docking can provide a high throughput way to test different combinations of ligand modifications and protein mutations. Foldx, a software suite by CRG, can be used to evaluate proposed mutations, providing energy readouts that give a good indication of the mutant protein\u2019s stability relative to the parent. Once a 3D coordinate file for the mutant protein has been generated, the \u2018bumped\u2019 ligand can be docked into it and their predicted binding energies evaluated.<\/p>\n\n\n\n<p>\u201cBut why bother redocking a small molecule whose binding conformation is, essentially, already known?\u201d I hear you ask. A very good point, and this is why constrained docking is so useful for this particular scenario. GOLD, by CCDC, allows users to provide a \u2018scaffold\u2019 when docking, which will bias its docking results towards those where a certain portion of your molecule is held in a user-defined position. In this case, a conserved portion of the ligand &#8211; i.e., a section that is exactly the same in both the original ligand and the bumped ligand &#8211; can be used to hold that section of the bumped ligand in the correct position and conformation during docking. This ensures that the computational power is used effectively, by exploring conformations of the modified section of the ligand, rather than trying to accurately redock a section whose binding conformation is already known to us.<\/p>\n\n\n\n<p>With this in mind, I find it most effective to build any \u2018bumped\u2019 ligand 3D coordinate files using the ligand from the original crystal structure as a template &#8211; e.g. using PyMOL\u2019s \u2018builder\u2019 function. This will not only preserve the coordinates of the conserved section (which, in theory, doesn\u2019t matter), but it will also preserve any ring conformations present (which docking software often has more trouble with). Contrast this with, for example, converting a SMILES string into a .mol2 files using RDKit, Openbabel, or otherwise: the conformation of the conserved section given by this may well be different to that seen in the crystal structure, and again it makes no sense the the docking software to attempt to recreate something for which we already have the answer.<\/p>\n\n\n\n<p>Another feature that GOLD offers is residue flexibility. How extensively you use this feature &#8211; if at all &#8211; will depend on your specific system. But allowing the mutated residue\u2019s side chain to explore multiple conformations along with the unconstrained half of the ligand will provide more insight than leaving it in the \u2018default\u2019 position assigned by Foldx.<\/p>\n\n\n\n<p>Once GOLD has returned its best docking attempts for a set of \u2018bumped\u2019 ligands and \u2018holed\u2019 proteins, Any suitable forcefield &#8211; such as the one built into Autodock 4 &#8211; can be used to predict their binding energies, and inspection of the full grid of ligands vs proteins will highlight any exciting avenues to be explored in the lab. From my own personal experience, I would encourage anyone attempting this method to be exhaustive in their docking, pitching all twenty potential mutants for any given position against a set of all the chemical modifications you can think of\/easily make 3D coordinate files for. Don\u2019t limit your testing to just testing the bumps and holes that you think will complement each other, but do the full exhaustive grid, as some predicted interactions with good binding energies may well end up being between combinations that wouldn\u2019t have appeared intuitive.<\/p>\n\n\n\n<p>Furthermore, whilst taking synthetic feasibility into account when proposing chemical modifications is obviously sensible, don\u2019t be frightened of testing ideas in silico that might be hard to make in the lab &#8211; after all, it won\u2019t add much extra time to your computational workflow. Besides, spending time making the right compound from the get go, however tricky the synthesis, may end up proving quicker than sticking to easier syntheses, but having to do loads because none of your compounds work properly.<\/p>\n","protected":false},"excerpt":{"rendered":"<p>Selectivity is an important trait to consider when designing small molecule probes for chemical biology. If you wish to use a small molecule to study a particular protein, but that small molecule is fairly promiscuous in its binding habits, there are risks that any effects you observe may be due to it binding other proteins [&hellip;]<\/p>\n","protected":false},"author":76,"featured_media":0,"comment_status":"closed","ping_status":"open","sticky":false,"template":"","format":"standard","meta":{"nf_dc_page":"","wikipediapreview_detectlinks":true,"_monsterinsights_skip_tracking":false,"_monsterinsights_sitenote_active":false,"_monsterinsights_sitenote_note":"","_monsterinsights_sitenote_category":0,"ngg_post_thumbnail":0,"_jetpack_memberships_contains_paid_content":false,"footnotes":""},"categories":[228,291],"tags":[339,335,22,338,336,340,287],"ppma_author":[551],"class_list":["post-6099","post","type-post","status-publish","format-standard","hentry","category-protein-structure","category-protein-ligand-docking","tag-bromodomains","tag-bump-and-hole-methodology","tag-docking","tag-foldx","tag-gold","tag-imids","tag-protein-ligand-docking"],"jetpack_featured_media_url":"","jetpack_sharing_enabled":true,"authors":[{"term_id":551,"user_id":76,"is_guest":0,"slug":"patrick","display_name":"Patrick Brennan","avatar_url":"https:\/\/secure.gravatar.com\/avatar\/3fce00c610506084a8ff25e90b2dcd13a78befd877255c5128c69e358d8962eb?s=96&d=mm&r=g","0":null,"1":"","2":"","3":"","4":"","5":"","6":"","7":"","8":""}],"_links":{"self":[{"href":"https:\/\/www.blopig.com\/blog\/wp-json\/wp\/v2\/posts\/6099","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/www.blopig.com\/blog\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/www.blopig.com\/blog\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/www.blopig.com\/blog\/wp-json\/wp\/v2\/users\/76"}],"replies":[{"embeddable":true,"href":"https:\/\/www.blopig.com\/blog\/wp-json\/wp\/v2\/comments?post=6099"}],"version-history":[{"count":2,"href":"https:\/\/www.blopig.com\/blog\/wp-json\/wp\/v2\/posts\/6099\/revisions"}],"predecessor-version":[{"id":6140,"href":"https:\/\/www.blopig.com\/blog\/wp-json\/wp\/v2\/posts\/6099\/revisions\/6140"}],"wp:attachment":[{"href":"https:\/\/www.blopig.com\/blog\/wp-json\/wp\/v2\/media?parent=6099"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/www.blopig.com\/blog\/wp-json\/wp\/v2\/categories?post=6099"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/www.blopig.com\/blog\/wp-json\/wp\/v2\/tags?post=6099"},{"taxonomy":"author","embeddable":true,"href":"https:\/\/www.blopig.com\/blog\/wp-json\/wp\/v2\/ppma_author?post=6099"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}