Structure-Function Relationships of Cellular Retinoic Acid-binding Proteins
QUANTITATIVE ANALYSIS OF THE LIGAND BINDING PROPERTIES OF THE WILD-TYPE PROTEINS AND SITE-DIRECTED MUTANTS*

(Received for publication, August 1, 1996, and in revised form, September 27, 1996)

Lincong Wang , Yue Li and Honggao Yan Dagger

From the Department of Biochemistry, Michigan State University, East Lansing, Michigan 48824

ABSTRACT
INTRODUCTION
EXPERIMENTAL PROCEDURES
RESULTS
DISCUSSION
FOOTNOTES
REFERENCES


ABSTRACT

It has been suggested that electrostatic interactions are critical for binding of retinoic acid by cellular retinoic acid-binding proteins (CRABP-I and CRABP-II). However, the roles of two conserved arginine residues (Arg-111 and Arg-131 in CRABP-I; Arg-111 and Arg-132 in CRABP-II) that interact with the carboxyl group of retinoic acid have not been evaluated. A novel competitive binding assay has been developed for measuring the relative dissociation constants of the site-directed mutants of CRABPs. Arg-111 and Arg-132 of CRABP-II were replaced with methionine by site-directed mutagenesis. The relative dissociation constants of R111M and R132M (Kd (R111M)/Kd (CRABP-II) and Kd (R132M)/Kd(CRABP-II)) were determined to be 40-45 and 6-8, respectively. The ring protons of the aromatic residues of the wild-type CRABP-II and the two mutants were sequentially assigned by two-dimensional homonuclear NMR in conjunction with three-dimensional heteronuclear NMR. Detailed analysis of the nuclear Overhauser effect spectroscopy spectra of the proteins indicated that the conformations of the two mutants are highly similar to that of the wild-type CRABP-II. These results taken together showed that Arg-111 and Arg-132 are important for binding retinoic acid but contribute to the binding energy only by ~2.2 and 1.2 kcal/mol, respectively. In addition, the relative dissociation constant of CRABP-II and CRABP-I (Kd (CRABP-II)/Kd (CRABP-I)) was determined to be 2-3, in close agreement with that calculated using the apparent Kd values determined under the same conditions by fluorometric titrations.


INTRODUCTION

Retinoic acid (RA),1 a hormonally active metabolite of vitamin A, has profound effects on cell growth, differentiation, and morphogenesis. Two types of proteins have been found to bind RA: nuclear retinoic acid receptors (RARs and RXRs) and cellular retinoic acid-binding proteins (CRABPs). RARs and RXRs are RA-activated transcriptional factors that regulate expression of target genes (1). Although the physiological roles of CRABPs are not clear at present, they are thought to be involved in cellular transport and metabolism of RA (2). Two isoforms (CRABP-I and CRABP-II) have been characterized. Both CRABP-I and CRABP-II bind specifically to all-trans-retinoic acid, but they differ in several respects: (i) CRABP-I has higher affinity for RA than CRABP-II (3-6); (ii) CRABP-I is expressed in many adult mouse and human tissues, but the expression of CRABP-II is limited to skin (7, 8); and (iii) RA stimulates expression of CRABP-II but not that of CRABP-I (8). It appears that the two isoforms may have distinct functions. The idea is supported by the fact that the sequence identity of human and mouse CRABP-I (99.3%) or human and mouse CRABP-II (93.5%) is much higher than the sequence identity (73.7%) between the two isoforms from the same source.

CRABPs are members of a family of intracellular lipid-binding proteins that bind small hydrophobic molecules such as retinoids and fatty acids (9). The family of proteins includes CRABPs, cellular retinol-binding proteins, fatty acid-binding proteins, P2 myelin protein, a mammary gland protein, and gastropin. The structures of 11 different intracellular lipid-binding proteins, including CRABP-I and CRABP-II, have been determined by x-ray crystallography (9-11). Although the sequence identity among intracellular lipid-binding proteins is rather low (~20%), these crystal structures are remarkably similar with respect to backbone folding. They are composed of two nearly orthogonally packed five-stranded beta -sheets and two short alpha -helices. The ligand binding pocket in each protein is deep inside the interior of the beta -barrel formed by the two beta -sheets. The helix-turn-helix motif is located at the ligand entrance. These proteins apparently lack a true hydrophobic core that is important for folding and stability of many other proteins (12).

Electrostatic interactions are thought to play major roles in binding of RA by CRABPs because the two proteins bind only retinoic acid but not retinol and retinal. Two arginine residues (Arg-111 and Arg-131 in CRABP-I; Arg-111 and Arg-132 in CRABP-II) have been identified by crystallography to interact with the carboxyl group of the bound RA (10). However, previous site-directed mutagenesis studies suggest that both Arg-111 and Arg-131 are critical for CRABP-I to bind RA (13) but that only Arg-132 is important for CRABP-II to bind RA (14). Because of the limitations of the current methods for measuring RA binding, the dissociation constants of the critical mutants have not been measured. Furthermore, the conformations of the CRABP-I mutants have been characterized only by CD, and the conformations of the CRABP-II mutants have not been characterized. Thus, the contributions of the electrostatic interactions between the arginine residues and RA to binding are still uncertain.

We are interested in studying the quantitative structure-function relationships of CRABPs with respect to binding of RA. In this paper, we report a novel competitive binding assay developed for measuring the dissociation constants of the site-directed mutants of CRABPs. We have used the novel method to evaluate the contributions of both Arg-111 and Arg-132 of CRABP-II to binding of RA in conjunction with site-directed mutagenesis and NMR. The results show that like CRABP-I, both Arg-111 and Arg-132 in CRABP-II are important for binding of RA, contrary to the results of Chen et al. (14). However, Arg-111 and Arg-132 contribute to the overall binding energy only by ~2.2 and 1.2 kcal/mol, respectively.


