<html><body><div style="color:#000; background-color:#fff; font-family:arial, helvetica, sans-serif;font-size:12pt"><div><span>OK, I have to build a .rtp file. But I'd like to know if I should make it corresponding to relation between.itp and .rtp files format?</span></div><div><br><span></span></div><div><span>Cheers,</span></div><div><span>Shima<br></span></div><div><br></div> <div style="font-family: arial, helvetica, sans-serif; font-size: 12pt;"> <div style="font-family: times new roman, new york, times, serif; font-size: 12pt;"> <div dir="ltr"> <font face="Arial" size="2"> <hr size="1"> <b><span style="font-weight:bold;">From:</span></b> Justin A. Lemkul <jalemkul@vt.edu><br> <b><span style="font-weight: bold;">To:</span></b> Discussion list for GROMACS users <gmx-users@gromacs.org> <br> <b><span style="font-weight: bold;">Sent:</span></b> Wednesday, May 9, 2012 3:58 PM<br> <b><span style="font-weight: bold;">Subject:</span></b> Re:
[gmx-users] Formyl parameters<br> </font> </div> <br>
<br><br>On 5/9/12 4:48 AM, Shima Arasteh wrote:<br>> <br>> Oh my god! :-)<br>> The situation seems going complicated! :-(<br>> Let me explain what I'm gonna do in detail:<br>> I have formyl residue as the N-terminus in my protein (.pdb file) ! As you know,<br>> the topology database in CHARMM27 doesn't have formyl. So I got its parameter as<br>> an .itp file through swissparam site. Now I need to include the parameters in<br>> the top file. But don't know it well! I just did such as this:<br>> <br>> 1. seperating the formyl residue from my .pdb file.<br>> 2. run pdb2gmx to get top file ( for the .pdb file which doesn't have formyl )<br>> 3. include the formyl .itp file in my top file<br>> <br>> <br>> Now please let me know if I've done correctly?<br>> <br><br>No, you haven't. What you've defined now are two separate molecules, when in fact the formyl group needs to be an added moiety to your
protein.<br><br>Dealing with this species can be done in one of two ways, both of which require an .rtp file. Whether or not making a formyl .rtp entry is appropriate is up to you based on the parameterization you have already done. I would be very surprised if the presence of the formyl group did not change some aspect (e.g. charges) of the amino acid to which it is connected. In this case, it is better in my mind to create an .rtp entry for the complete formylated amino acid, change residue and/or atom names as necessary in the .pdb file, and run the entire structure (with the formyl group) through pdb2gmx. The general workflow for doing all of this (in addition to information in the manual) can be found at:<br><br>http://www.gromacs.org/Documentation/How-tos/Adding_a_Residue_to_a_Force_Field<br><br>-Justin<br><br>-- ========================================<br><br>Justin A. Lemkul<br>Ph.D. Candidate<br>ICTAS Doctoral
Scholar<br>MILES-IGERT Trainee<br>Department of Biochemistry<br>Virginia Tech<br>Blacksburg, VA<br>jalemkul[at]<a target="_blank" href="http://vt.edu">vt.edu</a> | (540) 231-9080<br>http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin<br><br>========================================<br>-- gmx-users mailing list <a ymailto="mailto:gmx-users@gromacs.org" href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a><br>http://lists.gromacs.org/mailman/listinfo/gmx-users<br>Please search the archive at <a href="http://www.gromacs.org/Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/Support/Mailing_Lists/Search</a> before posting!<br>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>Can't post? Read
http://www.gromacs.org/Support/Mailing_Lists<br><br><br> </div> </div> </div></body></html>