In their 2003 review on force fields protein simulations (Adv. in Protein Chem. 2003, 66, 27) Ponder and Case concluded that “An increase in computer power of at least two orders of magnitude should occur over the next decade. Without further research into the accuracy of force field potentials, future macromolecular modeling may well be limited more by validity of the energy functions, particularly electrostatic terms, than by technical ability to perform the computations.”
We aim to address this issue by systematically scrutinizing the way that well-known force fields (AMBER, CHARMMm, GROMOS,…) are constructed. This is of course a long-term goal but important steps have already been taken. Although point charges are popular they are inherently limited and best replaced by multipole moments. In the study of crystal polymorphism for example multipole moments provide a much needed accurate description of the electrostatic interaction.
The topological approach divides molecules into atoms in a natural and unbiased way, following the construction of the gradient vector field. As such it generates unique atoms, carved out of their molecular environment. The first three-dimensional pictures of topological atoms were produced in 1994 , in a paper that stressed that the topological approach is a departure from the conventional way of defining molecular fragments in the Hilbert space of the basis functions. Because each topological atom is unique there are millions of atoms but transferability enables them to be grouped in atom types by similarity. We have generated a data set of all atoms occurring in one conformation of each natural amino acid and smaller derived molecules. Cluster analysis revealed atom types for carbon , hydrogen, nitrogen, oxygen and sulfur . These atom types are constructed in the opposite way to the construction of atom types in force fields. In force field design one starts from a single atom type per atomic number and gradually introduces new types according to the need for accuracy from the force field. We spent much attention to the question how much the integration error influences atomic properties . To our knowledge this is the first time that atom types are being computed rather than postulated. They may guide the design of future force fields .
To be updated
To be updated