Thanks Justin for the help,<br />
<br />
I have Manganese in my protein with other two ligands at the active site. <br />
I am using GROMOS96 43a1 force field.<br />
<br />
Mangnese may have topology similar to Magnesium 2+, What I found is that I need to do changes in ffG43a1.rtp and ffG43a1nb.itp. But how should I calculate the nonbond_params for ffG43a1nb.itp.<br />
<br />
For both the other two ligands which contains phosphate and enol groups, I was able to create topology using PRODRG server.<br />
<br />
Regards<br />
<br />
onetwo wrote:<br />
> Hello All,<br />
><br />
> I am new in this field, and I want to do simulation study on one of the protein conatining a metal ion and study its ability to form co-oridnation bond with other ligand, which is not defined in GROMOS force field which I have tried. If this choice of force field is correct for such type of study ?<br />
><br />
> I have read Gromacs Manual Chapter 5, also I am following this gromacs mailing list for quite some time to get a help on how to include a new metal ion or a new ligand in their simulation and they have been refered to refer to the http://www.gromacs.org/Documentation/How-tos/Parameterization,<br />
><br />
> but in this its not mentioned that how to actually do parameterization<br />
><br />
> neither in manual it has been told on how to generate it ( due apologies,,I know, I may be wrong )<br />
><br />
<br />
The manual will not cover every possible topic.  Parameterization is described in the primary literature for the force field you wish you use.  "GROMOS" is not very specific - there are numerous parameter sets within this class of force field.<br />
<br />
> and it is more difficult for a person like me who doesnt have much knowledge in this field, so if someone can please help me guiding how to do parameterisation, I know this is not a trivial task, but still help in this regard may help many other fellow users besides me.<br />
<br />
Aside from telling you to read all the literature you can, there's not much more help anyone can do to help you.  You haven't described very well what your system is.  What kind of functional groups are in your ligand?  If the functional groups are shared with existing parameters, then you'll have a rather easy time constructing a topology.  What type of metal ion is it?  If you're dealing with a transition metal, using a standard MM force field may not work very well.  Note the "Exotic Species" heading on the page you quoted above.<br />
<br />
-Justin<br />
<br />
> Thanks in advance<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]vt.edu | (540) 231-9080<br />
http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin<br />
<br />
========================================<br />
-- gmx-users mailing list    gmx-users@gromacs.org<br />
http://lists.gromacs.org/mailman/listinfo/gmx-users<br />
Please search the archive at http://www.gromacs.org/search before posting!<br />
Please don't post (un)subscribe requests to the list. Use the www interface or send it to gmx-users-request@gromacs.org.<br />
Can't post? Read http://www.gromacs.org/mailing_lists/users.php<br />
<br><Table border=0 Width=644 Height=57 cellspacing=0 cellpadding=0 style="font-family:Verdana;font-size:11px;line-height:15px;"><TR><td><A HREF="http://sigads.rediff.com/RealMedia/ads/click_nx.ads/www.rediffmail.com/signatureline.htm@Middle?" target="_blank"><IMG SRC="http://sigads.rediff.com/RealMedia/ads/adstream_nx.ads/www.rediffmail.com/signatureline.htm@Middle"></A></td></TR></Table>