<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
<html>
  <head>
    <meta content="text/html; charset=ISO-8859-1"
      http-equiv="Content-Type">
    <title></title>
  </head>
  <body text="#000000" bgcolor="#ffffff">
    On 11/11/2010 7:46 PM, NG HUI WEN wrote:
    <blockquote
cite="mid:EE6BA5FF4C277348A3A297EDC1D54DC5C330BE@wsexbe2.nottingham.edu.my"
      type="cite">
      <meta http-equiv="Content-Type" content="text/html;
        charset=ISO-8859-1">
      <meta name="Generator" content="Microsoft Word 12 (filtered
        medium)">
      <style>
<!--
 /* Font Definitions */
 @font-face
        {font-family:Calibri;
        panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
        {font-family:Tahoma;
        panose-1:2 11 6 4 3 5 4 4 2 4;}
@font-face
        {font-family:Consolas;
        panose-1:2 11 6 9 2 2 4 3 2 4;}
@font-face
        {font-family:"Courier New \;";
        panose-1:0 0 0 0 0 0 0 0 0 0;}
 /* Style Definitions */
 p.MsoNormal, li.MsoNormal, div.MsoNormal
        {margin:0in;
        margin-bottom:.0001pt;
        font-size:11.0pt;
        font-family:"Calibri","sans-serif";
        color:black;}
a:link, span.MsoHyperlink
        {mso-style-priority:99;
        color:blue;
        text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
        {mso-style-priority:99;
        color:purple;
        text-decoration:underline;}
p.MsoPlainText, li.MsoPlainText, div.MsoPlainText
        {mso-style-priority:99;
        mso-style-link:"Plain Text Char";
        margin:0in;
        margin-bottom:.0001pt;
        font-size:10.5pt;
        font-family:Consolas;
        color:black;}
p
        {mso-style-priority:99;
        mso-margin-top-alt:auto;
        margin-right:0in;
        mso-margin-bottom-alt:auto;
        margin-left:0in;
        font-size:12.0pt;
        font-family:"Times New Roman","serif";
        color:black;}
p.MsoListParagraph, li.MsoListParagraph, div.MsoListParagraph
        {mso-style-priority:34;
        margin-top:0in;
        margin-right:0in;
        margin-bottom:0in;
        margin-left:.5in;
        margin-bottom:.0001pt;
        font-size:11.0pt;
        font-family:"Calibri","sans-serif";
        color:black;}
span.PlainTextChar
        {mso-style-name:"Plain Text Char";
        mso-style-priority:99;
        mso-style-link:"Plain Text";
        font-family:Consolas;}
span.EmailStyle19
        {mso-style-type:personal;
        font-family:"Calibri","sans-serif";
        color:windowtext;}
span.EmailStyle21
        {mso-style-type:personal-reply;
        font-family:"Calibri","sans-serif";
        color:#1F497D;}
.MsoChpDefault
        {mso-style-type:export-only;
        font-size:10.0pt;}
@page WordSection1
        {size:8.5in 11.0in;
        margin:1.0in 1.0in 1.0in 1.0in;}
div.WordSection1
        {page:WordSection1;}
 /* List Definitions */
 @list l0
        {mso-list-id:1137138639;
        mso-list-type:hybrid;
        mso-list-template-ids:1018980054 67698705 67698713 67698715 67698703 67698713 67698715 67698703 67698713 67698715;}
@list l0:level1
        {mso-level-text:"%1\)";
        mso-level-tab-stop:none;
        mso-level-number-position:left;
        text-indent:-.25in;}
@list l1
        {mso-list-id:1605770817;
        mso-list-type:hybrid;
        mso-list-template-ids:-1509657810 67698705 67698713 67698715 67698703 67698713 67698715 67698703 67698713 67698715;}
@list l1:level1
        {mso-level-text:"%1\)";
        mso-level-tab-stop:none;
        mso-level-number-position:left;
        text-indent:-.25in;}
ol
        {margin-bottom:0in;}
ul
        {margin-bottom:0in;}
-->
</style><!--[if gte mso 9]><xml>
 <o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
 <o:shapelayout v:ext="edit">
  <o:idmap v:ext="edit" data="1" />
 </o:shapelayout></xml><![endif]-->
      <div class="WordSection1">
        <p class="MsoNormal"><span style="color: rgb(31, 73, 125);">Thanks
            a lot Justin and Mark for
            your useful input.<o:p></o:p></span></p>
        <p class="MsoNormal"><span style="color: rgb(31, 73, 125);"><o:p>&nbsp;</o:p></span></p>
        <p class="MsoNormal"><span style="color: rgb(31, 73, 125);">Indeed
            Justin was right, the
            quest to dissect the total energy of the system to get that
            contributed by the
            protein alone was not trivial at all! </span></p>
      </div>
    </blockquote>
    <br>
    Indeed. If I'd observed you were doing PME, I'd have made the same
    point.<span style="color: rgb(31, 73, 125);"><o:p> <br>
      </o:p></span><br>
    <blockquote
cite="mid:EE6BA5FF4C277348A3A297EDC1D54DC5C330BE@wsexbe2.nottingham.edu.my"
      type="cite">
      <div class="WordSection1">
        <p class="MsoNormal"><span style="color: rgb(31, 73, 125);">I
            missed out this thread
            yesterday <a moz-do-not-send="true"
              href="http://www.mail-archive.com/gmx-users@gromacs.org/msg34610.html">http://www.mail-archive.com/gmx-users@gromacs.org/msg34610.html</a>
            as well as Mark&#8217;s trick to decompose bonded terms by hacking
            the .top
            file. (However, I read from <a moz-do-not-send="true"
              href="http://www.gromacs.org/Documentation/Gromacs_Utilities/g_energy">http://www.gromacs.org/Documentation/Gromacs_Utilities/g_energy</a>
            that the Coul-recip term is not decomposable regardless of
            the tricks used&#8230;)<o:p></o:p></span></p>
        <p class="MsoNormal"><span style="color: rgb(31, 73, 125);"><o:p>&nbsp;</o:p></span></p>
        <p class="MsoNormal"><span style="color: rgb(31, 73, 125);">Mark,
            following your advice, I
            added &#8211;nsteps 0 (-1 didn&#8217;t work) to tpbconv. I managed to
            get the interactive
            prompt succesfully this time. </span></p>
      </div>
    </blockquote>
    <br>
    Yes, nsteps=-1 is a GROMACS 4.5-ism, sorry.<span style="font-family:
      &quot;Courier New&quot;;"><o:p> <br>
      </o:p></span><br>
    <blockquote
cite="mid:EE6BA5FF4C277348A3A297EDC1D54DC5C330BE@wsexbe2.nottingham.edu.my"
      type="cite">
      <div class="WordSection1">
        <p class="MsoPlainText"><span style="font-family: &quot;Courier
            New&quot;;">Reading toplogy
            and shit from rerun1.tpr<br>
            Reading file rerun1.tpr, VERSION 4.0.7 (single precision)<br>
            Setting nsteps to 0<br>
            Group&nbsp;&nbsp;&nbsp;&nbsp; 0 (&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; System) has
            100581 elements<br>
            Group&nbsp;&nbsp;&nbsp;&nbsp; 1 (&nbsp;&nbsp;&nbsp;&nbsp; Protein) has&nbsp;
            4002 elements<br>
            Group&nbsp;&nbsp;&nbsp;&nbsp; 2 (&nbsp;&nbsp; Protein-H) has&nbsp; 3145
            elements<br>
            Group&nbsp;&nbsp;&nbsp;&nbsp; 3 (&nbsp;&nbsp;&nbsp;&nbsp; C-alpha)
            has&nbsp;&nbsp; 412 elements<br>
            Group&nbsp;&nbsp;&nbsp;&nbsp; 4 (&nbsp;&nbsp;&nbsp; Backbone) has&nbsp; 1236
            elements<br>
            Group&nbsp;&nbsp;&nbsp;&nbsp; 5 (&nbsp;&nbsp; MainChain) has&nbsp; 1649
            elements<br>
            Group&nbsp;&nbsp;&nbsp;&nbsp; 6 (MainChain+Cb) has&nbsp; 2027 elements<br>
            Group&nbsp;&nbsp;&nbsp;&nbsp; 7 ( MainChain+H) has&nbsp; 2039 elements<br>
            Group&nbsp;&nbsp;&nbsp;&nbsp; 8 (&nbsp;&nbsp; SideChain) has&nbsp; 1963
            elements<br>
            Group&nbsp;&nbsp;&nbsp;&nbsp; 9 ( SideChain-H) has&nbsp; 1496 elements<br>
            Group&nbsp;&nbsp;&nbsp; 10 ( Prot-Masses) has&nbsp; 4002 elements<br>
            Group&nbsp;&nbsp;&nbsp; 11 ( Non-Protein) has 96579 elements<br>
            Group&nbsp;&nbsp;&nbsp; 12 (&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; POPE)
            has 13052 elements<br>
            Group&nbsp;&nbsp;&nbsp; 13 (&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
            SOL) has 83523 elements<br>
            Group&nbsp;&nbsp;&nbsp; 14 (&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
            CL-) has&nbsp;&nbsp;&nbsp;&nbsp; 4 elements<br>
            Group&nbsp;&nbsp;&nbsp; 15 (&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Other) has
            96579 elements<br>
            Group&nbsp;&nbsp;&nbsp; 16 (&nbsp;&nbsp;&nbsp;&nbsp; SOL_CL-) has 83527
            elements<br>
            Group&nbsp;&nbsp;&nbsp; 17 (Protein_POPE) has 17054 elements<br>
            Select a group: </span><o:p></o:p></p>
        <p class="MsoNormal"><span style="color: rgb(31, 73, 125);"><o:p>&nbsp;</o:p></span></p>
        <p class="MsoNormal"><span style="color: rgb(31, 73, 125);">As
            suggested, I have also tried
            to compare the (average) values of subset.tpr and
            fullsystem.tpr, these are the
            total energies of:<o:p></o:p></span></p>
        <p class="MsoListParagraph" style="text-indent: -0.25in;"><!--[if !supportLists]--><span
            style="color: rgb(31, 73, 125);"><span style="">1)<span
                style="font: 7pt &quot;Times New Roman&quot;;">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
              </span></span></span><!--[endif]--><span style="color:
            rgb(31, 73, 125);">System=&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; -1.28209 e+06<o:p></o:p></span></p>
        <p class="MsoListParagraph" style="text-indent: -0.25in;"><!--[if !supportLists]--><span
            style="color: rgb(31, 73, 125);"><span style="">2)<span
                style="font: 7pt &quot;Times New Roman&quot;;">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
              </span></span></span><!--[endif]--><span style="color:
            rgb(31, 73, 125);">Protein= &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; -9571.86<o:p></o:p></span></p>
        <p class="MsoListParagraph" style="text-indent: -0.25in;"><!--[if !supportLists]--><span
            style="color: rgb(31, 73, 125);"><span style="">3)<span
                style="font: 7pt &quot;Times New Roman&quot;;">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
              </span></span></span><!--[endif]--><span style="color:
            rgb(31, 73, 125);">Lipid= &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; -296121<o:p></o:p></span></p>
        <p class="MsoListParagraph" style="text-indent: -0.25in;"><!--[if !supportLists]--><span
            style="color: rgb(31, 73, 125);"><span style="">4)<span
                style="font: 7pt &quot;Times New Roman&quot;;">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
              </span></span></span><!--[endif]--><span style="color:
            rgb(31, 73, 125);">SOL_CL= &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; -821353<o:p></o:p></span></p>
        <p class="MsoListParagraph" style="text-indent: -0.25in;"><!--[if !supportLists]--><span
            style="color: rgb(31, 73, 125);"><span style="">5)<span
                style="font: 7pt &quot;Times New Roman&quot;;">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
              </span></span></span><!--[endif]--><span style="color:
            rgb(31, 73, 125);">Other= &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; -1.22684e+06<o:p></o:p></span></p>
        <p class="MsoNormal"><span style="color: rgb(31, 73, 125);">It
            was found that the: <o:p></o:p></span></p>
        <p class="MsoListParagraph" style="text-indent: -0.25in;"><!--[if !supportLists]--><span
            style="color: rgb(31, 73, 125);"><span style="">1)<span
                style="font: 7pt &quot;Times New Roman&quot;;">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
              </span></span></span><!--[endif]--><span style="color:
            rgb(31, 73, 125);">Total energies (Protein
            + lipid + SOL_CL-) not equal to total energy of the system<o:p></o:p></span></p>
        <p class="MsoListParagraph" style="text-indent: -0.25in;"><!--[if !supportLists]--><span
            style="color: rgb(31, 73, 125);"><span style="">2)<span
                style="font: 7pt &quot;Times New Roman&quot;;">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
              </span></span></span><!--[endif]--><span style="color:
            rgb(31, 73, 125);">Total energies of
            (Lipid + SOL_CL-) not equal to total energy of other<o:p></o:p></span></p>
        <p class="MsoListParagraph" style="text-indent: -0.25in;"><!--[if !supportLists]--><span
            style="color: rgb(31, 73, 125);"><span style="">3)<span
                style="font: 7pt &quot;Times New Roman&quot;;">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
              </span></span></span><!--[endif]--><span style="color:
            rgb(31, 73, 125);">Total energies of
            (others + protein) not equal to total energy of the system<o:p></o:p></span></p>
        <p class="MsoNormal"><span style="color: rgb(31, 73, 125);"><o:p>&nbsp;</o:p></span></p>
        <p class="MsoNormal"><span style="color: rgb(31, 73, 125);">I
            have yet to work out the
            reasons for the discrepancies, will spend some time
            pondering and trying to make
            sense of these (and other) values.</span></p>
      </div>
    </blockquote>
    <br>
    Aren't you omitting the inter-group non-bonded terms? Anyway, I was
    suggesting this comparison for seeing whether the bonded components
    can be decomposed group-wise. It is already known that the
    non-bonded decomposition works.<br>
    <br>
    Mark<br>
    <br>
    <blockquote
