<div>It did work. Thank you!</div>
<div> </div>
<div>Cheers,</div>
<div>Andrea<br><br></div>
<div class="gmail_quote">2008/10/23 Mark Abraham <span dir="ltr"><<a href="mailto:Mark.Abraham@anu.edu.au">Mark.Abraham@anu.edu.au</a>></span><br>
<blockquote class="gmail_quote" style="PADDING-LEFT: 1ex; MARGIN: 0px 0px 0px 0.8ex; BORDER-LEFT: #ccc 1px solid">
<div class="Ih2E3d">Andrea Muntean wrote:<br>
<blockquote class="gmail_quote" style="PADDING-LEFT: 1ex; MARGIN: 0px 0px 0px 0.8ex; BORDER-LEFT: #ccc 1px solid">Thank you for your answer.<br> If I am looking at the PDB standard format and my file, I observe that I have only one number between the type of residue and coordinates, and I suppose it is the residue current number. I am indeed missing the chain identifier code, which should be of A1 Fortran format, which I don't really understand (one carachter, but I have 32 chains, so it should be minimum A2 or I2). That I could introduce.<br>
</blockquote><br></div>Get a PDB file with just one polymer molecule and work on getting a topology for that. Walking before running :-)
<div class="Ih2E3d"><br><br>
<blockquote class="gmail_quote" style="PADDING-LEFT: 1ex; MARGIN: 0px 0px 0px 0.8ex; BORDER-LEFT: #ccc 1px solid">In my pdb entries I dont define any connectors. Should they be a dummy atom, with the corresponding bond, angle, dihedhral, which after that I redefine in the tdb file?<br>
</blockquote><br></div>First you need a functional *residue topology*. Check chapter 5 of the manual and look at the other force fields to see how (for example) protein force fields define the connectors with + and - prefixes. Once you've got that, then it's possible for pdb2gmx to join them head to tail (using the residue number field to indicate the breaks, and the name to match up with the topology you've defined), and you may then need either terminating residues of a different name or a .tdb entry.
<div class="Ih2E3d"><br><br>
<blockquote class="gmail_quote" style="PADDING-LEFT: 1ex; MARGIN: 0px 0px 0px 0.8ex; BORDER-LEFT: #ccc 1px solid">I dont think it is a visualisation artefact, because it really gives less number of bonds, corresponding to the missing intermonomeric bonds.<br>
</blockquote><br></div>Your visualization software is not reading a file format that is encoding any topology (since you haven't produced one yet). Thus it's inventing some bonds based on the coordinates. The output of pdb2gmx that contains topology information is the .top file, not the coordinate file.
<div class="Ih2E3d"><br><br>
<blockquote class="gmail_quote" style="PADDING-LEFT: 1ex; MARGIN: 0px 0px 0px 0.8ex; BORDER-LEFT: #ccc 1px solid">Should I do now both (introducing a chain identifier in the pdb, but how?, and defining the intermonomeric bonds via dummy in the pdb and tdb?)?<br>
</blockquote><br></div>I suppose there are some circumstances where chain identifiers might be useful, but since you haven't even told us what your monomer and system are, it's hard to help you.
<div>
<div></div>
<div class="Wj3C7c"><br><br>Mark<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>