Protein design automation

Abstract

We have conceived and implemented a cyclical protein design strategy that couples theory, computation, and experimental testing. The combinatorially large number of possible sequences and the incomplete understanding of the factors that control protein structure are the primary obstacles in protein design. Our protein design automation algorithm objectively predicts protein sequences likely to achieve a desired fold. Using a rotamer description of the side chains, we implemented a fast discrete search algorithm based on the Dead-End Elimination Theorem to rapidly find the globally optimal sequence in its optimal geometry from the vast number of possible solutions. Rotamer sequences were scored for steric complementarity using a van der Waals potential. A Monte Carlo search was then executed, starting at the optimal sequence, in order to find other high-scoring sequences. As a test of the design methodology, high-scoring sequences were found for the buried hydrophobic residues of a homodimeric coiled coil based on GCN4-p1. The corresponding peptides were synthesized and characterized by CD spectroscopy and size-exclusion chromatography. All peptides were dimeric and nearly 100% helical at 1 °C, with melting temperatures ranging from 24 °C to 57 °C. A quantitative structure activity relation analysis was performed on the designed peptides, and a significant correlation was found with surface area burial. Incorporation of a buried surface area potential in the scoring of sequences greatly improved the correlation between predicted and measured stabilities and demonstrated experimental feedback in a complete design cycle.