<html>
<head>
<style>
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
FONT-SIZE: 10pt;
FONT-FAMILY:Tahoma
}
</style>
</head>
<body class='hmmessage'>
<BR>Hi Luciano,<BR>
use the -label option to separately label the gro files (e.g A for the first peptide and B for the 2nd), before you cat them with an editor. pdb2gmx should now recognise them as separate molecules.<BR>
Rgds<BR>
John<BR>
<BR>&gt; Date: Sat, 8 Sep 2007 15:43:33 +0200<BR>&gt; Subject: Re: [gmx-users] Adding two peptides in a water box<BR>&gt; From: Lars.Schaefer@mpi-bpc.mpg.de<BR>&gt; To: gmx-users@gromacs.org<BR>&gt; <BR>&gt; &gt; Hi Everybory,<BR>&gt; &gt;<BR>&gt; &gt; I would like to simulate the interaction between two separated peptides in<BR>&gt; &gt; a water box. Every peptide needs to be N-and C-terminated. I have tried<BR>&gt; &gt; different alternatives in order to create the corresponding GROMACS files<BR>&gt; &gt; to start my simulation. I have used the GROMACS porograms tools as<BR>&gt; &gt; described in the manual. However, I still have some problems and need of<BR>&gt; &gt; your expeirences.<BR>&gt; &gt;<BR>&gt; &gt; This is what I do:<BR>&gt; &gt;<BR>&gt; &gt; 1.- I start from a pdb of one peptide.<BR>&gt; &gt; 2.- Create the gro and top files. In this step I terminate the peptide<BR>&gt; &gt; without any problem.<BR>&gt; &gt; 3.- Apply corresponding translations and rotation to the existing gro file<BR>&gt; &gt; (obtained in step 1) and create a new .gro file. This one is termed the<BR>&gt; &gt; second peptide.<BR>&gt; &gt; 4.- Cat them with an editor and save the result in a file named dimer.gro.<BR>&gt; &gt; I can visualize it without any problem. So far, everything is going good.<BR>&gt; &gt; 5.- HERE STARTS THE PROBLEM. I use dimer.gro as input to pdb2gmx to create<BR>&gt; &gt; the new, dimer.top file. The pdb2gmx removes the hydrogens of the N-term<BR>&gt; &gt; of one petide and the oxigen of the C-term in the second peptide.<BR>&gt; &gt; Although, they are far apart.<BR>&gt; &gt;<BR>&gt; &gt; Does anyone know how to overcome this porblem?<BR>&gt; <BR>&gt; Hi,<BR>&gt; why do you run pdb2gmx for the dimer? run it for the monomer, and then<BR>&gt; just put a "Protein 2" instead of a "Protein 1" (or however you name it)<BR>&gt; at the end of your topol.top file (just as is done for the water<BR>&gt; molecules, ions, etc).<BR>&gt; Lars<BR>&gt; &gt;<BR>&gt; <BR>&gt; &gt; Thank in advance for you time.<BR>&gt; &gt;<BR>&gt; &gt; Ciao,<BR>&gt; &gt;<BR>&gt; &gt; Luciano<BR>&gt; &gt;<BR>&gt; &gt; Dr. Luciano Triguero College of Art and Science Department of Physics and<BR>&gt; &gt; Chemistry Cox Science Building 1301 Memorial Drive, Room 146 P.O Box<BR>&gt; &gt; 249118 Coral Gables, FL 33124-0431 Cellular: 305-904-2419 Office:<BR>&gt; &gt; 305-284-3938<BR>&gt; &gt; _______________________________________________<BR>&gt; &gt; gmx-users mailing list gmx-users@gromacs.org<BR>&gt; &gt; http://www.gromacs.org/mailman/listinfo/gmx-users<BR>&gt; &gt; Please search the archive at http://www.gromacs.org/search before posting!<BR>&gt; &gt; Please don't post (un)subscribe requests to the list. Use the<BR>&gt; &gt; www interface or send it to gmx-users-request@gromacs.org.<BR>&gt; &gt; Can't post? Read http://www.gromacs.org/mailing_lists/users.php<BR>&gt; &gt;<BR>&gt; <BR>&gt; <BR>&gt; _______________________________________________<BR>&gt; gmx-users mailing list gmx-users@gromacs.org<BR>&gt; http://www.gromacs.org/mailman/listinfo/gmx-users<BR>&gt; Please search the archive at http://www.gromacs.org/search before posting!<BR>&gt; Please don't post (un)subscribe requests to the list. Use the <BR>&gt; www interface or send it to gmx-users-request@gromacs.org.<BR>&gt; Can't post? Read http://www.gromacs.org/mailing_lists/users.php<BR><BR><br /><hr />Make your little one a shining star! <a href='http://www.reallivemoms.com?ocid=TXT_TAGHM&loc=us' target='_new'>Shine on!</a></body>
</html>