Three-dimensional modeling of and ligand docking to vitamin D receptor ligand binding domain

The ligand binding domain of the human vitamin D receptor (VDR) was modeled based on the crystal structure of the retinoic acid receptor. The ligand binding pocket of our VDR model is spacious at the helix 11 site and confined at the β-turn site. The ligand 1α,25-dihydroxyvitamin D3 was assumed to be anchored in the ligand binding pocket with its side chain heading to helix 11 (site 2) and the A-ring toward the β-turn (site 1). Three residues forming hydrogen bonds with the functionally important 1α- and 25-hydroxyl groups of 1α,25-dihydroxyvitamin D3 were identified and confirmed by mutational analysis: the 1α-hydroxyl group is forming pincer-type hydrogen bonds with S237 and R274 and the 25-hydroxyl group is interacting with H397. Docking potential for various ligands to the VDR model was examined, and the results are in good agreement with our previous three-dimensional structure-function theory.

Small antibody-like proteins with prescribed ligand specificities derived from the lipocalin fold

We demonstrate that the ligand pocket of a lipocalin from Pieris brassicae, the bilin-binding protein (BBP), can be reshaped by combinatorial protein design such that it recognizes fluorescein, an established immunological hapten. For this purpose 16 residues at the center of the binding site, which is formed by four loops on top of an eight-stranded β-barrel, were subjected to random mutagenesis. Fluorescein-binding BBP variants were then selected from the mutant library by bacterial phage display. Three variants were identified that complex fluorescein with high affinity, exhibiting dissociation constants as low as 35.2 nM. Notably, one of these variants effects almost complete quenching of the ligand fluorescence, similarly as an anti-fluorescein antibody. Detailed ligand-binding studies and site-directed mutagenesis experiments indicated (i) that the molecular recognition of fluorescein is specific and (ii) that charged residues at the center of the pocket are responsible for tight complex formation. Sequence comparison of the BBP variants directed against fluorescein with the wild-type protein and with further variants that were selected against several other ligands revealed that all of the randomized amino acid positions are variable. Hence, a lipocalin can be used for generating molecular pockets with a diversity of shapes. We term this class of engineered proteins “anticalins.” Their one-domain scaffold makes them a promising alternative to antibodies to create a stable receptor protein for a ligand of choice.

A peroxisome proliferator-activated receptor γ ligand inhibits adipocyte differentiation

The peroxisome proliferator-activated receptors (PPARs) are nuclear hormone receptors that regulate glucose and lipid homeostasis. The PPARγ subtype plays a central role in the regulation of adipogenesis and is the molecular target for the 2,4-thiazolidinedione class of antidiabetic drugs. Structural studies have revealed that agonist ligands activate the PPARs through direct interactions with the C-terminal region of the ligand-binding domain, which includes the activation function 2 helix. GW0072 was identified as a high-affinity PPARγ ligand that was a weak partial agonist of PPARγ transactivation. X-ray crystallography revealed that GW0072 occupied the ligand-binding pocket by using different epitopes than the known PPAR agonists and did not interact with the activation function 2 helix. In cell culture, GW0072 was a potent antagonist of adipocyte differentiation. These results establish an approach to the design of PPAR ligands with modified biological activities.

Solution structure and peptide binding studies of the C-terminal Src homology 3-like domain of the diphtheria toxin repressor protein

The diphtheria toxin repressor (DtxR) is the best-characterized member of a family of homologous proteins that regulate iron uptake and virulence gene expression in the Gram-positive bacteria. DtxR contains two domains that are separated by a short, unstructured linker. The N-terminal domain is structurally well-defined and is responsible for Fe2+ binding, dimerization, and DNA binding. The C-terminal domain adopts a fold similar to eukaryotic Src homology 3 domains, but the functional role of the C-terminal domain in repressor activity is unknown. The solution structure of the C-terminal domain, consisting of residues N130-L226 plus a 13-residue N-terminal extension, has been determined by using NMR spectroscopy. Residues before A147 are highly mobile and adopt a random coil conformation, but residues A147-L226 form a single structured domain consisting of five β-strands and three helices arranged into a partially orthogonal, two-sheet β-barrel, similar to the structure observed in the crystalline Co2+ complex of full-length DtxR. Chemical shift perturbation studies demonstrate that a proline-rich peptide corresponding to residues R125-G139 of intact DtxR binds to the C-terminal domain in a pocket formed by residues in β-strands 2, 3, and 5, and helix 3. Binding of the proline-rich peptide by the C-terminal domain of DtxR presents an example of peptide binding by a prokaryotic Src homology 3-like protein. The results of this study, combined with previous x-ray studies of intact DtxR, provide insights into a possible biological function of the C-terminal domain in regulating repressor activity.

