On Wed, Dec 02, 2020, Alexander Izvorski wrote:
>
>I'm trying to figure out how to use GAFF2 in antechamber/tleap.
>
>antechamber -i lig_in.mol2 -fi mol2 -o lig.mol2 -fo mol2 -c bcc -s 2 -rn LIG
>
>I noticed that "-at gaff2" is a possible option here - is it needed?
Yes.
>what happens is a bunch of errors when running tleap later,
>"Error! Could not find bond parameter for: ns - ca" and similar
Can you provide a small example that shows this error?
>2. Run parmchk: "-s gaff2" enables GAFF2 here
>3. Run tleap: source leaprc.gaff2 loads GAFF2
>
Above are correct.
>This works using sustiva.mol2
Since it works for sustiva, we would need to see the case that fails. I
use gaff2 a lot, and have not seen the problem you describe.
...thanks...dave
cc-ing to junmei: can't we finally make gaff2 the default? I suspect
tons of people are still using gaff1 because one needs to make three
different changes.
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Thu Dec 03 2020 - 06:00:02 PST