Please ignore my previous mail. There was typo error.<br><br>Dear colleagues,<br>
<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; I want to do MD simulation of only the extracellular
domain of a pentamer, excluding the transmembrane domain. So, I want to
restrain the part of the extracellular domain in contact with the
transmembrane domain. I am using "position restraints" for few of those
residues as well as to the two ends of one of the subunits, in which
there are missing residues.<br>
<br>
[ position_restraints ]<br>
; atom&nbsp; type&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; fx&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; fy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; fz<br>
&nbsp;&nbsp; 457&nbsp;&nbsp;&nbsp;&nbsp; 1&nbsp; 1000&nbsp; 1000&nbsp; 1000<br>
&nbsp;&nbsp; 459&nbsp;&nbsp;&nbsp;&nbsp; 1&nbsp; 1000&nbsp; 1000&nbsp; 1000<br>
&nbsp;&nbsp; 460&nbsp;&nbsp;&nbsp;&nbsp; 1&nbsp; 1000&nbsp; 1000&nbsp; 1000<br>
&nbsp;&nbsp; 461&nbsp;&nbsp;&nbsp;&nbsp; 1&nbsp; 1000&nbsp; 1000&nbsp; 1000<br>
<br>
<br>
&nbsp;&nbsp;&nbsp;&nbsp; After few ns of MD, I observed that 3 of the subunits are distorted while the other 2 aren't. <br>
<br>
1) Is this because of imposing position restraints or is it because of the gap in between (due to missing residues) in one of the subunits? <br>
<br>
2)Is there any other way to restrain (apart from position restraints) the protein for this purpose? <br>
<br>
&nbsp;&nbsp;&nbsp; Coming to the missing residues part, I decided to use distance
restraints between the two ends. Please note that only one of the chain
is in covalently linked with the transmembrane domain. Rest of them
have vanderwaals interactions.<br>
<br>
Thanks for any kind of suggestion in this regard.<br>
<br>
Prema.<br><br>----- Original Message -----<br>From: &lt;plmallip@mail.uh.edu&gt;<br>Date: Friday, August 15, 2008 2:00 pm<br>Subject: Position restraints on a pentamer<br>To: gmx-users@gromacs.org<br><br>&gt; Dear colleagues,<br>&gt; <br>&gt; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; I want to do MD simulation of only the <br>&gt; extracellular domain of a pentamer, excluding the transmembrane <br>&gt; domain. So, I want to restrain the part of the extracellular <br>&gt; domain in contact with the transmembrane domain. I am using <br>&gt; "position restraints" for few of those residues as well as to <br>&gt; the two ends of one of the subunits, in which there are missing <br>&gt; residues.<br>&gt; [ position_restraints ]<br>&gt; ; atom&nbsp; type&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; fx&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; fy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; fz<br>&gt; &nbsp;&nbsp; 457&nbsp;&nbsp;&nbsp;&nbsp; 1&nbsp; 1000&nbsp; 1000&nbsp; 1000<br>&gt; &nbsp;&nbsp; 459&nbsp;&nbsp;&nbsp;&nbsp; 1&nbsp; 1000&nbsp; 1000&nbsp; 1000<br>&gt; &nbsp;&nbsp; 460&nbsp;&nbsp;&nbsp;&nbsp; 1&nbsp; 1000&nbsp; 1000&nbsp; 1000<br>&gt; &nbsp;&nbsp; 461&nbsp;&nbsp;&nbsp;&nbsp; 1&nbsp; 1000&nbsp; 1000&nbsp; 1000<br>&gt; <br>&gt; <br>&gt; &nbsp;&nbsp;&nbsp;&nbsp; After few ns of MD, I observed that 3 of the subunits are <br>&gt; distorted while the other 2 aren't. <br>&gt; <br>&gt; 1) Is this because of imposing position restraints or is it <br>&gt; because of the gap between two subunits in on the subunits? <br>&gt; <br>&gt; 2)Is there any other way to restrain (apart from position <br>&gt; restraints) the protein for this purpose? <br>&gt; <br>&gt; &nbsp;&nbsp;&nbsp; Coming to the missing residues part, I decided to use <br>&gt; distance restraints between the two ends. Please note that only <br>&gt; one of the chain is in covalently linked with the transmembrane <br>&gt; domain. Rest of them have vanderwaals interactions.<br>&gt; <br>&gt; Thanks for any kind of suggestion in this regard.<br>&gt; <br>&gt; Prema.<br>&gt; <br>&gt;