Molecular recognition and packing frustration in a helical protein

Biomolecular recognition entails attractive forces for the functional native states and discrimination against potential nonnative interactions that favor alternate stable configurations. The challenge posed by the competition of nonnative stabilization against native-centric forces is conceptualized as frustration. Experiment indicates that frustration is often minimal in evolved biological systems although nonnative possibilities are intuitively abundant. Much of the physical basis of minimal frustration in protein folding thus remains to be elucidated. Here we make progress by studying the colicin immunity protein Im9. To assess the energetic favorability of nonnative versus native interactions, we compute free energies of association of various combinations of the four helices in Im9 (referred to as H1, H2, H3, and H4) by extensive explicit-water molecular dynamics simulations (total simulated time > 300 μs), focusing primarily on the pairs with the largest native contact surfaces, H1-H2 and H1-H4. Frustration is detected in H1-H2 packing in that a nonnative packing orientation is significantly stabilized relative to native, whereas such a prominent nonnative effect is not observed for H1-H4 packing. However, in contrast to the favored nonnative H1-H2 packing in isolation, the native H1-H2 packing orientation is stabilized by H3 and loop residues surrounding H4. Taken together, these results showcase the contextual nature of molecular recognition, and suggest further that nonnative effects in H1-H2 packing may be largely avoided by the experimentally inferred Im9 folding transition state with native packing most developed at the H1-H4 rather than the H1-H2 interface.

[1]  S. Radford,et al.  Rapid folding with and without populated intermediates in the homologous four-helix proteins Im7 and Im9. , 1999, Journal of molecular biology.

[2]  Y. Sugita,et al.  Replica-exchange molecular dynamics method for protein folding , 1999 .

[3]  K. Hukushima,et al.  Exchange Monte Carlo Method and Application to Spin Glass Simulations , 1995, cond-mat/9512035.

[4]  T. Darden,et al.  A smooth particle mesh Ewald method , 1995 .

[5]  B. Roux The calculation of the potential of mean force using computer simulations , 1995 .

[6]  Sheena E Radford,et al.  Structural analysis of the rate-limiting transition states in the folding of Im7 and Im9: similarities and differences in the folding of homologous proteins. , 2003, Journal of molecular biology.

[7]  C L Brooks,et al.  Exploring the origins of topological frustration: design of a minimally frustrated model of fragment B of protein A. , 1999, Proceedings of the National Academy of Sciences of the United States of America.

[8]  M. Levitt,et al.  Computer simulation of protein folding , 1975, Nature.

[9]  Martin Zacharias,et al.  Role of tyrosine hot‐spot residues at the interface of colicin E9 and immunity protein 9: A comparative free energy simulation study , 2013, Proteins.

[10]  Gerhard Hummer,et al.  Native contacts determine protein folding mechanisms in atomistic simulations , 2013, Proceedings of the National Academy of Sciences.

[11]  Colin Kleanthous,et al.  Protein-Protein Interaction Specificity of Im9 for the Endonuclease Toxin Colicin E9 Defined by Homologue-scanning Mutagenesis* , 1997, The Journal of Biological Chemistry.

[12]  Cinque S. Soto,et al.  Evaluating conformational free energies: The colony energy and its application to the problem of loop prediction , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[13]  Hue Sun Chan,et al.  Effects of desolvation barriers and sidechains on local-nonlocal coupling and chevron behaviors in coarse-grained models of protein folding. , 2014, Physical chemistry chemical physics : PCCP.

[14]  Jane R. Allison,et al.  Current computer modeling cannot explain why two highly similar sequences fold into different structures. , 2011, Biochemistry.

[15]  Benoît Roux,et al.  Extension to the weighted histogram analysis method: combining umbrella sampling with free energy calculations , 2001 .

[16]  R. Biltonen,et al.  Validity of the “two‐state” hypothesis for conformational transitions of proteins , 1966, Biopolymers.

[17]  David A Agard,et al.  How general is the nucleation–condensation mechanism? , 2008, Proteins.

[18]  W. Kabsch,et al.  Dictionary of protein secondary structure: Pattern recognition of hydrogen‐bonded and geometrical features , 1983, Biopolymers.

[19]  Hue Sun Chan,et al.  Hydrophobic association of alpha-helices, steric dewetting, and enthalpic barriers to protein folding. , 2007, Proceedings of the National Academy of Sciences of the United States of America.

[20]  Chris Neale,et al.  Simulated Tempering Distributed Replica Sampling, Virtual Replica Exchange, and Other Generalized-Ensemble Methods for Conformational Sampling. , 2009, Journal of chemical theory and computation.

[21]  Michele Vendruscolo,et al.  The Mechanism of Folding of Im7 Reveals Competition between Functional and Kinetic Evolutionary Constraints , 2009, Nature Structural &Molecular Biology.

[22]  P. Wolynes,et al.  Spin glasses and the statistical mechanics of protein folding. , 1987, Proceedings of the National Academy of Sciences of the United States of America.

[23]  R. Best Computational and theoretical advances in studies of intrinsically disordered proteins. , 2017, Current opinion in structural biology.