EXPERIMENTAL PROCEDURES

Materials

Nonradioactive RA was purchased from Sigma. [11,12-3H]RA was purchased from DuPont NEN. A DNA sequencing kit was obtained from U.S. Biochemical Corp. Enzymes for recombinant DNA experiments were purchased from Life Technologies, Inc., or New England Biolabs. Other chemicals were analytical or reagent grade from commercial sources.

Cloning and Expression

The cDNAs encoding human CRABP-I and CRABP-II were kindly provided by Dr. Anders Åström (8). The genes were subcloned into the expression vector pET-17b (Novagen) by polymerase chain reaction. The primers for polymerase chain reaction were: 5'-GGAATTCCATATGCCCAACTTCGCCGGC-3' (forward, CRABP-I), 5'-CGGGATCCTCATTCCCGGACATAAATTC-3' (reverse, CRABP-I), 5'-CGGGATCCATATGCCAAACTTCTCTGGCAACTG-3' (forward, CRABP-II), and 5'-GGAATTCTCACTCTCGGACGTAGAC-3' (reverse, CRABP-II). To ensure that there were no undesired mutations in the amplified genes, they were sequenced by double-stranded DNA sequencing from both orientations. The expression constructs were then transformed into the bacterial strain BL21(DE3)pLysS (15). The transformed bacterial cells were grown at 37 °C in a LB agar plate containing both ampicillin (100 µg/ml) and chloramphenicol (25 µg/ml). Expression of CRABPs was verified by growing the colonies in 5 ml of LB medium followed by induction with isopropyl-1-thio-beta -D-galactopyranoside and SDS-polyacrylamide gel electrophoresis of the harvested cells.

Site-directed Mutagenesis

The oligonucleotides for making human CRABP-II mutants were 5'-TCGTGGACCATGGAACTGACCA-3' (R111M) and 5'-GTGTGCACCATGGTCTACGTC-3' (R132M). The mutants were generated by the method of Kunkel (16) and screened by DNA sequencing. In order to ensure that there were no unintended mutations in the mutants, the entire sequences of the mutated genes were determined.

Protein Purification

All proteins were purified by the same procedure as described below. 100 ml of LB medium containing 100 µg/ml of ampicillin and 25 µg/ml of chloramphenicol was inoculated by a small piece of a frozen seed culture and incubated at 37 °C with vigorous shaking (200 rpm) overnight. It was then used to inoculate 4 liters of LB medium containing both antibiotics and grown at 37 °C. When A600 of the culture reached 0.6-0.8, isopropyl-1-thio-beta -D-galactopyranoside was added to a final concentration of 0.4 mM. The culture was allowed to grow further at 30 °C for 3-5 h. The bacterial cells were then harvested by centrifugation and suspended in 100 ml of buffer A (10 mM Tris-HCl, 1 mM dithiothreitol, pH 8.2). The suspension was sonicated on ice and centrifuged (27,000 × g) at 4 °C for 30 min. The supernatant was applied to a DEAE-cellulose column equilibrated with the same buffer. The column was washed with buffer A until A280 of the eluent was less than 0.05. Elution of the column was achieved by a linear NaCl gradient (0-200 mM in buffer A) and monitored by A280 and 15% SDS-polyacrylamide gel electrophoresis. The fractions containing CRABP were pooled and concentrated by an Amicon ultrafiltration cell using a YM 10 membrane. The concentrated protein solution was centrifuged (27,000 × g) for 20 min. The supernatant was loaded onto a Sephadex G-50 column equilibrated with phosphate-buffered saline (4 mM NaH2PO4, 16 mM Na2HPO4, 150 mM NaCl, pH 7.3) and eluted with the same buffer. The fractions containing CRABP were located by A280 and 15% SDS-polyacrylamide gel electrophoresis. The fractions of greater than 99% purity were pooled and concentrated. The protein solution was dialyzed against double-distilled water and lyophilized.

Fluorometric Titration

Fluorescence binding assays were carried out by a procedure modified from that of Cogan et al. (17) using a Hitachi 4500 fluorometer. Briefly, CRABPs were dissolved in phosphate-buffered saline. The concentrations of the protein stock solutions were measured by A280. The absorption coefficients used were 19,868 M-1 cm-1 for CRABP-I and 19,480 M-1 cm-1 for CRABP-II determined by the method of Gill and Hippel (18). RA stock solutions were prepared in absolute ethanol. The concentrations of the RA stock solutions were determined by A336 using the absorption coefficient of 45,000 M-1 cm-1. The final ethanol content for each titration was kept less than 2%. The samples were excited at wavelength of 283 nm with a slit width of 2.5 nm. The emission wavelengths were 330 nm for CRABP-I and 338 nm for CRABP-II with a slit width of 20 nm. The excitation shutter was closed between measurements. The inner filter effects were not corrected because they were negligible at the protein and RA concentrations we used in the experiments. The data were analyzed by nonlinear least square fit to Equation 1.
<FR><NU>F</NU><DE>F<SUB>f</SUB></DE></FR>=1+<FENCE><FR><NU>F<SUB>b</SUB></NU><DE>F<SUB>f</SUB></DE></FR>−1</FENCE>×<FENCE><FR><NU>P<SUB>t</SUB>+R<SUB>t</SUB>+K<SUB>d</SUB>−<RAD><RCD>(P<SUB>t</SUB>+R<SUB>t</SUB>+K<SUB>d</SUB>)<SUP>2</SUP>−4P<SUB>t</SUB>R<SUB>t</SUB></RCD></RAD></NU><DE>2P<SUB>t</SUB></DE></FR></FENCE> (Eq. 1)
where F is the observed fluorescence, Fb is the fluorescence of the bound CRABP, Ff is the fluorescence of the free CRABP, Pt is the total protein concentration, and Rt is the total RA concentration.

