Article ID: 266 - Last Modified: December 4, 2010
Can I generate Glide grids with hydrophobic constraints from the command line, or via a Python script?
Unfortunately, it isn't possible to script the creation of a hydrophobic constraint for Grid Generation, because the hydrophobic constraint requires an auxiliary .phob file containing details of the hydrophobic boxes generated in Maestro for a given receptor structure.
One alternative would be to use a positional constraint instead of a hydrophobic constraint. During the docking job setup, you can specify the ligand feature to be Hydrophobic for this positional constraint. The relevant Grid Generation input file keywords for a positional constraint are
where the *POS keywords indicate the coordinates of the center of the positional constraint, and the RPOS keyword indicates the radius of the positional constraint. You will need some mechanism for computing the appropriate coordinates for each receptor, but this can be done via a Python script (e.g., by calculating the centroid of certain residues common to all receptor structures). You also will need to increment GLIDE_NTOTALCONS and add a constraint label to GLIDE_CONSNAMES when you add the positional constraint.
Keywords: grid, generation, hydrophobic, constraint, script, command line
Type the words or phrases on which you would like to search, or click here to view a list of all
Knowledge Base articles