Maybe, one could also think about calculating the RESP charges via QM,
e.g. with Gaussian...Further parameters, like force-constants and
dihedrals could also be evaluated...
But what's against an expensive random number generator?
I fear I'm doing this with FEP all day long...;)
Raja: We are (mostly) not doing LEGO here. You should inform yourself
about the theoretic parts in general. Then you would see, that LEGOing a
molecule can't be feasible....I suggest e.g. Leach, Molecular Modelling,
Section 3: Empirical force field models.
There you can find lots of information about what MD in general resembles.
Best regards
Maik Goette, Dipl. Biol.
Max Planck Institute for Biophysical Chemistry
Theoretical & computational biophysics department
Am Fassberg 11
37077 Goettingen
Germany
Tel. : ++49 551 201 2310
Fax : ++49 551 201 2302
Email : mgoette[at]mpi-bpc.mpg.de
mgoette2[at]gwdg.de
WWW : http://www.mpibpc.gwdg.de/groups/grubmueller/
Mark Abraham wrote:
raja wrote:
Dear gmxions,
I am writing a topology for a ligand using AMBER99 force field
convention. The ligand contains the hydroxamate group, for your quick
reference I have drawn it below
O
HO ||
\ / \
N
Your N has an unsatisfied valence... what's missing? (Don't bother to
tell me, it doesn't matter)
Please let me know which atoms type in amber can be used for these
atoms.... Especially the Nitrogen atom in this group. Since I am using
carbonyl group and hydroxyl group atoms types for oxygen (But I am not
sure of its validity). Please help me to chose appropriate atoms type.
Probably, none of them. Force fields are not magic numbers that can be
applied to all possible "types" of atoms. They are constructs designed
to be reliably applicable to fairly narrow sets of atom types under the
range of conditions they were parametrized with (also normally narrow).
You will see this when you read the paper that describes the development
process of this force field. Notwithstanding, people routinely abuse
them with gay abandon, probably on the bad advice of others, or their
own experience of the fairly good transferability of quantum-chemical
methods.
Your only vaguely reliable course is to find a set of experimental data
consistent with those used for the development of the force field, work
out what new values you need to derive and to follow a process to
optimize those values to reproduce the experimental data. This process
is not for the inexperienced, or for the faint of heart.
Otherwise, you'll have the usual expensive random number generator.
Mark
_______________________________________________
gmx-users mailing list gmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please don't post (un)subscribe requests to the list. Use the www
interface or send it to [EMAIL PROTECTED]
Can't post? Read http://www.gromacs.org/mailing_lists/users.php
.
_______________________________________________
gmx-users mailing list gmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please don't post (un)subscribe requests to the list. Use the
www interface or send it to [EMAIL PROTECTED]
Can't post? Read http://www.gromacs.org/mailing_lists/users.php