<html>
<head>
<style><!--
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 10pt;
font-family:Verdana
}
--></style>
</head>
<body class='hmmessage'>
<span class="Apple-style-span" style="font-family:Tahoma, Verdana, Arial, sans-serif;color:rgb(68, 68, 68)"><pre style="white-space:normal">Hi,</pre><pre style="white-space:normal">Thank you for your help.Now there is this question that I have just .pdb file and when use protonate command it is" protonate -f drg.pdb -o drg.gro", this is without hydrogens atoms too.</pre><pre style="white-space:normal">I think something is wrong, but I don't know what it is.In definition of "protonate" there is "<span class="Apple-style-span" style="color:rgb(0, 0, 0);font-family:Times, Verdana, Arial, sans-serif;font-size:16px"><tt style="font-family:courier, 'lucida console', serif">protonate</tt>&nbsp;reads (a) conformation(s) and adds all missing hydrogens as defined in&nbsp;<tt style="font-family:courier, 'lucida console', serif">ffgmx2.<a href="http://rocks5.vki.ac.be/gromacs/online/hdb.html" style="text-decoration:none">hdb</a></tt>."&nbsp;<span class="Apple-style-span" style="font-fami
 ly:monospace, Verdana, Arial, sans-serif;font-size:13px;color:rgb(68, 68, 68)">but I can't add hydrogens. What is my problem?</span></span></pre><pre style="white-space:normal">Thanks</pre><pre style="white-space:normal">you zou wrote:<br>&gt; Hi again,<br>&gt;&nbsp;<br>&gt; Sorry, I have one question now, what is the meaning of structure? I think<br>&gt; coordinates is structure, is it true?<br>&gt;&nbsp;<br><br>Yes, a coordinate file contains a structure.<br><br>&gt; If it is true, when I used "editconf -f drg.pdb -o drg.gro" number of atoms<br>&gt; are different from top file and editconf can not add hydrogens to drg.gro. If<br>&gt; Gromacs can handle .pdb, How can it do this, because number of atoms are<br>&gt; different(Which command I have to use?). If can't handle it how can I add<br>&gt; Hydrogens to drg.gro?<br>&gt;&nbsp;<br><br>The underlying assumption when running any simulation is that you have developed&nbsp;<br>the proper parameters for the ligand and that it 
 has an appropriate structure.&nbsp;<br>If you need additional hydrogens, the Gromacs "protonate" tool can generate an&nbsp;<br>all-atom structure.<br><br>-Justin<br><br>&gt; Thanks,<br>&gt;&nbsp;<br>&gt;&nbsp;<br>&gt;&nbsp;<br>&gt; you zou wrote:<br>&gt;&gt; Hi again,<br>&gt;&gt;&nbsp;<br>&gt;&gt; Sorry I confused you with my question. My question is How can I make .gro&nbsp;<br>&gt;&gt; file and .top file from<br>&gt; drug.pdb (that removed from drug-enzyme.pdb)?<br>&gt;&gt;&nbsp;<br>&gt;&gt; If I can use x2top command I will make .top file just, is it true? I think&nbsp;<br>&gt;&gt; .gro file is dependent on forcefiled too so If I use editconf command I<br>&gt;&gt; will miss something, is it true?<br>&gt;&nbsp;<br>&gt; If you want to use x2top, the assumption is that the structure is already&nbsp;<br>&gt; appropriate as is, that is it is properly protonated. The only tool that is&nbsp;<br>&gt; smart enough to add force field-specific hydrogens is pdb2gmx. If you're<br>&gt;
  using OPLS-AA, then you should have all hydrogens present, anyway. If that's<br>&gt; true, then you can use editconf to create a .gro file (which is not<br>&gt; absolutely necessary; Gromacs can handle .pdb files just fine). If you don't<br>&gt; have all the appropriate atoms present in your molecule's structure, then you<br>&gt; need to build a proper structure.<br>&gt;&nbsp;<br>&gt; -Justin<br>&gt;&nbsp;<br>&gt;&gt;&nbsp;<br>&gt;&gt; Thank you again<br>&gt;&gt;&nbsp;<br>&gt;&gt;&nbsp;<br>&gt;&gt; you<br>&gt; zou wrote:<br>&gt;&gt;&gt; Hi Justin,<br>&gt;&gt;&gt;&nbsp;<br>&gt;&gt;&gt; Thank you for your help, But when I run x2top command there is one error<br>&gt;&gt;&gt; that is: " Can not find forcefield for atom C1-1 with 2 bonds Can not<br>&gt;&gt;&gt; find forcefield for atom C4-4 with 2 bonds ...<br>&gt;&gt; &amp;g t; Program x2top, VERSION 4.0.5<br>&gt;&gt;&gt; Source code file: x2top.c, line: 207<br>&gt;&gt;&gt;&nbsp;<br>&gt;&gt;&gt; Fatal error: Could only find a f
 orcefield type for 6 out of 24 atoms"<br>&gt;&gt;&gt;&nbsp;<br>&gt;&gt;&nbsp;<br>&gt;&gt; Not all of your atom types are described by ffop
