Re: [AMBER] Generating a HETATM database to use with reduce

From: David A Case <>
Date: Thu, 14 Jul 2016 20:51:23 -0400

On Thu, Jul 14, 2016, Kucharski Jr., Amir N wrote:
> I am trying to do simulations on several protein-ligand complexes. I
> know from the SUS tutorial that reducing the ligand to add hydrogens
> is a key early step, but these ligands do not have entries in the
> database that reduce uses. To remedy this, we have tried making our
> own database from connectivity information in the mol2 files. We have
> followed the input format for the default library, but when we direct
> reduce to our database using the -DB tag, still no hydrogens are
> added. Is there another file we would need to make? Any help is much
> appreciated. Thanks!

What you describe should work, and you might have to put debugging information
into the reduce code to see why it is not doing what you expect.

However: if you already have a mol2 file with hydrogens, you can use tleap
to add hydrogens, bypassing the reduce step. Just load the mol2 file,
and a pdb file (missing hydrogens), and the output pdb file will have the
hydrogens added in their canonical positions.

...good luck...dac

AMBER mailing list
Received on Thu Jul 14 2016 - 18:00:03 PDT
Custom Search