First principles predictions of the structure and function of g-protein-coupled receptors: validation for bovine rhodopsin.

G-protein-coupled receptors (GPCRs) are involved in cell communication processes and with mediating such senses as vision, smell, taste, and pain. They constitute a prominent superfamily of drug targets, but an atomic-level structure is available for only one GPCR, bovine rhodopsin, making it difficult to use structure-based methods to design receptor-specific drugs. We have developed the MembStruk first principles computational method for predicting the three-dimensional structure of GPCRs. In this article we validate the MembStruk procedure by comparing its predictions with the high-resolution crystal structure of bovine rhodopsin. The crystal structure of bovine rhodopsin has the second extracellular (EC-II) loop closed over the transmembrane regions by making a disulfide linkage between Cys-110 and Cys-187, but we speculate that opening this loop may play a role in the activation process of the receptor through the cysteine linkage with helix 3. Consequently we predicted two structures for bovine rhodopsin from the primary sequence (with no input from the crystal structure)-one with the EC-II loop closed as in the crystal structure, and the other with the EC-II loop open. The MembStruk-predicted structure of bovine rhodopsin with the closed EC-II loop deviates from the crystal by 2.84 A coordinate root mean-square (CRMS) in the transmembrane region main-chain atoms. The predicted three-dimensional structures for other GPCRs can be validated only by predicting binding sites and energies for various ligands. For such predictions we developed the HierDock first principles computational method. We validate HierDock by predicting the binding site of 11-cis-retinal in the crystal structure of bovine rhodopsin. Scanning the whole protein without using any prior knowledge of the binding site, we find that the best scoring conformation in rhodopsin is 1.1 A CRMS from the crystal structure for the ligand atoms. This predicted conformation has the carbonyl O only 2.82 A from the N of Lys-296. Making this Schiff base bond and minimizing leads to a final conformation only 0.62 A CRMS from the crystal structure. We also used HierDock to predict the binding site of 11-cis-retinal in the MembStruk-predicted structure of bovine rhodopsin (closed loop). Scanning the whole protein structure leads to a structure in which the carbonyl O is only 2.85 A from the N of Lys-296. Making this Schiff base bond and minimizing leads to a final conformation only 2.92 A CRMS from the crystal structure. The good agreement of the ab initio-predicted protein structures and ligand binding site with experiment validates the use of the MembStruk and HierDock first principles' methods. Since these methods are generic and applicable to any GPCR, they should be useful in predicting the structures of other GPCRs and the binding site of ligands to these proteins.

[1]  Bernard Maigret,et al.  Rhodopsin crystal: new template yielding realistic models of G-protein-coupled receptors? , 2003, Trends in pharmacological sciences.

[2]  William A. Goddard,et al.  Fidelity of Phenylalanyl-tRNA Synthetase in Binding the Natural Amino Acids , 2003 .

[3]  C. Cowan,et al.  A comparison of the efficiency of G protein activation by ligand-free and light-activated forms of rhodopsin. , 1997, Biophysical journal.

[4]  S. L. Mayo,et al.  DREIDING: A generic force field for molecular simulations , 1990 .

[5]  A. Sali,et al.  Comparative protein structure modeling of genes and genomes. , 2000, Annual review of biophysics and biomolecular structure.

[6]  Andrei L. Lomize,et al.  Structural organization of G-protein-coupled receptors , 1999, J. Comput. Aided Mol. Des..

[7]  C. DeLisi,et al.  Hydrophobicity scales and computational techniques for detecting amphipathic structures in proteins. , 1987, Journal of molecular biology.

[8]  P Herzyk,et al.  Automated method for modeling seven-helix transmembrane receptors from experimental data. , 1995, Biophysical journal.

[9]  J. Thornton,et al.  PROCHECK: a program to check the stereochemical quality of protein structures , 1993 .

[10]  G. Heijne,et al.  Genome‐wide analysis of integral membrane proteins from eubacterial, archaean, and eukaryotic organisms , 1998, Protein science : a publication of the Protein Society.

