
<!DOCTYPE html PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">

<html>

<head>

  <title></title>

</head>

<body>

<div class="moz-text-flowed"

 style="font-family: -moz-fixed; font-size: 12px;" lang="x-western">Thanks

Maik, <br>

&nbsp;&nbsp;&nbsp;&nbsp; <br>

But as david sugged me that i have to include topology of my linker in

rtp file , so now i am not using PRODRG server. <br>

I also find the option for fixing the optimal no of Hydrogen in

PRODRG. <br>

It is DELHYD "atom name" flag above the PDB file. <br>

Thanks again for ur suggesion. <br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <br>

&nbsp;&nbsp;&nbsp;&nbsp; &nbsp; Alok Jain <br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <br>

Maik Goette wrote: <br>

<br>

<blockquote type="cite">Hi <br>

  <br>

I had similar problems longer time ago. <br>

PRODRG seems to guess the optimal no of hydrogens. <br>

But you are able to tell the PRODRG-Server to don't place specific

hydrogens. <br>

How it works is listed in a howto on the servers homepage. <br>

I preferred ASCII-Input, so I am not sure how to include this option

into PDB. This may also occur in the "manual". <br>

  <br>

Concerning the charges, one has to pay attention, that the charges may

change significantly, depending on what atom will be attached to the

linkers nitrogen. This might even influence the AA on the opposite side. <br>

  <br>

A possible solvation may be, to change the linker into glycins. <br>

  <br>

The other way would be a QM of the linker, I guess. <br>

  <br>

Regards <br>

  <br>

Maik Goette, Dipl. Biol. <br>

Max Planck Institute for Biophysical Chemistry <br>

Theoretical &amp; computational biophysics department <br>

Am Fassberg 11 <br>

37077 Goettingen <br>

Germany <br>

Tel.&nbsp; : ++49 551 201 2310 <br>

Fax&nbsp;&nbsp; : ++49 551 201 2302 <br>

Email : mgoette[at]mpi-bpc.mpg.de <br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; mgoette2[at]gwdg.de <br>

WWW&nbsp;&nbsp; : <a class="moz-txt-link-freetext"

 href="http://www.mpibpc.gwdg.de/groups/grubmueller/">http://www.mpibpc.gwdg.de/groups/grubmueller/</a> <br>

  <br>

  <br>

Alok wrote: <br>

  <br>

  <blockquote type="cite">hello gmx users, <br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; i want to simulate a peptide-linker-peptide system. The <br>

system is XXX-NH-CH2-CH2-NH-XXX where X represents the <br>

standard amino acids. <br>

&nbsp;To simulate this system i want to first generate a topology file for

the <br>

linker section of the system.I did so by using&nbsp; PRODRG server by

suppling <br>

the linker section as the input.. <br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; The input coordinate file to the server and the

output <br>

from the server are as follows: <br>

    <br>

Input: <br>

HETATM&nbsp;&nbsp; 14&nbsp; N&nbsp;&nbsp;&nbsp;&nbsp; 2&nbsp;&nbsp;&nbsp;&nbsp; 2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; -0.564&nbsp;&nbsp; 0.106&nbsp;&nbsp; 1.392 <br>

HETATM&nbsp;&nbsp; 15&nbsp; C&nbsp;&nbsp;&nbsp;&nbsp; 2&nbsp;&nbsp;&nbsp;&nbsp; 2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.779&nbsp;&nbsp; 0.380&nbsp;&nbsp; 0.903 <br>

HETATM&nbsp;&nbsp; 16&nbsp; C&nbsp;&nbsp;&nbsp;&nbsp; 2&nbsp;&nbsp;&nbsp;&nbsp; 2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1.220&nbsp;&nbsp; 1.824&nbsp;&nbsp; 1.201 <br>

HETATM&nbsp;&nbsp; 17&nbsp; N&nbsp;&nbsp;&nbsp;&nbsp; 2&nbsp;&nbsp;&nbsp;&nbsp; 2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 2.587&nbsp;&nbsp; 2.069&nbsp;&nbsp; 0.782 <br>

    <br>

Output: <br>

    <br>

[ moleculetype ] <br>

; Name nrexcl <br>

&nbsp;_2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 3 <br>

    <br>

[ atoms ] <br>

;&nbsp;&nbsp; nr&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; type&nbsp; resnr resid&nbsp; atom&nbsp; cgnr&nbsp;&nbsp; charge&nbsp;&nbsp;&nbsp;&nbsp; mass <br>

&nbsp;&nbsp;&nbsp;&nbsp; 1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; NL&nbsp;&nbsp;&nbsp;&nbsp; 1&nbsp;&nbsp; _2&nbsp;&nbsp;&nbsp;&nbsp; NAA&nbsp;&nbsp;&nbsp;&nbsp; 1&nbsp;&nbsp;&nbsp; 0.704&nbsp; 14.0067 <br>

&nbsp;&nbsp;&nbsp;&nbsp; 2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; H&nbsp;&nbsp;&nbsp;&nbsp; 1&nbsp;&nbsp; _2&nbsp;&nbsp;&nbsp;&nbsp; HAA&nbsp;&nbsp;&nbsp;&nbsp; 1&nbsp;&nbsp;&nbsp; 0.059&nbsp;&nbsp; 1.0080 <br>