cite="mid:EE6BA5FF4C277348A3A297EDC1D54DC5C330BE@wsexbe2.nottingham.edu.my"
      type="cite">
      <div class="WordSection1">
        <p class="MsoNormal"><span style="color: rgb(31, 73, 125);"><o:p></o:p></span></p>
        <p class="MsoNormal"><span style="color: rgb(31, 73, 125);"><o:p>&nbsp;</o:p></span></p>
        <p class="MsoNormal"><span style="color: rgb(31, 73, 125);">Thanks
            a bunch again!<o:p></o:p></span></p>
        <p class="MsoNormal"><span style="color: rgb(31, 73, 125);"><o:p>&nbsp;</o:p></span></p>
        <p class="MsoNormal"><span style="color: rgb(31, 73, 125);">HW<o:p></o:p></span></p>
        <p class="MsoNormal"><span style="color: rgb(31, 73, 125);"><o:p>&nbsp;</o:p></span></p>
        <p class="MsoNormal"><span style="color: rgb(31, 73, 125);"><o:p>&nbsp;</o:p></span></p>
        <div>
          <div style="border-right: medium none; border-width: 1pt
            medium medium; border-style: solid none none; border-color:
            rgb(181, 196, 223) -moz-use-text-color -moz-use-text-color;
            padding: 3pt 0in 0in;">
            <p class="MsoNormal"><b><span style="font-size: 10pt;
                  font-family:
                  &quot;Tahoma&quot;,&quot;sans-serif&quot;; color:
                  windowtext;">From:</span></b><span style="font-size:
                10pt; font-family:
                &quot;Tahoma&quot;,&quot;sans-serif&quot;; color:
                windowtext;"> <a class="moz-txt-link-abbreviated" href="mailto:gmx-users-bounces@gromacs.org">gmx-users-bounces@gromacs.org</a>
                [<a class="moz-txt-link-freetext" href="mailto:gmx-users-bounces@gromacs.org">mailto:gmx-users-bounces@gromacs.org</a>] <b>On Behalf Of
                </b>Mark Abraham<br>
                <b>Sent:</b> Wednesday, November 10, 2010 9:08 PM<br>
                <b>To:</b> Discussion list for GROMACS users<br>
                <b>Subject:</b> Re: [gmx-users] mdrun -rerun: bonded
                interactions of protein<o:p></o:p></span></p>
          </div>
        </div>
        <p class="MsoNormal"><o:p>&nbsp;</o:p></p>
        <p class="MsoNormal">On 10/11/2010 11:29 PM, NG HUI WEN wrote: <o:p></o:p></p>
        <p class="MsoNormal">Hi Gmxusers,<o:p></o:p></p>
        <p class="MsoNormal">&nbsp;<o:p></o:p></p>
        <p class="MsoNormal">I have been trying to run mdrun &#8211;rerun to
          get the
          energy of the protein in my protein-lipid system. I know
          similar questions have
          been raised on this topic before, I have tried to glean useful
          information from
          them to solve my problem but unfortunately to no avail. Thanks
          for your
          patience! &nbsp;<o:p></o:p></p>
        <p class="MsoNormal">&nbsp;<o:p></o:p></p>
        <p class="MsoNormal">As I did not have energygrps in the initial
          .mdp file, I
          included it this time &nbsp;<o:p></o:p></p>
        <p class="MsoNormal">&nbsp;<o:p></o:p></p>
        <p class="MsoPlainText"><span style="font-size: 10pt;
            font-family: &quot;Courier New&quot;;">integrator&nbsp;
            = md&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <br>
            nsteps&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0&nbsp;&nbsp; <br>
            dt&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
            0.002&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <br>
            nstxout&nbsp;&nbsp;&nbsp;&nbsp; =
            50000&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <br>
            nstvout&nbsp;&nbsp;&nbsp;&nbsp; = 50000&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
            <br>
            nstenergy&nbsp;&nbsp; = 500&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <br>
            nstlog&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 500&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <br>
            Continuation&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
            yes&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <br>
            constraint_algorithm = lincs&nbsp; <br>
            constraints = all-bonds <br>
            lincs_iter&nbsp; = 1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <br>
            lincs_order = 4&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <br>
            ns_type&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
            grid&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <br>
            nstlist&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
            5&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; </span><o:p></o:p></p>
        <p class="MsoPlainText"><span style="font-size: 10pt;
            font-family: &quot;Courier New&quot;;">rlist&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
            = 1.2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <br>
            rcoulomb&nbsp;&nbsp;&nbsp; = 1.2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <br>
            rvdw&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
            1.2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <br>
            coulombtype = PME&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <br>
            pme_order&nbsp;&nbsp; = 4&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <br>
            fourierspacing&nbsp;&nbsp;&nbsp; =
            0.16&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <br>
            tcoupl&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
            Nose-Hoover&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <br>
            tc-grps&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = Protein
            POPE&nbsp;&nbsp;&nbsp; SOL_CL-&nbsp;&nbsp;&nbsp;&nbsp; </span><o:p></o:p></p>
        <p class="MsoPlainText"><span style="font-size: 10pt;
            font-family: &quot;Courier New&quot;;">tau_t&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
            = 0.1 0.1&nbsp;&nbsp; 0.1&nbsp;&nbsp; <br>
            ref_t&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 323 &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
            323&nbsp;&nbsp; 323&nbsp;&nbsp; <br>
            pcoupl&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
            Parrinello-Rahman&nbsp;&nbsp;&nbsp;&nbsp; <br>
            pcoupltype&nbsp; =
            semiisotropic&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; </span><o:p></o:p></p>
        <p class="MsoPlainText"><span style="font-size: 10pt;
            font-family: &quot;Courier New&quot;;">tau_p&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
            = 5.0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <br>
            ref_p&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 1.0
            1.0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <br>
            compressibility = 4.5e-5&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
            4.5e-5&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <br>
            pbc&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
            xyz&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <br>
            DispCorr&nbsp;&nbsp;&nbsp; = EnerPres&nbsp; <br>
            gen_vel&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
            no&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <br>
            nstcomm&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 1<br>
            comm-mode&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = Linear<br>
            comm-grps&nbsp; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;= Protein_POPE SOL_CL-<br>
            <b>energygrps&nbsp; = Protein SOL POPE</b></span><o:p></o:p></p>
        <p class="MsoPlainText"><span style="font-size: 10pt;
            font-family: &quot;Courier New&quot;;">&nbsp;</span><o:p></o:p></p>
        <p class="MsoPlainText"><span style="font-size: 10pt;
            font-family: &quot;Courier New&quot;;">I
            then did </span><o:p></o:p></p>
        <p class="MsoPlainText"><span style="font-size: 10pt;
            font-family: &quot;Courier New&quot;;">1)grompp
            &#8211;f new.mdp &#8211;n index.ndx &#8211;c old.tpr &#8211;o rerun.tpr
            &#8211;p topol.top</span><o:p></o:p></p>
        <p class="MsoPlainText"><span style="font-size: 10pt;
            font-family: &quot;Courier New&quot;;">2)trjconv
            &#8211;f old.trr &#8211;n index.ndx &#8211;s rerun.tpr &#8211;o rerun.trr (when
            prompted, I selected &#8220;0&#8221; system)</span><o:p></o:p></p>
        <p class="MsoPlainText"><span style="font-size: 10pt;
            font-family: &quot;Courier New&quot;;">3)mdrun
            &#8211;s rerun.tpr &#8211;rerun rerun.trr</span><o:p></o:p></p>
        <p class="MsoPlainText"><span style="font-size: 10pt;
            font-family: &quot;Courier New&quot;;">&nbsp;</span><o:p></o:p></p>
        <p class="MsoPlainText"><span style="font-size: 10pt;
            font-family: &quot;Courier New&quot;;">I
            notice that the previous post <a moz-do-not-send="true"
