Hi Guillem,<br><br>First regarding the force field. You should not use the ffgmx one. That is basically the old GROMOS87 force field, which has lots of errors in it. Rather take one of the newer GROMOS96 force fields, or maybe an all atom one, such as OPLS. That does pose you the problem of getting the right parameters for the ligand.., sorry.
<br><br>As for the rectangular rhombic dodecahedron.., that's explained in the manual in Chapter 3 (periodic boundary conditions). Every box can be represented as a triclinic box or a rectangular box with exactly equal content. This also seems the cause of your molecule being split. Every part is oriented such that the first atom is placed inside the rectangular box you see. That has all to do with efficiency in the MD code. It does seem to suggest though, that you have more subunits in the protein. Should it really be one piece? How many *_?.itp files did you obtain after running pdb2gmx? If there's more than one, while it is actually one chain, run pdb2gmx with the option -merge
<br><br>If some parts of this don't make sense, feel free to ask further... (but also consider the manual and the user list archive).<br><br>Cheers,<br><br>Tsjerk<br><br><br><div><span class="gmail_quote">On 6/19/06, <b class="gmail_sendername">
Olli Lehtonen</b> <<a href="mailto:olehtone@focus.hut.fi">olehtone@focus.hut.fi</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;">
<br><br>On Mon, 19 Jun 2006, Guillem Plasencia wrote:<br><br>> Then i was curious and did trjconv my system (before running mdrun energy<br>> minimization) into a pdb, and had a look at it, and i saw my protein was<br>
> splitted in three different molecules, one was inside the water box (but far<br>> away from box center, lying in touch of one of the boxes' sides) and the<br>> other two parts were completely outside the box. Worst of all, they were no
<br>> more one single peptidic chain, but separated one.<br><br>I was having somewhat similar problem recently when I had multiple<br>peptides in my system (I used genconf to generate them) but in the final<br>structure they all had the same chain identifier, so pdb2gmx thought that
<br>I have only one long peptide in the system (with some 1 nm long bonds).<br><br>If this is the reason, one way to check how many chains there are is to<br>interactively choose termini groups with pdb2gmx, you should be prompted
<br>for each peptide unit.<br><br>By the way, if it is somehow possible to tell genconf to generate repeated<br>units with different chain identifiers, I would be interested to know.<br><br>regards,<br><br>Olli<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><br clear="all"><br>-- <br><br>Tsjerk A. Wassenaar, M.Sc.<br>Groningen Biomolecular Sciences and Biotechnology Institute (GBB)<br>Dept. of Biophysical Chemistry<br>University of Groningen<br>
Nijenborgh 4<br>9747AG Groningen, The Netherlands<br>+31 50 363 4336<br>