Competitive Binding Assay

The assay was carried out in a Spectra/Por equilibrium dialyzer at room temperature. The two compartments of each dialysis cell were separated by a semipermeable membrane with a molecular mass cut-off of 6-8 kDa. One compartment was filled with the wild-type CRABP-II (1 ml), and the other was filled with a CRABP-II mutant or the wild-type CRABP-I. The buffer used was phosphate-buffered saline plus 5 mM dithiothreitol. An equal amount of [3H]RA (100 nM) was added to each compartment. The protein concentrations in both compartments were much greater (at least 20-fold) than the RA concentration. 100 µl of samples was taken from the two compartments after various times of incubation at room temperature, mixed with 5 ml of scintillation fluids, and counted by a liquid scintillation counter. The equilibria in the two compartments that contained the wild-type CRABP-II and mutant proteins can be described by Equations 2 and 3 (similar equations can be written for CRABP-I and CRABP-II),
K<SUB>d(<UP>WT</UP>)</SUB>=<FR><NU>[<UP>WT</UP>][<UP>RA</UP>]</NU><DE>[<UP>WT</UP> · <UP>RA</UP>]</DE></FR> (Eq. 2)
K<SUB>d(<UP>MT</UP>)</SUB>=<FR><NU>[<UP>MT</UP>][<UP>RA</UP>]</NU><DE>[<UP>MT</UP> · <UP>RA</UP>]</DE></FR> (Eq. 3)
where WT, MT, WT·RA, and MT·RA represent the wild-type CRABP-II, a CRABP-II mutant, and their RA complexes, respectively. Therefore,
<FR><NU>K<SUB>d(<UP>MT</UP>)</SUB></NU><DE>K<SUB>d(<UP>WT</UP>)</SUB></DE></FR>=<FR><NU>[<UP>MT</UP>][<UP>WT</UP> · <UP>RA</UP>]</NU><DE>[<UP>WT</UP>][<UP>MT</UP> · <UP>RA</UP>]</DE></FR> (Eq. 4)
Since the concentrations of the proteins were much greater than their respective dissociation constants and the concentration of RA, [WT] approx  [WT]total, [MT] approx  [MT]total, [WT·RA] >>  [RA], and [MT·RA] >>  [RA]. Then the relative dissociation constant can be calculated by Equation 5,
<FR><NU>K<SUB>d(<UP>MT</UP>)</SUB></NU><DE>K<SUB>d(<UP>WT</UP>)</SUB></DE></FR>=<FR><NU>[<UP>MT</UP>]<SUB><UP>total</UP></SUB>C<SUB><UP>WT</UP></SUB></NU><DE>[<UP>WT</UP>]<SUB><UP>total</UP></SUB>C<SUB><UP>MT</UP></SUB></DE></FR> (Eq. 5)
where CWT and CMT are the measured radioactivities of the two compartments containing the wild-type CRABP-II and the mutant, respectively. It turned out that the system could not reach equilibrium in 2 days, presumably because of few free RAs in solution to diffuse across the membrane. Since RA is not stable even in the dark, the assay was redesigned to match the equilibrium conditions by varying the ratio of the protein concentrations of the wild-type and the mutant ([MT]total/[WT]total). Thus, the concentration of the mutant was varied while keeping the concentration of the wild type at ~2 µM. Initially, the concentration of the mutant was increased in an exponential manner (e.g. 2, 20, 200 µM). Then it was varied in a small range. Since an equal amount of RA was added to the two compartments of the dialysis cell, the two compartments should have the same RA concentration and radioactivity at the beginning of the assay. If [MT]total/[WT]total not equal  Kd (MT)/Kd (WT), there would have a net transfer of RA across the semipermeable membrane separating the two compartments. Thus, the radioactivity counts of the two compartments (CWT and CMT) would differ after incubation for a certain period. When [MT]total/[WT]total < Kd (MT)/Kd (WT), then CWT - CMT > 0. When [MT]total/[WT]total > Kd (MT)/Kd (WT), then CWT - CMT < 0. When [MT]total/[WT]total = Kd (MT)/Kd (WT), then CWT - CMT = 0.

NMR Spectroscopy

All NMR measurements were performed at 32 °C on a Varian VXR-500 spectrometer operating at a proton frequency of 500 MHz. The proteins were dissolved in 20 mM sodium phosphate, pH 7.5 (direct pH meter reading), 100 mM NaCl, 5 mM dithiothreitol in D2O. The protein concentration was ~2 mM. The spectral width was 7200 Hz. For wild-type CRABP-II, one phase-sensitive DQF-COSY spectrum (19), two MLEV-17 clean TOCSY spectra (20-22) with mixing times of 20 and 40 ms and one NOESY spectrum (23, 24) with a mixing time of 150 ms were acquired. For R111M, one DQF-COSY spectrum, three TOCSY spectra with mixing times of 20, 40, and 75 ms and one NOESY spectrum with a mixing time of 150 ms were acquired. For R132M, one TOCSY spectrum with a mixing time of 40 ms and one NOESY spectrum with a mixing time of 150 ms were acquired. All of the spectra were acquired in the hypercomplex mode with standard phase cycling schemes. The data were usually acquired with 2048 complex points in the t2 dimension and 256 complex points in the t1 dimension. 96 transients were collected for each FID. Data processing was performed on a Sun Sparc 10 station using VNMR software from Varian. The time domain data were zero-filled once and multiplied by shifted sine bell or Gaussian functions before Fourier transformation in both dimensions. Base lines were corrected in t2 dimensions using a 5-order polynomial. Chemical shifts were referenced to internal sodium 3-(trimethyl silyl)-propionate-2,2,3,3-d4.


