{"id":936,"date":"2019-04-18T19:10:40","date_gmt":"2019-04-18T19:10:40","guid":{"rendered":"https:\/\/ccsb.scripps.edu\/adfr\/?page_id=936"},"modified":"2019-04-18T21:21:06","modified_gmt":"2019-04-18T21:21:06","slug":"tutorial-covalent","status":"publish","type":"page","link":"https:\/\/ccsb.scripps.edu\/adfr\/tutorial-covalent\/","title":{"rendered":"tutorial covalent"},"content":{"rendered":"<div class=\"fusion-fullwidth fullwidth-box fusion-builder-row-1 nonhundred-percent-fullwidth non-hundred-percent-height-scrolling\" style=\"--awb-border-radius-top-left:0px;--awb-border-radius-top-right:0px;--awb-border-radius-bottom-right:0px;--awb-border-radius-bottom-left:0px;--awb-flex-wrap:wrap;\" ><div class=\"fusion-builder-row fusion-row\"><div class=\"fusion-layout-column fusion_builder_column fusion-builder-column-0 fusion_builder_column_1_1 1_1 fusion-one-full fusion-column-first fusion-column-last\" style=\"--awb-bg-size:cover;\"><div class=\"fusion-column-wrapper fusion-flex-column-wrapper-legacy\"><div class=\"fusion-text fusion-text-1\"><h3>Docking a covalently bound ligand<\/h3>\n<p>In this tutorial we illustrate how to re-dock the native covalent ligand of PDB id 3c9w.\u00a0Both, the receptor (3cw9.pdbqt) and ligand (3c9w_ligandWithSideChain_<span class=\"GramE\">random.pdbqt<\/span>) are available in the data file associate with this tutorial.<\/p>\n<p>For covalent docking, the receptor and covalent ligand need to share 3 atoms as sown below. 2 atoms (green) for the covalent bond and 1 atom (yellow) is an anchor atom on the receptor side. These atoms are used to transform the ligand into place with respect to the receptor to create the covalent bond.<\/p>\n<p>NOTE: for covalent docking no translational points are needed as the ligand is positioned in the docking box using the covalent bond.<\/p>\n<p><img decoding=\"async\" class=\" wp-image-938 alignnone\" src=\"https:\/\/ccsb.scripps.edu\/adfr\/wp-content\/uploads\/sites\/36\/2019\/04\/covalentFig.png\" alt=\"\" width=\"582\" height=\"265\" srcset=\"https:\/\/ccsb.scripps.edu\/adfr\/wp-content\/uploads\/sites\/36\/2019\/04\/covalentFig-200x91.png 200w, https:\/\/ccsb.scripps.edu\/adfr\/wp-content\/uploads\/sites\/36\/2019\/04\/covalentFig-300x137.png 300w, https:\/\/ccsb.scripps.edu\/adfr\/wp-content\/uploads\/sites\/36\/2019\/04\/covalentFig-400x182.png 400w, https:\/\/ccsb.scripps.edu\/adfr\/wp-content\/uploads\/sites\/36\/2019\/04\/covalentFig.png 468w\" sizes=\"(max-width: 582px) 100vw, 582px\" \/><\/p>\n<p><strong>In this tutorial you will learn<\/strong>:<\/p>\n<ul>\n<li>to generate a target file for covalent docking<\/li>\n<li>to run <em>ADFR<\/em> to dock the ligand<\/li>\n<li>to understand the output of an <em>ADFR<\/em> docking run<\/li>\n<\/ul>\n<\/div><div class=\"fusion-sep-clear\"><\/div><div class=\"fusion-separator fusion-full-width-sep\" style=\"margin-left: auto;margin-right: auto;margin-top:15px;margin-bottom:25px;width:100%;\"><div class=\"fusion-separator-border sep-shadow\" style=\"--awb-height:20px;--awb-amount:20px;background:radial-gradient(ellipse at 50% -50% , #e0dede 0px, rgba(255, 255, 255, 0) 80%) repeat scroll 0 0 rgba(0, 0, 0, 0);background:-webkit-radial-gradient(ellipse at 50% -50% , #e0dede 0px, rgba(255, 255, 255, 0) 80%) repeat scroll 0 0 rgba(0, 0, 0, 0);background:-moz-radial-gradient(ellipse at 50% -50% , #e0dede 0px, rgba(255, 255, 255, 0) 80%) repeat scroll 0 0 rgba(0, 0, 0, 0);background:-o-radial-gradient(ellipse at 50% -50% , #e0dede 0px, rgba(255, 255, 255, 0) 80%) repeat scroll 0 0 rgba(0, 0, 0, 0);\"><\/div><\/div><div class=\"fusion-sep-clear\"><\/div><div class=\"fusion-tabs fusion-tabs-1 clean vertical-tabs icon-position-left mobile-mode-accordion\" style=\"--awb-title-border-radius-top-left:0px;--awb-title-border-radius-top-right:0px;--awb-title-border-radius-bottom-right:0px;--awb-title-border-radius-bottom-left:0px;--awb-inactive-color:#9acae0;--awb-background-color:#ffffff;--awb-border-color:#9acae0;--awb-active-border-color:#d65169;\"><div class=\"nav\"><ul class=\"nav-tabs\" role=\"tablist\"><li class=\"active\" role=\"presentation\"><a class=\"tab-link\" data-toggle=\"tab\" role=\"tab\" aria-controls=\"tab-74538996c2b940c7d33\" aria-selected=\"true\" id=\"fusion-tab-generatetargetfile\" href=\"#tab-74538996c2b940c7d33\"><h4 class=\"fusion-tab-heading\"><i class=\"fontawesome-icon fa-cogs fas\" aria-hidden=\"true\" style=\"font-size:13px;\"><\/i>Generate target file<\/h4><\/a><\/li><li role=\"presentation\"><a class=\"tab-link\" data-toggle=\"tab\" role=\"tab\" aria-controls=\"tab-00c3697ba7ca4dd7edf\" aria-selected=\"false\" tabindex=\"-1\" id=\"fusion-tab-inspecttargetfile\" href=\"#tab-00c3697ba7ca4dd7edf\"><h4 class=\"fusion-tab-heading\"><i class=\"fontawesome-icon fa-search fas\" aria-hidden=\"true\" style=\"font-size:13px;\"><\/i>Inspect target file<\/h4><\/a><\/li><li role=\"presentation\"><a class=\"tab-link\" data-toggle=\"tab\" role=\"tab\" aria-controls=\"tab-4d64c9c233eda02d1a1\" aria-selected=\"false\" tabindex=\"-1\" id=\"fusion-tab-dockingligand\" href=\"#tab-4d64c9c233eda02d1a1\"><h4 class=\"fusion-tab-heading\"><i class=\"fontawesome-icon fa-cogs fas\" aria-hidden=\"true\" style=\"font-size:13px;\"><\/i>Docking ligand<\/h4><\/a><\/li><li role=\"presentation\"><a class=\"tab-link\" data-toggle=\"tab\" role=\"tab\" aria-controls=\"tab-630b551d9b728b54dea\" aria-selected=\"false\" tabindex=\"-1\" id=\"fusion-tab-understandingtheoutput\" href=\"#tab-630b551d9b728b54dea\"><h4 class=\"fusion-tab-heading\"><i class=\"fontawesome-icon fa-lightbulb far\" aria-hidden=\"true\" style=\"font-size:13px;\"><\/i>Understanding the output<\/h4><\/a><\/li><li role=\"presentation\"><a class=\"tab-link\" data-toggle=\"tab\" role=\"tab\" aria-controls=\"tab-69ff3988bccc2c93296\" aria-selected=\"false\" tabindex=\"-1\" id=\"fusion-tab-dockingmeta-data\" href=\"#tab-69ff3988bccc2c93296\"><h4 class=\"fusion-tab-heading\"><i class=\"fontawesome-icon fa-search fas\" aria-hidden=\"true\" style=\"font-size:13px;\"><\/i>Docking meta-data<\/h4><\/a><\/li><\/ul><\/div><div class=\"tab-content\"><div class=\"nav fusion-mobile-tab-nav\"><ul class=\"nav-tabs\" role=\"tablist\"><li class=\"active\" role=\"presentation\"><a class=\"tab-link\" data-toggle=\"tab\" role=\"tab\" aria-controls=\"tab-74538996c2b940c7d33\" aria-selected=\"true\" id=\"mobile-fusion-tab-generatetargetfile\" href=\"#tab-74538996c2b940c7d33\"><h4 class=\"fusion-tab-heading\"><i class=\"fontawesome-icon fa-cogs fas\" aria-hidden=\"true\" style=\"font-size:13px;\"><\/i>Generate target file<\/h4><\/a><\/li><\/ul><\/div><div class=\"tab-pane fade fusion-clearfix in active\" role=\"tabpanel\" tabindex=\"0\" aria-labelledby=\"fusion-tab-generatetargetfile\" id=\"tab-74538996c2b940c7d33\">\n<h4>Generate the target file containing the affinity maps.<\/h4>\n<style type=\"text\/css\" scopped=\"scopped\">.fusion-syntax-highlighter-1 > .CodeMirror, .fusion-syntax-highlighter-1 > .CodeMirror .CodeMirror-gutters {background-color:var(--awb-color1);}<\/style><div class=\"fusion-syntax-highlighter-container fusion-syntax-highlighter-1 fusion-syntax-highlighter-theme-light\" style=\"opacity:0;margin-top:0px;margin-right:0px;margin-bottom:0px;margin-left:0px;font-size:14px;border-width:1px;border-style:solid;border-color:var(--awb-color3);\"><div class=\"syntax-highlighter-copy-code\"><span class=\"syntax-highlighter-copy-code-title\" data-id=\"fusion_syntax_highlighter_1\" style=\"font-size:14px;\">Copy to Clipboard<\/span><\/div><label for=\"fusion_syntax_highlighter_1\" class=\"screen-reader-text\">Syntax Highlighter<\/label><textarea class=\"fusion-syntax-highlighter-textarea\" id=\"fusion_syntax_highlighter_1\" data-readOnly=\"nocursor\" data-lineNumbers=\"\" data-lineWrapping=\"1\" data-theme=\"default\" data-mode=\"text\/x-sh\">agfr -r data\/3c9w.pdbqt -b user 28.565 6.329 6.985 22.5 22.5 22.5 -c 1593 1596 -t 1591 -x A:CYS164 -o 3c9w_cov_cmdline<\/textarea><\/div>\n<div class=\"fusion-sep-clear\"><\/div><div class=\"fusion-separator fusion-full-width-sep\" style=\"margin-left: auto;margin-right: auto;margin-top:15px;margin-bottom:15px;width:100%;\"><div class=\"fusion-separator-border sep-double sep-solid\" style=\"--awb-height:20px;--awb-amount:20px;border-color:#e0dede;border-top-width:1px;border-bottom-width:1px;\"><\/div><\/div><div class=\"fusion-sep-clear\"><\/div>\n<p><strong>Details<\/strong>:<\/p>\n<ul>\n<li>In this example we specify the box placement and size manually to illustrate this capability. Note that this is the only <strong>-b\/-boxMode<\/strong> option for which padding is ignored.<\/li>\n<li>The PDB serial numbers of the 2 receptor atoms forming the covalent bond are specified using the <strong>-c\/&#8211;covalentBond<\/strong> option. The numbers are the serial codes appearing in the pdbqt file.<\/li>\n<li>The third atom defining the covalent attachment is used to compute the torsion angle of the covalent bond. It is specified using the <strong>-t\/&#8211;covalentBondTorsionAtom<\/strong> option.<\/li>\n<li>The <strong>-x\/&#8211;covalentResidues<\/strong> option allows to limit the traversal of the receptor to a list or residues. This is needed sometimes as covalent residues can create bonds with the receptor other than the covalent attachment and this is the case with the native ligand of 3cw9. When <em>agfr<\/em> identifies the sub-tree beyond the covalent bond to find which atoms to cut out of the receptor for calculating affinity maps, it would include a large part of the receptor because of the spurious bond the ligand makes with the receptor and cut out in excess of 1600 atoms. The -x option prevents this from happening.<\/li>\n<li>the target file is called\u00a03c9w_cov_cmdline.trg (<strong>-o<\/strong> option)<\/li>\n<\/ul>\n<p>The output of this commands is shown below:<\/p>\n<style type=\"text\/css\" scopped=\"scopped\">.fusion-syntax-highlighter-2 > .CodeMirror, .fusion-syntax-highlighter-2 > .CodeMirror .CodeMirror-gutters {background-color:var(--awb-color1);}.fusion-syntax-highlighter-2 > .CodeMirror .CodeMirror-gutters { background-color: var(--awb-color2); }.fusion-syntax-highlighter-2 > .CodeMirror .CodeMirror-linenumber { color: var(--awb-color8); }<\/style><div class=\"fusion-syntax-highlighter-container fusion-syntax-highlighter-2 fusion-syntax-highlighter-theme-light\" style=\"opacity:0;margin-top:0px;margin-right:0px;margin-bottom:0px;margin-left:0px;font-size:14px;border-width:1px;border-style:solid;border-color:var(--awb-color3);\"><label for=\"fusion_syntax_highlighter_2\" class=\"screen-reader-text\">Syntax Highlighter<\/label><textarea class=\"fusion-syntax-highlighter-textarea\" id=\"fusion_syntax_highlighter_2\" data-readOnly=\"nocursor\" data-lineNumbers=\"1\" data-lineWrapping=\"\" data-theme=\"default\" data-mode=\"text\/x-sh\">#################################################################\n# If you used AGFR in your work, please cite:                   #\n#                                                               #\n# P.A. Ravindranath S. Forli, D.S. Goodsell, A.J. Olson and     #\n# M.F. Sanner                                                   #\n# AutoDockFR: Advances in Protein-Ligand Docking with           #\n# Explicitly Specified Binding Site Flexibility                 #\n# PLoS Comput Biol 11(12): e1004586                             #\n# DOI:10.1371\/journal.pcbi.1004586                              #\n#                                                               #\n# P. Ananad Ravindranath and M.F. Sanner                        #\n# AutoSite: an automated approach for pseudoligands prediction  #\n# - From ligand binding sites identification to predicting key  #\n# ligand atoms                                                  #\n# Bioinformatics (2016)                                         #\n# DOI:10.1093\/bioinformatics\/btw367                             #\n#                                                               #\n# Please see http:\/\/adfr.scripps.edu for more information.      #\n#################################################################\n\nComputing grids on fiji a Darwin-17.7.0-x86_64-i386-64bit computer\nDate Thu Apr 18 12:31:04 2019\n\nloading receptor: data\/3c9w.pdbqt\nset box using user\n    Box center:    28.565     6.329     6.985\n    Box length:    22.500    22.500    22.500\n    Box size  :        60        60        60\n    padding   :     0.000\n    spacing   :     0.375\n\nsetting map types using: all to ['HS', 'Mg', 'HD', 'NA', 'Fe', 'Br', 'NS', 'A', 'C', 'Mn', 'G', 'F', 'I', 'H', 'J', 'N', 'Q', 'P', 'S', 'GA', 'Z', 'Zn', 'Cl', 'Ca', 'OA', 'SA', 'OS']\n\ncomputing maps for center=(28.565 6.329 6.985) size=(22.500 22.500 22.500) dims=(60 60 60) ...\n    maps computed in 9.86 (sec)\nthe following 33 covalent ligand atoms did not contribute to the grid calculation:\n  A:CYS164:CA,CB,SG,O23,C7,O8,C9,C25,C10,C11,C12,C2,C1,O20,H20,C6,C5,O21,C22,C4,C3,C18,O31,C17,C16,C15,O29,H29,C14,O28,H28,C13,O24,\nAdding gradient to maps ...\nprocessing maps ... done 26.2598040104\nwriting maps ... done 4.58263802528\ndone adding gradient to maps 33.02 (sec)\nmaking target file 3c9w_cov_cmdline.trg ...done.\n    done. 45.70 (sec)\n<\/textarea><\/div>\n<\/div><div class=\"nav fusion-mobile-tab-nav\"><ul class=\"nav-tabs\" role=\"tablist\"><li role=\"presentation\"><a class=\"tab-link\" data-toggle=\"tab\" role=\"tab\" aria-controls=\"tab-00c3697ba7ca4dd7edf\" aria-selected=\"false\" tabindex=\"-1\" id=\"mobile-fusion-tab-inspecttargetfile\" href=\"#tab-00c3697ba7ca4dd7edf\"><h4 class=\"fusion-tab-heading\"><i class=\"fontawesome-icon fa-search fas\" aria-hidden=\"true\" style=\"font-size:13px;\"><\/i>Inspect target file<\/h4><\/a><\/li><\/ul><\/div><div class=\"tab-pane fade fusion-clearfix\" role=\"tabpanel\" tabindex=\"0\" aria-labelledby=\"fusion-tab-inspecttargetfile\" id=\"tab-00c3697ba7ca4dd7edf\">\n<h4>Display information about a target file<\/h4>\n<style type=\"text\/css\" scopped=\"scopped\">.fusion-syntax-highlighter-3 > .CodeMirror, .fusion-syntax-highlighter-3 > .CodeMirror .CodeMirror-gutters {background-color:var(--awb-color1);}<\/style><div class=\"fusion-syntax-highlighter-container fusion-syntax-highlighter-3 fusion-syntax-highlighter-theme-light\" style=\"opacity:0;margin-top:0px;margin-right:0px;margin-bottom:0px;margin-left:0px;font-size:14px;border-width:1px;border-style:solid;border-color:var(--awb-color3);\"><div class=\"syntax-highlighter-copy-code\"><span class=\"syntax-highlighter-copy-code-title\" data-id=\"fusion_syntax_highlighter_3\" style=\"font-size:14px;\">Copy to Clipboard<\/span><\/div><label for=\"fusion_syntax_highlighter_3\" class=\"screen-reader-text\">Syntax Highlighter<\/label><textarea class=\"fusion-syntax-highlighter-textarea\" id=\"fusion_syntax_highlighter_3\" data-readOnly=\"nocursor\" data-lineNumbers=\"\" data-lineWrapping=\"\" data-theme=\"default\" data-mode=\"text\/x-sh\">about 3c9w_cov_cmdline.trg<\/textarea><\/div>\n<div class=\"fusion-sep-clear\"><\/div><div class=\"fusion-separator fusion-full-width-sep\" style=\"margin-left: auto;margin-right: auto;margin-top:15px;margin-bottom:15px;width:100%;\"><div class=\"fusion-separator-border sep-double sep-solid\" style=\"--awb-height:20px;--awb-amount:20px;border-color:#e0dede;border-top-width:1px;border-bottom-width:1px;\"><\/div><\/div><div class=\"fusion-sep-clear\"><\/div>\n<p>Details: the target file meta data is read and displayed.<\/p>\n<p>The command produces the following output:<\/p>\n<style type=\"text\/css\" scopped=\"scopped\">.fusion-syntax-highlighter-4 > .CodeMirror, .fusion-syntax-highlighter-4 > .CodeMirror .CodeMirror-gutters {background-color:var(--awb-color1);}.fusion-syntax-highlighter-4 > .CodeMirror .CodeMirror-gutters { background-color: var(--awb-color2); }.fusion-syntax-highlighter-4 > .CodeMirror .CodeMirror-linenumber { color: var(--awb-color8); }<\/style><div class=\"fusion-syntax-highlighter-container fusion-syntax-highlighter-4 fusion-syntax-highlighter-theme-light\" style=\"opacity:0;margin-top:0px;margin-right:0px;margin-bottom:0px;margin-left:0px;font-size:14px;border-width:1px;border-style:solid;border-color:var(--awb-color3);\"><label for=\"fusion_syntax_highlighter_4\" class=\"screen-reader-text\">Syntax Highlighter<\/label><textarea class=\"fusion-syntax-highlighter-textarea\" id=\"fusion_syntax_highlighter_4\" data-readOnly=\"nocursor\" data-lineNumbers=\"1\" data-lineWrapping=\"1\" data-theme=\"default\" data-mode=\"text\/x-sh\">docking target file\n  date       : Thu Apr 18 13:41:22 2019\n  node       : fiji\n  AGFR       : v1.2\n\n  receptor   : 3c9w.pdbqt\n    FlexRec  : None\n    covBond  : [3c9w:_A:CYS164:N (1591)] 3c9w:_A:CYS164:CA (1593) -- 3c9w:_A:CYS164:CB (1596)\n    coords   : [23.521, 5.081, 12.263] [24.700, 4.867, 11.312] [23.696, 4.237, 13.432]\n    covRes   : A:CYS164\n    covFile  : covalenLig.pdbqt\n    ignAtms  : [1592, 1595, 1596, 3349, 3350, 3351, 3352, 3353, 3354, 3355, 3356, 3357, 3358, 3359, 3360, 3361, 3362, 3363, 3364, 3365, 3366, 3367, 3368, 3369, 3370, 3371, 3372, 3373, 3374, 3375, 3376, 3377, 3378]\n\n  box        :  mode ['user', '28.565', '6.329', '6.985', '22.5', '22.5', '22.5'], padding 4.00\n    center   : 28.565 6.329 6.985\n    length   : 22.500 22.500 22.500\n    size     : 0060 0060 0060\n    spacing  : 0.375\n\n  maps       : \n    types    :  C  GA  Fe   A   e   d   S   Q   F  SA   G   P  Mg  Cl  Zn  OA  NS  Ca  Br  Mn  OS  \n               NA   J   I  HS  HD   H   Z   N   W \n    W map    : weight 0.60 entropy -0.20\n    gradients: Yes,  kept largest negative cluster\n\n  pocketMode : ['user', '28.565', '6.329', '6.985', '22.5', '22.5', '22.5']\n    #fillpts : 0 points\n<\/textarea><\/div>\n<\/div><div class=\"nav fusion-mobile-tab-nav\"><ul class=\"nav-tabs\" role=\"tablist\"><li role=\"presentation\"><a class=\"tab-link\" data-toggle=\"tab\" role=\"tab\" aria-controls=\"tab-4d64c9c233eda02d1a1\" aria-selected=\"false\" tabindex=\"-1\" id=\"mobile-fusion-tab-dockingligand\" href=\"#tab-4d64c9c233eda02d1a1\"><h4 class=\"fusion-tab-heading\"><i class=\"fontawesome-icon fa-cogs fas\" aria-hidden=\"true\" style=\"font-size:13px;\"><\/i>Docking ligand<\/h4><\/a><\/li><\/ul><\/div><div class=\"tab-pane fade fusion-clearfix\" role=\"tabpanel\" tabindex=\"0\" aria-labelledby=\"fusion-tab-dockingligand\" id=\"tab-4d64c9c233eda02d1a1\">\n<h3>Dock the randomized ligand using the generated target file<\/h3>\n<style type=\"text\/css\" scopped=\"scopped\">.fusion-syntax-highlighter-5 > .CodeMirror, .fusion-syntax-highlighter-5 > .CodeMirror .CodeMirror-gutters {background-color:var(--awb-color1);}<\/style><div class=\"fusion-syntax-highlighter-container fusion-syntax-highlighter-5 fusion-syntax-highlighter-theme-light\" style=\"opacity:0;margin-top:0px;margin-right:0px;margin-bottom:0px;margin-left:0px;font-size:14px;border-width:1px;border-style:solid;border-color:var(--awb-color3);\"><div class=\"syntax-highlighter-copy-code\"><span class=\"syntax-highlighter-copy-code-title\" data-id=\"fusion_syntax_highlighter_5\" style=\"font-size:14px;\">Copy to Clipboard<\/span><\/div><label for=\"fusion_syntax_highlighter_5\" class=\"screen-reader-text\">Syntax Highlighter<\/label><textarea class=\"fusion-syntax-highlighter-textarea\" id=\"fusion_syntax_highlighter_5\" data-readOnly=\"nocursor\" data-lineNumbers=\"\" data-lineWrapping=\"1\" data-theme=\"default\" data-mode=\"text\/x-sh\">adfr -l data\/3c9w_ligandWithSideChain_random.pdbqt -t 3c9w_cov_cmdline.trg --jobName covalent -C 1 2 3 --nbRuns 8 --maxEvals 100000 -O --seed 1<\/textarea><\/div>\n<div class=\"fusion-sep-clear\"><\/div><div class=\"fusion-separator fusion-full-width-sep\" style=\"margin-left: auto;margin-right: auto;margin-top:15px;margin-bottom:15px;width:100%;\"><div class=\"fusion-separator-border sep-double sep-solid\" style=\"--awb-height:20px;--awb-amount:20px;border-color:#e0dede;border-top-width:1px;border-bottom-width:1px;\"><\/div><\/div><div class=\"fusion-sep-clear\"><\/div>\n<p>Details: Here we\u00a0re-dock the native ligand, that has been randomized (i.e. its conformation as well as it positions and orientation in the crystal structure have been randomly modified).\u00a0<em>adfr<\/em> detects the number of cores available and by default will use them all to perform 8 independent searches (&#8211;<span class=\"SpellE\">nbRuns<\/span> 8) each using up to 100&#8217;000 evaluations of the scoring function (&#8211;<span class=\"SpellE\">maxEvals<\/span> 100&#8217;00). By default <em>adfr<\/em> performs 50 searches, each allotted 2.5 million evaluations. Typically, more complex docking problems require more searches to be performed to increase the chances to find the best possible docked pose (i.e. global minimum of the scoring function). Here we set these parameters to lower values to perform a quick run that is sufficient to illustrate the docking principles.<\/p>\n<p>This calculation generates the following files:<\/p>\n<ul>\n<li>\n<pre>3c9w_ligandWithSideChain_random_covalent_summary.dlg\u00a0 \u00a0# the docking log file that captures most of what is displayed on the terminal<\/pre>\n<\/li>\n<li>\n<pre>3c9w_ligandWithSideChain_random_covalent_out.pdbqt\u00a0 \u00a0 \u00a0# the docking pose file containing the docking solutions<\/pre>\n<\/li>\n<li>\n<pre>3c9w_ligandWithSideChain_random_covalent.dro\u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0# the docking object file that contains input, output, and meta-data for this docking run<\/pre>\n<\/li>\n<\/ul>\n<p>NOTES:<\/p>\n<ul>\n<li>The output files are named using\u00a0the ligand name followed by the job name (&#8211;<span class=\"SpellE\">jobName<\/span>\u00a0if specified)<\/li>\n<li>ADFR&#8217;s search procedure is stochastic, meaning that docking the same ligand into the same target twice can produce different results if\u00a0 different random number generator seeds are used. However, the energy landscape for this receptor and ligand is the same in both runs. If both docking runs find the global minimum of this energy landscape, the solutions produced by both runs will be the same, independently of the paths taken by the search to get there. On the other hand, searches that get trapped in a local minima, yield docking poses that differ from each other. Specifying the seeds used by the random number generator (&#8211;seed) allows reproducing a docking calculation, for a given version of the code.<\/li>\n<\/ul>\n<p>The output of the command is listed below:<\/p>\n<style type=\"text\/css\" scopped=\"scopped\">.fusion-syntax-highlighter-6 > .CodeMirror, .fusion-syntax-highlighter-6 > .CodeMirror .CodeMirror-gutters {background-color:var(--awb-color1);}.fusion-syntax-highlighter-6 > .CodeMirror .CodeMirror-gutters { background-color: var(--awb-color2); }.fusion-syntax-highlighter-6 > .CodeMirror .CodeMirror-linenumber { color: var(--awb-color8); }<\/style><div class=\"fusion-syntax-highlighter-container fusion-syntax-highlighter-6 fusion-syntax-highlighter-theme-light\" style=\"opacity:0;margin-top:0px;margin-right:0px;margin-bottom:0px;margin-left:0px;font-size:14px;border-width:1px;border-style:solid;border-color:var(--awb-color3);\"><label for=\"fusion_syntax_highlighter_6\" class=\"screen-reader-text\">Syntax Highlighter<\/label><textarea class=\"fusion-syntax-highlighter-textarea\" id=\"fusion_syntax_highlighter_6\" data-readOnly=\"nocursor\" data-lineNumbers=\"1\" data-lineWrapping=\"1\" data-theme=\"default\" data-mode=\"text\/x-sh\">#################################################################\n# If you used ADFR in your work, please cite:                   #\n#                                                               #\n# P.A. Ravindranath S. Forli, D.S. Goodsell, A.J. Olson and     #\n# M.F. Sanner                                                   #\n#                                                               #\n# AutoDockFR: Advances in Protein-Ligand Docking with           #\n# Explicitly Specified Binding Site Flexibility                 #\n# PLoS Comput Biol 11(12): e1004586                             #\n#                                                               #\n# DOI:10.1371\/journal.pcbi.1004586                              #\n#                                                               #\n# Please see http:\/\/adfr.scripps.edu for more information.      #\n#################################################################\n\nDocking on fiji a Darwin-17.7.0-x86_64-i386-64bit computer\nDate Thu Apr 18 13:46:27 2019\nreading ligand \/Users\/sanner\/test1.1\/data\/3c9w_ligandWithSideChain_random.pdbqt\nDetected 4 cores, using 4 cores\nUnpacking maps \/Users\/sanner\/test1.1\/3c9w_cov_cmdline.trg\nPerforming search (8 GA evolutions with 100000 maxEvals each) ...\n0%   10   20   30   40   50   60   70   80   90   100%\n|----|----|----|----|----|----|----|----|----|----|\n***************************************************\nTermination status\n    0\/   8  0.0% runs failed\n    8\/   8 100.0% runs exhausted their evaluations\n    0\/   8  0.0% runs stopped converged 1 or 2 clusters\n    0\/   8  0.0% runs stopped after no improvement in clusters\n    0\/   8  0.0% runs stopped because GA ran out of choices\n    0\/   8  0.0% runs stopped because GA population converged\n\nRefining results ...\ndone.\n\nmode |  affinity  | clust. | ref. | clust. | rmsd | energy | best |\n     | (kcal\/mol) | rmsd   | rmsd |  size  | stdv |  stdv  | run  |\n-----+------------+--------+------+--------+------+--------+------+\n   1         -6.3     0.0    -1.0       8     0.2     0.1    006\nWriting poses to 3c9w_ligandWithSideChain_random_covalent_out.pdbqt\n\nDocking performed in 13.17 seconds, i.e. 0 hours 00 minutes 13.165729 seconds \n\n*************** first GA command ***************************\n\"\/Users\/sanner\/test1.1\/mgltools2_x86_64Darwin_1.2\/bin\/pythonsh\" \"\/Users\/sanner\/test1.1\/ADFR\/bin\/runOneGA.py\" -F \"\/var\/folders\/gm\/frz8gxj57hzbyyzwtt5bq8m40000gn\/T\/tmpaAcjqR\/3c9w_cov_cmdline\" -M rigidReceptor -R \"\/var\/folders\/gm\/frz8gxj57hzbyyzwtt5bq8m40000gn\/T\/tmpaAcjqR\/3c9w_cov_cmdline\/3c9w.pdbqt\" \"-l\" \"data\/3c9w_ligandWithSideChain_random.pdbqt\" \"--jobName\" \"covalent\" \"-C\" \"1\" \"2\" \"3\" \"--maxEvals\" \"100000\" \"-O\" -V \"1591\" \"1593\" \"1596\" -S 1 -j 1 -o \"3c9w_ligandWithSideChain_random_covalent\/covalent0001.dlg\"\n<\/textarea><\/div>\n<\/div><div class=\"nav fusion-mobile-tab-nav\"><ul class=\"nav-tabs\" role=\"tablist\"><li role=\"presentation\"><a class=\"tab-link\" data-toggle=\"tab\" role=\"tab\" aria-controls=\"tab-630b551d9b728b54dea\" aria-selected=\"false\" tabindex=\"-1\" id=\"mobile-fusion-tab-understandingtheoutput\" href=\"#tab-630b551d9b728b54dea\"><h4 class=\"fusion-tab-heading\"><i class=\"fontawesome-icon fa-lightbulb far\" aria-hidden=\"true\" style=\"font-size:13px;\"><\/i>Understanding the output<\/h4><\/a><\/li><\/ul><\/div><div class=\"tab-pane fade fusion-clearfix\" role=\"tabpanel\" tabindex=\"0\" aria-labelledby=\"fusion-tab-understandingtheoutput\" id=\"tab-630b551d9b728b54dea\">\n<p>Here we describe line by line the messages output during the docking procedure.<\/p>\n<style type=\"text\/css\" scopped=\"scopped\">.fusion-syntax-highlighter-7 > .CodeMirror, .fusion-syntax-highlighter-7 > .CodeMirror .CodeMirror-gutters {background-color:var(--awb-color1);}.fusion-syntax-highlighter-7 > .CodeMirror .CodeMirror-gutters { background-color: var(--awb-color2); }.fusion-syntax-highlighter-7 > .CodeMirror .CodeMirror-linenumber { color: var(--awb-color8); }<\/style><div class=\"fusion-syntax-highlighter-container fusion-syntax-highlighter-7 fusion-syntax-highlighter-theme-light\" style=\"opacity:0;margin-top:0px;margin-right:0px;margin-bottom:15px;margin-left:0px;font-size:14px;border-width:1px;border-style:solid;border-color:var(--awb-color3);\"><label for=\"fusion_syntax_highlighter_7\" class=\"screen-reader-text\">Syntax Highlighter<\/label><textarea class=\"fusion-syntax-highlighter-textarea\" id=\"fusion_syntax_highlighter_7\" data-readOnly=\"nocursor\" data-lineNumbers=\"1\" data-lineWrapping=\"\" data-theme=\"default\" data-mode=\"text\/x-sh\">Docking on fiji a Darwin-17.7.0-x86_64-i386-64bit computer\nDate Thu Apr 18 13:46:27 2019\nreading ligand \/Users\/sanner\/test1.1\/data\/3c9w_ligandWithSideChain_random.pdbqt\nDetected 4 cores, using 4 cores\nUnpacking maps \/Users\/sanner\/test1.1\/3c9w_cov_cmdline.trg<\/textarea><\/div>\n<ol>\n<li>Hostname and platform architecture on which the program is running<\/li>\n<li>Date and time of execution<\/li>\n<li>ligand docked<\/li>\n<li>number of detected and used cores.<\/li>\n<li>target files used<\/li>\n<\/ol>\n<p>NOTES:<\/p>\n<ul>\n<li>Number of cores. By default ADFR will use all cores available to parallelize the search threads comprised in a run. Use the &#8220;-c&#8221; command line option to limit the number of used cores.<\/li>\n<\/ul>\n<style type=\"text\/css\" scopped=\"scopped\">.fusion-syntax-highlighter-8 > .CodeMirror, .fusion-syntax-highlighter-8 > .CodeMirror .CodeMirror-gutters {background-color:var(--awb-color1);}.fusion-syntax-highlighter-8 > .CodeMirror .CodeMirror-gutters { background-color: var(--awb-color2); }.fusion-syntax-highlighter-8 > .CodeMirror .CodeMirror-linenumber { color: var(--awb-color8); }<\/style><div class=\"fusion-syntax-highlighter-container fusion-syntax-highlighter-8 fusion-syntax-highlighter-theme-light\" style=\"opacity:0;margin-top:0px;margin-right:0px;margin-bottom:15px;margin-left:0px;font-size:14px;border-width:1px;border-style:solid;border-color:var(--awb-color3);\"><label for=\"fusion_syntax_highlighter_8\" class=\"screen-reader-text\">Syntax Highlighter<\/label><textarea class=\"fusion-syntax-highlighter-textarea\" id=\"fusion_syntax_highlighter_8\" data-readOnly=\"nocursor\" data-lineNumbers=\"1\" data-lineWrapping=\"\" data-theme=\"default\" data-mode=\"text\/x-sh\">Performing search (8 GA evolutions with 100000 maxEvals each) ...\n0%   10   20   30   40   50   60   70   80   90   100%\n|----|----|----|----|----|----|----|----|----|----|\n***************************************************\nTermination status\n    0\/   8  0.0% runs failed\n    8\/   8 100.0% runs exhausted their evaluations\n    0\/   8  0.0% runs stopped converged 1 or 2 clusters\n    0\/   8  0.0% runs stopped after no improvement in clusters\n    0\/   8  0.0% runs stopped because GA ran out of choices\n    0\/   8  0.0% runs stopped because GA population converged\n<\/textarea><\/div>\n<p>By default, ADFR performs 50 independent searches, i.e. 50 evolutions of a population of 100 individuals using a Genetic Algorithm (GA). In this example we intentionally reduced this number to 8 very short runs.<\/p>\n<p>Lines 1-4 display a progress bar indicating the percentage of these runs that completed.<\/p>\n<p>The lines below provide statistics over the termination status of these searches. <em>ADFR<\/em>\u00a0implements several termination criteria in its search method. In this example all search terminated because they reached their maximum number of evaluations. The default number of evaluations is 2.5 millions and is usually never reached because of other termination criteria such as convergence of the population, meaning that there is no more diversity in the population and the chances to discover new solutions has become small, or the population still has diversity (i.e. it contains multiples competitive solutions) but none of these solution has improves over a user-defined number of generations (default 5).<\/p>\n<p>Typically, you want searches to end because the population converged or there was no improvement. A result like the one shown here is a clear indication that this docking problem needs more evaluations per search (i.e. increased &#8211;maxEvals).<\/p>\n<p>The next section lists the results:<\/p>\n<style type=\"text\/css\" scopped=\"scopped\">.fusion-syntax-highlighter-9 > .CodeMirror, .fusion-syntax-highlighter-9 > .CodeMirror .CodeMirror-gutters {background-color:var(--awb-color1);}.fusion-syntax-highlighter-9 > .CodeMirror .CodeMirror-gutters { background-color: var(--awb-color2); }.fusion-syntax-highlighter-9 > .CodeMirror .CodeMirror-linenumber { color: var(--awb-color8); }<\/style><div class=\"fusion-syntax-highlighter-container fusion-syntax-highlighter-9 fusion-syntax-highlighter-theme-light\" style=\"opacity:0;margin-top:0px;margin-right:0px;margin-bottom:15px;margin-left:0px;font-size:14px;border-width:1px;border-style:solid;border-color:var(--awb-color3);\"><label for=\"fusion_syntax_highlighter_9\" class=\"screen-reader-text\">Syntax Highlighter<\/label><textarea class=\"fusion-syntax-highlighter-textarea\" id=\"fusion_syntax_highlighter_9\" data-readOnly=\"nocursor\" data-lineNumbers=\"1\" data-lineWrapping=\"\" data-theme=\"default\" data-mode=\"text\/x-sh\">mode |  affinity  | clust. | ref. | clust. | rmsd | energy | best |\n     | (kcal\/mol) | rmsd   | rmsd |  size  | stdv |  stdv  | run  |\n-----+------------+--------+------+--------+------+--------+------+\n   1         -6.3     0.0    -1.0       8     0.2     0.1    006<\/textarea><\/div>\n<p>In this docking run, the 8 searches lead to the same solution as indicated by the cluster size 8 on the first and only solution (<em>clust. size<\/em> column).<\/p>\n<\/div><div class=\"nav fusion-mobile-tab-nav\"><ul class=\"nav-tabs\" role=\"tablist\"><li role=\"presentation\"><a class=\"tab-link\" data-toggle=\"tab\" role=\"tab\" aria-controls=\"tab-69ff3988bccc2c93296\" aria-selected=\"false\" tabindex=\"-1\" id=\"mobile-fusion-tab-dockingmeta-data\" href=\"#tab-69ff3988bccc2c93296\"><h4 class=\"fusion-tab-heading\"><i class=\"fontawesome-icon fa-search fas\" aria-hidden=\"true\" style=\"font-size:13px;\"><\/i>Docking meta-data<\/h4><\/a><\/li><\/ul><\/div><div class=\"tab-pane fade fusion-clearfix\" role=\"tabpanel\" tabindex=\"0\" aria-labelledby=\"fusion-tab-dockingmeta-data\" id=\"tab-69ff3988bccc2c93296\">\n<h4>Display information about a docking result<\/h4>\n<style type=\"text\/css\" scopped=\"scopped\">.fusion-syntax-highlighter-10 > .CodeMirror, .fusion-syntax-highlighter-10 > .CodeMirror .CodeMirror-gutters {background-color:var(--awb-color1);}<\/style><div class=\"fusion-syntax-highlighter-container fusion-syntax-highlighter-10 fusion-syntax-highlighter-theme-light\" style=\"opacity:0;margin-top:0px;margin-right:0px;margin-bottom:0px;margin-left:0px;font-size:14px;border-width:1px;border-style:solid;border-color:var(--awb-color3);\"><div class=\"syntax-highlighter-copy-code\"><span class=\"syntax-highlighter-copy-code-title\" data-id=\"fusion_syntax_highlighter_10\" style=\"font-size:14px;\">Copy to Clipboard<\/span><\/div><label for=\"fusion_syntax_highlighter_10\" class=\"screen-reader-text\">Syntax Highlighter<\/label><textarea class=\"fusion-syntax-highlighter-textarea\" id=\"fusion_syntax_highlighter_10\" data-readOnly=\"nocursor\" data-lineNumbers=\"\" data-lineWrapping=\"\" data-theme=\"default\" data-mode=\"text\/x-sh\">about 3c9w_ligandWithSideChain_random_covalent.dro<\/textarea><\/div>\n<div class=\"fusion-sep-clear\"><\/div><div class=\"fusion-separator fusion-full-width-sep\" style=\"margin-left: auto;margin-right: auto;margin-top:15px;margin-bottom:15px;width:100%;\"><div class=\"fusion-separator-border sep-double sep-solid\" style=\"--awb-height:20px;--awb-amount:20px;border-color:#e0dede;border-top-width:1px;border-bottom-width:1px;\"><\/div><\/div><div class=\"fusion-sep-clear\"><\/div>\n<p>Details: the meta data about this docking run is displayed<\/p>\n<style type=\"text\/css\" scopped=\"scopped\">.fusion-syntax-highlighter-11 > .CodeMirror, .fusion-syntax-highlighter-11 > .CodeMirror .CodeMirror-gutters {background-color:var(--awb-color1);}.fusion-syntax-highlighter-11 > .CodeMirror .CodeMirror-gutters { background-color: var(--awb-color2); }.fusion-syntax-highlighter-11 > .CodeMirror .CodeMirror-linenumber { color: var(--awb-color8); }<\/style><div class=\"fusion-syntax-highlighter-container fusion-syntax-highlighter-11 fusion-syntax-highlighter-theme-light\" style=\"opacity:0;margin-top:0px;margin-right:0px;margin-bottom:0px;margin-left:0px;font-size:14px;border-width:1px;border-style:solid;border-color:var(--awb-color3);\"><label for=\"fusion_syntax_highlighter_11\" class=\"screen-reader-text\">Syntax Highlighter<\/label><textarea class=\"fusion-syntax-highlighter-textarea\" id=\"fusion_syntax_highlighter_11\" data-readOnly=\"nocursor\" data-lineNumbers=\"1\" data-lineWrapping=\"1\" data-theme=\"default\" data-mode=\"text\/x-sh\">docking result file\n  date       : Thu Apr 18 13:46:27 2019\n  node       : fiji\n  platform   : Darwin-17.7.0-x86_64-i386-64bit\n  ncores     : 4\n\n  receptor   : 3c9w_cov_cmdline.trg\n   docking target file\n     date       : Thu Apr 18 13:41:22 2019\n     node       : fiji\n     AGFR       : v1.2\n\n     receptor   : 3c9w.pdbqt\n       FlexRec  : None\n       covBond  : [3c9w:_A:CYS164:N (1591)] 3c9w:_A:CYS164:CA (1593) -- 3c9w:_A:CYS164:CB (1596)\n       coords   : [23.521, 5.081, 12.263] [24.700, 4.867, 11.312] [23.696, 4.237, 13.432]\n       covRes   : A:CYS164\n       covFile  : covalenLig.pdbqt\n       ignAtms  : [1592, 1595, 1596, 3349, 3350, 3351, 3352, 3353, 3354, 3355, 3356, 3357, 3358, 3359, 3360, 3361, 3362, 3363, 3364, 3365, 3366, 3367, 3368, 3369, 3370, 3371, 3372, 3373, 3374, 3375, 3376, 3377, 3378]\n\n     box        :  mode ['user', '28.565', '6.329', '6.985', '22.5', '22.5', '22.5'], padding 4.00\n       center   : 28.565 6.329 6.985\n       length   : 22.500 22.500 22.500\n       size     : 0060 0060 0060\n       spacing  : 0.375\n\n     maps       : \n       types    :  C  GA  Fe   A   e   d   S   Q   F  SA   G   P  Mg  Cl  Zn  OA  NS  Ca  Br  Mn  OS  \n                  NA   J   I  HS  HD   H   Z   N   W \n       W map    : weight 0.60 entropy -0.20\n       gradients: Yes,  kept largest negative cluster\n\n     pocketMode : ['user', '28.565', '6.329', '6.985', '22.5', '22.5', '22.5']\n       #fillpts : 0 points\n\n  ligand     : 3c9w_ligandWithSideChain_random.pdbqt\n     nbAtoms : 34\n     types   : A:6 C:16 OA:8 HD:3 SA:1 \n  lig ref    : None\n  summary    : 3c9w_ligandWithSideChain_random_covalent_summary.dlg\n  GA params  : \n    maxEvals : 100000\n    nbRuns   : 50\n    jobName  : covalent\n    seed     : 1\n    covLig   : [1, 2, 3]\n<\/textarea><\/div>\n<\/div><\/div><\/div><div class=\"fusion-clearfix\"><\/div><\/div><\/div><\/div><\/div>\n","protected":false},"excerpt":{"rendered":"","protected":false},"author":1,"featured_media":0,"parent":0,"menu_order":0,"comment_status":"closed","ping_status":"closed","template":"","meta":{"footnotes":""},"class_list":["post-936","page","type-page","status-publish","hentry"],"_links":{"self":[{"href":"https:\/\/ccsb.scripps.edu\/adfr\/wp-json\/wp\/v2\/pages\/936","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/ccsb.scripps.edu\/adfr\/wp-json\/wp\/v2\/pages"}],"about":[{"href":"https:\/\/ccsb.scripps.edu\/adfr\/wp-json\/wp\/v2\/types\/page"}],"author":[{"embeddable":true,"href":"https:\/\/ccsb.scripps.edu\/adfr\/wp-json\/wp\/v2\/users\/1"}],"replies":[{"embeddable":true,"href":"https:\/\/ccsb.scripps.edu\/adfr\/wp-json\/wp\/v2\/comments?post=936"}],"version-history":[{"count":9,"href":"https:\/\/ccsb.scripps.edu\/adfr\/wp-json\/wp\/v2\/pages\/936\/revisions"}],"predecessor-version":[{"id":940,"href":"https:\/\/ccsb.scripps.edu\/adfr\/wp-json\/wp\/v2\/pages\/936\/revisions\/940"}],"wp:attachment":[{"href":"https:\/\/ccsb.scripps.edu\/adfr\/wp-json\/wp\/v2\/media?parent=936"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}