Modified genetic algorithm resolves ambiguous NOE restraints and reduces unsightly NOE violations.

Abstract

In an ideal world, every NOE cross peak would have a unique assignment. However, the interpretation of NOE peaks is frequently complicated by overlapping resonances. In theory, ambiguous assignments could be resolved by performing separate structure calculations with each possible interpretation. Unfortunately, this would require an astronomical amount of computing time. A modified genetic algorithm has been developed that efficiently resolves hundreds of ambiguous restraints in parallel. Each NOE assignment becomes a gene that can be passed on to a new generation. New individuals are constructed by making a constraint lists from a subset of the genes. The constraint lists are then tested for self-consistency by using molecular dynamics to generate new structures for each list. To a first-degree approximation, there is enough information retained in each list to determine the global fold of the protein. Self-consistent constraint lists receive higher scores and their genes (or NOEs) stand a better chance of surviving into the next generation. The process selects NOEs that are consistent with the global fold. Under normal conditions, the program converges in 3 to 8 generations using 70 structures per generation. The final constraints are self-consistent and contain almost no residual NOE violations.