RESULTS

Fluorometric Titration

Probably because of the simplicity of the method, fluorometric titration has been frequently used for measuring the affinities of CRABPs for RA. However, there are large discrepancies in the Kd values obtained by the method. These could be due to variations in the assay methodology and/or conditions. Some problems with the method have been recently discussed (2, 4, 6). In order to compare the ligand binding properties of CRABPs under the same conditions, we carried out fluorometric titration experiments. The results are shown in Fig. 1. The apparent Kd values obtained by nonlinear fitting were 0.63 nM for CRABP-I and 1.9 nM for CRABP-II. While our work was in progress, Norris et al. (6) reported the apparent Kd values of mouse CRABPs measured by an improved fluorescence titration method. The major changes were the lowering of protein and RA concentrations and the use of nonlinear least square fitting for data analysis. These modifications were in line with our approach. Our fluorometric titration data were analyzed by nonlinear least square fitting. The protein concentrations we used in the assays (~50 nM) were the lowest among those used in previous studies except those of Norris et al. (6). We did not include the additive gelatin used by Norris et al. (6) in the assays. The apparent Kd values determined by us are in close agreement with those (0.4 nM for CRABP-I and 2 nM for CRABP-II) reported by Norris et al. (6).


Fig. 1. Fluorometric titration of wild-type CRABP-I (A) and CRABP-II (B). The concentrations of CRABP-I and CRABP-II were 40 and 50 nM, respectively.
[View Larger Version of this Image (15K GIF file)]


Competitive Binding Assay

In order to measure the affinities of the mutants of CRABPs for RA, we developed a novel competitive binding assay. We first used the method to measure the relative affinity of CRABP-I and CRABP-II for RA. The result is shown in Fig. 2A. When the ratio of the concentrations of the two proteins ([CRABP-II]/[CRABP-I]) was <= 2, there was a net transfer of RA from the compartment containing CRABP-II to the compartment containing CRABP-I. When the ratio of the concentrations of the two proteins was >= 3, there was a net transfer of RA from the compartment containing CRABP-I to the compartment containing CRABP-II. As described under "Experimental Procedures," the relative Kd of the two proteins lies between the points with opposite net transfers. Thus the Kd of CRABP-II relative to that of CRABP-I (Kd(CRABP-II)/Kd(CRABP-I)) is 2-3. It is in close agreement with the results of the fluorometric titration studies. We then used the competitive binding assay to measure the Kd values of R111M and R132M relative to that of the wild-type CRABP-II. The results are shown in Fig. 2 (B and C). In the case of R111M (Fig. 2B), there was a net transfer of RA from the compartment containing R111M to the compartment containing the wild-type CRABP-II when [R111M]/[WT] <=  40. When [R111M]/[WT] >=  45, there was a net transfer of RA from the compartment containing the wild-type protein to the compartment containing R111M. Therefore, the relative affinity of R111M for RA (Kd (R111M)/Kd (WT)) is 40-45. As shown in Fig. 2C, the relative affinity of R132M for RA (Kd (R132M)/Kd (WT)) is 6-8.


Fig. 2. Competitive binding assays for measuring the relative dissociation constants of CRABPs. The relative radioactivity in A is the radioactivity count of the compartment containing CRABP-I minus that of the compartment containing CRABP-II. The relative radioactivity in B is the radioactivity count of the compartment containing CRABP-II minus that of the compartment containing R111M. The relative radioactivity in C is the radioactivity count of the compartment containing CRABP-II minus that of the compartment containing R132M.
[View Larger Version of this Image (17K GIF file)]


Sequential Resonance Assignment of the Ring Protons of the Aromatic Residues of the Wild-type CRABP-II and Site-directed Mutants

Since perturbations in the ligand binding property of a mutant may be due to conformational changes caused by the mutation (25), we characterized the conformations of the mutants by NMR. A prerequisite for NMR structural analysis is sequential resonance assignment. Currently, we are in the process of making total sequential resonance assignments of the wild-type CRABP-II and R111M by multidimensional multinuclear NMR. We have sequentially assigned the ring protons of aromatic residues of the wild-type and mutants based on the homonuclear two-dimensional NMR (DQF-COSY, TOCSY, and NOESY) and heteronuclear three-dimensional NMR.2 The chemical shifts of the ring protons are listed in Table I. The sequential resonance assignment is labeled in Fig. 3 for R132M. In comparison with the chemical shifts of the wild-type ring protons, there are quite a few chemical shift changes in the ring protons of the mutants. Some chemical shift changes can be easily rationalized, while the causes for other chemical shift changes are not obvious. For example, the guanidino group of Arg-132 stacks on the aromatic ring of Phe-15 and is hydrogen-bonded to the hydroxyl group of Tyr-134 in the crystal structure of the wild-type CRABP-II (10). Removal of the guanidino group is likely to be the cause for the chemical shift changes in Phe-15 and Tyr-134 of R132M. The crystal structure also reveals that the guanidino group of Arg-111 interacts with Trp-109 through an amido/aromatic-ring hydrogen bond (10). Thus, removal of the guanidino group of Arg-111 may cause the chemical shift changes in Trp-109 of R111M. Some small chemical shift differences may result from changes in pH, temperature, and other conditions. They are not necessarily indicative of conformational changes. Conformational changes are best characterized by analysis of NOESY spectra as described below. Interestingly, 4H and 7H of the Trp-87 ring protons in R132M have two sets of NMR signals (Fig. 3), indicating that there is localized conformational heterogeneity and that the Trp ring is in two distinct conformational states. The lifetimes of the two conformational states are at least 50 ms. The ratio of the populations in the two states (major and minor forms) is about 4:1. 1H of the Trp-87 ring protons in the wild-type protein and R111M also show two distinct resonance signals (data not shown), indicating that the two proteins also have localized conformational heterogeneities. The lifetimes of the two conformational states are at least 40 ms in the wild-type protein and at least 33 ms in R111M. The ratio of the major and minor forms is about 3:1 in the wild-type protein and about 4:1 in R111M. Conformational dynamics may be important for RA to get into the deep binding pocket of the protein. However, Trp-87 is not at the RA-binding pocket and is far away from the ligand entrance. The functional significance of the conformational heterogeneity is not clear.