href="http://oldwww.gromacs.org/pipermail/gmx-users/2009-January/038968.html">http://oldwww.gromacs.org/pipermail/gmx-users/2009-January/038968.html</a>
            suggested to use tpbconv (on rerun.tpr) and trjconv (on
            rerun.trr) to extract
            the protein only. While it was possible to do so with
            trjconv, it wasn&#8217;t
            feasible with tpbconv (I&#8217;m using gromacs 4.0.7) &#8211; I might
            have
            missed out something as I did not get any prompt/output (see
            below)</span><o:p></o:p></p>
        <p class="MsoPlainText"><span style="font-size: 10pt;
            font-family: &quot;Courier New&quot;;">&nbsp;</span><o:p></o:p></p>
        <p class="MsoPlainText"><span style="font-size: 9pt;
            font-family: &quot;Courier New&quot;;">tpbconv
            -s topol.tpr -n index_P.ndx -o rerun2.tpr</span><o:p></o:p></p>
        <p class="MsoPlainText"><span style="font-size: 10pt;
            font-family: &quot;Courier New&quot;;">&nbsp;</span><o:p></o:p></p>
        <p class="MsoPlainText"><span style="font-size: 9pt;
            font-family: &quot;Courier New&quot;;">Reading
            toplogy and shit from topol.tpr<br>
            Reading file topol.tpr, VERSION 4.0.7 (single precision)<br>
            0 steps (0 ps) remaining from first run.<br>
            You've simulated long enough. Not writing tpr file</span><o:p></o:p></p>
        <p class="MsoNormal"><span style="font-size: 12pt; font-family:
            &quot;Times New Roman&quot;,&quot;serif&quot;;"><br>
            Looking at the code, tpbconv checks for whether there are
            any more steps to
            simulate before even considering letting you use it in the
            "create subset"
            mode. You could argue that this is buggy, because the way it
            computes whether
            there are any more steps will always indicate no more steps
            in such cases.
            However, the workaround is to use tpbconv -nsteps -1 (as
            well as the other
            stuff). Let me know how this goes and I'll update the
            documentation.<br>
            <br>
            <br>
            <o:p></o:p></span></p>
        <p class="MsoPlainText"><span style="font-size: 10pt;
            font-family: &quot;Courier New&quot;;">&nbsp;</span><o:p></o:p></p>
        <p class="MsoNormal">Using the g_energy command on the output
          energy.edr file, I
          got among others, these options to choose<o:p></o:p></p>
        <p class="MsoPlainText"><span style="font-size: 9pt;
            font-family: &quot;Courier New&quot;;">49&nbsp;
            Coul-SR:Protein-Protein&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
            50&nbsp;
            LJ-SR:Protein-Protein&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
          </span><o:p></o:p></p>
        <p class="MsoPlainText"><span style="font-size: 9pt;
            font-family: &quot;Courier New&quot;;">51&nbsp;
            Coul-14:Protein-Protein&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
            52&nbsp;
            LJ-14:Protein-Protein&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
          </span><o:p></o:p></p>
        <p class="MsoNormal">&nbsp;<o:p></o:p></p>
        <p class="MsoNormal">In order to get the energy of the protein,
          I reckon I have
          to add 49,50,51,52 (to account for the nonbonded components)
          and<o:p></o:p></p>
        <p class="MsoNormal"><span style="font-size: 12pt; font-family:
            &quot;Times New Roman&quot;,&quot;serif&quot;;"><br>
            I think that you can make a case either way for 1-4
            interactions - they're
            algorithmically similar to the other non-bonded
            interactions, but their
            parameter values should be tightly coupled to some other of
            the bonded
            parameters, but then in several forcefields those values are
            just scaled
            versions of the normal ones... I'd guess most people call
            them non-bonded.<br>
            <br>
            <br>
            <o:p></o:p></span></p>
        <p class="MsoPlainText"><span style="font-size: 9pt;
            font-family: &quot;Courier New&quot;;">1&nbsp;
            Angle&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 2&nbsp;
            G96Angle&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 3&nbsp;
            Proper-Dih.&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 4&nbsp; Ryckaert-Bell.&nbsp; 5&nbsp;
            Improper-Dih&nbsp; </span><o:p></o:p></p>
        <p class="MsoPlainText"><span style="font-size: 9pt;
            font-family: &quot;Courier New&quot;;">f</span>or
          the bonded components. However, I think 1-5 is the bonded
          terms for the system
          and not the protein alone. Can anyone help me with this?<o:p></o:p></p>
        <p class="MsoNormal"><span style="font-size: 12pt; font-family:
            &quot;Times New Roman&quot;,&quot;serif&quot;;"><br>
            Compare values from reruns on the subset-tpr and the
            full-tpr. I don't know
            whether/how well that works. In extremis, you should be able
            to reduce the .tpr
            to a single interaction.<br>
            <br>
            <br>
            <o:p></o:p></span></p>
        <p class="MsoPlainText">&nbsp;<o:p></o:p></p>
        <p class="MsoPlainText">&nbsp;<o:p></o:p></p>
        <p class="MsoPlainText">Also, on a slightly different note, this
          post <span style="font-size: 10pt; font-family: &quot;Courier
            New ;&quot;,&quot;serif&quot;;"><a moz-do-not-send="true"
href="http://oldwww.gromacs.org/pipermail/gmx-users/2005-July/016307.html">http://oldwww.gromacs.org/pipermail/gmx-users/2005-July/016307.html</a></span>
          suggested that the force constant of the solvent (DMSO) to be
          adjusted to zero.<span style="font-size: 10pt; font-family:
            &quot;Courier New&quot;;"> Am I right to think that it
            does not apply to my case as my protein-lipid system </span>is
          solvated with
          SPC? (I read that SPC is rigid water. I did not add &#8211;DFLEXIBLE
          in .mdp)<o:p></o:p></p>
        <p class="MsoNormal"><span style="font-size: 12pt; font-family:
            &quot;Times New Roman&quot;,&quot;serif&quot;;"><br>
            That was in the context of a full-tpr rerun. The point is
            that atoms that have
            been constrained don't contribute to relevant sums. Because
            you are using
            constraints, there's no "Bonds" energy sum for any atom
            pair. Because
            the DMSO was a flexible model, its bonded-interactions
            contributions need to be
            subtracted (i.e. parameters set to zero) to get a group-wise
            bonded-interaction
            value.<br>
            <br>
            Mark<br>
            <br>
            <br>
            <o:p></o:p></span></p>
        <p class="MsoPlainText">&nbsp;<o:p></o:p></p>
        <p class="MsoPlainText">Any help would be highly
          appreciated.Thanks!<o:p></o:p></p>
        <p class="MsoPlainText">&nbsp;<o:p></o:p></p>
        <p class="MsoPlainText">HW<o:p></o:p></p>
        <p class="MsoNormal"><span style="font-size: 12pt; font-family:
            &quot;Times New Roman&quot;,&quot;serif&quot;;">&lt;&lt;
            <o:p></o:p></span></p>
        <p><span style="font-size: 10pt; font-family:
            &quot;Arial&quot;,&quot;sans-serif&quot;;">This message
            and any attachment are intended solely for the addressee and
            may contain
            confidential information. If you have received this message
            in error, please
            send it back to me, and immediately delete it. Please do not
            use, copy or
            disclose the information contained in this message or in any
            attachment. Any
            views or opinions expressed by the author of this email do
            not necessarily
            reflect the views of the University of Nottingham. <br>
            <br>
            This message has been checked for viruses but the contents
            of an attachment may
            still contain software viruses which could damage your
            computer system: you are
            advised to perform your own checks. Email communications
            with the University of
            Nottingham may be monitored as permitted by UK &amp;
            Malaysia legislation.</span><o:p></o:p></p>
        <p class="MsoNormal"><span style="font-size: 12pt; font-family:
            &quot;Times New Roman&quot;,&quot;serif&quot;;">&gt;&gt;
            <o:p></o:p></span></p>
        <p class="MsoNormal"><span style="font-size: 12pt; font-family:
            &quot;Times New Roman&quot;,&quot;serif&quot;;"><o:p>&nbsp;</o:p></span></p>
      </div>
      &lt;&lt;
      <p><font face="Arial" size="2">
          This message and any attachment are intended solely for the
          addressee and may contain confidential information. If you
          have received this message in error, please send it back to
          me, and immediately delete it. Please do not use, copy or
          disclose the information contained in this message or in any
          attachment. Any views or opinions expressed by the author of
          this email do not necessarily reflect the views of the
          University of Nottingham. <br>
          <br>
          This message has been checked for viruses but the contents of
          an attachment may still contain software viruses which could
          damage your computer system: you are advised to perform your
          own checks. Email communications with the University of
          Nottingham may be monitored as permitted by UK &amp; Malaysia
          legislation.</font></p>
      &gt;&gt;
    </blockquote>
    <br>
  </body>
</html>