Chih-Ying Lin wrote:
Hi
The system = one protein molecule + 10 ligands + water molecules in 6
x 6 x 6 box.
after 100 ns simulatin, 10 ligands have attached on the protein.
Now, I tried to put another 10 ligands into the system.
The steps are as follows.
1. Take one coordinate file and remove all the coordinates of the
water molecules.
2. Center the protein with the attached 10 ligands
3. With "genbox" command, to put another 10 ligands into the simulation box.
4. Solvate the system, with "genbox" command
5. Add ions, with "genion" command
6. Energy minimization of the solvated system
pbc = xyz (minim.mdp)
=> the system is crystallized with visualization
What's the energy and maximum force?
7. Relaxation of solvent and hydrogen atom positions
Run => Position restrained MD
=> simulation break
Check the box size at each step makes sense. Load the structure at each
step into a visualization program and check that it makes sense with the
box. I'm expecting you've got some molecule clashing with another. If
everything's OK, trying a smaller integration time step can sometimes help.
Otherwise, the procedure looks good.
Mark
--
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