Importance of receptor flexibility in binding of cyclam compounds to the chemokine receptor CXCR4.

We have elucidated the binding sites of four moncyclam and one bicyclam antagonist AMD3100, in the human chemokine receptor CXCR4. Using the predicted structural models of CXCR4, we have further predicted the binding sites of these cyclam compounds. We used the computational method LITiCon to map the differences in receptor structure stabilized by the mono and bicyclam compounds. Accounting for the receptor flexibility lead to a single binding mode for the cyclam compounds, that has not been possible previously using a single receptor structural model and fixed receptor docking algorithms. There are several notable differences in the receptor conformations stabilized by monocyclam antagonist compared to a bicylam antagonist. The loading of the Cu(2+) ions in the cyclam compounds, shrinks the size of the cyclam rings and the residue D262(6.58) plays an important role in bonding to the copper ion in the monocylam compounds while residue E288(7.39) is important for the bicyclam compound.

Elucidation of binding sites of dual antagonists in the human chemokine receptors CCR2 and CCR5.

Design of dual antagonists for the chemokine receptors CCR2 and CCR5 will be greatly facilitated by knowledge of the structural differences of their binding sites. Thus, we computationally predicted the binding site of the dual CCR2/CCR5 antagonist N-dimethyl-N-[4-[[[2-(4-methylphenyl)-6,7-dihydro-5H-benzohepten-8-yl] carbonyl]amino]benzyl]tetrahydro-2H-pyran-4-aminium (TAK-779), and a CCR2-specific antagonist N-(carbamoylmethyl)-3-trifluoromethyl benzamido-parachlorobenzyl 3-aminopyrrolidine (Teijin compound 1) in an ensemble of predicted structures of human CCR2 and CCR5. Based on our predictions of the protein-ligand interactions, we examined the activity of the antagonists for cells expressing thirteen mutants of CCR2 and five mutants of CCR5. The results show that residues Trp98(2.60) and Thr292(7.40) contribute significantly to the efficacy of both TAK-779 and Teijin compound 1, whereas His121(3.33) and Ile263(6.55) contribute significantly only to the antagonistic effect of Teijin compound 1 at CCR2. Mutation of residues Trp86(2.60) and Tyr108(3.32) adversely affected the efficacy of TAK-779 in antagonizing CCR5-mediated chemotaxis. Y49A(1.39) and E291A(7.39) mutants of CCR2 showed a complete loss of CCL2 binding and chemotaxis, despite robust cell surface expression, suggesting that these residues are critical in maintaining the correct receptor architecture. Modeling studies support the hypothesis that the residues Tyr49(1.39), Trp98(2.60), Tyr120(3.32), and Glu291(7.39) of CCR2 form a tight network of aromatic cluster and polar contacts between transmembrane helices 1, 2, 3, and 7.

Ligand-stabilized conformational states of human beta(2) adrenergic receptor: insight into G-protein-coupled receptor activation.

G-protein-coupled receptors (GPCRs) are known to exist in dynamic equilibrium between inactive- and several active-state conformations, even in the absence of a ligand. Recent experimental studies on the beta(2) adrenergic receptor (beta(2)AR) indicate that structurally different ligands with varying efficacies trigger distinct conformational changes and stabilize different receptor conformations. We have developed a computational method to study the ligand-induced rotational orientation changes in the transmembrane helices of GPCRs. This method involves a systematic spanning of the rotational orientation of the transmembrane helices (TMs) that are in the vicinity of the ligand for predicting the helical rotations that occur on ligand binding. The predicted ligand-stabilized receptor conformations are characterized by a simultaneous lowering of the ligand binding energy and a significant gain in interhelical and receptor-ligand hydrogen bonds. Using the beta(2)AR as a model, we show that the receptor conformational state depends on the structure and efficacy of the ligand for a given signaling pathway. We have studied the ligand-stabilized receptor conformations of five different ligands, a full agonist, norepinephrine; a partial agonist, salbutamol; a weak partial agonist, dopamine; a very weak agonist, catechol; and an inverse agonist, ICI-115881. The predicted ligand-stabilized receptor models correlate well with the experimentally observed conformational switches in beta(2)AR, namely, the breaking of the ionic lock between R131(3.50) at the intracellular end of TM3 (part of the DRY motif) and E268(6.30) on TM6, and the rotamer toggle switch on W286(6.48) on TM6. In agreement with trp-bimane quenching experiments, we found that norepinephrine and dopamine break the ionic lock and engage the rotamer toggle switch, whereas salbutamol, a noncatechol partial agonist only breaks the ionic lock, and the weak agonist catechol only engages the rotamer toggle switch. Norepinephrine and dopamine occupy the same binding region, between TM3, TM5, and TM6, whereas the binding site of salbutamol is shifted toward TM4. Catechol binds deeper into the protein cavity compared to the other ligands, making contact with TM5 and TM6. A part of the catechol binding site overlaps with those of dopamine and norepinephrine but not with that of salbutamol. Virtual ligand screening on 10,060 ligands on the norepinephrine-stabilized receptor conformation shows an enrichment of 38% compared to ligand unbound receptor conformation. These results show that ligand-induced conformational changes are important for developing functionally specific drugs that will stabilize a particular receptor conformation. These studies represent the first step toward a more universally applicable computational method for studying ligand efficacy and GPCR activation.

