[gmx-users] Adding Hydrogens by pdb2gmx

[James Starlight]

Justin,

Could you tell me more about adding caps by amber tools?

1) I've found in that software ACE as well as other groups but they are not
standart

e.g ACE is below instead of simple CH3 group
             949 HH31 ACE    29       2.000   1.000  -0.000  1.00  0.00
ATOM    950  CH3 ACE    29       2.000   2.090   0.000  1.00  0.00
ATOM    951 HH32 ACE    29       1.486   2.454   0.890  1.00  0.00
ATOM    952 HH33 ACE    29       1.486   2.454  -0.890  1.00  0.00
ATOM    953  C   ACE    29       3.427   2.641  -0.000  1.00  0.00
ATOM    954  O   ACE    29       4.391   1.877  -0.000  1.00  0.00
ATOM    955  H7  ACE    29       3.582   3.629  -0.000  1.00  0.00

have you edited that group to standart form in amber tools or used
something else ?

2) How you place CAP groups to desired coordinates?
E.g if I use Edit- Unit Import Unit I can just import desired groups to my
pdb coordinates but how make linckage with N and C termi of my molecule ?
Or that linkage would be build by pdb2gmx -ter option ?

3) What another CAP groups in amber tools should I use for capping of termi
wich would be compatible with gromacs ff?

James

2011/10/31 Justin A. Lemkul <jalemkul at vt.edu>

>
>
> James Starlight wrote:
>
>> exactly
>>
>> all hydrogens were represented as a part of the groupd in wich they were.
>> So is there possible way to explicit it?
>>
>>
> I'm not quite sure I follow.  A united atom force field, by definition,
> does not have aliphatic hydrogens.  I suggest you consult the literature
> and understand the intrinsics of the chosen force field before you attempt
> to use it.
>
>
>  Another question about pdb2gmx. about -term
>> as I understood this must be used only if atoms for CAP groups are
>> presented in the PDB file mustnt it?
>>
>
> The -ter option can be useful in changing the protonation state of the
> termini or in adding caps.
>
>
>  So what term exactly do? Doest it make a connection beetween CAP and C
>> and N termi or anothing else?
>>
>>
> By default, pdb2gmx builds charged termini (NH3+ and COO-), per the
> typical protonation states of these groups in aqueous solution.  If a cap
> is present, then you don't want pdb2gmx doing that, so by choosing "None"
> for the termini in conjunction with the caps, only one proton is added to
> the N-terminus (as in an amide), and only one carbonyl oxygen remains on
> the C-terminus, such that the caps can be added in a sensible way.
>
>
> -Justin
>
> --
> ==============================**==========
>
> Justin A. Lemkul
> Ph.D. Candidate
> ICTAS Doctoral Scholar
> MILES-IGERT Trainee
> Department of Biochemistry
> Virginia Tech
> Blacksburg, VA
> jalemkul[at]vt.edu | (540) 231-9080
> http://www.bevanlab.biochem.**vt.edu/Pages/Personal/justin<http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin>
>
> ==============================**==========
> --
> gmx-users mailing list    gmx-users at gromacs.org
> http://lists.gromacs.org/**mailman/listinfo/gmx-users<http://lists.gromacs.org/mailman/listinfo/gmx-users>
> Please search the archive at http://www.gromacs.org/**
> Support/Mailing_Lists/Search<http://www.gromacs.org/Support/Mailing_Lists/Search>before posting!
> Please don't post (un)subscribe requests to the list. Use the www
> interface or send it to gmx-users-request at gromacs.org.
> Can't post? Read http://www.gromacs.org/**Support/Mailing_Lists<http://www.gromacs.org/Support/Mailing_Lists>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://maillist.sys.kth.se/pipermail/gromacs.org_gmx-users/attachments/20111102/6b93193e/attachment.html>