nikhil damle wrote:
Hi All,
I am trying to simulate kinase domain (chain A) in complex with cyclin
(chain B). During solvent dynamics, PBC on and PME treatment of
electrostatics, i need to keep only kinase domain flexible but cyclin fixed.
i put POSRES_B in define variable and tried to simulated the entire
solvated system in GROMACS 4 package. But time and again it failed with
LINCS warnings. On the contrary, if entire system is kept either
flexible or position restrained, run is obtained without any errors. I
tried 1ns MD. Unable to understand possible reasons.
Probably by restraining whatever configuration you're starting with, you're
preserving some unresolved clashes in the system. By either restraining both
chains or allowing both to move, you are overcoming the problem.
Could anybody help in this regard ?
Watch the trajectory to see what's going on. A few other general tips:
http://www.gromacs.org/Documentation/Terminology/Blowing_Up#Diagnosing_an_Unstable_System
-Justin
Regards,
Nikhil
--
========================================
Justin A. Lemkul
Ph.D. Candidate
ICTAS Doctoral Scholar
MILES-IGERT Trainee
Department of Biochemistry
Virginia Tech
Blacksburg, VA
jalemkul[at]vt.edu | (540) 231-9080
http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin
========================================
--
gmx-users mailing list [email protected]
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 [email protected].
Can't post? Read http://www.gromacs.org/mailing_lists/users.php