Table I.

Chemical shifts of the ring protons of the aromatic residues of wild-type CRABP II, R111M mutant, and R132M mutant

The underlined values are the mutant resonances that differ by >0.02 ppm from the corresponding resonances of the wild-type CRABP-II.
Residue Wild-type CRABP-II R111M R132M

Phe-3 7.60a 7.34b 7.52c 7.61a  <UNL>7.38</UNL>b  <UNL>7.45</UNL>c 7.60a  <UNL>7.34</UNL>b 7.52c
Phe-15 7.42a 7.02b 7.11c  <UNL>7.39</UNL>a  <UNL>7.12</UNL>b  <UNL>7.20</UNL>c 7.40a  <UNL>7.11</UNL>b 7.11c
Phe-50 6.66a 7.26b 6.99c 6.68a 7.24b 6.99c 6.64a 7.24b 6.99c
Phe-65 6.42a 6.70b 6.48c  <UNL>6.38</UNL>a 6.70b  <UNL>6.44</UNL>c 6.42a 6.69b  <UNL>6.42</UNL>c
Phe-71 7.48a 7.00b 6.70c  <UNL>7.44</UNL>a 6.99b  <UNL>6.73</UNL>c 7.48a 7.00b 6.69c
Tyr-51 7.12d 6.79e 7.12d  <UNL>6.74</UNL>e 7.13d 6.78e
Tyr-134 6.98d 6.65e 7.00d  <UNL>6.86</UNL>e  <UNL>7.03</UNL>d  <UNL>6.75</UNL>e
Trp-7 7.25f 7.00g 7.41h 7.30i  <UNL>7.22</UNL>f 6.99g 7.39h 7.29i 7.25f 6.99g 7.43h 7.30i
Trp-87 8.16f 7.13g 7.39h 7.80i 8.16f  <UNL>7.04</UNL>g 7.39h  <UNL>7.76</UNL>i 8.17f  <UNL>7.05</UNL>g 7.39h 7.80i
Trp-87'j 8.07f 7.05g 7.38h 7.76i
Trp-109 7.08f 6.53g 6.61h 6.87i  <UNL>7.01</UNL>f  <UNL>6.59</UNL>g 6.59h 6.89i 7.09f 6.52g 6.61h 6.87i

a-c 2,6H, 3,5H, and 4H of phenylalanine, respectively.
d,e 2,6H and 3,5H of tyrosine, respectively.
f-i 4H, 5H, 6H, and 7H of tryptophan, respectively.
j  Trp-87 has two sets of resonances because of conformational heterogeneity.


Fig. 3. TOCSY spectrum of the aromatic protons of R132M at 32 °C showing two sets of resonances for Trp-87. The mixing time was 40 ms. All of the cross-peaks have been sequentially assigned. The chemical shifts of the resonances are listed in Table I.
[View Larger Version of this Image (44K GIF file)]


Conformational Comparison of the Wild-type CRABP-II, R111M, and R132M

The aromatic-aromatic NOEs are shown in the lower panels of Fig. 4 (wild-type CRABP-II), Fig. 5 (R111M), and Fig. 6 (R132M). Some aromatic-aliphatic NOEs are shown in the upper panels of the figures. Interresidue NOEs are labeled, because they are indicative of the conformation of a protein. The interresidue NOE peaks of the three proteins are the same except for some minor variations in peak intensity. The observed NOEs are also consistent with the atomic distances measured from the coordinate of the crystal structure of CRABP-II complexed with RA. The identities of the aliphatic protons shown in the upper panels of these figures are not known at present. However, qualitatively, the three proteins have very similar aromatic-aliphatic NOE patterns. The results of the NOESY analysis suggest that the conformations of R111M and R132M are highly similar to that of the wild-type protein. Thus, the decreases in affinity for RA of the two mutants are most likely due to the disruption of the interactions of the arginine residues with RA by mutagenesis. Arg-111 and Arg-132 contribute to the binding energy by ~2.2 and 1.2 kcal/mol, respectively.


Fig. 4. Parts of the 500-MHz NOESY spectrum of the wild-type CRABP-II at 32 °C. The mixing time was 150 ms. Only the interresidue NOEs are labeled. The identities of the NOEs are as follows. A, Phe-65 2,6H to Phe-71 2,6H; B, Phe-65 2,6H to Phe-71 3,5H; C, Phe-65 4H to Trp-109 6H; D, Phe-65 3,5H to Trp-109 6H; E, Phe-50 2,6H to Trp-87 7H; F, Phe-71 4H to Trp-109 7H; G, Phe-65 4H to Trp-109 7H; H, Phe-65 2,6H to Trp-109 7H; I, Phe-50 3,5H to Phe-65 4H; J, Trp-87 5H to Phe-3, 4H; K, Phe-50 3,5H to Phe-3 2,6H; L, Phe-50 3,5H to Trp-87 7H; M, Phe-50 4H to Trp-87 7H; N, Phe-50 4H to Phe-3 2,6H; O, Phe-50 2,6H to Trp-87 6H; P, Phe-50 2,6H to F3-2,6H; Q, Trp-87 5H to Phe-3 2,6H; R, Phe-3 4H to Trp-87 4H.
[View Larger Version of this Image (34K GIF file)]



