Dear Eudes,

You need to technically make it 1 "big molecule". As long as there are no bonds specified between molecule A and B. The big molecule will behave as if there are two molecules in it. I do it to dock small molecules to proteins.

Steps running Gromacs 3.3.3:

1) Combine molecules A and B into a big structure file (i.e. pdb or gro).
2) Edit this file, so that they are the same chain.
3) run pdb2gmx
4) Change all the residue names to MOL and change Protein_A to MOL in the top and pdb files. 5) Make sure that there is no bonds between molecule A and B in the top file. There probably won't be, but it is good to check anyway.
5) run grompp to check this file to see if it is acceptable to gromacs.
6) add distance restraints after the improper dihedrals. If you put it anywhere else, it will not work. 7) run grompp again to check this file to see if it is acceptable to gromacs

Other advise:
1) I would roughly dock the molecules together, so that you don't have to have such a large water box. 2) I would experiment with force constants and time constants for the distance restraints. A useful article is:

Torda, A. E., Scheek, R. M., and Gunsteren, W. F. v. (1989) Time- dependent distance restraints in molecular dynamics simulations, Chem. Phys. Lett. 157, 289-294.

I hope that helps.


On Nov 15, 2008, at 3:58 AM, Eudes Fileti wrote:

Dear gmx-users, I need to use constraints between two
polyatomic molecules in order to calculate PMF using free energy code.
Can anybody help me to use these constraints over two different molecules?
I ve already read something in the forum but not helped me enough.
(see http://www.gromacs.org/pipermail/gmx-users/2006-February/020049.html) .
bests
eef
_______________________________________
Eudes Eterno Fileti
Centro de Ciências Naturais e Humanas
Universidade Federal do ABC
Rua Santa Adélia, 166 - Bloco B, Sala 1048
09210-170  Santo André - SP Brasil
+55.11.4437-8408
skype: eefileti
http://cromo.ufabc.edu.br/~fileti/
_______________________________________________
gmx-users mailing list    [email protected]
http://www.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the
www interface or send it to [EMAIL PROTECTED]
Can't post? Read http://www.gromacs.org/mailing_lists/users.php

_______________________________________________
gmx-users mailing list    [email protected]
http://www.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to [EMAIL PROTECTED]
Can't post? Read http://www.gromacs.org/mailing_lists/users.php

Reply via email to