<html><head><style type="text/css"><!-- DIV {margin:0px;} --></style></head><body><div style="font-family:times new roman,new york,times,serif;font-size:12pt">What are the suitable options for&nbsp; Chirality&nbsp; , Full charges ,&nbsp; Energy
Minimization in PRODRG&nbsp; to get input files&nbsp; for hydrocarbon molecules ?
And If I choose Full charges : YES , I will get the files that include
charge info. Can I erase the charge info and&nbsp; write the charge info I
got from Gaussian ( net charge is zero ) according to the order of
atoms that I paste to the PRODRG screen ? Could you give me some
advices about this ?<br>
<br>
Thanks in advance<br>
<br>
Adrien<div style="font-family: times new roman,new york,times,serif; font-size: 12pt;"><br><div style="font-family: arial,helvetica,sans-serif; font-size: 13px;"><font face="Tahoma" size="2"><hr size="1"><b><span style="font-weight: bold;">From:</span></b> Justin A. Lemkul &lt;jalemkul@vt.edu&gt;<br><b><span style="font-weight: bold;">To:</span></b> adriendelmont@yahoo.com; Discussion list for GROMACS users &lt;gmx-users@gromacs.org&gt;<br><b><span style="font-weight: bold;">Sent:</span></b> Thursday, March 12, 2009 12:17:57 AM<br><b><span style="font-weight: bold;">Subject:</span></b> Re: [gmx-users] get a box of 100 molecules<br></font><br>
<br><br>Adrien Delmont wrote:<br>&gt; If I create an .rtp file for my molecule how can I use this file when I producing gro and top files and Can I use this .rtp file with oplsaa and gromos 96 force fields for two different md run successfuly ? Because I plan to run md simulations for all atom and united atom force fields and compare the results. Secondly If I use PRODRG , the only force field I can use is oplsaa. Is it true ? Can't I use a gromos 96 force field ?<br>&gt; <br><br>I think you should start by reading the manual about .rtp files and how they are used.&nbsp; Each force field has its own .rtp, so no, you cannot use the same .rtp file with different force fields.&nbsp; Make a local copy (in your working directory) if you choose to modify your .rtp files; if you are unsure of what you are doing you can seriously screw things up.<br><br>PRODRG produces Gromos-compatible force fields, not OPLS; this is answered in the FAQ at the PRODRG
 site.<br><br>I reiterate: read the manual, be comfortable with how the file types are used, then proceed.<br><br>-Justin<br><br>&gt; --- On *Wed, 3/11/09, Justin A. Lemkul /&lt;<a ymailto="mailto:jalemkul@vt.edu" href="mailto:jalemkul@vt.edu">jalemkul@vt.edu</a>&gt;/* wrote:<br>&gt; <br>&gt;&nbsp; &nbsp;  From: Justin A. Lemkul &lt;<a ymailto="mailto:jalemkul@vt.edu" href="mailto:jalemkul@vt.edu">jalemkul@vt.edu</a>&gt;<br>&gt;&nbsp; &nbsp;  Subject: Re: [gmx-users] get a box of 100 molecules<br>&gt;&nbsp; &nbsp;  To: "Gromacs Users' List" &lt;<a ymailto="mailto:gmx-users@gromacs.org" href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a>&gt;<br>&gt;&nbsp; &nbsp;  Date: Wednesday, March 11, 2009, 9:53 PM<br>&gt; <br>&gt;&nbsp; &nbsp;  Adrien Delmont wrote:<br>&gt;&nbsp; &nbsp;  &gt; The molecule I use is a simple hydrocarbon ( n-heptane ) and I think that<br>&gt;&nbsp; &nbsp;  I should add something like&nbsp; dihedral, bond, angle parameters
 after I search the<br>&gt;&nbsp; &nbsp;  literature. I found some information about that. I think I will add parameters<br>&gt;&nbsp; &nbsp;  into the related force fields .itp file. Is it true ? And then I will create a<br>&gt;&nbsp; &nbsp;  box of 100 molecules without any error.<br>&gt;&nbsp; &nbsp;  &gt; <br>&gt;&nbsp; &nbsp;  Most (if not all) of the parameters you will need will already be present in<br>&gt;&nbsp; &nbsp;  the force field.&nbsp; No modification of any of these files should be necessary..&nbsp; If<br>&gt;&nbsp; &nbsp;  you want to use pdb2gmx to create your topology, you need to create an .rtp<br>&gt;&nbsp; &nbsp;  entry. For a simple alkane (especially with a UA representation), it should be<br>&gt;&nbsp; &nbsp;  trivial to write the topology by hand with a text editor and a thorough<br>&gt;&nbsp; &nbsp;  understanding of Chapter 5.<br>&gt; <br>&gt;&nbsp; &nbsp;  Alternate ideas:<br>&gt; <br>&gt;&nbsp; &nbsp;  1. Use PRODRG to
 generate a rough topology. Verify its contents before<br>&gt;&nbsp; &nbsp;  proceeding, make any modifications you need, and validate the resulting<br>&gt;&nbsp; &nbsp;  topology.<br>&gt;&nbsp; &nbsp;  2. Don't use Gromos96.&nbsp; AFAIK, none of the Gromos variants deal very well<br>&gt;&nbsp; &nbsp;  with hydrocarbon chains.&nbsp; OPLS might be better, but I think there are some<br>&gt;&nbsp; &nbsp;  limitations there, too.<br>&gt; <br>&gt;&nbsp; &nbsp;  -Justin<br>&gt; <br>&gt;&nbsp; &nbsp;  &gt; --- On *Wed, 3/11/09, Justin A. Lemkul /&lt;<a ymailto="mailto:jalemkul@vt.edu" href="mailto:jalemkul@vt.edu">jalemkul@vt.edu</a>&gt;/* wrote:<br>&gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  From: Justin A. Lemkul &lt;<a ymailto="mailto:jalemkul@vt.edu" href="mailto:jalemkul@vt.edu">jalemkul@vt.edu</a>&gt;<br>&gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  Subject: Re: [gmx-users] get a box of 100 molecules<br>&gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  To: <a
 ymailto="mailto:adriendelmont@yahoo..com" href="mailto:adriendelmont@yahoo..com">adriendelmont@yahoo..com</a>, "Discussion list for GROMACS<br>&gt;&nbsp; &nbsp;  users"<br>&gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  &lt;<a ymailto="mailto:gmx-users@gromacs.org" href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a>&gt;<br>&gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  Date: Wednesday, March 11, 2009, 9:11 PM<br>&gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  Adrien Delmont wrote:<br>&gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  &gt;&nbsp; I just want to add missing parameters for my molecule in order<br>&gt;&nbsp; &nbsp;  to get<br>&gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  ...gro and top files for example GROMOS96 43a2 force field . How can I<br>&gt;&nbsp; &nbsp;  add<br>&gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  missing parameters ? I didn't find any information about this<br>&gt;&nbsp; &nbsp;  procedure.<br>&gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  I'm waiting for your
 help.<br>&gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp; &nbsp; &gt;&nbsp; &nbsp;  You'll have to define what you mean by "parameters" - is<br>&gt;&nbsp; &nbsp;  this<br>&gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  some new small molecule that has not been parameterized under this<br>&gt;&nbsp; &nbsp;  force field?&nbsp; &nbsp;  Are you starting from scratch?&nbsp; If so, this is not a trivial<br>&gt;&nbsp; &nbsp;  problem for a<br>&gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  beginner. Read here:<br>&gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  <a href="http://wiki.gromacs.org/index.php/Parameterization" target="_blank">http://wiki.gromacs.org/index.php/Parameterization</a><br>&gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  If you have adequate parameters for most of your molecule and are<br>&gt;&nbsp; &nbsp;  missing<br>&gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  something like a dihedral, etc. then read the primary literature and<br>&gt;&nbsp;
 &nbsp;  derive it<br>&gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  yourself.&nbsp; This process is also complicated, but is probably less work<br>&gt;&nbsp; &nbsp;  than<br>&gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  coming up with parameters for a brand new molecule.<br>&gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  -Justin<br>&gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  &gt; Thanks in advance<br>&gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  &gt; --- On *Wed, 3/11/09, Mark Abraham<br>&gt;&nbsp; &nbsp;  /&lt;<a ymailto="mailto:Mark.Abraham@anu.edu.au" href="mailto:Mark.Abraham@anu.edu.au">Mark.Abraham@anu.edu.au</a>&gt;/*<br>&gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  wrote:<br>&gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  From: Mark Abraham &lt;<a ymailto="mailto:Mark.Abraham@anu.edu.au" href="mailto:Mark.Abraham@anu.edu.au">Mark.Abraham@anu.edu.au</a>&gt;<br>&gt;&nbsp; &nbsp;  &gt;&nbsp;
 &nbsp;  &gt;&nbsp; &nbsp;  Subject: Re: [gmx-users] get a box of 100 molecules<br>&gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  To: "Discussion list for GROMACS users"<br>&gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  &lt;<a ymailto="mailto:gmx-users@gromacs.org" href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a>&gt;<br>&gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  Date: Wednesday, March 11, 2009, 4:58 PM<br>&gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  Adrien Delmont wrote:<br>&gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  &gt;<br>&gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp; &nbsp; &nbsp; ------------------------------------------------------------------------<br>&gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  &gt; *From:* Adrien Delmont &lt;<a ymailto="mailto:adriendelmont@yahoo.com"
 href="mailto:adriendelmont@yahoo.com">adriendelmont@yahoo.com</a>&gt;<br>&gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  &gt; *To:* <a ymailto="mailto:jalemkul@vt.edu" href="mailto:jalemkul@vt.edu">jalemkul@vt.edu</a><br>&gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  &gt; *Sent:* Wednesday, March 11, 2009 6:09:53 PM<br>&gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  &gt; *Subject:* Re: [gmx-users] get a box of 100 molecules<br>&gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  &gt; I didn't get gro and top files with any<br>&gt;&nbsp; &nbsp;  force field<br>&gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  in Gromacs. I<br>&gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  basicly want to use oplsaa and GROMOS96 43a2 force field<br>&gt;&nbsp; &nbsp;  (improved<br>&gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  alkane<br>&gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  dihedrals) . I think I should add all the
 parameters but&nbsp; how<br>&gt;&nbsp; &nbsp;  can I<br>&gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  add my<br>&gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  parametes both into oplsaa and GROMOS96 43a2 force field ?<br>&gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  You can't just mix force fields. See<br>&gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  <a href="http://wiki.gromacs.org/index.php/force_field" target="_blank">http://wiki.gromacs.org/index.php/force_field</a><br>&gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  Mark<br>&gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  _______________________________________________<br>&gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  gmx-users mailing list&nbsp; &nbsp; <a ymailto="mailto:gmx-users@gromacs.org" href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a><br>&gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  &gt;&nbsp;
 &nbsp;  <a href="http://www.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://www.gromacs.org/mailman/listinfo/gmx-users</a><br>&gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  Please search the archive at <a href="http://www.gromacs.org/search" target="_blank">http://www.gromacs.org/search</a><br>&gt;&nbsp; &nbsp;  before<br>&gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  posting!<br>&gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  Please don't post (un)subscribe requests to the list. Use<br>&gt;&nbsp; &nbsp;  the www<br>&gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  interface<br>&gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  or send it to <a ymailto="mailto:gmx-users-request@gromacs.org" href="mailto:gmx-users-request@gromacs.org">gmx-users-request@gromacs.org</a>.<br>&gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  Can't post? Read<br>&gt;&nbsp; &nbsp;  <a href="http://www.gromacs.org/mailing_lists/users.php"
 target="_blank">http://www.gromacs.org/mailing_lists/users.php</a><br>&gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  &gt;<br>&gt;&nbsp; &nbsp;  ------------------------------------------------------------------------<br>&gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  &gt; _______________________________________________<br>&gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  &gt; gmx-users mailing list&nbsp; &nbsp; <a ymailto="mailto:gmx-users@gromacs.org" href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a><br>&gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  &gt; <a href="http://www.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://www.gromacs.org/mailman/listinfo/gmx-users</a><br>&gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  &gt; Please search the archive at <a href="http://www.gromacs.org/search" target="_blank">http://www.gromacs.org/search</a> before<br>&gt;&nbsp; &nbsp;  posting!<br>&gt;&nbsp; &nbsp; 
 &gt;&nbsp; &nbsp;  &gt; Please don't post (un)subscribe requests to the list. Use the<br>&gt;&nbsp; &nbsp;  www<br>&gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  interface or send it to <a ymailto="mailto:gmx-users-request@gromacs.org" href="mailto:gmx-users-request@gromacs.org">gmx-users-request@gromacs.org</a>.<br>&gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  &gt; Can't post? Read<br>&gt;&nbsp; &nbsp;  <a href="http://www.gromacs.org/mailing_lists/users.php" target="_blank">http://www.gromacs.org/mailing_lists/users.php</a><br>&gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  -- ========================================<br>&gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  Justin A. Lemkul<br>&gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  Graduate Research Assistant<br>&gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  ICTAS Doctoral Scholar<br>&gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  Department of Biochemistry<br>&gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  Virginia Tech<br>&gt;&nbsp;
 &nbsp;  &gt;&nbsp; &nbsp;  Blacksburg, VA<br>&gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  jalemkul[at]<a target="_blank" href="http://vt.edu">vt.edu</a> | (540) 231-9080<br>&gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  <a href="http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin" target="_blank">http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin</a><br>&gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  ========================================<br>&gt;&nbsp; &nbsp;  &gt;&nbsp; &nbsp;  &gt; <br>&gt;&nbsp; &nbsp;  -- ========================================<br>&gt; <br>&gt;&nbsp; &nbsp;  Justin A. Lemkul<br>&gt;&nbsp; &nbsp;  Graduate Research Assistant<br>&gt;&nbsp; &nbsp;  ICTAS Doctoral Scholar<br>&gt;&nbsp; &nbsp;  Department of Biochemistry<br>&gt;&nbsp; &nbsp;  Virginia Tech<br>&gt;&nbsp; &nbsp;  Blacksburg, VA<br>&gt;&nbsp; &nbsp;  jalemkul[at]vt.edu | (540) 231-9080<br>&gt;&nbsp; &nbsp;  <a
 href="http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin" target="_blank">http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin</a><br>&gt; <br>&gt;&nbsp; &nbsp;  ========================================<br>&gt;&nbsp; &nbsp;  _______________________________________________<br>&gt;&nbsp; &nbsp;  gmx-users mailing list&nbsp; &nbsp; <a ymailto="mailto:gmx-users@gromacs.org" href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a><br>&gt;&nbsp; &nbsp;  <a href="http://www.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://www.gromacs.org/mailman/listinfo/gmx-users</a><br>&gt;&nbsp; &nbsp;  Please search the archive at <a href="http://www.gromacs.org/search" target="_blank">http://www.gromacs.org/search</a> before posting!<br>&gt;&nbsp; &nbsp;  Please don't post (un)subscribe requests to the list. Use the www interface<br>&gt;&nbsp; &nbsp;  or send it to <a ymailto="mailto:gmx-users-request@gromacs.org"
 href="mailto:gmx-users-request@gromacs.org">gmx-users-request@gromacs.org</a>.<br>&gt;&nbsp; &nbsp;  Can't post? Read <a href="http://www.gromacs.org/mailing_lists/users.php" target="_blank">http://www.gromacs.org/mailing_lists/users.php</a><br>&gt; <br>&gt; <br>&gt; <br>&gt; ------------------------------------------------------------------------<br>&gt; <br>&gt; _______________________________________________<br>&gt; gmx-users mailing list&nbsp; &nbsp; <a ymailto="mailto:gmx-users@gromacs.org" href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a><br>&gt; <a href="http://www.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://www.gromacs.org/mailman/listinfo/gmx-users</a><br>&gt; Please search the archive at <a href="http://www.gromacs.org/search" target="_blank">http://www.gromacs.org/search</a> before posting!<br>&gt; Please don't post (un)subscribe requests to the list. Use the www interface or send it to <a
 ymailto="mailto:gmx-users-request@gromacs.org" href="mailto:gmx-users-request@gromacs.org">gmx-users-request@gromacs.org</a>.<br>&gt; Can't post? Read <a href="http://www.gromacs.org/mailing_lists/users.php" target="_blank">http://www.gromacs.org/mailing_lists/users.php</a><br><br>-- ========================================<br><br>Justin A. Lemkul<br>Graduate Research Assistant<br>ICTAS Doctoral Scholar<br>Department of Biochemistry<br>Virginia Tech<br>Blacksburg, VA<br>jalemkul[at]vt.edu | (540) 231-9080<br><a href="http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin" target="_blank">http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin</a><br><br>========================================<br></div></div></div><br>

      </body></html>