Coarse-grained model for colloidal protein interactions, B(22), and protein cluster formation.

Reversible protein cluster formation is an important initial step in the processes of native and non-native protein aggregation, but involves relatively long time and length scales for detailed atomistic simulations and extensive mapping of free energy landscapes. A coarse-grained (CG) model is presented to semiquantitatively characterize the thermodynamics and key configurations involved in the landscape for protein oligomerization, as well as experimental measures of interactions such as the osmotic second virial coefficient (B22). Based on earlier work (Grüenberger et al., J. Phys. Chem. B 2013, 117, 763), this CG model treats proteins as rigid bodies composed of one bead per amino acid, with each amino acid having specific parameters for its size, hydrophobicity, and charge. The net interactions are a combination of steric repulsions, short-range attractions, and screened long-range charge-charge interactions. Model parametrization was done by fitting simulation results against experimental value of B22 as a function of solution ionic strength for α-chymotrypsinogen A and γD-Crystallin (gD-Crys). The CG model is applied to characterize the pairwise interactions and dimerization of gD-Crys and the dependence on temperature, protein concentration, and ionic strength. The results illustrate that at experimentally relevant conditions where stable dimers do not form, the entropic contributions are predominant in the free-energy of protein cluster formation and colloidal protein interactions, arguing against interpretations that treat B22 primarily from energetic considerations alone. Additionally, the results suggest that electrostatic interactions help to modulate the population of the different stable configurations for protein nearest-neighbor pairs, while short-range attractions determine the relative orientations of proteins within these configurations. Finally, simulation results are combined with Principal Component Analysis to identify those amino-acids/surface patches that form interprotein contacts at conditions that favor dimerization of gD-Crys. The resulting regions agree with previously found aggregation-prone sites, as well as suggesting new ones that may be important.

[1]  Charles F. Zukoski,et al.  Protein interactions and crystallization , 1996 .

[2]  N. Mousseau,et al.  Structures and thermodynamics of Alzheimer's amyloid-beta Abeta(16-35) monomer and dimer by replica exchange molecular dynamics simulations: implication for full-length Abeta fibrillation. , 2009, The journal of physical chemistry. B.

[3]  Jing Zhang,et al.  Effect of protein–protein interactions on protein aggregation kinetics , 2003 .

[4]  Babatunde A. Ogunnaike,et al.  Multi-variate approach to global protein aggregation behavior and kinetics: effects of pH, NaCl, and temperature for alpha-chymotrypsinogen A. , 2010, Journal of pharmaceutical sciences.

[5]  Valentina Tozzini,et al.  Coarse-grained models for proteins. , 2005, Current opinion in structural biology.

[6]  D. Frenkel,et al.  Fluid-fluid coexistence in colloidal systems with short-ranged strongly directional attraction , 2003 .

[7]  Arup K Chakraborty,et al.  Atomistic understanding of kinetic pathways for single base-pair binding and unbinding in DNA , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[8]  Luis Fernández Pacios,et al.  Distinct Molecular Surfaces and Hydrophobicity of Amino Acid Residues in Proteins , 2001, J. Chem. Inf. Comput. Sci..

[9]  Flavio Seno,et al.  Insight into the Structure of Amyloid Fibrils from the Analysis of Globular Proteins , 2006, PLoS Comput. Biol..

[10]  C. Pace,et al.  Toward a molecular understanding of protein solubility: increased negative surface charge correlates with increased solubility. , 2012, Biophysical journal.

[11]  John F. Carpenter,et al.  Physical Stability of Proteins in Aqueous Solution: Mechanism and Driving Forces in Nonnative Protein Aggregation , 2003, Pharmaceutical Research.

[12]  R. Nussinov,et al.  Principles of protein-protein interactions: what are the preferred ways for proteins to interact? , 2008, Chemical reviews.

[13]  John E Straub,et al.  Role of water in protein aggregation and amyloid polymorphism. , 2011, Accounts of chemical research.

[14]  R. Nussinov,et al.  Protein binding versus protein folding: the role of hydrophilic bridges in protein associations. , 1997, Journal of molecular biology.

[15]  J. King,et al.  Protein misfolding and aggregation in cataract disease and prospects for prevention. , 2012, Trends in molecular medicine.

[16]  Wilfred F van Gunsteren,et al.  On developing coarse-grained models for biomolecular simulation: a review. , 2012, Physical chemistry chemical physics : PCCP.

[17]  P. Cremer,et al.  Interactions between macromolecules and ions: The Hofmeister series. , 2006, Current opinion in chemical biology.

[18]  R. Jernigan,et al.  Residue-residue potentials with a favorable contact pair term and an unfavorable high packing density term, for simulation and threading. , 1996, Journal of molecular biology.

