Hi, I'm trying to do some surface area calculations on MARTINI models using APBS, which requires supplying vdw radii via a pqr file. Given how MARTINI treats LJ interactions (0.47 nm radius for all particle types), how do I calculate reasonable effective vdw radii? Obviously the PC headgroup (NC3) will have a larger vdw radius than the PE headgroup (NH3), and g_sas seems to account for this (reported ASA for -NC3 group is always higher than for -NH3 group) but I can't figure out how. Indeed, using "editconf -mead" to produce a pqr file from a MARTINI structure shows 2.35 A as the vdw radius for all lipid groups, which I fear will cause inaccurate calculations in APBS. Has anyone had this problem before?
Mike -- ==================================== Michael D. Daily Postdoctoral research associate Pacific Northwest National Lab (PNNL) 509-375-4581 (formerly Qiang Cui group, University of Wisconsin-Madison)
-- 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