NMR reveals hydrogen bonds between oxygen and distal histidines in oxyhemoglobin

Compared with free heme, the proteins hemoglobin (Hb) and myoglobin (Mb) exhibit greatly enhanced affinity for oxygen relative to carbon monoxide. This physiologically vital property has been attributed to either steric hindrance of CO or stabilization of O2 binding by a hydrogen bond with the distal histidine. We report here the first direct evidence of such a hydrogen bond in both α- and β-chains of oxyhemoglobin, as revealed by heteronuclear NMR spectra of chain-selectively labeled samples. Using these spectra, we have assigned the imidazole ring 1H and 15N chemical shifts of the proximal and distal histidines in both carbonmonoxy- and oxy-Hb. Because of their proximity to the heme, these chemical shifts are extremely sensitive to the heme pocket conformation. Comparison of the measured chemical shifts with values predicted from x-ray structures suggests differences between the solution and crystal structures of oxy-Hb. The chemical shift discrepancies could be accounted for by very small displacements of the proximal and distal histidines. This suggests that NMR could be used to obtain very high-resolution heme pocket structures of Hb, Mb, and other heme proteins.

Molecular mechanisms underlying differential odor responses of a mouse olfactory receptor

The prevailing paradigm for G protein-coupled receptors is that each receptor is narrowly tuned to its ligand and closely related agonists. An outstanding problem is whether this paradigm applies to olfactory receptor (ORs), which is the largest gene family in the genome, in which each of 1,000 different G protein-coupled receptors is believed to interact with a range of different odor molecules from the many thousands that comprise “odor space.” Insights into how these interactions occur are essential for understanding the sense of smell. Key questions are: (i) Is there a binding pocket? (ii) Which amino acid residues in the binding pocket contribute to peak affinities? (iii) How do affinities change with changes in agonist structure? To approach these questions, we have combined single-cell PCR results [Malnic, B., Hirono, J., Sato, T. & Buck, L. B. (1999) Cell 96, 713–723] and well-established molecular dynamics methods to model the structure of a specific OR (OR S25) and its interactions with 24 odor compounds. This receptor structure not only points to a likely odor-binding site but also independently predicts the two compounds that experimentally best activate OR S25. The results provide a mechanistic model for olfactory transduction at the molecular level and show how the basic G protein-coupled receptor template is adapted for encoding the enormous odor space. This combined approach can significantly enhance the identification of ligands for the many members of the OR family and also may shed light on other protein families that exhibit broad specificities, such as chemokine receptors and P450 oxidases.

Two amino acid substitutions convert a guanylyl cyclase, RetGC-1, into an adenylyl cyclase

Guanylyl cyclases (GCs) and adenylyl cyclases (ACs) have fundamental roles in a wide range of cellular processes. Whereas GCs use GTP as a substrate to form cGMP, ACs catalyze the analogous conversion of ATP to cAMP. Previously, a model based on the structure of adenylate cyclase was used to predict the structure of the nucleotide-binding pocket of a membrane guanylyl cyclase, RetGC-1. Based on this model, we replaced specific amino acids in the guanine-binding pocket of GC with their counterparts from AC. A change of two amino acids, E925K together with C995D, is sufficient to completely alter the nucleotide specificity from GTP to ATP. These experiments strongly validate the AC-derived RetGC-1 structural model and functionally confirm the role of these residues in nucleotide discrimination.

Structural requirements for antigen presentation by mouse CD1

