[gmx-users] Help with position restrains for heterodimer system.

Xavier Periole x.periole at rug.nl
Tue May 24 15:04:14 CEST 2005


itai bloch wrote:

>Hi everyone,
>I am trying to simulate a bilayer system with a heterodimer inside it.
>
>I have a problem with the position restrained stage:
>
>For some reason, only one of the two peptide chains remains fixed.
>I have noticed that the number of POSRES in the grompp output is
>exactly twice the number of my restrained heavy atoms (the total of
>both peptide chains).
>I have searched the GMX mailing list and tried the suggested solutions
>(i.e. ordering the molecules differently and several different
>'#include' locations), but it failed.
>
>Can anyone help me with it???
>
I checked something for you. I had a dimer in DPPC same stuff as you.
By using the shuffle option in grompp I completely lost the position 
restrains
definition that was ok without it ... the first atom (using shuffle) was 
a DPPC
atom.

might be something top fix in next version.


-- 
----------------------------------
 Xavier Periole - Ph.D.

 Dept. of Biophysical Chemistry / MD Group   
 Univ. of Groningen
 Nijenborgh 4
 9747 AG Groningen
 The Netherlands
 
 Tel: +31-503634329
 Fax: +31-503634800
 email: x.periole at rug.nl
 web-page: http://md.chem.rug.nl/~periole
----------------------------------




More information about the gromacs.org_gmx-users mailing list