<div>Hi!</div>
<div>Thanks for the reply. So would Gromos96 53a6 be available in the latest version of GROMACS? I am currently using Gromacs version 3.3.3.</div>
<div>Thanks</div>
<div>JJ<br><br></div>
<div class="gmail_quote">On Wed, Jan 7, 2009 at 4:28 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="PADDING-LEFT: 1ex; MARGIN: 0px 0px 0px 0.8ex; BORDER-LEFT: #ccc 1px solid"><br><br>jayant james wrote:<br>
<blockquote class="gmail_quote" style="PADDING-LEFT: 1ex; MARGIN: 0px 0px 0px 0.8ex; BORDER-LEFT: #ccc 1px solid">Hi!
<div class="Ih2E3d"><br>Thanks for the reply.<br>I am using the option 0 of pdb2gmx which corresponds to GROMOS96 43a1 force field. But earlier today I was linked to a thesis which used GROMACS force field so I am still wondering if it would be right for me to go to the *.top file and change the charges to amino acids for simulations in GROMOS96 force field found in GROMACS package.<br>
</div></blockquote><br>I tend to opt for the most recent revision of a modern force field, something like Gromos96 53a6, unless there is a specific basis for a choice like 43a1 (i.e., consistency with previous work, parameters that have been specifically developed for use with that parameter set, etc.) Otherwise, blindly choosing a force field with no basis for the choice is not a good idea.<br>
<br>If you are considering making ad hoc changes to your .top to force a certain behavior, I think you are asking for trouble, especially if it comes to peer review. Force fields aren't designed to be tinkered with according to the whim of a user. If you are willing to expend a great deal of effort in developing new parameters to accurately replicate real-life behavior, then by all means do so. Otherwise, I would be very suspicious of results that come from untested modifications.<br>
<br>-Justin<br><br>
<blockquote class="gmail_quote" style="PADDING-LEFT: 1ex; MARGIN: 0px 0px 0px 0.8ex; BORDER-LEFT: #ccc 1px solid">Thanks<br>JJ
<div class="Ih2E3d"><br><br>On Wed, Jan 7, 2009 at 3:07 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> jayant james wrote:<br><br> Hi!<br> In the protein of my interest I find that oxygen atoms of Thr,<br> Tyr and Ser are close to the bound Ca2+. So is there any<br> information on charge modifications that need to be made for the<br>
specified amino acids in the gromacs force field, such that the<br> Ca2+ would stay and not stray!<br><br><br> I should hope you're not using the "Gromacs force field" (ffgmx), it<br> has long been deprecated. You are better off with a more modern<br>
Gromos96 parameter set.<br><br> I think the best procedure would first to refer to the citation you<br> were provided earlier today, then proceed with a further literature<br> search. Has anyone simulated a similar protein with any success?<br>
You may be in for some parameterization if you find the existing<br> parameters inadequate. Beware, such work is very difficult.<br><br> -Justin<br><br> Thanks<br> Jayant James<br><br><br><br> On Wed, Jan 7, 2009 at 10:36 AM, 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="Wj3C7c"> <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>
jayant james wrote:<br><br> Hi!<br> Thanks for your suggestions. In my protein there are two<br> phosphorylated serines (-vely charged) adjacent to each<br> other on<br>
chain A and these serines start to move towards the Ca2+<br> ion in<br> chain B. It is only after a while of them being close to each<br> other that the Ca2+ starts popping off. Also to modify the<br>
charges of the residues adjoining the Ca2+ ion, what file<br> should<br> I be looking for to make these changes? would it be the<br> topology<br> file?<br><br><br> Well, that behavior may not be entirely unexpected. Two adjacent<br>
phosphates would likely repel each other, making it difficult to<br> bind ions at that location. Is there some basis to suggest that<br> Ca2+ should even bind there?<br><br> Yes, charges are in your topology. If you generated your<br>
.top with<br> pdb2gmx, then obviously there are entries in the appropriate .rtp<br> file that can be modified to generate future topologies with the<br> same charges.<br><br><br> As per your suggestions I would be going with the PME so<br>
should<br> the mdp file modified to incorporate the parameters below be<br> sufficient?!!<br> <br> Well, assuming those are the appropriate cutoff's for your force<br>
field, and the rest of your .mdp file is satisfactory (and you<br> actually specify PME), then it is probably reasonable enough.<br><br> -Justin<br><br> ns_type = grid<br> rlist = 1.0<br>
rcoulomb = 1.0<br> rvdw = 1.4<br> vdwtype = cut-off<br> fourierspacing = 0.12<br> fourier_nx = 0<br> fourier_ny = 0<br>
fourier_nx = 0<br> pme_order = 4<br> ewald_rtol = 1e-5<br> optimize_fft = yes<br> disre = simple<br> disre_weighting = equal<br>
Thanks<br> Jayant James<br><br><br> ------------------------------------------------------------------------<br><br><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>
<mailto:<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><br>
<a href="http://www.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://www.gromacs.org/mailman/listinfo/gmx-users</a><br> Please search the archive at <a href="http://www.gromacs.org/search" target="_blank">http://www.gromacs.org/search</a><br>
before posting!<br> Please don't post (un)subscribe requests to the list. Use the<br> www 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<br> <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><br> Justin A. Lemkul<br> Graduate Research Assistant<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="Ih2E3d"><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>
<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 class="Ih2E3d"><br><br> <a href="http://www.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://www.gromacs.org/mailman/listinfo/gmx-users</a><br> Please search the archive at <a href="http://www.gromacs.org/search" target="_blank">http://www.gromacs.org/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>
</div>
<div>
<div></div>
<div class="Wj3C7c"> <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/mailing_lists/users.php" target="_blank">http://www.gromacs.org/mailing_lists/users.php</a><br>
<br><br><br><br> -- Jayasundar Jayant James<br><br> <a href="http://www.chick.com/reading/tracts/0096/0096_01.asp" target="_blank">www.chick.com/reading/tracts/0096/0096_01.asp</a><br> <<a href="http://www.chick.com/reading/tracts/0096/0096_01.asp" target="_blank">http://www.chick.com/reading/tracts/0096/0096_01.asp</a>><br>
<<a href="http://www.chick.com/reading/tracts/0096/0096_01.asp" target="_blank">http://www.chick.com/reading/tracts/0096/0096_01.asp</a>>)<br><br><br><br> ------------------------------------------------------------------------<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://www.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://www.gromacs.org/mailman/listinfo/gmx-users</a><br> Please search the archive at <a href="http://www.gromacs.org/search" target="_blank">http://www.gromacs.org/search</a><br>
before posting!<br> Please don't post (un)subscribe requests to the list. Use the<br> www 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/mailing_lists/users.php" target="_blank">http://www.gromacs.org/mailing_lists/users.php</a><br>
<br><br> -- ========================================<br><br> Justin A. Lemkul<br> Graduate Research Assistant<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> _______________________________________________<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://www.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://www.gromacs.org/mailman/listinfo/gmx-users</a><br>
Please search the archive at <a href="http://www.gromacs.org/search" target="_blank">http://www.gromacs.org/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>
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><br>-- <br>Jayasundar Jayant James<br><br><a href="http://www.chick.com/reading/tracts/0096/0096_01.asp" target="_blank">www.chick.com/reading/tracts/0096/0096_01.asp</a> <<a href="http://www.chick.com/reading/tracts/0096/0096_01.asp" target="_blank">http://www.chick.com/reading/tracts/0096/0096_01.asp</a>>)<br>
<br></div></div></blockquote><br>-- <br>
<div>
<div></div>
<div class="Wj3C7c">========================================<br><br>Justin A. Lemkul<br>Graduate Research Assistant<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://www.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://www.gromacs.org/mailman/listinfo/gmx-users</a><br>
Please search the archive at <a href="http://www.gromacs.org/search" target="_blank">http://www.gromacs.org/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/mailing_lists/users.php" target="_blank">http://www.gromacs.org/mailing_lists/users.php</a><br></div></div></blockquote></div><br><br clear="all"><br>-- <br>Jayasundar Jayant James<br>
<br><a href="http://www.chick.com/reading/tracts/0096/0096_01.asp">www.chick.com/reading/tracts/0096/0096_01.asp</a>) <br><br>