Thanks for your guidance.<br>I looked that file,But I think the name of functional groups are different in .rtp file because I can't find no one of them in this file.<br><br>please let me know how can I know the correct name f or functional groups<br>
for example:HYDROXYL,CARBOXYL,HALO,AMINO and ...<br><br>Thanks in advance for your help<br><br><div class="gmail_quote">On Mon, Jan 24, 2011 at 4:01 PM, Justin A. Lemkul <span dir="ltr"><<a href="mailto:jalemkul@vt.edu">jalemkul@vt.edu</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;"><div class="im"><br>
<br>
mohsen ramezanpour wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Dear Dr.justin<br>
Actually by doing this we are using two different force fields in one simulation.<br>
I had done it before and the result was that I discussed before in gmx-users(LINCS Error,Exploding system,Bad contacts between atoms)<br>
Then,this approch seems to doesn't work about my system.<br>
<br>
Then I want to find charges and charge groups for gromos 43A1 and replace them for my drug(to edit PRODRG file manually) and work totally in gromos 43A1.<br>
</blockquote>
<br></div>
Please keep your story consistent. In the last message, you said you wanted to work completely within 53A6, so I advised you on how to do that, now you say that you're trying to work completely within 43A1.<div class="im">
<br>
<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Unfortunately I can't obtain these parameter.<br>
</blockquote>
<br></div>
You certainly do have these parameters. 43A1 is part of the Gromacs installation; in the .rtp file you'll find all of the functional groups that were derived in 43A1, as applied to amino acids and a few other groups.<div class="im">
<br>
<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Please let me have if you have it.<br>
Can I use some Ab Initio software for determining partial charges of my drug?<br>
for example ABINIT or Gaussian!<br>
<br>
</blockquote>
<br></div>
My paper that you said you read has discussion and recommendations on this point. But be very clear: none of the QM methods we tested were able to reproduce the charges that are assigned to known functional groups since the Gromos parameterization methodology calls for "empirical refinement." Thus, manual modification and thorough validation are always necessary.<br>
<br>
-Justin<br>
<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Thanks in advance<div class="im"><br>
<br>
On Sat, Jan 22, 2011 at 8:03 PM, Justin A. Lemkul <<a href="mailto:jalemkul@vt.edu" target="_blank">jalemkul@vt.edu</a> <mailto:<a href="mailto:jalemkul@vt.edu" target="_blank">jalemkul@vt.edu</a>>> wrote:<br>
<br>
<br>
<br>
mohsen ramezanpour wrote:<br>
<br>
Ok<br>
then,I can use PRODRG server to generate .top and .gro files<br>
for drug.<br>
since it's reported charges are not very accurate ,we can<br>
replace all charges completely with them in 53A6(if was present).<br>
But it means we are working in 53A6 force field.<br>
then,we must generate .top and .gro files for our protein with<br>
53A6 too.<br>
and work completely with 53A6.<br>
Am i right?<br>
thanks in advance<br>
<br>
<br>
That sounds like a reasonable approach. Be sure to validate the<br>
drug topology. In my experience, this procedure is pretty good, but<br>
you always have to convince reviewers...<br>
<br>
-Justin<br>
<br>
<br>
On Sat, Jan 22, 2011 at 4:43 PM, Justin A. Lemkul<br>
<<a href="mailto:jalemkul@vt.edu" target="_blank">jalemkul@vt.edu</a> <mailto:<a href="mailto:jalemkul@vt.edu" target="_blank">jalemkul@vt.edu</a>><br></div><div><div></div><div class="h5">
<mailto:<a href="mailto:jalemkul@vt.edu" target="_blank">jalemkul@vt.edu</a> <mailto:<a href="mailto:jalemkul@vt.edu" target="_blank">jalemkul@vt.edu</a>>>> wrote:<br>
<br>
<br>
<br>
mohsen ramezanpour wrote:<br>
<br>
Dear Justin<br>
<br>
I read your articles about PRODRG server,they were very<br>
useful.<br>
But I have a question:<br>
are charges of functional groups and generally other atom<br>
groups<br>
the same in all force fields?<br>
Because you have modified charges of your molecules by<br>
Gromos96<br>
53A6 while prodrg server is generating topology files in<br>
43A1.<br>
I want to know can I replace charges from gromos 53A6 or<br>
other<br>
forcefields?<br>
thanks in advance<br>
<br>
<br>
Charges are not the same between force fields. We did our study<br>
with 43A1 since that is what PRODRG purports to produce. I would<br>
say that our recommendations carry to other Gromos force<br>
fields, as<br>
well, but don't take charges from 43A1 and apply them to<br>
53A6. Be<br>
consistent within the force field.<br>
<br>
The atom types produced by PRODRG are largely shared between 43A1<br>
and 53A6, so if you *completely* replace all charges with<br>
those from<br>
53A6, you should have a topology that is compatible with 53A6.<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></div></div>
jalemkul[at]<a href="http://vt.edu" target="_blank">vt.edu</a> <<a href="http://vt.edu" target="_blank">http://vt.edu</a>> <<a href="http://vt.edu" target="_blank">http://vt.edu</a>> | (540)<div class="im">
<br>
231-9080<br>
<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>
-- gmx-users mailing list <a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br>
<mailto:<a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a>><br></div>
<mailto:<a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a> <mailto:<a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a>>><div><div></div><div class="h5">
<br>
<br>
<a href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/mailman/listinfo/gmx-users</a><br>
Please search the archive at<br>
<a href="http://www.gromacs.org/Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/Support/Mailing_Lists/Search</a> before<br>
posting!<br>
Please don't post (un)subscribe requests to the list. Use the www<br>
interface or send it to <a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a><br>
<mailto:<a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a>><br>
<mailto:<a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a><br>
<mailto:<a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a>>>.<br>
<br>
Can't post? Read <a href="http://www.gromacs.org/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/Support/Mailing_Lists</a><br>
<br>
<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 href="http://vt.edu" target="_blank">vt.edu</a> <<a href="http://vt.edu" target="_blank">http://vt.edu</a>> | (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>
-- gmx-users mailing list <a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br>
<mailto:<a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a>><br>
<a href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/mailman/listinfo/gmx-users</a><br>
Please search the archive at<br>
<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<br>
interface or send it to <a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a><br>
<mailto:<a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a>>.<br>
Can't post? Read <a href="http://www.gromacs.org/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/Support/Mailing_Lists</a><br>
<br>
<br>
</div></div></blockquote><div><div></div><div class="h5">
<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 href="http://vt.edu" target="_blank">vt.edu</a> | (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>
-- <br>
gmx-users mailing list <a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br>
<a href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/mailman/listinfo/gmx-users</a><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 href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a>.<br>
Can't post? Read <a href="http://www.gromacs.org/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/Support/Mailing_Lists</a><br>
</div></div></blockquote></div><br>