A Fast and Precise Approach for Computational Saturation Mutagenesis and its Experimental Validation by Using an Artificial (βα)8‐Barrel Protein

We present a computational saturation mutagenesis protocol (CoSM) that predicts the impact on stability of all possible amino acid substitutions for a given site at an internal protein interface. CoSM is an efficient algorithm that uses a combination of rotamer libraries, side‐chain flips, energy minimization, and molecular dynamics equilibration. Because CoSM considers full side‐chain and backbone flexibility in the local environment of the mutated position, amino acids larger than the wild‐type residue are also modeled in a proper manner. To assess the performance of CoSM, the effect of point mutations on the stability of an artificial (βα)8‐barrel protein that has been designed from identical (βα)4‐half barrels, was studied. In this protein, position 234N is a previously identified stability hot‐spot that is located at the interface of the two half barrels. By using CoSM, changes in protein stability were predicted for all possible single point mutations replacing wild‐type Val234N. In parallel, the stabilities of 14 representative mutants covering all amino acid classes were experimentally determined. A linear correlation of computationally and experimentally determined energy values yielded an R2 value of 0.90, which is statistically significant. This degree of coherence is stronger than the ones we obtained for established computational methods of mutational analysis.

[1]  R. Wierenga,et al.  The TIM‐barrel fold: a versatile framework for efficient enzymes , 2001, FEBS letters.

[2]  Feng Ding,et al.  Modeling backbone flexibility improves protein stability estimation. , 2007, Structure.

[3]  Haruki Nakamura,et al.  Announcing the worldwide Protein Data Bank , 2003, Nature Structural Biology.

[4]  N. Pokala,et al.  Energy functions for protein design: adjustment with protein-protein complex affinities, models for the unfolded state, and negative design of solubility and specificity. , 2005, Journal of molecular biology.

[5]  Yuval Inbar,et al.  Protein structure modelling and evaluation based on a 4-distance description of side-chain interactions , 2010, BMC Bioinformatics.

[6]  T. Blundell,et al.  Comparative protein modelling by satisfaction of spatial restraints. , 1993, Journal of molecular biology.

[7]  Shuangye Yin,et al.  Eris: an automated estimator of protein stability , 2007, Nature Methods.

[8]  D. W. Bolen,et al.  Unfolding free energy changes determined by the linear extrapolation method. 1. Unfolding of phenylmethanesulfonyl alpha-chymotrypsin using different denaturants. , 1988, Biochemistry.

[9]  Birte Höcker,et al.  High-resolution crystal structure of an artificial (betaalpha)(8)-barrel protein designed from identical half-barrels. , 2009, Biochemistry.

[10]  Johannes Söding,et al.  The HHpred interactive server for protein homology detection and structure prediction , 2005, Nucleic Acids Res..

[11]  T. Tsukihara,et al.  Conformational changes in the α-subunit coupled to binding of the β2-Subunit of tryptophan synthase from Escherichia coli: Crystal structure of the tryptophan synthase α-subunit alone , 2005 .

[12]  Paul R. Gerber,et al.  MAB, a generally applicable molecular force field for structure modelling in medicinal chemistry , 1995, J. Comput. Aided Mol. Des..

[13]  David Baker,et al.  Protein Structure Prediction Using Rosetta , 2004, Numerical Computer Methods, Part D.

[14]  A. Fersht,et al.  Dissection of the structure and activity of the tyrosyl-tRNA synthetase by site-directed mutagenesis. , 1987, Biochemistry.

[15]  B. Höcker,et al.  Mimicking enzyme evolution by generating new (betaalpha)8-barrels from (betaalpha)4-half-barrels. , 2004, Proceedings of the National Academy of Sciences of the United States of America.

[16]  Marco Bocola,et al.  Stabilisation of a (βα) 8 -Barrel Protein Designed from Identical Half Barrels , 2007 .

[17]  Piero Fariselli,et al.  I-Mutant2.0: predicting stability changes upon mutation from the protein sequence or structure , 2005, Nucleic Acids Res..

[18]  J. Kirsch,et al.  Energetic analysis of an antigen/antibody interface: Alanine scanning mutagenesis and double mutant cycles on the hyhel‐10/lysozyme interaction , 1999, Protein science : a publication of the Protein Society.

[19]  Birte Höcker,et al.  Catalytic versatility, stability, and evolution of the (betaalpha)8-barrel enzyme fold. , 2005, Chemical reviews.

[20]  T Tsujita,et al.  Dependence of conformational stability on hydrophobicity of the amino acid residue in a series of variant proteins substituted at a unique position of tryptophan synthase alpha subunit. , 1987, Proceedings of the National Academy of Sciences of the United States of America.

[21]  P. Kollman,et al.  Calculating structures and free energies of complex molecules: combining molecular mechanics and continuum models. , 2000, Accounts of chemical research.

[22]  Lei Wang,et al.  Accuracy of structure-derived properties in simple comparative models of protein structures , 2005, Nucleic acids research.

[23]  Wilfred F. van Gunsteren,et al.  Thermodynamic cycle integration by computer simulation as a tool for obtaining free energy differences in molecular chemistry , 1987, J. Comput. Aided Mol. Des..

[24]  N. Guex,et al.  SWISS‐MODEL and the Swiss‐Pdb Viewer: An environment for comparative protein modeling , 1997, Electrophoresis.

[25]  G. Schreiber,et al.  Assessing computational methods for predicting protein stability upon mutation: good on average but not in the details. , 2009, Protein engineering, design & selection : PEDS.

[26]  B. L. de Groot,et al.  Predicting free energy changes using structural ensembles. , 2009, Nature methods.

[27]  Paul R. Gerber,et al.  Charge distribution from a simple molecular orbital type calculation and non-bonding interaction terms in the force field MAB , 1998, J. Comput. Aided Mol. Des..

[28]  M Wilmanns,et al.  Structural evidence for evolution of the beta/alpha barrel scaffold by gene duplication and fusion. , 2000, Science.

[29]  Birte Höcker,et al.  Dissection of a (βα)8-barrel enzyme into two folded halves , 2001, Nature Structural Biology.

[30]  Michael A. Johnston,et al.  Integrated prediction of the effect of mutations on multiple protein characteristics , 2011, Proteins.

[31]  B. Matthews,et al.  Structural and genetic analysis of protein stability. , 1993, Annual review of biochemistry.

[32]  Jan H. Jensen,et al.  Very fast empirical prediction and rationalization of protein pKa values , 2005, Proteins.

[33]  L. Serrano,et al.  Predicting changes in the stability of proteins and protein complexes: a study of more than 1000 mutations. , 2002, Journal of molecular biology.

[34]  Wilfred F van Gunsteren,et al.  Biomolecular modeling: Goals, problems, perspectives. , 2006, Angewandte Chemie.

[35]  W. F. V. Gunsteren,et al.  Biomolekulare Modellierung: Ziele, Probleme, Perspektiven , 2006 .

[36]  Adrian A Canutescu,et al.  Access the most recent version at doi: 10.1110/ps.03154503 References , 2003 .