<br><br>
<div><span class="gmail_quote">2008/1/10, Mark Abraham &lt;<a href="mailto:Mark.Abraham@anu.edu.au">Mark.Abraham@anu.edu.au</a>&gt;:</span></div>
<blockquote class="gmail_quote" style="PADDING-LEFT: 1ex; MARGIN: 0px 0px 0px 0.8ex; BORDER-LEFT: #ccc 1px solid">Stanislav Bobritsky wrote:<br><br>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&gt; Well, but (as I know from Manual, Ch. 6, p. 120) there is
<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp; &quot;approach to do targeted MD&quot; in GMX.<br>&gt;<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp; Did you add position restraint to your protein<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp; topology?<br>&gt;<br>&gt;<br>&gt; No, I did not.<br>&gt; What group must be restrained therefore? Whole protein or whole protein
<br>&gt; except loop?<br><br>From &quot;man grompp&quot;,<br><br>&gt; When using position restraints a file with restraint coordinates can be<br>&gt; supplied with -r, otherwise restraining will be done with respect to the
<br>&gt; conformation from the -c option. For free energy calculation the the<br>&gt; coordinates for the B topology can be supplied with -rb, otherwise they will<br>&gt; be equal to those of the A topology.<br><br>Thus I expect your earlier grompp command line didn&#39;t do the job you
<br>thought it did.<br><br>You could restrain just the loop if you knew the protein was stable, but<br>easiest is to present conformation A as -c and B as -r. pdb2gmx will<br>have supplied you with a posre.itp file, and you need -DPOSRES in the
<br>include .mdp directive so that the former works.<br><br>Your .mdp file as supplied earlier should not T-couple only your solute,<br>and neither should it couple your counterions separately from the<br>solvent. See <a href="http://wiki.gromacs.org/index.php/Thermostats">
http://wiki.gromacs.org/index.php/Thermostats</a></blockquote>
<div>&nbsp;</div>
<div>&nbsp;</div>
<div>
<p><font style="BACKGROUND-COLOR: #ff6666">Some hints on practices that generally <b>not a good idea</b> to use: </font></p>
<ul>
<li><font style="BACKGROUND-COLOR: #ff6666">Do not use separate thermostats for different components of your system. Some molecular dynamics thermostats only work well in the thermodynamic limit. If you use one thermostat for, say, a small molecule, another for protein, and another for water, you are likely introducing errors and artifacts that are hard to predict. In particular, do not couple ions in aqueous solvent differently from that 
</font></li>
<li><font style="BACKGROUND-COLOR: #ff6666">solvent.</font> </li></ul></div>
<div><font style="BACKGROUND-COLOR: #33ffff">Sorry that I do not actually understand here. The link I copied above shows that better not to &quot;couple ions in aqueous solvent differently from that solvent&quot;. Maybe not separately but differently (mean different temperature)?
</font></div>
<div><font style="BACKGROUND-COLOR: #33ffff"></font>&nbsp;</div>
<div><font style="BACKGROUND-COLOR: #33ffff">I am actually a novice in Gromacs.</font></div>
<div>&nbsp;</div>
<div>&nbsp;</div>
<div>&nbsp;</div>
<div>&nbsp;</div>
<div><br>&nbsp;</div>
<blockquote class="gmail_quote" style="PADDING-LEFT: 1ex; MARGIN: 0px 0px 0px 0.8ex; BORDER-LEFT: #ccc 1px solid">Based on the foregoing, I&#39;d guess that you&#39;re fairly new to MD and/or<br>GROMACS. I&#39;d strongly recommend doing some (more) tutorials and
<br>background reading so that you have more idea what &quot;normal&quot; MD is, and<br>the consequences of changes you might make. Otherwise the odds are<br>excellent you&#39;re just using an complicated quasi-random number generator.
<br><br>Mark<br>_______________________________________________<br>gmx-users mailing list&nbsp;&nbsp;&nbsp;&nbsp;<a href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a><br><a href="http://www.gromacs.org/mailman/listinfo/gmx-users">
http://www.gromacs.org/mailman/listinfo/gmx-users</a><br>Please search the archive at <a href="http://www.gromacs.org/search">http://www.gromacs.org/search</a> before posting!<br>Please don&#39;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">gmx-users-request@gromacs.org</a>.<br>Can&#39;t post? Read <a href="http://www.gromacs.org/mailing_lists/users.php">http://www.gromacs.org/mailing_lists/users.php
</a><br></blockquote><br>