Dear Mark,<br><br>Sorry, perhaps I should have clarified more. I did not just make a blind substitution of a choline by a H-atom.<br><br>First, I did conserve charge.<br>
Second, I did introduce reasonable non-bonded and bonded interaction parameters
for the new H-atom. The non-bonded interaction parameters were similar to those suggested by
the mono-anionic phosphate group on a threonine or serine.<br><br>I guess you will agree force fields are often set up that entire functional groups from one molecule can be transferred to others. For example, one can use the same head group for both POPC and DPPC, or conversely, one can use different types of hydrocarbon chains for the same head group, from 12 carbons to 20 carbons, as has often been done. The only parameters that might change in getting POPA from POPC are near the H-atom itself, and those are the ones that I had posted earlier, just to see if anyone could notice a mistake.<br>
<br>Similar methods have been used in literature earlier, except that a DFT calculation was made for the mono-anionic phosphate group to get new charges. The DFT did reveal that charges on DOPA were similar to DOPC.<br><br>
The main issue I have is that I get a somewhat skewed phosphate group that leads to unfavorable positions for the oxygen atoms on phosphate. I am guessing that this is because of the bonded-interactions of the new H-atom.<br>
<br><br><div class="gmail_quote">On Mon, May 5, 2008 at 3:10 PM, Mark Abraham <<a href="mailto:Mark.Abraham@anu.edu.au">Mark.Abraham@anu.edu.au</a>> 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 class="Ih2E3d">himanshu khandelia wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
Hi,<br>
<br>
I am trying to run a simulation of a POPC bilayer mixed with some<br>
mono-anionic Phosphatidic acid (POPA) , where the choline group is<br>
replaced by a hydrogen atom. However, the energy of my simulation box<br>
diverges, and I am trying to fix the problem. I also tried running a<br>
simulation of PA solvated in water and a single sodium ion, but even<br>
that does not work. I will try to provide as much detail as possible,<br>
I hope someone can help point out some obvious error which I have not<br>
been able to debug over 3 days. In the simulation of a single PA<br>
molecule in water, MDRUN keeps complaining that the distance between<br>
atoms H1 and O4 is very high (> 4 nm or so). The problem probably lies<br>
in my topology, possibly in the bonded interactions?<br>
</blockquote>
<br></div>
Merely replacing a functional group with a hydrogen doesn't make for a sensible model physics. First, you probably didn't conserve charge, and second the parameters don't necessarily make sense for a particle of mass 1 rather than mass 12 (plus secondary effect of substituent group), and third there's no way to fix both the foregoing and the other issues.<br>
<br>
You need at least to find a proper model of POPA and to implement it, or if need be, develop one yourself. Search the wiki for Parameterization to get a start on the latter, but it is an expert topic.<br><font color="#888888">
<br>
Mark</font><div><div></div><div class="Wj3C7c"><br>
_______________________________________________<br>
gmx-users mailing list <a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br>
<a href="http://www.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://www.gromacs.org/mailman/listinfo/gmx-users</a><br>
Please search the archive at <a href="http://www.gromacs.org/search" target="_blank">http://www.gromacs.org/search</a> before posting!<br>
Please don't post (un)subscribe requests to the list. Use the www interface or send it to <a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a>.<br>
Can't post? Read <a href="http://www.gromacs.org/mailing_lists/users.php" target="_blank">http://www.gromacs.org/mailing_lists/users.php</a><br>
</div></div></blockquote></div><br><br clear="all"><br>-- <br>Maria G.<br>Technical University of Denmark<br>Copenhagen