On 08/10/17 04:12, atanu das wrote:
Hi,
I am intending to run a simulation of a protein in 100 nM NaCl solution by
turning off protein's electrostatic interactions, but keeping the LJ alive, so
that water and salt behave as usual and the protein only behaves as a
hydrophobic polymer chain. The most obvious way that I could think of is
editing the topol.top file and setting the partial charges of all the protein
atoms to be zero. However, I am slightly skeptical about using such a crude
method. Is there any better way to turn off the electrostatics of the protein
only?
ThanksAtanu
You may want to slowly turn of the charges using free energy code
although it will be slower. Of course even a hydrophobic chain has
charges, and doing what you plan to will mess up the balance between
Coulomb and Lennard-Jones interaction. You may get some insight but
nothing of physical relevance that you can compare to any experiment.
--
David van der Spoel, Ph.D., Professor of Biology
Head of Department, Cell & Molecular Biology, Uppsala University.
Box 596, SE-75124 Uppsala, Sweden. Phone: +46184714205.
http://www.icm.uu.se
--
Gromacs Users mailing list
* Please search the archive at
http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!
* Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a
mail to gmx-users-requ...@gromacs.org.