Test of the Binding Threshold Hypothesis for olfactory receptors: explanation of the differential binding of ketones to the mouse and human orthologs of olfactory receptor 912-93.

We tested the Binding Threshold Hypothesis (BTH) for activation of olfactory receptors (ORs): To activate an OR, the odorant must bind to the OR with binding energy above some threshold value. The olfactory receptor (OR) 912-93 is known experimentally to be activated by ketones in mouse, but is inactive to ketones in human, despite an amino acid sequence identity of approximately 66%. To investigate the origins of this difference, we used the MembStruk first-principles method to predict the tertiary structure of the mouse OR 912-93 (mOR912-93), and the HierDock first-principles method to predict the binding site for ketones to this receptor. We found that the strong binding of ketones to mOR912-93 is dominated by a hydrogen bond of the ketone carbonyl group to Ser105. All ketones predicted to have a binding energy stronger than EBindThresh = 26 kcal/mol were observed experimentally to activate this OR, while the two ketones predicted to bind more weakly do not. In addition, we predict that 2-undecanone and 2-dodecanone both bind sufficiently strongly to activate mOR912-93. A similar binding site for ketones was predicted in hOR912-93, but the binding is much weaker because the human ortholog has a Gly at the position of Ser105. We predict that mutating this Gly to Ser in human should lead to activation of hOR912-93 by these ketones. Experimental substantiations of the above predictions would provide further tests of the validity of the BTH, our predicted 3D structures, and our predicted binding sites for these ORs.

Position of helical kinks in membrane protein crystal structures and the accuracy of computational prediction.

The structural features of helical transmembrane (TM) proteins, such as helical kinks, tilts, and rotational orientations are important in modulation of their function and these structural features give rise to functional diversity in membrane proteins with similar topology. In particular, the helical kinks caused by breaking of the backbone hydrogen bonds lead to hinge bending flexibility in these helices. Therefore it is important to understand the nature of the helical kinks and to be able to reproduce these kinks in structural models of membrane proteins. We have analyzed the position and extent of helical kinks in the transmembrane helices of all the crystal structures of membrane proteins taken from the MPtopo database, which are about 405 individual helices of length between 19 and 35 residues. 44% of the crystal structures of TM helices showed a significant helical kink, and 35% of these kinks are caused by prolines. Many of the non-proline helical kinks are caused by other residues like Ser and Gly that are located at the center of helical kinks. The side chain of Ser makes a hydrogen bond with the main chain carbonyl of the i - 4th or i + 4th residue thus making a kink. We have also studied how well molecular dynamics (MD) simulations on isolated helices can reproduce the position of the helical kinks in TM helices. Such a method is useful for structure prediction of membrane proteins. We performed MD simulations, starting from a canonical helix for the 405 TM helices. 1 ns of MD simulation results show that we can reproduce about 79% of the proline kinks, only 59% of the vestigial proline kinks and 18% of the non-proline helical kinks. We found that similar results can be obtained from choosing the lowest potential energy structure from the MD simulation. 4-14% more of the vestigial prolines were reproduced by replacing them with prolines before performing MD simulations, and changing the amino acid back to proline after the MD simulations. From these results we conclude that the position of the helical kinks is dependent on the TM sequence. However the extent of helical kinking may depend on the packing of the rest of the protein and the lipid bilayer.

Agonist-induced conformational changes in bovine rhodopsin: insight into activation of G-protein-coupled receptors.

