[gmx-users] How to restrain the center of gravity
kazunori.hirayama at astellas.com
Tue Nov 10 11:02:53 CET 2015
Hello, this is the first post for me. My apology for any rudeness.
My molecule of interest is a protein of about 12kD.
When I elevated temperature from 400K to 500K with .mdp file according to a tutorial, it seems that the proteins are shaking or parallel shifting.
Is it necessary for us to restrain the center of gravity? If so, could you show us which parameter to modify in .mdp file.
My parameters relevant to restraint are shown below:
define = -DPOSRES ; position restrain the protein
constraint_algorithm = lincs ; holonomic constraints
constraints = all-bonds ; all bonds (even heavy atom-H bonds) constrained
lincs_iter = 1 ; accuracy of LINCS
lincs_order = 4 ; also related to accuracy
Thank you in advance.
Bioscience Research Labs. Drug Discovery Research
Astellas Pharma Inc.
Tel: +81-29-863-7058 E-mail: kazunori.hirayama at astellas.com<mailto:kazunori.hirayama at astellas.com>
More information about the gromacs.org_gmx-users