<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=us-ascii">
<meta name="Generator" content="Microsoft Word 12 (filtered medium)">
<style><!--
/* Font Definitions */
@font-face
{font-family:Helvetica;
panose-1:2 11 6 4 2 2 2 2 2 4;}
@font-face
{font-family:Helvetica;
panose-1:2 11 6 4 2 2 2 2 2 4;}
@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;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
margin-bottom:.0001pt;
font-size:12.0pt;
font-family:"Times New Roman","serif";}
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;}
span.EmailStyle17
{mso-style-type:personal-reply;
font-family:"Calibri","sans-serif";
color:#1F497D;}
.MsoChpDefault
{mso-style-type:export-only;}
@page WordSection1
{size:612.0pt 792.0pt;
margin:72.0pt 72.0pt 72.0pt 72.0pt;}
div.WordSection1
{page:WordSection1;}
--></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]-->
</head>
<body lang="EN-US" link="blue" vlink="purple">
<div class="WordSection1">
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">Sounds like the ligand is near the edge of the box, so only has to move short distance then does the jumping you are observing.<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><br>
The suggestion that has already been make, -pbc nojump, should stop that.<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">Another option is to center the box on a residue within the protein which is near the spot at which the ligand is interacting.<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Calibri","sans-serif";color:#1F497D">Catch ya,<br>
<br>
Dr. Dallas Warren<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Calibri","sans-serif";color:#1F497D">Medicinal Chemistry and Drug Action<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Calibri","sans-serif";color:#1F497D">Monash Institute of Pharmaceutical Sciences</span><span style="font-size:10.0pt;font-family:"Calibri","sans-serif";color:#1F497D">, Monash University<br>
381 Royal Parade, Parkville VIC 3010<br>
dallas.warren@monash.edu<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Calibri","sans-serif";color:#1F497D">+61 3 9903 9304<br>
---------------------------------<br>
When the only tool you own is a hammer, every problem begins to resemble a nail.</span><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">
<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p>
<div style="border:none;border-left:solid blue 1.5pt;padding:0cm 0cm 0cm 4.0pt">
<div>
<div style="border:none;border-top:solid #B5C4DF 1.0pt;padding:3.0pt 0cm 0cm 0cm">
<p class="MsoNormal"><b><span style="font-size:10.0pt;font-family:"Tahoma","sans-serif"">From:</span></b><span style="font-size:10.0pt;font-family:"Tahoma","sans-serif""> gmx-users-bounces@gromacs.org [mailto:gmx-users-bounces@gromacs.org]
<b>On Behalf Of </b>Steven Neumann<br>
<b>Sent:</b> Tuesday, 8 November 2011 8:30 PM<br>
<b>To:</b> jalemkul@vt.edu; Discussion list for GROMACS users<br>
<b>Subject:</b> Re: [gmx-users] Re: PBC - Protein - ligand<o:p></o:p></span></p>
</div>
</div>
<p class="MsoNormal"><o:p> </o:p></p>
<div>
<p class="MsoNormal">Thank you Justin, Mark and Tsjerk.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">I used the following workflow<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> <o:p></o:p></p>
</div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:18.0pt">
<span lang="PL" style="font-size:10.0pt;font-family:"Helvetica","sans-serif"">trjconv -s md.tpr -f md.xtc -o pbc_fix.xtc -pbc mol</span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:18.0pt">
<span lang="PL" style="font-size:10.0pt;font-family:"Helvetica","sans-serif";color:blue">trjconv -s md.tpr -f pbc_fix.xtc -n index.ndx -pbc cluster -o pbcfixcluster.xtc (Protein+ligand group)</span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:18.0pt">
<span style="font-size:10.0pt;font-family:"Helvetica","sans-serif"">trjconv -s md.tpr -f pbcfixcluster.xtc -o center.xtc -center (center on the protein)</span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:18.0pt">
<span style="font-size:10.0pt;font-family:"Helvetica","sans-serif""><br>
trjconv -s md.tpr -f center.xtc -o Cluster1.xtc -fit rot+trans (choose protein for
<br>
fitting)</span><o:p></o:p></p>
<div>
<p class="MsoNormal"><br>
Trajectory looks very good from the time when ligand stacked to the protein (90% of trajectory) but at the begining of the trajectory (when it is away from protein) it still jumps. I think that is the best solution I have found. If you know how to fix begining
please let me know.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">Steven<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><br>
<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">On Tue, Nov 8, 2011 at 8:53 AM, Steven Neumann <<a href="mailto:s.neumann08@gmail.com">s.neumann08@gmail.com</a>> wrote:<o:p></o:p></p>
<p class="MsoNormal" style="margin-bottom:12.0pt"><o:p> </o:p></p>
<div>
<div>
<p class="MsoNormal">On Mon, Nov 7, 2011 at 9:47 PM, Justin A. Lemkul <<a href="mailto:jalemkul@vt.edu" target="_blank">jalemkul@vt.edu</a>> wrote:<o:p></o:p></p>
<div>
<p class="MsoNormal"><br>
<br>
Steven Neumann wrote:<o:p></o:p></p>
<p class="MsoNormal" style="margin-bottom:12.0pt">Hi Tsjerk,<br>
<br>
Thank you. Unfortunately my ligand is not with protein. I put my ligand around my protein (in water) running separate simulations to see where can it bind. It is close to protein but not within. Any other suggestion?<br>
I used also pbc -res so I observe my ligand close to protein but sometimes still changing its position rapidly... No clue for now how to solve it...<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<p class="MsoNormal" style="margin-bottom:12.0pt">I have no idea why the proposed protocol isn't working, but I know that one should be able to do something very simple, along the lines of the following, for this to work:<br>
<br>
1. trjconv -s md.tpr -f md.xtc -o pbc_fix.xtc -pbc mol<br>
2. trjconv -s md.tpr -f pbc_fix.xtc -o center.xtc -center (center on the protein)<br>
3. trjconv -s md.tpr -f center.xtc -o fit.xtc -fit rot+trans (choose protein for fitting)<br>
<br>
-Justin<o:p></o:p></p>
<div>
<p class="MsoNormal"> <o:p></o:p></p>
</div>
</div>
<div>
<p class="MsoNormal">Thank you Justin. From this workflow my ligand is binding the protein most of the frames but sometimes it rapidly jumps to different part of the box and come back again. Then remains with protein and situation is repeated: in one frame
it changes its position and come back to protein remaining. :((((((((( no clue...<o:p></o:p></p>
</div>
<div>
<div>
<div>
<p class="MsoNormal"> <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> <o:p></o:p></p>
</div>
<blockquote style="border:none;border-left:solid #CCCCCC 1.0pt;padding:0cm 0cm 0cm 6.0pt;margin-left:4.8pt;margin-right:0cm">
<blockquote style="border:none;border-left:solid #CCCCCC 1.0pt;padding:0cm 0cm 0cm 6.0pt;margin-left:4.8pt;margin-right:0cm">
<p class="MsoNormal">Steven <o:p></o:p></p>
<div>
<p class="MsoNormal"><br>
<br>
On Monday, November 7, 2011, Tsjerk Wassenaar <<a href="mailto:tsjerkw@gmail.com" target="_blank">tsjerkw@gmail.com</a> <mailto:<a href="mailto:tsjerkw@gmail.com" target="_blank">tsjerkw@gmail.com</a>>> wrote:<br>
> Hi Steven,<br>
><br>
> Step 2: Cluster your molecules.<br>
> This is where you have to forge a reference frame that you can use to<br>
> remove jumps from your trajectory. If the ligand is not with the<br>
> protein at the start, you'll have to shift it so that it is. Maybe<br>
> -pbc cluster is your friend there. I do assume that the ligand is<br>
> really with the protein and not in the solvent...<br>
><br>
> Cheers,<br>
><br>
> Tsjerk<br>
><o:p></o:p></p>
</div>
<p class="MsoNormal"> > On Mon, Nov 7, 2011 at 5:17 PM, Steven Neumann <<a href="mailto:s.neumann08@gmail.com" target="_blank">s.neumann08@gmail.com</a> <mailto:<a href="mailto:s.neumann08@gmail.com" target="_blank">s.neumann08@gmail.com</a>>> wrote:
<o:p></o:p></p>
<div>
<div>
<p class="MsoNormal"><br>
>><br>
>><br>
>> On Mon, Nov 7, 2011 at 2:26 PM, Justin A. Lemkul <<a href="mailto:jalemkul@vt.edu" target="_blank">jalemkul@vt.edu</a> <mailto:<a href="mailto:jalemkul@vt.edu" target="_blank">jalemkul@vt.edu</a>>> wrote:<br>
>>><br>
>>><br>
>>> Steven Neumann wrote:<br>
>>>><br>
>>>> Dear Gmx Users,<br>
>>>> I know that this problem has been discussed may times but I cannot find<br>
>>>> the solution to get rid of pbc in my system: protein and ligand. I followed<br>
>>>> the workflow:<br>
>>>><br>
>>>> 1. First make your molecules whole if you want them whole<br>
>>>><br>
>>>> trjconv -f md.trr -s md.tpr -pbc whole -ur compact -o mdwhole.xtc<br>
>>>><br>
>>>> 2. Cluster your molecules/particles if you want them clustered<br>
>>>><br>
>>>> 3. Extract the first frame from the trajectory as reference for<br>
>>>> removing jumps if you want to remove jumps.<br>
>>>><br>
>>>> trjconv -f mdwhole.xtc -s md.tpr -dump 0 -o 1stframe.pdb<br>
>>>><br>
>>>> 4. Remove jumps if you want to have them removed using the first<br>
>>>> frame<br>
>>>><br>
>>>> trjconv -f mdwhole.xtc -s 1stframe.pdb -pbc nojump -o mdwholeNOjump.xtc<br>
>>>><br>
>>>> 5. Center your system using some criterion. Doing so shifts the<br>
>>>> system, so don't use |trjconv -|pbc| nojump| after this step.<br>
>>>><br>
>>>> trjconv -f mdwholeNOjump.xtc -center -o mdwholeNOjumpCENTER.xtc<br>
>>>><br>
>>>> 6. Put everything in some box.<br>
>>>><br>
>>>> trjconv -f mdwholeNOjumpCENTER.xtc -box 6 6 6 -o<br>
>>>> mdwholeNOjumpCENTERbox.xtc<br>
>>>><br>
>>>> 7. Fit if desired and don't use any PBC related option afterwards.<br>
>>>><br>
>>>> trjconv -f mdwholeNOjumpCENTERbox.xtc -s 1stframe.pdb -fit rot+trans -o<br>
>>>> mdfinal.xtc<br>
>>>><br>
>>>><br>
>>>> I used SYSTEM everywhere as output orinput. However, my ligand is still<br>
>>>> jumping like a fly around the stable protein. Do you have any suggestions?<br>
>>>><br>
>>>><br>
>>><br>
>>> Center on either the protein, the ligand, or some custom index group of<br>
>>> residues surrounding the ligand. Centering on the whole system usually<br>
>>> doesn't do anything useful.<br>
>>><br>
>>> -Justin<br>
>>><br>
>><br>
>> Thank you guys but...<br>
>><br>
>> I am trying and it does not work... my ligand is jumping like an idiot<br>
>> outside the box changing its position even two dimensions of box in one<br>
>> frame. I removed -ur compact from the first line and I tried centering on<br>
>> ligand or protein (centering group: LIG or Protein, output: SYSTEM). No<br>
>> results...<br>
>> My ligand at the begining of the simualtion is not within the protein.<br>
>> Please, help :(((( I tried this workflow with many ligands and same protein<br>
>> - it worked! Now it does not...<br>
>> Here is my workflow:<br>
>><br>
>><br>
>> 1. First make your molecules whole if you want them whole.<br>
>><br>
>> trjconv -f md.trr -s md.tpr -pbc whole -o mdwhole.xtc<br>
>><br>
>> 2. Cluster your molecules/particles if you want them clustered<br>
>><br>
>> 3. Extract the first frame from the trajectory as reference for<br>
>> removing jumps if you want to remove jumps.<br>
>><br>
>> trjconv -f mdwhole.xtc -s md.tpr -dump 0 -o 1stframe.pdb<br>
>><br>
>> 4. Remove jumps if you want to have them removed using the first frame<br>
>> (system)<br>
>><br>
>> trjconv -f mdwhole.xtc -s 1stframe.pdb -pbc nojump -o mdwholeNOjump.xtc<br>
>><br>
>> 5. Center your system using some criterion. Doing so shifts the system,<br>
>> so don't use trjconv -pbc nojump after this step (tried centering on LIG or<br>
>> PROTEIN)<br>
>><br>
>> trjconv -f mdwholeNOjump.xtc -n Ligand.ndx -center -o<br>
>> mdwholeNOjumpCENTER.xtc<br>
>><br>
>> 6. Put everything in some box.<br>
>><br>
>> trjconv -f mdwholeNOjumpCENTER.xtc -box 6 6 6 -o mdwholeNOjumpCENTERbox.xtc<br>
>><br>
>> 7. Fit if desired and don't use any PBC related option afterwards.<br>
>><br>
>> trjconv -f mdwholeNOjumpCENTERbox.xtc -s 1stframe.pdb -fit rot+trans -o<br>
> --<br>
> Tsjerk A. Wassenaar, Ph.D.<br>
><br>
> post-doctoral researcher<br>
> Molecular Dynamics Group<br>
> * Groningen Institute for Biomolecular Research and Biotechnology<br>
> * Zernike Institute for Advanced Materials<br>
> University of Groningen<br>
> The Netherlands<br>
> --<o:p></o:p></p>
</div>
</div>
<p class="MsoNormal"> > gmx-users mailing list <a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a> <mailto:<a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a>>
<o:p></o:p></p>
<div>
<p class="MsoNormal"><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<o:p></o:p></p>
</div>
<p class="MsoNormal"> > www interface or send it to <a href="mailto:gmx-users-request@gromacs.org" target="_blank">
gmx-users-request@gromacs.org</a> <mailto:<a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a>>.
<o:p></o:p></p>
<div>
<p class="MsoNormal" style="margin-bottom:12.0pt"><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>
><o:p></o:p></p>
</div>
</blockquote>
<p class="MsoNormal"><br>
-- <o:p></o:p></p>
<div>
<p class="MsoNormal">========================================<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]<a href="http://vt.edu/" target="_blank">vt.edu</a> | <a href="tel:%28540%29%20231-9080" target="_blank">
(540) 231-9080</a><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>
========================================<o:p></o:p></p>
</div>
<div>
<div>
<p class="MsoNormal">-- <br>
gmx-users mailing list <a href="mailto:gmx-users@gromacs.org" target="_blank">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" target="_blank">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><o:p></o:p></p>
</div>
</div>
</blockquote>
</div>
</div>
</div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
</div>
</body>
</html>