The structural basis for the T cell response to glycolipid antigens (Ags) remains poorly understood. T lymphocytes autoreactive for mouse CD1 (mCD1.1) or reactive for the glycosphingolipid αgalactosylceramide (α-GalCer) presented by mCD1.1 have been described previously. In this paper it is shown that mutations at the top of the α helices and in the bottom of the Ag-binding groove can disrupt both mCD1.1 autoreactivity and α-GalCer recognition. The locations of the positions that affect T cell responses indicate that recognition of mCD1.1 is not likely to be unconventional or superantigen-like. Furthermore, the effects of the bottom of the pocket mutation suggest that the autoreactive response could require an autologous ligand, and they indicate that α-GalCer binds to the groove of mCD1.1, most likely with the shorter 18-carbon hydrophobic chain in the A′ pocket. Natural killer T cell hybridomas with identical T cell antigen receptor (TCR) α chains and different β chains respond differently to α-GalCer presented by mCD1.1 mutants. This finding indicates a role for TCR β in defining natural killer T cell specificity, despite the more restricted diversity of the α chains in these cells. Overall, the data are consistent with a mode of lipoglycan recognition similar to that proposed for glycopeptides, in which the TCR α and β chains survey a surface composed of both mCD1.1 and the carbohydrate portion of α-GalCer.

Redesigning an FKBP–ligand interface to generate chemical dimerizers with novel specificity

FKBP ligand homodimers can be used to activate signaling events inside cells and animals that have been engineered to express fusions between appropriate signaling domains and FKBP. However, use of these dimerizers in vivo is potentially limited by ligand binding to endogenous FKBP. We have designed ligands that bind specifically to a mutated FKBP over the wild-type protein by remodeling an FKBP-ligand interface to introduce a specificity binding pocket. A compound bearing an ethyl substituent in place of a carbonyl group exhibited sub-nanomolar affinity and 1,000-fold selectivity for a mutant FKBP with a compensating truncation of a phenylalanine residue. Structural and functional analysis of the new pocket showed that recognition is surprisingly relaxed, with the modified ligand only partially filling the engineered cavity. We incorporated the specificity pocket into a fusion protein containing FKBP and the intracellular domain of the Fas receptor. Cells expressing this modified chimeric protein potently underwent apoptosis in response to AP1903, a homodimer of the modified ligand, both in culture and when implanted into mice. Remodeled dimerizers such as AP1903 are ideal reagents for controlling the activities of cells that have been modified by gene therapy procedures, without interference from endogenous FKBP.

Interaction of the Hsp70 molecular chaperone, DnaK, with its cochaperone DnaJ

Chaperones of the Hsp70 family bind to unfolded or partially folded polypeptides to facilitate many cellular processes. ATP hydrolysis and substrate binding, the two key molecular activities of this chaperone, are modulated by the cochaperone DnaJ. By using both genetic and biochemical approaches, we provide evidence that DnaJ binds to at least two sites on the Escherichia coli Hsp70 family member DnaK: under the ATPase domain in a cleft between its two subdomains and at or near the pocket of substrate binding. The lower cleft of the ATPase domain is defined as a binding pocket for the J-domain because (i) a DnaK mutation located in this cleft (R167H) is an allele-specific suppressor of the binding defect of the DnaJ mutation, D35N and (ii) alanine substitution of two residues close to R167 in the crystal structure, N170A and T173A, significantly decrease DnaJ binding. A second binding determinant is likely to be in the substrate-binding domain because some DnaK mutations in the vicinity of the substrate-binding pocket are defective in either the affinity (G400D, G539D) or rate (D526N) of both peptide and DnaJ binding to DnaK. Binding of DnaJ may propagate conformational changes to the nearby ATPase catalytic center and substrate-binding sites as well as facilitate communication between these two domains to alter the molecular properties of Hsp70.

Crystal structure of the Holliday junction migration motor protein RuvB from Thermus thermophilus HB8

We report here the crystal structure of the RuvB motor protein from Thermus thermophilus HB8, which drives branch migration of the Holliday junction during homologous recombination. RuvB has a crescent-like architecture consisting of three consecutive domains, the first two of which are involved in ATP binding and hydrolysis. DNA is likely to interact with a large basic cleft, which encompasses the ATP-binding pocket and domain boundaries, whereas the junction-recognition protein RuvA may bind a flexible β-hairpin protruding from the N-terminal domain. The structures of two subunits, related by a noncrystallographic pseudo-2-fold axis, imply that conformational changes of motor protein coupled with ATP hydrolysis may reflect motility essential for its translocation around double-stranded DNA.

In vitro assembly of phytochrome B apoprotein with synthetic analogs of the phytochrome chromophore