&nbsp;&nbsp;&nbsp;&nbsp; 3&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; H&nbsp;&nbsp;&nbsp;&nbsp; 1&nbsp;&nbsp; _2&nbsp;&nbsp;&nbsp;&nbsp; HAB&nbsp;&nbsp;&nbsp;&nbsp; 1&nbsp;&nbsp;&nbsp; 0.059&nbsp;&nbsp; 1.0080 <br>

&nbsp;&nbsp;&nbsp;&nbsp; 4&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; H&nbsp;&nbsp;&nbsp;&nbsp; 1&nbsp;&nbsp; _2&nbsp;&nbsp;&nbsp;&nbsp; HAC&nbsp;&nbsp;&nbsp;&nbsp; 1&nbsp;&nbsp;&nbsp; 0.059&nbsp;&nbsp; 1.0080 <br>

&nbsp;&nbsp;&nbsp;&nbsp; 5&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; CH2&nbsp;&nbsp;&nbsp;&nbsp; 1&nbsp;&nbsp; _2&nbsp;&nbsp;&nbsp;&nbsp; CAB&nbsp;&nbsp;&nbsp;&nbsp; 1&nbsp;&nbsp;&nbsp; 0.119&nbsp; 14.0270 <br>

&nbsp;&nbsp;&nbsp;&nbsp; 6&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; CH2&nbsp;&nbsp;&nbsp;&nbsp; 1&nbsp;&nbsp; _2&nbsp;&nbsp;&nbsp;&nbsp; CAC&nbsp;&nbsp;&nbsp;&nbsp; 2&nbsp;&nbsp;&nbsp; 0.119&nbsp; 14.0270 <br>

&nbsp;&nbsp;&nbsp;&nbsp; 7&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; NL&nbsp;&nbsp;&nbsp;&nbsp; 1&nbsp;&nbsp; _2&nbsp;&nbsp;&nbsp;&nbsp; NAD&nbsp;&nbsp;&nbsp;&nbsp; 2&nbsp;&nbsp;&nbsp; 0.704&nbsp; 14.0067 <br>

&nbsp;&nbsp;&nbsp;&nbsp; 8&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; H&nbsp;&nbsp;&nbsp;&nbsp; 1&nbsp;&nbsp; _2&nbsp;&nbsp;&nbsp;&nbsp; HAE&nbsp;&nbsp;&nbsp;&nbsp; 2&nbsp;&nbsp;&nbsp; 0.059&nbsp;&nbsp; 1.0080 <br>

&nbsp;&nbsp;&nbsp;&nbsp; 9&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; H&nbsp;&nbsp;&nbsp;&nbsp; 1&nbsp;&nbsp; _2&nbsp;&nbsp;&nbsp;&nbsp; HAF&nbsp;&nbsp;&nbsp;&nbsp; 2&nbsp;&nbsp;&nbsp; 0.059&nbsp;&nbsp; 1.0080 <br>

&nbsp;&nbsp;&nbsp; 10&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; H&nbsp;&nbsp;&nbsp;&nbsp; 1&nbsp;&nbsp; _2&nbsp;&nbsp;&nbsp;&nbsp; HAD&nbsp;&nbsp;&nbsp;&nbsp; 2&nbsp;&nbsp;&nbsp; 0.059&nbsp;&nbsp; 1.0080 <br>

    <br>

The queries to the output are : <br>

    <br>

1)&nbsp;&nbsp; The linker now has three hydrogens attached to nitrogen whereas

the linker with the peptide will have only one hydrogen attached to <br>

nitrogen.So can we restrict the number of hydrogens to one or should i

now <br>

delete two hydrogens from the topology file before using it for my <br>

system???? <br>

    <br>

2)&nbsp;&nbsp; Should I use the same parameters including the charges directly as <br>

given in the output in the topology file for my simulation or are there <br>

other issues to look into before using this topology file????? <br>

    <br>

3)&nbsp; These parameters were developed by the server for an isolated

linker <br>

(N-C-C-N).So can i use the same parameters for the linker when it is

bound <br>

to my peptide (XXX-Linker-XXX)???? <br>

    <br>

Please shed some light into this issue. <br>

    <br>

Thanking all in advance. <br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Alok Jain <br>

    <br>

    <br>

    <br>

_______________________________________________ <br>

gmx-users mailing list <br>

    <a class="moz-txt-link-abbreviated"

 href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a> <br>

    <a class="moz-txt-link-freetext"

 href="http://www.gromacs.org/mailman/listinfo/gmx-users">http://www.gromacs.org/mailman/listinfo/gmx-users</a> <br>

Please don't post (un)subscribe requests to the list. Use the www

interface or send it to <a class="moz-txt-link-abbreviated"

 href="mailto:gmx-users-request@gromacs.org">gmx-users-request@gromacs.org</a>.<br>

    <br>

. <br>

    <br>

  </blockquote>

_______________________________________________ <br>

gmx-users mailing list <br>

  <a class="moz-txt-link-abbreviated"

 href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a> <br>

  <a class="moz-txt-link-freetext"

 href="http://www.gromacs.org/mailman/listinfo/gmx-users">http://www.gromacs.org/mailman/listinfo/gmx-users</a> <br>

Please don't post (un)subscribe requests to the list. Use the www

interface or send it to <a class="moz-txt-link-abbreviated"

 href="mailto:gmx-users-request@gromacs.org">gmx-users-request@gromacs.org</a>.<br>

</blockquote>

<br>

<br>

</div>

</body>

</html>