Fig. 5. Parts of the 500-MHz NOESY spectrum of R111M at 32 °C. The mixing time was 150 ms. Only the interresidue NOEs are labeled. The identities of the NOEs are given in the legend to Fig. 4.
[View Larger Version of this Image (35K GIF file)]



Fig. 6. Parts of the 500-MHz NOESY spectrum of R132M at 32 °C. The mixing time was 150 ms. Only the interresidue NOEs are labeled. The identities of the NOEs are given in the legend to Fig. 4.
[View Larger Version of this Image (32K GIF file)]



DISCUSSION

Competitive Binding Assay

Two types of methods have been in general use for measuring binding of RA to CRABPs: fluorometry and radiometry. The radiometric method involves separation of bound from free RA by dextran-coated charcoal, gel filtration, and other means. Substantial loss of bound ligand during the separation process makes the method unsuitable for measuring the dissociation constants of site-directed mutants with greatly decreased affinity for RA. The very limited solubility of RA in water (~200 nM) (26) also makes the fluorometric method inapplicable for determining the dissociation constants of these mutants. Studies of the quantitative structure-function relationships of CRABPs have been hampered by the lack of methods for measuring the affinities of site-directed mutants for RA (13, 14). We have developed a novel competitive binding assay for measuring the affinities of site-directed mutants for RA relative to that of the wild-type CRABP. The essence of the method is to monitor the competition between a mutant and the wild-type protein for binding of limited RA. The method uses an equilibrium dialyzer. The two compartments of each dialysis cell are filled with the wild-type and mutant proteins, respectively. The absolute concentration of RA is not important as long as the concentration of free RA is much smaller than that of bound RA. There is no need to separate bound from free RA by dextran-coated charcoal, gel filtration, and other means. The transfer of RA from one compartment to the other is determined by measuring the radioactivities of the samples taken from the two compartments. The direction of the net transfer is dependent on the relative affinity of the proteins and the ratio of the protein concentrations of the two compartments. If the solubility of a mutant is the same as that of the wild-type CRABP-II (~2 mM), then the ratio of the concentrations of the wild-type and mutant proteins can be varied by more than 1,000-fold, and a relative Kd of as large as 1000 (Kd (MT)/Kd (WT) = 1000) can still be measured by the method. More importantly, determination of the relative dissociation constant of a point mutant is sufficient for estimating the energetic contribution of the amino acid residue to ligand binding (Delta Delta G = RTln(Kd (MT)/Kd (WT))).

In addition to vitamin A, many other hydrophobic molecules such as vitamin D and steroid hormones play vital roles in a variety of cellular processes. Studies of the structure-function relationships of the receptor proteins for these bioactive hydrophobic molecules are of fundamental interest and biomedical significance. The method can also be used for measuring the relative dissociation constants of the site-directed mutants of these important proteins and may facilitate these studies.

Relative Affinity of CRABP-I and CRABP-II for RA

CRABP-I and CRABP-II differ in affinity for RA, expression pattern, and regulation. Although it is generally accepted that CRABP-I has a higher affinity for RA than CRABP-II, it has been difficult to determine precisely the relative affinity of the two proteins. A wide range of apparent Kd values have been reported. The apparent Kd values for purified CRABP-I are in the range of 0.4-39 nM (4-6, 13, 27-29). The apparent Kd values for purified CRABP-II are in the range of 2-64 nM (3-6, 14). The best estimation of the relative affinity of the two proteins for RA was made by Fiorella et al. (5). In their study, RA was first incubated with CRABP-II and then mixed with CRABP-I. The protein·RA complexes (CRABP-I·RA and CRABP-II·RA) were separated by chromatofocusing. The result indicated that CRABP-II has about 3-fold lower affinity for RA than CRABP-I. More recently, Norris et al. (6) compared the affinities of CRABP-I and CRABP-II for RA by fluorometric titration under the same conditions. The results showed that CRABP-I has 4-fold higher affinity for RA than CRABP-II (although a 10-fold stronger affinity of CRABP-I for RA, compared with CRABP-II, was stated in the paper, apparently due to an error in calculation). The apparent dissociation constant of CRABP-II relative to that of CRABP-I (Kd (CRABP-II)/Kd(CRABP-I)) is 3, based on our fluorometric measurements. The results of our competitive binding assays indicated that the relative apparent dissociation constant of CRABP-II (Kd (CRABP-II)/Kd(CRABP-I)) is 2-3. The results of our fluorometric measurements and competitive binding studies are in close agreement. They are also in line with the results of Fiorella et al. (5). The relative affinity of CRABP-I and CRABP-II measured by Fiorella et al. (5) may be the upper limit, because some RA bound to CRABP-II could be preferentially lost during separation of the protein·RA complexes by chromatofocusing. The difference in affinity for RA between CRABP-I and CRABP-II is rather small (1-2-fold).

Both Arg-111 and Arg-132 of CRABP-II Are Important for Binding of RA

