<html>

<head>

<style><!--

.hmmessage P

{

margin:0px;

padding:0px

}

body.hmmessage

{

font-size: 10pt;

font-family:Verdana

}

--></style>

</head>

<body class='hmmessage'>

Dear GROMACS users, <br><br>I have a protein structure with two oppositely charged residues facing each other. They are close enough for an electrostatic interaction to occur between them. However, this does not occur. After minimization they are in the exact same positions as before. <br><br>Is there some way to confirm that an electrostatic interaction is not 

being made in the structure with two charged residues facing each other?

 And if this is the case, do you have any idea as to why this is 

happening?<br><br>Perhaps the following information will prove helpful: <br><br><ul><li>I have another structure where there is only one charge (GLU) and one uncharged residue in the same positions as the charged residues in the previous structure. The coulomb energy values of the two structures after minimization do not show much difference (-3.9e+04 compared to -4.0e+4, where -4.0e+04 corresponds to the structure with two charges facing each other). Theoretically, this difference should be greater as the unsatisfied single charge is not as stable. To me the similar coloumb energies suggest that the electrostaic interactions are not being made in the structure with two opposite charges. </li></ul><br><ul><li>My .mdp file is as follows:</li></ul>&nbsp;&nbsp;&nbsp; &nbsp; ;<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ;&nbsp;&nbsp;&nbsp; User spoel (236)<br>&nbsp;&nbsp;&nbsp; &nbsp; ;&nbsp;&nbsp;&nbsp; Wed Nov&nbsp; 3 17:12:44 1993<br>&nbsp;&nbsp;&nbsp; &nbsp; ;&nbsp;&nbsp;&nbsp; Input file<br>&nbsp;&nbsp; &nbsp;&nbsp; ;<br>&nbsp;&nbsp; &nbsp;&nbsp; cpp&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; /usr/bin/cpp<br>&nbsp;&nbsp;&nbsp; &nbsp; define&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; -DFLEXIBLE <br>&nbsp;&nbsp; &nbsp;&nbsp; constraints&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; none<br>&nbsp;&nbsp; &nbsp;&nbsp; integrator&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; steep<br>&nbsp;&nbsp;&nbsp; &nbsp; nsteps&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; 1000<br>&nbsp;&nbsp; &nbsp; ;<br>&nbsp;&nbsp;&nbsp;&nbsp; ;&nbsp;&nbsp;&nbsp; Energy minimizing stuff<br>&nbsp;&nbsp; &nbsp; ;<br>&nbsp;&nbsp; &nbsp; emtol&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; 2000<br>&nbsp;&nbsp; &nbsp; emstep&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; 0.01<br><br>&nbsp; &nbsp; nstcomm&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; 1<br>&nbsp;&nbsp;&nbsp; ns_type&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; grid<br>&nbsp; &nbsp; rlist&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; 1<br>&nbsp; &nbsp; rcoulomb&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; 1.0<br>&nbsp; &nbsp; rvdw&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; 1.0<br>&nbsp; &nbsp; Tcoupl&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; no<br>&nbsp; &nbsp; Pcoupl&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; no<br>&nbsp; &nbsp; gen_vel&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; no<br><br><br>Thank you in advance!<br>Abdullah&nbsp; <br><br>&nbsp;<br>                                               <br /><hr />Hotmail: Powerful Free email with security by Microsoft. <a href='https://signup.live.com/signup.aspx?id=60969' target='_new'>Get it now.</a></body>

</html>