Matt,<br>
<br>
I agree with Berk -- I'm doing free energy calculations with the Pande
group AMBER port and get no problems. For electrostatic calculations I
generally get good agreement with what I get in the AMBER simulation
package itself (things are different for the vdW component since AMBER
doesn't have soft core) so it seems to be working properly. So I'd look
closely at your topologies and run input. What have you turned on
during the minimization you're running? Try it with free_energy=no;
also, have you tried using soft core or something? (It seems possible
that you might occasionally get strange things in minimization if you
use soft core with odd parameters). Also, make sure your torsion
parameters, etc. are set properly, and check for warnings. Again, I've
never seen anything like this, so I strongly suspect there's something
funny going on on your end.<br>
<br>
David<br>
<br><br><div><span class="gmail_quote">On 1/31/06, <b class="gmail_sendername"><a href="mailto:mernst@tricity.wsu.edu">mernst@tricity.wsu.edu</a></b> <<a href="mailto:mernst@tricity.wsu.edu">mernst@tricity.wsu.edu</a>
> wrote:</span><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">I am attempting to perform a free energy simulation that involves the<br>transformation of the lesioned DNA base 8-oxoguanine to guanine. Before I
<br>set up the B state in my .itp file, I was able to simulate 8-oxoguanine<br>with MD and its behavior looked to be in reasonable agreement with<br>simulations I have performed in another package.<br><br>Now that I have defined B-state parameters, MD and even energy
<br>minimization give strange and undesirable results, even before I begin the<br>free energy calculations or set lambda to a non-zero value.<br><br>I've narrowed it down to a simple test case where I supply a B state for
<br>just one atom.<br><br>If I use this line in my .itp file, which changes nothing between A and B,<br>the energy minimization gives a reasonable-looking structure:<br><br>46
amber99_35 2 8oG N7 46
-0.5129 14.01 amber99_35 -0.5129
14.01<br><br>This is what the output looks like:<br><br><a href="http://www.sciencemadness.org/scipics/goodtrimer.png">http://www.sciencemadness.org/scipics/goodtrimer.png</a><br><br>If I use this line, using a different atom type in the B state,
<br><br>46
amber99_35 2 8oG N7 46
-0.5129 14.01 amber99_36 -0.5129
14.01<br><br>the minimization gives me a hydrogen that sticks out at a strange angle:<br><br><a href="http://www.sciencemadness.org/scipics/badtrimer.png">http://www.sciencemadness.org/scipics/badtrimer.png</a><br><br>The amber99_35 atom is AMBER type NA, corresponding to sp2 N in 5
<br>memb.ring w/H atom (HIS). The amber99_36 atom is AMBER type NB,<br>corresponding to sp2 N in 5 memb.ring w/LP (HIS,ADE,GUA). In the A state<br>there is a hydrogen that becomes a dummy in the B state, so I would like<br>
to go from NA to NB. The AMBER force field port to Gromacs is from the<br>Pande group.<br><br>I cannot figure out why this effect shows up when lambda=0, and how I<br>might set up the free energy simulation so that it is not horribly
<br>distorted and prone to crashing in the equilibrium period.<br><br>Matt Ernst<br>Washington State University<br><br>_______________________________________________<br>gmx-users mailing list <a href="mailto:gmx-users@gromacs.org">
gmx-users@gromacs.org</a><br><a href="http://www.gromacs.org/mailman/listinfo/gmx-users">http://www.gromacs.org/mailman/listinfo/gmx-users</a><br>Please don't post (un)subscribe requests to the list. Use the<br>www interface or send it to
<a href="mailto:gmx-users-request@gromacs.org">gmx-users-request@gromacs.org</a>.<br>Can't post? Read <a href="http://www.gromacs.org/mailing_lists/users.php">http://www.gromacs.org/mailing_lists/users.php</a><br></blockquote>
</div><br>