Activation of G-protein-coupled receptors (GPCRs) is initiated by conformational changes in the transmembrane (TM) helices and the intra- and extracellular loops induced by ligand binding. Understanding the conformational changes in GPCRs leading to activation is imperative in deciphering the role of these receptors in the pathology of diseases. Since the crystal structures of activated GPCRs are not yet available, computational methods and biophysical techniques have been used to predict the structures of GPCR active states. We have recently applied the computational method LITiCon to understand the ligand-induced conformational changes in beta(2)-adrenergic receptor by ligands of varied efficacies. Here we report a study of the conformational changes associated with the activation of bovine rhodopsin for which the crystal structure of the inactive state is known. Starting from the inactive (dark) state, we have predicted the TM conformational changes that are induced by the isomerization of 11-cis retinal to all-trans retinal leading to the fully activated state, metarhodopsin II. The predicted active state of rhodopsin satisfies all of the 30 known experimental distance constraints. The predicted model also correlates well with the experimentally observed conformational switches in rhodopsin and other class A GPCRs, namely, the breaking of the ionic lock between R135(3.50) at the intracellular end of TM3 (part of the DRY motif) and E247(6.30) on TM6, and the rotamer toggle switch on W265(6.48) on TM6. We observe that the toggling of the W265(6.48) rotamer modulates the bend angle of TM6 around the conserved proline. The rotamer toggling is facilitated by the formation of a water wire connecting S298(7.45), W265(6.48) and H211(5.46). As a result, the intracellular ends of TMs 5 and 6 move outward from the protein core, causing large conformational changes at the cytoplasmic interface. The predicted outward movements of TM5 and TM6 are in agreement with the recently published crystal structure of opsin, which is proposed to be close to the active-state structure. In the predicted active state, several residues in the intracellular loops, such as R69, V139(3.54), T229, Q237, Q239, S240, T243 and V250(6.33), become more water exposed compared to the inactive state. These residues may be involved in mediating the conformational signal from the receptor to the G protein. From mutagenesis studies, some of these residues, such as V139(3.54), T229 and V250(6.33), are already implicated in G-protein activation. The predicted active state also leads to the formation of new stabilizing interhelical hydrogen-bond contacts, such as those between W265(6.48) and H211(5.46) and E122(3.37) and C167(4.56). These hydrogen-bond contacts serve as potential conformational switches offering new opportunities for future experimental investigations. The calculated retinal binding energy surface shows that binding of an agonist makes the receptor dynamic and flexible and accessible to many conformations, while binding of an inverse agonist traps the receptor in the inactive state and makes the other conformations inaccessible.

The predicted 3D structures of the human M1 muscarinic acetylcholine receptor with agonist or antagonist bound.

The muscarinic acetylcholine G-protein-coupled receptors are implicated in diseases ranging from cognitive dysfunctions to smooth-muscle disorders. To provide a structural basis for drug design, we used the MembStruk computational method to predict the 3D structure of the human M1 muscarinic receptor. We validated this structure by using the HierDock method to predict the binding sites for three agonists and four antagonists. The intermolecular ligand-receptor contacts at the predicted binding sites agree well with deductions from available mutagenesis experiments, and the calculated relative binding energies correlate with measured binding affinities. The predicted binding site of all four antagonists is located between transmembrane (TM) helices 3, 4, 5, 6, and 7, whereas the three agonists prefer a site involving residues from TM3, TM6, and TM7. We find that Trp 157(4) contributes directly to antagonist binding, whereas Pro 159(4) provides an indirect conformational switch to position Trp 157(4) in the binding site (the number in parentheses indicates the TM helix). This explains the large decrease in ligand binding affinity and signaling efficacy by mutations of Trp 157(4) and Pro 159(4) not previously explained by homology models. We also found that Asp 105(3) and aromatic residues Tyr 381(6), Tyr 404(7), and Tyr 408(7) are critical for binding the quaternary ammonium head group of the ligand through cation-pi interactions. For ligands with a charged tertiary amine head group, we suggest that proton transfer from the ligand to Asp 105(3) occurs upon binding. Furthermore, we found that an extensive aromatic network involving Tyr 106(3), Trp 157(4), Phe 197(5), Trp 378(6), and Tyr 381(6) is important in stabilizing antagonist binding. For antagonists with two terminal phenyl rings, this aromatic network extends to Trp 164(4), Tyr 179(extracellular loop 2), and Phe 390(6) located at the extracellular end of the TMs. We find that Asn 382(6) forms hydrogen bonds with selected antagonists. Tyr381(6) and Ser 109(3) form hydrogen bonds with the ester moiety of acetylcholine, which binds in the gauche conformation.

Predicted 3-D structures for mouse I7 and rat I7 olfactory receptors and comparison of predicted odor recognition profiles with experiment.

