I was looking at the g_membed paper recently, and noticed the authors "only" spent 10ns to re-equilibrate the membedded-system after g_membeding. This they also did after bumping the temperature from 300 to 323K "to prevent ordering of the bilayer" even though the "typical" simulation temperatures of POPC is 300 to 310K...Anyone know why 323K was used? Was it to compare it to the DOPC environment or something?
I noticed for my particular system, I had to run a round of EM on the post-membedded system in order to resolve some clashes before I could mdrun. Then I ran 10ns in NPT with gen_vel=yes at 310K and position restraints on the helix CA atoms. I am inferring that because the membrane patch I used had been previously equilibrated after I constructed it and ran it under NPT at 300K for 100ns, 10ns should be sufficient time to re-equilibrate with the protein in the middle and 310K temperature? After the 10ns with protein restrained, I ran Grid_MAT on it and got a reasonable APL (61-62A^2/lipid when taking protein atoms into account, comparable to pre-g_membed patch) and the box dimensions look stable. Is this sufficient to answer the abvove question? Now I am going to be introducing a ligand to the protein. Is there a way to preserve any state, like velocities of all the previous atoms? I was thinking of the doing the EM while freezing everything but the ligand atoms and allowing the ligand to change conformation (like a docking refinement). After this, can I have the thermostat and barostat rescale the original velocities (and gradually heat the ligand) instead of reinitialising everything? -- ================================================================== Peter C. Lai | University of Alabama-Birmingham Programmer/Analyst | BEC 257 Genetics, Div. of Research | 1150 10th Avenue South p...@uab.edu | Birmingham AL 35294-4461 (205) 690-0808 | ================================================================== -- 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/Support/Mailing_Lists/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/Support/Mailing_Lists