Of course I can, but the only way to do that is to capture the pdb2gmx or grompp output as it seems...<br><br><div><span class="gmail_quote">On 12/21/05, <b class="gmail_sendername">Florian Haberl</b> &lt;<a href="mailto:Florian.Haberl@chemie.uni-erlangen.de">
Florian.Haberl@chemie.uni-erlangen.de</a>&gt; wrote:</span><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">hi,<br><br>On Wednesday 21 December 2005 17:53, Tamas Horvath wrote:
<br>&gt; The thing is, I'm working with literally thousands of proteins, and running<br>&gt; some basic simulations on them. Therefore I need to run these simulations<br>&gt; using scripts, so I guess u can't expect me to know every charge
<br>&gt; distribution... These simulations may well be worthless and bogus, but I try<br>&gt; to do my best to make them meaningful...<br><br>you can change your script&nbsp;&nbsp;so it can automatically neutralise your systems.<br>
<br><br><br>&gt;<br>&gt; On 12/21/05, Mark Abraham &lt;<a href="mailto:Mark.Abraham@anu.edu.au">Mark.Abraham@anu.edu.au</a>&gt; wrote:<br>&gt; &gt;<br>&gt; &gt; Tamas Horvath wrote:<br>&gt; &gt; &gt; I have a protein in water system. The protein has some charge, which I
<br>&gt; &gt; &gt; want to neutralise, introducing stochastic ammount of K+ or Cl-. As far<br>&gt; &gt; &gt; as I figured out, the best way to do it is to use genion. Only I don't<br>&gt; &gt; &gt; know beforehand what the protein's charge would be, so I'm wondering if
<br>&gt; &gt; &gt; it's possible to tell genion to neutralise the system?<br>&gt; &gt;<br>&gt; &gt; Not knowing the pattern of charge distribution in your protein, and thus<br>&gt; &gt; the net charge, sounds to me like a great way to run the risk of doing
<br>&gt; &gt; some worthless simulations and/or bogus science for want of a few<br>&gt; &gt; minutes of observation and thought. I'd call the gromacs utility<br>&gt; &gt; programs' inability to automatically neutralise a system a feature, not
<br>&gt; &gt; a bug!<br>&gt; &gt;<br>&gt; &gt; Cheers,<br>&gt; &gt;<br>&gt; &gt; Mark Abraham<br>&gt; &gt; _______________________________________________<br>&gt; &gt; gmx-users mailing list<br>&gt; &gt; <a href="mailto:gmx-users@gromacs.org">
gmx-users@gromacs.org</a><br>&gt; &gt; <a href="http://www.gromacs.org/mailman/listinfo/gmx-users">http://www.gromacs.org/mailman/listinfo/gmx-users</a><br>&gt; &gt; Please don't post (un)subscribe requests to the list. Use the
<br>&gt; &gt; www interface or send it to <a href="mailto:gmx-users-request@gromacs.org">gmx-users-request@gromacs.org</a>.<br>&gt; &gt;<br>&gt;<br><br>Greetings,<br><br>Florian<br><br>--<br>-------------------------------------------------------------------------------
<br> Florian Haberl<br> Computer-Chemie-Centrum<br> Universitaet Erlangen/ Nuernberg<br> Naegelsbachstr 25<br> D-91052 Erlangen<br> Mailto: florian.haberl AT <a href="http://chemie.uni-erlangen.de">chemie.uni-erlangen.de</a>
<br>-------------------------------------------------------------------------------<br>_______________________________________________<br>gmx-users mailing list<br><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></blockquote></div><br>