[Dock-fans] ligand preparation
pett at cgl.ucsf.edu
Wed Nov 30 10:52:07 PST 2011
> Is there a way to prepare the ligands (if more than 1000s) in a stand-alone
> program, basically to calculate the charges and write a mol2 file. I have a
> ligprep format ligands that i would like to convert to mol2 format with
> proper charges.
> I assume antechamber tools included in chimera can do this on a command
> prompt mode. Any ideas.
Can you convert your "ligprep format" files to some format that Chimera can read (PDB/Mol2/XYZ)?
If so, you can use the simple scripting approach described here:
to loop through your files, adding charges and writing them out as Mol2 files (with the "addcharge" and "write" commands respectively; check the Chimera documentation for the command arguments you would want to use).
Antechamber/addcharge will probably fail for a (hopefully small) fraction of your structures. Assuming you run in "command prompt" (nogui) mode the failures will be noted in the output that Chimera sends to the terminal. You could look through that to find them (though that output will be voluminous if you have thousands of structures) or you could write some kind of post-processing script that identifies Mol2 files where the charges are all zero (i.e. addcharge failed).
For the failures where Chimera estimated the wrong formal charge, you could rerun those structures using the graphical version of Chimera and specifying the correct formal charge when prompted by the AddCharge tool, then write the Mol2 files "by hand" with File->Save Mol2. For the failures where antechamber simply could not handle the ligand structure, you will have to use some other method for getting the charges for that structure.
If Chimera estimates the wrong formal charge for a structure, feel free to report it as a bug (Help->Report a Bug in Chimera), attaching the structure to the bug report.
UCSF Computer Graphics Lab
More information about the Dock-fans