Arg-111, Arg-132, and Tyr-134 of CRABP-II have been predicted to be involved in binding of RA on the basis of homologous modeling (13, 30, 31). The crystal structure of CRABP-II in complex with RA reveals that one of the carboxylate oxygens of RA is hydrogen-bonded to the guanidino group Arg-132 and the hydroxyl group of Tyr-134, and the other carboxylate oxygen of RA interacts with Arg-111 via a water molecule (10). We have substituted Arg-111 and Arg-132 with methionine. Our biochemical and structural characterizations show that guanidino groups of Arg-111 and Arg-132 contribute to the binding energy by ~2.2 and 1.2 kcal/mol, respectively. The results indicate that both Arg-111 and Arg-132 are indeed involved in binding of RA and that Arg-111 contributes more to the binding energy than Arg-132. Fluorometric titration also suggests that R132M has higher affinity for RA than R111M (data not shown). Zhang et al. (13) have replaced the corresponding residues in CRABP-I with glutamine (R111Q and R131Q). Quantitative comparison of our results with those of Zhang et al. (13) is not possible, because the dissociation constants of their mutants have not been determined. Qualitatively, the ligand binding properties of our mutants are similar to those of the corresponding CRABP-I mutants. Zhang et al. (13) observed that both the wild-type CRABP-I and R131Q can induce CD of RA in the wavelength range between 320 and 360 nm, but R111Q cannot induce the CD signal. The results of the CD experiments with CRABP-I are consistent with our measurements of the energetic contributions of Arg-111 and Arg-132 to the binding of RA. Chen et al. (14) have mutated Arg-111 of CRABP-II to alanine (R111A) and Arg-132 of CRABP-II to alanine and glutamine (R132A and R132Q). In contrast to our R111M mutant, the apparent dissociation constant of their R111A mutant is very similar to that of the wild-type CRABP-II (only 50% higher than that of the wild type). The causes of the discrepancy are not clear, but several possibilities can be ruled out. (i) It is unlikely that there are mutations other than substitution of Arg-111 with a methionine in our R111M mutant. We have sequenced the entire gene after mutagenesis and subcloning into the expression vector. (ii) Steric hindrance is unlikely to be the cause, because the side chain of methionine is smaller than that of arginine. (iii) Although small conformational changes cannot be ruled out, NMR characterizations suggest that R111M is properly folded, and its conformation is very similar to that of the wild-type protein. Furthermore, R111M has been crystallized. The preliminary crystal structure at 2.2-Å resolution shows that the structure of R111M is indeed very similar to that of the wild-type protein.3 The ligand binding property of our R132M mutant is qualitatively similar to those of the R132A and R132Q generated by Chen et al. (14). However, quantitative comparison is again not possible, because the apparent dissociation constants of the two mutants have not been determined.

The results of our mutagenesis studies appear to be partly inconsistent with the crystallographic data. The crystal structure of holo-CRABP-II suggests that Arg-132 interacts directly with the carboxyl group of RA, whereas the interaction between Arg-111 and the carboxyl group of RA is mediated by a water molecule (10). Our results agree with the crystallographic observation that both Arg-111 and Arg-132 are involved in binding of RA. However, in contrast to the crystallographic data, our mutagenesis results indicate that Arg-111 is more important than Arg-132 for RA binding. But this is not necessarily inconsistent with the crystallographic studies. Since Arg-132 and Tyr-134 are hydrogen-bonded to the same carboxylate oxygen of RA, the loss of the hydrogen bond between Arg-132 and RA in R132M may be partially compensated by the interaction between Tyr-134 and RA. On the other hand, the other carboxylate oxygen of RA only interacts with Arg-111 (albeit via a water molecule), the loss of the interaction in R111M may not be compensated. Alternatively, in solution, the position of the carboxyl group of the bound RA may be slightly different from that observed in the crystals. It is noted that the RA binding cavity is much larger than necessary to accommodate the ligand, and a large portion of the deep binding cavity is not occupied by RA. It appears that the beta -ionone ring of the bound RA is well fixed, but the isoprene tail and carboxyl group have room to move. The guanidino group of Arg-111 could be closer to the carboxyl group of RA than to that of Arg-132 in solution. Solution NMR study of holo-CRABP-II is in progress.

CRABPs have very stringent retinoid specificity (2). The proteins bind only RA and reject both retinol and retinal. On the other hand, cellular retinol-binding proteins bind both retinol and retinal and reject RA. Amino acid sequence analysis reveals that the two conserved arginine residues in CRABPs are replaced with glutamine in cellular retinol-binding proteins. It has long been thought that the electrostatic interactions between the guanidino groups of the two conserved arginine residues and the carboxyl group of RA play major roles in binding of RA (31). Our mutagenesis studies indicate that Arg-111 and Arg-132 are indeed involved in RA binding, but the electrostatic interactions between the guanidino groups of Arg-111 and Arg-132 and the carboxyl group of RA contribute to the overall binding energy only by ~2.2 and 1.2 kcal/mol. Thus, amino acid residues that interact with the hydrophobic moiety of RA may be more important for binding of RA.


FOOTNOTES

*   This work was supported by funds from the REF Center of Protein Structure and Design and the Cancer Center at Michigan State University. The costs of publication of this article were defrayed in part by the payment of page charges. The article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
Dagger    To whom correspondence should be addressed. Tel.: 517-353-8786; Fax: 517-353-9334; E-mail: yan{at}nmr1.bch.msu.edu.
1    The abbreviations used are: RA, all-trans-retinoic acid; CRABP, cellular retinoic acid-binding protein; DQF-COSY, double quantum-filtered correlation spectroscopy; NOE, nuclear Overhauser effect; NOESY, nuclear Overhauser effect spectroscopy; RAR, retinoic acid receptor; RXR, retinoid X receptor; TOCSY, total correlation spectroscopy; WT, wild-type CRABP-II; MT, CRABP-II mutant.
2    L. Wan, Y. Li, and H. Yan, unpublished results.
3    Ji, X., L. Wang, Y. Li, and H. Yan, unpublished results.

