<html dir="ltr">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style id="owaParaStyle" type="text/css">P {margin-top:0;margin-bottom:0;}</style>
</head>
<body ocsi="0" fpstyle="1">
<div style="direction: ltr;font-family: Tahoma;color: #000000;font-size: 10pt;">Hi everyone,
<br>
<br>
<br>
I'm a new Gromacs user and I bumped into a problem with the parameter assignment to the constraint bonds in COH and CSH groups in serines and cysteines when flag
<i>-vsite hydrogens</i> is being used in the topology preparation. Namely, after running an energy minimisation on the system built with virtual sites, I noticed that CSH angles in cysteines ended up with a strange pointy form corresponding to 75 degrees, instead
 of the equilibrium value of 96. After a bit of a digging, it turned out that the problem is created during the grompp run, where only a single function is created to represent the constraint bond that prevents angle bending for both COH and CSH groups. There
 parameter used in that particular function corresponds to the COH group (0.198 nm), which is too short for CSH (should be 0.237 nm) and results with the distorted CSH angle. The reason why this happens seems to be the vague definition of the constraint bond
 given through the atom types (CH2 and H), which are unfortunately the same for both CYS and SER. In addition,
<font size="2" color="black" face="Tahoma"><span style="font-size:10pt;" dir="ltr">the CH2-H distance parameter for CSH is actually completely omitted from the current ffG54a7bon.file!
</span></font>There is no problem with other COH groups, which are defined through different atom types, as shown below in the excerpt from ffG54a7bon.itp file with parts of interest highlighted:<br>
<br>
&nbsp;<font size="2" color="black" face="Tahoma"><span style="font-size:10pt;" dir="ltr">*************************************************************************<br>
; get the constraint distances for dummy atom constructions<br>
<br>
#include &quot;ff_dum.itp&quot;<br>
<br>
[ constrainttypes ]<br>
; now the constraints for the rigid NH3 groups<br>
&nbsp;MNH3&nbsp;&nbsp;&nbsp; C&nbsp;&nbsp;&nbsp; 2&nbsp;&nbsp; DC_MNC1<br>
&nbsp;MNH3&nbsp; CH1&nbsp;&nbsp;&nbsp; 2&nbsp;&nbsp; DC_MNC2<br>
&nbsp;MNH3&nbsp; CH2&nbsp;&nbsp;&nbsp; 2&nbsp;&nbsp; DC_MNC2<br>
&nbsp;MNH3 MNH3&nbsp;&nbsp;&nbsp; 2&nbsp;&nbsp; DC_MNMN<br>
; and the angle-constraints for OH and SH groups in proteins:<br>
<span style="background-color:#FFCC99;">&nbsp; CH2&nbsp;&nbsp;&nbsp; H&nbsp;&nbsp;&nbsp; 2&nbsp;&nbsp; DC_CO</span><br>
&nbsp; CH1&nbsp;&nbsp;&nbsp; H&nbsp;&nbsp;&nbsp; 2&nbsp;&nbsp; DC_CO<br>
&nbsp;&nbsp;&nbsp; C&nbsp;&nbsp;&nbsp; H&nbsp;&nbsp;&nbsp; 2&nbsp;&nbsp; DC_CO<br>
&nbsp;&nbsp;&nbsp; P&nbsp;&nbsp;&nbsp; H&nbsp;&nbsp;&nbsp; 2&nbsp;&nbsp; DC_PO<br>
&nbsp;<span style="background-color:#FFCC99;"> ????????????????????? missing DC_CS value</span><br>
**************************************************************************</span></font><br>
<br>
I'm not really sure if this belongs to the developers or to the users mailing list, but I gues that depends on the solution of the problem. It can be probably done on a more general level by redefining the constraints using more unique atom descriptors or just
 by calculating the atom distance for each constraint of interest by looking at the atoms directly and not reading in a fixed value. As I have no experience in coding in C, this goes way beyond my programming skills so I came up with an easier solution - introducing
 additional atom type for hydrogen bound to sulphur named HS. With this, the grompp was able to override this issue and this seems to be working:<br>
<br>
<font size="2" color="black" face="Tahoma"><span style="font-size:10pt;" dir="ltr">**************************************************************************<br>
; get the constraint distances for dummy atom constructions<br>
<br>
#include &quot;ff_dum.itp&quot;<br>
<br>
[ constrainttypes ]<br>
; now the constraints for the rigid NH3 groups<br>
&nbsp;MNH3&nbsp;&nbsp;&nbsp; C&nbsp;&nbsp;&nbsp; 2&nbsp;&nbsp; DC_MNC1<br>
&nbsp;MNH3&nbsp; CH1&nbsp;&nbsp;&nbsp; 2&nbsp;&nbsp; DC_MNC2<br>
&nbsp;MNH3&nbsp; CH2&nbsp;&nbsp;&nbsp; 2&nbsp;&nbsp; DC_MNC2<br>
&nbsp;MNH3 MNH3&nbsp;&nbsp;&nbsp; 2&nbsp;&nbsp; DC_MNMN<br>
; and the angle-constraints for OH and SH groups in proteins:<br>
<span style="background-color:#FFCC99;">&nbsp; CH2&nbsp;&nbsp;&nbsp; H&nbsp;&nbsp;&nbsp; 2&nbsp;&nbsp; DC_CO</span><br>
&nbsp; CH1&nbsp;&nbsp;&nbsp; H&nbsp;&nbsp;&nbsp; 2&nbsp;&nbsp; DC_CO<br>
&nbsp;&nbsp;&nbsp; C&nbsp;&nbsp;&nbsp; H&nbsp;&nbsp;&nbsp; 2&nbsp;&nbsp; DC_CO<br>
&nbsp;&nbsp;&nbsp; P&nbsp;&nbsp;&nbsp; H&nbsp;&nbsp;&nbsp; 2&nbsp;&nbsp; DC_PO<br>
<span style="background-color:#FFCC00;">&nbsp; CH2&nbsp;&nbsp; HS&nbsp;&nbsp;&nbsp; 2&nbsp;&nbsp; DC_CS</span><br>
<br>
************************************************************************<br>
<br>
<br>
I duplicated all the parameters given for atom type H and reassigned them to HS type. The modified force field files (version G54a7) are in the attachment and have suffix _HS.<br>
<br>
I was trying to find if this problem was already known, but I couldn't find anything related to this issue. In any case, it would be good to fix this as many proteins contain cysteines, but in case it has already been fixed, I apologize for spamming.<br>
<br>
All the best, <br>
<br>
Karmen<br>
<br>
</span></font>
<div><br>
<div style="font-family:Tahoma; font-size:13px">
<div style="font-family:Tahoma; font-size:13px">
<div style="font-family:Tahoma; font-size:13px">
<div style="font-family:Tahoma; font-size:13px">
<div style="font-family:Tahoma; font-size:13px">
<div><font size="2" face="Tahoma"><br>
-----------------------------<br>
Dr. Karmen Condic-Jurkic</font><font size="2" face="Tahoma"><span class="signature"><br>
Postdoctoral Researcher <br>
School of Chemistry and Molecular Biosciences<br>
The University of Queensland</span></font><font size="2" face="Tahoma"><span class="signature"><font size="2" face="Tahoma"><span class="signature"><br>
</span></font>Brisbane<font size="2"> St. Luc<font size="2">ia, </font>QLD</font> 407<font size="2">2<br>
Aus<font size="2">tralia</font></font> </span></font></div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</body>
</html>