[24]  Hue Sun Chan,et al.  Desolvation is a likely origin of robust enthalpic barriers to protein folding. , 2005, Journal of molecular biology.

[25]  T. Darden,et al.  Particle mesh Ewald: An N⋅log(N) method for Ewald sums in large systems , 1993 .

[26]  D. Baker,et al.  A surprising simplicity to protein folding , 2000, Nature.

[27]  Sheena E. Radford,et al.  Im7 folding mechanism: misfolding on a path to the native state , 2002, Nature Structural Biology.

[28]  P. Faccioli,et al.  Folding Mechanism of Proteins Im7 and Im9: Insight from All-Atom Simulations in Implicit and Explicit Solvent. , 2016, The journal of physical chemistry. B.

[29]  Shankar Kumar,et al.  Multidimensional free‐energy calculations using the weighted histogram analysis method , 1995, J. Comput. Chem..

[30]  H S Chan,et al.  Modeling protein density of states: Additive hydrophobic effects are insufficient for calorimetric two‐state cooperativity , 2000, Proteins.

[31]  Sheena E. Radford,et al.  Ultrarapid mixing experiments reveal that Im7 folds via an on-pathway intermediate , 2001, Nature Structural Biology.

[32]  Sheena E Radford,et al.  Switching two-state to three-state kinetics in the helical protein Im9 via the optimisation of stabilising non-native interactions by design. , 2004, Journal of molecular biology.

[33]  Peter G. Wolynes,et al.  Protein Frustratometer 2: a tool to localize energetic frustration in protein molecules, now with electrostatics , 2016, Nucleic Acids Res..

[34]  Paul Robustelli,et al.  Water dispersion interactions strongly influence simulated structural properties of disordered protein states. , 2015, The journal of physical chemistry. B.

[35]  S. Whittington,et al.  Monte carlo study of the interacting self-avoiding walk model in three dimensions , 1996 .

[36]  Chris Sander,et al.  The double cubic lattice method: Efficient approaches to numerical integration of surface area and volume and to dot surface contouring of molecular assemblies , 1995, J. Comput. Chem..

[37]  Peter G Wolynes,et al.  Consequences of localized frustration for the folding mechanism of the IM7 protein , 2007, Proceedings of the National Academy of Sciences.

[38]  Dengming Ming,et al.  Energetic Frustrations in Protein Folding at Residue Resolution: A Homologous Simulation Study of Im9 Proteins , 2014, PloS one.

[39]  Hue Sun Chan,et al.  Charge pattern matching as a ‘fuzzy’ mode of molecular recognition for the functional phase separations of intrinsically disordered proteins , 2017, 1707.08990.

[40]  Hue Sun Chan,et al.  Solvation and desolvation effects in protein folding: native flexibility, kinetic cooperativity and enthalpic barriers under isostability conditions , 2005, Physical biology.

[41]  Stefano Piana,et al.  Assessing the accuracy of physical models used in protein-folding simulations: quantitative evidence from long molecular dynamics simulations. , 2014, Current opinion in structural biology.

[42]  W. L. Jorgensen,et al.  Comparison of simple potential functions for simulating liquid water , 1983 .

[43]  J. Onuchic,et al.  Topological and energetic factors: what determines the structural details of the transition state ensemble and "en-route" intermediates for protein folding? An investigation for small globular proteins. , 2000, Journal of molecular biology.

[44]  Peter G. Wolynes,et al.  As simple as can be? , 1997, Nature Structural Biology.

[45]  N. Go,et al.  Studies on protein folding, unfolding and fluctuations by computer simulation. I. The effect of specific amino acid sequence represented by specific inter-unit interactions. , 2009 .

[46]  Hue Sun Chan,et al.  Theoretical Insights into the Biophysics of Protein Bi-stability and Evolutionary Switches , 2016, PLoS Comput. Biol..

[47]  William A. Eaton,et al.  Structural origin of slow diffusion in protein folding , 2015, Science.

[48]  Carsten Kutzner,et al.  GROMACS 4:  Algorithms for Highly Efficient, Load-Balanced, and Scalable Molecular Simulation. , 2008, Journal of chemical theory and computation.

[49]  R. Swendsen,et al.  THE weighted histogram analysis method for free‐energy calculations on biomolecules. I. The method , 1992 .

[50]  P. Kollman,et al.  Settle: An analytical version of the SHAKE and RATTLE algorithm for rigid water models , 1992 .

[51]  J. Onuchic,et al.  Protein folding mediated by solvation: Water expulsion and formation of the hydrophobic core occur after the structural collapse , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[52]  Hue Sun Chan,et al.  Pressure-Dependent Properties of Elementary Hydrophobic Interactions: Ramifications for Activation Properties of Protein Folding. , 2014, The journal of physical chemistry. B.

[53]  H. G. Petersen,et al.  Error estimates on averages of correlated data , 1989 .

[54]  R. Friesner,et al.  Evaluation and Reparametrization of the OPLS-AA Force Field for Proteins via Comparison with Accurate Quantum Chemical Calculations on Peptides† , 2001 .

[55]  M. Karplus,et al.  Comparison of the transition state ensembles for folding of Im7 and Im9 determined using all‐atom molecular dynamics simulations with ϕ value restraints , 2003, Proteins.