REFERENCES

  1. Mangelsdorf, D. J., Umesono, K., and Evans, R. M. (1994) in The Retinoids: Biology, Chemistry, and Medicine (Sporn, M. B., Rorberts, A. B., and Goodman, D. S., eds), pp. 319-349, Raven Press, New York
  2. Ong, D. E., Newcomer, M. E., and Chytil, F. (1994) in The Retinoids: Biology, Chemistry, and Medicine (Sporn, M. B., Rorberts, A. B., and Goodman, D. S., eds), pp. 283-317, Raven Press, New York
  3. Bailey, J. S., and Siu, C.-H. (1988) J. Biol. Chem. 263, 9326-9332 [Abstract/Free Full Text]
  4. Fogh, K., Voorhees, J. J., and Åström (1993) Arch. Biochem. Biophys. 300, 751-755 [CrossRef][Medline] [Order article via Infotrieve]
  5. Fiorella, P. D., Giguère, V., and Napoli, J. L. (1993) J. Biol. Chem. 268, 21545-21552 [Abstract/Free Full Text]
  6. Norris, A. W., Cheng, L., Giguère, V., Rosenberger, M., and Li, E. (1994) Biochim. Biophys. Acta 1209, 10-18 [Medline] [Order article via Infotrieve]
  7. Giguère, V., Lyn, S., Yip, P., Siu, C.-H., and Amin, S. (1990) Proc. Natl. Acad. Sci. U. S. A. 87, 6233-6237 [Abstract]
  8. Åström, A., Tavakkol, A., Pettersson, U., Cromie, M., Elder, J. T., and Voorhees, J. J. (1991) J. Biol. Chem. 266, 17662-17666 [Abstract/Free Full Text]
  9. Banaszak, L., Winter, N., Xu, Z., Bernlohr, D. A., Cowan, S. W., and Jones, T. A. (1994) Adv. Protein Chem. 45, 89-151 [Medline] [Order article via Infotrieve]
  10. Kleywegt, G. J., Bergfors, T., Senn, H., Le Motte, P., Gsell, B., Shudo, K., and Jones, T. A. (1994) Structure 2, 1241-1258 [Medline] [Order article via Infotrieve]
  11. Thompson, J. R., Bratt, J. M., and Banaszak, J. (1995) J. Mol. Biol. 252, 433-446 [CrossRef][Medline] [Order article via Infotrieve]
  12. Scapin, G., Gordon, J. I., and Sacchettini, J. C. (1992) J. Biol. Chem. 267, 4253-4269 [Abstract/Free Full Text]
  13. Zhang, J., Liu, Z.-P., Jones, T. A., Gierasch, L. M., and Sambrook, J. F. (1992) Proteins Struct. Funct. Genet. 13, 87-89 [Medline] [Order article via Infotrieve]
  14. Chen, L. X., Zhang, Z.-P., Scafonas, A., Cavalli, R. C., Gabriel, J. L., Soprano, K. J., and Soprano, D. R. (1995) J. Biol. Chem. 270, 4518-4525 [Abstract/Free Full Text]
  15. Studier, F. W., Rosenberg, A. H., Dunn, J. J., and Dubendorff, J. W. (1987) Methods Enzymol. 1985, 60-89
  16. Kunkel, T. A. (1985) Proc. Natl. Acad. Sci. U. S. A. 82, 488-492 [Abstract]
  17. Cogan, U., Kopelman, M., Mokady, S., and Shinitzky, M. (1976) Eur. J. Biochem. 65, 71-78 [Abstract]
  18. Gill, S. C., and Hippel, P. T. (1989) Anal. Biochem. 182, 319-326 [Medline] [Order article via Infotrieve]
  19. Piantini, U., Sorensen, O. W., and Ernst, R. P. (1982) J. Am. Chem. Soc. 104, 6800-6801
  20. Braunschweiler, L., and Ernst, R. P. (1983) J. Magn. Reson. 53, 5521-5528
  21. Bax, A., and Davis, D. V. (1985) J. Magn. Reson. 65, 355-360
  22. Griesinger, C., Otting, G., Wuthrich, K., and Ernst, R. P. (1988) J. Am. Chem. Soc. 110, 7870-7872
  23. Jenneer, J., and Ernst, R. P. (1979) J. Chem. Phys. 71, 4546-4553 [CrossRef]
  24. Macura, S., and Ernst, R. P. (1980) Mol. Phys. 41, 95-117
  25. Tsai, M.-D., and Yan, H. (1991) Biochemistry 30, 6806-6818 [Medline] [Order article via Infotrieve]
  26. Szuts, E. Z., and Harosi, F. I. (1991) Arch. Biochem. Biophys. 287, 297-304 [Medline] [Order article via Infotrieve]
  27. Daly, A. K., and Redfern, C. P. F. (1988) Biochem. Biophys. Acta 965, 118-126 [Medline] [Order article via Infotrieve]
  28. Fiorella, P. D., and Napoli, J. L. (1991) J. Biol. Chem. 266, 16572-16579 [Abstract/Free Full Text]
  29. Ong, D. E., and Chytil, F. (1978) J. Biol. Chem. 253, 4551-4554 [Abstract]
  30. Cowan, S. W., Newcomer, M. E., and Jones, T. A. (1993) J. Mol. Biol. 230, 1225-1246 [CrossRef][Medline] [Order article via Infotrieve]
  31. Jones, T. A., Bergfors, T., Sedzik, J., and Unge, T. (1988) EMBO J. 7, 1597-1604 [Abstract]

©1997 by The American Society for Biochemistry and Molecular Biology, Inc.