Phytochrome B (PhyB), one of the major photosensory chromoproteins in plants, mediates a variety of light-responsive developmental processes in a photoreversible manner. To analyze the structural requirements of the chromophore for the spectral properties of PhyB, we have designed and chemically synthesized 20 analogs of the linear tetrapyrrole (bilin) chromophore and reconstituted them with PhyB apoprotein (PHYB). The A-ring acts mainly as the anchor for ligation to PHYB, because the modification of the side chains at the C2 and C3 positions did not significantly influence the formation or difference spectra of adducts. In contrast, the side chains of the B- and C-rings are crucial to position the chromophore properly in the chromophore pocket of PHYB and for photoreversible spectral changes. The side-chain structure of the D-ring is required for the photoreversible spectral change of the adducts. When methyl and ethyl groups at the C17 and C18 positions are replaced with an n-propyl, n-pentyl, or n-octyl group, respectively, the photoreversible spectral change of the adducts depends on the length of the side chains. From these studies, we conclude that each pyrrole ring of the linear tetrapyrrole chromophore plays a different role in chromophore assembly and the photochromic properties of PhyB.

Crystallographic structures of the ligand-binding domains of the androgen receptor and its T877A mutant complexed with the natural agonist dihydrotestosterone

The structures of the ligand-binding domains (LBD) of the wild-type androgen receptor (AR) and the T877A mutant corresponding to that in LNCaP cells, both bound to dihydrotestosterone, have been refined at 2.0 Å resolution. In contrast to the homodimer seen in the retinoid-X receptor and estrogen receptor LBD structures, the AR LBD is monomeric, possibly because of the extended C terminus of AR, which lies in a groove at the dimerization interface. Binding of the natural ligand dihydrotestosterone by the mutant LBD involves interactions with the same residues as in the wild-type receptor, with the exception of the side chain of threonine 877, which is an alanine residue in the mutant. This structural difference in the binding pocket can explain the ability of the mutant AR found in LNCaP cells (T877A) to accommodate progesterone and other ligands that the wild-type receptor cannot.

Crystal structures of the vitamin D receptor complexed to superagonist 20-epi ligands

The crystal structures of the ligand-binding domain (LBD) of the vitamin D receptor complexed to 1α,25(OH)2D3 and the 20-epi analogs, MC1288 and KH1060, show that the protein conformation is identical, conferring a general character to the observation first made for retinoic acid receptor (RAR) that, for a given LBD, the agonist conformation is unique, the ligands adapting to the binding pocket. In all complexes, the A- to D-ring moieties of the ligands adopt the same conformation and form identical contacts with the protein. Differences are observed only for the 17β-aliphatic chains that adapt their conformation to anchor the 25-hydroxyl group to His-305 and His-397. The inverted geometry of the C20 methyl group induces different paths of the aliphatic chains. The ligands exhibit a low-energy conformation for MC1288 and a more strained conformation for the two others. KH1060 compensates this energy cost by additional contacts. Based on the present data, the explanation of the superagonist effect is to be found in higher stability and longer half-life of the active complex, thereby excluding different conformations of the ligand binding domain.

Mutant G protein α subunit activated by Gβγ: A model for receptor activation?

How receptors catalyze exchange of GTP for GDP bound to the Gα subunit of trimeric G proteins is not known. One proposal is that the receptor uses the G protein's βγ heterodimer as a lever, tilting it to pull open the guanine nucleotide binding pocket of Gα. To test this possibility, we designed a mutant Gα that would bind to βγ in the tilted conformation. To do so, we excised a helical turn (four residues) from the N-terminal region of αs, the α subunit of GS, the stimulatory regulator of adenylyl cyclase. In the presence, but not in the absence, of transiently expressed β1 and γ2, this mutant (αsΔ), markedly stimulated cAMP accumulation. This effect depended on the ability of the coexpressed β protein to interact normally with the lip of the nucleotide binding pocket of αsΔ. We substituted alanine for an aspartate in β1 that binds to a lysine (K206) in the lip of the α subunit's nucleotide binding pocket. Coexpressed with αsΔ and γ2, this mutant, β1-D228A, elevated cAMP much less than did β1-wild type; it did bind to αsΔ normally, however, as indicated by its unimpaired ability to target αsΔ to the plasma membrane. We conclude that βγ can activate αs and that this effect probably involves both a tilt of βγ relative to αs and interaction of β with the lip of the nucleotide binding pocket. We speculate that receptors use a similar mechanism to activate trimeric G proteins.