On 2011-05-27 10:50:14AM -0500, simon sham wrote: > Hi, > I have recently done two simulations on a protein at high temperature near > its melting temperature. At the beginning I used CHARMM forcefield, and then > OPLSAA to double check the results. There are some differences in the > structures between the forcefield used. I know different forcefields can > give different results. All parameters in the simulations were the same > except the forcefield. Is there anyway I can tell which forcefield gives more > reliable results? >
Also you never stated what RMSD/RMSFs between the two and if they are significant to what you care to see (i.e. if you expect, say, alpha helix to random coil transition, what changes in the unfolding do you see that make you suspect one forcefield is better than the other?). Also if this is in a solvated system, the rest of the forcefield has been parameterized with the implicit assumption that TIPS3P (Charmm-specific) water will be used so make sure you have switched your water type when using CHARMM. -- ================================================================== 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