Hello Everyone, I was willing to conduct umbrella sampling on the dihedrals derived from the sidechains of a ion channel. So i tried dihedral restraint to achieve the goal. However I encountered an awkward problem, hope someone can help me with that. The dihedrals of the initial conformation I used in my MD is -60 and 153 degree, respectively for phi and psi. And then I restricted both of them to 180 degree and conducted the dihedral restraint MD for 2ns. However the dihedral analysis of my final result was quite awkward, since they kept on fluctuating just around their initial value during this 2ns MD process. At first, I considered maybe it was caused by the small restraint force constant, which was set as 100 in my .mdp file. Nevertheless, after I changed it to 500, no better result was given out. I am wondering whether the energy barrier is too high for them to cross, so in the 2ns short MD process, they have to move around their initial configuration. If this is the key to my problem, how can I solve this problem? Many thanks in advance.
-- View this message in context: http://gromacs.5086.n6.nabble.com/dihedral-restraint-MD-tp4998407.html Sent from the GROMACS Users Forum mailing list archive at Nabble.com. -- 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