[56]  Emanuele Paci,et al.  The effect of increasing the stability of non-native interactions on the folding landscape of the bacterial immunity protein Im9. , 2007, Journal of molecular biology.

[57]  S. Radford,et al.  Acidic conditions stabilise intermediates populated during the folding of Im7 and Im9. , 2001, Journal of molecular biology.

[58]  C Kleanthous,et al.  Three-dimensional solution structure and 13C nuclear magnetic resonance assignments of the colicin E9 immunity protein Im9. , 1996, Biochemistry.

[59]  G. Torrie,et al.  Nonphysical sampling distributions in Monte Carlo free-energy estimation: Umbrella sampling , 1977 .

[60]  R. Dror,et al.  How Fast-Folding Proteins Fold , 2011, Science.

[61]  Alan M. Ferrenberg,et al.  New Monte Carlo technique for studying phase transitions. , 1988, Physical review letters.

[62]  H Videler,et al.  A structural comparison of the colicin immunity proteins Im7 and Im9 gives new insights into the molecular determinants of immunity-protein specificity. , 1998, The Biochemical journal.

[63]  Sheena E. Radford,et al.  Conformational Properties of the Unfolded State of Im7 in Nondenaturing Conditions , 2012, Journal of molecular biology.

[64]  Collin M. Stultz,et al.  Perturbing the folding energy landscape of the bacterial immunity protein Im7 by site-specific N-linked glycosylation , 2010, Proceedings of the National Academy of Sciences.

[65]  Berk Hess,et al.  P-LINCS:  A Parallel Linear Constraint Solver for Molecular Simulation. , 2008, Journal of chemical theory and computation.

[66]  W. L. Jorgensen,et al.  Development and Testing of the OPLS All-Atom Force Field on Conformational Energetics and Properties of Organic Liquids , 1996 .

[67]  Wookyung Yu,et al.  Benchmarking all-atom simulations using hydrogen exchange , 2014, Proceedings of the National Academy of Sciences.

[68]  H. Berendsen,et al.  A LEAP-FROG ALGORITHM FOR STOCHASTIC DYNAMICS , 1988 .

[69]  Peter G Wolynes,et al.  Frustration in biomolecules , 2013, Quarterly Reviews of Biophysics.

[70]  R. Best,et al.  Balanced Protein–Water Interactions Improve Properties of Disordered Proteins and Non-Specific Protein Association , 2014, Journal of chemical theory and computation.

[71]  U. Hansmann Parallel tempering algorithm for conformational studies of biological molecules , 1997, physics/9710041.

[72]  Joan-Emma Shea,et al.  Simulations of disordered proteins and systems with conformational heterogeneity. , 2017, Current opinion in structural biology.

[73]  Justin L. MacCallum,et al.  Hydrophobic association of α-helices, steric dewetting, and enthalpic barriers to protein folding , 2007, Proceedings of the National Academy of Sciences.

[74]  J. Onuchic,et al.  Water mediation in protein folding and molecular recognition. , 2006, Annual review of biophysics and biomolecular structure.

[75]  N. Go Theoretical studies of protein folding. , 1983, Annual review of biophysics and bioengineering.

[76]  Monika Fuxreiter,et al.  Fuzzy complexes: a more stochastic view of protein function. , 2012, Advances in experimental medicine and biology.

[77]  W F Drew Bennett,et al.  Statistical Convergence of Equilibrium Properties in Simulations of Molecular Solutes Embedded in Lipid Bilayers. , 2011, Journal of chemical theory and computation.

[78]  Régis Pomès,et al.  Accelerating Convergence in Molecular Dynamics Simulations of Solutes in Lipid Membranes by Conducting a Random Walk along the Bilayer Normal. , 2013, Journal of chemical theory and computation.

[79]  Hue Sun Chan,et al.  Cooperativity, local-nonlocal coupling, and nonnative interactions: principles of protein folding from coarse-grained models. , 2011, Annual review of physical chemistry.

[80]  Tao Chen,et al.  Native Contact Density and Nonnative Hydrophobic Effects in the Folding of Bacterial Immunity Proteins , 2015, PLoS Comput. Biol..

[81]  H. Berendsen,et al.  Molecular dynamics with coupling to an external bath , 1984 .

[82]  Wenwei Zheng,et al.  Empirical Optimization of Interactions between Proteins and Chemical Denaturants in Molecular Simulations. , 2015, Journal of chemical theory and computation.

[83]  Hue Sun Chan,et al.  Competition between native topology and nonnative interactions in simple and complex folding kinetics of natural and designed proteins , 2010, Proceedings of the National Academy of Sciences.

[84]  Shachi Gosavi,et al.  Using the folding landscapes of proteins to understand protein function. , 2016, Current opinion in structural biology.

[85]  H. Chan,et al.  Theoretical perspectives on nonnative interactions and intrinsic disorder in protein folding and binding. , 2015, Current opinion in structural biology.

[86]  Hue Sun Chan,et al.  Volumetric Physics of Polypeptide Coil-Helix Transitions. , 2016, Biochemistry.