[11]  F. Glaser,et al.  Modeling the 3D structure of GPCRs from sequence , 2001, Medicinal research reviews.

[12]  G. Schertler,et al.  Structure of rhodopsin. , 2007, Eye.

[13]  Peter L. Freddolino,et al.  The predicted 3D structure of the human D2 dopamine receptor and the binding site and binding affinities for agonists and antagonists. , 2004, Proceedings of the National Academy of Sciences of the United States of America.

[14]  D. Donnelly Modelling alpha-helical transmembrane domains. , 1993, Biochemical Society transactions.

[15]  Abhinandan Jain,et al.  Constant temperature constrained molecular dynamics: The Newton-Euler inverse mass operator method , 1996 .

[16]  Francesca Fanelli,et al.  Mutagenesis and modelling of the alpha(1b)-adrenergic receptor highlight the role of the helix 3/helix 6 interface in receptor activation. , 2002, Molecular pharmacology.

[17]  G M Shepherd,et al.  Molecular mechanisms underlying differential odor responses of a mouse olfactory receptor. , 2000, Proceedings of the National Academy of Sciences of the United States of America.

[18]  W. Goddard,et al.  Ab Initio Quantum Mechanical Study of the Structures and Energies for the Pseudorotation of 5‘-Dehydroxy Analogues of 2‘-Deoxyribose and Ribose Sugars , 1999 .

[19]  W. Goddard,et al.  Mechanism for antibody catalysis of the oxidation of water by singlet dioxygen , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[20]  L. Buck,et al.  Combinatorial Receptor Codes for Odors , 1999, Cell.

[21]  Peter L. Freddolino,et al.  Predicted 3D structure for the human beta 2 adrenergic receptor and its binding site for agonists and antagonists. , 2004, Proceedings of the National Academy of Sciences of the United States of America.

[22]  John P. Overington,et al.  The prediction and orientation of alpha-helices from sequence alignments: the combined use of environment-dependent substitution tables, Fourier transform methods and helix capping rules. , 1994, Protein engineering.

[23]  F Glaser,et al.  Modeling the 3D structure of GPCRs from sequence , 2001, Medicinal research reviews.

[24]  William A. Goddard,et al.  Predicted 3 D structure for the human 2 adrenergic receptor and its binding site for agonists and antagonists , 2004 .

[25]  Nagarajan Vaidehi,et al.  Virtual screening for binding of phenylalanine analogues to phenylalanyl-tRNA synthetase. , 2002, Journal of the American Chemical Society.

[26]  J. Klein-Seetharaman,et al.  Structure and function in rhodopsin: mapping light-dependent changes in distance between residue 65 in helix TM1 and residues in the sequence 306-319 at the cytoplasmic end of helix TM7 and in helix H8. , 2001, Biochemistry.

[27]  W. Goddard,et al.  Charge equilibration for molecular dynamics simulations , 1991 .

[28]  Abhinandan Jain,et al.  Protein simulations using techniques suitable for very large systems: The cell multipole method for nonbond interactions and the Newton‐Euler inverse mass operator method for internal coordinate dynamics , 1994, Proteins.

[29]  G. Schertler,et al.  Structure of rhodopsin , 1998, Eye.

[30]  S. W. Lin,et al.  Specific tryptophan UV-absorbance changes are probes of the transition of rhodopsin to its active state. , 1996, Biochemistry.

[31]  D. Donnelly Modelling α-helical transmembrane domains , 1993 .

[32]  K. Palczewski,et al.  Crystal Structure of Rhodopsin: A G‐Protein‐Coupled Receptor , 2000, Science.

[33]  D Bergsma,et al.  Orphan G-protein coupled receptors: novel drug targets for the pharmaceutical industry. , 2000, Drug design and discovery.

[34]  K. Nakanishi,et al.  Movement of retinal along the visual transduction path. , 2000, Science.