[19]  Anthony W. Fitzpatrick,et al.  Inversion of the Balance between Hydrophobic and Hydrogen Bonding Interactions in Protein Folding and Aggregation , 2011, PLoS Comput. Biol..

[20]  Valerie L. Anderson,et al.  Identification of a helical intermediate in trifluoroethanol-induced alpha-synuclein aggregation , 2010, Proceedings of the National Academy of Sciences.

[21]  Changbong Hyeon,et al.  Capturing the essence of folding and functions of biomolecules using coarse-grained models. , 2011, Nature communications.

[22]  Bernhardt L Trout,et al.  Aggregation in protein-based biotherapeutics: computational studies and tools to identify aggregation-prone regions. , 2011, Journal of pharmaceutical sciences.

[23]  Z. H. Wang,et al.  Origin of the native driving force for protein folding. , 1999, Physical review letters.

[24]  E. Kaler,et al.  Effects of pH on protein-protein interactions and implications for protein phase behavior. , 2008, Biochimica et biophysica acta.

[25]  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.

[26]  Yuji Sugita,et al.  Replica-exchange multicanonical algorithm and multicanonical replica-exchange method for simulating systems with rough energy landscape , 2000, cond-mat/0009119.

[27]  H. M. Schaink,et al.  Determination of the osmotic second virial coefficient and the dimerization of β-lactoglobulin in aqueous solutions with added salt at the isoelectric point , 2000 .

[28]  Giuseppe Pellicane,et al.  Microscopic determination of the phase diagrams of lysozyme and γ-crystallin solutions , 2004 .

[29]  Gerhard Hummer,et al.  Coarse-grained models for simulations of multiprotein complexes: application to ubiquitin binding. , 2008, Journal of molecular biology.

[30]  J. King,et al.  Hydrophobic Core Mutations Associated with Cataract Development in Mice Destabilize Human γD-Crystallin* , 2009, The Journal of Biological Chemistry.

[31]  J. King,et al.  Inhibition of unfolding and aggregation of lens protein human gamma D crystallin by sodium citrate. , 2011, Experimental eye research.

[32]  Ken A Dill,et al.  Use of the Weighted Histogram Analysis Method for the Analysis of Simulated and Parallel Tempering Simulations. , 2007, Journal of chemical theory and computation.

[33]  Charles S Henry,et al.  Colloidal behavior of proteins: effects of the second virial coefficient on solubility, crystallization and aggregation of proteins in aqueous solution. , 2005, Current pharmaceutical biotechnology.

[34]  William F Weiss,et al.  Principles, approaches, and challenges for predicting protein aggregation rates and shelf life. , 2009, Journal of pharmaceutical sciences.

[35]  D. M. Ford.  Enthalpy-entropy compensation is not a general feature of weak association. , 2005, Journal of the American Chemical Society.

[36]  A. Lenhoff,et al.  A consistent experimental and modeling approach to light-scattering studies of protein-protein interactions in solution. , 2005, Biophysical journal.

[37]  Alexander Grünberger,et al.  Coarse-grained modeling of protein second osmotic virial coefficients: sterics and short-ranged attractions. , 2013, The journal of physical chemistry. B.

[38]  Rudesh D. Toofanny,et al.  A comprehensive multidimensional-embedded, one-dimensional reaction coordinate for protein unfolding/folding. , 2010, Biophysical journal.

[39]  Stephen L Mayo,et al.  A de novo designed protein–protein interface , 2007, Protein science : a publication of the Protein Society.

[40]  Panagiotis Katsonis,et al.  Metastable liquid clusters in super- and undersaturated protein solutions. , 2007, The journal of physical chemistry. B.

[41]  J. Thornton,et al.  Diversity of protein–protein interactions , 2003, The EMBO journal.

[42]  Cecilia Clementi,et al.  Coarse-grained models of protein folding: toy models or predictive tools? , 2008, Current opinion in structural biology.

[43]  Melissa S Kosinski-Collins,et al.  Interdomain side‐chain interactions in human γD crystallin influencing folding and stability , 2005, Protein science : a publication of the Protein Society.

[44]  D. Kofke,et al.  Mayer sampling: calculation of cluster integrals using free-energy perturbation methods. , 2004, Physical review letters.

[45]  Holger Gohlke,et al.  A natural coarse graining for simulating large biomolecular motion. , 2006, Biophysical journal.

[46]  Valentina Tozzini,et al.  Minimalist models for proteins: a comparative analysis , 2010, Quarterly Reviews of Biophysics.

[47]  C. Roberts,et al.  A quasichemical approach for protein-cluster free energies in dilute solution. , 2007, The Journal of chemical physics.

[48]  William F. Weiss,et al.  Computational design and biophysical characterization of aggregation-resistant point mutations for γD crystallin illustrate a balance of conformational stability and intrinsic aggregation propensity. , 2011, Biochemistry.

