论文信息 - Ab initio computational modeling of loops in G‐protein‐coupled receptors: Lessons from the crystal structure of rhodopsin

Ab initio computational modeling of loops in G‐protein‐coupled receptors: Lessons from the crystal structure of rhodopsin

With the help of the crystal structure of rhodopsin an ab initio method has been developed to calculate the three‐dimensional structure of the loops that connect the transmembrane helices (TMHs). The goal of this procedure is to calculate the loop structures in other G‐protein coupled receptors (GPCRs) for which only model coordinates of the TMHs are available. To mimic this situation a construct of rhodopsin was used that only includes the experimental coordinates of the TMHs while the rest of the structure, including the terminal domains, has been removed. To calculate the structure of the loops a method was designed based on Monte Carlo (MC) simulations which use a temperature annealing protocol, and a scaled collective variables (SCV) technique with proper structural constraints. Because only part of the protein is used in the calculations the usual approach of modeling loops, which consists of finding a single, lowest energy conformation of the system, is abandoned because such a single structure may not be a representative member of the native ensemble. Instead, the method was designed to generate structural ensembles from which the single lowest free energy ensemble is identified as representative of the native folding of the loop. To find the native ensemble a successive series of SCV‐MC simulations are carried out to allow the loops to undergo structural changes in a controlled manner. To increase the chances of finding the native funnel for the loop, some of the SCV‐MC simulations are carried out at elevated temperatures. The native ensemble can be identified by an MC search starting from any conformation already in the native funnel. The hypothesis is that native structures are trapped in the conformational space because of the high‐energy barriers that surround the native funnel. The existence of such ensembles is demonstrated by generating multiple copies of the loops from their crystal structures in rhodopsin and carrying out an extended SCV‐MC search. For the extracellular loops e1 and e3, and the intracellular loop i1 that were used in this work, the procedure resulted in dense clusters of structures with Cα‐RMSD ∼0.5 Å. To test the predictive power of the method the crystal structure of each loop was replaced by its extended conformations. For e1 and i1 the procedure identifies native clusters with Cα‐RMSD ∼0.5 Å and good structural overlap of the side chains; for e3, two clusters were found with Cα‐RMSD ∼1.1 Å each, but with poor overlap of the side chains. Further searching led to a single cluster with lower Cα‐RMSD but higher energy than the two previous clusters. This discrepancy was found to be due to the missing elements in the constructs available from experiment for use in the calculations. Because this problem will likely appear whenever parts of the structural information are missing, possible solutions are discussed. Proteins 2006. © 2006 Wiley‐Liss, Inc.

[1] M. DePristo,et al. Ab initio construction of polypeptide fragments: Accuracy of loop decoy discrimination by an all‐atom statistical potential and the AMBER force field with the Generalized Born solvation model , 2003, Proteins.

[2] J. Ballesteros,et al. Structural mimicry in G protein-coupled receptors: implications of the high-resolution structure of rhodopsin for structure-function analysis of rhodopsin-like receptors. , 2001, Molecular pharmacology.

[3] K. Palczewski,et al. Crystal Structure of Rhodopsin: A G‐Protein‐Coupled Receptor , 2002, Chembiochem : a European journal of chemical biology.

[4] David Chandler,et al. Hydrophobicity: Two faces of water , 2002, Nature.

[5] U H Hansmann,et al. New Monte Carlo algorithms for protein folding. , 1999, Current opinion in structural biology.

[6] S. Hassan. Amino acid side chain interactions in the presence of salts. , 2005, The journal of physical chemistry. B.

[7] Bernard R. Brooks,et al. Solvent-Induced Forces between Two Hydrophilic Groups , 1994 .

[8] R. Abagyan,et al. Biased probability Monte Carlo conformational searches and electrostatic calculations for peptides and proteins. , 1994, Journal of molecular biology.

[9] Luhua Lai,et al. A fast and efficient program for modeling protein loops , 1997 .

[10] Dmitri I Svergun,et al. Addition of missing loops and domains to protein models by x-ray solution scattering. , 2002, Biophysical journal.

[11] B. Honig,et al. A hierarchical approach to all‐atom protein loop prediction , 2004, Proteins.

[12] Michael A. Patton,et al. Mutations in PTPN11, encoding the protein tyrosine phosphatase SHP-2, cause Noonan syndrome , 2001, Nature Genetics.

[13] Lei Shi,et al. The binding site of aminergic G protein-coupled receptors: the transmembrane segments and second extracellular loop. , 2002, Annual review of pharmacology and toxicology.

[14] Krzysztof Palczewski,et al. The crystallographic model of rhodopsin and its use in studies of other G protein-coupled receptors. , 2003, Annual review of biophysics and biomolecular structure.

[15] D. Baker,et al. Clustering of low-energy conformations near the native structures of small proteins. , 1998, Proceedings of the National Academy of Sciences of the United States of America.

[16] D. Baker,et al. Modeling structurally variable regions in homologous proteins with rosetta , 2004, Proteins.

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

[18] Pieter Rein ten Wolde,et al. Drying-induced hydrophobic polymer collapse , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[19] Bernhard Pfeiffer,et al. Over one hundred peptide-activated G protein-coupled receptors recognize ligands with turn structure. , 2005, Chemical reviews.

[20] R A Friesner,et al. Prediction of loop geometries using a generalized born model of solvation effects , 1999, Proteins.

[21] Christopher J. Woods,et al. Enhanced configurational sampling in binding free-energy calculations , 2003 .

[22] Yoshinori Shichida,et al. Functional role of internal water molecules in rhodopsin revealed by x-ray crystallography , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[23] H. Weinstein,et al. Mechanistic Insights from a Refined Three-dimensional Model of Integrin αIIbβ3* , 2004, Journal of Biological Chemistry.

[24] Sergio A Hassan,et al. Molecular dynamics simulations of peptides and proteins with a continuum electrostatic model based on screened Coulomb potentials , 2003, Proteins.

[25] Peter A. Kollman,et al. A general analysis of noncovalent intermolecular interactions , 1977 .

[26] Ernest L. Mehler,et al. A general screened Coulomb potential based implicit solvent model: Calculation of secondary structure of small peptides , 2001 .

[27] Francesca Fanelli,et al. Computational modeling approaches to structure-function analysis of G protein-coupled receptors. , 2005, Chemical reviews.

[28] Henry A. Gabb,et al. Collective‐variable Monte Carlo simulation of DNA , 1997 .

[29] G. Némethy,et al. Hydrophobic interactions. , 1967, Angewandte Chemie.

[30] Sergio A Hassan,et al. Long dynamics simulations of proteins using atomistic force fields and a continuum representation of solvent effects: Calculation of structural and dynamic properties , 2005, Proteins.