On 15 Apr 2010, at 18:06, Trang wrote:
My target system is a protein with lipid molecules added randomly (using GENBOX). Running MD, I expect to

I hope you're using a larger van der Waals distance (0.24nm or so) when inserting the lipids.

broke down the problem, that is, to run md simulation for the protein molecule only, and in vacuum. Still no improvement. Although all the distances in the minimized structure are visually proper, the system exploded.

If you can't run the protein even in vacuum, then the problem is either in your MD parameters, starting co-ordinates or some simple mistake somewhere. Since the .mdp files you posted seem ok, my bet is on the starting co-ordinates, just as Xavier proposed.

Could you post the system topology (.top) and the protein topology (.itp), so we can rule out any mistakes in those? You should also double check whether you have close contacts between some atoms in the protein; e.g. in VMD this is easily done using the "dynamic bonds" representation.

-martti-
--
Post-doctoral research fellow
Moleculaire Dynamica
University of Groningen
Nijenborgh 4, 9747AG Groningen, the Netherlands
tel. +(31) 50 363 4339 | fax. +(31) 50 363 4398

--
gmx-users mailing list    gmx-users@gromacs.org
http://lists.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 gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/mailing_lists/users.php

Reply via email to