[49]  Tristan Bereau,et al.  Generic coarse-grained model for protein folding and aggregation. , 2009, The Journal of chemical physics.

[50]  J. King,et al.  Crystal cataracts: Human genetic cataract caused by protein crystallization , 2001, Proceedings of the National Academy of Sciences of the United States of America.

[51]  Erinc Sahin,et al.  Predicting solution aggregation rates for therapeutic proteins: approaches and challenges. , 2011, International journal of pharmaceutics.

[52]  O. Velev,et al.  Why is the osmotic second virial coefficient related to protein crystallization , 1999 .

[53]  Mitsunori Ikeguchi,et al.  Partial rigid‐body dynamics in NPT, NPAT and NPγT ensembles for proteins and membranes , 2004, J. Comput. Chem..

[54]  Vincent K Shen,et al.  Osmotic virial coefficients for model protein and colloidal solutions: importance of ensemble constraints in the analysis of light scattering data. , 2012, The Journal of chemical physics.

[55]  A. Minton.  Static light scattering from concentrated protein solutions, I: General theory for protein mixtures and application to self-associating proteins. , 2007, Biophysical journal.

[56]  Structure and thermodynamics of colloidal protein cluster formation: comparison of square-well and simple dipolar models. , 2009, The Journal of chemical physics.

[57]  Steven J Shire,et al.  Challenges in the development of high protein concentration formulations. , 2004, Journal of pharmaceutical sciences.

[58]  Alan M. Ferrenberg,et al.  Optimized Monte Carlo data analysis. , 1989, Physical Review Letters.

[59]  A. Lenhoff,et al.  Light-scattering studies of protein solutions: role of hydration in weak protein-protein interactions. , 2005, Biophysical journal.

[60]  G. Benedek,et al.  Phase separation in solutions of monoclonal antibodies and the effect of human serum albumin , 2011, Proceedings of the National Academy of Sciences.

[61]  Thomas F. Miller,et al.  Symplectic quaternion scheme for biophysical molecular dynamics , 2002 .

[62]  L. Serrano,et al.  Prediction of sequence-dependent and mutational effects on the aggregation of peptides and proteins , 2004, Nature Biotechnology.

[63]  Amedeo Caflisch,et al.  Interpreting the aggregation kinetics of amyloid peptides. , 2006, Journal of molecular biology.

[64]  U. Andley.  Crystallins in the eye: Function and pathology , 2007, Progress in Retinal and Eye Research.

[65]  David J Wales,et al.  Potential energy and free energy landscapes. , 2006, The journal of physical chemistry. B.

[66]  O. Velev,et al.  Protein interactions in solution characterized by light and neutron scattering: comparison of lysozyme and chymotrypsinogen. , 1998, Biophysical journal.

[67]  Francesc X. Avilés,et al.  AGGRESCAN: a server for the prediction and evaluation of "hot spots" of aggregation in polypeptides , 2007, BMC Bioinform..

[68]  M. Basan,et al.  Altered phase diagram due to a single point mutation in human γD-crystallin , 2007, Proceedings of the National Academy of Sciences.

[69]  Fractionation effects in phase equilibria of polydisperse hard-sphere colloids. , 2004, Physical review. E, Statistical, nonlinear, and soft matter physics.

[70]  D. Kalonia,et al.  Nature and consequences of protein-protein interactions in high protein concentration solutions. , 2008, International journal of pharmaceutics.

[71]  Patrick Garidel,et al.  Correlation of protein-protein interactions as assessed by affinity chromatography with colloidal protein stability: A case study with lysozyme , 2009, Pharmaceutical development and technology.

[72]  D. Raleigh,et al.  Ionic strength effects on amyloid formation by amylin are a complicated interplay among Debye screening, ion selectivity, and Hofmeister effects. , 2012, Biochemistry.

[73]  Frédéric Cardinaux,et al.  Equilibrium cluster formation in concentrated protein solutions and colloids , 2004, Nature.

[74]  Yi Li,et al.  Reexamining protein-protein and protein-solvent interactions from Kirkwood-Buff analysis of light scattering in multi-component solutions. , 2011, The Journal of chemical physics.

[75]  Huan‐Xiang Zhou.  Influence of crowded cellular environments on protein folding, binding, and oligomerization: Biological consequences and potentials of atomistic modeling , 2013, FEBS letters.

[76]  D. Kofke,et al.  Higher-order virial coefficients of water models. , 2007, The journal of physical chemistry. B.

[77]  Paolo Carloni,et al.  Hybrid Molecular Mechanics/Coarse-Grained Simulations for Structural Prediction of G-Protein Coupled Receptor/Ligand Complexes , 2012, PloS one.

[78]  C. Pang,et al.  A novel γD-crystallin mutation causes mild changes in protein properties but leads to congenital coralliform cataract , 2009, Molecular vision.