Rapid Computational Prediction of Thermostabilizing Mutations for G Protein-Coupled Receptors

G protein-coupled receptors (GPCRs) are highly dynamic and often denature when extracted in detergents. Deriving thermostable mutants has been a successful strategy to stabilize GPCRs in detergents, but this process is experimentally tedious. We have developed a computational method to predict the position of the thermostabilizing mutations for a given GPCR sequence. We have validated the method against experimentally measured thermostability data for single mutants of the β1-adrenergic receptor (β1AR), adenosine A2A receptor (A2AR) and neurotensin receptor 1 (NTSR1). To make these predictions we started from homology models of these receptors of varying accuracies and generated an ensemble of conformations by sampling the rigid body degrees of freedom of transmembrane helices. Then, an all-atom force field function was used to calculate the enthalpy gain, known as the “stability score” upon mutation of every residue, in these receptor structures, to alanine. For all three receptors, β1AR, A2AR, and NTSR1, we observed that mutations of hydrophobic residues in the transmembrane domain to alanine that have high stability scores correlate with high experimental thermostability. The prediction using the stability score improves when using an ensemble of receptor conformations compared to a single structure, showing that receptor flexibility is important. We also find that our previously developed LITiCon method for generating conformation ensembles is similar in performance to predictions using ensembles obtained from microseconds of molecular dynamics simulations (which is computationally hundred times slower than LITiCon). We improved the thermostability prediction by including other properties such as residue-based stress and the extent of allosteric communication by each residue in the stability score. Our method is the first step toward a computational method for rapid prediction of thermostable mutants of GPCRs.

[1]  A. Plückthun,et al.  Stabilizing membrane proteins through protein engineering. , 2013, Current opinion in chemical biology.

[2]  Christopher G Tate,et al.  A crystal clear solution for determining G-protein-coupled receptor structures. , 2012, Trends in biochemical sciences.

[3]  Yoko Shibata,et al.  Thermostabilisation of the neurotensin receptor NTS1 , 2009, Journal of molecular biology.

[4]  M. Congreve,et al.  Structure of the adenosine A(2A) receptor in complex with ZM241385 and the xanthines XAC and caffeine. , 2011, Structure.

[5]  C. Tate,et al.  Transferability of thermostabilizing mutations between β-adrenergic receptors , 2009, Molecular membrane biology.

[6]  N. Vaidehi,et al.  Structural insights into conformational stability of wild-type and mutant beta1-adrenergic receptor. , 2010, Biophysical journal.

[7]  Christopher G. Tate,et al.  Crystal Structures of a Stabilized β1-Adrenoceptor Bound to the Biased Agonists Bucindolol and Carvedilol , 2012, Structure.

[8]  B. Stoddard,et al.  Computational Thermostabilization of an Enzyme , 2005, Science.

[9]  N. Vaidehi,et al.  Dynamic Behavior of the Active and Inactive States of the Adenosine A2A Receptor , 2014, The journal of physical chemistry. B.

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

[11]  Gebhard F. X. Schertler,et al.  Structure of a β1-adrenergic G-protein-coupled receptor , 2008, Nature.

[12]  C. Tate,et al.  Optimising the combination of thermostabilising mutations in the neurotensin receptor for structure determination. , 2013, Biochimica et biophysica acta.

[13]  Edsger W. Dijkstra,et al.  A note on two problems in connexion with graphs , 1959, Numerische Mathematik.

[14]  Jonathan Bronson,et al.  Computational design and selections for an engineered, thermostable terpene synthase , 2011, Protein science : a publication of the Protein Society.

[15]  N. Vaidehi,et al.  Differences in allosteric communication pipelines in the inactive and active states of a GPCR. , 2014, Biophysical journal.

[16]  P. Kitabgi,et al.  Constitutive activation of the neurotensin receptor 1 by mutation of Phe358 in Helix seven , 2002, British journal of pharmacology.

[17]  Nagarajan Vaidehi,et al.  LITiCon: a discrete conformational sampling computational method for mapping various functionally selective conformational states of transmembrane helical proteins. , 2012, Methods in molecular biology.

[18]  Jens Meiler,et al.  Structure of a Class C GPCR Metabotropic Glutamate Receptor 1 Bound to an Allosteric Modulator , 2014, Science.

[19]  Nagarajan Vaidehi,et al.  Ligand-stabilized conformational states of human beta(2) adrenergic receptor: insight into G-protein-coupled receptor activation. , 2008, Biophysical journal.

[20]  J. Shiloach,et al.  Structure of the agonist-bound neurotensin receptor , 2012, Nature.

[21]  Peter M. Kasson,et al.  GROMACS 4.5: a high-throughput and highly parallel open source molecular simulation toolkit , 2013, Bioinform..

[22]  Roland L. Dunbrack,et al.  proteins STRUCTURE O FUNCTION O BIOINFORMATICS Improved prediction of protein side-chain conformations with SCWRL4 , 2022 .

[23]  C. Tate,et al.  Thermostabilisation of an Agonist-Bound Conformation of the Human Adenosine A2A Receptor , 2011, Journal of molecular biology.

[24]  Jan Steyaert,et al.  Nanobody stabilization of G protein-coupled receptor conformational states. , 2011, Current opinion in structural biology.

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

[26]  Yoko Shibata,et al.  Conformational thermostabilization of the β1-adrenergic receptor in a detergent-resistant form , 2008, Proceedings of the National Academy of Sciences.

[27]  Tong Li,et al.  Engineering a More Thermostable Blue Light Photo Receptor Bacillus subtilis YtvA LOV Domain by a Computer Aided Rational Design Method , 2013, PLoS Comput. Biol..

[28]  Andrew E. Torda,et al.  The GROMOS biomolecular simulation program package , 1999 .

[29]  Ali Jazayeri,et al.  Structure of class B GPCR corticotropin-releasing factor receptor 1 , 2013, Nature.

[30]  Joshua M. Kunken,et al.  Fusion partner toolchest for the stabilization and crystallization of G protein-coupled receptors. , 2012, Structure.

[31]  R. Stevens,et al.  GPCR Engineering Yields High-Resolution Structural Insights into β2-Adrenergic Receptor Function , 2007, Science.

[32]  C. Tate,et al.  Engineering G protein-coupled receptors to facilitate their structure determination. , 2009, Current opinion in structural biology.

[33]  Jianpeng Ma,et al.  CHARMM: The biomolecular simulation program , 2009, J. Comput. Chem..

[34]  N. Vaidehi,et al.  Thermostabilization of the β1-adrenergic receptor correlates with increased entropy of the inactive state. , 2013, The journal of physical chemistry. B.

[35]  Guocheng Du,et al.  Enhanced thermostability of keratinase by computational design and empirical mutation , 2013, Journal of Industrial Microbiology & Biotechnology.

[36]  Frauke Gräter,et al.  Implementation of force distribution analysis for molecular dynamics simulations , 2011, BMC Bioinformatics.

[37]  Patrick Barth,et al.  Naturally evolved G protein-coupled receptors adopt metastable conformations , 2012, Proceedings of the National Academy of Sciences.

[38]  Laxmikant V. Kalé,et al.  Scalable molecular dynamics with NAMD , 2005, J. Comput. Chem..

[39]  A. Doré,et al.  Structure of class C GPCR metabotropic glutamate receptor 5 transmembrane domain , 2014, Nature.

[40]  Yoko Shibata,et al.  Co-evolving stability and conformational homogeneity of the human adenosine A2a receptor , 2008, Proceedings of the National Academy of Sciences.