<html><body><div style="color:#000; background-color:#fff; font-family:times new roman, new york, times, serif;font-size:12pt"><div><span>Thanks.</span></div><div><span></span>&nbsp;</div><div><span>The original PDB file (2BEG) has 5 chains. In your tutorial you fix chain B. Thus if I am right, the energy calculated should be the disintegation energy of the whole protein to the 5 chain peptide. In addition the protein chain in the original pdb is from residue 17 to 42. </span></div><div><span></span>&nbsp;</div><div><span>This makes me confused when you define "1-27" to chain A, and "28-54" to chain B.</span></div><div><span></span>&nbsp;</div><div><span>It seems your tutorial protein is a 2 chain protein, not the 5 chain protein.</span></div><div><span></span>&nbsp;</div><div><span>Can you give me further
 clarification?</span></div><div><span></span>&nbsp;</div><div><span>Cheers,</span></div><div><span></span>&nbsp;</div><div><span>Acoot</span></div>&nbsp; <div style="font-family: times new roman, new york, times, serif; font-size: 12pt;"> <div style="font-family: times new roman, new york, times, serif; font-size: 12pt;"> <div dir="ltr"> <font size="2" face="Arial"> <div style="margin: 5px 0px; padding: 0px; border: 1px solid rgb(204, 204, 204); height: 0px; line-height: 0; font-size: 0px;" class="hr" contentEditable="false" readonly="true"></div>  <b><span style="font-weight: bold;">From:</span></b> Justin A. Lemkul &lt;jalemkul@vt.edu&gt;<br> <b><span style="font-weight: bold;">To:</span></b> Discussion list for GROMACS users &lt;gmx-users@gromacs.org&gt; <br> <b><span style="font-weight: bold;">Sent:</span></b> Sunday, 1 April 2012 10:49 PM<br> <b><span style="font-weight: bold;">Subject:</span></b> Re: [gmx-users] on Umbrella Sampling<br> </font>
 </div> <br><br><br>Acoot Brett wrote:<br>&gt; Here I also mean how can I get the index.npt in the "grompp -f md_pull.mdp -c npt.gro -p topol.top -n index.ndx -t npt.cpt -o pull.tpr".<br>&gt;&nbsp; Looking forward to getting your advice.<br>&gt;&nbsp; <br><br>That file is produced by make_ndx.&nbsp; It is written when you type 'q' and Enter to quit the program.&nbsp; If you have not specified any other output name, the default is "index.ndx."<br><br>-Justin<br><br>&gt; Cheers,<br>&gt;&nbsp; Acoot<br>&gt; <br>&gt; *From:* Justin A. Lemkul &lt;<a href="mailto:jalemkul@vt.edu" ymailto="mailto:jalemkul@vt.edu">jalemkul@vt.edu</a>&gt;<br>&gt; *To:* Acoot Brett &lt;<a href="mailto:acootbrett@yahoo.com" ymailto="mailto:acootbrett@yahoo.com">acootbrett@yahoo.com</a>&gt;; Discussion list for GROMACS users &lt;<a href="mailto:gmx-users@gromacs.org" ymailto="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a>&gt;<br>&gt; *Sent:* Sunday, 1 April 2012 10:32
 PM<br>&gt; *Subject:* Re: [gmx-users] on Umbrella Sampling<br>&gt; <br>&gt; <br>&gt; <br>&gt; Acoot Brett wrote:<br>&gt;&nbsp; &gt; Dear All,<br>&gt;&nbsp; &gt;&nbsp; In the on-line tutorial "Umbrella Sampling " by*Justin Lemkul in "<a href="http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin/gmx-tutorials/umbrella/05_pull.html" target="_blank">http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin/gmx-tutorials/umbrella/05_pull.html</a>"* , there is a command "grompp -f md_pull.mdp -c npt.gro -p topol.top -n index.ndx -t npt.cpt -o pull.tpr". Here is the index.ndx exactly same as following?<br>&gt;&nbsp; &gt;&nbsp; make_ndx -f npt.gro<br>&gt;&nbsp; &gt; (&gt; indicates the make_ndx prompt)<br>&gt;&nbsp; &gt;&nbsp; &gt; r 1-27<br>&gt;&nbsp; &gt;&nbsp; &gt; name 19 Chain_A<br>&gt;&nbsp; &gt;&nbsp; &gt; r 28-54<br>&gt;&nbsp; &gt;&nbsp; &gt; name 20 Chain_B<br>&gt;&nbsp; &gt;&nbsp; &gt; q<br>&gt;&nbsp; &gt;&nbsp; Then what is the meaning of "r1-27",
 what is the meaning of " name 19 Chain_A", and what is the meaning of "q"?<br>&gt;&nbsp; &gt; <br>&gt; Type "help" at the make_ndx prompt for an explanation.<br>&gt; <br>&gt;&nbsp; &gt; Suppose I have a protein-peptide complex A, the protein is of 400 aas,the peptide is 20 aas. Suppose I have another protein-peptide complx B, and the only difference beteen protein complex A and complex B is that in position 3 of peptide of complex A, it is Ala. In the coresponsing position of the peptide of complex B, it is Arg.<br>&gt;&nbsp; &gt;&nbsp; Can the results of this Umbrella Sampling calculation be accurate enough to show that the binding energy between the protein and peptide in complex A and B is different and cuurently reflects the fact that the difference is caused by the substitution of Ala by Arg?<br>&gt;&nbsp; &gt; <br>&gt; Probably.<br>&gt; <br>&gt; -Justin<br>&gt; <br>&gt; -- ========================================<br>&gt; <br>&gt; Justin A.
 Lemkul<br>&gt; Ph.D. Candidate<br>&gt; ICTAS Doctoral Scholar<br>&gt; MILES-IGERT Trainee<br>&gt; Department of Biochemistry<br>&gt; Virginia Tech<br>&gt; Blacksburg, VA<br>&gt; jalemkul[at]vt.edu | (540) 231-9080<br>&gt; <a href="http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin" target="_blank">http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin</a><br>&gt; <br>&gt; ========================================<br>&gt; <br>&gt; <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><a href="http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin" target="_blank">http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin</a><br><br>========================================<br>-- gmx-users mailing list&nbsp; &nbsp; <a href="mailto:gmx-users@gromacs.org"
 ymailto="mailto: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/Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/Support/Mailing_Lists/Search</a> before posting!<br>Please don't post (un)subscribe requests to the list. Use the www interface or send it to <a href="mailto:gmx-users-request@gromacs.org" ymailto="mailto:gmx-users-request@gromacs.org">gmx-users-request@gromacs.org</a>.<br>Can't post? Read <a href="http://www.gromacs.org/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/Support/Mailing_Lists</a><br><br><br> </div> </div>  </div></body></html>