[35]  W. Goddard,et al.  Atomic level simulations on a million particles: The cell multipole method for Coulomb and London nonbond interactions , 1992 .

[36]  Peter L. Freddolino,et al.  Prediction of structure and function of G protein-coupled receptors , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[37]  J. Thompson,et al.  CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. , 1994, Nucleic acids research.

[38]  Klaus Schulten,et al.  Molecular dynamics investigation of primary photoinduced events in the activation of rhodopsin. , 2002, Biophysical journal.

[39]  E. Meng,et al.  Rhodopsin Sees the Light , 2000, Science.

[40]  Sharon Brunett,et al.  Molecular dynamics for very large systems on massively parallel computers: The MPSim program , 1997 .

[41]  Abhinandan Jain,et al.  A fast recursive algorithm for molecular dynamics simulation , 1993 .

[42]  H. Khorana,et al.  Structural features and light-dependent changes in the cytoplasmic interhelical E-F loop region of rhodopsin: a site-directed spin-labeling study. , 1996, Biochemistry.

[43]  Wely B. Floriano,et al.  Interaction of e. coli outer‐membrane protein A with sugars on the receptors of the brain microvascular endothelial cells , 2002, Proteins.

[44]  D. Eisenberg,et al.  The hydrophobic moment detects periodicity in protein hydrophobicity. , 1984, Proceedings of the National Academy of Sciences of the United States of America.

[45]  C. Strader,et al.  Structure and function of G protein-coupled receptors. , 1994, Annual review of biochemistry.

[46]  K. Palczewski,et al.  Activation of rhodopsin: new insights from structural and biochemical studies. , 2001, Trends in biochemical sciences.

[47]  J. Klein-Seetharaman,et al.  Structure and function in rhodopsin: mapping light-dependent changes in distance between residue 316 in helix 8 and residues in the sequence 60-75, covering the cytoplasmic end of helices TM1 and TM2 and their connection loop CL1. , 2001, Biochemistry.

[48]  E. Myers,et al.  Basic local alignment search tool. , 1990, Journal of molecular biology.

[49]  Structure. Rhodopsin sees the light. , 2000, Science.

[50]  Nagarajan Vaidehi,et al.  Making sense of olfaction through predictions of the 3-D structure and function of olfactory receptors. , 2004, Chemical senses.

[51]  H. Khorana,et al.  Requirement of Rigid-Body Motion of Transmembrane Helices for Light Activation of Rhodopsin , 1996, Science.

[52]  Nagarajan Vaidehi,et al.  HierVLS hierarchical docking protocol for virtual ligand screening of large-molecule databases. , 2004, Journal of medicinal chemistry.

[53]  Velin Z. Spassov,et al.  Introducing an Implicit Membrane in Generalized Born/Solvent Accessibility Continuum Solvent Models , 2002 .

[54]  G. Zamanakos A fast and accurate analytical method for the computation of solvent effects in molecular simulations , 2002 .

[55]  Todd J. A. Ewing,et al.  Critical evaluation of search algorithms for automated molecular docking and database screening , 1997, J. Comput. Chem..

[56]  Roland L. Dunbrack,et al.  Prediction of protein side-chain rotamers from a backbone-dependent rotamer library: a new homology modeling tool. , 1997, Journal of molecular biology.

[57]  T. Gudermann,et al.  The structural basis of G-protein-coupled receptor function and dysfunction in human diseases. , 2002, Reviews of physiology, biochemistry and pharmacology.

[58]  P. Strange Three-state and two-state models. , 1998, Trends in pharmacological sciences.

[59]  D. Eisenberg,et al.  Hydrophobic moments and protein structure , 1982 .

[60]  G Vriend,et al.  WHAT IF: a molecular modeling and drug design program. , 1990, Journal of molecular graphics.

[61]  Alexander D. MacKerell,et al.  All-atom empirical potential for molecular modeling and dynamics studies of proteins. , 1998, The journal of physical chemistry. B.