<!DOCTYPE html PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
<html>
<head>
<meta content="text/html;charset=ISO-8859-1" http-equiv="Content-Type">
</head>
<body bgcolor="#ffffff" text="#000000">
<div align="justify">Dear Sonali,<br>
<br>
Unfortunately development of force-field parameters is difficult, and
even more so for bivalent ions and transition metals. This is certainly
not suitable as a project to begin with. <br>
<br>
Having said that, a careful search in the literature will reveal that
many sets of parameters for metal ions do exist. Try to find a set that
works best with the water model that you're using or your particular
model, and carry out some simulations to validate that the results are
correct, reproducible and robust. The following references discuss some
aspects of simulations of ions and proteins:<br>
</div>
<br>
@article{Project2006,<br>
Author = {Project, E. and Friedman, R. and Nachliel, E. and Gutman,
M.},<br>
Title = {A Molecular Dynamics Study of the Effect of Ca{2+} Removal
on Calmodulin Structure},<br>
Journal = {Biophys. J.},<br>
Volume= {90},<br>
Pages = {3842-3850},<br>
Year = {2006} <br>
}<br>
<br>
<br>
@Article{Fyta2010,<br>
author = "Fyta, M and Kalcher, I and Dzubiella, J and Vrbka, L and
Netz, R R",<br>
title = {Ionic force field optimization based on single-ion and
ion-pair solvation properties},<br>
journal = "J Chem Phys",<br>
year = "2010",<br>
volume = "132",<br>
pages = "024911-024911"<br>
}<br>
<br>
Ran.<br>
<br>
sonali dhindwal wrote:
<blockquote cite="mid787770.25100.qm@web94613.mail.in2.yahoo.com"
type="cite">
<table border="0" cellpadding="0" cellspacing="0">
<tbody>
<tr>
<td
style="font-family: inherit; font-style: inherit; font-variant: inherit; font-weight: inherit; font-size: inherit; line-height: inherit; font-size-adjust: inherit; font-stretch: inherit;"
valign="top">Hello All ,<br>
I am also trying to simulate my protein with Mn ion present in it.<br>
So can I create the topology entry for Mn ion similar to MG2+ ion ?
and how can I get the values of C6 and C12 leonard jones potential in
[atom type] entry and will it be required to add [nonbond_params] also
? how can I get them ?<br>
<br>
Also I want to use Fe2+ also, but it is also not included in gromos
force field. This problem of adding other ions and molecule in the
system always remains. And being new to this field, can someone suggest
in simple terms how to include these considering the person is not an
expert in this field.<br>
Please help.<br>
<br>
<div>--<br>
Sonali Dhindwal</div>
<br>
<br>
--- On <b>Tue, 8/6/10, Vitaly Chaban <i><a class="moz-txt-link-rfc2396E" href="mailto:vvchaban@gmail.com"><vvchaban@gmail.com></a></i></b>
wrote:<br>
<blockquote
style="border-left: 2px solid rgb(16, 16, 255); margin-left: 5px; padding-left: 5px;"><br>
From: Vitaly Chaban <a class="moz-txt-link-rfc2396E" href="mailto:vvchaban@gmail.com"><vvchaban@gmail.com></a><br>
Subject: [gmx-users] Re: Simulation with CsCl<br>
To: <a class="moz-txt-link-abbreviated" href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a><br>
Date: Tuesday, 8 June, 2010, 9:03 PM<br>
<br>
<div class="plainMail">> Hi all:<br>
> I am trying to simulate a polysaccharide in solution of water and
CsCl, but<br>
> cesium is not parametrized in the gromacs 4. I am using the force<br>
> field GROMOS 96.<br>
> I have looked for the parameters of Cs+ in the OPLS Force field
and I have<br>
> created 4 files: Cs.atp, Cs.itp, Csnb.itp and Cs.rtp. Also I have
include<br>
> the parameters for Cs+ in the "ions.itp" file, but it didn't work,
Can<br>
> someone help me?<br>
> Thanks in advance.<br>
> Cecilia.<br>
<br>
Cecilia -<br>
<br>
You do not need so many files so cesium. Just copy your parameters to<br>
the force field used and create the topology entry for cesium similar<br>
to<br>
<br>
[ moleculetype ]<br>
Ar 1<br>
<br>
[ atoms ]<br>
1 Ar 1 Ar Ar 0 0.0<br>
<br>
Do not forget, cesium+ is an ion...<br>
<br>
-- <br>
Dr. Vitaly Chaban<br>
--<br>
gmx-users mailing list <a ymailto="mailto:gmx-users@gromacs.org"
href="/mc/compose?to=gmx-users@gromacs.org">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/search"
target="_blank">http://www.gromacs.org/search</a> before posting!<br>
Please don't post (un)subscribe requests to the list. Use the<br>
www interface or send it to <a
ymailto="mailto:gmx-users-request@gromacs.org"
href="/mc/compose?to=gmx-users-request@gromacs.org">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>
</blockquote>
</td>
</tr>
</tbody>
</table>
<br>
</blockquote>
<br>
<br>
<pre class="moz-signature" cols="72">--
------------------------------------------------------
Ran Friedman
Postdoctoral Fellow
Computational Structural Biology Group (A. Caflisch)
Department of Biochemistry
University of Zurich
Winterthurerstrasse 190
CH-8057 Zurich, Switzerland
Tel. +41-44-6355559
Email: <a class="moz-txt-link-abbreviated" href="mailto:r.friedman@bioc.uzh.ch">r.friedman@bioc.uzh.ch</a>
Skype: ran.friedman
------------------------------------------------------
</pre>
</body>
</html>