lsaa.n2t so you will have<br>&gt;&gt; to add them. There are only a limited number of types that are covered by&nbsp;<br>&gt;&gt; default.<br>&gt;&gt;&nbsp;<br>&gt;&gt;&nbsp;<a href="http://www.gromacs.org/Documentation/File_Formats/.n2t_File" style="font-weight:inherit;text-decoration:underline;color:rgb(0, 104, 207);cursor:pointer">http://www.gromacs.org/Documentation/File_Formats/.n2t_File</a><br>&gt;&gt;&nbsp;<br>&gt;&gt;&gt; I don't know how<br>&gt; can I adjust this error. I have one more question again,<br>&gt;&gt;&gt; this command give me a top file, if I want gro file of this pdb (drug<br>&gt;&gt;&gt; that has removed from drug-enzyme complex) how can I do that?<br>&gt;&gt;&gt;&nbsp;<br>&gt;&gt;&nbsp;<br>&gt;&gt; Do you just need a .gro file, and not a .top? My understanding from your&nbsp;<br>&gt;&gt; first message was that you needed a topo logy. If you just need a .gro,<br>&gt;&gt; then simply pass your .pdb file to editconf.<br>&gt;&gt;&nbsp;<br>&gt;&gt; -Justin<
 br>&gt;&gt;&nbsp;<br>&gt;&gt;&gt; you zou wrote:<br>&gt;&gt;&gt;&gt; Dear Users,<br>&gt;&gt;&gt;&gt;&nbsp;<br>&gt;&gt;&gt;&gt; I have one question about Drug-Enzyme Complex,Similar to tutorial If I<br>&gt;&gt;&gt;&gt;&nbsp;<br>&gt;&gt;&gt; want to use GROMOS96 43a1, I can use "Prodrg Beta version" for drug<br>&gt;&gt;&gt;&gt; but If I want to use OPLS-AA/L all-atom force field I can use "Prodrg&nbsp;<br>&gt;&gt;&gt;&gt; Beta version" server too, or not?<br>&gt;&gt;&gt;&nbsp;<br>&gt;&gt;&gt; No. You can't use two different force fields in<br>&gt; one simulation system.<br>&gt;&gt;&gt;&nbsp;<br>&gt;&gt;&gt;&gt; If I can't use this server, how can I make .gro file and .itp file for<br>&gt;&gt;&gt;&gt;&nbsp;<br>&gt; &amp; gt;&gt;&gt; drug that remove from initial .pdb file?<br>&gt;&gt;&gt;&gt;&nbsp;<br>&gt;&gt;&gt;&nbsp;<br>&gt;&gt;&gt; There are several programs in the User Contributions from the website,&nbsp;<br>&gt;&gt;&gt; x2top (which is distributed with Gromacs), or you c
 an build the topology<br>&gt;&gt;&gt; by hand. No matter what you choose, you ne<br>&gt;&gt; ed a thorough understanding of the mechanics of<br>&gt;&gt;&gt; your chosen force field, methods of validation, and of course Chapter 5<br>&gt;&gt;&gt; in the<br>&gt;&gt;&gt;&nbsp;<br>&gt;&gt;&gt; Gromacs manual.<br>&gt;&gt;&gt;&nbsp;</pre></span>
                                               <br /><hr />Your E-mail and More On-the-Go. Get Windows Live Hotmail Free. <a href='https://signup.live.com/signup.aspx?id=60969' target='_new'>Sign up now.</a></body>
</html>