Add hydrogen atoms on protein and water for Autodock.

I was asked to do docking with explicit water. I am confused whether I should add water manually or just use the crystal water. By searching for a while, it seemed the docking with water only use the water inside the crystal structure. So I just do that for the job now.

Autodock can not add hydrogen atoms on the lonely oxygen atoms in x-ray crystal structure. I tried some other programs, VMD, Schrodinger, HAAD server, WHAT IF. And finally, I found USCF chimera has a module to add hydrogen atoms very well. Although I am not sure whether the hydrogen atoms are optimized or not, at least they did not orient to one direction.

Autodock seems have trouble to deal with the structure got from chimera. I used VMD to extract water and add hydrogen atoms by chimera and save to PDB. Then I used VMD to extract protein and glycan of the protein crystal PDB file and save to another PDB file.

cat water.pdb protein.pdb > total.pdb

I use Autodock to open total.pdb, then add the polar hydrogen atoms. When I choose it for Grid process. The protein and water all be treated as the receptor for the following docking.

Maybe I can choose water as flexible residue for better docking result.


This afternoon, I got reply from Stephanie Nicole Lewis from Virginia Polytechnic Institute and State University. She is a great mentor and expert in Molecular Docking field. The docking with explicit water just use the crystal water. But in some condition, only water in specific region like docking pocket or caveat are meaningful. And she also use Chimera to add hydrogen atoms.

I read some papers and files from Google about FlexX. It’s a famous docking software, commercial, a part of LeadIT from BioSolveIt. This docking program have better solution for crystal water. It has three method to deal with water. One, add hydrogen atoms on the oxygen, and the orientation of hydrogen atoms should allow to form more hydrogen bonds like Chimera does. Two, do not add hydrogen atoms but treat it as a oxygen of water with uncertain orientations. This is better for ligands to form hydrogen bonds with the oxygen atoms of water. Three, the oxygen still does not have hydrogen atoms, and could bond with ligand or receptor in any direction. They can also be replaced by ligands if the water is an obstacle for the docking process. So the third way may be the best one. I may learn the software some day when I have time.


Leave a Reply

Fill in your details below or click an icon to log in: Logo

You are commenting using your account. Log Out / Change )

Twitter picture

You are commenting using your Twitter account. Log Out / Change )

Facebook photo

You are commenting using your Facebook account. Log Out / Change )

Google+ photo

You are commenting using your Google+ account. Log Out / Change )

Connecting to %s