The structure in the original PDB file of ethylene glycol was in the gauche conformation.  When I minimised the solvated system, the molecule remained in the gauche conformation.  However, during equilibration of the system, the conformation of the molecule changed to trans within less than one picosecond.  It remained in the trans conformation for the rest of the equilibration.<br>
<br><br>Thanks,<br>Nancy<br><br><br><br><div class="gmail_quote">On Sun, Aug 9, 2009 at 5:49 PM, Nagy, Peter I. <span dir="ltr">&lt;<a href="mailto:PNAGY@utnet.utoledo.edu">PNAGY@utnet.utoledo.edu</a>&gt;</span> wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">






<div>

<br>

<p><font size="2">Dear Nancy,<br>
<br>
  It was useful to find charges for the solute... I am not a GROMACS user, so I do not know what checks have to be done for a non-protein or non-DNA solute. Anyway, this case helped develop data management for a non-standard solute.<br>

  It is normal to use united atom model for CH, CH2 and CH3 in GROMACS. As far as the charges are concerned, it is always a matter of speculation. Nobody knows what charges are &quot;correct&quot;. Only a quantum MD (for the whole system) helps, where you need not derive charges since you do not use empirical pair potentials. Charges of -0.574 and 0.398 for the OH atoms are reasonable, although somewhat larger absolute values could also be used. Nonetheless, they already triggered the solute-solvent H-bonds, which was the Nr. 1 problem.<br>

  I do not know what conformation you derived from the web. Was it ORIGINALLY trans or gauche? As I wrote, the experiment found only 12% trans in aqueous solution. Try to start the structure from the so-called tGg&#39; conformation or some other gauche(you may find them e.g. in our papers JACS,1991, 113, 6719; 1992, 114, 4752 and a complete set of conformations in the paper of Cramer and Truhlar in the JACS about 1995, I do not have the reference handy). If the system stays in gauche, then the simulation has been short starting from trans and/or the OCCO torsion barrier is too high for equilibration. I would start from the tGg&#39; and would conduct at least 500 ps long simulation, check if trans appears and what is the gauche/trans ratio. If you do not find a trans starting from gauche or no gauche starting from trans, then I would reduce the OCCO torisonal. Also check whether the intramolcular hydrogen bonds remain maintained in the gauche structures. Hopefully these above help.<br>

  The ratio is a sensitive problem, and could depend on the applied charge sets, too. Thus you may need to conduct a more complicated test finding combination(s) of charges and torsion parameters, which produce good results. Maybe there are different combinations working equally well. I did not meet this problem when I performed the free energy perturbation calculations via Monte Carlo simulations with rigid solute geometries.<br>

<br>
Best regards<br><font color="#888888">
Peter Nagy  </font><div><div></div><div class="h5"><br>
<br>
</div></div></font>
</p>

</div>
</blockquote></div><br><br><br><br>