[gmx-users] Use pdb file generated with Maestro, in Gromacs
bio hpc
biohpc2015 at gmail.com
Thu Mar 31 12:55:52 CEST 2016
Thanks.
Apart from this, is is there some script that keeps my original hydrogens from Maestro?
> El 30/3/2016, a las 19:43, Tsjerk Wassenaar <tsjerkw at gmail.com> escribió:
>
> Hey :)
>
> -ignh does ignore the hydrogens in the input file. It builds those
> hydrogens that are specified in the force field. For histidines, the
> protonation state is determined from the possible hydrogen-bonded network,
> but it is possible to assign specific states interactively, using the
> option -his, or by specifying the residue name in the input file.
>
> Hope it helps,
>
> Tsjerk
>
> On Wed, Mar 30, 2016 at 4:42 PM, sun <sun.iba2 at gmail.com> wrote:
>
>> I think using -ignh does not "remove" hydrogens. Hence, you can use it.
>>
>> Sent from my iPhone
>>
>>> On 30-Mar-2016, at 6:04 pm, rajendra kumar <rjdkmr at gmail.com> wrote:
>>>
>>> Hi,
>>>
>>> As suggested above, use -ignh to ignore hydrogen atoms. To use a specific
>>> Histidine, you may change residue name HIS to
>> HIP(Amber)/HSP(Charmm)/HISH(
>>> OPLS
>>> ) or HID(Amber)/HSD(Charmm)/HISD(
>>> OPLS
>>> ) or HIE(Amber)/HSE(Charmm)/HISE(
>>> OPLS
>>> ) either in PDB file or interactively through pdb2gmx usng -his command.
>>>
>>>
>>>> On Wed, Mar 30, 2016 at 12:27 PM, bio hpc <biohpc2015 at gmail.com> wrote:
>>>>
>>>> Hi,
>>>>
>>>> we have created some protein pdb files with Maestro. When we try to un
>> an
>>>> MD simulation with gromacs, we get errors like:
>>>>
>>>>>> Atom HD11 in residue ASN 5 was not found in rtp entry ASN with 14
>> atoms
>>>>
>>>> I tried to find a solution and it seems that internal gromacs dictionary
>>>> for hydrogens is different from the one from Maestro and other formats.
>> It
>>>> seems, from what I read, that one must convert all hydrogen types to
>>>> Gromacs format. Is is there any script that can do this automatically?
>>>>
>>>> Of course, there is the option “ignore hydrogens”, but this would be bad
>>>> for us, since we need some histidines in the binding site, and so on.
>>>>
>>>> Thanks,
>>>> BIO-HPC
>>>>
>>>>
>>>>
>>>> --
>>>> Gromacs Users mailing list
>>>>
>>>> * Please search the archive at
>>>> http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before
>>>> posting!
>>>>
>>>> * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
>>>>
>>>> * For (un)subscribe requests visit
>>>> https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or
>>>> send a mail to gmx-users-request at gromacs.org.
>>> --
>>> Gromacs Users mailing list
>>>
>>> * Please search the archive at
>> http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before
>> posting!
>>>
>>> * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
>>>
>>> * For (un)subscribe requests visit
>>> https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or
>> send a mail to gmx-users-request at gromacs.org.
>> --
>> Gromacs Users mailing list
>>
>> * Please search the archive at
>> http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before
>> posting!
>>
>> * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
>>
>> * For (un)subscribe requests visit
>> https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or
>> send a mail to gmx-users-request at gromacs.org.
>>
>
>
>
> --
> Tsjerk A. Wassenaar, Ph.D.
> --
> Gromacs Users mailing list
>
> * Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!
>
> * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
>
> * For (un)subscribe requests visit
> https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a mail to gmx-users-request at gromacs.org.
More information about the gromacs.org_gmx-users
mailing list