<br><br>----- Original Message -----<br>From: Roland Schulz &lt;roland@utk.edu&gt;<br>Date: Thursday, August 26, 2010 3:06<br>Subject: Re: [gmx-users] How to tell if the molecule is correctly reduced?<br>To: Discussion list for GROMACS users &lt;gmx-users@gromacs.org&gt;<br><br><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font>Hi,<div><br></div><div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font>if the gro file written by pdb2gmx contains the same number of hydrogens then before, than you&nbsp;disulfide&nbsp;bonds haven't changed. pdb2gmx automatically forms disulfide bridges if the atoms are within some distance (see specbond). Look at the output of pdb2gmx and make sure it is doing what you expect/want.</div>  <div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font><br></div><div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font>Roland<br><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font><br><div class="gmail_quote"><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font>On Wed, Aug 25, 2010 at 12:44 PM, Rabab Toubar <span dir="ltr">&lt;<a href="javascript:main.compose('new', 't=rtoubar@yahoo.com')">rtoubar@yahoo.com</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;"><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font>  Hi<br> <font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font><br><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font> I did simulations for a protein with disulfide bonds using opls. Then I reduced the ssbonds in vmd and did simulations for the reduced. When loading gro(or pdb) and trr files on vmd I do not see reduced bonds. I also compared the two pdb files and they looked the same in terms of number of hydrogens in cysteines. So how can I tell if the molecule is properly reduced before going further.<br></blockquote></div></div><br>It can't be "properly reduced" by VMD. The coordinates and the topology both have to change, and VMD neither reads or writes .tpr and .top.<br><br>Instead, use pdb2gmx -ss to allow you to choose which disulfides are made in the first place.<br><br>Mark