The first step in the perception of an odor is the activation of one or more olfactory receptors (ORs) following binding of the odorant molecule to the OR. In order to initiate the process of determining how the molecular level receptor-odorant interactions are related to odor perception, we used the MembStruk computational method to predict the three-dimensional (3-D) structure of the I7 OR for both mouse and rat. We then used the HierDock ligand docking computational method to predict the binding site and binding energy for the library of 56 odorants to these receptors for which experiment response data are now available. We find that the predicted 3-D structures of the mouse and rat I7 OR lead to predictions of odorant binding that are in good agreement with the experimental results, thus validating the accuracy of both the 3-D structure and the predicted binding site. In particular we predict that heptanal and octanal both bind strongly to both mouse and rat I7 ORs, which conflicts with the older literature but agrees with recent experiments. To provide the basis of additional validations of our 3-D structures, we also report the odorant binding site for a new odorant (8-hydroxy-octanal) with a novel functionality designed to bind strongly to mouse I7. Such validated computational methods should be very useful in predicting the structure and function of many other ORs.

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.

Predicted 3D structure for the human beta 2 adrenergic receptor and its binding site for agonists and antagonists.

We report the 3D structure of human beta2 adrenergic receptor (AR) predicted by using the MembStruk first principles method. To validate this structure, we use the HierDock first principles method to predict the ligand-binding sites for epinephrine and norepinephrine and for eight other ligands, including agonists and antagonists to beta 2 AR and ligands not observed to bind to beta 2 AR. The binding sites agree well with available mutagenesis data, and the calculated relative binding energies correlate reasonably with measured binding affinities. In addition, we find characteristic differences in the predicted binding sites of known agonists and antagonists that allow us to infer the likely activity of other ligands. The predicted ligand-binding properties validate the methods used to predict the 3D structure and function. This validation is a successful step toward applying these procedures to predict the 3D structures and function of the other eight subtypes of ARs, which should enable the development of subtype-specific antagonists and agonists with reduced side effects.

Prediction of structure and function of G protein-coupled receptors.

G protein-coupled receptors (GPCRs) mediate our sense of vision, smell, taste, and pain. They are also involved in cell recognition and communication processes, and hence have emerged as a prominent superfamily for drug targets. Unfortunately, the atomic-level structure is available for only one GPCR (bovine rhodopsin), making it difficult to use structure-based methods to design drugs and mutation experiments. We have recently developed first principles methods (MembStruk and HierDock) for predicting structure of GPCRs, and for predicting the ligand binding sites and relative binding affinities. Comparing to the one case with structural data, bovine rhodopsin, we find good accuracy in both the structure of the protein and of the bound ligand. We report here the application of MembStruk and HierDock to beta1-adrenergic receptor, endothelial differential gene 6, mouse and rat I7 olfactory receptors, and human sweet receptor. We find that the predicted structure of beta1-adrenergic receptor leads to a binding site for epinephrine that agrees well with the mutation experiments. Similarly the predicted binding sites and affinities for endothelial differential gene 6, mouse and rat I7 olfactory receptors, and human sweet receptor are consistent with the available experimental data. These predicted structures and binding sites allow the design of mutation experiments to validate and improve the structure and function prediction methods. As these structures are validated they can be used as targets for the design of new receptor-selective antagonists or agonists for GPCRs.

The predicted 3D structure of the human D2 dopamine receptor and the binding site and binding affinities for agonists and antagonists.

Dopamine neurotransmitter and its receptors play a critical role in the cell signaling process responsible for information transfer in neurons functioning in the nervous system. Development of improved therapeutics for such disorders as Parkinson's disease and schizophrenia would be significantly enhanced with the availability of the 3D structure for the dopamine receptors and of the binding site for dopamine and other agonists and antagonists. We report here the 3D structure of the long isoform of the human D2 dopamine receptor, predicted from primary sequence using first-principles theoretical and computational techniques (i.e., we did not use bioinformatic or experimental 3D structural information in predicting structures). The predicted 3D structure is validated by comparison of the predicted binding site and the relative binding affinities of dopamine, three known dopamine agonists (antiparkinsonian), and seven known antagonists (antipsychotic) in the D2 receptor to experimentally determined values. These structures correctly predict the critical residues for binding dopamine and several antagonists, identified by mutation studies, and give relative binding affinities that correlate well with experiments. The predicted binding site for dopamine and agonists is located between transmembrane (TM) helices 3, 4, 5, and 6, whereas the best antagonists bind to a site involving TM helices 2, 3, 4, 6, and 7 with minimal contacts to TM helix 5. We identify characteristic differences between the binding sites of agonists and antagonists.