I am using the Amber force field in Gromacs on a compound that I do not have a topology for (I do NOT have a copy of AMBER itself). There are no metals, simply O, C, N and H (it is a derivative of glucose). I can derive the partial charges quite easily using R.E.D. tools, but I am unsure of how to derive the bond and angle equilibrium data to go into my .itp file.
The parmchk2 in the AMBER toolkit just takes a guess, but I don't have access to this tool. Would this have been the equivalent of taking a "guess" from the ffbonded.itp file to construct the structure in the .itp?
Or ought I use Gaussian to derive distances and angles and then also fit the potential energy surface of displaced atoms to derive the force constants?